Abstract:Accurate estimation of solvation free energy (SFE) lays the foundation for accurate prediction of binding free energy. The Poisson-Boltzmann (PB) or generalized Born (GB) combined with surface area (SA) continuum solvation method (PBSA and GBSA) have been widely used in SFE calculations because they can achieve good balance between accuracy and efficiency. However, the accuracy of these methods can be affected by several factors such as the charge models, polar and nonpolar SFE calculation methods and the atom radii used in the calculation. In this work, the performance of the ABCG2 (AM1-BCC-GAFF2) charge model as well as other two charge models, that is, RESP (Restrained Electrostatic Potential) and AM1-BCC (Austin Model 1-bond charge corrections), on the SFE prediction of 544 small molecules in water by PBSA/GBSA was evaluated. In order to improve the performance of the PBSA prediction based on the ABCG2 charge, we further explored the influence of atom radii on the prediction accuracy and yielded a set of atom radius parameters for more accurate SFE prediction using PBSA based on the ABCG2/GAFF2 by reproducing the thermodynamic integration (TI) calculation results. The PB radius parameters of carbon, oxygen, sulfur, phosphorus, chloride, bromide and iodine, were adjusted. New atom types, on, oi, hn1, hn2, hn3, were introduced to further improve the fitting performance. Then, we tuned the parameters in the nonpolar SFE model using the experimental SFE data and the PB calculation results. By adopting the new radius parameters and new nonpolar SFE model, the root mean square error (RMSE) of the SFE calculation for the 544 molecules decreased from 2.38 to 1.05 kcal/mol. Finally, the new radius parameters were applied in the prediction of protein-ligand binding free energies using the MM-PBSA method. For the eight systems tested, we could observe higher correlation between the experiment data and calculation results and smaller prediction errors for the absolute binding free energies, demonstrating that our new radius parameters can improve the free energy calculation using the MM-PBSA method.

Accurate, Robust, And Reliable Calculations Of Poisson-Boltzmann Binding Energies

Accurate, robust and reliable calculations of Poisson-Boltzmann solvation energies

Poisson-Boltzmann model for protein-surface electrostatic interactions and grid-convergence study using the PyGBe code

Poisson-Boltzmann based machine learning (PBML) model for electrostatic analysis

Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Accurate boundary-integral formulations for the calculation of electrostatic forces with an implicit-solvent model

Incorporating Dipolar Solvents with Variable Density in Poisson-Boltzmann Electrostatics

Development and test of highly accurate endpoint free energy methods. 1: Evaluation of ABCG2 charge model on solvation free energy prediction and optimization of atom radii suitable for more accurate solvation free energy prediction by the PBSA method

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

How Well Does the Extended Linear Interaction Energy Method Perform in Accurate Binding Free Energy Calculations?

Improving the Performance of MM/PBSA in Protein-Protein Interactions Via the Screening Electrostatic Energy.

Assessing the performance of the MM/PBSA and MM/GBSA methods. 10. Impacts of enhanced sampling and variable dielectric model on protein-protein Interactions

On delivering polar solvation free energy of proteins from energy minimized structures using a regularized super‐Gaussian Poisson–Boltzmann model

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

Optimized Parallelization of Boundary Integral Poisson-Boltzmann Solvers

MM/PB(GB)SA benchmarks on soluble proteins and membrane proteins

How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis.

End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design

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