The role of backbone conformational heat capacity in protein stability: temperature dependent dynamics of the B1 domain of Streptococcal protein G

Specific anion effect on properties of HRV 3C protease

Specific salts effect is intensively studied from the prospective of modification of different physico-chemical properties of biomacromolecules. Limited knowledge of the specific salts effect on enzymes led us to address the influence of five sodium anions: sulfate, phosphate, chloride, bromide, and perchlorate, on catalytic and conformational properties of human rhinovirus-14 (HRV) 3C protease. The enzyme conformation was monitored by circular dichroism spectrum (CD) and by tyrosines fluorescence. Stability and flexibility of the enzyme have been analyzed by CD in the far-UV region, differential scanning calorimetry and molecular dynamics simulations, respectively. We showed significant influence of the anions on the enzyme properties in accordance with the Hofmeister effect. The HRV 3C protease in the presence of kosmotropic anions, in contrast with chaotropic anions, exhibits increased stability, rigidity. Correlations of stabilization effect of anions on the enzyme with their charge density and the rate constant of the enzyme with the viscosity B-coefficients of anions suggest direct interaction of the anions with HRV 3C protease. The role of stabilization and decreased fluctuation of the polypeptide chain of HRV 3C protease on its activation in the presence of kosmotropic anions is discussed within the frame of the macromolecular rate theory.

Structure and activity of a thermally stable mutant of Acanthamoeba actophorin

Actophorin, which was recently tested for crystallization under microgravity on the International Space Station, was subjected to mutagenesis to identify a construct with improved biophysical properties that were expected to improve the extent of diffraction. First, 20 mutations, including one C-terminal deletion of three residues, were introduced individually into actophorin, resulting in modest increases in thermal stability of between +0.5°C and +2.2°C. All but two of the stabilizing mutants increased both the rates of severing F-actin filaments and of spontaneous polymerization of pyrenyl G-actin in vitro. When the individual mutations were combined into a single actophorin variant, Acto-2, the overall thermal stability was 22°C higher than that of wild-type actophorin. When an inactivating S2P mutation in Acto-2 was restored, Acto-2/P2S was more stable by 20°C but was notably more active than the wild-type protein. The inactivating S2P mutation reaffirms the importance that Ser2 plays in the F-actin-severing reaction. The crystal structure of Acto-2 was solved to 1.7 Å resolution in a monoclinic space group, a first for actophorin. Surprisingly, despite the increase in thermal stability, the extended β-turn region, which is intimately involved in interactions with F-actin, is disordered in one copy of Acto-2 in the asymmetric unit. These observations emphasize the complex interplay among protein thermal stability, function and dynamics.

Conformational Entropy from Mobile Bond Vectors in Proteins: A Viewpoint that Unifies NMR Relaxation Theory and Molecular Dynamics Simulation Approaches

A new method for determining conformational entropy in proteins is reported. Proteins prevail as conformational ensembles, p ∝ exp(-u). By selecting a bond vector (e.g., N-H) as a conformation representative, molecular dynamics simulations can provide (relative to a reference structure) p as approximate Boltzmann probability density and u as N-H potential of mean force (POMF). The latter is as accurate as implied by the force field but statistical in character; this limits the insights it can provide and its utilization. Conformational entropy is given exclusively by u. Deriving it from POMFs renders it accurate but statistical in character. Previously, we devised explicit (i.e., analytical but not exact) potentials made of Wigner functions, D_0K^L, with L ≤ 4, which closely resemble the corresponding POMFs in form; hence, they also approach the latter in accuracy. Such potentials can be beneficially characterized/compared in terms of composition, symmetry, and associated order parameters. In this study, we develop a method for deriving conformational entropy from them, which also features the benefits specified above. The method developed is applied to the dimerization of the Rho GTPase-binding domain of plexin-B1. Insights into local ordering, entropy compensation, and features of allostery are gained. In previous work, we developed the slowly relaxing local structure (SRLS) approach for the analysis of NMR relaxation from restricted bond vector motion in proteins. SRLS comprises explicit (restricting) potentials of the kind developed here. It also comprises diffusion tensors describing the local motion and related features of local geometry. The complete model fits experimental data. In future work, the explicit potentials developed here will be inserted unchanged in SRLS-based data fitting, thereby improving the picture of structural dynamics. Given that SRLS is unique in featuring potentials that can closely approach the corresponding POMFs in accuracy, the present study is an important step toward generally improving protein dynamics by NMR relaxation.

Structural and dynamic origins of ESR lineshapes in spin-labeled GB1 domain: the insights from spin dynamics simulations based on long MD trajectories

Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 μs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1 and χ2 are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.

Insulin fibrillation: toward strategies for attenuating the process

The environmental factors affecting the rate of insulin fibrillation. The factors are representative.

