Backbone dynamics of free barnase and its complex with barstar determined by 15N NMR relaxation study

Backbone dynamics of uniformly ¹⁵N-labeled free barnase and its complex with unlabelled barstar have been studied at 40°C, pH 6.6, using ¹⁵N relaxation data obtained from proton-detected 2D {¹H}-¹⁵N NMR spectroscopy. ¹⁵N spin-lattice relaxation rate constants (R₁), spin-spin relaxation rate constants (R₂), and steady-state heteronuclear {¹H}-¹⁵N NOEs have been measured at a magnetic field strength of 14.1 Tesla for 91 residues of free barnase and for 90 residues out of a total of 106 in the complex (excluding three prolines and the N-terminal residue) backbone amide ¹⁵N sites of barnase. The primary relaxation data for both the cases have been analyzed in the framework of the model-free formalism using both isotropic and axially symmetric models of the rotational diffusion tensor. As per the latter, the overall rotational correlation times (_m) are 5.0 and 9.5 ns for the free and complexed barnase, respectively. The average order parameter is found to be 0.80 for free barnase and 0.86 for the complex. However, the changes are not uniform along the backbone and for about 5 residues near the binding interface there is actually a significant decrease in the order parameters on complex formation. These residues are not involved in the actual binding. For the residues where the order parameter increases, the magnitudes vary significantly. It is observed that the complex has much less internal mobility, compared to free barnase. From the changes in the order parameters, the entropic contribution of NH bond vector motion to the free energy of complex formation has been calculated. It is apparent that these motions cause significant unfavorable contributions and therefore must be compensated by many other favorable contributions to effect tight complex formation. The observed variations in the motion and their different locations with regard to the binding interface may have important implications for remote effects and regulation of the enzyme action.     

Backbone dynamics of the Bacillus subtilis glucose permease IIA domain determined from 15N NMR relaxation measurements.

The backbone dynamics of the uniformly 15N-labeled IIA domain of the glucose permease of Bacillus subtilis have been characterized using inverse-detected two-dimensional 1H-15N NMR spectroscopy. Longitudinal (T1) and transverse (T2) 15N relaxation time constants and steady-state (1H)-15N NOEs were measured, at a spectrometer proton frequency of 500 MHz, for 137 (91%) of the 151 protonated backbone nitrogens. These data were analyzed by using a model-free dynamics formalism to determine the generalized order parameter (S2), the effective correlation time for internal motions (tau e), and 15N exchange broadening contributions (Rex) for each residue, as well as the overall molecular rotational correlation time (tau m). The T1 and T2 values for most residues were in the ranges 0.45-0.55 and 0.11-0.15 s, respectively; however, a small number of residues exhibited significantly slower relaxation. Similarly, (1H)-15N NOE values for most residues were in the range 0.72-0.80, but a few residues had much smaller positive NOEs and some exhibited negative NOEs. The molecular rotational correlation time was 6.24 +/- 0.01 ns; most residues had order parameters in the range 0.75-0.90 and tau e values of less than ca. 25 ps. Residues found to be more mobile than the average were concentrated in three areas: the N-terminal residues (1-13), which were observed to be highly disordered; the loop from P25 to D41, the apex of which is situated adjacent to the active site and may have a role in binding to other proteins; and the region from A146 to S149. All mobile residues occurred in regions close to termini, in loops, or in irregular secondary structure.     

Backbone dynamics of a bacterially expressed peptide from the receptor binding domain of Pseudomonas aeruginosa pilin strain PAK from heteronuclear 1H-15N NMR spectroscopy

