1
|
Maghsoud Y, Vázquez-Montelongo EA, Yang X, Liu C, Jing Z, Lee J, Harger M, Smith AK, Espinoza M, Guo HF, Kurie JM, Dalby KN, Ren P, Cisneros GA. Computational Investigation of a Series of Small Molecules as Potential Compounds for Lysyl Hydroxylase-2 (LH2) Inhibition. J Chem Inf Model 2023; 63:986-1001. [PMID: 36779232 PMCID: PMC10233724 DOI: 10.1021/acs.jcim.2c01448] [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] [Indexed: 02/01/2023]
Abstract
The catalytic function of lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required cosubstrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of 44 small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intramolecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than those in the apoenzyme, suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all of the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.
Collapse
Affiliation(s)
- Yazdan Maghsoud
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
| | - Erik Antonio Vázquez-Montelongo
- Department of Physical Medicine and Rehabilitation, University of Texas Southwestern Medical Center, Dallas, Texas 75390, United States
| | - Xudong Yang
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Juhoon Lee
- Division of Chemical Biology and Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Matthew Harger
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Ally K Smith
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Miguel Espinoza
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Hou-Fu Guo
- Department of Molecular and Cellular Biochemistry, College of Medicine, The University of Kentucky, Lexington, Kentucky 40536, United States
| | - Jonathan M Kurie
- Department of Thoracic/Head and Neck Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas 77005, United States
| | - Kevin N Dalby
- Division of Chemical Biology and Medicinal Chemistry, College of Pharmacy, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - G Andrés Cisneros
- Department of Chemistry and Biochemistry, The University of Texas at Dallas, Richardson, Texas 75080, United States
- Department of Physics, The University of Texas at Dallas, Richardson, Texas 75080, United States
| |
Collapse
|
2
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
3
|
Lambros E, Lipparini F, Cisneros GA, Paesani F. A Many-Body, Fully Polarizable Approach to QM/MM Simulations. J Chem Theory Comput 2020; 16:7462-7472. [PMID: 33213149 PMCID: PMC8131112 DOI: 10.1021/acs.jctc.0c00932] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
We present a new development in quantum mechanics/molecular mechanics (QM/MM) methods by replacing conventional MM models with data-driven many-body (MB) representations rigorously derived from high-level QM calculations. The new QM/MM approach builds on top of mutually polarizable QM/MM schemes developed for polarizable force fields with inducible dipoles and uses permutationally invariant polynomials to effectively account for quantum-mechanical contributions (e.g., exchange-repulsion and charge transfer and penetration) that are difficult to describe by classical expressions adopted by conventional MM models. Using the many-body MB-pol and MB-DFT potential energy functions for water, which include explicit two-body and three-body terms fitted to reproduce the corresponding CCSD(T) and PBE0 two-body and three-body energies for water, we demonstrate a smooth energetic transition as molecules are transferred between QM and MM regions, without the need of a transition layer. By effectively elevating the accuracy of both the MM region and the QM/MM interface to that of the QM region, the new QM/MB-MM approach achieves an accuracy comparable to that obtained with a fully QM treatment of the entire system.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | | | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
4
|
Vázquez-Montelongo EA, Vázquez-Cervantes JE, Cisneros GA. Polarizable ab initio QM/MM Study of the Reaction Mechanism of N- tert-Butyloxycarbonylation of Aniline in [EMIm][BF₄]. Molecules 2018; 23:E2830. [PMID: 30384470 PMCID: PMC6278528 DOI: 10.3390/molecules23112830] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2018] [Revised: 10/27/2018] [Accepted: 10/29/2018] [Indexed: 12/25/2022] Open
Abstract
N-t e r t-butoxycarbonylation of amines in solution (water, organic solvents, or ionic liquids) is a common reaction for the preparation of drug molecules. To understand the reaction mechanism and the role of the solvent, quantum mechanical/molecular mechanical simulations using a polarizable multipolar force field with long⁻range electrostatic corrections were used to optimize the minimum energy paths (MEPs) associated with various possible reaction mechanisms employing the nudged elastic band (NEB) and the quadratic string method (QSM). The calculated reaction energies and energy barriers were compared with the corresponding gas-phase and dichloromethane results. Complementary Electron Localization Function (ELF)/NCI analyses provide insights on the critical structures along the MEP. The calculated results suggest the most likely path involves a sequential mechanism with the rate⁻limiting step corresponding to the nucleophilic attack of the aniline, followed by proton transfer and the release of CO 2 without the direct involvement of imidazolium cations as catalysts.
Collapse
Affiliation(s)
| | | | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, TX 76201, USA.
- The Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, Denton, TX 76201, USA.
| |
Collapse
|
5
|
Ganguly A, Boulanger E, Thiel W. Importance of MM Polarization in QM/MM Studies of Enzymatic Reactions: Assessment of the QM/MM Drude Oscillator Model. J Chem Theory Comput 2017; 13:2954-2961. [DOI: 10.1021/acs.jctc.7b00016] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Abir Ganguly
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Eliot Boulanger
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
6
|
Boulanger E, Thiel W. Toward QM/MM Simulation of Enzymatic Reactions with the Drude Oscillator Polarizable Force Field. J Chem Theory Comput 2014; 10:1795-809. [DOI: 10.1021/ct401095k] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Eliot Boulanger
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| |
Collapse
|
7
|
Nakano H, Yamamoto T. Variational calculation of quantum mechanical/molecular mechanical free energy with electronic polarization of solvent. J Chem Phys 2012; 136:134107. [DOI: 10.1063/1.3699234] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
|
8
|
Wang Q, Bryce RA. Accounting for non-optimal interactions in molecular recognition: a study of ion-π complexes using a QM/MM model with a dipole-polarisable MM region. Phys Chem Chem Phys 2011; 13:19401-8. [PMID: 21960295 DOI: 10.1039/c1cp21944h] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
For a quantitative understanding of molecular structure, interaction and dynamics, accurate modelling of the energetics of both near-equilibrium and less optimal contacts is important. In this work, we explore the potential energy surfaces of representative ion-π complexes. We examine the performance of a semi-empirical QM/MM approach and the corresponding QM/MMpol model, where inducible point dipoles are additionally employed in the MM region. The predicted potential energy surfaces of cation-benzene complexes are improved by inclusion of explicit MM polarisation of the π-molecule. For cation-formamide complexes, inducible dipoles appreciably improve energetic estimates at geometries forming non-optimal interactions. Energetic component analysis suggests that the implicit MM polarisation of the fixed charge QM/MM model mirrors the behaviour of the QM/MMpol dipole model for the energetics of near-equilibrium conformations. However, for complexes at less optimal orientations, the QM/MM model exhibits higher errors than the QM/MMpol approach, being unable to capture orientation-dependent variations in polarisation energy.
Collapse
Affiliation(s)
- Qiantao Wang
- School of Pharmacy and Pharmaceutical Sciences, University of Manchester, Oxford Road, Manchester, M13 9PT, UK
| | | |
Collapse
|
9
|
Mazanetz MP, Ichihara O, Law RJ, Whittaker M. Prediction of cyclin-dependent kinase 2 inhibitor potency using the fragment molecular orbital method. J Cheminform 2011; 3:2. [PMID: 21219630 PMCID: PMC3032746 DOI: 10.1186/1758-2946-3-2] [Citation(s) in RCA: 85] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2010] [Accepted: 01/10/2011] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The reliable and robust estimation of ligand binding affinity continues to be a challenge in drug design. Many current methods rely on molecular mechanics (MM) calculations which do not fully explain complex molecular interactions. Full quantum mechanical (QM) computation of the electronic state of protein-ligand complexes has recently become possible by the latest advances in the development of linear-scaling QM methods such as the ab initio fragment molecular orbital (FMO) method. This approximate molecular orbital method is sufficiently fast that it can be incorporated into the development cycle during structure-based drug design for the reliable estimation of ligand binding affinity. Additionally, the FMO method can be combined with approximations for entropy and solvation to make it applicable for binding affinity prediction for a broad range of target and chemotypes. RESULTS We applied this method to examine the binding affinity for a series of published cyclin-dependent kinase 2 (CDK2) inhibitors. We calculated the binding affinity for 28 CDK2 inhibitors using the ab initio FMO method based on a number of X-ray crystal structures. The sum of the pair interaction energies (PIE) was calculated and used to explain the gas-phase enthalpic contribution to binding. The correlation of the ligand potencies to the protein-ligand interaction energies gained from FMO was examined and was seen to give a good correlation which outperformed three MM force field based scoring functions used to appoximate the free energy of binding. Although the FMO calculation allows for the enthalpic component of binding interactions to be understood at the quantum level, as it is an in vacuo single point calculation, the entropic component and solvation terms are neglected. For this reason a more accurate and predictive estimate for binding free energy was desired. Therefore, additional terms used to describe the protein-ligand interactions were then calculated to improve the correlation of the FMO derived values to experimental free energies of binding. These terms were used to account for the polar and non-polar solvation of the molecule estimated by the Poisson-Boltzmann equation and the solvent accessible surface area (SASA), respectively, as well as a correction term for ligand entropy. A quantitative structure-activity relationship (QSAR) model obtained by Partial Least Squares projection to latent structures (PLS) analysis of the ligand potencies and the calculated terms showed a strong correlation (r2 = 0.939, q2 = 0.896) for the 14 molecule test set which had a Pearson rank order correlation of 0.97. A training set of a further 14 molecules was well predicted (r2 = 0.842), and could be used to obtain meaningful estimations of the binding free energy. CONCLUSIONS Our results show that binding energies calculated with the FMO method correlate well with published data. Analysis of the terms used to derive the FMO energies adds greater understanding to the binding interactions than can be gained by MM methods. Combining this information with additional terms and creating a scaled model to describe the data results in more accurate predictions of ligand potencies than the absolute values obtained by FMO alone.
Collapse
Affiliation(s)
- Michael P Mazanetz
- Evotec (UK) limited, 114 Milton Park, Abingdon, Oxfordshire, OX14 4SA, UK
| | - Osamu Ichihara
- Evotec (UK) limited, 114 Milton Park, Abingdon, Oxfordshire, OX14 4SA, UK
| | - Richard J Law
- Evotec (UK) limited, 114 Milton Park, Abingdon, Oxfordshire, OX14 4SA, UK
| | - Mark Whittaker
- Evotec (UK) limited, 114 Milton Park, Abingdon, Oxfordshire, OX14 4SA, UK
| |
Collapse
|
10
|
Illingworth CJR, Furini S, Domene C. Computational Studies on Polarization Effects and Selectivity in K+ Channels. J Chem Theory Comput 2010. [DOI: 10.1021/ct100276c] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Christopher J. R. Illingworth
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Simone Furini
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Carmen Domene
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
11
|
Hayik SA, Dunbrack R, Merz KM. A Mixed QM/MM Scoring Function to Predict Protein-Ligand Binding Affinity. J Chem Theory Comput 2010; 6:3079-3091. [PMID: 21221417 DOI: 10.1021/ct100315g] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Computational methods for predicting protein-ligand binding free energy continue to be popular as a potential cost-cutting method in the drug discovery process. However, accurate predictions are often difficult to make as estimates must be made for certain electronic and entropic terms in conventional force field based scoring functions. Mixed quantum mechanics/molecular mechanics (QM/MM) methods allow electronic effects for a small region of the protein to be calculated, treating the remaining atoms as a fixed charge background for the active site. Such a semi-empirical QM/MM scoring function has been implemented in AMBER using DivCon and tested on a set of 23 metalloprotein-ligand complexes, where QM/MM methods provide a particular advantage in the modeling of the metal ion. The binding affinity of this set of proteins can be calculated with an R(2) of 0.64 and a standard deviation of 1.88 kcal/mol without fitting and 0.71 and a standard deviation of 1.69 kcal/mol with fitted weighting of the individual scoring terms. In this study we explore using various methods to calculate terms in the binding free energy equation, including entropy estimates and minimization standards. From these studies we found that using the rotational bond estimate to ligand entropy results in a reasonable R(2) of 0.63 without fitting. We also found that using the ESCF energy of the proteins without minimization resulted in an R(2) of 0.57, when using the rotatable bond entropy estimate.
Collapse
Affiliation(s)
- Seth A Hayik
- Institute for Cancer Research, Fox Chase Cancer Center, 333 Cottman Avenue, Philadelphia, PA 19111, USA
| | | | | |
Collapse
|
12
|
Ranaghan KE, Mulholland AJ. Investigations of enzyme-catalysed reactions with combined quantum mechanics/molecular mechanics (QM/MM) methods. INT REV PHYS CHEM 2010. [DOI: 10.1080/01442350903495417] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
13
|
Illingworth CJ, Domene C. Many-body effects and simulations of potassium channels. Proc Math Phys Eng Sci 2009. [DOI: 10.1098/rspa.2009.0014] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The electronic polarizability of an ion or a molecule is a measure of the relative tendency of its electron cloud to be distorted from its normal shape by an electric field. On the molecular scale, in a condensed phase, any species sits in an electric field due to its neighbours, and the resulting polarization is an important contribution to the total interaction energy. Electrostatic interactions are crucial for determining the majority of chemical–physical properties of the system and electronic polarization is a fundamental component of these interactions. Thus, polarization effects should be taken into account if accurate descriptions are desired. In classical computer simulations, the forces required to drive the system are typically based on interatomic interaction potentials derived in part from electronic structure calculations or from experimental data. Owing to the difficulties in including polarization effects in classical force fields, most of them are based just on pairwise additive interaction potentials. At present, major efforts are underway to develop polarizable interaction potentials for biomolecular simulations. In this review, various ways of introducing explicit polarizability into biomolecular models and force fields are reviewed, and the progress that might be achieved in applying such methods to study potassium channels is described.
Collapse
Affiliation(s)
- Christopher J. Illingworth
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of OxfordOxford OX1 3QZ, UK
| | - Carmen Domene
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of OxfordOxford OX1 3QZ, UK
| |
Collapse
|