Cooperative energetic effects elicited by the yeast Shwachman-Diamond syndrome protein (Sdo1) and guanine nucleotides modulate the complex conformational landscape of the elongation factor-like 1 (Efl1) GTPase

One of the final maturation steps of the large ribosomal subunit requires the joint action of the elongation factor-like 1 (human EFL1, yeast Efl1) GTPase and the Shwachman-Diamond syndrome protein (human SBDS, yeast Sdo1) to release the eukaryotic translation initiation factor 6 (human eIF6, yeast Tif6) and allow the assembly of mature ribosomes. EFL1 function is driven by conformational changes. However, the nature of such conformational changes or the mechanism by which they are prompted are still largely unknown. In previous studies, it has been established that this GTPase interacts with its cofactor in solution in an inverted orientation with respect to the binding mode derived from 60S ribosome subunit cryo-EM data. To shed new light on this conundrum, we characterized calorimetrically the energetic basis describing the recognition of Efl1 to GT(D)P, Sdo1 and their intercommunication in solution. A structural-based analysis of the binding signatures indicates that Efl1 has a large structural flexibility. The mutual effects of Sdo1 and nucleotides on Efl1 modulate in a very specific and robust way the complex conformational landscape of Efl1, resembling the behavior observed with other GTPases and their cofactors.

Entropic contribution to enhanced thermal stability in the thermostable P450 CYP119

The enhanced thermostability of thermophilic proteins with respect to their mesophilic counterparts is often attributed to the enthalpy effect, arising from strong interactions between protein residues. Intuitively, these strong interresidue interactions will rigidify the biomolecules. However, the present work utilizing neutron scattering and solution NMR spectroscopy measurements demonstrates a contrary example that the thermophilic cytochrome P450, CYP119, is much more flexible than its mesophilic counterpart, CYP101A1, something which is not apparent just from structural comparison of the two proteins. A mechanism to explain this apparent contradiction is that higher flexibility in the folded state of CYP119 increases its conformational entropy and thereby reduces the entropy gain during denaturation, which will increase the free energy needed for unfolding and thus stabilize the protein. This scenario is supported by thermodynamic data on the temperature dependence of unfolding free energy, which shows a significant entropic contribution to the thermostability of CYP119 and lends an added dimension to enhanced stability, previously attributed only to presence of aromatic stacking interactions and salt bridge networks. Our experimental data also support the notion that highly thermophilic P450s such as CYP119 may use a mechanism that partitions flexibility differently from mesophilic P450s between ligand binding and thermal stability.

Heat Capacity Changes for Transition-State Analogue Binding and Catalysis with Human 5'-Methylthioadenosine Phosphorylase

Human 5'-methylthioadenosine phosphorylase (MTAP) catalyzes the phosphorolysis of 5'-methylthioadenosine (MTA). Its action regulates cellular MTA and links polyamine synthesis to S-adenosylmethionine (AdoMet) salvage. Transition state analogues with picomolar dissociation constants bind to MTAP in an entropically driven process at physiological temperatures, suggesting increased hydrophobic character or dynamic structure for the complexes. Inhibitor binding exhibits a negative heat capacity change (-ΔC_p), and thus the changes in enthalpy and entropy upon binding are strongly temperature-dependent. The ΔC_p of inhibitor binding by isothermal titration calorimetry does not follow conventional trends and is contrary to that expected from the hydrophobic effect. Thus, ligands of increasing hydrophobicity bind with increasing values of ΔC_p. Crystal structures of MTAP complexed to transition-state analogues MT-DADMe-ImmA, BT-DADMe-ImmA, PrT-ImmA, and a substrate analogue, MT-tubercidin, reveal similar active site contacts and overall protein structural parameters, despite large differences in ΔC_p for binding. In addition, ΔC_p values are not correlated with K_d values. Temperature dependence of presteady state kinetics revealed the chemical step for the MTAP reaction to have a negative heat capacity for transition state formation (-ΔC_p^‡). A comparison of the ΔC_p^‡ for MTAP presteady state chemistry and ΔC_p for inhibitor binding revealed those transition-state analogues most structurally and thermodynamically similar to the transition state. Molecular dynamics simulations of MTAP apoenzyme and complexes with MT-DADMe-ImmA and MT-tubercidin show small, but increased dynamic motion in the inhibited complexes. Variable temperature CD spectroscopy studies for MTAP-inhibitor complexes indicate remarkable protein thermal stability (to T_m = 99 °C) in complexes with transition-state analogues.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Site-Specific Internal Motions in GB1 Protein Microcrystals Revealed by 3D ²H-¹³C-¹³C Solid-State NMR Spectroscopy

