Dynamic NMR line-shape analysis demonstrates that the villin headpiece subdomain folds on the microsecond time scale

Experimental Evidence for Millisecond-Timescale Structural Evolution Following the Microsecond-Timescale Folding of a Small Protein

Prior work has shown that small proteins can fold (i.e., convert from unstructured to structured states) within 10  μs. Here we use time-resolved solid state nuclear magnetic resonance (ssNMR) methods to show that full folding of the 35-residue villin headpiece subdomain (HP35) requires a slow annealing process that has not been previously detected. ^{13}C ssNMR spectra of frozen HP35 solutions, acquired with a variable time τ_{e} at 30 °C after rapid cooling from 95 °C and before rapid freezing, show changes on the 3-10 ms timescale, attributable to slow rearrangements of protein sidechains during τ_{e}.

Line Shape Analysis of 19F NMR-Monitored Chemical Denaturation of a Fold-Switching Protein RfaH Reveals Its Slow Folding Dynamics

The recent discovery of metamorphic proteins, which can switch between multiple conformations under native conditions, has challenged the well-established one sequence-one structure paradigm of protein folding. This is exemplified in the C-terminal domain of the multidomain transcription factor RfaH, which converts from an α-helical coiled-coil conformation in its autoinhibited state to a β-barrel conformation upon activation. Here, we use multisite line shape analysis of 19F NMR-monitored equilibrium chemical denaturation measurements of two 19F-labeled variants of full-length RfaH, to show that it folds/unfolds slowly on the NMR time scale, in an apparent all-or-none fashion at physiological pH and room temperature in the 3.3-4.8 M urea concentration range. The significant thermodynamic stability and slow unfolding rate (kinetic stability) are likely responsible for maintaining the closed autoinhibited state of RfaH, preventing spurious interactions with RNA polymerase (RNAP) when not functional. Our results provide a quantitative understanding of the folding-function relationship in the model fold-switching protein RfaH.

Monitoring protein unfolding transitions by NMR-spectroscopy

NMR-spectroscopy has certain unique advantages for recording unfolding transitions of proteins compared e.g. to optical methods. It enables per-residue monitoring and separate detection of the folded and unfolded state as well as possible equilibrium intermediates. This allows a detailed view on the state and cooperativity of folding of the protein of interest and the correct interpretation of subsequent experiments. Here we summarize in detail practical and theoretical aspects of such experiments. Certain pitfalls can be avoided, and meaningful simplification can be made during the analysis. Especially a good understanding of the NMR exchange regime and relaxation properties of the system of interest is beneficial. We show by a global analysis of signals of the folded and unfolded state of GB1 how accurate values of unfolding can be extracted and what limits different NMR detection and unfolding methods. E.g. commonly used exchangeable amides can lead to a systematic under determination of the thermodynamic protein stability. We give several perspectives of how to deal with more complex proteins and how the knowledge about protein stability at residue resolution helps to understand protein properties under crowding conditions, during phase separation and under high pressure.

Structure and Carbohydrate Recognition by the Nonmitogenic Lectin Horcolin

Lectins are sugar-binding proteins that have shown considerable promise as antiviral agents because of their ability to interact with envelope glycoproteins present on the surface of viruses such as HIV-1. However, their therapeutic potential has been compromised by their mitogenicity that stimulates uncontrolled division of T-lymphocytes. Horcolin, a member of the jacalin family of lectins, tightly binds the HIV-1 envelope glycoprotein gp120 and neutralizes HIV-1 particles but is nonmitogenic. In this report, we combine X-ray crystallography and NMR spectroscopy to obtain atomic-resolution insights into the structure of horcolin and the molecular basis for its carbohydrate recognition. Each protomer of the horcolin dimer adopts a canonical β-prism I fold with three Greek key motifs and carries two carbohydrate-binding sites. The carbohydrate molecule binds in a negatively charged pocket and is stabilized by backbone and side chain hydrogen bonds to conserved residues in the ligand-binding loop. NMR titrations reveal a two-site binding mode and equilibrium dissociation constants for the two binding sites determined from two-dimensional (2D) lineshape modeling are 4-fold different. Single-binding-site variants of horcolin confirm the dichotomy in binding sites and suggest that there is allosteric communication between the two sites. An analysis of the horcolin structure shows a network of hydrogen bonds linking the two carbohydrate-binding sites directly and through a secondary binding site, and this coupling between the two sites is expected to assume importance in the interaction of horcolin with high-mannose glycans found on viral envelope glycoproteins.

Quantitative Analysis of Protein Unfolded State Energetics: Experimental and Computational Studies Demonstrate That Non-Native Side-Chain Interactions Stabilize Local Native Backbone Structure

Proteins fold on relatively smooth free energy landscapes which are biased toward the native state, but even simple topologies which fold rapidly can experience roughness on their free energy landscape. The details of these interactions are difficult to elucidate experimentally. Closely related to the problem of deciphering the details of the free energy landscape is the problem of defining the interactions in the denatured state ensemble (DSE) which is populated under native conditions, that is, under conditions where the native state is stable. The DSE of many proteins deviates from random coil models, but quantifying and defining the energetics of the transiently populated interactions in this ensemble is extremely challenging. Characterization of the DSE of proteins which fold to compact structures is also relevant to studies of intrinsically disordered proteins (IDPs) since interactions in the dynamic ensemble populated by IDPs can modulate their behavior. Here we show how experimental thermodynamic and pK_a measurements can be combined with computational thermodynamic integration to quantify interactions in the DSE. We show that non-native side chain interactions can stabilize native backbone structure in the DSE and demonstrate that that even rapidly folding proteins can form energetically significant non-native interactions in their DSE. As an example, we characterize a non-native salt bridge that stabilizes local native backbone structure in the DSE of a widely studied fast-folding protein, the villin headpiece helical domain. The combined computational experimental approach is applicable to other protein unfolded states and provides insight that is impossible to obtain with either method alone.

