Abstract:In the current work, an empirical solvation model based on accessible surface areas was reported, which can be used to predict the solvation free energies for both organic and biological molecules very fast. This solvation model is based on atom-weighted solvent accessible surface area (SAWSA). The parameterization procedure for different kinds of atoms was performed as follows: first, the atoms in a molecule were defined to different atom types based on SMARTS language; then the solvent accessible surface area for each atom (or charged group) was calculated; finally, 4 genetic algorithm (GA) was applied to optimize the solvation parameters for different atom types in order to reproduce the experimental solvation free energies. The derived model possessed promising predictive ability as indicated by the good statistical significance and good prediction on the external test set. Using the solvation model based on all 377 neutral molecules, we have achieved an average unsigned error of 0:51 kcal/mol and standard deviation of 0.46 kcal/mol, which was better than the model proposed by Wang et al. The solvation model developed in the current work was applied to predict the solvation free energies of small organic molecules and proteins. For the,51 small organic molecules, the SAWSA model could give consistent results with the AM1/SM2.1 model in addition to several molecules with large conjugate systems. Moreover, the predictions from SAWSA were much better than those from SM5.0R, a solvation model based on geometry-dependent atom surface tensions. For the IS proteins randomly selected from the Brookhaven PDB database, the solvation free energies predicted by the SAWSA model showed high linear correlation (r = 0.99) than those predicted by PBSA, which were much better than those given by the Ooi model and the Vial model. Finally, we have successfully applied this model to predict the relative binding free energies for four binding modes of EGFR/quinazoline. The most favorable binding mode identified by MM-PBSA could also be correctly recognized by MM-SAWSA. The relative solvation free energies calculated by SAWSA show obvious correlation with those calculated by PBSA. The SAWSA should have potential applications in QSAR, molecular docking, protein folding, free energy calculations, and so forth.

Electrostatic Interactions in Protein Molecules: A Novel Model for the Calculation of the Intrinsic Pk a of Lysozyme

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

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

Development of molecular simulation methods to accurately represent protein-surface interactions: Method assessment for the calculation of electrostatic effects

Accurate p K a Calculations in Proteins with Reactive Molecular Dynamics Provide Physical Insight Into the Electrostatic Origins of Their Values

Calculating Ph-Dependent Free Energy of Proteins by Using Monte Carlo Protonation Probabilities of Ionizable Residues

A Statistical Approach To The Prediction Of Pk(A) Values In Proteins

Continuum Electrostatics Approaches to Calculating p$K_a$s and $E_m$s in Proteins

Monte Carlo Protein Folding: Simulations of Met-Enkephalin with Solvent-Accessible Area Parameterizations

Electrostatic interactions in concentrated protein solutions

Charge Interactions in a Highly Charge-depleted Protein

Protein electrostatics: a review of the equations and methods used to model electrostatic equations in biomolecules--applications in biotechnology

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

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

Advances in Molecular Modeling of Ion-Protein Interaction Systems Towards Accurate Electrostatics: Methods and Applications

An effective method of electrostatic potential calculation on protein binding site

Exploring the pH Sensitivity of Ion-Pair Interactions on a Self-Assembled Monolayer by Scanning Electrochemical Microscopy

Empirical Aqueous Solvation Model for Protein Based on Accessible Surface Areas

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

Electrostatics, structure prediction, and the energy landscapes for protein folding and binding

Electrostatic solvation energy for two oppositely charged ions in a solvated protein system: salt bridges can stabilize proteins.