(2)H quadrupolar line shapes deliver rich information about protein dynamics. A newly designed 3D (2)H-(13)C-(13)C solid-state NMR magic angle spinning (MAS) experiment is presented and demonstrated on the microcrystalline β1 immunoglobulin binding domain of protein G (GB1). The implementation of (2)H-(13)C adiabatic rotor-echo-short-pulse-irradiation cross-polarization (RESPIRATION CP) ensures the accuracy of the extracted line shapes and provides enhanced sensitivity relative to conventional CP methods. The 3D (2)H-(13)C-(13)C spectrum reveals (2)H line shapes for 140 resolved aliphatic deuterium sites. Motional-averaged (2)H quadrupolar parameters obtained from the line-shape fitting identify side-chain motions. Restricted side-chain dynamics are observed for a number of polar residues including K13, D22, E27, K31, D36, N37, D46, D47, K50, and E56, which we attribute to the effects of salt bridges and hydrogen bonds. In contrast, we observe significantly enhanced side-chain flexibility for Q2, K4, K10, E15, E19, N35, N40, and E42, due to solvent exposure and low packing density. T11, T16, and T17 side chains exhibit motions with larger amplitudes than other Thr residues due to solvent interactions. The side chains of L5, V54, and V29 are highly rigid because they are packed in the core of the protein. High correlations were demonstrated between GB1 side-chain dynamics and its biological function. Large-amplitude side-chain motions are observed for regions contacting and interacting with immunoglobulin G (IgG). In contrast, rigid side chains are primarily found for residues in the structural core of the protein that are absent from protein binding and interactions.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R² = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.

Chimeric Avidin--NMR structure and dynamics of a 56 kDa homotetrameric thermostable protein

Chimeric avidin (ChiAVD) is a product of rational protein engineering remarkably resistant to heat and harsh conditions. In quest of the fundamentals behind factors affecting stability we have elucidated the solution NMR spectroscopic structure of the biotin–bound form of ChiAVD and characterized the protein dynamics through ¹⁵N relaxation and hydrogen/deuterium (H/D) exchange of this and the biotin–free form. To surmount the challenges arising from the very large size of the protein for NMR spectroscopy, we took advantage of its high thermostability. Conventional triple resonance experiments for fully protonated proteins combined with methyl–detection optimized experiments acquired at 58°C were adequate for the structure determination of this 56 kDa protein. The model–free parameters derived from the ¹⁵N relaxation data reveal a remarkably rigid protein at 58°C in both the biotin–bound and the free forms. The H/D exchange experiments indicate a notable increase in hydrogen protection upon biotin binding.

Proton-decoupled CPMG: a better experiment for measuring (15)N R2 relaxation in disordered proteins

(15)N R2 relaxation is one of the most informative experiments for characterization of intrinsically disordered proteins (IDPs). Small changes in nitrogen R2 rates are often used to determine how IDPs respond to various biologically relevant perturbations such as point mutations, posttranslational modifications and weak ligand interactions. However collecting high-quality (15)N relaxation data can be difficult. Of necessity, the samples of IDPs are often prepared with low protein concentration and the measurement time can be limited because of rapid sample degradation. Furthermore, due to hardware limitations standard experiments such as (15)N spin-lock and CPMG can sample the relaxation decay only to ca. 150ms. This is much shorter than (15)N T2 times in disordered proteins at or near physiological temperature. As a result, the sampling of relaxation decay profiles in these experiments is suboptimal, which further lowers the precision of the measurements. Here we report a new implementation of the proton-decoupled (PD) CPMG experiment which allows one to sample (15)N R2 relaxation decay up to ca. 0.5-1s. The new experiment has been validated through comparison with the well-established spin-lock measurement. Using dilute samples of denatured ubiquitin, we have demonstrated that PD-CPMG produces up to 3-fold improvement in the precision of the data. It is expected that for intrinsically disordered proteins the gains may be even more substantial. We have also shown that this sequence has a number of favorable properties: (i) the spectra are recorded with narrow linewidth in nitrogen dimension; (ii) (15)N offset correction is small and easy to calculate; (iii) the experiment is immune to various spurious effects arising from solvent exchange; (iv) the results are stable with respect to pulse miscalibration and rf field inhomogeneity; (v) with minimal change, the pulse sequence can also be used to measure R2 relaxation of (15)N(ε) spins in arginine side chains. We anticipate that the new experiment will be a valuable addition to the NMR toolbox for studies of IDPs.

Change in heat capacity for enzyme catalysis determines temperature dependence of enzyme catalyzed rates

The increase in enzymatic rates with temperature up to an optimum temperature (Topt) is widely attributed to classical Arrhenius behavior, with the decrease in enzymatic rates above Topt ascribed to protein denaturation and/or aggregation. This account persists despite many investigators noting that denaturation is insufficient to explain the decline in enzymatic rates above Topt. Here we show that it is the change in heat capacity associated with enzyme catalysis (ΔC(‡)p) and its effect on the temperature dependence of ΔG(‡) that determines the temperature dependence of enzyme activity. Through mutagenesis, we demonstrate that the Topt of an enzyme is correlated with ΔC(‡)p and that changes to ΔC(‡)p are sufficient to change Topt without affecting the catalytic rate. Furthermore, using X-ray crystallography and molecular dynamics simulations we reveal the molecular details underpinning these changes in ΔC(‡)p. The influence of ΔC(‡)p on enzymatic rates has implications for the temperature dependence of biological rates from enzymes to ecosystems.

