Characterization of threonine side chain dynamics in an antifreeze protein using natural abundance 13C NMR spectroscopy

The dynamics of threonine side chains of the Tenebrio molitor antifreeze protein (TmAFP) were investigated using natural abundance (13)C NMR. In TmAFP, the array of threonine residues on one face of the protein is responsible for conferring its ability to bind crystalline ice and inhibit its growth. Heteronuclear longitudinal and transverse relaxation rates and the [(1)H]-(13)C NOE were determined in this study. The C alpha H relaxation measurements were compared to the previously measured (15)N backbone parameters and these are found to be in agreement. For the analysis of the threonine side chain motions, the model of restricted rotational diffusion about the chi(1) dihedral angle was employed [London and Avitabile (1978) J. Am. Chem. Soc., 100, 7159-7165]. We demonstrate that the motion experienced by the ice binding threonine side chains is highly restricted, with an approximate upper limit of less than +/-25 degrees.     

The role of side chain conformational flexibility in surface recognition by Tenebrio molitor antifreeze protein

Two-dimensional nuclear magnetic resonance spectroscopy was used to investigate the flexibility of the threonine side chains in the beta-helical Tenebrio molitor antifreeze protein (TmAFP) at low temperatures. From measurement of the (3)J(alphabeta) (1)H-(1)H scalar coupling constants, the chi(1) angles and preferred rotamer populations can be calculated. It was determined that the threonines on the ice-binding face of the protein adopt a preferred rotameric conformation at near freezing temperatures, whereas the threonines not on the ice-binding face sample many rotameric states. This suggests that TmAFP maintains a preformed ice-binding conformation in solution, wherein the rigid array of threonines that form the AFP-ice interface matches the ice crystal lattice. A key factor in binding to the ice surface and inhibition of ice crystal growth appears to be the close surface-to-surface complementarity between the AFP and crystalline ice, and the lack of an entropic penalty associated with freezing out motions in a flexible ligand.     

Crystal structure of beta-helical antifreeze protein points to a general ice binding model

Reported here is the 2.3 A resolution crystal structure of spruce budworm (Choristoneura fumiferana) antifreeze protein (CfAFP), solved by single anomalous scattering. The structure reveals an extremely regular left-handed beta-helical platform consisting of 15-amino acid loops with a repetitive Thr-X-Thr motif displayed on one of the helix's three faces. This motif results in a two-dimensional array of threonine residues in an identical orientation to those in the nonhomologous, right-handed beta-helical beetle AFP from Tenebrio molitor (TmAFP). The CfAFP structure led us to reevaluate our ice binding model, and the analysis of three possible modes of docking gives rise to a binding mechanism based on surface complementarity. This general mechanism is applicable to both fish and insect AFPs.     

Computational study on the function of water within a beta-helix antifreeze protein dimer and in the process of ice-protein binding

Antifreeze proteins (AFPs) help many organisms protect themselves from freezing in subzero temperatures. The most active AFPs found to date are those from insects, which possess exceptionally regular beta-helical structures. On the ice-binding surface of these proteins, regularly arrayed water molecules are observed within the repeating Thr-Xxx-Thr motif, but the exact role of these water molecules remains unknown. In this work, we have employed a number of computational methods to examine the role of these water molecules in an AFP from Tenebrio molitor (TmAFP). Our investigation involved a combination of molecular and quantum mechanical approaches. Properties such as stability, interaction energy, orbital overlap, and conformational analysis of various systems, including TmAFP-water, TmAFP-water-ice, and TmAFP-ice, were systematically evaluated and compared. The regularly arrayed water molecules were found to remain associated with TmAFP before ice binding, demonstrating that they are an intrinsic part of the protein. These water molecules may assist TmAFP in the process of ice recognition and binding. However, after facilitating the initial stages of ice recognition and binding, these water molecules are excluded in the final formation of the AFP-ice complex. The departure of these water molecules enables a better two-dimensional match between TmAFP and ice. These results agree with experimental observations showing that although these water molecules are aligned with the ice-binding hydroxyl groups of Thr residues in one dimension, they are in fact positioned slightly off in the second dimension, making a good two-dimensional match impossible.     

甲虫抗冻蛋白热滞活性的二维吸附结合模型

