1
|
Scheerer D, Levy D, Casier R, Riven I, Mazal H, Haran G. Enzyme activation by urea reveals the interplay between conformational dynamics and substrate binding: a single-molecule FRET study. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.09.01.610662. [PMID: 39257823 PMCID: PMC11384010 DOI: 10.1101/2024.09.01.610662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
Proteins often harness extensive motions of domains and subunits to promote their function. Deciphering how these movements impact activity is key for understanding life's molecular machinery. The enzyme adenylate kinase is an intriguing example for this relationship; it ensures efficient catalysis by large-scale domain motions that lead to the enclosure of the bound substrates ATP and AMP. At high concentrations, AMP also operates as an allosteric inhibitor of the protein. Surprisingly, the enzyme is activated by urea, a compound commonly acting as a denaturant. Combining single-molecule FRET spectroscopy and enzymatic activity studies, we find that urea interferes with two key mechanisms that contribute to enzyme efficacy. First, urea promotes the open conformation of the enzyme, aiding the proper positioning of the substrates. Second, urea decreases AMP affinity, paradoxically facilitating a more efficient progression towards the catalytically active complex. These results signify the important interplay between conformational dynamics and chemical steps, including binding, in the activity of enzymes. State-of-the-art tools, such as single-molecule fluorescence spectroscopy, offer new insights into how enzymes balance different conformations to regulate activity.
Collapse
Affiliation(s)
- David Scheerer
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Dorit Levy
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Remi Casier
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Inbal Riven
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Hisham Mazal
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
- Max Planck Institute for the Science of Light, Erlangen 91058, Germany
| | - Gilad Haran
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| |
Collapse
|
2
|
Nam K, Thodika ARA, Tischlik S, Phoeurk C, Nagy TM, Schierholz L, Ådén J, Rogne P, Drescher M, Sauer-Eriksson AE, Wolf-Watz M. Magnesium induced structural reorganization in the active site of adenylate kinase. SCIENCE ADVANCES 2024; 10:eado5504. [PMID: 39121211 PMCID: PMC11313852 DOI: 10.1126/sciadv.ado5504] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Accepted: 07/08/2024] [Indexed: 08/11/2024]
Abstract
Phosphoryl transfer is a fundamental reaction in cellular signaling and metabolism that requires Mg2+ as an essential cofactor. While the primary function of Mg2+ is electrostatic activation of substrates, such as ATP, the full spectrum of catalytic mechanisms exerted by Mg2+ is not known. In this study, we integrate structural biology methods, molecular dynamic (MD) simulations, phylogeny, and enzymology assays to provide molecular insights into Mg2+-dependent structural reorganization in the active site of the metabolic enzyme adenylate kinase. Our results demonstrate that Mg2+ induces a conformational rearrangement of the substrates (ATP and ADP), resulting in a 30° adjustment of the angle essential for reversible phosphoryl transfer, thereby optimizing it for catalysis. MD simulations revealed transitions between conformational substates that link the fluctuation of the angle to large-scale enzyme dynamics. The findings contribute detailed insight into Mg2+ activation of enzymes and may be relevant for reversible and irreversible phosphoryl transfer reactions.
Collapse
Affiliation(s)
- Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, TX 76019, USA
| | | | - Sonja Tischlik
- Department of Chemistry, Konstanz Research School Chemical Biology, University of Konstanz, 78464 Konstanz, Germany
| | - Chanrith Phoeurk
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
- Department of Bio-Engineering, Royal University of Phnom Penh, Phnom Penh, Cambodia
| | | | - Léon Schierholz
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
- Department of Molecular Biology, Umeå University, Umeå, 901 87, Sweden
| | - Jörgen Ådén
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | - Per Rogne
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | - Malte Drescher
- Department of Chemistry, Konstanz Research School Chemical Biology, University of Konstanz, 78464 Konstanz, Germany
| | | | | |
Collapse
|
3
|
Alavi Z, Casanova-Morales N, Quiroga-Roger D, Wilson CAM. Towards the understanding of molecular motors and its relationship with local unfolding. Q Rev Biophys 2024; 57:e7. [PMID: 38715547 DOI: 10.1017/s0033583524000052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2024]
Abstract
Molecular motors are machines essential for life since they convert chemical energy into mechanical work. However, the precise mechanism by which nucleotide binding, catalysis, or release of products is coupled to the work performed by the molecular motor is still not entirely clear. This is due, in part, to a lack of understanding of the role of force in the mechanical-structural processes involved in enzyme catalysis. From a mechanical perspective, one promising hypothesis is the Haldane-Pauling hypothesis which considers the idea that part of the enzymatic catalysis is strain-induced. It suggests that enzymes cannot be efficient catalysts if they are fully complementary to the substrates. Instead, they must exert strain on the substrate upon binding, using enzyme-substrate energy interaction (binding energy) to accelerate the reaction rate. A novel idea suggests that during catalysis, significant strain energy is built up, which is then released by a local unfolding/refolding event known as 'cracking'. Recent evidence has also shown that in catalytic reactions involving conformational changes, part of the heat released results in a center-of-mass acceleration of the enzyme, raising the possibility that the heat released by the reaction itself could affect the enzyme's integrity. Thus, it has been suggested that this released heat could promote or be linked to the cracking seen in proteins such as adenylate kinase (AK). We propose that the energy released as a consequence of ligand binding/catalysis is associated with the local unfolding/refolding events (cracking), and that this energy is capable of driving the mechanical work.
Collapse
Affiliation(s)
- Zahra Alavi
- Department of Physics, Loyola Marymount University, Los Angeles, CA, USA
| | | | - Diego Quiroga-Roger
- Biochemistry and Molecular Biology Department, Faculty of Chemistry and Pharmaceutical Sciences, Universidad de Chile, Santiago, Chile
| | - Christian A M Wilson
- Biochemistry and Molecular Biology Department, Faculty of Chemistry and Pharmaceutical Sciences, Universidad de Chile, Santiago, Chile
| |
Collapse
|
4
|
Nam K, Arattu Thodika AR, Grundström C, Sauer UH, Wolf-Watz M. Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release. J Chem Inf Model 2024; 64:150-163. [PMID: 38117131 PMCID: PMC10778088 DOI: 10.1021/acs.jcim.3c01618] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 12/01/2023] [Accepted: 12/07/2023] [Indexed: 12/21/2023]
Abstract
This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.
Collapse
Affiliation(s)
- Kwangho Nam
- Department
of Chemistry and Biochemistry, University
of Texas at Arlington, Arlington, Texas 76019, United States
| | - Abdul Raafik Arattu Thodika
- Department
of Chemistry and Biochemistry, University
of Texas at Arlington, Arlington, Texas 76019, United States
| | | | - Uwe H. Sauer
- Department
of Chemistry, Umeå University, Umeå 90187, SE, Sweden
| | - Magnus Wolf-Watz
- Department
of Chemistry, Umeå University, Umeå 90187, SE, Sweden
| |
Collapse
|
5
|
Tischlik S, Oelker M, Rogne P, Sauer-Eriksson AE, Drescher M, Wolf-Watz M. Insights into Enzymatic Catalysis from Binding and Hydrolysis of Diadenosine Tetraphosphate by E. coli Adenylate Kinase. Biochemistry 2023; 62:2238-2243. [PMID: 37418448 PMCID: PMC10399197 DOI: 10.1021/acs.biochem.3c00189] [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] [Received: 04/12/2023] [Revised: 07/04/2023] [Indexed: 07/09/2023]
Abstract
Adenylate kinases play a crucial role in cellular energy homeostasis through the interconversion of ATP, AMP, and ADP in all living organisms. Here, we explore how adenylate kinase (AdK) from Escherichia coli interacts with diadenosine tetraphosphate (AP4A), a putative alarmone associated with transcriptional regulation, stress, and DNA damage response. From a combination of EPR and NMR spectroscopy together with X-ray crystallography, we found that AdK interacts with AP4A with two distinct modes that occur on disparate time scales. First, AdK dynamically interconverts between open and closed states with equal weights in the presence of AP4A. On a much slower time scale, AdK hydrolyses AP4A, and we suggest that the dynamically accessed substrate-bound open AdK conformation enables this hydrolytic activity. The partitioning of the enzyme into open and closed states is discussed in relation to a recently proposed linkage between active site dynamics and collective conformational dynamics.
Collapse
Affiliation(s)
- Sonja Tischlik
- Department
of Chemistry, Konstanz Research School Chemical Biology, University of Konstanz, 78457 Konstanz, Germany
| | - Melanie Oelker
- Department
of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - Per Rogne
- Department
of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - A. Elisabeth Sauer-Eriksson
- Department
of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
- Centre
of Microbial Research (UCMR), Umeå
University, SE-901 87 Umeå, Sweden
| | - Malte Drescher
- Department
of Chemistry, Konstanz Research School Chemical Biology, University of Konstanz, 78457 Konstanz, Germany
| | - Magnus Wolf-Watz
- Department
of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
- Centre
of Microbial Research (UCMR), Umeå
University, SE-901 87 Umeå, Sweden
| |
Collapse
|
6
|
Nag R, Joshi S, Rathore AS, Majumder S. Profiling Enzyme Activity of l-Asparaginase II by NMR-Based Methyl Fingerprinting at Natural Abundance. J Am Chem Soc 2023; 145:10826-10838. [PMID: 37154467 DOI: 10.1021/jacs.3c02154] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
l-asparaginase II (MW 135 kDa) from E. coli is an FDA-approved protein drug used for the treatment of childhood leukemia. Despite its long history as a chemotherapeutic, the structural basis of enzyme action, in solution, remains widely contested. In this work, methyl-based 2D [1H-13C]-heteronuclear single-quantum correlation (HSQC) NMR, at natural abundance, has been used to profile the enzymatic activity of the commercially available enzyme drug. The [1H-13C]-HSQC NMR spectra of the protein reveal the role of a flexible loop segment in the activity of the enzyme, in solution. Addition of asparagine to the protein results in distinct conformational changes of the loop that could be signatures of intermediates formed in the catalytic reaction. To this end, an isothermal titration calorimetry (ITC)-based assay has been developed to measure the enzymatic reaction enthalpy, as a marker for its activity. Combining both ITC and NMR, it was shown that the disruption of the protein conformation can result in the loss of function. The scope, robustness, and validity of the loop fingerprints in relation to enzyme activity have been tested under different solution conditions. Overall, our results indicate that 2D NMR can be used reliably to gauge the structure-function of this enzyme, bypassing the need to label the protein. Such natural abundant NMR methods can be potentially extended to probe the structure-function aspects of high-molecular-weight protein therapeutics (glycosylated protein drugs, enzymes, therapeutic monoclonal antibodies, antibody-drug conjugates, and Fc-fusion proteins), where (a) flexible loops are required for their function and (b) isotope labeling may not be straightforward.
Collapse
Affiliation(s)
- Rachayita Nag
- Biophysics & Structural Genomics, Saha Institute of Nuclear Physics, Kolkata 700064, India
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India
| | - Srishti Joshi
- Department of Chemical Engineering, Indian Institute of Technology, Hauz Khas, New Delhi 110016, India
| | - Anurag Singh Rathore
- Department of Chemical Engineering, Indian Institute of Technology, Hauz Khas, New Delhi 110016, India
| | - Subhabrata Majumder
- Biophysics & Structural Genomics, Saha Institute of Nuclear Physics, Kolkata 700064, India
- Homi Bhabha National Institute, Anushaktinagar, Mumbai 400094, India
| |
Collapse
|
7
|
Scheerer D, Adkar BV, Bhattacharyya S, Levy D, Iljina M, Riven I, Dym O, Haran G, Shakhnovich EI. Allosteric communication between ligand binding domains modulates substrate inhibition in adenylate kinase. Proc Natl Acad Sci U S A 2023; 120:e2219855120. [PMID: 37094144 PMCID: PMC10160949 DOI: 10.1073/pnas.2219855120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2022] [Accepted: 03/22/2023] [Indexed: 04/26/2023] Open
Abstract
Enzymes play a vital role in life processes; they control chemical reactions and allow functional cycles to be synchronized. Many enzymes harness large-scale motions of their domains to achieve tremendous catalytic prowess and high selectivity for specific substrates. One outstanding example is provided by the three-domain enzyme adenylate kinase (AK), which catalyzes phosphotransfer between ATP to AMP. Here we study the phenomenon of substrate inhibition by AMP and its correlation with domain motions. Using single-molecule FRET spectroscopy, we show that AMP does not block access to the ATP binding site, neither by competitive binding to the ATP cognate site nor by directly closing the LID domain. Instead, inhibitory concentrations of AMP lead to a faster and more cooperative domain closure by ATP, leading in turn to an increased population of the closed state. The effect of AMP binding can be modulated through mutations throughout the structure of the enzyme, as shown by the screening of an extensive AK mutant library. The mutation of multiple conserved residues reduces substrate inhibition, suggesting that substrate inhibition is an evolutionary well conserved feature in AK. Combining these insights, we developed a model that explains the complex activity of AK, particularly substrate inhibition, based on the experimentally observed opening and closing rates. Notably, the model indicates that the catalytic power is affected by the microsecond balance between the open and closed states of the enzyme. Our findings highlight the crucial role of protein motions in enzymatic activity.
Collapse
Affiliation(s)
- David Scheerer
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Bharat V Adkar
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
| | | | - Dorit Levy
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Marija Iljina
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Inbal Riven
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Orly Dym
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Gilad Haran
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 761001, Israel
| | - Eugene I Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
| |
Collapse
|
8
|
Irukuvajjula SS, Jithender Reddy G, Rao K, Vadrevu LR. Contrasting effect of ficoll on apo and holo forms of bacterial chemotaxis protein Y: Selective destabilization of the conformationally altered holo form. Int J Biol Macromol 2023; 232:123505. [PMID: 36736516 DOI: 10.1016/j.ijbiomac.2023.123505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 01/13/2023] [Accepted: 01/28/2023] [Indexed: 02/04/2023]
Abstract
Chemotaxis Y (CheY), upon metal binding, displays a drastic alteration in its structure and stability. This premise prompted us to study the effect of crowding on the two conformationally distinct states of the same test protein. A comparative analysis on the structure and thermal stability in the presence and absence of the macromolecular crowder, ficoll, and its monomeric unit, sucrose, revealed a contrasting effect of ficoll on the apo and holo forms. In the presence of ficoll while the thermal stability (Tm) of the apo form is enhanced, the thermal stability of the holo form is reduced. The selective lowering of Tm for the holo form in the combined presence of ficoll and sucrose and not in sucrose alone suggests that the contrasting effect is due to the macromolecular nature of ficoll. Since metal-protein interaction remains unperturbed in the presence of ficoll and Mg2+ sequestration is ruled out in a systematic manner the alternative possibility for the exclusive reduction in the thermal stability of the holo form is the ficoll-induced modulation of the relative population of apo and holo forms of CheY.
Collapse
Affiliation(s)
- Shivkumar Sharma Irukuvajjula
- Department of Biological Sciences, Birla Institute of Science and Technology - Pilani, Hyderabad Campus, Shamirpet, Hyderabad 500078, India.
| | - G Jithender Reddy
- NMR Division, Department of Analytical & Structural Chemistry, CSIR-Indian Institute of Chemical Technology, Ministry of Science and Technology, Uppal Road, Tarnaka, Hyderabad 500007, India
| | - Krishna Rao
- Tata Institute of Fundamental Research, 36/P, Gopanpally Mandal, Ranga Reddy District, Hyderabad, Telangana State 500107, India
| | - Late Ramakrishna Vadrevu
- Department of Biological Sciences, Birla Institute of Science and Technology - Pilani, Hyderabad Campus, Shamirpet, Hyderabad 500078, India
| |
Collapse
|
9
|
Verma A, Åberg-Zingmark E, Sparrman T, Mushtaq AU, Rogne P, Grundström C, Berntsson R, Sauer UH, Backman L, Nam K, Sauer-Eriksson E, Wolf-Watz M. Insights into the evolution of enzymatic specificity and catalysis: From Asgard archaea to human adenylate kinases. SCIENCE ADVANCES 2022; 8:eabm4089. [PMID: 36332013 PMCID: PMC9635829 DOI: 10.1126/sciadv.abm4089] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Accepted: 09/15/2022] [Indexed: 06/16/2023]
Abstract
Enzymatic catalysis is critically dependent on selectivity, active site architecture, and dynamics. To contribute insights into the interplay of these properties, we established an approach with NMR, crystallography, and MD simulations focused on the ubiquitous phosphotransferase adenylate kinase (AK) isolated from Odinarchaeota (OdinAK). Odinarchaeota belongs to the Asgard archaeal phylum that is believed to be the closest known ancestor to eukaryotes. We show that OdinAK is a hyperthermophilic trimer that, contrary to other AK family members, can use all NTPs for its phosphorylation reaction. Crystallographic structures of OdinAK-NTP complexes revealed a universal NTP-binding motif, while 19F NMR experiments uncovered a conserved and rate-limiting dynamic signature. As a consequence of trimerization, the active site of OdinAK was found to be lacking a critical catalytic residue and is therefore considered to be "atypical." On the basis of discovered relationships with human monomeric homologs, our findings are discussed in terms of evolution of enzymatic substrate specificity and cold adaptation.
Collapse
Affiliation(s)
- Apoorv Verma
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | | | - Tobias Sparrman
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | | | - Per Rogne
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | | | - Ronnie Berntsson
- Department of Medical Biochemistry and Biophysics, Umeå University, 901 87 Umeå, Sweden
- Wallenberg Centre for Molecular Medicine, Umeå University, 901 87 Umeå, Sweden
| | - Uwe H. Sauer
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | - Lars Backman
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, TX 76019, USA
| | | | | |
Collapse
|
10
|
Lu J, Scheerer D, Haran G, Li W, Wang W. Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase. J Phys Chem B 2022; 126:8188-8201. [PMID: 36222098 PMCID: PMC9589722 DOI: 10.1021/acs.jpcb.2c05497] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The catalytic cycle of the enzyme adenylate kinase involves large conformational motions between open and closed states. A previous single-molecule experiment showed that substrate binding tends to accelerate both the opening and the closing rates and that a single turnover event often involves multiple rounds of conformational switching. In this work, we showed that the repeated conformational transitions of adenylate kinase are essential for the relaxation of incorrectly bound substrates into the catalytically competent conformation by combining all-atom and coarse-grained molecular simulations. In addition, free energy calculations based on all-atom and coarse-grained models demonstrated that the enzyme with incorrectly bound substrates has much a lower free energy barrier for domain opening compared to that with the correct substrate conformation, which may explain the the acceleration of the domain opening rate by substrate binding. The results of this work provide mechanistic understanding to previous experimental observations and shed light onto the interplay between conformational dynamics and enzyme catalysis.
Collapse
Affiliation(s)
- Jiajun Lu
- Department
of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing210093, China,Wenzhou
Key Laboratory of Biophysics, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang325000, China
| | - David Scheerer
- Department
of Chemical and Biological Physics, Weizmann
Institute of Science, Rehovot761001, Israel
| | - Gilad Haran
- Department
of Chemical and Biological Physics, Weizmann
Institute of Science, Rehovot761001, Israel,
| | - Wenfei Li
- Department
of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing210093, China,Wenzhou
Key Laboratory of Biophysics, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang325000, China,
| | - Wei Wang
- Department
of Physics, National Laboratory of Solid State Microstructure, Nanjing University, Nanjing210093, China,
| |
Collapse
|
11
|
Zhang Y, Chen M, Lu J, Li W, Wolynes PG, Wang W. Frustration and the Kinetic Repartitioning Mechanism of Substrate Inhibition in Enzyme Catalysis. J Phys Chem B 2022; 126:6792-6801. [PMID: 36044985 PMCID: PMC9483917 DOI: 10.1021/acs.jpcb.2c03832] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
![]()
Substrate inhibition, whereby enzymatic activity decreases
with
excess substrate after reaching a maximum turnover rate, is among
the most elusive phenomena in enzymatic catalysis. Here, based on
a dynamic energy landscape model, we investigate the underlying mechanism
by performing molecular simulations and frustration analysis for a
model enzyme adenylate kinase (AdK), which catalyzes the phosphoryl
transfer reaction ATP + AMP ⇋ ADP + ADP. Intriguingly, these
reveal a kinetic repartitioning mechanism of substrate inhibition,
whereby excess substrate AMP suppresses the population of an energetically
frustrated, but kinetically activated, catalytic pathway going through
a substrate (ATP)-product (ADP) cobound complex with steric incompatibility.
Such a frustrated pathway plays a crucial role in facilitating the
bottleneck product ADP release, and its suppression by excess substrate
AMP leads to a slow down of product release and overall turnover.
The simulation results directly demonstrate that substrate inhibition
arises from the rate-limiting product-release step, instead of the
steps for populating the catalytically competent complex as often
suggested in previous works. Furthermore, there is a tight interplay
between the enzyme conformational equilibrium and the extent of substrate
inhibition. Mutations biasing to more closed conformations tend to
enhance substrate inhibition. We also characterized the key features
of single-molecule enzyme kinetics with substrate inhibition effect.
We propose that the above molecular mechanism of substrate inhibition
may be relevant to other multisubstrate enzymes in which product release
is the bottleneck step.
Collapse
Affiliation(s)
- Yangyang Zhang
- Department of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China.,Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
| | - Mingchen Chen
- Department of Research and Development, neoX Biotech, Beijing 102206, China.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Jiajun Lu
- Department of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| | - Wenfei Li
- Department of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China.,Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Wei Wang
- Department of Physics, National Laboratory of Solid State Microstructure, and Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China
| |
Collapse
|
12
|
Machine learning/molecular dynamic protein structure prediction approach to investigate the protein conformational ensemble. Sci Rep 2022; 12:10018. [PMID: 35705565 PMCID: PMC9200820 DOI: 10.1038/s41598-022-13714-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2021] [Accepted: 05/11/2022] [Indexed: 11/25/2022] Open
Abstract
Proteins exist in several different conformations. These structural changes are often associated with fluctuations at the residue level. Recent findings show that co-evolutionary analysis coupled with machine-learning techniques improves the precision by providing quantitative distance predictions between pairs of residues. The predicted statistical distance distribution from Multi Sequence Analysis reveals the presence of different local maxima suggesting the flexibility of key residue pairs. Here we investigate the ability of the residue-residue distance prediction to provide insights into the protein conformational ensemble. We combine deep learning approaches with mechanistic modeling to a set of proteins that experimentally showed conformational changes. The predicted protein models were filtered based on energy scores, RMSD clustering, and the centroids selected as the lowest energy structure per cluster. These models were compared to the experimental-Molecular Dynamics (MD) relaxed structure by analyzing the backbone residue torsional distribution and the sidechain orientations. Our pipeline allows to retrieve the experimental structural dynamics experimentally represented by different X-ray conformations for the same sequence as well the conformational space observed with the MD simulations. We show the potential correlation between the experimental structure dynamics and the predicted model ensemble demonstrating the susceptibility of the current state-of-the-art methods in protein folding and dynamics prediction and pointing out the areas of improvement.
Collapse
|
13
|
Stiller JB, Otten R, Häussinger D, Rieder PS, Theobald DL, Kern D. Structure determination of high-energy states in a dynamic protein ensemble. Nature 2022; 603:528-535. [PMID: 35236984 PMCID: PMC9126080 DOI: 10.1038/s41586-022-04468-9] [Citation(s) in RCA: 39] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Accepted: 01/25/2022] [Indexed: 01/24/2023]
Abstract
Macromolecular function frequently requires that proteins change conformation into high-energy states1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts5 (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion6 (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.
Collapse
Affiliation(s)
- John B Stiller
- Department of Biochemistry and Howard Hughes Medical Institute, Brandeis University, Waltham, MA, USA
| | - Renee Otten
- Department of Biochemistry and Howard Hughes Medical Institute, Brandeis University, Waltham, MA, USA
| | | | - Pascal S Rieder
- Department of Chemistry, University of Basel, Basel, Switzerland
| | | | - Dorothee Kern
- Department of Biochemistry and Howard Hughes Medical Institute, Brandeis University, Waltham, MA, USA.
| |
Collapse
|
14
|
Orädd F, Ravishankar H, Goodman J, Rogne P, Backman L, Duelli A, Nors Pedersen M, Levantino M, Wulff M, Wolf-Watz M, Andersson M. Tracking the ATP-binding response in adenylate kinase in real time. SCIENCE ADVANCES 2021; 7:eabi5514. [PMID: 34788091 PMCID: PMC8597995 DOI: 10.1126/sciadv.abi5514] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Accepted: 09/27/2021] [Indexed: 05/25/2023]
Abstract
The biological function of proteins is critically dependent on dynamics inherent to the native structure. Such structural dynamics obey a predefined order and temporal timing to execute the specific reaction. Determination of the cooperativity of key structural rearrangements requires monitoring protein reactions in real time. In this work, we used time-resolved x-ray solution scattering (TR-XSS) to visualize structural changes in the Escherichia coli adenylate kinase (AdK) enzyme upon laser-induced activation of a protected ATP substrate. A 4.3-ms transient intermediate showed partial closing of both the ATP- and AMP-binding domains, which indicates a cooperative closing mechanism. The ATP-binding domain also showed local unfolding and breaking of an Arg131-Asp146 salt bridge. Nuclear magnetic resonance spectroscopy data identified similar unfolding in an Arg131Ala AdK mutant, which refolded in a closed, substrate-binding conformation. The observed structural dynamics agree with a “cracking mechanism” proposed to underlie global structural transformation, such as allostery, in proteins.
Collapse
Affiliation(s)
- Fredrik Orädd
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| | - Harsha Ravishankar
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| | - Jack Goodman
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| | - Per Rogne
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| | - Lars Backman
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| | - Annette Duelli
- Department of Biomedical Sciences, University of Copenhagen, Blegdamsvej 3, 2200 Copenhagen, Denmark
| | - Martin Nors Pedersen
- ESRF—The European Synchrotron, 71 Avenue des Martyrs, CS40220, 38043 Grenoble, Cedex 9, France
| | - Matteo Levantino
- ESRF—The European Synchrotron, 71 Avenue des Martyrs, CS40220, 38043 Grenoble, Cedex 9, France
| | - Michael Wulff
- ESRF—The European Synchrotron, 71 Avenue des Martyrs, CS40220, 38043 Grenoble, Cedex 9, France
| | - Magnus Wolf-Watz
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| | - Magnus Andersson
- Department of Chemistry, Umeå University, Linnaeus Väg 10, 901 87 Umeå, Sweden
| |
Collapse
|
15
|
Ojeda-May P, Mushtaq AU, Rogne P, Verma A, Ovchinnikov V, Grundström C, Dulko-Smith B, Sauer UH, Wolf-Watz M, Nam K. Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions. Biochemistry 2021; 60:2246-2258. [PMID: 34250801 PMCID: PMC8297476 DOI: 10.1021/acs.biochem.1c00221] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
![]()
Enzymes employ a wide range of protein motions to achieve efficient catalysis of
chemical reactions. While the role of collective protein motions in substrate binding,
product release, and regulation of enzymatic activity is generally understood, their
roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations,
enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are
combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the
roles of catalytic residues in catalysis and the conformational change in the enzyme.
This study reveals that the motions in the active site, which occur on a time scale of
picoseconds to nanoseconds, link the catalytic reaction to the slow conformational
dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In
particular, substantial conformational rearrangement occurs in the active site following
the catalytic reaction. This rearrangement not only affects the reaction barrier but
also promotes a more open conformation of the enzyme after the reaction, which then
results in an accelerated opening of the enzyme compared to that of the reactant state.
The results illustrate a linkage between enzymatic catalysis and collective protein
motions, whereby the disparate time scales between the two processes are bridged by a
cascade of intermediate-scale motion of catalytic residues modulating the free energy
landscapes of the catalytic and conformational change processes.
Collapse
Affiliation(s)
- Pedro Ojeda-May
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden.,High Performance Computing Centre North (HPC2N), Umeå University, Umeå SE-90187, Sweden
| | | | - Per Rogne
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden
| | - Apoorv Verma
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden
| | - Victor Ovchinnikov
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | | | - Beata Dulko-Smith
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Uwe H Sauer
- Department of Chemistry, Umeå University, Umeå SE-90187, Sweden
| | | | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| |
Collapse
|
16
|
Wang L, Zhang J, Han M, Zhang L, Chen C, Huang A, Xie R, Wang G, Zhu J, Wang Y, Liu X, Zhuang W, Li Y, Wang J. A Genetically Encoded Two‐Dimensional Infrared Probe for Enzyme Active‐Site Dynamics. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202016880] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Li Wang
- School of Life Sciences University of Chinese Academy of Sciences Yuquan Road, Shijingshan District Beijing 100049 China
- Institute of Biophysics Chinese Academy of Sciences Datun Road, Chaoyang District Beijing 100101 China
| | - Jia Zhang
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics Institute of Physics Chinese Academy of Sciences Beijing 100190 China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing 100049 China
| | - Ming‐Jie Han
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences West 7th Avenue, Tianjin Airport Economic Area Tianjin 300308 China
- Shenzhen Institute of Transfusion Medicine Shenzhen Blood Center Futian District Shenzhen 518052 China
| | - Lu Zhang
- State Key Laboratory of Structural Chemistry Fujian Institute of Research on the Structure of Matter Chinese Academy of Sciences Fuzhou 350002 China
| | - Chao Chen
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences West 7th Avenue, Tianjin Airport Economic Area Tianjin 300308 China
| | - Aiping Huang
- Tianjin Institute of Industrial Biotechnology Chinese Academy of Sciences West 7th Avenue, Tianjin Airport Economic Area Tianjin 300308 China
| | - Ruipei Xie
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics Institute of Physics Chinese Academy of Sciences Beijing 100190 China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing 100049 China
| | - Guosheng Wang
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics Institute of Physics Chinese Academy of Sciences Beijing 100190 China
| | - Jiangrui Zhu
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics Institute of Physics Chinese Academy of Sciences Beijing 100190 China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing 100049 China
| | - Yuchuan Wang
- Shenzhen Institute of Transfusion Medicine Shenzhen Blood Center Futian District Shenzhen 518052 China
| | - Xiaohong Liu
- Institute of Biophysics Chinese Academy of Sciences Datun Road, Chaoyang District Beijing 100101 China
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry Fujian Institute of Research on the Structure of Matter Chinese Academy of Sciences Fuzhou 350002 China
- Institute of urban environment Chinese Academy of Sciences Xiamen Fujian 361021 China
| | - Yunliang Li
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics Institute of Physics Chinese Academy of Sciences Beijing 100190 China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing 100049 China
- Songshan Lake Materials Laboratory Dongguan Guangdong 523808 China
| | - Jiangyun Wang
- Institute of Biophysics Chinese Academy of Sciences Datun Road, Chaoyang District Beijing 100101 China
- Shenzhen Institute of Transfusion Medicine Shenzhen Blood Center Futian District Shenzhen 518052 China
- School of Life Sciences University of Chinese Academy of Sciences Yuquan Road, Shijingshan District Beijing 100049 China
| |
Collapse
|
17
|
Shinobu A, Kobayashi C, Matsunaga Y, Sugita Y. Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins. J Chem Inf Model 2021; 61:2427-2443. [PMID: 33956432 DOI: 10.1021/acs.jcim.1c00286] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.
Collapse
Affiliation(s)
- Ai Shinobu
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Yasuhiro Matsunaga
- Graduate School of Science and Engineering, Saitama University, Saitama 338-8570, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan.,Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan.,Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
| |
Collapse
|
18
|
Wang L, Zhang J, Han MJ, Zhang L, Chen C, Huang A, Xie R, Wang G, Zhu J, Wang Y, Liu X, Zhuang W, Li Y, Wang J. A Genetically Encoded Two-Dimensional Infrared Probe for Enzyme Active-Site Dynamics. Angew Chem Int Ed Engl 2021; 60:11143-11147. [PMID: 33644946 DOI: 10.1002/anie.202016880] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2020] [Revised: 02/10/2021] [Indexed: 11/08/2022]
Abstract
While two-dimensional infrared (2D-IR) spectroscopy is uniquely suitable for monitoring femtosecond (fs) to picosecond (ps) water dynamics around static protein structures, its utility for probing enzyme active-site dynamics is limited due to the lack of site-specific 2D-IR probes. We demonstrate the genetic incorporation of a novel 2D-IR probe, m-azido-L-tyrosine (N3Y) in the active-site of DddK, an iron-dependent enzyme that catalyzes the conversion of dimethylsulfoniopropionate to dimethylsulphide. Our results show that both the oxidation of active-site iron to FeIII , and the addition of denaturation reagents, result in significant decrease in enzyme activity and active-site water motion confinement. As tyrosine residues play important roles, including as general acids and bases, and electron transfer agents in many key enzymes, the genetically encoded 2D-IR probe N3Y should be broadly applicable to investigate how the enzyme active-site motions at the fs-ps time scale direct reaction pathways to facilitating specific chemical reactions.
Collapse
Affiliation(s)
- Li Wang
- School of Life Sciences, University of Chinese Academy of Sciences, Yuquan Road, Shijingshan District, Beijing, 100049, China.,Institute of Biophysics, Chinese Academy of Sciences, Datun Road, Chaoyang District, Beijing, 100101, China
| | - Jia Zhang
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China.,School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Ming-Jie Han
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, West 7th Avenue, Tianjin Airport Economic Area, Tianjin, 300308, China.,Shenzhen Institute of Transfusion Medicine, Shenzhen Blood Center, Futian District, Shenzhen, 518052, China
| | - Lu Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, 350002, China
| | - Chao Chen
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, West 7th Avenue, Tianjin Airport Economic Area, Tianjin, 300308, China
| | - Aiping Huang
- Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, West 7th Avenue, Tianjin Airport Economic Area, Tianjin, 300308, China
| | - Ruipei Xie
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China.,School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Guosheng Wang
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Jiangrui Zhu
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China.,School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Yuchuan Wang
- Shenzhen Institute of Transfusion Medicine, Shenzhen Blood Center, Futian District, Shenzhen, 518052, China
| | - Xiaohong Liu
- Institute of Biophysics, Chinese Academy of Sciences, Datun Road, Chaoyang District, Beijing, 100101, China
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, 350002, China.,Institute of urban environment, Chinese Academy of Sciences, Xiamen, Fujian, 361021, China
| | - Yunliang Li
- Beijing National Laboratory for Condensed Matter Physics and Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China.,School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China.,Songshan Lake Materials Laboratory, Dongguan, Guangdong, 523808, China
| | - Jiangyun Wang
- Institute of Biophysics, Chinese Academy of Sciences, Datun Road, Chaoyang District, Beijing, 100101, China.,Shenzhen Institute of Transfusion Medicine, Shenzhen Blood Center, Futian District, Shenzhen, 518052, China.,School of Life Sciences, University of Chinese Academy of Sciences, Yuquan Road, Shijingshan District, Beijing, 100049, China
| |
Collapse
|
19
|
Zhang Y, Cao Z, Zhang JZ, Xia F. Double-Well Ultra-Coarse-Grained Model to Describe Protein Conformational Transitions. J Chem Theory Comput 2020; 16:6678-6689. [PMID: 32926616 DOI: 10.1021/acs.jctc.0c00551] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The double-well model is usually used to describe the conformational transition between two states of a protein. Since conformational changes usually occur within a relatively large time scale, coarse-grained models are often used to accelerate the dynamic process due to their inexpensive computational cost. In this work, we develop a double-well ultra-coarse-grained (DW-UCG) model to describe the conformational transitions of the adenylate kinase, glutamine-binding protein, and lactoferrin. The coarse-grained simulation results show that the DW-UCG model of adenylate kinase captures the crucial intermediate states in the LID-closing and NMP-closing pathways, reflecting the key secondary structural changes in the conformational transition. A comparison of the different DW-UCG models of adenylate kinase indicates that an appropriate choice of bead resolution could generate the free energy landscape that is comparable to that from the residue-based model. The coarse-grained simulations for the glutamine-binding protein and lactoferrin also demonstrate that the DW-UCG model is valid in reproducing the correct two-state behavior for their functional study, which indicates the potential application of the DW-UCG model in investigating the mechanism of conformational changes of large proteins.
Collapse
Affiliation(s)
- Yuwei Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemistry Engineering, Xiamen University, Xiamen 361005, China.,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Zexing Cao
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemistry Engineering, Xiamen University, Xiamen 361005, China
| | - John Zenghui Zhang
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Fei Xia
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
20
|
Rogne P, Dulko-Smith B, Goodman J, Rosselin M, Grundström C, Hedberg C, Nam K, Sauer-Eriksson AE, Wolf-Watz M. Structural Basis for GTP versus ATP Selectivity in the NMP Kinase AK3. Biochemistry 2020; 59:3570-3581. [PMID: 32822537 DOI: 10.1021/acs.biochem.0c00549] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
ATP and GTP are exceptionally important molecules in biology with multiple, and often discrete, functions. Therefore, enzymes that bind to either of them must develop robust mechanisms to selectively utilize one or the other. Here, this specific problem is addressed by molecular studies of the human NMP kinase AK3, which uses GTP to phosphorylate AMP. AK3 plays an important role in the citric acid cycle, where it is responsible for GTP/GDP recycling. By combining a structural biology approach with functional experiments, we present a comprehensive structural and mechanistic understanding of the enzyme. We discovered that AK3 functions by recruitment of GTP to the active site, while ATP is rejected and nonproductively bound to the AMP binding site. Consequently, ATP acts as an inhibitor with respect to GTP and AMP. The overall features with specific recognition of the correct substrate and nonproductive binding by the incorrect substrate bear a strong similarity to previous findings for the ATP specific NMP kinase adenylate kinase. Taken together, we are now able to provide the fundamental principles for GTP and ATP selectivity in the large NMP kinase family. As a side-result originating from nonlinearity of chemical shifts in GTP and ATP titrations, we find that protein surfaces offer a general and weak binding affinity for both GTP and ATP. These nonspecific interactions likely act to lower the available intracellular GTP and ATP concentrations and may have driven evolution of the Michaelis constants of NMP kinases accordingly.
Collapse
Affiliation(s)
- Per Rogne
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - Beata Dulko-Smith
- Department of Chemistry and Biochemistry, University of Texas at Alington, Arlington, Texas 76019-0065, United States
| | - Jack Goodman
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - Marie Rosselin
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | | | | | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Alington, Arlington, Texas 76019-0065, United States
| | | | - Magnus Wolf-Watz
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| |
Collapse
|
21
|
Yuan Y, Zhu Q, Song R, Ma J, Dong H. A Two-Ended Data-Driven Accelerated Sampling Method for Exploring the Transition Pathways between Two Known States of Protein. J Chem Theory Comput 2020; 16:4631-4640. [PMID: 32320614 DOI: 10.1021/acs.jctc.9b01184] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Conformational transitions of protein between different states are often associated with their biological functions. These dynamic processes, however, are usually not easy to be well characterized by experimental measurements, mainly because of inadequate temporal and spatial resolution. Meantime, sampling of configuration space with molecular dynamics (MD) simulations is still a challenge. Here we proposed a robust two-ended data-driven accelerated (teDA2) conformational sampling method, which drives the structural change in an adaptively updated feature space without introducing a bias potential. teDA2 was applied to explore adenylate kinase (ADK), a model with well characterized "open" and "closed" states. A single conformational transition event of ADK could be achieved within only a few or tens of nanoseconds sampled with teDA2. By analyzing hundreds of transition events, we reproduced different mechanisms and the associated pathways for domain motion of ADK reported in the literature. The multiroute characteristic of ADK was confirmed by the fact that some metastable states identified with teDA2 resemble available crystal structures determined at different conditions. This feature was further validated with Markov state modeling with independent MD simulations. Therefore, our work provides strong evidence for the conformational plasticity of protein, which is mainly due to the inherent degree of flexibility. As a reliable and efficient enhanced sampling protocol, teDA2 could be used to study the dynamics between functional states of various biomolecular machines.
Collapse
Affiliation(s)
- Yigao Yuan
- Kuang Yaming Honors School, Nanjing University, 210023 Nanjing, China
| | - Qiang Zhu
- Kuang Yaming Honors School, Nanjing University, 210023 Nanjing, China.,Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, 210023 Nanjing China
| | - Ruiheng Song
- Kuang Yaming Honors School, Nanjing University, 210023 Nanjing, China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, 210023 Nanjing China
| | - Hao Dong
- Kuang Yaming Honors School, Nanjing University, 210023 Nanjing, China.,Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| |
Collapse
|
22
|
Wang J, Peng C, Yu Y, Chen Z, Xu Z, Cai T, Shao Q, Shi J, Zhu W. Exploring Conformational Change of Adenylate Kinase by Replica Exchange Molecular Dynamic Simulation. Biophys J 2020; 118:1009-1018. [PMID: 31995738 DOI: 10.1016/j.bpj.2020.01.001] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Revised: 12/28/2019] [Accepted: 01/03/2020] [Indexed: 12/19/2022] Open
Abstract
Replica exchange molecular dynamics (REMD) simulation is a popular enhanced sampling method that is widely used for exploring the atomic mechanism of protein conformational change. However, the requirement of huge computational resources for REMD, especially with the explicit solvent model, largely limits its application. In this study, the availability and efficiency of a variant of velocity-scaling REMD (vsREMD) was assessed with adenylate kinase as an example. Although vsREMD achieved results consistent with those from conventional REMD and experimental studies, the number of replicas required for vsREMD (30) was much less than that for conventional REMD (80) to achieve a similar acceptance rate (∼0.2), demonstrating high efficiency of vsREMD to characterize the protein conformational change and associated free-energy profile. Thus, vsREMD is a highly efficient approach for studying the large-scale conformational change of protein systems.
Collapse
Affiliation(s)
- Jinan Wang
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.
| | - Cheng Peng
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China; University of Chinese Academy of Sciences, Beijing, China
| | - Yuqu Yu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Zhaoqiang Chen
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Zhijian Xu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China; University of Chinese Academy of Sciences, Beijing, China
| | - Tingting Cai
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Qiang Shao
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China; University of Chinese Academy of Sciences, Beijing, China
| | - Jiye Shi
- UCB Biopharma SPRL, Braine-l'Alleud, Belgium
| | - Weiliang Zhu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China; Open Studio for Druggability Research of Marine Lead Compounds, Qingdao National Laboratory for Marine Science and Technology, Jimo, Qingdao, China; University of Chinese Academy of Sciences, Beijing, China.
| |
Collapse
|
23
|
Rogne P, Andersson D, Grundström C, Sauer-Eriksson E, Linusson A, Wolf-Watz M. Nucleation of an Activating Conformational Change by a Cation-π Interaction. Biochemistry 2019; 58:3408-3412. [PMID: 31339702 DOI: 10.1021/acs.biochem.9b00538] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
As a key molecule in biology, adenosine triphosphate (ATP) has numerous crucial functions in, for instance, energetics, post-translational modifications, nucleotide biosynthesis, and cofactor metabolism. Here, we have discovered an intricate interplay between the enzyme adenylate kinase and its substrate ATP. The side chain of an arginine residue was found to be an efficient sensor of the aromatic moiety of ATP through the formation of a strong cation-π interaction. In addition to recognition, the interaction was found to have dual functionality. First, it nucleates the activating conformational transition of the ATP binding domain and also affects the specificity in the distant AMP binding domain. In light of the functional consequences resulting from the cation-π interaction, it is possible that the mode of ATP recognition may be a useful tool in enzyme design.
Collapse
Affiliation(s)
- Per Rogne
- Department of Chemistry , Umeå University , SE-901 87 Umeå , Sweden
| | - David Andersson
- Department of Chemistry , Umeå University , SE-901 87 Umeå , Sweden
| | | | | | - Anna Linusson
- Department of Chemistry , Umeå University , SE-901 87 Umeå , Sweden
| | - Magnus Wolf-Watz
- Department of Chemistry , Umeå University , SE-901 87 Umeå , Sweden
| |
Collapse
|
24
|
Inomata K, Kamoshida H, Ikari M, Ito Y, Kigawa T. Impact of cellular health conditions on the protein folding state in mammalian cells. Chem Commun (Camb) 2018; 53:11245-11248. [PMID: 28960222 DOI: 10.1039/c7cc06004a] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
By using in-cell NMR experiments, we have demonstrated that the protein folding state in cells is significantly influenced by the cellular health conditions. hAK1 was denatured in cells under stressful culture conditions, while it remained functional and properly folded in cells continuously supplied with a fresh medium.
Collapse
Affiliation(s)
- Kohsuke Inomata
- RIKEN Quantitative Biology Center, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan. and PRESTO/Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Hajime Kamoshida
- Department of Chemistry, Graduate School of Science and Engineering, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397, Japan
| | - Masaomi Ikari
- RIKEN Quantitative Biology Center, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.
| | - Yutaka Ito
- Department of Chemistry, Graduate School of Science and Engineering, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji, Tokyo 192-0397, Japan
| | - Takanori Kigawa
- RIKEN Quantitative Biology Center, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan. and Department of Computer Science, School of Computing, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8503, Japan
| |
Collapse
|
25
|
Direct observation of ultrafast large-scale dynamics of an enzyme under turnover conditions. Proc Natl Acad Sci U S A 2018. [PMID: 29531052 PMCID: PMC5879700 DOI: 10.1073/pnas.1720448115] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The potential effect of conformational dynamics of enzymes on their chemical steps has been intensely debated recently. We use single-molecule FRET experiments on adenylate kinase (AK) to shed new light on this question. AK closes its domains to bring its two substrate close together for reaction. We show that domain closure takes only microseconds to complete, which is two orders of magnitude faster than the chemical reaction. Nevertheless, active-site mutants that reduce the rate of domain closure also reduce the reaction rate, suggesting a connection between the two phenomena. We propose that ultrafast domain closure is used by enzymes as a mechanism to optimize mutual orientation of substrates, a novel mode of coupling between conformational dynamics and catalysis. The functional cycle of many proteins involves large-scale motions of domains and subunits. The relation between conformational dynamics and the chemical steps of enzymes remains under debate. Here we show that in the presence of substrates, domain motions of an enzyme can take place on the microsecond time scale, yet exert influence on the much-slower chemical step. We study the domain closure reaction of the enzyme adenylate kinase from Escherichia coli while in action (i.e., under turnover conditions), using single-molecule FRET spectroscopy. We find that substrate binding increases dramatically domain closing and opening times, making them as short as ∼15 and ∼45 µs, respectively. These large-scale conformational dynamics are likely the fastest measured to date, and are ∼100–200 times faster than the enzymatic turnover rate. Some active-site mutants are shown to fully or partially prevent the substrate-induced increase in domain closure times, while at the same time they also reduce enzymatic activity, establishing a clear connection between the two phenomena, despite their disparate time scales. Based on these surprising observations, we propose a paradigm for the mode of action of enzymes, in which numerous cycles of conformational rearrangement are required to find a mutual orientation of substrates that is optimal for the chemical reaction.
Collapse
|
26
|
Molecular mechanism of ATP versus GTP selectivity of adenylate kinase. Proc Natl Acad Sci U S A 2018; 115:3012-3017. [PMID: 29507216 DOI: 10.1073/pnas.1721508115] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Enzymatic substrate selectivity is critical for the precise control of metabolic pathways. In cases where chemically related substrates are present inside cells, robust mechanisms of substrate selectivity are required. Here, we report the mechanism utilized for catalytic ATP versus GTP selectivity during adenylate kinase (Adk) -mediated phosphorylation of AMP. Using NMR spectroscopy we found that while Adk adopts a catalytically competent and closed structural state in complex with ATP, the enzyme is arrested in a catalytically inhibited and open state in complex with GTP. X-ray crystallography experiments revealed that the interaction interfaces supporting ATP and GTP recognition, in part, are mediated by coinciding residues. The mechanism provides an atomic view on how the cellular GTP pool is protected from Adk turnover, which is important because GTP has many specialized cellular functions. In further support of this mechanism, a structure-function analysis enabled by synthesis of ATP analogs suggests that a hydrogen bond between the adenine moiety and the backbone of the enzyme is vital for ATP selectivity. The importance of the hydrogen bond for substrate selectivity is likely general given the conservation of its location and orientation across the family of eukaryotic protein kinases.
Collapse
|
27
|
Mazal H, Aviram H, Riven I, Haran G. Effect of ligand binding on a protein with a complex folding landscape. Phys Chem Chem Phys 2018; 20:3054-3062. [PMID: 28721412 DOI: 10.1039/c7cp03327c] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Ligand binding to a protein can stabilize it significantly against unfolding. The variation of the folding free energy, ΔΔG0, due to ligand binding can be derived from a simple reaction scheme involving exclusive binding to the native state. One obtains the following expression: , where Kd is the ligand dissociation constant and L is its concentration, R is the universal gas constant and T is the temperature. This expression has been shown to correctly describe experimental results on multiple proteins. In the current work we studied the effect of ligand binding on the stability of the multi-domain protein adenylate kinase from E. coli (AKE). Unfolding experiments were conducted using single-molecule FRET spectroscopy, which allowed us to directly obtain the fraction of unfolded protein in a model-free way from FRET efficiency histograms. Surprisingly, it was found that the effect of two inhibitors (Ap5A and AMPPNP) and a substrate (AMP) on the stability of AKE was much smaller than expected based on Kd values obtained independently using microscale thermophoresis. To shed light on this issue, we measured the Kd for Ap5A over a range of chemical denaturant concentrations where the protein is still folded. It was found that Kd increases dramatically over this range, likely due to the population of folding intermediates, whose binding to the ligand is much weaker than that of the native state. We propose that binding to folding intermediates may dominate the effect of ligands on the stability of multi-domain proteins, and could therefore have a strong impact on protein homeostasis in vivo.
Collapse
Affiliation(s)
- Hisham Mazal
- Department of Chemical Physics, Weizmann Institute of Science, Rehovot 7610001, Israel.
| | | | | | | |
Collapse
|
28
|
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: 24] [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
|
29
|
Boswell ZK, Rahman S, Canny MD, Latham MP. A dynamic allosteric pathway underlies Rad50 ABC ATPase function in DNA repair. Sci Rep 2018; 8:1639. [PMID: 29374232 PMCID: PMC5786021 DOI: 10.1038/s41598-018-19908-8] [Citation(s) in RCA: 14] [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: 11/22/2017] [Accepted: 01/10/2018] [Indexed: 11/23/2022] Open
Abstract
The Mre11-Rad50 protein complex is an initial responder to sites of DNA double strand breaks. Many studies have shown that ATP binding to Rad50 causes global changes to the Mre11-Rad50 structure, which are important for DNA repair functions. Here we used methyl-based NMR spectroscopy on a series of mutants to describe a dynamic allosteric pathway within Rad50. Mutations result in changes in the side chain methyl group chemical environment that are correlated with altered nanosecond timescale dynamics. We also observe striking relationships between the magnitude of chemical shift perturbations and Rad50 and Mre11 activities. Together, these data suggest an equilibrium between a ground state and an "active" dimerization competent state of Rad50 that has locally altered structure and dynamics and is poised for ATP-induced dimerization and eventual ATP hydrolysis. Thus, this sparsely populated intermediate is critical for Mre11-Rad50-directed DNA double strand break repair.
Collapse
Affiliation(s)
- Zachary K Boswell
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX, 79423, USA
| | - Samiur Rahman
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX, 79423, USA
| | - Marella D Canny
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX, 79423, USA
| | - Michael P Latham
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, TX, 79423, USA.
| |
Collapse
|
30
|
Mehaffey MR, Cammarata MB, Brodbelt JS. Tracking the Catalytic Cycle of Adenylate Kinase by Ultraviolet Photodissociation Mass Spectrometry. Anal Chem 2018; 90:839-846. [PMID: 29188992 PMCID: PMC5750083 DOI: 10.1021/acs.analchem.7b03591] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The complex interplay of dynamic protein plasticity and specific side-chain interactions with substrate molecules that allows enzymes to catalyze reactions has yet to be fully unraveled. Top-down ultraviolet photodissociation (UVPD) mass spectrometry is used to track snapshots of conformational fluctuations in the phosphotransferase adenylate kinase (AK) throughout its active reaction cycle by characterization of complexes containing AK and each of four different adenosine phosphate ligands. Variations in efficiencies of UVPD backbone cleavages were consistently observed for three α-helices and the adenosine binding regions for AK complexes representing different steps of the catalytic cycle, implying that these stretches of the protein sample various structural microstates as the enzyme undergoes global open-to-closed transitions. Focusing on the conformational impact of recruiting or releasing the Mg2+ cofactor highlights two loop regions for which fragmentation increases upon UVPD, signaling an increase in loop flexibility as the metal cation disrupts the loop interactions with the substrate ligands. Additionally, the observation of holo ions and variations in UVPD backbone cleavage efficiency at R138 implicate this conserved active site residue in stabilizing the donor phosphoryl group during catalysis. This study showcases the utility of UVPD-MS to provide insight into conformational fluctuations of single residues for active enzymes.
Collapse
Affiliation(s)
- M. Rachel Mehaffey
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712
| | | | | |
Collapse
|
31
|
Kar P, Feig M. Hybrid All-Atom/Coarse-Grained Simulations of Proteins by Direct Coupling of CHARMM and PRIMO Force Fields. J Chem Theory Comput 2017; 13:5753-5765. [PMID: 28992696 DOI: 10.1021/acs.jctc.7b00840] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Hybrid all-atom/coarse-grained (AA/CG) simulations of proteins offer a computationally efficient compromise where atomistic details are only applied to biologically relevant regions while benefiting from the speedup of treating the remaining parts of a given system at the CG level. The recently developed CG model, PRIMO, allows a direct coupling with an atomistic force field with no additional modifications or coupling terms and the ability to carry out dynamic simulations without any restraints on secondary or tertiary structures. A hybrid AA/CG scheme based on combining all-atom CHARMM and coarse-grained PRIMO representations was validated via molecular dynamics and replica exchange simulations of soluble and membrane proteins. The AA/CG scheme was also tested in the calculation of the free energy profile for the transition from the closed to the open state of adenylate kinase via umbrella sampling molecular dynamics method. The overall finding is that the AA/CG scheme generates dynamics and energetics that are qualitatively and quantitatively comparable to AA simulations while offering the computational advantages of coarse-graining. This model opens the door to challenging applications where high accuracy is required only in parts of large biomolecular complexes.
Collapse
Affiliation(s)
- Parimal Kar
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|
32
|
Onuk E, Badger J, Wang YJ, Bardhan J, Chishti Y, Akcakaya M, Brooks DH, Erdogmus D, Minh DDL, Makowski L. Effects of Catalytic Action and Ligand Binding on Conformational Ensembles of Adenylate Kinase. Biochemistry 2017; 56:4559-4567. [PMID: 28767234 DOI: 10.1021/acs.biochem.7b00351] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Crystal structures of adenylate kinase (AdK) from Escherichia coli capture two states: an "open" conformation (apo) obtained in the absence of ligands and a "closed" conformation in which ligands are bound. Other AdK crystal structures suggest intermediate conformations that may lie on the transition pathway between these two states. To characterize the transition from open to closed states in solution, X-ray solution scattering data were collected from AdK in the apo form and with progressively increasing concentrations of five different ligands. Scattering data from apo AdK are consistent with scattering predicted from the crystal structure of AdK in the open conformation. In contrast, data from AdK samples saturated with Ap5A do not agree with that calculated from AdK in the closed conformation. Using cluster analysis of available structures, we selected representative structures in five conformational states: open, partially open, intermediate, partially closed, and closed. We used these structures to estimate the relative abundances of these states for each experimental condition. X-ray solution scattering data obtained from AdK with AMP are dominated by scattering from AdK in the open conformation. For AdK in the presence of high concentrations of ATP and ADP, the conformational ensemble shifts to a mixture of partially open and closed states. Even when AdK is saturated with Ap5A, a significant proportion of AdK remains in a partially open conformation. These results are consistent with an induced-fit model in which the transition of AdK from an open state to a closed state is initiated by ATP binding.
Collapse
Affiliation(s)
- Emre Onuk
- Radiation Oncology Department, University of California , Los Angeles, California 90095, United States
| | - John Badger
- DeltaG Technologies , San Diego, California 92122, United States
| | - Yu Jing Wang
- Department of Bioengineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Jaydeep Bardhan
- Department of Mechanical and Industrial Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Yasmin Chishti
- Department of Bioengineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Murat Akcakaya
- Department of Electrical and Computer Engineering, University of Pittsburgh , Pittsburgh, Pennsylvania 15261, United States
| | - Dana H Brooks
- Department of Electrical and Computer Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - Deniz Erdogmus
- Department of Electrical and Computer Engineering, Northeastern University , Boston, Massachusetts 02115, United States
| | - David D L Minh
- Department of Chemistry, Illinois Institute of Technology , Chicago, Illinois 60616, United States
| | - Lee Makowski
- Department of Bioengineering, Northeastern University , Boston, Massachusetts 02115, United States.,Department of Chemistry and Chemical Biology, Northeastern University , Boston, Massachusetts 02115, United States
| |
Collapse
|
33
|
Rogne P, Wolf-Watz M. Urea-Dependent Adenylate Kinase Activation following Redistribution of Structural States. Biophys J 2017; 111:1385-1395. [PMID: 27705762 DOI: 10.1016/j.bpj.2016.08.028] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Revised: 08/22/2016] [Accepted: 08/26/2016] [Indexed: 11/29/2022] Open
Abstract
Proteins are often functionally dependent on conformational changes that allow them to sample structural states that are sparsely populated in the absence of a substrate or binding partner. The distribution of such structural microstates is governed by their relative stability, and the kinetics of their interconversion is governed by the magnitude of associated activation barriers. Here, we have explored the interplay among structure, stability, and function of a selected enzyme, adenylate kinase (Adk), by monitoring changes in its enzymatic activity in response to additions of urea. For this purpose we used a 31P NMR assay that was found useful for heterogeneous sample compositions such as presence of urea. It was found that Adk is activated at low urea concentrations whereas higher urea concentrations unfolds and thereby deactivates the enzyme. From a quantitative analysis of chemical shifts, it was found that urea redistributes preexisting structural microstates, stabilizing a substrate-bound open state at the expense of a substrate-bound closed state. Adk is rate-limited by slow opening of substrate binding domains and the urea-dependent redistribution of structural states is consistent with a model where the increased activity results from an increased rate-constant for domain opening. In addition, we also detected a strong correlation between the catalytic free energy and free energy of substrate (ATP) binding, which is also consistent with the catalytic model for Adk. From a general perspective, it appears that urea can be used to modulate conformational equilibria of folded proteins toward more expanded states for cases where a sizeable difference in solvent-accessible surface area exists between the states involved. This effect complements the action of osmolytes, such as trimethylamine N-oxide, that favor more compact protein states.
Collapse
Affiliation(s)
- Per Rogne
- Department of Chemistry, Umeå University, Umeå, Sweden
| | | |
Collapse
|
34
|
Halder R, Manna RN, Chakraborty S, Jana B. Modulation of the Conformational Dynamics of Apo-Adenylate Kinase through a π–Cation Interaction. J Phys Chem B 2017; 121:5699-5708. [DOI: 10.1021/acs.jpcb.7b01736] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Ritaban Halder
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Rabindra Nath Manna
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Sandipan Chakraborty
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Biman Jana
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
35
|
Ozgur B, Ozdemir ES, Gursoy A, Keskin O. Relation between Protein Intrinsic Normal Mode Weights and Pre-Existing Conformer Populations. J Phys Chem B 2017; 121:3686-3700. [DOI: 10.1021/acs.jpcb.6b10401] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Beytullah Ozgur
- Center for Computational Biology and Bioinformatics, ‡Chemical and Biological
Engineering, and §Computer Engineering,
College of Engineering, Koc University, 34450 Istanbul, Turkey
| | - E. Sila Ozdemir
- Center for Computational Biology and Bioinformatics, ‡Chemical and Biological
Engineering, and §Computer Engineering,
College of Engineering, Koc University, 34450 Istanbul, Turkey
| | - Attila Gursoy
- Center for Computational Biology and Bioinformatics, ‡Chemical and Biological
Engineering, and §Computer Engineering,
College of Engineering, Koc University, 34450 Istanbul, Turkey
| | - Ozlem Keskin
- Center for Computational Biology and Bioinformatics, ‡Chemical and Biological
Engineering, and §Computer Engineering,
College of Engineering, Koc University, 34450 Istanbul, Turkey
| |
Collapse
|
36
|
Abstract
It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.
Collapse
|
37
|
Tükenmez H, Magnussen HM, Kovermann M, Byström A, Wolf-Watz M. Linkage between Fitness of Yeast Cells and Adenylate Kinase Catalysis. PLoS One 2016; 11:e0163115. [PMID: 27642758 PMCID: PMC5028032 DOI: 10.1371/journal.pone.0163115] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Accepted: 09/03/2016] [Indexed: 01/08/2023] Open
Abstract
Enzymes have evolved with highly specific values of their catalytic parameters kcat and KM. This poses fundamental biological questions about the selection pressures responsible for evolutionary tuning of these parameters. Here we are address these questions for the enzyme adenylate kinase (Adk) in eukaryotic yeast cells. A plasmid shuffling system was developed to allow quantification of relative fitness (calculated from growth rates) of yeast in response to perturbations of Adk activity introduced through mutations. Biophysical characterization verified that all variants studied were properly folded and that the mutations did not cause any substantial differences to thermal stability. We found that cytosolic Adk is essential for yeast viability in our strain background and that viability could not be restored with a catalytically dead, although properly folded Adk variant. There exist a massive overcapacity of Adk catalytic activity and only 12% of the wild type kcat is required for optimal growth at the stress condition 20°C. In summary, the approach developed here has provided new insights into the evolutionary tuning of kcat for Adk in a eukaryotic organism. The developed methodology may also become useful for uncovering new aspects of active site dynamics and also in enzyme design since a large library of enzyme variants can be screened rapidly by identifying viable colonies.
Collapse
Affiliation(s)
- Hasan Tükenmez
- Department of Molecular Biology, Umeå University, SE-901 87, Umeå, Sweden
| | | | | | - Anders Byström
- Department of Molecular Biology, Umeå University, SE-901 87, Umeå, Sweden
- * E-mail: (MWW); (AB)
| | - Magnus Wolf-Watz
- Department of Chemistry, Umeå University, SE-901 87, Umeå, Sweden
- * E-mail: (MWW); (AB)
| |
Collapse
|
38
|
Whitney DS, Volkman BF, Prehoda KE. Evolution of a Protein Interaction Domain Family by Tuning Conformational Flexibility. J Am Chem Soc 2016; 138:15150-15156. [PMID: 27502157 DOI: 10.1021/jacs.6b05954] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Conformational flexibility allows proteins to adopt multiple functionally important conformations but can also lead to nonfunctional structures. We analyzed the dynamic behavior of the enzyme guanylate kinase as it evolved into the GK protein interaction domain (GKPID) to investigate the role of flexibility in the evolution of new protein functions. We found that the ancestral enzyme is very flexible, allowing it to adopt open conformations that can bind nucleotide and closed ones that enable catalysis of phosphotransfer from ATP to GMP. Historical mutations that converted the GK from an enzyme to a protein interaction domain dramatically reduce flexibility, predominantly by inhibiting rotations of the protein backbone that are coupled to the global closing motion. Removing flexibility prevents adoption of conformations that cannot fit the protein partner in the binding site. Our results highlight the importance of mutations that optimize protein conformational flexibility with function during evolution.
Collapse
Affiliation(s)
- Dustin S Whitney
- Department of Biochemistry, Medical College of Wisconsin , Milwaukee, Wisconsin 53226, United States
| | - Brian F Volkman
- Department of Biochemistry, Medical College of Wisconsin , Milwaukee, Wisconsin 53226, United States
| | - Kenneth E Prehoda
- Department of Chemistry and Biochemistry, Institute of Molecular Biology, University of Oregon , Eugene, Oregon 97403, United States
| |
Collapse
|
39
|
Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase. Biophys J 2016; 109:1978-85. [PMID: 26536274 DOI: 10.1016/j.bpj.2015.08.049] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2015] [Accepted: 08/27/2015] [Indexed: 11/21/2022] Open
Abstract
The enzyme adenylate kinase (ADK) features two substrate binding domains that undergo large-scale motions during catalysis. In the apo state, the enzyme preferentially adopts a globally open state with accessible binding sites. Binding of two substrate molecules (AMP + ATP or ADP + ADP) results in a closed domain conformation, allowing efficient phosphoryl-transfer catalysis. We employed molecular dynamics simulations to systematically investigate how the individual domain motions are modulated by the binding of substrates. Two-dimensional free-energy landscapes were calculated along the opening of the two flexible lid domains for apo and holo ADK as well as for all single natural substrates bound to one of the two binding sites of ADK. The simulations reveal a strong dependence of the conformational ensembles on type and binding position of the bound substrates and a nonsymmetric behavior of the lid domains. Altogether, the ensembles suggest that, upon initial substrate binding to the corresponding lid site, the opposing lid is maintained open and accessible for subsequent substrate binding. In contrast, ATP binding to the AMP-lid induces global domain closing, preventing further substrate binding to the ATP-lid site. This might constitute a mechanism by which the enzyme avoids the formation of a stable but enzymatically unproductive state.
Collapse
|
40
|
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: 6.0] [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
|
41
|
Maximova T, Moffatt R, Ma B, Nussinov R, Shehu A. Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics. PLoS Comput Biol 2016; 12:e1004619. [PMID: 27124275 PMCID: PMC4849799 DOI: 10.1371/journal.pcbi.1004619] [Citation(s) in RCA: 132] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.
Collapse
Affiliation(s)
- Tatiana Maximova
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Ryan Moffatt
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
- Department of Biongineering, George Mason University, Fairfax, Virginia, United States of America
- School of Systems Biology, George Mason University, Manassas, Virginia, United States of America
| |
Collapse
|
42
|
Subnanometre enzyme mechanics probed by single-molecule force spectroscopy. Nat Commun 2016; 7:10848. [PMID: 26906294 PMCID: PMC4770092 DOI: 10.1038/ncomms10848] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Accepted: 01/26/2016] [Indexed: 01/01/2023] Open
Abstract
Enzymes are molecular machines that bind substrates specifically, provide an adequate chemical environment for catalysis and exchange products rapidly, to ensure fast turnover rates. Direct information about the energetics that drive conformational changes is difficult to obtain. We used subnanometre single-molecule force spectroscopy to study the energetic drive of substrate-dependent lid closing in the enzyme adenylate kinase. Here we show that in the presence of the bisubstrate inhibitor diadenosine pentaphosphate (AP5A), closing and opening of both lids is cooperative and tightly coupled to inhibitor binding. Surprisingly, binding of the substrates ADP and ATP exhibits a much smaller energetic drive towards the fully closed state. Instead, we observe a new dominant energetic minimum with both lids half closed. Our results, combining experiment and molecular dynamics simulations, give detailed mechanical insights into how an enzyme can cope with the seemingly contradictory requirements of rapid substrate exchange and tight closing, to ensure efficient catalysis. Adenylate kinase catalyses the interconversion of adenosine phosphates, and plays a crucial role in maintaining cellular energy homeostasis. Here, the authors use single molecule optical tweezers to understand how the enzyme's conformation dynamics modulates catalysis.
Collapse
|
43
|
Shao Q. Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions. Phys Chem Chem Phys 2016; 18:29170-29182. [DOI: 10.1039/c6cp05634b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.
Collapse
Affiliation(s)
- Qiang Shao
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| |
Collapse
|
44
|
Rogne P, Sparrman T, Anugwom I, Mikkola JP, Wolf-Watz M. Realtime (31)P NMR Investigation on the Catalytic Behavior of the Enzyme Adenylate kinase in the Matrix of a Switchable Ionic Liquid. CHEMSUSCHEM 2015; 8:3764-3768. [PMID: 26494201 DOI: 10.1002/cssc.201501104] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Indexed: 06/05/2023]
Abstract
The integration of highly efficient enzymatic catalysis with the solvation properties of ionic liquids for an environmentally friendly and efficient use of raw materials such as wood requires fundamental knowledge about the influence of relevant ionic liquids on enzymes. Switchable ionic liquids (SIL) are promising candidates for implementation of enzymatic treatments of raw materials. One industrially interesting SIL is constituted by monoethanol amine (MEA) and 1,8-diazabicyclo-[5.4.0]-undec-7-ene (DBU) formed with sulfur dioxide (SO2) as the coupling media (DBU-SO2-MEASIL). It has the ability to solubilize the matrix of lignocellulosic biomass while leaving the cellulose backbone intact. Using a novel (31)P NMR-based real-time assay we show that this SIL is compatible with enzymatic catalysis because a model enzyme, adenylate kinase, retains its activity in up to at least 25 wt % of DBU-SO2-MEASIL. Thus this SIL appears suitable for, for example, enzymatic degradation of hemicellulose.
Collapse
Affiliation(s)
- Per Rogne
- Department of Chemistry, Umeå University, Linnaeus väg 10, 901 87, Umeå, Sweden
| | - Tobias Sparrman
- Department of Chemistry, Umeå University, Linnaeus väg 10, 901 87, Umeå, Sweden
| | - Ikenna Anugwom
- Department of Chemistry, Umeå University, Linnaeus väg 10, 901 87, Umeå, Sweden
| | - Jyri-Pekka Mikkola
- Department of Chemistry, Umeå University, Linnaeus väg 10, 901 87, Umeå, Sweden
- Laboratory of Industrial Chemistry and Reaction Engineering, Johan Gadolin Process Chemistry Centre, Åbo Akademi University, Åbo-Turku, 20500, Finland
| | - Magnus Wolf-Watz
- Department of Chemistry, Umeå University, Linnaeus väg 10, 901 87, Umeå, Sweden.
| |
Collapse
|
45
|
Seyler SL, Kumar A, Thorpe MF, Beckstein O. Path Similarity Analysis: A Method for Quantifying Macromolecular Pathways. PLoS Comput Biol 2015; 11:e1004568. [PMID: 26488417 PMCID: PMC4619321 DOI: 10.1371/journal.pcbi.1004568] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2015] [Accepted: 09/23/2015] [Indexed: 01/03/2023] Open
Abstract
Diverse classes of proteins function through large-scale conformational changes and various sophisticated computational algorithms have been proposed to enhance sampling of these macromolecular transition paths. Because such paths are curves in a high-dimensional space, it has been difficult to quantitatively compare multiple paths, a necessary prerequisite to, for instance, assess the quality of different algorithms. We introduce a method named Path Similarity Analysis (PSA) that enables us to quantify the similarity between two arbitrary paths and extract the atomic-scale determinants responsible for their differences. PSA utilizes the full information available in 3N-dimensional configuration space trajectories by employing the Hausdorff or Fréchet metrics (adopted from computational geometry) to quantify the degree of similarity between piecewise-linear curves. It thus completely avoids relying on projections into low dimensional spaces, as used in traditional approaches. To elucidate the principles of PSA, we quantified the effect of path roughness induced by thermal fluctuations using a toy model system. Using, as an example, the closed-to-open transitions of the enzyme adenylate kinase (AdK) in its substrate-free form, we compared a range of protein transition path-generating algorithms. Molecular dynamics-based dynamic importance sampling (DIMS) MD and targeted MD (TMD) and the purely geometric FRODA (Framework Rigidity Optimized Dynamics Algorithm) were tested along with seven other methods publicly available on servers, including several based on the popular elastic network model (ENM). PSA with clustering revealed that paths produced by a given method are more similar to each other than to those from another method and, for instance, that the ENM-based methods produced relatively similar paths. PSA was applied to ensembles of DIMS MD and FRODA trajectories of the conformational transition of diphtheria toxin, a particularly challenging example. For the AdK transition, the new concept of a Hausdorff-pair map enabled us to extract the molecular structural determinants responsible for differences in pathways, namely a set of conserved salt bridges whose charge-charge interactions are fully modelled in DIMS MD but not in FRODA. PSA has the potential to enhance our understanding of transition path sampling methods, validate them, and to provide a new approach to analyzing conformational transitions. Many proteins are nanomachines that perform mechanical or chemical work by changing their three-dimensional shape and cycle between multiple conformational states. Computer simulations of such conformational transitions provide mechanistic insights into protein function but such simulations have been challenging. In particular, it is not clear how to quantitatively compare current simulation methods or to assess their accuracy. To that end, we present a general and flexible computational framework for quantifying transition paths—by measuring mutual geometric similarity—that, compared with existing approaches, requires minimal a-priori assumptions and can take advantage of full atomic detail alongside heuristic information derived from intuition. Using our Path Similarity Analysis (PSA) framework in parallel with several existing quantitative approaches, we examine transitions generated for a toy model of a transition and two biological systems, the enzyme adenylate kinase and diphtheria toxin. Our results show that PSA enables the quantitative comparison of different path sampling methods and aids the identification of potentially important atomistic motions by exploiting geometric information in transition paths. The method has the potential to enhance our understanding of transition path sampling methods, validate them, and to provide a new approach to analyzing macromolecular conformational transitions.
Collapse
Affiliation(s)
- Sean L. Seyler
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| | - Avishek Kumar
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
| | - M. F. Thorpe
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, United Kingdom
| | - Oliver Beckstein
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, Arizona, United States of America
- * E-mail:
| |
Collapse
|
46
|
Onuk AE, Akcakaya M, Bardhan JP, Erdogmus D, Brooks DH, Makowski L. Constrained Maximum Likelihood Estimation of Relative Abundances of Protein Conformation in a Heterogeneous Mixture from Small Angle X-Ray Scattering Intensity Measurements. IEEE TRANSACTIONS ON SIGNAL PROCESSING : A PUBLICATION OF THE IEEE SIGNAL PROCESSING SOCIETY 2015; 63:5383-5394. [PMID: 26924916 PMCID: PMC4767180 DOI: 10.1109/tsp.2015.2455515] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
In this paper, we describe a model for maximum likelihood estimation (MLE) of the relative abundances of different conformations of a protein in a heterogeneous mixture from small angle X-ray scattering (SAXS) intensities. To consider cases where the solution includes intermediate or unknown conformations, we develop a subset selection method based on k-means clustering and the Cramér-Rao bound on the mixture coefficient estimation error to find a sparse basis set that represents the space spanned by the measured SAXS intensities of the known conformations of a protein. Then, using the selected basis set and the assumptions on the model for the intensity measurements, we show that the MLE model can be expressed as a constrained convex optimization problem. Employing the adenylate kinase (ADK) protein and its known conformations as an example, and using Monte Carlo simulations, we demonstrate the performance of the proposed estimation scheme. Here, although we use 45 crystallographically determined experimental structures and we could generate many more using, for instance, molecular dynamics calculations, the clustering technique indicates that the data cannot support the determination of relative abundances for more than 5 conformations. The estimation of this maximum number of conformations is intrinsic to the methodology we have used here.
Collapse
Affiliation(s)
- A Emre Onuk
- Electrical and Computer Engineering Department, Northeastern University, Boston, MA
| | - Murat Akcakaya
- Electrical and Computer Engineering Department, University of Pittsburgh, Pittsburgh, PA
| | - Jaydeep P Bardhan
- Department of Mechanical and Industrial Engineering, Northeastern University, Boston, MA
| | - Deniz Erdogmus
- Electrical and Computer Engineering Department, Northeastern University, Boston, MA
| | - Dana H Brooks
- Electrical and Computer Engineering Department, Northeastern University, Boston, MA
| | - Lee Makowski
- Bioengineering Department, Northeastern University, Boston, MA; Department of Chemistry and Chemical Biology, Northeastern University, Boston, MA
| |
Collapse
|
47
|
Unan H, Yildirim A, Tekpinar M. Opening mechanism of adenylate kinase can vary according to selected molecular dynamics force field. J Comput Aided Mol Des 2015; 29:655-65. [DOI: 10.1007/s10822-015-9849-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2015] [Accepted: 05/22/2015] [Indexed: 11/29/2022]
|
48
|
Formoso E, Limongelli V, Parrinello M. Energetics and structural characterization of the large-scale functional motion of adenylate kinase. Sci Rep 2015; 5:8425. [PMID: 25672826 PMCID: PMC4325324 DOI: 10.1038/srep08425] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2014] [Accepted: 01/08/2015] [Indexed: 12/22/2022] Open
Abstract
Adenylate Kinase (AK) is a signal transducing protein that regulates cellular energy homeostasis balancing between different conformations. An alteration of its activity can lead to severe pathologies such as heart failure, cancer and neurodegenerative diseases. A comprehensive elucidation of the large-scale conformational motions that rule the functional mechanism of this enzyme is of great value to guide rationally the development of new medications. Here using a metadynamics-based computational protocol we elucidate the thermodynamics and structural properties underlying the AK functional transitions. The free energy estimation of the conformational motions of the enzyme allows characterizing the sequence of events that regulate its action. We reveal the atomistic details of the most relevant enzyme states, identifying residues such as Arg119 and Lys13, which play a key role during the conformational transitions and represent druggable spots to design enzyme inhibitors. Our study offers tools that open new areas of investigation on large-scale motion in proteins.
Collapse
Affiliation(s)
- Elena Formoso
- 1] Department of Chemistry and Applied Biosciences, ETH Zurich, and Faculty of Informatics, Institute of Computational Science, Università della Svizzera Italiana, via G. Buffi 13, CH-6900 Lugano, Switzerland [2] Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU) and Donostia International Physics Center (DIPC), PK 1072, 20080 Donostia, Euskadi, Spain
| | - Vittorio Limongelli
- 1] Università della Svizzera Italiana (USI), Faculty of Informatics, Institute of Computational Science, via G. Buffi 13, CH-6900 Lugano, Switzerland [2] Department of Pharmacy, University of Naples "Federico II", via D. Montesano 49, I-80131 Naples, Italy
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, and Faculty of Informatics, Institute of Computational Science, Università della Svizzera Italiana, via G. Buffi 13, CH-6900 Lugano, Switzerland
| |
Collapse
|
49
|
Esteban-Martín S, Fenwick RB, Ådén J, Cossins B, Bertoncini CW, Guallar V, Wolf-Watz M, Salvatella X. Correlated inter-domain motions in adenylate kinase. PLoS Comput Biol 2014; 10:e1003721. [PMID: 25078441 PMCID: PMC4117416 DOI: 10.1371/journal.pcbi.1003721] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2014] [Accepted: 05/28/2014] [Indexed: 11/18/2022] Open
Abstract
Correlated inter-domain motions in proteins can mediate fundamental biochemical processes such as signal transduction and allostery. Here we characterize at structural level the inter-domain coupling in a multidomain enzyme, Adenylate Kinase (AK), using computational methods that exploit the shape information encoded in residual dipolar couplings (RDCs) measured under steric alignment by nuclear magnetic resonance (NMR). We find experimental evidence for a multi-state equilibrium distribution along the opening/closing pathway of Adenylate Kinase, previously proposed from computational work, in which inter-domain interactions disfavour states where only the AMP binding domain is closed. In summary, we provide a robust experimental technique for study of allosteric regulation in AK and other enzymes.
Collapse
Affiliation(s)
- Santiago Esteban-Martín
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Barcelona Supercomputing Center - BSC, Barcelona, Spain
- * E-mail: (SEM); (XS)
| | - Robert Bryn Fenwick
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine – IRB Barcelona, Barcelona, Spain
| | - Jörgen Ådén
- Department of Chemistry, Chemical Biological Centre, Umeå University, Umeå, Sweden
| | - Benjamin Cossins
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Barcelona Supercomputing Center - BSC, Barcelona, Spain
| | - Carlos W. Bertoncini
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine – IRB Barcelona, Barcelona, Spain
| | - Victor Guallar
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Barcelona Supercomputing Center - BSC, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats - ICREA, Barcelona, Spain
| | - Magnus Wolf-Watz
- Department of Chemistry, Chemical Biological Centre, Umeå University, Umeå, Sweden
| | - Xavier Salvatella
- Joint BSC-CRG-IRB Research Programme in Computational Biology, Institute for Research in Biomedicine – IRB Barcelona, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats - ICREA, Barcelona, Spain
- * E-mail: (SEM); (XS)
| |
Collapse
|
50
|
Seyler SL, Beckstein O. Sampling large conformational transitions: adenylate kinase as a testing ground. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.919497] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|