The role of backbone motions in ligand binding to the c-Src SH3 domain.

The Src homology 3 (SH3) domain of pp60(c-src) (Src) plays dual roles in signal transduction, through stabilizing the repressed form of the Src kinase and through mediating the formation of activated signaling complexes. Transition of the Src SH3 domain between a variety of binding partners during progression through the cell cycle requires adjustment of a delicate free energy balance. Although numerous structural and functional studies of SH3 have provided an in-depth understanding of structural determinants for binding, the origins of binding energy in SH3-ligand interactions are not fully understood. Considering only the protein-ligand interface, the observed favorable change in standard enthalpy (DeltaH=-9.1 kcal/mol) and unfavorable change in standard entropy (TDeltaS=-2.7 kcal/mol) upon binding the proline-rich ligand RLP2 (RALPPLPRY) are inconsistent with the predominantly hydrophobic interaction surface. To investigate possible origins of ligand binding energy, backbone dynamics of free and RLP2-bound SH3 were performed via (15)N NMR relaxation and hydrogen-deuterium (H/(2)H) exchange measurements. On the ps-ns time scale, assuming uncorrelated motions, ligand binding results in a significant reduction in backbone entropy (-1.5(+/-0.6) kcal/mol). Binding also suppresses motions on the micros-ms time scale, which may additionally contribute to an unfavorable change in entropy. A large increase in protection from H/(2)H exchange is observed upon ligand binding, providing evidence for entropy loss due to motions on longer time scales, and supporting the notion that stabilization of pre-existing conformations within a native state ensemble is a fundamental paradigm for ligand binding. Observed changes in motion on all three time scales occur at locations both near and remote from the protein-ligand interface. The propagation of ligand binding interactions across the SH3 domain has potential consequences in target selection through altering both free energy and geometry in intact Src, and suggests that looking beyond the protein-ligand interface is essential in understanding ligand binding energetics.

An improved method for distinguishing between anisotropic tumbling and chemical exchange in analysis of 15N relaxation parameters.

Although an accurate description of global tumbling of a protein is essential for correct analysis of internal motions. proper distinction between the effects of anisotropic rotational diffusion and conformational exchange has remained a challenge. We present a novel two-part filtering procedure designed specifically to distinguish between the effects of anisotropy and conformational exchange. The efficacy of this method is assessed using synthetic data sets. The method is then applied to two proteins of dramatically different size and shape, OspA and SH3. The large size and extreme anisotropy of OspA provide a challenging case, where conformational exchange is a small perturbation of the effects of anisotropy on transverse relaxation rates. Conversely, in the chicken c-Src SH3 domain, with its small size and nearly spherical shape, anisotropy is a small perturbation of the effects of conformational exchange on transverse relaxation rates. Accurate extraction of the global tumbling parameters for each protein allows optimal characterization of conformational exchange processes, as well as ps-ns time scale motions.

Importance of the N terminus of rous sarcoma virus protease for structure and enzymatic function.

All retrovirus proteases (PRs) are homodimers, and dimerization is essential for enzymatic function. The dimer is held together largely by a short four-stranded antiparallel beta sheet composed of the four or five N-terminal amino acid residues and a similar stretch of residues from the C terminus. We have found that the enzymatic and structural properties of Rous sarcoma virus (RSV) PR are exquisitely sensitive to mutations at the N terminus. Deletion of one or three residues, addition of one residue, or substitution of alanine for the N-terminal leucine reduced enzymatic activity on peptide and protein substrates 100- to 1,000-fold. The purified mutant proteins remained monomeric up to a concentration of about 2 mg/ml, as determined by dynamic light scattering. At higher concentrations, dimerization was observed, but the dimer lacked or was deficient in enzymatic activity and thus was inferred to be structurally distinct from a wild-type dimer. The mutant protein lacking three N-terminal residues (DeltaLAM), a form of PR occurring naturally in virions, was examined by nuclear magnetic resonance spectroscopy and found to be folded at concentrations where it was monomeric. This result stands in contrast to the report that a similarly engineered monomeric PR of human immunodeficiency virus type 1 is unstructured. Heteronuclear single quantum coherence spectra of the mutant at concentrations where either monomers or dimers prevail were nearly identical. However, these spectra differed from that of the dimeric wild-type RSV PR. These results imply that the chemical environment of many of the amide protons differed and thus that the three-dimensional structure of the DeltaLAM PR mutant is different from that of the wild-type PR. The structure of this mutant protein may serve as a model for the structure of the PR domain of the Gag polyprotein and may thus give clues to the initiation of proteolytic maturation in retroviruses.