Heterogeneity in the Folding of Villin Headpiece Subdomain HP36

Small single domain proteins that fold on the microsecond time scale have been the subject of intense interest as models for probing the complexity of folding energy landscapes. The villin headpiece subdomain (HP36) has been extensively studied because of its simple three helix structure, ultrafast folding lifetime of a few microseconds, and stable native fold. We have previously shown that folding as measured by a single ¹³C═¹⁸O isotopic label on residue A57 in helix 2 occurs at a different rate than that measured by global probes of folding, indicating noncooperative complexity in the folding of HP36. In order to determine whether this complexity reflects intermediates or parallel pathways over a small activation barrier, ¹³C═¹⁸O labels were individually incorporated at six different positions in HP36, including into all 3 helices. The equilibrium thermal unfolding transitions and the folding/unfolding dynamics were monitored using the unique IR signature of the ¹³C═¹⁸O label by temperature dependent FTIR and temperature jump IR spectroscopy, respectively. Equilibrium experiments reveal that the ¹³C═¹⁸O labels at different positions in HP36 show drastic differences in the midpoint of their transitions ( T_m), ranging from 45 to 67 °C. Heterogeneity is also observed in the relaxation kinetics; there are differences in the microsecond phase when different labeled positions are probed. At a final temperature of 45 °C, the relaxation rate for ¹³C═¹⁸O A57 is 2.4e + 05 s^-1 whereas for ¹³C═¹⁸O L69 HP36 the relaxation rate is 5.1e + 05 s^-1, two times faster. The observation of site-dependent midpoints for the equilibrium unfolding transitions and differences in the relaxation rates of the labeled positions enables us to probe the progressive accumulation of the folded structure, providing insight into the microscopic details of the folding mechanism.

Slow-Down in Diffusion in Crowded Protein Solutions Correlates with Transient Cluster Formation

For a long time, the effect of a crowded cellular environment on protein dynamics has been largely ignored. Recent experiments indicate that proteins diffuse more slowly in a living cell than in a diluted solution, and further studies suggest that the diffusion depends on the local surroundings. Here, detailed insight into how diffusion depends on protein-protein contacts is presented based on extensive all-atom molecular dynamics simulations of concentrated villin headpiece solutions. After force field adjustments in the form of increased protein-water interactions to reproduce experimental data, translational and rotational diffusion was analyzed in detail. Although internal protein dynamics remained largely unaltered, rotational diffusion was found to slow down more significantly than translational diffusion as the protein concentration increased. The decrease in diffusion is interpreted in terms of a transient formation of protein clusters. These clusters persist on sub-microsecond time scales and follow distributions that increasingly shift toward larger cluster size with increasing protein concentrations. Weighting diffusion coefficients estimated for different clusters extracted from the simulations with the distribution of clusters largely reproduces the overall observed diffusion rates, suggesting that transient cluster formation is a primary cause for a slow-down in diffusion upon crowding with other proteins.

Optimized purification of a fusion protein by reversed-phase high performance liquid chromatography informed by the linear solvent strength model

Fusion protein systems are commonly used for expression of small proteins and peptides. An important criterion for a fusion protein system to be useful is the ability to separate the protein of interest from the tag. Additionally, because no protease cleaves fusion proteins with 100% efficiency, the ability to separate the desired peptide from any remaining uncleaved protein is also necessary. This is likely to be the more difficult task as at least a portion of the sequence of the fusion protein is identical to that of the protein of interest. When a high level of purity is required, gradient elution reversed-phase HPLC is frequently used as a final purification step. Shallow gradients are often advantageous for maximizing both the purity and yield of the final product; however, the relationship between relative retention times at shallow gradients and those at steeper gradients typically used for analytical HPLC are not always straightforward. In this work, we report reversed-phase HPLC results for the fusion protein system consisting of the N-terminal domain of ribosomal protein L9 (NTL9) and the 36-residue villin headpiece subdomain (HP36) linked by a recognition sequence for the protease factor Xa. This system represents an excellent example of the difficulties in purification that may arise from this unexpected elution behavior at shallow gradients. Additionally, we report on the sensitivity of this elution behavior to the concentration of the additive trifluoroacetic acid in the mobile phase and present optimized conditions for separating HP36 from the full fusion protein by reversed-phase HPLC using a shallow gradient. Finally, we suggest that these findings are relevant to the purification of other fusion protein systems, for which similar problems may arise, and support this suggestion using insights from the linear solvent strength model of gradient elution liquid chromatography.

Studying the role of cooperative hydration in stabilizing folded protein states

