Chimeric spider silk proteins mediated by intein result in artificial hybrid silks

Hybrid silks hold a great potential as specific biomaterials due to its controlled mechanical properties. To produce fibers with tunable properties, here we firstly made chimeric proteins in vitro, called W2C4CT and W2C8CT, with ligation of MaSp repetitive modules (C) with AcSp modules (W) by intein trans splicing technology from smaller precursors without final yield reduction. Intein mediated chimeric proteins form fibers at a low concentration of 0.4 mg/mL in 50 mM K₃PO₄ pH 7.5 just drawn by hand. Hybrid fibers show smoother surface, and also have stronger chemical resistance as compared with fibers from W2CT (W fibers) and mixture of W2CT/C8CT (MHF8 fibers). Fibers from chimeric protein W2C4CT (HFH4) have improved mechanical properties than W fibers; however, with more C modules W2C8CT fibers (HFH8) properties decreased, indicates the length proportion of various modules is very important and should be optimized for fibers with specific properties. Generally, hybrid silks generated via chimeric proteins, which can be simplified by intein trans splicing, has greater potential to produce fibers with tunable properties. Our research shows that intein mediated directional protein ligation is a novel way to make large chimeric spider silk proteins and hybrid silks. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 385–392, 2016.     

A proposed model for dragline spider silk self‐assembly: Insights from the effect of the repetitive domain size on fiber properties

Dragline spider silk has been intensively studied for its superior qualities as a biomaterial. In previous studies, we made use of the baculovirus mediated expression system for the production of a recombinant Araneus diadematus spider silk dragline ADF4 protein and its self‐assembly into intricate fibers in host insect cells. In this study, our aim was to explore the function of the major repetitive domain of the dragline spider silk. Thus, we generated an array of synthetic proteins, each containing a different number of identical repeats up to the largest recombinantly expressed spider silk to date. Study of the self‐assembly properties of these proteins showed that depending on the increasing number of repeats they give rise to different assembly phenotypes, from a fully soluble protein to bona fide fibers with superior qualities. The different assembly forms, the corresponding chemical resistance properties obtained as well as ultrastructural studies, revealed novel insights concerning the structure and intermolecular interactions of the repetitive and nonrepetitive domains. Based on these observations and current knowledge in the field, we hereby present a comprehensive hypothetical model for the mechanism of dragline silk self‐assembly and fiber formation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 458–468, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Mapping long-range contacts in a highly unfolded protein

Insights into the earliest events in protein folding can be obtained by analysis of the conformational propensities of unfolded or partly folded states. The structure of the acid-unfolded state of apomyoglobin has been characterized using paramagnetic spin labeling and NMR. Nitroxide side-chains, introduced by coupling to mutant cysteine residues at positions 18, 77, and 133, were used as probes of chain compaction and long-range tertiary contacts. Significant interactions are observed within and between the N and C termini, while the central region of the polypeptide chain behaves as a random polymer. Even in this highly denatured form, the protein samples transient compact states in which there are native-like contacts between the N and C-terminal regions.     

Reducing the computational complexity of protein folding via fragment folding and assembly

Understanding, and ultimately predicting, how a 1-D protein chain reaches its native 3-D fold has been one of the most challenging problems during the last few decades. Data increasingly indicate that protein folding is a hierarchical process. Hence, the question arises as to whether we can use the hierarchical concept to reduce the practically intractable computational times. For such a scheme to work, the first step is to cut the protein sequence into fragments that form local minima on the polypeptide chain. The conformations of such fragments in solution are likely to be similar to those when the fragments are embedded in the native fold, although alternate conformations may be favored during the mutual stabilization in the combinatorial assembly process. Two elements are needed for such cutting: (1) a library of (clustered) fragments derived from known protein structures and (2) an assignment algorithm that selects optimal combinations to "cover" the protein sequence. The next two steps in hierarchical folding schemes, not addressed here, are the combinatorial assembly of the fragments and finally, optimization of the obtained conformations. Here, we address the first step in a hierarchical protein-folding scheme. The input is a target protein sequence and a library of fragments created by clustering building blocks that were generated by cutting all protein structures. The output is a set of cutout fragments. We briefly outline a graph theoretic algorithm that automatically assigns building blocks to the target sequence, and we describe a sample of the results we have obtained.     