An increase in side chain entropy facilitates effector binding: NMR characterization of the side chain methyl group dynamics in Cdc42Hs.

Cdc42Hs is a signal transduction protein that is involved in cytoskeletal growth and organization. We describe here the methyl side chain dynamics of three forms of (2)H,(13)C,(15)N-Cdc42Hs [GDP-bound (inactive), GMPPCP-bound (active), and GMPPCP/PBD46-bound (effector-bound)] from (13)C-(1)H NMR measurements of deuterium T(1) and T(1 rho) relaxation times. A wide variation in flexibility was observed throughout the protein, with methyl axis order parameters (S(2)(axis)) ranging from 0.2 to 0.4 (highly disordered) in regions near the PBD46 binding site to 0.8--1.0 (highly ordered) in some helices. The side chain dynamics of the GDP and GMPPCP forms are similar, with methyl groups on the PBD46 binding surface experiencing significantly greater mobility (lower S(2)(axis)) than those not on the binding surface. Binding of PBD46 results in a significant increase in the disorder and a corresponding increase in entropy for the majority of methyl groups. Many of the methyl groups that experience an increase in mobility are found in residues that are not part of the PBD46 binding interface. This entropy gain represents a favorable contribution to the overall entropy of effector binding and partially offsets unfavorable entropy losses such as those that occur in the backbone.

Phosphorylation-induced structural changes in the amyloid precursor protein cytoplasmic tail detected by NMR.

The cytoplasmic tail of the amyloid precursor protein (APPc) interacts with several cellular factors implicated in intracellular signaling or proteolytic production of amyloid beta peptide found in senile plaques of Alzheimer's disease patients. APPc contains two threonine residues (654 and 668 relative to APP695, or 6 and 20 relative to APPc) and a serine residue (655 or 7, respectively) that are known to be phosphorylated in vivo and may play regulatory roles in these events. We show by solution NMR spectroscopy of a 49 residue cytoplasmic tail peptide (APP-C) that in all three cases, phosphorylation induces changes in backbone dihedral angles that can be attributed to formation of local hydrogen bonds between the phosphate group and nearby amide protons. Phosphorylation of S7 also induces chemical shift changes in the hydrophobic cluster (residues I8-V13), indicating additional medium-range effects. The most pronounced changes occur upon phosphorylation of T20, a neuron-specific phosphorylation site, where the N-terminal helix capping box previously characterized for this region is altered. Characterization of torsion angles and transient hydrogen bonds indicates that prolyl isomerization of the pThr-Pro peptide bond results from both destabilization of the N-terminal helix capping box and stabilization of the cis isomer by transient hydrogen bonds. The significant population of the cis isomer (9 %) present after phosphorylation of T20 suggests a potential role of selective recognition of cis versus trans isomers in response to phosphorylation of APP. Together, these structural changes indicate that phosphorylation may act as a conformational switch in the cytoplasmic tail of APP to alter specificity and affinity of binding to cytosolic partners, particularly in response to the abnormal phosphorylation events associated with Alzheimer's disease.

Solution structure of ThiS and implications for the evolutionary roots of ubiquitin.

ThiS is a sulfur carrier protein that plays a central role in thiamin biosynthesis in Escherichia coli. Here we report the solution NMR structure of ThiS, the first for this class of sulfur carrier proteins. Although ThiS shares only 14% sequence identity with ubiquitin, it possesses the ubiquitin fold. This structural homology, combined with established functional similarities involving sulfur chemistry, demonstrates that the eukaryotic ubiquitin and the prokaryotic ThiS evolved from a common ancestor. This illustrates how structure determination is essential in establishing evolutionary links between proteins in which structure and function have been conserved through eons of evolution despite loss of sequence identity. The ThiS structure reveals both hydrophobic and electrostatic surface features that are likely determinants for interactions with binding partners. Comparison with surface features of ubiquitin and ubiquitin homologs SUMO-1, RUB-1 and NEDD8 suggest how Nature has utilized this single fold to incorporate similar chemistry into a broad array of highly specific biological processes.

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR.

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.