Understanding and modelling protein folding remains a key scientific and engineering challenge. Two key questions in protein folding are (1) why many proteins adopt a folded state and (2) how these proteins transition from the random coil ensemble to a folded state. In this paper we employ molecular dynamics simulations to address the first of these questions. Computational methods are well-placed to address this issue due to their ability to analyze systems at atomic-level resolution. Traditionally, the stability of folded proteins has been ascribed to the balance of two types of intermolecular interactions: hydrogen-bonding interactions and hydrophobic contacts. In this study, we explore a third type of intermolecular interaction: cooperative hydration of protein surface residues. To achieve this, we consider multiple independent simulations of the villin headpiece domain to quantify the contributions of different interactions to the energy of the native and fully extended states. In addition, we consider whether these findings are robust with respect to the protein forcefield, the water model, and the presence of salt. In all cases, we identify many cooperatively hydrated interactions that are transient but energetically favor the native state. Whilst further work on additional protein structures, forcefields, and water models is necessary, these results suggest a role for cooperative hydration in protein folding that should be explored further. Rational design of cooperative hydration on the protein surface could be a viable strategy for increasing protein stability.

Selenomethionine, p-cyanophenylalanine pairs provide a convenient, sensitive, non-perturbing fluorescent probe of local helical structure

The use of selenomethionine (MSe)-p-cyanophenylalanine (FCN) pairs to probe protein structure is demonstrated. MSe quenches FCN fluorescence via electron transfer. Both residues can be incorporated recombinantly or by peptide synthesis. Time-resolved and steady-state fluorescence measurements demonstrate that MSe-FCN pairs provide specific local probes of helical structure.

Folding stability modulation of the villin headpiece helical subdomain by 4-fluorophenylalanine and 4-methylphenylalanine

HP36, the helical subdomain of villin headpiece, contains a hydrophobic core composed of three phenylalanine residues (Phe47, Phe51, and Phe58). Hydrophobic effects and electrostatic interactions were shown to be the critical factors in stabilizing this core and the global structure. To assess the interactions among Phe47, Phe51, and Phe58 residues and investigate how they affect the folding stability, we implanted 4-fluorophenylalanine (Z) and 4-methylphenylalanine (X) into the hydrophobic core of HP36. We chemically synthesized HP36 and its seven variants including four single mutants whose Phe51 or Phe58 was replaced with Z or X, and three double mutants whose Phe51 and Phe58 were both substituted. Circular dichroism and nuclear magnetic resonance measurements show that the variants exhibit a native HP36 like fold, of which F51Z and three double mutants are more stable than the wild type. Molecular modeling provided detailed interaction energy within the phenylalanine residues, revealing that electrostatic interactions dominate the stability modulation upon the introduction of 4-fluorophenylalanine and 4-methylphenylalanine. Our results show that these two non-natural amino acids can successfully tune the interactions in a relatively compact hydrophobic core and the folding stability without inducing dramatic steric effects. Such an approach may be applied to other folded motifs or proteins.

Two-Dimensional NMR Lineshape Analysis

NMR titration experiments are a rich source of structural, mechanistic, thermodynamic and kinetic information on biomolecular interactions, which can be extracted through the quantitative analysis of resonance lineshapes. However, applications of such analyses are frequently limited by peak overlap inherent to complex biomolecular systems. Moreover, systematic errors may arise due to the analysis of two-dimensional data using theoretical frameworks developed for one-dimensional experiments. Here we introduce a more accurate and convenient method for the analysis of such data, based on the direct quantum mechanical simulation and fitting of entire two-dimensional experiments, which we implement in a new software tool, TITAN (TITration ANalysis). We expect the approach, which we demonstrate for a variety of protein-protein and protein-ligand interactions, to be particularly useful in providing information on multi-step or multi-component interactions.

Characterizing a partially ordered miniprotein through folding molecular dynamics simulations: Comparison with the experimental data

The villin headpiece helical subdomain (HP36) is one of the best known model systems for computational studies of fast-folding all-α miniproteins. HP21 is a peptide fragment-derived from HP36-comprising only the first and second helices of the full domain. Experimental studies showed that although HP21 is mostly unfolded in solution, it does maintain some persistent native-like structure as indicated by the analysis of NMR-derived chemical shifts. Here we compare the experimental data for HP21 with the results obtained from a 15-μs long folding molecular dynamics simulation performed in explicit water and with full electrostatics. We find that the simulation is in good agreement with the experiment and faithfully reproduces the major experimental findings, namely that (a) HP21 is disordered in solution with <10% of the trajectory corresponding to transiently stable structures, (b) the most highly populated conformer is a native-like structure with an RMSD from the corresponding portion of the HP36 crystal structure of <1 Å, (c) the simulation-derived chemical shifts-over the whole length of the trajectory-are in reasonable agreement with the experiment giving reduced χ(2) values of 1.6, 1.4, and 0.8 for the Δδ(13) C(α) , Δδ(13) CO, and Δδ(13) C(β) secondary shifts, respectively (becoming 0.8, 0.7, and 0.3 when only the major peptide conformer is considered), and finally, (d) the secondary structure propensity scores are in very good agreement with the experiment and clearly indicate the higher stability of the first helix. We conclude that folding molecular dynamics simulations can be a useful tool for the structural characterization of even marginally stable peptides.

Equilibrium kinetic network of the villin headpiece in implicit solvent

We applied the single-replica multiple-state transition-interface sampling method to elucidate the equilibrium kinetic network of the 35-residue-fragment (HP-35) villin headpiece in implicit water at room temperature. Starting from the native Protein Data Bank structure, nine (meta)stable states of the system were identified, from which the kinetic network was built by sampling pathways between these states. Application of transition path theory allowed analysis of the (un)folding mechanism. The resulting (un)folding rates agree well with experiments. This work demonstrates that high (un)folding barriers can now be studied.

Effects of J couplings and unobservable minor states on kinetics parameters extracted from CEST data

