Amide backbone and water-related H/D isotope effects on the dynamics of a protein folding reaction

Probing the role of protein conformational changes in the mechanism of prenylated-FMN-dependent phenazine-1-carboxylic acid decarboxylase

Phenazine-1-carboxylic acid decarboxylase (PhdA) is a prenylated-FMN-dependent (prFMN) enzyme belonging to the UbiD family of decarboxylases. Many UbiD-like enzymes catalyze (de)carboxylation reactions on aromatic rings and conjugated double bonds and are potentially valuable industrial catalysts. We have investigated the mechanism of PhdA using a slow turnover substrate, 2,3-dimethylquinoxaline-5-carboxylic acid (DQCA). Detailed analysis of the pH dependence and solvent deuterium isotope effects associated with the reaction uncovered unusual kinetic behavior. At low substrate concentrations, a substantial inverse solvent isotope effect (SIE) is observed on Vmax/KM of ∼ 0.5 when reaction rates of DQCA in H2O and D2O are compared. Under the same conditions, a normal SIE of 4.15 is measured by internal competition for proton transfer to the product. These apparently contradictory results indicate that the SIE values report on different steps in the mechanism. A proton inventory analysis of the reaction under Vmax/KM and Vmax conditions points to a "medium effect" as the source of the inverse SIE. Molecular dynamics simulations of the effect of D2O on PhdA structure support that D2O reduces the conformational lability of the enzyme and results in a more compact structure, akin to the active, "closed" conformer observed in crystal structures of some UbiD-like enzymes. Consistent with the simulations, PhdA was found to be more stable in D2O and to bind DQCA more tightly, leading to the observed rate enhancement under Vmax/KM conditions.

Effects of Hydrogen/Deuterium Exchange on Protein Stability in Solution and in the Gas Phase

Mass spectrometry (MS)-based techniques are widely used for probing protein structure and dynamics in solution. H/D exchange (HDX)-MS is one of the most common approaches in this context. HDX is often considered to be a "benign" labeling method, in that it does not perturb protein behavior in solution. However, several studies have reported that D₂O pushes unfolding equilibria toward the native state. The origin, and even the existence of this protein stabilization remain controversial. Here we conducted thermal unfolding assays in solution to confirm that deuterated proteins in D₂O are more stable, with 2-4 K higher melting temperatures than unlabeled proteins in H₂O. Previous studies tentatively attributed this phenomenon to strengthened H-bonds after deuteration, an effect that may arise from the lower zero-point vibrational energy of the deuterated species. Specifically, it was proposed that strengthened water-water bonds (W···W) in D₂O lower the solubility of nonpolar side chains. The current work takes a broader view by noting that protein stability in solution also depends on water-protein (W···P) and protein-protein (P···P) H-bonds. To help unravel these contributions, we performed collision-induced unfolding (CIU) experiments on gaseous proteins generated by native electrospray ionization. CIU profiles of deuterated and unlabeled proteins were indistinguishable, implying that P···P contacts are insensitive to deuteration. Thus, protein stabilization in D₂O is attributable to solvent effects, rather than alterations of intraprotein H-bonds. Strengthening of W···W contacts represents one possible explanation, but the stabilizing effect of D₂O can also originate from weakened W···P bonds. Future work will be required to elucidate which of these two scenarios is correct, or if both contribute to protein stabilization in D₂O. In any case, the often-repeated adage that "D-bonds are more stable than H-bonds" does not apply to intramolecular contacts in native proteins.

A methyl-TROSY based 13C relaxation dispersion NMR experiment for studies of chemical exchange in proteins

A methyl Transverse Relaxation Optimized Spectroscopy (methyl-TROSY) based, multiple quantum (MQ) ¹³C Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiment is described. The experiment is derived from the previously developed MQ ¹³C-¹H CPMG scheme (Korzhnev in J Am Chem Soc 126: 3964-73, 2004) supplemented with a CPMG train of refocusing ¹H pulses applied with constant frequency and synchronized with the ¹³C CPMG pulse train. The optimal ¹H 'decoupling' scheme that minimizes the amount of fast-relaxing methyl MQ magnetization present during CPMG intervals, makes use of an XY-4 phase cycling of the refocusing composite ¹H pulses. For small-to-medium sized proteins, the MQ ¹³C CPMG experiment has the advantage over its single quantum (SQ) ¹³C counterpart of significantly reducing intrinsic, exchange-free relaxation rates of methyl coherences. For high molecular weight proteins, the MQ ¹³C CPMG experiment eliminates complications in the interpretation of MQ ¹³C-¹H CPMG relaxation dispersion profiles arising from contributions to exchange from differences in methyl ¹H chemical shifts between ground and excited states. The MQ ¹³C CPMG experiment is tested on two protein systems: (1) a triple mutant of the Fyn SH3 domain that interconverts slowly on the chemical shift time scale between the major folded state and an excited state folding intermediate; and (2) the 82-kDa enzyme Malate Synthase G (MSG), where chemical exchange at individual Ile δ1 methyl positions occurs on a much faster time-scale.

