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

HP35 Protein in the Mesopore of MIL-101(Cr) MOF: A Model to Study Cotranslocational Unfolding

The immobilization of enzymes in metal-organic framework (MOF) cages is important in biotechnology. In this context, the mechanism of translocation of proteins through the cavities of the MOF and the roles played by confinement and MOF chemistry in giving rise to stable protein intermediates that are otherwise transiently populated in the physiological environment are important questions to be addressed. These unexplored aspects are examined with villin headpiece (HP35) as a model protein confined within a mesopore of MIL-101(Cr) using molecular dynamics simulations. At equilibrium, the protein is located farther from the center of the cavity and closer to the MOF surface. Molecular interactions with the MOF partially unfold helix-1 at its N-terminus. Umbrella sampling simulations inform the range of conformations that HP35 undertakes during translocation from one cavity to another and associated changes in free energy. Relative to its equilibrium state within the cavity, the free energy barrier for the unfolded protein at the cage window is estimated to be 16 kcal/mol. This study of MOF-based protein conformation can also be a general approach to observing intermediates in folding-unfolding pathways.

Evolutionary conservation of amino acids contributing to the protein folding transition state

The question of whether amino acids critical to protein folding kinetics are evolutionarily conserved has been investigated intensively in the past, but no consensus has yet been reached. Recently, we have demonstrated that the transition state, dictating folding kinetics, is characterized as the state of maximum dynamic cooperativity, i.e., the state of maximum correlations between amino acid contact formations. Here, we investigate the evolutionary conservation of those amino acids contributing significantly to the dynamic cooperativity. We find a strong indication of a new kind of relationship-necessary but not sufficient causality-between the evolutionary conservation and the dynamic cooperativity: larger contributions to the dynamic cooperativity arise from more conserved residues, but not vice versa. This holds for all the protein systems for which long folding simulation trajectories are available. To our knowledge, this is the first systematic demonstration of any kind of evolutionary conservation of amino acids relevant to folding kinetics.

Targeting AVIL, a New Cytoskeleton Regulator in Glioblastoma

Glioblastoma (GBM) is the most common adult neural malignancy and the deadliest. The standard of care is optimal, safe, cytoreductive surgery followed by combined radiation therapy and alkylating chemotherapy with temozolomide. Recurrence is common and therapeutic options in the recurrent setting are limited. The dismal prognosis of GBM has led to novel treatments being a serious roadblock in the field, with most new treatments failing to show efficacy. Targeted therapies have shown some success in many cancers, but GBM remains one of the most difficult to treat, especially in recurrence. New chemotherapeutic directions need to be explored, possibly expanding the targeted chemotherapy spectrum in previously unforeseen ways. In this perspective paper, we will explain why AVIL, an actin-binding protein recently found to be overexpressed in GBM and a driving force for GBM, could prove versatile in the fight against cancer. By looking at AVIL and its potential to regulate FOXM1 and LIN28B, we will be able to highlight a way to improve outcomes for GBM patients who normally have very little hope.

F-actin-bundling sites are conserved in proteins with villin-type headpiece domains

Villin is a major actin-bundling protein that assembles the brush border of intestinal and renal epithelial cells. The villin "headpiece" domain and the actin-binding residues within it regulate its actin-bundling function. Substantial experimental and theoretical information about the three-dimensional structure of the isolated villin headpiece, including a description of the actin-binding residues within the headpiece, is available. Despite that, the actin-bundling site in the full-length (FL) villin protein remains unidentified. We used this existing villin headpiece nuclear magnetic resonance data and performed mutational analysis and functional assays to identify the actin-bundling site in FL human villin protein. By careful evaluation of these conserved actin-binding residues in human advillin protein, we demonstrate their functional significance in the over 30 proteins that contain a villin-type headpiece domain. Our study is the first that combines the available structural data on villin headpiece with functional assays to identify the actin-binding residues in FL villin that regulate its filament-bundling activity. Our findings could have wider implications for other actin-bundling proteins that contain a villin-type headpiece domain.

Steroidal glycoalkaloids from Solanum nigrum target cytoskeletal proteins: an in silico analysis

Background

Solanum nigrum (black nightshade; S. nigrum), a member of family Solanaceae, has been endowed with a heterogeneous array of secondary metabolites of which the steroidal glycoalkaloids (SGAs) and steroidal saponins (SS) have vast potential to serve as anticancer agents. Since there has been much controversy regarding safety of use of glycoalkaloids as anticancer agents, this area has remained more or less unexplored. Cytoskeletal proteins like actin play an important role in maintaining cell shape, synchronizing cell division, cell motility, etc. and along with their accessory proteins may also serve as important therapeutic targets for potential anticancer candidates. In the present study, glycoalkaloids and saponins from S. nigrum were screened for their interaction and binding affinity to cytoskeletal proteins, using molecular docking.