The backbone dynamics of a ¹⁵N-labeled recombinant PAK pilin peptide spanning residues 128–144 in the C-terminal receptor binding domain of Pseudomonas aeruginosa pilin protein strain PAK (Lys¹²⁸-Cys-Thr-Ser-Asp-Gln-Asp-Glu-Gln-Phe-Ile-Pro-Lys-Gly-Cys-Ser-Lys¹⁴⁴) were probed by measurements of ¹⁵N NMR relaxation. This PAK(128–144) sequence is a target for the design of a synthetic peptide vaccine effective against multiple strains of P. aeruginosa infection. The ¹⁵N longitudinal (T₁) and transverse (T₂) relaxation rates and the steady-state heteronuclear {¹H}-¹⁵N NOE were measured at three fields (7.04, 11.74 and 14.1 Tesla), five temperatures (5, 10, 15, 20, and 25°C ) and at pH 4.5 and 7.2. Relaxation data was analyzed using both the `model-free' formalism [Lipari, G. and Szabo, A. (1982) J. Am. Chem. Soc., 104, 4546–4559 and 4559–4570] and the reduced spectral density mapping approach [Farrow, N.A., Szabo, A., Torchia, D.A. and Kay, L.E. (1995) J. Biomol. NMR, 6, 153–162]. The relaxation data, spectral densities and order parameters suggest that the type I and type II -turns spanning residues Asp¹³⁴-Glu-Gln-Phe¹³⁷ and Pro¹³⁹-Lys-Gly-Cys¹⁴², respectively, are the most ordered and structured regions of the peptide. The biological implications of these results will be discussed in relation to the role that backbone motions play in PAK pilin peptide immunogenicity, and within the framework of developing a pilin peptide vaccine capable of conferring broad immunity across P. aeruginosa strains.     

A novel view of domain flexibility in E. coli adenylate kinase based on structural mode-coupling (15)N NMR relaxation.

Adenylate kinase from Escherichia coli (AKeco), consisting of a single 23.6 kDa polypeptide chain folded into domains CORE, AMPbd and LID, catalyzes the reaction AMP+ATP-->2ADP. In the ligand-free enzyme the domains AMPbd and LID execute large-amplitude movements controlling substrate binding and product release during catalysis. Domain flexibility is investigated herein with the slowly relaxing local structure (SRLS) model for (15)N relaxation. SRLS accounts rigorously for coupling between the global and local N-H motions through a local ordering potential exerted by the protein structure at the N-H bond. The latter reorients with respect to its protein surroundings, which reorient on the slower time scale associated with the global protein tumbling. AKeco diffuses globally with correlation time tau(m)=15.1 ns, while locally two different dynamic cases prevail. The domain CORE features ordering about the equilibrium N-H bond orientation with order parameters, S(2), of 0.8-0.9 and local motional correlation times, tau, mainly between 5-130 ps. This represents a conventional rigid protein structure with rapid small-amplitude N-H fluctuations. The domains AMPbd and LID feature small parallel (Z(M)) ordering of S(2)=0.2-0.5 which can be reinterpreted as high perpendicular (Y(M)) ordering. M denotes the local ordering/local diffusion frame. Local motion about Z(M) is given by tau( parallel) approximately 5 ps and local motion of the effective Z(M) axis about Y(M) by tau( perpendicular)=6-11 ns. Z(M) is tilted at approximately 20 degrees from the N-H bond. The orientation of the Y(M) axis may be considered parallel to the C(alpha)(i-1)-C(alpha)(i) axis. The tau( perpendicular) mode reflects collective nanosecond peptide-plane motions, interpretable as domain motion. A powerful new model of protein flexibility/domain motion has been established. Conformational exchange (R(ex)) processes accompany the tau( perpendicular) mode. The SRLS analysis is compared with the conventional model-free analysis.     

Backbone dynamics of the natively unfolded pro-peptide of subtilisin by heteronuclear NMR relaxation studies

The dynamics of the natively unfolded form of the pro-peptide of subtilisin (PPS) have been characterized at two different pHs (6.0 and 3.0) by ¹⁵N relaxation experiments. ¹⁵N relaxation data is obtained at multiple field strengths and a detailed comparison of spectral density mapping, the model free approach and the recently proposed Cole–Cole model free (CC-MF) analysis is presented. The CC-MF analysis provides a better fit to the observed magnetic field dependence of ¹⁵N relaxation data of unfolded PPS than conventional model free approaches and shows that fluctuations in R₂ may be accounted for by a distribution of correlation times on the nanosecond timescale. A new parameter derives from the analysis and represents the width of the distribution function and the heterogeneity of the dynamics on the nanosecond timescale at a particular site. Particularly interesting is the observation that is sensitive to pH changes and that PPS samples a wider distribution of nanosecond time scale motions at less acidic pHs than at more acidic pHs. These results suggest that PPS experiences a higher degree of correlated motion at pH 6.0 and that electrostatic interactions may be important for inducing correlated motions on the nanosecond timescale in unfolded PPS.     

Backbone dynamics of the cytotoxic ribonuclease alpha-sarcin by 15N NMR relaxation methods

The cytotoxic ribonuclease -sarcin is a 150-residue protein that inactivates ribosomes by selectively cleaving a single phosphodiester bond in a strictly conserved rRNA loop. In order to gain insights on the molecular basis of its highly specific activity, we have previously determined its solution structure and studied its electrostatics properties. Here, we complement those studies by analysing the backbone dynamics of -sarcin through measurement of longitudinal relaxation rates R₁, off resonance rotating frame relaxation rates R₁, and the ¹⁵N¹HNOE of the backbone amide ¹⁵N nuclei at two different magnetic field strengths (11.7 and 17.6 T). The two sets of relaxation parameters have been analysed in terms of the reduced spectral density mapping formalism, as well as by the model-free approach. -Sarcin behaves as an axial symmetric rotor of the prolate type (D/D=1.16 ± 0.02) which tumbles with a correlation time _m of 7.54 ± 0.02 ns. The rotational diffusion properties have been also independently evaluated by hydrodynamic calculations and are in good agreement with the experimental results. The analysis of the internal dynamics reveals that -sarcin is composed of a rigid hydrophobic core and some exposed segments which undergo fast (ps to ns) internal motions. Slower motions in the s to ms time scale are less abundant and in some cases can be assigned to specific motional processes. All dynamic data are discussed in relation to the role of some particular residues of -sarcin in the process of recognition of its ribosomal target.     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.     

Structural and dynamic characterization of the urea denatured state of the immunoglobulin binding domain of streptococcal protein G by multidimensional heteronuclear NMR spectroscopy.

The structure and dynamics of the urea-denatured B1 immunoglobulin binding domain of streptococcal protein G (GB1) has been investigated by multidimensional heteronuclear NMR spectroscopy. Complete 1H, 15N, and 13C assignments are obtained by means of sequential through-bond correlations. The nuclear Overhauser enhancement, chemical shift, and 3JHN alpha coupling constant data provide no evidence for the existence of any significant population of residual native or nonnative ordered structure. 15N relaxation measurements at 500 and 600 MHz, however, provide evidence for conformationally restricted motions in three regions of the polypeptide that correspond to the second beta-hairpin, the N-terminus of the alpha-helix, and the middle of the alpha-helix in the native protein. The time scale of these motions is longer than the apparent overall correlation time (approximately 3 ns) and could range from about 6 ns in the case of one model to between 4 microseconds and 2 ms in another; it is not possible to distinguish between these two cases with certainty because the dynamics are highly complex and hence the analysis of the time scale of this slower motion is highly model dependent. It is suggested that these three regions may correspond to nucleation sites for the folding of the GB1 domain. With the exception of the N- and C-termini, where end effects predominate, the amplitude of the subnanosecond motions, on the other hand, are fairly uniform and model independent, with an overall order parameter S2 ranging from 0.4 to 0.5.     

Redox-controlled backbone dynamics of human cytochrome c revealed by N NMR relaxation measurements

Redox-controlled backbone dynamics in cytochrome c (Cyt c) were revealed by 2D ¹⁵N NMR relaxation experiments. ¹⁵N T₁ and T₂ values and ¹H-¹⁵N NOEs of uniformly ¹⁵N-labeled reduced and oxidized Cyt c were measured, and the generalized order parameters (S²), the effective correlation time for internal motion (τ_e), the ¹⁵N exchange broadening contributions (R_ex) for each residue, and the overall correlation time (τ_m) were estimated by model-free dynamics formalism. These dynamic parameters clearly showed that the backbone dynamics of Cyt c are highly restricted due to the covalently bound heme that functions as the stable hydrophobic core. Upon oxidation of the heme iron in Cyt c, the average S² value was increased from 0.88 ± 0.01 to 0.92 ± 0.01, demonstrating that the mobility of the backbone is further restricted in the oxidized form. Such increases in the S² values were more prominent in the loop regions, including amino acid residues near the thioether bonds to the heme moiety and positively charged region around Lys87. Both of the regions are supposed to form the interaction site for cytochrome c oxidase (CcO) and the electron pathway from Cyt c to CcO. The redox-dependent mobility of the backbone in the interaction site for the electron transfer to CcO suggests an electron transfer mechanism regulated by the backbone dynamics in the Cyt c-CcO system.     

Off-resonance effects in 15N T2 CPMG measurements

The systematic difference between T ₂ values obtained from CPMG and T ₁ experiments was observed for backbone ¹⁵N nuclei of bacterial ribonuclease barnase. Theoretical consideration suggests that the observed difference is caused by off-resonance effects of 180° pulses of the CPMG pulse train. Namely, at off-resonance conditions T ₁-dependent secondary echo coherence pathways considerably contribute to the signal decay in the CPMG experiment and result in systematic (up to 10%) offset-dependent overestimation of ¹⁵N T ₂ measured by the CPMG technique. Under certain circumstances off-resonance effects result in dependence of ¹⁵N T ₂ on CPMG frequency, which might be erroneously interpreted as conformational exchange on the millisecond time-scale. A procedure for numerical correction of ¹⁵N T ₂ (CPMG) data is proposed.     

Backbone dynamics of the regulatory domain of calcium vector protein, studied by 15N relaxation at four fields, reveals unique mobility characteristics of the intermotif linker

CaVP is a calcium-binding protein from amphioxus. It has a modular composition with two domains, but only the two EF-hand motifs localized in the C-terminal domain are functional. We recently determined the solution structure of this regulatory half (C-CaVP) in the Ca(2+)-saturated form and characterized the stepwise ion binding. This paper reports the (15)N nuclear relaxation rates of the Ca(2+)-saturated C-CaVP, measured at four different NMR fields (9.39, 11.74, 14.1, and 18.7 T), which were used to map the spectral density function for the majority of the amide H(N)-N vectors. Fitting the spectral density values at eight frequencies by a model-free approach, we obtained the microdynamic parameters characterizing the global and internal movements of the polypeptide backbone. The two EF-hand motifs, including the ion binding loops, behave like compact structural units with restricted mobility as reflected in the quite uniform order parameter and short internal correlation time (< 20 nsec). Comparative analysis of the two Ca(2+) binding sites shows that site III, having a larger affinity for the metal ion, is generally more rigid, and the amide vector in the second residue of each loop is significantly less restricted. The linker fragment is animated simultaneously by a larger amplitude fast motion and a slow conformational exchange on a microsecond to millisecond time scale. The backbone dynamics of C-CaVP characterized here is discussed in relation with other well-characterized Ca(2+)-binding proteins. Supplemental material: See www.proteinscience.org     

A suite of Mathematica notebooks for the analysis of protein main chain 15N NMR relaxation data

A suite of Mathematica notebooks has been designed to ease the analysis of protein main chain ¹⁵N NMR relaxation data collected at a single magnetic field strength. Individual notebooks were developed to perform the following tasks: nonlinear fitting of ¹⁵N-T ₁ and -T ₂ relaxation decays to a two parameter exponential decay, calculation of the principal components of the inertia tensor from protein structural coordinates, nonlinear optimization of the principal components and orientation of the axially symmetric rotational diffusion tensor, model-free analysis of ¹⁵N-T ₁, -T ₂, and {¹H}–¹⁵N NOE data, and reduced spectral density analysis of the relaxation data. The principle features of the notebooks include use of a minimal number of input files, integrated notebook data management, ease of use, cross-platform compatibility, automatic visualization of results and generation of high-quality graphics, and output of analyses in text format.L. Spyracopoulos is an AHFMR Medical Research Senior Scholar     

Backbone dynamics of the calcium-signaling protein apo-S100B as determined by 15N NMR relaxation

Backbone dynamics of homodimeric apo-S100B were studied by (15)N nuclear magnetic resonance relaxation at 9.4 and 14.1 T. Longitudinal relaxation (T(1)), transverse relaxation (T(2)), and the (15)N-[(1)H] NOE were measured for 80 of 91 backbone amide groups. Internal motional parameters were determined from the relaxation data using the model-free formalism while accounting for diffusion anisotropy. Rotational diffusion of the symmetric homodimer has moderate but statistically significant prolate axial anisotropy (D( parallel)/D( perpendicular) = 1.15 +/- 0.02), a global correlation time of tau(m) = 7.80 +/- 0.03 ns, and a unique axis in the plane normal to the molecular symmetry axis. Of 29 residues at the dimer interface (helices 1 and 4), only one has measurable internal motion (Q71), and the order parameters of the remaining 28 were the highest in the protein (S(2) = 0.80 to 0.91). Order parameters in the typical EF hand calcium-binding loop (S(2) = 0.73 to 0.87) were slightly lower than in the pseudo-EF hand (S(2) = 0.75 to 0.89), and effective internal correlation times, tau(e), distinct from global tumbling, were detected in the calcium-binding loops. Helix 3, which undergoes a large, calcium-induced conformational change necessary for target-protein binding, does not show evidence of interchanging between the apo and Ca(2+)-bound orientations in the absence of calcium but has rapid motion in several residues throughout the helix (S(2) = 0.78 to 0.88; 10 < or = tau(e) < or = 30 ps). The lowest order parameters were found in the C-terminal tail (S(2) = 0.62 to 0.83). Large values for chemical exchange also occur in this loop and in regions nearby in space to the highly mobile C-terminal loop, consistent with exchange broadening effects observed.     

A comparison of 15N NMR relaxation measurements with a molecular dynamics simulation: Backbone dynamics of the glucocorticoid receptor DNA-binding domain

The rapid motions of the backbone of the DNA-binding domain of the glucocorticoid receptor (GR DBD) have been investigated using proton-detected heteronuclear NMR experiments on ¹⁵N-labeled protein at pH 6.0 and with a 200 psec molecular dynamics simulation of hydrated GR DBD. The experimental data were interpreted in terms of a generalized order parameter (S²) and an effective correlation time (τ_e) for the internal motion of each amide bond. A back calculation, using the same model, yielded the {¹H}-¹⁵N nuclear Overhauser effects (NOEs) and the ¹⁵N spin-lattice relaxation times (T₁) from the simulated data. The rapid motions of the backbone turned out to be rather limited and uniform throughout the protein, with a somewhat reduced mobility in the two major α-helical regions and a slightly enhanced flexibility for some residues in the first zinc coordinating region. The agreement between the experimental and simulated S²-values was as good as quantitative for most of the residues, except for some residues that were subject to a more large-scale, and in the simulation thus poorly sampled, motion. Examples of such motions that were found in the simulation include jumps of the amide bond of Ile-487 between the charged oxygens of the side chain of Asp-485 and less distinct large scale motions for some of the residues in the extended regions, that were shown to give rise to noisy and/or fast decaying internal reorientational correlation functions. For these residues large differences in the simulated and experimental τ_e-values were found, indicating that motions on different time scales were dominating in the experimental and simulated values. The lower (<0.7) experimental NOEs for these residues could not be reproduced in the simulation and were shown to be a consequence of the lower τ_e-values estimated in the simulation. By combining information from the simulation and the experiment a more complete picture of the motions for these residues can be obtained as is illustrated with an estimation of the jump angle and jump frequency for the amide bond of Ile-487. © 1993 Wiley-Liss, Inc.     

Backbone dynamics of oxidized and reduced D. vulgaris flavodoxin in solution

Recombinant Desulfovibrio vulgaris flavodoxin was produced inEscherichia coli. A complete backbone NMR assignment for the two-electronreduced protein revealed significant changes of chemical shift valuescompared to the oxidized protein, in particular for the flavinemononucleotide (FMN)-binding site. A comparison of homo- and heteronuclearNOESY spectra for the two redox states led to the assumption that reductionis not accompanied by significant changes of the global fold of the protein.The backbone dynamics of both the oxidized and reduced forms of D. vulgarisflavodoxin were investigated using two-dimensional¹⁵N-¹H correlation NMR spectroscopy.T₁, T₂ and NOE data are obtained for 95%of the backbone amide groups in both redox states. These values wereanalysed in terms of the model-free approach introduced by Lipari andSzabo [(1982) J. Am. Chem. Soc., 104, 4546-;4559, 4559-;4570]. Acomparison of the two redox states indicates that in the reduced speciessignificantly more flexibility occurs in the two loop regions enclosing FMN.Also, a higher amplitude of local motion could be found for the N(3)H groupof FMN bound to the reduced protein compared to the oxidized state.     

Backbone dynamics measurements on leukemia inhibitory factor, a rigid four-helical bundle cytokine

The backbone dynamics of the four-helical bundle cytokine leukemia inhibitory factor (LIF) have been investigated using 15N NMR relaxation and amide proton exchange measurements on a murine-human chimera, MH35-LIF. For rapid backbone motions (on a time scale of 10 ps to 100 ns), as probed by 15N relaxation measurements, the dynamics parameters were calculated using the model-free formalism incorporating the model selection approach. The principal components of the inertia tensor of MH35-LIF, as calculated from its NMR structure, were 1:0.98:0.38. The global rotational motion of the molecule was, therefore, assumed to be axially symmetric in the analysis of its relaxation data. This yielded a diffusion anisotropy D(parallel)/D(perpendicular) of 1.31 and an effective correlation time (4D(perpendicular) + 2D(parallel))(-1) of 8.9 ns. The average values of the order parameters (S2) for the four helices, the long interhelical loops, and the N-terminus were 0.91, 0.84, and 0.65, respectively, indicating that LIF is fairly rigid in solution, except at the N-terminus. The S2 values for the long interhelical loops of MH35-LIF were higher than those of their counterparts in short-chain members of the four-helical bundle cytokine family. Residues involved in LIF receptor binding showed no consistent pattern of backbone mobilities, with S2 values ranging from 0.71 to 0.95, but residues contributing to receptor binding site III had relatively lower S2 values, implying higher amplitude motions than for the backbone of sites I and II. In the relatively slow motion regime, backbone amide exchange measurements showed that a number of amides from the helical bundle exchanged extremely slowly, persisting for several months in 2H2O at 37 degrees C. Evidence for local unfolding was considered, and correlations among various structure-related parameters and the backbone amide exchange rates were examined. Both sets of data concur in showing that LIF is one of the most rigid four-helical bundle cytokines.     

Backbone dynamics of the human CC-chemokine eotaxin

Eotaxin is a CC chemokine with potent chemoattractant activity towards eosinophils. ¹⁵N NMR relaxation data have been used to characterize the backbone dynamics of recombinant human eotaxin. ¹⁵N longitudinal (R₁) and transverse (R₂) auto relaxation rates, heteronuclear ¹H-¹⁵N steady-state NOEs, and transverse cross-relaxation rates (_xy) were obtained at 30 °C for all resolved backbone secondary amide groups using ¹ H-detected two-dimensional NMR experiments. Ratios of transverse auto and cross relaxation rates were used to identify NH groups influenced by slow conformational rearrangement. Relaxation data were fit to the extended model free dynamics formalism, yielding parameters describing axially symmetric molecular rotational diffusion and the internal dynamics of each NH group. The molecular rotational correlation time (_m) is 5.09±0.02 ns, indicating that eotaxin exists predominantly as a monomer under the conditions of the NMR study. The ratio of diffusion rates about unique and perpendicular axes (D/D) is 0.81±0.02. Residues with large amplitudes of subnanosecond motion are clustered in the N-terminal region (residues 1–19), the C-terminus (residues 68–73) and the loop connecting the first two -strands (residues 30–37). N-terminal flexibility appears to be conserved throughout the chemokine family and may have implications for the mechanism of chemokine receptor activation. Residues exhibiting significant dynamics on the microsecond–millisecond time scale are located close to the two conserved disulfide bonds, suggesting that these motions may be coupled to disulfide bond isomerization.