Dataset Information

Optimized Lennard-Jones Parameters for Druglike Small Molecules.

ABSTRACT: Meaningful efforts in computer-aided drug design (CADD) require accurate molecular mechanical force fields to quantitatively characterize protein-ligand interactions, ligand hydration free energies, and other ligand physical properties. Atomic models of new compounds are commonly generated by analogy from the predefined tabulated parameters of a given force field. Two widely used approaches following this strategy are the General Amber Force Field (GAFF) and the CHARMM General Force Field (CGenFF). An important limitation of using pretabulated parameter values is that they may be inadequate in the context of a specific molecule. To resolve this issue, we previously introduced the General Automated Atomic Model Parameterization (GAAMP) for automatically generating the parameters of atomic models of small molecules, using the results from ab initio quantum mechanical (QM) calculations as target data. The GAAMP protocol uses QM data to optimize the bond, valence angle, and dihedral angle internal parameters, and atomic partial charges. However, since the treatment of van der Waals interactions based on QM is challenging and may often be unreliable, the Lennard-Jones 6-12 parameters are kept unchanged from the initial atom types assignments (GAFF or CGenFF), which limits the accuracy that can be achieved by these models. To address this issue, a new set of Lennard-Jones 6-12 parameters was systematically optimized to reproduce experimental neat liquid densities and enthalpies of vaporization for a large set of 430 compounds, covering a wide range of chemical functionalities. Calculations of the hydration free energy indicate that optimal accuracy for these models is achieved when the molecule-water van der Waals dispersion is rescaled by a factor of 1.115. The final optimized model yields an average unsigned error of 0.79 kcal/mol in the hydration free energies.

SUBMITTER: Boulanger E

PROVIDER: S-EPMC5997559 | biostudies-literature | 2018 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Optimized Lennard-Jones Parameters for Druglike Small Molecules.

Boulanger Eliot E Huang Lei L Rupakheti Chetan C MacKerell Alexander D AD Roux Benoît B

Journal of chemical theory and computation 20180507 6

Meaningful efforts in computer-aided drug design (CADD) require accurate molecular mechanical force fields to quantitatively characterize protein-ligand interactions, ligand hydration free energies, and other ligand physical properties. Atomic models of new compounds are commonly generated by analogy from the predefined tabulated parameters of a given force field. Two widely used approaches following this strategy are the General Amber Force Field (GAFF) and the CHARMM General Force Field (CGenF ...[more]

PMID: 29694035

Similar Datasets

Project description:Molecular dynamics (MD) simulations based on atomic models play an important role in the drug-discovery process to screen molecules, estimate binding free energies, and optimize lead compounds in chemical space. Accurate computations of thermodynamic and kinetic properties using MD simulations are highly dependent on the accuracy of the underlying atomic force field. In this context, going beyond the nonpolarizable fixed-charge model by accounting explicitly for induced polarization is highly desirable. The CHARMM polarizable force field based on classical Drude oscillators, in which an auxiliary charged particle is attached via a harmonic spring to its parent nucleus, offers both a computationally convenient and rigorous framework to model explicitly induced electronic polarization in MD simulations. For any molecule of interest, electrostatic partial charges, atomic polarizabilities, and Thole shielding factors, as well as bonded parameters can either be determined from ab initio calculations or ascribed from the knowledge-based library of the CHARMM Generalized force field. While this approach is fairly reliable in general, it is well understood that the overall accuracy of the models with respect to thermodynamic properties such as bulk density, enthalpies, and solvation free energies is particularly sensitive to the nonbonded Lennard-Jones (LJ) parameters. In the present study, we systematically refined the set of LJ parameters for the atom types available in the Drude force field to best match the experimental thermodynamic properties for 416 small druglike organic molecules. To further test the transferability of the optimized parameters, the hydration free energy of 372 molecules was computed. The calculations resulted in a small average error of 0.46 kcal/mol and a Pearson R of 0.9, representing a significant improvement over the additive GAFF force field in our previous study, where an average error of ∼2 kcal/mol was obtained. Such an improvement is consistent with the ability of the polarizable Drude model to more accurately model interactions in different environments. The effort provides a roadmap for the global optimization of force field parameters using experimental data. It is hoped that the present effort will further the application of the Drude polarizable force field in molecular simulations including drug design and discovery.

Project description:Metal ions play significant roles in biological systems. Accurate molecular dynamics (MD) simulations on these systems require a validated set of parameters. Although there are more detailed ways to model metal ions, the nonbonded model, which employs a 12-6 Lennard-Jones (LJ) term plus an electrostatic potential is still widely used in MD simulations today due to its simple form. However, LJ parameters have limited transferability due to different combining rules, various water models and diverse simulation methods. Recently, simulations employing a Particle Mesh Ewald (PME) treatment for long-range electrostatics have become more and more popular owing to their speed and accuracy. In the present work we have systematically designed LJ parameters for 24 +2 metal (M(II)) cations to reproduce different experimental properties appropriate for the Lorentz-Berthelot combining rules and PME simulations. We began by testing the transferability of currently available M(II) ion LJ parameters. The results showed that there are differences between simulations employing Ewald summation with other simulation methods and that it was necessary to design new parameters specific for PME based simulations. Employing the thermodynamic integration (TI) method and performing periodic boundary MD simulations employing PME, allowed for the systematic investigation of the LJ parameter space. Hydration free energies (HFEs), the ion-oxygen distance in the first solvation shell (IOD) and coordination numbers (CNs) were obtained for various combinations of the parameters of the LJ potential for four widely used water models (TIP3P, SPC/E, TIP4P and TIP4PEW). Results showed that the three simulated properties were highly correlated. Meanwhile, M(II) ions with the same parameters in different water models produce remarkably different HFEs but similar structural properties. It is difficult to reproduce various experimental values simultaneously because the nonbonded model underestimates the interaction between the metal ions and water molecules at short range. Moreover, the extent of underestimation increases successively for the TIP3P, SPC/E, TIP4PEW and TIP4P water models. Nonetheless, we fitted a curve to describe the relationship between ? (the well depth) and radius (Rmin/2) from experimental data on noble gases to facilitate the generation of the best possible compromise models. Hence, by targeting different experimental values, we developed three sets of parameters for M(II) cations for three different water models (TIP3P, SPC/E and TIP4PEW). These parameters we feel represent the best possible compromise that can be achieved using the nonbonded model for the ions in combination with simple water models. From a computational uncertainty analysis we estimate that the uncertainty in our computed HFEs is on the order of ±1kcal/mol. Further improvements will require more advanced non-bonded models likely with inclusion of polarization.

Dataset Information

Optimized Lennard-Jones Parameters for Druglike Small Molecules.

Publications

Optimized Lennard-Jones Parameters for Druglike Small Molecules.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets