Hydrogen exchange kinetics of RNase A and the urea:TMAO paradigm

A key paradigm in the biology of adaptation holds that urea affects protein function by increasing the fluctuations of the native state, while trimethylamine N-oxide (TMAO) affects function in the opposite direction by decreasing the normal fluctuations of the native ensemble. Using urea and TMAO separately and together, hydrogen exchange (HX) studies on RNase A at pH* 6.35 were used to investigate the basic tenets of the urea:TMAO paradigm. TMAO (1 M) alone decreases HX rate constants of a select number of sites exchanging from the native ensemble, and low urea alone increases the rate constants of some of the same sites. Addition of TMAO to urea solutions containing RNase A also suppresses HX rate constants. The data show that urea and TMAO independently or in combination affect the dynamics of the native ensemble in opposing ways. The results provide evidence in support of the counteraction aspect of the urea:TMAO paradigm linking structural dynamics with protein function in urea-rich organs and organisms. RNase A is so resistant to urea denaturation at pH* 6.35 that even in the presence of 4.8 M urea, the native ensemble accounts for >99.5% of the protein. An essential test, devised to determine the HX mechanism of exchangeable protons, shows that over the 0-4.8 M urea concentration range nearly 80% of all observed sites convert from EX2 to EX1. The slow exchange sites are all EX1; they do not exhibit global exchange even at urea concentrations (5.8 M) well into the denaturation transition zone, and their energetically distinct activated complexes leading to exchange gives evidence of residual structure. Under these experimental conditions, the use of DeltaG(HX) as a basis for HX analysis of RNase A urea denaturation is invalid.     

No effect of trimethylamine N-oxide on the internal dynamics of the protein native fold

Trimethylamine N-oxide (TMAO) is a natural osmolyte accumulated in cells of organisms as they adapt to environmental stresses. In vitro, TMAO increases protein stability and forces partially unfolded structures to refold. Its effects on the native fold are unknown. To investigate the interrelationship between protein stability, internal dynamics and function, the influence of TMAO on the flexibility of the native fold was examined with four different proteins by Trp phosphorescence spectroscopy. Its influence on conformational dynamics was assessed by both the intrinsic phosphorescence lifetime, which reports on the local structure about the triplet probe, and the acrylamide bimolecular quenching rate constant that is a measure of the average acrylamide diffusion coefficient through the macromolecule. The results demonstrate that for apoazurin, alcohol dehydrogenase, alkaline phosphatase and glyceraldehydes-3-phosphate dehydrogenase 1.8 M TMAO does not perturb the flexibility of these macromolecules in a temperature range between - 10 degreesC and up to near the melting temperature. This unexpected finding contrasts with the dampening effect observed with polyols as well as with the expectations based on the preferential exclusion of the osmolyte from the protein surface.     

Covariation of backbone motion throughout a small protein domain

The synchronization (correlation) of conformational fluctuations in folded proteins may influence the rates of enzyme catalysis and ligand binding as well as the stabilities of native proteins and their complexes. However, experimental detection of correlated motions remains difficult. Herein, we present an analysis of the covariation of NMR-derived backbone dynamical parameters among a family of ten mutants of a small protein. Both the spatial restriction and the time scales of backbone motions exhibit a higher degree of covariation than would be expected if the internal motions of each group were independent, providing experimental support for correlated dynamics. Application of this approach to other proteins may reveal dynamical correlations that influence catalysis, ligand-binding and/or protein stability.     

Concordance of residual dipolar couplings, backbone order parameters and crystallographic B-factors for a small alpha/beta protein: a unified picture of high probability, fast atomic motions in proteins

