An atomic radii set and Generalized Born implicit solvation model trained using explicit water solvation free energy data

Carlos Simmerling,Yuzhang Wang,Chuan Tian,Jorge Pincay

DOI: https://doi.org/10.26434/chemrxiv-2024-mhv7n

2024-06-10

Abstract:Compared with the other common implicit and explicit water models, Generalized Born (GB) models can provide a fast approximation of solvation free energy that is reasonably accurate but fast enough to use in molecular dynamics (MD) simulations. This enhances conformational sampling of the solute molecules, and also can be a valuable component of multi-scale simulations. We previously developed the GB-Neck2 model, which exhibited improved secondary structure balance and was used to successfully fold a series of small proteins. More recent simulations using GB-Neck2 with updated protein models suggest that α-helices remain somewhat over-stabilized. Here, we develop a more self-consistent model, retraining both the intrinsic solvation radii as well as the GB model parameters, using the solvation free energies of an explicit water model as training references. The new radii set, named MIRO, when used with the GBNSR6 implicit solvent model leads to improved reproduction of solvation free energies calculated in explicit water. The new GB-Neck3 model shows a good balance of secondary structures: the stability of β-sheets is improved, while the previously over-stabilized α-helices became less favorable, as expected. GB-Neck3 and MIRO radii should extend the range of problems accessible to biomolecular simulation.

Chemistry

What problem does this paper attempt to address?

The problems that this paper attempts to solve mainly focus on improving the Generalized Born (GB) model in the implicit solvent model to enhance its accuracy and applicability in biomolecular simulations. Specifically, the researchers are concerned with: 1. **Improving the equilibrium of α - helices and β - sheets**: Early GB models (such as GB - Neck2) had the problem of over - stabilizing α - helices when simulating proteins, which limited the accuracy of the model in predicting protein structures and dynamics. One of the goals of this paper is to improve the balance between α - helices and β - sheets by retraining the GB model parameters and the set of atomic radii. 2. **Optimizing the set of atomic radii**: Traditional sets of atomic radii (such as the Bondi radius set) may not be entirely suitable for all types of biomolecules, especially when dealing with complex systems. Therefore, the researchers developed a new set of atomic radii - MIRO, aiming to more accurately describe the behavior of different atoms in aqueous solutions. 3. **Enhancing the universality of the model**: By using the solvent - free - energy data of explicit water models (such as OPC) as a reference to train new GB model parameters and the set of atomic radii, the researchers hope to develop a general GB model that can calculate quickly and maintain high accuracy, so as to be applied to a wider range of biomolecular simulations. 4. **Reducing error compensation**: Traditional GB models often rely on error compensation with specific force fields to obtain reasonable results, but this error compensation has poor transferability between different systems. The method proposed in this paper aims to reduce the dependence on error compensation through a more self - consistent training process, thereby improving the reliability and universality of the model. In summary, the main goal of this paper is to develop a more accurate, more efficient, and widely applicable implicit solvent model by retraining the GB model parameters and the set of atomic radii, in order to improve the quality of biomolecular simulations.

An atomic radii set and Generalized Born implicit solvation model trained using explicit water solvation free energy data

A combined steepest descent and genetic algorithm (SD/GA) approach for the optimization of solvation parameters

Implicit Solvent with Explicit Ions Generalized Born Model in Molecular Dynamics: Application to DNA

Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations.

The AGBNP2 Implicit Solvation Model

Improving the Efficiency of Non-equilibrium Sampling in the Aqueous Environment via Implicit-Solvent Simulations.

Accurate and efficient generalized born model based on solvent accessibility: Derivation and application for LogP octanol/water prediction and flexible peptide docking

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains

Comparison of Implicit and Explicit Solvent Models for the Calculation of Solvation Free Energy in Organic Solvents.

Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?

Comparison of Protein Solution Structures Refined by Molecular Dynamics Simulation in Vacuum, with a Generalized Born Model, and with Explicit Water

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Optimal Dielectric Boundary for Binding Free Energy Estimates in the Implicit Solvent

A Test of Implicit Solvent Models on the Folding Simulation of the GB1 Peptide.

Parameterization and validation of solvation corrected atomic radii.

General Graph Neural Network Based Implicit Solvation Model for Organic Molecules in Water

Model for Aqueous Solvation Based on Class IV Atomic Charges and First Solvation Shell Effects

Empirical Aqueous Solvation Models Based on Accessible Surface Areas with Implicit Electrostatics

A general graph neural network based implicit solvation model for organic molecules in water

GCMe: Efficient implementation of the Gaussian core model with smeared electrostatic interactions for molecular dynamics simulations of soft matter systems

Integral equation models for solvent in macromolecular crystals