Methods

Bioactivity score and Prediction of Activity Spectra for Substances (PASS) analysis were performed using softwares Molinspiration and Osiris Data Explorer respectively, to assess the feasibility of selected phytoconstituents as potential drug candidates. The results were compared with two standard reference drugs doxorubicin hydrochloride (anticancer) and tetracycline (antibiotic). Multivariate data obtained were analyzed using principal component analysis (PCA).

Results

Docking analysis revealed that the binding affinities of the phytoconstituents towards the target cytoskeletal proteins decreased in the order coronin>villin>ezrin>vimentin>gelsolin>thymosin>cofilin. Glycoalkaloid solasonine displayed the greatest binding affinity towards the target proteins followed by alpha-solanine whereas amongst the saponins, nigrumnin-I showed maximum binding affinity. PASS Analysis of the selected phytoconstituents revealed 1 to 3 violations of Lipinski’s parameters indicating the need for modification of their structure-activity relationship (SAR) for improvement of their bioactivity and bioavailability. Glycoalkaloids and saponins all had bioactivity scores between −5.0 and 0.0 with respect to various receptor proteins and target enzymes. Solanidine, solasodine and solamargine had positive values of druglikeness which indicated that these compounds have the potential for development into future anticancer drugs. Toxicity potential evaluation revealed that glycoalkaloids and saponins had no toxicity, tumorigenicity or irritant effect(s). SAR analysis revealed that the number, type and location of sugar or the substitution of hydroxyl group on alkaloid backbone had an effect on the activity and that the presence of α-L-rhamnopyranose sugar at C-2 was critical for a compound to exhibit anticancer activity.

Conclusion

The present study revealed some cytoskeletal target(s) for S. nigrum phytoconstituents by docking analysis that have not been previously reported and thus warrant further investigations both in vitro and in vivo.

Hydrolytic zinc metallopeptides using a computational multi-state design approach

Combination of multi-state design and long-timescale conformational dynamics as a powerful strategy to obtain metalloenzymes.

Advillin acts upstream of phospholipase C ϵ1 in steroid-resistant nephrotic syndrome

Steroid-resistant nephrotic syndrome (SRNS) is a frequent cause of chronic kidney disease. Here, we identified recessive mutations in the gene encoding the actin-binding protein advillin (AVIL) in 3 unrelated families with SRNS. While all AVIL mutations resulted in a marked loss of its actin-bundling ability, truncation of AVIL also disrupted colocalization with F-actin, thereby leading to impaired actin binding and severing. Additionally, AVIL colocalized and interacted with the phospholipase enzyme PLCE1 and with the ARP2/3 actin-modulating complex. Knockdown of AVIL in human podocytes reduced actin stress fibers at the cell periphery, prevented recruitment of PLCE1 to the ARP3-rich lamellipodia, blocked EGF-induced generation of diacylglycerol (DAG) by PLCE1, and attenuated the podocyte migration rate (PMR). These effects were reversed by overexpression of WT AVIL but not by overexpression of any of the 3 patient-derived AVIL mutants. The PMR was increased by overexpression of WT Avil or PLCE1, or by EGF stimulation; however, this increased PMR was ameliorated by inhibition of the ARP2/3 complex, indicating that ARP-dependent lamellipodia formation occurs downstream of AVIL and PLCE1 function. Together, these results delineate a comprehensive pathogenic axis of SRNS that integrates loss of AVIL function with alterations in the action of PLCE1, an established SRNS protein.

The F-actin-binding RapGEF GflB is required for efficient macropinocytosis in Dictyostelium

Macropinocytosis involves the uptake of large volumes of fluid, which is regulated by various small GTPases. The Dictyostelium discoideum protein GflB is a guanine nucleotide exchange factor (GEF) of Rap1, and is involved in chemotaxis. Here, we studied the role of GflB in macropinocytosis, phagocytosis and cytokinesis. In plate culture of vegetative cells, compared with the parental strain AX2, gflB-knockout (KO) cells were flatter and more polarized, whereas GflB-overproducing cells were rounder. The gflB-KO cells exhibited impaired crown formation and retraction, particularly retraction, resulting in more crowns (macropinocytic cups) per cell and longer crown lifetimes. Accordingly, gflB-KO cells showed defects in macropinocytosis and also in phagocytosis and cytokinesis. F-actin levels were elevated in gflB-KO cells. GflB localized to the actin cortex most prominently at crowns and phagocytic cups. The villin headpiece domain (VHP)-like N-terminal domain of GflB directly interacted with F-actin in vitro Furthermore, a domain enriched in basic amino acids interacted with specific membrane cortex structures such as the cleavage furrow. In conclusion, GflB acts as a key local regulator of actin-driven membrane protrusion possibly by modulating Rap1 signaling pathways.

