The temperature dependences of the nanosecond dynamics of different chemical classes of amino acid residue have been analyzed by combining elastic incoherent neutron scattering experiments with molecular dynamics simulations on cytochrome P450cam. At T = 100-160 K, anharmonic motion in hydrophobic and aromatic residues is activated, whereas hydrophilic residue motions are suppressed because of hydrogen-bonding interactions. In contrast, at T = 180-220 K, water-activated jumps of hydrophilic side chains, which are strongly coupled to the relaxation rates of the hydrogen bonds they form with hydration water, become apparent. Thus, with increasing temperature, first the hydrophobic core awakens, followed by the hydrophilic surface.

The derivation of mean-square displacements from elastic incoherent neutron scattering (EINS) of proteins is examined, with the aid of experiments on camphor-bound cytochrome P450cam and complementary molecular dynamics simulations. It is shown that a q(4) correction to the elastic incoherent structure factor (where q is the scattering vector) can be simply used to reliably estimate from the experiment both the average mean-square atomic displacement, <Δr(2)> of the nonexchanged hydrogen atoms in the protein and its variance, σ(2). The molecular dynamics simulation results are in broad agreement with the experimentally derived <Δr(2)> and σ(2) derived from EINS on instruments at two different energy resolutions, corresponding to dynamics on the ∼100 ps and ∼1 ns time scales. Significant dynamical heterogeneity is found to arise from methyl-group rotations. The easy-to-apply q(4) correction extends the information extracted from elastic incoherent neutron scattering experiments and should be of wide applicability.

A refined model for the solution structure of oxidized putidaredoxin (Pdxo), a Cys4Fe2S2 ferredoxin, has been determined. A previous structure (Pochapsky et al. (1994) Biochemistry 33, 6424-6432; PDB entry ) was calculated using the results of homonuclear two-dimensional NMR experiments. New data has made it possible to calculate a refinement of the original Pdxo solution structure. First, essentially complete assignments for diamagnetic 15N and 13C resonances of Pdxo have been made using multidimensional NMR methods, and 15N- and 13C-resolved NOESY experiments have permitted the identification of many new NOE restraints for structural calculations. Stereospecific assignments for leucine and valine CH3 resonances were made using biosynthetically directed fractional 13C labeling, improving the precision of NOE restraints involving these residues. Backbone dihedral angle restraints have been obtained using a combination of two-dimensional J-modulated 15N,1H HSQC and 3D (HN)CO(CO)NH experiments. Second, the solution structure of a diamagnetic form of Pdx, that of the C85S variant of gallium putidaredoxin, in which a nonligand Cys is replaced by Ser, has been determined (Pochapsky et al. (1998) J. Biomol. NMR 12, 407-415), providing information concerning structural features not observable in the native ferredoxin due to paramagnetism. Third, a crystal structure of a closely related ferredoxin, bovine adrenodoxin, has been published (Müller et al. (1998) Structure 6, 269-280). This structure has been used to model the metal binding site structure in Pdx. A family of fourteen structures is presented that exhibits an rmsd of 0.51 A for backbone heavy atoms and 0.83 A for all heavy atoms. Exclusion of the modeled metal binding loop region reduces overall the rmsd to 0.30 A for backbone atoms and 0.71 A for all heavy atoms.

