A THEORETICAL-STUDY OF POLYELECTROLYTE EFFECTS IN PROTEIN DNA INTERACTIONS - MONTE-CARLO FREE-ENERGY SIMULATIONS ON THE ION ATMOSPHERE CONTRIBUTION TO THE THERMODYNAMICS OF LAMBDA REPRESSOR OPERATOR COMPLEX-FORMATION

We report herein a theoretical calculation of the ion atmosphere contribution to the free energy of association for a protein-DNA complex based on Monte Carlo computer simulations and thermodynamic perturbation theory. The system considered is the dimer of the amino-terminal fragment of the lambda-cI repressor in a complex with a 17 base pair oligonucleotide of DNA, based on the crystal structure of Pabo and Sauer. Only the movements of the small ions (sodium and chloride ions) are considered explicitly, with solvent water modeled as a dielectric continuum (a ''primitive model''). The free energies are determined as a function of both distance of separation between the protein and the DNA, each of which is fixed in its respective crystal geometry, and ionic strength at a temperature of 298 K. Results of our simulations indicate that the ion atmosphere contribution to the free energy of association is favorable only at short distances of separation and is at a maximum when the protein approaches the DNA from a distance of approximately 7 angstrom. This distance corresponds to the radius of the shroud of condensed counterions around B-DNA. At larger distances of separation between the protein and the DNA, the uncondensed diffuse ionic cloud opposes complexation. The effect known as ''counterion release'' in the context of DNA-ligand association appears to be short-ranged and a property of the condensed counterions only.