Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.     

The regulation of purine ribonucleotide biosynthesis by glucocorticoid hormones

The influence of the glucocorticoid hormones (cortisone, cortisol, corticosterone) on the biosynthesis of purine nucleotides (inosinic acid, guanylic acid and adenylic acid) in different organs was investigated in vivo, by following the incorporation of formate-14C into the acid-soluble nucleotides, after administration of the hormones to adrenalectomized rats. Cortisone and corticosterone show a remarkable and comparable increase of the incorporation of formate-14C only in the purine bases of the liver: cortisol is much more effective, increasing the incorporation of formate-14C into the purine bases even ten times over the basal values. No specific effect is evident either in the kidney or in the heart after glucocorticoid administration. Results are interpreted considering that the action of an individual hormone is specifically restricted to the purine nucleotide synthesis in the liver, and that cortisol seems to be the most efficient from this point of view.     

Deprotonated carboxylic group of amino acids transforms adenine into its rare prototropic tautomers

By UV spectroscopic data for anhydrous DMSO solutions and ab initio HF/6-31G** calculations in vacuum it was shown for the first time that deprotonated amino acid carboxylic group is able to change tautomeric state of a nucleotide base, exactly to convert the N9H ground-state prototropic tautomer of adenine into the N7H and N1H rare ones.     

A theoretical study of the cis-syn pyrimidine dimers in the gas phase and water cluster and a tautomer-bypass mechanism for the origin of UV-induced mutations

A quantum mechanical study of all cis-syn cyclobutane pyrimidine photodimers including the normal and rare tautomeric forms of bases has been performed using the ab initio method at HF/6-31G(d.p), MP2(fc)//HF/6-31G(d,p) and MP2(fc)/6-31G(d,p) levels. A puckering angle of the cyclobutyl ring and twist angle of pyrimidine rings with respect to each other is well described by these calculations. It is predicted that in the gas phase all photodimers containing the rare imino form of cytosine are more stable than those containing its normal form. The Monte Carlo simulations show that the dimer containing the imino form of cytosine is more stabilized by water cluster than that containing its amino forms. The possible biological significance stems from the fact that the cytosine in the dimer directs the incorporation of adenine in the complementary strand during replicative bypass. Data obtained point to the cytosine tautomerism as a possible mechanism for the origin of UV-induced mutation.     

Specific interactions of isoguanine with the neutral and deprotonated carboxylic group of amino acids: results of model quantum chemical calculations

The MNDO/H quantum chemical calculations performed in order to estimate energetic features of the isoguanine (isoGua) prototropic tautomers complexes with acetic acid and its carboxylate-ion (models of neutral and deprotonated forms of amino acid carboxylic group) demonstrate ability of the latter to induce the N9H-->N7H tautomeric transition in the base, being characteristic to other purine bases as well. By contrast, the neutral carboxylic group forms the most stable complex with the ground-state isoGua tautomer N3HN9H.     

外显子与内含子序列分维的双碱基研究

引入碱基间的关联,研究了外显子和内含子序列以双碱基为单位的分维,我们发现在这种情况下,外显子和内显子序列在短程和中程存在自相似性并分别定义了这两个区域的分维。结果表明,短程的分维值Ｄｇ一般比中程的Ｄｍ大,外显子的两个分维值比内含子大。我们改变双联体的位相而分维却不变,这反映出在双联体基础上,外显子的不规则性大于内含子,短程的不规则性大于中程,外显子和内含子序列对以２为周期的结构没有位相的特异性。     

Consequences of one-electron oxidation and one-electron reduction for 4-aminopyrimidine—DFT studies

The consequences of one-electron oxidation and one-electron reduction were studied for 4-aminopyrimidine (4APM), which displays prototropic tautomerism. Since experimental techniques are incapable of detecting less than 0.1% of minor tautomers, quantum-chemical calculations [DFT(B3LYP)/6-311+G(d,p)] were carried out for all possible tautomers of neutral 4AMP and its redox forms, 4APM (+ ?) and 4APM (- ?). Four tautomers were considered: one amine and three imine tautomers (two NH and one CH form). Geometric isomerism of the exo?=?NH group was also taken into account. One-electron oxidation (4APM - e → 4APM (+??)) has no significant effect on the tautomeric preferences; it influences solely the composition of the tautomeric mixture. The amine tautomer is favored for both 4APM (+??) and 4APM. An interesting change in the tautomeric preference occurs for 4APM (- ?). One-electron reduction (4APM?+?e → 4APM (- ?)) favors the C5 atom for the labile proton. The preference of the imine CH tautomer in the tautomeric mixture of 4APM (- ?) may partially explain the origin of CH tautomers in nucleobases.     