Silkworm, Bombyx mori larvae expressed the spider silk protein through a novel Bac-to-Bac/BmNPV baculovirus

Abstract:  The silkworm has become an ideal multicellular eukaryotic model system for basic research. The major advantages of expressing foreign genes in silkworm larvae are the low cost of feeding, the extremely high levels of expression achievable compared with expression in cell lines and increased safety because the baculovirus is noninfectious to vertebrates. In this study, we used a recently developed Bombyx mori Nucleopolyhedrovirus (BmNPV) bacmid to express the spider flagelliform silk gene in silkworm larvae. The recombinant bacmid baculoviruses (rBacmid/BmNPV/Flag) were introduced into the first-day larvae of the fifth instar by subcutaneous injection. The worms presented symptoms typical of NPV infection from 72 h after injection compared with control. The haemolymph was collected from the infected larvae 120 h post-infection and the recombinant 6× His-tagged Flag protein was purified by the Ni-NTA spin kit under denaturing conditions with 8  m urea. A 37.0-kDa protein was visualized both in rBacmid/BmNPV/Flag-infected haemolymph and eluting fraction. The results showed that the Bac-to-Bac/BmNPV baculovirus expression system is an efficient tool to express the target gene in silkworm larvae, which takes only 7–10 days for generating recombinant baculovirus, compared with the traditional homologous recombination method, which needs at least 40 days for multiple rounds of purification and amplification of viruses.     

Assembly of protein tertiary structures from secondary structures using optimized potentials

We present a simulated annealing-based method for the prediction of the tertiary structures of proteins given knowledge of the secondary structure associated with each amino acid in the sequence. The backbone is represented in a detailed fashion whereas the sidechains and pairwise interactions are modeled in a simplified way, following the LINUS model of Srinivasan and Rose. A perceptron-based technique is used to optimize the interaction potentials for a training set of three proteins. For these proteins, the procedure is able to reproduce the tertiary structures to below 3 A in root mean square deviation (rmsd) from the PDB targets. We present the results of tests on twelve other proteins. For half of these, the lowest energy decoy has a rmsd from the native state below 6 A and, in 9 out of 12 cases, we obtain decoys whose rmsd from the native states are also well below 5 A.     

A structural view on spider silk proteins and their role in fiber assembly

Spider silk is the toughest known biomaterial and even outrivals modern synthetic high‐performance materials. The question of understanding fiber formation is how the spider can prevent premature and fatal aggregation processes inside its own body and how the chemical and mechanical stimuli used to induce the fiber formation process translate into structural changes of the silk material, finally leading to controlled and irreversible aggregation. Here, the focus will be on the structure and function of the highly conserved N‐domains and C‐terminal domains of spider dragline silk which, unlike the very long repetitive sequence elements, adopt a folded conformation in solution and are therefore able to control intermolecular interactions and aggregation between other spider silk molecules. The structures of these domains add valuable details for the construction of a molecular picture of the complicated and highly optimized silk assembly process that might be beneficial for large‐scale in vitro fiber formation attempts with recombinant silk material. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

Long- and short-range interactions in native protein structures are consistent/minimally frustrated in sequence space

