Influence of the hypermodified Y base on the A·U pairing in the codon–anticodon interaction

Structural significance of hypermodified nucleic acid base hydroxywybutine (OHyW) which occur at 37th position in the anticodon loop of yeast tRNA(Phe)

Conformational preferences of hypermodified nucleic acid base hydroxywybutine (OHyW) have been studied using quantum chemical single point semi-empirical PM3 method. Automated geometry optimization using semi-empirical RM1, molecular mechanics force field (MMFF) along with ab-initio HF-SCF (6-31G** basis set) and DFT (B3LYP/6-31G** basis set) calculations have also been made to compare the salient features. Molecular electrostatic potentials (MEPs) depict the polarities of hydroxywybutine (OHyW) side chain. Another conformational study showed that hydroxywybutosine side chain interacts with adjacent bases within the anticodon loop of tRNA(Phe). The solvent accessible surface area (SASA) calculations revealed the structural role of hydroxywybutine in anticodon loop. Explicit molecular dynamics (MD) simulation has been done over the PM3 most stable structure of OHyW. The hydroxywybutine side chain prefers 'distal' conformation i.e. spreads away from the cyclic five membered imidazole moiety of modified tricyclic guanine base. The predicted preferred conformation of hydroxywybutine may prevent extended Watson-Crick base pairing during protein biosynthesis process. This conformation of OHyW stabilized by intramolecular interactions between O(6)⋯HO(16), O(6)⋯HC(15) and O(20)⋯HC(17). Further stabilization is also expected from interactions between O(22)⋯HC(16) and O(23)⋯HC(15). Explicit molecular dynamics (MD) simulation over the PM3 most stable structure of OHyW support the preferred geometry by preserving the 'distal' orientation of hydroxywybutine side chain and intramolecular hydrogen bonding interactions. MD simulation study revealed the role of hydroxyl group of OHyW to avoid fluctuations and prevent multiple iso-energetic conformations of hydroxywybutine side chain as compared to wybutine (yW).

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.

Energy and structural analysis of double nucleic acid triplets

In order to investigate the energy and structural character of RNA-RNA triplets and RNA-DNA duplex base triplets, 64 sets of three-dimensional models of RNA-DNA duplex base triplets and mRNA-tRNA triplex base triplets were constructed and optimized by homologous modeling method using the software InsightII. The comparative statistical method and cluster analysis were adopted to study these features. The result showed: (i) all energy parameters of monomer RNA-DNA hybrid triplets and ternary complexes appeared significantly different; and some parameters related with overall molecules such as overall energy, bond energy and coulomb energy have statistically significant correlations between the structures in vacuum and aquatic solutions while other parameters, including theta energy, phi energy, hydrogen bond energy and non-bond energy, changed significantly, but not continuously. (ii) However, the case of mRNA-tRNA triplets was much more complicated in that only the bond energy's correlation coefficient is -0.8. Typically, the main contribution of GC pairs and G/A/U bases were interesting. The models of RNA-DNA hybrid triplets and mRNA-tRNA triplet should be helpful for the study of base pairing in codons and the biological effectiveness of antisense nucleic acids.

Molecular dynamics of the anticodon domain of yeast tRNA(Phe): codon-anticodon interaction

We have studied the effect of codon-anticodon interaction on the structure and dynamics of transfer RNAs using molecular dynamics simulations over a nanosecond time scale. From our molecular dynamical investigations of the solvated anticodon domain of yeast tRNA(Phe) in the presence and absence of the codon trinucleotides UUC and UUU, we find that, although at a gross level the structures are quite similar for the free and the bound domains, there are small but distinct differences in certain parts of the molecule, notably near the Y37 base. Comparison of the dynamics in terms of interatomic or inter-residual distance fluctuation for the free and the bound domains showed regions of enhanced rigidity in the loop region in the presence of codons. Because fluorescence experiments suggested the existence of multiple conformers of the anticodon domain, which interconvert on a much larger time scale than our simulations, we probed the conformational space using five independent trajectories of 500 ps duration. A generalized ergodic measure analysis of the trajectories revealed that at least for this time scale, all the trajectories populated separate parts of the conformational space, indicating a need for even longer simulations or enhanced sampling of the conformational space to give an unequivocal answer to this question.