Conformational preferences of modified nucleic acid bases N6-methyl-N6-(N-threonylcarbonyl) adenine and 2-methylthio-N6-(N-threonylcarbonyl) adenine by the quantum chemical PCILO calculations

Conformational preferences of modified nucleoside N(4)-acetylcytidine, ac4C occur at "wobble" 34th position in the anticodon loop of tRNA

Conformational preferences of modified nucleoside, N(4)-acetylcytidine, ac(4)C have been investigated using quantum chemical semi-empirical RM1 method. Automated geometry optimization using PM3 method along with ab initio methods HF SCF (6-31G**), and density functional theory (DFT; B3LYP/6-31G**) have also been made to compare the salient features. The most stable conformation of N(4)-acetyl group of ac(4)C prefers "proximal" orientation. This conformation is stabilized by intramolecular hydrogen bonding between O(7)···HC(5), O(2)···HC2', and O4'···HC(6). The "proximal" conformation of N(4)-acetyl group has also been observed in another conformational study of anticodon loop of E. coli elongator tRNA(Met). The solvent accessible surface area (SASA) calculations revealed the role of ac(4)C in anticodon loop. The explicit molecular dynamics simulation study also shows the "proximal" orientation of N(4)-acetyl group. The predicted "proximal" conformation would allow ac(4)C to interact with third base of codon AUG/AUA whereas the 'distal' orientation of N(4)-acetyl cytidine side-chain prevents such interactions. Single point energy calculation studies of various models of anticodon-codon bases revealed that the models ac(4)C(34)(Proximal):G3, and ac(4)C(34)(Proximal):A3 are energetically more stable as compared to models ac(4)C(34)(Distal):G3, and ac(4)C(34)(Distal):A3, respectively. MEPs calculations showed the unique potential tunnels between the hydrogen bond donor-acceptor atoms of ac(4)C(34)(Proximal):G3/A3 base pairs suggesting role of ac(4)C in recognition of third letter of codons AUG/AUA. The "distal" conformation of ac(4)C might prevent misreading of AUA codon. Hence, this study could be useful to understand the role of ac(4)C in the tertiary structure folding of tRNA as well as in the proper recognition of codons during protein biosynthesis process.

Conformational preferences of modified nucleoside N(2)-methylguanosine (m(2)G) and its derivative N(2), N(2)-dimethylguanosine (m(2)(2)G) occur at 26th position (hinge region) in tRNA

Conformational preferences of the modified nucleosides N(2)-methylguanosine (m(2)G) and N(2), N(2)-dimethylguanosine (m(2)(2)G) have been studied theoretically by using quantum chemical perturbative configuration interaction with localized orbitals (PCILO) method. Automated complete geometry optimization using semiempirical quantum chemical RM1, along with ab initio molecular orbital Hartree-Fock (HF-SCF), and density functional theory (DFT) calculations has also been made to compare the salient features. Single-point energy calculation studies have been made on various models of m(2)G26:C/A/U44 and m(2)(2)G26:C/A/U44. The glycosyl torsion angle prefers "syn" (χ = 286°) conformation for m(2)G and m(2)(2)G molecules. These conformations are stabilized by N(3)-HC2' and N(3)-HC3' by replacing weak interaction between O5'-HC(8). The N(2)-methyl substituent of (m(2)G26) prefers "proximal" or s-trans conformation. It may also prefer "distal" or s-cis conformation that allows base pairing with A/U44 instead of C at the hinge region. Thus, N(2)-methyl group of m(2)G may have energetically two stable s-trans m(2)G:C/A/U or s-cis m(2)G:A/U rotamers. This could be because of free rotations around C-N bond. Similarly, N(2), N(2)-dimethyl substituent of (m(2)(2)G) prefers "distal" conformation that may allow base pairing with A/U instead of C at 44th position. Such orientations of m(2)G and m(2)(2)G could play an important role in base-stacking interactions at the hinge region of tRNA during protein biosynthesis process.

Conformation effects of base modification on the anticodon stem-loop of Bacillus subtilis tRNA(Tyr)

