Solution structures of the C-terminal headpiece subdomains of human villin and advillin, evaluation of headpiece F-actin-binding requirements

Headpiece (HP) is a 76-residue F-actin-binding module at the C terminus of many cytoskeletal proteins. Its 35-residue C-terminal subdomain is one of the smallest known motifs capable of autonomously adopting a stable, folded structure in the absence of any disulfide bridges, metal ligands, or unnatural amino acids. We report the three-dimensional solution structures of the C-terminal headpiece subdomains of human villin (HVcHP) and human advillin (HAcHP), determined by two-dimensional 1H-NMR. They represent the second and third structures of such C-terminal headpiece subdomains to be elucidated so far. A comparison with the structure of the chicken villin C-terminal subdomain reveals a high structural conservation. Both C-terminal subdomains bind specifically to F-actin. Mutagenesis is used to demonstrate the involvement of Trp 64 in the F-actin-binding surface. The latter residue is part of a conserved structural feature, in which the surface-exposed indole ring is stacked on the proline and lysine side chain embedded in a PXWK sequence motif. On the basis of the structural and mutational data concerning Trp 64 reported here, the results of a cysteine-scanning mutagenesis study of full headpiece, and a phage display mutational study of the 69-74 fragment, we propose a modification of the model, elaborated by Vardar and coworkers, for the binding of headpiece to F-actin.     

Crystal structure of a pH-stabilized mutant of villin headpiece

Villin-type headpiece domains are compact F-actin-binding motifs that have been used extensively as a model system to investigate protein folding by both experimental and computational methods. Villin headpiece (HP67) harbors a highly helical, thermostable, and autonomously folding subdomain in the C terminus (HP35), and because of this feature, HP67 is usually considered to be composed of a N- and C-terminal subdomain. Unlike the C-terminal subdomain, the N-terminal subdomain consists mainly of loops and turns, and the folding is dependent upon the presence of the C-terminal subdomain. The pH sensitivity of this subdomain is thought to arise from, at least partially, protonation of H41 buried in the hydrophobic core. Substitution of this histidine with tyrosine, another permissive residue at this position for naturally occurring sequences, increases not only the pH stability of HP67 but also the thermal stability and the cooperativity of thermal unfolding over a wide pH range (0.9-7.5). The crystal structures of wild-type HP67 and the H41Y mutant, determined under the same conditions, indicate that the H41Y substitution causes only localized rearrangement around the mutated residue. The F-actin-binding motif remains essentially the same after the mutation, accounting for the negligible effect of the mutation on F-actin affinity. The hydrogen bond formed between the imidazole ring of H41 and the backbone carbonyl of E14 of HP67 is eliminated by the H41Y mutation, which renders the extreme N terminus of H41Y more mobile; the hydrogen bond formed between the imidazole ring of H41 and the backbone nitrogen of D34 is replaced with that between the hydroxyl group of Y41 and the backbone nitrogen of D34 after the H41Y substitution. The increased hydrophobicity of tyrosine compensates for the loss of hydrogen bonds in the extreme N terminus and accounts for the increased stability and cooperativity of the H41Y mutant.     

Folding network of villin headpiece subdomain

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 μs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.     

Multistate folding of the villin headpiece domain

The villin headpiece (HP67) is a 67 residue, monomeric protein derived from the C-terminal domain of villin. Wild-type HP67 (WT HP67) is the smallest fragment of villin that retains strong in vitro actin-binding activity. WT HP67 is made up of two subdomains, which form a tightly packed interface. The C-terminal subdomain of WT HP67, denoted HP35, is rich in helical structure, folds in isolation, and has been widely used as a model system for folding studies. In contrast, very little is known about the folding of the intact villin headpiece domain. Here, NMR, CD and H/2H amide exchange measurements are used to follow the pH, thermal and urea-induced unfolding of WT HP67 and a mutant (HP67 H41Y) in which a buried conserved histidine in the N-terminal subdomain, His41, has been mutated to Tyr. Although most small proteins display two-state equilibrium unfolding, the results presented here demonstrate that unfolding of the villin headpiece is a multistate process. The presence of a folded N-terminal subdomain is shown to stabilize the C-terminal subdomain, increasing the midpoints of the thermal and urea-induced unfolding transitions and increasing protection factors for H/2H exchange. Histidine 41 has been shown to act as a pH-dependent switch in wild-type HP67: the N-terminal subdomain is unfolded when His41 is protonated, while the C-terminal subdomain remains folded irrespective of the protonation state of His41. Mutation of His41 to Tyr eliminates the segmental pH-dependent unfolding of the headpiece. The mutation stabilizes both domains, but folding is still multistate, indicating that His41 is not solely responsible for the unusual equilibrium unfolding behavior of villin headpiece domain.     