NMR- and circular dichroism-monitored lipid binding studies suggest a general role for the FATC domain as membrane anchor of phosphatidylinositol 3-kinase-related kinases (PIKK)

The FATC domain is shared by all members of the family of phosphatidylinositol-3 kinase-related kinases (PIKKs). It has been shown that the FATC domain plays an important role for the regulation of each PIKK. However, other than an involvement in protein-protein interactions, a common principle for the action of the FATC domain has not been detected. A detailed characterization of the structure and lipid binding properties of the FATC domain of the Ser/Thr kinase target of rapamycin (TOR) revealed that it contains a redox-sensitive membrane anchor in its C terminus. Because the C-terminal regions of the FATC domains of all known PIKKs are rather hydrophobic and especially rich in aromatic residues, we examined whether the ability to interact with lipids and membranes might be a general property. Here, we present the characterization of the interactions with lipids and different membrane mimetics for the FATC domains of human DNA-PKcs, human ATM, human ATR, human SMG-1, and human TRRAP by NMR and CD spectroscopy. The data indicate that all of these can interact with different membrane mimetics and may have different preferences only for membrane properties such as surface charge, curvature, and lipid packing. The oxidized form of the TOR FATC domain is well structured overall and forms an α-helix that is followed by a disulfide-bonded loop. In contrast, the FATC domains of the other PIKKs are rather unstructured in the isolated form and only significantly populate α-helical secondary structure upon interaction with membrane mimetics.

Thermal coefficients of the methyl groups within ubiquitin

Physiological processes such as protein folding and molecular recognition are intricately linked to their dynamic signature, which is reflected in their thermal coefficient. In addition, the local conformational entropy is directly related to the degrees of freedom, which each residue possesses within its conformational space. Therefore, the temperature dependence of the local conformational entropy may provide insight into understanding how local dynamics may affect the stability of proteins. Here, we analyze the temperature dependence of internal methyl group dynamics derived from the cross-correlated relaxation between dipolar couplings of two CH bonds within ubiquitin. Spanning a temperature range from 275 to 308 K, internal methyl group dynamics tend to increase with increasing temperature, which translates to a general increase in local conformational entropy. With this data measured over multiple temperatures, the thermal coefficient of the methyl group order parameter, the characteristic thermal coefficient, and the local heat capacity were obtained. By analyzing the distribution of methyl group thermal coefficients within ubiquitin, we found that the N-terminal region has relatively high thermostability. These results indicate that methyl groups contribute quite appreciably to the total heat capacity of ubiquitin through the regulation of local conformational entropy.

Temperature dependence of molecular interactions involved in defining stability of glutamine binding protein and its complex with L-glutamine

The temperature dependence of dynamic parameters derived from nuclear magnetic resonance (NMR) relaxation data is related to conformational entropy of the system under study. This provides information such as macromolecules stability and thermodynamics of ligand binding. We studied the temperature dependence of NMR order parameter of glutamine binding protein (GlnBP), a periplasmic binding protein (PBP) highly specific to L-glutamine associated with its ABC transporter, with the goal of elucidating the dynamical differences between the respective ligand bound and free forms. We found that the protein-ligand interaction, which is stabilized at higher temperature, has a striking effect on the stability of the hydrophobic core of the large domain of GlnBP. Moreover, in contrast to what was found for less specific PBPs, the decreasing backbone motion of the hinge region at increasing temperature supports the idea that the likelihood that GlnBP can adopt a ligand free closed conformation in solution diminishes at higher temperatures. Our results support the induced-fit model as mode of action for GlnBP. In addition, we found that the backbones of residues involved in a salt bridge do not necessarily become more rigid as the temperature rises as it was previously suggested [Vinther, J. M., et al. (2011) J. Am. Chem. Soc., 133, 271-278]. Our results show that for this to happen these residues have to also directly interact with a region of the protein that is becoming more rigid as the temperature increases.

SRLS analysis of 15N spin relaxation from E. coli ribonuclease HI: the tensorial perspective

