A rapid solvent accessible surface area estimator for coarse grained molecular simulations

The rapid and accurate calculation of solvent accessible surface area (SASA) is extremely useful in the energetic analysis of biomolecules. For example, SASA models can be used to estimate the transfer free energy associated with biophysical processes, and when combined with coarse‐grained simulations, can be particularly useful for accounting for solvation effects within the framework of implicit solvent models. In such cases, a fast and accurate, residue‐wise SASA predictor is highly desirable. Here, we develop a predictive model that estimates SASAs based on Cα ‐only protein structures. Through an extensive comparison between this method and a comparable method, POPS‐R, we demonstrate that our new method, Protein‐C α Solvent Accessibilities or PCASA, shows better performance, especially for unfolded conformations of proteins. We anticipate that this model will be quite useful in the efficient inclusion of SASA‐based solvent free energy estimations in coarse‐grained protein folding simulations. PCASA is made freely available to the academic community at https://github.com/atfrank/PCASA. © 2017 Wiley Periodicals, Inc.

[1]  C. Brooks,et al.  Recent advances in the development and application of implicit solvent models in biomolecule simulations. , 2004, Current opinion in structural biology.

[2]  D Fischer,et al.  Molecular surface representations by sparse critical points , 1994, Proteins.

[3]  C. Tanford.  Protein denaturation. , 1968, Advances in protein chemistry.

[4]  Franca Fraternali,et al.  POPS: a fast algorithm for solvent accessible surface areas at atomic and residue level , 2003, Nucleic Acids Res..

[5]  A. Shrake,et al.  Environment and exposure to solvent of protein atoms. Lysozyme and insulin. , 1973, Journal of molecular biology.

[6]  Ján Plávka,et al.  ARVO: A Fortran package for computing the solvent accessible surface area and the excluded volume of overlapping spheres via analytic equations , 2005, Comput. Phys. Commun..

[7]  J. Janin,et al.  Analytical approximation to the accessible surface area of proteins. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[8]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[9]  Michael L. Connolly,et al.  Computation of molecular volume , 1985 .

[10]  Vinayak Kulkarni,et al.  Exhaustive prediction of disease susceptibility to coding base changes in the human genome , 2008, BMC Bioinformatics.

[11]  T. Richmond,et al.  Solvent accessible surface area and excluded volume in proteins. Analytical equations for overlapping spheres and implications for the hydrophobic effect. , 1984, Journal of molecular biology.

[12]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[13]  Franca Fraternali,et al.  Parameter optimized surfaces (POPS): analysis of key interactions and conformational changes in the ribosome , 2002, Nucleic Acids Res..

[14]  M. L. Connolly.  Analytical molecular surface calculation , 1983 .

[15]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[16]  Andrew D. Martin,et al.  MCMCpack: Markov chain Monte Carlo in R , 2011 .

[17]  D. Thirumalai,et al.  Molecular origin of constant m-values, denatured state collapse, and residue-dependent transition midpoints in globular proteins. , 2009, Biochemistry.

[18]  John Karanicolas,et al.  Improved Gō-like models demonstrate the robustness of protein folding mechanisms towards non-native interactions. , 2003, Journal of molecular biology.

[19]  W. C. Still,et al.  A rapid approximation to the solvent accessible surface areas of atoms , 1988 .

[20]  Hamid Pezeshk,et al.  Impact of residue accessible surface area on the prediction of protein secondary structures , 2008, BMC Bioinformatics.

[21]  W. C. Still,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999 .

[22]  Guy Ziv,et al.  Protein folding, protein collapse, and tanford's transfer model: lessons from single-molecule FRET. , 2009, Journal of the American Chemical Society.

[23]  Michel Petitjean,et al.  On the analytical calculation of van der Waals surfaces and volumes: Some numerical aspects , 1994, J. Comput. Chem..

[24]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[25]  Wolfgang Wenzel,et al.  Derivatives of molecular surface area and volume: Simple and exact analytical formulas , 2011, J. Comput. Chem..

[26]  Bernard R Brooks,et al.  Effects of denaturants and osmolytes on proteins are accurately predicted by the molecular transfer model , 2008, Proceedings of the National Academy of Sciences.