G Br 6 NL : A generalized Born method for accurately reproducing solvation energy of the nonlinear Poisson-Boltzmann equation

Abstract

The nonlinear Poisson-Boltzmann (NLPB) equation can provide accurate modeling of electrostatic effects for nucleic acids and highly charged proteins. Generalized Born methods have been developed to mimic the linearized Poisson-Boltzmann (LPB) equation at substantially reduced cost. The computer time for solving the NLPB equation is approximately fivefold longer than for the LPB equation, thus presenting an even greater obstacle. Here we present the first generalized Born method, GBr(6)NL, for mimicking the NLPB equation. GBr(6)NL is adapted from GBr(6), a generalized Born method recently developed to reproduce the solvation energy of the LPB equation [Tjong and Zhou, J. Phys. Chem. B 111, 3055 (2007)]. Salt effects predicted by GBr(6)NL on 55 proteins overall deviate from NLPB counterparts by 0.5 kcal/mol from ionic strengths from 10 to 1000 mM, which is approximately 10% of the average magnitudes of the salt effects. GBr(6)NL predictions for the salts effects on the electrostatic interaction energies of two protein:RNA complexes are very promising.