Free-energy landscape of the villin headpiece in an all-atom force field

We investigate the landscape of the internal free-energy of the 36 amino acid villin headpiece with a modified basin hopping method in the all-atom force field PFF01, which was previously used to predictively fold several helical proteins with atomic resolution. We identify near native conformations of the protein as the global optimum of the force field. More than half of the twenty best simulations started from random initial conditions converge to the folding funnel of the native conformation, but several competing low-energy metastable conformations were observed. From 76,000 independently generated conformations we derived a decoy tree which illustrates the topological structure of the entire low-energy part of the free-energy landscape and characterizes the ensemble of metastable conformations. These emerge as similar in secondary content, but differ in tertiary arrangement.     

Villin-type headpiece domains show a wide range of F-actin-binding affinities

The villin-type "headpiece" domain is a modular motif found at the extreme C-terminus of larger "core" domains in over 25 cytoskeletal proteins in plants and animals. Although headpiece is classified as an F-actin-binding domain, it has been suggested that some expressed fusion-proteins containing headpiece may lack F-actin-binding in vivo. To determine the intrinsic F-actin affinity of headpiece domains, we quantified the F-actin affinity of seven headpiece domains and three N-terminal truncations, under identical in vitro conditions. The constructs are folded and adopt the native headpiece structure. However, they show a wide range of affinities that can be grouped into high, low, and nonspecific-binding categories. Computer models of the structure and charged surface potential of these headpiece domains suggest features important for high F-actin affinity. We conclude that not all headpiece domains are intrinsically F-actin-binding motifs, and suggest that the surface charge distribution may be an important element for F-actin recognition.     

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.     

NMR characterization of a peptide model provides evidence for significant structure in the unfolded state of the villin headpiece helical subdomain

The villin headpiece subdomain (HP36) is the smallest naturally occurring protein that folds cooperatively. The protein folds on a microsecond time scale. Its small size and very rapid folding have made it a popular target for biophysical studies of protein folding. Temperature-dependent one-dimensional (1D) NMR studies of the full-length protein together with CD and 1D NMR studies of the 21-residue peptide fragment (HP21) derived from HP36 have shown that there is significant structure in the unfolded state of HP36 and have demonstrated that HP21 is a good model of these interactions. Here, we characterized the model peptide HP21 in detail by two-dimensional NMR. Strongly upfield shifted C(alpha) protons, the magnitude of the 3J(NH,alpha) coupling constants, and the pattern of backbone-backbone and backbone-side chain NOEs indicate that the ensemble of structures populated by HP21 contains alpha-helical structure and native as well as non-native hydrophobic contacts. The hydrogen-bonded secondary structure inferred from the NOEs is, however, not sufficient to confer significant protection against amide H-D exchange. These studies indicate that there is significant secondary structure and hydrophobic clustering in the unfolded state of HP36. The implications for the folding of HP36 are discussed.     

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.     

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

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

On unsatisfied hydrogen bonds in the N-terminal subdomain of villin headpiece

Villin headpiece is a small autonomously folding protein that has emerged as a model system for understanding the fundamental tenets governing protein folding. In this communication, we employ NMR and X-ray crystallography to characterize a point mutant, H41F, which retains actin-binding activity, is more thermostable but, interestingly, does not exhibit the partially folded intermediate observed of either wild-type or other similar point mutants.     

Slow motions in chicken villin headpiece subdomain probed by cross-correlated NMR relaxation of amide NH bonds in successive residues

The villin headpiece subdomain (HP36) is a widely used system for protein-folding studies. Nuclear magnetic resonance cross-correlated relaxation rates arising from correlated fluctuations of two N-H^N dipole-dipole interactions involving successive residues were measured at two temperatures at which HP36 is at least 99% folded. The experiment revealed the presence of motions slower than overall tumbling of the molecule. Based on the theoretical analysis of the spectral densities we show that the structural and dynamic contributions to the experimental cross-correlated relaxation rate can be separated under certain conditions. As a result, dynamic cross-correlated order parameters describing slow microsecond-to-millisecond motions of N-H bonds in neighboring residues can be introduced for any extent of correlations in the fluctuations of the two bond vectors. These dynamic cross-correlated order parameters have been extracted for HP36. The comparison of their values at two different temperatures indicates that when the temperature is raised, slow motions increase in amplitude. The increased amplitude of these fluctuations may reflect the presence of processes directly preceding the unfolding of the protein.     

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: ¹³C′ longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, ¹³C′/¹³C′-¹³C^α CSA/dipolar and ¹³C′/¹³C′–¹⁵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.     

The unfolded state of the villin headpiece helical subdomain: computational studies of the role of locally stabilized structure

