1
|
Chehrazi E. Molecular Dynamics Simulations of Gas Transport Properties in Cross-Linked Polyamide Membranes: Tracing the Morphology and Addition of Silicate Nanotubes. ACS OMEGA 2024; 9:33425-33436. [PMID: 39130576 PMCID: PMC11307296 DOI: 10.1021/acsomega.3c10108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Revised: 05/12/2024] [Accepted: 07/17/2024] [Indexed: 08/13/2024]
Abstract
This study employs molecular dynamics (MD) simulations to fundamentally provide insight into the role of cross-link density in the CO2 separation properties of interfacially polymerized polyamide (PA) membranes. For this purpose, two atomistic models of pure polyamide membranes with different cross-link densities are constructed by MD simulations to conceptually determine how the fractional free volume of polyamide affects the gas separation performance of the membrane. The PA membrane with a lower cross-link density (LCPA) shows a higher gas diffusion coefficient, a lower gas solubility coefficient, and a higher gas permeability than the PA membrane with a higher cross-link density (HCPA). Moreover, the pristine and modified silicate nanotubes (SNTs) as the fast gas transport channels are incorporated into the polyamide membranes to assess the effect of the SNT/PA interface chemistry on the CO2 separation properties of the membranes. SNTs are systematically modified by three modifying agents with different CO2-philic groups and different interfacial interaction energies with the polyamide matrix. The results of MD simulations demonstrate that the incorporation of silicate nanotubes into the PA matrix increases the gas diffusivity and permeability and decreases the CO2/gas selectivity. Moreover, the membranes containing modified SNTs possessing high CO2-philicity and high SNTs/PA interfacial interactions show a high CO2 separation performance.
Collapse
Affiliation(s)
- Ehsan Chehrazi
- Department of Polymer Chemistry
and Materials, Faculty of Chemistry and Petroleum Sciences, Shahid Beheshti University, Tehran 1983969411, Iran
| |
Collapse
|
2
|
Tkaczyk S, Karwounopoulos J, Schöller A, Woodcock HL, Langer T, Boresch S, Wieder M. Reweighting from Molecular Mechanics Force Fields to the ANI-2x Neural Network Potential. J Chem Theory Comput 2024; 20:2719-2728. [PMID: 38527958 DOI: 10.1021/acs.jctc.3c01274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/27/2024]
Abstract
To achieve chemical accuracy in free energy calculations, it is necessary to accurately describe the system's potential energy surface and efficiently sample configurations from its Boltzmann distribution. While neural network potentials (NNPs) have shown significantly higher accuracy than classical molecular mechanics (MM) force fields, they have a limited range of applicability and are considerably slower than MM potentials, often by orders of magnitude. To address this challenge, Rufa et al. [Rufa et al. bioRxiv 2020, 10.1101/2020.07.29.227959.] suggested a two-stage approach that uses a fast and established MM alchemical energy protocol, followed by reweighting the results using NNPs, known as endstate correction or indirect free energy calculation. This study systematically investigates the accuracy and robustness of reweighting from an MM reference to a neural network target potential (ANI-2x) for an established data set in vacuum, using single-step free-energy perturbation (FEP) and nonequilibrium (NEQ) switching simulation. We assess the influence of longer switching lengths and the impact of slow degrees of freedom on outliers in the work distribution and compare the results to those of multistate equilibrium free energy simulations. Our results demonstrate that free energy calculations between NNPs and MM potentials should be preferably performed using NEQ switching simulations to obtain accurate free energy estimates. NEQ switching simulations between the MM potentials and NNPs are efficient, robust, and trivial to implement.
Collapse
Affiliation(s)
- Sara Tkaczyk
- Department of Pharmaceutical Sciences, Pharmaceutical Chemistry Division, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
- Vienna Doctoral School of Pharmaceutical, Nutritional and Sport Sciences (PhaNuSpo), University of Vienna, 1090 Vienna, Austria
| | - Johannes Karwounopoulos
- Faculty of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Währingerstrasse 17, 1090 Vienna, Austria
- Vienna Doctoral School of Chemistry (DoSChem), University of Vienna, Währingerstrasse 42, 1090 Vienna, Austria
| | - Andreas Schöller
- Faculty of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Währingerstrasse 17, 1090 Vienna, Austria
- Vienna Doctoral School of Chemistry (DoSChem), University of Vienna, Währingerstrasse 42, 1090 Vienna, Austria
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, 4202 E. Fowler Ave., CHE205, Tampa, Florida 33620-5250, United States
| | - Thierry Langer
- Department of Pharmaceutical Sciences, Pharmaceutical Chemistry Division, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Stefan Boresch
- Faculty of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Währingerstrasse 17, 1090 Vienna, Austria
| | - Marcus Wieder
- Faculty of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Währingerstrasse 17, 1090 Vienna, Austria
| |
Collapse
|
3
|
Rizzi A, Carloni P, Parrinello M. Free energies at QM accuracy from force fields via multimap targeted estimation. Proc Natl Acad Sci U S A 2023; 120:e2304308120. [PMID: 37931103 PMCID: PMC10655219 DOI: 10.1073/pnas.2304308120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Accepted: 09/25/2023] [Indexed: 11/08/2023] Open
Abstract
Accurate predictions of ligand binding affinities would greatly accelerate the first stages of drug discovery campaigns. However, using highly accurate interatomic potentials based on quantum mechanics (QM) in free energy methods has been so far largely unfeasible due to their prohibitive computational cost. Here, we present an efficient method to compute QM free energies from simulations using cheap reference potentials, such as force fields (FFs). This task has traditionally been out of reach due to the slow convergence of computing the correction from the FF to the QM potential. To overcome this bottleneck, we generalize targeted free energy methods to employ multiple maps-implemented with normalizing flow neural networks (NNs)-that maximize the overlap between the distributions. Critically, the method requires neither a separate expensive training phase for the NNs nor samples from the QM potential. We further propose a one-epoch learning policy to efficiently avoid overfitting, and we combine our approach with enhanced sampling strategies to overcome the pervasive problem of poor convergence due to slow degrees of freedom. On the drug-like molecules in the HiPen dataset, the method accelerates the calculation of the free energy difference of switching from an FF to a DFTB3 potential by three orders of magnitude compared to standard free energy perturbation and by a factor of eight compared to previously published nonequilibrium calculations. Our results suggest that our method, in combination with efficient QM/MM calculations, may be used in lead optimization campaigns in drug discovery and to study protein-ligand molecular recognition processes.
Collapse
Affiliation(s)
- Andrea Rizzi
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich52428, Germany
- Atomistic Simulations, Italian Institute of Technology, Genova16163, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich52428, Germany
- Department of Physics and Universitätsklinikum, RWTH Aachen University, Aachen52074, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Genova16163, Italy
| |
Collapse
|
4
|
Khadka D, Jayasinghe-Arachchige VM, Prabhakar R, Ramamurthy V. Application of molecular dynamic simulations in modeling the excited state behavior of confined molecules. Photochem Photobiol Sci 2023:10.1007/s43630-023-00486-2. [PMID: 37843722 DOI: 10.1007/s43630-023-00486-2] [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: 06/27/2023] [Accepted: 09/20/2023] [Indexed: 10/17/2023]
Abstract
Relative to isotropic organic solvent medium, the structure and conformation of a reactant molecule in an organized and confining medium are often different. In addition, because of the rigidity of the immediate environment, the reacting molecule have a little freedom to undergo large changes even upon gaining energy or modifications in the electronic structure. These alterations give rise to differences in the photochemistry of a molecular and supramolecular species. In this study, one such example is presented. α-Alkyl dibenzylketones upon excitation in isotropic solvents give products via Norrish type I and type II reactions that are independent of the chain length of the alkyl substituent. On the other hand, when these molecules are enclosed within an organic capsule of volume ~ 550 Å3, they give products that are strikingly dependent on the length of the α-alkyl substitution. These previously reported experimental observations are rationalized based on the structures generated by molecular modeling (docking and molecular dynamics (MD) simulations). It is shown that MD simulations that are utilized extensively in biologically important macromolecules can also be useful to understand the excited state behavior of reactive molecules that are part of supramolecular assemblies. These simulations can provide structural information of the reactant molecule and the surroundings complementing that with the one obtained from 1 and 2D NMR experiments. MD simulated structures of seven α-alkyl dibenzylketones encapsulated within the octa acid capsule provide a clear understanding of their unique behavior in this restricted medium. Because of the rigidity of the medium, these structures although generated in the ground state can rationalize the photochemical behavior of the molecules in the excited state.
Collapse
Affiliation(s)
- Dipendra Khadka
- Department of Chemistry, University of Miami, Coral Gables, FL, 33124, USA
| | | | - Rajeev Prabhakar
- Department of Chemistry, University of Miami, Coral Gables, FL, 33124, USA.
| | | |
Collapse
|
5
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
6
|
Barros EP, Ries B, Böselt L, Champion C, Riniker S. Recent developments in multiscale free energy simulations. Curr Opin Struct Biol 2021; 72:55-62. [PMID: 34534706 DOI: 10.1016/j.sbi.2021.08.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 08/06/2021] [Accepted: 08/16/2021] [Indexed: 11/26/2022]
Abstract
Physics-based free energy simulations enable the rigorous calculation of properties, such as conformational equilibria, solvation or binding free energies. While historically most applications have occurred at the atomistic level of resolution, a range of advances in the past years make it possible now to reliably cross the temporal, spatial and theory scales for the modeling of complex systems or the efficient prediction of results at the accuracy level of expensive quantum-mechanical calculations. In this mini-review, we discuss recent methodological advances as well as opportunities opened up by the introduction of machine learning approaches, which tackle the diverse challenges across the different scales, improve the accuracy and feasibility, and push the boundaries of multiscale free energy simulations.
Collapse
Affiliation(s)
- Emilia P Barros
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Benjamin Ries
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Lennard Böselt
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Candide Champion
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093, Zurich, Switzerland.
| |
Collapse
|
7
|
Galimberti DR, Sauer J. Chemically Accurate Vibrational Free Energies of Adsorption from Density Functional Theory Molecular Dynamics: Alkanes in Zeolites. J Chem Theory Comput 2021; 17:5849-5862. [PMID: 34459582 PMCID: PMC8444336 DOI: 10.1021/acs.jctc.1c00519] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a methodology
to compute, at reduced computational cost,
Gibbs free energies, enthalpies, and entropies of adsorption from
molecular dynamics. We calculate vibrational partition functions from
vibrational energies, which we obtain from the vibrational density
of states by projection on the normal modes. The use of a set of well-chosen
reference structures along the trajectories accounts for the anharmonicities
of the modes. For the adsorption of methane, ethane, and propane in
the H-CHA zeolite, we limit our treatment to a set of vibrational
modes localized at the adsorption site (zeolitic OH group) and the
alkane molecule interacting with it. Only two short trajectories (1–20
ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic
functions. The mean absolute deviations from the experimentally measured
values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the
enthalpy, and the entropy term (−TΔS),
respectively. In particular, the entropy terms show a major improvement
compared to the harmonic approximation and almost reach the accuracy
of the previous use of anharmonic frequencies obtained with curvilinear
distortions of individual modes. The thermodynamic functions so obtained
follow the trend of the experimental values for methane, ethane, and
propane, and the Gibbs free energy of adsorption at experimental conditions
is correctly predicted to change from positive for methane (5.9 kJ/mol)
to negative for ethane (−4.8 kJ/mol) and propane (−7.1
kJ/mol).
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany.,Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität, Unter den Linden 6, 10117 Berlin, Germany
| |
Collapse
|
8
|
Chen LD, Lawniczak JJ, Ding F, Bygrave PJ, Riahi S, Manby FR, Mukhopadhyay S, Miller TF. Embedded Mean-Field Theory for Solution-Phase Transition-Metal Polyolefin Catalysis. J Chem Theory Comput 2020; 16:4226-4237. [PMID: 32441933 DOI: 10.1021/acs.jctc.0c00169] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Decreasing the wall-clock time of quantum mechanics/molecular mechanics (QM/MM) calculations without sacrificing accuracy is a crucial prerequisite for widespread simulation of solution-phase dynamical processes. In this work, we demonstrate the use of embedded mean-field theory (EMFT) as the QM engine in QM/MM molecular dynamics (MD) simulations to examine polyolefin catalysts in solution. We show that employing EMFT in this mode preserves the accuracy of hybrid-functional DFT in the QM region, while providing up to 20-fold reductions in the cost per SCF cycle, thereby increasing the accessible simulation time-scales. We find that EMFT reproduces DFT-computed binding energies and optimized bond lengths to within chemical accuracy, as well as consistently ranking conformer stability. Furthermore, solution-phase EMFT/MM simulations provide insight into the interaction strength of strongly coordinating and bulky counterions.
Collapse
Affiliation(s)
- Leanne D Chen
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - James J Lawniczak
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Feizhi Ding
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Peter J Bygrave
- Centre for Computational Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Saleh Riahi
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Frederick R Manby
- Centre for Computational Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | | | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
9
|
Pecina A, Eyrilmez SM, Köprülüoğlu C, Miriyala VM, Lepšík M, Fanfrlík J, Řezáč J, Hobza P. SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design. Chempluschem 2020; 85:2362-2371. [PMID: 32609421 DOI: 10.1002/cplu.202000120] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 05/27/2020] [Indexed: 12/17/2022]
Abstract
Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Saltuk M Eyrilmez
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Vijay Madhav Miriyala
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| |
Collapse
|
10
|
Dong TG, Peng H, He XF, Wang X, Gao J. Hybrid Molecular Dynamics for Elucidating Cooperativity Between Halogen Bond and Water Molecules During the Interaction of p53-Y220C and the PhiKan5196 Complex. Front Chem 2020; 8:344. [PMID: 32457871 PMCID: PMC7221198 DOI: 10.3389/fchem.2020.00344] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 04/02/2020] [Indexed: 12/15/2022] Open
Abstract
The cooperativity between hydrogen and halogen bonds plays an important role in rational drug design. However, mimicking the dynamic cooperation between these bonds is a challenging issue, which has impeded the development of the halogen bond force field. In this study, the Y220C–PhiKan5196 complex of p53 protein was adopted as a model, and the functions of three water molecules that formed hydrogen bonds with halogen atoms were analyzed by the simulation method governed by the hybrid quantum mechanical/molecular mechanical molecular dynamics. A comparison with the water-free model revealed that the strength of the halogen bond in the complex was consistently stronger. This confirmed that the water molecules formed weak hydrogen bonds with the halogen atom and cooperated with the halogen atom to enhance the halogen bond. Further, it was discovered that the roles of the three water molecules were not the same. Therefore, the results obtained herein can facilitate a rational drug design. Further, this work emphasizes on the fact that, in addition to protein pockets and ligands, the role of voids should also be considered with regard to the water molecules surrounding them.
Collapse
Affiliation(s)
- Tian-Ge Dong
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Hui Peng
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Xue-Feng He
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Xiaocong Wang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| |
Collapse
|
11
|
Anumalla B, Prabhu NP. Surface hydration and preferential interaction directs the charged amino acids-induced changes in protein stability. J Mol Graph Model 2020; 98:107602. [PMID: 32251994 DOI: 10.1016/j.jmgm.2020.107602] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Revised: 03/20/2020] [Accepted: 03/23/2020] [Indexed: 01/23/2023]
Abstract
In the present study, we investigate the interaction of amino acid osmolytes, Arg, Lys, Asp and Glu, and a denaturant, guanidinium chloride (Gdm) with proteins. To achieve this, molecular dynamics (MD) simulation of RNase A and α-lactalbumin was performed in the presence of three charged amino acids Arg, Lys, and Asp and the molecular mechanism of amino acid-induced (de)stabilization of the proteins was examined by combining with our earlier report on Glu. As Arg has the side chain similar to that of Gdm and destabilizes the proteins, MD simulation was carried out in the presence of Gdm as well. Radial distribution function and hydration fraction around the protein surface reveals that preferential hydration increases upon the addition of any of the cosolvent; however, the extent of increase is more in the presence of stabilizing cosolvents (stAAs: Lys, Asp and Glu) compared to destabilizing cosolvents (Arg and Gdm). Moreover, the preferential interaction of Arg and Gdm with the proteins is higher than that of stAAs. Residue-level interaction analysis suggests that stAAs preferably interacts with charged amino acids of the proteins whereas Arg and Gdm interactions could be found on almost all the surface exposed residues which might provide higher preferential interaction for these residues. From the results, we propose that the net outcome of preferential hydration versus preferential interaction of the amino acids might determine their effect on the stability of proteins.
Collapse
Affiliation(s)
- Bramhini Anumalla
- Department of Biotechnology and Bioinformatics, School of Life Sciences, University of Hyderabad, Hyderabad, 500 046, India
| | - N Prakash Prabhu
- Department of Biotechnology and Bioinformatics, School of Life Sciences, University of Hyderabad, Hyderabad, 500 046, India.
| |
Collapse
|
12
|
Abstract
Quantum mechanics (QM) methods provide a fine description of receptor-ligand interactions and of chemical reactions. Their use in drug design and drug discovery is increasing, especially for complex systems including metal ions in the binding sites, for the design of highly selective inhibitors, for the optimization of bi-specific compounds, to understand enzymatic reactions, and for the study of covalent ligands and prodrugs. They are also used for generating molecular descriptors for predictive QSAR/QSPR models and for the parameterization of force fields. Thanks to the continuous increase of computational power offered by GPUs and to the development of sophisticated algorithms, QM methods are becoming part of the standard tools used in computer-aided drug design (CADD). We present the most used QM methods and software packages, and we discuss recent representative applications in drug design and drug discovery.
Collapse
Affiliation(s)
- Martin Kotev
- Global Research Informatics/Cheminformatics and Drug Design, Evotec (France) SAS, Toulouse, France
| | - Laurie Sarrat
- Global Research Informatics/Cheminformatics and Drug Design, Evotec (France) SAS, Toulouse, France
| | | |
Collapse
|
13
|
Giese TJ, York DM. Development of a Robust Indirect Approach for MM → QM Free Energy Calculations That Combines Force-Matched Reference Potential and Bennett's Acceptance Ratio Methods. J Chem Theory Comput 2019; 15:5543-5562. [PMID: 31507179 DOI: 10.1021/acs.jctc.9b00401] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We use the PBE0/6-31G* density functional method to perform ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations under periodic boundary conditions with rigorous electrostatics using the ambient potential composite Ewald method in order to test the convergence of MM → QM/MM free energy corrections for the prediction of 17 small-molecule solvation free energies and eight ligand binding free energies to T4 lysozyme. The "indirect" thermodynamic cycle for calculating free energies is used to explore whether a series of reference potentials improve the statistical quality of the predictions. Specifically, we construct a series of reference potentials that optimize a molecular mechanical (MM) force field's parameters to reproduce the ab initio QM/MM forces from a QM/MM simulation. The optimizations form a systematic progression of successively expanded parameters that include bond, angle, dihedral, and charge parameters. For each reference potential, we calculate benchmark quality reference values for the MM → QM/MM correction by performing the mixed MM and QM/MM Hamiltonians at 11 intermediate states, each for 200 ps. We then compare forward and reverse application of Zwanzig's relation, thermodynamic integration (TI), and Bennett's acceptance ratio (BAR) methods as a function of reference potential, simulation time, and the number of simulated intermediate states. We find that Zwanzig's equation is inadequate unless a large number of intermediate states are explicitly simulated. The TI and BAR mean signed errors are very small even when only the end-state simulations are considered, and the standard deviations of the TI and BAR errors are decreased by choosing a reference potential that optimizes the bond and angle parameters. We find a robust approach for the data sets of fairly rigid molecules considered here is to use bond + angle reference potential together with the end-state-only BAR analysis. This requires QM/MM simulations to be performed in order to generate reference data to parametrize the bond + angle reference potential, and then this same simulation serves a dual purpose as the full QM/MM end state. The convergence of the results with respect to time suggests that computational resources may be used more efficiently by running multiple simulations for no more than 50 ps, rather than running one long simulation.
Collapse
Affiliation(s)
- Timothy J Giese
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| | - 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-8087 , United States
| |
Collapse
|
14
|
Molecular mechanisms of the protein-protein interaction-regulated binding specificity of basic-region leucine zipper transcription factors. J Mol Model 2019; 25:246. [PMID: 31342181 DOI: 10.1007/s00894-019-4138-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Accepted: 07/14/2019] [Indexed: 10/26/2022]
Abstract
It is well known that the DNA-binding specificity of transcription factors (TFs) is influenced by protein-protein interactions (PPIs). However, the underlying molecular mechanisms remain largely unknown. In this work, we adopted the cAMP-response element-binding protein (CREB) of the basic leucine zipper (bZIP) TF family as a model system, and a workflow of combined bioinformatics and molecular modeling analysis of protein-DNA interaction was tested. First, the multiple sequence alignment and SDPsite method were used to find potential bZIP family binding specificity determining positions (SDPs) within the protein-protein interaction region. Second, the mutation system was analyzed using molecular dynamics simulation. Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) free energy calculations confirmed the enhancement of the binding affinity of the mutation, which was in agreement with experimental results. The root mean square fluctuation (RMSF) and hydrogen bonding changes suggested an open and close protein dimerization process after the system was mutated, which resulted in the change of the hydrogen bonding of the protein-DNA interface and a slight conformational change. We believe that this work will contribute to understanding the protein-protein interaction-regulated binding specificity of bZIP transcription factors.
Collapse
|
15
|
Hudson PS, Woodcock HL, Boresch S. Use of Interaction Energies in QM/MM Free Energy Simulations. J Chem Theory Comput 2019; 15:4632-4645. [PMID: 31142113 DOI: 10.1021/acs.jctc.9b00084] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The use of the most accurate (i.e., QM or QM/MM) levels of theory for free energy simulations (FES) is typically not possible. Primarily, this is because the computational cost associated with the extensive configurational sampling needed for converging FES is prohibitive. To ensure the feasibility of QM-based FES, the "indirect" approach is generally taken, necessitating a free energy calculation between the MM and QM/MM potential energy surfaces. Ideally, this step is performed with standard free energy perturbation (Zwanzig's equation) as it only requires simulations be carried out at the low level of theory; however, work from several groups over the past few years has conclusively shown that Zwanzig's equation is ill-suited to this task. As such, many approximations have arisen to mitigate difficulties with Zwanzig's equation. One particularly popular notion is that the convergence of Zwanzig's equation can be improved by using interaction energy differences instead of total energy differences. Although problematic numerical fluctuations (a major problem when using Zwanzig's equation) are indeed reduced, our results and analysis demonstrate that this "interaction energy approximation" (IEA) is theoretically incorrect, and the implicit approximation invoked is spurious at best. Herein, we demonstrate this via solvation free energy calculations using IEA from two different low levels of theory to the same target high level. Results from this proof-of-concept consistently yield the wrong results, deviating by ∼1.5 kcal/mol from the rigorously obtained value.
Collapse
Affiliation(s)
- Phillip S Hudson
- Department of Chemistry , University of South Florida , 4202 East Fowler Avenue, CHE205 , Tampa , Florida 33620-5250 , United States.,Laboratory of Computational Biology , National Institutes of Health, National Heart, Lung and Blood Institute , 12 South Drive, Rm 3053 , Bethesda , Maryland 20892-5690 , United States
| | - H Lee Woodcock
- Department of Chemistry , University of South Florida , 4202 East Fowler Avenue, CHE205 , Tampa , Florida 33620-5250 , United States
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry , University of Vienna , Währingerstraße 17 , Vienna A-1090 , Austria
| |
Collapse
|
16
|
Sokka IK, Ekholm FS, Johansson MP. Increasing the Potential of the Auristatin Cancer-Drug Family by Shifting the Conformational Equilibrium. Mol Pharm 2019; 16:3600-3608. [PMID: 31199662 PMCID: PMC6750905 DOI: 10.1021/acs.molpharmaceut.9b00437] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
![]()
Monomethyl auristatin E and monomethyl
auristatin F are widely
used cytotoxic agents in antibody–drug conjugates (ADCs), a
group of promising cancer drugs. The ADCs specifically target cancer
cells, releasing the auristatins inside, which results in the prevention
of mitosis. The auristatins suffer from a potentially serious flaw,
however. In solution, the molecules exist in an equal mixture of two
conformers, cis and trans. Only the trans-isomer is biologically active
and the isomerization process, i.e., the conversion of cis to trans
is slow. This significantly diminishes the efficiency of the drugs
and their corresponding ADCs, and perhaps more importantly, raises
concerns over drug safety. The potency of the auristatins would be
enhanced by decreasing the amount of the biologically inactive isomer,
either by stabilizing the trans-isomer or destabilizing the cis-isomer.
Here, we follow the computer-aided design strategy of shifting the
conformational equilibrium and employ high-level quantum chemical
modeling to identify promising candidates for improved auristatins.
Coupled cluster calculations predict that a simple halogenation in
the norephedrine/phenylalanine residues shifts the isomer equilibrium
almost completely toward the active trans-conformation, due to enhanced
intramolecular interactions specific to the active isomer.
Collapse
Affiliation(s)
- Iris K Sokka
- Department of Chemistry , University of Helsinki , P.O. Box 55, FI-00014 Helsinki , Finland
| | - Filip S Ekholm
- Department of Chemistry , University of Helsinki , P.O. Box 55, FI-00014 Helsinki , Finland
| | - Mikael P Johansson
- Department of Chemistry , University of Helsinki , P.O. Box 55, FI-00014 Helsinki , Finland.,Helsinki Institute of Sustainability Science, HELSUS , FI-00014 Helsinki , Finland
| |
Collapse
|
17
|
Li J, Fu A, Zhang L. An Overview of Scoring Functions Used for Protein-Ligand Interactions in Molecular Docking. Interdiscip Sci 2019; 11:320-328. [PMID: 30877639 DOI: 10.1007/s12539-019-00327-w] [Citation(s) in RCA: 185] [Impact Index Per Article: 37.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2018] [Revised: 02/06/2019] [Accepted: 03/06/2019] [Indexed: 12/17/2022]
Abstract
Currently, molecular docking is becoming a key tool in drug discovery and molecular modeling applications. The reliability of molecular docking depends on the accuracy of the adopted scoring function, which can guide and determine the ligand poses when thousands of possible poses of ligand are generated. The scoring function can be used to determine the binding mode and site of a ligand, predict binding affinity and identify the potential drug leads for a given protein target. Despite intensive research over the years, accurate and rapid prediction of protein-ligand interactions is still a challenge in molecular docking. For this reason, this study reviews four basic types of scoring functions, physics-based, empirical, knowledge-based, and machine learning-based scoring functions, based on an up-to-date classification scheme. We not only discuss the foundations of the four types scoring functions, suitable application areas and shortcomings, but also discuss challenges and potential future study directions.
Collapse
Affiliation(s)
- Jin Li
- College of Computer and Information Science, Southwest University, Chongqing, 400715, China.,School of Medical Information and Engineering, Southwest Medical University, Luzhou, 646000, China
| | - Ailing Fu
- College of Pharmaceutical Sciences, Southwest University, Chongqing, 400715, China
| | - Le Zhang
- College of Computer and Information Science, Southwest University, Chongqing, 400715, China. .,College of Computer Science, Sichuan University, Chengdu, 610065, China. .,Medical Big Data Center, Sichuan University, Chengdu, 610065, China. .,Zdmedical, Information Polytron Technologies Inc Chongqing, Chongqing, 401320, China.
| |
Collapse
|
18
|
Wang M, Mei Y, Ryde U. Host-Guest Relative Binding Affinities at Density-Functional Theory Level from Semiempirical Molecular Dynamics Simulations. J Chem Theory Comput 2019; 15:2659-2671. [PMID: 30811192 DOI: 10.1021/acs.jctc.8b01280] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Relative free energies for the binding of nine cyclic carboxylate ligands to the octa-acid deep-cavity host were calculated at the combined density-functional theory and molecular mechanics (DFT/MM) level of theory. The DFT calculations employed the BLYP functional and the 6-31G* basis set for the ligand. We employed free-energy perturbations (FEP) with the reference-potential approach and used molecular dynamics (MD) simulations with the semiempirical quantum mechanical (SQM) PM6-DH+ method for the ligand as an intermediate level between MM and DFT/MM to improve the convergence. Thus, the relative binding free energy of two ligands was first calculated at the MM level by an alchemical transformation from one ligand to another in both the bound and unbound states. Then, for each ligand the free-energy correction for going from the MM to the SQM/MM potentials was calculated using explicit SQM/MM MD simulations. Finally, the free-energy correction for going from the SQM/MM to the DFT/MM potentials was estimated with FEP without running any DFT/MM simulations. Instead, the free energy was calculated by single-step exponential averaging (ssEA) or employing the cumulant approximation to the second order (CA). The results show that CA converges much better than ssEA, and with 500-4500 DFT/MM single-point energy calculations, converged free energies with a precision of 0.3 kJ/mol can be obtained. These free energies reproduce the experimental binding free energy differences with a mean absolute deviation of 3.4 kJ/mol, a correlation ( R2) of 0.97, and correct signs for all of the eight free-energy differences. This is appreciably better than the results obtained at the SQM/MM level of theory and also slightly better than those obtained with MM. We show that the convergence of the SQM/MM → DFT/MM perturbations can be monitored by the use of Wu and Kofke's bias metric Π and by the standard deviation of the difference between the SQM/MM and DFT/MM energies. Finally, we show that the use of the intermediate SQM/MM MD simulations improves the convergence of the free energies by a factor of at least two, compared to doing direct MM → DFT/MM perturbations.
Collapse
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science , East China Normal University , Shanghai 200062 , China.,Department of Theoretical Chemistry , Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science , East China Normal University , Shanghai 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062 , China.,Collaborative Innovation Center of Extreme Optics , Shanxi University , Taiyuan , Shanxi 030006 , China
| | - Ulf Ryde
- Department of Theoretical Chemistry , Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
19
|
Wang M, Mei Y, Ryde U. Predicting Relative Binding Affinity Using Nonequilibrium QM/MM Simulations. J Chem Theory Comput 2018; 14:6613-6622. [DOI: 10.1021/acs.jctc.8b00685] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
20
|
Liu M, Youmans KN, Gao J. Dual QM and MM Approach for Computing Equilibrium Isotope Fractionation Factor of Organic Species in Solution. Molecules 2018; 23:E2644. [PMID: 30326599 PMCID: PMC6222756 DOI: 10.3390/molecules23102644] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Revised: 10/10/2018] [Accepted: 10/11/2018] [Indexed: 11/17/2022] Open
Abstract
A dual QM and MM approach for computing equilibrium isotope effects has been described. In the first partition, the potential energy surface is represented by a combined quantum mechanical and molecular mechanical (QM/MM) method, in which a solute molecule is treated quantum mechanically, and the remaining solvent molecules are approximated classically by molecular mechanics. In the second QM/MM partition, differential nuclear quantum effects responsible for the isotope effect are determined by a statistical mechanical double-averaging formalism, in which the nuclear centroid distribution is sampled classically by Newtonian molecular dynamics and the quantum mechanical spread of quantized particles about the centroid positions is treated using the path integral (PI) method. These partitions allow the potential energy surface to be properly represented such that the solute part is free of nuclear quantum effects for nuclear quantum mechanical simulations, and the double-averaging approach has the advantage of sampling efficiency for solvent configuration and for path integral convergence. Importantly, computational precision is achieved through free energy perturbation (FEP) theory to alchemically mutate one isotope into another. The PI-FEP approach is applied to model systems for the 18O enrichment found in cellulose of trees to determine the isotope enrichment factor of carbonyl compounds in water. The present method may be useful as a general tool for studying isotope fractionation in biological and geochemical systems.
Collapse
Affiliation(s)
- Meiyi Liu
- Laboratory of Theoretical and Computational Chemistry, Theoretical Chemistry Institute, Jilin University, Changchun 130023, China.
| | - Katelyn N Youmans
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455, USA.
| | - Jiali Gao
- Laboratory of Theoretical and Computational Chemistry, Theoretical Chemistry Institute, Jilin University, Changchun 130023, China.
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455, USA.
| |
Collapse
|
21
|
Huggins DJ, Biggin PC, Dämgen MA, Essex JW, Harris SA, Henchman RH, Khalid S, Kuzmanic A, Laughton CA, Michel J, Mulholland AJ, Rosta E, Sansom MSP, van der Kamp MW. Biomolecular simulations: From dynamics and mechanisms to computational assays of biological activity. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1393] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- David J. Huggins
- TCM Group, Cavendish Laboratory University of Cambridge Cambridge UK
- Unilever Centre, Department of Chemistry University of Cambridge Cambridge UK
- Department of Physiology and Biophysics Weill Cornell Medical College New York NY
| | | | - Marc A. Dämgen
- Department of Biochemistry University of Oxford Oxford UK
| | - Jonathan W. Essex
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | - Sarah A. Harris
- School of Physics and Astronomy University of Leeds Leeds UK
- Astbury Centre for Structural and Molecular Biology University of Leeds Leeds UK
| | - Richard H. Henchman
- Manchester Institute of Biotechnology The University of Manchester Manchester UK
- School of Chemistry The University of Manchester Oxford UK
| | - Syma Khalid
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | | | - Charles A. Laughton
- School of Pharmacy University of Nottingham Nottingham UK
- Centre for Biomolecular Sciences University of Nottingham Nottingham UK
| | - Julien Michel
- EaStCHEM school of Chemistry University of Edinburgh Edinburgh UK
| | - Adrian J. Mulholland
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
| | - Edina Rosta
- Department of Chemistry King's College London London UK
| | | | - Marc W. van der Kamp
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
- School of Biochemistry, Biomedical Sciences Building University of Bristol Bristol UK
| |
Collapse
|
22
|
Eken Y, Patel P, Díaz T, Jones MR, Wilson AK. SAMPL6 host-guest challenge: binding free energies via a multistep approach. J Comput Aided Mol Des 2018; 32:1097-1115. [PMID: 30225724 DOI: 10.1007/s10822-018-0159-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2018] [Accepted: 08/31/2018] [Indexed: 12/14/2022]
Abstract
In this effort in the SAMPL6 host-guest binding challenge, a combination of molecular dynamics and quantum mechanical methods were used to blindly predict the host-guest binding free energies of a series of cucurbit[8]uril (CB8), octa-acid (OA), and tetramethyl octa-acid (TEMOA) hosts bound to various guest molecules in aqueous solution. Poses for host-guest systems were generated via molecular dynamics (MD) simulations and clustering analyses. The binding free energies for the structures obtained via cluster analyses of MD trajectories were calculated using the MMPBSA method and density functional theory (DFT) with the inclusion of Grimme's dispersion correction, an implicit solvation model to model the aqueous solution, and the resolution-of-the-identity (RI) approximation (MMPBSA, RI-B3PW91-D3, and RI-B3PW91, respectively). Among these three methods tested, the results for OA and TEMOA systems showed MMPBSA and RI-B3PW91-D3 methods can be used to qualitatively rank binding energies of small molecules with an overbinding by 7 and 37 kcal/mol respectively, and RI-B3PW91 gave the poorest quality results, indicating the importance of dispersion correction for the binding free energy calculations. Due to the complexity of the CB8 systems, all of the methods tested show poor correlation with the experimental results. Other quantum mechanical approaches used for the calculation of binding free energies included DFT without the RI approximation, utilizing truncated basis sets to reduce the computational cost (memory, disk space, CPU time), and a corrected dielectric constant to account for ionic strength within the implicit solvation model.
Collapse
Affiliation(s)
- Yiğitcan Eken
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824, USA
| | - Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824, USA
| | - Thomas Díaz
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824, USA
| | - Michael R Jones
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, 48824, USA.
| |
Collapse
|
23
|
Caldararu O, Olsson MA, Misini Ignjatović M, Wang M, Ryde U. Binding free energies in the SAMPL6 octa-acid host-guest challenge calculated with MM and QM methods. J Comput Aided Mol Des 2018; 32:1027-1046. [PMID: 30203229 DOI: 10.1007/s10822-018-0158-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 08/31/2018] [Indexed: 01/15/2023]
Abstract
We have estimated free energies for the binding of eight carboxylate ligands to two variants of the octa-acid deep-cavity host in the SAMPL6 blind-test challenge (with or without endo methyl groups on the four upper-rim benzoate groups, OAM and OAH, respectively). We employed free-energy perturbation (FEP) for relative binding energies at the molecular mechanics (MM) and the combined quantum mechanical (QM) and MM (QM/MM) levels, the latter obtained with the reference-potential approach with QM/MM sampling for the MM → QM/MM FEP. The semiempirical QM method PM6-DH+ was employed for the ligand in the latter calculations. Moreover, binding free energies were also estimated from QM/MM optimised structures, combined with COSMO-RS estimates of the solvation energy and thermostatistical corrections from MM frequencies. They were performed at the PM6-DH+ level of theory with the full host and guest molecule in the QM system (and also four water molecules in the geometry optimisations) for 10-20 snapshots from molecular dynamics simulations of the complex. Finally, the structure with the lowest free energy was recalculated using the dispersion-corrected density-functional theory method TPSS-D3, for both the structure and the energy. The two FEP approaches gave similar results (PM6-DH+/MM slightly better for OAM), which were among the five submissions with the best performance in the challenge and gave the best results without any fit to data from the SAMPL5 challenge, with mean absolute deviations (MAD) of 2.4-5.2 kJ/mol and a correlation coefficient (R2) of 0.77-0.93. This is the first time QM/MM approaches give binding free energies that are competitive to those obtained with MM for the octa-acid host. The QM/MM-optimised structures gave somewhat worse performance (MAD = 3-8 kJ/mol and R2 = 0.1-0.9), but the results were improved compared to previous studies of this system with similar methods.
Collapse
Affiliation(s)
- Octav Caldararu
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Martin A Olsson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Majda Misini Ignjatović
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Meiting Wang
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|