Protein dynamics measurements by TROSY-based NMR experiments.

The described TROSY-based experiments for investigating backbone dynamics of proteins make it possible to elucidate internal motions in large proteins via measurements of T(1), T(2), and NOE of backbone (15)N nuclei. In our proposed sequences, the INEPT sequence is eliminated and the PEP sequence is replaced by the ST2-PT sequence from the HSQC-based experiments. This has the benefit of shortening the pulse sequences by 5.4 ms (=1/2J) and results in an increase in the intrinsic sensitivity of the proposed TROSY-based experiments. The TROSY-based experiments are on average of 13% more sensitive than the corresponding HSQC-based experiments on a uniformly (15)N-labeled Xenopus laevis calcium-bound calmodulin sample on a 750-MHz spectrometer at 5 degrees C. The amide proton linewidths of the TROSY-based experiments are 2-13 Hz narrower than those of the HSQC experiments. More sensitivity gain and higher resolution are expected if the protein sample is deuterated.

Transient structure of the amyloid precursor protein cytoplasmic tail indicates preordering of structure for binding to cytosolic factors.

The cytoplasmic tail of the amyloid precursor protein (APP) appears to play two important roles in the cell through participation in intracellular signaling and proteolytic processing of APP. Hence, knowledge of the structure of the 47 residue cytoplasmic tail of APP is important for understanding the molecular interactions involved in normal cell function as well as in the pathogenesis of Alzheimer's disease. Multidimensional solution NMR spectroscopy has been applied to examine the structural features of a 49-residue peptide (APP-C) containing two N-terminal residues (GS) and the APP cytoplasmic tail, over the pH range of 4.2-7.1. Although the peptide does not adopt a stable folded structure, regions of unstable structure exist over the pH range examined and have been characterized by a combination of H(alpha) chemical shifts, NOE analysis, and (3)J(HNH)(alpha) coupling constants and by identification of transient hydrogen bonds between amide protons and titrating carboxylate groups. These studies extend the work of others [Kroenke et al. (1997) Biochemistry 36, 8145-8152] by identifying an additional nascent helix and a hydrophobic cluster within the N-terminal 20 amino acid residues and by further characterizing the TPEE turn as a helix capping box. The transient structure of APP-C provides insight into the importance of preordering of this cytoplasmic tail in governing specificity and affinity for cytosolic binding partners.

Backbone dynamics of inactive, active, and effector-bound Cdc42Hs from measurements of (15)N relaxation parameters at multiple field strengths.

Cdc42Hs, a member of the Ras superfamily of GTP-binding proteins, initiates a cascade that begins with the activation of several kinases, including p21-activated kinase (PAK). We have previously determined the structure of Cdc42Hs and found that the regions involved in effector (Switch I) and regulator (Switch II) actions are partially disordered [Feltham, J. L., et al. (1997) Biochemistry 36, 8755-8766]. Recently, we used a 46-amino acid fragment of PAK (PBD46) to define the binding surface on Cdc42Hs, which includes the beta2 strand and a portion of Switch I [Guo, W., et al. (1998) Biochemistry 37, 14030-14037]. Here we describe the backbone dynamics of three constructs of [(15)N]Cdc42Hs (GDP-, GMPPCP-, and GMPPCP- and PBD46-bound) using (15)N-(1)H NMR measurements of T(1), T(1)(rho), and the steady-state NOE at three magnetic field strengths. Residue-specific values of the generalized order parameters (S(s)(2) and S(f)(2)), local correlation time (tau(e)), and exchange rate (R(ex)) were obtained using the Lipari-Szabo model-free formalism. Residues in Switch I were found to exhibit high-amplitude (low-order) motions on a nanosecond time scale, whereas those in Switch II experience low-amplitude motion on the nanosecond time scale and chemical (conformational) exchange on a millisecond time scale. The Insert region of Cdc42Hs-GDP exhibits high-order, nanosecond motions; the time scale of motion in the Insert is reduced in Cdc42Hs-GMPPCP and Cdc42Hs-PBD46. Overall, significant flexibility was observed mainly in the regions of Cdc42Hs that are involved in protein-protein interactions (Switch I, Switch II, and Insert), and flexibility was reduced upon interaction with a protein ligand. These results suggest that protein flexibility is important for high-affinity binding interactions.