Chemical exchange saturation transfer (CEST) experiments have emerged as a powerful tool for characterizing dynamics and sparse populated conformers of protein in slow exchanging systems. We show that J couplings and 'invisible' minor states can cause systematic errors in kinetics parameters and chemical shifts extracted from CEST data. For weakly coupled spin systems, the J coupling effect can be removed using an approximation method. This method is warranted through detailed theoretical derivation, supported by results from simulations and experiments on an acyl carrier protein domain. Simulations demonstrate that the effect of 'invisible' minor states on the extracted kinetics parameters depends on the chemical shifts, populations, exchange rates of the 'invisible' states to the observed major or minor state and exchange models. Moreover, the extracted chemical shifts of the observed minor state can also be influenced by the "invisible" minor states. The presence of an off-pathway folding intermediate in the acyl carrier protein domain explains why the exchange rates obtained with a two-state model from individual residues that displayed only two obvious CEST dips varied significantly and the extracted exchange rates for ¹⁵N and ¹³CO spins located in the same peptide bond could be very different. The approximation method described here simplifies CEST data analysis in many situations where the coupling effect cannot be ignored and decoupling techniques are not desirable. In addition, this study also raises alerts for 'invisible' minor states which can cause errors in not only kinetics parameters but also chemical shifts of the observed minor state.

Balanced Protein-Water Interactions Improve Properties of Disordered Proteins and Non-Specific Protein Association

Some frequently encountered deficiencies in all-atom molecular simulations, such as nonspecific protein-protein interactions being too strong, and unfolded or disordered states being too collapsed, suggest that proteins are insufficiently well solvated in simulations using current state-of-the-art force fields. To address these issues, we make the simplest possible change, by modifying the short-range protein-water pair interactions, and leaving all the water-water and protein-protein parameters unchanged. We find that a modest strengthening of protein-water interactions is sufficient to recover the correct dimensions of intrinsically disordered or unfolded proteins, as determined by direct comparison with small-angle X-ray scattering (SAXS) and Förster resonance energy transfer (FRET) data. The modification also results in more realistic protein-protein affinities, and average solvation free energies of model compounds which are more consistent with experiment. Most importantly, we show that this scaling is small enough not to affect adversely the stability of the folded state, with only a modest effect on the stability of model peptides forming α-helix and β-sheet structures. The proposed adjustment opens the way to more accurate atomistic simulations of proteins, particularly for intrinsically disordered proteins, protein-protein association, and crowded cellular environments.

Statistical analyses of protein folding rates from the view of quantum transition

Understanding protein folding rate is the primary key to unlock the fundamental physics underlying protein structure and its folding mechanism. Especially, the temperature dependence of the folding rate remains unsolved in the literature. Starting from the assumption that protein folding is an event of quantum transition between molecular conformations, we calculated the folding rate for all two-state proteins in a database and studied their temperature dependencies. The non-Arrhenius temperature relation for 16 proteins, whose experimental data had previously been available, was successfully interpreted by comparing the Arrhenius plot with the first-principle calculation. A statistical formula for the prediction of two-state protein folding rate was proposed based on quantum folding theory. The statistical comparisons of the folding rates for 65 two-state proteins were carried out, and the theoretical vs. experimental correlation coefficient was 0.73. Moreover, the maximum and the minimum folding rates given by the theory were consistent with the experimental results.

Probing the origins of two-state folding

Many protein systems fold in a two-state manner. Random models, however, rarely display two-state kinetics and thus such behavior should not be accepted as a default. While theories for the prevalence of two-state kinetics have been presented, none sufficiently explain the breadth of experimental observations. A model, making minimal assumptions, is introduced that suggests two-state behavior is likely for any system with an overwhelmingly populated native state. We show two-state folding is a natural consequence of such two-state thermodynamics, and is strengthened by increasing the population of the native state. Further, the model exhibits hub-like behavior, with slow interconversions between unfolded states. Despite this, the unfolded state equilibrates quickly relative to the folding time. This apparent paradox is readily understood through this model. Finally, our results compare favorable with measurements of folding rates as a function of chain length and Keq, providing new insight into these relations.

Electronic polarization stabilizes tertiary structure prediction of HP-36

Molecular dynamic (MD) simulations with both implicit and explicit solvent models have been carried out to study the folding dynamics of HP-36 protein. Starting from the extended conformation, the secondary structure of all three helices in HP-36 was formed in about 50 ns and remained stable in the remaining simulation. However, the formation of the tertiary structure was difficult. Although some intermediates were close to the native structure, the overall conformation was not stable. Further analysis revealed that the large structure fluctuation of loop and hydrophobic core regions was devoted mostly to the instability of the structure during MD simulation. The backbone root-mean-square deviation (RMSD) of the loop and hydrophobic core regions showed strong correlation with the backbone RMSD of the whole protein. The free energy landscape indicated that the distribution of main chain torsions in loop and turn regions was far away from the native state. Starting from an intermediate structure extracted from the initial AMBER simulation, HP-36 was found to generally fold to the native state under the dynamically adjusted polarized protein-specific charge (DPPC) simulation, while the peptide did not fold into the native structure when AMBER force filed was used. The two best folded structures were extracted and taken into further simulations in water employing AMBER03 charge and DPPC for 25 ns. Result showed that introducing polarization effect into interacting potential could stabilize the near-native protein structure.

Fast protein folding kinetics

