Subject(s)
Betacoronavirus , Coronavirus Infections , Pandemics , Pneumonia, Viral , Autopsy , COVID-19 , Humans , SARS-CoV-2ABSTRACT
ABSTRACT: Measurement of corpse temperature is mainly used for estimation of early postmortem interval, and rectal temperature is often used as a representative of body's core temperature in actual work because it is simple, quick and non-invasive. At present, the rectal temperature postmortem interval estimation method internationally accepted and widely used is HENSSGE's nomogram method, while many domestic scholars also deduced their own regression equations through a large number of case data. Estimation of postmortem interval based on rectal temperature still needs further study. The nomogram method needs to be optimized and extended, and quantification of its influencing factors needs to be dealt with more scientifically. There is still a lack of consensus on the probability and duration of the temperature plateau. There is no clear understanding of the probability and extent of the change in initial temperature caused by various causes. New methods and ideas enrich methodological research, but it still lacks systemicity and practicality. This article reviews the researches on estimation of postmortem interval based on rectal temperature in order to summarize the current situation of previous researches and seek new breakthrough points. Because the decline of body temperature can be easily influenced by many factors in vitro and vivo, and the influencing factors in different regions vary greatly, regionalization research and application may be a practical exploration to improve the accuracy of postmortem interval determination.
Subject(s)
Body Temperature , Postmortem Changes , Temperature , Autopsy , Cadaver , Humans , Probability , Time FactorsSubject(s)
Myelodysplastic Syndromes/blood , Myelodysplastic Syndromes/complications , Platelet Count , Receptors, Fc/therapeutic use , Recombinant Fusion Proteins/therapeutic use , Thrombocytopenia/drug therapy , Thrombocytopenia/etiology , Thrombopoietin/therapeutic use , Humans , Myelodysplastic Syndromes/mortality , Thrombocytopenia/mortality , Treatment OutcomeABSTRACT
The INK4A locus encodes two tumor suppressor genes, p16(INK4A) and p14(ARF), transcribed using alternative exons 1alpha or 1beta spliced onto the same exons 2 and 3. Both p16(INK4A) and p14(ARF) are capable of inhibiting the cell-cycle progression, albeit in different manner; p16(INK4A) is phosphorylation of retinoblastoma (pRB) dependent while p14(ARF) is p53-dependent. In this study, we report the discovery of a novel variant of p16(INK4A), termed p16gamma, in a primary T-cell acute lymphoblastic leukemia (T-ALL) patient sample and a neuroblastoma cell line, which was expressed at both the transcriptional and translational levels. Cloning and sequencing of the p16gamma cDNA revealed that p16gamma was identical to p16(INK4A), except that it contained an in-frame insertion of 197 bp between exons 2 and 3. p16gamma expression was detected in the majority of p16(INK4A)-expressing primary T-ALL and B-ALL patient samples and other p16(INK4A)-expressing tumor samples, but was only barely detectable in some normal mononuclear cells and other non-tumor samples. Structural analysis by nuclear magnetic resonance and circular dichroism confirmed that p16gamma, like p16(INK4A), is also an ankyrin-repeat protein. Functional analysis of p16gamma revealed that p16gamma protein interacted with cyclin D-dependent kinase4 and inhibited its kinase activity. Using a luciferase reporter assay, the transfection of p16gamma repressed the E2F response, the downstream target of pRB, with an efficacy equivalent to that of p16(INK4A). Moreover p16gamma, like p16(INK4A), induced cell-cycle arrest at G(0)/G(1), and inhibited cell growth in colony formation assay.
Subject(s)
Burkitt Lymphoma/metabolism , Cyclin-Dependent Kinase Inhibitor p16/genetics , Cyclin-Dependent Kinase Inhibitor p16/metabolism , G1 Phase , Leukemia-Lymphoma, Adult T-Cell/metabolism , Neuroblastoma/metabolism , Resting Phase, Cell Cycle , Alternative Splicing , Blotting, Western , Burkitt Lymphoma/genetics , Circular Dichroism , Colony-Forming Units Assay , Cyclin-Dependent Kinase 4/metabolism , E2F Transcription Factors/metabolism , Humans , Immunoprecipitation , Leukemia-Lymphoma, Adult T-Cell/genetics , Luciferases/metabolism , Neuroblastoma/genetics , Two-Hybrid System TechniquesABSTRACT
Computational methods that predict three-dimensional structures from amino acid sequences have become increasingly accurate and have provided insights into structure-function relationships for proteins in the absence of structural data. However, the accuracy of computational structural models requires experimental approaches for validation. Here we report direct testing of the predictions of a previously reported structural model of the C-terminus of the human heart Na(+) channel. We focused on understanding the structural basis for the unique effects of an inherited C-terminal mutation (Y1795C), associated with long QT syndrome variant 3 (LQT-3), that has pronounced effects on Na(+) channel inactivation. Here we provide evidence that this mutation, in which a cysteine replaces a tyrosine at position 1795 (Y1795C), enables the formation of disulfide bonds with a partner cysteine in the channel. Using the predictions of the model, we identify the cysteine and show that three-dimensional information contained in the sequence for the channel protein is necessary to understand the structural basis for some of the effects of the mutation. The experimental evidence supports the accuracy of the predicted structural model of the human heart Na(+) channel C-terminal domain and provides insight into a structural basis for some of the mutation-induced altered channel function underlying the disease phenotype.
Subject(s)
Ion Channel Gating/physiology , Models, Cardiovascular , Models, Chemical , Myocytes, Cardiac/chemistry , Myocytes, Cardiac/physiology , Sodium Channels/chemistry , Sodium Channels/physiology , Amino Acid Substitution , Cell Line , Computer Simulation , Humans , Membrane Potentials/physiology , Models, Molecular , Mutagenesis, Site-Directed , Mutation , NAV1.5 Voltage-Gated Sodium Channel , Recombinant Proteins/chemistry , Recombinant Proteins/metabolism , Sequence Analysis, Protein , Structure-Activity RelationshipABSTRACT
Despite recent attention given to the concept of modularity and its potential contribution to the evolvability of organisms, there has been little mention of how such a contribution may affect rates of diversification or how this would be assessed. A first key prediction is that lineages with relatively greater degrees of modularity in given traits should exhibit higher rates of diversification. Four general conditions for testing this prediction of the modular evolvability hypothesis are outlined here. The potential role of modularity as a deterministic factor in adaptive radiations is best examined by looking at historic patterns of diversification rather than just levels of extant diversity, the focus of most analyses of key innovations. Recent developmental evidence supports the notion that phenotypes of juvenile and adult stages of insects with "complete" metamorphosis (Holometabola) are distinct developmental and evolvable modules compared to the highly correlated life stages of insects with "incomplete" metamorphosis (Hemimetabola). Family-level rates of diversification for these two groups were calculated from the fossil record. The Holometabola was found to have a significantly and characteristically higher rate of diversification compared to the less modular Hemimetabola, consistent with the idea that intrinsic differences in modularity can influence the long-term evolvability of organisms. The modular evolvability hypothesis also makes a second key prediction: that characters in more modular clades will exhibit greater levels of variation due to their independence. This provides an independent, phenotypically based test of the hypothesis. We discuss here how this second prediction may be tested in the case of the Hemi- and Holometabola.
Subject(s)
Adaptation, Physiological , Biological Evolution , Insecta/anatomy & histology , Animals , Fossils , Insecta/classification , Insecta/genetics , Insecta/physiology , Phenotype , PhylogenyABSTRACT
We have identified a novel isoform of rat caspase-9 in which the C terminus of full-length caspase-9 is replaced with an alternative peptide sequence. Casp-9-CTD (where CTD is carboxyl-terminal divergent) is expressed in multiple tissues, with the relative highest expression observed in ovary and heart. Casp-9-CTD was found primarily in the cytoplasm and was not detected in the nucleus. Structural predictions suggest that in contrast to full-length caspase-9, casp-9-CTD will not be processed. Our model is supported by reduced protease activity of casp-9-CTD preparations in vitro and by the lack of detectable processing of casp-9-CTD proenzyme or the induction of cell death following transfection into cells. Both neuronal and non-neuronal cell types transfected with casp-9-CTD were resistant to death evoked by trophic factor deprivation or DNA damage. In addition, cytosolic lysates prepared from cells permanently expressing exogenous casp-9-CTD were resistant to caspase induction by cytochrome c in reconstitution assays. Taken together, our observations indicate that casp-9-CTD acts as a dominant-negative variant. Its expression in various tissues indicates a physiological role in regulating cell death.
Subject(s)
Apoptosis/physiology , Caspases/physiology , Amino Acid Sequence , Animals , Base Sequence , Caspase 3 , Caspase 9 , Caspases/chemistry , Caspases/genetics , Caspases/metabolism , Cell Line , Cloning, Molecular , DNA , DNA Primers , Enzyme Activation , Humans , Molecular Sequence Data , PC12 Cells , Rats , Recombinant Proteins/metabolism , Substrate SpecificityABSTRACT
We have devised and implemented in PrISM (protein informatics system for modeling) a new measure of protein structural relationships, the protein structural distance (PSD). The PSD is designed to describe relationships between protein structures in quantitative rather than descriptive terms and is applicable both when two structures are very similar, and when they are very different. It is calculated with a structural alignment procedure that uses double dynamic programming to align secondary structure elements and an iterative rigid body superposition that minimizes the root-mean-square deviation of C(alpha) atoms. The alignment algorithm, as implemented on a modest workstation, is computationally efficient, allowing for large-scale structural comparisons. PSD scores for more than one and a half million pairs of proteins were calculated and compared to the discrete classification of proteins in the SCOP database. The PSD scores, which were obtained automatically, are in large part consistent with the manually derived classifications in SCOP. Discrepancies do arise, however, due, in part, to the fact that SCOP uses criteria other than structural similarity to derive classifications while the PrISM procedure is exclusively structure based. Analysis of PSD scores suggests that there is a continuous aspect of protein conformation space, even though various classification schemes are extremely useful. The use of a continuous measure for structural distance between all pairs of proteins allows us, as described in the two accompanying papers to derive sequence/structure relationships in a more quantitative way than has previously been possible. An important strength of the approach implemented in PrISM is its ability to address many different kinds of queries interactively, making its structural comparison procedure a convenient computational tool that complements structural classification databases such as SCOP and CATH.
Subject(s)
Protein Conformation , Sequence Alignment/methods , Software , Algorithms , Amino Acids/chemistry , Computer Simulation , Databases, Factual , Models, Molecular , Models, Statistical , Protein Folding , Protein Structure, SecondaryABSTRACT
Here, we discuss the relationship between protein sequence and protein structural similarity. It is established that a protein structural distance (PSD) of 2.0 is a threshold above which two proteins are unlikely to have a detectable pairwise sequence relationship. A precise correlation is established between the level of sequence similarity, defined by a normalized Smith-Waterman score, and the probability that two proteins will have a similar structure (defined by pairwise PSD<2). This correlation can be used in evaluating the likelihood for success in a comparative modeling procedure. We establish the existence of a correlation between sequence and structural similarity for pairs of proteins that are related in structure but whose sequence relationship is not detectable using standard pairwise sequence alignments. Although it is well known that there is a close relationship between sequence and structural similarity for pairwise sequence identities greater than about 30 %, there has been little discussion as to the possible existence of such a relationship for pairs of proteins in or below the twilight zone of sequence similarity (<25 % pairwise sequence identity). Possible implications of our results for the evolution of protein structure are discussed.
Subject(s)
Protein Conformation , Sequence Alignment/methods , Software , Algorithms , Amino Acids/chemistry , Computer Simulation , Databases, Factual , Models, Statistical , Protein Folding , Protein Structure, Secondary , Sensitivity and SpecificityABSTRACT
The information required to generate a protein structure is contained in its amino acid sequence, but how three-dimensional information is mapped onto a linear sequence is still incompletely understood. Multiple structure alignments of similar protein structures have been used to investigate conserved sequence features but contradictory results have been obtained, due, in large part, to the absence of subjective criteria to be used in the construction of sequence profiles and in the quantitative comparison of alignment results. Here, we report a new procedure for multiple structure alignment and use it to construct structure-based sequence profiles for similar proteins. The definition of "similar" is based on the structural alignment procedure and on the protein structural distance (PSD) described in paper I of this series, which offers an objective measure for protein structure relationships. Our approach is tested in two well-studied groups of proteins; serine proteases and Ig-like proteins. It is demonstrated that the quality of a sequence profile generated by a multiple structure alignment is quite sensitive to the PSD used as a threshold for the inclusion of proteins in the alignment. Specifically, if the proteins included in the aligned set are too distant in structure from one another, there will be a dilution of information and patterns that are relevant to a subset of the proteins are likely to be lost. In order to understand better how the same three-dimensional information can be encoded in seemingly unrelated sequences, structure-based sequence profiles are constructed for subsets of proteins belonging to nine superfolds. We identify patterns of relatively conserved residues in each subset of proteins. It is demonstrated that the most conserved residues are generally located in the regions where tertiary interactions occur and that are relatively conserved in structure. Nevertheless, the conservation patterns are relatively weak in all cases studied, indicating that structure-determining factors that do not require a particular sequential arrangement of amino acids, such as secondary structure propensities and hydrophobic interactions, are important in encoding protein fold information. In general, we find that similar structures can fold without having a set of highly conserved residue clusters or a well-conserved sequence profile; indeed, in some cases there is no apparent conservation pattern common to structures with the same fold. Thus, when a group of proteins exhibits a common and well-defined sequence pattern, it is more likely that these sequences have a close evolutionary relationship rather than the similarities having arisen from the structural requirements of a given fold.
Subject(s)
Conserved Sequence , Sequence Alignment/methods , Software , Algorithms , Amino Acid Sequence , Amino Acids/chemistry , Computer Simulation , Databases, Factual , Immunoglobulins/chemistry , Models, Molecular , Models, Statistical , Molecular Sequence Data , Protein Folding , Protein Structure, Secondary , Sequence Homology, Amino Acid , Serine Endopeptidases/chemistryABSTRACT
PrISM (Protein Informatics System for Modeling) is a protein analysis and modeling system in which informatics, alignment, modeling, and assessment modules are integrated in a computational environment where protein analysis and modeling protocols can be designed and assessed interactively. It can then be used automatically and repetitively in response to a variety of protein analysis and modeling problems. PrISM was used to predict a single model for each of the 43 targets in the CASP3 experiment. In this paper, we present results for 13 target sequences, which we consider to be comparative modeling targets with clearly related structural templates. We emphasize the problem of aligning a target sequence to a template structure with various alignment methods. When more than one alignment method and/or parameter set are applied, the final alignment is chosen on the basis of a model ranking system also used in PrISM's fold recognition module. Advanced sequence-template alignment procedures in PrISM are useful in some cases when standard pairwise dynamic programming algorithm fail to make any reasonable global alignment. The same procedures, however, failed in other cases, corresponding to remotely related query-template pairs that involved extensive insertions and deletions.
Subject(s)
Models, Molecular , Proteins/chemistry , Algorithms , Amino Acid Sequence , Databases as Topic , Molecular Sequence Data , Sequence AlignmentABSTRACT
Wilson disease (WD) and Menkes disease (MNK) are inherited disorders of copper metabolism. The genes that mutate to give rise to these disorders encode highly homologous copper transporting ATPases. We use yeast and mammalian two-hybrid systems, along with an in vitro assay to demonstrate a specific, copper-dependent interaction between the six metal-binding domains of the WD and MNK ATPases and the cytoplasmic copper chaperone HAH1. We demonstrate that several metal-binding domains interact independently or in combination with HAH1p, although notably domains five and six of WDp do not. Alteration of either the Met or Thr residue of the HAH1p MTCXXC motif has no observable effect on the copper-dependent interaction, whereas alteration of either of the two Cys residues abolishes the interaction. Mutation of any one of the HAH1p C-terminal Lys residues (Lys(56), Lys(57), or Lys(60)) to Gly does not affect the interaction, although deletion of the 15 C-terminal residues abolishes the interaction. We show that apo-HAH1p can bind in vitro to copper-loaded WDp, suggesting reversibility of copper transfer from HAH1p to WD/MNKp. The in vitro HAH1/WDp interaction is metalospecific; HAH1 preincubated with Cu(2+) or Hg(+) but not with Zn(2+), Cd(2+), Co(2+), Ni(3+), Fe(3+), or Cr(3+) interacted with WDp. Finally, we model the protein-protein interaction and present a theoretical representation of the HAH1p.Cu.WD/MNKp complex.
Subject(s)
Adenosine Triphosphatases/metabolism , Carrier Proteins/metabolism , Cation Transport Proteins , Molecular Chaperones , Recombinant Fusion Proteins , Adenosine Triphosphatases/chemistry , Amino Acid Motifs , Carrier Proteins/chemistry , Cell Line , Copper Transport Proteins , Copper-Transporting ATPases , Cytoplasm/metabolism , Humans , Metallochaperones , Models, Molecular , Protein Binding , Two-Hybrid System TechniquesABSTRACT
A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.
Subject(s)
DNA, Viral/metabolism , DNA-Binding Proteins/metabolism , Repressor Proteins/metabolism , Hydrogen-Ion Concentration , Models, Biological , Operator Regions, Genetic/genetics , Static Electricity , Thermodynamics , Viral Proteins , Viral Regulatory and Accessory ProteinsABSTRACT
The effect of unilateral transection of the anterior cruciate ligament on the confined compression and swelling properties of the distal femoral articular cartilage of skeletally mature rabbits at 9 weeks after surgery was determined. Gross morphological grading of the transected and contralateral control distal femora stained with India ink confirmed that cartilage degeneration had been induced by ligament transection. Osteochondral cores, 1.8 mm in diameter, were harvested from the medial femoral condyles. The modulus, permeability, and electrokinetic (streaming potential) coefficient of the articular cartilage of the osteochondral cores were assessed by confined compression creep experiments. The properties (mean +/- SD) of control cartilage were: confined compression modulus, 0.75 +/- 0.28 MPa; hydraulic permeability, 0.63 +/- 0.28 x 10(-15) m2/Pa*sec; and electrokinetic coefficient, 0.16 +/- 0.31 x 10(-9) V/Pa. In transected knees, the modulus was reduced by 18% (p = 0.04), while the permeability and electrokinetic coefficient were not detectably altered. The change in modulus was accompanied by a trend (p = 0.07) toward a decrease (-11%) in the glycosaminoglycan density within the tissue, a significant increase (p < 0.001) in the water content of the cartilage after equilibration in 1 x phosphate buffered saline from 70.3 +/- 4.1% in control knees to 75.2 +/- 4.0% in transected knees, and little further swelling after tissue equilibration in hypotonic saline. The compressive modulus of the cartilage from both control and transected knees was positively correlated with the density of tissue glycosaminoglycan. The alterations in the physical properties of the articular cartilage after transection of the anterior cruciate ligament in the rabbit show trends similar to those observed in human and other animal models of osteoarthritis and provide further support for the use of this model in the study of cartilage degeneration.
Subject(s)
Anterior Cruciate Ligament/surgery , Cartilage, Articular/physiopathology , Animals , Biomechanical Phenomena , Cartilage, Articular/metabolism , Cartilage, Articular/pathology , Glycosaminoglycans/metabolism , Knee Joint/physiopathology , Postoperative Period , RabbitsABSTRACT
The target cytosines of (cytosine-5)-DNA methyltransferases in prokaryotic and eukaryotic DNA show increased rates of C-->T transition mutations compared to non-target cytosines. These mutations are induced either by the spontaneous deamination of 5-mC-->T generating inefficiently repaired G:T rather than G:U mismatches, or by the enzyme-induced C-->U deamination which occurs under conditions of reduced levels of S-adenosylmethionine (AdoMet) and S-adenosylhomocysteine (AdoHcy). We tested whether various inhibitors of (cytosine-5)-DNA methyltransferases analogous to AdoMet and AdoHcy would affect the rate of enzyme-induced deamination of the target cytosine by M.HpaII and M.SssI. Interestingly, we found two compounds, sinefungin and 5'-amino-5'-deoxyadenosine, that increased the rate of deamination 10(3)-fold in the presence and 10(4)-fold in the absence of AdoMet and AdoHcy. We have therefore identified the first mutagenic compounds specific for the target sites of (cytosine-5)-DNA methyltransferases. A number of analogs of AdoMet and AdoHcy have been considered as possible antiviral, anticancer, antifungal and antiparasitic agents. Our findings show that chemotherapeutic agents with affinities to the cofactor binding pocket of (cytosine-5)-DNA methyltransferase should be tested for their potential mutagenic effects.
Subject(s)
DNA-Cytosine Methylases/antagonists & inhibitors , Enzyme Inhibitors/pharmacology , S-Adenosylhomocysteine/analogs & derivatives , S-Adenosylmethionine/analogs & derivatives , Bacteria/enzymology , Binding, Competitive , DNA/metabolism , DNA-Cytosine Methylases/metabolism , Deamination , Dose-Response Relationship, Drug , Methylation , Protein Binding , Structure-Activity RelationshipABSTRACT
The stability of beta-turns is calculated as a function of sequence and turn type with a Monte Carlo sampling technique. The conformational energy of four internal hydrogen-bonded turn types, I, I', II and II', is obtained by evaluating their gas phase energy with the CHARMM force field and accounting for solvation effects with the Finite Difference Poisson-Boltzmann (FDPB) method. All four turn types are found to be less stable than the coil state, independent of the sequence in the turn. The free-energy penalties associated with turn formation vary between 1.6 kcal/mol and 7.7 kcal/mol, depending on the sequence and turn type. Differences in turn stability arise mainly from intraresidue interactions within the two central residues of the turn. For each combination of the two central residues, except for -Gly-Gly-, the most stable beta-turn type is always found to occur most commonly in native proteins. The fact that a model based on local interactions accounts for the observed preference of specific sequences suggests that long-range tertiary interactions tend to play a secondary role in determining turn conformation. In contrast, for beta-hairpins, long-range interactions appear to dominate. Specifically, due to the right-handed twist of beta-strands, type I' turns for -Gly-Gly- are found to occur with high frequency, even when local energetics would dictate otherwise. The fact that any combination of two residues is found able to adopt a relatively low-energy turn structure explains why the amino acid sequence in turns is highly variable. The calculated free-energy cost of turn formation, when combined with related numbers obtained for alpha-helices and beta-sheets, suggests a model for the initiation of protein folding based on metastable fragments of secondary structure.
Subject(s)
Protein Folding , Protein Structure, Secondary , Data Interpretation, Statistical , Databases, Factual , Dipeptides/chemistry , Hydrogen Bonding , Monte Carlo Method , Protein Conformation , Software , ThermodynamicsABSTRACT
Cytosine to thymine transition mutations at the CpG dinucleotide are the most common point mutations in cancer and genetic disease. We calculated the in vivo rate of CpG mutation in the primate germline by deriving a primordial consensus sequence for an Alu repetitive element which inserted into intron 6 of the primate p53 gene 35 to 55 million years ago. Comparison of this primordial sequence to the Alu sequence in intron 6 of present-day primates was used to determine the nature and rate of mutations which occurred during evolution. We estimate the half-life of a CpG nucleotide to be 24 to 60 million years, and the rate constant for mutation at this dinucleotide to be 1.2 x 1O(-8) to 2.9 x 1O(-8) years(-1). These results were confirmed by the analysis of a second Alu sequence in intron 10 of the p53 gene. The in vivo mutation rate is at least 1250-fold slower than the in vitro chemical rate of 5-methylcytosine deamination in double-stranded DNA, showing that current estimates of CpG mutation repair have been significantly underestimated. Furthermore, the mutability of the CpG dinucleotide has led to the depletion of this dinucleotide from the vertebrate genome, and calculations in this study suggest that current levels of the CpG dinucleotide in the primate genome are very close to a steady state equilibrium in which the rate of CpG mutation is equal to the rate of CpG formation by random mutation.
Subject(s)
Dinucleoside Phosphates/genetics , Genes, p53 , Germ-Line Mutation , Primates/genetics , Repetitive Sequences, Nucleic Acid , Animals , Base Sequence , Biological Evolution , Cell Line , DNA , DNA-Cytosine Methylases/metabolism , Humans , Introns , Methylation , Molecular Sequence Data , Polymorphism, Single-Stranded Conformational , Sequence Homology, Nucleic Acid , Tumor Cells, CulturedABSTRACT
C --> T transitions at CpG sites are the most prevalent mutations found in the p53 tumor suppressor gene in human colon tumors and in the germline (Li-Fraumeni syndrome). All of the mutational hot spots are methylated to 5-methylcytosine, and it has been hypothesized that the majority of these mutations are caused by spontaneous hydrolytic deamination of this base to thymine. We have previously reported that bacterial methyltransferases induce transition mutations at CpG sites by increasing the deamination rate of C --> U when the concentration of the methyl group donor S-adenosylmethionine (AdoMet) drops below its Km, suggesting an alternative mechanism to create these mutations. Unrepaired uracil pairs with adenine during replication, completing the C --> T transition mutation. To determine whether this mechanism could contribute to the development of human colon cancer, we examined the level of DNA (cytosine-5)-methyltransferase (MTase) expression, the concentration of AdoMet, and the activity of uracil-DNA glycosylase in human colon tissues, and searched for the presence of mutations in the MTase gene. Using reverse transcription-PCR methods, we found that average MTase mRNA expression levels were only 3.7-fold elevated in tumor tissues compared with surrounding normal mucosa from the same patient. Also, no mutations were found in conserved regions of the gene in 10 tumors sequenced. High-performance liquid chromatographic analysis of extracts from the same tissues showed that AdoMet concentrations were not reduced below the Km value for the mammalian enzyme, and the concentration ratio of AdoMet:S-adenosylhomocysteine, the breakdown product of AdoMet and the competitive MTase inhibitor, did not differ significantly. Finally, extracts from the tumor tissue efficiently removed uracil from DNA. Therefore, biochemical conditions favoring a mutagenic pathway of C --> U --> T were not found in a target tissue known to undergo a high rate of C --> T transitions at CpG sites.