No abstract

Kinetic data on a number of protein-protein associations have provided evidence for the initial formation of a pre-equilibrium encounter complex that subsequently relaxes to the final stereospecific complex. Site-directed mutagenesis and brownian dynamics simulations have suggested that the rate of association can be modulated by perturbations in charge distribution outside the direct interaction surfaces. Furthermore, rate enhancement through non-specific binding may occur by either a reduction in dimensionality or the presence of a short-range, non-specific attractive potential. Here, using paramagnetic relaxation enhancement, we directly demonstrate the existence and visualize the distribution of an ensemble of transient, non-specific encounter complexes under equilibrium conditions for a relatively weak protein-protein complex between the amino-terminal domain of enzyme I and the phosphocarrier protein HPr. Neither the stereospecific complex alone nor any single alternative conformation can account fully for the intermolecular paramagnetic relaxation enhancement data. Restrained rigid-body simulated annealing refinement against the paramagnetic relaxation enhancement data enables us to obtain an atomic probability distribution map of the non-specific encounter complex ensemble that qualitatively correlates with the electrostatic surface potentials on the interacting proteins. Qualitatively similar results are presented for two other protein-protein complexes.

Paramagnetic relaxation enhancement (PRE) measurements on (1)H nuclei have the potential to play an important role in NMR structure determination of macromolecules by providing unique long-range (10-35 A) distance information. Recent methodological advances for covalently attaching paramagnetic groups at specific sites on both proteins and nucleic acids have permitted the application of the PRE to various biological macromolecules. However, because artificially introduced paramagnetic groups are exposed to solvent and linked to the macromolecule by several freely rotatable bonds, they are intrinsically flexible. This renders conventional back-calculation of the (1)H-PRE using a single-point representation inaccurate, thereby severely limiting the utility of the (1)H-PRE as a tool for structure refinement. To circumvent these limitations, we have developed a theoretical framework and computational strategy with which to accurately back-calculate (1)H-PREs arising from flexible paramagnetic groups attached to macromolecules. In this scheme, the (1)H-PRE is calculated using a modified Solomon-Bloembergen equation incorporating a "model-free" formalism, based on a multiple-structure representation of the paramagnetic group in simulated annealing calculations. The ensemble approach for (1)H-PRE back-calculation was examined using several SRY/DNA complexes incorporating dT-EDTA-Mn(2+) at three distinct sites in the DNA, permitting a large data set comprising 435 experimental backbone and side-chain (1)H-PREs to be obtained in a straightforward manner from 2D through-bond correlation experiments. Calculations employing complete cross-validation demonstrate that the ensemble representation provides a means to accurately utilize backbone and side-chain (1)H-PRE data arising from a flexible paramagnetic group in structure refinement. The results of (1)H-PRE based refinement, in conjunction with previously obtained NMR restraints, indicate that significant gains in accuracy can be readily obtained. This is particularly significant in the case of macromolecular complexes where intermolecular translational restraints derived from nuclear Overhauser enhancement data may be limited.