1
|
Durgam L, Guruprasad L. Molecular mechanism of ATP and RNA binding to Zika virus NS3 helicase and identification of repurposed drugs using molecular dynamics simulations. J Biomol Struct Dyn 2022; 40:12642-12659. [PMID: 34516356 DOI: 10.1080/07391102.2021.1973909] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Congenital Zika virus syndrome has caused a public health emergency of international concern. So far, there are no drugs available to prevent or treat the infection caused by Zika virus. The Zika virus NS3 helicase is a potential protein target for drug discovery due to its vital role in viral genome replication. NS3 helicase unwinds the viral RNA to enable the reproduction of the viral genome by the NS5 protein. NS3 helicase has two crucial binding sites; the ATP binding site and the RNA binding site. We used molecular docking and molecular dynamics (MD) simulations to study the structural behavior of Zika virus NS3 helicase in its apo form and in the presence of ATP, single-stranded RNA, and both ATP-RNA to understand their potential implications in NS3 helicase activity. Further, we have carried out virtual screening of FDA approved drugs, followed by molecular docking to identify the ATP-competitive hit molecules as probable Zika virus NS3 helicase inhibitors. The MD simulations trajectories were analyzed using normal mode analysis and principal component analysis that reveals fluctuations in the R-loop. These findings aid in understanding the molecular mechanisms of the simultaneous binding of ATP and RNA, and guide the design and discovery of new inhibitors of the Zika virus NS3 helicase as a promising drug target to treat the Zika virus infection. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Laxman Durgam
- School of Chemistry, University of Hyderabad, Hyderabad, India
| | | |
Collapse
|
2
|
Nishikawa T, Wojciak JM, Dyson HJ, Wright PE. RNA Binding by the KTS Splice Variants of Wilms' Tumor Suppressor Protein WT1. Biochemistry 2020; 59:3889-3901. [PMID: 32955251 DOI: 10.1021/acs.biochem.0c00602] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Wilms' tumor suppressor protein WT1 regulates the expression of multiple genes through binding of the Cys2-His2 zinc finger domain to promoter sites. WT1 has also been proposed to be involved in post-transcriptional regulation, by binding to RNA using the same set of zinc fingers. WT1 has two major splice variants, where the Lys-Thr-Ser (KTS) tripeptide is inserted into the linker between the third and fourth zinc fingers. To obtain insights into the mechanism by which the different WT1 splice variants recognize both DNA and RNA, we have determined the solution structure of the WT1 (-KTS) zinc finger domain in complex with a 29mer stem-loop RNA. Zinc fingers 1-3 bind in a widened major groove favored by the presence of a bulge nucleotide in the double-stranded helical stem. Fingers 2 and 3 make specific contacts with the nucleobases in a conserved AUGG sequence in the helical stem. Nuclear magnetic resonance chemical shift mapping and relaxation analysis show that fingers 1-3 of the two splice variants (-KTS and +KTS) of WT1 form similar complexes with RNA. Finger 4 of the -KTS isoform interacts weakly with the RNA loop, an interaction that is abrogated in the +KTS isoform, and both isoforms bind with similar affinity to the RNA. In contrast, finger 4 is required for high-affinity binding to DNA and insertion of KTS into the linker of fingers 3 and 4 abrogates DNA binding. While finger 1 is required for RNA binding, it is dispensable for binding to consensus DNA sites.
Collapse
Affiliation(s)
- Tadateru Nishikawa
- Department of Integrative Structural and Computational Biology and Skaggs Institute of Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Jonathan M Wojciak
- Department of Integrative Structural and Computational Biology and Skaggs Institute of Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - H Jane Dyson
- Department of Integrative Structural and Computational Biology and Skaggs Institute of Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Peter E Wright
- Department of Integrative Structural and Computational Biology and Skaggs Institute of Chemical Biology, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, California 92037, United States
| |
Collapse
|
3
|
Ren H, Cheyne CG, Fleming AM, Burrows CJ, White HS. Single-Molecule Titration in a Protein Nanoreactor Reveals the Protonation/Deprotonation Mechanism of a C:C Mismatch in DNA. J Am Chem Soc 2018; 140:5153-5160. [PMID: 29562130 DOI: 10.1021/jacs.8b00593] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Measurement of single-molecule reactions can elucidate microscopic mechanisms that are often hidden from ensemble analysis. Herein, we report the acid-base titration of a single DNA duplex confined within the wild-type α-hemolysin (α-HL) nanopore for up to 3 h, while monitoring the ionic current through the nanopore. Modulation between two states in the current-time trace for duplexes containing the C:C mismatch in proximity to the latch constriction of α-HL is attributed to the base flipping of the C:C mismatch. As the pH is lowered, the rate for the C:C mismatch to flip from the intra-helical state to the extra-helical state ( kintra-extra) decreases, while the rate for base flipping from the extra-helical state to the intra-helical state ( kextra-intra) remains unchanged. Both kintra-extra and kextra-intra are on the order of 1 × 10-2 s-1 to 1 × 10-1 s-1 and remain stable over the time scale of the measurement (several hours). Analysis of the pH-dependent kinetics of base flipping using a hidden Markov kinetic model demonstrates that protonation/deprotonation occurs while the base pair is in the intra-helical state. We also demonstrate that the rate of protonation is limited by transport of H+ into the α-HL nanopore. Single-molecule kinetic isotope experiments exhibit a large kinetic isotope effect (KIE) for kintra-extra ( kH/ kD ≈ 5) but a limited KIE for kextra-intra ( kH/ kD ≈ 1.3), supporting our model. Our experiments correspond to the longest single-molecule measurements performed using a nanopore, and demonstrate its application in interrogating mechanisms of single-molecule reactions in confined geometries.
Collapse
Affiliation(s)
- Hang Ren
- Department of Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| | - Cameron G Cheyne
- Department of Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| | - Aaron M Fleming
- Department of Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| | - Cynthia J Burrows
- Department of Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| | - Henry S White
- Department of Chemistry , University of Utah , 315 South 1400 East , Salt Lake City , Utah 84112 , United States
| |
Collapse
|
4
|
The human mitochondrial tRNAMet: structure/function relationship of a unique modification in the decoding of unconventional codons. J Mol Biol 2010; 406:257-74. [PMID: 21168417 DOI: 10.1016/j.jmb.2010.11.042] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2010] [Revised: 11/18/2010] [Accepted: 11/22/2010] [Indexed: 11/21/2022]
Abstract
Human mitochondrial mRNAs utilize the universal AUG and the unconventional isoleucine AUA codons for methionine. In contrast to translation in the cytoplasm, human mitochondria use one tRNA, hmtRNA(Met)(CAU), to read AUG and AUA codons at both the peptidyl- (P-), and aminoacyl- (A-) sites of the ribosome. The hmtRNA(Met)(CAU) has a unique post-transcriptional modification, 5-formylcytidine, at the wobble position 34 (f(5)C(34)), and a cytidine substituting for the invariant uridine at position 33 of the canonical U-turn in tRNAs. The structure of the tRNA anticodon stem and loop domain (hmtASL(Met)(CAU)), determined by NMR restrained molecular modeling, revealed how the f(5)C(34) modification facilitates the decoding of AUA at the P- and the A-sites. The f(5)C(34) defined a reduced conformational space for the nucleoside, in what appears to have restricted the conformational dynamics of the anticodon bases of the modified hmtASL(Met)(CAU). The hmtASL(Met)(CAU) exhibited a C-turn conformation that has some characteristics of the U-turn motif. Codon binding studies with both Escherichia coli and bovine mitochondrial ribosomes revealed that the f(5)C(34) facilitates AUA binding in the A-site and suggested that the modification favorably alters the ASL binding kinetics. Mitochondrial translation by many organisms, including humans, sometimes initiates with the universal isoleucine codons AUU and AUC. The f(5)C(34) enabled P-site codon binding to these normally isoleucine codons. Thus, the physicochemical properties of this one modification, f(5)C(34), expand codon recognition from the traditional AUG to the non-traditional, synonymous codons AUU and AUC as well as AUA, in the reassignment of universal codons in the mitochondria.
Collapse
|
5
|
Sarzyńska J, Réblová K, Šponer J, Kuliński T. Conformational transitions of flanking purines in HIV-1 RNA dimerization initiation site kissing complexes studied by CHARMM explicit solvent molecular dynamics. Biopolymers 2008; 89:732-46. [DOI: 10.1002/bip.21001] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
6
|
Niemann M, Brecht M, Schlüter E, Weitzel K, Zacharias M, Göringer HU. TbMP42 is a structure-sensitive ribonuclease that likely follows a metal ion catalysis mechanism. Nucleic Acids Res 2008; 36:4465-73. [PMID: 18603593 PMCID: PMC2490751 DOI: 10.1093/nar/gkn410] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2008] [Revised: 06/10/2008] [Accepted: 06/10/2008] [Indexed: 12/22/2022] Open
Abstract
RNA editing in African trypanosomes is characterized by a uridylate-specific insertion and/or deletion reaction that generates functional mitochondrial transcripts. The process is catalyzed by a multi-enzyme complex, the editosome, which consists of approximately 20 proteins. While for some of the polypeptides a contribution to the editing reaction can be deduced from their domain structure, the involvement of other proteins remains elusive. TbMP42, is a component of the editosome that is characterized by two C(2)H(2)-type zinc-finger domains and a putative oligosaccharide/oligonucleotide-binding fold. Recombinant TbMP42 has been shown to possess endo/exoribonuclease activity in vitro; however, the protein lacks canonical nuclease motifs. Using a set of synthetic gRNA/pre-mRNA substrate RNAs, we demonstrate that TbMP42 acts as a topology-dependent ribonuclease that is sensitive to base stacking. We further show that the chelation of Zn(2+) cations is inhibitory to the enzyme activity and that the chemical modification of amino acids known to coordinate Zn(2+) inactivates rTbMP42. Together, the data are suggestive of a Zn(2+)-dependent metal ion catalysis mechanism for the ribonucleolytic activity of rTbMP42.
Collapse
Affiliation(s)
- Moritz Niemann
- Genetics, Darmstadt University of Technology, Schnittspahnstraße 10, 64287 Darmstadt and Computational Biology, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Michael Brecht
- Genetics, Darmstadt University of Technology, Schnittspahnstraße 10, 64287 Darmstadt and Computational Biology, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Elke Schlüter
- Genetics, Darmstadt University of Technology, Schnittspahnstraße 10, 64287 Darmstadt and Computational Biology, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Kerstin Weitzel
- Genetics, Darmstadt University of Technology, Schnittspahnstraße 10, 64287 Darmstadt and Computational Biology, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Martin Zacharias
- Genetics, Darmstadt University of Technology, Schnittspahnstraße 10, 64287 Darmstadt and Computational Biology, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - H. Ulrich Göringer
- Genetics, Darmstadt University of Technology, Schnittspahnstraße 10, 64287 Darmstadt and Computational Biology, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
7
|
Popenda L, Adamiak RW, Gdaniec Z. Bulged Adenosine Influence on the RNA Duplex Conformation in Solution. Biochemistry 2008; 47:5059-67. [DOI: 10.1021/bi7024904] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Lukasz Popenda
- Institute of Bioorganic Chemistry, Polish Academy of Science, Noskowskiego 12/14, 61-704 Poznań, Poland
| | - Ryszard W. Adamiak
- Institute of Bioorganic Chemistry, Polish Academy of Science, Noskowskiego 12/14, 61-704 Poznań, Poland
| | - Zofia Gdaniec
- Institute of Bioorganic Chemistry, Polish Academy of Science, Noskowskiego 12/14, 61-704 Poznań, Poland
| |
Collapse
|
8
|
McCrate NE, Varner ME, Kim KI, Nagan MC. Molecular dynamics simulations of human tRNA Lys,3 UUU: the role of modified bases in mRNA recognition. Nucleic Acids Res 2006; 34:5361-8. [PMID: 17012271 PMCID: PMC1636460 DOI: 10.1093/nar/gkl580] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accuracy in translation of the genetic code into proteins depends upon correct tRNA–mRNA recognition in the context of the ribosome. In human tRNAUUULys,3 three modified bases are present in the anticodon stem–loop—2-methylthio-N6-threonylcarbamoyladenosine at position 37 (ms2t6A37), 5-methoxycarbonylmethyl-2-thiouridine at position 34 (mcm5s2U34) and pseudouridine (ψ) at position 39—two of which, ms2t6A37 and mcm5s2U34, are required to achieve wild-type binding activity of wild-type human tRNAUUULys,3 [C. Yarian, M. Marszalek, E. Sochacka, A. Malkiewicz, R. Guenther, A. Miskiewicz and P. F. Agris (2000) Biochemistry, 39, 13390–13395]. Molecular dynamics simulations of nine tRNA anticodon stem–loops with different combinations of nonstandard bases were performed. The wild-type simulation exhibited a canonical anticodon stair-stepped conformation. The ms2t6 modification at position 37 is required for maintenance of this structure and reduces solvent accessibility of U36. Ms2t6A37 generally hydrogen bonds across the loop and may prevent U36 from rotating into solution. A water molecule does coordinate to ψ39 most of the simulation time but weakly, as most of the residence lifetimes are <40 ps.
Collapse
Affiliation(s)
| | | | | | - Maria C. Nagan
- To whom correspondence should be addressed. Tel: +1 660 785 4084; Fax: +1 660 785 4045;
| |
Collapse
|
9
|
|
10
|
Lescrinier E, Nauwelaerts K, Zanier K, Poesen K, Sattler M, Herdewijn P. The naturally occurring N6-threonyl adenine in anticodon loop of Schizosaccharomyces pombe tRNAi causes formation of a unique U-turn motif. Nucleic Acids Res 2006; 34:2878-86. [PMID: 16738127 PMCID: PMC1474066 DOI: 10.1093/nar/gkl081] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Modified nucleosides play an important role in structure and function of tRNA. We have determined the solution structure of the anticodon stem-loop (ASL) of initiator tRNA of Schizosaccharomyces pombe. The incorporation of N6-threonylcarbamoyladenosine at the position 3' to the anticodon triplet (t6A37) results in the formation of a U-turn motif and enhances stacking interactions within the loop and stem regions (i.e. between A35 and t6A37) by bulging out U36. This conformation was not observed in a crystal structure of tRNAi including the same modification in its anticodon loop, nor in the solution structure of the unmodified ASL. A t6A modification also occurs in the well studied anti-stem-loop of lys-tRNA(UUU). A comparison of this stem-loop with our structure demonstrates different effects of the modification depending on the loop sequence.
Collapse
Affiliation(s)
| | | | - Katia Zanier
- EMBL, Structural & Computational Biology and Gene ExpressionMeyerhofstrasse 1, D-69117 Heidelberg, Germany
| | - Koen Poesen
- EMBL, Structural & Computational Biology and Gene ExpressionMeyerhofstrasse 1, D-69117 Heidelberg, Germany
| | - Michael Sattler
- EMBL, Structural & Computational Biology and Gene ExpressionMeyerhofstrasse 1, D-69117 Heidelberg, Germany
| | - Piet Herdewijn
- To whom correspondence should be addressed. Tel: +32 0 16 337387; Fax: +32 0 16 337340;
| |
Collapse
|
11
|
Barthel A, Zacharias M. Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study. Biophys J 2006; 90:2450-62. [PMID: 16399833 PMCID: PMC1403169 DOI: 10.1529/biophysj.105.076158] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.
Collapse
Affiliation(s)
- André Barthel
- School of Engineering and Science, International University Bremen, D-28759 Bremen, Germany
| | | |
Collapse
|
12
|
Koplin J, Mu Y, Richter C, Schwalbe H, Stock G. Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study. Structure 2005; 13:1255-67. [PMID: 16154083 DOI: 10.1016/j.str.2005.05.015] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2005] [Revised: 03/31/2005] [Accepted: 05/13/2005] [Indexed: 10/25/2022]
Abstract
Molecular dynamics simulations of the RNA tetraloop 5'-CGCUUUUGCG-3' with high melting temperature and significant conformational heterogeneity in explicit water solvent are presented and compared to NMR studies. The NMR data allow for a detailed test of the theoretical model, including the quality of the force field and the conformational sampling. Due to the conformational heterogeneity of the tetraloop, high temperature (350 K) and locally enhanced sampling simulations need to be invoked. The Amber98 force field leads to a good overall agreement with experimental data. Based on NMR data and a principal component analysis of the 350 K trajectory, the dynamic structure of the tetraloop is revealed. The principal component free energy surface exhibits four minima, which correspond to well-defined conformational structures that differ mainly by their base stacking in the loop region. No correlation between the motion of the sugar rings and the stacking dynamics of the loop bases is found.
Collapse
Affiliation(s)
- Jessica Koplin
- Institute of Physical and Theoretical Chemistry, Center for Biomolecular Magnetic Resonance, Johann Wolfgang Goethe-University, Frankfurt, Marie-Curie-Str 11, D-60439 Frankfurt/Main, Germany
| | | | | | | | | |
Collapse
|
13
|
Mu Y, Stock G. Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study. Biophys J 2005; 90:391-9. [PMID: 16239331 PMCID: PMC1367046 DOI: 10.1529/biophysj.105.069559] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.
Collapse
Affiliation(s)
- Yuguang Mu
- School of Biological Sciences, Nanyang Technological University, Singapore and School of Physics and Microelectronics, Shandong University, Jinan, China
| | | |
Collapse
|
14
|
Durant PC, Bajji AC, Sundaram M, Kumar RK, Davis DR. Structural effects of hypermodified nucleosides in the Escherichia coli and human tRNALys anticodon loop: the effect of nucleosides s2U, mcm5U, mcm5s2U, mnm5s2U, t6A, and ms2t6A. Biochemistry 2005; 44:8078-89. [PMID: 15924427 DOI: 10.1021/bi050343f] [Citation(s) in RCA: 104] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Previous nuclear magnetic resonance (NMR) studies of unmodified and pseudouridine39-modified tRNA(Lys) anticodon stem loops (ASLs) show that significant structural rearrangements must occur to attain a canonical anticodon loop conformation. The Escherichia coli tRNA(Lys) modifications mnm(5)s(2)U34 and t(6)A37 have indeed been shown to remodel the anticodon loop, although significant dynamic flexibility remains within the weakly stacked U35 and U36 anticodon residues. The present study examines the individual effects of mnm(5)s(2)U34, s(2)U34, t(6)A37, and Mg(2+) on tRNA(Lys) ASLs to decipher how the E. coli modifications accomplish the noncanonical to canonical structural transition. We also investigated the effects of the corresponding human tRNA(Lys,3) versions of the E. coli modifications, using NMR to analyze tRNA ASLs containing the nucleosides mcm(5)U34, mcm(5)s(2)U34, and ms(2)t(6)A37. The human wobble modification has a less dramatic loop remodeling effect, presumably because of the absence of a positive charge on the mcm(5) side chain. Nonspecific magnesium effects appear to play an important role in promoting anticodon stacking. Paradoxically, both t(6)A37 and ms(2)t(6)A37 actually decrease anticodon stacking compared to A37 by promoting U36 bulging. Rather than stack with U36, the t(6)A37 nucleotide in the free tRNAs is prepositioned to form a cross-strand stack with the first codon nucleotide as seen in the recent crystal structures of tRNA(Lys) ASLs bound to the 30S ribosomal subunit. Wobble modifications, t(6)A37, and magnesium each make unique contributions toward promoting canonical tRNA structure in the fundamentally dynamic tRNA(Lys)(UUU) anticodon.
Collapse
Affiliation(s)
- Philippe C Durant
- Department of Medicinal Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | | | | | | | | |
Collapse
|
15
|
Villescas-Diaz G, Zacharias M. Efficient search on energy minima for structure prediction of nucleic acid motifs. J Biomol Struct Dyn 2005; 22:355-64. [PMID: 15473709 DOI: 10.1080/07391102.2004.10507007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.
Collapse
Affiliation(s)
- G Villescas-Diaz
- Institut fur Molekulare Biotechnologie, Beutenbergstr. 11, D-07745 Jena, Germany
| | | |
Collapse
|
16
|
Golebiowski J, Antonczak S, Fernandez-Carmona J, Condom R, Cabrol-Bass D. Closing loop base pairs in RNA loop-loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations. J Mol Model 2004; 10:408-17. [PMID: 15597210 DOI: 10.1007/s00894-004-0216-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2004] [Accepted: 09/10/2004] [Indexed: 10/26/2022]
Abstract
Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.
Collapse
Affiliation(s)
- Jérôme Golebiowski
- Laboratoire Arômes, Synthèses, Interactions, Faculté des sciences, Université de Nice-Sophia Antipolis, Nice Cedex 2, 06108, France.
| | | | | | | | | |
Collapse
|
17
|
Rohs R, Sklenar H. Methylene Blue Binding to DNA with Alternating AT Base Sequence: Minor Groove Binding is Favored over Intercalation. J Biomol Struct Dyn 2004; 21:699-711. [PMID: 14769063 DOI: 10.1080/07391102.2004.10506960] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
The results presented in this paper on methylene blue (MB) binding to DNA with AT alternating base sequence complement the data obtained in two former modeling studies of MB binding to GC alternating DNA. In the light of the large amount of experimental data for both systems, this theoretical study is focused on a detailed energetic analysis and comparison in order to understand their different behavior. Since experimental high-resolution structures of the complexes are not available, the analysis is based on energy minimized structural models of the complexes in different binding modes. For both sequences, four different intercalation structures and two models for MB binding in the minor and major groove have been proposed. Solvent electrostatic effects were included in the energetic analysis by using electrostatic continuum theory, and the dependence of MB binding on salt concentration was investigated by solving the non-linear Poisson-Boltzmann equation. We find that the relative stability of the different complexes is similar for the two sequences, in agreement with the interpretation of spectroscopic data. Subtle differences, however, are seen in energy decompositions and can be attributed to the change from symmetric 5'-YpR-3' intercalation to minor groove binding with increasing salt concentration, which is experimentally observed for the AT sequence at lower salt concentration than for the GC sequence. According to our results, this difference is due to the significantly lower non-electrostatic energy for the minor groove complex with AT alternating DNA, whereas the slightly lower binding energy to this sequence is caused by a higher deformation energy of DNA. The energetic data are in agreement with the conclusions derived from different spectroscopic studies and can also be structurally interpreted on the basis of the modeled complexes. The simple static modeling technique and the neglect of entropy terms and of non-electrostatic solute-solvent interactions, which are assumed to be nearly constant for the compared complexes of MB with DNA, seem to be justified by the results.
Collapse
Affiliation(s)
- Remo Rohs
- Theoretical Biophysics Group, Max Delbrück Center for Molecular Medicine, Robert-Rössle-Str 10, 13092 Berlin, Germany.
| | | |
Collapse
|
18
|
Villescas-Diaz G, Zacharias M. Sequence context dependence of tandem guanine:adenine mismatch conformations in RNA: a continuum solvent analysis. Biophys J 2003; 85:416-25. [PMID: 12829496 PMCID: PMC1303097 DOI: 10.1016/s0006-3495(03)74486-5] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Guanine:adenine (G:A) mismatches and in particular tandem G:A (tG:A) mismatches are frequently observed in biological RNA molecules and can serve as sites for tertiary interaction, metal binding and protein recognition. Depending on the surrounding sequence tG:A mismatches can adopt different basepairing topologies. In the sequence context (5'-) GGAC (tandem G:A in bold) a face-to-face (imino or Watson-Crick-like) pairing is preferred whereas in the CGAG context, G and A adopt a sheared arrangement. Systematic conformational searches with a generalized Born continuum model and molecular dynamics simulations including explicit water molecules and ions have been used to generate face-to-face and sheared tG:A mismatches in both CGAG and GGAC sequence contexts. Conformations from both approaches were evaluated using the same force field and a Poisson-Boltzmann continuum solvent model. Although the substate analysis predicted the sheared arrangement to be energetically preferred in both sequence contexts, a significantly greater preference of the sheared form was found for the CGAG context. In agreement with the experimental observation, the analysis of molecular dynamics trajectories indicated a preference of the sheared form in the case of the CGAG-context and a favorization of the face-to-face form in the case of the GGAC context. The computational studies allowed to identify energetic contributions that stabilize or destabilize the face-to-face and sheared tandem mismatch topologies. The calculated nonpolar solvation and Lennard-Jones packing interaction were found to stabilize the sheared topology independent of the sequence context. Electrostatic contributions are predicted to make the most significant contribution to the sequence context dependence on the structural preference of tG:A mismatches.
Collapse
|
19
|
Shi K, Pan B, Sundaralingam M. The crystal structure of an alternating RNA heptamer r(GUAUACA) forming a six base-paired duplex with 3'-end adenine overhangs. Nucleic Acids Res 2003; 31:1392-7. [PMID: 12595546 PMCID: PMC149819 DOI: 10.1093/nar/gkg226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The crystal structure of an alternating RNA heptamer r(GUAUACA) has been determined to 2.0 A resolution and refined to an R(work) of 17.1% and R(free) of 18.5% using 2797 reflections. The heptamer crystallized in the space group C222 with a unit cell of a = 25.74, b = 106.58, c = 30.26 A and two independent strands in the asymmetric unit. Each heptamer forms a duplex with its symmetry-related strand and each duplex contains six Watson-Crick base pairs and 3'-end adenosine overhangs. Therefore, two kinds of duplex (duplex 1 and duplex 2) are formed. Duplexes 1 stack on each other forming a pseudo-continuous column, which is typical of the RNA packing mode, while duplex 2 is typical of A-DNA packing with its termini in abutting interactions. Overhang adenine residues stack within the duplexes with C3'-endo sugar pucker and C2'-endo sugar pucker in duplexes 1 and 2, respectively. A Na+ ion in the crystal lattice is water bridged to two N1 atoms of symmetry-related A7 bases.
Collapse
Affiliation(s)
- Ke Shi
- Department of Chemistry, The Ohio State University, 200 Johnston Laboratory, Columbus, OH 43210, USA
| | | | | |
Collapse
|
20
|
Bui CT, Rees K, Lambrinakos A, Bedir A, Cotton RGH. Site-selective reactions of imperfectly matched DNA with small chemical molecules: applications in mutation detection. Bioorg Chem 2002; 30:216-32. [PMID: 12392702 DOI: 10.1016/s0045-2068(02)00019-6] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
The last decade has witnessed many exciting scientific publications associated with site-selective reactions of small chemical molecules with imperfectly matched DNA. Typical examples are carbodiimide, hydroxylamine, potassium permanganate, osmium tetroxide, chemical tagging probes, biotinylated, chemiluminescent and fluorescent probes, and all of them selectively react with imperfectly matched DNA. More recently, some therapeutic agents including DNA intercalating drugs and groove binders have been found to promote the in vivo repair system to recognize and repair the mismatch more effectively. The results have established a novel method for detection of mismatches. Development of new chemical reactions for detection of imperfectly matched DNA and mutations is a rapidly growing field and has attracted significant interest of scientists from both chemical and biological fields and it is the main focus of this review.
Collapse
Affiliation(s)
- Chinh T Bui
- Genomic Disorders Research Centre, The University of Melbourne, 7th Floor, Daly Wing, 35 Victoria Parade, Fitzroy, Vic 3065, Melbourne, Australia.
| | | | | | | | | |
Collapse
|
21
|
|
22
|
Burkhardt C, Zacharias M. Modelling ion binding to AA platform motifs in RNA: a continuum solvent study including conformational adaptation. Nucleic Acids Res 2001; 29:3910-8. [PMID: 11574672 PMCID: PMC60250 DOI: 10.1093/nar/29.19.3910] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Binding of monovalent and divalent cations to two adenine-adenine platform structures from the Tetrahymena group I intron ribozyme has been studied using continuum solvent models based on the generalised Born and the finite-difference Poisson-Boltzmann approaches. The adenine-adenine platform RNA motif forms an experimentally characterised monovalent ion binding site important for ribozyme folding and function. Qualitative agreement between calculated and experimental ion placements and binding selectivity was obtained. The inclusion of solvation effects turned out to be important to obtain low energy structures and ion binding placements in agreement with the experiment. The calculations indicate that differences in solvation of the isolated ions contribute to the calculated ion binding preference. However, Coulomb attraction and van der Waals interactions due to ion size differences and RNA conformational adaptation also influence the calculated ion binding affinity. The calculated alkali ion binding selectivity for both platforms followed the order K(+) > Na(+) > Rb(+) > Cs(+) > Li(+) (Eisenman series VI) in the case of allowing RNA conformational relaxation during docking. With rigid RNA an Eisenman series V was obtained (K(+) > Rb(+) > Na(+) > Cs(+) > Li(+)). Systematic energy minimisation docking simulations starting from several hundred initial placements of potassium ions on the surface of platform containing RNA fragments identified a coordination geometry in agreement with the experiment as the lowest energy binding site. The approach could be helpful to identify putative ion binding sites in nucleic acid structures determined at low resolution or with experimental methods that do not allow identification of ion binding sites.
Collapse
Affiliation(s)
- C Burkhardt
- AG Theoretische Biophysik, Institut für Molekulare Biotechnologie, Beutenbergstrasse 11, D-07745 Jena, Germany
| | | |
Collapse
|
23
|
Abstract
During the past year, major improvements have been made in methods used to solve RNA structures from crystals, find RNA patterns in sequence data and determine RNA secondary structure. Computational methods for assisting an interactive computer graphics human modeler, searching the conformational space of RNA tertiary structure, studying the dynamics of complexes involving RNA and simulating RNA catalytic activities have also been advanced.
Collapse
Affiliation(s)
- F Major
- Département d'Informatique et de Recherche Opérationnelle, Université de Montréal, CP 6128, Succ Centre-Ville, Montréal, Québec, H3C 3J7, Canada.
| | | |
Collapse
|
24
|
Zacharias M. Conformational analysis of DNA-trinucleotide-hairpin-loop structures using a continuum solvent model. Biophys J 2001; 80:2350-63. [PMID: 11325735 PMCID: PMC1301424 DOI: 10.1016/s0006-3495(01)76205-4] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
A number of trinucleotide sequences in DNA can form compact and stable hairpin loops that may have significance for DNA replication and transcription. The conformational analysis of these motifs is important for an understanding of the function and design of nucleic acid structures. Extensive conformational searches have been performed on three experimentally known trinucleotide hairpin loops (AGC, AAA, and GCA) closed by a four-base-pair stem. An implicit solvation model based on the generalized Born method has been employed during energy minimization and conformational search. In addition, energy-minimized conformers were evaluated using a finite-difference Poisson-Boltzmann approach. For all three loop sequences, conformations close to experiment were found as lowest-energy structures among several thousand alternative energy minima. The inclusion of reaction-field contributions was found to be important for a realistic conformer ranking. Most generated hairpin loop structures within approximately 5 kcal x mol(-1) of the lowest-energy structure have a similar topology. Structures within approximately 10 kcal x mol(-1) could be classified into about five structural families representing distinct arrangements of loop nucleotides. Although a large number of backbone torsion angle combinations were compatible with each structural class, some specific patterns could be identified. Harmonic mode analysis was used to account for differences in conformational flexibility of low-energy sub-states. Class-specific differences in the pattern of atomic fluctuations along the sequence were observed; however, inclusion of conformational entropy contributions did not change ranking of structural classes. For an additional loop sequence (AAG) with no available experimental structure, the approach suggests a lowest-energy loop topology overall similar to the other three loop sequences but closed by a different non-canonical base-pairing scheme.
Collapse
Affiliation(s)
- M Zacharias
- AG Theoretische Biophysik, Institut für Molekulare Biotechnologie, 07745 Jena, Germany.
| |
Collapse
|
25
|
Sarzynska J, Kulinski T, Nilsson L. Conformational dynamics of a 5S rRNA hairpin domain containing loop D and a single nucleotide bulge. Biophys J 2000; 79:1213-27. [PMID: 10968986 PMCID: PMC1301018 DOI: 10.1016/s0006-3495(00)76376-4] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022] Open
Abstract
Molecular modeling and molecular dynamics have been employed to study the conformation and flexibility of a 15-nucleotide fragment of the plant 5S rRNA containing loop D and a single uridine bulge. Two different model built initial structures were used: one with the bulge localized inside the helical stem and another with the bulge pointing out from the helix. Several independent 700-ps-long trajectories in aqueous solution with Na(+) conterions were produced for each starting structure. The bulge nucleotide inside the helix stayed in two main conformations, both of which affected the geometry of the stem part opposite the bulge. When the bulge nucleotide was located outside the helix, we found high base mobility and local backbone flexibility. The dynamics of the hydrogen bond network and conformational changes from a direct to a water mediated hydrogen bond in the sheared G-A basepair in the tetraloop was described. Our results correlate with lead ion induced cleavage patterns in 5S rRNA. Sites resistant to nonspecific lead cleavage appeared in all our simulations as the most rigid fragments independent of the localization of the bulge nucleotide.
Collapse
Affiliation(s)
- J Sarzynska
- Center for Structural Biochemistry, Department of Bioscience at NOVUM, Karolinska Institutet, SE-141 57 Huddinge, Sweden
| | | | | |
Collapse
|
26
|
Thiviyanathan V, Guliaev AB, Leontis NB, Gorenstein DG. Solution conformation of a bulged adenosine base in an RNA duplex by relaxation matrix refinement. J Mol Biol 2000; 300:1143-54. [PMID: 10903860 DOI: 10.1006/jmbi.2000.3931] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Bulges are common structural motifs in RNA secondary structure and are thought to play important roles in RNA-protein and RNA-drug interactions. Adenosine bases are the most commonly occurring unpaired base in double helical RNA secondary structures. The solution conformation and dynamics of a 25-nucleotide RNA duplex containing an unpaired adenosine, r(GGCAGAGUGCCGC): r(GCGGCACCUGCC) have been studied by NMR spectroscopy and MORASS iterative relaxation matrix structural refinement. The results show that the bulged adenosine residue stacks into the RNA duplex with little perturbation around the bulged region. Most of the bases in the RNA duplex adopt C(3)'-endo conformation, exhibiting the N-type sugar pucker as found in the A form helices. The sugars of the bulged residue and the 5' flanking residue to it are found to exhibit C(2)'-endo conformation. None of the residues are in syn conformation.
Collapse
Affiliation(s)
- V Thiviyanathan
- Sealy Center for Structural Biology and Department of Human Biological Chemistry and Genetics, University of Texas Medical Branch, Galveston, TX 77555-1157, USA
| | | | | | | |
Collapse
|
27
|
Abstract
Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.
Collapse
Affiliation(s)
- M Zacharias
- AG Theoretische Biophysik, Institut für Molekulare Biotechnologie, Jena, D-07745, Germany.
| |
Collapse
|
28
|
Sudarsanakumar C, Xiong Y, Sundaralingam M. Crystal structure of an adenine bulge in the RNA chain of a DNA.RNA hybrid, d(CTCCTCTTC).r(gaagagagag). J Mol Biol 2000; 299:103-12. [PMID: 10860725 DOI: 10.1006/jmbi.2000.3730] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Crystal structure of a DNA.RNA hybrid, d(CTCCTCTTC).r(gaagagagag), with an adenine bulge in the polypurine RNA strand was determined at 2.3 A resolution. The structure was solved by the molecular replacement method and refined to a final R-factor of 19.9% (Rfree 22.2%). The hybrid duplex crystallized in the space group I222 with unit cell dimensions, a = 46.66 A, b = 47.61 A and c = 54.05 A, and adopts the A-form conformation. All RNA and DNA sugars are in the C3'-endo conformation, the glycosyl angles in anti conformation and the majority of the C4'-C5' torsion angles in g+ except two trans angles, in conformity with the C3'-endo rigid nucleotide hypothesis. The adenine bulge is looped out and it is also in the anti C3'-endo conformation. The bulge is involved in a base-triple (C.g)*a interaction with the end base-pair (C9.g10) in the minor groove of a symmetry-related molecule. The 2' hydroxyl group of g15 is hydrogen bonded to O2P and O5' of g17, skipping the bulged adenine a16 and stabilizing the sugar-phosphate backbone of the hybrid. The hydrogen bonding and the backbone conformation at the bulged adenine site is very similar to that found in the crystal structure of a protein-RNA complex.
Collapse
Affiliation(s)
- C Sudarsanakumar
- Department of Chemistry, Biochemistry, Ohio State University, Columbus 43210, USA
| | | | | |
Collapse
|
29
|
Abstract
The harmonic mode analysis method was used to characterize the conformational deformability of regular Watson-Crick paired, mismatch- and bulge-containing RNA. Good agreement between atomic Debye-Waller factors derived from x-ray crystallography of a regular RNA oligonucleotide and calculated atomic fluctuations was obtained. Calculated helical coordinate fluctuations showed a small sequence dependence of up to approximately 30-50%. A negative correlation between motions at a given base pair step and neighboring steps was found for most helical coordinates. Only very few calculated modes contribute significantly to global motions such as bending, twisting, and stretching of the RNA molecules. With respect to a local helical description of the RNA helix our calculations suggest that RNA bending is mostly due to a periodic change in the base pair step descriptors slide and roll. The presence of single guanine:uridine or guanine:adenine mismatches had little influence on the calculated RNA flexibility. In contrast, for tandem guanine:adenine base pairs the harmonic mode approach predicts a significantly reduced conformational flexibility in the case of a sheared arrangement and slightly enhanced flexibility for a face-to-face (imino proton) pairing relative to regular RNA. The presence of a single extra adenine bulge nucleotide stacked between flanking sequences resulted in an increased local atomic mobility around the bulge site (approximately 40%) and a slightly enhanced global bending flexibility. For an adenine bulge nucleotide in a looped-out conformation a strongly enhanced bulge nucleotide mobility but no increased bending flexibility compared to regular RNA was found.
Collapse
Affiliation(s)
- M Zacharias
- Institute for Molecular Biotechnology, 07745 Jena, Germany.
| | | |
Collapse
|
30
|
Rohs R, Sklenar H, Lavery R, Röder B. Methylene Blue Binding to DNA with Alternating GC Base Sequence: A Modeling Study. J Am Chem Soc 2000. [DOI: 10.1021/ja992966k] [Citation(s) in RCA: 134] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Remo Rohs
- Contribution from the Arbeitsgruppe Theoretische Biophysik, Max-Delbrück-Centrum für Molekulare Medizin, Robert-Rössle-Strasse 10, D-13125 Berlin, Germany, Laboratoire de Biochimie Théorique CNRS UPR 9080, Institut de Biologie Physico-Chimique, 13, rue Pierre et Marie Curie, F-75005 Paris, France, and Arbeitsgruppe Photobiophysik, Humboldt-Universität zu Berlin, Invalidenstrasse 110, D-10115 Berlin, Germany
| | - Heinz Sklenar
- Contribution from the Arbeitsgruppe Theoretische Biophysik, Max-Delbrück-Centrum für Molekulare Medizin, Robert-Rössle-Strasse 10, D-13125 Berlin, Germany, Laboratoire de Biochimie Théorique CNRS UPR 9080, Institut de Biologie Physico-Chimique, 13, rue Pierre et Marie Curie, F-75005 Paris, France, and Arbeitsgruppe Photobiophysik, Humboldt-Universität zu Berlin, Invalidenstrasse 110, D-10115 Berlin, Germany
| | - Richard Lavery
- Contribution from the Arbeitsgruppe Theoretische Biophysik, Max-Delbrück-Centrum für Molekulare Medizin, Robert-Rössle-Strasse 10, D-13125 Berlin, Germany, Laboratoire de Biochimie Théorique CNRS UPR 9080, Institut de Biologie Physico-Chimique, 13, rue Pierre et Marie Curie, F-75005 Paris, France, and Arbeitsgruppe Photobiophysik, Humboldt-Universität zu Berlin, Invalidenstrasse 110, D-10115 Berlin, Germany
| | - Beate Röder
- Contribution from the Arbeitsgruppe Theoretische Biophysik, Max-Delbrück-Centrum für Molekulare Medizin, Robert-Rössle-Strasse 10, D-13125 Berlin, Germany, Laboratoire de Biochimie Théorique CNRS UPR 9080, Institut de Biologie Physico-Chimique, 13, rue Pierre et Marie Curie, F-75005 Paris, France, and Arbeitsgruppe Photobiophysik, Humboldt-Universität zu Berlin, Invalidenstrasse 110, D-10115 Berlin, Germany
| |
Collapse
|
31
|
Ayadi L, Coulombeau C, Lavery R. Abasic sites in duplex DNA: molecular modeling of sequence-dependent effects on conformation. Biophys J 1999; 77:3218-26. [PMID: 10585943 PMCID: PMC1300592 DOI: 10.1016/s0006-3495(99)77152-3] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Molecular modeling calculations using JUnction Minimization of Nucleic Acids (JUMNA) have been used to study sequence effects on the conformation of abasic sites within duplex DNA. We have considered lesions leading to all possible unpaired bases (X), adenine, guanine, cytosine, or thymine contained within two distinct sequence contexts, CXC and GXG. Calculations were carried out on DNA 11-mers using extensive conformational search techniques to locate the most stable abasic conformations and using Poisson-Boltzmann corrected electrostatics to account for solvation effects. The results, which are in very good agreement with available experimental data, point to strong sequence effects on both the position of the unpaired base (intra or extrahelical) and on the overall curvature induced by the abasic lesion. For CXC, unpaired purines are found to lie within the helix, while unpaired pyrimidines are either extrahelical or in equilibrium between the intra and extrahelical forms. For GXG, all unpaired bases lead to intrahelical forms, but with marked, sequence-dependent differences in induced curvature.
Collapse
Affiliation(s)
- L Ayadi
- Laboratoire d'Etudes Dynamiques et Structurales de la Sélectivité, Laboratoire de Chimie Bioorganique, Université Joseph Fourier, F-38041 Grenoble, Cedex 9, France
| | | | | |
Collapse
|
32
|
Ennifar E, Yusupov M, Walter P, Marquet R, Ehresmann B, Ehresmann C, Dumas P. The crystal structure of the dimerization initiation site of genomic HIV-1 RNA reveals an extended duplex with two adenine bulges. Structure 1999; 7:1439-49. [PMID: 10574792 DOI: 10.1016/s0969-2126(00)80033-7] [Citation(s) in RCA: 122] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
BACKGROUND An important step in retroviral replication is dimerization of the genomic RNA prior to encapsidation. Dimerization is initiated by the formation of a transient 'kissing-loop complex' that is thought to be subsequently matured into an extended duplex by the nucleocapsid protein (NCp). Although chemical probing and nuclear magnetic resonance spectroscopy have provided insight into the structure of the kissing-loop structure, no structural information concerning the extended-duplex state is available so far. RESULTS The structure of a minimal HIV-1 RNA dimerization initiation site has been solved at 2.3 A resolution in two different space groups. It reveals a 22 base pair extended duplex with two noncanonical Watson-Crick-like G-A mismatches, each adjacent to a bulged-out adenine. The structure shows significant asymmetry in deep groove width and G-A base-pair conformations. A network of eight magnesium cations was clearly identified, one being unusually chelated by the 3' phosphate of each bulge across an extremely narrowed deep major groove. CONCLUSIONS These crystal structures represent the putative matured form of the initial kissing-loop complex. They show the ability of this self-complementary RNA hairpin loop to acquire a more stable extended duplex structure. Both bulged adenines form a striking 'base grip' that could be a recognition signal, either in cis for another viral RNA sequence, or in trans for a protein, possibly the NCp. Magnesium binding might be important to promote and stabilize the observed extrahelical conformation of these bulges.
Collapse
Affiliation(s)
- E Ennifar
- UPR9002 - Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg cedex, F-67084, France
| | | | | | | | | | | | | |
Collapse
|