Collin M. Stultz

Constructing Ensembles for Intrinsically Disordered Proteins

The relatively flat energy landscapes associated with intrinsically disordered proteins makes modeling these systems especially problematic. A comprehensive model for these proteins requires one to build an ensemble consisting of a finite collection of structures, and their corresponding relative stabilities, which adequately capture the range of accessible states of the protein. In this regard, methods that use computational techniques to interpret experimental data in terms of such ensembles are an essential part of the modeling process. In this review, we critically assess the advantages and limitations of current techniques and discuss new methods for the validation of these ensembles.

Strong underwater adhesives made by self-assembling multi-protein nanofibres

Many natural underwater adhesives harness hierarchically assembled amyloid nanostructures to achieve strong and robust interfacial adhesion under dynamic and turbulent environments. Despite recent advances, our understanding of the molecular design, self-assembly and structure-function relationships of these natural amyloid fibres remains limited. Thus, designing biomimetic amyloid-based adhesives remains challenging. Here, we report strong and multi-functional underwater adhesives obtained from fusing mussel foot proteins (Mfps) of Mytilus galloprovincialis with CsgA proteins, the major subunit of Escherichia coli amyloid curli fibres. These hybrid molecular materials hierarchically self-assemble into higher-order structures, in which, according to molecular dynamics simulations, disordered adhesive Mfp domains are exposed on the exterior of amyloid cores formed by CsgA. Our fibres have an underwater adhesion energy approaching 20.9 mJ m(-2), which is 1.5 times greater than the maximum of bio-inspired and bio-derived protein-based underwater adhesives reported thus far. Moreover, they outperform Mfps or curli fibres taken on their own and exhibit better tolerance to auto-oxidation than Mfps at pH ≥ 7.0.Many natural underwater adhesives harness hierarchically assembled amyloid nanostructures to achieve strong and robust interfacial adhesion under dynamic and turbulent environments. Despite recent advances, our understanding of the molecular design, self-assembly, and structure-function relationship of those natural amyloid fibers remains limited. Thus, designing biomimetic amyloid-based adhesives remains challenging. Here, we report strong and multi-functional underwater adhesives obtained from fusing mussel foot proteins (Mfps) of Mytilus galloprovincialis with CsgA proteins, the major subunit of Escherichia coli amyloid curli fibers. These hybrid molecular materials hierarchically self-assemble into higher-order structures, in which, according to molecular dynamics simulations, disordered adhesive Mfp domains are exposed on the exterior of amyloid cores formed by CsgA. Our fibers have an underwater adhesion energy approaching 20.9 mJ/m2, which is 1.5 times greater than the maximum of bio-inspired and bio-derived protein-based underwater adhesives reported thus far. Moreover, they outperform Mfps or curli fibers taken on their own at all pHs and exhibit better tolerance to auto-oxidation than Mfps at pH ≥7.0. This work establishes a platform for engineering multi-component self-assembling materials inspired by nature.

Localized unfolding of collagen explains collagenase cleavage near imino-poor sites.

Collagenases cleave all three chains of type III collagen at specific sites characterized by a Gly-Leu or a Gly-Ile bond that is upstream from an imino acid-poor region. Molecular dynamics trajectories were used to calculate the free energy of unfolding for collagen-like model peptides. The free energy profiles suggest that such imino-poor regions can adopt a low-energy, partially unfolded state where one of the peptide chains forms a solvent-exposed loop. The results are consistent with a model for collagenase cleavage where partial unfolding of imino-poor regions enables collagenases to gain access to their cleavage sites.

Explaining the Structural Plasticity of α-Synuclein

Given that α-synuclein has been implicated in the pathogenesis of several neurodegenerative disorders, deciphering the structure of this protein is of particular importance. While monomeric α-synuclein is disordered in solution, it can form aggregates rich in cross-β structure, relatively long helical segments when bound to micelles or lipid vesicles, and a relatively ordered helical tetramer within the native cell environment. To understand the physical basis underlying this structural plasticity, we generated an ensemble for monomeric α-synuclein using a Bayesian formalism that combines data from NMR chemical shifts, RDCs, and SAXS with molecular simulations. An analysis of the resulting ensemble suggests that a non-negligible fraction of the ensemble (0.08, 95% confidence interval 0.03–0.12) places the minimal toxic aggregation-prone segment in α-synuclein, NAC(8–18), in a solvent exposed and extended conformation that can form cross-β structure. Our data also suggest that a sizable fraction of structures in the ensemble (0.14, 95% confidence interval 0.04–0.23) contains long-range contacts between the N- and C-termini. Moreover, a significant fraction of structures that contain these long-range contacts also place the NAC(8–18) segment in a solvent exposed orientation, a finding in contrast to the theory that such long-range contacts help to prevent aggregation. Lastly, our data suggest that α-synuclein samples structures with amphipathic helices that can self-associate via hydrophobic contacts to form tetrameric structures. Overall, these observations represent a comprehensive view of the unfolded ensemble of monomeric α-synuclein and explain how different conformations can arise from the monomeric protein.

The Dynamic Structure of α‑Synuclein Multimers

