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

Toward a Benchmark for Markov State Models: The Folding of HP35

Adopting a 300 μs long MD trajectory of the folding of villin headpiece (HP35) by D. E. Shaw Research, we recently constructed a Markov state model (MSM) based on inter-residue contacts. The model reproduces the folding time and predicts that the native basin and unfolded region consist of metastable substates that are structurally well-characterized. Recognizing the need to establish well-defined benchmark problems, we study to what extent and in what sense this MSM can be employed as a reference model. Hence, we test the robustness of the MSM by comparing it to models that use alternative combinations of features, dimensionality reduction methods, and clustering schemes. The study suggests some main characteristics of the folding of HP35 that should be reproduced by other competitive models. Moreover, the discussion reveals which parts of the MSM workflow matter most for the considered problem and illustrates the promises and pitfalls of state-based models for the interpretation of biomolecular simulations.

Size-and-Shape Space Gaussian Mixture Models for Structural Clustering of Molecular Dynamics Trajectories

Determining the optimal number and identity of structural clusters from an ensemble of molecular configurations continues to be a challenge. Recent structural clustering methods have focused on the use of internal coordinates due to the innate rotational and translational invariance of these features. The vast number of possible internal coordinates necessitates a feature space supervision step to make clustering tractable but yields a protocol that can be system type-specific. Particle positions offer an appealing alternative to internal coordinates but suffer from a lack of rotational and translational invariance, as well as a perceived insensitivity to regions of structural dissimilarity. Here, we present a method, denoted shape-GMM, that overcomes the shortcomings of particle positions using a weighted maximum likelihood alignment procedure. This alignment strategy is then built into an expectation maximization Gaussian mixture model (GMM) procedure to capture metastable states in the free-energy landscape. The resulting algorithm distinguishes between a variety of different structures, including those indistinguishable by root-mean-square displacement and pairwise distances, as demonstrated on several model systems. Shape-GMM results on an extensive simulation of the fast-folding HP35 Nle/Nle mutant protein support a four-state folding/unfolding mechanism, which is consistent with previous experimental results and provides kinetic details comparable to previous state-of-the art clustering approaches, as measured by the VAMP-2 score. Currently, training of shape-GMMs is recommended for systems (or subsystems) that can be represented by ≲200 particles and ≲100k configurations to estimate high-dimensional covariance matrices and balance computational expense. Once a shape-GMM is trained, it can be used to predict the cluster identities of millions of configurations.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Exposing the distinctive modular behavior of β-strands and α-helices in folded proteins

Although folded proteins are commonly depicted as simplistic combinations of β-strands and α-helices, the actual properties and functions of these secondary-structure elements in their native contexts are just partly understood. The principal reason is that the behavior of individual β- and α-elements is obscured by the global folding cooperativity. In this study, we have circumvented this problem by designing frustrated variants of the mixed α/β-protein S6, which allow the structural behavior of individual β-strands and α-helices to be targeted selectively by stopped-flow kinetics, X-ray crystallography, and solution-state NMR. Essentially, our approach is based on provoking intramolecular "domain swap." The results show that the α- and β-elements have quite different characteristics: The swaps of β-strands proceed via global unfolding, whereas the α-helices are free to swap locally in the native basin. Moreover, the α-helices tend to hybridize and to promote protein association by gliding over to neighboring molecules. This difference in structural behavior follows directly from hydrogen-bonding restrictions and suggests that the protein secondary structure defines not only tertiary geometry, but also maintains control in function and structural evolution. Finally, our alternative approach to protein folding and native-state dynamics presents a generally applicable strategy for in silico design of protein models that are computationally testable in the microsecond-millisecond regime.

PET and FRET utility of an amino acid pair: tryptophan and 4-cyanotryptophan

Methods based on fluorescence resonance energy transfer (FRET) and photo-induced electron transfer (PET) are widely used in the biological sciences, employing mostly dye-based FRET and PET pairs. While very useful and important, dye-based reporters are not always applicable without concern, for example, in cases where the fluorophore size needs to be minimized. Therefore, development and characterization of smaller, ideally amino acid-based PET and FRET pairs will expand the biological spectroscopy toolbox to enable new applications. Herein, we show that, depending on the excitation wavelength, tryptophan and 4-cyanotrptophan can interact with each other via the mechanism of either energy or electron transfer, hence constituting a dual FRET and PET pair. The biological utility of this amino acid pair is further demonstrated by applying it to study the end-to-end collision rate of a short peptide, the mode of interaction between a ligand and BSA, and the activity of a protease.

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.

Comparison of biological chromophores: photophysical properties of cyanophenylalanine derivatives

