Native like structure in the unfolded state of the villin headpiece helical subdomain, an ultrafast folding protein

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔG_u°'(H₂O), of UBA(1) is 2.4 kcal mol^-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, k_f, in the absence of a denaturant of 13,000 s^-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone ¹⁵NH, ¹³CO, and ¹³Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Scaffold Hopping Transformations Using Auxiliary Restraints for Calculating Accurate Relative Binding Free Energies

In silico screening of drug-target interactions is a key part of the drug discovery process. Changes in the drug scaffold via contraction or expansion of rings, the breaking of rings, and the introduction of cyclic structures from acyclic structures are commonly applied by medicinal chemists to improve binding affinity and enhance favorable properties of candidate compounds. These processes, commonly referred to as scaffold hopping, are challenging to model computationally. Although relative binding free energy (RBFE) calculations have shown success in predicting binding affinity changes caused by perturbing R-groups attached to a common scaffold, applications of RBFE calculations to modeling scaffold hopping are relatively limited. Scaffold hopping inevitably involves breaking and forming bond interactions of quadratic functional forms, which is highly challenging. A novel method for handling ring opening/closure/contraction/expansion and linker contraction/expansion is presented here. To the best of our knowledge, RBFE calculations on linker contraction/expansion have not been previously reported. The method uses auxiliary restraints to hold the atoms at the ends of a bond in place during the breaking and forming of the bonds. The broad applicability of the method was demonstrated by examining perturbations involving small-molecule macrocycles and mutations of proline in proteins. High accuracy was obtained using the method for most of the perturbations studied. The rigor of the method was isolated from the force field by validating the method using relative and absolute hydration free energy calculations compared to standard simulation results. Unlike other methods that rely on λ-dependent functional forms for bond interactions, the method presented here can be employed using modern molecular dynamics software without modification of codes or force field functions.

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.

Fast Pairwise Approximation of Solvent Accessible Surface Area for Implicit Solvent Simulations of Proteins on CPUs and GPUs

We propose a pairwise and readily parallelizable SASA-based nonpolar solvation approach for protein simulations, inspired by our previous pairwise GB polar solvation model development. In this work, we developed a novel function to estimate the atomic and molecular SASAs of proteins, which results in comparable accuracy as the LCPO algorithm in reproducing numerical icosahedral-based SASA values. Implemented in Amber software and tested on consumer GPUs, our pwSASA method reasonably reproduces LCPO simulation results, but accelerates MD simulations up to 30 times compared to the LCPO implementation, which is greatly desirable for protein simulations facing sampling challenges. The value of incorporating the nonpolar term in implicit solvent simulations is explored on a peptide fragment containing the hydrophobic core of HP36 and evaluating thermal stability profiles of four small proteins.

CHARMM36m: an improved force field for folded and intrinsically disordered proteins

