Partition coefficients of methylated DNA bases obtained from free energy calculations with molecular electron density derived atomic charges

2018 ◽  
Vol 39 (22) ◽  
pp. 1728-1737 ◽  
Author(s):  
A. Lara ◽  
M. Riquelme ◽  
E. Vöhringer-Martinez
Keyword(s):  
Free Energy ◽  
Electron Density ◽  
Partition Coefficients ◽  
Free Energy Calculations ◽  
Dna Bases ◽  
Atomic Charges ◽  
Energy Calculations ◽  
Methylated Dna ◽  
Molecular Electron

Partition Coefficients of Methylated DNA Bases Obtained from Free Energy Calculations with Molecular Electron Density Derived Atomic Charges.

2018 ◽  
Author(s):  
Alejandro Lara ◽  
Maximiliano Riquelme ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Electron Density ◽  
Partition Coefficients ◽  
Dna Bases ◽  
Atomic Charges ◽  
Free Energies ◽  
Solvation Model ◽  
Minimal Basis ◽  
Methylated Dna ◽  
Implicit Solvation Model ◽  
Good Agreement

<div> <div> <div> <p>Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have been also used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in vac- uum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform-water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER-99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account deriving the atomic charges of polar DNA bases and when the energy needed to polarize the electron den- sity of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogues. Comparison of the two partitioning methods Hirsheld-I and Minimal Basis Iterative Stockholder (MBIS) revealed some deficiencies in the Hirshfeld-I method related to nonexistent isolated anionic nitrogen pro-atoms used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. </p> </div> </div> </div>

SAMPL6 Octanol-Water Partition Coefficients from Alchemical Free Energy Calculations with MBIS Atomic Charges

2019 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Free Energy ◽  
Electron Density ◽  
Free Energy Calculations ◽  
Basis Set ◽  
Energetic Cost ◽  
Atomic Charges ◽  
Free Energies ◽  
Energy Calculations ◽  
Alchemical Free Energy ◽  
Alchemical Free Energy Calculations

In molecular modeling the description of the interactions between molecules forms the basis for a correct prediction of macroscopic observables. Here, we derive atomic charges from the implicitly polarized electron density of eleven molecules in the SAMPL6 challenge using the Hirshfeld-I and Minimal Basis Set Iterative Stockholder(MBIS) partitioning method. These atomic charges combined with other parameters in the GAFF force field and different water/octanol models were then used in alchemical free energy calculations to obtain hydration and solvation free energies, which after correction for the polarization cost, result in the blind prediction of the partition coefficient. From the tested partitioning methods and water models the S-MBIS atomic charges with the TIP3P water model presented the smallest deviation from the experiment. Conformational dependence of the free energies and the energetic cost associated with the polarization of the electron density are discussed.

SAMPL6 Octanol-Water Partition Coefficients from Alchemical Free Energy Calculations with MBIS Atomic Charges

2019 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Free Energy ◽  
Electron Density ◽  
Free Energy Calculations ◽  
Basis Set ◽  
Energetic Cost ◽  
Atomic Charges ◽  
Free Energies ◽  
Energy Calculations ◽  
Alchemical Free Energy ◽  
Alchemical Free Energy Calculations

In molecular modeling the description of the interactions between molecules forms the basis for a correct prediction of macroscopic observables. Here, we derive atomic charges from the implicitly polarized electron density of eleven molecules in the SAMPL6 challenge using the Hirshfeld-I and Minimal Basis Set Iterative Stockholder(MBIS) partitioning method. These atomic charges combined with other parameters in the GAFF force field and different water/octanol models were then used in alchemical free energy calculations to obtain hydration and solvation free energies, which after correction for the polarization cost, result in the blind prediction of the partition coefficient. From the tested partitioning methods and water models the S-MBIS atomic charges with the TIP3P water model presented the smallest deviation from the experiment. Conformational dependence of the free energies and the energetic cost associated with the polarization of the electron density are discussed.