The Influence of Deuterium Isotope Effects on Structural Rearrangements, Ensemble Equilibria, and Hydrogen Bonding in Protic Ionic Liquids

We report strong isotope effects for the protic ionic liquid triethylammonium methanesulfonate [TEA][OMs] by means of deuterium solid-state NMR spectroscopy covering broad temperature ranges from 65 K to 313 K. Both isotopically labelled PILs differ in non-deuterated and fully deuterated ethyl groups of the triethyl ammonium cations. The N-D bond of both cations is used as sensitive probe for hydrogen bonding and structural ordering. The ² H NMR line shape analysis provides the deuteron quadrupole coupling constants and the characteristics of a broad heterogeneous phase with simultaneously present static and mobile states indicating plastic crystal behavior. The temperatures where both states are equally populated differ by about 80 K for the two PILs, showing that deuteration of the ethyl groups in the trialkylammonium cations tremendously shifts the equilibrium towards the static state. In addition, it leads to a significant less cooperative transition, associated with a significantly reduced standard molar transition entropy.

Bulk phase behaviour vs interface adsorption: Effects of anions and isotopes on β-lactoglobulin (BLG) interactions

HYPOTHESIS

Protein adsorption is highly relevant in numerous applications ranging from food processing to medical implants. In this context, it is important to gain a deeper understanding of protein-protein and protein-surface interactions. Thus, the focus of this investigation is on the interplay of bulk properties and surface properties on protein adsorption. It was hypothesised that the type of solvent and ions in solution should significantly influence the protein's bulk and interface behaviour, which has been observed in literature and previous work for other net negatively charged, globular proteins such as bovine serum albumin (BSA).

EXPERIMENTS

The phase behaviour of β-lactoglobulin (BLG) with lanthanum chloride (LaCl₃) and iodide (LaI₃) in normal water H₂O(l) and heavy water (D₂O(l)) was established via optical microscopy and ultraviolet-visible spectroscopy. The formation of an adsorption layer and its properties such as thickness, density, structure, and hydration was investigated via neutron reflectivity, quartz-crystal microbalance with dissipation, and infra-red measurements.

FINDINGS

β-lactoglobulin does not show significant anion-induced or isotope-induced effects - neither in bulk nor at the solid-liquid interface, which deviates strongly from the behaviour of bovine serum albumin. We also provide a comprehensive discussion and comparison of protein-specific bulk and interface behaviour between bovine serum albumin and β-lactoglobulin dependent on anion, cation, solvent, and substrate properties. These findings pave the way for understanding the transition from adsorption to crystallisation.

Low amounts of heavy water increase the phase separation propensity of a fragment of the androgen receptor activation domain

The phase equilibria of intrinsically disordered proteins are exquisitely sensitive to changes in solution conditions and this can be used to investigate the driving forces of phase separation in vitro as well as the biological roles of phase transitions in live cells. Here we investigate how using D₂ O as co-solvent in an aqueous buffer changes the phase equilibrium of a fragment of the activation domain of the androgen receptor, a transcription factor that plays a role in the development of the male phenotype and is a therapeutic target for castration resistant prostate cancer. We show how replacing even small fractions of H₂ O with D₂ O increases the propensity of this fragment to undergo liquid-liquid phase separation, likely reflecting a stabilization of the hydrophobic interactions that drive condensation. Our results indicate that it is necessary to take this effect into consideration when studying phase separation phenomena with biophysical methods that require using D₂ O as a co-solvent. In addition, they suggest that additions of D₂ O may be used to enhance phase separation phenomena in cells, facilitating their observation.

Nano-second protein dynamics of key residue at Position 38 in catechol-O-methyltransferase system: a time-resolved fluorescence study

Human catechol-O-methyltransferase, a key enzyme related to neurotransmitter metabolism, catalyses a methyl transfer from S-adenosylmethionine to catechol. Although extensive studies aim to understand the enzyme mechanisms, the connection of protein dynamics and enzyme catalysis is still not clear. Here, W38in (Trp143Phe) and W38in/Y68A (Trp143Phe with Tyr68Ala) mutants were carried out to study the relationship of dynamics and catalysis in nano-second timescale using time-resolved fluorescence lifetimes and Stokes shifts in various solvents. The comprehensive data implied the mutant W38in/Y68A with lower activity is more rigid than the ‘WT’−W38in, suggesting the importance of flexibility at residue 38 to maintain the optimal catalysis.

X-ray crystallographic studies on the hydrogen isotope effects of green fluorescent protein at sub-ångström resolutions

