1
|
Lazzeri G, Micheletti C, Pasquali S, Faccioli P. RNA folding pathways from all-atom simulations with a variationally improved history-dependent bias. Biophys J 2023; 122:3089-3098. [PMID: 37355771 PMCID: PMC10432211 DOI: 10.1016/j.bpj.2023.06.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 05/03/2023] [Accepted: 06/15/2023] [Indexed: 06/26/2023] Open
Abstract
Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. As a step forward in tackling these issues, we extend to RNA an enhanced path-sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force promoting the generation of productive folding trajectories in an all-atom force field with explicit solvent. A rigorous variational principle is then applied to minimize the effect of the bias. Here, we report on an application of this method to RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. By comparison with analog simulations performed on small proteins with similar size and architecture, we show that the RNA folding landscape is significantly more frustrated, even for relatively small chains with a simple topology. The predicted RNA folding mechanisms are found to be consistent with the available experiments and some of the existing coarse-grained models. Due to its computational performance, this scheme provides a promising platform to efficiently gather atomistic RNA folding trajectories, thus retain the information about the chemical composition of the sequence.
Collapse
Affiliation(s)
- Gianmarco Lazzeri
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany; Physics Department of Trento University, Povo (Trento), Italy
| | | | - Samuela Pasquali
- Laboratoire Cibles Thérapeutiques et Conception de Médicaments, Université Paris Cité, Paris, France; Laboratoire Biologie Fonctionnelle et Adaptative, Université Paris Cité, Paris, France.
| | - Pietro Faccioli
- Physics Department of Trento University, Povo (Trento), Italy; INFN-TIFPA, Povo (Trento), Italy.
| |
Collapse
|
2
|
Wang X, Tan YL, Yu S, Shi YZ, Tan ZJ. Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions. Biophys J 2023; 122:1503-1516. [PMID: 36924021 PMCID: PMC10147842 DOI: 10.1016/j.bpj.2023.03.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 03/05/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023] Open
Abstract
RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.
Collapse
Affiliation(s)
- Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
3
|
Interdomain dynamics in human Replication Protein A regulates kinetics and thermodynamics of its binding to ssDNA. PLoS One 2023; 18:e0278396. [PMID: 36656834 PMCID: PMC9851514 DOI: 10.1371/journal.pone.0278396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Accepted: 11/15/2022] [Indexed: 01/20/2023] Open
Abstract
Human Replication Protein A (hRPA) is a multidomain protein that interacts with ssDNA intermediates to provide the latter much-needed stability during DNA metabolism and maintain genomic integrity. Although the ssDNA organization with hRPA was studied recently through experimental means, characterizing the underlying mechanism at the atomic level remains challenging because of the dynamic domain architecture of hRPA and poorly understood heterogeneity of ssDNA-protein interactions. Here, we used a computational framework, precisely tailored to capture protein-ssDNA interactions, and investigated the binding of hRPA with a 60 nt ssDNA. Two distinct binding mechanisms are realized based on the hRPA domain flexibility. For a rigid domain architecture of hRPA, ssDNA binds sequentially with hRPA domains, resulting in slow association kinetics. The binding pathway involves the formation of stable and distinct intermediate states. On contrary, for a flexible domain architecture of hRPA, ssDNA binds synergistically to the A and B domains followed by the rest of hRPA. The domain dynamics in hRPA alleviates the free energy cost of domain orientation necessary for specific binding with ssDNA, leading to fast association kinetics along a downhill binding free energy landscape. An ensemble of free energetically degenerate intermediate states is encountered that makes it arduous to characterize them structurally. An excellent match between our results with the available experimental observations provides new insights into the rich dynamics of hRPA binding to ssDNA and in general paves the way to investigate intricate details of ssDNA-protein interactions, crucial for cellular functioning.
Collapse
|
4
|
Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods. QRB DISCOVERY 2022. [PMID: 37529277 PMCID: PMC10392686 DOI: 10.1017/qrd.2022.19] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Abstract
While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.
Collapse
|
5
|
Dey P, Bhattacherjee A. Structural Basis of Enhanced Facilitated Diffusion of DNA-Binding Protein in Crowded Cellular Milieu. Biophys J 2020; 118:505-517. [PMID: 31862109 PMCID: PMC6976804 DOI: 10.1016/j.bpj.2019.11.3388] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2019] [Revised: 11/03/2019] [Accepted: 11/18/2019] [Indexed: 02/06/2023] Open
Abstract
Although the fast association between DNA-binding proteins (DBPs) and DNA is explained by a facilitated diffusion mechanism, in which DBPs adopt a weighted combination of three-dimensional diffusion and one-dimensional (1D) sliding and hopping modes of transportation, the role of cellular environment that contains many nonspecifically interacting proteins and other biomolecules is mostly overlooked. By performing large-scale computational simulations with an appropriately tuned model of protein and DNA in the presence of nonspecifically interacting bulk and DNA-bound crowders (genomic crowders), we demonstrate the structural basis of the enhanced facilitated diffusion of DBPs inside a crowded cellular milieu through, to our knowledge, novel 1D scanning mechanisms. In this one-dimensional scanning mode, the protein can float along the DNA under the influence of nonspecific interactions of bulk crowder molecules. The search mode is distinctly different compared to usual 1D sliding and hopping dynamics in which protein diffusion is regulated by the DNA electrostatics. In contrast, the presence of genomic crowders expedites the target search process by transporting the protein over DNA segments through the formation of a transient protein-crowder bridged complex. By analyzing the ruggedness of the associated potential energy landscape, we underpin the molecular origin of the kinetic advantages of these search modes and show that they successfully explain the experimentally observed acceleration of facilitated diffusion of DBPs by molecular crowding agents and crowder-concentration-dependent enzymatic activity of transcription factors. Our findings provide crucial insights into gene regulation kinetics inside the crowded cellular milieu.
Collapse
Affiliation(s)
- Pinki Dey
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Arnab Bhattacherjee
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India.
| |
Collapse
|
6
|
Dey P, Bhattacherjee A. Mechanism of Facilitated Diffusion of DNA Repair Proteins in Crowded Environment: Case Study with Human Uracil DNA Glycosylase. J Phys Chem B 2019; 123:10354-10364. [DOI: 10.1021/acs.jpcb.9b07342] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Pinki Dey
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India 110067
| | - Arnab Bhattacherjee
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India 110067
| |
Collapse
|
7
|
Hori N, Denesyuk NA, Thirumalai D. Frictional Effects on RNA Folding: Speed Limit and Kramers Turnover. J Phys Chem B 2018; 122:11279-11288. [PMID: 30179471 DOI: 10.1021/acs.jpcb.8b07129] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
We investigated frictional effects on the folding rates of a human telomerase hairpin (hTR HP) and H-type pseudoknot from the Beet Western Yellow Virus (BWYV PK) using simulations of the Three Interaction Site (TIS) model for RNA. The heat capacity from TIS model simulations, calculated using temperature replica exchange simulations, reproduces nearly quantitatively the available experimental data for the hTR HP. The corresponding results for BWYV PK serve as predictions. We calculated the folding rates ( kF) from more than 100 folding trajectories for each value of the solvent viscosity (η) at a fixed salt concentration of 200 mM. By using the theoretical estimate (∝ √N where N is the number of nucleotides) for folding free energy barrier, kF data for both the RNAs are quantitatively fit using one-dimensional Kramers's theory with two parameters specifying the curvatures in the unfolded basin and the barrier top. In the high-friction regime (η ≳ 10-5 Pa·s), for both HP and PK, kF values decrease as 1/η, whereas in the low friction regime, kF values increase as η increases, leading to a maximum folding rate at a moderate viscosity (∼10-6 Pa·s), which is the Kramers turnover. From the fits, we find that the speed limit to RNA folding at water viscosity is between 1 and 4 μs, which is in accord with our previous theoretical prediction as well as results from several single molecule experiments. Both the RNA constructs fold by parallel pathways. Surprisingly, we find that the flux through the pathways could be altered by changing solvent viscosity, a prediction that is more easily testable in RNA than in proteins.
Collapse
Affiliation(s)
- Naoto Hori
- Department of Chemistry , University of Texas , Austin , Texas 78712 , United States
| | - Natalia A Denesyuk
- Biophysics program, Institute for Physical Science and Technology , University of Maryland , College Park , Maryland 20742-2431 , United States
| | - D Thirumalai
- Department of Chemistry , University of Texas , Austin , Texas 78712 , United States
| |
Collapse
|
8
|
Wang Y, Wang Z, Liu T, Gong S, Zhang W. Effects of flanking regions on HDV cotranscriptional folding kinetics. RNA (NEW YORK, N.Y.) 2018; 24:1229-1240. [PMID: 29954950 PMCID: PMC6097654 DOI: 10.1261/rna.065961.118] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/06/2018] [Accepted: 06/25/2018] [Indexed: 05/20/2023]
Abstract
Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.
Collapse
Affiliation(s)
- Yanli Wang
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, P.R. China
| | - Zhen Wang
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, P.R. China
| | - Taigang Liu
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, P.R. China
| | - Sha Gong
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, P.R. China
| | - Wenbing Zhang
- Department of Physics, Wuhan University, Wuhan, Hubei 430072, P.R. China
| |
Collapse
|
9
|
Monovalent ions modulate the flux through multiple folding pathways of an RNA pseudoknot. Proc Natl Acad Sci U S A 2018; 115:E7313-E7322. [PMID: 30012621 PMCID: PMC6077692 DOI: 10.1073/pnas.1717582115] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The assembly mechanism of RNA, vital to describing its functions, depends on both the sequence and the metal ion concentration. How the latter influences the folding trajectories remains an important unsolved problem. Here, we examine the folding pathways of an RNA pseudoknot (PK) with key functional roles in transcription and translation, using a combination of experiments and simulations. We demonstrate that the PK, consisting of two hairpins with differing stabilities, folds by parallel pathways. Surprisingly, the flux between them is modulated by monovalent salt concentration. Our work shows that the order of assembly of PKs is determined by the relative stability of the hairpins, implying that the folding landscape can be controlled by sequence and ion concentration. The functions of RNA pseudoknots (PKs), which are minimal tertiary structural motifs and an integral part of several ribozymes and ribonucleoprotein complexes, are determined by their structure, stability, and dynamics. Therefore, it is important to elucidate the general principles governing their thermodynamics/folding mechanisms. Here, we combine laser temperature-jump experiments and coarse-grained simulations to determine the folding/unfolding pathways of VPK, a variant of the mouse mammary tumor virus (MMTV) PK involved in ribosomal frameshifting. Fluorescent nucleotide analogs (2-aminopurine and pyrrolocytidine) placed at different stem/loop positions in the PK serve as local probes allowing us to monitor the order of assembly of VPK that has two constituent hairpins with different intrinsic stabilities. We show that at 50 mM KCl, the dominant folding pathway populates only the more stable hairpin intermediate; as the salt concentration is increased, a parallel folding pathway emerges involving the less stable hairpin as an alternate intermediate. Notably, the flux between the pathways is modulated by the ionic strength. Our findings support the principle that the order of PK structure formation is determined by the relative stabilities of the hairpins, which can be altered by sequence variations or salt concentrations. The experimental results of salt effects on the partitioning between the two folding pathways are in remarkable agreement with simulations that were performed with no adjustable parameters. Our study not only unambiguously demonstrates that VPK folds by parallel pathways but also showcases the power of combining experiments and simulations for a more enriched description of RNA self-assembly.
Collapse
|
10
|
Villada-Balbuena M, Carbajal-Tinoco MD. One-bead coarse-grained model for RNA dynamics. J Chem Phys 2018; 146:045101. [PMID: 28147510 DOI: 10.1063/1.4974899] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
We present a revised version of a coarse-grained model for RNA dynamics. In such approach, the description of nucleotides is reduced to single points that interact between them through a series of effective pair potentials that were obtained from an improved analysis of RNA structures from the Protein Data Bank. These interaction potentials are the main constituents of a Brownian dynamics simulation algorithm that allows to perform a variety of tasks by taking advantage of the reduced number of variables. Such tasks include the prediction of the three-dimensional configuration of a series of test molecules. Moreover, the model permits the inclusion of effective magnesium ions and the ends of the RNA chains can be pulled with an external force to study the process of unfolding. In spite of the simplicity of the model, we obtain a good agreement with the experimental results.
Collapse
Affiliation(s)
- Mario Villada-Balbuena
- Departamento de Física, Centro de Investigación y de Estudios Avanzados del IPN, Av. Instituto Politécnico Nacional No. 2508, Colonia San Pedro Zacatenco, CP 07360 Ciudad de México, Mexico
| | - Mauricio D Carbajal-Tinoco
- Departamento de Física, Centro de Investigación y de Estudios Avanzados del IPN, Av. Instituto Politécnico Nacional No. 2508, Colonia San Pedro Zacatenco, CP 07360 Ciudad de México, Mexico
| |
Collapse
|
11
|
Shi YZ, Jin L, Feng CJ, Tan YL, Tan ZJ. Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions. PLoS Comput Biol 2018; 14:e1006222. [PMID: 29879103 PMCID: PMC6007934 DOI: 10.1371/journal.pcbi.1006222] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Revised: 06/19/2018] [Accepted: 05/22/2018] [Indexed: 01/30/2023] Open
Abstract
RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.
Collapse
Affiliation(s)
- Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematics and Computer Science, Wuhan Textile University, Wuhan, China
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Lei Jin
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Chen-Jie Feng
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Lan Tan
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| |
Collapse
|
12
|
Role of Macromolecular Crowding on the Intracellular Diffusion of DNA Binding Proteins. Sci Rep 2018; 8:844. [PMID: 29339733 PMCID: PMC5770392 DOI: 10.1038/s41598-017-18933-3] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2017] [Accepted: 12/20/2017] [Indexed: 11/08/2022] Open
Abstract
Recent experiments suggest that cellular crowding facilitates the target search dynamics of proteins on DNA, the mechanism of which is not yet known. By using large scale computer simulations, we show that two competing factors, namely the width of the depletion layer that separates the crowder cloud from the DNA molecule and the degree of protein-crowder crosstalk, act in harmony to affect the target search dynamics of proteins. The impacts vary from nonspecific to specific target search regime. During a nonspecific search, dynamics of a protein is only minimally affected, whereas, a significantly different behaviour is observed when the protein starts forming a specific protein-DNA complex. We also find that the severity of impacts largely depends upon physiological crowder concentration and deviation from it leads to attenuation in the binding kinetics. Based on extensive kinetic study and binding energy landscape analysis, we further present a comprehensive molecular description of the search process that allows us to interpret the experimental findings.
Collapse
|
13
|
Markegard CB, Gallivan CP, Cheng DD, Nguyen HD. Effects of Concentration and Temperature on DNA Hybridization by Two Closely Related Sequences via Large-Scale Coarse-Grained Simulations. J Phys Chem B 2016; 120:7795-806. [PMID: 27447850 DOI: 10.1021/acs.jpcb.6b03937] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
A newly developed coarse-grained model called BioModi is utilized to elucidate the effects of temperature and concentration on DNA hybridization in self-assembly. Large-scale simulations demonstrate that complementary strands of either the tetrablock sequence or randomized sequence with equivalent number of cytosine or guanine nucleotides can form completely hybridized double helices. Even though the end states are the same for the two sequences, there exist multiple kinetic pathways that are populated with a wider range of transient aggregates of different sizes in the system of random sequences compared to that of the tetrablock sequence. The ability of these aggregates to undergo the strand displacement mechanism to form only double helices depends upon the temperature and DNA concentration. On one hand, low temperatures and high concentrations drive the formation and enhance stability of large aggregating species. On the other hand, high temperatures destabilize base-pair interactions and large aggregates. There exists an optimal range of moderate temperatures and low concentrations that allow minimization of large aggregate formation and maximization of fully hybridized dimers. Such investigation on structural dynamics of aggregating species by two closely related sequences during the self-assembly process demonstrates the importance of sequence design in avoiding the formation of metastable species. Finally, from kinetic modeling of self-assembly dynamics, the activation energy for the formation of double helices was found to be in agreement with experimental results. The framework developed in this work can be applied to the future design of DNA nanostructures in both fields of structural DNA nanotechnology and dynamic DNA nanotechnology wherein equilibrium end states and nonequilibrium dynamics are equally important requiring investigation in cooperation.
Collapse
Affiliation(s)
- Cade B Markegard
- Department of Chemical Engineering and Materials Science, University of California, Irvine , Irvine, California 92697, United States
| | - Cameron P Gallivan
- Department of Chemical Engineering and Materials Science, University of California, Irvine , Irvine, California 92697, United States
| | - Darrell D Cheng
- Department of Chemical Engineering and Materials Science, University of California, Irvine , Irvine, California 92697, United States
| | - Hung D Nguyen
- Department of Chemical Engineering and Materials Science, University of California, Irvine , Irvine, California 92697, United States
| |
Collapse
|
14
|
Hori N, Denesyuk NA, Thirumalai D. Salt Effects on the Thermodynamics of a Frameshifting RNA Pseudoknot under Tension. J Mol Biol 2016; 428:2847-59. [PMID: 27315694 DOI: 10.1016/j.jmb.2016.06.002] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Revised: 06/06/2016] [Accepted: 06/07/2016] [Indexed: 12/24/2022]
Abstract
Because of the potential link between -1 programmed ribosomal frameshifting and response of a pseudoknot (PK) RNA to force, a number of single-molecule pulling experiments have been performed on PKs to decipher the mechanism of programmed ribosomal frameshifting. Motivated in part by these experiments, we performed simulations using a coarse-grained model of RNA to describe the response of a PK over a range of mechanical forces (fs) and monovalent salt concentrations (Cs). The coarse-grained simulations quantitatively reproduce the multistep thermal melting observed in experiments, thus validating our model. The free energy changes obtained in simulations are in excellent agreement with experiments. By varying f and C, we calculated the phase diagram that shows a sequence of structural transitions, populating distinct intermediate states. As f and C are changed, the stem-loop tertiary interactions rupture first, followed by unfolding of the 3'-end hairpin (I⇌F). Finally, the 5'-end hairpin unravels, producing an extended state (E⇌I). A theoretical analysis of the phase boundaries shows that the critical force for rupture scales as (logCm)(α) with α=1(0.5) for E⇌I (I⇌F) transition. This relation is used to obtain the preferential ion-RNA interaction coefficient, which can be quantitatively measured in single-molecule experiments, as done previously for DNA hairpins. A by-product of our work is the suggestion that the frameshift efficiency is likely determined by the stability of the 5'-end hairpin that the ribosome first encounters during translation.
Collapse
Affiliation(s)
- Naoto Hori
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| | - Natalia A Denesyuk
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| | - D Thirumalai
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
15
|
Coarse-grained models for studying protein diffusion along DNA. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1262] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
16
|
Zhang X, Xu X, Yang Z, Burcke AJ, Gates KS, Chen SJ, Gu LQ. Mimicking Ribosomal Unfolding of RNA Pseudoknot in a Protein Channel. J Am Chem Soc 2015; 137:15742-52. [PMID: 26595106 PMCID: PMC4886178 DOI: 10.1021/jacs.5b07910] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Pseudoknots are a fundamental RNA tertiary structure with important roles in regulation of mRNA translation. Molecular force spectroscopic approaches such as optical tweezers can track the pseudoknot's unfolding intermediate states by pulling the RNA chain from both ends, but the kinetic unfolding pathway induced by this method may be different from that in vivo, which occurs during translation and proceeds from the 5' to 3' end. Here we developed a ribosome-mimicking, nanopore pulling assay for dissecting the vectorial unfolding mechanism of pseudoknots. The pseudoknot unfolding pathway in the nanopore, either from the 5' to 3' end or in the reverse direction, can be controlled by a DNA leader that is attached to the pseudoknot at the 5' or 3' ends. The different nanopore conductance between DNA and RNA translocation serves as a marker for the position and structure of the unfolding RNA in the pore. With this design, we provided evidence that the pseudoknot unfolding is a two-step, multistate, metal ion-regulated process depending on the pulling direction. Most notably, unfolding in both directions is rate-limited by the unzipping of the first helix domain (first step), which is Helix-1 in the 5' → 3' direction and Helix-2 in the 3' → 5' direction, suggesting that the initial unfolding step in either pulling direction needs to overcome an energy barrier contributed by the noncanonical triplex base-pairs and coaxial stacking interactions for the tertiary structure stabilization. These findings provide new insights into RNA vectorial unfolding mechanisms, which play an important role in biological functions including frameshifting.
Collapse
Affiliation(s)
- Xinyue Zhang
- Department of Bioengineering and Dalton Cardiovascular Research Center, University of Missouri, Columbia, Missouri 65211, United States
| | - Xiaojun Xu
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Zhiyu Yang
- Department of Chemistry and Department of Biochemistry, University of Missouri, Columbia, Missouri 65211, United States
| | - Andrew J. Burcke
- Department of Bioengineering and Dalton Cardiovascular Research Center, University of Missouri, Columbia, Missouri 65211, United States
| | - Kent S. Gates
- Department of Chemistry and Department of Biochemistry, University of Missouri, Columbia, Missouri 65211, United States
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Li-Qun Gu
- Department of Bioengineering and Dalton Cardiovascular Research Center, University of Missouri, Columbia, Missouri 65211, United States
| |
Collapse
|
17
|
Mondal A, Bhattacherjee A. Searching target sites on DNA by proteins: Role of DNA dynamics under confinement. Nucleic Acids Res 2015; 43:9176-86. [PMID: 26400158 PMCID: PMC4627088 DOI: 10.1093/nar/gkv931] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2015] [Revised: 08/15/2015] [Accepted: 09/07/2015] [Indexed: 02/07/2023] Open
Abstract
DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.
Collapse
Affiliation(s)
- Anupam Mondal
- Center for Computational Biology, Indraprastha Institute of Information Technology (IIIT) Delhi, New Delhi-110020, India
| | - Arnab Bhattacherjee
- Center for Computational Biology, Indraprastha Institute of Information Technology (IIIT) Delhi, New Delhi-110020, India
| |
Collapse
|
18
|
MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery. Int J Mol Sci 2015; 16:15872-902. [PMID: 26184179 PMCID: PMC4519929 DOI: 10.3390/ijms160715872] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2015] [Revised: 07/07/2015] [Accepted: 07/08/2015] [Indexed: 12/21/2022] Open
Abstract
While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.
Collapse
|
19
|
Cragnolini T, Laurin Y, Derreumaux P, Pasquali S. Coarse-Grained HiRE-RNA Model for ab Initio RNA Folding beyond Simple Molecules, Including Noncanonical and Multiple Base Pairings. J Chem Theory Comput 2015; 11:3510-22. [PMID: 26575783 DOI: 10.1021/acs.jctc.5b00200] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
HiRE-RNA is a coarse-grained model for RNA structure prediction and the dynamical study of RNA folding. Using a reduced set of particles and detailed interactions accounting for base-pairing and stacking, we show that noncanonical and multiple base interactions are necessary to capture the full physical behavior of complex RNAs. In this paper, we give a full account of the model and present results on the folding, stability, and free energy surfaces of 16 systems with 12 to 76 nucleotides of increasingly complex architectures, ranging from monomers to dimers, using a total of 850 μs of simulation time.
Collapse
Affiliation(s)
- Tristan Cragnolini
- Laboratoire de Biochimie Théorique UPR 9080 CNRS, Université Paris Diderot , Sorbonne, Paris Cité, IBPC 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Yoann Laurin
- Laboratoire de Biochimie Théorique UPR 9080 CNRS, Université Paris Diderot , Sorbonne, Paris Cité, IBPC 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique UPR 9080 CNRS, Université Paris Diderot , Sorbonne, Paris Cité, IBPC 13 rue Pierre et Marie Curie, 75005 Paris, France.,Institut Universitaire de France , Boulevard Saint-Michel, 75005 Paris, France
| | - Samuela Pasquali
- Laboratoire de Biochimie Théorique UPR 9080 CNRS, Université Paris Diderot , Sorbonne, Paris Cité, IBPC 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
20
|
Markegard CB, Fu IW, Reddy KA, Nguyen HD. Coarse-grained simulation study of sequence effects on DNA hybridization in a concentrated environment. J Phys Chem B 2015; 119:1823-34. [PMID: 25581253 DOI: 10.1021/jp509857k] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
A novel coarse-grained model is developed to elucidate thermodynamics and kinetic mechanisms of DNA self-assembly. It accounts for sequence and solvent conditions to capture key experimental results such as sequence-dependent thermal property and salt-dependent persistence length of ssDNA and dsDNA. Moreover, constant-temperature simulations on two single strands of a homogeneous sequence show two main mechanisms of hybridization: a slow slithering mechanism and a one-order faster zippering mechanism. Furthermore, large-scale simulations at a high DNA strand concentration demonstrate that DNA self-assembly is a robust and enthalpically driven process in which the formation of double helices is deciphered to occur via multiple self-assembly pathways including the strand displacement mechanism. However, sequence plays an important role in shifting the majority of one pathway over the others and controlling size distribution of self-assembled aggregates. This study yields a complex picture on the role of sequence on programmable self-assembly and demonstrates a promising simulation tool that is suitable for studies in DNA nanotechnology.
Collapse
Affiliation(s)
- Cade B Markegard
- Department of Chemical Engineering and Materials Science, University of California-Irvine , Irvine, California 92697-2575, United States
| | | | | | | |
Collapse
|
21
|
Chakraborty D, Collepardo-Guevara R, Wales DJ. Energy Landscapes, Folding Mechanisms, and Kinetics of RNA Tetraloop Hairpins. J Am Chem Soc 2014; 136:18052-61. [DOI: 10.1021/ja5100756] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Debayan Chakraborty
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | | | - David J. Wales
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
22
|
Bhattacherjee A, Levy Y. Search by proteins for their DNA target site: 1. The effect of DNA conformation on protein sliding. Nucleic Acids Res 2014; 42:12404-14. [PMID: 25324308 PMCID: PMC4227778 DOI: 10.1093/nar/gku932] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The recognition of DNA-binding proteins (DBPs) to their specific site often precedes by a search technique in which proteins slide, hop along the DNA contour or perform inter-segment transfer and 3D diffusion to dissociate and re-associate to distant DNA sites. In this study, we demonstrated that the strength and nature of the non-specific electrostatic interactions, which govern the search dynamics of DBPs, are strongly correlated with the conformation of the DNA. We tuned two structural parameters, namely curvature and the extent of helical twisting in circular DNA. These two factors are mutually independent of each other and can modulate the electrostatic potential through changing the geometry of the circular DNA conformation. The search dynamics for DBPs on circular DNA is therefore markedly different compared with linear B-DNA. Our results suggest that, for a given DBP, the rotation-coupled sliding dynamics is precluded in highly curved DNA (as well as for over-twisted DNA) because of the large electrostatic energy barrier between the inside and outside of the DNA molecule. Under such circumstances, proteins prefer to hop in order to explore interior DNA sites. The change in the balance between sliding and hopping propensities as a function of DNA curvature or twisting may result in different search efficiency and speed.
Collapse
Affiliation(s)
- Arnab Bhattacherjee
- Department of Structural Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Yaakov Levy
- Department of Structural Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| |
Collapse
|
23
|
Shi YZ, Wang FH, Wu YY, Tan ZJ. A coarse-grained model with implicit salt for RNAs: Predicting 3D structure, stability and salt effect. J Chem Phys 2014; 141:105102. [DOI: 10.1063/1.4894752] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Affiliation(s)
- Ya-Zhou Shi
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Feng-Hua Wang
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Yuan-Yan Wu
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| |
Collapse
|
24
|
Single-molecule measurements of the CCR5 mRNA unfolding pathways. Biophys J 2014; 106:244-52. [PMID: 24411256 DOI: 10.1016/j.bpj.2013.09.036] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2013] [Revised: 08/13/2013] [Accepted: 09/04/2013] [Indexed: 01/07/2023] Open
Abstract
Secondary or tertiary structure in an mRNA, such as a pseudoknot, can create a physical barrier that requires the ribosome to generate additional force to translocate. The presence of such a barrier can dramatically increase the probability that the ribosome will shift into an alternate reading frame, in which a different set of codons is recognized. The detailed biophysical mechanism by which frameshifting is induced remains unknown. Here we employ optical trapping techniques to investigate the structure of a -1 programmed ribosomal frameshift (-1 PRF) sequence element located in the CCR5 mRNA, which encodes a coreceptor for HIV-1 and is, to our knowledge, the first known human -1 PRF signal of nonviral origin. We begin by presenting a set of computationally predicted structures that include pseudoknots. We then employ what we believe to be new analytical techniques for measuring the effective free energy landscapes of biomolecules. We find that the -1 PRF element manifests several distinct unfolding pathways when subject to end-to-end force, one of which is consistent with a proposed pseudoknot conformation, and another of which we have identified as a folding intermediate. The dynamic ensemble of conformations that CCR5 mRNA exhibits in the single-molecule experiments may be a significant feature of the frameshifting mechanism.
Collapse
|
25
|
Hyeon C, Denesyuk NA, Thirumalai D. Development and Applications of Coarse-Grained Models for RNA. Isr J Chem 2014. [DOI: 10.1002/ijch.201400029] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
26
|
Leuchter JD, Green AT, Gilyard J, Rambarat CG, Cho SS. Coarse-Grained and Atomistic MD Simulations of RNA and DNA Folding. Isr J Chem 2014. [DOI: 10.1002/ijch.201400022] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
27
|
Atomistic picture for the folding pathway of a hybrid-1 type human telomeric DNA G-quadruplex. PLoS Comput Biol 2014; 10:e1003562. [PMID: 24722458 PMCID: PMC3983051 DOI: 10.1371/journal.pcbi.1003562] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Accepted: 02/26/2014] [Indexed: 12/20/2022] Open
Abstract
In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and K(+) ions explicitly modeled. Enabled by the powerful bias-exchange metadynamics and large-scale conventional molecular dynamic simulations, the free energy landscape of this G-DNA was obtained for the first time and four folding intermediates were identified, including a triplex and a basically formed quadruplex. The simulations also provided atomistic pictures for the structures and cation binding patterns of the intermediates. The results showed that the structure formation and cation binding are cooperative and mutually supporting each other. The syn/anti reorientation dynamics of the intermediates was also investigated. It was found that the nucleotides usually take correct syn/anti configurations when they form native and stable hydrogen bonds with the others, while fluctuating between two configurations when they do not. Misfolded intermediates with wrong syn/anti configurations were observed in the early intermediates but not in the later ones. Based on the simulations, we also discussed the roles of the non-native interactions. Besides, the formation process of the parallel conformation in the first two G-repeats and the associated reversal loop were studied. Based on the above results, we proposed a folding pathway for the hybrid-1 type G-quadruplex with atomistic details, which is new and more complete compared with previous ones. The knowledge gained for this type of G-DNA may provide a general insight for the folding of the other G-quadruplexes.
Collapse
|
28
|
Chen J, Gong S, Wang Y, Zhang W. Kinetic partitioning mechanism of HDV ribozyme folding. J Chem Phys 2014; 140:025102. [DOI: 10.1063/1.4861037] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
29
|
Harrison JG, Zheng YB, Beal PA, Tantillo DJ. Computational approaches to predicting the impact of novel bases on RNA structure and stability. ACS Chem Biol 2013; 8:2354-9. [PMID: 24063428 DOI: 10.1021/cb4006062] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The use of computational modeling techniques to gain insight into nucleobase interactions has been a challenging endeavor to date. Accurate treatment requires the tackling of many challenges but also holds the promise of great rewards. The development of effective computational approaches to predict the binding affinities of nucleobases and analogues can, for example, streamline the process of developing novel nucleobase modifications, which should facilitate the development of new RNAi-based therapeutics. This brief review focuses on available computational approaches to predicting base pairing affinity in RNA-based contexts such as nucleobase-nucleobase interactions in duplexes and nucleobase-protein interactions. The challenges associated with such modeling along with potential future directions for the field are highlighted.
Collapse
Affiliation(s)
- Jason G. Harrison
- Department of Chemistry, University of California−Davis, Davis, California 95616, United States
| | - Yvonne B. Zheng
- Department of Chemistry, University of California−Davis, Davis, California 95616, United States
| | - Peter A. Beal
- Department of Chemistry, University of California−Davis, Davis, California 95616, United States
| | - Dean J. Tantillo
- Department of Chemistry, University of California−Davis, Davis, California 95616, United States
| |
Collapse
|
30
|
Li R, Chen R, Chen P, Wen Y, Ke PC, Cho SS. Computational and experimental characterizations of silver nanoparticle-apolipoprotein biocorona. J Phys Chem B 2013; 117:13451-6. [PMID: 24073791 DOI: 10.1021/jp4061158] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
With the advancement of nanotoxicology and nanomedicine, it has been realized that nanoparticles (NPs) interact readily with biomolecular species and other chemical and organic matter to result in biocorona formation. The field of the environmental health and safety of nanotechnology, or NanoEHS, is currently lacking significant molecular-resolution data, and we set out to characterize biocorona formation through electron microscopy imaging and circular dichroism spectroscopy that inspired a novel approach for molecular dynamics (MD) simulations of protein-NP interactions. In our present study, we developed a novel GPU-optimized coarse-grained MD simulation methodology for the study of biocorona formation, a first in the field. Specifically, we performed MD simulations of a spherical, negatively charged citrate-covered silver nanoparticle (AgNP) interacting with 15 apolipoproteins. At low ion concentrations, we observed the formation of an AgNP-apolipoprotein biocorona. Consistent with the circular dichroism (CD) spectra, we observed a decrease in α-helices coupled with an increase in β-sheets in apolipoprotein upon biocorona formation.
Collapse
Affiliation(s)
- Rongzhong Li
- Department of Physics, Wake Forest University , Winston-Salem, North Carolina 27106, United States
| | | | | | | | | | | |
Collapse
|
31
|
Li R, Ge HW, Cho SS. Sequence-dependent base-stacking stabilities guide tRNA folding energy landscapes. J Phys Chem B 2013; 117:12943-52. [PMID: 23841777 DOI: 10.1021/jp402114p] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The folding of bacterial tRNAs with disparate sequences has been observed to proceed in distinct folding mechanisms despite their structural similarity. To explore the folding landscapes of tRNA, we performed ion concentration-dependent coarse-grained TIS model MD simulations of several E. coli tRNAs to compare their thermodynamic melting profiles to the classical absorbance spectra of Crothers and co-workers. To independently validate our findings, we also performed atomistic empirical force field MD simulations of tRNAs, and we compared the base-to-base distances from coarse-grained and atomistic MD simulations to empirical base-stacking free energies. We then projected the free energies to the secondary structural elements of tRNA, and we observe distinct, parallel folding mechanisms whose differences can be inferred on the basis of their sequence-dependent base-stacking stabilities. In some cases, a premature, nonproductive folding intermediate corresponding to the Ψ hairpin loop must backtrack to the unfolded state before proceeding to the folded state. This observation suggests a possible explanation for the fast and slow phases observed in tRNA folding kinetics.
Collapse
Affiliation(s)
- Rongzhong Li
- Department of Physics, Wake Forest University , Winston-Salem, North Carolina 27106, United States
| | | | | |
Collapse
|
32
|
Narayanan R, Zhu L, Velmurugu Y, Roca J, Kuznetsov SV, Prehna G, Lapidus LJ, Ansari A. Exploring the Energy Landscape of Nucleic Acid Hairpins Using Laser Temperature-Jump and Microfluidic Mixing. J Am Chem Soc 2012; 134:18952-63. [DOI: 10.1021/ja301218e] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
| | - Li Zhu
- Department of Physics
and Astronomy, Michigan State University, East Lansing, Michigan 48824,
United States
- Advanced
Photonics Center, Southeast University,
Nanjing 210096, China
| | | | | | | | | | - Lisa J. Lapidus
- Department of Physics
and Astronomy, Michigan State University, East Lansing, Michigan 48824,
United States
| | | |
Collapse
|
33
|
Abstract
Transitions between the different conformational states play a critical role in many RNA catalytic and regulatory functions. In this study, we use the Kinetic Monte Carlo method to investigate the kinetic mechanism for the conformational switches between bistable RNA hairpins. We find three types of conformational switch pathways for RNA hairpins: refolding after complete unfolding, folding through basepair-exchange pathways and through pseudoknot-assisted pathways, respectively. The result of the competition between the three types of pathways depends mainly on the location of the rate-limiting base stacks (such as the GC base stacks) in the structures. Depending on the structural relationships between the two bistable hairpins, the conformational switch can follow single or multiple dominant pathways. The predicted folding pathways are supported by the activation energy results derived from the Arrhenius plot as well as the NMR spectroscopy data.
Collapse
Affiliation(s)
- Xiaojun XU
- Department of Physics and Department of Biochemistry University of Missouri, Columbia, MO 65211
| | - Shi-Jie CHEN
- Department of Physics and Department of Biochemistry University of Missouri, Columbia, MO 65211
| |
Collapse
|
34
|
Nonnative interactions regulate folding and switching of myristoylated protein. Proc Natl Acad Sci U S A 2012; 109:17839-44. [PMID: 22847411 DOI: 10.1073/pnas.1201803109] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We present an integrated experimental and computational study of the molecular mechanisms by which myristoylation affects protein folding and function, which has been little characterized to date. Myristoylation, the covalent linkage of a hydrophobic C14 fatty acyl chain to the N-terminal glycine in a protein, is a common modification that plays a critical role in vital regulated cellular processes by undergoing reversible energetic and conformational switching. Coarse-grained folding simulations for the model pH-dependent actin- and membrane-binding protein hisactophilin reveal that nonnative hydrophobic interactions of the myristoyl with the protein as well as nonnative electrostatic interactions have a pronounced effect on folding rates and thermodynamic stability. Folding measurements for hydrophobic residue mutations of hisactophilin and atomistic simulations indicate that the nonnative interactions of the myristoyl group in the folding transition state are nonspecific and robust, and so smooth the energy landscape for folding. In contrast, myristoyl interactions in the native state are highly specific and tuned for sensitive control of switching functionality. Simulations and amide hydrogen exchange measurements provide evidence for increases as well as decreases in stability localized on one side of the myristoyl binding pocket in the protein, implicating strain and altered dynamics in switching. The effects of folding and function arising from myristoylation are profoundly different from the effects of other post-translational modifications.
Collapse
|