Drosophila Dunc-115 mediates axon projection through actin binding

A central step in organizing the central nervous system development is the growth cone of an axon navigating through guidance cues to reach its specific target. While a great deal of this process has been understood especially in identifying the extracellular guidance cues and their membrane receptors, much less is known about how guidance signals are further relayed to the actin filaments that are central to the mobility of the growth cone. The previous results from our laboratory have shown that Drosophila gene dunc-115 regulates axon projection in the eye and the central nervous system. Furthermore, Dunc-115 has a villin-headpiece (VHD) domain, implying the possibility of binding to actin. To further characterize Dunc-115's functions, we have identified the isoform Dunc-115L as a possible downstream target in relaying guidance cues further down to the cytoskeleton. Specifically, we have shown that Dunc-115 regulates neural connections in both the eye and the central nervous system in Drosophila and that Dunc-115 contains an actin-binding domain potentially capable of binding to actin filaments. In this report, we show that Dunc-115 binds to actin via its VHD domain directly, suggesting a possible mechanism for how Dunc-115 relays guidance signals.

Following easy slope paths on a free energy landscape: the case study of the Trp-cage folding mechanism

In this work a new method for the automatic exploration and calculation of multidimensional free energy landscapes is proposed. Inspired by metadynamics, it uses several collective variables that are relevant for the investigated process and a bias potential that discourages the sampling of already visited configurations. The latter potential allows escaping a local free energy minimum following the direction of slow motions. This is different from metadynamics in which there is no specific direction of the biasing force and the computational effort increases significantly with the number of collective variables. The method is tested on the Ace-Ala3-Nme peptide, and then it is applied to investigate the Trp-cage folding mechanism. For this protein, within a few hundreds of nanoseconds, a broad range of conformations is explored, including nearly native ones, initiating the simulation from a completely unfolded conformation. Finally, several folding/unfolding trajectories give a systematic description of the Trp-cage folding pathways, leading to a unified view for the folding mechanisms of this protein. The proposed mechanism is consistent with NMR chemical shift data at increasing temperature and recent experimental observations pointing to a pivotal role of secondary structure elements in directing the folding process toward the native state.

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.

A Simple Mechanism Underlying the Effect of Protecting Osmolytes on Protein Folding

Osmolytes are small organic compounds that confer to the cell an enhanced adaptability to external conditions. Many osmolytes not only protect the cell from osmotic stress but also stabilize the native structure of proteins. While simplified models able to predict changes to protein stability are available, a general physicochemical explanation of the underlying microscopic mechanism is still missing. Here, we address this issue by performing very long all-atom MD simulations, free energy calculations, and experiments on a well-characterized mini-protein, the villin headpiece. Comparisons between the folding free energy landscapes in pure water and osmolyte solutions, together with experimental validation by means of circular dichroism, unfolding experiments, and NMR, led us to formulate a simple hypothesis for the protecting mechanism. Taken together, our results support a novel mechanistic explanation according to which the main driving force behind native state protection is a change in the solvent rotational diffusion.

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.

Competition between intradomain and interdomain interactions: a buried salt bridge is essential for villin headpiece folding and actin binding

Villin-type headpiece domains are ∼70 residue motifs that reside at the C-terminus of a variety of actin-associated proteins. Villin headpiece (HP67) is a commonly used model system for both experimental and computational studies of protein folding. HP67 is made up of two subdomains that form a tightly packed interface. The isolated C-terminal subdomain of HP67 (HP35) is one of the smallest autonomously folding proteins known. The N-terminal subdomain requires the presence of the C-terminal subdomain to fold. In the structure of HP67, a conserved salt bridge connects N- and C-terminal subdomains. This buried salt bridge between residues E39 and K70 is unusual in a small protein domain. We used mutational analysis, monitored by CD and NMR, and functional assays to determine the role of this buried salt bridge. First, the two residues in the salt bridge were replaced with strictly hydrophobic amino acids, E39M/K70M. Second, the two residues in the salt bridge were swapped, E39K/K70E. Any change from the wild-type salt bridge residues results in unfolding of the N-terminal subdomain, even when the mutations were made in a stabilized variant of HP67. The C-terminal subdomain remains folded in all mutants and is stabilized by some of the mutations. Using actin sedimentation assays, we find that a folded N-terminal domain is essential for specific actin binding. Therefore, the buried salt bridge is required for the specific folding of the N-terminal domain which confers actin-binding activity to villin-type headpiece domains, even though the residues required for this specific interaction destabilize the C-terminal subdomain.

Structural analysis of the conserved ubiquitin-binding motifs (UBMs) of the translesion polymerase iota in complex with ubiquitin