We show that long- and short-range interactions in almost all protein native structures are actually consistent with each other for coarse-grained energy scales; specifically we mean the long-range inter-residue contact energies and the short-range secondary structure energies based on peptide dihedral angles, which are potentials of mean force evaluated from residue distributions observed in protein native structures. This consistency is observed at equilibrium in sequence space rather than in conformational space. Statistical ensembles of sequences are generated by exchanging residues for each of 797 protein native structures with the Metropolis method. It is shown that adding the other category of interaction to either the short- or long-range interactions decreases the means and variances of those energies for essentially all protein native structures, indicating that both interactions consistently work by more-or-less restricting sequence spaces available to one of the interactions. In addition to this consistency, independence by these interaction classes is also indicated by the fact that there are almost no correlations between them when equilibrated using both interactions and significant but small, positive correlations at equilibrium using only one of the interactions. Evidence is provided that protein native sequences can be regarded approximately as samples from the statistical ensembles of sequences with these energy scales and that all proteins have the same effective conformational temperature. Designing protein structures and sequences to be consistent and minimally frustrated among the various interactions is a most effective way to increase protein stability and foldability.     

Building native protein conformation from NMR backbone chemical shifts using Monte Carlo fragment assembly

We have been analyzing the extent to which protein secondary structure determines protein tertiary structure in simple protein folds. An earlier paper demonstrated that three-dimensional structure can be obtained successfully using only highly approximate backbone torsion angles for every residue. Here, the initial information is further diluted by introducing a realistic degree of experimental uncertainty into this process. In particular, we tackle the practical problem of determining three-dimensional structure solely from backbone chemical shifts, which can be measured directly by NMR and are known to be correlated with a protein's backbone torsion angles. Extending our previous algorithm to incorporate these experimentally determined data, clusters of structures compatible with the experimentally determined chemical shifts were generated by fragment assembly Monte Carlo. The cluster that corresponds to the native conformation was then identified based on four energy terms: steric clash, solvent-squeezing, hydrogen-bonding, and hydrophobic contact. Currently, the method has been applied successfully to five small proteins with simple topology. Although still under development, this approach offers promise for high-throughput NMR structure determination.     

pH-dependent stability and folding kinetics of a protein with an unusual alpha-beta topology: the C-terminal domain of the ribosomal protein L9

The folding kinetics and thermodynamics of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) have been studied as a function of pH. CTL9 is an alpha-beta protein that contains a single beta-sheet with an unusual mixed parallel, anti-parallel topology. The folding is fully reversible and two-state over the entire pH range. Stopped-flow fluorescence and CD experiments yield the same folding rate, and the chevron plots have the characteristic V-shape expected for two-state folding. The values of DeltaG*(H2O) and the m value calculated from the kinetic experiments are in excellent agreement with the equilibrium measurements. The extrapolated initial amplitudes of both the stopped-flow fluorescence and CD measurements show that there is no detectable burst phase intermediate. The domain contains three histidine residues, two of which are largely buried in the native state. They do not participate in salt-bridges or take part in a hydrogen bonded network. NMR measurements reveal that the buried histidine residues have significantly perturbed pK(a) values in the native state. The equilibrium stability and the folding rate are found to be strongly dependent upon their ionization state. There is a linear relationship between the log of the folding rate and DeltaG* (H2O) . The protein is much more stable and folds noticeably faster at pH values above the native state pK(a) values. DeltaG*(H2O) of unfolding increases from 2.90 kcal mol(-1) at pH 5.0 to 6.40 kcal mol(-1) at pH 8.0 while the folding rate increases from 0.60 to 18.7 s(-1). Tanford linkage analysis revealed that the interactions involving the two histidine residues are largely developed in the transition state. The results are compared to other studies of the pH-dependence of folding.     

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.     

Intra- versus intermolecular interactions in monellin: contribution of surface charges to protein assembly

The relative significance of weak non-covalent interactions in biological context has been much debated. Here, we have addressed the contribution of Coulombic interactions to protein stability and assembly experimentally. The sweet protein monellin, a non-covalently linked heterodimeric protein, was chosen for this study because of its ability to spontaneously reconstitute from separated fragments. The reconstitution of monellin mutants containing large surface charge perturbations was compared to the thermostability of structurally equivalent single-chain monellin containing the same sets of mutations under varying salt concentrations. The affinity between monellin fragments is found to correlate with the thermostability of single chain monellin, indicating the involvement of the same underlying Coulombic interactions. This confirms that there are no principal differences in the interactions involved in folding and binding. Based on comparison with a previous mutational study involving hydrophobic core residues, the relative contribution of Coulombic interactions to stability and affinity is modest. However, the Coulombic perturbations only affect the association rates of reconstitution in contrast to perturbations involving hydrophobic residues, which affect primarily the dissociation rates. These results indicate that Coulombic interactions are likely to be of main importance for the association of protein assembly, relevant for functions of proteins.     