tRNA molecules contain 93 chemically unique nucleotide base modifications that expand the chemical and biophysical diversity of RNA and contribute to the overall fitness of the cell. Nucleotide modifications of tRNA confer fidelity and efficiency to translation and are important in tRNA-dependent RNA-mediated regulatory processes. The three-dimensional structure of the anticodon is crucial to tRNA-mRNA specificity, and the diverse modifications of nucleotide bases in the anticodon region modulate this specificity. We have determined the solution structures and thermodynamic properties of Bacillus subtilis tRNA(Tyr) anticodon arms containing the natural base modifications N(6)-dimethylallyl adenine (i(6)A(37)) and pseudouridine (ψ(39)). UV melting and differential scanning calorimetry indicate that the modifications stabilize the stem and may enhance base stacking in the loop. The i(6)A(37) modification disrupts the hydrogen bond network of the unmodified anticodon loop including a C(32)-A(38)(+) base pair and an A(37)-U(33) base-base interaction. Although the i(6)A(37) modification increases the dynamic nature of the loop nucleotides, metal ion coordination reestablishes conformational homogeneity. Interestingly, the i(6)A(37) modification and Mg(2+) are sufficient to promote the U-turn fold of the anticodon loop of Escherichia coli tRNA(Phe), but these elements do not result in this signature feature of the anticodon loop in tRNA(Tyr).

Iso-energetic multiple conformations of hypermodified nucleic acid base wybutine (yW) which occur at 37(th) position in anticodon loop of tRNA(Phe)

Conformational preferences of wybutine (yW) have been studied by quantum chemical semi-empirical Perturbative Configuration Interaction with Localized Orbitals (PCILO) and PM3 methods. Automated full geometry optimization by using RM1 along with ab-initio Hartree-Fock (HF-SCF) and Density Functional Theory (B3LYP/6-31G**) calculations have also been made to compare the salient features. Molecular dynamics (MD) simulation has been performed to see the solvation effect on wybutine side chain. The preferred conformations of wybutine side chain spreads away 'distal' from five membered imidazole moiety of tricyclic base. The intramolecular interactions provide stability to the preferred wybutine structure. The most stable and alternative stable structures obtained by PCILO and PM3 methods reveal that wybutine side chain may have multiple iso-energetic conformations. Molecular dynamics (MD) simulation study also confirms multiple conformations of wybutine side chain by showing regular periodical fluctuations over the 2 ns time period. These fluctuations occur when torsion angle α takes value ±90° and ±120° as observed in the most stable and alternative stable structures resulted by PCILO and PM3 methods. Such conformational behavior of wybutine may have certain implications on frameshifting to prevent extended Watson-Crick base pairing by maintaining proper codon-anticodon interactions during protein biosynthesis process.

Conformational preferences of the base substituent in hypermodified nucleotide queuosine 5'-monophosphate 'pQ' and protonated variant 'pQH+'

Conformational preferences of the base substituent in hypermodified nucleotide queuosine 5'-monophosphate 'pQ' and its protonated form 'pQH+' have been studied using quantum chemical Perturbative Configuration Interaction with Localized Orbitals PCILO method. The salient points have also been examined using molecular mechanics force field MMFF, parameterized modified neglect of differential overlap PM3 and Hartree Fock-Density Functional Theory HF DFT (pBP/DN*) approaches. Aqueous solvation of pQ and pQH+ has also been studied using molecular dynamics simulations. Consistent with the observed crystal structure, in isolated protonated form pQH+, the quaternary amine HN(13)(+)H, of the sidechain having 7-aminomethyl linkage, hydrogen bonds with the carbonyl oxygen O(10) of the base. However, N(13)H-O(10) hydrogen bonding is not preferred for unprotonated pQ, whether isolated or hydrated. Interaction between the 5'-phosphate and the 7-aminomethyl group is more likely for isolated pQ. The cyclopentenediol hydroxyl group O4"H may hydrogen bond with the O(10) in isolated pQ as well as in pQH+. The O4"H may hydrogen bond with the 5'-phosphate as well. The presence of -CH2-NH- and O"H groups in pQ and pQH+ allows interesting possibilities for intranucleotide hydrogen bonds and interactions across the anticodon loop. Simultaneous hydrogen bonds O2P-HN(13)+H-O(10) are indicated for hydrated pQH+. Unlike weak involvement of O4"H, these interactions also persist in hydrated pQH+ and may much reduce backbone flexibility. Resulting sub-optimal Q:C base pairing leads to unbiased reading of U or C as the third codon letter. Cyclopentenediol hydroxyl groups may interact with other biomolecules, allowing specific recognition. Prospective pQ(34) and pQ(34)H+ sites for codon-anticodon base pairing remain unhindered, but non canonical Q:G base pairing (amber-suppression) is ruled out.