15N–H relaxation parameters from ribonuclease HI (RNase H), acquired in previous work at magnetic fields of 14.1 and 18.8 T, and at 300 K, are analyzed with the mode-coupling slowly relaxing local structure (SRLS) approach. In accordance with standard theoretical treatments of restricted motions, SRLS approaches N-H bond dynamics from a tensorial perspective. As shown previously, a physically adequate description of this phenomenon has to account for the asymmetry of the local spatial restrictions. So far, we used rhombic local ordering tensors; this is straightforward but computationally demanding. Here, we propose substantiating the asymmetry of the local spatial restrictions in terms of tilted axial local ordering (S) and local diffusion (D2) tensors. Although less straightforward, this description provides physically sound structural and dynamic information and is efficient computationally. We find that the local order parameter, S(0)2, is on average 0.89 (0.84, and may be as small as 0.6) for the secondary structure elements (loops). The main local ordering axis deviates from the C(i-1)α-C(i)α axis by less than 6°. At 300 K, D(2,perpendicular) is virtually the same as the global diffusion rate, D1 = 1.8 × 10(7) s(-1). The correlation time 1/6D(2,parallel) ranges from 3-125 (208-344) ps for the secondary structure elements (loops) and is on average 125 ps for the C-terminal segment. The main local diffusion axis deviates from the N-H bond by less than 2° (10°) for the secondary structure elements (loops). An effective data-fitting protocol, which leads in most cases to unambiguous results with limited uncertainty, has been devised. A physically sound and computationally effective methodology for analyzing 15N relaxation in proteins, that provides a new picture of N–H bond structural dynamics in proteins, has been set forth.

Water proton spin saturation affects measured protein backbone 15N spin relaxation rates

Protein backbone 15N NMR spin relaxation rates are useful in characterizing the protein dynamics and structures. To observe the protein nuclear-spin resonances a pulse sequence has to include a water suppression scheme. There are two commonly employed methods, saturating or dephasing the water spins with pulse field gradients and keeping them unperturbed with flip-back pulses. Here different water suppression methods were incorporated into pulse sequences to measure 15N longitudinal T1 and transversal rotating-frame T1ρ spin relaxation. Unexpectedly the 15N T1 relaxation time constants varied significantly with the choice of water suppression method. For a 25-kDa Escherichiacoli. glutamine binding protein (GlnBP) the T1 values acquired with the pulse sequence containing a water dephasing gradient are on average 20% longer than the ones obtained using a pulse sequence containing the water flip-back pulse. In contrast the two T1ρ data sets are correlated without an apparent offset. The average T1 difference was reduced to 12% when the experimental recycle delay was doubled, while the average T1 values from the flip-back measurements were nearly unchanged. Analysis of spectral signal to noise ratios (s/n) showed the apparent slower 15N relaxation obtained with the water dephasing experiment originated from the differences in 1HN recovery for each relaxation time point. This in turn offset signal reduction from 15N relaxation decay. The artifact becomes noticeable when the measured 15N relaxation time constant is comparable to recycle delay, e.g., the 15N T1 of medium to large proteins. The 15N relaxation rates measured with either water suppression schemes yield reasonable fits to the structure. However, data from the saturated scheme results in significantly lower Model-Free order parameters (<S2>=0.81) than the non-saturated ones (<S2>=0.88), indicating such order parameters may be previously underestimated.

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2-16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.

Enhanced stability of a protein with increasing temperature

The unusual stability of a structured but locally flexible protein, human growth hormone (hGH) at pH 2.7, was investigated using the temperature dependence of the nanosecond-picosecond dynamics of the backbone amide groups obtained from (15)N NMR relaxation data. It is found that the flexibility of the backbone of the helices decreases with temperature in the range from 24 °C to ∼40 °C, corresponding to an increasing stability. A concomitant increase with temperature of the electrostatic interactions between charged residues forming an interhelical network of salt bridges at the center of the four-helix core suggests that these interactions give rise to the decreasing flexibility and increasing stability of the protein. However, numerous hydrophobic interactions in the interior of the four-helix core may also contribute. Above ∼40 °C, where the thermal energy overcomes the electrostatic and hydrophobic interactions, a substantial increase in the flexibility of the helix backbones results in a highly positive contribution from the local conformational heat capacity, C(p, conf), of the helix backbones to the total heat capacity, C(p), of the protein. This reduces the change in heat capacity upon unfolding, ΔC(p), increases the change in the Gibbs free energy, ΔG(unfold), and stabilizes the protein at high temperatures. A similar decrease in flexibility is found near other salt bridges in hGH and in Calmodulin and may be of general importance for the thermostability of proteins and, in particular, of the salt bridge intensive thermophilic proteins.

Dimerization of protein G B1 domain at low pH: a conformational switch caused by loss of a single hydrogen bond