Hydrogen atoms are critical to the nature and properties of proteins, and thus deuteration has the potential to influence protein function. In fact, it has been reported that some deuterated proteins show different physical and chemical properties to their protiated counterparts. Consequently, it is important to investigate protonation states around the active site when using deuterated proteins. Here, hydrogen isotope effects on the S65T/F99S/M153T/V163A variant of green fluorescent protein (GFP), in which the deprotonated B form is dominant at pH 8.5, were investigated. The pH/pD dependence of the absorption and fluorescence spectra indicates that the protonation state of the chromophore is the same in protiated GFP in H₂O and protiated GFP in D₂O at pH/pD 8.5, while the pK_a of the chromophore became higher in D₂O. Indeed, X-ray crystallographic analyses at sub-ångström resolution revealed no apparent changes in the protonation state of the chromophore between the two samples. However, detailed comparisons of the hydrogen OMIT maps revealed that the protonation state of His148 in the vicinity of the chromophore differed between the two samples. This indicates that protonation states around the active site should be carefully adjusted to be the same as those of the protiated protein when neutron crystallographic analyses of proteins are performed.

Opposing Electronic and Nuclear Quantum Effects on Hydrogen Bonds in H2 O and D2 O

The effect of extending the O-H bond length(s) in water on the hydrogen-bonding strength has been investigated using static ab initio molecular orbital calculations. The "polar flattening" effect that causes a slight σ-hole to form on hydrogen atoms is strengthened when the bond is stretched, so that the σ-hole becomes more positive and hydrogen bonding stronger. In opposition to this electronic effect, path-integral ab initio molecular-dynamics simulations show that the nuclear quantum effect weakens the hydrogen bond in the water dimer. Thus, static electronic effects strengthen the hydrogen bond in H2 O relative to D2 O, whereas nuclear quantum effects weaken it. These quantum fluctuations are stronger for the water dimer than in bulk water.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Hydrogen isotope replacement changes hydration and large scale structure, but not small scale structure, of agarose hydrogel networks

Agarose samples of low (Ag1) and high (Ag2) O -methyl content on position 6 of the galactose residue were studied in H₂O and D₂O. Differential scanning calorimetry, turbidity and rheological measurements showed a [Formula: see text] 2 ^° C shift in the coil-to-helix transition temperature, indicating higher helix stability in D₂O. The differential scanning calorimetry data could be superimposed using a temperature shift factor, suggesting similar extents of helix aggregation in both solvents. Small angle X-ray scattering of H₂O and D₂O gels were essentially identical, indicating no change in the small scale ( [Formula: see text] 0.05-20 nm) network structure on isotopic exchange. Larger ([Formula: see text] 1 μm) scale heterogeneities were more pronounced in deuterium gels. The ¹HT₂ relaxation times were measured at different H/D ratios. These relaxation times were analyzed using a model assuming regular solution mixing of H₂O, HDO and D₂O between the solvent and gel phases. The fit results suggested that H₂O has higher affinity for the agarose network than HDO and D₂O. The difference, however, was much larger for the Ag2 sample. This finding implies that the higher hydrophobic effect observed in D₂O affects the hydration state much more strongly for the more hydrophobic (and more polarizable) agarose sample Ag2. As a consequence, Ag2 (but not Ag1) gels retained more H₂O than D₂O. In contrast, the bulk rheology of either hydrogel was not affected by the isotopic exchange.

Enthalpic stabilization of an SH3 domain by D2 O

The stability of a protein is vital for its biological function, and proper folding is partially driven by intermolecular interactions between protein and water. In many studies, H2 O is replaced by D2 O because H2 O interferes with the protein signal. Even this small perturbation, however, affects protein stability. Studies in isotopic waters also might provide insight into the role of solvation and hydrogen bonding in protein folding. Here, we report a complete thermodynamic analysis of the reversible, two-state, thermal unfolding of the metastable, 7-kDa N-terminal src-homology 3 domain of the Drosophila signal transduction protein drk in H2 O and D2 O using one-dimensional 19 F NMR spectroscopy. The stabilizing effect of D2 O compared with H2 O is enthalpic and has a small to insignificant effect on the temperature of maximum stability, the entropy, and the heat capacity of unfolding. We also provide a concise summary of the literature about the effects of D2 O on protein stability and integrate our results into this body of data.

Site-Specific Hydrogen Exchange in a Membrane Environment Analyzed by Infrared Spectroscopy

Hydrogen exchange is a powerful method to examine macromolecules. In membrane proteins, exchange can distinguish between solvent-accessible and -inaccessible residues due to shielding by the hydrophobic environment of the lipid bilayer. Herein, rather than examining which residues undergo hydrogen exchange, we employ a protocol that enables the full deuteration of all polar hydrogens in a membrane protein. We then measure the impact of hydrogen exchange on the shift of the amide I vibrational mode of individually labeled sites. The results enable us to correlate polarity with vibrational shifts, thereby providing a powerful tool to examine specific locations within a membrane protein in its native membrane environment.

Equilibrium folding dynamics of meACP in water, heavy water, and low concentration of urea