Fast-folding proteins have been a major focus of computational and experimental study because they are accessible to both techniques: they are small and fast enough to be reasonably simulated with current computational power, but have dynamics slow enough to be observed with specially developed experimental techniques. This coupled study of fast-folding proteins has provided insight into the mechanisms, which allow some proteins to find their native conformation well <1 ms and has uncovered examples of theoretically predicted phenomena such as downhill folding. The study of fast folders also informs our understanding of even 'slow' folding processes: fast folders are small; relatively simple protein domains and the principles that govern their folding also govern the folding of more complex systems. This review summarizes the major theoretical and experimental techniques used to study fast-folding proteins and provides an overview of the major findings of fast-folding research. Finally, we examine the themes that have emerged from studying fast folders and briefly summarize their application to protein folding in general, as well as some work that is left to do.

Salt effects on hydrophobic-core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica-exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt-bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near-native conformations in pure water. Our analyses also show that such salt-bridge formation hinders the fully native hydrophobic-core packing at the final stages of folding.

NMR at low and ultralow temperatures

Solid state nuclear magnetic resonance (NMR) measurements at low temperatures have been common in physical sciences for many years and are becoming increasingly important in studies of biomolecular systems. This Account reviews a diverse set of projects from my laboratory, dating back to the early 1990s, that illustrate the motivations for low-temperature solid state NMR, the types of information that are available from the measurements, and likely directions for future research. These projects include NMR studies of both physical and biological systems, performed at low (cooled with nitrogen, down to 77 K) and ultralow (cooled with helium, below 77 K) temperatures, and performed with and without magic-angle spinning (MAS). NMR studies of physical systems often focus on phenomena that occur only at low temperatures. Two examples from my laboratory are studies of molecular rotation and orientational ordering in solid C60 at low temperatures and studies of unusual electronic states, called skyrmions, in two-dimensionally confined electron systems within semiconductor quantum wells. To study quantum wells, we used optical pumping of nuclear spin polarizations to enhance their NMR signals. The optical pumping phenomenon exists only at ultralow temperatures. In studies of biomolecular systems, low-temperature NMR has several motivations. In some cases, low temperatures suppress molecular tumbling, thereby permitting solid state NMR measurements on soluble proteins. Studies of AIDS-related peptide/antibody complexes illustrate this effect. In other cases, low temperatures suppress conformational exchange, thereby permitting quantitation of conformational distributions. Studies of chemically denatured states of the model protein HP35 illustrate this effect. Low temperatures and rapid freeze-quenching can also be used to trap transient intermediate states in a non-equilibrium kinetic process, as shown in studies of a transient intermediate in the rapid folding pathway of HP35. NMR sensitivity generally increases with decreasing sample temperature. Therefore, it can be useful to carry out experiments at the lowest possible temperatures, particularly in studies of biomolecular systems in frozen solutions. However, solid state NMR studies of biomolecular systems generally require rapid MAS. A novel MAS NMR probe design that uses nitrogen gas for sample spinning and cold helium only for sample cooling allows a wide variety of solid state NMR measurements to be performed on biomolecular systems at 20-25 K, where signals are enhanced by factors of 12-15 relative to measurements at room temperature. MAS NMR at ultralow temperatures also facilitates dynamic nuclear polarization (DNP), allowing sizeable additional signal enhancements and large absolute NMR signal amplitudes with relatively low microwave powers. Current research in my laboratory seeks to develop and exploit DNP-enhanced MAS NMR at ultralow temperatures, for example, in studies of transient intermediates in protein folding and aggregation processes and studies of peptide/protein complexes that can be prepared only at low concentrations.

Rational modification of protein stability by targeting surface sites leads to complicated results

The rational modification of protein stability is an important goal of protein design. Protein surface electrostatic interactions are not evolutionarily optimized for stability and are an attractive target for the rational redesign of proteins. We show that surface charge mutants can exert stabilizing effects in distinct and unanticipated ways, including ones that are not predicted by existing methods, even when only solvent-exposed sites are targeted. Individual mutation of three solvent-exposed lysines in the villin headpiece subdomain significantly stabilizes the protein, but the mechanism of stabilization is very different in each case. One mutation destabilizes native-state electrostatic interactions but has a larger destabilizing effect on the denatured state, a second removes the desolvation penalty paid by the charged residue, whereas the third introduces unanticipated native-state interactions but does not alter electrostatics. Our results show that even seemingly intuitive mutations can exert their effects through unforeseen and complex interactions.

Origin of abrupt rise in deuteron NMR longitudinal relaxation times of protein methyl groups below 90 K

In order to examine the origin of the abrupt change in the temperature dependence of (2)H NMR longitudinal relaxation times observed previously for methyl groups of L69 in the hydrophobic core of villin headpiece protein at around 90 K (Vugmeyster et al. J. Am. Chem. Soc. 2010, 132, 4038-4039), we extended the measurements to several other methyl groups in the hydrophobic core. We show that, for all methyl groups, relaxation times experience a dramatic jump several orders of magnitude around this temperature. Theoretical modeling supports the conclusion that the origin of the apparent transition in the relaxation times is due to the existence of the distribution of conformers distinguished by their activation energy for methyl three-site hops. It is also crucial to take into account the differential contribution of individual conformers into overall signal intensity. When a particular conformer approaches the regime at which its three-site hop rate constant is on the order of the quadrupolar coupling interaction constant, the intensity of the signal due to this conformer experiences a sharp drop, thus changing the balance of the contributions of different conformers into the overall signal. As a result, the observed apparent transition in the relaxation rates can be explained without the assumption of an underlying transition in the rate constants. This work in combination with earlier results also shows that the model based on the distribution of conformers explains the relaxation behavior in the entire temperature range between 300 and 70 K.

