1
|
Zhang S, Wang Z, Qiao J, Yu T, Zhang W. The effect of the loop on the thermodynamic and kinetic of single base pair in pseudoknot. J Chem Phys 2024; 161:085105. [PMID: 39212209 DOI: 10.1063/5.0216593] [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: 04/30/2024] [Accepted: 08/11/2024] [Indexed: 09/04/2024] Open
Abstract
RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases.
Collapse
Affiliation(s)
- Shuhao Zhang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Zhen Wang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Jie Qiao
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Ting Yu
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Wenbing Zhang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| |
Collapse
|
2
|
Sinha S, Pindi C, Ahsan M, Arantes PR, Palermo G. Machines on Genes through the Computational Microscope. J Chem Theory Comput 2023; 19:1945-1964. [PMID: 36947696 PMCID: PMC10104023 DOI: 10.1021/acs.jctc.2c01313] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/24/2023]
Abstract
Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.
Collapse
Affiliation(s)
- Souvik Sinha
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 52512, United States
| | - Chinmai Pindi
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 52512, United States
| | - Mohd Ahsan
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 52512, United States
| | - Pablo R. Arantes
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 52512, United States
| | - Giulia Palermo
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 52512, United States
- Department of Chemistry, University of California Riverside, 900 University Avenue, Riverside, CA 52512, United States
| |
Collapse
|
3
|
Bandyopadhyay D, Mishra PP. Revealing the DNA Unwinding Activity and Mechanism of Fork Reversal by RecG While Exposed to Variants of Stalled Replication-fork at Single-Molecular Resolution. J Mol Biol 2022; 434:167822. [PMID: 36108776 DOI: 10.1016/j.jmb.2022.167822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 08/23/2022] [Accepted: 09/06/2022] [Indexed: 11/25/2022]
Abstract
RecG, belonging to the category of Superfamily-2 plays a vital role in rescuing different kinds of stalled fork. The elemental mechanism of the helicase activity of RecG with several non-homologous stalled fork structures resembling intermediates formed during the process of DNA repair has been investigated in the present study to capture the dynamic stages of genetic rearrangement. The functional characterization has been exemplified through quantifying the response of the substrate in terms of their molecular heterogeneity and dynamical response by employing single-molecule fluorescence methods. An elevated processivity of RecG is observed for the stalled fork where progression of lagging daughter strand is ahead as compared to that of the leading strand. Through precise alteration of its function in terms of unwinding, depending upon the substrate DNA, RecG catalyzes the formation of Holliday junction from a stalled fork DNA. RecG is found to adopt an asymmetric mode of locomotion to unwind the lagging daughter strand for facilitating formation of Holliday junction that acts as a suitable intermediate for recombinational repair pathway. Our results emphasize the mechanism adopted by RecG during its 'sliding back' mode along the lagging daughter strand to be 'active translocation and passive unwinding'. This also provide clues as to how this helicase decides and controls the mode of translocation along the DNA to unwind.
Collapse
Affiliation(s)
- Debolina Bandyopadhyay
- Single Molecule Biophysics Lab, Chemical Sciences Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 700064, India; Homi Bhaba National Institute, Mumbai, India. https://twitter.com/DebolinaBandyo2
| | - Padmaja Prasad Mishra
- Single Molecule Biophysics Lab, Chemical Sciences Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata 700064, India; Homi Bhaba National Institute, Mumbai, India.
| |
Collapse
|
4
|
Measuring thermodynamic preferences to form non-native conformations in nucleic acids using ultraviolet melting. Proc Natl Acad Sci U S A 2022; 119:e2112496119. [PMID: 35671421 PMCID: PMC9214542 DOI: 10.1073/pnas.2112496119] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Thermodynamic preferences to form non-native conformations are crucial for understanding how nucleic acids fold and function. However, they are difficult to measure experimentally because this requires accurately determining the population of minor low-abundance (<10%) conformations in a sea of other conformations. Here, we show that melting experiments enable facile measurements of thermodynamic preferences to adopt nonnative conformations in DNA and RNA. The key to this "delta-melt" approach is to use chemical modifications to render specific minor non-native conformations the major state. The validity and robustness of delta-melt is established for four different non-native conformations under various physiological conditions and sequence contexts through independent measurements of thermodynamic preferences using NMR. Delta-melt is faster relative to NMR, simple, and cost-effective and enables thermodynamic preferences to be measured for exceptionally low-populated conformations. Using delta-melt, we obtained rare insights into conformational cooperativity, obtaining evidence for significant cooperativity (1.0 to 2.5 kcal/mol) when simultaneously forming two adjacent Hoogsteen base pairs. We also measured the thermodynamic preferences to form G-C+ and A-T Hoogsteen and A-T base open states for nearly all 16 trinucleotide sequence contexts and found distinct sequence-specific variations on the order of 2 to 3 kcal/mol. This rich landscape of sequence-specific non-native minor conformations in the DNA double helix may help shape the sequence specificity of DNA biochemistry. Thus, melting experiments can now be used to access thermodynamic information regarding regions of the free energy landscape of biomolecules beyond the native folded and unfolded conformations.
Collapse
|
5
|
Dohnalová H, Lankaš F. Deciphering the mechanical properties of
B‐DNA
duplex. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1575] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| | - Filip Lankaš
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| |
Collapse
|
6
|
Grinevich AA, Masulis IS, Yakushevich LV. Mathematical Modeling of Transcription Bubble Behavior in the pPF1 Plasmid and its Modified Versions: The Link between the Plasmid Energy Profile and the Direction of Transcription. Biophysics (Nagoya-shi) 2021. [DOI: 10.1134/s000635092102007x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
|
7
|
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
|
8
|
Duan Q, Tao P, Wang J, Xiao Y. Molecular dynamics study of ways of RNA base-pair formation. Phys Rev E 2020; 102:032403. [PMID: 33076020 DOI: 10.1103/physreve.102.032403] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2019] [Accepted: 08/21/2020] [Indexed: 06/11/2023]
Abstract
Base pairing is a fundamental phenomenon in RNA structure and function. However, although there have been considerable recent advances, some important aspects of base-pair formation are still unknown, including the ways of base-pair formation and detailed roles of metal ions. Here we show that base pairs can form through four different ways: stabilizing, bridging, rotating, and shifting. Among them the stabilizing and bridging ways involve direct binding of metal ions while the rotating and shifting ways do not in most cases. Furthermore, we find that the formations of base pairs in different positions of the hairpin stem may adopt different ways.
Collapse
Affiliation(s)
- Qiangqiang Duan
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Peng Tao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Jun Wang
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
9
|
Suma A, Coronel L, Bussi G, Micheletti C. Directional translocation resistance of Zika xrRNA. Nat Commun 2020; 11:3749. [PMID: 32719310 PMCID: PMC7385498 DOI: 10.1038/s41467-020-17508-7] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Accepted: 06/29/2020] [Indexed: 12/14/2022] Open
Abstract
xrRNAs from flaviviruses survive in host cells because of their exceptional dichotomic response to the unfolding action of different enzymes. They can be unwound, and hence copied, by replicases, and yet can resist degradation by exonucleases. How the same stretch of xrRNA can encode such diverse responses is an open question. Here, by using atomistic models and translocation simulations, we uncover an elaborate and directional mechanism for how stress propagates when the two xrRNA ends, [Formula: see text] and [Formula: see text], are driven through a pore. Pulling the [Formula: see text] end, as done by replicases, elicits a progressive unfolding; pulling the [Formula: see text] end, as done by exonucleases, triggers a counterintuitive molecular tightening. Thus, in what appears to be a remarkable instance of intra-molecular tensegrity, the very pulling of the [Formula: see text] end is what boosts resistance to translocation and consequently to degradation. The uncovered mechanistic principle might be co-opted to design molecular meta-materials.
Collapse
Affiliation(s)
- Antonio Suma
- Dipartimento di Fisica, Università di Bari and INFN Sezione di Bari, via Amendola 173, 70126, Bari, Italy
- Institute for Computational Molecular Science (ICMS), Temple University, 19122, Philadelphia, PA, Italy
| | - Lucia Coronel
- Physics Area, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136, Trieste, Italy
| | - Giovanni Bussi
- Physics Area, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136, Trieste, Italy
| | - Cristian Micheletti
- Physics Area, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Via Bonomea 265, 34136, Trieste, Italy.
| |
Collapse
|
10
|
Liebl K, Zacharias M. How global DNA unwinding causes non-uniform stress distribution and melting of DNA. PLoS One 2020; 15:e0232976. [PMID: 32413048 PMCID: PMC7228070 DOI: 10.1371/journal.pone.0232976] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 04/24/2020] [Indexed: 12/31/2022] Open
Abstract
DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.
Collapse
Affiliation(s)
- Korbinian Liebl
- Physics Department T38, Technical University of Munich, Garching, Germany
| | - Martin Zacharias
- Physics Department T38, Technical University of Munich, Garching, Germany
- * E-mail:
| |
Collapse
|