Many proteins fold in apparent two-state behavior, as partially folded intermediates only transiently accumulate and easily escape detection. Besides a native form and a mainly unfolded form, we captured a partially unfolded form of an acyl carrier protein from Micromonospora echinospora (meACP) in the folding/unfolding equilibrium using chemical exchange saturation transfer NMR experiments. The C-terminal region of the partially unfolded form is mainly folded and the N-terminal is unfolded. Furthermore, to understand how the folding process of meACP is influenced by solvent environments, we compared the folding dynamics of meACP in D₂O, H₂O and low concentration of urea. As the environment becomes more denaturing from D₂O to H₂O and then to urea, the unfolded state becomes increasingly populated, and the folding rate decreases. Adding a small amount of urea, which does not change solvent viscosity, has little effects on the unfolding rates, while changing H₂O to D₂O reduces the unfolding rates possibly due to the increase of solvent viscosity. The quantified solvent effects on the protein folding Gibbs energy and activation energy suggest that the transition state of folding may have a similar structure to the native state of the protein.

Extended Law of Corresponding States Applied to Solvent Isotope Effect on a Globular Protein

Investigating proteins with techniques such as NMR or neutron scattering frequently requires the partial or complete substitution of D2O for H2O as a solvent, often tacitly assuming that such a solvent substitution does not significantly alter the properties of the protein. Here, we report a systematic investigation of the solvent isotope effect on the phase diagram of the lens protein γB-crystallin in aqueous solution as a model system exhibiting liquid-liquid phase separation. We demonstrate that the observed strong variation of the critical temperature Tc can be described by the extended law of corresponding states for all H2O/D2O ratios, where scaling of the temperature by Tc or the reduced second virial coefficient accurately reproduces the binodal, spinodal, and osmotic compressibility. These findings highlight the impact of H2O/D2O substitution on γB-crystallin properties and warrant further investigations into the universality of this phenomenon and its underlying mechanisms.

Heavy water: a simple solution to increasing the brightness of fluorescent proteins in super-resolution imaging

Photoactivatable fluorescent proteins (PA-FPs) are widely used in live single-molecule super-resolution imaging but emit substantially fewer photons than organic dyes do. Herein, we show that in heavy water (D2O) instead of H2O, common PA-FPs emit 26-54% more photons, effectively improving the localization precision in super-resolution imaging.

Rational design of an efficient, genetically encodable, protein-encased singlet oxygen photosensitizer

Singlet oxygen, O(2)(a(1)Δ(g)), plays a key role in many processes of cell signaling. Limitations in mechanistic studies of such processes are generally associated with the difficulty of controlling the amount and location of O(2)(a(1)Δ(g)) production in or on a cell. As such, there is great need for a system that (a) selectively produces O(2)(a(1)Δ(g)) in appreciable and accurately quantifiable yields and (b) can be localized in a specific place at the suborganelle level. A genetically encodable, protein-encased photosensitizer is one way to achieve this goal. Through a systematic and rational approach involving mutations to a LOV2 protein that binds the chromophore flavin mononucleotide (FMN), we have developed a promising photosensitizer that overcomes many of the problems that affect related systems currently in use. Specifically, by decreasing the extent of hydrogen bonding between FMN and a specific amino acid residue in the local protein environment, we decrease the susceptibility of FMN to undesired photoinitiated electron-transfer reactions that kinetically compete with O(2)(a(1)Δ(g)) production. As a consequence, our protein-encased FMN system produces O(2)(a(1)Δ(g)) with the uniquely large quantum efficiency of 0.25 ± 0.03. We have also quantified other key photophysical parameters that characterize this sensitizer system, including unprecedented H(2)O/D(2)O solvent isotope effects on the O(2)(a(1)Δ(g)) formation kinetics and yields. As such, our results facilitate future systematic developments in this field.

A peptide's perspective of water dynamics

This Perspective is focused on amide groups of peptides interacting with water. The 2D IR spectroscopy has already enabled structural aspects of the peptide backbone to be determined through its ability to measure the coupling between different amide-I modes. Here we describe why nonlinear IR is emerging as the method of choice to examine the fast components of the water dynamics near peptides and how isotopically edited peptide links can be used to probe the local water at a residue level in proteins. This type of research necessarily involves an intimate mix of theory and experiment. The description of the results is underpinned by relatively well established quantum-statistical theories that describe the important manifestations of peptide vibrational frequency fluctuations.

Hydrogen bond dynamics in heavy water studied with quantum dynamical simulations

The structure and dynamics of the hydrogen-bond network in heavy water (D(2)O) is studied as a function of the temperature using quantum dynamical simulations. Our approach combines an ab initio-based representation of the water interactions with an explicit quantum treatment of the molecular motion. A direct connection between the calculated linear and nonlinear vibrational spectra and the underlying molecular dynamics is made, which provides new insights into the rearrangement of the hydrogen-bond network in heavy water. A comparison with previous calculations on liquid H(2)O suggests that tunneling does not effectively contribute to the dynamics of the water hydrogen-bond network above the melting point. However, the effects of nuclear quantization are not negligible at all temperatures and become increasingly important near the melting point, which is in agreement with recent experimental analysis of the structural properties of liquid water as well as of the proton momentum distribution in supercooled water.