Ubiquitin-binding domains (UBDs) provide specificity to the ubiquitin system, which is also involved in translesion synthesis (TLS) in eukaryotic cells. Upon DNA damage, the UBDs (UBM domains) of polymerase iota (Pol ι) interact with ubiquitinated proliferating cell nuclear antigen to regulate the interchange between processive DNA polymerases and TLS. We report a biophysical analysis and solution structures of the two conserved UBM domains located in the C-terminal tail of murine Pol ι in complex with ubiquitin. The 35-amino acid core folds into a helix-turn-helix motif, which belongs to a novel domain fold. Similar to other UBDs, UBMs bind to ubiquitin on the hydrophobic surface delineated by Leu-8, Ile-44, and Val-70, however, slightly shifted toward the C terminus. In addition, UBMs also use electrostatic interactions to stabilize binding. NMR and fluorescence spectroscopy measurements revealed that UBMs bind monoubiquitin, and Lys-63- but not Lys-48-linked chains. Importantly, these biophysical data are supported by functional studies. Indeed, yeast cells expressing ubiquitin mutants specifically defective for UBM binding are viable but sensitive to DNA damaging conditions that require TLS for repair.

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.

Crystallization of small proteins assisted by green fluorescent protein

The generation of crystal lattice contacts by proteinaceous tags fused to target proteins is an attractive approach to aid in the crystallization of otherwise intractable proteins. Here, the use of green fluorescent protein (GFP) fusions for this purpose is demonstrated, using ubiquitin and the ubiquitin-binding motif (UBM) of Y-family polymerase ι as examples. The structure of the GFP-ubiquitin fusion protein revealed that the crystal lattice was formed by GFP moieties. Ubiquitin was accommodated in the lattice through interactions with the peripheral loops of GFP. However, in the GFP-UBM fusion crystal UBM formed extensive interactions with GFP and these interactions, together with UBM dimerization, mediated the crystal packing. Interestingly, the tyrosine residues that are involved in mediating crystal contacts in both GFP-ubiquitin and GFP-UBM crystals are arranged in a belt on the surface of the β-barrel structure of GFP. Therefore, it is likely that GFP can assist in the crystallization of small proteins and of protein domains in general.

Probing the folding transition state structure of the villin headpiece subdomain via side chain and backbone mutagenesis

Backbone-backbone hydrogen bonds are a common feature of native protein structures, yet their thermodynamic and kinetic influence on folding has long been debated. This is reflected by the disparity between current protein folding models, which place hydrogen bond formation at different stages along the folding trajectory. For example, previous studies have suggested that the denatured state of the villin headpiece subdomain contains a residual helical structure that may provide a bias toward the folded state by confining the conformational search associated with its folding. Although helical hydrogen bonds clearly stabilize the folded state, here we show, using an amide-to-ester mutation strategy, that the formation of backbone hydrogen bonds within helices is not rate-limiting in the folding of the subdomain, thereby suggesting that such hydrogen bonds are unlikely to be formed en route from the denatured to the transition state. On the other hand, elimination of hydrogen bonds within the turn region elicits a slower folding rate, consistent with the hypothesis that these residues are involved in the formation of a folding nucleus. While illustrating a potentially conserved aspect of helix-turn-helix folding, our results further underscore the inherent importance of turns in protein supersecondary structure formation.

Nuclear actin and actin-binding proteins in the regulation of transcription and gene expression

Nuclear actin is involved in the transcription of all three RNA polymerases, in chromatin remodeling and in the formation of heterogeneous nuclear ribonucleoprotein complexes, as well as in recruitment of the histone modifier to the active gene. In addition, actin-binding proteins (ABPs) control actin nucleation, bundling, filament capping, fragmentation and monomer availability in the cytoplasm. In recent years, more and more attention has focused on the role of actin and ABPs in the modulation of the subcellular localization of transcriptional regulators. This review focuses on recent developments in the study of transcription and transcriptional regulation by nuclear actin, and the regulation of muscle-specific gene expression, nuclear receptor and transcription complexes by ABPs. Among the ABPs, striated muscle activator of Rho signaling and actin-binding LIM protein regulate actin dynamics and serum response factor-dependent muscle-specific gene expression. Functionally and structurally unrelated cytoplasmic ABPs interact cooperatively with nuclear receptor and regulate its transactivation. Furthermore, ABPs also participate in the formation of transcription complexes.

The 3D structure of villin as an unusual F-Actin crosslinker