The all-atom additive CHARMM36 protein force field is widely used in molecular modeling and simulations. We present its refinement, CHARMM36m (http://mackerell.umaryland.edu/charmm_ff.shtml), with improved accuracy in generating polypeptide backbone conformational ensembles for intrinsically disordered peptides and proteins.

CHARMM36m: an improved force field for folded and intrinsically disordered proteins

The all-atom additive CHARMM36 protein force field is widely used in molecular modeling and simulations. We present its refinement, CHARMM36m (http://mackerell.umaryland.edu/charmm_ff.shtml), with improved accuracy in generating polypeptide backbone conformational ensembles for intrinsically disordered peptides and proteins.

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.

The unusual internal motion of the villin headpiece subdomain

The thermostable 36-residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side-chain motion in this protein. Here, the fast motion of methyl-bearing amino acid side chains is explored over a range of temperatures using site-resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid-state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.

A compact native 24-residue supersecondary structure derived from the villin headpiece subdomain

Many small proteins fold highly cooperatively in an all-or-none fashion and thus their native states are well protected from thermal fluctuations by an extensive network of interactions across the folded structure. Because protein structures are stabilized by local and nonlocal interactions among distal residues, dissecting individual substructures from the context of folded proteins results in large destabilization and loss of unique three-dimensional structure. Thus, mini-protein substructures can only rarely be derived from natural templates. Here, we describe a compact native 24-residues-long supersecondary structure derived from the hyperstable villin headpiece subdomain consisting of helices 2 and 3 (HP24). Using a combination of experimental techniques, including NMR and small-angle x-ray scattering, as well as all-atom replica exchange molecular-dynamics simulations, we show that a variant with stabilizing substitutions (HP24stab) forms a densely packed and compact conformation. In HP24stab, interactions between helices 2 and 3 are similar to those observed in native folded HP35, and the two helices cooperatively stabilize each other by completing the hydrophobic core lining the central part of HP35. Interestingly, even though the HP24wt fragment shows a more expanded and less structured conformation, NMR and simulations demonstrate a preference for a native-like topology. Thus, the two stabilizing residues in HP24stab shift the energy balance toward the native state, leading to a minimal folding motif.

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.

Low folding cooperativity of HP35 revealed by single-molecule force spectroscopy and molecular dynamics simulation

Some small proteins, such as HP35, fold at submicrosecond timescale with low folding cooperativity. Although these proteins have been extensively investigated, still relatively little is known about their folding mechanism. Here, using single-molecule force spectroscopy and steered molecule dynamics simulation, we study the unfolding of HP35 under external force. Our results show that HP35 unfolds at extremely low forces without a well-defined unfolding transition state. Subsequently, we probe the structure of unfolded HP35 using the persistence length obtained in the force spectroscopy. We found that the persistence length of unfolded HP35 is around 0.72 nm, >40% longer than typical unstructured proteins, suggesting that there are a significant amount of residual secondary structures in the unfolded HP35. Molecular dynamics simulation further confirmed this finding and revealed that many native contacts are preserved in HP35, even its two ends have been extended up to 8 nm. Our results therefore suggest that retaining a significant amount of secondary structures in the unfolded state of HP35 may be an efficient way to reduce the entropic cost for the formation of tertiary structure and increase the folding speed, although the folding cooperativity is compromised. Moreover, we anticipate that the methods we used in this work can be extended to the study of other proteins with complex folding behaviors and even intrinsically disordered ones.

Native state conformational heterogeneity of HP35 revealed by time-resolved FRET

The villin headpiece subdomain (HP35) has become one of the most widely used model systems in protein folding studies, due to its small size and ultrafast folding kinetics. Here, we use HP35 as a test bed to show that the fluorescence decay kinetics of an unnatural amino acid, p-cyanophenylalanine (Phe(CN)), which are modulated by a nearby quencher (e.g., tryptophan or 7-azatryptophan) through the mechanism of fluorescence resonance energy transfer (FRET), can be used to detect protein conformational heterogeneity. This method is based on the notion that protein conformations having different donor-acceptor distances and interconverting slowly compared to the fluorescence lifetime of the donor (Phe(CN)) would exhibit different donor fluorescence lifetimes. Our results provide strong evidence suggesting that the native free energy basin of HP35 is populated with conformations that differ mostly in the position and mean helicity of the C-terminal helix. This finding is consistent with several previous experimental and computational studies. Moreover, this result holds strong implications for computational investigation of the folding mechanism of HP35.

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.

Highly perturbed pKa values in the unfolded state of hen egg white lysozyme

The majority of pK(a) values in protein unfolded states are close to the amino acid model pK(a) values, thus reflecting the weak intramolecular interactions present in the unfolded ensemble of most proteins. We have carried out thermal denaturation measurements on the WT and eight mutants of HEWL from pH 1.5 to pH 11.0 to examine the unfolded state pK(a) values and the pH dependence of protein stability for this enzyme. The availability of accurate pK(a) values for the folded state of HEWL and separate measurements of mutant-induced effects on the folded state pK(a) values, allows us to estimate the pK(a) values of seven acidic residues in the unfolded state of HEWL. Asp-48 and Asp-66 display pK(a) values of 2.9 and 3.1 in our analysis, thus representing the most depressed unfolded state pK(a) values observed to date. We observe a strong correlation between the folded state pK(a) values and the unfolded state pK(a) values of HEWL, thus suggesting that the unfolded state of HEWL possesses a large degree of native state characteristics.

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.

Residual structure in unfolded proteins

The denatured state ensemble (DSE) of unfolded proteins, once considered to be well-modeled by an energetically featureless random coil, is now well-known to contain flickering elements of residual structure. The position and nature of DSE residual structure may provide clues toward deciphering the protein folding code. This review focuses on recent advances in our understanding of the nature of DSE collapse under folding conditions, the quantification of the stability of residual structure in the DSE, the determination of the location and types of residues involved in thermodynamically significant residual structure and advances in detection of long-range interactions in the DSE.

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.

Quantitative comparison of villin headpiece subdomain simulations and triplet-triplet energy transfer experiments

As the fastest folding protein, the villin headpiece (HP35) serves as an important bridge between simulation and experimental studies of protein folding. Despite the simplicity of this system, experiments continue to reveal a number of surprises, including structure in the unfolded state and complex equilibrium dynamics near the native state. Using 2.5 ms of molecular dynamics and Markov state models, we connect to current experimental results in three ways. First, we present and validate a novel method for the quantitative prediction of triplet-triplet energy transfer experiments. Second, we construct a many-state model for HP35 that is consistent with previous experiments. Finally, we predict contact-formation time traces for all 1,225 possible triplet-triplet energy transfer experiments on HP35.

Dry molten globule intermediates and the mechanism of protein unfolding

New experimental results show that either gain or loss of close packing can be observed as a discrete step in protein folding or unfolding reactions. This finding poses a significant challenge to the conventional two-state model of protein folding. Results of interest involve dry molten globule (DMG) intermediates, an expanded form of the protein that lacks appreciable solvent. When an unfolding protein expands to the DMG state, side chains unlock and gain conformational entropy, while liquid-like van der Waals interactions persist. Four unrelated proteins are now known to form DMGs as the first step of unfolding, suggesting that such an intermediate may well be commonplace in both folding and unfolding. Data from the literature show that peptide amide protons are protected in the DMG, indicating that backbone structure is intact despite loss of side-chain close packing. Other complementary evidence shows that secondary structure formation provides a major source of compaction during folding. In our model, the major free-energy barrier separating unfolded from native states usually occurs during the transition between the unfolded state and the DMG. The absence of close packing at this barrier provides an explanation for why phi-values, derived from a Brønsted-Leffler plot, depend primarily on structure at the mutational site and not on specific side-chain interactions. The conventional two-state folding model breaks down when there are DMG intermediates, a realization that has major implications for future experimental work on the mechanism of protein folding.

The folding of single domain proteins--have we reached a consensus?

Rather than stressing the most recent advances in the field, this review highlights the fundamental topics where disagreement remains and where adequate experimental data are lacking. These topics include properties of the denatured state and the role of residual structure, the nature of the fundamental steps and barriers, the extent of pathway heterogeneity and non-native interactions, recent comparisons between theory and experiment, and finally, dynamical properties of the folding reaction.

A critical assessment of putative gatekeeper interactions in the villin headpiece helical subdomain

The helical subdomain of the villin headpiece (HP36) is one of the smallest naturally occurring proteins that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extraordinary popular model system for computational, theoretical, and experimental studies of protein folding. Aromatic-proline interactions involving Trp64 and Pro62 have been proposed to play a critical role in specifying the subdomain fold by acting as gatekeeper residues. Note that the numbering corresponds to full-length headpiece. Mutation of Pro62 has been shown to lead to a protein that does not fold, but this may arise for two different reasons: The residue may make interactions that are critical for the specificity of the fold or the mutation may simply destabilize the domain. In the first case, the protein cannot fold, while in the second, the small fraction of molecules that do fold adopt the correct structure. The modest stability of the wild type prevents a critical analysis of these interactions because even moderately destabilizing mutations lead to a very small folded state population. Using a hyperstable variant of HP36, denoted DM HP36, as our new wild type, we characterized a set of mutants designed to assess the role of the putative gatekeeper interactions. Four single mutants, DM Pro62Ala, DM Trp64Leu, DM Trp64Lys, and DM Trp64Ala, and a double mutant, DM Pro62Ala Trp64Leu, were prepared. All mutants are less stable than DM HP36, but all are well folded as judged by CD and (1)H NMR. All of the mutants display sigmoidal thermal unfolding and urea-induced unfolding curves. Double-mutant cycle analysis shows that the interactions between Pro62 and Trp64 are weak but favorable. Interactions involving Pro62 and proline-aromatic interactions are, thus, not required for specifying the subdomain fold. The implications for the design and thermodynamics of miniature proteins are discussed.

An unlocking/relocking barrier in conformational fluctuations of villin headpiece subdomain

A reversible structural unlocking reaction, in which the close-packed van der Waals interactions break cooperatively, has been found for the villin headpiece subdomain (HP35) using triplet-triplet-energy transfer to monitor conformational fluctuations from equilibrium. Unlocking is associated with an unfavorable enthalpy change (DeltaH(0) = 35 +/- 4 kJ/mol) which is nearly compensated in free energy by the entropy change (DeltaS(0) = 112 +/- 20 Jxmol(-1)xK(-1)). The unlocking reaction has a time constant of about 1 mus at 5 degrees C and is enthalpy-limited with an activation energy of 32 +/- 1 kJ/mol and a large Arrhenius preexponential factor of A = 7.5 x 10(11) s(-1). In the unlocked state a fast local conformational fluctuation with a time constant of 170 ns and a low activation barrier of 17 +/- 1 kJ/mol leads to unfolding of the C-terminal helix and to its undocking from the core. On a much slower time scale, global unfolding occurs from the unlocked state. These results suggest that native protein structures are locked into conformations with low amplitude motions. Large scale motions and global unfolding require an initial structural unlocking step leading to a state with properties of a dry molten globule. The experiments additionally yielded information on the dynamics of loop formation between different positions in unfolded HP35. Comparison of the results with dynamics in unstructured model peptides indicates slightly decelerated kinetics of local loop formation in the C-terminal region which points at residual, nonrandom structure. Dynamics of long-range loop formation, in contrast, are not influenced by residual structure, which argues against unfolded state properties as molecular origin for ultrafast folding of HP35.