Mutational effects on the folding dynamics of a minimized hairpin

The fold stabilities and folding dynamics of a series of mutants of a model hairpin, KTW-NPATGK-WTE (HP7), are reported. The parent system and the corresponding DPATGK loop species display submicrosecond folding time constants. The mutational studies revealed that ultrafast folding requires both some prestructuring of the loop and a favorable interaction between the chain termini in the transition state. In the case of YY-DPETGT-WY, another submicrosecond folding species [Davis, C. M., Xiao, S., Raleigh, D. P., and Dyer, R. B. (2012) J. Am. Chem. Soc. 134, 14476-14482], a hydrophobic cluster provides the latter. In the case of HP7, the Coulombic interaction between the terminal NH3(+) and CO2(-) units provides this; a C-terminal Glu to amidated Ala mutation results in a 5-fold retardation of the folding rate. The effects of mutations within the reversing loop indicate the balance between loop flexibility (favoring fast conformational searching) and turn formation in the unfolded state is a major factor in determining the folding dynamics. The -NAAAKX- loops examined display no detectable turn formation propensity in other hairpin constructs but do result in stable analogues of HP7. Peptide KTW-NAAAKK-WTE displays the same fold stability as HP7, but both the folding and unfolding time constants are greater by a factor of 20.

Reduced native state stability in crowded cellular environment due to protein-protein interactions

The effect of cellular crowding environments on protein structure and stability is a key issue in molecular and cellular biology. The classical view of crowding emphasizes the volume exclusion effect that generally favors compact, native states. Here, results from molecular dynamics simulations and NMR experiments show that protein crowders may destabilize native states via protein-protein interactions. In the model system considered here, mixtures of villin head piece and protein G at high concentrations, villin structures become increasingly destabilized upon increasing crowder concentrations. The denatured states observed in the simulation involve partial unfolding as well as more subtle conformational shifts. The unfolded states remain overall compact and only partially overlap with unfolded ensembles at high temperature and in the presence of urea. NMR measurements on the same systems confirm structural changes upon crowding based on changes of chemical shifts relative to dilute conditions. An analysis of protein-protein interactions and energetic aspects suggests the importance of enthalpic and solvation contributions to the crowding free energies that challenge an entropic-centered view of crowding effects.

Glassy dynamics of protein methyl groups revealed by deuteron NMR

We investigated site-specific dynamics of key methyl groups in the hydrophobic core of chicken villin headpiece subdomain (HP36) over the temperature range between 298 and 140 K using deuteron solid-state NMR longitudinal relaxation measurements. The relaxation of the longitudinal magnetization is weakly nonexponential (glassy) at high temperatures and exhibits a stronger degree of nonexponentiality below about 175 K. In addition, the characteristic relaxation times deviate from the simple Arrhenius law. We interpret this behavior via the existence of distribution of activation energy barriers for the three-site methyl jumps, which originates from somewhat different methyl environments within the local energy landscape. The width of the distribution of the activation barriers for methyl jumps is rather significant, about 1.4 kJ/mol. Our experimental results and modeling allow for the description of the apparent change at about 175 K without invoking a specific transition temperature. For most residues in the core, the relaxation behavior at high temperatures points to the existence of conformational exchange between the substates of the landscape, and our model takes into account the kinetics of this process. The observed dynamics are the same for dry and hydrated protein. We also looked at the effect of F58L mutation inside the hydrophobic core on the dynamics of one of the residues and observed a significant increase in its conformational exchange rate constant at high temperatures.

Thermally induced protein unfolding probed by isotope-edited IR spectroscopy

Infrared (IR) spectroscopy has been widely utilized for the study of protein folding, unfolding, and misfolding processes. We have previously developed a theoretical method for calculating IR spectra of proteins in the amide I region. In this work, we apply this method, in combination with replica-exchange molecular dynamics simulations, to study the equilibrium thermal unfolding transition of the villin headpiece subdomain (HP36). Temperature-dependent IR spectra and spectral densities are calculated. The spectral densities correctly reflect the unfolding conformational changes in the simulation. With the help of isotope labeling, we are able to capture the feature that helix 2 of HP36 loses its secondary structure before global unfolding occurs, in agreement with experiment.

Conformational dynamics and stability of HP35 studied with 2D IR vibrational echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.

Simple few-state models reveal hidden complexity in protein folding

Markov state models constructed from molecular dynamics simulations have recently shown success at modeling protein folding kinetics. Here we introduce two methods, flux PCCA+ (FPCCA+) and sliding constraint rate estimation (SCRE), that allow accurate rate models from protein folding simulations. We apply these techniques to fourteen massive simulation datasets generated by Anton and Folding@home. Our protocol quantitatively identifies the suitability of describing each system using two-state kinetics and predicts experimentally detectable deviations from two-state behavior. An analysis of the villin headpiece and FiP35 WW domain detects multiple native substates that are consistent with experimental data. Applying the same protocol to GTT, NTL9, and protein G suggests that some beta containing proteins can form long-lived native-like states with small register shifts. Even the simplest protein systems show folding and functional dynamics involving three or more states.

Temperature dependence of water interactions with the amide carbonyls of α-helices