Villin is an F-actin nucleating, crosslinking, severing, and capping protein within the gelsolin superfamily. We have used electron tomography of 2D arrays of villin-crosslinked F-actin to generate 3D images revealing villin's crosslinking structure. In these polar arrays, neighboring filaments are spaced 125.9 +/- 7.1 A apart, offset axially by 17 A, with one villin crosslink per actin crossover. More than 6500 subvolumes containing a single villin crosslink and the neighboring actin filaments were aligned and classified to produce 3D subvolume averages. Placement of a complete villin homology model into the average density reveals that full-length villin binds to different sites on F-actin from those used by other actin-binding proteins and villin's close homolog gelsolin.

Advillin folding takes place on a hypersurface of small dimensionality

All-atom explicit-solvent molecular dynamics simulations have been used to investigate the topological structure of the space explored during folding by the c-terminal fragment of the Advillin headpiece, a 36 amino-acid protein. A fractal dimension analysis shows that the hypersurface explored during the folding process has an approximate dimensionality of only three. It is shown that this low dimensionality persists well above the unfolding temperature and is not present in simple coarse-grained models.

Using an amino acid fluorescence resonance energy transfer pair to probe protein unfolding: application to the villin headpiece subdomain and the LysM domain

Previously, we have shown that p-cyanophenylalanine (Phe CN) and tryptophan (Trp) constitute an efficient fluorescence resonance energy transfer (FRET) pair that has several advantages over commonly used dye pairs. Here, we aim to examine the general applicability of this FRET pair in protein folding-unfolding studies by applying it to the urea-induced unfolding transitions of two small proteins, the villin headpiece subdomain (HP35) and the lysin motif (LysM) domain. Depending on whether Phe CN is exposed to solvent, we are able to extract either qualitative information about the folding pathway, as demonstrated by HP35, which has been suggested to unfold in a stepwise manner, or quantitative thermodynamic and structural information, as demonstrated by LysM, which has been shown to be an ideal two-state folder. Our results show that the unfolding transition of HP35 reported by FRET occurs at a denaturant concentration lower than that measured by circular dichroism (CD) and that the loop linking helix 2 and helix 3 remains compact in the denatured state, which are consistent with the notion that HP35 unfolds in discrete steps and that its unfolded state contains residual structures. On the other hand, our FRET results on the LysM domain allow us to develop a model for extracting structural and thermodynamic parameters about its unfolding, and we find that our results are in agreement with those obtained by other methods. Given the fact that Phe CN is a non-natural amino acid and, thus, amenable to incorporation into peptides and proteins via existing peptide synthesis and protein expression methods, we believe that the FRET method demonstrated here is widely applicable to protein conformational studies, especially to the study of relatively small proteins.

Regulation of cell structure and function by actin-binding proteins: villin's perspective

Villin is a tissue-specific actin modifying protein that is associated with actin filaments in the microvilli and terminal web of epithelial cells. It belongs to a large family of actin-binding proteins which includes actin-capping, -nucleating and/or -severing proteins such as gelsolin, severin, fragmin, adseverin/scinderin and actin crosslinking proteins such as dematin and supervillin. Studies done in epithelial cell lines and villin knock-out mice have demonstrated the function of villin in regulating actin dynamics, cell morphology, epithelial-to-mesenchymal transition, cell migration and cell survival. In addition, the ligand-binding properties of villin (F-actin, G-actin, calcium, phospholipids and phospholipase C-gamma1) are mechanistically important for the crosstalk between signaling pathways and actin reorganization in epithelial cells.

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.

Evolutionary Pareto-optimization of stably folding peptides

BACKGROUND

As a rule, peptides are more flexible and unstructured than proteins with their substantial stabilizing hydrophobic cores. Nevertheless, a few stably folding peptides have been discovered. This raises the question whether there may be more such peptides that are unknown as yet. These molecules could be helpful in basic research and medicine.

RESULTS

As a method to explore the space of conformationally stable peptides, we have developed an evolutionary algorithm that allows optimization of sequences with respect to several criteria simultaneously, for instance stability, accessibility of arbitrary parts of the peptide, etc. In a proof-of-concept experiment we have perturbed the sequence of the peptide Villin Headpiece, known to be stable in vitro. Starting from the perturbed sequence we applied our algorithm to optimize peptide stability and accessibility of a loop. Unexpectedly, two clusters of sequences were generated in this way that, according to our criteria, should form structures with higher stability than the wild-type. The structures in one of the clusters possess a fold that markedly differs from the native fold of Villin Headpiece. One of the mutants predicted to be stable was selected for synthesis, its molecular 3D-structure was characterized by nuclear magnetic resonance spectroscopy, and its stability was measured by circular dichroism. Predicted structure and stability were in good agreement with experiment. Eight other sequences and structures, including five with a non-native fold are provided as bona fide predictions.

CONCLUSION

The results suggest that much more conformationally stable peptides may exist than are known so far, and that small fold classes could comprise well-separated sub-folds.

