1
|
Kim J, Moon S, Romo TD, Yang Y, Bae E, Phillips GN. Conformational dynamics of adenylate kinase in crystals. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2024; 11:014702. [PMID: 38389978 PMCID: PMC10883716 DOI: 10.1063/4.0000205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 12/14/2023] [Indexed: 02/24/2024]
Abstract
Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.
Collapse
Affiliation(s)
- Junhyung Kim
- Department of Agricultural Biotechnology, Seoul National University, Seoul 08826, South Korea
| | - Sojin Moon
- Department of Agricultural Biotechnology, Seoul National University, Seoul 08826, South Korea
| | - Tod D Romo
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, New York 14642, USA
| | - Yifei Yang
- Departments of BioSciences, Rice University, Houston, Texas 77005, USA
| | | | | |
Collapse
|
2
|
Curado‐Carballada C, Feixas F, Osuna S. Molecular Dynamics Simulations on
Aspergillus niger
Monoamine Oxidase: Conformational Dynamics and Inter‐monomer Communication Essential for Its Efficient Catalysis. Adv Synth Catal 2019. [DOI: 10.1002/adsc.201900158] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Christian Curado‐Carballada
- CompBioLab group, Institut de Química Computacional i Catàlisi (IQCC)Departament de Química Carrer Maria Aurèlia Capmany 69 17003 Girona, Catalonia Spain
| | - Ferran Feixas
- CompBioLab group, Institut de Química Computacional i Catàlisi (IQCC)Departament de Química Carrer Maria Aurèlia Capmany 69 17003 Girona, Catalonia Spain
| | - Sílvia Osuna
- CompBioLab group, Institut de Química Computacional i Catàlisi (IQCC)Departament de Química Carrer Maria Aurèlia Capmany 69 17003 Girona, Catalonia Spain
- ICREAPg. Lluís Companys 23 08010 Barcelona Spain
| |
Collapse
|
3
|
Ye C, Ding C, Ma R, Wang J, Zhang Z. Electrostatic interactions determine entrance/release order of substrates in the catalytic cycle of adenylate kinase. Proteins 2019; 87:337-347. [PMID: 30615212 DOI: 10.1002/prot.25655] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2018] [Revised: 01/01/2019] [Accepted: 01/02/2019] [Indexed: 12/29/2022]
Abstract
Adenylate kinase is a monomeric phosphotransferase with important biological function in regulating concentration of adenosine triphosphate (ATP) in cells, by transferring the terminal phosphate group from ATP to adenosine monophosphate (AMP) and forming two adenosine diphosphate (ADP) molecules. During this reaction, the kinase may undergo a large conformational transition, forming different states with its substrates. Although many structures of the protein are available, atomic details of the whole process remain unclear. In this article, we use both conventional molecular dynamics (MD) simulation and an enhanced sampling technique called parallel cascade selection MD simulation to explore different conformational states of the Escherichia coli adenylate kinase. Based on the simulation results, we propose a possible entrance/release order of substrates during the catalytic cycle. The substrate-free protein prefers an open conformation, but changes to a closed state once ATP·Mg enters into its binding pocket first and then AMP does. After the reaction of ATP transferring the terminal phosphate group to AMP, ADP·Mg and ADP are released sequentially, and finally the whole catalyze cycle is completed. Detailed contact and distance analysis reveals that the entrance/release order of substrates may be largely controlled by electrostatic interactions between the protein and the substrates.
Collapse
Affiliation(s)
- Chun Ye
- High Magnetic Field Laboratory, Chinese Academy of Sciences, Hefei, Anhui, China.,Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China
| | - Chengtao Ding
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China
| | - Rongsheng Ma
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China
| | - Junfeng Wang
- High Magnetic Field Laboratory, Chinese Academy of Sciences, Hefei, Anhui, China
| | - Zhiyong Zhang
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, China
| |
Collapse
|
4
|
Zheng Y, Cui Q. Multiple Pathways and Time Scales for Conformational Transitions in apo-Adenylate Kinase. J Chem Theory Comput 2018; 14:1716-1726. [PMID: 29378407 DOI: 10.1021/acs.jctc.7b01064] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The open/close transition in adenylate kinase (AK) is regarded as a representative example for large-scale conformational transition in proteins, yet its mechanism remains unclear despite numerous experimental and computational studies. Using extensive (∼50 μs) explicit solvent atomistic simulations and Markov state analysis, we shed new lights on the mechanism of this transition in the apo form of AK. The closed basin of apo AK features an open NMP domain while the LID domain closes and rotates toward it. Therefore, although the computed structural properties of the closed ensemble are consistent with previously reported FRET and PRE measurements, our simulations suggest that NMP closure is likely to follow AMP binding, in contrast to the previous interpretation of FRET and PRE data that the apo state was able to sample the fully closed conformation for "ligand selection". The closed state ensemble is found to be kinetically heterogeneous; multiple pathways and time scales are associated with the open/close transition, providing new clues to the disparate time scales observed in different experiments. Besides interdomain interactions, a novel mutual information analysis identifies specific intradomain interactions that correlate strongly to transition kinetics, supporting observations from previous chimera experiments. While our results underscore the role of internal domain properties in determining the kinetics of open/close transition in apo AK, no evidence is observed for any significant degree of local unfolding during the transition. These observations about AK have general implications to our view of conformational states, transition pathways, and time scales of conformational changes in proteins. The key features and time scales of observed transition pathways are robust and similar from simulations using two popular fixed charge force fields.
Collapse
Affiliation(s)
- Yuqing Zheng
- Graduate Program in Biophysics and Department of Chemistry , University of Wisconsin-Madison , 1101 University Avenue , Madison , Wisconsin 53706 , United States
| | - Qiang Cui
- Graduate Program in Biophysics and Department of Chemistry , University of Wisconsin-Madison , 1101 University Avenue , Madison , Wisconsin 53706 , United States
| |
Collapse
|
5
|
Hashem S, Tiberti M, Fornili A. Allosteric modulation of cardiac myosin dynamics by omecamtiv mecarbil. PLoS Comput Biol 2017; 13:e1005826. [PMID: 29108014 PMCID: PMC5690683 DOI: 10.1371/journal.pcbi.1005826] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2017] [Revised: 11/16/2017] [Accepted: 10/16/2017] [Indexed: 01/10/2023] Open
Abstract
New promising avenues for the pharmacological treatment of skeletal and heart muscle diseases rely on direct sarcomeric modulators, which are molecules that can directly bind to sarcomeric proteins and either inhibit or enhance their activity. A recent breakthrough has been the discovery of the myosin activator omecamtiv mecarbil (OM), which has been shown to increase the power output of the cardiac muscle and is currently in clinical trials for the treatment of heart failure. While the overall effect of OM on the mechano-chemical cycle of myosin is to increase the fraction of myosin molecules in the sarcomere that are strongly bound to actin, the molecular basis of its action is still not completely clear. We present here a Molecular Dynamics study of the motor domain of human cardiac myosin bound to OM, where the effects of the drug on the dynamical properties of the protein are investigated for the first time with atomistic resolution. We found that OM has a double effect on myosin dynamics, inducing a) an increased coupling of the motions of the converter and lever arm subdomains to the rest of the protein and b) a rewiring of the network of dynamic correlations, which produces preferential communication pathways between the OM binding site and distant functional regions. The location of the residues responsible for these effects suggests possible strategies for the future development of improved drugs and the targeting of specific cardiomyopathy-related mutations.
Collapse
Affiliation(s)
- Shaima Hashem
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
| | - Matteo Tiberti
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
| | - Arianna Fornili
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
- The Thomas Young Centre for Theory and Simulation of Materials, London, United Kingdom
- * E-mail:
| |
Collapse
|
6
|
Lu X, Ovchinnikov V, Demapan D, Roston D, Cui Q. Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin. Biochemistry 2017; 56:1482-1497. [PMID: 28225609 DOI: 10.1021/acs.biochem.7b00016] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis. Comparison of the observed hydrolysis pathways and the corresponding free energy profiles leads to detailed models for the mechanism of ATP hydrolysis in the pre-powerstroke state and proposes factors that regulate the hydrolysis activity in different conformational states. In the pre-powerstroke state, the scissile Pγ-O3β bond breaks early in the reaction. Proton transfer from the lytic water to the γ-phosphate through active site residues is an important part of the kinetic bottleneck; several hydrolysis pathways that feature distinct proton transfer routes are found to have similar free energy barriers, suggesting a significant degree of plasticity in the hydrolysis mechanism. Comparison of hydrolysis in the pre-powerstroke state and the closed post-rigor model suggests that optimization of residues beyond the active site for electrostatic stabilization and preorganization is likely important to enzyme design.
Collapse
Affiliation(s)
- Xiya Lu
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Victor Ovchinnikov
- Department of Chemistry and Chemical Biology, Harvard University , 12 Oxford Street, Boston, Massachusetts 02138, United States
| | - Darren Demapan
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Daniel Roston
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
7
|
Highly selective inhibition of myosin motors provides the basis of potential therapeutic application. Proc Natl Acad Sci U S A 2016; 113:E7448-E7455. [PMID: 27815532 DOI: 10.1073/pnas.1609342113] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Direct inhibition of smooth muscle myosin (SMM) is a potential means to treat hypercontractile smooth muscle diseases. The selective inhibitor CK-2018571 prevents strong binding to actin and promotes muscle relaxation in vitro and in vivo. The crystal structure of the SMM/drug complex reveals that CK-2018571 binds to a novel allosteric pocket that opens up during the "recovery stroke" transition necessary to reprime the motor. Trapped in an intermediate of this fast transition, SMM is inhibited with high selectivity compared with skeletal muscle myosin (IC50 = 9 nM and 11,300 nM, respectively), although all of the binding site residues are identical in these motors. This structure provides a starting point from which to design highly specific myosin modulators to treat several human diseases. It further illustrates the potential of targeting transition intermediates of molecular machines to develop exquisitely selective pharmacological agents.
Collapse
|
8
|
Li D, Liu MS, Ji B. Mapping the Dynamics Landscape of Conformational Transitions in Enzyme: The Adenylate Kinase Case. Biophys J 2016; 109:647-60. [PMID: 26244746 DOI: 10.1016/j.bpj.2015.06.059] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2014] [Revised: 06/19/2015] [Accepted: 06/29/2015] [Indexed: 12/22/2022] Open
Abstract
Conformational transition describes the essential dynamics and mechanism of enzymes in pursuing their various functions. The fundamental and practical challenge to researchers is to quantitatively describe the roles of large-scale dynamic transitions for regulating the catalytic processes. In this study, we tackled this challenge by exploring the pathways and free energy landscape of conformational changes in adenylate kinase (AdK), a key ubiquitous enzyme for cellular energy homeostasis. Using explicit long-timescale (up to microseconds) molecular dynamics and bias-exchange metadynamics simulations, we determined at the atomistic level the intermediate conformational states and mapped the transition pathways of AdK in the presence and absence of ligands. There is clearly chronological operation of the functional domains of AdK. Specifically in the ligand-free AdK, there is no significant energy barrier in the free energy landscape separating the open and closed states. Instead there are multiple intermediate conformational states, which facilitate the rapid transitions of AdK. In the ligand-bound AdK, the closed conformation is energetically most favored with a large energy barrier to open it up, and the conformational population prefers to shift to the closed form coupled with transitions. The results suggest a perspective for a hybrid of conformational selection and induced fit operations of ligand binding to AdK. These observations, depicted in the most comprehensive and quantitative way to date, to our knowledge, emphasize the underlying intrinsic dynamics of AdK and reveal the sophisticated conformational transitions of AdK in fulfilling its enzymatic functions. The developed methodology can also apply to other proteins and biomolecular systems.
Collapse
Affiliation(s)
- Dechang Li
- Biomechanics and Biomaterials Laboratory, Department of Applied Mechanics, Beijing Institute of Technology, Beijing, China.
| | - Ming S Liu
- CSIRO - Digital Productivity Flagship, Clayton South, Victoria, Australia; Monash Institute of Medical Research, Clayton, Victoria, Australia.
| | - Baohua Ji
- Biomechanics and Biomaterials Laboratory, Department of Applied Mechanics, Beijing Institute of Technology, Beijing, China.
| |
Collapse
|
9
|
Ma W, Schulten K. Mechanism of substrate translocation by a ring-shaped ATPase motor at millisecond resolution. J Am Chem Soc 2015; 137:3031-40. [PMID: 25646698 PMCID: PMC4393844 DOI: 10.1021/ja512605w] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Ring-shaped, hexameric ATPase motors fulfill key functions in cellular processes, such as genome replication, transcription, or protein degradation, by translocating a long substrate through their central pore powered by ATP hydrolysis. Despite intense research efforts, the atomic-level mechanism transmitting chemical energy from hydrolysis into mechanical force that translocates the substrate is still unclear. Here we employ all-atom molecular dynamics simulations combined with advanced path sampling techniques and milestoning analysis to characterize how mRNA substrate is translocated by an exemplary homohexameric motor, the transcription termination factor Rho. We find that the release of hydrolysis product (ADP + Pi) triggers the force-generating process of Rho through a 0.1 millisecond-long conformational transition, the time scale seen also in experiment. The calculated free energy profiles and kinetics show that Rho unidirectionally translocates the single-stranded RNA substrate via a population shift of the conformational states of Rho; upon hydrolysis product release, the most favorable conformation shifts from the pretranslocation state to the post-translocation state. Via two previously unidentified intermediate states, the RNA chain is seen to be pulled by six K326 side chains, whose motions are induced by highly coordinated relative translation and rotation of Rho's six subunits. The present study not only reveals in new detail the mechanism employed by ring-shaped ATPase motors, for example the use of loosely bound and tightly bound hydrolysis reactant and product states to coordinate motor action, but also provides an effective approach to identify allosteric sites of multimeric enzymes in general.
Collapse
|
10
|
Ruvinsky AM, Vakser IA, Rivera M. Local packing modulates diversity of iron pathways and cooperative behavior in eukaryotic and prokaryotic ferritins. J Chem Phys 2014; 140:115104. [PMID: 24655206 DOI: 10.1063/1.4868229] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Ferritin-like molecules show a remarkable combination of the evolutionary conserved activity of iron uptake and release that engage different pores in the conserved ferritin shell. It was hypothesized that pore selection and iron traffic depend on dynamic allostery with no conformational changes in the backbone. In this study, we detect the allosteric networks in Pseudomonas aeruginosa bacterioferritin (BfrB), bacterial ferritin (FtnA), and bullfrog M and L ferritins (Ftns) by a network-weaving algorithm (NWA) that passes threads of an allosteric network through highly correlated residues using hierarchical clustering. The residue-residue correlations are calculated in the packing-on elastic network model that introduces atom packing into the common packing-off model. Applying NWA revealed that each of the molecules has an extended allosteric network mostly buried inside the ferritin shell. The structure of the networks is consistent with experimental observations of iron transport: The allosteric networks in BfrB and FtnA connect the ferroxidase center with the 4-fold pores and B-pores, leaving the 3-fold pores unengaged. In contrast, the allosteric network directly links the 3-fold pores with the 4-fold pores in M and L Ftns. The majority of the network residues are either on the inner surface or buried inside the subunit fold or at the subunit interfaces. We hypothesize that the ferritin structures evolved in a way to limit the influence of functionally unrelated events in the cytoplasm on the allosteric network to maintain stability of the translocation mechanisms. We showed that the residue-residue correlations and the resultant long-range cooperativity depend on the ferritin shell packing, which, in turn, depends on protein sequence composition. Switching from the packing-on to the packing-off model reduces correlations by 35%-38% so that no allosteric network can be found. The influence of the side-chain packing on the allosteric networks explains the diversity in mechanisms of iron traffic suggested by experimental approaches.
Collapse
Affiliation(s)
- Anatoly M Ruvinsky
- Infection Innovative Medicine, AstraZeneca R&D Boston, 35 Gatehouse Drive, Waltham, Massachusetts 02451, USA
| | - Ilya A Vakser
- Center for Bioinformatics, The University of Kansas, Lawrence, Kansas 66047, USA
| | - Mario Rivera
- Department of Chemistry, The University of Kansas, Lawrence, Kansas 66047, USA
| |
Collapse
|
11
|
Shen L, Hu H. Resolution-Adapted All-Atomic and Coarse-Grained Model for Biomolecular Simulations. J Chem Theory Comput 2014; 10:2528-36. [DOI: 10.1021/ct401029k] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Lin Shen
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
| | - Hao Hu
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
- The University of Hong Kong-Shenzhen Institute of Research and Innovation, Kejizhong Second Road, Shenzhen, China
| |
Collapse
|
12
|
Gur M, Madura JD, Bahar I. Global transitions of proteins explored by a multiscale hybrid methodology: application to adenylate kinase. Biophys J 2014; 105:1643-52. [PMID: 24094405 DOI: 10.1016/j.bpj.2013.07.058] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Revised: 06/25/2013] [Accepted: 07/11/2013] [Indexed: 12/19/2022] Open
Abstract
Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.
Collapse
Affiliation(s)
- Mert Gur
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania
| | | | | |
Collapse
|