Using ensemble refinement of the third immunoglobulin binding domain (GB3) of streptococcal protein G (a small alpha/beta protein of 56 residues), we demonstrate that backbone (N-H, N-C', Calpha-Halpha, Calpha-C') residual dipolar coupling data in five independent alignment media, generalized order parameters from 15N relaxation data, and B-factors from a high-resolution (1.1A), room temperature crystal structure are entirely consistent with one another within experimental error. The optimal ensemble size representation is between four and eight, as assessed by complete cross-validation of the residual dipolar couplings. Thus, in the case of GB3, all three observables reflect the same low-amplitude anisotropic motions arising from fluctuations in backbone phi/psi torsion angles in the picosecond to nanosecond regime in both solution and crystalline environments, yielding a unified picture of fast, high-probability atomic motions in proteins. An understanding of these motions is crucial for understanding the impact of protein dynamics on protein function, since they provide part of the driving force for triggered conformational changes that occur, for example, upon ligand binding, signal transduction and enzyme catalysis.     

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

Osmolyte effect on the stability and folding of a hyperthermophilic protein

Proteins are known to be stabilized by naturally occurring osmolytes such as amino acids, sugars, and methylamines. Here, we examine the effect of trimethylamine-N-oxide (TMAO) on the conformational stability of ribonuclease HII from a hyperthermophile, Thermococcus kodakaraensis (Tk-RNase HII), which inherently possesses high conformational stability. Heat- and guanidine hydrochloride-induced unfolding experiments demonstrated that the conformational stability of Tk-RNase HII in the presence of 0.5M TMAO was higher than that in the absence of TMAO at all examined temperatures. TMAO affected the unfolding and refolding kinetics of Tk-RNase HII to a similar extent. These results indicate that proteins are universally stabilized by osmolytes, regardless of their robustness, and suggest a stabilization mechanism by osmolytes, caused by the unfavorable interaction of osmolytes with protein backbones in the denatured state. Our results also imply that the basic protein folding principle is not dependent on protein stability and evolution.     

Temperature-mediated switching of protectant-denaturant behavior of trimethylamine-N-oxide and consequences on protein stability from a replica exchange molecular dynamics simulation study

The detailed mechanism of protein folding–unfolding processes with the aid of osmolytes has been a leading topic of discussion over many decades. We have used replica-exchange molecular dynamics simulation to propose the molecular mechanism of interaction of a 20-residue mini-protein with urea and trimethylamine N-oxide (TMAO) that act as denaturing and protecting osmolyte, respectively, in binary osmolyte solutions. Urea is found to exert its action by interacting directly with the protein residues. Temperature tolerance of TMAO’s action is particularly emphasised in this study. At lower range of temperature, TMAO acts as a successful protein protectant. Interestingly, the study discloses the tendency of TMAO molecules to prefer self-association at the protein surface at elevated temperature. A greater number of TMAO molecules in the protein hydration shell at higher temperature is also observed. Dihedral angle principal component analysis and free energy landscape plots sampled all possible conformations adopted by the protein that reveal highly folded behaviour of the protein in pure water and binary TMAO solutions and highly unfolded behaviour in presence of urea.     

Measuring the stability of partly folded proteins using TMAO

Standard methods for measuring free energy of protein unfolding by chemical denaturation require complete folding at low concentrations of denaturant so that a native baseline can be observed. Alternatively, proteins that are completely unfolded in the absence of denaturant can be folded by addition of the osmolyte trimethylamine N-oxide (TMAO), and the unfolding free energy can then be calculated through analysis of the refolding transition. However, neither chemical denaturation nor osmolyte-induced refolding alone is sufficient to yield accurate thermodynamic unfolding parameters for partly folded proteins, because neither method produces both native and denatured baselines in a single transition. Here we combine urea denaturation and TMAO stabilization as a means to bring about baseline-resolved structural transitions in partly folded proteins. For Barnase and the Notch ankyrin domain, which both show two-state equilibrium unfolding, we found that DeltaG degrees for unfolding depends linearly on TMAO concentration, and that the sensitivity of DeltaG degrees to urea (the m-value) is TMAO independent. This second observation confirms that urea and TMAO exert independent effects on stability over the range of cosolvent concentrations required to bring about baseline-resolved structural transitions. Thermodynamic parameters calculated using a global fit that assumes additive, linear dependence of DeltaG degrees on each cosolvent are similar to those obtained by standard urea-induced unfolding in the absence of TMAO. Finally, we demonstrate the applicability of this method to measurement of the free energy of unfolding of a partly folded protein, a fragment of the full-length Notch ankyrin domain.     

