Modeling ion-ion interaction in proteins: A molecular dynamics free energycalculation of the guanidinium-acetate association

The twin nitrogen-twin oxygen, C-2v association of guanidinium and acetate in water is examined by means of molecular dynamics free energy calculation s, using three approaches for handling Coulomb electrostatic interactions: (i) an Ewald lattice summation, (ii) a generalized reaction field correctio n, and, (iii) a smoothed spherical truncation. The potential of mean force obtained from the Ewald simulation exhibits the expected minima characteris tic of a contact and a solvent-separated ion pair. In contrast, the repulsi ve and, therefore, physically unrealistic profile resulting from a spherica l truncation of the electrostatic interactions at 12 Angstrom calls into qu estion the validity of such a protocol for simulating charged proteins. The introduction of a generalized reaction field improves the description of t he solution significantly, in spite of an artifactual behavior towards the edge of the cutoff sphere. Whereas this approach, compared to a conventiona l spherical truncation, implies virtually no additional computational effor t, it rapidly becomes costly for large ion-water and water-water cutoff rad ii. At 12 Angstrom, the computational investment is comparable to that of a n Ewald simulation with an appropriately chosen number of k-vectors and sep aration parameter, making the latter more cost-effective. (C) 2000 American Institute of Physics. [S0021-9606(00)51422-6].