Partition Coefficients of Methylated DNA Bases Obtained from Free Energy Calculations with Molecular Electron Density Derived Atomic Charges.

2018 ◽  
Author(s):  
Alejandro Lara ◽  
Maximiliano Riquelme ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Electron Density ◽  
Partition Coefficients ◽  
Dna Bases ◽  
Atomic Charges ◽  
Free Energies ◽  
Solvation Model ◽  
Minimal Basis ◽  
Methylated Dna ◽  
Implicit Solvation Model ◽  
Good Agreement

<div> <div> <div> <p>Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have been also used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in vac- uum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform-water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER-99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account deriving the atomic charges of polar DNA bases and when the energy needed to polarize the electron den- sity of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogues. Comparison of the two partitioning methods Hirsheld-I and Minimal Basis Iterative Stockholder (MBIS) revealed some deficiencies in the Hirshfeld-I method related to nonexistent isolated anionic nitrogen pro-atoms used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. </p> </div> </div> </div>

Free-Energy Calculations Reveal the Subtle Differences in the Interactions of DNA Bases with α-Hemolysin

2015 ◽  
Vol 11 (2) ◽  
pp. 810-816 ◽  
Author(s):  
Richard M. A. Manara ◽  
Andrew T. Guy ◽  
E. Jayne Wallace ◽  
Syma Khalid
Keyword(s):  
Free Energy ◽  
Free Energy Calculations ◽  
Dna Bases ◽  
Energy Calculations

scholarly journals Molecular Environment Specific Atomic Charges Improve Binding Affinity Predictions of SAMPL5 Host-Guest Systems

2021 ◽  
Author(s):  
Duvan Gonzalez ◽  
Luis Macaya ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Free Energy ◽  
Force Field ◽  
Binding Affinity ◽  
Free Energy Calculations ◽  
Atomic Charges ◽  
Binding Affinity Prediction ◽  
Energy Calculations ◽  
Affinity Prediction ◽  
Alchemical Free Energy ◽  
Alchemical Free Energy Calculations

<div> <div> <div> <p>Host-guest systems are widely used in benchmarks as model systems to improve computational methods for absolute binding free energy predictions. Recent advances in sampling algorithms for alchemical free energy calculations and the increase in computational power have made their binding affinity prediction primarily dependent on the quality of the force field. Here, we propose a new methodology to derive the atomic charges of host-guest systems based on QM/MM calculations and the MBIS partitioning of the polarized electron density. A newly developed interface between the OpenMM and ORCA software package provides D-MBIS charges that best represent the guest’s average electrostatic interactions in the hosts or the solvent. The simulation workflow also calculates the average energy required to polarize the guest in the bound and unbound state. Alchemical free energy calculations using the GAFF force field parameters with D-MBIS charges improve the binding affinity prediction of six guests bound to two octa-acid hosts compared to the AM1-BCC charge set after correction with the average energetic polarization cost. This correction results from the difference in the energetic polarization cost between the bound and unbound state and contributes significantly to the binding affinity of anionic guests. </p></div></div></div><div><div><div> </div> </div> </div>

scholarly journals Molecular Environment Specific Atomic Charges Improve Binding Affinity Predictions of SAMPL5 Host-Guest Systems

2021 ◽  
Author(s):  
Duvan Gonzalez ◽  
Luis Macaya ◽  
Esteban Vöhringer-Martinez
Keyword(s):  
Free Energy ◽  
Force Field ◽  
Binding Affinity ◽  
Free Energy Calculations ◽  
Atomic Charges ◽  
Binding Affinity Prediction ◽  
Energy Calculations ◽  
Affinity Prediction ◽  
Alchemical Free Energy ◽  
Alchemical Free Energy Calculations