甲虫抗冻蛋白是一种具有规则结构的昆虫抗冻蛋白。在相同浓度条件下,甲虫抗冻蛋白比鱼类抗冻蛋白有更高的热滞活性,目前已成为人们重点研究的一类抗冻蛋白。根据甲虫抗冻蛋白的结构特点及其在冰晶表面的吸附模式,应用二维吸附结合模型计算分析了具有6￣11个β-螺旋(β-helix)结构片段的甲虫抗冻蛋白变体分子,得到了它们的热滞活性随溶液浓度变化的规律,特别是热滞活性与甲虫抗冻蛋白的β-螺旋结构片段数的关系。结果显示,抗冻蛋白在冰晶表面的覆盖度是一个影响其热滞活性的重要因素。     

Temperature dependence of dynamics and thermodynamics of the regulatory domain of human cardiac troponin C

Binding of Ca(2+) to the regulatory domain of troponin C (TnC) in cardiac muscle initiates a series of protein conformational changes and modified protein-protein interactions that initiate contraction. Cardiac TnC contains two Ca(2+) binding sites, with one site being naturally defunct. Previously, binding of Ca(2+) to the functional site in the regulatory domain of TnC was shown to lead to a decrease in conformational entropy (TDeltaS) of 2 and 0.5 kcal mol(-1) for the functional and nonfunctional sites, respectively, using (15)N nuclear magnetic resonance (NMR) relaxation studies [Spyracopoulos, L., et al. (1998) Biochemistry 37, 18032-18044]. In this study, backbone dynamics of the Ca(2+)-free regulatory domain are investigated by backbone amide (15)N relaxation measurements at eight temperatures from 5 to 45 degrees C. Analysis of the relaxation measurements yields an order parameter (S(2)) indicating the degree of spatial restriction for a backbone amide H-N vector. The temperature dependence of S(2) allows estimation of the contribution to protein heat capacity from pico- to nanosecond time scale conformational fluctuations on a per residue basis. The average heat capacity contribution (C(p,j)) from backbone conformational fluctuations for regions of secondary structure for the regulatory domain of cardiac apo-TnC is 6 cal mol(-1) K(-1). The average heat capacity for Ca(2+) binding site 1 is larger than that for site 2 by 1.3 +/- 0.8 cal mol(-1) K(-1), and likely represents a mechanism where differences in affinity between Ca(2+) binding sites for EF hand proteins can be modulated.     

The physico-chemical characterization of a boiling stable antifreeze protein from a perennial grass (Lolium perenne)

We have characterized a cold-induced, boiling stable antifreeze protein. This highly active ice recrystallization inhibition protein shows a much lower thermal hysteresis effect and displays binding behavior that is uncharacteristic of any AFP from fish or insects. Ice-binding studies show it binds to the (1 0 1 0) plane of ice and FTIR studies reveal that it has an unusual type of highly beta-sheeted secondary structure. Ice-binding studies of both glycosylated and nonglycosylated expressed forms indicate that it adsorbs to ice through the protein backbone. These results are discussed in light of the currently proposed mechanisms of AFP action.     

Stabilization of protein by replacement of a fluctuating loop: structural analysis of a chimera of bovine alpha-lactalbumin and equine lysozyme

Equine lysozyme is a calcium-binding lysozyme and an evolutional intermediate between non-calcium binding c-type lysozyme and alpha-lactalbumin. We constructed a chimeric protein by substituting the fluctuating loop of bovine alpha-lactalbumin with the D-helix of equine lysozyme. The substitution affects the protection factors not only in the fluctuating loop but also in the antiparallel beta-sheet, the A- and B-helices, and the loop between the B-helix and the beta-sheet. Amide protons in these regions of the chimera are more protected from exchange than are those of bovine alpha-lactalbumin. We used model-free analysis based on 15N nuclear magnetic resonance relaxation measurements to investigate the dynamics of the main chain of the chimera and showed that the fluctuating loop of the chimera is as rigid as three major helices. When we analyzed the chemical shift deviations and backbone HN-H(alpha) scalar coupling constants, we found that the chimera showed an alpha-helical tendency in residues around the fluctuating loop. Our results suggest that the replacement of a highly fluctuating loop in a protein with a rigid structural element in a homologous one may be useful to stabilize the protein structure.     