The isolated sixth gelsolin repeat and headpiece domain of villin bundle F-actin in the presence of calcium and are linked by a 40-residue unstructured sequence

Villin is an F-actin regulating, modular protein with a gelsolin-like core and a distinct C-terminal "headpiece" domain. Localized in the microvilli of the absorptive epithelium, villin can bundle F-actin and, at higher calcium concentrations, is capable of a gelsolin-like F-actin severing. The headpiece domain can, in isolation, bind F-actin and is crucial for F-actin bundling by villin. While the three-dimensional structure of the isolated headpiece is known, its conformation in the context of attachment to the villin core remains unexplored. Furthermore, the dynamics of the linkage of the headpiece to the core has not been determined. To address these issues, we employ a 208-residue modular fragment of villin, D6-HP, which consists of the sixth gelsolin-like domain of villin (D6) and the headpiece (HP). We demonstrate that this protein fragment requires calcium for structural stability and, surprisingly, is capable of Ca2+-dependent F-actin bundling, suggesting that D6 contains a cryptic F-actin binding site. NMR resonance assignments and 15N relaxation measurements of D6-HP in 5 mM Ca2+ demonstrate that D6-HP consists of two independent structural domains (D6 and HP) connected by an unfolded 40-residue linker sequence. The headpiece domain in D6-HP retains its structure and interacts with D6 only through the linker sequence without engaging in other interactions. Chemical shift values indicate essentially the same secondary structure elements for D6 in D6-HP as in the highly homologous gelsolin domain 6. Thus, the headpiece domain of villin is structurally and functionally independent of the core domain.

Role of Drosophila gene dunc-115 in nervous system

Axonal guidance signals are transduced through growth cone surface receptors to the interior leading to changes of actin dynamics and actin binding proteins, which are critical in determining the outcome of actin cytoskeleton reorganization. We report here the characterization of the Drosophila actin binding protein abLIM/Unc-115 homolog Dunc-115 and its role in the nervous system. Three Dunc-115 isoforms are identified as Dunc-115L, M and S, respectively. While Dunc-115L is a canonical homolog of Unc-115 with four LIM domains and one villin headpiece domain, Dunc-115M and S are novel isoforms without counterparts in other species. Our molecular modeling shows Dunc-115L is likely to bind to actin. Mutant analysis reveals that Dunc-115 is involved in axonal projection in both the visual and central nervous system.

Implementation of π-π interactions in molecular dynamics simulation

No explicit pi-pi interaction term has been incorporated in the conventional molecular dynamics (MD) simulation programs in spite of its significant role in the folding of biomolecules and the clustering of organic chemicals. In this article, we propose a technique to emphasize the effect of pi-pi interactions using a function of energy and implement it into an MD simulation program. Several trial calculations show that the pi-pi incorporated program gives improved results consistent with experimental data on atom geometry and has no unfavorable interference with the conventional computational framework. This indicates an importance of the explicit consideration of pi-pi interactions in MD simulation.

Reconciling the solution and X-ray structures of the villin headpiece helical subdomain: molecular dynamics simulations and double mutant cycles reveal a stabilizing cation-pi interaction

The 36-residue helical subdomain of the villin headpiece, HP36, is one of the smallest cooperatively folded proteins, folding on the microsecond time scale. The domain is an extraordinarily popular model system for both experimental and computational studies of protein folding. The structure of HP36 has been determined using X-ray crystallography and NMR spectroscopy, with the resulting structures exhibiting differences in helix packing, van der Waals contacts, and hydrogen bonding. It is important to determine the solution structure of HP36 with as much accuracy as possible since this structure is widely used as a reference for simulations and experiments. We complement the existing data by using all-atom molecular dynamics simulations with explicit solvent to evaluate which of the experimental models is the better representation of HP36 in solution. After simulation for 50 ns initiated with the NMR structure, we observed that the protein spontaneously adopts structures with a backbone conformation, core packing, and C-capping motif on the third helix that are more consistent with the crystal structure. We also examined hydrogen bonding and side chain packing interactions between D44 and R55 and between F47 and R55, respectively, which were observed in the crystal structure but not in the NMR-based solution structure. Simulations showed large fluctuations in the distance between D44 and R55, while the distance between F47 and R55 remained stable, suggesting the formation of a cation-pi interaction between those residues. Experimental double mutant cycles confirmed that the F47-R55 pair has a larger energetic coupling than the D44-R55 interaction. Overall, these combined experimental and computational studies show that the X-ray crystal structure is the better reference structure for HP36 in solution at neutral pH. Our analysis also shows how detailed molecular dynamics simulations combined with experimental validation can help bridge the gap between NMR and crystallographic methods.

Gene expression patterns in a novel animal appendage: the sea urchin pluteus arm