Force spectroscopy of an elastic peptide: effect of D₂O and temperature on persistence length

This study explores the mechanical unfolding of elastic protein analogues as a function of temperature, in both H₂O and D₂O, using atomic force microscopy (AFM) force spectroscopy in a specially constructed AFM liquid cell. This represents the first time that the effect of D₂O on protein flexibility has been investigated at the single molecule level by this technique. Model elastic peptides, R6, were encoded from synthetic genes expressed in Escherichia coli. The peptides possess short N- and C-terminal domains with central repetitive domains containing 13 repeats of the motif PGQGQQ-plus-GYYPTSLQQ. These sequences mimic those in native high molecular weight subunit glutenin proteins which confer elasticity to bread dough. Fitting single molecule stretching events to the worm-like chain model, allows determination of the molecular persistence length under various experimental conditions. The effect of increasing the temperature is to increase the persistence length of the peptide in both H₂O and D₂O, consistent with the expected "thermal softening" effect. However, the effect is significantly enhanced in D₂O, in which the persistence length at 45°C is ∼25% greater than the value measured in H₂O at the same temperature. Stronger intrapeptide H-bonding due to isotopic substitution of hydrogen for deuterium is the most likely cause of the enhanced backbone rigidity.

Single-molecule spectroscopy of protein folding in a chaperonin cage

Molecular chaperones are known to be essential for avoiding protein aggregation in vivo, but it is still unclear how they affect protein folding mechanisms. We use single-molecule Förster resonance energy transfer to follow the folding of a protein inside the GroEL/GroES chaperonin cavity over a time range from milliseconds to hours. Our results show that confinement in the chaperonin decelerates the folding of the C-terminal domain in the substrate protein rhodanese, but leaves the folding rate of the N-terminal domain unaffected. Microfluidic mixing experiments indicate that strong interactions of the substrate with the cavity walls impede the folding process, but the folding hierarchy is preserved. Our results imply that no universal chaperonin mechanism exists. Rather, a competition between intra- and intermolecular interactions determines the folding rates and mechanisms of a substrate inside the GroEL/GroES cage.

The properties of water: insights from quantum simulations

The properties of water play a central role in many phenomena of relevance to different areas of science, including physics, chemistry, biology, geology, and climate research. Although well studied for decades, the behavior of water under different conditions and in different environments still remains mysterious and often surprising. In this article, various efforts aimed at providing a comprehensive representation of the water properties at a molecular level through computer modeling and simulation will be described. In particular, the unique role played by the hydrogen-bond network will be examined, first in liquid water, then in the solvation of model biological compounds, and finally in ice, especially highlighting the important effects related to the quantization of the nuclear motion.

A single-molecule perspective on the role of solvent hydrogen bonds in protein folding and chemical reactions

We present an array of force spectroscopy experiments that aim to identify the role of solvent hydrogen bonds in protein folding and chemical reactions at the single-molecule level. In our experiments we control the strength of hydrogen bonds in the solvent environment by substituting water (H(2)O) with deuterium oxide (D(2)O). Using a combination of force protocols, we demonstrate that protein unfolding, protein collapse, protein folding and a chemical reaction are affected in different ways by substituting H(2)O with D(2)O. We find that D(2)O molecules form an integral part of the unfolding transition structure of the immunoglobulin module of human cardiac titin, I27. Strikingly, we find that D(2)O is a worse solvent than H(2)O for the protein I27, in direct contrast with the behaviour of simple hydrocarbons. We measure the effect of substituting H(2)O with D(2)O on the force dependent rate of reduction of a disulphide bond engineered within a single protein. Altogether, these experiments provide new information on the nature of the underlying interactions in protein folding and chemical reactions and demonstrate the power of single-molecule techniques to identify the changes induced by a small change in hydrogen bond strength.

A molecular dynamics simulation study of the solvent isotope effect on copper plastocyanin

The effect of heavy water on the structure and dynamics of copper plastocyanin as well as on some aspects of the solvent dynamics at the protein-solvent interfacial region have been investigated by molecular dynamics simulation. The simulated system has been analyzed in terms of the atomic root mean square deviation and fluctuations, intraprotein H-bond pattern, dynamical cross-correlation map and the results have been compared with those previously obtained for plastocyanin in H2O (Ciocchetti et al. Biophys. Chem. 69 (1997), 185-198). The simulated plastocyanin structure in the two solvents, averaging 1 ns, is very similar along the beta-structure regions, while the most significant differences are registered, analogous to the turns and the regions likely involved in the electron transfer pathway. Moreover, plastocyanin in D2O shows an increase in the number of both the intraprotein H-bonds and the residues involved in correlated motions. An analysis of the protein-solvent coupling evidenced that D2O makes the H-bond formation more difficult with the solvent molecules for positively charged and polar residues, while an opposite trend is observed for negatively charged residues. On the other hand, the frequency of exchange of the solvent molecules involved in the protein-solvent H-bond formation is significantly depressed in D2O. The results are discussed also in connection with protein functionality and briefly with some experimental results connected with the thermostability of proteins in D2O.

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.