Synthesis of the tRNA(Lys,3) anticodon stem-loop domain containing the hypermodified ms2t6A nucleoside

The synthesis of a protected form of the hypermodified nucleoside, N-[(9-beta-D-ribofuranosyl-2-methylthiopurin-6-yl)carbamoyl]threonine, (ms2t6A) is reported. The hypermodified nucleoside was subsequently elaborated to the protected nucleoside phosphophoramidite using a protecting group strategy compatible with standard RNA oligonucleotide chemistry. The phosphoramidite reagent was then used to synthesize the 17-nucleotide RNA hairpin having the sequence of the anticodon stem-loop (ASL) domain of human tRNA(Lys,3), the primer for HIV-1 reverse transcriptase. Introduction of the modification at position 37 of the tRNA ASL modestly decreases the thermodynamic stability of the RNA hairpin as has been seen previously for the prokaryotic t6A nucleoside lacking the 2-methylthio substituent. 2D NOESY NMR spectra of the ms2t6A containing tRNA ASL indicate that the threonyl side chain adopts a conformation similar to that seen in the solution structure of the analogous t6A containing E. coli tRNA(Lys), despite the presence of the bulky methylthio group. This synthetic approach allows for site-specific incorporation of the hypermodified nucleoside and should facilitate future studies directed at understanding the roles of nucleoside modification in modulating the stability and specificity of biologically important RNA-RNA interactions. Our synthesis of the ms2t6A containing RNAs demonstrates that this methodology is suitable for obtaining quantities of RNA required for structural studies of the HIV primer tRNA.

Conformational preferences of anticodon 3'-adjacent hypermodified nucleic acid base cis-or trans-zeatin and its 2-methylthio derivative, cis-or trans- ms(2)zeatin

Conformational preferences of the hypermodified nucleic acid bases N6-(Delta(2)-cis-hydroxyisopentenyl)adenine, cis-io(6)Ade also known as cis-zeatin, and N(6)-(Delta(2)-trans-hydroxyisopentenyl)adenine, trans-io(6)ade or trans-zeatin, and 2-methylthio derivatives of these cis-ms(2)io(6)Ade or cis-ms(2)zeatin, and trans-ms(2)io6Ade or trans-ms(2)zeatin have been investigated theoretically by the quantum chemical Perturbative Configuration Interaction with Localized Orbitals (PCILO) method. Automated geometry optimization using quantum chemical MNDO, AM1 and PM3 methods has also been made to compare the salient features. The predicted most stable conformation of cis-io(6)Ade, trans-io(6)Ade, cis-ms(2)io(6)Ade and trans-ms(2)io(6)Ade are such that in each of these molecules the isopentenyl substituent spreads away (has "dista" conformation) from the five membered ring imidazole moiety of the adenine. The atoms N(6), C(10) and C(11) remain coplanar with the adenine ring in the predicted preferred conformation for each of these molecules. In cis-io(6)Ade as well as cis-ms(2)io(6)Ade the hydroxyl oxygen may participate in intramolecular hydrogen bonding with the H-C(10)-H group. In trans-io(6)Ade the hydroxyl group is oriented towards the H-C(2) instead. This orientation is retained in trans-ms(2)io(6)Ade, possible O-H...S hydrogen bonding may be a stabilizing factor. In all these four modified adenines C(11)-H is favourably placed to participate in intramolecular hydrogen bonding with N(1). In cis-ms(2)io(6)Ade as well as trans-ms(2)io(6)Ade the 2-methylthio group preferentially orients on the same side as C(2)-N(3) bond, due to this non-obstrusive placing, orientation of the hydroxyisopentenyl substituent remains unaffected by 2-methylthiolation. Thus the N(1) site remains shielded irrespective of the 2-methylthiolation status in these various cis-and trans-zeatin analogs alike. Firmly held orientation of hydroxyisopentenyl substituent in zeatin isomers and derivatives, in contrast to adaptable orientation of isopentenyl substituent in i(6)Ade and ms(2)i(6)Ade, may account for the increased efficiency of suppressor tRNA and reduced codon context sensitivity accompanied with the occurrence of ms(2)-zeatin (ms(2)io(6)Ade) modification.