Hydration is a key determinant of the folding, dynamics, and function of proteins. In this study, temperature-dependent Fourier transform infrared (FTIR) spectroscopy combined with singular value decomposition (SVD) and global fitting were used to investigate both the interaction of water with α-helical proteins and the cooperative thermal unfolding of these proteins. This methodology has been applied to an isolated α-helix (Fs peptide) and to globular α-helical proteins including the helical subdomain and full-length villin headpiece (HP36 and HP67). The results suggest a unique IR signature for the interaction of water with the helical amide carbonyl groups of the peptide backbone. The IR spectra indicate a weakening of the net hydrogen bond strength of water to the backbone carbonyls with increasing temperature. This weakening of the backbone solvation occurs as a discrete transition near the maximum of the temperature-dependent hydrophobic effect, not a continuous change with increasing temperature. Possible molecular origins of this effect are discussed with respect to previous molecular dynamics simulations of the temperature-dependent solvation of the helix backbone.

Microsecond folding dynamics of apomyoglobin at acidic pH

Apomyolgobin (apoMb) is an important model for understanding the folding mechanism of helical proteins. This study focuses on a partially structured state of sperm whale apoMb populated at pH 4.2 (M-state), which structurally resembles a late kinetic intermediate in the formation of the native state (N) at higher pH. The thermodynamics and cooperativity of apoMb folding at pH 4.2 and 6.2 were studied by global analysis of the urea-induced unfolding transitions monitored by tryptophan fluorescence and circular dichroism. The kinetics of folding and unfolding of apoMb at pH 4.2 was measured over a time window from 40 to 850 μs, using fluorescence-detected continuous-flow measurements. Our observation of biphasic kinetics provides clear evidence for rapid (<100 μs) accumulation of previously unresolved intermediate states in both refolding and unfolding experiments. Quantitative kinetic modeling of the results, using a four-state mechanism with two intermediates on a direct route between the unfolded and folded states (U↔I↔L↔M), gave new insight into the conformational states and barriers that precede the rate-limiting step in the formation of the N-state of apoMb.

Slow motions in the hydrophobic core of chicken villin headpiece subdomain and their contributions to configurational entropy and heat capacity from solid-state deuteron NMR measurements

We have investigated microsecond to millisecond time scale dynamics in several key hydrophobic core methyl groups of chicken villin headpiece subdomain protein (HP36) using a combination of single-site labeling, deuteron solid-state NMR line shape analysis, and computational modeling. Deuteron line shapes of hydrated powder samples are dominated by rotameric jumps and show a large variability of rate constants, activation energies, and rotameric populations. Site-specific activation energies vary from 6 to 38 kJ/mol. An additional mode of diffusion on a restricted arc is significant for some sites. In dry samples, the dynamics is quenched. Parameters of the motional models allow for calculations of configurational entropy and heat capacity, which, together with the rate constants, allow for observation of interplay between thermodynamic and kinetic picture of the landscape. Mutations at key phenylalanine residues at both distal (F47L&F51L) and proximal (F58L) locations to a relatively rigid side chain of L69 have a pronounced effect on alleviating the rigidity of this side chain at room temperature and demonstrate the sensitivity of the hydrophobic core environment to such perturbations.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

Bridging the Gap between Folding Simulations and Experiments: The Case of the Villin Headpiece

The increasing accuracy of molecular dynamics force fields parameters and the increasing resolution of experimental results allow one to carefully compare and complement in silico data with experimental observations. Here, we study the human villin headpiece C-terminal helical subdomain (HP35) with the recent highly optimized Amber99SB*-ILDN force field and compare the results with recent high resolution triplet-triplet energy transfer (TTET) experiments. The correct reproduction of the main structural features reveals a good agreement between experimental data and simulations.

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.

Nonequilibrium NMR Methods for Monitoring Protein and RNA Folding

The review introduces to time-resolved NMR spectroscopic investigations of the kinetics of protein and RNA folding. The description of the experimental investigations is discussed in the context of possible kinetic folding pathways showing the extent of information that can be gained from the various kinetic experiments. The review introduces to four different methods to initiate folding reactions in connection with time-resolved NMR experiments and discusses examples of refolding of the model proteinα-lactalbumin and of bistable RNAs.

Making connections between ultrafast protein folding kinetics and molecular dynamics simulations

Determining the rate of forming the truly folded conformation of ultrafast folding proteins is an important issue for both experiments and simulations. The double-norleucine mutant of the 35-residue villin subdomain is the focus of recent computer simulations with atomistic molecular dynamics because it is currently the fastest folding protein. The folding kinetics of this protein have been measured in laser temperature-jump experiments using tryptophan fluorescence as a probe of overall folding. The conclusion from the simulations, however, is that the rate determined by fluorescence is significantly larger than the rate of overall folding. We have therefore employed an independent experimental method to determine the folding rate. The decay of the tryptophan triplet-state in photoselection experiments was used to monitor the change in the unfolded population for a sequence of the villin subdomain with one amino acid difference from that of the laser temperature-jump experiments, but with almost identical equilibrium properties. Folding times obtained in a two-state analysis of the results from the two methods at denaturant concentrations varying from 1.5-6.0 M guanidinium chloride are in excellent agreement, with an average difference of only 20%. Polynomial extrapolation of all the data to zero denaturant yields a folding time of 220 (+100,-70) ns at 283 K, suggesting that under these conditions the barrier between folded and unfolded states has effectively disappeared--the so-called "downhill scenario."

Triplet-triplet energy transfer studies on conformational dynamics in peptides and a protein