Within this work, the family of cyanophenylalanine spectroscopic reporters is extended by showing the ortho and meta derivatives have intrinsic photophysical properties that are useful for studies of protein structure and dynamics. The molar absorptivities of 2-cyanophenylalanine and 3-cyanophenylalanine are shown to be comparable to that of 4-cyanophenylalanine with similar spectral features in their absorbance and emission profiles, demonstrating that these probes can be utilized interchangeably. The fluorescence quantum yields are also on the same scale as commonly used fluorophores in peptides and proteins, tyrosine and tryptophan. These new cyano-fluorophores can be paired with either 4-cyanophenylalanine or tryptophan to capture distances in peptide structure through Förster resonance energy transfer. Additionally, the spectroscopic properties of these chromophores can report the local solvent environment via changes in fluorescence emission intensity as a result of hydrogen bonding and/or hydration. A decrease in the quantum yield is also observed in basic environments due to photoinduced electron transfer from a deprotonated amine in the free PheCN species and at the N-terminus of a short peptide, providing an avenue to detect pH in biological systems. Our results show the potential of these probes, 2-cyanophenylalanine and 3-cyanophenylalanine, to be incorporated into a single peptide chain, either individually or in tandem with 4-cyanophenylalanine, tryptophan, or tyrosine, in order to obtain information about peptide structure and dynamics.

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.

Minimalist IR and fluorescence probes of protein function

Spectroscopic studies of small proteins and peptides, especially those requiring fine spatial and/or temporal resolution, demand synthetic probes that confer the minimal possible steric and functional change on the native properties. Here we review the recent progress in development of minimally disruptive probes for fluorescence and infrared spectroscopies, as well as the methods to efficiently incorporate them into proteins. Advances in spectroscopy on the one hand result in high specialization of synthetic probes for a particular purpose, but on the other hand allow for the same probes be used for different techniques to gather complementary biochemical information.

Infrared and Fluorescence Assessment of Protein Dynamics: From Folding to Function

While folding or performing functions, a protein can sample a rich set of conformational space. However, experimentally capturing all of the important motions with sufficient detail to allow a mechanistic description of their dynamics is nontrivial since such conformational events often occur over a wide range of time and length scales. Therefore, many methods have been employed to assess protein conformational dynamics, and depending on the nature of the conformational transition in question, some may be more advantageous than others. Herein, we describe our recent efforts, and also those of others, wherever appropriate, to use infrared- and fluorescence-based techniques to interrogate protein folding and functional dynamics. Specifically, we focus on discussing how to use extrinsic spectroscopic probes to enhance the structural resolution of these techniques and how to exploit various cross-linking strategies to acquire dynamic and mechanistic information that was previously difficult to attain.

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.

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.

Sensing pH via p-cyanophenylalanine fluorescence: Application to determine peptide pKa and membrane penetration kinetics

We expand the spectroscopic utility of a well-known infrared and fluorescence probe, p-cyanophenylalanine, by showing that it can also serve as a pH sensor. This new application is based on the notion that the fluorescence quantum yield of this unnatural amino acid, when placed at or near the N-terminal end of a polypeptide, depends on the protonation status of the N-terminal amino group of the peptide. Using this pH sensor, we are able to determine the N-terminal pKa values of nine tripeptides and also the membrane penetration kinetics of a cell-penetrating peptide. Taken together, these examples demonstrate the applicability of using this unnatural amino acid fluorophore to study pH-dependent biological processes or events that accompany a pH change.

New aspects of the structure and mode of action of the human cathelicidin LL-37 revealed by the intrinsic probe p-cyanophenylalanine

The human cathelicidin peptide LL-37 is an important effector of our innate immune system and contributes to host defence with direct antimicrobial activity and immunomodulatory properties, and by stimulating wound healing. Its sequence has evolved to confer specific structural characteristics that strongly affect these biological activities, and differentiate it from orthologues of other primate species. In the present paper we report a detailed study of the folding and self-assembly of this peptide in comparison with rhesus monkey peptide RL-37, taking into account the different stages of its trajectory from bulk solution to contact with, and insertion into, biological membranes. Phenylalanine residues in different positions throughout the native sequences of LL-37 and RL-37 were systematically replaced with the non-invasive fluorescent and IR probe p-cyanophenylalanine. Steady-state and time-resolved fluorescence studies showed that LL-37, in contrast to RL-37, forms oligomers with a loose hydrophobic core in physiological solutions, which persist in the presence of biological membranes. Fourier transform IR and surface plasmon resonance studies also indicated different modes of interaction for LL-37 and RL-37 with anionic and neutral membranes. This correlated with a distinctly different mode of bacterial membrane permeabilization, as determined using a flow cytometric method involving impermeant fluorescent dyes linked to polymers of defined sizes.

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.

Site-specific dynamics of amyloid formation and fibrillar configuration of Aβ1–23 using an unnatural amino acid

Distinct local dynamics of Aβ_1–23 amyloid formation are characterized using an unnatural amino acid p-cyanophenylalanine as a spectroscopic probe.

Molecular dynamics simulation of configurational ensembles compatible with experimental FRET efficiency data through a restraint on instantaneous FRET efficiencies