Low Energy Pathways and Non-native Interactions

Four versions of a beta-sheet protein (CD2.d1) have been made, each with a single artificial disulfide bond inserted into hairpin structures. Folding kinetics of reduced and oxidized forms shows bridge position strongly influences its effect on the folding reaction. Bridging residues 58 and 62 does not affect the rapidly formed intermediate (I) or rate-limiting transition (t) state, whereas bridging 33 and 38, or 31 and 41, lowers the t-state energy, with the latter having the stronger influence. Bridging residues 79 and 90 stabilizes both I- and t-states. To assess additivity in the energetic effects of these bridges, four double-bridge variants have also been made. All show precise additivity of overall stability, with two showing additivity when ground states and the rate-limiting t-state are assessed, i.e. no measurable change in the folding mechanism occurs. However, combining 31-41 and 79-90 bridges produces a molecule that folds through a different pathway, with a much more stable intermediate than expected and a much higher t-state barrier. This is explained by the artificial introduction of stabilizing, non-native contacts in the I-state. More surprisingly, for another double-bridge version (58-62 and 79-90) both I- and t-states are less stable than expected, showing that conformational constraints introduced by the two bridges prevent formation of non-native contacts that would otherwise stabilize the I- and t-states, thereby lowering the energy of the folding landscape in the wild-type (unbridged) molecule. We conclude that the lowest energy path for folding has I- and t-state structures that are stabilized by non-native interactions.

Isotopic Quantum Effects in Liquid Methanol

Density functional calculations (B3 LYP/6-31+G*) on molecular clusters and a quantum cluster equilibrium (QCE) model were used to calculate thermodynamic and structural properties of four isotopically labeled methanol species. The method allowed the reproduction of the characteristic differences in boiling points and heats of vaporization. Structural changes were also detected and related to recent experimental findings. It was shown that isotopic effects clearly have a quantum-mechanical origin.

Beyond the EX1 limit: probing the structure of high-energy states in protein unfolding