Heteronuclear NMR studies of the combined Src homology domains 2 and 3 of pp60 c-Src: effects of phosphopeptide binding.

The results of heteronuclear NMR studies on the combined Src homology domains 2 and 3 (SH3-SH2) of pp60 c-Src are presented. Resonance assignments were obtained using heteronuclear triple-resonance experiments in conjunction with 15N-separated nuclear Overhauser effect spectroscopy (NOESY) data. A modified three-dimensional 13CO-15N-1H spectral correlation experiment [(HACA)CO(CA)-NH] with improved sensitivity is presented that provided additional sequential information and resolved several ambiguities. Chemical shifts and sequential- and medium-range NOE cross peaks indicate that the structures of both the SH3 and SH2 portions of the polypeptide are very similar to those of the isolated SH3 and SH2 domains. Binding of a high-affinity phosphopeptide, EPQpYEEIPIYL, induces large chemical shift changes at several locations in the SH2 domain. Comparison with known results for peptide binding to SH2 domains shows that the residues displaying the largest effects are all involved in peptide binding or undergo significant conformational changes upon binding. However, subtle changes of both 1H and 15N chemical shifts are observed for residues within the SH3 domain and the connecting linker region, indicating possible cross-domain communication.

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.

Flexibility and function in HIV-1 protease.

HIV protease is a homodimeric protein whose activity is essential to viral function. We have investigated the molecular dynamics of the HIV protease, thought to be important for proteinase function, bound to high affinity inhibitors using NMR techniques. Analysis of 15N spin relaxation parameters, of all but 13 backbone amide sites, reveals the presence of significant internal motions of the protein backbone. In particular, the flaps that cover the proteins active site of the protein have terminal loops that undergo large amplitude motions on the ps to ns time scale, while the tips of the flaps undergo a conformational exchange on the microsecond time scale. This enforces the idea that the flaps of the proteinase are flexible structures that facilitate function by permitting substrate access to and product release from the active site of the enzyme.

Structural analysis of highly oriented poly(p-phenylene-terephthalamide) by 15N solid-state nuclear magnetic resonance.

The structure of an aromatic polyamide, poly(p-phenylene-terephthalamide) (PPTA), was studied in the solid state using 15N nuclear magnetic resonance (NMR) spectroscopy. Spectra of uniaxially aligned molecules placed with the axis of alignment both parallel with and perpendicular to the applied magnetic field were analyzed to yield the orientations of specific molecular bonds with respect to the fiber axis. The 15N chemical shift tensor was characterized by simulating powder pattern spectra of both PPTA and a model compound, benzanilide. Chemical shift and dipolar coupled chemical shift line shapes were calculated through Euler angle transformations from the principal axis system (PAS) reference frame to the fiber axis system (FAS) frame. The orientations of NH and NC' bonds in PPTA are determined as well as the orientational distribution of the PPTA fiber axis. The structural parameters determined for PPTA are compared with those obtained by X-ray diffraction.

Secondary structure and signal assignments of human-immunodeficiency-virus-1 protease complexed to a novel, structure-based inhibitor.

We report comprehensive NMR studies in solution of the human-immunodeficiency-virus (HIV)-1 protease. Stable solutions of the protease were obtained by complexing the protein to a designed cyclic urea inhibitor DMP 323. A variety of triple-resonance experiments provided essentially complete 1H, 13C and 15N NMR signal assignments of the protease. These assignments, together with short-range NOE constraints, coupling constants and hydrogen-exchange data, yielded the secondary structure of the protease in solution. The results reported herein open the way to the determination of the high-resolution three-dimensional solution structures of protease/inhibitor complexes, as well as to studies of protease dynamics and solvent interactions.

A method for studying the structure of uniaxially aligned biopolymers using solid state 15N-nmr: application to Bombyx mori silk fibroin fibers.