Förster resonance energy transfer (FRET) measurements are widely used to investigate (bio)molecular interactions or/and association. FRET efficiencies, the primary data obtained from this method, give, in combination with the common assumption of isotropic chromophore orientation, detailed insight into the lengthscale of molecular phenomena. This study illustrates the application of a FRET efficiency restraint during classical atomistic molecular dynamics simulations of a mutant mastoparan X peptide in either water or 7 M aqueous urea. The restraint forces acting on the donor and acceptor chromophores ensure that the sampled peptide configurational ensemble satisfies the experimental primary data by modifying interchromophore separation and chromophore transition dipole moment orientations. By means of a conformational cluster analysis, it is seen that indeed different configurational ensembles may be sampled without and with application of the restraint. In particular, while the FRET efficiency and interchromophore distances monitored in an unrestrained simulation may differ from the experimentally-determined values, they can be brought in agreement with experimental data through usage of the FRET efficiency restraining potential. Furthermore, the present results suggest that the assumption of isotropic chromophore orientation is not always justified. The FRET efficiency restraint allows the generation of configurational ensembles that may not be accessible with unrestrained simulations, and thereby supports a meaningful interpretation of experimental FRET results in terms of the underlying molecular degrees of freedom. Thus, it offers an additional tool to connect the realms of computer and wet-lab experimentation.

Ultrafast folding kinetics and cooperativity of villin headpiece in single-molecule force spectroscopy

In this study we expand the accessible dynamic range of single-molecule force spectroscopy by optical tweezers to the microsecond range by fast sampling. We are able to investigate a single molecule for up to 15 min and with 300-kHz bandwidth as the protein undergoes tens of millions of folding/unfolding transitions. Using equilibrium analysis and autocorrelation analysis of the time traces, the full energetics as well as real-time kinetics of the ultrafast folding of villin headpiece 35 and a stable asparagine 68 alanine/lysine 70 methionine variant can be measured directly. We also performed Brownian dynamics simulations of the response of the bead-DNA system to protein-folding fluctuations. All key features of the force-dependent deflection fluctuations could be reproduced: SD, skewness, and autocorrelation function. Our measurements reveal a difference in folding pathway and cooperativity between wild-type and stable variant of headpiece 35. Autocorrelation force spectroscopy pushes the time resolution of single-molecule force spectroscopy to ∼10 µs thus approaching the timescales accessible for all atom molecular dynamics simulations.

High-Resolution Free-Energy Landscape Analysis of α-Helical Protein Folding: HP35 and Its Double Mutant

The free-energy landscape can provide a quantitative description of folding dynamics, if determined as a function of an optimally chosen reaction coordinate. Here, we construct the optimal coordinate and the associated free-energy profile for all-helical proteins HP35 and its norleucine (Nle/Nle) double mutant, based on realistic equilibrium folding simulations [Piana et al. Proc. Natl. Acad. Sci. U.S.A.2012, 109, 17845]. From the obtained profiles, we directly determine such basic properties of folding dynamics as the configurations of the minima and transition states (TS), the formation of secondary structure and hydrophobic core during the folding process, the value of the pre-exponential factor and its relation to the transition path times, the relation between the autocorrelation times in TS and minima. We also present an investigation of the accuracy of the pre-exponential factor estimation based on the transition-path times. Four different estimations of the pre-exponential factor for both proteins give k₀^–1 values of approximately a few tens of nanoseconds. Our analysis gives detailed information about folding of the proteins and can serve as a rigorous common language for extensive comparison between experiment and simulation.

Allosteric activation of the Par-6 PDZ via a partial unfolding transition

Proteins exist in a delicate balance between the native and unfolded states, where thermodynamic stability may be sacrificed to attain the flexibility required for efficient catalysis, binding, or allosteric control. Partition-defective 6 (Par-6) regulates the Par polarity complex by transmitting a GTPase signal through the Cdc42/Rac interaction binding PSD-95/Dlg/ZO-1 (CRIB-PDZ) module that alters PDZ ligand binding. Allosteric activation of the PDZ is achieved by local rearrangement of the L164 and K165 side chains to stabilize the interdomain CRIB:PDZ interface and reposition a conserved element of the ligand binding pocket. However, microsecond to millisecond dynamics measurements revealed that L164/K165 exchange requires a larger rearrangement than expected. The margin of thermodynamic stability for the PDZ domain is modest (∼3 kcal/mol) and further reduced by transient interactions with the disordered CRIB domain. Measurements of local structural stability revealed that tertiary contacts within the PDZ are disrupted by a partial unfolding transition that enables interconversion of the L/K switch. The unexpected participation of partial PDZ unfolding in the allosteric mechanism of Par-6 suggests that native-state unfolding may be essential for the function of other marginally stable proteins.

Quenching of p-Cyanophenylalanine Fluorescence by Various Anions

To expand the spectroscopic utility of the non-natural amino acid p-cyanophenylalanine (Phe_CN), we examine the quenching efficiencies of a series of commonly encountered anions toward its fluorescence. We find that iodide exhibits an unusually large Stern-Volmer quenching constant, making it a convenient choice in Phe_CN fluorescence quenching studies. Indeed, using the villin headpiece subdomain as a testbed we demonstrate that iodide quenching of Phe_CN fluorescence offers a convenient means to reveal protein conformational heterogeneity. Furthermore, we show that the amino group of Phe_CN strongly quenches its fluorescence, suggesting that Phe_CN could be used as a local pH sensor.