Solution structure and internal dynamics of NSCP, a compact calcium-binding protein

The solution structure of Nereis diversicolor sarcoplasmic calcium-binding protein (NSCP) in the calcium-bound form was determined by NMR spectroscopy, distance geometry and simulated annealing. Based on 1859 NOE restraints and 262 angular restraints, 17 structures were generated with a rmsd of 0.87 A from the mean structure. The solution structure, which is highly similar to the structure obtained by X-ray crystallography, includes two open EF-hand domains, which are in close contact through their hydrophobic surfaces. The internal dynamics of the protein backbone were determined by studying amide hydrogen/deuterium exchange rates and 15N nuclear relaxation. The two methods revealed a highly compact and rigid structure, with greatly restricted mobility at the two termini. For most of the amide protons, the free energy of exchange-compatible structural opening is similar to the free energy of structural stability, suggesting that isotope exchange of these protons takes place through global unfolding of the protein. Enhanced conformational flexibility was noted in the unoccupied Ca2+-binding site II, as well as the neighbouring helices. Analysis of the experimental nuclear relaxation and the molecular dynamics simulations give very similar profiles for the backbone generalized order parameter (S2), a parameter related to the amplitude of fast (picosecond to nanosecond) movements of N(H)-H vectors. We also noted a significant correlation between this parameter, the exchange rate, and the crystallographic B factor along the sequence.     

Simultaneous processing of solid-state NMR relaxation and 1D-MAS exchange data: the backbone dynamics of free vs. binase-bound barstar

Two types of dynamic solid-state NMR experiments-relaxation and 1D-MAS exchange-were combined for the investigation of the backbone dynamics of a 15% randomly 15N-enriched protein barstar in both free and binase-bound states. The main novelty of this work is a simultaneous quantitative processing of the results of these two types of experiments that we call Simultaneous Relaxation and Exchange Data Analysis (SREDA) approach. It extends the well-known model-free approach such that it permits to discriminate between various motional models (jumps between different sites, wobbling in a cone, etc.). This objective cannot be achieved by analyzing the relaxation or exchange data separately. The SREDA approach was applied to probe a modification of the average backbone dynamics of barstar upon forming a complex with another protein binase. T(1) and off-resonance T(1rho) relaxation times of 15N backbone nuclei were measured at three temperatures between 0 and 45 degrees C, 1D-MAS exchange (CODEX) data were obtained at room temperature within the mixing time range from 0.3 to 200 ms. It has been found that the barstar backbone participates in two molecular processes with correlation times in the 10(-9)-10(-7) and 10(-3)-10(-2) s ranges. Forming the complex with binase results in a significant decrease of the amplitudes of both motions, suggesting that the complex is a more rigid and stable structure than free barstar.     

Flap opening and dimer-interface flexibility in the free and inhibitor-bound HIV protease, and their implications for function.

BACKGROUND: (1)H and (15)N transverse relaxation measurements on perdeuterated proteins are ideally suited for detecting backbone conformational fluctuations on the millisecond-microsecond timescale. The identification of conformational exchange on this timescale by measuring the relaxation of both (1)H and (15)N holds great promise for the elucidation of functionally relevant conformational changes in proteins. RESULTS: We measured the transverse (1)H and (15)N relaxation rates of backbone amides of HIV-1 protease in its free and inhibitor-bound forms. An analysis of these rates, obtained as a function of the effective rotating frame field, provided information about the timescale of structural fluctuations in several regions of the protein. The flaps that cover the active site of the inhibitor-bound protein undergo significant changes of backbone (φ,psi) angles, on the 100 micros timescale, in the free protein. In addition, the intermonomer beta-sheet interface of the bound form, which from protease structure studies appears to be rigid, was found to fluctuate on the millisecond timescale. CONCLUSIONS: We present a working model of the flap-opening mechanism in free HIV-1 protease which involves a transition from a semi-open to an open conformation that is facilitated by interaction of the Phe53 ring with the substrate. We also identify a surprising fluctuation of the beta-sheet intermonomer interface that suggests a structural requirement for maturation of the protease. Thus, slow conformational fluctuations identified by (1)H and (15)N transverse relaxation measurements can be related to the biological functions of proteins.     