Classification of pseudo pairs between nucleotide bases and amino acids by analysis of nucleotide-protein complexes

Nucleotide bases are recognized by amino acid residues in a variety of DNA/RNA binding and nucleotide binding proteins. In this study, a total of 446 crystal structures of nucleotide-protein complexes are analyzed manually and pseudo pairs together with single and bifurcated hydrogen bonds observed between bases and amino acids are classified and annotated. Only 5 of the 20 usual amino acid residues, Asn, Gln, Asp, Glu and Arg, are able to orient in a coplanar fashion in order to form pseudo pairs with nucleotide bases through two hydrogen bonds. The peptide backbone can also form pseudo pairs with nucleotide bases and presents a strong bias for binding to the adenine base. The Watson-Crick side of the nucleotide bases is the major interaction edge participating in such pseudo pairs. Pseudo pairs between the Watson-Crick edge of guanine and Asp are frequently observed. The Hoogsteen edge of the purine bases is a good discriminatory element in recognition of nucleotide bases by protein side chains through the pseudo pairing: the Hoogsteen edge of adenine is recognized by various amino acids while the Hoogsteen edge of guanine is only recognized by Arg. The sugar edge is rarely recognized by either the side-chain or peptide backbone of amino acid residues.     

Fixations of the HIV-1 env gene refute neutralism: New evidence for pan-selective evolution

We examined 103 nucleotide sequences of the HIV-1 env gene, sampled from 35 countries and tested: I) the random (neutral) distribution of the number of nucleotide changes; II) the proportion of bases at molecular equilibrium; III) the neutral expected homogeneity of the distribution of new fixated bases; IV) the hypothesis of the neighbor influence on the mutation rates in a site. The expected random number of fixations per site was estimated by Bose-Einstein statistics, and the expected frequencies of bases by matrices of mutation-fixation rates. The homogeneity of new fixations was analyzed using χ2 and trinomial tests for homogeneity. Fixations of the central base in trinucleotides were used to test the neighbor influence on base substitutions. Neither the number of fixations nor the frequencies of bases fitted the expected neutral distribution. There was a highly significant heterogeneity in the distribution of new fixations, and several sites showed more transversions than transitions, showing that each nucleotide site has its own pattern of change. These three independent results make the neutral theory, the nearly neutral and the neighbor influence hypotheses untenable and indicate that evolution of env is rather highly selective.     

recA protein filaments bind two molecules of single-stranded DNA with off rates regulated by nucleotide cofactor

To probe the role of nucleotide cofactor in the binding of single-stranded DNA to recA protein, we have developed a sedimentation assay using 5'-labeled 32P-poly(dT).recA.poly(dT) complexes sediment quantitatively when centrifuged at 100,000 x g for 45 min, whereas free poly(dT) remains in the supernatant. In the presence of ATP, between 6 and 7 bases cosediment per recA monomer; but when ADP is present or in the absence of added nucleotide cofactor, only 3-3.5 bases/recA monomer cosediment. In competition experiments in which recA.32P-poly(dT) complexes are incubated with unlabeled poly(dT), we again find 3-3.5 bases of labeled poly(dT) cosedimenting per recA monomer when no nucleotide cofactor is present. However, when the same experiment is performed with ATP, only half of the expected 6-7 bases of labeled poly(dT) remain bound to the DNA, demonstrating that half of the poly(dT) in the complex exchanges rapidly with free poly(dT), whereas the other half equilibrates slowly, like poly(dT) in the absence of nucleotide. The rate of exchange of the second more tightly bound poly(dT) is accelerated when ADP is present. Our observations are rationalized by a model in which each recA protein helical filament binds two strands of poly(dT) with a stoichiometry of 3-3.5 bases/recA monomer/strand.     

Solution structure and dynamics of DNA duplexes containing the universal base analogues 5-nitroindole and 5-nitroindole 3-carboxamide

Universal bases hybridize with all other natural DNA or RNA bases, and have applications in PCR and sequencing. We have analysed by nuclear magnetic resonance spectroscopy the structure and dynamics of three DNA oligonucleotides containing the universal base analogues 5-nitroindole and 5-nitroindole-3-carboxamide. In all systems studied, both the 5-nitroindole nucleotide and the opposing nucleotide adopt a standard anti conformation and are fully stacked within the DNA duplex. The 5-nitroindole bases do not base pair with the nucleotide opposite them, but intercalate between this base and an adjacent Watson–Crick pair. In spite of their smooth accommodation within the DNA double-helix, the 5-nitroindole-containing duplexes exist as a dynamic mixture of two different stacking configurations exchanging fast on the chemical shift timescale. These configurations depend on the relative intercalating positions of the universal base and the opposing base, and their exchange implies nucleotide opening motions on the millisecond time range. The structure of these nitroindole-containing duplexes explains the mechanism by which these artificial moieties behave as universal bases.     

