1
|
Wu Y, Cao S, Qiu Y, Huang X. Tutorial on how to build non-Markovian dynamic models from molecular dynamics simulations for studying protein conformational changes. J Chem Phys 2024; 160:121501. [PMID: 38516972 DOI: 10.1063/5.0189429] [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: 11/28/2023] [Accepted: 02/20/2024] [Indexed: 03/23/2024] Open
Abstract
Protein conformational changes play crucial roles in their biological functions. In recent years, the Markov State Model (MSM) constructed from extensive Molecular Dynamics (MD) simulations has emerged as a powerful tool for modeling complex protein conformational changes. In MSMs, dynamics are modeled as a sequence of Markovian transitions among metastable conformational states at discrete time intervals (called lag time). A major challenge for MSMs is that the lag time must be long enough to allow transitions among states to become memoryless (or Markovian). However, this lag time is constrained by the length of individual MD simulations available to track these transitions. To address this challenge, we have recently developed Generalized Master Equation (GME)-based approaches, encoding non-Markovian dynamics using a time-dependent memory kernel. In this Tutorial, we introduce the theory behind two recently developed GME-based non-Markovian dynamic models: the quasi-Markov State Model (qMSM) and the Integrative Generalized Master Equation (IGME). We subsequently outline the procedures for constructing these models and provide a step-by-step tutorial on applying qMSM and IGME to study two peptide systems: alanine dipeptide and villin headpiece. This Tutorial is available at https://github.com/xuhuihuang/GME_tutorials. The protocols detailed in this Tutorial aim to be accessible for non-experts interested in studying the biomolecular dynamics using these non-Markovian dynamic models.
Collapse
Affiliation(s)
- Yue Wu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
- Data Science Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
2
|
Wang X, Xu T, Yao Y, Cheung PPH, Gao X, Zhang L. SARS-CoV-2 RNA-Dependent RNA Polymerase Follows Asynchronous Translocation Pathway for Viral Transcription and Replication. J Phys Chem Lett 2023; 14:10119-10128. [PMID: 37922192 DOI: 10.1021/acs.jpclett.3c01249] [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: 11/05/2023]
Abstract
Translocation is one essential step for the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) to exert viral replication and transcription. Although cryo-EM structures of SARS-CoV-2 RdRp are available, the molecular mechanisms of dynamic translocation remain elusive. Herein, we constructed a Markov state model based on extensive molecular dynamics simulations to elucidate the translocation dynamics of the SARS-CoV-2 RdRp. We identified two intermediates that pinpoint the rate-limiting step of translocation and characterize the asynchronous movement of the template-primer duplex. The 3'-terminal nucleotide in the primer strand lags behind due to the uneven distribution of protein-RNA interactions, while the translocation of the template strand is delayed by the hurdle residue K500. Even so, the two strands share the same "ratchet" to stabilize the polymerase in the post-translocation state, suggesting a Brownian-ratchet model. Overall, our study provides intriguing insights into SARS-CoV-2 replication and transcription, which would open a new avenue for drug discoveries.
Collapse
Affiliation(s)
- Xiaowei Wang
- Department of Chemical and Biological Engineering and Department of Mathematics, Hong Kong University of Science and Technology Kowloon, Clear Water Bay, Hong Kong
| | - Tiantian Xu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, Fujian 350002, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yuan Yao
- Department of Chemical and Biological Engineering and Department of Mathematics, Hong Kong University of Science and Technology Kowloon, Clear Water Bay, Hong Kong
| | - Peter Pak-Hang Cheung
- Li Ka Shing Institute of Health Sciences, Department of Chemical Pathology, Chinese University of Hong Kong, 999077, Hong Kong
| | - Xin Gao
- Computational Bioscience Research Center, King Abdullah University of Science and Technology (KAUST), Thuwal 23955, Saudi Arabia
- Computer, Electrical and Mathematical Sciences and Engineering Division, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955, Saudi Arabia
| | - Lu Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, Fujian 350002, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Fuzhou, Fujian 361005, China
| |
Collapse
|
3
|
Unarta IC, Goonetilleke EC, Wang D, Huang X. Nucleotide addition and cleavage by RNA polymerase II: Coordination of two catalytic reactions using a single active site. J Biol Chem 2022; 299:102844. [PMID: 36581202 PMCID: PMC9860460 DOI: 10.1016/j.jbc.2022.102844] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2022] [Revised: 12/19/2022] [Accepted: 12/22/2022] [Indexed: 12/28/2022] Open
Abstract
RNA polymerase II (Pol II) incorporates complementary ribonucleotides into the growing RNA chain one at a time via the nucleotide addition cycle. The nucleotide addition cycle, however, is prone to misincorporation of noncomplementary nucleotides. Thus, to ensure transcriptional fidelity, Pol II backtracks and then cleaves the misincorporated nucleotides. These two reverse reactions, nucleotide addition and cleavage, are catalyzed in the same active site of Pol II, which is different from DNA polymerases or other endonucleases. Recently, substantial progress has been made to understand how Pol II effectively performs its dual role in the same active site. Our review highlights these recent studies and provides an overall model of the catalytic mechanisms of Pol II. In particular, RNA extension follows the two-metal-ion mechanism, and several Pol II residues play important roles to facilitate the catalysis. In sharp contrast, the cleavage reaction is independent of any Pol II residues. Interestingly, Pol II relies on its residues to recognize the misincorporated nucleotides during the backtracking process, prior to cleavage. In this way, Pol II efficiently compartmentalizes its two distinct catalytic functions using the same active site. Lastly, we also discuss a new perspective on the potential third Mg2+ in the nucleotide addition and intrinsic cleavage reactions.
Collapse
Affiliation(s)
- Ilona Christy Unarta
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Eshani C Goonetilleke
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Dong Wang
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, California, USA; Department of Cellular and Molecular Medicine, School of Medicine, University of California, San Diego, La Jolla, California, USA; Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California, USA.
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA.
| |
Collapse
|
4
|
Wang L, Song K, Yu J, Da LT. Computational investigations on target-site searching and recognition mechanisms by thymine DNA glycosylase during DNA repair process. Acta Biochim Biophys Sin (Shanghai) 2022; 54:796-806. [PMID: 35593467 PMCID: PMC9828053 DOI: 10.3724/abbs.2022050] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
DNA glycosylase, as one member of DNA repair machineries, plays an essential role in correcting mismatched/damaged DNA nucleotides by cleaving the N-glycosidic bond between the sugar and target nucleobase through the base excision repair (BER) pathways. Efficient corrections of these DNA lesions are critical for maintaining genome integrity and preventing premature aging and cancers. The target-site searching/recognition mechanisms and the subsequent conformational dynamics of DNA glycosylase, however, remain challenging to be characterized using experimental techniques. In this review, we summarize our recent studies of sequential structural changes of thymine DNA glycosylase (TDG) during the DNA repair process, achieved mostly by molecular dynamics (MD) simulations. Computational simulations allow us to reveal atomic-level structural dynamics of TDG as it approaches the target-site, and pinpoint the key structural elements responsible for regulating the translocation of TDG along DNA. Subsequently, upon locating the lesions, TDG adopts a base-flipping mechanism to extrude the mispaired nucleobase into the enzyme active-site. The constructed kinetic network model elucidates six metastable states during the base-extrusion process and suggests an active role of TDG in flipping the intrahelical nucleobase. Finally, the molecular mechanism of product release dynamics after catalysis is also summarized. Taken together, we highlight to what extent the computational simulations advance our knowledge and understanding of the molecular mechanism underlying the conformational dynamics of TDG, as well as the limitations of current theoretical work.
Collapse
Affiliation(s)
- Lingyan Wang
- Key Laboratory of Systems Biomedicine (Ministry of Education)Shanghai Center for Systems BiomedicineShanghai Jiao Tong UniversityShanghai200240China
| | - Kaiyuan Song
- Key Laboratory of Systems Biomedicine (Ministry of Education)Shanghai Center for Systems BiomedicineShanghai Jiao Tong UniversityShanghai200240China
| | - Jin Yu
- Department of Physics and AstronomyDepartment of ChemistryNSF-Simons Center for Multiscale Cell Fate ResearchUniversity of CaliforniaIrvineCA92697USA
| | - Lin-Tai Da
- Key Laboratory of Systems Biomedicine (Ministry of Education)Shanghai Center for Systems BiomedicineShanghai Jiao Tong UniversityShanghai200240China,Correspondence address. Tel: +86-21-34207348; E-mail:
| |
Collapse
|
5
|
Yuan C, Goonetilleke EC, Unarta IC, Huang X. Incorporation efficiency and inhibition mechanism of 2'-substituted nucleotide analogs against SARS-CoV-2 RNA-dependent RNA polymerase. Phys Chem Chem Phys 2021; 23:20117-20128. [PMID: 34514487 DOI: 10.1039/d1cp03049c] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The ongoing pandemic caused by SARS-CoV-2 emphasizes the need for effective therapeutics. Inhibition of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) by nucleotide analogs provides a promising antiviral strategy. One common group of RdRp inhibitors, 2'-modified nucleotides, are reported to exhibit different behaviors in the SARS-CoV-2 RdRp transcription assay. Three of these analogs, 2'-O-methyl UTP, Sofosbuvir, and 2'-methyl CTP, act as effective inhibitors in previous biochemical experiments, while Gemcitabine and ara-UTP show no inhibitory activity. To understand the impact of the 2'-modification on their inhibitory effects, we conducted extensive molecular dynamics simulations and relative binding free energy calculations using the free energy perturbation method on SARS-CoV-2 replication-transcription complex (RTC) with these five nucleotide analogs. Our results reveal that the five nucleotide analogs display comparable binding affinities to SARS-CoV-2 RdRp and they can all be added to the nascent RNA chain. Moreover, we examine how the incorporation of these nucleotide triphosphate (NTP) analogs will impact the addition of the next nucleotide. Our results indicate that 2'-O-methyl UTP can weaken the binding of the subsequent NTP and consequently lead to partial chain termination. Additionally, Sofosbuvir and 2'-methyl CTP can cause immediate termination due to the strong steric hindrance introduced by their bulky 2'-methyl groups. In contrast, nucleotide analogs with smaller substitutions, such as the fluorine atoms and the ara-hydroxyl group in Gemcitabine and ara-UTP, have a marginal impact on the polymerization process. Our findings are consistent with experimental observations, and more importantly, shed light on the detailed molecular mechanism of SARS-CoV-2 RdRp inhibition by 2'-substituted nucleotide analogs, and may facilitate the rational design of antiviral agents to inhibit SARS-CoV-2 RdRp.
Collapse
Affiliation(s)
- Congmin Yuan
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Eshani C Goonetilleke
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Ilona Christy Unarta
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| |
Collapse
|
6
|
Manigrasso J, De Vivo M, Palermo G. Controlled Trafficking of Multiple and Diverse Cations Prompts Nucleic Acid Hydrolysis. ACS Catal 2021; 11:8786-8797. [PMID: 35145762 DOI: 10.1021/acscatal.1c01825] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Recent in crystallo reaction intermediates have detailed how nucleic acid hydrolysis occurs in the RNA ribonuclease H1 (RNase H1), a fundamental metalloenzyme involved in maintaining the human genome. At odds with the previous characterization, these in crystallo structures unexpectedly captured multiple metal ions (K+ and Mg2+) transiently bound in the vicinity of the two-metal-ion active site. Using multi-microsecond all-atom molecular dynamics and free-energy simulations, we investigated the functional implications of the dynamic exchange of multiple K+ and Mg2+ ions at the RNase H1 reaction center. We found that such ions are timely positioned at non-overlapping locations near the active site, at different stages of catalysis, being crucial for both reactants' alignment and leaving group departure. We also found that this cation trafficking is tightly regulated by variations of the solution's ionic strength and is aided by two conserved second-shell residues, E188 and K196, suggesting a mechanism for the cations' recruitment during catalysis. These results indicate that controlled trafficking of multi-cation dynamics, opportunely prompted by second-shell residues, is functionally essential to the complex enzymatic machinery of the RNase H1. These findings revise the current knowledge on the RNase H1 catalysis and open new catalytic possibilities for other similar metalloenzymes including, but not limited to, CRISPR-Cas9, group II intron ribozyme and the human spliceosome.
Collapse
Affiliation(s)
- Jacopo Manigrasso
- Laboratory of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Genoa, 16163, Italy.,Department of Bioengineering, University of California Riverside, Riverside, CA 52512, United States
| | - Marco De Vivo
- Laboratory of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Genoa, 16163, Italy
| | - Giulia Palermo
- Department of Bioengineering, University of California Riverside, Riverside, CA 52512, United States.,Department of Chemistry, University of California Riverside, Riverside, CA 52512, United States
| |
Collapse
|
7
|
A comprehensive mechanism for 5-carboxylcytosine-induced transcriptional pausing revealed by Markov state models. J Biol Chem 2021; 296:100735. [PMID: 33991521 PMCID: PMC8191312 DOI: 10.1016/j.jbc.2021.100735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Revised: 04/27/2021] [Accepted: 04/28/2021] [Indexed: 11/23/2022] Open
Abstract
RNA polymerase II (Pol II) surveils the genome, pausing as it encounters DNA lesions and base modifications and initiating signals for DNA repair among other important regulatory events. Recent work suggests that Pol II pauses at 5-carboxycytosine (5caC), an epigenetic modification of cytosine, because of a specific hydrogen bond between the carboxyl group of 5caC and a specific residue in fork loop 3 of Pol II. This hydrogen bond compromises productive NTP binding and slows down elongation. Apart from this specific interaction, the carboxyl group of 5caC can potentially interact with numerous charged residues in the cleft of Pol II. However, it is not clear how other interactions between Pol II and 5caC contribute to pausing. In this study, we use Markov state models (a type of kinetic network models) built from extensive molecular dynamics simulations to comprehensively study the impact of 5caC on Pol II translocation. We describe two translocation intermediates with specific interactions that prevent the template base from loading into the Pol II active site. In addition to the previously observed state with 5caC constrained by fork loop 3, we discovered a new intermediate state with a hydrogen bond between 5caC and fork loop 2. Surprisingly, we find that 5caC may curb translocation by suppressing kinking of the helix bordering the active site (the bridge helix) because its high flexibility is critical to translocation. Our work provides new insights into how epigenetic modifications of genomic DNA can modulate Pol II translocation, inducing pauses in transcription.
Collapse
|
8
|
Geronimo I, Vidossich P, Donati E, Vivo M. Computational investigations of polymerase enzymes: Structure, function, inhibition, and biotechnology. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1534] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Inacrist Geronimo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| | - Pietro Vidossich
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| | - Elisa Donati
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| | - Marco Vivo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia Genoa Italy
| |
Collapse
|
9
|
Role of bacterial RNA polymerase gate opening dynamics in DNA loading and antibiotics inhibition elucidated by quasi-Markov State Model. Proc Natl Acad Sci U S A 2021; 118:2024324118. [PMID: 33883282 DOI: 10.1073/pnas.2024324118] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
To initiate transcription, the holoenzyme (RNA polymerase [RNAP] in complex with σ factor) loads the promoter DNA via the flexible loading gate created by the clamp and β-lobe, yet their roles in DNA loading have not been characterized. We used a quasi-Markov State Model (qMSM) built from extensive molecular dynamics simulations to elucidate the dynamics of Thermus aquaticus holoenzyme's gate opening. We showed that during gate opening, β-lobe oscillates four orders of magnitude faster than the clamp, whose opening depends on the Switch 2's structure. Myxopyronin, an antibiotic that binds to Switch 2, was shown to undergo a conformational selection mechanism to inhibit clamp opening. Importantly, we reveal a critical but undiscovered role of β-lobe, whose opening is sufficient for DNA loading even when the clamp is partially closed. These findings open the opportunity for the development of antibiotics targeting β-lobe of RNAP. Finally, we have shown that our qMSMs, which encode non-Markovian dynamics based on the generalized master equation formalism, hold great potential to be widely applied to study biomolecular dynamics.
Collapse
|
10
|
Wang X, Unarta IC, Cheung PPH, Huang X. Elucidating molecular mechanisms of functional conformational changes of proteins via Markov state models. Curr Opin Struct Biol 2020; 67:69-77. [PMID: 33126140 DOI: 10.1016/j.sbi.2020.10.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Revised: 09/28/2020] [Accepted: 10/07/2020] [Indexed: 01/01/2023]
Abstract
Functional conformational changes of proteins can facilitate numerous biological events in cells. The Markov state model (MSM) built from molecular dynamics simulations provide a powerful approach to study them. We here introduce a protocol that is tailor-made for constructing MSMs to study the functional conformational changes of proteins. In this protocol, one of the important steps is to select proper molecular features that can collectively describe the slowest timescales of conformational changes of interest. We recommend spectral oASIS, the modified version of oASIS, as a promising approach for automatic feature selection. Recently developed deep learning methods could also serve efficient approaches for selecting features and finding collective variables. Using DNA repair enzymes and RNA polymerases as examples, we review recent applications of MSMs to elucidate molecular mechanisms of functional conformational changes. Finally, we discuss remaining challenges and future perspectives for constructing MSMs to study functional conformational changes of proteins.
Collapse
Affiliation(s)
- Xiaowei Wang
- The Hong Kong University of Science and Technology-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China; Department of Chemistry, Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Ilona Christy Unarta
- Bioengineering Graduate Program, The Hong Kong University of Science and Technology, Kowloon, 4Hong Kong Center for Neurodegenerative Diseases, Hong Kong
| | - Peter Pak-Hang Cheung
- The Hong Kong University of Science and Technology-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China; Department of Chemistry, Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Xuhui Huang
- The Hong Kong University of Science and Technology-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China; Department of Chemistry, Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong; Bioengineering Graduate Program, The Hong Kong University of Science and Technology, Kowloon, 4Hong Kong Center for Neurodegenerative Diseases, Hong Kong.
| |
Collapse
|
11
|
Belogurov GA, Artsimovitch I. The Mechanisms of Substrate Selection, Catalysis, and Translocation by the Elongating RNA Polymerase. J Mol Biol 2019; 431:3975-4006. [PMID: 31153902 DOI: 10.1016/j.jmb.2019.05.042] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2019] [Revised: 05/24/2019] [Accepted: 05/24/2019] [Indexed: 11/15/2022]
Abstract
Multi-subunit DNA-dependent RNA polymerases synthesize all classes of cellular RNAs, ranging from short regulatory transcripts to gigantic messenger RNAs. RNA polymerase has to make each RNA product in just one try, even if it takes millions of successive nucleotide addition steps. During each step, RNA polymerase selects a correct substrate, adds it to a growing chain, and moves one nucleotide forward before repeating the cycle. However, RNA synthesis is anything but monotonous: RNA polymerase frequently pauses upon encountering mechanical, chemical and torsional barriers, sometimes stepping back and cleaving off nucleotides from the growing RNA chain. A picture in which these intermittent dynamics enable processive, accurate, and controllable RNA synthesis is emerging from complementary structural, biochemical, computational, and single-molecule studies. Here, we summarize our current understanding of the mechanism and regulation of the on-pathway transcription elongation. We review the details of substrate selection, catalysis, proofreading, and translocation, focusing on rate-limiting steps, structural elements that modulate them, and accessory proteins that appear to control RNA polymerase translocation.
Collapse
Affiliation(s)
| | - Irina Artsimovitch
- Department of Microbiology and The Center for RNA Biology, The Ohio State University, Columbus, OH, USA.
| |
Collapse
|
12
|
Long C, E. C, Da LT, Yu J. A Viral T7 RNA Polymerase Ratcheting Along DNA With Fidelity Control. Comput Struct Biotechnol J 2019; 17:638-644. [PMID: 31193497 PMCID: PMC6535458 DOI: 10.1016/j.csbj.2019.05.001] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 04/25/2019] [Accepted: 05/04/2019] [Indexed: 12/02/2022] Open
Abstract
RNA polymerase (RNAP) from bacteriophage T7 is a representative single-subunit viral RNAP that can transcribe with high promoter activities without assistances from transcription factors. We accordingly studied this small transcription machine computationally as a model system to understand underlying mechanisms of mechano-chemical coupling and fidelity control in the RNAP transcription elongation. Here we summarize our computational work from several recent publications to demonstrate first how T7 RNAP translocates via Brownian alike motions along DNA right after the catalytic product release. Then we show how the backward translocation motions are prevented at post-translocation upon successful nucleotide incorporation, which is also subject to stepwise nucleotide selection and acts as a pawl for "selective ratcheting". The structural dynamics and energetics features revealed from our atomistic molecular dynamics (MD) simulations and related analyses on the single-subunit T7 RNAP thus provided detailed and quantitative characterizations on the Brownian-ratchet working scenario of a prototypical transcription machine with sophisticated nucleotide selectivity for fidelity control. The presented mechanisms can be more or less general for structurally similar viral or mitochondrial RNAPs and some of DNA polymerases, or even for the RNAP engine of the more complicated transcription machinery in higher organisms.
Collapse
Affiliation(s)
- Chunhong Long
- Beijing Computational Science Research Center, Beijing, 100193, China
| | - Chao E.
- Beijing Computational Science Research Center, Beijing, 100193, China
| | - Lin-Tai Da
- Shanghai Center for Systems Biomedicine, Shanghai JiaoTong University, Shanghai 200240, China
| | - Jin Yu
- Beijing Computational Science Research Center, Beijing, 100193, China
| |
Collapse
|
13
|
Zhu L, Sheong FK, Cao S, Liu S, Unarta IC, Huang X. TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromolecules. J Chem Phys 2019; 150:124105. [DOI: 10.1063/1.5082633] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Lizhe Zhu
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen, Guangdong 518172, China
| | - Fu Kit Sheong
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Siqin Cao
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Song Liu
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Ilona C. Unarta
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Bioengineering Program, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- HKUST-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China
| |
Collapse
|
14
|
Vithani N, Ankush Jagtap PK, Verma SK, Tripathi R, Awasthi S, Nair NN, Prakash B. Mechanism of Mg 2+-Accompanied Product Release in Sugar Nucleotidyltransferases. Structure 2019. [PMID: 29514078 DOI: 10.1016/j.str.2018.02.002] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The nucleotidyl transfer reaction, catalyzed by sugar nucleotidyltransferases (SNTs), is assisted by two active site Mg2+ ions. While studying this reaction using X-ray crystallography, we captured snapshots of the pyrophosphate (product) as it exits along a pocket. Surprisingly, one of the active site Mg2+ ions remains coordinated to the exiting pyrophosphate. This hints at the participation of Mg2+ in the process of product release, besides its role in catalyzing nucleotidyl transfer. These observations are further supported by enhanced sampling molecular dynamics simulations. Free energy computations suggest that the product release is likely to be rate limiting in SNTs, and the origin of the high free energy barrier for product release could be traced back to the "slow" conformational change of an Arg residue at the exit end of the pocket. These results establish a dual role for Mg2+, and propose a general mechanism of product release during the nucleotidyl transfer by SNTs.
Collapse
Affiliation(s)
- Neha Vithani
- Department of Biological Sciences and Bioengineering, Indian Institute of Technology, Kanpur 208016, India
| | - Pravin Kumar Ankush Jagtap
- Department of Biological Sciences and Bioengineering, Indian Institute of Technology, Kanpur 208016, India
| | - Sunil Kumar Verma
- Department of Biological Sciences and Bioengineering, Indian Institute of Technology, Kanpur 208016, India
| | - Ravi Tripathi
- Department of Chemistry, Indian Institute of Technology, Kanpur 208016, India
| | - Shalini Awasthi
- Department of Chemistry, Indian Institute of Technology, Kanpur 208016, India
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology, Kanpur 208016, India.
| | - Balaji Prakash
- Department of Molecular Nutrition, CSIR-Central Food Technological Research Institute, Mysore 570020, India.
| |
Collapse
|
15
|
Affiliation(s)
- Vito Genna
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Elisa Donati
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| |
Collapse
|
16
|
Š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
|
17
|
Zeng X, Li ZW, Zheng X, Zhu L, Sun ZY, Lu ZY, Huang X. Improving the productivity of monodisperse polyhedral cages by the rational design of kinetic self-assembly pathways. Phys Chem Chem Phys 2018; 20:10030-10037. [PMID: 29620122 DOI: 10.1039/c8cp00522b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Hollow polyhedral cages hold great potential for application in nanotechnological and biomedical fields. Understanding the formation mechanism of these self-assembled structures could provide guidance for the rational design of the desired polyhedral cages. Here, by constructing kinetic network models from extensive coarse-grained molecular dynamics simulations, we elucidated the formation mechanism of the dodecahedral cage, which is formed by the self-assembly of patchy particles. We found that the dodecahedral cage is formed through increasing the aggregate size followed by structure rearrangement. Based on this mechanistic understanding, we improved the productivity of the dodecahedral cage through the rational design of the patch arrangement of patchy particles, which promotes the structural rearrangement process. Our results demonstrate that it should be a feasible strategy to achieve the rational design of the desired nanostructures via the kinetic analysis. We anticipate that this methodology could be extended to other self-assembly systems for the fabrication of functional nanomaterials.
Collapse
Affiliation(s)
- Xiangze Zeng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | | | | | | | | | | | | |
Collapse
|
18
|
Zeng X, Zhu L, Zheng X, Cecchini M, Huang X. Harnessing complexity in molecular self-assembly using computer simulations. Phys Chem Chem Phys 2018; 20:6767-6776. [PMID: 29479585 DOI: 10.1039/c7cp06181a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In molecular self-assembly, hundreds of thousands of freely-diffusing molecules associate to form ordered and functional architectures in the absence of an actuator. This intriguing phenomenon plays a critical role in biology and has become a powerful tool for the fabrication of advanced nanomaterials. Due to the limited spatial and temporal resolutions of current experimental techniques, computer simulations offer a complementary strategy to explore self-assembly with atomic resolution. Here, we review recent computational studies focusing on both thermodynamic and kinetic aspects. As we shall see, thermodynamic approaches based on modeling and statistical mechanics offer initial guidelines to design nanostructures with modest computational effort. Computationally more intensive analyses based on molecular dynamics simulations and kinetic network models (KNMs) reach beyond it, opening the door to the rational design of self-assembly pathways. Current limitations of these methodologies are discussed. We anticipate that the synergistic use of thermodynamic and kinetic analyses based on computer simulations will provide an important contribution to the de novo design of self-assembly.
Collapse
Affiliation(s)
- Xiangze Zeng
- Department of Chemistry, Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | | | | | | | | |
Collapse
|
19
|
Zhu L, Sheong FK, Zeng X, Huang X. Elucidation of the conformational dynamics of multi-body systems by construction of Markov state models. Phys Chem Chem Phys 2018; 18:30228-30235. [PMID: 27314275 DOI: 10.1039/c6cp02545e] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Constructing Markov State Models (MSMs) based on short molecular dynamics simulations is a powerful computational technique to complement experiments in predicting long-time kinetics of biomolecular processes at atomic resolution. Even though the MSM approach has been widely applied to study one-body processes such as protein folding and enzyme conformational changes, the majority of biological processes, e.g. protein-ligand recognition, signal transduction, and protein aggregation, essentially involve multiple entities. Here we review the attempts at constructing MSMs for multi-body systems, point out the challenges therein and discuss recent algorithmic progresses that alleviate these challenges. In particular, we describe an automatic kinetics based partitioning method that achieves optimal definition of the conformational states in a multi-body system, and discuss a novel maximum-likelihood approach that efficiently estimates the slow uphill kinetics utilizing pre-computed equilibrium populations of all states. We expect that these new algorithms and their combinations may boost investigations of important multi-body biological processes via the efficient construction of MSMs.
Collapse
Affiliation(s)
- Lizhe Zhu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Fu Kit Sheong
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China.
| | - Xiangze Zeng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| |
Collapse
|
20
|
Wang W, Cao S, Zhu L, Huang X. Constructing Markov State Models to elucidate the functional conformational changes of complex biomolecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1343] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Wei Wang
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Center of Systems Biology and Human HealthThe Hong Kong University of Science and Technology Kowloon Hong Kong
| | - Siqin Cao
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
| | - Lizhe Zhu
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Center of Systems Biology and Human HealthThe Hong Kong University of Science and Technology Kowloon Hong Kong
| | - Xuhui Huang
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Center of Systems Biology and Human HealthThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & ReconstructionThe Hong Kong University of Science and Technology Kowloon Hong Kong
- HKUST‐Shenzhen Research Institute Shenzhen China
| |
Collapse
|
21
|
Wu S, Li L, Li Q. Mechanism of NTP Binding to the Active Site of T7 RNA Polymerase Revealed by Free-Energy Simulation. Biophys J 2017; 112:2253-2260. [PMID: 28591598 PMCID: PMC5474740 DOI: 10.1016/j.bpj.2017.04.039] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Revised: 04/11/2017] [Accepted: 04/24/2017] [Indexed: 11/26/2022] Open
Abstract
In genetic transcription, molecular dynamic details and energetics of NTP binding to the active site of RNA polymerase (RNAP) are poorly understood. In this article, we investigated the NTP binding process in T7 RNAP using all-atom MD simulation combined with the umbrella sampling technique. Based on our simulations, a two-step mechanism was proposed to explain NTP binding: first, substrate NTP in aqueous solution, which carries a magnesium ion, diffuses through a secondary channel of RNAP to attain a pore region, where it undergoes conformational changes to give a correct orientation; next, the NTP establishes initial basepairing contacts with the template nucleoside (TN). Our free-energy calculations suggest that both steps are spontaneous. This mechanism can easily explain the problem of NTP binding with different orientations. Moreover, it is found that the nascent NTP:TN basepair is fragile and easily broken by thermal disturbance. Therefore, we speculate that the fingers domain will be triggered to close, so as to create a steady environment for the next chemical step. The observations from the work provide valuable information for comprehensively understanding the mechanism of the basic step in genetic transcription.
Collapse
Affiliation(s)
- Shaogui Wu
- College of Chemistry and Materials Science, Sichuan Normal University, Chengdu, China; State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China.
| | - Laicai Li
- College of Chemistry and Materials Science, Sichuan Normal University, Chengdu, China
| | - Quan Li
- College of Chemistry and Materials Science, Sichuan Normal University, Chengdu, China
| |
Collapse
|
22
|
Wang B, Sexton RE, Feig M. Kinetics of nucleotide entry into RNA polymerase active site provides mechanism for efficiency and fidelity. BIOCHIMICA ET BIOPHYSICA ACTA. GENE REGULATORY MECHANISMS 2017; 1860:482-490. [PMID: 28242207 PMCID: PMC5393355 DOI: 10.1016/j.bbagrm.2017.02.008] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Revised: 01/29/2017] [Accepted: 02/15/2017] [Indexed: 10/20/2022]
Abstract
During transcription, RNA polymerase II elongates RNA by adding nucleotide triphosphates (NTPs) complementary to a DNA template. Structural studies have suggested that NTPs enter and exit the active site via the narrow secondary pore but details have remained unclear. A kinetic model is presented that integrates molecular dynamics simulations with experimental data. Previous simulations of trigger loop dynamics and the dynamics of matched and mismatched NTPs in and near the active site were combined with new simulations describing NTP exit from the active site via the secondary pore. Markov state analysis was applied to identify major states and estimate kinetic rates for transitions between those states. The kinetic model predicts elongation and misincorporation rates in close agreement with experiment and provides mechanistic hypotheses for how NTP entry and exit via the secondary pore is feasible and a key feature for achieving high elongation and low misincorporation rates during RNA elongation.
Collapse
Affiliation(s)
- Beibei Wang
- Department of Biochemistry & Molecular Biology, 603 Wilson Rd., Room 218 BCH, Michigan State University, East Lansing, MI 48824, USA.
| | - Rachel E Sexton
- Department of Biochemistry & Molecular Biology, 603 Wilson Rd., Room 218 BCH, Michigan State University, East Lansing, MI 48824, USA.
| | - Michael Feig
- Department of Biochemistry & Molecular Biology, 603 Wilson Rd., Room 218 BCH, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
23
|
Elucidating Mechanisms of Molecular Recognition Between Human Argonaute and miRNA Using Computational Approaches. Methods Mol Biol 2016. [PMID: 27924488 DOI: 10.1007/978-1-4939-6563-2_18] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
MicroRNA (miRNA) and Argonaute (AGO) protein together form the RNA-induced silencing complex (RISC) that plays an essential role in the regulation of gene expression. Elucidating the underlying mechanism of AGO-miRNA recognition is thus of great importance not only for the in-depth understanding of miRNA function but also for inspiring new drugs targeting miRNAs. In this chapter we introduce a combined computational approach of molecular dynamics (MD) simulations, Markov state models (MSMs), and protein-RNA docking to investigate AGO-miRNA recognition. Constructed from MD simulations, MSMs can elucidate the conformational dynamics of AGO at biologically relevant timescales. Protein-RNA docking can then efficiently identify the AGO conformations that are geometrically accessible to miRNA. Using our recent work on human AGO2 as an example, we explain the rationale and the workflow of our method in details. This combined approach holds great promise to complement experiments in unraveling the mechanisms of molecular recognition between large, flexible, and complex biomolecules.
Collapse
|
24
|
Liu S, Zhu L, Sheong FK, Wang W, Huang X. Adaptive partitioning by local density-peaks: An efficient density-based clustering algorithm for analyzing molecular dynamics trajectories. J Comput Chem 2016; 38:152-160. [PMID: 27868222 DOI: 10.1002/jcc.24664] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2016] [Revised: 10/09/2016] [Accepted: 10/26/2016] [Indexed: 12/11/2022]
Abstract
We present an efficient density-based adaptive-resolution clustering method APLoD for analyzing large-scale molecular dynamics (MD) trajectories. APLoD performs the k-nearest-neighbors search to estimate the density of MD conformations in a local fashion, which can group MD conformations in the same high-density region into a cluster. APLoD greatly improves the popular density peaks algorithm by reducing the running time and the memory usage by 2-3 orders of magnitude for systems ranging from alanine dipeptide to a 370-residue Maltose-binding protein. In addition, we demonstrate that APLoD can produce clusters with various sizes that are adaptive to the underlying density (i.e., larger clusters at low-density regions, while smaller clusters at high-density regions), which is a clear advantage over other popular clustering algorithms including k-centers and k-medoids. We anticipate that APLoD can be widely applied to split ultra-large MD datasets containing millions of conformations for subsequent construction of Markov State Models. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Song Liu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Lizhe Zhu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.,Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Fu Kit Sheong
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Wei Wang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.,Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.,Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| |
Collapse
|
25
|
Zhu L, Jiang H, Sheong FK, Cui X, Wang Y, Gao X, Huang X. Understanding the core of RNA interference: The dynamic aspects of Argonaute-mediated processes. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2016; 128:39-46. [PMID: 27697475 DOI: 10.1016/j.pbiomolbio.2016.09.008] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2016] [Revised: 09/04/2016] [Accepted: 09/26/2016] [Indexed: 12/14/2022]
Abstract
At the core of RNA interference, the Argonaute proteins (Ago) load and utilize small guide nucleic acids to silence mRNAs or cleave foreign nucleic acids in a sequence specific manner. In recent years, based on extensive structural studies of Ago and its interaction with the nucleic acids, considerable progress has been made to reveal the dynamic aspects of various Ago-mediated processes. Here we review these novel insights into the guide-strand loading, duplex unwinding, and effects of seed mismatch, with a focus on two representative Agos, the human Ago 2 (hAgo2) and the bacterial Thermus thermophilus Ago (TtAgo). In particular, comprehensive molecular simulation studies revealed that although sharing similar overall structures, the two Agos have vastly different conformational landscapes and guide-strand loading mechanisms because of the distinct rigidity of their L1-PAZ hinge. Given the central role of the PAZ motions in regulating the exposure of the nucleic acid binding channel, these findings exemplify the importance of protein motions in distinguishing the overlapping, yet distinct, mechanisms of Ago-mediated processes in different organisms.
Collapse
Affiliation(s)
- Lizhe Zhu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong; Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Hanlun Jiang
- Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong; Bioengineering Graduate Program, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Fu Kit Sheong
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuefeng Cui
- King Abdullah University of Science and Technology (KAUST), Computational Bioscience Research Center (CBRC), Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, Thuwal, 23955, Saudi Arabia
| | - Yanli Wang
- Laboratory of Non-Coding RNA, Institute of Biophysics, Chinese Academy of Sciences, Beijing, 100101, China
| | - Xin Gao
- King Abdullah University of Science and Technology (KAUST), Computational Bioscience Research Center (CBRC), Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, Thuwal, 23955, Saudi Arabia
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong; Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong; Bioengineering Graduate Program, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| |
Collapse
|
26
|
Gunnoo M, Cazade PA, Galera-Prat A, Nash MA, Czjzek M, Cieplak M, Alvarez B, Aguilar M, Karpol A, Gaub H, Carrión-Vázquez M, Bayer EA, Thompson D. Nanoscale Engineering of Designer Cellulosomes. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2016; 28:5619-47. [PMID: 26748482 DOI: 10.1002/adma.201503948] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Revised: 10/01/2015] [Indexed: 05/27/2023]
Abstract
Biocatalysts showcase the upper limit obtainable for high-speed molecular processing and transformation. Efforts to engineer functionality in synthetic nanostructured materials are guided by the increasing knowledge of evolving architectures, which enable controlled molecular motion and precise molecular recognition. The cellulosome is a biological nanomachine, which, as a fundamental component of the plant-digestion machinery from bacterial cells, has a key potential role in the successful development of environmentally-friendly processes to produce biofuels and fine chemicals from the breakdown of biomass waste. Here, the progress toward so-called "designer cellulosomes", which provide an elegant alternative to enzyme cocktails for lignocellulose breakdown, is reviewed. Particular attention is paid to rational design via computational modeling coupled with nanoscale characterization and engineering tools. Remaining challenges and potential routes to industrial application are put forward.
Collapse
Affiliation(s)
- Melissabye Gunnoo
- Materials and Surface Science Institute and Department of Physics and Energy, University of Limerick, Limerick, Ireland
| | - Pierre-André Cazade
- Materials and Surface Science Institute and Department of Physics and Energy, University of Limerick, Limerick, Ireland
| | - Albert Galera-Prat
- Instituto Cajal, Consejo Superior de Investigaciones Cientificas (CSIC), IMDEA Nanociencias and CIBERNED, Madrid, Spain
| | - Michael A Nash
- Lehrstuhl für Angewandte Physik and Center for Nanoscience, Ludwig-Maximilians-University, 80799, Munich, Germany
| | - Mirjam Czjzek
- Sorbonne Universités, UPMC, Université Paris 06, and Centre National de la Recherche Scientifique, UMR 8227, Integrative Biology of Marine Models, Station Biologique, de Roscoff, CS 90074, F-29688, Roscoff cedex, Bretagne, France
| | - Marek Cieplak
- Laboratory of Biological Physics, Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| | - Beatriz Alvarez
- Biopolis S.L., Parc Científic de la Universitat de Valencia, Edificio 2, C/Catedrático Agustín Escardino 9, 46980, Paterna (Valencia), Spain
| | - Marina Aguilar
- Abengoa, S.A., Palmas Altas, Calle Energía Solar nº 1, 41014, Seville, Spain
| | - Alon Karpol
- Designer Energy Ltd., 2 Bergman St., Tamar Science Park, Rehovot, 7670504, Israel
| | - Hermann Gaub
- Lehrstuhl für Angewandte Physik and Center for Nanoscience, Ludwig-Maximilians-University, 80799, Munich, Germany
| | - Mariano Carrión-Vázquez
- Instituto Cajal, Consejo Superior de Investigaciones Cientificas (CSIC), IMDEA Nanociencias and CIBERNED, Madrid, Spain
| | - Edward A Bayer
- Department of Biological Chemistry, Weizmann Institute of Science, Rehovot, 76100, Israel
| | - Damien Thompson
- Materials and Surface Science Institute and Department of Physics and Energy, University of Limerick, Limerick, Ireland
| |
Collapse
|
27
|
Zhang L, Pardo-Avila F, Unarta IC, Cheung PPH, Wang G, Wang D, Huang X. Elucidation of the Dynamics of Transcription Elongation by RNA Polymerase II using Kinetic Network Models. Acc Chem Res 2016; 49:687-94. [PMID: 26991064 DOI: 10.1021/acs.accounts.5b00536] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
RNA polymerase II (Pol II) is an essential enzyme that catalyzes transcription with high efficiency and fidelity in eukaryotic cells. During transcription elongation, Pol II catalyzes the nucleotide addition cycle (NAC) to synthesize mRNA using DNA as the template. The transitions between the states of the NAC require conformational changes of both the protein and nucleotides. Although X-ray structures are available for most of these states, the dynamics of the transitions between states are largely unknown. Molecular dynamics (MD) simulations can predict structure-based molecular details and shed light on the mechanisms of these dynamic transitions. However, the employment of MD simulations on a macromolecule (tens to hundreds of nanoseconds) such as Pol II is challenging due to the difficulty of reaching biologically relevant timescales (tens of microseconds or even longer). For this challenge to be overcome, kinetic network models (KNMs), such as Markov State Models (MSMs), have become a popular approach to access long-timescale conformational changes using many short MD simulations. We describe here our application of KNMs to characterize the molecular mechanisms of the NAC of Pol II. First, we introduce the general background of MSMs and further explain procedures for the construction and validation of MSMs by providing some technical details. Next, we review our previous studies in which we applied MSMs to investigate the individual steps of the NAC, including translocation and pyrophosphate ion release. In particular, we describe in detail how we prepared the initial conformations of Pol II elongation complex, performed MD simulations, extracted MD conformations to construct MSMs, and further validated them. We also summarize our major findings on molecular mechanisms of Pol II elongation based on these MSMs. In addition, we have included discussions regarding various key points and challenges for applications of MSMs to systems as large as the Pol II elongation complex. Finally, to study the overall NAC, we combine the individual steps of the NAC into a five-state KNM based on a nonbranched Brownian ratchet scheme to explain the single-molecule optical tweezers experimental data. The studies complement experimental observations and provide molecular mechanisms for the transcription elongation cycle. In the long term, incorporation of sequence-dependent kinetic parameters into KNMs has great potential for identifying error-prone sequences and predicting transcription dynamics in genome-wide transcriptomes.
Collapse
Affiliation(s)
- Lu Zhang
- Department
of Chemistry and State Key Laboratory of Molecular Neuroscience, Center
for System Biology and Human Health, School of Science, and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Fátima Pardo-Avila
- Department
of Chemistry and State Key Laboratory of Molecular Neuroscience, Center
for System Biology and Human Health, School of Science, and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Ilona Christy Unarta
- Department
of Chemistry and State Key Laboratory of Molecular Neuroscience, Center
for System Biology and Human Health, School of Science, and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Peter Pak-Hang Cheung
- Department
of Chemistry and State Key Laboratory of Molecular Neuroscience, Center
for System Biology and Human Health, School of Science, and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Guo Wang
- Department
of Chemistry and State Key Laboratory of Molecular Neuroscience, Center
for System Biology and Human Health, School of Science, and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Dong Wang
- Department
of Cellular and Molecular Medicine, Skaggs School of Pharmacy and
Pharmaceutical Sciences, University of California, San Diego, La Jolla, California 92093, United States
| | - Xuhui Huang
- Department
of Chemistry and State Key Laboratory of Molecular Neuroscience, Center
for System Biology and Human Health, School of Science, and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| |
Collapse
|
28
|
Zeng X, Zhang L, Xiao X, Jiang Y, Guo Y, Yu X, Pu X, Li M. Unfolding mechanism of thrombin-binding aptamer revealed by molecular dynamics simulation and Markov State Model. Sci Rep 2016; 6:24065. [PMID: 27045335 PMCID: PMC4820715 DOI: 10.1038/srep24065] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2016] [Accepted: 03/16/2016] [Indexed: 12/27/2022] Open
Abstract
Thrombin-binding aptamer (TBA) with the sequence 5′GGTTGGTGTGGTTGG3′ could fold into G-quadruplex, which correlates with functionally important genomic regionsis. However, unfolding mechanism involved in the structural stability of G-quadruplex has not been satisfactorily elucidated on experiments so far. Herein, we studied the unfolding pathway of TBA by a combination of molecular dynamics simulation (MD) and Markov State Model (MSM). Our results revealed that the unfolding of TBA is not a simple two-state process but proceeds along multiple pathways with multistate intermediates. One high flux confirms some observations from NMR experiment. Another high flux exhibits a different and simpler unfolding pathway with less intermediates. Two important intermediate states were identified. One is similar to the G-triplex reported in the folding of G-quadruplex, but lack of H-bonding between guanines in the upper plane. More importantly, another intermediate state acting as a connector to link the folding region and the unfolding one, was the first time identified, which exhibits higher population and stability than the G-triplex-like intermediate. These results will provide valuable information for extending our understanding the folding landscape of G-quadruplex formation.
Collapse
Affiliation(s)
- Xiaojun Zeng
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Liyun Zhang
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Xiuchan Xiao
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Yuanyuan Jiang
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Yanzhi Guo
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Xinyan Yu
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Xuemei Pu
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| | - Menglong Li
- Faculty of Chemistry, Sichuan University, Chengdu 610064, People's Republic of China
| |
Collapse
|
29
|
Zhang L, Jiang H, Sheong F, Pardo-Avila F, Cheung PH, Huang X. Constructing Kinetic Network Models to Elucidate Mechanisms of Functional Conformational Changes of Enzymes and Their Recognition with Ligands. Methods Enzymol 2016; 578:343-71. [DOI: 10.1016/bs.mie.2016.05.026] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
30
|
Shen H, Li G. Bridging the Missing Link between Structure and Fidelity of the RNA-Dependent RNA Polymerase from Poliovirus through Free Energy Simulations. J Chem Theory Comput 2015; 10:5195-205. [PMID: 26584391 DOI: 10.1021/ct5006449] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
RNA-dependent RNA polymerases (RdRps) are enzymes catalyzing RNA replication from a RNA template. Active-site closure in RdRps, normally induced by correct nucleotide triphosphate (NTP) binding, is a prerequisite for the cycle of nucleotide incorporation. So, a complete understanding of polymerase function (in particular polymerase fidelity) of a RdRp requires more complete knowledge of active-site closure in the RdRp. In this work, based on solved crystal structures, we have built different models for the RNA-dependent RNA polymerase from poliovirus (termed PV 3D(pol)). Through MD simulations and free energy calculations of these PV 3D(pol) models, we have revealed the dynamic correlation between motif A and motif D and between motif A and incoming NTP, have deepened our understanding of polymerase fidelity from dynamic aspects, and have provided an explanation to the puzzle that arises from different observations based on kinetic studies and structural data.
Collapse
Affiliation(s)
- Hujun Shen
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, China
| |
Collapse
|
31
|
Da LT, E C, Duan B, Zhang C, Zhou X, Yu J. A Jump-from-Cavity Pyrophosphate Ion Release Assisted by a Key Lysine Residue in T7 RNA Polymerase Transcription Elongation. PLoS Comput Biol 2015; 11:e1004624. [PMID: 26599007 PMCID: PMC4658072 DOI: 10.1371/journal.pcbi.1004624] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2015] [Accepted: 10/26/2015] [Indexed: 11/18/2022] Open
Abstract
Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucleotide addition cycle. The kinetics and energetics of the process as well as how it proceeds with substantial conformational changes of the polymerase complex determine the mechano-chemical coupling mechanism of the transcription elongation. Here we investigated detailed dynamics of the PPi release process in a single-subunit RNA polymerase (RNAP) from bacteriophage T7, implementing all-atom molecular dynamics (MD) simulations. We obtained a jump-from-cavity kinetic model of the PPi release utilizing extensive nanosecond MD simulations. We found that the PPi release in T7 RNAP is initiated by the PPi dissociation from two catalytic aspartic acids, followed by a comparatively slow jump-from-cavity activation process. Combining with a number of microsecond long MD simulations, we also found that the activation process is hindered by charged residue associations as well as by local steric and hydrogen bond interactions. On the other hand, the activation is greatly assisted by a highly flexible lysine residue Lys472 that swings its side chain to pull PPi out. The mechanism can apply in general to single subunit RNA and DNA polymerases with similar molecular structures and conserved key residues. Remarkably, the flexible lysine or arginine residue appears to be a universal module that assists the PPi release even in multi-subunit RNAPs with charge facilitated hopping mechanisms. We also noticed that the PPi release is not tightly coupled to opening motions of an O-helix on the fingers domain of T7 RNAP according to the microsecond MD simulations. Our study thus supports the Brownian ratchet scenario of the mechano-chemical coupling in the transcription elongation of the single-subunit polymerase.
Collapse
Affiliation(s)
- Lin-Tai Da
- Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida, United States of America
| | - Chao E
- Beijing Computational Science Research Center, Beijing, China
| | - Baogen Duan
- Beijing Computational Science Research Center, Beijing, China
| | - Chuanbiao Zhang
- School of Physics, University of the Chinese Academy of Sciences, Beijing, China
| | - Xin Zhou
- School of Physics, University of the Chinese Academy of Sciences, Beijing, China
| | - Jin Yu
- Beijing Computational Science Research Center, Beijing, China
- * E-mail:
| |
Collapse
|
32
|
Duan B, Wu S, Da LT, Yu J. A critical residue selectively recruits nucleotides for t7 RNA polymerase transcription fidelity control. Biophys J 2015; 107:2130-40. [PMID: 25418098 PMCID: PMC4223216 DOI: 10.1016/j.bpj.2014.09.038] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2014] [Revised: 09/23/2014] [Accepted: 09/30/2014] [Indexed: 01/21/2023] Open
Abstract
Nucleotide selection is essential for fidelity control in gene replication and transcription. Recent work on T7 RNA polymerase suggested that a small posttranslocation free energy bias stabilizes Tyr(639) in the active site to aid nucleotide selection. However, it was not clear exactly how Tyr(639) assists the selection. Here we report a molecular-dynamics simulation study revealing atomistic detail of this critical selectivity. The study shows first that Tyr(639) blocks the active site at posttranslocation by marginally stacking to the end basepair of the DNA-RNA hybrid. The study then demonstrates that at the nucleotide preinsertion state, a cognate RNA nucleotide does not affect the local Tyr(639) stabilization, whereas a noncognate nucleotide substantially stabilizes Tyr(639) so that Tyr(639) keeps blocking the active site. As a result, further nucleotide insertion into the active site, which requires moving Tyr(639) out of the site, would be hindered for the noncognate nucleotide, but not for the cognate nucleotide. In particular, we note that water molecules assist the ribose recognition in the RNA nucleotide preinsertion, and help Tyr(639) stacking to the end basepair in the case of a DNA nucleotide. It was also seen that a base-mismatched nucleotide at preinsertion directly grabs Tyr(639) for the active site stabilization. We also find that in a mutant polymerase Y639F the strong stabilization of residue 639 in the active site cannot establish upon the DNA nucleotide preinsertion. The finding explains the reduced differentiation between ribo- and deoxyribonucleotides that has been recorded experimentally for the mutant polymerase.
Collapse
Affiliation(s)
- Baogen Duan
- Beijing Computational Science Research Center, Beijing, P. R. China
| | - Shaogui Wu
- Beijing Computational Science Research Center, Beijing, P. R. China
| | - Lin-Tai Da
- Department of Physics and Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida
| | - Jin Yu
- Beijing Computational Science Research Center, Beijing, P. R. China.
| |
Collapse
|
33
|
Zhang L, Silva DA, Pardo-Avila F, Wang D, Huang X. Structural Model of RNA Polymerase II Elongation Complex with Complete Transcription Bubble Reveals NTP Entry Routes. PLoS Comput Biol 2015; 11:e1004354. [PMID: 26134169 PMCID: PMC4489626 DOI: 10.1371/journal.pcbi.1004354] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2015] [Accepted: 05/22/2015] [Indexed: 12/27/2022] Open
Abstract
The RNA polymerase II (Pol II) is a eukaryotic enzyme that catalyzes the synthesis of the messenger RNA using a DNA template. Despite numerous biochemical and biophysical studies, it remains elusive whether the “secondary channel” is the only route for NTP to reach the active site of the enzyme or if the “main channel” could be an alternative. On this regard, crystallographic structures of Pol II have been extremely useful to understand the structural basis of transcription, however, the conformation of the unpaired non-template DNA part of the full transcription bubble (TB) is still unknown. Since diffusion routes of the nucleoside triphosphate (NTP) substrate through the main channel might overlap with the TB region, gaining structural information of the full TB is critical for a complete understanding of Pol II transcription process. In this study, we have built a structural model of Pol II with a complete transcription bubble based on multiple sources of existing structural data and used Molecular Dynamics (MD) simulations together with structural analysis to shed light on NTP entry pathways. Interestingly, we found that although both channels have enough space to allow NTP loading, the percentage of MD conformations containing enough space for NTP loading through the secondary channel is twice higher than that of the main channel. Further energetic study based on MD simulations with NTP loaded in the channels has revealed that the diffusion of the NTP through the main channel is greatly disfavored by electrostatic repulsion between the NTP and the highly negatively charged backbones of nucleotides in the non-template DNA strand. Taken together, our results suggest that the secondary channel is the major route for NTP entry during Pol II transcription. In eukaryotic cells, the RNA polymerase II (Pol II) is a central enzyme that reads the genetic information encoded in the DNA template to synthetize a messenger RNA. To perform its function, Pol II needs to have the substrate nucleoside triphosphate (NTP) diffuse into its deeply buried active site. Despite numerous efforts, the NTP entry routes remain elusive: NTP could diffuse only through the secondary channel, or also via the main channel. The structural information of the transcription bubble is essential to study this process, however, the unpaired non-template DNA of the transcription bubble is absent in the available X-ray crystal structures. In this regard, we have built a structural model of the Pol II elongation complex with reconstructed transcription bubble using existing experimental data. We then performed Molecular Dynamics (MD) simulations and applied structural analysis to study the routes of NTP diffusion. We found that sterically the probability of NTP loading through the secondary channel is more than twice that of the main channel. Further analysis of the non-bonded energetic contributions to NTP diffusion suggests that NTP diffusion through the main channel is greatly disfavored by the electrostatic repulsion between the substrate and negatively charged backbones of nucleotides in the non-template strand of the transcription bubble. Altogether, our findings suggest that the secondary channel is the more favorable NTP diffusion route for Pol II transcription elongation.
Collapse
Affiliation(s)
- Lu Zhang
- Department of Chemistry and State Key Laboratory of Molecular Neuroscience, Center for System Biology and Human Health, School of Science and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Daniel-Adriano Silva
- Department of Chemistry and State Key Laboratory of Molecular Neuroscience, Center for System Biology and Human Health, School of Science and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Fátima Pardo-Avila
- Department of Chemistry and State Key Laboratory of Molecular Neuroscience, Center for System Biology and Human Health, School of Science and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
| | - Dong Wang
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, California, United States of America
| | - Xuhui Huang
- Department of Chemistry and State Key Laboratory of Molecular Neuroscience, Center for System Biology and Human Health, School of Science and IAS, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- * E-mail:
| |
Collapse
|
34
|
Free energy landscape of activation in a signalling protein at atomic resolution. Nat Commun 2015; 6:7284. [PMID: 26073309 PMCID: PMC4470301 DOI: 10.1038/ncomms8284] [Citation(s) in RCA: 69] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2014] [Accepted: 04/26/2015] [Indexed: 11/24/2022] Open
Abstract
The interconversion between inactive and active protein states, traditionally described by two static structures, is at the heart of signaling. However, how folded states interconvert is largely unknown due to the inability to experimentally observe transition pathways. Here we explore the free energy landscape of the bacterial response regulator NtrC by combining computation and NMR, and discover unexpected features underlying efficient signaling. We find that functional states are defined purely in kinetic and not structural terms. The need of a well-defined conformer, crucial to the active state, is absent in the inactive state, which comprises a heterogeneous collection of conformers. The transition between active and inactive states occurs through multiple pathways, facilitated by a number of nonnative transient hydrogen bonds, thus lowering the transition barrier through both entropic and enthalpic contributions. These findings may represent general features for functional conformational transitions within the folded state.
Collapse
|
35
|
Bai Q, Pérez-Sánchez H, Zhang Y, Shao Y, Shi D, Liu H, Yao X. Ligand induced change of β2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis. Phys Chem Chem Phys 2015; 16:15874-85. [PMID: 24962153 DOI: 10.1039/c4cp01185f] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The reported crystal structures of β2 adrenergic receptor (β2AR) reveal that the open and closed states of the water channel are correlated with the inactive and active conformations of β2AR. However, more details about the process by which the water channel states are affected by the active to inactive conformational change of β2AR remain illusive. In this work, molecular dynamics simulations are performed to study the dynamical inactive and active conformational change of β2AR induced by inverse agonist ICI 118,551. Markov state model analysis and free energy calculation are employed to explore the open and close states of the water channel. The simulation results show that inverse agonist ICI 118,551 can induce water channel opening during the conformational transition of β2AR. Markov state model (MSM) analysis proves that the energy contour can be divided into seven states. States S1, S2 and S5, which represent the active conformation of β2AR, show that the water channel is in the closed state, while states S4 and S6, which correspond to the intermediate state conformation of β2AR, indicate the water channel opens gradually. State S7, which represents the inactive structure of β2AR, corresponds to the full open state of the water channel. The opening mechanism of the water channel is involved in the ligand-induced conformational change of β2AR. These results can provide useful information for understanding the opening mechanism of the water channel and will be useful for the rational design of potent inverse agonists of β2AR.
Collapse
Affiliation(s)
- Qifeng Bai
- Department of Chemistry, Lanzhou University, Lanzhou 730000, China.
| | | | | | | | | | | | | |
Collapse
|
36
|
Wang B, Opron K, Burton ZF, Cukier RI, Feig M. Five checkpoints maintaining the fidelity of transcription by RNA polymerases in structural and energetic details. Nucleic Acids Res 2014; 43:1133-46. [PMID: 25550432 PMCID: PMC4333413 DOI: 10.1093/nar/gku1370] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Transcriptional fidelity, which prevents the misincorporation of incorrect nucleoside monophosphates in RNA, is essential for life. Results from molecular dynamics (MD) simulations of eukaryotic RNA polymerase (RNAP) II and bacterial RNAP with experimental data suggest that fidelity may involve as many as five checkpoints. Using MD simulations, the effects of different active site NTPs in both open and closed trigger loop (TL) structures of RNAPs are compared. Unfavorable initial binding of mismatched substrates in the active site with an open TL is proposed to be the first fidelity checkpoint. The leaving of an incorrect substrate is much easier than a correct one energetically from the umbrella sampling simulations. Then, the closing motion of the TL, required for catalysis, is hindered by the presence of mismatched NTPs. Mismatched NTPs also lead to conformational changes in the active site, which perturb the coordination of magnesium ions and likely affect the ability to proceed with catalysis. This step appears to be the most important checkpoint for deoxy-NTP discrimination. Finally, structural perturbations in the template DNA and the nascent RNA in the presence of mismatches likely hinder nucleotide addition and provide the structural foundation for backtracking followed by removing erroneously incorporated nucleotides during proofreading.
Collapse
Affiliation(s)
- Beibei Wang
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Kristopher Opron
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Zachary F Burton
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
| | - Robert I Cukier
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA
| | - Michael Feig
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI 48824, USA Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
37
|
Yu J, Da LT, Huang X. Constructing kinetic models to elucidate structural dynamics of a complete RNA polymerase II elongation cycle. Phys Biol 2014; 12:016004. [DOI: 10.1088/1478-3975/12/1/016004] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
38
|
Gu S, Silva DA, Meng L, Yue A, Huang X. Quantitatively characterizing the ligand binding mechanisms of choline binding protein using Markov state model analysis. PLoS Comput Biol 2014; 10:e1003767. [PMID: 25101697 PMCID: PMC4125059 DOI: 10.1371/journal.pcbi.1003767] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2013] [Accepted: 06/22/2014] [Indexed: 01/05/2023] Open
Abstract
Protein-ligand recognition plays key roles in many biological processes. One of the most fascinating questions about protein-ligand recognition is to understand its underlying mechanism, which often results from a combination of induced fit and conformational selection. In this study, we have developed a three-pronged approach of Markov State Models, Molecular Dynamics simulations, and flux analysis to determine the contribution of each model. Using this approach, we have quantified the recognition mechanism of the choline binding protein (ChoX) to be ∼90% conformational selection dominant under experimental conditions. This is achieved by recovering all the necessary parameters for the flux analysis in combination with available experimental data. Our results also suggest that ChoX has several metastable conformational states, of which an apo-closed state is dominant, consistent with previous experimental findings. Our methodology holds great potential to be widely applied to understand recognition mechanisms underlining many fundamental biological processes.
Collapse
Affiliation(s)
- Shuo Gu
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Daniel-Adriano Silva
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
| | - Luming Meng
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Alexander Yue
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Division of Biomedical Engineering, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- * E-mail:
| |
Collapse
|
39
|
Chodera JD, Noé F. Markov state models of biomolecular conformational dynamics. Curr Opin Struct Biol 2014; 25:135-44. [PMID: 24836551 DOI: 10.1016/j.sbi.2014.04.002] [Citation(s) in RCA: 502] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 04/08/2014] [Accepted: 04/12/2014] [Indexed: 10/25/2022]
Abstract
It has recently become practical to construct Markov state models (MSMs) that reproduce the long-time statistical conformational dynamics of biomolecules using data from molecular dynamics simulations. MSMs can predict both stationary and kinetic quantities on long timescales (e.g. milliseconds) using a set of atomistic molecular dynamics simulations that are individually much shorter, thus addressing the well-known sampling problem in molecular dynamics simulation. In addition to providing predictive quantitative models, MSMs greatly facilitate both the extraction of insight into biomolecular mechanism (such as folding and functional dynamics) and quantitative comparison with single-molecule and ensemble kinetics experiments. A variety of methodological advances and software packages now bring the construction of these models closer to routine practice. Here, we review recent progress in this field, considering theoretical and methodological advances, new software tools, and recent applications of these approaches in several domains of biochemistry and biophysics, commenting on remaining challenges.
Collapse
Affiliation(s)
- John D Chodera
- Computational Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA.
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany.
| |
Collapse
|
40
|
Bowman GR, Meng L, Huang X. Quantitative comparison of alternative methods for coarse-graining biological networks. J Chem Phys 2014; 139:121905. [PMID: 24089717 DOI: 10.1063/1.4812768] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Markov models and master equations are a powerful means of modeling dynamic processes like protein conformational changes. However, these models are often difficult to understand because of the enormous number of components and connections between them. Therefore, a variety of methods have been developed to facilitate understanding by coarse-graining these complex models. Here, we employ Bayesian model comparison to determine which of these coarse-graining methods provides the models that are most faithful to the original set of states. We find that the Bayesian agglomerative clustering engine and the hierarchical Nyström expansion graph (HNEG) typically provide the best performance. Surprisingly, the original Perron cluster cluster analysis (PCCA) method often provides the next best results, outperforming the newer PCCA+ method and the most probable paths algorithm. We also show that the differences between the models are qualitatively significant, rather than being minor shifts in the boundaries between states. The performance of the methods correlates well with the entropy of the resulting coarse-grainings, suggesting that finding states with more similar populations (i.e., avoiding low population states that may just be noise) gives better results.
Collapse
Affiliation(s)
- Gregory R Bowman
- Departments of Chemistry and Molecular and Cell Biology, University of California, Berkeley, California 94720, USA
| | | | | |
Collapse
|
41
|
Application of Markov State Models to Simulate Long Timescale Dynamics of Biological Macromolecules. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 805:29-66. [DOI: 10.1007/978-3-319-02970-2_2] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
42
|
PARDO-AVILA FÁTIMA, DA LINTAI, WANG YING, HUANG XUHUI. THEORETICAL INVESTIGATIONS ON ELUCIDATING FUNDAMENTAL MECHANISMS OF CATALYSIS AND DYNAMICS INVOLVED IN TRANSCRIPTION BY RNA POLYMERASE. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2013. [DOI: 10.1142/s0219633613410058] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
RNA polymerase is the enzyme that synthesizes RNA during the transcription process. To understand its mechanism, structural studies have provided us pictures of the series of steps necessary to add a new nucleotide to the nascent RNA chain, the steps altogether known as the nucleotide addition cycle (NAC). However, these static snapshots do not provide dynamic information of these processes involved in NAC, such as the conformational changes of the protein and the atomistic details of the catalysis. Computational studies have made efforts to fill these knowledge gaps. In this review, we provide examples of different computational approaches that have improved our understanding of the transcription elongation process for RNA polymerase, such as normal mode analysis, molecular dynamic (MD) simulations, Markov state models (MSMs). We also point out some unsolved questions that could be addressed using computational tools in the future.
Collapse
Affiliation(s)
- FÁTIMA PARDO-AVILA
- Department of Chemistry, Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, Hong Kong University of Science and Technology, Clear Water Bay Road, Kowloon, Hong Kong
| | - LIN-TAI DA
- Department of Chemistry, Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, Hong Kong University of Science and Technology, Clear Water Bay Road, Kowloon, Hong Kong
| | - YING WANG
- Department of Chemistry, Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, Hong Kong University of Science and Technology, Clear Water Bay Road, Kowloon, Hong Kong
| | - XUHUI HUANG
- Department of Chemistry, Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, Hong Kong University of Science and Technology, Clear Water Bay Road, Kowloon, Hong Kong
- Division of Biomedical Engineering, Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, Hong Kong University of Science and Technology, Clear Water Bay Road, Kowloon, Hong Kong
| |
Collapse
|
43
|
Tan C, Li W, Wang W. Localized frustration and binding-induced conformational change in recognition of 5S RNA by TFIIIA zinc finger. J Phys Chem B 2013; 117:15917-25. [PMID: 24266699 DOI: 10.1021/jp4052165] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein TFIIIA is composed of nine tandemly arranged Cys2His2 zinc fingers. It can bind either to the 5S RNA gene as a transcription factor or to the 5S RNA transcript as a chaperone. Although structural and biochemical data provided valuable information on the recognition between the TFIIIIA and the 5S DNA/RNA, the involved conformational motions and energetic factors contributing to the binding affinity and specificity remain unclear. In this work, we conducted MD simulations and MM/GBSA calculations to investigate the binding-induced conformational changes in the recognition of the 5S RNA by the central three zinc fingers of TFIIIA and the energetic factors that influence the binding affinity and specificity at an atomistic level. Our results revealed drastic interdomain conformational changes between these three zinc fingers, involving the exposure/burial of several crucial DNA/RNA binding residues, which can be related to the competition between DNA and RNA for the binding of TFIIIA. We also showed that the specific recognition between finger 4/finger 6 and the 5S RNA introduces frustrations to the nonspecific interactions between finger 5 and the 5S RNA, which may be important to achieve optimal binding affinity and specificity.
Collapse
Affiliation(s)
- Cheng Tan
- National Laboratory of Solid State Microstructure and Department of Physics, Nanjing University , Nanjing, Jiangsu 210093, China
| | | | | |
Collapse
|