Investigation of ligand binding and protein dynamics in Bacillus subtilis chorismate mutase by transverse relaxation optimized spectroscopy-nuclear magnetic resonance

The structural and dynamical consequences of ligand binding to a monofunctional chorismate mutase from Bacillus subtilis have been investigated by solution NMR spectroscopy. TROSY methods were employed to assign 98% of the backbone (1)H(N), (1)H(alpha), (15)N, (13)C', and (13)C(alpha) resonances as well as 86% of the side chain (13)C resonances of the 44 kDa trimeric enzyme at 20 degrees C. This information was used to map chemical shift perturbations and changes in intramolecular mobility caused by binding of prephenate or a transition state analogue to the X-ray structure. Model-free interpretation of backbone dynamics for the free enzyme and its complexes based on (15)N relaxation data measured at 600 and 900 MHz showed significant structural consolidation of the protein in the presence of a bound ligand. In agreement with earlier structural and biochemical studies, substantial ordering of 10 otherwise highly flexible residues at the C-terminus is particularly notable. The observed changes suggest direct contact between this protein segment and the bound ligand, providing support for the proposal that the C-terminus can serve as a lid for the active site, limiting diffusion into and out of the pocket and possibly imposing conformational control over substrate once bound. Other regions of the protein that experience substantial ligand-induced changes also border the active site or lie along the subunit interfaces, indicating that the enzyme adapts dynamically to ligands by a sort of induced fit mechanism. It is believed that the mutase-catalyzed chorismate-to-prephenate rearrangement is partially encounter controlled, and backbone motions on the millisecond time scale, as seen here, may contribute to the reaction barrier.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.     

Backbone dynamics of receptor binding and antigenic regions of a Pseudomonas aeruginosa pilin monomer

Pilin is the major structural protein that forms type IV pili of various pathogenic bacteria, including Pseudomonas aeruginosa. Pilin is involved in attachment of the bacterium to host cells during infection, in the initiation of immune response, and serves as a receptor for a variety of bacteriophage. We have used (15)N nuclear magnetic resonance relaxation measurements to probe the backbone dynamics of an N-terminally truncated monomeric pilin from P. aeruginosa strain K122-4. (15)N-T(1), -T(2), and [(1)H]-(15)N nuclear Overhauser enhancement measurements were carried out at three magnetic field strengths. The measurements were interpreted using the Lipari-Szabo model-free analysis, which reveals the amplitude of spatial restriction for backbone N-NH bond vectors with respect to nano- to picosecond time-scale motions. Regions of well-defined secondary structure exhibited consistently low-amplitude spatial fluctuations, while the terminal and loop regions showed larger amplitude motions in the subnano- to picosecond time-scale. Interestingly, the C-terminal disulfide loop region that contains the receptor binding domain was found to be relatively rigid on the pico- to nanosecond time-scale but exhibited motion in the micro- to millisecond time-scale. It is notable that this disulfide loop displays a conserved antigenic epitope and mediates binding to the asialo-GM(1) cell surface receptor. The present study suggests that a rigid backbone scaffold mediates attachment to the host cell receptor, and also maintains the conformation of the conserved antigenic epitope for antibody recognition. In addition, slower millisecond time-scale motions are likely to be crucial for conferring a range of specificity for these interactions. Characterization of pilin dynamics will aid in developing a detailed understanding of infection, and will facilitate the design of more efficient anti-adhesin synthetic vaccines and therapeutics against pathogenic bacteria containing type IV pili.     

15N NMR relaxation studies of free and ligand-bound human acidic fibroblast growth factor