Comparison of the nucleotide sequences of a yeast gene family. II. Analysis of spontaneous deletions and insertions

We compared the nucleotide sequences of 3 yeast invertase genes in regions where the homology is better than 90%. In the noncoding region 40 gaps of 1-61 bases were found. This is about half as much as the nucleotide substitutions in the same sequences. We grouped the gaps into 5 categories by their length and the characteristics of their sequences. Group I gaps are about 20 nucleotides long and are flanked by repeated sequence of 6 bases which may trigger the deletion of one of the repeats and the sequence between the repeats. Group II gaps are characterized by a small repeated sequence which is missing in one of the invertase genes. Gaps which occur in sequences exclusively made up of one of the 4 bases are summarized in group III. The 4 gaps in group IV do not show any of these sequence characteristics and they are all just one base long. A 61 nucleotide sequence found in only one of the invertase genes seems to be of complex origin. We conclude that small repeated sequences or monotonous sequences are prone to deletion or insertion mutations.     

Identification of the functional ionic groups of papain by pH/rate profile analysis.

The pH dependence of papain catalysis was analyzed by a scheme which evaluates the kinetic contribution of both protonated and unprotonated species of functional groups involved in catalysis. Kinetic measurements were made at constant pH, without buffers, by automatic titration. The rate-determining step for papain-catalyzed hydrolysis of alpha-N-benzoyl-L-arginine ethyl ester, determined by nucleophile competition, changed from acylation below pH 6.5 to mixed acylation-deacylation above pH 6.5. Kinetic analysis indicated that three prototropic groups governed the pH-specificity of alpha-N-benzoyl-L-arginine ethyl ester hydrolysis. These prototropic groups had pKa values of 4.8, 6.5 to 6.7, and 8.7. Theoretical treatment of the kinetics provided an excellent fit with the experimentally found profile when the contribution of all three prototropic groups was considered. Analysis showed that, in acid, the pathways of papain catalysis were functional with either two or three active-site protons. In base, a single functional ionic pathway is associated with an active site with only one proton. Pathways involving an unprotonated active site are catalytically inoperative in both acid and base. These results indicate that papain exhibits several catalytically functional ionic pathways. The results are discussed in terms of pKa assignments, and the mechanism of papain catalysis.     

A comparative study of the hydrogen-bonding patterns and prototropism in solid 2-thiocytosine (potential antileukemic agent) and cytosine, as studied by 1H-14N NQDR and QTAIM/ DFT

A potential antileukemic and anticancer agent, 2-thiocytosine (2-TC), has been studied experimentally in the solid state by ¹H-¹⁴N NMR-NQR double resonance (NQDR) and theoretically by the quantum theory of atoms in molecules (QTAIM)/density functional theory (DFT). Eighteen resonance frequencies on ¹⁴N were detected at 180 K and assigned to particular nitrogen sites (−NH₂, –N=, and –NH–) in 2-thiocytosine. Factors such as the nonequivalence of molecules (connected to the duplication of sites) and possible prototropic tautomerism (capable of modifying the type of site due to proton transfer) were taken into account during frequency assignment. The result of replacing oxygen with sulfur, which leads to changes in the intermolecular interaction pattern and molecular aggregation, is discussed. This study demonstrates the advantages of combining NQDR and DFT to extract detailed information on the H-bonding properties of crystals with complex H-bonding networks. Solid-state properties were found to have a profound impact on the stabilities and reactivities of both compounds.     

A Theoretical Study of the cis-syn Pyrimidine Dimers in the Gas Phase and Water Cluster and a Tautomer-Bypass Mechanism for the Origin of UV-induced Mutations

Abstract

A quantum mechanical study of all cis-syn cyclobutane pyrimidine photodimers including the normal and rare tautomeric forms of bases has been performed using the ab initio method at HF/6–31G(d,p), MP2(fc)//HF/6–31G(d,p) and MP2(fc)/6–31G(d,p) levels. A puckering angle of the cyclobutyl ring and twist angle of pyrimidine rings with respect to each other is well described by these calculations. It is predicted that in the gas phase all photodimers containing the rare imino form of cytosine are more stable than those containing its normal form. The Monte Carlo simulations show that the dimer containing the imino form of cytosine is more stabilized by water cluster than that containing its amino forms. The possible biological significance stems from the fact that the cytosine in the dimer directs the incorporation of adenine in the complementary strand during replicative bypass. Data obtained point to the cytosine tautomerism as a possible mechanism for the origin of UV-induced mutation.     