A number of signals in the NMR spectrum of the B1 domain of staphylococcal protein G (GB1) show a chemical shift dependence on the concentration of the protein at pH 3 but not at neutral pH, implying the existence of self-association at low pH. NMR backbone relaxation experiments show that GB1 undergoes a slow conformational exchange at pH 3, which is not seen at higher pH. Analysis of relaxation dispersion experiments yields a self-association constant of 50 mM, and shows that (15)N chemical shift changes in the dimer interface are up to 3 ppm. The shift changes measured from concentration-dependent HSQC spectra and from relaxation dispersion show good consistency. Measurements of chemical shifts as a function of pH show that a hydrogen bond between the sidechains of Asp44 and Gln40 is broken when Asp44 is protonated, and that loss of this hydrogen bond leads to the breaking of the (i, i + 4) backbone helical hydrogen bond from Asp44 HN to Gln40 O, and therefore to a loss of two residues from the C-terminal end of the helix. This weakens the helix structure and facilitates the loss of further helical structure thus permitting dimerization, which is suggested to occur in the same way as observed for the A42F mutant of GB1 (Jee et al., Proteins 2007;71:1420-1431), by formation of an antiparallel beta-sheet between the edge strands 2 in two monomers. The monomer/dimer ratio is thus a finely balanced equilibrium even in the wild type protein.

Structural basis for the association of the redox-sensitive target of rapamycin FATC domain with membrane-mimetic micelles

The target of rapamycin (TOR) is a conserved eukaryotic Ser/Thr kinase that regulates cellular growth in response to the nutrient and energy state. TOR signaling plays an important role in the development of diseases such as cancer, obesity, and diabetes and in different redox-sensitive processes (hypoxia, apoptosis, and aging). Because TOR has been detected at different cellular membranes and in the nucleus, its localization may influence the specific signaling readout. To better understand how TOR can associate with different membranes, the lipid-binding properties of the redox-sensitive yeast TOR1 FATC domain (y1fatc) have been characterized by solution NMR spectroscopy. Binding studies with different lipids indicate that y1fatc interacts specifically with a membrane-mimetic environment but appears not to recognize a specific lipid headgroup. In both, the structures of oxidized and reduced micelle-bound y1fatc, residues Ile-2456 to Trp-2470 of the lipid-binding motif form a hydrophobic bulb that has a rim of charged residues. The diffusion constants for both micelle-bound states are consistent with the rotational correlation times from the analysis of the (15)N relaxation data. Based on the K(d) values, the oxidized form (K(d) approximately 0.31 mm) binds dodecyl phosphocholine micelles slightly tighter than the reduced form (K(d) approximately 1.86 mM). Binding studies with y1fatc in which one or both tryptophans (Trp-2466 and Trp-2470) were replaced by alanine suggest that these residues are important for the exact positioning in the membrane and that the other aromatic (His-2462, Tyr-2463, and Phe-2469) and aliphatic residues (Ile-2456, Leu-2459, Ile-2464, and Pro-2468) in the lipid-binding motif contribute significantly to the affinity.

Evidence for domain motion in proteins affecting global diffusion properties: a nuclear magnetic resonance study

The rotational diffusion of proteins is an important hydrodynamic property. Compact protein structures were found previously to exhibit hydration layer viscosity, etaloc, higher than the viscosity of bulk water, eta. This implies an apparent activation energy for rotational diffusion higher than the activation energy of water viscosity, Eeta=15.4+/-0.3 kJ/mol. In this study we examine etaloc of internally mobile proteins using 15N spin relaxation methods. We also examine the activation enthalpy, DeltaH#, and activation entropy, DeltaS#, for rotational diffusion. Of particular relevance are internally mobile ligand-free forms and compact ligand-bound forms of multidomain proteins. Adenylate kinase (AKeco) and Ca2+-calmodulin (Ca2+-CaM) are typical examples. For AKeco (Ca2+-CaM) we find that DeltaH# is 14.5+/-0.5 (15.7+/-0.4) kJ/mol. For the complex of AKeco with the inhibitor AP5A (the complex of Ca2+-CaM with the peptide smMLCKp), we find that DeltaH# is 18.1+/-0.7 (18.2+/-0.5) kJ/mol. The internally mobile outer surface protein A has DeltaH#=12.6+/-0.8 kJ/mol, and the compact protein Staphylococcal nuclease has DeltaH#=18.8+/-0.6 kJ/mol. For the internally mobile and compact proteins studied, <|DeltaS(|> equals 62+/-7 J/(mol K) and 44+/-5 J/(mol K), respectively. The fact is that etaloc>eta (DeltaH#>Eeta) for compact proteins was ascribed previously to electrostatic interactions between surface sites and water rigidifying the hydration layer. We find herein that obliteration of these interactions by domain motion leads to etaloc approximately eta, DeltaH# approximately Eeta, and large activation entropy for internally mobile protein structures.

Structural and motional contributions of the Bacillus subtilis ClpC N-domain to adaptor protein interactions

