1
|
Champion C, Gall R, Ries B, Rieder SR, Barros EP, Riniker S. Accelerating Alchemical Free Energy Prediction Using a Multistate Method: Application to Multiple Kinases. J Chem Inf Model 2023; 63:7133-7147. [PMID: 37948537 PMCID: PMC10685456 DOI: 10.1021/acs.jcim.3c01469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/23/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023]
Abstract
Alchemical free-energy methods based on molecular dynamics (MD) simulations have become important tools to identify modifications of small organic molecules that improve their protein binding affinity during lead optimization. The routine application of pairwise free-energy methods to rank potential binders from best to worst is impacted by the combinatorial increase in calculations to perform when the number of molecules to assess grows. To address this fundamental limitation, our group has developed replica-exchange enveloping distribution sampling (RE-EDS), a pathway-independent multistate method, enabling the calculation of alchemical free-energy differences between multiple ligands (N > 2) from a single MD simulation. In this work, we apply the method to a set of four kinases with diverse binding pockets and their corresponding inhibitors (42 in total), chosen to showcase the general applicability of RE-EDS in prospective drug design campaigns. We show that for the targets studied, RE-EDS is able to model up to 13 ligands simultaneously with high sampling efficiency, leading to a substantial decrease in computational cost when compared to pairwise methods.
Collapse
Affiliation(s)
- Candide Champion
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - René Gall
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | | | - Salomé R. Rieder
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Emilia P. Barros
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
2
|
Wade AD, Huggins DJ. Identification of Optimal Ligand Growth Vectors Using an Alchemical Free-Energy Method. J Chem Inf Model 2020; 60:5580-5594. [PMID: 32810401 DOI: 10.1021/acs.jcim.0c00610] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, a novel method to rationally design inhibitors with improved steric contacts and enhanced binding free energies is presented. This new method uses alchemical single step perturbation calculations to rapidly optimize the van der Waals interactions of a small molecule in a protein-ligand complex in order to maximize its binding affinity. The results of the optimizer are used to predict beneficial growth vectors on the ligand, and good agreement is found between the predictions from the optimizer and a more rigorous free energy calculation, with a Spearman's rank order correlation of 0.59. The advantage of the method presented here is the significant speed up of over 10-fold compared to traditional free energy calculations and sublinear scaling with the number of growth vectors assessed. Where experimental data were available, mutations from hydrogen to a methyl group at sites highlighted by the optimizer were calculated with MBAR, and the mean unsigned error between experimental and calculated values of the binding free energy was 0.83 kcal/mol.
Collapse
Affiliation(s)
- Alexander D Wade
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David J Huggins
- Tri-Institutional Therapeutics Discovery Institute, Belfer Research Building, 413 East 69th Street, 16th Floor, Box 300, New York, United States.,Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| |
Collapse
|
3
|
Perthold JW, Petrov D, Oostenbrink C. Toward Automated Free Energy Calculation with Accelerated Enveloping Distribution Sampling (A-EDS). J Chem Inf Model 2020; 60:5395-5406. [PMID: 32492343 PMCID: PMC7686955 DOI: 10.1021/acs.jcim.0c00456] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Free-energy
perturbation (FEP) methods are commonly used in drug
design to calculate relative binding free energies of different ligands
to a common host protein. Alchemical ligand transformations are usually
performed in multiple steps which need to be chosen carefully to ensure
sufficient phase-space overlap between neighboring states. With one-step
or single-step FEP techniques, a single reference state is designed
that samples phase-space not only representative of a full transformation
but also ideally resembles multiple ligand end states and hence allows
for efficient multistate perturbations. Enveloping distribution sampling
(EDS) is one example for such a method in which the reference state
is created by a mathematical combination of the different ligand end
states based on solid statistical mechanics. We have recently proposed
a novel approach to EDS which enables efficient barrier crossing between
the different end states, termed accelerated EDS (A-EDS). In this
work, we further simplify the parametrization of the A-EDS reference
state and demonstrate the automated calculation of multiple free-energy
differences between different ligands from a single simulation in
three different well-described drug design model systems.
Collapse
Affiliation(s)
- Jan Walther Perthold
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Dražen Petrov
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| |
Collapse
|
4
|
Li G, Zhang W, Xie Y, Li Y, Cao R, Zheng G, Huang N, Zhou Y. Structure-Based Optimization of 10-DEBC Derivatives as Potent and Selective Pim-1 Kinase Inhibitors. J Chem Inf Model 2020; 60:3287-3294. [PMID: 32407627 DOI: 10.1021/acs.jcim.0c00245] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Pim-1 kinase has been widely regarded as an attractive target for anticancer drugs. Here, we reported our continued efforts in structure-based optimization of compound 10-DEBC, a previously identified micromolar Pim-1 inhibitor. Guided by the Site Identification by Ligand Competitive Saturation (SILCS) method, we quickly obtained a series of 10-DEBC derivatives with significantly improved activity and selectivity. In particular, compound 26 exhibited an IC50 value of 0.9 nM, as well as 220- and 8-fold selectivity over Pim-2 and Pim-3 kinases, respectively.
Collapse
Affiliation(s)
- Gudong Li
- State Key Laboratory of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Wei Zhang
- National Institute of Biological Sciences, Beijing, Beijing 102206, China
| | - Yuting Xie
- National Institute of Biological Sciences, Beijing, Beijing 102206, China
| | - Yang Li
- National Institute of Biological Sciences, Beijing, Beijing 102206, China
| | - Rao Cao
- National Institute of Biological Sciences, Beijing, Beijing 102206, China
| | - Guojun Zheng
- State Key Laboratory of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Niu Huang
- National Institute of Biological Sciences, Beijing, Beijing 102206, China.,Tsinghua Institute of Multidisciplinary Biomedical Research, Tsinghua University, Beijing 102206, China
| | - Yu Zhou
- National Institute of Biological Sciences, Beijing, Beijing 102206, China
| |
Collapse
|
5
|
Miao M, Fu H, Zhang H, Shao X, Chipot C, Cai W. Avoiding non-equilibrium effects in adaptive biasing force calculations. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1775222] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Mengyao Miao
- Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin, People’s Republic of China
| | - Haohao Fu
- Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin, People’s Republic of China
| | - Hong Zhang
- Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin, People’s Republic of China
| | - Xueguang Shao
- Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin, People’s Republic of China
- State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin, People’s Republic of China
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana−Champaign, Vandœuvre-lès-Nancy cedex, France
- Department of Physics, University of Illinois at Urbana−Champaign, 1110 West Green Street, Urbana, IL 61801, USA
| | - Wensheng Cai
- Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin, People’s Republic of China
| |
Collapse
|
6
|
Wade AD, Huggins DJ. Optimization of Protein-Ligand Electrostatic Interactions Using an Alchemical Free-Energy Method. J Chem Theory Comput 2019; 15:6504-6512. [PMID: 31584802 PMCID: PMC7007198 DOI: 10.1021/acs.jctc.9b00976] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
![]()
We present an explicit solvent alchemical free-energy method for
optimizing the partial charges of a ligand to maximize the binding
affinity with a receptor. This methodology can be applied to known
ligand–protein complexes to determine an optimized set of ligand
partial atomic changes. Three protein–ligand complexes have
been optimized in this work: FXa, P38, and the androgen receptor.
The sets of optimized charges can be used to identify design principles
for chemical changes to the ligands which improve the binding affinity
for all three systems. In this work, beneficial chemical mutations
are generated from these principles and the resulting molecules tested
using free-energy perturbation calculations. We show that three quarters
of our chemical changes are predicted to improve the binding affinity,
with an average improvement for the beneficial mutations of approximately
1 kcal/mol. In the cases where experimental data are available, the
agreement between prediction and experiment is also good. The results
demonstrate that charge optimization in explicit solvent is a useful
tool for predicting beneficial chemical changes such as pyridinations,
fluorinations, and oxygen to sulfur mutations.
Collapse
Affiliation(s)
- Alexander D Wade
- TCM Group, Cavendish Laboratory , University of Cambridge , 19 J J Thomson Avenue , Cambridge CB3 0HE , U.K
| | - David J Huggins
- TCM Group, Cavendish Laboratory , University of Cambridge , 19 J J Thomson Avenue , Cambridge CB3 0HE , U.K.,Tri-Institutional Therapeutics Discovery Institute , Belfer Research Building, 413 East 69th Street, 16th Floor, Box 300 , New York 10021 , United States.,Department of Physiology and Biophysics , Weill Cornell Medicine , 1300 York Avenue , New York , New York 10065 , United States
| |
Collapse
|
7
|
Song LF, Lee TS, Chun-Zhu, York DM, Merz KM. Using AMBER18 for Relative Free Energy Calculations. J Chem Inf Model 2019; 59:3128-3135. [PMID: 31244091 PMCID: PMC7371000 DOI: 10.1021/acs.jcim.9b00105] [Citation(s) in RCA: 129] [Impact Index Per Article: 25.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
With renewed interest in free energy methods in contemporary structure-based drug design, there is a pressing need to validate against multiple targets and force fields to assess the overall ability of these methods to accurately predict relative binding free energies. We computed relative binding free energies using graphics processing unit accelerated thermodynamic integration (GPU-TI) on a data set originally assembled by Schrödinger, Inc. Using their GPU free energy code (FEP+) and the OPLS2.1 force field combined with the REST2 enhanced sampling approach, these authors obtained an overall MUE of 0.9 kcal/mol and an overall RMSD of 1.14 kcal/mol. In our study using GPU-TI from AMBER with the AMBER14SB/GAFF1.8 force field but without enhanced sampling, we obtained an overall MUE of 1.17 kcal/mol and an overall RMSD of 1.50 kcal/mol for the 330 perturbations contained in this data set. A more detailed analysis of our results suggested that the observed differences between the two studies arise from differences in sampling protocols along with differences in the force fields employed. Future work should address the problem of establishing benchmark quality results with robust statistical error bars obtained through multiple independent runs and enhanced sampling, which is possible with the GPU-accelerated features in AMBER.
Collapse
Affiliation(s)
- Lin Frank Song
- Department of Chemistry and the Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, MI 48824
| | - Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Chun-Zhu
- Department of Chemistry and the Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, MI 48824
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Kenneth M. Merz
- Department of Chemistry and the Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, MI 48824
- Institute for Cyber Enabled Research, Michigan State University, 567 Wilson Road, Room 1440, East Lansing, MI 48824
| |
Collapse
|
8
|
Wade AD, Rizzi A, Wang Y, Huggins DJ. Computational Fluorine Scanning Using Free-Energy Perturbation. J Chem Inf Model 2019; 59:2776-2784. [PMID: 31046267 DOI: 10.1021/acs.jcim.9b00228] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
We present perturbative fluorine scanning, a computational fluorine scanning approach using free-energy perturbation. This method can be applied to molecular dynamics simulations of a single compound and make predictions for the best binders out of numerous fluorinated analogues. We tested the method on nine test systems: renin, DPP4, menin, P38, factor Xa, CDK2, AKT, JAK2, and androgen receptor. The predictions were in excellent agreement with more rigorous alchemical free-energy calculations and in good agreement with experimental data for most of the test systems. However, the agreement with experiment was very poor in some of the test systems, and this highlights the need for improved force fields in addition to accurate treatment of tautomeric and protonation states. The method is of particular interest due to the wide use of fluorine in medicinal chemistry to improve binding affinity and ADME properties. The promising results on this test case suggest that perturbative fluorine scanning will be a useful addition to the available arsenal of free-energy methods.
Collapse
Affiliation(s)
- Alexander D Wade
- TCM Group, Cavendish Laboratory , University of Cambridge , 19 J J Thomson Avenue , Cambridge CB3 0HE , United Kingdom
| | - Andrea Rizzi
- Tri-Institutional Training Program in Computational Biology and Medicine , New York , New York 10065 , United States.,Computational and Systems Biology Program , Sloan Kettering Institute, Memorial Sloan-Kettering Cancer Center , New York , New York 10065 , United States
| | - Yuanqing Wang
- Physiology, Biophysics, and System Biology Program , Weill Cornell Medicine , 1300 York Avenue , New York , New York 10065 , United States
| | - David J Huggins
- TCM Group, Cavendish Laboratory , University of Cambridge , 19 J J Thomson Avenue , Cambridge CB3 0HE , United Kingdom.,Tri-Institutional Therapeutics Discovery Institute , Belfer Research Building , 413 East 69th Street, 16th Floor , Box 300, New York , New York 10021 , United States.,Department of Physiology and Biophysics , Weill Cornell Medicine , 1300 York Avenue , New York , New York 10065 , United States
| |
Collapse
|
9
|
Hahn DF, Hünenberger PH. Alchemical Free-Energy Calculations by Multiple-Replica λ-Dynamics: The Conveyor Belt Thermodynamic Integration Scheme. J Chem Theory Comput 2019; 15:2392-2419. [PMID: 30821973 DOI: 10.1021/acs.jctc.8b00782] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
A new method is proposed to calculate alchemical free-energy differences based on molecular dynamics (MD) simulations, called the conveyor belt thermodynamic integration (CBTI) scheme. As in thermodynamic integration (TI), K replicas of the system are simulated at different values of the alchemical coupling parameter λ. The number K is taken to be even, and the replicas are equally spaced on a forward-turn-backward-turn path, akin to a conveyor belt (CB) between the two physical end-states; and as in λ-dynamics (λD), the λ-values associated with the individual systems evolve in time along the simulation. However, they do so in a concerted fashion, determined by the evolution of a single dynamical variable Λ of period 2π controlling the advance of the entire CB. Thus, a change of Λ is always associated with K/2 equispaced replicas moving forward and K/2 equispaced replicas moving backward along λ. As a result, the effective free-energy profile of the replica system along Λ is periodic of period 2 πK-1, and the magnitude of its variations decreases rapidly upon increasing K, at least as K-1 in the limit of large K. When a sufficient number of replicas is used, these variations become small, which enables a complete and quasi-homogeneous coverage of the λ-range by the replica system, without application of any biasing potential. If desired, a memory-based biasing potential can still be added to further homogenize the sampling, the preoptimization of which is computationally inexpensive. The final free-energy profile along λ is calculated similarly to TI, by binning of the Hamiltonian λ-derivative as a function of λ considering all replicas simultaneously, followed by quadrature integration. The associated quadrature error can be kept very low owing to the continuous and quasi-homogeneous λ-sampling. The CBTI scheme can be viewed as a continuous/deterministic/dynamical analog of the Hamiltonian replica-exchange/permutation (HRE/HRP) schemes or as a correlated multiple-replica analog of the λD or λ-local elevation umbrella sampling (λ-LEUS) schemes. Compared to TI, it shares the advantage of the latter schemes in terms of enhanced orthogonal sampling, i.e. the availability of variable-λ paths to circumvent conformational barriers present at specific λ-values. Compared to HRE/HRP, it permits a deterministic and continuous sampling of the λ-range, is expected to be less sensitive to possible artifacts of the thermo- and barostating schemes, and bypasses the need to carefully preselect a λ-ladder and a swapping-attempt frequency. Compared to λ-LEUS, it eliminates (or drastically reduces) the dead time associated with the preoptimization of a biasing potential. The goal of this article is to provide the mathematical/physical formulation of the proposed CBTI scheme, along with an initial application of the method to the calculation of the hydration free energy of methanol.
Collapse
Affiliation(s)
- David F Hahn
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Philippe H Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| |
Collapse
|
10
|
Perthold JW, Oostenbrink C. Accelerated Enveloping Distribution Sampling: Enabling Sampling of Multiple End States while Preserving Local Energy Minima. J Phys Chem B 2018; 122:5030-5037. [DOI: 10.1021/acs.jpcb.8b02725] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Walther Perthold
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, Department for Material Sciences and Process Engineering, University of Natural Resources and Life Sciences (BOKU), Vienna, Muthgasse 18, 1190 Vienna, Austria
| |
Collapse
|
11
|
Jandova Z, Fast D, Setz M, Pechlaner M, Oostenbrink C. Saturation Mutagenesis by Efficient Free-Energy Calculation. J Chem Theory Comput 2018; 14:894-904. [PMID: 29262673 PMCID: PMC5813279 DOI: 10.1021/acs.jctc.7b01099] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Single-point mutations
in proteins can greatly influence protein
stability, binding affinity, protein function or its expression per
se. Here, we present accurate and efficient predictions of the free
energy of mutation of amino acids. We divided the complete mutational
free energy into an uncharging step, which we approximate by a third-power
fitting (TPF) approach, and an annihilation step, which we approximate
using the one-step perturbation (OSP) method. As a diverse set of
test systems, we computed the solvation free energy of all amino acid
side chain analogues and obtained an excellent agreement with thermodynamic
integration (TI) data. Moreover, we calculated mutational free energies
in model tripeptides and established an efficient protocol involving
a single reference state. Again, the approximate methods agreed excellently
with the TI references, with a root-mean-square error of only 3.6
kJ/mol over 17 mutations. Our combined TPF+OSP approach does show
not only a very good agreement but also a 2-fold higher efficiency
than full blown TI calculations.
Collapse
Affiliation(s)
- Zuzana Jandova
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Daniel Fast
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Martina Setz
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Maria Pechlaner
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| |
Collapse
|
12
|
Nguyen TH, Zhou HX, Minh DDL. Using the fast fourier transform in binding free energy calculations. J Comput Chem 2017; 39:621-636. [PMID: 29270990 DOI: 10.1002/jcc.25139] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2017] [Revised: 09/23/2017] [Accepted: 11/27/2017] [Indexed: 12/21/2022]
Abstract
According to implicit ligand theory, the standard binding free energy is an exponential average of the binding potential of mean force (BPMF), an exponential average of the interaction energy between the unbound ligand ensemble and a rigid receptor. Here, we use the fast Fourier transform (FFT) to efficiently evaluate BPMFs by calculating interaction energies when rigid ligand configurations from the unbound ensemble are discretely translated across rigid receptor conformations. Results for standard binding free energies between T4 lysozyme and 141 small organic molecules are in good agreement with previous alchemical calculations based on (1) a flexible complex ( R≈0.9 for 24 systems) and (2) flexible ligand with multiple rigid receptor configurations ( R≈0.8 for 141 systems). While the FFT is routinely used for molecular docking, to our knowledge this is the first time that the algorithm has been used for rigorous binding free energy calculations. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Trung Hai Nguyen
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois, 60616
| | - Huan-Xiang Zhou
- Departments of Chemistry and Physics, University of Illinois at Chicago, Chicago, Illinois, 60607
| | - David D L Minh
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois, 60616
| |
Collapse
|
13
|
Cournia Z, Allen B, Sherman W. Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations. J Chem Inf Model 2017; 57:2911-2937. [PMID: 29243483 DOI: 10.1021/acs.jcim.7b00564] [Citation(s) in RCA: 395] [Impact Index Per Article: 56.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.
Collapse
Affiliation(s)
- Zoe Cournia
- Biomedical Research Foundation, Academy of Athens , 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Bryce Allen
- Silicon Therapeutics , 300 A Street, Boston, Massachusetts 02210, United States
| | - Woody Sherman
- Silicon Therapeutics , 300 A Street, Boston, Massachusetts 02210, United States
| |
Collapse
|
14
|
Raman EP, Lakkaraju SK, Denny RA, MacKerell AD. Estimation of relative free energies of binding using pre-computed ensembles based on the single-step free energy perturbation and the site-identification by Ligand competitive saturation approaches. J Comput Chem 2017; 38:1238-1251. [PMID: 27782307 PMCID: PMC5403604 DOI: 10.1002/jcc.24522] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2016] [Revised: 09/17/2016] [Accepted: 10/04/2016] [Indexed: 12/19/2022]
Abstract
Accurate and rapid estimation of relative binding affinities of ligand-protein complexes is a requirement of computational methods for their effective use in rational ligand design. Of the approaches commonly used, free energy perturbation (FEP) methods are considered one of the most accurate, although they require significant computational resources. Accordingly, it is desirable to have alternative methods of similar accuracy but greater computational efficiency to facilitate ligand design. In the present study relative free energies of binding are estimated for one or two non-hydrogen atom changes in compounds targeting the proteins ACK1 and p38 MAP kinase using three methods. The methods include standard FEP, single-step free energy perturbation (SSFEP) and the site-identification by ligand competitive saturation (SILCS) ligand grid free energy (LGFE) approach. Results show the SSFEP and SILCS LGFE methods to be competitive with or better than the FEP results for the studied systems, with SILCS LGFE giving the best agreement with experimental results. This is supported by additional comparisons with published FEP data on p38 MAP kinase inhibitors. While both the SSFEP and SILCS LGFE approaches require a significant upfront computational investment, they offer a 1000-fold computational savings over FEP for calculating the relative affinities of ligand modifications once those pre-computations are complete. An illustrative example of the potential application of these methods in the context of screening large numbers of transformations is presented. Thus, the SSFEP and SILCS LGFE approaches represent viable alternatives for actively driving ligand design during drug discovery and development. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- E. Prabhu Raman
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn Street HSF II, Baltimore MD 21201
| | - Sirish Kaushik Lakkaraju
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn Street HSF II, Baltimore MD 21201
| | - Rajiah Aldrin Denny
- Medicine Design, Worldwide Research & Development, Pfizer Inc, 610 Main Street, Cambridge, MA 02139, USA
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn Street HSF II, Baltimore MD 21201
| |
Collapse
|
15
|
Abstract
Computational approaches are useful tools to interpret and guide experiments to expedite the antibiotic drug design process. Structure-based drug design (SBDD) and ligand-based drug design (LBDD) are the two general types of computer-aided drug design (CADD) approaches in existence. SBDD methods analyze macromolecular target 3-dimensional structural information, typically of proteins or RNA, to identify key sites and interactions that are important for their respective biological functions. Such information can then be utilized to design antibiotic drugs that can compete with essential interactions involving the target and thus interrupt the biological pathways essential for survival of the microorganism(s). LBDD methods focus on known antibiotic ligands for a target to establish a relationship between their physiochemical properties and antibiotic activities, referred to as a structure-activity relationship (SAR), information that can be used for optimization of known drugs or guide the design of new drugs with improved activity. In this chapter, standard CADD protocols for both SBDD and LBDD will be presented with a special focus on methodologies and targets routinely studied in our laboratory for antibiotic drug discoveries.
Collapse
|
16
|
Chiang YC, Pang YT, Wang Y. The role of intramolecular nonbonded interaction and angle sampling in single-step free energy perturbation. J Chem Phys 2016; 145:234109. [PMID: 28010084 DOI: 10.1063/1.4972230] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Single-step free energy perturbation (sFEP) has often been proposed as an efficient tool for a quick free energy scan due to its straightforward protocol and the ability to recycle an existing molecular dynamics trajectory for free energy calculations. Although sFEP is expected to fail when the sampling of a system is inefficient, it is often expected to hold for an alchemical transformation between ligands with a moderate difference in their sizes, e.g., transforming a benzene into an ethylbenzene. Yet, exceptions were observed in calculations for anisole and methylaniline, which have similar physical sizes as ethylbenzene. In this study, we show that such exceptions arise from the sampling inefficiency on an unexpected rigid degree of freedom, namely, the bond angle θ. The distributions of θ differ dramatically between two end states of a sFEP calculation, i.e., the conformation of the ligand changes significantly during the alchemical transformation process. Our investigation also reveals the interrelation between the ligand conformation and the intramolecular nonbonded interactions. This knowledge suggests a best combination of the ghost ligand potential and the dual topology setting, which improves the accuracy in a single reference sFEP calculation by bringing down its error from around 5kBT to kBT.
Collapse
Affiliation(s)
- Ying-Chih Chiang
- Department of Physics, Chinese University of Hong Kong, Shatin, NT, Hong Kong
| | - Yui Tik Pang
- Department of Physics, Chinese University of Hong Kong, Shatin, NT, Hong Kong
| | - Yi Wang
- Department of Physics, Chinese University of Hong Kong, Shatin, NT, Hong Kong
| |
Collapse
|
17
|
Leelananda SP, Lindert S. Computational methods in drug discovery. Beilstein J Org Chem 2016; 12:2694-2718. [PMID: 28144341 PMCID: PMC5238551 DOI: 10.3762/bjoc.12.267] [Citation(s) in RCA: 285] [Impact Index Per Article: 35.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2016] [Accepted: 11/22/2016] [Indexed: 12/11/2022] Open
Abstract
The process for drug discovery and development is challenging, time consuming and expensive. Computer-aided drug discovery (CADD) tools can act as a virtual shortcut, assisting in the expedition of this long process and potentially reducing the cost of research and development. Today CADD has become an effective and indispensable tool in therapeutic development. The human genome project has made available a substantial amount of sequence data that can be used in various drug discovery projects. Additionally, increasing knowledge of biological structures, as well as increasing computer power have made it possible to use computational methods effectively in various phases of the drug discovery and development pipeline. The importance of in silico tools is greater than ever before and has advanced pharmaceutical research. Here we present an overview of computational methods used in different facets of drug discovery and highlight some of the recent successes. In this review, both structure-based and ligand-based drug discovery methods are discussed. Advances in virtual high-throughput screening, protein structure prediction methods, protein-ligand docking, pharmacophore modeling and QSAR techniques are reviewed.
Collapse
Affiliation(s)
- Sumudu P Leelananda
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH 43210, USA
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH 43210, USA
| |
Collapse
|
18
|
Ghanakota P, Carlson HA. Driving Structure-Based Drug Discovery through Cosolvent Molecular Dynamics. J Med Chem 2016; 59:10383-10399. [PMID: 27486927 DOI: 10.1021/acs.jmedchem.6b00399] [Citation(s) in RCA: 65] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Identifying binding hotspots on protein surfaces is of prime interest in structure-based drug discovery, either to assess the tractability of pursuing a protein target or to drive improved potency of lead compounds. Computational approaches to detect such regions have traditionally relied on energy minimization of probe molecules onto static protein conformations in the absence of the natural aqueous environment. Advances in high performance computing now allow us to assess hotspots using molecular dynamics (MD) simulations. MD simulations integrate protein flexibility and the complicated role of water, thereby providing a more realistic assessment of the complex kinetics and thermodynamics at play. In this review, we describe the evolution of various cosolvent-based MD techniques and highlight a myriad of potential applications for such technologies in computational drug development.
Collapse
Affiliation(s)
- Phani Ghanakota
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church Street, Ann Arbor, Michigan 48109-1065, United States
| | - Heather A Carlson
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church Street, Ann Arbor, Michigan 48109-1065, United States
| |
Collapse
|
19
|
Ung PMU, Ghanakota P, Graham SE, Lexa KW, Carlson HA. Identifying binding hot spots on protein surfaces by mixed-solvent molecular dynamics: HIV-1 protease as a test case. Biopolymers 2016; 105:21-34. [PMID: 26385317 DOI: 10.1002/bip.22742] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2015] [Revised: 09/14/2015] [Accepted: 09/14/2015] [Indexed: 12/16/2022]
Abstract
Mixed-solvent molecular dynamics (MixMD) simulations use full protein flexibility and competition between water and small organic probes to achieve accurate hot-spot mapping on protein surfaces. In this study, we improved MixMD using human immunodeficiency virus type-1 protease (HIVp) as the test case. We used three probe-water solutions (acetonitrile-water, isopropanol-water, and pyrimidine-water), first at 50% w/w concentration and later at 5% v/v. Paradoxically, better mapping was achieved by using fewer probes; 5% simulations gave a superior signal-to-noise ratio and far fewer spurious hot spots than 50% MixMD. Furthermore, very intense and well-defined probe occupancies were observed in the catalytic site and potential allosteric sites that have been confirmed experimentally. The Eye site, an allosteric site underneath the flap of HIVp, has been confirmed by the presence of a 5-nitroindole fragment in a crystal structure. MixMD also mapped two additional hot spots: the Exo site (between the Gly16-Gly17 and Cys67-Gly68 loops) and the Face site (between Glu21-Ala22 and Val84-Ile85 loops). The Exo site was observed to overlap with crystallographic additives such as acetate and dimethyl sulfoxide that are present in different crystal forms of the protein. Analysis of crystal structures of HIVp in different symmetry groups has shown that some surface sites are common interfaces for crystal contacts, which means that they are surfaces that are relatively easy to desolvate and complement with organic molecules. MixMD should identify these sites; in fact, their occupancy values help establish a solid cut-off where "druggable" sites are required to have higher occupancies than the crystal-packing faces.
Collapse
Affiliation(s)
- Peter M U Ung
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, MI, 48109-1065
| | - Phani Ghanakota
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, MI, 48109-1065
| | - Sarah E Graham
- Department of Biophysics, College of LSA, University of Michigan, 930 N. University St., Ann Arbor, MI, 48109-1055
| | - Katrina W Lexa
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, MI, 48109-1065
| | - Heather A Carlson
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, MI, 48109-1065.,Department of Biophysics, College of LSA, University of Michigan, 930 N. University St., Ann Arbor, MI, 48109-1055
| |
Collapse
|
20
|
Guvench O. Computational functional group mapping for drug discovery. Drug Discov Today 2016; 21:1928-1931. [PMID: 27393487 DOI: 10.1016/j.drudis.2016.06.030] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Revised: 06/23/2016] [Accepted: 06/29/2016] [Indexed: 01/05/2023]
Abstract
Computational functional group mapping (cFGM) is emerging as a high-impact complement to existing widely used experimental and computational structure-based drug discovery methods. cFGM provides comprehensive atomic-resolution 3D maps of the affinity of functional groups that can constitute drug-like molecules for a given target, typically a protein. These 3D maps can be intuitively and interactively visualized by medicinal chemists to rapidly design synthetically accessible ligands. Given that the maps can inform selection of functional groups for affinity, specificity, and pharmacokinetic properties, they are of utility for both the optimization of existing drug candidates and creating novel ones. Here, I review recent advances in cFGM with emphasis on the unique information content in the approach that offers the potential of broadly facilitating structure-based ligand design.
Collapse
Affiliation(s)
- Olgun Guvench
- SilcsBio, LLC, 8 Market Street, Suite 300, Baltimore, MD 21202, USA.
| |
Collapse
|
21
|
Chiang YC, Wang Y. Virtual substitution scan via single-step free energy perturbation. Biopolymers 2016; 105:324-336. [DOI: 10.1002/bip.22820] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2015] [Revised: 02/01/2016] [Accepted: 02/03/2016] [Indexed: 11/10/2022]
Affiliation(s)
- Ying-Chih Chiang
- Department of Physics; the Chinese University of Hong Kong; Shatin N.T. Hong Kong
| | - Yi Wang
- Department of Physics; the Chinese University of Hong Kong; Shatin N.T. Hong Kong
| |
Collapse
|
22
|
Raman EP, MacKerell AD. Spatial analysis and quantification of the thermodynamic driving forces in protein-ligand binding: binding site variability. J Am Chem Soc 2015; 137:2608-21. [PMID: 25625202 DOI: 10.1021/ja512054f] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The thermodynamic driving forces behind small molecule-protein binding are still not well-understood, including the variability of those forces associated with different types of ligands in different binding pockets. To better understand these phenomena we calculate spatially resolved thermodynamic contributions of the different molecular degrees of freedom for the binding of propane and methanol to multiple pockets on the proteins Factor Xa and p38 MAP kinase. Binding thermodynamics are computed using a statistical thermodynamics based end-point method applied on a canonical ensemble comprising the protein-ligand complexes and the corresponding free states in an explicit solvent environment. Energetic and entropic contributions of water and ligand degrees of freedom computed from the configurational ensemble provide an unprecedented level of detail into the mechanisms of binding. Direct protein-ligand interaction energies play a significant role in both nonpolar and polar binding, which is comparable to water reorganization energy. Loss of interactions with water upon binding strongly compensates these contributions leading to relatively small binding enthalpies. For both solutes, the entropy of water reorganization is found to favor binding in agreement with the classical view of the "hydrophobic effect". Depending on the specifics of the binding pocket, both energy-entropy compensation and reinforcement mechanisms are observed. It is notable to have the ability to visualize the spatial distribution of the thermodynamic contributions to binding at atomic resolution showing significant differences in the thermodynamic contributions of water to the binding of propane versus methanol.
Collapse
Affiliation(s)
- E Prabhu Raman
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , 20 Penn Street HSF II, Baltimore, Maryland 21201, United States
| | | |
Collapse
|
23
|
Yu W, Lakkaraju SK, Raman EP, Fang L, MacKerell AD. Pharmacophore modeling using site-identification by ligand competitive saturation (SILCS) with multiple probe molecules. J Chem Inf Model 2015; 55:407-20. [PMID: 25622696 DOI: 10.1021/ci500691p] [Citation(s) in RCA: 56] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Receptor-based pharmacophore modeling is an efficient computer-aided drug design technique that uses the structure of the target protein to identify novel leads. However, most methods consider protein flexibility and desolvation effects in a very approximate way, which may limit their use in practice. The Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling protocol (SILCS-Pharm) was introduced recently to address these issues, as SILCS naturally takes both protein flexibility and desolvation effects into account by using full molecular dynamics simulations to determine 3D maps of the functional group-affinity patterns on a target receptor. In the present work, the SILCS-Pharm protocol is extended to use a wider range of probe molecules including benzene, propane, methanol, formamide, acetaldehyde, methylammonium, acetate and water. This approach removes the previous ambiguity brought by using water as both the hydrogen-bond donor and acceptor probe molecule. The new SILCS-Pharm protocol is shown to yield improved screening results, as compared to the previous approach based on three target proteins. Further validation of the new protocol using five additional protein targets showed improved screening compared to those using common docking methods, further indicating improvements brought by the explicit inclusion of additional feature types associated with the wider collection of probe molecules in the SILCS simulations. The advantage of using complementary features and volume constraints, based on exclusion maps of the protein defined from the SILCS simulations, is presented. In addition, reranking using SILCS-based ligand grid free energies is shown to enhance the diversity of identified ligands for the majority of targets. These results suggest that the SILCS-Pharm protocol will be of utility in rational drug design.
Collapse
Affiliation(s)
- Wenbo Yu
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | | | | | | | | |
Collapse
|
24
|
Raman EP, Yu W, Lakkaraju SK, MacKerell AD. Inclusion of multiple fragment types in the site identification by ligand competitive saturation (SILCS) approach. J Chem Inf Model 2013; 53:3384-98. [PMID: 24245913 DOI: 10.1021/ci4005628] [Citation(s) in RCA: 87] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The site identification by ligand competitive saturation (SILCS) method identifies the location and approximate affinities of small molecular fragments on a target macromolecular surface by performing molecular dynamics (MD) simulations of the target in an aqueous solution of small molecules representative of different chemical functional groups. In this study, we introduce a set of small molecules to map potential interactions made by neutral hydrogen bond donors and acceptors and charged donor and acceptor fragments in addition to nonpolar fragments. The affinity pattern is obtained in the form of discretized probability or, equivalently, free energy maps, called FragMaps, which can be visualized with the target surface. We performed SILCS simulations for four proteins for which structural and thermodynamic data is available for multiple diverse ligands. Good overlap is shown between high affinity regions identified by the FragMaps and the crystallographic positions of ligand functional groups with similar chemical functionality, thus demonstrating the validity of the qualitative information obtained from the simulations. To test the ability of FragMaps in providing quantitative predictions, we calculate the previously introduced ligand grid free energy (LGFE) metric and observe its correspondence with experimentally measured binding affinity. LGFE is computed for different conformational ensembles and improvement in prediction is shown with increasing ligand conformational sampling. Ensemble generation includes a Monte Carlo sampling approach that uses the GFE FragMaps directly as the energy function. The results show that some but not all experimental trends are predicted and warrant improvements in the scoring methodology. In addition, the potential utility of atom-based free energy contributions to the LGFE scores and the use of multiple ligands in SILCS to identify displaceable water molecules during ligand design are discussed.
Collapse
Affiliation(s)
- E Prabhu Raman
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , 20 Penn Street HSF II, Baltimore, Maryland 21201 United States
| | | | | | | |
Collapse
|
25
|
Cao X, Yap JL, Newell-Rogers MK, Peddaboina C, Jiang W, Papaconstantinou HT, Jupitor D, Rai A, Jung KY, Tubin RP, Yu W, Vanommeslaeghe K, Wilder PT, MacKerell AD, Fletcher S, Smythe RW. The novel BH3 α-helix mimetic JY-1-106 induces apoptosis in a subset of cancer cells (lung cancer, colon cancer and mesothelioma) by disrupting Bcl-xL and Mcl-1 protein-protein interactions with Bak. Mol Cancer 2013; 12:42. [PMID: 23680104 PMCID: PMC3663763 DOI: 10.1186/1476-4598-12-42] [Citation(s) in RCA: 74] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Accepted: 05/02/2013] [Indexed: 11/10/2022] Open
Abstract
Background It has been shown in many solid tumors that the overexpression of the pro-survival Bcl-2 family members Bcl-2/Bcl-xL and Mcl-1 confers resistance to a variety of chemotherapeutic agents. We designed the BH3 α-helix mimetic JY-1-106 to engage the hydrophobic BH3-binding grooves on the surfaces of both Bcl-xL and Mcl-1. Methods JY-1-106–protein complexes were studied using molecular dynamics (MD) simulations and the SILCS methodology. We have evaluated the in vitro effects of JY-1-106 by using a fluorescence polarization (FP) assay, an XTT assay, apoptosis assays, and immunoprecipitation and western-blot assays. A preclinical human cancer xenograft model was used to test the efficacy of JY-1-106 in vivo. Results MD and SILCS simulations of the JY-1-106–protein complexes indicated the importance of the aliphatic side chains of JY-1-106 to binding and successfully predicted the improved affinity of the ligand for Bcl-xL over Mcl-1. Ligand binding affinities were measured via an FP assay using a fluorescently labeled Bak-BH3 peptide in vitro. Apoptosis induction via JY-1-106 was evidenced by TUNEL assay and PARP cleavage as well as by Bax–Bax dimerization. Release of multi-domain Bak from its inhibitory binding to Bcl-2/Bcl-xL and Mcl-1 using JY-1-106 was detected via immunoprecipitation (IP) western blotting. At the cellular level, we compared the growth proliferation IC50s of JY-1-106 and ABT-737 in multiple cancer cell lines with various Bcl-xL and Mcl-1 expression levels. JY-1-106 effectively induced cell death regardless of the Mcl-1 expression level in ABT-737 resistant solid tumor cells, whilst toxicity toward normal human endothelial cells was limited. Furthermore, synergistic effects were observed in A549 cells using a combination of JY-1-106 and multiple chemotherapeutic agents. We also observed that JY-1-106 was a very effective agent in inducing apoptosis in metabolically stressed tumors. Finally, JY-1-106 was evaluated in a tumor-bearing nude mouse model, and was found to effectively repress tumor growth. Strong TUNEL signals in the tumor cells demonstrated the effectiveness of JY-1-106 in this animal model. No significant side effects were observed in mouse organs after multiple injections. Conclusions Taken together, these observations demonstrate that JY-1-106 is an effective pan-Bcl-2 inhibitor with very promising clinical potential.
Collapse
Affiliation(s)
- Xiaobo Cao
- Department of Surgery, Scott & White Memorial Hospital and Clinic, The Texas A&M University System, Health Science Center, College of Medicine, 702 SW HK Dodgen Loop, Temple, Texas 76504, USA
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|