A new spin on protein dynamics

Site-directed spin labeling is a general method for investigating structure and conformational switching in soluble and membrane proteins. It will also be an important tool for exploring protein backbone dynamics. A semi-empirical analysis of nitroxide sidechain dynamics in spin-labeled proteins reveals contributions from fluctuations in backbone dihedral angles and rigid-body (collective) motions of alpha helices. Quantitative analysis of sidechain dynamics is sometimes possible, and contributions from backbone modes can be expressed in terms of relative order parameters and rates. Dynamic sequences identified by site-directed spin labeling correlate with functional domains, and so nitroxide scanning could provide an efficient strategy for identifying such domains in high-molecular weight proteins, supramolecular complexes and membrane proteins.     

Non-redox-active small-molecules can accelerate oxidative protein folding by novel mechanisms

Multi-disulfide-bond-containing proteins acquire their native structures through an oxidative folding reaction which involves formation of native disulfide bonds through thiol-disulfide exchange reactions between cysteines and disulfides coupled to a conformational folding event. Oxidative folding rates of the four-disulfide-bond-containing protein bovine pancreatic ribonuclease A (RNase A) in the presence of the synthetic redox-active molecule, (+/-)-trans-1,2-bis(2-mercaptoacetamido)cyclohexane (BMC), and in combination with non-redox-active trimethylamine-N-oxide (TMAO), and trifluorethanol were determined by HPLC analysis. The data indicate that regeneration of RNase A is enhanced 2-fold by BMC (50 microM) and 3-fold upon addition of TMAO (0.2 M) and TFE (3% v/v) relative to control experiments performed in the absence of small-molecules. Examination of the native tendency of the fully-reduced polypeptide and the stability of key folding intermediates suggests that the increased oxidative folding rate can be attributed to native-like elements induced within the fully-reduced polypeptide and the stabilization of native-like species by added non-redox-active molecules.     

Protein dynamics in a family of laboratory evolved thermophilic enzymes

Molecular dynamics simulations were employed to study how protein solution structure and dynamics are affected by adaptation to high temperature. Simulations were carried out on a para-nitrobenzyl esterase (484 residues) and two thermostable variants that were generated by laboratory evolution. Although these variants display much higher melting temperatures than wild-type (up to 18 degrees C higher) they are both >97% identical in sequence to the wild-type. In simulations at 300 K the thermostable variants remain closer to their crystal structures than wild-type. However, they also display increased fluctuations about their time-averaged structures. Additionally, both variants show a small but significant increase in radius of gyration relative to wild-type. The vibrational density of states was calculated for each of the esterases. While the density of states profiles are similar overall, both thermostable mutants show increased populations of the very lowest frequency modes (<10 cm(-1)), with the more stable mutant showing the larger increase. This indicates that the thermally stable variants experience increased concerted motions relative to wild-type. Taken together, these data suggest that adaptation for high temperature stability has resulted in a restriction of large deviations from the native state and a corresponding increase in smaller scale fluctuations about the native state. These fluctuations contribute to entropy and hence to the stability of the native state. The largest changes in localized dynamics occur in surface loops, while other regions, particularly the active site residues, remain essentially unchanged. Several mutations, most notably L313F and H322Y in variant 8G8, are in the region showing the largest increase in fluctuations, suggesting that these mutations confer more flexibility to the loops. As a validation of our simulations, the fluctuations of Trp102 were examined in detail, and compared with Trp102 phosphorescence lifetimes that were previously measured. Consistent with expectations from the theory of phosphorescence, an inverse correlation between out-of-plane fluctuations on the picosecond time scale and phosphorescence lifetime was observed.     

Functional protection of human haemoglobin against protein dissociation