The AAA(+) (ATPases associated with a variety of cellular activities) superfamily protein ClpC is a key regulator of cell development in Bacillus subtilis. As part of a large oligomeric complex, ClpC controls an array of cellular processes by recognizing, unfolding, and providing misfolded and aggregated proteins as substrates for the ClpP peptidase. ClpC is unique compared to other HSP100/Clp proteins, as it requires an adaptor protein for all fundamental activities. The NMR solution structure of the N-terminal repeat domain of ClpC (N-ClpCR) comprises two structural repeats of a four-helix motif. NMR experiments used to map the MecA adaptor protein interaction surface of N-ClpCR reveal that regions involved in the interaction possess conformational flexibility and conformational exchange on the microsecond-to-millisecond timescale. The electrostatic surface of N-ClpCR differs substantially from the N-domain of Escherichia coli ClpA and ClpB, suggesting that the electrostatic surface characteristics of HSP100/Clp N-domains may play a role in adaptor protein and substrate interaction specificity, and perhaps contribute to the unique adaptor protein requirement of ClpC.

Remote changes in the dynamics of the phosphotyrosine-binding domain of insulin receptor substrate-1 induced by phosphopeptide binding

Protein structures consist of cooperative networks of interactions that can extend substantial distances across the protein structure and have significant consequences for stability and function. To characterize such interaction networks within the phosphotyrosine-binding domain of insulin receptor substrate-1 (IRS-PTB), we have used NMR relaxation to determine the dynamics of backbone amide groups, side chain methyl groups, and tryptophan side chain indole groups in IRS-PTB, both in the absence and in the presence of a bound phosphotyrosine-containing peptide. Although there are minimal differences between the structures of apo and peptide-bound forms, primarily localized close to the peptide binding site, we observe significant changes in dynamics for residues located remotely from the binding site. These residues are clustered together and appear to constitute a network of interactions that is coupled to the peptide binding site through intervening residues.

Conformational studies of the alpha-helical 28-43 fragment of the B3 domain of the immunoglobulin binding protein G from Streptococcus

To determine whether the alpha-helix in the B3 immunoglobulin binding domain of protein G from group G Streptococcus has conformational stability as an isolated fragment, we carried out a CD and NMR study of the 16-residue peptide in solution corresponding to this alpha-helix. Based on two-dimensional H-NMR spectra recorded at three different temperatures (283, 305, and 313 K), it was found that this peptide is mostly unstructured in water at these temperatures. Weak signals corresponding to i,i+3 or i,i+4 interactions, which are characteristic of formation of turn-like structures, were observed in the ROE spectra at all temperatures. The absence of a stable three-dimensional structure of the investigated peptide supports an earlier study (Blanco and Serrano, Eur J Biochem 1995, 230, 634-649) of a possible mechanism for folding of other (B1 and B2) immunoglobulin binding domains of Protein G.

Comparison of the backbone dynamics of a natural and a consensus designed 3-TPR domain

The tetratricopeptide repeat (TPR) is a 34-amino acid helix-turn-helix motif that occurs in tandem arrays in numerous proteins. Here we compare the backbone dynamics of a natural 3-repeat TPR domain, from the protein UBP, with the behavior of a designed protein CTPR3, which consists of three identical consensus TPR units. Although the three tandem TPR repeats in both CTPR3 and UBP behave as a single unit, with no evidence of independent repeat motions, the data indicate that certain positions in UBP are significantly more flexible than are the corresponding positions in CTPR3. Most of the dynamical changes occur at or adjacent to positions that are involved in intra-repeat packing interactions. These observations lead us to suggest that the three-TPR domain of UBP does not incorporate optimized packing, compared to that seen in the idealized CTPR. The natural TPR domain is not only less stable overall than CTPR3, but also presents increased local flexibility at the positions where the sequences differs from the conserved consensus.

Manipulating biopolymer dynamics by anisotropic nanoconfinement

How the geometry of nanosized confinement affects dynamics of biomaterials is interesting yet poorly understood. An elucidation of structural details upon nanosized confinement may benefit manufacturing pharmaceuticals in biomaterial sciences and medicine. The behavior of biopolymers in nanosized confinement is investigated using coarse-grained models and molecular simulations. Particularly, we address the effects of shapes of a confinement on protein-folding dynamics by measuring folding rates and dissecting structural properties of the transition states in nanosized spheres and ellipsoids. We find that when the form of a confinement resembles the geometrical properties of the transition states, the rates of folding kinetics are most enhanced. This knowledge of shape selectivity in identifying optimal conditions for reactions will have a broad impact in nanotechnology and pharmaceutical sciences.

Temperature dependence of the NMR generalized order parameter