Peptides and proteins are highly dynamic systems, which can adopt more or less stable conformations. The dynamics of these molecules, particularly those on the nanosecond to tens of microsecond time scale, are difficult to assess with conventional techniques. This review summarizes experiments using TTET, a technique that reports on van der Waals contact formation between a triplet donor and acceptor group, and which is sensitive in this time range. TTET allows to directly measure the chain dynamics of unstructured model peptides, i.e. large-amplitude fluctuations on the nanosecond time scale. Furthermore, contact formation can be used as irreversible probing reaction to study the kinetics of conformational equilibria. This approach enabled us to measure local α-helix folding and unfolding in helical peptides, which gave new insight into the equilibrium dynamics of this fundamental secondary structure element. TTET has also been applied to study the dynamics both in the native and unfolded state of a protein, the villin headpiece subdomain. The contact formation kinetics between different positions revealed an unlocking and local unfolding reaction in the native state of this model protein, and gave information about the chain dynamics in the unfolded state ensemble.

Dynamics of the folded and unfolded villin headpiece (HP35) measured with ultrafast 2D IR vibrational echo spectroscopy

A series of two-dimensional infrared vibrational echo experiments performed on nitrile-labeled villin headpiece [HP35-(CN)(2)] is described. HP35 is a small peptide composed of three alpha helices in the folded configuration. The dynamics of the folded HP35-(CN)(2) are compared to that of the guanidine-induced unfolded peptide, as well as the nitrile-functionalized phenylalanine (PheCN), which is used to differentiate the peptide dynamic contributions to the observables from those of the water solvent. Because the viscosity of solvent has a significant effect on fast dynamics, the viscosity of the solvent is held constant by adding glycerol. For the folded peptide, the addition of glycerol to the water solvent causes observable slowing of the peptide's dynamics. Holding the viscosity constant as GuHCl is added, the dynamics of unfolded peptide are much faster than those of the folded peptide, and they are very similar to that of PheCN. These observations indicate that the local environment of the nitrile in the unfolded peptide resembles that of PheCN, and the dynamics probed by the CN are dominated by the fluctuations of the solvent molecules, in contrast to the observations on the folded peptide.

Hydrophobic core formation and dehydration in protein folding studied by generalized-ensemble simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 A. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental phi-value analysis.

Multicolor tuning of aggregation-induced emission through substituent variation of diphenyl-o-carborane

We demonstrate multicolor tuning of aggregation-induced emission (AIE) derived from o-carborane. Both electron-donating and electron-accepting arylacetylenes underwent efficient palladium coupling reaction with bis(4-bromophenyl)-o-carborane, resulting in moderate yields. The emission spectra of these compounds span almost the entire visible spectrum (λ(max) = 452-662 nm). Study on AIE mechanism indicated that CT-based emission of o-carborane derivatives was enhanced by the restriction of molecular motions. The computational study also suggests the possibility of precise color-control of AIE through substituent variation.

Challenges in protein folding simulations: Timescale, representation, and analysis

Experimental studies of protein folding processes are frequently hampered by the fact that only low resolution structural data can be obtained with sufficient temporal resolution. Molecular dynamics simulations offer a complementary approach, providing extremely high resolution spatial and temporal data on folding processes. The effectiveness of such simulations is currently hampered by continuing questions regarding the ability of molecular dynamics force fields to reproduce the true potential energy surfaces of proteins, and ongoing difficulties with obtaining sufficient sampling to meaningfully comment on folding mechanisms. We review recent progress in the simulation of three common model systems for protein folding, and discuss how recent advances in technology and theory are allowing protein folding simulations to address their current shortcomings.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Modelling proteins: conformational sampling and reconstruction of folding kinetics

In the last decades biomolecular simulation has made tremendous inroads to help elucidate biomolecular processes in-silico. Despite enormous advances in molecular dynamics techniques and the available computational power, many problems involve long time scales and large-scale molecular rearrangements that are still difficult to sample adequately. In this review we therefore summarise recent efforts to fundamentally improve this situation by decoupling the sampling of the energy landscape from the description of the kinetics of the process. Recent years have seen the emergence of many advanced sampling techniques, which permit efficient characterisation of the relevant family of molecular conformations by dispensing with the details of the short-term kinetics of the process. Because these methods generate thermodynamic information at best, they must be complemented by techniques to reconstruct the kinetics of the process using the ensemble of relevant conformations. Here we review recent advances for both types of methods and discuss their perspectives to permit efficient and accurate modelling of large-scale conformational changes in biomolecules. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

What can solid state NMR contribute to our understanding of protein folding?

Complete understanding of the folding process that connects a structurally disordered state of a protein to an ordered, biochemically functional state requires detailed characterization of intermediate structural states with high resolution and site specificity. While the intrinsically inhomogeneous and dynamic nature of unfolded and partially folded states limits the efficacy of traditional X-ray diffraction and solution NMR in structural studies, solid state NMR methods applied to frozen solutions can circumvent the complications due to molecular motions and conformational exchange encountered in unfolded and partially folded states. Moreover, solid state NMR methods can provide both qualitative and quantitative structural information at the site-specific level, even in the presence of structural inhomogeneity. This article reviews relevant solid state NMR methods and their initial applications to protein folding studies. Using either chemical denaturation to prepare unfolded states at equilibrium or a rapid freezing apparatus to trap non-equilibrium, transient structural states on a sub-millisecond time scale, recent results demonstrate that solid state NMR can contribute essential information about folding processes that is not available from more familiar biophysical methods.