The ability to control pre-mRNA splicing with small molecules could facilitate the development of... more The ability to control pre-mRNA splicing with small molecules could facilitate the development of therapeutics or cell-based circuits that control gene function. Myotonic dystrophy type 1 (DM1) is caused by the dysregulation of alternative pre-mRNA splicing due to sequestration of muscleblind-like 1 protein (MBNL1) by expanded, non-coding r(CUG) repeats (r(CUG) exp ). Here we report two small molecules that induce or ameliorate alternative splicing dysregulation. The thiophene-containing small molecule (1) inhibits the interaction of MBNL1 with its natural pre-mRNA substrates. Compound (2), a substituted naphthyridine, binds r(CUG) exp and displaces MBNL1. Structural models show that 1 binds MBNL1 in the Zn-finger domain and that 2 interacts with UU loops in r(CUG) exp . This study provides a structural framework for small molecules that target MBNL1 by mimicking r(CUG) exp and shows that targeting MBNL1 causes dysregulation of alternative splicing, suggesting that MBNL1 is thus not a suitable therapeutic target for the treatment of DM1.
RNA is an immensely important target for small molecule therapeutics or chemical probes of functi... more RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.
Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expans... more Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode…
Myotonic Dystrophy 1 (DM1) is a genetic disease caused by expansion of CTG repeats in DNA. Once t... more Myotonic Dystrophy 1 (DM1) is a genetic disease caused by expansion of CTG repeats in DNA. Once transcribed, these repeats form RNA hairpins with repeating 1×1 nucleotide UU internal loop motifs, r(CUG)n, which attract muscleblind-like 1 (MBNL1) protein leading to the disease. In DM1 CUG can be repeated thousands of times, so these structures are intractable to characterization using structural biology. However, inhibition of MBNL1-r(CUG)n binding requires a detailed analysis of the 1×1 UU internal loops. In this contribution we employ regular and umbrella sampling molecular dynamics (MD) simulations to describe the structural and thermodynamic properties of 1×1 UU internal loops. Calculations were run on a reported crystal structure and a designed system, which mimics an infinitely long RNA molecule with continuous CUG repeats. Two-dimensional (2D) potential of mean force (PMF) surfaces were created by umbrella sampling, and the discrete path sampling (DPS) method was utilized to investigate the energy landscape of 1×1 UU RNA internal loops, revealing that 1×1 UU base pairs are dynamic and strongly prefer the anti-anti conformation. Two 2D PMF surfaces were calculated for the 1×1 UU base pairs, revealing several local minima and three syn-anti ↔ anti-anti transformation pathways. Although at room temperature the syn-anti ↔ anti-anti transformation is not observed on the MD time scale, one of these pathways dominates the dynamics of the 1×1 UU base pairs in temperature jump MD simulations. This mechanism has now been treated successfully using the DPS approach. Our results suggest that local minima predicted by umbrella sampling calculations could be stabilized by small molecules, which is of great interest for future drug design. Furthermore, distorted GC/CG conformations may be important in understanding how MBNL1 binds to RNA CUG repeats. Hence we provide new insight into the dynamic roles of RNA loops and their contributions to presently incurable diseases.
Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are ... more Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against (1)H-(1)H NOESY distances and (3) J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.
The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relati... more The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relative to the Watson-Crick GC pair, provides an expanded alphabet for understanding interactions that shape nucleic acid structure. Here, thermodynamic stabilities of tandem GA pairs flanked by iGiC pairs are reported along with the NMR structures of the RNA self-complementary duplexes (GCiGGAiCGCA) 2 and (GGiCGAiGCCA) 2 . A sheared GA pairing forms in (GCiGGAiCGCA) 2 , and an imino GA pairing forms in (GGiCGAiGCCA) 2 . The structures contrast with the formation of tandem imino and sheared GA pairs flanked by GC pairs in the RNA self-complementary duplexes (GCGGACGC) 2 and (GGCGAGCC) 2 , respectively. In both iGiC duplexes, Watson-Crick-like hydrogen bonds are formed between iG and iC, and iGiC substitutions result in less favorable loop stability. The results provide benchmarks for testing computations of molecular interactions that shape RNA three-dimensional structure.
Proceedings of the National Academy of Sciences, 2008
FRET between the zinc porphyrin (ZnP) chromophore in zincsubstituted cytochrome c (Zn-cyt c) and ... more FRET between the zinc porphyrin (ZnP) chromophore in zincsubstituted cytochrome c (Zn-cyt c) and an Alexa Fluor dye attached to specific surface sites was used to characterize Zn-cyt c unfolding. The use of ZnP as a fluorescent acceptor eliminates the need to doubly label the protein with exogenous dyes to perform FRET experiments in which both donor and acceptor fluorescence is monitored. The requirement for attachment of only one dye also minimizes perturbation to the protein. This sensitive technique allowed for the determination of distances between the label placed at six different sites and ZnP through a range of denaturant concentrations. Fitting of the data to a three-state model provides distances in the unfolding intermediate. The use of ZnP as a fluorescent acceptor of energy in FRET has a significant potential for application to a range of other systems including heme-binding proteins and proteins to which a covalently attached heme tag may be added.
RNA is an immensely important target for small molecule therapeutics or chemical probes of functi... more RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.
One class of functionally important RNA is repeating transcripts that cause disease through vario... more One class of functionally important RNA is repeating transcripts that cause disease through various mechanisms. For example, expanded CAG repeats can cause Huntington's and other disease through translation of toxic proteins. Herein, a crystal structure of r[5'UUGGGC(CAG)3GUCC]2, a model of CAG expanded transcripts, refined to 1.65 Å resolution is disclosed that shows both anti-anti and syn-anti orientations for 1 × 1 nucleotide AA internal loops. Molecular dynamics (MD) simulations using AMBER force field in explicit solvent were run for over 500 ns on the model systems r(5'GCGCAGCGC)2 (MS1) and r(5'CCGCAGCGG)2 (MS2). In these MD simulations, both anti-anti and syn-anti AA base pairs appear to be stable. While anti-anti AA base pairs were dynamic and sampled multiple anti-anti conformations, no syn-anti ↔ anti-anti transformations were observed. Umbrella sampling simulations were run on MS2, and a 2D free energy surface was created to extract transformation pathways. In addition, an explicit solvent MD simulation over 800 ns was run on r[5'GGGC(CAG)3GUCC]2, which closely represents the refined crystal structure. One of the terminal AA base pairs (syn-anti conformation), transformed to anti-anti conformation. The pathway followed in this transformation was the one predicted by umbrella sampling simulations. Further analysis showed a binding pocket near AA base pairs in syn-anti conformations. Computational results combined with the refined crystal structure show that global minimum conformation of 1 × 1 nucleotide AA internal loops in r(CAG) repeats is anti-anti but can adopt syn-anti depending on the environment. These results are important to understand RNA dynamic-function relationships and to develop small molecules that target RNA dynamic ensembles.
All-atom force fields are important for predicting thermodynamic, structural, and dynamic propert... more All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG) 2 , r(GGCC) 2 , r(GCGC) 2 , and r(CGCG) 2 are replaced by isocytidineÀisoguanosine (iCiG) pairs. Agreement with experiment was improved when ε/ζ, α/γ, β, and χ torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG f iCiCiGiG, GGCC f iGiGiCiC, GCGC f iGiCiGiC, and CGCG f iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC f iGiGiCiC and GCGC f iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG) 2 retained A-form conformation for g82% of the time. This is consistent with NMR spectra of r(iGiGiCiC) 2 , which reveal an A-form conformation. In MD simulations, r(CCGG) 2 retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods.
Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loop... more Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair and finally by a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and basephosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 µs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by MM-PBSA free energy calculations, locally enhanced sampling (LES) runs, targeted MD (TMD) and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and netneutralizing Na + vs. excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the Corresponding authors: Kamila Réblová, . SUPPORTING INFORMATION Detailed description of standard MD simulations run with parm99, description of standard MD simulations run with parmbsc0, and description of LES simulations of ribosomal H40. Simulated annealing protocol used for minimization in the NEB trials (in ). Calculation of free energy basin around H40 X-ray structure based on parm99 and parmbsc0 simulations using potential of mean force coordinate method. . This material is available free of charge via the Internet at major groove, loss of base-phosphate interaction and other instabilities. The alternative X-ray geometry even undergoes conformational transition towards the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.
A reparameterization of the torsional parameters for the glycosidic dihedral angle, , for the AMB... more A reparameterization of the torsional parameters for the glycosidic dihedral angle, , for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99 . Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99 force fields are compared with NMR results. For each nucleoside and force field, 10 individual molecular dynamics simulations of 30 ns each were run. For cytidine with AMBER99 force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) 1 H, steady-state 1D 1 H nuclear Overhauser effect (NOE), and transient 1D 1 H NOE, was used to determine the sugar puckering and preferred base orientation with respect to the ribose of cytidine and uridine. The AMBER99 force field overestimates the population of syn conformations of the base orientation and of C2′-endo sugar puckering of the pyrimidines, while the AMBER99 force field's predictions are more consistent with NMR results. Moreover, the AMBER99 force field prefers high anti conformations with glycosidic dihedral angles around 310°for the base orientation of purines. The AMBER99 force field prefers anti conformations around 185°, which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99 force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.
Locked Nucleic Acids (LNA) considerably enhance the thermodynamic stability of DNA and RNA duplex... more Locked Nucleic Acids (LNA) considerably enhance the thermodynamic stability of DNA and RNA duplexes. We report the thermodynamic stabilities of LNA-2'-O-methyl RNA/RNA duplexes designed to provide insight into the contributions of stacking and hydrogen bonding interactions to the enhanced stability. The results show that hydrogen bonding of LNA nucleotides is similar to that of 2'-O-methyl RNA nucleotides, whereas the 3'-stacking interactions are on average about 0.7 kcal/mol more favorable at 37 °C than for 2'-O-methyl or RNA nucleotides. Moreover, NMR spectra suggest helical pre-organization of the single stranded tetramer, C L A M A L U M , probably due to restriction of some torsion angles. Thus enhanced stacking interactions and helical pre-organization of single stranded oligonucleotides contribute to the extraordinary stabilization of duplexes by LNA nucleotides.
Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pai... more Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pairwise hydrogen-bonding models. Such effects include the stabilities of isoguanosine-isocytidine (iG-iC) base pairs and of various 2 x 2 nucleotide internal loops. Presumably, these results can be explained by base stacking effects, which can be partitioned into Coulombic and overlap effects. We review experimental measurements that provide benchmarks for testing the approximations and theories used for modeling nucleic acids. Quantitative agreement between experiment and theory will indicate understanding of the interactions determining RNA stability and structure.
Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expans... more Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode occurs via the major groove, in which one C residue is unstacked and the cross-strand…
Many transcriptional activators act at a distance from core promoter elements and work by recruit... more Many transcriptional activators act at a distance from core promoter elements and work by recruiting RNA polymerase through protein-protein interactions. We show here how the prokaryotic regulatory protein CueR both represses and activates transcription by differentially modulating local DNA structure within the promoter. Structural studies reveal that the repressor state slightly bends the promoter DNA, precluding optimal RNA polymerase-promoter recognition. Upon binding a metal ion in the allosteric site, CueR switches into an activator conformation. It maintains all protein-DNA contacts but introduces torsional stresses that kink and undertwist the promoter, stabilizing an A-form DNA-like conformation. These factors switch on and off transcription by exerting dynamic control of DNA stereochemistry, reshaping the core promoter and making it a better or worse substrate for polymerase.
Spinocerebellar ataxia type 10 (SCA10) is caused by a pentanucleotide repeat expansion of r(AUUCU... more Spinocerebellar ataxia type 10 (SCA10) is caused by a pentanucleotide repeat expansion of r(AUUCU) within intron 9 of the ATXN10 pre-mRNA. The RNA causes disease by a gain-of-function mechanism in which it inactivates proteins involved in RNA biogenesis. Spectroscopic studies showed that r(AUUCU) repeats form a hairpin structure; however, there were no high-resolution structural models prior to this work. Herein, we report the first crystal structure of model r(AUUCU) repeats refined to 2.8 Å and analysis of the structure via molecular dynamics simulations. The r(AUUCU) tracts adopt an overall A-form geometry in which 3 × 3 nucleotide (5')UCU(3')/(3')UCU(5') internal loops are closed by AU pairs. Helical parameters of the refined structure as well as the corresponding electron density map on the crystallographic model reflect dynamic features of the internal loop. The computational analyses captured dynamic motion of the loop closing pairs, which can form single-stra...
Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are ... more Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against (1)H-(1)H NOESY distances and (3) J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.
The ability to control pre-mRNA splicing with small molecules could facilitate the development of... more The ability to control pre-mRNA splicing with small molecules could facilitate the development of therapeutics or cell-based circuits that control gene function. Myotonic dystrophy type 1 (DM1) is caused by the dysregulation of alternative pre-mRNA splicing due to sequestration of muscleblind-like 1 protein (MBNL1) by expanded, non-coding r(CUG) repeats (r(CUG) exp ). Here we report two small molecules that induce or ameliorate alternative splicing dysregulation. The thiophene-containing small molecule (1) inhibits the interaction of MBNL1 with its natural pre-mRNA substrates. Compound (2), a substituted naphthyridine, binds r(CUG) exp and displaces MBNL1. Structural models show that 1 binds MBNL1 in the Zn-finger domain and that 2 interacts with UU loops in r(CUG) exp . This study provides a structural framework for small molecules that target MBNL1 by mimicking r(CUG) exp and shows that targeting MBNL1 causes dysregulation of alternative splicing, suggesting that MBNL1 is thus not a suitable therapeutic target for the treatment of DM1.
RNA is an immensely important target for small molecule therapeutics or chemical probes of functi... more RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.
Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expans... more Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode…
Myotonic Dystrophy 1 (DM1) is a genetic disease caused by expansion of CTG repeats in DNA. Once t... more Myotonic Dystrophy 1 (DM1) is a genetic disease caused by expansion of CTG repeats in DNA. Once transcribed, these repeats form RNA hairpins with repeating 1×1 nucleotide UU internal loop motifs, r(CUG)n, which attract muscleblind-like 1 (MBNL1) protein leading to the disease. In DM1 CUG can be repeated thousands of times, so these structures are intractable to characterization using structural biology. However, inhibition of MBNL1-r(CUG)n binding requires a detailed analysis of the 1×1 UU internal loops. In this contribution we employ regular and umbrella sampling molecular dynamics (MD) simulations to describe the structural and thermodynamic properties of 1×1 UU internal loops. Calculations were run on a reported crystal structure and a designed system, which mimics an infinitely long RNA molecule with continuous CUG repeats. Two-dimensional (2D) potential of mean force (PMF) surfaces were created by umbrella sampling, and the discrete path sampling (DPS) method was utilized to investigate the energy landscape of 1×1 UU RNA internal loops, revealing that 1×1 UU base pairs are dynamic and strongly prefer the anti-anti conformation. Two 2D PMF surfaces were calculated for the 1×1 UU base pairs, revealing several local minima and three syn-anti ↔ anti-anti transformation pathways. Although at room temperature the syn-anti ↔ anti-anti transformation is not observed on the MD time scale, one of these pathways dominates the dynamics of the 1×1 UU base pairs in temperature jump MD simulations. This mechanism has now been treated successfully using the DPS approach. Our results suggest that local minima predicted by umbrella sampling calculations could be stabilized by small molecules, which is of great interest for future drug design. Furthermore, distorted GC/CG conformations may be important in understanding how MBNL1 binds to RNA CUG repeats. Hence we provide new insight into the dynamic roles of RNA loops and their contributions to presently incurable diseases.
Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are ... more Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against (1)H-(1)H NOESY distances and (3) J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.
The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relati... more The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relative to the Watson-Crick GC pair, provides an expanded alphabet for understanding interactions that shape nucleic acid structure. Here, thermodynamic stabilities of tandem GA pairs flanked by iGiC pairs are reported along with the NMR structures of the RNA self-complementary duplexes (GCiGGAiCGCA) 2 and (GGiCGAiGCCA) 2 . A sheared GA pairing forms in (GCiGGAiCGCA) 2 , and an imino GA pairing forms in (GGiCGAiGCCA) 2 . The structures contrast with the formation of tandem imino and sheared GA pairs flanked by GC pairs in the RNA self-complementary duplexes (GCGGACGC) 2 and (GGCGAGCC) 2 , respectively. In both iGiC duplexes, Watson-Crick-like hydrogen bonds are formed between iG and iC, and iGiC substitutions result in less favorable loop stability. The results provide benchmarks for testing computations of molecular interactions that shape RNA three-dimensional structure.
Proceedings of the National Academy of Sciences, 2008
FRET between the zinc porphyrin (ZnP) chromophore in zincsubstituted cytochrome c (Zn-cyt c) and ... more FRET between the zinc porphyrin (ZnP) chromophore in zincsubstituted cytochrome c (Zn-cyt c) and an Alexa Fluor dye attached to specific surface sites was used to characterize Zn-cyt c unfolding. The use of ZnP as a fluorescent acceptor eliminates the need to doubly label the protein with exogenous dyes to perform FRET experiments in which both donor and acceptor fluorescence is monitored. The requirement for attachment of only one dye also minimizes perturbation to the protein. This sensitive technique allowed for the determination of distances between the label placed at six different sites and ZnP through a range of denaturant concentrations. Fitting of the data to a three-state model provides distances in the unfolding intermediate. The use of ZnP as a fluorescent acceptor of energy in FRET has a significant potential for application to a range of other systems including heme-binding proteins and proteins to which a covalently attached heme tag may be added.
RNA is an immensely important target for small molecule therapeutics or chemical probes of functi... more RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.
One class of functionally important RNA is repeating transcripts that cause disease through vario... more One class of functionally important RNA is repeating transcripts that cause disease through various mechanisms. For example, expanded CAG repeats can cause Huntington's and other disease through translation of toxic proteins. Herein, a crystal structure of r[5'UUGGGC(CAG)3GUCC]2, a model of CAG expanded transcripts, refined to 1.65 Å resolution is disclosed that shows both anti-anti and syn-anti orientations for 1 × 1 nucleotide AA internal loops. Molecular dynamics (MD) simulations using AMBER force field in explicit solvent were run for over 500 ns on the model systems r(5'GCGCAGCGC)2 (MS1) and r(5'CCGCAGCGG)2 (MS2). In these MD simulations, both anti-anti and syn-anti AA base pairs appear to be stable. While anti-anti AA base pairs were dynamic and sampled multiple anti-anti conformations, no syn-anti ↔ anti-anti transformations were observed. Umbrella sampling simulations were run on MS2, and a 2D free energy surface was created to extract transformation pathways. In addition, an explicit solvent MD simulation over 800 ns was run on r[5'GGGC(CAG)3GUCC]2, which closely represents the refined crystal structure. One of the terminal AA base pairs (syn-anti conformation), transformed to anti-anti conformation. The pathway followed in this transformation was the one predicted by umbrella sampling simulations. Further analysis showed a binding pocket near AA base pairs in syn-anti conformations. Computational results combined with the refined crystal structure show that global minimum conformation of 1 × 1 nucleotide AA internal loops in r(CAG) repeats is anti-anti but can adopt syn-anti depending on the environment. These results are important to understand RNA dynamic-function relationships and to develop small molecules that target RNA dynamic ensembles.
All-atom force fields are important for predicting thermodynamic, structural, and dynamic propert... more All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG) 2 , r(GGCC) 2 , r(GCGC) 2 , and r(CGCG) 2 are replaced by isocytidineÀisoguanosine (iCiG) pairs. Agreement with experiment was improved when ε/ζ, α/γ, β, and χ torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG f iCiCiGiG, GGCC f iGiGiCiC, GCGC f iGiCiGiC, and CGCG f iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC f iGiGiCiC and GCGC f iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG) 2 retained A-form conformation for g82% of the time. This is consistent with NMR spectra of r(iGiGiCiC) 2 , which reveal an A-form conformation. In MD simulations, r(CCGG) 2 retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods.
Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loop... more Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair and finally by a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and basephosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 µs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by MM-PBSA free energy calculations, locally enhanced sampling (LES) runs, targeted MD (TMD) and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and netneutralizing Na + vs. excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the Corresponding authors: Kamila Réblová, . SUPPORTING INFORMATION Detailed description of standard MD simulations run with parm99, description of standard MD simulations run with parmbsc0, and description of LES simulations of ribosomal H40. Simulated annealing protocol used for minimization in the NEB trials (in ). Calculation of free energy basin around H40 X-ray structure based on parm99 and parmbsc0 simulations using potential of mean force coordinate method. . This material is available free of charge via the Internet at major groove, loss of base-phosphate interaction and other instabilities. The alternative X-ray geometry even undergoes conformational transition towards the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.
A reparameterization of the torsional parameters for the glycosidic dihedral angle, , for the AMB... more A reparameterization of the torsional parameters for the glycosidic dihedral angle, , for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99 . Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99 force fields are compared with NMR results. For each nucleoside and force field, 10 individual molecular dynamics simulations of 30 ns each were run. For cytidine with AMBER99 force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) 1 H, steady-state 1D 1 H nuclear Overhauser effect (NOE), and transient 1D 1 H NOE, was used to determine the sugar puckering and preferred base orientation with respect to the ribose of cytidine and uridine. The AMBER99 force field overestimates the population of syn conformations of the base orientation and of C2′-endo sugar puckering of the pyrimidines, while the AMBER99 force field's predictions are more consistent with NMR results. Moreover, the AMBER99 force field prefers high anti conformations with glycosidic dihedral angles around 310°for the base orientation of purines. The AMBER99 force field prefers anti conformations around 185°, which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99 force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.
Locked Nucleic Acids (LNA) considerably enhance the thermodynamic stability of DNA and RNA duplex... more Locked Nucleic Acids (LNA) considerably enhance the thermodynamic stability of DNA and RNA duplexes. We report the thermodynamic stabilities of LNA-2'-O-methyl RNA/RNA duplexes designed to provide insight into the contributions of stacking and hydrogen bonding interactions to the enhanced stability. The results show that hydrogen bonding of LNA nucleotides is similar to that of 2'-O-methyl RNA nucleotides, whereas the 3'-stacking interactions are on average about 0.7 kcal/mol more favorable at 37 °C than for 2'-O-methyl or RNA nucleotides. Moreover, NMR spectra suggest helical pre-organization of the single stranded tetramer, C L A M A L U M , probably due to restriction of some torsion angles. Thus enhanced stacking interactions and helical pre-organization of single stranded oligonucleotides contribute to the extraordinary stabilization of duplexes by LNA nucleotides.
Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pai... more Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pairwise hydrogen-bonding models. Such effects include the stabilities of isoguanosine-isocytidine (iG-iC) base pairs and of various 2 x 2 nucleotide internal loops. Presumably, these results can be explained by base stacking effects, which can be partitioned into Coulombic and overlap effects. We review experimental measurements that provide benchmarks for testing the approximations and theories used for modeling nucleic acids. Quantitative agreement between experiment and theory will indicate understanding of the interactions determining RNA stability and structure.
Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expans... more Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disorder caused by a r(CCUG) expansion (r(CCUG)(exp)) that folds into an extended hairpin with periodically repeating 2×2 nucleotide internal loops (5'CCUG/3'GUCC). We designed multivalent compounds that improve DM2-associated defects using information about RNA-small molecule interactions. We also report the first crystal structure of r(CCUG) repeats refined to 2.35 Å. Structural analysis of the three 5'CCUG/3'GUCC repeat internal loops (L) reveals that the CU pairs in L1 are each stabilized by one hydrogen bond and a water-mediated hydrogen bond, while CU pairs in L2 and L3 are stabilized by two hydrogen bonds. Molecular dynamics (MD) simulations reveal that the CU pairs are dynamic and stabilized by Na(+) and water molecules. MD simulations of the binding of the small molecule to r(CCUG) repeats reveal that the lowest free energy binding mode occurs via the major groove, in which one C residue is unstacked and the cross-strand…
Many transcriptional activators act at a distance from core promoter elements and work by recruit... more Many transcriptional activators act at a distance from core promoter elements and work by recruiting RNA polymerase through protein-protein interactions. We show here how the prokaryotic regulatory protein CueR both represses and activates transcription by differentially modulating local DNA structure within the promoter. Structural studies reveal that the repressor state slightly bends the promoter DNA, precluding optimal RNA polymerase-promoter recognition. Upon binding a metal ion in the allosteric site, CueR switches into an activator conformation. It maintains all protein-DNA contacts but introduces torsional stresses that kink and undertwist the promoter, stabilizing an A-form DNA-like conformation. These factors switch on and off transcription by exerting dynamic control of DNA stereochemistry, reshaping the core promoter and making it a better or worse substrate for polymerase.
Spinocerebellar ataxia type 10 (SCA10) is caused by a pentanucleotide repeat expansion of r(AUUCU... more Spinocerebellar ataxia type 10 (SCA10) is caused by a pentanucleotide repeat expansion of r(AUUCU) within intron 9 of the ATXN10 pre-mRNA. The RNA causes disease by a gain-of-function mechanism in which it inactivates proteins involved in RNA biogenesis. Spectroscopic studies showed that r(AUUCU) repeats form a hairpin structure; however, there were no high-resolution structural models prior to this work. Herein, we report the first crystal structure of model r(AUUCU) repeats refined to 2.8 Å and analysis of the structure via molecular dynamics simulations. The r(AUUCU) tracts adopt an overall A-form geometry in which 3 × 3 nucleotide (5')UCU(3')/(3')UCU(5') internal loops are closed by AU pairs. Helical parameters of the refined structure as well as the corresponding electron density map on the crystallographic model reflect dynamic features of the internal loop. The computational analyses captured dynamic motion of the loop closing pairs, which can form single-stra...
Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are ... more Molecular dynamics (MD) simulations for RNA tetramers r(AAAA), r(CAAU), r(GACC), and r(UUUU) are benchmarked against (1)H-(1)H NOESY distances and (3) J scalar couplings to test effects of RNA torsion parametrizations. Four different starting structures were used for r(AAAA), r(CAAU), and r(GACC), while five starting structures were used for r(UUUU). On the basis of X-ray structures, criteria are reported for quantifying stacking. The force fields, AMBER ff99, parmbsc0, parm99χ_Yil, ff10, and parmTor, all predict experimentally unobserved stacks and intercalations, e.g., base 1 stacked between bases 3 and 4, and incorrect χ, ϵ, and sugar pucker populations. The intercalated structures are particularly stable, often lasting several microseconds. Parmbsc0, parm99χ_Yil, and ff10 give similar agreement with NMR, but the best agreement is only 46%. Experimentally unobserved intercalations typically are associated with reduced solvent accessible surface area along with amino and hydroxyl hydrogen bonds to phosphate nonbridging oxygens. Results from an extensive set of MD simulations suggest that recent force field parametrizations improve predictions, but further improvements are necessary to provide reasonable agreement with NMR. In particular, intramolecular stacking and hydrogen bonding interactions may not be well balanced with the TIP3P water model. NMR data and the scoring method presented here provide rigorous benchmarks for future changes in force fields and MD methods.
Uploads
Papers by Ilyas Yildirim