1
|
Dound AS, Fandilolu PM, Sonawane KD. Structural Significance of Conformational Preferences and Ribose-Ring-Puckering of Hyper Modified Nucleotide 5'-Monophosphate 2-Methylthio Cyclic N 6-Threonylcarbamoyladenosine (p-ms 2ct 6A) Present at 37th Position in Anticodon Loop of tRNA Lys. Cell Biochem Biophys 2022; 80:665-680. [PMID: 35965304 DOI: 10.1007/s12013-022-01086-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2021] [Accepted: 07/28/2022] [Indexed: 12/01/2022]
Abstract
Structural significance of conformational preferences and ribose ring puckering of newly discovered hyper modified nucleotide, 5'-monophosphate 2-methylthio cyclic N6-threonylcarbamoyladenosine (p-ms2ct6A) have been investigated using quantum chemical semi-empirical RM1 and molecular dynamics simulation techniques. Automated geometry optimization of most stable structure of p-ms2ct6A has also been carried out with the help of abinitio (HF SCF, DFT) as well as semi empirical quantum chemical (RM1, AM1, PM3, and PM6) methods. Most stable structure of p-ms2ct6A is stabilized by intramolecular interactions between N(3)…HC(2'), N(1)…HC(16), O(13)…HC(15), and O(13)…HO(14). The torsion angles alpha (α) and beta (β) show the significant characteristic patterns with the involvement of intramolecular hydrogen bonding to provide stability to the p-ms2ct6A. Further, molecular dynamics simulations of p-ms2ct6A revealed the role of ribose sugar ring puckering i.e. C2'-endo and C3'-endo on the structural dynamics of ms2ct6A side chain. The modified nucleotide p-ms2ct6A periodically prefers both the C2'-endo and C3'-endo sugar with 'anti' and 'syn' conformations. This property of p-ms2ct6A could be useful to recognize the starting ANN codons. All atom explicit MD simulation of anticodon loop (ACL) of tRNALys of Bacillus subtilis containing ms2ct6A at 37th position showed the U-turn feature, base stacking ability with other adjacent bases and hydrogen bonding interactions similar to the isolated base p-ms2ct6A. The ribose sugar puckering contributes to the orientation of the side chain conformation of p-ms2ct6A. Thus, the present study could be helpful to understand the structure-function relationship of the hypermodified nucleoside, ms2ct6A in recognition of the proper codons AAA/AAG during protein biosynthesis.
Collapse
Affiliation(s)
- Ambika S Dound
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416 004, Maharashtra, India
| | - Prayagraj M Fandilolu
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416 004, Maharashtra, India
| | - Kailas D Sonawane
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416 004, Maharashtra, India.
- Department of Microbiology, Shivaji University, Kolhapur, 416004, Maharashtra, India.
| |
Collapse
|
2
|
Seelam Prabhakar P, Takyi NA, Wetmore SD. Posttranscriptional modifications at the 37th position in the anticodon stem-loop of tRNA: structural insights from MD simulations. RNA (NEW YORK, N.Y.) 2021; 27:202-220. [PMID: 33214333 PMCID: PMC7812866 DOI: 10.1261/rna.078097.120] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Accepted: 11/16/2020] [Indexed: 06/11/2023]
Abstract
Transfer RNA (tRNA) is the most diversely modified RNA. Although the strictly conserved purine position 37 in the anticodon stem-loop undergoes modifications that are phylogenetically distributed, we do not yet fully understand the roles of these modifications. Therefore, molecular dynamics simulations are used to provide molecular-level details for how such modifications impact the structure and function of tRNA. A focus is placed on three hypermodified base families that include the parent i6A, t6A, and yW modifications, as well as derivatives. Our data reveal that the hypermodifications exhibit significant conformational flexibility in tRNA, which can be modulated by additional chemical functionalization. Although the overall structure of the tRNA anticodon stem remains intact regardless of the modification considered, the anticodon loop must rearrange to accommodate the bulky, dynamic hypermodifications, which includes changes in the nucleotide glycosidic and backbone conformations, and enhanced or completely new nucleobase-nucleobase interactions compared to unmodified tRNA or tRNA containing smaller (m1G) modifications at the 37th position. Importantly, the extent of the changes in the anticodon loop is influenced by the addition of small functional groups to parent modifications, implying each substituent can further fine-tune tRNA structure. Although the dominant conformation of the ASL is achieved in different ways for each modification, the molecular features of all modified tRNA drive the ASL domain to adopt the functional open-loop conformation. Importantly, the impact of the hypermodifications is preserved in different sequence contexts. These findings highlight the likely role of regulating mRNA structure and translation.
Collapse
MESH Headings
- Adenosine/analogs & derivatives
- Adenosine/metabolism
- Anticodon/chemistry
- Anticodon/genetics
- Anticodon/metabolism
- Base Pairing
- Base Sequence
- Escherichia coli/genetics
- Escherichia coli/metabolism
- Isopentenyladenosine/chemistry
- Isopentenyladenosine/metabolism
- Molecular Dynamics Simulation
- Nucleic Acid Conformation
- Nucleosides/chemistry
- Nucleosides/metabolism
- RNA Processing, Post-Transcriptional
- RNA, Transfer, Lys/chemistry
- RNA, Transfer, Lys/genetics
- RNA, Transfer, Lys/metabolism
- RNA, Transfer, Phe/chemistry
- RNA, Transfer, Phe/genetics
- RNA, Transfer, Phe/metabolism
Collapse
Affiliation(s)
- Preethi Seelam Prabhakar
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| | - Nathania A Takyi
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| |
Collapse
|
3
|
Fandilolu P, Kamble AS, Dound AS, Sonawane KD. Role of Wybutosine and Mg 2+ Ions in Modulating the Structure and Function of tRNA Phe: A Molecular Dynamics Study. ACS OMEGA 2019; 4:21327-21339. [PMID: 31867527 PMCID: PMC6921629 DOI: 10.1021/acsomega.9b02238] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Accepted: 10/31/2019] [Indexed: 06/10/2023]
Abstract
Transfer RNA remains to be a mysterious molecule of the cell repertoire. With its modified bases and selectivity of codon recognition, it remains to be flexible inside the ribosomal machinery for smooth and hassle-free protein biosynthesis. Structural changes occurring in tRNA due to the presence or absence of wybutosine, with and without Mg2+ ions, have remained a point of interest for structural biologists. Very few studies have come to a conclusion correlating the changes either with the structure and flexibility or with the codon recognition. Considering the above facts, we have implemented molecular modeling methods to address these problems using multiple molecular dynamics (MD) simulations of tRNAPhe along with codons. Our results highlight some of the earlier findings and also shed light on some novel structural and functional aspects. Changes in the stability of tRNAPhe in native or codon-bound states result from the conformations of constituent nucleotides with respect to each other. A smaller change in their conformations leads to structural distortions in the base-pairing geometry and eventually in the ribose-phosphate backbone. MD simulation studies highlight the preference of UUC codons over UUU by tRNAPhe in the presence of wybutosine and Mg2+ ions. This study also suggests that magnesium ions are required by tRNAPhe for proper recognition of UUC/UUU codons during ribosomal interactions with tRNA.
Collapse
Affiliation(s)
- Prayagraj
M. Fandilolu
- Structural
Bioinformatics Unit, Department of Biochemistry and Department of
Microbiology, Shivaji University, Kolhapur 416004, Maharashtra, India
| | - Asmita S. Kamble
- Structural
Bioinformatics Unit, Department of Biochemistry and Department of
Microbiology, Shivaji University, Kolhapur 416004, Maharashtra, India
| | - Ambika S. Dound
- Structural
Bioinformatics Unit, Department of Biochemistry and Department of
Microbiology, Shivaji University, Kolhapur 416004, Maharashtra, India
| | - Kailas D. Sonawane
- Structural
Bioinformatics Unit, Department of Biochemistry and Department of
Microbiology, Shivaji University, Kolhapur 416004, Maharashtra, India
| |
Collapse
|
4
|
Kumbhar NM, Gopal JS. Structural significance of hypermodified nucleoside 5-carboxymethylaminomethyluridine (cmnm 5U) from 'wobble' (34th) position of mitochondrial tRNAs: Molecular modeling and Markov state model studies. J Mol Graph Model 2018; 86:66-83. [PMID: 30336453 DOI: 10.1016/j.jmgm.2018.10.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Revised: 10/02/2018] [Accepted: 10/04/2018] [Indexed: 11/28/2022]
Abstract
A quantum chemical semi-empirical RM1 approach was used to deduce the structural role of hypermodified nucleoside 5-carboxymethylaminomethyluridine 5'-monophosphate (pcmnm5U) from 'wobble' (34th) position of mitochondrial tRNAs. The energetically preferred pcmnm5U(34) adopted a 'skew' conformation for C5-substituted side chain (-CH2-NH2+-CH2-COO-) moiety that orient towards the 5'-ribose-phosphate backbone, which support 'anti' orientation of glycosyl (χ34) torsion angle. Preferred conformation of pcmnm5U(34) was stabilized by O(4) … HC(10), O1P⋯HN(11), O(15) … HN(11), O(15) … HC(10), O4' … HC(6) and O(2) … HC2' hydrogen bonding interactions. The high flexibility of side chain moiety displayed different structural properties for pcmnm5U(34). Three different conformations of pcmnm5U(34) were observed in molecular dynamics simulations and Markov state model studies. The unmodified uracil revealed 'syn' and 'anti' orientations for glycosyl (χ34) torsion angle that substantiate the role of "-CH2-NH2+-CH2-COO-" moiety in maintaining the 'anti' orientation of pcmnm5U(34). The preferred conformation of pcmnm5U(34) helps to recognize Guanosine more proficiently than Adenosine from the third position of codons. The role of pcmnm5U(34) in tRNA biogenesis paves the way to understand its structural significance in usual mitochondrial metabolism and respiration.
Collapse
Affiliation(s)
- Navanath M Kumbhar
- Garware Research Centre, Department of Chemistry, Savitribai Phule Pune University (Formerly University of Pune), Pune, 411007, India.
| | - Janhavi S Gopal
- Institute of Bioinformatics and Biotechnology, Savitribai Phule Pune University (Formerly University of Pune), Pune, 411007, India
| |
Collapse
|
5
|
Fandilolu PM, Kamble AS, Sambhare SB, Sonawane KD. Conformational preferences and structural analysis of hypermodified nucleoside, peroxywybutosine (o2yW) found at 37 th position in anticodon loop of tRNA Phe and its role in modulating UUC codon-anticodon interactions. Gene 2017; 641:310-325. [PMID: 29107006 DOI: 10.1016/j.gene.2017.10.072] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Revised: 09/02/2017] [Accepted: 10/25/2017] [Indexed: 12/17/2022]
Abstract
Hypermodified bases present at 3'-adjacent (37th) position in anticodon loop of tRNAPhe are well known for their contribution in modulating codon-anticodon interactions. Peroxywybutosine (o2yW), a wyosine family member, is one of such tricyclic modified bases observed at the 37th position in tRNAPhe. Conformational preferences and three-dimensional structural analysis of peroxywybutosine have not been investigated in detail at atomic level. Hence, in the present study quantum chemical semi-empirical RM1 and multiple molecular dynamics (MD) simulations have been used to study structural significance of peroxywybutosine in tRNAPhe. Full geometry optimizations over the peroxywybutosine base have also been performed using ab-initio HF-SCF (6-31G**), DFT (B3LYP/6-31G**) and semi-empirical PM6 method to compare the salient properties. RM1 predicted most stable structure shows that the amino-carboxy-propyl side chain of o2yW remains 'distal' to the five membered imidazole ring of tricyclic guanosine. MD simulation trajectory of the isolated peroxy base showed restricted periodical fluctuations of peroxywybutosine side chain which might be helpful to maintain proper anticodon loop structure and mRNA reading frame during protein biosynthesis process. Another comparative MD simulation study of the anticodon stem loop with codon UUC showed various properties, which justify the functional implications of peroxywybutosine at 37th position along with other modified bases present in ASL of tRNAPhe. Thus, this study presents an atomic view into the structural properties of peroxywybutosine, which can be useful to determine its role in the anticodon stem loop in context of codon-anticodon interactions and frame shift mutations.
Collapse
Affiliation(s)
- Prayagraj M Fandilolu
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416004, Maharashtra (M.S.), India
| | - Asmita S Kamble
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416004, Maharashtra (M.S.), India
| | - Susmit B Sambhare
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416004, Maharashtra (M.S.), India
| | - Kailas D Sonawane
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur, 416004, Maharashtra (M.S.), India; Department of Microbiology, Shivaji University, Kolhapur, 416004, Maharashtra (MS), India.
| |
Collapse
|
6
|
Conformational Preferences of Modified Nucleoside 5-Taurinomethyluridine, τm(5)U Occur at 'wobble' 34th Position in the Anticodon Loop of tRNA. Cell Biochem Biophys 2016; 71:1589-603. [PMID: 25388845 DOI: 10.1007/s12013-014-0382-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Abstract
Conformational preferences of hypermodified nucleoside 5-taurinomethyluridine 5'-monophoshate 'p-τm(5)U' (-CH2-NH2(+)-CH2-CH2-SO3(-)) have been investigated using semi-empirical RM1 method. Automated geometry optimization using ab initio molecular orbital HF-SCF (6-31G**) and DFT (B3LYP/6-31G**) calculations have also been made to compare the salient features. The RM1 preferred most stable conformation of 'p-τm(5)U' has been stabilized by hydrogen bonding interactions between O(11a)…HN(8), O1P(34)…HN(8), and O1P(34)…HC(10). Another conformational study of 5-taurinomethyluridine side chain has also been performed in context of anticodon loop bases of E. coli tRNA(Leu). The atom O(11a) of τm(5)U(34) side chain interacts with adenosine (A35) as well as ribose-phosphate backbone which might provide structural stability to the anticodon loop. The glycosyl torsion angle of τm(5)U retains 'anti'-conformation. The solvent accessible surface area calculations revealed the role of τm(5)U in tRNA(Leu) anticodon loop. MD simulation results are found in agreement with RM1 preferred stable structure. The MEPs calculations of τm(5)U(34):G3 model show unique potential tunnels between the hydrogen bond donor and acceptor atoms as compared to τm(5)U(34):A3 model. Thus, these results could pave the way to understand the role of τm(5)U(34) to recognize UUG/UUA codons at atomic level in the mitochondrial disease, MELAS.
Collapse
|
7
|
Sonawane KD, Sambhare SB. The influence of hypermodified nucleosides lysidine and t(6)A to recognize the AUA codon instead of AUG: a molecular dynamics simulation study. Integr Biol (Camb) 2016. [PMID: 26215455 DOI: 10.1039/c5ib00058k] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hypermodified nucleosides lysidine (L) and N(6)-threonylcarbamoyladenosine (t(6)A) influence codon-anticodon interactions during the protein biosynthesis process. Lysidine prevents the misrecognition of the AUG codon as isoleucine and that of AUA as methionine. The structural significance of these modified bases has not been studied in detail at the atomic level. Hence, in the present study we performed multiple molecular dynamics (MD) simulations of anticodon stem loop (ASL) of tRNA(Ile) in the presence and absence of modified bases 'L' and 't(6)A' at the 34th and 37th positions respectively along with trinucleotide 'AUA' and 'AUG' codons. Hydrogen bonding interactions formed by the tautomeric form of lysidine may assist in reading the third base adenine of the 'AUA' codon, unlike the guanine of the 'AUG' codon. Such interactions might be useful to restrict codon specificity to recognize isoleucine tRNA instead of methionine tRNA. The t(6)A side chain interacts with the purine ring of the first codon nucleotide adenine, which might provide base stacking interactions and could be responsible for restricting extended codon-anticodon recognition. We found that ASL tRNA(Ile) in the absence of modifications at the 34th and 37th positions cannot establish proper hydrogen bonding interactions to recognize the isoleucine codon 'AUA' and subsequently disturbs the anticodon loop structure. The binding free energy calculations revealed that tRNA(Ile) ASL with modified nucleosides prefers the codon AUA over AUG. Thus, these findings might be useful to understand the role of modified bases L and t(6)A to recognize the AUA codon instead of AUG.
Collapse
Affiliation(s)
- Kailas D Sonawane
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur-416 004, Maharashtra, India. and Department of Microbiology, Shivaji University, Kolhapur- 416 004, Maharashtra, India
| | - Susmit B Sambhare
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur-416 004, Maharashtra, India.
| |
Collapse
|
8
|
Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach. Cell Biochem Biophys 2016; 74:157-73. [PMID: 27216172 DOI: 10.1007/s12013-016-0731-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2014] [Accepted: 05/01/2016] [Indexed: 12/13/2022]
Abstract
Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.
Collapse
|
9
|
Structural significance of modified nucleoside 5-taurinomethyl-2-thiouridine, τm5s2U, found at ‘wobble’ position in anticodon loop of human mitochondrial tRNALys. Struct Chem 2015. [DOI: 10.1007/s11224-015-0642-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
10
|
Labidi NS. Comparative study of kinetics isomerization of substituted polyacetylene (Cl, F, Br and I): Semi empirical RM1 study. JOURNAL OF SAUDI CHEMICAL SOCIETY 2015. [DOI: 10.1016/j.jscs.2012.01.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
11
|
An investigation of in vitro cytotoxicity and apoptotic potential of aromatic diselenides. Bioorg Med Chem Lett 2014; 24:3440-6. [DOI: 10.1016/j.bmcl.2014.05.075] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2014] [Revised: 05/11/2014] [Accepted: 05/22/2014] [Indexed: 11/18/2022]
|
12
|
Conformational preferences of modified nucleoside N(4)-acetylcytidine, ac4C occur at "wobble" 34th position in the anticodon loop of tRNA. Cell Biochem Biophys 2014; 66:797-816. [PMID: 23408308 DOI: 10.1007/s12013-013-9525-8] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
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.
Collapse
|
13
|
Sambhare SB, Kumbhar BV, Kamble AD, Bavi RS, Kumbhar NM, Sonawane KD. Structural significance of modified nucleosides k2C and t6A present in the anticodon loop of tRNAIle. RSC Adv 2014. [DOI: 10.1039/c3ra47335j] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
14
|
Bavi RS, Sambhare SB, Sonawane KD. MD simulation studies to investigate iso-energetic conformational behaviour of modified nucleosides m(2)G and m(2) 2G present in tRNA. Comput Struct Biotechnol J 2013; 5:e201302015. [PMID: 24688708 PMCID: PMC3962230 DOI: 10.5936/csbj.201302015] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2013] [Revised: 05/16/2013] [Accepted: 05/19/2013] [Indexed: 11/28/2022] Open
Abstract
Modified nucleic acid bases are most commonly found in tRNA. These may contain modifications from simple methylation to addition of bulky groups. Methylation of the four canonical nucleotide bases at a wide variety of positions is particularly prominent among the known modification. Methylation of N2 group of guanine is a relatively common modification in tRNA and rRNA. N2-methylguanosine (m2G) is the second most often encountered nucleoside in E. coli tRNAs. N2, N2- dimethylguanosine (m22G) is found in the majority of eukaryotic tRNAs and involved in forming base pair interactions with adjacent bases. Hence, in order to understand the structural significance of these methylated nucleic acid bases we have carried out molecular dynamics simulation to see the salvation effect. The results obtained shows iso-energetic conformational behaviors for m2G and m22G. The simulation trajectory of m2G shows regular periodical fluctuations suggesting that m2G is equally stable as either s-cis or s-trans rotamers. The two rotamers of m2G may interact canonically or non-canonically with opposite base as s-trans m2G26:C/A/U44 and s-cis m2G26:A/U44. The free rotations around the C-N bond could be the possible reason for these iso-energetic conformations. Dimethylation of G has almost no influence on base pairing with either A or U. Thus, these results reveal that modified nucleosides m2G and m22G may play an important role to prevent tRNA from adopting the unusual mitochondrial like conformation.
Collapse
Affiliation(s)
- Rohit S Bavi
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur 416 004, Maharashtra, India
| | - Susmit B Sambhare
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur 416 004, Maharashtra, India
| | - Kailas D Sonawane
- Structural Bioinformatics Unit, Department of Biochemistry, Shivaji University, Kolhapur 416 004, Maharashtra, India ; Department of Microbiology, Shivaji University, Kolhapur 416 004, Maharashtra, India
| |
Collapse
|
15
|
Kumbhar NM, Kumbhar BV, Sonawane KD. Structural significance of hypermodified nucleic acid base hydroxywybutine (OHyW) which occur at 37th position in the anticodon loop of yeast tRNA(Phe). J Mol Graph Model 2012; 38:174-85. [PMID: 23073221 DOI: 10.1016/j.jmgm.2012.07.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2012] [Revised: 07/18/2012] [Accepted: 07/20/2012] [Indexed: 10/28/2022]
Abstract
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).
Collapse
Affiliation(s)
- Navanath M Kumbhar
- Department of Biotechnology, Shivaji University, Kolhapur 416004, Maharashtra-MS, India
| | | | | |
Collapse
|
16
|
Molecular dynamics simulation and quantum mechanical calculations on α-d-N-acetylneuraminic acid. Carbohydr Res 2012; 351:93-7. [DOI: 10.1016/j.carres.2012.01.015] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2011] [Revised: 01/18/2012] [Accepted: 01/19/2012] [Indexed: 11/19/2022]
|
17
|
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. Cell Biochem Biophys 2012; 61:507-21. [PMID: 21735129 DOI: 10.1007/s12013-011-9233-1] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
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.
Collapse
|