Rapid protein domain assignment from amino acid sequence using predicted secondary structure

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (+/-20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed.     

A folded and functional protein domain in an amyloid-like fibril

The effect of the polypeptide environment on polyalanine-induced fibril formation was investigated with amyloidogenic fragments from PAPBN1, a nuclear protein controlling polyadenylation. Mutation-caused extensions of the natural 10 alanine sequence up to maximally 17 alanines result in fibril formation of PABPN1 and the development of the disease oculopharyngeal muscular dystrophy (OPMD). We explored the influence of fibril formation on the structure and function of a one-domain protein linked to the fibril-forming part of PABPN1. The well-characterized, stably folded, one-domain protein, cold-shock protein CspB from Bacillus subtilis, was fused either to the C terminus of the entire N-terminal domain of PABPN1 or directly to peptides consisting of 10 or 17 alanine residues. The fusion protein between the N-terminal domain of PABPN1 and CspB formed fibrils in which the structure and activity of CspB were retained. In the fibrils formed by fusions in which the polyalanine sequence was directly linked to CspB, CspB was unfolded. These results indicate that the folded conformation and the function of a protein domain can be maintained in amyloid-like fibrils, and that the distance between this domain and the fibril plays an important role.     

Functionalized silk assembled from a recombinant spider silk fusion protein (Z‐4RepCT) produced in the methylotrophic yeast Pichia pastoris

Functional biological materials are a growing research area with potential applicability in medicine and biotechnology. Using genetic engineering, the possibility to introduce additional functions into spider silk‐based materials has been realized. Recently, a recombinant spider silk fusion protein, Z‐4RepCT, was produced intracellularly in Escherichia coli and could after purification self‐assemble into silk‐like fibers with ability to bind antibodies via the IgG‐binding Z domain. In this study, the use of the methylotrophic yeast Pichia pastoris for production of Z‐4RepCT has been investigated. Temperature, pH and production time were influencing the amount of soluble Z‐4RepCT retrieved from the extracellular fraction. Purification of secreted Z‐4RepCT resulted in a mixture of full‐length and degraded silk proteins that failed to self‐assemble into fibers. A position in the C‐terminal domain of 4RepCT was identified as being subjected to proteolytic cleavage by proteases in the Pichia culture supernatant. Moreover, the C‐terminal domain was subjected to glycosylation during production in P. pastoris. These observed alterations of the CT domain are suggested to contribute to the failure in fiber assembly. As alternative approach, Z‐4RepCT retrieved from the intracellular fraction, which was less degraded, was used and shown to retain ability to assemble into silk‐like fibers after enzymatic deglycosylation.     

Very fast folding and association of a trimerization domain from bacteriophage T4 fibritin

The foldon domain constitutes the C-terminal 30 amino acid residues of the trimeric protein fibritin from bacteriophage T4. Its function is to promote folding and trimerization of fibritin. We investigated structure, stability and folding mechanism of the isolated foldon domain. The domain folds into the same trimeric beta-propeller structure as in fibritin and undergoes a two-state unfolding transition from folded trimer to unfolded monomers. The folding kinetics involve several consecutive reactions. Structure formation in the region of the single beta-hairpin of each monomer occurs on the submillisecond timescale. This reaction is followed by two consecutive association steps with rate constants of 1.9(+/-0.5)x10(6)M(-1)s(-1) and 5.4(+/-0.3)x10(6)M(-1)s(-1) at 0.58 M GdmCl, respectively. This is similar to the fastest reported bimolecular association reactions for folding of dimeric proteins. At low concentrations of protein, folding shows apparent third-order kinetics. At high concentrations of protein, the reaction becomes almost independent of protein concentrations with a half-time of about 3 ms, indicating that a first-order folding step from a partially folded trimer to the native protein (k=210 +/- 20 s(-1)) becomes rate-limiting. Our results suggest that all steps on the folding/trimerization pathway of the foldon domain are evolutionarily optimized for rapid and specific initiation of trimer formation during fibritin assembly. The results further show that beta-hairpins allow efficient and rapid protein-protein interactions during folding.     

Polytopic membrane protein folding and assembly in vitro and in vivo (Review)

The insertion and folding of proteins in biological membranes during protein synthesis in vivo is fundamental to membrane biogenesis. At present, however, certain molecular aspects of this process can only be understood by complementary studies in vitro. We bring together in vitro and in vivo results, highlighting how the studies inform each other and increase our knowledge of the folding and assembly of polytopic membrane proteins. A notable recent advance is the high-resolution crystal structure of the protein machinery responsible for membrane protein insertion into the endoplasmic reticulum. This provides an opportunity to combine in vitro and in vivo studies at a more sophisticated level and address mechanistic aspects of polytopic protein insertion and folding. Quality control is another important aspect of membrane biogenesis, and we give an overview of the current understanding of this process, focusing on cystic fibrosis as a well-studied paradigm. Mutations in the associated membrane protein, the cystic fibrosis transmembrane conductance regulator (CFTR), can cause the quality control mechanisms to prevent the mutant protein reaching its normal site of action, the cell surface. In vitro studies of CFTR shed light on the possible origins of other clinically relevant folding mutants and highlight the potential synergy between in vitro and in vivo approaches.     

A career pathway in protein folding: from model peptides to postreductionist protein science

This review discusses the inherent challenge of linking "reductionist" approaches to decipher the information encoded in protein sequences with burgeoning efforts to explore protein folding in native environments-"postreductionist" approaches. Because the invitation to write this article came as a result of my selection to receive the 2010 Dorothy Hodgkin Award of the Protein Society, I use examples from my own work to illustrate the evolution from the reductionist to the postreductionist perspective. I am incredibly honored to receive the Hodgkin Award, but I want to emphasize that it is the combined effort, creativity, and talent of many students, postdoctoral fellows, and collaborators over several years that has led to any accomplishments on which this selection is based. Moreover, I do not claim to have unique insight into the topics discussed here; but this writing opportunity allows me to illustrate some threads in the evolution of protein folding research with my own experiences and to point out to those embarking on careers how the twists and turns in anyone's scientific path are influenced and enriched by the scientific context of our research. The path my own career has taken thus far has been shaped by the timing of discoveries in the field of protein science; together with our contemporaries, we become part of a knowledge evolution. In my own case, this has been an epoch of great discovery in protein folding and I feel very fortunate to have participated in it.     

Phase separation and mechanical properties of an elastomeric biomaterial from spider wrapping silk and elastin block copolymers

Elastin and silk spidroins are fibrous, structural proteins with elastomeric properties of extension and recoil. While elastin is highly extensible and has excellent recovery of elastic energy, silks are particularly strong and tough. This study describes the biophysical characterization of recombinant polypeptides designed by combining spider wrapping silk and elastin‐like sequences as a strategy to rationally increase the strength of elastin‐based materials while maintaining extensibility. We demonstrate a thermo‐responsive phase separation and spontaneous colloid‐like droplet formation from silk‐elastin block copolymers, and from a 34 residue disordered region of Argiope trifasciata wrapping silk alone, and measure a comprehensive suite of tensile mechanical properties from cross‐linked materials. Silk‐elastin materials exhibited significantly increased strength, toughness, and stiffness compared to an elastin‐only material, while retaining high failure strains and low energy loss upon recoil. These data demonstrate the mechanical tunability of protein polymer biomaterials through modular, chimeric recombination, and provide structural insights into mechanical design. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 693–703, 2016.