Theory and applications of the generalized Born solvation model in macromolecular simulations

Generalized Born (GB) models provide an attractive way to include some thermodynamic aspects of aqueous solvation into simulations that do not explicitly model the solvent molecules. Here we discuss our recent experience with this model, presenting in detail the way it is implemented and parallelized in the AMBER molecular modeling code. We compare results using the GB model (or GB plus a surface-area based "hydrophobic" term) to explicit solvent simulations for a 10 base-pair DNA oligomer, and for the 108-residue protein thioredoxin. A slight modification of our earlier suggested parameters makes the GB results more like those found in explicit solvent, primarily by slightly increasing the strength of NH [bond] O and NH [bond] N internal hydrogen bonds. Timing and energy stability results are reported, with an eye toward using these model for simulations of larger macromolecular systems and longer time scales.