1
|
Joshi DC, Gosse C, Huang SY, Lin JH. A Curvilinear-Path Umbrella Sampling Approach to Characterizing the Interactions Between Rapamycin and Three FKBP12 Variants. Front Mol Biosci 2022; 9:879000. [PMID: 35874613 PMCID: PMC9304761 DOI: 10.3389/fmolb.2022.879000] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 05/12/2022] [Indexed: 11/13/2022] Open
Abstract
Rapamycin is an immunosuppressant macrolide that exhibits anti-proliferative properties through inhibiting the mTOR kinase. In fact, the drug first associates with the FKBP12 enzyme before interacting with the FRB domain of its target. Despite the availability of structural and thermodynamic information on the interaction of FKBP12 with rapamycin, the energetic and mechanistic understanding of this process is still incomplete. We recently reported a multiple-walker umbrella sampling simulation approach to characterizing the protein–protein interaction energetics along curvilinear paths. In the present paper, we extend our investigations to a protein-small molecule duo, the FKBP12•rapamycin complex. We estimate the binding free energies of rapamycin with wild-type FKBP12 and two mutants in which a hydrogen bond has been removed, D37V and Y82F. Furthermore, the underlying mechanistic details are analyzed. The calculated standard free energies of binding agree well with the experimental data, and the roles of the hydrogen bonds are shown to be quite different for each of these two mutated residues. On one hand, removing the carboxylate group of D37 strongly destabilizes the association; on the other hand, the hydroxyl group of Y82 is nearly unnecessary for the stability of the complex because some nonconventional, cryptic, indirect interaction mechanisms seem to be at work.
Collapse
Affiliation(s)
| | - Charlie Gosse
- Institut de Biologie de l’Ecole Normale Supérieure, ENS, CNRS, INSERM, PSL Research University, Paris, France
| | - Shu-Yu Huang
- Research Center for Applied Sciences, Academia Sinica, Taipei, Taiwan
| | - Jung-Hsin Lin
- Research Center for Applied Sciences, Academia Sinica, Taipei, Taiwan
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
- Biomedical Translation Research Center, National Biotechnology Research Park, Academia Sinica, Taipei, Taiwan
- School of Pharmacy, College of Medicine, National Taiwan University, Taipei, Taiwan
- College of Engineering Sciences, Chang Gung University, Taoyuan, Taiwan
- *Correspondence: Jung-Hsin Lin,
| |
Collapse
|
2
|
Solomentsev G, Diehl C, Akke M. Conformational Entropy of FK506 Binding to FKBP12 Determined by Nuclear Magnetic Resonance Relaxation and Molecular Dynamics Simulations. Biochemistry 2018; 57:1451-1461. [DOI: 10.1021/acs.biochem.7b01256] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Gleb Solomentsev
- Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Carl Diehl
- Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mikael Akke
- Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
3
|
Weininger U, Modig K, Geitner AJ, Schmidpeter PAM, Koch JR, Akke M. Dynamics of Aromatic Side Chains in the Active Site of FKBP12. Biochemistry 2016; 56:334-343. [DOI: 10.1021/acs.biochem.6b01157] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Ulrich Weininger
- Department
of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O.
Box 124, SE-22100 Lund, Sweden
- Institute
of Physics, Biophysics, Martin-Luther-University Halle-Wittenberg, D-06120 Halle (Saale), Germany
| | - Kristofer Modig
- Department
of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O.
Box 124, SE-22100 Lund, Sweden
| | - Anne-Juliane Geitner
- Laboratorium
für Biochemie, Bayreuther Zentrum für Molekulare Biowissenschaften, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Philipp A. M. Schmidpeter
- Laboratorium
für Biochemie, Bayreuther Zentrum für Molekulare Biowissenschaften, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Johanna R. Koch
- Laboratorium
für Biochemie, Bayreuther Zentrum für Molekulare Biowissenschaften, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Mikael Akke
- Department
of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, P.O.
Box 124, SE-22100 Lund, Sweden
| |
Collapse
|
4
|
Structure-affinity properties of a high-affinity ligand of FKBP12 studied by molecular simulations of a binding intermediate. PLoS One 2014; 9:e114610. [PMID: 25502559 PMCID: PMC4264844 DOI: 10.1371/journal.pone.0114610] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2014] [Accepted: 11/10/2014] [Indexed: 11/19/2022] Open
Abstract
With a view to explaining the structure-affinity properties of the ligands of the protein FKBP12, we characterized a binding intermediate state between this protein and a high-affinity ligand. Indeed, the nature and extent of the intermolecular contacts developed in such a species may play a role on its stability and, hence, on the overall association rate. To find the binding intermediate, a molecular simulation protocol was used to unbind the ligand by gradually decreasing the biasing forces introduced. The intermediate was subsequently refined with 17 independent stochastic boundary molecular dynamics simulations that provide a consistent picture of the intermediate state. In this state, the core region of the ligand remains stable, notably because of the two anchoring oxygen atoms that correspond to recurrent motifs found in all FKBP12 ligand core structures. Besides, the non-core regions participate in numerous transient intermolecular and intramolecular contacts. The dynamic aspect of most of the contacts seems important both for the ligand to retain at least a part of its configurational entropy and for avoiding a trapped state along the binding pathway. Since the transient and anchoring contacts contribute to increasing the stability of the intermediate, as a corollary, the dissociation rate constant of this intermediate should be decreased, resulting in an increase of the affinity constant . The present results support our previous conclusions and provide a coherent rationale for explaining the prevalence in high-affinity ligands of (i) the two oxygen atoms found in carbonyl or sulfonyl groups of dissimilar core structures and of (ii) symmetric or pseudo-symmetric mobile groups of atoms found as non-core moieties. Another interesting aspect of the intermediate is the distortion of the flexible 80 s loop of the protein, mainly in its tip region, that promotes the accessibility to the bound state.
Collapse
|
5
|
Olivieri L, Gardebien F. Classical force field parameters for two high-affinity ligands of FKBP12. J Mol Graph Model 2014; 49:118-28. [PMID: 24657432 DOI: 10.1016/j.jmgm.2014.02.003] [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] [Received: 06/17/2013] [Revised: 02/16/2014] [Accepted: 02/18/2014] [Indexed: 11/25/2022]
Abstract
FKBP12 is an important target in the treatment of transplant rejection and is also a promising target for cancer and neurodegenerative diseases. We determined for two ligands of nanomolar affinity the set of parameters in the CHARMM force field. The fitting procedure was based on reproducing the quantum chemistry data (distances, angles, and energies). Since the dynamical behavior of such ligands strongly depends on the dihedral angles, care was taken to derive the corresponding parameters. Moreover, since each of the central core region of these two ligands is similar to other known ligands or drugs of other proteins, part at least of these parameters could also be useful for these other ligands.
Collapse
Affiliation(s)
- Lilian Olivieri
- DSIMB, INSERM, U1134, Paris F-75015, France; Université de la Réunion, UMR_S 1134, Faculté des Sciences et Technologies, 15, avenue René Cassin, BP 7151, 97715 Saint Denis Messag Cedex 09, Réunion; Institut National de la Transfusion Sanguine, F-75015 Paris, France; Laboratory of Excellence GR-Ex
| | - Fabrice Gardebien
- DSIMB, INSERM, U1134, Paris F-75015, France; Université de la Réunion, UMR_S 1134, Faculté des Sciences et Technologies, 15, avenue René Cassin, BP 7151, 97715 Saint Denis Messag Cedex 09, Réunion; Institut National de la Transfusion Sanguine, F-75015 Paris, France; Laboratory of Excellence GR-Ex.
| |
Collapse
|
6
|
Lapelosa M, Gallicchio E, Levy RM. Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. J Chem Theory Comput 2011; 8:47-60. [PMID: 22368530 DOI: 10.1021/ct200684b] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The Binding Energy Distribution Analysis Method (BEDAM) is employed to compute the standard binding free energies of a series of ligands to a FK506 binding protein (FKBP12) with implicit solvation. Binding free energy estimates are in reasonably good agreement with experimental affinities. The conformations of the complexes identified by the simulations are in good agreement with crystallographic data, which was not used to restrain ligand orientations. The BEDAM method is based on λ -hopping Hamiltonian parallel Replica Exchange (HREM) molecular dynamics conformational sampling, the OPLS-AA/AGBNP2 effective potential, and multi-state free energy estimators (MBAR). Achieving converged and accurate results depends on all of these elements of the calculation. Convergence of the binding free energy is tied to the level of convergence of binding energy distributions at critical intermediate states where bound and unbound states are at equilibrium, and where the rate of binding/unbinding conformational transitions is maximal. This finding mirrors similar observations in the context of order/disorder transitions as for example in protein folding. Insights concerning the physical mechanism of ligand binding and unbinding are obtained. Convergence for the largest FK506 ligand is achieved only after imposing strict conformational restraints, which however require accurate prior structural knowledge of the structure of the complex. The analytical AGBNP2 model is found to underestimate the magnitude of the hydrophobic driving force towards binding in these systems characterized by loosely packed protein-ligand binding interfaces. Rescoring of the binding energies using a numerical surface area model corrects this deficiency. This study illustrates the complex interplay between energy models, exploration of conformational space, and free energy estimators needed to obtain robust estimates from binding free energy calculations.
Collapse
Affiliation(s)
- Mauro Lapelosa
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers the State University of New Jersey, Piscataway, NJ 08854
| | | | | |
Collapse
|