The 36 residue villin headpiece helical subdomain (HP36) is one of the fastest cooperatively folding proteins, folding on the microsecond timescale. HP36's simple three helix topology, fast folding and small size have made it an attractive model system for computational and experimental studies of protein folding. Recent experimental studies have explored the denatured state of HP36 using fragment analysis coupled with relatively low-resolution spectroscopic techniques. These studies have shown that there is apparently only a small tendency to form locally stabilized secondary structure. Here, we complement the experimental studies by using replica exchange molecular dynamics with explicit solvent to investigate the structural features of these peptide models of unfolded HP36. To ensure convergence, two sets of simulations for each fragment were performed with different initial structures, and simulations were continued until these generated very similar final ensembles. These simulations reveal low populations of native-like structure and early folding events that cannot be resolved by experiment. For each fragment, calculated J-coupling constants and helical propensities are in good agreement with experimental trends. HP-1, corresponding to residues 41 to 53 and including the first alpha-helix, contains the highest helical population. HP-3, corresponding to residues 62 through 75 and including the third alpha-helix, contains a small population of helical turn residing at the N terminus while HP-2, corresponding to residues 52 through 61 and including the second alpha-helix, formed little to no structure in isolation. Overall, HP-1 was the only fragment to adopt a native-like conformation, but the low population suggests that formation of significant structure only occurs after formation of specific tertiary interactions.     

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.     

High-resolution crystal structures of villin headpiece and mutants with reduced F-actin binding activity

Villin-type headpiece domains are approximately 70 amino acid modular motifs found at the C terminus of a variety of actin cytoskeleton-associated proteins. The headpiece domain of villin, a protein found in the actin bundles of the brush border epithelium, is of interest both as a compact F-actin binding domain and as a model folded protein. We have determined the high-resolution crystal structures of chicken villin headpiece (HP67) at 1.4 A resolution as well as two mutants, R37A and W64Y, at 1.45 and 1.5 A resolution, respectively. Replacement of R37 causes a 5-fold reduction in F-actin binding affinity in sedimentation assays. Replacement of W64 results in a much more drastic reduction in F-actin binding affinity without significant changes in headpiece structure or stability. The detailed comparison of these crystal structures with each other and to our previously determined NMR structures of HP67 and the 35-residue autonomously folding subdomain in villin headpiece, HP35, provides the details of the headpiece fold and further defines the F-actin binding site of villin-type headpiece domains.     

31P NMR spectra of an oligodeoxyribonucleotide duplex lac operator-repressor headpiece complex

The interaction of a symmetric lac operator duplex, d(TGTGAGCGCTCACA)2, with the N-terminal 56-residue headpiece fragment of the lac repressor protein was monitored by 31P NMR spectroscopy. The changes in the 31P chemical shifts upon addition of the headpiece demonstrated an end point of two headpiece fragments per symmetric 14-mer duplex with each headpiece binding to the T1pG2pT3pG4pA5 ends of the duplex. The specific phosphate 31P perturbations observed are consistent with those residues implicated in protein binding by previous NMR, molecular biological, and biochemical techniques. Upon complexation, the 31P signals of phosphates G2-A5 showed upfield or downfield shifts (less than 0.2 ppm) while most other residues were unperturbed. The interactions were dependent on ionic strength. The 31P NMR data provide direct evidence for predominant recognition of the 5' strand of the 5'-TGTGA/3'-ACACT binding site.     

The role of aromatic residues in the hydrophobic core of the villin headpiece subdomain

Small autonomously folding proteins are of interest as model systems to study protein folding, as the same molecule can be used for both experimental and computational approaches. The question remains as to how well these minimized peptide model systems represent larger native proteins. For example, is the core of a minimized protein tolerant to mutation like larger proteins are? Also, do minimized proteins use special strategies for specifying and stabilizing their folded structure? Here we examine these questions in the 35‐residue autonomously folding villin headpiece subdomain (VHP subdomain). Specifically, we focus on a cluster of three conserved phenylalanine (F) residues F47, F51, and F58, that form most of the hydrophobic core. These three residues are oriented such that they may provide stabilizing aromatic–aromatic interactions that could be critical for specifying the fold. Circular dichroism and 1D‐NMR spectroscopy show that point mutations that individually replace any of these three residues with leucine were destabilized, but retained the native VHP subdomain fold. In pair‐wise replacements, the double mutant that retains F58 can adopt the native fold, while the two double mutants that lack F58 cannot. The folding of the double mutant that retains F58 demonstrates that aromatic–aromatic interactions within the aromatic cluster are not essential for specifying the VHP subdomain fold. The ability of the VHP subdomain to tolerate mutations within its hydrophobic core indicates that the information specifying the three dimensional structure is distributed throughout the sequence, as observed in larger proteins. Thus, the VHP subdomain is a legitimate model for larger, native proteins.