Re: [AMBER] setting ionic strength in leap

From: George M. Giambasu <>
Date: Sat, 7 Jul 2018 00:34:45 -0400

Here are few further points to consider when trying to reach a desired
salt concentration in nucleic acids MD simulations:

- the resulting salt concentration after equilibration will be different
from any initial guess, and so, in order to reach a desired
concentration you will want to iteratively change the composition of
your system and equilibrate (here is an example:

- equilibration of ions and water takes a long time: tens of nanoseconds
for simple helical nucleic acids and at least one order of magnitude
more for cases where ions might bind in cavities. is a good starting point.

- imposing specific solution salt conditions refers to maintaining a
specific *bulk* salt concentration, that is the ion concentration in
regions remote from the nucleic acid solute where the local density of
ions is not perturbed by the presence of the nucleic acid.

- typical simulation boxes (having clearances of 10-15 Ang.) are not
able to model a solution bulk, as one has to go *at least* 25 Angstroms
away from the nucleic acid surface to reach the bulk. So, if you want
your your simulation to reflect specific salt conditions (that can be
replicated experimentally) you better choose a large enough box (see
more here: ,



Dr. George M. Giambasu, Ph.D.
Rutgers, The State University of New Jersey
174 Frelinghuysen Road, Room 308E Piscataway,
NJ 08854, USA

AMBER mailing list
Received on Fri Jul 06 2018 - 22:00:01 PDT
Custom Search