15N NMR relaxation data have been used to characterize the backbone dynamics of the human acidic fibroblast growth factor (hFGF-1) in its free and sucrose octasulfate (SOS)-bound states. (15)N longitudinal (R(1)), transverse (R(2)) relaxation rates and (1H)-(15)N steady-state nuclear Overhauser effects were obtained at 500 and 600 MHz (at 25 degrees C) for all resolved backbone amide groups using (1)H- detected two-dimensional NMR experiments. Relaxation data were fit to the extended model free dynamics for each NH group. The overall correlation time (tau(m)) for the free and SOS-bound forms were estimated to be 10.4 +/- 1.07 and 11.1 +/- 1.35 ns, respectively. Titration experiments with SOS reveals that the ligand binds specifically to the C-terminal domain of the protein in a 1:1 ratio. Binding of SOS to hFGF-1 is found to induce a subtle conformational change in the protein. Significant conformational exchange (R(ex)) is observed for several residues in the free form of the protein. However, in the SOS-bound form only three residues exhibit significant R(ex) values, suggesting that the dynamics on the micro- to millisecond time scale in the free form is coupled to the cis-trans-proline isomerization. hFGF-1 is a rigid molecule with an average generalized parameter (S(2)) value of 0.89 +/- 0.03. Upon binding to SOS, there is a marked decrease in the overall flexibility (S(2) = 0.94 +/- 0.02) of the hFGF-1 molecule. However, the segment comprising residues 103-111 shows increased flexibility in the presence of SOS. Significant correlation is found between residues that show high flexibility and the putative receptor binding sites on the protein.     

Backbone dynamics of reduced plastocyanin from the cyanobacterium Anabaena variabilis: regions involved in electron transfer have enhanced mobility

The dynamics of the backbone of the electron-transfer protein plastocyanin from the cyanobacterium Anabaena variabilis were determined from the (15)N and (13)C(alpha) R(1) and R(2) relaxation rates and steady-state [(1)H]-(15)N and [(1)H]-(13)C nuclear Overhauser effects (NOEs) using the model-free approach. The (13)C relaxation studies were performed using (13)C in natural abundance. Overall, it is found that the protein backbone is rigid. However, the regions that are important for the function of the protein show moderate mobility primarily on the microsecond to millisecond time scale. These regions are the "northern" hydrophobic site close to the metal site, the metal site itself, and the "eastern" face of the molecule. In particular, the mobility of the latter region is interesting in light of recent findings indicating that residues also on the eastern face of plastocyanins from prokaryotes are important for the function of the protein. The study also demonstrates that relaxation rates and NOEs of the (13)C(alpha) nuclei of proteins are valuable supplements to the conventional (15)N relaxation measurements in studies of protein backbone dynamics.     

Structure, dynamics, and thermodynamics of the structural domain of troponin C in complex with the regulatory peptide 1-40 of troponin I

The structure of the calcium-saturated C-domain of skeletal troponin C (CTnC) in complex with a regulatory peptide comprising residues 1-40 (Rp40) of troponin I (TnI) was determined using nuclear magnetic resonance (NMR) spectroscopy. The solution structure determined by NMR is similar to the structure of the C-domain from intact TnC in complex with TnI(1)(-)(47) determined by X-ray crystallography [Vassylyev, D. G., Takeda, S., Wakatsuki, S., Maeda, K., and Maeda, Y. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 4847-4852]. Changes in the dynamic properties of CTnC.2Ca2+ induced by Rp40 binding were investigated using backbone amide (15)N NMR relaxation measurements. Analysis of NMR relaxation data allows for extraction of motional order parameters on a per residue basis, from which the contribution of changes in picosecond to nanosecond time scale motions to the conformational entropy associated with complex formation can be estimated. The results indicate that binding of Rp40 decreases backbone flexibility in CTnC, particularly at the end of the C-terminal helix. The backbone conformational entropy change (-TDeltaS) associated with binding of Rp40 to CTnC.2Ca2+ determined from (15)N relaxation data is 9.6 +/- 0.7 kcal mol(-1) at 30 degrees C. However, estimation of thermodynamic quantities using a structural approach [Lavigne, P., Bagu, J. R., Boyko, R., Willard, L., Holmes, C. F., and Sykes, B. D. (2000) Protein Sci. 9, 252-264] reveals that the change in solvation entropy upon complex formation is dominant and overcomes the thermodynamic "cost" associated with "stiffening" of the protein backbone upon Rp40 binding. Additionally, backbone amide (15)N relaxation data measured at different concentrations of CTnC.2Ca2+.Rp40 reveal that the complex dimerizes in solution. Fitting of the apparent global rotational correlation time as a function of concentration to a monomer-dimer equilibrium yields a dimerization constant of approximately 8.3 mM.     

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.