Recent advances in the application of solid state nmr spectroscopy to uniformly aligned biopolymers have opened a window through which to view the detailed structure of biological macromolecules that are unable to be seen with standard techniques for structure determination such as x-ray diffraction. Atomic resolution structural details are obtained from solid state nmr data in the form of bond orientations, which yield the relative positions of specific atoms within the molecule. For static aligned systems such as fibers, in which rapid reorientation about the axis of alignment does not occur, it has generally been necessary to perform trial and error line-shape simulations to extract structural details from nmr spectra arising from a single type of nuclear spin interaction. In the present work, a new method is developed in which solid state 15N-nmr spectra obtained from uniaxially aligned molecules placed with the axis of alignment both parallel and perpendicular to the magnetic field are analyzed to yield the orientations of specific molecular bonds. Analytical expressions are derived that utilize spectral features read from 15N chemical shift anisotropy line shapes to calculate a discrete number of possible orientations for a specific site. The 15N-1H dipolar interaction is employed to further narrow the number of unique orientations possible for a given site. With this method, a neighborhood of possible orientations is quickly determined, and full line-shape simulations within this region of allowed space can be performed to refine the limits of orientation. This technique demonstrates the use of a single type of isotopic label to determine the orientation of a specific molecular group such as a peptide plane within a protein. Results from the application of this method to the Bombyx mori silk fibroin protein provide structural detail that is consistent with currently accepted structural models based on fiber diffraction studies.

Dynamics of methyl groups in proteins as studied by proton-detected 13C NMR spectroscopy. Application to the leucine residues of staphylococcal nuclease.

This paper describes the application of recently developed nuclear magnetic resonance (NMR) pulse sequences to obtain information about the internal dynamics of isotopically enriched hydrophobic side chains in proteins. The two-dimensional spectra provided by the pulse sequences enable one to make accurate measurements of nuclear Overhauser effects (NOE) and longitudinal (T1) and transverse (T2) relaxation times of enriched methyl carbons in proteins. Herein, these techniques are used to investigate the internal dynamics of the 11 leucine side chains of staphylococcal nuclease (SNase), a small enzyme having Mr = 16.8K, in the absence and presence of ligands thymidine 3',5'-bisphosphate (pdTp) and Ca2+. We report the synthesis of [5,5'-13C2]leucine, the preparation of SNase containing the labeled leucine, the sequential assignment of the leucine methyl carbons and protons in the liganded and unliganded proteins, and the measurement of the 13C T1, T2, and NOE values for the SNase leucine methyl carbons. Analysis of the relaxation parameters using the formalism of Lipari and Szabo shows that the internal motions of the leucine methyl carbons are characterized by effective correlation times tau f (5-80 ps) and tau s (less than 2 ns). The fast motion is identified with the rapid rotation of the methyl group about the C gamma-C delta bond axis, while the slow motion is associated with reorientation of the C gamma-C delta bond axis itself. The mean squared order parameters associated with the latter motion, Ss2, lie in the range 0.34-0.92. The values of Ss2 correlate reasonably well with the temperature factors of the leucine methyl carbons obtained from the crystal structures, but some are smaller than anticipated on the basis of the fact that nearly all leucine methyl carbons are buried and have temperature factors no larger than that of the leucine backbone atoms. Five leucine residues in liganded SNase and eight in unliganded SNase have values of Ss2 less than 0.71. These order parameters correspond to large amplitude motions (angular excursions of 27-67 degrees) of the C gamma-C delta bond axis. These results indicate that, in solution, the internal motions of the leucine side chains of SNase are significantly larger than suggested by the X-ray structures or by qualitative analysis of NOESY spectra. Comparison of Ss2 values obtained from liganded and unliganded SNase reveals a strong correlation between delta Ss2 and distance between the leucine methyl carbon and the ligands.(ABSTRACT TRUNCATED AT 400 WORDS)

Molecular dynamics computations and solid state nuclear magnetic resonance of the gramicidin cation channel.

This paper reports on a coupled approach to determining the structure of the gramicidin A ion channel, utilizing solid state nuclear magnetic resonance (NMR) of isotopically labeled gramicidin channels aligned parallel to the magnetic field direction, and molecular dynamics (MD). MD computations using an idealized right-handed beta-helix as a starting point produce a refined molecular structure that is in excellent agreement with atomic resolution solid state NMR data. The data provided by NMR and MD are complementary to each other. When applied in a coordinated manner they provide a powerful approach to structure determination in molecular systems not readily amenable to x-ray diffraction.