The spectroscopic and functional properties of human adult haemoglobin are clearly disrupted by concentrations of urea > 0.4 M. This disruption of structure and function is completely obviated by the presence of 0.2 M trimethylamine N-oxide (TMAO). Spectroscopic data suggest that TMAO prevents urea-induced production of high-spin haem. Functional analysis shows that TMAO exerts its influence by counteracting urea-induced destabilisation of the T state of the haemoglobin protein. Further studies show, however, that TMAO is not able to exert any such stabilising influences in the presence of high concentrations of typical organic solvent denaturants.     

Global changes in local protein dynamics reduce the entropic cost of carbohydrate binding in the arabinose-binding protein

Protein dynamics make important but poorly understood contributions to molecular recognition phenomena. To address this, we measure changes in fast protein dynamics that accompany the interaction of the arabinose-binding protein (ABP) with its ligand, d-galactose, using NMR relaxation and molecular dynamics simulation. These two approaches present an entirely consistent view of the dynamic changes that occur in the protein backbone upon ligand binding. Increases in the amplitude of motions are observed throughout the protein, with the exception of a few residues in the binding site, which show restriction of dynamics. These counter-intuitive results imply that a localised binding event causes a global increase in the extent of protein dynamics on the pico- to nanosecond timescale. This global dynamic change constitutes a substantial favourable entropic contribution to the free energy of ligand binding. These results suggest that the structure and dynamics of ABP may be adapted to exploit dynamic changes to reduce the entropic costs of binding.     

Preventing misfolding of the prion protein by trimethylamine N-oxide

Transmissible spongiform encephalopathies are a class of fatal neurodegenerative diseases linked to the prion protein. The prion protein normally exists in a soluble, globular state (PrP(C)) that appears to participate in copper metabolism in the central nervous system and/or signal transduction. Infection or disease occurs when an alternatively folded form of the prion protein (PrP(Sc)) converts soluble and predominantly alpha-helical PrP(C) into aggregates rich in beta-structure. The structurally disordered N-terminus adopts beta-structure upon conversion to PrP(Sc) at low pH. Chemical chaperones, such as trimethylamine N-oxide (TMAO), can prevent formation of PrP(Sc) in scrapie-infected mouse neuroblastoma cells [Tatzelt, J., et al. (1996) EMBO J. 15, 6363-6373]. To explore the mechanism of TMAO protection of PrP(C) at the atomic level, molecular dynamics simulations were performed under conditions normally leading to conversion (low pH) with and without 1 M TMAO. In PrP(C) simulations at low pH, the helix content drops and the N-terminus is brought into the small native beta-sheet, yielding a PrP(Sc)-like state. Addition of 1 M TMAO leads to a decreased radius of gyration, a greater number of protein-protein hydrogen bonds, and a greater number of tertiary contacts due to the N-terminus forming an Omega-loop and packing against the structured core of the protein, not due to an increase in the level of extended structure as with the PrP(C) to PrP(Sc) simulation. In simulations beginning with the "PrP(Sc)-like" structure (derived from PrP(C) simulated at low pH in pure water) in 1 M TMAO, similar structural reorganization at the N-terminus occurred, disrupting the extended sheet. The mechanism of protection by TMAO appears to be exclusionary in nature, consistent with previous theoretical and experimental studies. The TMAO-induced N-terminal conformational change prevents residues that are important in the conversion of PrP(C) to PrP(Sc) from assuming extended sheet structure at low pH.     

Accuracy and precision of NMR relaxation experiments and MD simulations for characterizing protein dynamics