<div> <div> <div> <p>Host-guest systems are widely used in benchmarks as model systems to improve computational methods for absolute binding free energy predictions. Recent advances in sampling algorithms for alchemical free energy calculations and the increase in computational power have made their binding affinity prediction primarily dependent on the quality of the force field. Here, we propose a new methodology to derive the atomic charges of host-guest systems based on QM/MM calculations and the MBIS partitioning of the polarized electron density. A newly developed interface between the OpenMM and ORCA software package provides D-MBIS charges that best represent the guest’s average electrostatic interactions in the hosts or the solvent. The simulation workflow also calculates the average energy required to polarize the guest in the bound and unbound state. Alchemical free energy calculations using the GAFF force field parameters with D-MBIS charges improve the binding affinity prediction of six guests bound to two octa-acid hosts compared to the AM1-BCC charge set after correction with the average energetic polarization cost. This correction results from the difference in the energetic polarization cost between the bound and unbound state and contributes significantly to the binding affinity of anionic guests. </p></div></div></div><div><div><div> </div> </div> </div>

Automated Assessment of Binding Affinity via Alchemical Free Energy Calculations

2020 ◽  
Author(s):  
Maximilian Kuhn ◽  
Stuart Firth-Clark ◽  
Paolo Tosco ◽  
Antonia S. J. S. Mey ◽  
Mark Mackey ◽  
...  
Keyword(s):  
Free Energy ◽  
Protein Structures ◽  
Free Energy Calculations ◽  
Distinct Advantage ◽  
Binding Affinities ◽  
Structure Based Drug Design ◽  
Energy Calculations ◽  
Kendall’S Τ ◽  
Experimental Values ◽  
Better Than

Free energy calculations have seen increased usage in structure-based drug design. Despite the rising interest, automation of the complex calculations and subsequent analysis of their results are still hampered by the restricted choice of available tools. In this work, an application for automated setup and processing of free energy calculations is presented. Several sanity checks for assessing the reliability of the calculations were implemented, constituting a distinct advantage over existing open-source tools. The underlying workflow is built on top of the software Sire, SOMD, BioSimSpace and OpenMM and uses the AMBER14SB and GAFF2.1 force fields. It was validated on two datasets originally composed by Schrödinger, consisting of 14 protein structures and 220 ligands. Predicted binding affinities were in good agreement with experimental values. For the larger dataset the average correlation coefficient Rp was 0.70 ± 0.05 and average Kendall’s τ was 0.53 ± 0.05 which is broadly comparable to or better than previously reported results using other methods. <br>

Reaction-based Enumeration, Active Learning, and Free Energy Calculations to Rapidly Explore Synthetically Tractable Chemical Space and Optimize Potency of Cyclin Dependent Kinase 2 Inhibitors

2019 ◽  
Author(s):  
Kyle Konze ◽  
Pieter Bos ◽  
Markus Dahlgren ◽  
Karl Leswing ◽  
Ivan Tubert-Brohman ◽  
...  
Keyword(s):  
Free Energy ◽  
Drug Discovery ◽  
Active Learning ◽  
Large Scale ◽  
Chemical Space ◽  
Population Based ◽  
Free Energy Calculations ◽  
Computational Technique ◽  
Cyclin Dependent Kinase ◽  
Energy Calculations

We report a new computational technique, PathFinder, that uses retrosynthetic analysis followed by combinatorial synthesis to generate novel compounds in synthetically accessible chemical space. Coupling PathFinder with active learning and cloud-based free energy calculations allows for large-scale potency predictions of compounds on a timescale that impacts drug discovery. The process is further accelerated by using a combination of population-based statistics and active learning techniques. Using this approach, we rapidly optimized R-groups and core hops for inhibitors of cyclin-dependent kinase 2. We explored greater than 300 thousand ideas and identified 35 ligands with diverse commercially available R-groups and a predicted IC<sub>50</sub> < 100 nM, and four unique cores with a predicted IC<sub>50</sub> < 100 nM. The rapid turnaround time, and scale of chemical exploration, suggests that this is a useful approach to accelerate the discovery of novel chemical matter in drug discovery campaigns.

Sign in / Sign up

Export Citation Format

Share Document