The larval arms of echinoid plutei are used for locomotion and feeding. They are composed of internal calcite skeletal rods covered by an ectoderm layer bearing a ciliary band. Skeletogenesis includes an autonomous molecular differentiation program in primary mesenchyme cells (PMCs), initiated when PMCs leave the vegetal plate for the blastocoel, and a patterning of the differentiated skeletal units that requires molecular cues from the overlaying ectoderm. The arms represent a larval feature that arose in the echinoid lineage during the Paleozoic and offers a subject for the study of gene co-option in the evolution of novel larval features. We isolated new molecular markers in two closely related but differently developing species, Heliocidaris tuberculata and Heliocidaris erythrogramma. We report the expression of a larval arm-associated ectoderm gene tetraspanin, as well as two new PMC markers, advillin and carbonic anhydrase. Tetraspanin localizes to the animal half of blastula stage H. tuberculata and then undergoes a restriction into the putative oral ectoderm and future location of the postoral arms, where it continues to be expressed at the leading edge of both the postoral and anterolateral arms. In H. erythrogramma, its expression initiates in the animal half of blastulae and expands over the entire ectoderm from gastrulation onward. Advillin and carbonic anhydrase are upregulated in the PMCs postgastrulation and localized to the leading edge of the growing larval arms of H. tuberculata but do not exhibit coordinated expression in H. erythrogramma larvae. The tight spatiotemporal regulation of these genes in H. tuberculata along with other ontogenetic and phylogenetic evidence suggest that pluteus arms are novel larval organs, distinguishable from the processes of skeletogenesis per se. The dissociation of expression control in H. erythrogramma suggest that coordinate gene expression in H. tuberculata evolved as part of the evolution of pluteus arms, and is not required for larval or adult development.

Solution structure of a small protein containing a fluorinated side chain in the core

We report the first high-resolution structure for a protein containing a fluorinated side chain. Recently we carried out a systematic evaluation of phenylalanine to pentafluorophenylalanine (Phe --> F(5)-Phe) mutants for the 35-residue chicken villin headpiece subdomain (c-VHP), the hydrophobic core of which features a cluster of three Phe side chains (residues 6, 10, and 17). Phe --> F(5)-Phe mutations are interesting because aryl-perfluoroaryl interactions of optimal geometry are intrinsically more favorable than either aryl-aryl or perfluoroaryl-perfluoroaryl interactions, and because perfluoroaryl units are more hydrophobic than are analogous aryl units. Only one mutation, Phe10 --> F(5)-Phe, was found to provide enhanced tertiary structural stability relative to the native core (by approximately 1 kcal/mol, according to guanidinium chloride denaturation studies). The NMR structure of this mutant, described here, reveals very little variation in backbone conformation or side chain packing relative to the wild type. Thus, although Phe --> F(5)-Phe mutations offer the possibility of greater tertiary structural stability from side chain-side chain attraction and/or side chain desolvation, the constraints associated with the native c-VHP fold apparently prevent the modified polypeptide from taking advantage of this possibility. Our findings are important because they complement several studies that have shown that fluorination of saturated side chain carbon atoms can provide enhanced conformational stability.

Two novel members of the ABLIM protein family, ABLIM-2 and -3, associate with STARS and directly bind F-actin

In addition to regulating cell motility, contractility, and cytokinesis, the actin cytoskeleton plays a critical role in the regulation of transcription and gene expression. We have previously identified a novel muscle-specific actin-binding protein, STARS (striated muscle activator of Rho signaling), which directly binds actin and stimulates serum-response factor (SRF)-dependent transcription. To further dissect the STARS/SRF pathway, we performed a yeast two-hybrid screen of a skeletal muscle cDNA library using STARS as bait, and we identified two novel members of the ABLIM protein family, ABLIM-2 and -3, as STARS-interacting proteins. ABLIM-1, which is expressed in retina, brain, and muscle tissue, has been postulated to function as a tumor suppressor. ABLIM-2 and -3 display distinct tissue-specific expression patterns with the highest expression levels in muscle and neuronal tissue. Moreover, these novel ABLIM proteins strongly bind F-actin, are localized to actin stress fibers, and synergistically enhance STARS-dependent activation of SRF. Conversely, knockdown of endogenous ABLIM expression utilizing small interfering RNA significantly blunted SRF-dependent transcription in C2C12 skeletal muscle cells. These findings suggest that the members of the novel ABLIM protein family may serve as a scaffold for signaling modules of the actin cytoskeleton and thereby modulate transcription.

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.

Identification of the PXW sequence as a structural gatekeeper of the headpiece C-terminal subdomain fold