Putidaredoxin (Pdx), a Cys4Fe2S2 ferredoxin from Pseudomonas putida, exhibits redox-dependent binding to its physiological redox partner, cytochrome P450(cam) (CYP101), with the reduced form of Pdx (Pdx(r)) binding with greater affinity to oxidized camphor-bound CYP101 than the oxidized form, Pdx(o). It has been previously shown that Pdx(o) is more dynamic than Pdx(r) on all accessible time scales, and it has been proposed that Pdx(r) samples only a fraction of the conformational substates populated by Pdx(o) on a time average. It is postulated that the ensemble subset populated by Pdx(r) is the same subset that binds CYP101, providing a mechanism for coupling the Pdx oxidation state to binding affinity for CYP101. Evidence from a variety of sources, including redox-dependent shifts of 15N and 13C resonances, indicates that the metal cluster binding loop of Pdx is the primary determinant of redox-dependent conformational selection. Patterns of paramagnetic effects suggest that the metal cluster binding loop contracts around the metal cluster upon reduction, possibly due to the strengthening of hydrogen bonds between the sulfur atoms of the metal cluster and the surrounding polypeptide NH and OH groups. Effects of this perturbation are then transmitted mechanically to other affected regions of the protein. A specific mutation has been introduced into the metal binding loop of Pdx, G40N, that slows conformational exchange sufficiently that the ensemble of conformational substates in Pdx(o) are directly observable as severe broadenings or splittings in affected NMR resonances. Many of the residues most affected by the mutation also show significant exchange contributions to 15N T(2) relaxation in wild-type Pdx(o). As predicted, G40N Pdx(r) shows a collapse of many of these multiplets and broadened lines to form much sharper resonances that are essentially identical to those observed in wild-type Pdx(r), indicating that Pdx(r) occupies fewer conformational substates than does Pdx(o). This is the first direct observation of such redox-dependent ensembles at slow exchange on the chemical shift time scale. These results confirm that conformational selection within the Fe2S2 cluster binding loop is the primary source of redox-dependent changes in protein dynamics in Pdx.

The binding of a nitroxide spin-labeled analog of N-acetyllactosamine to galectin-3, a mammalian lectin of 26 kD size, is studied to map the binding sites of this small oligosaccharide on the protein surface. Perturbation of intensities of cross-peaks in the 15 N heteronuclear single quantum coherence (HSQC) spectrum of full-length galectin-3 owing to the bound spin label is used qualitatively to identify protein residues proximate to the binding site for N-acetyllactosamine. A protocol for converting intensity measurements to a more quantitative determination of distances between discrete protein amide protons and the bound spin label is then described. This protocol is discussed as part of a drug design strategy in which subsequent perturbation of chemical shifts of distance mapped amide cross-peaks can be used effectively to screen a library of compounds for other ligands that bind to the target protein at distances suitable for chemical linkage to the primary ligand. This approach is novel in that it bypasses the need for structure determination and resonance assignment of the target protein.Keywords: Nitroxide spin label; galectin-3; drug design; drug screening; distance mapping Knowing the geometric relationship of various small molecules that bind to protein surfaces can be an important point of reference in the design of effective inhibitors of protein function (Shuker et al. 1996;Hadjuk et al. 1997). Here, we illustrate a method for providing intermolecular distance information based on the use of spin-labeled analogs of known protein ligands to perturb cross-peaks in a 15 N-1 H heteronuclear single quantum coherence (HSQC) spectrum in a distance-dependent fashion. Labeling crosspeaks in terms of distance from a primary ligand binding site allows use of those cross-peaks in subsequent chemical shift perturbation screens for secondary ligands that bind at appropriate distances for chemical linkage to the primary ligand.There are precedents in the literature for using a spinlabeled ligand to gain information about protein-ligand interactions (Kosen 1989;Johnson et al. 1999). The most common application employs nitroxide spin labels incorporated in chemically stable molecules such as 2,2,6,6-tetramethyl-1-piperidine-1-oxyl (TEMPO). The route to distance information makes use of the paramagnetic relaxation of NMR resonances caused by these spin-labeled compounds. The relaxation rate enhancement depends on the distance of the nuclei from the unpaired electron, a property that can be used to map distances of nuclei from the spin label. The unpaired electron-nucleus interaction is fairly long-range, unlike the nuclear Overhauser effect (NOE), which is short range in nature (<5 Å). Hence, distances up to 20 Å can be mapped. This general procedure has been used recently to characterize the nature of the interaction between a cellulose-derived ligand and the cellulose binding domains of ␤-1,4 glucanase CenC from Cellulomonas fimi (Johnson et al. 1999). Also, there is some related work published recently t...