The study of protein conformational dynamics is motivated in large part by a desire to understand the forces present at different sites throughout the molecular structure. The generalized order parameter determined by NMR spectroscopy has played a central role in the study of protein dynamics on the picosecond-nanosecond time scale. A modeling procedure is presented for analysis of the temperature dependence of the generalized order parameter that extends a previous analysis (Massi and Palmer, J Am Chem Soc 2003;125:11158-11159). As part of this procedure, the potential of mean force is characterized for the N-H bond vectors of the protein backbone. This procedure accounts for the observed temperature dependence of the generalized order parameter in a representative data set from the B1 domain of Streptococcal protein G (Seewald, Pichumani, Stowell, Tibbals, Regan, and Stone, Protein Sci 2000;9:1177-1193). The results indicate a general trend, in which the force constants associated with the potential of mean force decrease with increasing temperature. The analysis also provides evidence for variations in the potential of mean force for different secondary structural elements.

NMR characterization of the energy landscape of SUMO-1 in the native-state ensemble

Characterizing the low energy excited states in the energy landscape of a protein is one of the exciting and demanding problems in structural biology at the present time. These describe the adaptability of the protein structure to external perturbations. In this context, we used here non-linear dependence of amide proton chemical shifts on temperature to identify residues accessing alternative conformations in SUMO-1 in the native state as well as in the near-native states created by sub-denaturing concentrations of urea. The number of residues accessing alternative conformations increases and the profiles of curved temperature dependence also change with increasing urea concentration. In every case these alternative conformations lie within 2 kcal/mol from the ground state, and are separated from it by low energy barriers. The residues that access alternative conformations span the length of the protein chain but are located at particular regions on the protein structure. These include many of the loops, beta2 and beta5 strands, and some edges of the helices. We observed that some of the regions of the protein structure that exhibit such fluctuations coincide with the protein's binding surfaces with different substrate like GTPase effector domain (GED) of dynamin, SUMO binding motifs (SBM), E1 (activating enzyme, SAE1/SAE2) and E2 (conjugating enzyme, UBC9) enzymes of sumoylation machinery, reported earlier. We speculate that this would have significant implications for the binding of diversity of targets by SUMO-1 for the variety of functions it is involved in.

Cross-Correlated Relaxation between H1‘ Chemical Shift Anisotropy and H1‘−H2‘ Dipolar Relaxation Mechanisms in Ribonucleosides: Application to the Characterization of Their Anomeric Configuration

Cross-correlated nuclear spin relaxation between 1H chemical shift anisotropy (CSA) and 1H-1H dipolar relaxation mechanisms in ribonucleosides in solution phase are observed and used to identify their anomeric configuration. Only alpha-ribonucleosides showed the presence of cross-correlated spin relaxation through differential spin-lattice relaxation (T1) of the H1' doublet. Dependence of the magnitude and the orientation of the H1' CSA tensor values on the glycosidic torsion angle and the fast time-scale internal motions present in the ribose moiety play a significant role in the characterization of the anomeric configuration of the nucleosides via cross-correlated relaxation.

An Experimental Investigation of Conformational Fluctuations in Proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.

Analysis of side chain mobility among protein G B1 domain mutants with widely varying stabilities

"Host-guest" studies of the B1 domain from Streptococcal protein G have been used previously to establish a thermodynamic scale for the beta-sheet-forming propensities of the 20 common amino acids. To investigate the contribution of side chain conformational entropy to the relative stabilities of B1 domain mutants, we have determined the dynamics of side chain methyl groups in 10 of the 20 mutants used in a previous study. Deuterium relaxation rates were measured using two-dimensional NMR techniques for 13CH2D groups. Analysis of the relaxation data using the Lipari-Szabo model-free formalism showed that mutations introduced at the guest position caused small but statistically significant changes in the methyl group dynamics. In addition, there was a low level of covariation of the Lipari-Szabo order parameters among the 10 mutants. The variations in conformational free energy estimated from the order parameters were comparable in magnitude to the variations in global stability of the 10 mutants but did not correlate with the global stability of the domain or with the structural properties of the guest amino acids. The data support the view that conformational entropy in the folded state is one of many factors that can influence the folding thermodynamics of proteins.

Gated electron transfers and electron pathways in azurin: a NMR dynamic study at multiple fields and temperatures

Dynamic properties of electron transfer pathways in a small blue copper cupredoxin are explored using an extensive 15N NMR relaxation study of reduced Pseudomonas aeruginosa azurin at four magnetic fields (500-900 MHz) and at two temperatures chosen well below the melting point of the protein. Following a careful model-free analysis, several protein regions with different dynamic regimes are identified. Nanosecond time-scale mobility characterizes various residues of the hydrophobic surface patch believed to mark the natural entry point for electrons, notably the surface-exposed copper-ligand His117. These findings are consistent with a gated electron transfer process according to the "dynamic docking" model. Residues 47-49 along intramolecular pathways of electrons show rigidity that is remarkably conserved when increasing the temperature. Three different conformational exchange processes were observed in the millisecond range, one near the only disulfide bridge in the molecule and two near the copper ion. The latter two processes are consistent with previous data such as crystal structures at various pH values and NMR relaxation dispersion experiments; they may indicate an additional gated electron transfer mechanism at slower time-scales.

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.