Hydrogen exchange kinetics in native solvent conditions have been used to explore the conformational fluctuations of an immunoglobulin domain (CD2.domain1). The global folding/unfolding kinetics of the protein are unaltered between pH 4.5 and pH 9.5, allowing us to use the pH-dependence of amide hydrogen/deuterium exchange to characterise conformational states with energies up to 7.2kcal/mol higher than the folded ground state. The study was intended to search for discreet unfolding intermediates in this region of the energy spectrum, their presence being revealed by the concerted exchange behaviour of subsets of amide groups that become accessible at a given free energy, i.e. the spectrum would contain discreet groupings. Protection factors for 58 amide groups were measured across the pH range and the hydrogen-exchange energy profile is described. More interestingly, exchange behaviour could be grouped into three categories; the first two unremarkable, the third unexpected. (1) In 33 cases, amide exchange was dominated by rapid fluctuation, i.e. the free energy difference between the ground state and the rapidly accessed open state is sufficiently low that the contribution from crossing the unfolding barrier is negligible. (2) In 18 cases exchange is dominated by the global folding transition barrier across the whole pH range measured. The relationship between hydroxyl ion concentration and observed exchange rate is hyperbolic, with the limiting rate being that for global unfolding; the so-called EX1 limit. For these, the free energy difference between the folded ground state and any rapidly-accessed open state is too great for the proton to be exchanged through such fluctuations, even at the highest pH employed in this study. (3) For the third group, comprising five cases, we observe a behaviour that has not been described. In this group, as in category 2, the rate of exchange reaches a plateau; the EX1 limit. However, as the intrinsic exchange rate (k(int)) is increased, this limit is breached and the rate begins to rise again. This unintuitive behaviour does not result from pH instability, rather it is a consequence of amide groups experiencing two processes; rapid fluctuation of structure and crossing the global barrier for unfolding. The boundary at which the EX1 limit is overcome is determined by the equilibrium distribution of the fluctuating open and closed states (K(O/C)) and the rate constant for unfolding (k(u)). This critical boundary is reached when k(int)K(O/C)=k(u). Given that, in a simple transition state formalism: k(u)=K(#)k' (where K(#) describes the equilibrium distribution between the transition and ground state and k' describes the rate of a barrierless rearrangement), it follows that if the pH is raised to a level where k(int)=k', then the entire free energy spectrum from ground state to transition state could be sampled.

Effect of heavy water on protein flexibility

The effects of heavy water (D(2)O) on internal dynamics of proteins were assessed by both the intrinsic phosphorescence lifetime of deeply buried Trp residues, which reports on the local structure about the triplet probe, and the bimolecular acrylamide phosphorescence quenching rate constant that is a measure of the average acrylamide diffusion coefficient through the macromolecule. The results obtained with several protein systems (ribonuclease T1, superoxide dismutase, beta-lactoglobulin, liver alcohol dehydrogenase, alkaline phosphatase, and apo- and Cd-azurin) demonstrate that in most cases D(2)O does significantly increase the rigidity the native structure. With the exception of alkaline phosphatase, the kinetics of the structure tightening effect of deuteration are rapid compared with the rate of H/D exchange of internal protons, which would then assign the dampening of structural fluctuations in D(2)O to a solvent effect, rather than to stronger intramolecular D bonding. Structure tightening by heavy water is generally amplified at higher temperatures, supporting a mostly hydrophobic nature of the underlying interaction, and under conditions that destabilize the globular fold.

A kinetic folding intermediate probed by native state hydrogen exchange

Stopped-flow fluorescence studies on the N-terminal domain of rat CD2 (CD2.d1) have demonstrated that folding from the fully denatured state (U) proceeds via the transient accumulation of an apparent intermediate (I) in a so-called burst phase that precedes the rate-limiting transition leading to the native state (N). A previous pH-dependent equilibrium hydrogen exchange (HX) study identified a subset of amides in CD2.d1 which, under EX2 conditions, exchange from N with free energies greater than or equal to the free energy difference between the N and I states calculated from the stopped-flow data. Under EX1 conditions the rates of HX for these amides tend towards an asymptote that matches the global unfolding rate calculated from the stopped-flow data, suggesting that exchange for these amides requires traversing the N-to-I transition state barrier. Exchange for these amides presumably occurs from exchange-competent forms comprising the kinetic burst phase therefore. To explore this idea further, native state HX (NHX) data have been collected for CD2.d1 under EX2 conditions using denaturant concentrations which span either side of the denaturant concentration where, according to the stopped-flow data, the apparent U and I states are iso-energetic. The data fit to a two-component, sub-global (sg)/global (g) NHX mechanism, yielding Delta G and m value parameters (where the m value is a measure of hydrocarbon solvation). Regression analysis demonstrates that the (m(sg), Delta G(sg)) and (m(g), Delta G(g)) values calculated for this subset of amides correspond with those describing the kinetic burst phase transition. This result confirms the ability of the NHX technique to explore the structural and energetic properties of kinetic folding intermediates.

Progressive stabilization of intermediate and transition states in protein folding reactions by introducing surface hydrophobic residues

It can be argued from the principle of solvent exclusion that the introduction of hydrophobic residues onto the surface of a protein will not destabilize the folded state because the nonpolar side chain will be at least as exposed in the unfolded state as it is when the protein chain is folded. A comparison of the folding pathway of wild type and 11 site-directed mutants of CD2.d1 shows this to be true. In fact, owing to partial burial of nonpolar groups as folding proceeds, we find that the rapidly formed intermediate state and, to a greater extent, the transition state are generally stabilized by hydrophobic surface mutations. This effect is slightly moderated in the folded state presumably by the perturbation of van der Waals' contacts and/or local electrostatic interactions that have a greater influence in this fully compact structure. The fact that in all but one case we find that stabilization of the rapidly collapsed intermediate is accompanied by a faster acquisition of the folded state refutes the argument that I states are generally "off pathway" conformations or ensembles that lead to the inhibition of otherwise more rapid folding trajectories.

A statistical appraisal of native state hydrogen exchange data: evidence for a burst phase continuum?

For a number of proteins, folding occurs via the rapid accumulation of secondary and tertiary structural features in a so-called burst phase, preceding the relatively slow, highly activated transition leading to the native state. A fundamental question is: do these burst phase reactions comprise two phase-separated thermodynamic states or a continuum of states? Ribonuclease HI (RNase H) from Escherichia coli and phage T4 lysozyme (T4L) both exhibit such a phenomenon. Native-state hydrogen exchange (NHX) data have been collected for these proteins, providing residue-specific free energies and m-values (a measure of hydrocarbon solvation) for the manifold of partially unfolded, exchange-competent forms that are accessible from the native state (DeltaG(sg) and m(sg), where the sg subscript denotes sub-global). There is good evidence that these parameters pertain to exchange-competent species comprising the burst phase observed in the global folding kinetics. We combine the results from the global folding kinetics of these proteins with a statistical analysis of their NHX parameters to determine if the distribution of experimental (m(sg), DeltaG(sg)) values derive from a mechanism where the burst phase is two-state. For RNase H, this analysis demonstrates that the burst phase of this protein is not two-state; the results imply a distribution of states, m and DeltaG exhibiting a linear functional relationship consistent with the global folding parameters. For T4L, it is difficult to distinguish the observed distribution of m(sg), DeltaG(sg) values from that expected for a mechanism where the burst phase is two-state. The results for RNase H* lend support for the idea that the burst phase reaction of this protein comprises a continuum of states. This has important implications for how we model the process of structural acquisition in protein folding reactions.

Conservation of folding and stability within a protein family: the tyrosine corner as an evolutionary cul-de-sac

What are the selective pressures on protein sequences during evolution? Amino acid residues may be highly conserved for functional or structural (stability) reasons. Theoretical studies have proposed that residues involved in the folding nucleus may also be highly conserved. To test this we are using an experimental "fold approach" to the study of protein folding. This compares the folding and stability of a number of proteins that share the same fold, but have no common amino acid sequence or biological activity. The fold selected for this study is the immunoglobulin-like beta-sandwich fold, which is a fold that has no specifically conserved function. Four model proteins are used from two distinct superfamilies that share the immunoglobulin-like fold, the fibronectin type III and immunoglobulin superfamilies. Here, the fold approach and protein engineering are used to question the role of a highly conserved tyrosine in the "tyrosine corner" motif that is found ubiquitously and exclusively in Greek key proteins. In the four model beta-sandwich proteins characterised here, the tyrosine is the only residue that is absolutely conserved at equivalent sites. By mutating this position to phenylalanine, we show that the tyrosine hydroxyl is not required to nucleate folding in the immunoglobulin superfamily, whereas it is involved to some extent in early structure formation in the fibronectin type III superfamily. The tyrosine corner is important for stability, mutation to phenylalanine costs between 1.5 and 3 kcal mol(-1). We propose that the high level of conservation of the tyrosine is related to the structural restraints of the loop connecting the beta-sheets, representing an evolutionary "cul-de-sac".

Folding of beta-sandwich proteins: three-state transition of a fibronectin type III module

An analysis of the folding of the 94 residue tenth fibronectin type III (fnIII) domain of human fibronectin (FNfn10) is presented. Use of guanidine isothiocyanate as a denaturant allows us to obtain equilibrium and kinetic data across a broad range of denaturant concentrations that are unavailable in guanidine hydrochloride. Equilibrium unfolding experiments show that FNfn10 is significantly more stable than has been reported previously. Comparison of equilibrium and kinetic parameters reveals the presence of an intermediate that accumulates at low denaturant concentrations. This is the first demonstration of three-state folding kinetics for a fnIII domain. We have previously shown that a homologous domain from human tenascin (TNfn3) folds by a two-state mechanism, but this does not necessarily indicate that the two proteins fold by different folding pathways.

Role of hydration water in protein unfolding

In this paper, following our work on the two-state outer neighbor mixed bonding model of water, it is proposed that polar groups promote the formation of the low density ice Ih-type bonding in their neighborhood, whereas nonpolar groups tend to promote the higher density ice II-type structure. In a protein, because of the large numbers of exposed polar and nonpolar groups, large changes in the neighboring water structure can occur. These changes, of course, depend on whether the protein is in its native or its unfolded state and will be shown here to have a direct impact on the thermodynamics of protein unfolding at both high and low temperatures. For example, it is known that the polar hydration entropies become rapidly more negative with increasing temperature. This very unusual behavior can be directly related to the promotion in the outer bulk liquid of the more stable Ih-type bonding at the expense of II-type bonding by polar groups of the protein. In contrast, nonpolar groups have an opposite effect on the thermodynamics. It is the delicate balance created by these outer hydration contributions, mixed with ordinary thermodynamic contributions from the inner hydration shell and those from hydrogen-bond and van der Waals forces within the protein molecule itself that is responsible for both heat and cold denaturation of proteins.

The cooperativity of burst phase reactions explored

The denaturant-dependence of the major, observable relaxation rates for folding (kobs) of ribonuclease HI from Escherichia coli (RNase H) and phage T4 lysozyme (T4L) reveal that, for both proteins, folding begins with the rapid and transient accumulation of intermediate species in a "burst phase" which precedes the rate-limiting formation of the native state; this is evidenced by a "rollover" in the folding limb of the rate profiles (kobs versus denaturant, or chevron plot). These rate profiles are most simply described by a three-state mechanism (unfolded-to-intermediate-to-native), which implies that the burst phase represents a transition between two distinct thermodynamic states. It is shown here that the equilibrium properties of these burst phase reactions can be equally well modeled by a mechanism involving a continuum of states where the free energy of each state is linearly related to its m-value (the parameter describing the linear relationship between free energy and denaturant). A numerical model is also developed to describe the time evolution of such a system, which exhibits nearly perfect exponential behavior. Both models emphasize how a continuum of states operating under a linear free energy relationship may behave like a two state system. Such a scheme finds experimental justification from an interpretation of recent native state hydrogen exchange data. The analytical model described for a continuum can account for the observed kinetic profiles of several other model proteins. The results, however, appear context specific, suggesting that burst phase reactions are not entirely random and non-specific. The results reported in this study have important implications for the concept of cooperativity in protein folding reactions.