1
|
Auerbach A. Dynamics of receptor activation by agonists. Biophys J 2024:S0006-3495(24)00003-1. [PMID: 38178577 DOI: 10.1016/j.bpj.2024.01.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Revised: 11/18/2023] [Accepted: 01/02/2024] [Indexed: 01/06/2024] Open
Abstract
How do agonists turn on receptors? The model system we have used to address this question is the adult-type skeletal muscle nicotinic acetylcholine receptor. This ligand-gated ion channel has two orthosteric sites (for neurotransmitters) in the extracellular domain linked to an allosteric site (a gate) in the transmembrane domain. The goal of this perspective is to summarize how measurements of agonist binding energy reveal the dynamics of the neurotransmitter sites and the fundamental link between binding and gating.
Collapse
Affiliation(s)
- Anthony Auerbach
- Department of Physiology and Biophysics, State University of New York at Buffalo, Buffalo, New York.
| |
Collapse
|
2
|
Grebenkov DS. Diffusion-Controlled Reactions: An Overview. Molecules 2023; 28:7570. [PMID: 38005291 PMCID: PMC10674959 DOI: 10.3390/molecules28227570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 10/20/2023] [Accepted: 10/21/2023] [Indexed: 11/26/2023] Open
Abstract
We review the milestones in the century-long development of the theory of diffusion-controlled reactions. Starting from the seminal work by von Smoluchowski, who recognized the importance of diffusion in chemical reactions, we discuss perfect and imperfect surface reactions, their microscopic origins, and the underlying mathematical framework. Single-molecule reaction schemes, anomalous bulk diffusions, reversible binding/unbinding kinetics, and many other extensions are presented. An alternative encounter-based approach to diffusion-controlled reactions is introduced, with emphasis on its advantages and potential applications. Some open problems and future perspectives are outlined.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée, CNRS-Ecole Polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France
| |
Collapse
|
3
|
Indurthi DC, Auerbach A. Agonist efficiency links binding and gating in a nicotinic receptor. eLife 2023; 12:e86496. [PMID: 37399234 DOI: 10.7554/elife.86496] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2023] [Accepted: 06/15/2023] [Indexed: 07/05/2023] Open
Abstract
Receptors signal by switching between resting (C) and active (O) shapes ('gating') under the influence of agonists. The receptor's maximum response depends on the difference in agonist binding energy, O minus C. In nicotinic receptors, efficiency (η) represents the fraction of agonist binding energy applied to a local rearrangement (an induced fit) that initiates gating. In this receptor, free energy changes in gating and binding can be interchanged by the conversion factor η. Efficiencies estimated from concentration-response curves (23 agonists, 53 mutations) sort into five discrete classes (%): 0.56 (17), 0.51(32), 0.45(13), 0.41(26), and 0.31(12), implying that there are 5 C versus O binding site structural pairs. Within each class efficacy and affinity are corelated linearly, but multiple classes hide this relationship. η unites agonist binding with receptor gating and calibrates one link in a chain of coupled domain rearrangements that comprises the allosteric transition of the protein.
Collapse
Affiliation(s)
- Dinesh C Indurthi
- Department of Physiology and Biophysics, University at Buffalo, State University of New York, Buffalo, United States
| | - Anthony Auerbach
- Department of Physiology and Biophysics, University at Buffalo, State University of New York, Buffalo, United States
| |
Collapse
|
4
|
Sakyi PO, Broni E, Amewu RK, Miller WA, Wilson MD, Kwofie SK. Homology Modeling, de Novo Design of Ligands, and Molecular Docking Identify Potential Inhibitors of Leishmania donovani 24-Sterol Methyltransferase. Front Cell Infect Microbiol 2022; 12:859981. [PMID: 35719359 PMCID: PMC9201040 DOI: 10.3389/fcimb.2022.859981] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 04/28/2022] [Indexed: 11/13/2022] Open
Abstract
The therapeutic challenges pertaining to leishmaniasis due to reported chemoresistance and toxicity necessitate the need to explore novel pathways to identify plausible inhibitory molecules. Leishmania donovani 24-sterol methyltransferase (LdSMT) is vital for the synthesis of ergosterols, the main constituents of Leishmania cellular membranes. So far, mammals have not been shown to possess SMT or ergosterols, making the pathway a prime candidate for drug discovery. The structural model of LdSMT was elucidated using homology modeling to identify potential novel 24-SMT inhibitors via virtual screening, scaffold hopping, and de-novo fragment-based design. Altogether, six potential novel inhibitors were identified with binding energies ranging from −7.0 to −8.4 kcal/mol with e-LEA3D using 22,26-azasterol and S1–S4 obtained from scaffold hopping via the ChEMBL, DrugBank, PubChem, ChemSpider, and ZINC15 databases. These ligands showed comparable binding energy to 22,26-azasterol (−7.6 kcal/mol), the main inhibitor of LdSMT. Moreover, all the compounds had plausible ligand efficiency-dependent lipophilicity (LELP) scores above 3. The binding mechanism identified Tyr92 to be critical for binding, and this was corroborated via molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The ligand A1 was predicted to possess antileishmanial properties with a probability of activity (Pa) of 0.362 and a probability of inactivity (Pi) of 0.066, while A5 and A6 possessed dermatological properties with Pa values of 0.205 and 0.249 and Pi values of 0.162 and 0.120, respectively. Structural similarity search via DrugBank identified vabicaserin, daledalin, zanapezil, imipramine, and cefradine with antileishmanial properties suggesting that the de-novo compounds could be explored as potential antileishmanial agents.
Collapse
Affiliation(s)
- Patrick O. Sakyi
- Department of Chemistry, School of Physical and Mathematical Sciences, College of Basic and Applied Sciences, University of Ghana, Accra, Ghana
- Department of Chemical Sciences, School of Sciences, University of Energy and Natural Resources, Sunyani, Ghana
| | - Emmanuel Broni
- Department of Biomedical Engineering, School of Engineering Sciences, College of Basic & Applied Sciences, University of Ghana, Accra, Ghana
- Department of Parasitology, Noguchi Memorial Institute for Medical Research (NMIMR), College of Health Sciences (CHS), University of Ghana, Accra, Ghana
| | - Richard K. Amewu
- Department of Chemistry, School of Physical and Mathematical Sciences, College of Basic and Applied Sciences, University of Ghana, Accra, Ghana
| | - Whelton A. Miller
- Department of Medicine, Loyola University Medical Center, Maywood, IL, United States
- Department of Molecular Pharmacology and Neuroscience, Loyola University Medical Center, Maywood, IL, United States
- Department of Chemical and Biomolecular Engineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA, United States
| | - Michael D. Wilson
- Department of Parasitology, Noguchi Memorial Institute for Medical Research (NMIMR), College of Health Sciences (CHS), University of Ghana, Accra, Ghana
- Department of Medicine, Loyola University Medical Center, Maywood, IL, United States
| | - Samuel Kojo Kwofie
- Department of Biomedical Engineering, School of Engineering Sciences, College of Basic & Applied Sciences, University of Ghana, Accra, Ghana
- Department of Biochemistry, Cell and Molecular Biology, West African Centre for Cell Biology of Infectious Pathogens, College of Basic and Applied Sciences, University of Ghana, Accra, Ghana
- *Correspondence: Samuel Kojo Kwofie,
| |
Collapse
|
5
|
Punia R, Goel G. Computation of the Protein Conformational Transition Pathway on Ligand Binding by Linear Response-Driven Molecular Dynamics. J Chem Theory Comput 2022; 18:3268-3283. [PMID: 35484642 DOI: 10.1021/acs.jctc.1c01243] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
While extremely important for relating the protein structure to its biological function, determination of the protein conformational transition pathway upon ligand binding is made difficult due to the transient nature of intermediates, a large and rugged conformational space, and coupling between protein dynamics and ligand-protein interactions. Existing methods that rely on prior knowledge of the bound (holo) state structure are restrictive. A second concern relates to the correspondence of intermediates obtained to the metastable states on the apo → holo transition pathway. Here, we have taken the protein apo structure and ligand-binding site as only inputs and combined an elastic network model (ENM) representation of the protein Hamiltonian with linear response theory (LRT) for protein-ligand interactions to identify the set of slow normal modes of protein vibrations that have a high overlap with the direction of the protein conformational change. The structural displacement along the chosen direction was performed using excited normal modes molecular dynamics (MDeNM) simulations rather than by the direct use of LRT. Herein, the MDeNM excitation velocity was optimized on-the-fly on the basis of its coupling to protein dynamics and ligand-protein interactions. Thus, a determined set of structures was validated against crystallographic and simulation data on four protein-ligand systems, namely, adenylate kinase-di(adenosine-5')pentaphosphate, ribose binding protein-β-d-ribopyranose, DNA β-glucosyltransferase-uridine-5'-diphosphate, and G-protein α subunit-guanosine-5'-triphosphate, which present important differences in protein conformational heterogeneity, ligand binding mechanism, viz. induced-fit or conformational selection, extent, and nonlinearity in protein conformational changes upon ligand binding, and presence of allosteric effects. The obtained set of intermediates was used as an input to path metadynamics simulations to obtain the free energy profile for the apo → holo transition.
Collapse
Affiliation(s)
- Rajat Punia
- Department of Chemical Engineering, Indian Institute of Technology Delhi, Hauz Khas, Delhi 110016, India
| | - Gaurav Goel
- Department of Chemical Engineering, Indian Institute of Technology Delhi, Hauz Khas, Delhi 110016, India
| |
Collapse
|
6
|
Ghorbani M, Prasad S, Klauda J, Brooks B. GraphVAMPNet, using graph neural networks and variational approach to Markov processes for dynamical modeling of biomolecules. J Chem Phys 2022; 156:184103. [PMID: 35568532 PMCID: PMC9094994 DOI: 10.1063/5.0085607] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Finding low dimensional representation of data from long-timescale trajectories of biomolecular processes such as protein-folding or ligand-receptor binding is of fundamental importance and kinetic models such as Markov modeling have proven useful in describing the kinetics of these systems. Recently, an unsupervised machine learning technique called VAMPNet was introduced to learn the low dimensional representation and linear dynamical model in an end-to-end manner. VAMPNet is based on variational approach to Markov processes (VAMP) and relies on neural networks to learn the coarse-grained dynamics. In this contribution, we combine VAMPNet and graph neural networks to generate an end-to-end framework to efficiently learn high-level dynamics and metastable states from the long-timescale molecular dynamics trajectories. This method bears the advantages of graph representation learning and uses graph message passing operations to generate an embedding for each datapoint which is used in the VAMPNet to generate a coarse-grained representation. This type of molecular representation results in a higher resolution and more interpretable Markov model than the standard VAMPNet enabling a more detailed kinetic study of the biomolecular processes. Our GraphVAMPNet approach is also enhanced with an attention mechanism to find the important residues for classification into different metastable states.
Collapse
Affiliation(s)
- Mahdi Ghorbani
- University of Maryland at College Park, United States of America
| | - Samarjeet Prasad
- National Heart Lung and Blood Institute, United States of America
| | - Jeffery Klauda
- Chemical and Biomolecular Engineering, University of Maryland at College Park, United States of America
| | - Bernard Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, United States of America
| |
Collapse
|
7
|
Köhn B, Schwarz P, Wittung-Stafshede P, Kovermann M. Impact of crowded environments on binding between protein and single-stranded DNA. Sci Rep 2021; 11:17682. [PMID: 34480058 PMCID: PMC8417293 DOI: 10.1038/s41598-021-97219-1] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Accepted: 08/20/2021] [Indexed: 11/09/2022] Open
Abstract
The concept of Molecular Crowding depicts the high density of diverse molecules present in the cellular interior. Here, we determine the impact of low molecular weight and larger molecules on binding capacity of single-stranded DNA (ssDNA) to the cold shock protein B (CspB). Whereas structural features of ssDNA-bound CspB are fully conserved in crowded environments as probed by high-resolution NMR spectroscopy, intrinsic fluorescence quenching experiments reveal subtle changes in equilibrium affinity. Kinetic stopped-flow data showed that DNA-to-protein association is significantly retarded independent of choice of the molecule that is added to the solution, but dissociation depends in a nontrivial way on its size and chemical characteristics. Thus, for this DNA-protein interaction, excluded volume effect does not play the dominant role but instead observed effects are dictated by the chemical properties of the crowder. We propose that surrounding molecules are capable of specific modification of the protein's hydration shell via soft interactions that, in turn, tune protein-ligand binding dynamics and affinity.
Collapse
Affiliation(s)
- Birgit Köhn
- Department of Chemistry, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany.,Konstanz Research School Chemical Biology KoRS-CB, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany
| | - Patricia Schwarz
- Department of Chemistry, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany
| | - Pernilla Wittung-Stafshede
- Department of Biology and Biological Engineering, Chalmers University of Technology, 41296, Gothenburg, Sweden
| | - Michael Kovermann
- Department of Chemistry, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany. .,Konstanz Research School Chemical Biology KoRS-CB, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany.
| |
Collapse
|
8
|
Ahn SH, Jagger BR, Amaro RE. Ranking of Ligand Binding Kinetics Using a Weighted Ensemble Approach and Comparison with a Multiscale Milestoning Approach. J Chem Inf Model 2020; 60:5340-5352. [PMID: 32315175 DOI: 10.1021/acs.jcim.9b00968] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
To improve lead optimization efforts in finding the right ligand, pharmaceutical industries need to know the ligand's binding kinetics, such as binding and unbinding rate constants, which often correlate with the ligand's efficacy in vivo. To predict binding kinetics efficiently, enhanced sampling methods, such as milestoning and the weighted ensemble (WE) method, have been used in molecular dynamics (MD) simulations of these systems. However, a comparison of these enhanced sampling methods in ranking ligands has not been done. Hence, a WE approach called the concurrent adaptive sampling (CAS) algorithm that uses MD simulations was used to rank seven ligands for β-cyclodextrin, a system in which a multiscale milestoning approach called simulation enabled estimation of kinetic rates (SEEKR) was also used, which uses both MD and Brownian dynamics simulations. Overall, the CAS algorithm can successfully rank ligands using the unbinding rate constant koff values and binding free energy ΔG values, as SEEKR did, with reduced computational cost that is about the same as SEEKR. We compare the CAS algorithm simulations with different parameters and discuss the impact of parameters in ranking ligands and obtaining rate constant and binding free energy estimates. We also discuss similarities and differences and advantages and disadvantages of SEEKR and the CAS algorithm for future use.
Collapse
Affiliation(s)
- Surl-Hee Ahn
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Benjamin R Jagger
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
9
|
Ray D, Gokey T, Mobley DL, Andricioaei I. Kinetics and free energy of ligand dissociation using weighted ensemble milestoning. J Chem Phys 2020; 153:154117. [PMID: 33092382 DOI: 10.1063/5.0021953] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider the recently developed weighted ensemble milestoning (WEM) scheme [D. Ray and I. Andricioaei, J. Chem. Phys. 152, 234114 (2020)] and test its capability of simulating ligand-receptor dissociation dynamics. We performed WEM simulations on the following host-guest systems: Na+/Cl- ion pair and 4-hydroxy-2-butanone ligand with FK506 binding protein. As a proof of principle, we show that the WEM formalism reproduces the Na+/Cl- ion pair dissociation timescale and the free energy profile obtained from long conventional MD simulation. To increase the accuracy of WEM calculations applied to kinetics and thermodynamics in protein-ligand binding, we introduced a modified WEM scheme called weighted ensemble milestoning with restraint release (WEM-RR), which can increase the number of starting points per milestone without adding additional computational cost. WEM-RR calculations obtained a ligand residence time and binding free energy in agreement with experimental and previous computational results. Moreover, using the milestoning framework, the binding time and rate constants, dissociation constants, and committor probabilities could also be calculated at a low computational cost. We also present an analytical approach for estimating the association rate constant (kon) when binding is primarily diffusion driven. We show that the WEM method can efficiently calculate multiple experimental observables describing ligand-receptor binding/unbinding and is a promising candidate for computer-aided inhibitor design.
Collapse
Affiliation(s)
- Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| | - Trevor Gokey
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| | - David L Mobley
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| | - Ioan Andricioaei
- Department of Chemistry, University of California Irvine, Irvine, California 92697, USA
| |
Collapse
|
10
|
Thomas T, Yuriev E, Chalmers DK. Markov State Model Analysis of Haloperidol Binding to the D 3 Dopamine Receptor. J Chem Theory Comput 2020; 16:3879-3888. [PMID: 32324998 DOI: 10.1021/acs.jctc.0c00013] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
We have developed Markov state models (MSMs) and hidden Markov models (HMMs) that describe the binding of haloperidol to the D3 dopamine receptor. Haloperidol is an antipsychotic drug that binds with nanomolar affinity to the D3 dopamine receptor, where it functions as an inverse agonist. The models were constructed using an adaptive sampling approach from 519 individual molecular dynamics simulations totaling 122 μs of simulated time and encompass the entire drug binding process. They reveal short-lived metastable bound states and two distinct long-lived bound conformations that cannot be separated in affinity using our current methodology. This work extends the use of MSMs and HMMs to study ligand binding, which thus far has been limited to simpler systems.
Collapse
Affiliation(s)
- Trayder Thomas
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Pde, Parkville, Victoria 3052, Australia
| | - Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Pde, Parkville, Victoria 3052, Australia
| | - David K Chalmers
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Pde, Parkville, Victoria 3052, Australia
| |
Collapse
|
11
|
Dreßler C, Kabbe G, Brehm M, Sebastiani D. Exploring non-equilibrium molecular dynamics of mobile protons in the solid acid CsH2PO4 at the micrometer and microsecond scale. J Chem Phys 2020; 152:164110. [DOI: 10.1063/5.0002167] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Christian Dreßler
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Gabriel Kabbe
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Martin Brehm
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Daniel Sebastiani
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| |
Collapse
|
12
|
Dreßler C, Kabbe G, Brehm M, Sebastiani D. Dynamical matrix propagator scheme for large-scale proton dynamics simulations. J Chem Phys 2020; 152:114114. [PMID: 32199428 DOI: 10.1063/1.5140635] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
We derive a matrix formalism for the simulation of long range proton dynamics for extended systems and timescales. On the basis of an ab initio molecular dynamics simulation, we construct a Markov chain, which allows us to store the entire proton dynamics in an M × M transition matrix (where M is the number of oxygen atoms). In this article, we start from common topology features of the hydrogen bond network of good proton conductors and utilize them as constituent constraints of our dynamic model. We present a thorough mathematical derivation of our approach and verify its uniqueness and correct asymptotic behavior. We propagate the proton distribution by means of transition matrices, which contain kinetic data from both ultra-short (sub-ps) and intermediate (ps) timescales. This concept allows us to keep the most relevant features from the microscopic level while effectively reaching larger time and length scales. We demonstrate the applicability of the transition matrices for the description of proton conduction trends in proton exchange membrane materials.
Collapse
Affiliation(s)
- Christian Dreßler
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Gabriel Kabbe
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Martin Brehm
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Daniel Sebastiani
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| |
Collapse
|
13
|
Rifai EA, Ferrario V, Pleiss J, Geerke DP. Combined Linear Interaction Energy and Alchemical Solvation Free-Energy Approach for Protein-Binding Affinity Computation. J Chem Theory Comput 2020; 16:1300-1310. [PMID: 31894691 PMCID: PMC7017367 DOI: 10.1021/acs.jctc.9b00890] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Calculating free energies of binding (ΔGbind) between ligands and their target protein is of major interest to drug discovery and safety, yet it is still associated with several challenges and difficulties. Linear interaction energy (LIE) is an efficient in silico method for ΔGbind computation. LIE models can be trained and used to directly calculate binding affinities from interaction energies involving ligands in the bound and unbound states only, and LIE can be combined with statistical weighting to calculate ΔGbind for flexible proteins that may bind their ligands in multiple orientations. Here, we investigate if LIE predictions can be effectively improved by explicitly including the entropy of (de)solvation into our free-energy calculations. For that purpose, we combine LIE calculations for the protein-ligand-bound state with explicit free-energy perturbation to rigorously compute the unbound ligand's solvation free energy. We show that for 28 Cytochrome P450 2A6 (CYP2A6) ligands, coupling LIE with alchemical solvation free-energy calculation helps to improve obtained correlation between computed and reference (experimental) binding data.
Collapse
Affiliation(s)
- Eko Aditya Rifai
- AIMMS Division of Molecular and Computational Toxicology, Department of Chemistry and Pharmaceutical Sciences , Vrije Universiteit Amsterdam , De Boelelaan 1108 , 1081 HZ Amsterdam , The Netherlands
| | - Valerio Ferrario
- Institute of Biochemistry and Technical Biochemistry , Universität Stuttgart , Allmandring 31 , 70569 Stuttgart , Germany
| | - Jürgen Pleiss
- Institute of Biochemistry and Technical Biochemistry , Universität Stuttgart , Allmandring 31 , 70569 Stuttgart , Germany
| | - Daan P Geerke
- AIMMS Division of Molecular and Computational Toxicology, Department of Chemistry and Pharmaceutical Sciences , Vrije Universiteit Amsterdam , De Boelelaan 1108 , 1081 HZ Amsterdam , The Netherlands
| |
Collapse
|
14
|
Kaur G, Dogra N, Singh S. Health Risk Assessment of Occupationally Pesticide-Exposed Population of Cancer Prone Area of Punjab. Toxicol Sci 2019; 165:157-169. [PMID: 29893964 DOI: 10.1093/toxsci/kfy140] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The alarming health issues especially the unusually high number of cancer cases in agriculture community of Bathinda district of Punjab (India) is a serious concern. There is limited knowledge about the role of gene-environment interactions in oncogenesis prevalent in this area. The aim of this study was to evaluate the association of oxidative stress with CYP1A2, CYP2B6, CYP2C9, CYP3A4, and PON1 genetic variation in the pesticide-exposed (occupationally) population of Bathinda district of Punjab (India). This study demonstrated significantly elevated relative risk (RR) of lower antioxidant defense mechanism (Glutathione, Catalase, Superoxide Dismutase, Glutathione peroxidases, and Glutathione Reductase) in occupationally pesticide-exposed group (n = 120) as compared with unexposed group (n = 84) from Bathinda district of Punjab (India). Our data shows pesticide exposure to be a major risk factor leading to increased oxidative stress inside the body. Gas chromatographic analysis revealed the residues of organophosphates (chlorpyriphos, dichlorvos, ethoprophos) and herbicides (atrazine, butachlor, alachlor, metolachlor) in the blood samples of the exposed population. In vitro results showed a dose dependent decrease in cell viability following treatment of pesticides detected in blood samples in hPBMCs and A549 cell line. Genetic variation analysis revealed missense mutations in CYP2B6 (2 mutations), CY3A4 (1 mutation), and CYP2C9 (2 mutations). The observed mutations have been predicted to cause structural and conformation change in protein structure which could result in altered stability. In first of its kind of study, our data reveal oxidative stress and pesticide residue accumulation inside the body to be the major reasons for health concerns in Bathinda district.
Collapse
Affiliation(s)
- Gurpreet Kaur
- Department of Environmental Science and Technology, Central University of Punjab, Bathinda, Punjab 151001, India
| | - Nilambra Dogra
- Department of Experimental Medicine and Biotechnology, Post Graduate Institute of Medical Education and Research, Chandigarh 160012, India
| | - Sandeep Singh
- Department of Human Genetics and Molecular Medicine, Central University of Punjab, Bathinda, Punjab 151001, India
| |
Collapse
|
15
|
Bernetti M, Masetti M, Recanatini M, Amaro RE, Cavalli A. An Integrated Markov State Model and Path Metadynamics Approach To Characterize Drug Binding Processes. J Chem Theory Comput 2019; 15:5689-5702. [DOI: 10.1021/acs.jctc.9b00450] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Mattia Bernetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Maurizio Recanatini
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093-0340, United States
| | - Andrea Cavalli
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum, Università di Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
- Computational & Chemical Biology, Istituto Italiano di Tecnologia, Via Morego 30, I-16163 Genova, Italy
| |
Collapse
|
16
|
Khan MT, Malik SI. Structural dynamics behind variants in pyrazinamidase and pyrazinamide resistance. J Biomol Struct Dyn 2019; 38:3003-3017. [PMID: 31357912 DOI: 10.1080/07391102.2019.1650113] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Pyrazinamide (PZA) is an important component of first-line anti-tuberculosis (anti-TB) drugs. The anti-TB agent is activated into an active form, pyrazinoic acid (POA), by Mycobacterium tuberculosis (MTB) pncA gene encoding pyrazinamidase (PZase). The major cause of PZA-resistance has been associated with mutations in the pncA gene. We have detected several novel mutations including V131F, Q141P, R154T, A170P, and V180F (GeneBank Accession No. MH461111) in the pncA gene of PZA-resistant isolates during PZA drug susceptibility testing followed by pncA gene sequencing. Here, we investigated molecular mechanism of PZA-resistance by comparing the results of experimental and molecular dynamics. The mutants (MTs) and wild type (WT) PZase structures in apo and complex with PZA were subjected to molecular dynamic simulations (MD) at the 40 ns. Multiple factors, including root mean square deviations (RMSD), binding pocket, total energy, dynamic cross correlation, and root mean square fluctuations (RMSF) of MTs and WT were compared. The MTs attained a high deviation and fluctuation compared to WT. Binding pocket volumes of the MTs, were found, lower than the WT, and the docking scores were high than WT while shape complementarity scores were lower than that of the WT. Residual motion in MTs are seemed to be dominant in anti-correlated motion. Mutations at locations, V131F, Q141P, R154T, A170P, and V180F, might be involved in the structural changes, possibly affecting the catalytic property of PZase to convert PZA into POA. Our study provides useful information that will enhance the understanding for better management of TB. AbbreviationsDSTdrug susceptibility testingΔelecelectrostatic energyLJLowenstein-Jensen mediumMGITmycobacterium growth indicator tubesMTsmutantsMDmolecular dynamic simulationsMTBMycobacterium tuberculosisNALC-NaOHN-acetyl-l-cysteine-sodium hydroxideNIHNational Institutes of HealthNPTamount of substance (N), pressure (P) temperature (T)NVTmoles (N), volume (V) temperature (T)PZasepyrazinamidaseΔpspolar solvation energyPTRLProvincial Tuberculosis Reference LaboratoryRMSDroot mean square deviationsRMSFroot mean square fluctuationsΔSASAsolvent accessible surface area energyTBtuberculosisGTotaltotal binding free energyΔvdWVan der Waals energyWTwild typeCommunicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Shaukat Iqbal Malik
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| |
Collapse
|
17
|
Capelli R, Carloni P, Parrinello M. Exhaustive Search of Ligand Binding Pathways via Volume-Based Metadynamics. J Phys Chem Lett 2019; 10:3495-3499. [PMID: 31188006 DOI: 10.1021/acs.jpclett.9b01183] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Determining the complete set of ligands' binding-unbinding pathways is important for drug discovery and for rational interpretation of mutation data. Here we have developed a metadynamics-based technique that addresses this issue and allows estimating affinities in the presence of multiple escape pathways. Our approach is shown on a lysozyme T4 variant in complex with a benzene molecule. The calculated binding free energy is in agreement with experimental data. Remarkably, not only were we able to find all the previously identified ligand binding pathways, but also we identified three pathways previously not identified as such. These results were obtained at a small computational cost, making this approach valuable for practical applications, such as screening of small compound libraries.
Collapse
Affiliation(s)
- Riccardo Capelli
- Computational Biomedicine (INM-9/IAS-5) , Forschungszentrum Jülich , Wilhelm-Johnen-Straße , D-52425 Jülich , Germany
- JARA-HPC, Forschungszentrum Jülich , D-54245 Jülich , Germany
| | - Paolo Carloni
- Computational Biomedicine (INM-9/IAS-5) , Forschungszentrum Jülich , Wilhelm-Johnen-Straße , D-52425 Jülich , Germany
- Department of Physics , RWTH Aachen University , D-52078 Aachen , Germany
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences , ETH Zürich , c/o USI Campus, Via Giuseppe Buffi 13 , CH-6900 Lugano , Ticino Switzerland
- Facoltà di Informatica, Istituto di Scienze Computazionali , Università della Svizzera Italiana (USI) , Via Giuseppe Buffi 13 , CH-6900 Lugano , Ticino Switzerland
- Istituto Italiano di Tecnologia , Via Morego 30 , I-16163 Genova , Italy
| |
Collapse
|
18
|
Gonzalez TL, Rae JM, Colacino JA, Richardson RJ. Homology models of mouse and rat estrogen receptor- α ligand-binding domain created by in silico mutagenesis of a human template: molecular docking with 17ß-estradiol, diethylstilbestrol, and paraben analogs. COMPUTATIONAL TOXICOLOGY (AMSTERDAM, NETHERLANDS) 2019; 10:1-16. [PMID: 30740556 PMCID: PMC6363358 DOI: 10.1016/j.comtox.2018.11.003] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Crystal structures exist for human, but not rodent, estrogen receptor-α ligand-binding domain (ERα-LBD). Consequently, rodent studies involving binding of compounds to ERα-LBD are limited in their molecular-level interpretation and extrapolation to humans. Because the sequences of rodent and human ERα-LBDs are > 95% identical, we expected their 3D structures and ligand binding to be highly similar. To test this hypothesis, we used the human ERα-LBD structure (PDB 3UUD) as a template to produce rat and mouse homology models. Employing the rodent models and human structure, we generated docking poses of 23 Group A ligands (17ß-estradiol, diethylstilbestrol, and 21 paraben analogs) in AutoDock Vina for interspecies comparisons. Ligand RMSDs (Å) (median, 95% CI) were 0.49 (0.21-1.82) (human-mouse) and 1.19 (0.22-1.82) (human-rat), well below the 2.0-2.5 Å range for equivalent docking poses. Numbers of interspecies ligand-receptor residue contacts were highly similar, with Sorensen Sc (%) = 96.8 (90.0-100) (human-mouse) and 97.7 (89.5-100) (human-rat). Likewise, numbers of interspecies ligand-receptor residue contacts were highly correlated: Pearson r = 0.913 (human-mouse) and 0.925 (human-rat). Numbers of interspecies ligand-receptor atom contacts were even more tightly correlated: r = 0.979 (human-mouse) and 0.986 (human-rat). Pyramid plots of numbers of ligand-receptor atom contacts by residue exhibited high interspecies symmetry and had Spearman r s = 0.977 (human-mouse) and 0.966 (human-rat). Group B ligands included 15 ring-substituted parabens recently shown experimentally to exhibit decreased binding to human ERα and to exert increased antimicrobial activity. Ligand efficiencies calculated from docking ligands into human ERα-LBD were well correlated with those derived from published experimental data (Pearson partial r p = 0.894 and 0.918; Groups A and B, respectively). Overall, the results indicate that our constructed rodent ERα-LBDs interact with ligands in like manner to the human receptor, thus providing a high level of confidence in extrapolations of rodent to human ligand-receptor interactions.
Collapse
Affiliation(s)
- Thomas L. Gonzalez
- Department of Environmental Health Sciences, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
| | - James M. Rae
- Division of Hematology and Oncology, Department of Internal Medicine, University of Michigan Medical School, Ann Arbor, Michigan 48109, USA
| | - Justin A. Colacino
- Department of Environmental Health Sciences, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
- Department of Nutritional Sciences, University of Michigan School of Public Health, Ann Arbor, MI 48109 USA
- Center for Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109 USA
| | - Rudy J. Richardson
- Department of Environmental Health Sciences, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
- Center for Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109 USA
- Department of Neurology, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| |
Collapse
|
19
|
Röder K, Joseph JA, Husic BE, Wales DJ. Energy Landscapes for Proteins: From Single Funnels to Multifunctional Systems. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201800175] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Konstantin Röder
- Department of ChemistryUniversity of CambridgeLensfield Road CB2 1EW Cambridge UK
| | - Jerelle A. Joseph
- Department of ChemistryUniversity of CambridgeLensfield Road CB2 1EW Cambridge UK
| | - Brooke E. Husic
- Department of ChemistryUniversity of CambridgeLensfield Road CB2 1EW Cambridge UK
| | - David J. Wales
- Department of ChemistryUniversity of CambridgeLensfield Road CB2 1EW Cambridge UK
| |
Collapse
|
20
|
Yenjai S, Kumar CV, Kuno M, Liwporncharoenvong T, Samosorn S, Buranaprapuk A. Tuning the chain length of new pyrene derivatives for site-selective photocleavage of avidin. JOURNAL OF PHOTOCHEMISTRY AND PHOTOBIOLOGY. B, BIOLOGY 2018; 186:23-30. [PMID: 29990670 DOI: 10.1016/j.jphotobiol.2018.07.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/10/2018] [Revised: 06/08/2018] [Accepted: 07/02/2018] [Indexed: 06/08/2023]
Abstract
Rational design of photoreagents with systematic modifications of their structures can provide valuable information for a better understanding of the protein photocleavage mechanism by these reagents. Variation of the length of the linker connecting the photoactive moiety with the protein anchoring-group allowed us to investigate the control of the protein photocleavage site. A series of new photochemical reagents (PMA-1A, PMA-2A and PMA-3A) with increasing chain lengths is examined in the current study. Using avidin as a model system, we examined the interaction of these probes by UV-Vis, fluorescence spectroscopic methods, photocleavage and computational docking studies. Hypochromism of the absorption spectrum was observed for the binding of these new photochemical reagents with estimated binding constants (Kb) of 6.2 × 105, 6.7 × 105 and 4.6 × 105 M-1, respectively. No significant changes of Stern-Volmer quenching constant (Ksv) with Co(NH3)6Cl3 has been noted and the data indicated that the probes bind near the surface of the protein with sufficient exposure to the solvent. Photoexcitation of the probe-avidin complex, in the presence of Co(NH3)6Cl3, resulted in protein fragmentation, and the cleavage yield decreased with the increase in the linker length, and paralleled with the observed Ksv values. Amino acid sequencing of the photofragments indicated that avidin is cleaved between Thr77 and Val78, as a major cleavage site for all the three photoreagents. This site is proximate to the biotin binding site on avidin, and molecular docking studies indicated that the H-bonding interactions between the polar end-group of the photoreagents and hydrophilic amino acids of avidin were important in positioning the reagent on the protein. The major cleavage site, at residues 77-78, was within 5 Å of the pyrenyl moiety of the probe, and hence, molecular tuning of the linker provided a simple approach to position the photoreagent along the potential photocleavage site.
Collapse
Affiliation(s)
- Sudarat Yenjai
- Department of Chemistry, Faculty of Science, Srinakharinwirot University, Sukhumvit 23, Bangkok 10110, Thailand
| | - Challa V Kumar
- Department of Chemistry, 55 N. Eagleville Road, University of Connecticut, Storrs, CT 06269-3060, USA; Department of Molecular and Cellular Biology, 91 N. Eagleville Road, U-3125, University of Connecticut, Storrs, CT 06269-3125, USA
| | - Mayuso Kuno
- Department of Chemistry, Faculty of Science, Srinakharinwirot University, Sukhumvit 23, Bangkok 10110, Thailand
| | | | - Siritron Samosorn
- Department of Chemistry, Faculty of Science, Srinakharinwirot University, Sukhumvit 23, Bangkok 10110, Thailand
| | - Apinya Buranaprapuk
- Department of Chemistry, Faculty of Science, Srinakharinwirot University, Sukhumvit 23, Bangkok 10110, Thailand.
| |
Collapse
|
21
|
Zeng X, Zhu L, Zheng X, Cecchini M, Huang X. Harnessing complexity in molecular self-assembly using computer simulations. Phys Chem Chem Phys 2018; 20:6767-6776. [PMID: 29479585 DOI: 10.1039/c7cp06181a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In molecular self-assembly, hundreds of thousands of freely-diffusing molecules associate to form ordered and functional architectures in the absence of an actuator. This intriguing phenomenon plays a critical role in biology and has become a powerful tool for the fabrication of advanced nanomaterials. Due to the limited spatial and temporal resolutions of current experimental techniques, computer simulations offer a complementary strategy to explore self-assembly with atomic resolution. Here, we review recent computational studies focusing on both thermodynamic and kinetic aspects. As we shall see, thermodynamic approaches based on modeling and statistical mechanics offer initial guidelines to design nanostructures with modest computational effort. Computationally more intensive analyses based on molecular dynamics simulations and kinetic network models (KNMs) reach beyond it, opening the door to the rational design of self-assembly pathways. Current limitations of these methodologies are discussed. We anticipate that the synergistic use of thermodynamic and kinetic analyses based on computer simulations will provide an important contribution to the de novo design of self-assembly.
Collapse
Affiliation(s)
- Xiangze Zeng
- Department of Chemistry, Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | | | | | | | | |
Collapse
|
22
|
Zhu L, Sheong FK, Zeng X, Huang X. Elucidation of the conformational dynamics of multi-body systems by construction of Markov state models. Phys Chem Chem Phys 2018; 18:30228-30235. [PMID: 27314275 DOI: 10.1039/c6cp02545e] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Constructing Markov State Models (MSMs) based on short molecular dynamics simulations is a powerful computational technique to complement experiments in predicting long-time kinetics of biomolecular processes at atomic resolution. Even though the MSM approach has been widely applied to study one-body processes such as protein folding and enzyme conformational changes, the majority of biological processes, e.g. protein-ligand recognition, signal transduction, and protein aggregation, essentially involve multiple entities. Here we review the attempts at constructing MSMs for multi-body systems, point out the challenges therein and discuss recent algorithmic progresses that alleviate these challenges. In particular, we describe an automatic kinetics based partitioning method that achieves optimal definition of the conformational states in a multi-body system, and discuss a novel maximum-likelihood approach that efficiently estimates the slow uphill kinetics utilizing pre-computed equilibrium populations of all states. We expect that these new algorithms and their combinations may boost investigations of important multi-body biological processes via the efficient construction of MSMs.
Collapse
Affiliation(s)
- Lizhe Zhu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Fu Kit Sheong
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China.
| | - Xiangze Zeng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| |
Collapse
|
23
|
Abstract
Whereas protein-ligand binding affinities have long-established prominence, binding rate constants and binding mechanisms have gained increasing attention in recent years. Both new computational methods and new experimental techniques have been developed to characterize the latter properties. It is now realized that binding mechanisms, like binding rate constants, can and should be quantitatively determined. In this review, we summarize studies and synthesize ideas on several topics in the hope of providing a coherent picture of and physical insight into binding kinetics. The topics include microscopic formulation of the kinetic problem and its reduction to simple rate equations; computation of binding rate constants; quantitative determination of binding mechanisms; and elucidation of physical factors that control binding rate constants and mechanisms.
Collapse
Affiliation(s)
- Xiaodong Pang
- Department of Physics, Florida State University, Tallahassee, Florida 32306; .,Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306
| | - Huan-Xiang Zhou
- Department of Physics, Florida State University, Tallahassee, Florida 32306; .,Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida 32306
| |
Collapse
|
24
|
Bernetti M, Cavalli A, Mollica L. Protein-ligand (un)binding kinetics as a new paradigm for drug discovery at the crossroad between experiments and modelling. MEDCHEMCOMM 2017; 8:534-550. [PMID: 30108770 PMCID: PMC6072069 DOI: 10.1039/c6md00581k] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Accepted: 01/25/2017] [Indexed: 12/14/2022]
Abstract
In the last three decades, protein and nucleic acid structure determination and comprehension of the mechanisms, leading to their physiological and pathological functions, have become a cornerstone of biomedical sciences. A deep understanding of the principles governing the fates of cells and tissue at the molecular level has been gained over the years, offering a solid basis for the rational design of drugs aimed at the pharmacological treatment of numerous diseases. Historically, affinity indicators (i.e. Kd and IC50/EC50) have been assumed to be valid indicators of the in vivo efficacy of a drug. However, recent studies pointed out that the kinetics of the drug-receptor binding process could be as important or even more important than affinity in determining the drug efficacy. This eventually led to a growing interest in the characterisation and prediction of the rate constants of protein-ligand association and dissociation. For instance, a drug with a longer residence time can kinetically select a given receptor over another, even if the affinity for both receptors is comparable, thus increasing its therapeutic index. Therefore, understanding the molecular features underlying binding and unbinding processes is of central interest towards the rational control of drug binding kinetics. In this review, we report the theoretical framework behind protein-ligand association and highlight the latest advances in the experimental and computational approaches exploited to investigate the binding kinetics.
Collapse
Affiliation(s)
- M Bernetti
- Department of Pharmacy and Biotechnology , University of Bologna , via Belmeloro 6 , 40126 Bologna , Italy
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , 16163 Genova , Italy .
| | - A Cavalli
- Department of Pharmacy and Biotechnology , University of Bologna , via Belmeloro 6 , 40126 Bologna , Italy
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , 16163 Genova , Italy .
| | - L Mollica
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , 16163 Genova , Italy .
| |
Collapse
|
25
|
Moritsugu K, Terada T, Kidera A. Free-Energy Landscape of Protein–Ligand Interactions Coupled with Protein Structural Changes. J Phys Chem B 2017; 121:731-740. [DOI: 10.1021/acs.jpcb.6b11696] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kei Moritsugu
- Graduate
School of Medical Life Science, Yokohama City University, 1-7-29
Suehirocho, Tsurumi, Yokohama 230-0045, Japan
| | - Tohru Terada
- Graduate
School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo, Tokyo 113-8657, Japan
| | - Akinori Kidera
- Graduate
School of Medical Life Science, Yokohama City University, 1-7-29
Suehirocho, Tsurumi, Yokohama 230-0045, Japan
| |
Collapse
|
26
|
Yuan S, Chan HS, Hu Z. Using
PyMOL
as a platform for computational drug design. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1298] [Citation(s) in RCA: 113] [Impact Index Per Article: 16.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Shuguang Yuan
- Laboratory of Physical Chemistry of Polymers and MembranesEcole Polytechnique Fédérale de Lausanne (EPFL) Lausanne Switzerland
| | | | - Zhenquan Hu
- High Magnetic Field LaboratoryChinese Academy of Science Hefei China
| |
Collapse
|
27
|
Elucidating Mechanisms of Molecular Recognition Between Human Argonaute and miRNA Using Computational Approaches. Methods Mol Biol 2016. [PMID: 27924488 DOI: 10.1007/978-1-4939-6563-2_18] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
MicroRNA (miRNA) and Argonaute (AGO) protein together form the RNA-induced silencing complex (RISC) that plays an essential role in the regulation of gene expression. Elucidating the underlying mechanism of AGO-miRNA recognition is thus of great importance not only for the in-depth understanding of miRNA function but also for inspiring new drugs targeting miRNAs. In this chapter we introduce a combined computational approach of molecular dynamics (MD) simulations, Markov state models (MSMs), and protein-RNA docking to investigate AGO-miRNA recognition. Constructed from MD simulations, MSMs can elucidate the conformational dynamics of AGO at biologically relevant timescales. Protein-RNA docking can then efficiently identify the AGO conformations that are geometrically accessible to miRNA. Using our recent work on human AGO2 as an example, we explain the rationale and the workflow of our method in details. This combined approach holds great promise to complement experiments in unraveling the mechanisms of molecular recognition between large, flexible, and complex biomolecules.
Collapse
|
28
|
Li Y, Li X, Dong Z. Exploration of gated ligand binding recognizes an allosteric site for blocking FABP4-protein interaction. Phys Chem Chem Phys 2016; 17:32257-67. [PMID: 26580122 DOI: 10.1039/c5cp04784f] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Fatty acid binding protein 4 (FABP4), reversibly binding to fatty acids and other lipids with high affinities, is a potential target for treatment of cancers. The binding site of FABP4 is buried in an interior cavity and thereby ligand binding/unbinding is coupled with opening/closing of FABP4. It is a difficult task both experimentally and computationally to illuminate the entry or exit pathway, especially with the conformational gating. In this report we combine extensive computer simulations, clustering analysis, and the Markov state model to investigate the binding mechanism of FABP4 and troglitazone. Our simulations capture spontaneous binding and unbinding events as well as the conformational transition of FABP4 between the open and closed states. An allosteric binding site on the protein surface is recognized for the development of novel FABP4 inhibitors. The binding affinity is calculated and compared with the experimental value. The kinetic analysis suggests that ligand residence on the protein surface may delay the binding process. Overall, our results provide a comprehensive picture of ligand diffusion on the protein surface, ligand migration into the buried cavity, and the conformational change of FABP4 at an atomic level.
Collapse
Affiliation(s)
- Yan Li
- The Hormel Institute, University of Minnesota, Austin Minnesota 55912, USA.
| | - Xiang Li
- Department of Physiology and Pathophysiology, School of Basic Medical Sciences, Zhengzhou University, 450001 Zhengzhou, Henan, China
| | - Zigang Dong
- The Hormel Institute, University of Minnesota, Austin Minnesota 55912, USA.
| |
Collapse
|
29
|
Bucci A, Yu TQ, Vanden-Eijnden E, Abrams CF. Kinetics of O2 Entry and Exit in Monomeric Sarcosine Oxidase via Markovian Milestoning Molecular Dynamics. J Chem Theory Comput 2016; 12:2964-72. [PMID: 27168219 DOI: 10.1021/acs.jctc.6b00071] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The flavoenzyme monomeric sarcosine oxidase (MSOX) catalyzes a complex set of reactions currently lacking a consensus mechanism. A key question that arises in weighing competing mechanistic models of MSOX function is to what extent ingress of O2 from the solvent (and its egress after an unsuccessful oxidation attempt) limits the overall catalytic rate. To address this question, we have applied to the MSOX/O2 system the relatively new simulation method of Markovian milestoning molecular dynamics simulations, which, as we recently showed [ Yu et al. J. Am. Chem. Soc. 2015 , 137 , 3041 ], accurately predicted the entry and exit kinetics of CO in myoglobin. We show that the mechanism of O2 entry and exit, in terms of which possible solvent-to-active-site channels contribute to the flow of O2, is sensitive to the presence of the substrate-mimicking competitive inhibitor 2-furoate in the substrate site. The second-order O2 entry rate constants were computed to be 8.1 × 10(6) and 3.1 × 10(6) M(-1) s(-1) for bound and apo MSOX, respectively, both of which moderately exceed the experimentally determined second-order rate constant of (2.83 ± 0.07) × 10(5) M(-1) s(-1) for flavin oxidation by O2 in MSOX. This suggests that the rate of flavin oxidation by O2 is likely not strongly limited by diffusion from the solvent to the active site. The first-order exit rate constants were computed to be 10(7) s(-1) and 7.2 × 10(6) s(-1) for the apo and bound states, respectively. The predicted faster entry and slower exit of O2 for the bound state indicate a longer residence time within MSOX, increasing the likelihood of collisions with the flavin isoalloxazine ring, a step required for reduction of molecular O2 and subsequent reoxidation of the flavin. This is also indirectly supported by previous experimental evidence favoring the so-called modified ping-pong mechanism, the distinguishing feature of which is an intermediate complex involving O2, the flavin, and the oxidized substrate simultaneously in the cavity. These findings demonstrate the utility of the Markovian milestoning approach in contributing new understanding of complicated enyzmatic function.
Collapse
Affiliation(s)
- Anthony Bucci
- Department of Chemical and Biological Engineering, Drexel University , Philadelphia, Pennsylvania 19104, United States
| | - Tang-Qing Yu
- Courant Institute of Mathematical Sciences, New York University , New York, New York 10012, United States
| | - Eric Vanden-Eijnden
- Courant Institute of Mathematical Sciences, New York University , New York, New York 10012, United States
| | - Cameron F Abrams
- Department of Chemical and Biological Engineering, Drexel University , Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
30
|
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
|
31
|
Rudzinski JF, Kremer K, Bereau T. Communication: Consistent interpretation of molecular simulation kinetics using Markov state models biased with external information. J Chem Phys 2016; 144:051102. [DOI: 10.1063/1.4941455] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | - Tristan Bereau
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| |
Collapse
|
32
|
Du X, Li Y, Xia YL, Ai SM, Liang J, Sang P, Ji XL, Liu SQ. Insights into Protein-Ligand Interactions: Mechanisms, Models, and Methods. Int J Mol Sci 2016; 17:ijms17020144. [PMID: 26821017 PMCID: PMC4783878 DOI: 10.3390/ijms17020144] [Citation(s) in RCA: 706] [Impact Index Per Article: 88.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2015] [Revised: 01/13/2016] [Accepted: 01/18/2016] [Indexed: 01/16/2023] Open
Abstract
Molecular recognition, which is the process of biological macromolecules interacting with each other or various small molecules with a high specificity and affinity to form a specific complex, constitutes the basis of all processes in living organisms. Proteins, an important class of biological macromolecules, realize their functions through binding to themselves or other molecules. A detailed understanding of the protein–ligand interactions is therefore central to understanding biology at the molecular level. Moreover, knowledge of the mechanisms responsible for the protein-ligand recognition and binding will also facilitate the discovery, design, and development of drugs. In the present review, first, the physicochemical mechanisms underlying protein–ligand binding, including the binding kinetics, thermodynamic concepts and relationships, and binding driving forces, are introduced and rationalized. Next, three currently existing protein-ligand binding models—the “lock-and-key”, “induced fit”, and “conformational selection”—are described and their underlying thermodynamic mechanisms are discussed. Finally, the methods available for investigating protein–ligand binding affinity, including experimental and theoretical/computational approaches, are introduced, and their advantages, disadvantages, and challenges are discussed.
Collapse
Affiliation(s)
- Xing Du
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
| | - Yi Li
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
| | - Yuan-Ling Xia
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
| | - Shi-Meng Ai
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
- Department of Applied Mathematics, Yunnan Agricultural University, Kunming 650201, China.
| | - Jing Liang
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
| | - Peng Sang
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
- Laboratory of Molecular Cardiology, Department of Cardiology, The First Affiliated Hospital of Kunming Medical University, Kunming 650032, China.
| | - Xing-Lai Ji
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
- Key Laboratory for Tumor molecular biology of High Education in Yunnan Province, School of Life Sciences, Yunnan University, Kunming 650091, China.
| | - Shu-Qun Liu
- Laboratory for Conservation and Utilization of Bio-Resources, Yunnan University, Kunming 650091, China.
- Key Laboratory for Tumor molecular biology of High Education in Yunnan Province, School of Life Sciences, Yunnan University, Kunming 650091, China.
| |
Collapse
|
33
|
Zeng X, Li B, Qiao Q, Zhu L, Lu ZY, Huang X. Elucidating dominant pathways of the nano-particle self-assembly process. Phys Chem Chem Phys 2016; 18:23494-9. [DOI: 10.1039/c6cp01808d] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Self-assembly processes play a key role in the fabrication of functional nano-structures with wide application in drug delivery and micro-reactors.
Collapse
Affiliation(s)
- Xiangze Zeng
- Department of Chemistry
- The Hong Kong University of Science and Technology
- Kowloon
- Hong Kong
| | - Bin Li
- State Key Laboratory of Supramolecular Structure and Materials
- Institute of Theoretical Chemistry
- Jilin University
- Changchun
- China
| | - Qin Qiao
- Department of Chemistry
- The Hong Kong University of Science and Technology
- Kowloon
- Hong Kong
| | - Lizhe Zhu
- Department of Chemistry
- The Hong Kong University of Science and Technology
- Kowloon
- Hong Kong
| | - Zhong-Yuan Lu
- State Key Laboratory of Supramolecular Structure and Materials
- Institute of Theoretical Chemistry
- Jilin University
- Changchun
- China
| | - Xuhui Huang
- Department of Chemistry
- The Hong Kong University of Science and Technology
- Kowloon
- Hong Kong
- Division of Biomedical Engineering
| |
Collapse
|
34
|
Boninsegna L, Gobbo G, Noé F, Clementi C. Investigating Molecular Kinetics by Variationally Optimized Diffusion Maps. J Chem Theory Comput 2015; 11:5947-60. [DOI: 10.1021/acs.jctc.5b00749] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Lorenzo Boninsegna
- Center
for Theoretical Biological Physics and Department of Chemistry, Rice University, 6100 Main Street, Houston, Texas 77005, United States
| | - Gianpaolo Gobbo
- Maxwell
Institute for Mathematical Sciences and School of Mathematics, The University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom
| | - Frank Noé
- Department
of Mathematics, Computer Science and Bioinformatics, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Cecilia Clementi
- Center
for Theoretical Biological Physics and Department of Chemistry, Rice University, 6100 Main Street, Houston, Texas 77005, United States
| |
Collapse
|
35
|
Trendelkamp-Schroer B, Wu H, Paul F, Noé F. Estimation and uncertainty of reversible Markov models. J Chem Phys 2015; 143:174101. [DOI: 10.1063/1.4934536] [Citation(s) in RCA: 79] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
| | - Hao Wu
- Institut für Mathematik und Informatik, FU Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Fabian Paul
- Institut für Mathematik und Informatik, FU Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Frank Noé
- Institut für Mathematik und Informatik, FU Berlin, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
36
|
Scherer MK, Trendelkamp-Schroer B, Paul F, Pérez-Hernández G, Hoffmann M, Plattner N, Wehmeyer C, Prinz JH, Noé F. PyEMMA 2: A Software Package for Estimation, Validation, and Analysis of Markov Models. J Chem Theory Comput 2015; 11:5525-42. [PMID: 26574340 DOI: 10.1021/acs.jctc.5b00743] [Citation(s) in RCA: 678] [Impact Index Per Article: 75.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Markov (state) models (MSMs) and related models of molecular kinetics have recently received a surge of interest as they can systematically reconcile simulation data from either a few long or many short simulations and allow us to analyze the essential metastable structures, thermodynamics, and kinetics of the molecular system under investigation. However, the estimation, validation, and analysis of such models is far from trivial and involves sophisticated and often numerically sensitive methods. In this work we present the open-source Python package PyEMMA ( http://pyemma.org ) that provides accurate and efficient algorithms for kinetic model construction. PyEMMA can read all common molecular dynamics data formats, helps in the selection of input features, provides easy access to dimension reduction algorithms such as principal component analysis (PCA) and time-lagged independent component analysis (TICA) and clustering algorithms such as k-means, and contains estimators for MSMs, hidden Markov models, and several other models. Systematic model validation and error calculation methods are provided. PyEMMA offers a wealth of analysis functions such that the user can conveniently compute molecular observables of interest. We have derived a systematic and accurate way to coarse-grain MSMs to few states and to illustrate the structures of the metastable states of the system. Plotting functions to produce a manuscript-ready presentation of the results are available. In this work, we demonstrate the features of the software and show new methodological concepts and results produced by PyEMMA.
Collapse
Affiliation(s)
- Martin K Scherer
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | | | - Fabian Paul
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Guillermo Pérez-Hernández
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Moritz Hoffmann
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Nuria Plattner
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Christoph Wehmeyer
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Jan-Hendrik Prinz
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Frank Noé
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| |
Collapse
|
37
|
Zhou J, Li M, Chen N, Wang S, Luo HB, Zhang Y, Wu R. Computational design of a time-dependent histone deacetylase 2 selective inhibitor. ACS Chem Biol 2015; 10:687-92. [PMID: 25546141 PMCID: PMC4372102 DOI: 10.1021/cb500767c] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Development of isoform-selective
histone deacetylase (HDAC) inhibitors is of great biological and medical
interest. Among 11 zinc-dependent HDAC isoforms, it is particularly
challenging to achieve isoform inhibition selectivity between HDAC1
and HDAC2 due to their very high structural similarities. In this
work, by developing and applying a novel de novo reaction-mechanism-based
inhibitor design strategy to exploit the reactivity difference, we
have discovered the first HDAC2-selective inhibitor, β-hydroxymethyl
chalcone. Our bioassay experiments show that this new compound has
a unique time-dependent selective inhibition on HDAC2, leading to
about 20-fold isoform-selectivity against HDAC1. Furthermore, our
ab initio QM/MM molecular dynamics simulations, a state-of-the-art
approach to study reactions in biological systems, have elucidated
how the β-hydroxymethyl chalcone can achieve the distinct time-dependent
inhibition toward HDAC2.
Collapse
Affiliation(s)
- Jingwei Zhou
- School of Pharmaceutical
Sciences, Sun Yat-sen University, Guangzhou 510006, P.R. China
| | - Min Li
- School of Pharmaceutical
Sciences, Sun Yat-sen University, Guangzhou 510006, P.R. China
| | - Nanhao Chen
- School of Pharmaceutical
Sciences, Sun Yat-sen University, Guangzhou 510006, P.R. China
| | - Shenglong Wang
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Hai-Bin Luo
- School of Pharmaceutical
Sciences, Sun Yat-sen University, Guangzhou 510006, P.R. China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, New York 10003, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, P.R. China
| | - Ruibo Wu
- School of Pharmaceutical
Sciences, Sun Yat-sen University, Guangzhou 510006, P.R. China
| |
Collapse
|
38
|
Sheong FK, Silva DA, Meng L, Zhao Y, Huang X. Automatic state partitioning for multibody systems (APM): an efficient algorithm for constructing Markov state models to elucidate conformational dynamics of multibody systems. J Chem Theory Comput 2014; 11:17-27. [PMID: 26574199 DOI: 10.1021/ct5007168] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The conformational dynamics of multibody systems plays crucial roles in many important problems. Markov state models (MSMs) are powerful kinetic network models that can predict long-time-scale dynamics using many short molecular dynamics simulations. Although MSMs have been successfully applied to conformational changes of individual proteins, the analysis of multibody systems is still a challenge because of the complexity of the dynamics that occur on a mixture of drastically different time scales. In this work, we have developed a new algorithm, automatic state partitioning for multibody systems (APM), for constructing MSMs to elucidate the conformational dynamics of multibody systems. The APM algorithm effectively addresses different time scales in the multibody systems by directly incorporating dynamics into geometric clustering when identifying the metastable conformational states. We have applied the APM algorithm to a 2D potential that can mimic a protein-ligand binding system and the aggregation of two hydrophobic particles in water and have shown that it can yield tremendous enhancements in the computational efficiency of MSM construction and the accuracy of the models.
Collapse
Affiliation(s)
- Fu Kit Sheong
- HKUST Shenzhen Research Institute , Nanshan, Shenzhen 518057, China
| | | | - Luming Meng
- HKUST Shenzhen Research Institute , Nanshan, Shenzhen 518057, China
| | | | - Xuhui Huang
- HKUST Shenzhen Research Institute , Nanshan, Shenzhen 518057, China
| |
Collapse
|
39
|
Gu S, Silva DA, Meng L, Yue A, Huang X. Quantitatively characterizing the ligand binding mechanisms of choline binding protein using Markov state model analysis. PLoS Comput Biol 2014; 10:e1003767. [PMID: 25101697 PMCID: PMC4125059 DOI: 10.1371/journal.pcbi.1003767] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2013] [Accepted: 06/22/2014] [Indexed: 01/05/2023] Open
Abstract
Protein-ligand recognition plays key roles in many biological processes. One of the most fascinating questions about protein-ligand recognition is to understand its underlying mechanism, which often results from a combination of induced fit and conformational selection. In this study, we have developed a three-pronged approach of Markov State Models, Molecular Dynamics simulations, and flux analysis to determine the contribution of each model. Using this approach, we have quantified the recognition mechanism of the choline binding protein (ChoX) to be ∼90% conformational selection dominant under experimental conditions. This is achieved by recovering all the necessary parameters for the flux analysis in combination with available experimental data. Our results also suggest that ChoX has several metastable conformational states, of which an apo-closed state is dominant, consistent with previous experimental findings. Our methodology holds great potential to be widely applied to understand recognition mechanisms underlining many fundamental biological processes.
Collapse
Affiliation(s)
- Shuo Gu
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Daniel-Adriano Silva
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
| | - Luming Meng
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Alexander Yue
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Division of Biomedical Engineering, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Center of Systems Biology and Human Health, Institute for Advance Study and School of Science, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- * E-mail:
| |
Collapse
|
40
|
Chodera JD, Noé F. Markov state models of biomolecular conformational dynamics. Curr Opin Struct Biol 2014; 25:135-44. [PMID: 24836551 DOI: 10.1016/j.sbi.2014.04.002] [Citation(s) in RCA: 493] [Impact Index Per Article: 49.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 04/08/2014] [Accepted: 04/12/2014] [Indexed: 10/25/2022]
Abstract
It has recently become practical to construct Markov state models (MSMs) that reproduce the long-time statistical conformational dynamics of biomolecules using data from molecular dynamics simulations. MSMs can predict both stationary and kinetic quantities on long timescales (e.g. milliseconds) using a set of atomistic molecular dynamics simulations that are individually much shorter, thus addressing the well-known sampling problem in molecular dynamics simulation. In addition to providing predictive quantitative models, MSMs greatly facilitate both the extraction of insight into biomolecular mechanism (such as folding and functional dynamics) and quantitative comparison with single-molecule and ensemble kinetics experiments. A variety of methodological advances and software packages now bring the construction of these models closer to routine practice. Here, we review recent progress in this field, considering theoretical and methodological advances, new software tools, and recent applications of these approaches in several domains of biochemistry and biophysics, commenting on remaining challenges.
Collapse
Affiliation(s)
- John D Chodera
- Computational Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA.
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany.
| |
Collapse
|
41
|
Bowman GR, Meng L, Huang X. Quantitative comparison of alternative methods for coarse-graining biological networks. J Chem Phys 2014; 139:121905. [PMID: 24089717 DOI: 10.1063/1.4812768] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Markov models and master equations are a powerful means of modeling dynamic processes like protein conformational changes. However, these models are often difficult to understand because of the enormous number of components and connections between them. Therefore, a variety of methods have been developed to facilitate understanding by coarse-graining these complex models. Here, we employ Bayesian model comparison to determine which of these coarse-graining methods provides the models that are most faithful to the original set of states. We find that the Bayesian agglomerative clustering engine and the hierarchical Nyström expansion graph (HNEG) typically provide the best performance. Surprisingly, the original Perron cluster cluster analysis (PCCA) method often provides the next best results, outperforming the newer PCCA+ method and the most probable paths algorithm. We also show that the differences between the models are qualitatively significant, rather than being minor shifts in the boundaries between states. The performance of the methods correlates well with the entropy of the resulting coarse-grainings, suggesting that finding states with more similar populations (i.e., avoiding low population states that may just be noise) gives better results.
Collapse
Affiliation(s)
- Gregory R Bowman
- Departments of Chemistry and Molecular and Cell Biology, University of California, Berkeley, California 94720, USA
| | | | | |
Collapse
|
42
|
Ferraris DM, Spallek R, Oehlmann W, Singh M, Rizzi M. Crystal structure of the Mycobacterium tuberculosis phosphate binding protein PstS3. Proteins 2014; 82:2268-74. [PMID: 24615888 DOI: 10.1002/prot.24548] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2013] [Revised: 02/20/2014] [Accepted: 03/04/2014] [Indexed: 11/10/2022]
Abstract
Mycobacterium tuberculosis evades host immune responses by colonizing macrophages. Intraphagosomal M. tuberculosis is exposed to environmental stresses such as reactive oxygen and nitrogen intermediates as well as acid shock and inorganic phosphate (Pi) depletion. Experimental evidence suggests that expression levels of mycobacterial protein PstS3 (Rv0928) are significantly increased when M. tuberculosis bacilli are exposed to Pi starvation. Hence, PstS3 may be important for survival of Mtb in conditions where there is limited supply of Pi. We report here the structure of PstS3 from M. tuberculosis at 2.3-Å resolution. The protein presents a structure typical for ABC phosphate transfer receptors. Comparison with its cognate receptor PstS1 showed a different pattern distribution of surface charges in proximity to the Pi recognition site, suggesting complementary roles of the two proteins in Pi uptake.
Collapse
Affiliation(s)
- Davide M Ferraris
- Department of Pharmaceutical Sciences, Università del Piemonte Orientale "A. Avogadro,", Largo Donegani 2, 28100, Novara, Italy
| | | | | | | | | |
Collapse
|
43
|
Bisignano P, Doerr S, Harvey MJ, Favia AD, Cavalli A, De Fabritiis G. Kinetic characterization of fragment binding in AmpC β-lactamase by high-throughput molecular simulations. J Chem Inf Model 2014; 54:362-6. [PMID: 24444037 DOI: 10.1021/ci4006063] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Small molecules used in fragment-based drug discovery form multiple, promiscuous binding complexes difficult to capture experimentally. Here, we identify such binding poses and their associated energetics and kinetics using molecular dynamics simulations on AmpC β-lactamase. Only one of the crystallographic binding poses was found to be thermodynamically favorable; however, the ligand shows several binding poses within the pocket. This study demonstrates free-binding molecular simulations in the context of fragment-to-lead development and its potential application in drug design.
Collapse
Affiliation(s)
- P Bisignano
- Department of Drug Discovery and Development, Istituto Italiano di Tecnologia , via Morego, 30, 16163 Genova, Italy
| | | | | | | | | | | |
Collapse
|
44
|
Takahashi R, Gil VA, Guallar V. Monte Carlo Free Ligand Diffusion with Markov State Model Analysis and Absolute Binding Free Energy Calculations. J Chem Theory Comput 2013; 10:282-8. [PMID: 26579911 DOI: 10.1021/ct400678g] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Obtaining absolute binding free energies from unbiased ligand diffusion has attracted a significant amount of attention due to its implications in drug design. Several studies have used special purpose computers and software to achieve microsecond molecular dynamics which, combined with a Markov state model analysis, are capable of providing absolute binding free energies. We have recently developed a Monte Carlo based technique, PELE, capable of performing a dynamical exploration of the protein-ligand energy landscape including free ligand diffusion into the active site, at a fraction of the computational cost of molecular dynamics techniques. We demonstrate here the capabilities of our Monte Carlo technique in obtaining absolute binding free energies for a series of benzamidine like inhibitors into trypsin. Our results are in good agreement with experimental data and other molecular dynamics simulations, indicating that PELE can be a useful tool for quick estimates of binding free energies and mechanisms.
Collapse
Affiliation(s)
- Ryoji Takahashi
- Joint BSC-IRB Research Program in Computational Biology, Barcelona Supercomputing Center , c/Jordi Girona 29, 08034 Barcelona, Spain
| | - Víctor A Gil
- Joint BSC-IRB Research Program in Computational Biology, Barcelona Supercomputing Center , c/Jordi Girona 29, 08034 Barcelona, Spain
| | - Victor Guallar
- Joint BSC-IRB Research Program in Computational Biology, Barcelona Supercomputing Center , c/Jordi Girona 29, 08034 Barcelona, Spain.,Institució Catalana de Recerca i Estudis Avançats (ICREA) , Passeig Lluís Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
45
|
Yao Y, Cui RZ, Bowman GR, Silva DA, Sun J, Huang X. Hierarchical Nyström methods for constructing Markov state models for conformational dynamics. J Chem Phys 2013; 138:174106. [PMID: 23656113 DOI: 10.1063/1.4802007] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Markov state models (MSMs) have become a popular approach for investigating the conformational dynamics of proteins and other biomolecules. MSMs are typically built from numerous molecular dynamics simulations by dividing the sampled configurations into a large number of microstates based on geometric criteria. The resulting microstate model can then be coarse-grained into a more understandable macrostate model by lumping together rapidly mixing microstates into larger, metastable aggregates. However, finite sampling often results in the creation of many poorly sampled microstates. During coarse-graining, these states are mistakenly identified as being kinetically important because transitions to/from them appear to be slow. In this paper, we propose a formalism based on an algebraic principle for matrix approximation, i.e., the Nyström method, to deal with such poorly sampled microstates. Our scheme builds a hierarchy of microstates from high to low populations and progressively applies spectral clustering on sets of microstates within each level of the hierarchy. It helps spectral clustering identify metastable aggregates with highly populated microstates rather than being distracted by lowly populated states. We demonstrate the ability of this algorithm to discover the major metastable states on two model systems, the alanine dipeptide and trpzip2 peptide.
Collapse
Affiliation(s)
- Yuan Yao
- School of Mathematical Sciences, LMAM-LMEQF-LMPR, Peking University, Beijing 100871, China.
| | | | | | | | | | | |
Collapse
|
46
|
Lindner B, Yi Z, Prinz JH, Smith JC, Noé F. Dynamic neutron scattering from conformational dynamics. I. Theory and Markov models. J Chem Phys 2013; 139:175101. [DOI: 10.1063/1.4824070] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
|
47
|
Trendelkamp-Schroer B, Noé F. Efficient Bayesian estimation of Markov model transition matrices with given stationary distribution. J Chem Phys 2013; 138:164113. [DOI: 10.1063/1.4801325] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
|
48
|
Bai F, Xu Y, Chen J, Liu Q, Gu J, Wang X, Ma J, Li H, Onuchic JN, Jiang H. Free energy landscape for the binding process of Huperzine A to acetylcholinesterase. Proc Natl Acad Sci U S A 2013; 110:4273-8. [PMID: 23440190 PMCID: PMC3600462 DOI: 10.1073/pnas.1301814110] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
Drug-target residence time (t = 1/k(off), where k(off) is the dissociation rate constant) has become an important index in discovering better- or best-in-class drugs. However, little effort has been dedicated to developing computational methods that can accurately predict this kinetic parameter or related parameters, k(off) and activation free energy of dissociation (ΔG(off)≠). In this paper, energy landscape theory that has been developed to understand protein folding and function is extended to develop a generally applicable computational framework that is able to construct a complete ligand-target binding free energy landscape. This enables both the binding affinity and the binding kinetics to be accurately estimated. We applied this method to simulate the binding event of the anti-Alzheimer's disease drug (-)-Huperzine A to its target acetylcholinesterase (AChE). The computational results are in excellent agreement with our concurrent experimental measurements. All of the predicted values of binding free energy and activation free energies of association and dissociation deviate from the experimental data only by less than 1 kcal/mol. The method also provides atomic resolution information for the (-)-Huperzine A binding pathway, which may be useful in designing more potent AChE inhibitors. We expect this methodology to be widely applicable to drug discovery and development.
Collapse
Affiliation(s)
- Fang Bai
- Department of Engineering Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, and
- Faculty of Chemical, Environmental, and Biological Science and Technology, Dalian University of Technology, Dalian 116023, China
| | - Yechun Xu
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Jing Chen
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Qiufeng Liu
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Junfeng Gu
- Department of Engineering Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, and
| | - Xicheng Wang
- Department of Engineering Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, and
| | - Jianpeng Ma
- Department of Bioengineering, and
- Verna and Marrs McLean Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, TX 77030; and
| | - Honglin Li
- State Key Laboratory of Bioreactor Engineering and Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai 200237, China
| | - José N. Onuchic
- Center for Theoretical Biological Physics and Department of Physics, Rice University, Houston, TX 77005
| | - Hualiang Jiang
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
- State Key Laboratory of Bioreactor Engineering and Shanghai Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, Shanghai 200237, China
| |
Collapse
|
49
|
Peuker S, Cukkemane A, Held M, Noé F, Kaupp UB, Seifert R. Kinetics of ligand-receptor interaction reveals an induced-fit mode of binding in a cyclic nucleotide-activated protein. Biophys J 2013; 104:63-74. [PMID: 23332059 DOI: 10.1016/j.bpj.2012.11.3816] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2012] [Revised: 11/14/2012] [Accepted: 11/27/2012] [Indexed: 12/17/2022] Open
Abstract
Many receptors and ion channels are activated by ligands. One key question concerns the binding mechanism. Does the ligand induce conformational changes in the protein via the induced-fit mechanism? Or does the protein preexist as an ensemble of conformers and the ligand selects the most complementary one, via the conformational selection mechanism? Here, we study ligand binding of a tetrameric cyclic nucleotide-gated channel from Mesorhizobium loti and of its monomeric binding domain (CNBD) using rapid mixing, mutagenesis, and structure-based computational biology. Association rate constants of ∼10(7) M(-1) s(-1) are compatible with diffusion-limited binding. Ligand binding to the full-length CNG channel and the isolated CNBD differ, revealing allosteric control of the CNBD by the effector domain. Finally, mutagenesis of allosteric residues affects only the dissociation rate constant, suggesting that binding follows the induced-fit mechanism. This study illustrates the strength of combining mutational, kinetic, and computational approaches to unravel important mechanistic features of ligand binding.
Collapse
Affiliation(s)
- Sebastian Peuker
- Department of Molecular Sensory Systems, Center of Advanced European Studies and Research, Bonn, Germany
| | | | | | | | | | | |
Collapse
|
50
|
Held M, Imhof P, Keller BG, Noé F. Modulation of a ligand's energy landscape and kinetics by the chemical environment. J Phys Chem B 2012; 116:13597-607. [PMID: 23025812 DOI: 10.1021/jp3006684] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Understanding how the chemical environment modulates the predominant conformations and kinetics of flexible molecules is a core interest of biochemistry and a prerequisite for the rational design of synthetic catalysts. This study combines molecular dynamics simulation and Markov state models (MSMs) to a systematic computational strategy for investigating the effect of the chemical environment of a molecule on its conformations and kinetics. MSMs allow quantities to be computed that are otherwise difficult to access, such as the metastable sets, their free energies, and the relaxation time scales related to the rare transitions between metastable states. Additionally, MSMs are useful to identify observables that may act as sensors for the conformational or binding state of the molecule, thus guiding the design of experiments. In the present study, the conformation dynamics of UDP-GlcNAc are studied in vacuum, water, water + Mg(2+), and in the protein UDP-GlcNAc 2-epimerase. It is found that addition of Mg(2+) significantly affects the conformational stability, thermodynamics, and kinetics of UDP-GlcNAc. In particular, the slowest structural process, puckering of the GlcNAc sugar, depends on the overall conformation of UDP-GlcNAc and may thus act as a sensor of whether Mg(2+) is bound or not. Interestingly, transferring the molecule from vacuum to water makes the protein-binding conformations UDP-GlcNAc first accessible, while adding Mg(2+) further stabilizes them by specifically associating to binding-competent conformations. While Mg(2+) is not cocrystallized in the UDP-GlcNAc 2-epimerase complex, the selectively stabilized Mg(2+)/UDP-GlcNAc complex may be a template for the bound state, and Mg(2+) may accompany the binding-competent ligand conformation to the binding pocket. This serves as a possible explanation of the enhanced epimerization rate in the presence of Mg(2+). This role of Mg(2+) has previously not been described and opens the question whether "binding co-factors" may be a concept of general relevance for protein-ligand binding.
Collapse
Affiliation(s)
- Martin Held
- Institute of Mathematics, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | | | | | | |
Collapse
|