Model-free parameters obtained from nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations commonly are used to describe the intramolecular dynamical properties of proteins. To assess the relative accuracy and precision of experimental and simulated model-free parameters, three independent data sets derived from backbone ¹⁵N NMR relaxation experiments and two independent data sets derived from MD simulations of Escherichia coli ribonuclease HI are compared. The widths of the distributions of the differences between the order parameters for pairs of NMR data sets are congruent with the uncertainties derived from statistical analyses of individual data sets; thus, current protocols for analyzing NMR data encapsulate random uncertainties appropriately. Large differences in order parameters for certain residues are attributed to systematic differences between samples for intralaboratory comparisons and unknown, possibly magnetic field-dependent, experimental effects for interlaboratory comparisons. The widths of distributions of the differences between the order parameters for two NMR sets are similar to widths of distributions for an NMR and an MD set or for two MD sets. The linear correlations between the order parameters for an MD set and an NMR set are within the range of correlations observed between pairs of NMR sets. These comparisons suggest that the NMR and MD generalized order parameters for the backbone amide N—H bond vectors are of comparable accuracy for residues exhibiting motions on a fast time scale (<100 ps). Large discrepancies between NMR and MD order parameters for certain residues are attributed to the occurrence of “rare” motional events over the simulation trajectories, the disruption of an element of secondary structure in one of the simulations, and lack of consensus among the experimental data sets. Consequently, (easily detectable) severe distortions of local protein structure and infrequent motional events in MD simulations appear to be the most serious artifacts affecting the accuracy and precision, respectively, of MD order parameters relative to NMR values. In addition, MD order parameters for motions on a fast (<100 ps) timescale are more precisely determined than their NMR counterparts, thereby permitting more detailed dynamic characterization of biologically important residues by MD simulation than is sometimes possible by experimental methods. Proteins 28:481–493, 1997. © 1997 Wiley-Liss, Inc.     

Accurate and efficient description of protein vibrational dynamics: comparing molecular dynamics and Gaussian models

Current all-atom potential based molecular dynamics (MD) allows the identification of a protein's functional motions on a wide-range of timescales, up to few tens of nanoseconds. However, functional, large-scale motions of proteins may occur on a timescale currently not accessible by all-atom potential based MD. To avoid the massive computational effort required by this approach, several simplified schemes have been introduced. One of the most satisfactory is the Gaussian network approach based on the energy expansion in terms of the deviation of the protein backbone from its native configuration. Here, we consider an extension of this model that captures in a more realistic way the distribution of native interactions due to the introduction of effective side-chain centroids. Since their location is entirely determined by the protein backbone, the model is amenable to the same exact and computationally efficient treatment as previous simpler models. The ability of the model to describe the correlated motion of protein residues in thermodynamic equilibrium is established through a series of successful comparisons with an extensive (14 ns) MD simulation based on the AMBER potential of HIV-1 protease in complex with a peptide substrate. Thus, the model presented here emerges as a powerful tool to provide preliminary, fast yet accurate characterizations of protein near-native motion.     

Singular efficacy of trimethylamine N-oxide to counter protein destabilization in ice

This study reports the first quantitative estimate of the thermodynamic stability (Delta G degrees ) of a protein in low-temperature partly frozen aqueous solutions in the presence of the protective osmolytes trimethylamine N-oxide (TMAO), glycine betaine, and sarcosine. The method, based on guanidinium chloride denaturation of the azurin mutant C112S from Pseudomonas aeruginosa, distinguishes between the deleterious effects of subfreezing temperatures from those due specifically to the formation of a solid ice phase. The results point out that in the liquid state molar concentrations of these osmolytes stabilize significantly the native fold and that their effect is maintained on cooling to -15 degrees C. At this temperature, freezing of the solution in the absence of any additive causes a progressive destabilization of the protein, Delta G degrees decreasing up to 3-4 kcal/mol as the fraction of liquid water in equilibrium with ice ( V L) is reduced to less than 1%. The ability of the three osmolytes to prevent the decrease in protein stability at small V L varies significantly among them, ranging from the complete inertness of sarcosine to full protection by TMAO. The singular effectiveness of TMAO among the osmolytes tested until now is maintained high even at concentrations as low as 0.1 M when the additive stabilization of the protein in the liquid state is negligible. In all cases the reduction in Delta G degrees caused by the solidification of water correlates with the decrease in m-value entailing that protein-ice interactions generally conduct to partial unfolding of the native state. It is proposed that the remarkable effectiveness of TMAO to counter the ice perturbation is owed to binding of the osmolyte to ice, thereby inhibiting protein adsorption to the solid phase.     

