1
|
Walker-Gibbons R, Zhu X, Behjatian A, Bennett TJD, Krishnan M. Sensing the structural and conformational properties of single-stranded nucleic acids using electrometry and molecular simulations. Sci Rep 2024; 14:20582. [PMID: 39232063 PMCID: PMC11375218 DOI: 10.1038/s41598-024-70641-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Accepted: 08/20/2024] [Indexed: 09/06/2024] Open
Abstract
Inferring the 3D structure and conformation of disordered biomolecules, e.g., single stranded nucleic acids (ssNAs), remains challenging due to their conformational heterogeneity in solution. Here, we use escape-time electrometry (ETe) to measure with sub elementary-charge precision the effective electrical charge in solution of short to medium chain length ssNAs in the range of 5-60 bases. We compare measurements of molecular effective charge with theoretically calculated values for simulated molecular conformations obtained from Molecular Dynamics simulations using a variety of forcefield descriptions. We demonstrate that the measured effective charge captures subtle differences in molecular structure in various nucleic acid homopolymers of identical length, and also that the experimental measurements can find agreement with computed values derived from coarse-grained molecular structure descriptions such as oxDNA, as well next generation ssNA force fields. We further show that comparing the measured effective charge with calculations for a rigid, charged rod-the simplest model of a nucleic acid-yields estimates of molecular structural dimensions such as linear charge spacings that capture molecular structural trends observed using high resolution structural analysis methods such as X-ray scattering. By sensitively probing the effective charge of a molecule, electrometry provides a powerful dimension supporting inferences of molecular structural and conformational properties, as well as the validation of biomolecular structural models. The overall approach holds promise for a high throughput, microscopy-based biomolecular analytical approach offering rapid screening and inference of molecular 3D conformation, and operating at the single molecule level in solution.
Collapse
Affiliation(s)
- Rowan Walker-Gibbons
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK
| | - Xin Zhu
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK
| | - Ali Behjatian
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK
| | - Timothy J D Bennett
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK
| | - Madhavi Krishnan
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ, UK.
- The Kavli Institute for Nanoscience Discovery, Sherrington Road, Oxford, OX1 3QU, UK.
| |
Collapse
|
2
|
Mondal B, Chakraborty D, Hori N, Nguyen HT, Thirumalai D. Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences. J Chem Theory Comput 2024; 20:2934-2946. [PMID: 38498914 DOI: 10.1021/acs.jctc.3c01193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Interplay between divalent cations (Mg2+ and Ca2+) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT30 and dA30, with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg2+ and Ca2+, we find that Mg2+ binds to the minor grooves and phosphate groups, whereas Ca2+ binds specifically to the minor groove. Both Mg2+ and Ca2+ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA30 conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT30 populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.
Collapse
Affiliation(s)
- Balaka Mondal
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - Debayan Chakraborty
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - Naoto Hori
- School of Pharmacy, University of Nottingham, Nottingham NG7 2RD, U.K
| | - Hung T Nguyen
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
| | - D Thirumalai
- Department of Chemistry, The University of Texas, Austin, Texas 78712, United States
- Department of Physics, The University of Texas, Austin, Texas 78712, United States
| |
Collapse
|
3
|
Tan LH, Kwoh CK, Mu Y. RmsdXNA: RMSD prediction of nucleic acid-ligand docking poses using machine-learning method. Brief Bioinform 2024; 25:bbae166. [PMID: 38695120 PMCID: PMC11063749 DOI: 10.1093/bib/bbae166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 03/15/2024] [Accepted: 03/19/2024] [Indexed: 05/04/2024] Open
Abstract
Small molecule drugs can be used to target nucleic acids (NA) to regulate biological processes. Computational modeling methods, such as molecular docking or scoring functions, are commonly employed to facilitate drug design. However, the accuracy of the scoring function in predicting the closest-to-native docking pose is often suboptimal. To overcome this problem, a machine learning model, RmsdXNA, was developed to predict the root-mean-square-deviation (RMSD) of ligand docking poses in NA complexes. The versatility of RmsdXNA has been demonstrated by its successful application to various complexes involving different types of NA receptors and ligands, including metal complexes and short peptides. The predicted RMSD by RmsdXNA was strongly correlated with the actual RMSD of the docked poses. RmsdXNA also outperformed the rDock scoring function in ranking and identifying closest-to-native docking poses across different structural groups and on the testing dataset. Using experimental validated results conducted on polyadenylated nuclear element for nuclear expression triplex, RmsdXNA demonstrated better screening power for the RNA-small molecule complex compared to rDock. Molecular dynamics simulations were subsequently employed to validate the binding of top-scoring ligand candidates selected by RmsdXNA and rDock on MALAT1. The results showed that RmsdXNA has a higher success rate in identifying promising ligands that can bind well to the receptor. The development of an accurate docking score for a NA-ligand complex can aid in drug discovery and development advancements. The code to use RmsdXNA is available at the GitHub repository https://github.com/laiheng001/RmsdXNA.
Collapse
Affiliation(s)
- Lai Heng Tan
- Interdisciplinary Graduate School, Nanyang Technological University, 61 Nanyang Drive, 637335 Singapore, Singapore
| | - Chee Keong Kwoh
- School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, 639798 Singapore, Singapore
| | - Yuguang Mu
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, 637551 Singapore, Singapore
| |
Collapse
|
4
|
Love O, Winkler L, Cheatham TE. van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD. J Chem Theory Comput 2024; 20:625-643. [PMID: 38157247 PMCID: PMC10809421 DOI: 10.1021/acs.jctc.3c01164] [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] [Received: 10/20/2023] [Revised: 12/11/2023] [Accepted: 12/14/2023] [Indexed: 01/03/2024]
Abstract
Molecular dynamics simulations can be used in combination with experimental techniques to uncover the intricacies of biomolecular structure, dynamics, and the resulting interactions. However, many noncanonical nucleic acid structures have proven to be challenging to replicate in accurate agreement with experimental data, often attributed to known force field deficiencies. A common force field criticism is the handling of van der Waals (vdW) parameters, which have not been updated since the regular use of Ewald's methods became routine. This work dives into the effects of minute vdW radii shifts on RNA tetranucleotide, B-DNA, and Z-DNA model systems described by commonly used Amber force fields. Using multidimensional replica exchange molecular dynamics (M-REMD), the GACC RNA tetranucleotide demonstrated changes in the structural distribution between the NMR minor and anomalous structure populations based on the O2' vdW radii scanning. However, no significant change in the NMR Major conformation population was observed. There were minimal changes in the B-DNA structure but there were more substantial improvements in Z-DNA structural descriptions, specifically with the Tumuc1 force field. This occurred with both LJbb vdW radii adjustments and incorporation of the CUFIX nonbonded parameter modifications. Though the limited vdW modifications tested did not provide a universal fix to the challenge of simulating the various known nucleic acid structures, they do provide direction and a greater understanding for future force field development efforts.
Collapse
Affiliation(s)
| | | | - Thomas E. Cheatham
- Department of Medicinal Chemistry,
College of Pharmacy, University of Utah, 2000 East 30 South Skaggs 306, Salt Lake City, Utah 84112, United States
| |
Collapse
|
5
|
Ali Z, Kaur S, Kukhta T, Abu-Saleh AAAA, Jhunjhunwala A, Mitra A, Trant JF, Sharma P. Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA. J Phys Chem B 2023. [PMID: 37369074 DOI: 10.1021/acs.jpcb.3c02268] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.
Collapse
Affiliation(s)
- Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Sarabjeet Kaur
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
- Surface Chemistry and Catalysis: Characterisation and Application Team (COK-KAT), Leuven (Arenberg) Celestijnenlaan 200f─Box 2461, 3001 Leuven, Belgium
| | - Teagan Kukhta
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
| | - Abd Al-Aziz A Abu-Saleh
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
- Binary Star Research Services, LaSalle, Ontario N9J 3X8, Canada
| | - Ayush Jhunjhunwala
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - John F Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
- Binary Star Research Services, LaSalle, Ontario N9J 3X8, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
| |
Collapse
|
6
|
Li J, Chen SJ. RNAJP: enhanced RNA 3D structure predictions with non-canonical interactions and global topology sampling. Nucleic Acids Res 2023; 51:3341-3356. [PMID: 36864729 PMCID: PMC10123122 DOI: 10.1093/nar/gkad122] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 01/14/2023] [Accepted: 02/25/2023] [Indexed: 03/04/2023] Open
Abstract
RNA 3D structures are critical for understanding their functions. However, only a limited number of RNA structures have been experimentally solved, so computational prediction methods are highly desirable. Nevertheless, accurate prediction of RNA 3D structures, especially those containing multiway junctions, remains a significant challenge, mainly due to the complicated non-canonical base pairing and stacking interactions in the junction loops and the possible long-range interactions between loop structures. Here we present RNAJP ('RNA Junction Prediction'), a nucleotide- and helix-level coarse-grained model for the prediction of RNA 3D structures, particularly junction structures, from a given 2D structure. Through global sampling of the 3D arrangements of the helices in junctions using molecular dynamics simulations and in explicit consideration of non-canonical base pairing and base stacking interactions as well as long-range loop-loop interactions, the model can provide significantly improved predictions for multibranched junction structures than existing methods. Moreover, integrated with additional restraints from experiments, such as junction topology and long-range interactions, the model may serve as a useful structure generator for various applications.
Collapse
Affiliation(s)
- Jun Li
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, Columbia, MO 65211, USA
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry and Institute for Data Science and Informatics, University of Missouri, Columbia, MO 65211, USA
| |
Collapse
|
7
|
Liebl K, Zacharias M. Toward Force Fields with Improved Base Stacking Descriptions. J Chem Theory Comput 2023; 19:1529-1536. [PMID: 36795949 DOI: 10.1021/acs.jctc.2c01121] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.
Collapse
Affiliation(s)
- Korbinian Liebl
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Garching 85748, Germany
| |
Collapse
|
8
|
Shin JH, Bonilla SL, Denny SK, Greenleaf WJ, Herschlag D. Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements. Proc Natl Acad Sci U S A 2023; 120:e2220485120. [PMID: 36897989 PMCID: PMC10243134 DOI: 10.1073/pnas.2220485120] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Accepted: 02/01/2023] [Indexed: 03/12/2023] Open
Abstract
Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.
Collapse
Affiliation(s)
- John H. Shin
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA94305
- Department of Chemical Engineering, Stanford University, Stanford, CA94305
| | - Steve L. Bonilla
- Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, CO80045
| | - Sarah K. Denny
- Department of Genetics, Stanford University School of Medicine, Stanford, CA94305
- Scribe Therapeutics, Alameda, CA94501
| | - William J. Greenleaf
- Department of Genetics, Stanford University School of Medicine, Stanford, CA94305
- Department of Applied Physics, Stanford University, Stanford, CA94305
- Chan Zuckerberg Biohub, San Francisco, CA94158
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University School of Medicine, Stanford, CA94305
- Department of Chemical Engineering, Stanford University, Stanford, CA94305
- ChEM-H Institute, Stanford University, Stanford, CA94305
| |
Collapse
|
9
|
Ramasanoff RR, Sokolov PA. The binding model of adenosine-specific DNA aptamer: Umbrella sampling study. J Mol Graph Model 2023; 118:108338. [PMID: 36201878 DOI: 10.1016/j.jmgm.2022.108338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 09/16/2022] [Accepted: 09/17/2022] [Indexed: 11/25/2022]
Abstract
We report a novel model of the selective binding mechanism of adenosine-specific DNA aptamer. Our theoretical investigations of AMP (Adenosine monophosphate) dissociation from aptamer-AMP complexes reveals new details of aptamer molecular specificity and stabilisation factors. Umbrella sampling MD calculations using parmbsc1 force field shows that the disordered structure of the internal loop of the unbound aptamer hairpin has a characteristic packing of guanines, which prevents barrier-free penetration of ligands into the site cavity. Also, this disordered structure of the unbound aptamer has a network of hydrogen bonds stabilising the cavity near the target guanines within the binding sites during the whole binding process. We suggested that the first AMP molecule binds to the disordered structure of the site closest to the aptamer hairpin stem and spends some free energy on ordering of the internal loop. Then the second AMP molecule binds to the ordered site closest to the aptamer hairpin loop with a lower energy gain. As a result, the induced-fit binding model is the most applicable for this aptamer and does not contradict the modern experimental NMR and calorimetry data.
Collapse
Affiliation(s)
- Ruslan R Ramasanoff
- Sevastopol State University, Universitetskaya 33, 299053, Sevastopol, Russia.
| | - Petr A Sokolov
- Sevastopol State University, Universitetskaya 33, 299053, Sevastopol, Russia; Saint Petersburg State University, Universitetskaya Nab. 7/9, 199034, Saint Petersburg, Russia
| |
Collapse
|
10
|
Nosaki S, Mitsuda N, Sakamoto S, Kusubayashi K, Yamagami A, Xu Y, Bui TBC, Terada T, Miura K, Nakano T, Tanokura M, Miyakawa T. Brassinosteroid-induced gene repression requires specific and tight promoter binding of BIL1/BZR1 via DNA shape readout. NATURE PLANTS 2022; 8:1440-1452. [PMID: 36522451 DOI: 10.1038/s41477-022-01289-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Accepted: 10/26/2022] [Indexed: 05/12/2023]
Abstract
BRZ-INSENSITIVE-LONG 1 (BIL1)/BRASSINAZOLE-RESISTANT 1 (BZR1) and its homologues are plant-specific transcription factors that convert the signalling of the phytohormones brassinosteroids (BRs) to transcriptional responses, thus controlling various physiological processes in plants. Although BIL1/BZR1 upregulates some BR-responsive genes and downregulates others, the molecular mechanism underlying the dual roles of BIL1/BZR1 is still poorly understood. Here we show that BR-responsive transcriptional repression by BIL1/BZR1 requires the tight binding of BIL1/BZR1 alone to the 10 bp elements of DNA fragments containing the known 6 bp core-binding motifs at the centre. Furthermore, biochemical and structural evidence demonstrates that the selectivity for two nucleobases flanking the core motifs is realized by the DNA shape readout of BIL1/BZR1 without direct recognition of the nucleobases. These results elucidate the molecular and structural basis of transcriptional repression by BIL1/BZR1 and contribute to further understanding of the dual roles of BIL1/BZR1 in BR-responsive gene regulation.
Collapse
Affiliation(s)
- Shohei Nosaki
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
- Faculty of Life and Environmental Sciences, University of Tsukuba, Tsukuba, Ibaraki, Japan
- Tsukuba Plant-Innovation Research Center (T-PIRC), University of Tsukuba, Tsukuba, Ibaraki, Japan
| | - Nobutaka Mitsuda
- Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
- Global Zero Emission Research Center, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| | - Shingo Sakamoto
- Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
- Global Zero Emission Research Center, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| | - Kazuki Kusubayashi
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Ayumi Yamagami
- Graduate School of Biostudies, Kyoto University, Sakyo-ku, Kyoto, Japan
- Gene Discovery Research Group, RIKEN CSRS, Wako, Saitama, Japan
| | - Yuqun Xu
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Thi Bao Chau Bui
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Tohru Terada
- Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Kenji Miura
- Faculty of Life and Environmental Sciences, University of Tsukuba, Tsukuba, Ibaraki, Japan
- Tsukuba Plant-Innovation Research Center (T-PIRC), University of Tsukuba, Tsukuba, Ibaraki, Japan
| | - Takeshi Nakano
- Graduate School of Biostudies, Kyoto University, Sakyo-ku, Kyoto, Japan
- Gene Discovery Research Group, RIKEN CSRS, Wako, Saitama, Japan
| | - Masaru Tanokura
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan.
| | - Takuya Miyakawa
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo-ku, Tokyo, Japan.
- Graduate School of Biostudies, Kyoto University, Sakyo-ku, Kyoto, Japan.
| |
Collapse
|
11
|
Capobianco A, Landi A, Peluso A. Duplex DNA Retains the Conformational Features of Single Strands: Perspectives from MD Simulations and Quantum Chemical Computations. Int J Mol Sci 2022; 23:ijms232214452. [PMID: 36430930 PMCID: PMC9697240 DOI: 10.3390/ijms232214452] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2022] [Revised: 11/09/2022] [Accepted: 11/13/2022] [Indexed: 11/22/2022] Open
Abstract
Molecular dynamics simulations and geometry optimizations carried out at the quantum level as well as by quantum mechanical/molecular mechanics methods predict that short, single-stranded DNA oligonucleotides adopt conformations very similar to those observed in crystallographic double-stranded B-DNA, with rise coordinates close to ≈3.3 Å. In agreement with the experimental evidence, the computational results show that DNA single strands rich in adjacent purine nucleobases assume more regular arrangements than poly-thymine. The preliminary results suggest that single-stranded poly-cytosine DNA should also retain a substantial helical order in solution. A comparison of the structures of single and double helices confirms that the B-DNA motif is a favorable arrangement also for single strands. Indeed, the optimal geometry of the complementary single helices is changed to a very small extent in the formation of the duplex.
Collapse
|
12
|
Xiao L, Fang L, Kool ET. Acylation probing of "generic" RNA libraries reveals critical influence of loop constraints on reactivity. Cell Chem Biol 2022; 29:1341-1352.e8. [PMID: 35662395 PMCID: PMC9391288 DOI: 10.1016/j.chembiol.2022.05.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Revised: 03/16/2022] [Accepted: 05/12/2022] [Indexed: 01/08/2023]
Abstract
The reactivity of RNA 2'-OH acylation is broadly useful both in probing structure and in preparing conjugates. To date, this reactivity has been analyzed in limited sets of biological RNA sequences, leaving open questions of how reactivity varies inherently without regard to sequence in structured contexts. We constructed and probed "generic" structured RNA libraries using homogeneous loop sequences, employing deep sequencing to carry out a systematic survey of reactivity. We find a wide range of RNA reactivities among single-stranded sequences, with nearest neighbors playing substantial roles. Remarkably, certain small loops are found to be far more reactive on average (up to 4,000-fold) than single-stranded RNAs, due to conformational constraints that enhance reactivity. Among loops, we observe large variations in reactivity based on size, type, and position. The results lend insights into RNA designs for achieving high-efficiency local conjugation and provide new opportunities to refine structure analysis.
Collapse
Affiliation(s)
- Lu Xiao
- Department of Chemistry and ChEM-H Institute, Stanford University, Stanford, CA 94305, USA
| | - Linglan Fang
- Department of Chemistry and ChEM-H Institute, Stanford University, Stanford, CA 94305, USA
| | - Eric T Kool
- Department of Chemistry and ChEM-H Institute, Stanford University, Stanford, CA 94305, USA.
| |
Collapse
|
13
|
Taghavi A, Riveros I, Wales DJ, Yildirim I. Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations. J Chem Theory Comput 2022; 18:3637-3653. [PMID: 35652685 DOI: 10.1021/acs.jctc.2c00178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.
Collapse
Affiliation(s)
- Amirhossein Taghavi
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States.,Department of Chemistry, Scripps Research Institute Florida, Jupiter, Florida 33458, United States
| | - Ivan Riveros
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States
| | - David J Wales
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - Ilyas Yildirim
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States
| |
Collapse
|
14
|
Rieger M, Zacharias M. Nearest-Neighbor dsDNA Stability Analysis Using Alchemical Free-Energy Simulations. J Phys Chem B 2022; 126:3640-3647. [PMID: 35549273 DOI: 10.1021/acs.jpcb.2c01138] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The thermodynamic stability of double-stranded (ds)DNA depends on its sequence. It is influenced by the base pairing and stacking with neighboring bases along DNA molecules. Semiempirical schemes are available that allow us to predict the thermodynamic stability of DNA sequences based on empirically derived nearest-neighbor contributions of base pairs formed in the context of all possible nearest-neighbor base pairs. Current molecular dynamics (MD) simulations allow one to simulate the dynamics of DNA molecules in good agreement with experimentally obtained structures and available data on conformational flexibility. However, the suitability of current force field methods to reproduce dsDNA stability and its sequence dependence has been much less well tested. We have employed alchemical free-energy simulations of whole base pair transversions in dsDNA and in unbound single-stranded partner molecules. Such transversions change the sequence context but not the nucleotide content or base pairing in dsDNA and allow a direct comparison with the empirical nearest-neighbor dsDNA stability model. For the alchemical free-energy changes in the unbound single-stranded (ss)DNA partner molecules, we tested different setups assuming either complete unstacking or unrestrained simulations with partial stacking in the unbound ssDNA. The free-energy simulations predicted nearest-neighbor effects of similar magnitude, as observed experimentally but showed overall limited correlation with experimental data. An inaccurate description of stacking interactions and other possible reasons such as the neglect of electronic polarization effects are discussed. The results indicate the need to improve the realistic description of stacking interactions in current molecular mechanic force fields.
Collapse
Affiliation(s)
- Manuel Rieger
- Physics Department and Center of Protein Assemblies, Technical University of Munich, 85748 Garching, Germany
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, 85748 Garching, Germany
| |
Collapse
|
15
|
Jash B, Kool ET. Conjugation of RNA via 2'-OH acylation: Mechanisms determining nucleotide reactivity. Chem Commun (Camb) 2022; 58:3693-3696. [PMID: 35226025 PMCID: PMC9211027 DOI: 10.1039/d2cc00660j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The acylation reactivity of RNA 2'-OH groups has proven broadly useful for labeling and mapping RNA. Here we perform kinetics studies to test the mechanisms governing this reaction, and we find strong steric and inductive effects modulating reactivity. The results shed light on new strategies for improved conjugation and mapping.
Collapse
Affiliation(s)
- Biswarup Jash
- Department of Chemistry and ChEM-H Institute, Stanford University, Stanford, CA 94305, USA.
| | - Eric T Kool
- Department of Chemistry and ChEM-H Institute, Stanford University, Stanford, CA 94305, USA.
| |
Collapse
|
16
|
Abstract
Modifications are present on many classes of RNA, including tRNA, rRNA, and mRNA. These modifications modulate diverse biological processes such as genetic recoding and mRNA export and folding. In addition, modifications can be introduced to RNA molecules using chemical probing strategies that reveal RNA structure and dynamics. Many methods exist to detect RNA modifications by short-read sequencing; however, limitations on read length inherent to short-read-based methods dissociate modifications from their native context, preventing single-molecule modification analysis. Here, we demonstrate direct RNA nanopore sequencing to detect endogenous and exogenous RNA modifications on long RNAs at the single-molecule level. We detect endogenous 2'-O-methyl and base modifications across E. coli and S. cerevisiae ribosomal RNAs as shifts in current signal and dwell times distally through interactions with the helicase motor protein. We further use the 2'-hydroxyl reactive SHAPE reagent acetylimidazole to probe RNA structure at the single-molecule level with readout by direct nanopore sequencing. Stephenson et al. employ direct RNA nanopore sequencing to detect endogenous and exogenous modifications on single RNA molecules. The authors demonstrate detection of endogenous 2'-O-methylation (Nm) on native ribosomal RNAs, confirming known modification patterns. They describe the development of nanoSHAPE, a method that involves exogenously labeling RNA with a small-adduct-generating chemical probe that can reveal RNA structure using long-read sequencing.
Collapse
|
17
|
Sengar A, Ouldridge TE, Henrich O, Rovigatti L, Šulc P. A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results. Front Mol Biosci 2021; 8:693710. [PMID: 34235181 PMCID: PMC8256390 DOI: 10.3389/fmolb.2021.693710] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2021] [Accepted: 05/27/2021] [Indexed: 11/13/2022] Open
Abstract
The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.
Collapse
Affiliation(s)
- A. Sengar
- Centre for Synthetic Biology, Department of Bioengineering, Imperial College London, London, United Kingdom
| | - T. E. Ouldridge
- Centre for Synthetic Biology, Department of Bioengineering, Imperial College London, London, United Kingdom
| | - O. Henrich
- Department of Physics, SUPA, University of Strathclyde, Glasgow, United Kingdom
| | - L. Rovigatti
- Department of Physics, Sapienza University of Rome, Rome, Italy
- CNR Institute of Complex Systems, Sapienza University of Rome, Rome, Italy
| | - P. Šulc
- Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, AZ, United States
- School of Molecular Sciences, Arizona State University, Tempe, AZ, United States
| |
Collapse
|
18
|
Beiranvand N, Freindorf M, Kraka E. Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study. Molecules 2021; 26:2268. [PMID: 33919989 PMCID: PMC8071019 DOI: 10.3390/molecules26082268] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2021] [Revised: 04/08/2021] [Accepted: 04/09/2021] [Indexed: 12/13/2022] Open
Abstract
In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.
Collapse
Affiliation(s)
| | | | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Ave, Dallas, TX 75275-0314, USA; (N.B.); (M.F.)
| |
Collapse
|
19
|
Abstract
![]()
Coacervates are a
type of liquid–liquid phase separated
(LLPS) droplets that can serve as models of membraneless organelles
(MLOs) in living cells. Peptide–nucleotide coacervates have
been widely used to mimic properties of ribonucleoprotein (RNP) granules,
but the thermal stability and the role of base stacking is still poorly
understood. Here, we report a systematic investigation of coacervates
formed by five different nucleoside triphosphates (NTPs) with poly-l-lysine and poly-l-arginine as a function of temperature.
All studied combinations exhibit an upper critical solution temperature
(UCST), and a temperature-dependent critical salt concentration, originating
from a significant nonelectrostatic contribution to the mixing free
energy. Both the enthalpic and entropic parts of this nonelectrostatic
interaction decrease in the order G/A/U/C/T, in accordance with nucleobase
stacking free energies. Partitioning of two dyes proves that the local
hydrophobicity inside the peptide–nucleotide coacervates is
different for every nucleoside triphosphate. We derive a simple relation
between the temperature and salt concentration at the critical point
based on a mean-field model of phase separation. Finally, when different
NTPs are mixed with one common oppositely charged peptide, hybrid
coacervates were formed, characterized by a single intermediate UCST
and critical salt concentration. NTPs with lower critical salt concentrations
can remain condensed in mixed coacervates far beyond their original
critical salt concentration. Our results show that NTP-based coacervates
have a strong temperature sensitivity due to base stacking interactions
and that mixing NTPs can significantly influence the stability of
condensates and, by extension, their bioavailability.
Collapse
Affiliation(s)
- Tiemei Lu
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Karina K Nakashima
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Evan Spruijt
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| |
Collapse
|
20
|
Levintov L, Paul S, Vashisth H. Reaction Coordinate and Thermodynamics of Base Flipping in RNA. J Chem Theory Comput 2021; 17:1914-1921. [PMID: 33594886 DOI: 10.1021/acs.jctc.0c01199] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.
Collapse
Affiliation(s)
- Lev Levintov
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, United States
| | - Sanjib Paul
- Department of Chemistry, New York University, New York 10003, New York, United States
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, United States
| |
Collapse
|
21
|
Zacharias M. Base-Pairing and Base-Stacking Contributions to Double-Stranded DNA Formation. J Phys Chem B 2020; 124:10345-10352. [PMID: 33156627 DOI: 10.1021/acs.jpcb.0c07670] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Double-stranded (ds)DNA formation and dissociation are of fundamental biological importance. The negative DNA charge influences the dsDNA stability. However, the base pairing and the stacking between neighboring bases are responsible for the sequence-dependent stability of dsDNA. The stability of a dsDNA molecule can be estimated from empirical nearest-neighbor models based on contributions assigned to base-pair steps along the DNA and additional parameters because of DNA termini. In efforts to separate contributions, it has been concluded that base stacking dominates dsDNA stability, whereas base pairing contributes negligibly. Using a different model for dsDNA formation, we reanalyze dsDNA stability contributions and conclude that base stacking contributes already at the level of separate ssDNAs but that pairing contributions drive the dsDNA formation. The theoretical model also predicts that stability contributions of base-pair steps that contain only guanine/cytosine, mixed steps, and steps with only adenine/thymine follow the order 6:5:4, respectively, as expected based on the formed hydrogen bonds. The model is fully consistent with the available stacking data and the nearest-neighbor dsDNA parameters. It allows assigning a narrowly distributed value for the effective free energy contribution per formed hydrogen bond during dsDNA formation of -0.72 kcal·mol-1 based entirely on the experimental data.
Collapse
Affiliation(s)
- Martin Zacharias
- Physics Department T38, Technical University of Munich, 85748 Garching, Germany
| |
Collapse
|
22
|
Giraud T, Bouguet-Bonnet S, Marchal P, Pickaert G, Averlant-Petit MC, Stefan L. Improving and fine-tuning the properties of peptide-based hydrogels via incorporation of peptide nucleic acids. NANOSCALE 2020; 12:19905-19917. [PMID: 32985645 DOI: 10.1039/d0nr03483e] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Peptide self-assemblies have attracted intense research interest over the last few decades thanks to their implications in key biological processes (e.g., amyloid formation) and their use in biotechnological and (bio)material fields. In particular, peptide-based hydrogels have been highly considered as high potential supramolecular materials in the biomedical domain and open new horizons in terms of applications. To further understand their self-assembly mechanisms and to optimize their properties, several strategies have been proposed with the modification of the constituting amino acid chains via, per se, the introduction of d-amino acids, halogenated amino acids, pseudopeptide bonds, or other chemical moieties. In this context, we report herein on the incorporation of DNA-nucleobases into their peptide nucleic acid (PNA) forms to develop a new series of hybrid nucleopeptides. Thus, depending on the nature of the nucleobase (i.e., thymine, cytosine, adenine or guanine), the physicochemical and mechanical properties of the resulting hydrogels can be significantly improved and fine-tuned with, for instance, drastic enhancements of both the gel stiffness (up to 70-fold) and the gel resistance to external stress (up to 40-fold), and the generation of both thermo-reversible and uncommon red-edge excitation shift (REES) properties. To decipher the actual role of each PNA moiety in the self-assembly processes, the induced modifications from the molecular to the macroscopic scales are studied thanks to the multiscale approach based on a large panel of analytical techniques (i.e., rheology, NMR relaxometry, TEM, thioflavin T assays, FTIR, CD, fluorescence, NMR chemical shift index). Thus, such a strategy provides new opportunities to adapt and fit hydrogel properties to the intended ones and pushes back the limits of supramolecular materials.
Collapse
Affiliation(s)
- Tristan Giraud
- Université de Lorraine, CNRS, LCPM, F-54000 Nancy, France.
| | | | | | | | | | | |
Collapse
|
23
|
Tang TTL, Passmore LA. Recognition of Poly(A) RNA through Its Intrinsic Helical Structure. COLD SPRING HARBOR SYMPOSIA ON QUANTITATIVE BIOLOGY 2020; 84:21-30. [PMID: 32295929 PMCID: PMC7116106 DOI: 10.1101/sqb.2019.84.039818] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The polyadenosine (poly(A)) tail, which is found on the 3’ end of almost all eukaryotic messenger RNAs (mRNAs), plays an important role in the posttranscriptional regulation of gene expression. Shortening of the poly(A) tail, a process known as deadenylation, is thought to be the first and rate-limiting step of mRNA turnover. Deadenylation is performed by the Pan2–Pan3 and Ccr4–Not complexes that contain highly conserved exonuclease enzymes Pan2, and Ccr4 and Caf1, respectively. These complexes have been extensively studied, but the mechanisms of how the deadenylase enzymes recognize the poly(A) tail were poorly understood until recently. Here, we summarize recent work from our laboratory demonstrating that the highly conserved Pan2 exonuclease recognizes the poly(A) tail, not through adenine-specific functional groups, but through the conformation of poly(A) RNA. Our biochemical, biophysical, and structural investigations suggest that poly(A) forms an intrinsic base-stacked, single-stranded helical conformation that is recognized by Pan2, and that disruption of this structure inhibits both Pan2 and Caf1. This intrinsic structure has been shown to be important in poly(A) recognition in other biological processes, further underlining the importance of the unique conformation of poly(A).
Collapse
Affiliation(s)
- Terence T L Tang
- MRC Laboratory of Molecular Biology, Cambridge CB2 0QH, United Kingdom
| | - Lori A Passmore
- MRC Laboratory of Molecular Biology, Cambridge CB2 0QH, United Kingdom
| |
Collapse
|
24
|
Molecular Dynamics Simulations of DNA Adsorption on Graphene Oxide and Reduced Graphene Oxide-PEG-NH2 in the Presence of Mg2+ and Cl− ions. COATINGS 2020. [DOI: 10.3390/coatings10030289] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Graphene and its functionalised derivatives are transforming the development of biosensors that are capable of detecting nucleic acid hybridization. Using a Molecular Dynamics (MD) approach, we explored single-stranded or double-stranded deoxyribose nucleic acid (ssDNA or dsDNA) adsorption on two graphenic species: graphene oxide (GO) and reduced graphene oxide functionalized with aminated polyethylene glycol (rGO-PEG-NH2). Innovatively, we included chloride (Cl−) and magnesium (Mg2+) ions that influenced both the ssDNA and dsDNA adsorption on GO and rGO-PEG-NH2 surfaces. Unlike Cl−, divalent Mg2+ ions formed bridges between the GO surface and DNA molecules, promoting adsorption through electrostatic interactions. For rGO-PEG-NH2, the Mg2+ ions were repulsed from the graphenic surface. The subsequent ssDNA adsorption, mainly influenced by electrostatic forces and hydrogen bonds, could be supported by π–π stacking interactions that were absent in the case of dsDNA. We provide a novel insight for guiding biosensor development.
Collapse
|
25
|
Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction. Biochim Biophys Acta Gen Subj 2020; 1864:129600. [PMID: 32179130 DOI: 10.1016/j.bbagen.2020.129600] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2019] [Revised: 03/05/2020] [Accepted: 03/11/2020] [Indexed: 11/21/2022]
Abstract
BACKGROUND Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs. METHODS We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair. RESULTS We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear. CONCLUSIONS The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity. GENERAL SIGNIFICANCE We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.
Collapse
|
26
|
Naskar S, Saurabh S, Jang YH, Lansac Y, Maiti PK. Liquid crystal ordering of nucleic acids. SOFT MATTER 2020; 16:634-641. [PMID: 31840704 DOI: 10.1039/c9sm01816f] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Several analytical calculations and computer simulations propose that cylindrical monodispersive rods having an aspect ratio (ratio of length to diameter) greater than 4 can exhibit liquid crystal (LC) ordering. But, recent experiments demonstrated the signature of LC ordering in systems of 4- to 20-base pair (bp) long nucleic acids (NAs) that do not satisfy the shape anisotropy criterion. Mechanisms of end-to-end adhesion and stacking have been proposed to explain this phenomenon. In this study, using all-atom molecular dynamics (MD) simulation, we explicitly verify the end-to-end stacking of double-stranded RNA (dsRNA) and demonstrate the LC ordering at the microscopic level. Using umbrella sampling (US) calculation, we quantify the potential of mean force (PMF) between two dsRNAs for various reaction coordinates (RCs) and compare our results with previously reported PMFs for double-stranded DNA (dsDNA). The PMF profiles demonstrate the anisotropic nature of inter-NA interaction. We find that, like dsDNA, dsRNA also prefers to stack on top of each other while repelling sideways, leading to the formation of supra-molecular-columns that undergo LC ordering at high NA volume fraction (φ). We also demonstrate and quantify the nematic ordering of the RNAs using several hundred nanosecond-long MD simulations that remain almost invariant for different initial configurations and under different external physiological conditions.
Collapse
Affiliation(s)
- Supriyo Naskar
- Center for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India.
| | | | | | | | | |
Collapse
|
27
|
Ohmann A, Göpfrich K, Joshi H, Thompson RF, Sobota D, Ranson NA, Aksimentiev A, Keyser UF. Controlling aggregation of cholesterol-modified DNA nanostructures. Nucleic Acids Res 2019; 47:11441-11451. [PMID: 31642494 PMCID: PMC6868430 DOI: 10.1093/nar/gkz914] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Revised: 09/27/2019] [Accepted: 10/07/2019] [Indexed: 12/31/2022] Open
Abstract
DNA nanotechnology allows for the design of programmable DNA-built nanodevices which controllably interact with biological membranes and even mimic the function of natural membrane proteins. Hydrophobic modifications, covalently linked to the DNA, are essential for targeted interfacing of DNA nanostructures with lipid membranes. However, these hydrophobic tags typically induce undesired aggregation eliminating structural control, the primary advantage of DNA nanotechnology. Here, we study the aggregation of cholesterol-modified DNA nanostructures using a combined approach of non-denaturing polyacrylamide gel electrophoresis, dynamic light scattering, confocal microscopy and atomistic molecular dynamics simulations. We show that the aggregation of cholesterol-tagged ssDNA is sequence-dependent, while for assembled DNA constructs, the number and position of the cholesterol tags are the dominating factors. Molecular dynamics simulations of cholesterol-modified ssDNA reveal that the nucleotides wrap around the hydrophobic moiety, shielding it from the environment. Utilizing this behavior, we demonstrate experimentally that the aggregation of cholesterol-modified DNA nanostructures can be controlled by the length of ssDNA overhangs positioned adjacent to the cholesterol. Our easy-to-implement method for tuning cholesterol-mediated aggregation allows for increased control and a closer structure-function relationship of membrane-interfacing DNA constructs - a fundamental prerequisite for employing DNA nanodevices in research and biomedicine.
Collapse
Affiliation(s)
- Alexander Ohmann
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK
| | - Kerstin Göpfrich
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK
- Max Planck Institute for Medical Research, Department of Cellular Biophysics, Jahnstraße 29, 69120 Heidelberg, Germany
| | - Himanshu Joshi
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA
| | | | - Diana Sobota
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK
| | - Neil A Ranson
- Faculty of Biological Sciences, University of Leeds, Leeds LS2 9JT, UK
| | - Aleksei Aksimentiev
- Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA
| | - Ulrich F Keyser
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge CB3 0HE, UK
| |
Collapse
|
28
|
Rouquette M, Lepetre-Mouelhi S, Couvreur P. Adenosine and lipids: A forced marriage or a love match? Adv Drug Deliv Rev 2019; 151-152:233-244. [PMID: 30797954 DOI: 10.1016/j.addr.2019.02.005] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2019] [Revised: 02/12/2019] [Accepted: 02/13/2019] [Indexed: 12/21/2022]
Abstract
Adenosine is a fascinating compound, crucial in many biochemical processes: this ubiquitous nucleoside serves as an essential building block of RNA, is also a component of ATP and regulates numerous pathophysiological mechanisms via binding to four extracellular receptors. Due to its hydrophilic nature, it belongs to a different world than lipids, and has no affinity for them. Since the 1970's, however, new discoveries have emerged and prompted the scientific community to associate adenosine with the lipid family, especially via liposomal preparations and bioconjugation. This seems to be an arranged marriage, but could it turn into a true love match? This review considered all types of unions established between adenosine and lipids. Even though exciting supramolecular structures were observed with adenosine-lipid conjugates, as well as with liposomal preparations which resulted in promising pre-clinical results, the translation of these technologies to the clinic is still limited.
Collapse
|
29
|
Lee JY, Kim YJ, Lee C, Lee JG, Yagyu H, Tabata O, Kim DN. Investigating the sequence-dependent mechanical properties of DNA nicks for applications in twisted DNA nanostructure design. Nucleic Acids Res 2019; 47:93-102. [PMID: 30476210 PMCID: PMC6326809 DOI: 10.1093/nar/gky1189] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Accepted: 11/09/2018] [Indexed: 01/09/2023] Open
Abstract
DNA nick can be used as a design motif in programming the shape and reconfigurable deformation of synthetic DNA nanostructures, but its mechanical properties have rarely been systematically characterized at the level of base sequences. Here, we investigated sequence-dependent mechanical properties of DNA nicks through molecular dynamics simulation for a comprehensive set of distinct DNA oligomers constructed using all possible base-pair steps with and without a nick. We found that torsional rigidity was reduced by 28–82% at the nick depending on its sequence and location although bending and stretching rigidities remained similar to those of regular base-pair steps. No significant effect of a nick on mechanically coupled deformation such as the twist-stretch coupling was observed. These results suggest that the primary structural role of nick is the relaxation of torsional constraint by backbones known to be responsible for relatively high torsional rigidity of DNA. Moreover, we experimentally demonstrated the usefulness of quantified nick properties in self-assembling DNA nanostructure design by constructing twisted DNA origami structures to show that sequence design of nicks successfully controls the twist angle of structures. Our study illustrates the importance as well as the opportunities of considering sequence-dependent properties in structural DNA nanotechnology.
Collapse
Affiliation(s)
- Jae Young Lee
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Young-Joo Kim
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Chanseok Lee
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Jae Gyung Lee
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Hiromasa Yagyu
- Department of Mechanical Engineering, Kanto Gakuin University, Yokohama 236-8501, Japan
| | - Osamu Tabata
- Department of Micro Engineering, Kyoto University, Kyoto Daigaku-Katsura, Kyoto 615-8540, Japan
| | - Do-Nyun Kim
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea.,Institute of Advanced Machines and Design, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| |
Collapse
|
30
|
Carter CW, Wills PR. Hierarchical groove discrimination by Class I and II aminoacyl-tRNA synthetases reveals a palimpsest of the operational RNA code in the tRNA acceptor-stem bases. Nucleic Acids Res 2019; 46:9667-9683. [PMID: 30016476 PMCID: PMC6182185 DOI: 10.1093/nar/gky600] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Accepted: 07/12/2018] [Indexed: 01/01/2023] Open
Abstract
Class I and II aaRS recognition of opposite grooves was likely among the earliest determinants fixed in the tRNA acceptor stem bases. A new regression model identifies those determinants in bacterial tRNAs. Integral coefficients relate digital dependent to independent variables with perfect agreement between observed and calculated grooves for all twenty isoaccepting tRNAs. Recognition is mediated by the Discriminator base 73, the first base pair, and base 2 of the acceptor stem. Subsets of these coefficients also identically compute grooves recognized by smaller numbers of aaRS. Thus, the model is hierarchical, suggesting that new rules were added to pre-existing ones as new amino acids joined the coding alphabet. A thermodynamic rationale for the simplest model implies that Class-dependent aaRS secondary structures exploited differential tendencies of the acceptor stem to form the hairpin observed in Class I aaRS•tRNA complexes, enabling the earliest groove discrimination. Curiously, groove recognition also depends explicitly on the identity of base 2 in a manner consistent with the middle bases of the codon table, confirming a hidden ancestry of codon-anticodon pairing in the acceptor stem. That, and the lack of correlation with anticodon bases support prior productive coding interaction of tRNA minihelices with proto-mRNA.
Collapse
Affiliation(s)
- Charles W Carter
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599-7260, USA
| | - Peter R Wills
- Department of Physics, Centre for Computational Evolution, and Te Ao Marama Centre for Fundamental Enquiry, University of Auckland, PB 92109, Auckland 1142, New Zealand
| |
Collapse
|
31
|
Duchi M, O'Hagan MP, Kumar R, Bennie SJ, Galan MC, Curchod BFE, Oliver TAA. Exploring ultraviolet photoinduced charge-transfer dynamics in a model dinucleotide of guanine and thymine. Phys Chem Chem Phys 2019; 21:14407-14417. [PMID: 30869082 DOI: 10.1039/c8cp07864e] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An understanding of the initial photoexcited states of DNA is essential to unravelling deleterious photoinduced chemical reactions and the intrinsic ultrafast photoprotection of the genetic code for all life. In our combined experimental and theoretical study, we have elucidated the primary non-radiative relaxation dynamics of a model nucleotide of guanine and thymine (2'-deoxyguanosine 3'-monophosphate 5'-thymidine, d(GpT)) in buffered aqueous solution. Experimentally, we unequivocally demonstrate that the Franck-Condon excited states of d(GpT) are significantly delocalised across both nucleobases, and mediate d(G+pT-) exciplex product formation on an ultrafast (<350 fs) timescale. Theoretical studies show that the nature of the vertical excited states is very dependent on the specific geometry of the dinucleotide, and dictate the degree of delocalised, charge-transfer or localised character. Our mechanism for prompt exciplex formation involves a rapid change in electronic structure and includes a diabatic surface crossing very close to the Franck-Condon region mediating fast d(G+pT-) formation. Exciplexes are quickly converted back to neutral ground state molecules on a ∼10 ps timescale with a high quantum yield, ensuring the photostability of the nucleotide sequence.
Collapse
Affiliation(s)
- Marta Duchi
- School of Chemistry, Cantock's University of Bristol, Bristol, BS8 1TS, UK.
| | - Michael P O'Hagan
- School of Chemistry, Cantock's University of Bristol, Bristol, BS8 1TS, UK.
| | - Rhea Kumar
- School of Chemistry, Cantock's University of Bristol, Bristol, BS8 1TS, UK.
| | - Simon J Bennie
- School of Chemistry, Cantock's University of Bristol, Bristol, BS8 1TS, UK.
| | - M Carmen Galan
- School of Chemistry, Cantock's University of Bristol, Bristol, BS8 1TS, UK.
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham, DH1 3LE, UK
| | - Thomas A A Oliver
- School of Chemistry, Cantock's University of Bristol, Bristol, BS8 1TS, UK.
| |
Collapse
|
32
|
Carter CW, Wills PR. Class I and II aminoacyl-tRNA synthetase tRNA groove discrimination created the first synthetase-tRNA cognate pairs and was therefore essential to the origin of genetic coding. IUBMB Life 2019; 71:1088-1098. [PMID: 31190358 DOI: 10.1002/iub.2094] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Revised: 04/14/2019] [Accepted: 04/15/2019] [Indexed: 12/20/2022]
Abstract
The genetic code likely arose when a bidirectional gene replicating as a quasi-species began to produce ancestral aminoacyl-tRNA synthetases (aaRS) capable of distinguishing between two distinct sets of amino acids. The synthetase class division therefore necessarily implies a mechanism by which the two ancestral synthetases could also discriminate between two different kinds of tRNA substrates. We used regression methods to uncover the possible patterns of base sequences capable of such discrimination and find that they appear to be related to thermodynamic differences in the relative stabilities of a hairpin necessary for recognition of tRNA substrates by Class I aaRS. The thermodynamic differences appear to be exploited by secondary structural differences between models for the ancestral aaRS called synthetase Urzymes and reinforced by packing of aromatic amino acid side chains against the nonpolar face of the ribose of A76 if and only if the tRNA CCA sequence forms a hairpin. The patterns of bases 1, 2, and 73 and stabilization of the hairpin by structural complementarity with Class I, but not Class II, aaRS Urzymes appear to be necessary and sufficient to have enabled the generation of the first two aaRS-tRNA cognate pairs, and the launch of a rudimentary binary genetic coding related recognizably to contemporary cognate pairs. As a consequence, it seems likely that nonrandom aminoacylation of tRNAs preceded the advent of the tRNA anticodon stem-loop. Consistent with this suggestion, coding rules in the acceptor-stem bases also reveal a palimpsest of the codon-anticodon interaction, as previously proposed. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1088-1098, 2019.
Collapse
Affiliation(s)
- Charles W Carter
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Peter R Wills
- Department of Physics and Te Ao Marama Centre for Fundamental Inquiry, University of Auckland, Auckland, New Zealand
| |
Collapse
|
33
|
Kim HW, Rhee YM, Shin SK. Charge-dipole interactions in G-quadruplex thrombin-binding aptamer. Phys Chem Chem Phys 2019; 20:21068-21074. [PMID: 30074033 DOI: 10.1039/c8cp03050b] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
DNAs form various structures through hydrogen-bonding, base-stacking and electrostatic interactions. Although these noncovalent interactions are known to be cooperative in stabilizing a G-quadruplex (G4) structure of DNA, we find from all-atom molecular dynamics simulations that the electrostatic charge-dipole interaction is competitive with both hydrogen-bonding and base-stacking interactions. For the thrombin-binding aptamer (TBA) forming a chair-type antiparallel G4 structure, we have examined effects of an intercalating metal ion [K+, Sr2+, Mn+: an ion having a charge of n+ (n = 1-4) with the ionic radius of K+] on structural properties and noncovalent interactions. When K+ in the TBA·K+ complex is replaced with Sr2+, guanine dipoles in the two G-tetrads are realigned toward the central metal ion, thereby distorting the planar G4 geometry. Replacing K+ with Sr2+ significantly enhances the charge-dipole interaction but substantially reduces the number of hydrogen bonds in the G-tetrads. In the case of TBA·Mn+ complexes, as the charge n increases, the charge-dipole interaction increases but both of the hydrogen-bonding and base-stacking interactions decrease. These results suggest that the charge-dipole interaction realigning guanine dipoles in the G-tetrads is not cooperative but competitive with both hydrogen-bonding and base-stacking interactions favoring the planar G-tetrad geometry. Obviously, the charge state of an intercalating metal ion is as important as the ionic radius in forming a stable G4 structure. Thus, a delicate balance between these competing noncovalent interactions makes the chair-type antiparallel G4 structure of TBA selective for intercalating metal ions.
Collapse
Affiliation(s)
- Hyun Woo Kim
- Center for Molecular Modeling and Simulation, Korea Research Institute of Chemical Technology, Daejeon 34114, Korea
| | | | | |
Collapse
|
34
|
Baek K, Noblett AD, Ren P, Suggs LJ. Design and Characterization of Nucleopeptides for Hydrogel Self-Assembly. ACS APPLIED BIO MATERIALS 2019; 2:2812-2821. [DOI: 10.1021/acsabm.9b00229] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Kiheon Baek
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Alexander D. Noblett
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Laura J. Suggs
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
35
|
Remington JM, McCullagh M, Kohler B. Molecular Dynamics Simulations of 2-Aminopurine-Labeled Dinucleoside Monophosphates Reveal Multiscale Stacking Kinetics. J Phys Chem B 2019; 123:2291-2304. [PMID: 30767498 DOI: 10.1021/acs.jpcb.8b12172] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics (MD) simulations of 2-aminopurine (2Ap)-labeled DNA dinucleoside monophosphates (DNMPs) were performed to investigate the hypothesis that base stacking dynamics occur on timescales sufficiently rapid to influence the emission signals measured in time-resolved fluorescence experiments. Analysis of multiple microsecond-length trajectories shows that the DNMPs sample all four coplanar stacking motifs. In addition, three metastable unstacked conformations are detected. A hidden Markov-state model (HMSM) was applied to the simulations to estimate transition rates between the stacked and unstacked states. Transitions between different stacked states generally occur at higher rates when the number of nucleobase faces requiring desolvation is minimized. Time constants for structural relaxation range between 1.6 and 25 ns, suggesting that emission from photoexcited 2Ap, which has an excited-state lifetime of 10 ns, is sensitive to base stacking kinetics. A master equation model for the excited-state population of 2Ap predicts multiexponential emission decays that reproduce the sub-10 ns emission decay lifetimes and amplitudes seen in experiments. Combining MD simulations with HMSM analysis is a powerful way to understand the dynamics that influence 2Ap excited-state relaxation and represents an important step toward using observed emission signals to validate MD simulations.
Collapse
Affiliation(s)
- Jacob M Remington
- Department of Chemistry and Biochemistry , Montana State University , Bozeman , Montana 59717 , United States
| | - Martin McCullagh
- Department of Chemistry , Colorado State University , Fort Collins , Colorado 80523 , United States
| | - Bern Kohler
- Department of Chemistry and Biochemistry , The Ohio State University , 100 West 18th Avenue , Columbus , Ohio 43210 , United States
| |
Collapse
|
36
|
Erlenbach N, Grünewald C, Krstic B, Heckel A, Prisner TF. "End-to-end" stacking of small dsRNA. RNA (NEW YORK, N.Y.) 2019; 25:239-246. [PMID: 30404925 PMCID: PMC6348986 DOI: 10.1261/rna.068130.118] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Accepted: 11/06/2018] [Indexed: 05/08/2023]
Abstract
PELDOR (pulsed electron-electron double resonance) is an established method to study intramolecular distances and can give evidence for conformational changes and flexibilities. However, it can also be used to study intermolecular interactions as for example oligerimization. Here, we used PELDOR to study the "end-to-end" stacking of small double-stranded (ds) RNAs. For this study, the dsRNA molecules were only singly labeled with the spin label TPA to avoid multispin effects and to measure only the intermolecular stacking interactions. It can be shown that small dsRNAs tend to assemble to rod-like structures due to π-π interactions between the base pairs at the end of the strands. On the one hand, these interactions can influence or complicate measurements aimed at the determining of the structure and dynamics of the dsRNA molecule itself. On the other hand, it can be interesting to study such intermolecular stacking interactions in more detail, as for example their dependence on ion concentration. We quantitatively determined the stacking probability as a function of the monovalent NaCl salt and the dsRNA concentration. From these data, the dissociation constant Kd was deduced and found to depend on the ratio between the NaCl salt and dsRNA concentrations. Additionally, the distances and distance distributions obtained predict a model for the stacking geometry of dsRNAs. Introducing a nucleotide overhangs at one end of the dsRNA molecule restricts the stacking to the other end, leading only to dimer formations. Introducing such an overhang at both ends of the dsRNA molecule fully suppresses stacking, as we demonstrate by PELDOR experiments quantitatively.
Collapse
Affiliation(s)
- Nicole Erlenbach
- Institute of Physical and Theoretical Chemistry, Center of Biomolecular Magnetic Resonance, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Christian Grünewald
- Institute of Organic Chemistry and Chemical Biology, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Bisera Krstic
- Institute of Physical and Theoretical Chemistry, Center of Biomolecular Magnetic Resonance, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Alexander Heckel
- Institute of Organic Chemistry and Chemical Biology, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Thomas F Prisner
- Institute of Physical and Theoretical Chemistry, Center of Biomolecular Magnetic Resonance, Goethe University, D-60438 Frankfurt am Main, Germany
| |
Collapse
|
37
|
Schrodt MV, Andrews CT, Elcock AH. Correction to Large-Scale Analysis of 48 DNA and 48 RNA Tetranucleotides Studied by 1 μs Explicit-Solvent Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:6744-6745. [PMID: 30427670 DOI: 10.1021/acs.jctc.8b00595] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
38
|
Andrews CT, Campbell BA, Elcock AH. Correction to Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:6746-6747. [DOI: 10.1021/acs.jctc.8b00596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
39
|
Smith LG, Tan Z, Spasic A, Dutta D, Salas-Estrada LA, Grossfield A, Mathews DH. Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations. J Chem Theory Comput 2018; 14:6598-6612. [PMID: 30375860 DOI: 10.1021/acs.jctc.8b00633] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χOL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.
Collapse
Affiliation(s)
- Louis G Smith
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Zhen Tan
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Aleksandar Spasic
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Debapratim Dutta
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Leslie A Salas-Estrada
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States
| | - Alan Grossfield
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States
| | - David H Mathews
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Department of Biostatistics and Computational Biology , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| |
Collapse
|
40
|
Chen YL, Lee T, Elber R, Pollack L. Conformations of an RNA Helix-Junction-Helix Construct Revealed by SAXS Refinement of MD Simulations. Biophys J 2018; 116:19-30. [PMID: 30558889 DOI: 10.1016/j.bpj.2018.11.020] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2018] [Revised: 11/02/2018] [Accepted: 11/12/2018] [Indexed: 10/27/2022] Open
Abstract
RNA is involved in a broad range of biological processes that extend far beyond translation. Many of RNA's recently discovered functions rely on folding to a specific conformation or transitioning between conformations. The RNA structure contains rigid, short basepaired regions connected by more flexible linkers. Studies of model constructs such as small helix-junction-helix (HJH) motifs are useful in understanding how these elements work together to determine RNA conformation. Here, we reveal the full ensemble of solution structures assumed by a model RNA HJH. We apply small-angle x-ray scattering and an ensemble optimization method to selectively refine models generated by all-atom molecular dynamics simulations. The expectation of a broad distribution of helix orientations, at and above physiological ionic strength, is not met. Instead, this analysis shows that the HJH structures are dominated by two distinct conformations at moderate to high ionic strength. Atomic structures, selected from the molecular dynamics simulations, reveal strong base-base interactions in the junction that critically constrain the conformational space available to the HJH molecule and lead to a surprising re-extension at high salt. These results are corroborated by comparison with previous single-molecule fluorescence resonance energy transfer experiments on the same constructs.
Collapse
Affiliation(s)
- Yen-Lin Chen
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York
| | - Tongsik Lee
- Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas
| | - Ron Elber
- Department of Chemistry and Biochemistry, University of Texas at Austin, Austin, Texas; Institute of Computational Sciences and Engineering, University of Texas at Austin, Austin, Texas
| | - Lois Pollack
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York.
| |
Collapse
|
41
|
Brown RF, Andrews CT, Elcock AH. Correction to Stacking Free Energies of All DNA and RNA Nucleoside Pairs and Dinucleoside-Monophosphates Computed Using Recently Revised AMBER Parameters and Compared with Experiment. J Chem Theory Comput 2018; 14:6742-6743. [PMID: 30427669 DOI: 10.1021/acs.jctc.8b00594] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
42
|
Chen YL, Sutton JL, Pollack L. How the Conformations of an Internal Junction Contribute to Fold an RNA Domain. J Phys Chem B 2018; 122:11363-11372. [PMID: 30285445 DOI: 10.1021/acs.jpcb.8b07262] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Like proteins, some RNAs fold to compact structures. We can model functional RNAs as a series of short, rigid, base-paired elements, connected by non-base-paired nucleotides that serve as junctions. These connecting regions bend and twist, facilitating the formation of tertiary contacts that stabilize compact states. Here, we explore the roles of salt and junction sequence in determining the structures of a ubiquitous connector: an asymmetric internal loop. We focus on the J5/5a junction from the widely studied P4-P6 domain of the Tetrahymena ribozyme. Following the addition of magnesium ions to fold P4-P6, this junction bends dramatically, bringing the two halves of the RNA domain together for tertiary contact engagement. Using single-molecule fluorescence resonance energy transfer (smFRET), we examine the role of sequence and salt on model RNA constructs that contain these junction regions. We explore the wild-type J5/5a junction as well as two sequence variants. These junctions display distinct, salt-dependent conformations. Small-angle X-ray scattering (SAXS) measurements verify that these effects persist in the full-length P4-P6 domain. These measurements underscore the importance of junction sequence and interactions with ions in facilitating RNA folding.
Collapse
Affiliation(s)
- Yen-Lin Chen
- School of Applied and Engineering Physics , Cornell University , Ithaca , New York 14853 , United States
| | - Julie L Sutton
- School of Applied and Engineering Physics , Cornell University , Ithaca , New York 14853 , United States
| | - Lois Pollack
- School of Applied and Engineering Physics , Cornell University , Ithaca , New York 14853 , United States
| |
Collapse
|
43
|
Capobianco A, Velardo A, Peluso A. Single-Stranded DNA Oligonucleotides Retain Rise Coordinates Characteristic of Double Helices. J Phys Chem B 2018; 122:7978-7989. [PMID: 30070843 DOI: 10.1021/acs.jpcb.8b04542] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The structures of single-stranded DNA oligonucleotides from dimeric to hexameric sequences have been thoroughly investigated. Computations performed at the density functional level of theory including dispersion forces and solvation show that single-stranded helices adopt conformations very close to crystallographic B-DNA, with rise coordinates amounting up to 3.3 Å. Previous results, suggesting that single strands should be shorter than double helices, largely originated from the incompleteness of the adopted basis set. Although sensible deviations with respect to standard B-DNA are predicted, computations indicate that sequences rich in stacked adenines are the most ordered ones, favoring the B-DNA pattern and inducing regular arrangements also on flanking nucleobases. Several structural properties of double helices rich in adenine are indeed already reflected by the corresponding single strands.
Collapse
Affiliation(s)
- Amedeo Capobianco
- Dipartimento di Chimica e Biologia "A. Zambelli" , Università di Salerno , Via Giovanni Paolo II , I-84084 Fisciano (SA) , Italy
| | - Amalia Velardo
- Dipartimento di Chimica e Biologia "A. Zambelli" , Università di Salerno , Via Giovanni Paolo II , I-84084 Fisciano (SA) , Italy
| | - Andrea Peluso
- Dipartimento di Chimica e Biologia "A. Zambelli" , Università di Salerno , Via Giovanni Paolo II , I-84084 Fisciano (SA) , Italy
| |
Collapse
|
44
|
Chakraborty D, Hori N, Thirumalai D. Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA. J Chem Theory Comput 2018; 14:3763-3779. [PMID: 29870236 PMCID: PMC6423546 DOI: 10.1021/acs.jctc.8b00091] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.
Collapse
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Naoto Hori
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - D. Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
45
|
Svidritskiy E, Demo G, Korostelev AA. Mechanism of premature translation termination on a sense codon. J Biol Chem 2018; 293:12472-12479. [PMID: 29941456 DOI: 10.1074/jbc.aw118.003232] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Accurate translation termination by release factors (RFs) is critical for the integrity of cellular proteomes. Premature termination on sense codons, for example, results in truncated proteins, whose accumulation could be detrimental to the cell. Nevertheless, some sense codons are prone to triggering premature termination, but the structural basis for this is unclear. To investigate premature termination, we determined a cryo-EM structure of the Escherichia coli 70S ribosome bound with RF1 in response to a UAU (Tyr) sense codon. The structure reveals that RF1 recognizes a UAU codon similarly to a UAG stop codon, suggesting that sense codons induce premature termination because they structurally mimic a stop codon. Hydrophobic interaction between the nucleobase of U3 (the third position of the UAU codon) and conserved Ile-196 in RF1 is important for misreading the UAU codon. Analyses of RNA binding in ribonucleoprotein complexes or by amino acids reveal that Ile-U packing is a frequent protein-RNA-binding motif with key functional implications. We discuss parallels with eukaryotic translation termination by the release factor eRF1.
Collapse
Affiliation(s)
- Egor Svidritskiy
- From the RNA Therapeutics Institute, University of Massachusetts Medical School, Worcester, Massachusetts 01605
| | - Gabriel Demo
- From the RNA Therapeutics Institute, University of Massachusetts Medical School, Worcester, Massachusetts 01605
| | - Andrei A Korostelev
- From the RNA Therapeutics Institute, University of Massachusetts Medical School, Worcester, Massachusetts 01605
| |
Collapse
|
46
|
Boehm BJ, Whidborne C, Button AL, Pukala TL, Huang DM. DNA triplex structure, thermodynamics, and destabilisation: insight from molecular simulations. Phys Chem Chem Phys 2018; 20:14013-14023. [PMID: 29744501 DOI: 10.1039/c8cp02385a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Molecular dynamics simulations are used to elucidate the structure and thermodynamics of DNA triplexes associated with the neurodegenerative disease Friedreich's ataxia (FRDA), as well as complexes of these triplexes with the small molecule netropsin, which is known to destabilise triplexes. The ability of molecular simulations in explicit solvent to accurately capture triplex thermodynamics is verified for the first time, with the free energy to dissociate a 15-base antiparallel purine triplex-forming oligomer (TFO) from the duplex found to be slightly higher than reported experimentally. The presence of netropsin in the minor groove destabilises the triplex as expected, reducing the dissociation free energy by approximately 50%. Netropsin binding is associated with localised narrowing of the minor groove near netropsin, an effect that has previously been under contention. This leads to localised widening of the major groove, weakening hydrogen bonds between the TFO and duplex. Consequently, destabilisation is found to be highly localised, occurring only when netropsin is bound directly opposite the TFO. The simulations also suggest that near saturation of the minor groove with ligand is required for complete triplex dissociation. A structural analysis of the DNA triplexes that can form with the FRDA-related duplex sequence indicates that the triplex with a parallel homopyrimidine TFO is likely to be more stable than the antiparallel homopurine-TFO triplex, which may have implications for disease onset and treatment.
Collapse
Affiliation(s)
- Belinda J Boehm
- Department of Chemistry, School of Physical Sciences, The University of Adelaide, Adelaide, Australia.
| | | | | | | | | |
Collapse
|
47
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
48
|
Yoo J, Aksimentiev A. New tricks for old dogs: improving the accuracy of biomolecular force fields by pair-specific corrections to non-bonded interactions. Phys Chem Chem Phys 2018; 20:8432-8449. [PMID: 29547221 DOI: 10.1039/c7cp08185e] [Citation(s) in RCA: 162] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
In contrast to ordinary polymers, the vast majority of biological macromolecules adopt highly ordered three-dimensional structures that define their functions. The key to folding of a biopolymer into a unique 3D structure or to an assembly of several biopolymers into a functional unit is a delicate balance between the attractive and repulsive forces that also makes such self-assembly reversible under physiological conditions. The all-atom molecular dynamics (MD) method has emerged as a powerful tool for studies of individual biomolecules and their functional assemblies, encompassing systems of ever increasing complexity. However, advances in parallel computing technology have outpaced the development of the underlying theoretical models-the molecular force fields, pushing the MD method into an untested territory. Recent tests of the MD method have found the most commonly used molecular force fields to be out of balance, overestimating attractive interactions between charged and hydrophobic groups, which can promote artificial aggregation in MD simulations of multi-component protein, nucleic acid, and lipid systems. One route towards improving the force fields is through the NBFIX corrections method, in which the intermolecular forces are calibrated against experimentally measured quantities such as osmotic pressure by making atom pair-specific adjustments to the non-bonded interactions. In this article, we review development of the NBFIX (Non-Bonded FIX) corrections to the AMBER and CHARMM force fields and discuss their implications for MD simulations of electrolyte solutions, dense DNA systems, Holliday junctions, protein folding, and lipid bilayer membranes.
Collapse
Affiliation(s)
- Jejoong Yoo
- Center for the Physics of Living Cells, Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA. and Center for Self-assembly and Complexity, Institute for Basic Science, Pohang, 37363, Republic of Korea
| | - Aleksei Aksimentiev
- Center for the Physics of Living Cells, Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA.
| |
Collapse
|
49
|
Milovanović B, Kojić M, Petković M, Etinski M. New Insight into Uracil Stacking in Water from ab Initio Molecular Dynamics. J Chem Theory Comput 2018; 14:2621-2632. [DOI: 10.1021/acs.jctc.8b00139] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000 Belgrade, Serbia
| | - Marko Kojić
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000 Belgrade, Serbia
| |
Collapse
|
50
|
New developments in force fields for biomolecular simulations. Curr Opin Struct Biol 2018; 49:129-138. [DOI: 10.1016/j.sbi.2018.02.002] [Citation(s) in RCA: 136] [Impact Index Per Article: 22.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2017] [Revised: 01/30/2018] [Accepted: 02/04/2018] [Indexed: 11/18/2022]
|