Pyrimidine scavenging by Mycobacterium leprae

Mycobacterium leprae incorporated exogenously supplied pyrimidines as bases or nucleosides, but not as a nucleotide, into its nucleic acids. Notably, thymine was incorporated approximately 5 times more rapidly than thymidine by both suspensions of, or intracellular M. leprae. Thymine incorporation was significantly inhibited by clofazamine and dapsone at near-pharmacological levels. Therefore, incorporation of thymine is preferable as an activity for assessing viability of M. leprae. Nucleosides were converted to nucleotides through kinases, bases through phosphoribosyltransferases. Alternatively, thymine and uracil could first be converted to nucleosides. Cytosine and uracil bases were interconvertible, and uracil alone could supply all the pyrimidine requirements of M. leprae, though conversion to the thymine base was extremely slow. Overall, pyrimidine scavenging occurs at a slower rate than, and appears not to be so important as purine scavenging in M. leprae.     

Internucleotide correlations and nucleotide periodicity in Drosophila mtDNA: new evidence for panselective evolution

Analysis for the homogeneity of the distribution of the second base of dinucleotides in relation to the first, whose bases are separated by 0, 1, 2,... 21 nucleotide sites, was performed with the VIH-1 genome (cDNA), the Drosophila mtDNA, the Drosophila Torso gene and the human p-globin gene. These four DNA segments showed highly significant heterogeneities of base distributions that cannot be accounted for by neutral or nearly neutral evolution or by the "neighbor influence" of nucleotides on mutation rates. High correlations are found in the bases of dinucleotides separated by 0, 1 and more number of sites. A periodicity of three consecutive significance values (measured by the x29) was found only in Drosophila mtDNA. This periodicity may be due to an unknown structure or organization of mtDNA. This non-random distribution of the two bases of dinucleotides widespread throughout these DNA segments is rather compatible with panselective evolution and generalized internucleotide co-adaptation.     

Avoidance of base runs in switch regions of immune-system genes

Mouse immunoglobulin (Ig) switch-region sequences are anti-"runny"; that is, they have a smaller amount of their total bases in homonucleotide tracts ("runs") than would be expected if each nucleotide in the sequence were a random selection from a pool of the composition of the region. The switch sequences involve the first intron of rearranged Ig heavy-chain genes; this intron differs strikingly from the succeeding ones, which are "runny" (have more bases than expected in runs). Switch regions are the only category of sequences so far found to be antirunny by statistical test. This sequence characteristic is related to the presence in switch sequences of repeating heteronucleotides. We suggest that the resulting base dispersion and increased complexity favor more specific interactions between sequences, which may be advantageous in recombinational processes such as switching and translocation.      

Mechanisms of targeted frameshift mutations: Insertions arising during error-prone or SOS synthesis of DNA containing cis-syn cyclobutane thymine dimers

It is still unclear how frameshift mutations arise at cyclobutane pyrimidine dimers. The polymerase model is commonly used to explain the mechanisms of various mutations. An alternative polymerase-tautomer model was developed for UV-induced mutagenesis. A mechanism was proposed for targeted insertions caused by cis-syn cyclobutane thymine dimers. Targeted insertions are frameshift mutations due to addition of one or more nucleotides in a DNA sequence opposite to a lesion capable of stopping DNA synthesis. Among other factors, cyclobutane pyrimidine dimers can cause targeted insertions. UV irradiation can change the tautomeric form of DNA bases. Five rare tautomeric forms are possible for thymine, and they are stable when the thymine is a component of a cyclobutane dimer. A structural analysis showed that none of the canonical nucleotides can be added opposite to a specific rare thymine tautomer so that hydrogen bonds form between the two bases. A single nucleotide gap is consequently left in the corresponding site of the nascent strand when a specialized or modified DNA polymerase drives SOS or error-prone DNA synthesis on a template containing cis-syn cyclobutane thymine dimers with a base occurring in the rare tautomeric form. If the DNA composition is homogenous within the region, the end of the growing DNA strand may slip to form a complementary pair with the nucleotide adjacent to the dimer according to the Streisinger model, thus producing a loop. A targeted insertion is thereby generated to make the daughter strand longer. Targeted insertions were for the first time assumed to result from the cis-syn cyclobutane thymine dimers wherein one or both of the bases occur in the specific tautomeric form that does not allow the addition and hydrogen bonding of any canonical nucleotide in the opposite position. A model was developed to explain how targeted insertions of one or more nucleotides are caused by cis-syn cyclobutane thymine dimers. Thus, the polymerase-tautomer model can explain the nature and formation of targeted frameshift mutations in addition to hot and cold spots or targeted or untargeted nucleotide substitutions.