1
|
Vergara S, Zhou X, Santiago U, Conway JF, Sluis-Cremer N, Calero G. Structures of kinetic intermediate states of HIV-1 reverse transcriptase DNA synthesis. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.18.572243. [PMID: 38187617 PMCID: PMC10769260 DOI: 10.1101/2023.12.18.572243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/09/2024]
Abstract
Reverse transcription of the retroviral single-stranded RNA into double-stranded DNA is an integral step during HIV-1 replication, and reverse transcriptase (RT) is a primary target for antiviral therapy. Despite a wealth of structural information on RT, we lack critical insight into the intermediate kinetic states of DNA synthesis. Using catalytically active substrates, and a novel blot/diffusion cryo-electron microscopy approach, we captured 11 structures that define the substrate binding, reactant, transition and product states of dATP addition by RT at 1.9 to 2.4 Å resolution in the active site. Initial dATP binding to RT-template/primer complex involves a single Mg 2+ (site B), and promotes partial closure of the active site pocket by a large conformational change in the β3-β4 loop in the Fingers domain, and formation of a negatively charged pocket where a second "drifting" Mg 2+ can bind (site A). During the transition state, the α-phosphate oxygen from a previously unobserved dATP conformer aligns with the site A Mg 2+ and the primer 3'-OH for nucleophilic attack. In the product state, we captured two substrate conformations in the active site: 1) dATP that had yet to be incorporated into the nascent DNA, and 2) an incorporated dAMP with the pyrophosphate leaving group coordinated by metal B and stabilized through H- bonds in the active site of RT. This study provides insights into a fundamental chemical reaction that impacts polymerase fidelity, nucleoside inhibitor drug design, and mechanisms of drug resistance.
Collapse
|
2
|
Dangerfield TL, Johnson KA. Design and interpretation of experiments to establish enzyme pathway and define the role of conformational changes in enzyme specificity. Methods Enzymol 2023; 685:461-492. [PMID: 37245912 DOI: 10.1016/bs.mie.2023.03.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We describe the experimental methods and analysis to define the role of enzyme conformational changes in specificity based on published studies using DNA polymerases as an ideal model system. Rather than give details of how to perform transient-state and single-turnover kinetic experiments, we focus on the rationale of the experimental design and interpretation. We show how initial experiments to measure kcat and kcat/Km can accurately quantify specificity but do not define its underlying mechanistic basis. We describe methods to fluorescently label enzymes to monitor conformational changes and to correlate fluorescence signals with rapid-chemical-quench flow assays to define the steps in the pathway. Measurements of the rate of product release and of the kinetics of the reverse reaction complete the kinetic and thermodynamic description of the full reaction pathway. This analysis showed that the substrate-induced change in enzyme structure from an open to a closed state was much faster than rate-limiting chemical bond formation. However, because the reverse of the conformational change was much slower than chemistry, specificity is governed solely by the product of the binding constant for the initial weak substrate binding and the rate constant for the conformational change (kcat/Km=K1k2) so that the specificity constant does not include kcat. The enzyme conformational change leads to a closed complex in which the substrate is bound tightly and is committed to the forward reaction. In contrast, an incorrect substrate is bound weakly, and the rate of chemistry is slow, so the mismatch is released from the enzyme rapidly. Thus, the substrate-induced-fit is the major determinant of specificity. The methods outlined here should be applicable to other enzyme systems.
Collapse
Affiliation(s)
- Tyler L Dangerfield
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX, United States
| | - Kenneth A Johnson
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX, United States.
| |
Collapse
|
3
|
Cardenas AE, Hunter A, Wang H, Elber R. ScMiles2: A Script to Conduct and Analyze Milestoning Trajectories for Long Time Dynamics. J Chem Theory Comput 2022; 18:6952-6965. [PMID: 36191005 PMCID: PMC10336853 DOI: 10.1021/acs.jctc.2c00708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Milestoning is a theory and an algorithm that computes kinetics and thermodynamics at long time scales. It is based on partitioning the (phase) space into cells and running a large number of short trajectories between the boundaries of the cells. The termination points of the trajectories are analyzed with the Milestoning theory to obtain kinetic and thermodynamic information. Managing the tens to hundreds of thousands of Milestoning trajectories is a challenge, which we handle with a python script, ScMiles. Here, we introduce a new version of the python script ScMiles2 to conduct Milestoning simulations. Major enhancements are: (i) post analysis of Milestoning trajectories to obtain the free energy, mean first passage time, the committor function, and exit times; (ii) similar to (i) but the post analysis is for a single long trajectory; (iii) we support the use of the GROMACS software in addition to NAMD; (iv) a restart option; (v) the automated finding, sampling, and launching trajectories from new milestones that are found on the fly; and (vi) support Milestoning calculations with several coarse variables and for complex reaction coordinates. We also evaluate the simulation parameters and suggest new algorithmic features to enhance the rate of convergence of observables. We propose the use of an iteration-averaged kinetic matrix for a rapid approach to asymptotic values. Illustrations are provided for small systems and one large example.
Collapse
Affiliation(s)
- Alfredo E. Cardenas
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
| | - Allison Hunter
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
| | - Hao Wang
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, China
| | - Ron Elber
- The Oden Institute, University of Texas at Austin, Austin, Texas, 78712, USA
- Department of Chemistry, University of Texas at Austin, Austin, Texas, 78712, USA
| |
Collapse
|
4
|
Kumari A, Yadav A, Lahiri I. Transient State Kinetics of Plasmodium falciparum Apicoplast DNA Polymerase Suggests the Involvement of Accessory Factors for Efficient and Accurate DNA Synthesis. Biochemistry 2022; 61:2319-2333. [PMID: 36251801 DOI: 10.1021/acs.biochem.2c00446] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Plasmodium, the causative agent of malaria, belongs to the phylum Apicomplexa. Most apicomplexans, including Plasmodium, contain an essential nonphotosynthetic plastid called the apicoplast that harbors its own genome that is replicated by a dedicated organellar replisome. This replisome employs a single DNA polymerase (apPol), which is expected to perform both replicative and translesion synthesis. Unlike other replicative polymerases, no processivity factor for apPol has been identified. While preliminary structural and biochemical studies have provided an overall characterization of apPol, the kinetic mechanism of apPol's activity remains unknown. We have used transient state methods to determine the kinetics of replicative and translesion synthesis by apPol and show that apPol has low processivity and efficiency while copying undamaged DNA. Moreover, while apPol can bypass oxidatively damaged lesions, the bypass is error-prone. Taken together, our results raise the following question─how does a polymerase with low processivity, efficiency, and fidelity (for translesion synthesis) faithfully replicate the apicoplast organellar DNA within the hostile environment of the human host? We hypothesize that interactions with putative components of the apicoplast replisome and/or an as-yet-undiscovered processivity factor transform apPol into an efficient and accurate enzyme.
Collapse
Affiliation(s)
- Anamika Kumari
- Department of Biological Sciences, Indian Institute of Science Education and Research Mohali, Punjab 140306, India
| | - Anjali Yadav
- Department of Biological Sciences, Indian Institute of Science Education and Research Mohali, Punjab 140306, India
| | - Indrajit Lahiri
- Department of Biological Sciences, Indian Institute of Science Education and Research Mohali, Punjab 140306, India.,Molecular Microbiology, School of Biosciences, University of Sheffield, Sheffield S10 2TN, U.K
| |
Collapse
|
5
|
Jones PE, Pérez-Segura C, Bryer AJ, Perilla JR, Hadden-Perilla JA. Molecular dynamics of the viral life cycle: progress and prospects. Curr Opin Virol 2021; 50:128-138. [PMID: 34464843 PMCID: PMC8651149 DOI: 10.1016/j.coviro.2021.08.003] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Revised: 08/09/2021] [Accepted: 08/09/2021] [Indexed: 01/29/2023]
Abstract
Molecular dynamics (MD) simulations across spatiotemporal resolutions are widely applied to study viruses and represent the central technique uniting the field of computational virology. We discuss the progress of MD in elucidating the dynamics of the viral life cycle, including the status of modeling intact extracellular virions and leveraging advanced simulations to mimic active life cycle processes. We further remark on the prospects of MD for continued contributions to the basic science characterization of viruses, especially given the increasing availability of high-quality experimental data and supercomputing power. Overall, integrative computational methods that are closely guided by experiments are unmatched in the level of detail they provide, enabling-now and in the future-new discoveries relevant to thwarting viral infection.
Collapse
Affiliation(s)
- Peter Eugene Jones
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Carolina Pérez-Segura
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Alexander J Bryer
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Juan R Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Jodi A Hadden-Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States.
| |
Collapse
|
6
|
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
|
7
|
Abstract
The protein HIV Reverse Transcriptase (HIV RT) synthesizes a DNA strand according to a template. During the synthesis, the polymerase slides on the double stranded DNA to allow the entry of a new nucleotide to the active site. We use Molecular Dynamics simulations to estimate the free energy profile and the time scale of the DNA-protein's relative displacement in the complex's closed state. We illustrate that the presence of the catalytic magnesium slows down the process. Upon removing the catalytic magnesium ion, the process is rapid and significantly faster than reopening the active site in preparation for the new substrate. We speculate that magnesium regulates DNA translocation. The magnesium locks the DNA into a specific orientation during the chemical addition of the nucleotide. The release of Mg2+ eases DNA sliding and the acceptance of a new substrate.
Collapse
Affiliation(s)
- Hao Wang
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, United States.,Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
8
|
Alberini G, Benfenati F, Maragliano L. Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations. J Chem Inf Model 2021; 61:1354-1367. [PMID: 33570938 PMCID: PMC8023575 DOI: 10.1021/acs.jcim.0c01407] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Activation of voltage-gated
ion channels is regulated by conformational
changes of the voltage sensor domains (VSDs), four water- and ion-impermeable
modules peripheral to the central, permeable pore domain. Anomalous
currents, defined as ω-currents, have been recorded in response
to mutations of residues on the VSD S4 helix and associated with ion
fluxes through the VSDs. In humans, gene defects in the potassium
channel Kv7.2 result in a broad range of epileptic disorders, from
benign neonatal seizures to severe epileptic encephalopathies. Experimental
evidence suggests that the R207Q mutation in S4, associated with peripheral
nerve hyperexcitability, induces ω-currents at depolarized potentials,
but the fine structural details are still elusive. In this work, we
use atom-detailed molecular dynamics simulations and a refined model
structure of the Kv7.2 VSD in the active conformation in a membrane/water
environment to study the effect of R207Q and four additional mutations
of proven clinical importance. Our results demonstrate that the R207Q
mutant shows the most pronounced increase of hydration in the internal
VSD cavity, a feature favoring the occurrence of ω-currents.
Free energy and kinetics calculations of sodium permeation through
the native and mutated VSD indicate as more favorable the formation
of a cationic current in the latter. Overall, our simulations establish
a mechanistic linkage between genetic variations and their physiological
outcome, by providing a computational description that includes both
thermodynamic and kinetic features of ion permeation associated with
ω-currents.
Collapse
Affiliation(s)
- Giulio Alberini
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,Department of Experimental Medicine, Università degli Studi di Genova, Viale Benedetto XV, 3, 16132 Genova, Italy
| | - Fabio Benfenati
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi, 10, 16132 Genova, Italy
| | - Luca Maragliano
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| |
Collapse
|
9
|
Elber R, Fathizadeh A, Ma P, Wang H. Modeling molecular kinetics with Milestoning. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1512] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ron Elber
- Department of Chemistry, The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Arman Fathizadeh
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Piao Ma
- Department of Chemistry University of Texas at Austin Austin Texas USA
| | - Hao Wang
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| |
Collapse
|
10
|
Gong S, Kirmizialtin S, Chang A, Mayfield JE, Zhang YJ, Johnson KA. Kinetic and thermodynamic analysis defines roles for two metal ions in DNA polymerase specificity and catalysis. J Biol Chem 2020; 296:100184. [PMID: 33310704 PMCID: PMC7948414 DOI: 10.1074/jbc.ra120.016489] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Revised: 12/05/2020] [Accepted: 12/11/2020] [Indexed: 11/06/2022] Open
Abstract
Magnesium ions play a critical role in catalysis by many enzymes and contribute to the fidelity of DNA polymerases through a two-metal ion mechanism. However, specificity is a kinetic phenomenon and the roles of Mg2+ ions in each step in the catalysis have not been resolved. We first examined the roles of Mg2+ by kinetic analysis of single nucleotide incorporation catalyzed by HIV reverse transcriptase. We show that Mg.dNTP binding induces an enzyme conformational change at a rate that is independent of free Mg2+ concentration. Subsequently, the second Mg2+ binds to the closed state of the enzyme-DNA-Mg.dNTP complex (Kd = 3.7 mM) to facilitate catalysis. Weak binding of the catalytic Mg2+ contributes to fidelity by sampling the correctly aligned substrate without perturbing the equilibrium for nucleotide binding at physiological Mg2+ concentrations. An increase of the Mg2+ concentration from 0.25 to 10 mM increases nucleotide specificity (kcat/Km) 12-fold largely by increasing the rate of the chemistry relative to the rate of nucleotide release. Mg2+ binds very weakly (Kd ≤ 37 mM) to the open state of the enzyme. Analysis of published crystal structures showed that HIV reverse transcriptase binds only two metal ions prior to incorporation of a correct base pair. Molecular dynamics simulations support the two-metal ion mechanism and the kinetic data indicating weak binding of the catalytic Mg2+. Molecular dynamics simulations also revealed the importance of the divalent cation cloud surrounding exposed phosphates on the DNA. These results enlighten the roles of the two metal ions in the specificity of DNA polymerases.
Collapse
Affiliation(s)
- Shanzhong Gong
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas, USA
| | - Serdal Kirmizialtin
- Chemistry Program, Science Division, New York University Abu Dhabi, Abu Dhabi, United Arab Emirates
| | - Adrienne Chang
- Chemistry Program, Science Division, New York University Abu Dhabi, Abu Dhabi, United Arab Emirates
| | - Joshua E Mayfield
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas, USA
| | - Yan Jessie Zhang
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas, USA
| | - Kenneth A Johnson
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, Texas, USA.
| |
Collapse
|
11
|
Dangerfield TL, Johnson KA. Conformational dynamics during high-fidelity DNA replication and translocation defined using a DNA polymerase with a fluorescent artificial amino acid. J Biol Chem 2020; 296:100143. [PMID: 33273013 PMCID: PMC7857513 DOI: 10.1074/jbc.ra120.016617] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 12/02/2020] [Accepted: 12/03/2020] [Indexed: 12/21/2022] Open
Abstract
We address the role of enzyme conformational dynamics in specificity for a high-fidelity DNA polymerase responsible for genome replication. We present the complete characterization of the conformational dynamics during the correct nucleotide incorporation forward and reverse reactions using stopped-flow and rapid-quench methods with a T7 DNA polymerase variant containing a fluorescent unnatural amino acid, (7-hydroxy-4-coumarin-yl) ethylglycine, which provides a signal for enzyme conformational changes. We show that the forward conformational change (>6000 s−1) is much faster than chemistry (300 s−1) while the enzyme opening to allow release of bound nucleotide (1.7 s−1) is much slower than chemistry. These parameters show that the conformational change selects a correct nucleotide for incorporation through an induced-fit mechanism. We also measured conformational changes occurring after chemistry and during pyrophosphorolysis, providing new insights into processive polymerization. Pyrophosphorolysis occurs via a conformational selection mechanism as the pyrophosphate binds to a rare pretranslocation state of the enzyme–DNA complex. Global data fitting was achieved by including experiments in the forward and reverse directions to correlate conformational changes with chemical reaction steps. This analysis provided well-constrained values for nine rate constants to establish a complete free-energy profile including the rates of DNA translocation during processive synthesis. Translocation does not follow Brownian ratchet or power stroke models invoking nucleotide binding as the driving force. Rather, translocation is rapid and thermodynamically favorable after enzyme opening and pyrophosphate release, and it appears to limit the rate of processive synthesis at 4 °C.
Collapse
Affiliation(s)
- Tyler L Dangerfield
- Institute for Cellular and Molecular Biology, Department of Molecular Biosciences, University of Texas, Austin, Texas, USA
| | - Kenneth A Johnson
- Institute for Cellular and Molecular Biology, Department of Molecular Biosciences, University of Texas, Austin, Texas, USA.
| |
Collapse
|
12
|
Wang H, Huang N, Dangerfield T, Johnson KA, Gao J, Elber R. Exploring the Reaction Mechanism of HIV Reverse Transcriptase with a Nucleotide Substrate. J Phys Chem B 2020; 124:4270-4283. [PMID: 32364738 PMCID: PMC7260111 DOI: 10.1021/acs.jpcb.0c02632] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Enzymatic reactions consist of several steps: (i) a weak binding event of the substrate to the enzyme, (ii) an induced fit or a protein conformational transition upon ligand binding, (iii) the chemical reaction, and (iv) the release of the product. Here we focus on step iii of the reaction of a DNA polymerase, HIV RT, with a nucleotide. We determine the rate and the free energy profile for the addition of a nucleotide to a DNA strand using a combination of a QM/MM model, the string method, and exact Milestoning. The barrier height and the time scale of the reaction are consistent with experiment. We show that the observables (free energies and mean first passage time) converge rapidly, as a function of the Milestoning iteration number. We also consider the substitution of an oxygen of the incoming nucleotide by a nonbridging sulfur atom and its impact on the enzymatic reaction. This substitution has been suggested in the past as a tool to examine the influence of the chemical step on the overall rate. Our joint computational and experimental study suggests that the impact of the substitution is small. Computationally, the differences between the two are within the estimated error bars. Experiments suggest a small difference. Finally, we examine step i, the weak binding of the nucleotide to the protein surface. We suggest that this step has only a small contribution to the selectivity of the enzyme. Comments are made on the impact of these steps on the overall mechanism.
Collapse
Affiliation(s)
- Hao Wang
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin TX 78712
| | - Nathan Huang
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712
| | - Tyler Dangerfield
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712
| | - Kenneth A. Johnson
- Department of Molecular Biosciences, University of Texas at Austin, Austin, TX 78712
| | - Jiali Gao
- Department of Chemistry, University of Minnesota, Minneapolis, MN, 55455-0431
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin TX 78712
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712
| |
Collapse
|
13
|
Narayan B, Fathizadeh A, Templeton C, He P, Arasteh S, Elber R, Buchete NV, Levy RM. The transition between active and inactive conformations of Abl kinase studied by rock climbing and Milestoning. Biochim Biophys Acta Gen Subj 2020; 1864:129508. [PMID: 31884066 PMCID: PMC7012767 DOI: 10.1016/j.bbagen.2019.129508] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2019] [Revised: 12/11/2019] [Accepted: 12/19/2019] [Indexed: 12/12/2022]
Abstract
BACKGROUND Kinases are a family of enzymes that catalyze the transfer of the ɤ-phosphate group from ATP to a protein's residue. Malfunctioning kinases are involved in many health problems such as cardiovascular diseases, diabetes, and cancer. Kinases transitions between multiple conformations of inactive to active forms attracted considerable interest. METHOD A reaction coordinate is computed for the transition between the active to inactive conformation in Abl kinase with a focus on the DFG-in to DFG-out flip. The method of Rock Climbing is used to construct a path locally, which is subsequently optimized using a functional of the entire path. The discrete coordinate sets along the reaction path are used in a Milestoning calculation of the free energy landscape and the rate of the transition. RESULTS The estimated transition times are between a few milliseconds and seconds, consistent with simulations of the kinetics and with indirect experimental data. The activation requires the transient dissociation of the salt bridge between Lys271 and Glu286. The salt bridge reforms once the DFG motif is stabilized by a locked conformation of Phe382. About ten residues are identified that contribute significantly to the process and are included as part of the reaction space. CONCLUSIONS The transition from DFG-in to DFG-out in Abl kinase was simulated using atomic resolution of a fully solvated protein yielding detailed description of the kinetics and the mechanism of the DFG flip. The results are consistent with other computational methods that simulate the kinetics and with some indirect experimental measurements. GENERAL SIGNIFICANCE The activation of kinases includes a conformational transition of the DFG motif that is important for enzyme activity but is not accessible to conventional Molecular Dynamics. We propose a detailed mechanism for the transition, at a timescale longer than conventional MD, using a combination of reaction path and Milestoning algorithms. The mechanism includes local structural adjustments near the binding site as well as collective interactions with more remote residues.
Collapse
Affiliation(s)
- Brajesh Narayan
- School of Physics, University College Dublin, Belfield, Dublin 4, Ireland
| | - Arman Fathizadeh
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, 201 E. 24(th) Street, 1 University Station (C0200), Austin, TX 78712-1229, USA
| | - Clark Templeton
- McKetta Department of Chemical Engineering, University of Texas at Austin, 200 E. Dean Keaton St. Stop C0400, Austin, TX 78712-1589, USA
| | - Peng He
- Department of Chemistry, Temple University, 1801 N Broad Street, Philadelphia, PA 19122, USA
| | - Shima Arasteh
- Department of Chemistry, Temple University, 1801 N Broad Street, Philadelphia, PA 19122, USA
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences, University of Texas at Austin, 201 E. 24(th) Street, 1 University Station (C0200), Austin, TX 78712-1229, USA; Department of Chemistry, University of Texas at Austin, 2506 Speedway STOP A5300, Austin, TX 78712-1224, USA.
| | | | - Ron M Levy
- Department of Chemistry, Temple University, 1801 N Broad Street, Philadelphia, PA 19122, USA
| |
Collapse
|
14
|
Xiao BL, Ning YN, Niu NN, Li D, Moosavi-Movahedi AA, Sheibani N, Hong J. Steered molecular dynamic simulations of conformational lock of Cu, Zn-superoxide dismutase. Sci Rep 2019; 9:4353. [PMID: 30867507 PMCID: PMC6416402 DOI: 10.1038/s41598-019-40892-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2018] [Accepted: 02/25/2019] [Indexed: 01/01/2023] Open
Abstract
The conformational lock was a bio-thermodynamic theory to explain the characteristics of interfaces in oligomeric enzymes and their effects on catalytic activity. The previous studies on superoxide dismutases (Cu, Zn-SODs) showed that the dimeric structure contributed to the high catalytic efficiency and the stability. In this study, steered molecular dynamics simulations were used firstly to study the main interactions between two subunits of Cu, Zn-SODs. The decomposition process study showed that there were not only four pairs of hydrogen bonds but also twenty-five residue pairs participating hydrophobic interactions between A and B chains of SOD, and van der Waals interactions occupied a dominant position among these residue pairs. Moreover, the residue pairs of hydrogen bonds played a major role in maintaining the protein conformation. The analysis of the energy and conformational changes in the SMD simulation showed that there were two groups (two conformational locks) between A and B chains of SOD. The first group consisted of one hydrogen-bond residues pair and seven hydrophobic interactions residues pairs with a total average energy of −30.10 KJ/mol, and the second group of three hydrogen-bond residues pair and eighteen hydrophobic interactions residues pairs formed with a total average energy of −115.23 KJ/mol.
Collapse
Affiliation(s)
- Bao-Lin Xiao
- School of Life Sciences, Henan University, JinMing Road, Kaifeng, 475000, China
| | - Yan-Na Ning
- School of Life Sciences, Henan University, JinMing Road, Kaifeng, 475000, China
| | - Nan-Nan Niu
- School of Life Sciences, Henan University, JinMing Road, Kaifeng, 475000, China
| | - Di Li
- School of Life Sciences, Henan University, JinMing Road, Kaifeng, 475000, China
| | - Ali Akbar Moosavi-Movahedi
- Institute of Biochemistry and Biophysics, University of Tehran, Enquelab Avenue, P.O. Box 13145-1384, Tehran, Iran
| | - Nader Sheibani
- Department of Ophthalmology and Visual Sciences and Biomedical Engineering, University of Wisconsin, School of Medicine and Public Health, Madison, WI, 53726, USA
| | - Jun Hong
- School of Life Sciences, Henan University, JinMing Road, Kaifeng, 475000, China. .,Henan Engineering Laboratory for Mammary Bioreactor, School of Life Sciences, Henan University JinMing Road, Kaifeng, 475000, China.
| |
Collapse
|