Relationships between protein structure and dynamics from a database of NMR-derived backbone order parameters

The amplitude of protein backbone NH group motions on a time-scale faster than molecular tumbling may be determined by analysis of (15)N NMR relaxation data according to the Lipari-Szabo model free formalism. An internet-accessible database has been compiled containing 1855 order parameters from 20 independent NMR relaxation studies on proteins whose three-dimensional structures are known. A series of statistical analyses has been performed to identify relationships between the structural features and backbone dynamics of these proteins. Comparison of average order parameters for different amino acid types indicates that amino acids with small side-chains tend to have greater backbone flexibility than those with large side-chains. In addition, the motions of a given NH group are also related to the sizes of the neighboring amino acids in the primary sequence. The secondary structural environment appears to influence backbone dynamics relatively weakly, with only subtle differences between the order parameter distributions of loop structures and regular hydrogen bonded secondary structure elements. However, NH groups near helix termini are more mobile on average than those in the central regions of helices. Tertiary structure influences are also relatively weak but in the expected direction, with more exposed residues being more flexible on average than residues that are relatively inaccessible to solvent.     

Temperature-dependent dynamics of the villin headpiece helical subdomain,an unusually small thermostable protein

(15)N spin relaxation experiments were used to measure the temperature-dependence of protein backbone conformational fluctuations in the thermostable helical subdomain, HP36, of the F-actin-binding headpiece domain of chicken villin. HP36 is the smallest domain of a naturally occurring protein that folds cooperatively to a compact native state. Spin-lattice, spin-spin, and heteronuclear nuclear Overhauser effect relaxation data for backbone amide (15)N spins were collected at five temperatures in the range of 275-305 K. The data were analyzed using a model-free formalism to determine generalized order parameters, S, that describe the distribution of N-H bond vector orientations in a molecular reference frame. A novel parameter, Lambda=dln(1-S)/dln T is introduced to characterize the temperature-dependence of S. An average value of Lambda=4.5 is obtained for residues in helical conformations in HP36. This value of Lambda is not reproduced by model potential energy functions commonly used to parameterize S. The maximum entropy principle was used to derive a new model potential function that reproduces both S and Lambda. Contributions to the entropy, S(r), and heat capacity, C(r)(p), from reorientational conformational fluctuations were analyzed using this potential energy function. Values of S(r) show a qualitative dependence on S similar to that obtained for the diffusion-in-a-cone model; however, quantitative differences of up to 0.5k, in which k is the Boltzmann constant, are observed. Values of C(r)(p) approach zero for small values of S and approach k for large values of S; the largest values of C(r)(p) are predicted to occur for intermediate values of S. The results suggest that backbone dynamics, as probed by relaxation measurements, make very little contribution to the heat capacity difference between folded and unfolded states for HP36.     

Dynamic transition associated with the thermal denaturation of a small Beta protein

We studied the temperature dependence of the picosecond internal dynamics of an all-beta protein, neocarzinostatin, by incoherent quasielastic neutron scattering. Measurements were made between 20 degrees C and 71 degrees C in heavy water solution. At 20 degrees C, only 33% of the nonexchanged hydrogen atoms show detectable dynamics, a number very close to the fraction of protons involved in the side chains of random coil structures, therefore suggesting a rigid structure in which the only detectable diffusive movements are those involving the side chains of random coil structures. At 61.8 degrees C, although the protein structure is still native, slight dynamic changes are detected that could reflect enhanced backbone and beta-sheet side-chain motions at this higher temperature. Conversely, all internal dynamics parameters (amplitude of diffusive motions, fraction of immobile scatterers, mean-squared vibration amplitude) rapidly change during heat-induced unfolding, indicating a major loss of rigidity of the beta-sandwich structure. The number of protons with diffusive motion increases markedly, whereas the volume occupied by the diffusive motion of protons is reduced. At the half-transition temperature (T = 71 degrees C) most of backbone and beta-sheet side-chain hydrogen atoms are involved in picosecond dynamics.