Energetics of sequence-specific protein-DNA association: computational analysis of integrase Tn916 binding to its...

The N-terminal domain of the bacterial integrase Tn916 specifically recognizes the 11 bp DNA target site by positioning the face of a three-stranded beta-sheet into the major groove. Binding is linked to structural adaptation. We have characterized INT-DBD binding to DNA in detail by calorimetry [Milev, S., Gorfe, A., Karshikoff, A., Clubb, R. T., Bosshard, H. R., and Jelesarov, I. (2003) Biochemistry 42, 3481-3491]. Our thermodynamic analysis has indicated that the major driving force of association is the hydrophobic effect while polar interactions contribute less. To gain more comprehensive information about the binding process, we performed a computational analysis of the binding free energy and report here the results. A hybrid molecular mechanics/continuum approach was followed. The total binding free energy is predicted with reasonable accuracy. The calculations confirm that nonpolar effects stabilize the protein-DNA complex while electrostatics opposes binding. Structural changes optimizing surface complementarity are costly in terms of energy. The energetic consequences from the replacement of nine DNA-contacting residues by alanine were investigated. The calculations correctly predict the binding affinity decrease of eight mutations and the destabilizing effect of one wild-type residue. Bulky side chains stabilize the wild-type complex through packing interactions and favorable nonpolar dehydration, but the net nonpolar energy changes do not correlate with the relative affinity loss upon mutation. Discrete protein-DNA electrostatic interactions may be net stabilizing or net destabilizing depending on the local environment. In contrast to nonpolar energy changes, the magnitude of the electrostatic free energy ranks the mutations according to the experimentally measured DeltaDeltaG. Free energy decomposition analysis from a structural perspective leads to detailed information about the thermodynamic strategy used by INT-DBD for sequence-specific DNA binding.