α-Synuclein, a protein that forms ordered aggregates in the brains of patients with Parkinsons disease, is intrinsically disordered in the monomeric state. Several studies, however, suggest that it can form soluble multimers in vivo that have significant secondary structure content. A number of studies demonstrate that α-synuclein can form β-strand-rich oligomers that are neurotoxic, and recent observations argue for the existence of soluble helical tetrameric species in cellulo that do not form toxic aggregates. To gain further insight into the different types of multimeric states that this protein can adopt, we generated an ensemble for an α-synuclein construct that contains a 10-residue N-terminal extension, which forms multimers when isolated from Escherichia coli. Data from NMR chemical shifts and residual dipolar couplings were used to guide the construction of the ensemble. Our data suggest that the dominant state of this ensemble is a disordered monomer, complemented by a small fraction of helical trimers and tetramers. Interestingly, the ensemble also contains trimeric and tetrameric oligomers that are rich in β-strand content. These data help to reconcile seemingly contradictory observations that indicate the presence of a helical tetramer in cellulo on the one hand, and a disordered monomer on the other. Furthermore, our findings are consistent with the notion that the helical tetrameric state provides a mechanism for storing α-synuclein when the protein concentration is high, thereby preventing non-membrane-bound monomers from aggregating.

Perturbing the folding energy landscape of the bacterial immunity protein Im7 by site-specific N-linked glycosylation

N-linked glycosylation modulates protein folding and stability through a variety of mechanisms. As such there is considerable interest in the development of general rules to predict the structural consequences of site-specific glycosylation and to understand how these effects can be exploited in the design and development of modified proteins with advantageous properties. In this study, expressed protein ligation is used to create site-specifically glycosylated variants of the bacterial immunity protein Im7 modified with the chitobiose disaccharide (GlcNAc-GlcNAc). Glycans were introduced at seven solvent exposed sites within the Im7 sequence and the kinetic and thermodynamic consequences of N-linked glycosylation analyzed. The values for glycan incorporation were found to range from +5.2 to -3.8 kJ·mol-1. In several cases, glycosylation influences folding by modulating the local conformational preferences of the glycosylated sequence. These locally mediated effects are most prominent in the center of α-helices where glycosylation negatively effects folding and in compact turn motifs between segments of ordered secondary structure where glycosylation promotes folding and enhances the overall stability of the native protein. The studies also provide insight into why glycosylation is commonly identified at the transition between different types of secondary structure and when glycosylation may be used to elaborate protein structure to protect disordered sequences from proteolysis or immune system recognition.

Predicting Protein Structure With Probabilistic Models

Publisher Summary This chapter discusses a probabilistic approach to various algorithmic structure-prediction problems. Important information about protein structure that is difficult or impossible to handle using standard statistical prediction algorithms, such as circular-dichroism studies and disulfide-bridge mapping, are modeled in the chapter as part of the mathematical formulation. The resulting prediction algorithms compute probabilities about the protein structure, given the amino-acid sequence and other experimental results. These probabilities measure the support for modeled hypotheses and are based on the particular protein sequence being analyzed. The probabilistic formulation involves two kinds of quantities—observed and unobserved. The chapter discusses the use of observed quantities to draw inferences about unobserved quantities, such as secondary structure, active sites, or amino-acid interactions with each other, with the solvent, or with ligands. It associates a measure of confidence with each inference.

Fibronectin Unfolding Revisited: Modeling Cell Traction-Mediated Unfolding of the Tenth Type-III Repeat

Fibronectin polymerization is essential for the development and repair of the extracellular matrix. Consequently, deciphering the mechanism of fibronectin fibril formation is of immense interest. Fibronectin fibrillogenesis is driven by cell-traction forces that mechanically unfold particular modules within fibronectin. Previously, mechanical unfolding of fibronectin has been modeled by applying tensile forces at the N- and C-termini of fibronectin domains; however, physiological loading is likely focused on the solvent-exposed RGD loop in the 10th type-III repeat of fibronectin (10FNIII), which mediates binding to cell-surface integrin receptors. In this work we used steered molecular dynamics to study the mechanical unfolding of 10FNIII under tensile force applied at this RGD site. We demonstrate that mechanically unfolding 10FNIII by pulling at the RGD site requires less work than unfolding by pulling at the N- and C- termini. Moreover, pulling at the N- and C-termini leads to 10FNIII unfolding along several pathways while pulling on the RGD site leads to a single exclusive unfolding pathway that includes a partially unfolded intermediate with exposed hydrophobic N-terminal β-strands – residues that may facilitate fibronectin self-association. Additional mechanical unfolding triggers an essential arginine residue, which is required for high affinity binding to integrins, to move to a position far from the integrin binding site. This cell traction-induced conformational change may promote cell detachment after important partially unfolded kinetic intermediates are formed. These data suggest a novel mechanism that explains how cell-mediated forces promote fibronectin fibrillogenesis and how cell surface integrins detach from newly forming fibrils. This process enables cells to bind and unfold additional fibronectin modules – a method that propagates matrix assembly.

The Effect of a ΔK280 Mutation on the Unfolded State of a Microtubule-Binding Repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimers disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a ΔK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and ΔK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of β-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.

Cleavage site specificity and conformational selection in type I collagen degradation.

Excessive degradation of type I collagen is associated with a variety of human diseases such as arthritis, tumor metastasis, and atherosclerosis. Methods that further our understanding of collagenolysis may therefore provide insights into the mechanism of several important disorders. Prior experiments suggest that cleavage of collagen in vitro requires intact full-length collagenase, a multidomain protein containing both a catalytic and a hemopexin-like domain. In this work we demonstrate that type I collagen can be degraded at room temperature, a temperature well below the melting temperature of type I collagen, by collagenase deletion mutants that only contain the catalytic domain of the enzyme. Furthermore, these mutant enzymes hydrolyze the same peptide bond that is recognized by the corresponding full-length enzymes. Hence enzyme specificity at room temperature is achieved without the hemopexin-like domain. We demonstrate that these findings can be explained in light of a conformational selection mechanism that dictates that collagenases preferentially recognize and cleave preformed partially unfolded states of collagen.