The HeadPiece (HP) domain, present in several F-actin-binding multi-domain proteins, features a well-conserved, solvent-exposed PXWK motif in its C-terminal subdomain. The latter is an autonomously folding subunit comprised of three alpha-helices organised around a hydrophobic core, with the sequence motif preceding the last helix. We report the contributions of each conserved residue in the PXWK motif to human villin HP function and structure, as well as the structural implications of the naturally occurring Pro to Ala mutation in dematin HP. NMR shift perturbation mapping reveals that substitution of each residue by Ala induces only minor, local perturbations in the full villin HP structure. CD spectroscopic thermal analysis, however, shows that the Pro and Trp residues in the PXWK motif afford stabilising interactions. This indicates that, in addition to the residues in the hydrophobic core, the Trp-Pro stacking within the motif contributes to HP stability. This is reinforced by our data on isolated C-terminal HP subdomains where the Pro is also essential for structure formation, since the villin, but not the dematin, C-terminal subdomain is structured. Proper folding can be induced in the dematin C-terminal subdomain by exchanging the Ala for Pro. Conversely, the reverse substitution in the villin C-terminal subdomain leads to loss of structure. Thus, we demonstrate a crucial role for this proline residue in structural stability and folding potential of HP (sub)domains consistent with Pro-Trp stacking as a more general determinant of protein stability.

A Phosphorylation-Induced Conformation Change in Dematin Headpiece

Dematin is an actin binding protein from the junctional complex of the erythrocyte cytoskeleton. The protein has two actin binding sites and bundles actin filaments in vitro. This actin bundling activity is reversibly regulated by phosphorylation in the carboxyl terminal "headpiece" domain (DHP). DHP is a typical villin-type headpiece actin binding motif and contains a flexible N-terminal loop and an alpha-helical C-terminal subdomain that is phosphorylated at Ser74. The NMR structure of a Ser74-to-Glu mutant (DHPs74e) closely mimics the conformation of phosphorylated DHP. The negative charge at Ser74 does not alter the conformation of the C-terminal subdomain, but attracts the N-terminal loop toward the C terminus, changing the orientation of the N-terminal subdomain. NMR relaxation studies also indicate reduced mobility in the N-terminal loop in DHPs74e. Thus, phosphorylation in DHP serves as a switch controlling the conformational state of DHP and the actin bundling activity of dematin.

Efficient high level expression of peptides and proteins as fusion proteins with the N-terminal domain of L9: application to the villin headpiece helical subdomain

The efficient expression of small to midsize polypeptides and small marginally stable proteins can be difficult. A new protein fusion system is developed to allow the expression of peptides and small proteins. The polypeptide of interest is linked via a Factor Xa cleavage sequence to the C-terminus of the N-terminal domain of the ribosomal protein L9 (NTL9). NTL9 is a small (56 residue) basic protein. The C-terminus of the protein is part of an alpha-helix which extends away from the globular structure thus additional domains can be fused without altering the fold of NTL9. NTL9 expresses at high levels, is extremely soluble, and remains fully folded over a wide temperature and pH range. The protein has a high net positive charge, facilitating purification of fusion proteins by ion exchange chromatography. NTL9 fusions can also be easily purified by reverse phase HPLC. As a test case we demonstrate the high level expression of a small, 36 residue, three helix bundle, the villin headpiece subdomain. This protein is widely used as a model system for folding studies and the development of a simple expression system should facilitate experimental studies of the subdomain. The yield of purified fusion protein is 70 mg/L of culture and the yield of purified villin headpiece subdomain is 24 mg/L of culture. We also demonstrate the use of the fusion system to express a smaller marginally folded peptide fragment of the villin headpiece domain.

High-resolution x-ray crystal structures of the villin headpiece subdomain, an ultrafast folding protein

The 35-residue subdomain of the villin headpiece (HP35) is a small ultrafast folding protein that is being intensely studied by experiments, theory, and simulations. We have solved the x-ray structures of HP35 and its fastest folding mutant [K24 norleucine (nL)] to atomic resolution and compared their experimentally measured folding kinetics by using laser temperature jump. The structures, which are in different space groups, are almost identical to each other but differ significantly from previously solved NMR structures. Hence, the differences between the x-ray and NMR structures are probably not caused by lattice contacts or crystal/solution differences, but reflect the higher accuracy of the x-ray structures. The x-ray structures reveal important details of packing of the hydrophobic core and some additional features, such as cross-helical H bonds. Comparison of the x-ray structures indicates that the nL substitution produces only local perturbations. Consequently, the finding that the small stabilization by the mutation is completely reflected in an increased folding rate suggests that this region of the protein is as structured in the transition state as in the folded structure. It is therefore a target for engineering to increase the folding rate of the subdomain from approximately 0.5 micros(-1) for the nL mutant to the estimated theoretical speed limit of approximately 3 micros(-1).