Computational design of proteins stereochemically optimized in size, stability, and folding speed

Artificial proteins potentially barrier-free in the folding kinetics are approached computationally under the guidance of protein-folding theories. The smallest and fastest folding globular protein triple-helix-bundle (THB) is so modified as to minimize or eliminate its presumed barriers in folding speed. As the barriers may reside in the ordering of either secondary or tertiary structure, the elements of both secondary and tertiary structure in the protein are targeted for prenucleation with suitable stereochemically constrained amino acid residues. The required elements of topology and sequence for the THB are optimized independently; first the topology is optimized with simulated annealing in polypeptides of highly simplified alphabet; next, the sequence in side chains is optimized using the standard inverse design methods. The resultant three best-adapted THBs, variable in topology and distinctive in sequences, are assessed by comparing them with a few benchmark proteins. The results of mainly molecular dynamics (MD) comparisons, undertaken in explicit water at different temperatures, show that the designed sequences are favorably placed against the chosen benchmarks as THB proteins potentially thermostable in the native folds. Folding simulation experiments with MD establish that the designed sequences are rapid in the folding of individual helices, but not in the evolution of tertiary structure; energetic cum topological frustrations remain but could be the artifacts of the starting conformations that were chosen in the THBs in the folding simulations. Overall, a practical high-throughput approach for de novo protein design has been developed that may have fruitful application for any type of tertiary structure.     

Beta-sheet folding of 11-kDa fibrillogenic polypeptide is completely aggregation driven

A de novo polypeptide GH(6)[(GA)(3)GY(GA)(3)GE](8)GAH(6) (YE8) was designed and genetically engineered to form antiparallel beta-strands of GAGAGA repeats. Modulation of pH enables control of solubility, folding, and aggregation of YE8 by control of the overall polypeptide charge, a consequence of the protonation or deprotonation of the glutamic acid and histidine residues. YE8 exhibits all the major properties of a fibrillogenic protein providing an excellent model for detailed study of the fibrillation. At neutral pH, YE8 is soluble in disordered form, yet at pH 3.5 folds into a predominantly beta-sheet conformation that is fibrillogenic. Atomic force microscopy and transmission electron microscopy indicated the formation of fibrillar aggregates on well-defined, hydrophobic surfaces. The beta-sheet folding of YE8 exhibited a lag phase that could be eliminated by seeding or stirring. The strong dependence of lag time on polypeptide concentration established the limiting step in aggregation as initiation of beta-sheet folding.     

Reduced alphabet for protein folding prediction

What are the key building blocks that would have been needed to construct complex protein folds? This is an important issue for understanding protein folding mechanism and guiding de novo protein design. Twenty naturally occurring amino acids and eight secondary structures consist of a 28‐letter alphabet to determine folding kinetics and mechanism. Here we predict folding kinetic rates of proteins from many reduced alphabets. We find that a reduced alphabet of 10 letters achieves good correlation with folding rates, close to the one achieved by full 28‐letter alphabet. Many other reduced alphabets are not significantly correlated to folding rates. The finding suggests that not all amino acids and secondary structures are equally important for protein folding. The foldable sequence of a protein could be designed using at least 10 folding units, which can either promote or inhibit protein folding. Reducing alphabet cardinality without losing key folding kinetic information opens the door to potentially faster machine learning and data mining applications in protein structure prediction, sequence alignment and protein design. Proteins 2015; 83:631–639. © 2015 Wiley Periodicals, Inc.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

In silico chaperonin-like cycle helps folding of proteins for structure prediction

Currently, one of the most serious problems in protein-folding simulations for de novo structure prediction is conformational sampling of medium-to-large proteins. In vivo, folding of these proteins is mediated by molecular chaperones. Inspired by the functions of chaperonins, we designed a simple chaperonin-like simulation protocol within the framework of the standard fragment assembly method: in our protocol, the strength of the hydrophobic interaction is periodically modulated to help the protein escape from misfolded structures. We tested this protocol for 38 proteins and found that, using a certain defined criterion of success, our method could successfully predict the native structures of 14 targets, whereas only those of 10 targets were successfully predicted using the standard protocol. In particular, for non-α-helical proteins, our method yielded significantly better predictions than the standard approach. This chaperonin-inspired protocol that enhanced de novo structure prediction using folding simulations may, in turn, provide new insights into the working principles underlying the chaperonin system.     

Protein folding in the ER

The endoplasmic reticulum (ER) is a major protein folding compartment for secreted, plasma membrane and organelle proteins. Each of these newly-synthesized polypeptides folds in a deterministic process, affected by the unique conditions that exist in the ER. An understanding of protein folding in the ER is a fundamental biomolecular challenge at two levels. The first level addresses how the amino acid sequence programs that polypeptide to efficiently arrive at a particular fold out of a multitude of alternatives, and how different sequences obtain similar folds. At the second level are the issues introduced by folding not in the cytosol, but in the ER, including the risk of aggregation in a molecularly crowded environment, accommodation of post-translational modifications and the compatibility with subsequent intracellular trafficking. This review discusses both the physicochemical and cell biological constraints of folding, which are the challenges that the ER molecular chaperones help overcome.     

In vivo newly translated polypeptides are sequestered in a protected folding environment.

Molecular chaperones play a fundamental role in cellular protein folding. Using intact mammalian cells we examined the contribution of cytosolic chaperones to de novo folding. A large fraction of newly translated polypeptides associate transiently with Hsc70 and the chaperonin TRiC/CCT during their biogenesis. The substrate repertoire observed for Hsc70 and TRiC is not identical: Hsc70 interacts with a wide spectrum of polypeptides larger than 20 kDa, while TRiC associates with a diverse set of proteins between 30 and 60 kDa. Overexpression of a bacterial chaperonin 'trap' that irreversibly captures unfolded polypeptides did not interrupt the productive folding pathway. The trap was unable to bind newly translated polypeptides, indicating that folding in mammalian cells occurs without the release of non-native folding intermediates into the bulk cytosol. We conclude that de novo protein folding occurs in a protected environment created by a highly processive chaperone machinery and is directly coupled to translation.     

Force field influences in beta-hairpin folding simulations

All-atom force fields are now routinely used for more detailed understanding of protein folding mechanisms. However, it has been pointed out that use of all-atom force fields does not guarantee more accurate representations of proteins; in fact, sometimes it even leads to biased structural distributions. Indeed, several issues remain to be solved in force field developments, such as accurate treatment of implicit solvation for efficient conformational sampling and proper treatment of backbone interactions for secondary structure propensities. In this study, we first investigate the quality of several recently improved backbone interaction schemes in AMBER for folding simulations of a beta-hairpin peptide, and further study their influences on the peptide's folding mechanism. Due to the significant number of simulations needed for a thorough analysis of tested force fields, the implicit Poisson-Boltzmann solvent was used in all simulations. The chosen implicit solvent was found to be reasonable for studies of secondary structures based on a set of simulations of both alpha-helical and beta-hairpin peptides with the TIP3P explicit solvent as benchmark. Replica exchange molecular dynamics was also utilized for further efficient conformational sampling. Among the tested AMBER force fields, ff03 and a revised ff99 force field were found to produce structural and thermodynamic data in comparably good agreement with the experiment. However, detailed folding pathways, such as the order of backbone hydrogen bond zipping and the existence of intermediate states, are different between the two force fields, leading to force field-dependent folding mechanisms.     

Conformational subspace in simulation of early-stage protein folding

A probability calculus was used to simulate the early stages of protein folding in ab initio structure prediction. The probabilities of particular phi and psi angles for each of 20 amino acids as they occur in crystal forms of proteins were used to calculate the amount of information necessary for the occurrence of given phi and psi angles to be predicted. It was found that the amount of information needed to predict phi and psi angles with 5 degrees precision is much higher than the amount of information actually carried by individual amino acids in the polypeptide chain. To handle this problem, a limited conformational space for the preliminary search for optimal polypeptide structure is proposed based on a simplified geometrical model of the polypeptide chain and on the probability calculus. These two models, geometric and probabilistic, based on different sources, yield a common conclusion concerning how a limited conformational space can represent an early stage of polypeptide chain-folding simulation. The ribonuclease molecule was used to test the limited conformational space as a tool for modeling early-stage folding.     

A mixed-alpha,beta miniprotein stereochemically reprogrammed to high-binding affinity for acetylcholine

The protein-structure space is limited to L configuration in the asymmetric alpha-amino acid structures; the function space, on other hand, seems limitless because of the chemical diversity in the amino acid side chain structures. The chemical diversity in side chain structure may be multiplied beneficially with the stereochemical diversity in main chain structure; thus, de novo protein design may be explored for customizing molecular structures stereochemically and molecular functions chemically. Illustrating de novo design in the structure space of L and D alphabet, canonical all-beta folds of poly-L structure were reprogrammed as bracelet, boat, and canoe-shaped molecules-the "boat" as a receptor-like pocket and the "canoe" as a metal-ion receptor-simply by mutating specific L-amino acid residues to the corresponding D stereochemical structure. Demonstrating customization of molecular shape stereochemically and function chemically, a 15-residue mixed-alpha, beta miniprotein of canonical poly-L structure is now reprogrammed stereochemically as a cup-shaped receptor for acetylcholine via cation-pi interaction with a triad of aromatic side chains placed in mimicry of the acetylcholine-receptor sites both natural and artificial. Evidence from CD, fluorescence, NMR, DSC, ITC, MD, and molecular-docking studies is presented to show that a rationally designed 15-residue mixed-L, D peptide is a cooperatively ordered molecular fold in the stereochemically specified molecular morphology, submicromolar in affinity of acetylcholine and thus an acetylcholine receptor exceptionally small and simple. .     

Ligand-driven vectorial folding of ribosome-bound human CFTR NBD1

The mechanism by which protein folding is coupled to biosynthesis is a critical, but poorly understood, aspect of protein conformational diseases. Here we use fluorescence resonance energy transfer (FRET) to characterize tertiary structural transitions of nascent polypeptides and show that the first nucleotide-binding domain (NBD1) of human CFTR, whose folding is defective in cystic fibrosis, folds via a cotranslational multistep pathway as it is synthesized on the ribosome. Folding begins abruptly as NBD1 residues 389-500 emerge from the ribosome exit tunnel, initiating compaction of a small, N-terminal α/β-subdomain. Real-time kinetics of synchronized nascent chains revealed that subdomain folding is rapid, occurs coincident with synthesis, and is facilitated by direct ATP binding to the nascent polypeptide. These findings localize the major CF defect late in the NBD1 folding pathway and establish a paradigm wherein a cellular ligand promotes vectorial domain folding by facilitating an energetically favored local peptide conformation.     

Lattice simulations of cotranslational folding of single domain proteins

We use lattice protein models and Monte Carlo simulations to study cotranslational folding of small single domain proteins. We show that the assembly of native structure begins during late extrusion stages, but final formation of native state occurs during de novo folding, when all residues are extruded. There are three main results in our study. First, for the sequences displaying two-state refolding mechanism de novo cotranslational folding pathway differs from that sampled in in vitro refolding. The change in folding pathways is due to partial assembly of native interactions during extrusion that results in different starting conditions for in vitro refolding and for de novo cotranslational folding. For small single domain proteins cotranslational folding is slower than in vitro refolding, but is generally fast enough to be completed before the release from a ribosome. Second, we found that until final stages of biosynthesis cotranslational folding is essentially equilibrium. This observation is explained by low stability of structured states for partially extruded chains. Finally, our data suggest that the proteins, which refold in vitro slowly via intermediates, complete their de novo folding after the release from a ribosome. Comparison of our lattice cotranslational simulations with recent experimental and computational studies is discussed.     

Slow and Bimolecular Folding of a De Novo Designed Monomeric Protein DS119

De novo protein design offers a unique means to test and advance our understanding of how proteins fold. However, most current design methods are native structure eccentric and folding kinetics has rarely been considered in the design process. Here, we show that a de novo designed mini-protein DS119, which folds into a βαβ structure, exhibits unusually slow and concentration-dependent folding kinetics. For example, the folding time for 50 μM of DS119 was estimated to be ∼2 s. Stopped-flow fluorescence resonance energy transfer experiments further suggested that its folding was likely facilitated by a transient dimerization process. Taken together, these results highlight the need for consideration of the entire folding energy landscape in de novo protein design and provide evidence suggesting nonnative interactions can play a key role in protein folding.     

GroE-mediated folding of bacterial luciferases in vivo

In this study we present evidence indicating that GroE chaperonins mediate de novo protein folding of heterodimeric and monomeric luciferases under heat shock or sub-heat shock conditions in vivo. The effects of additional groESL and groEL genes on the bioluminescence of Escherichia coli cells expressing different bacterial luciferase genes at various temperatures were directly studied in cells growing in liquid culture. Data indicate that at 42° C GroESL chaperonins are required for the folding of the subunit polypeptide of the heterodimeric luciferase from the mesophilic bacterium Vibrio harveyi MAV (B392). In contrast, the small number of amino acid substitutions present in the luciferase subunit polypeptide from the thermotolerant V. harveyi CTP5 suppresses this requirement for GroE chaperonins, and greatly reduces interaction between the subunit polypeptide and GroEL chaperonin. In addition, GroESL are required for the de novo folding at 37° C of a MAV luciferase fusion polypeptide that is functional as a monomer. No such requirement for luciferase activity is observed at that temperature with a fusion of the CTP5 and subunit polypeptides, although GroE chaperonins can still mediate folding of the CTP5 fusion luciferase. Bacterial luciferases provide a unique system for direct observation of the effects of GroE chaperonins on protein folding and enzyme assembly in living cells. Furthermore, they offer a sensitive and simple assay system for the identification of polypeptide domains required for GroEL protein binding.     

A de novo redesign of the WW domain

We have used a sequence prediction algorithm and a novel sampling method to design protein sequences for the WW domain, a small beta-sheet motif. The procedure, referred to as SPANS, designs sequences to be compatible with an ensemble of closely related polypeptide backbones, mimicking the inherent flexibility of proteins. Two designed sequences (termed SPANS-WW1 and SPANS-WW2), using only naturally occurring L-amino acids, were selected for study and the corresponding polypeptides were prepared in Escherichia coli. Circular dichroism data suggested that both purified polypeptides adopted secondary structure features related to those of the target without the aid of disulfide bridges or bound cofactors. The structure exhibited by SPANS-WW2 melted cooperatively by raising the temperature of the solution. Further analysis of this polypeptide by proton nuclear magnetic resonance spectroscopy demonstrated that at 5 degrees C, it folds into a structure closely resembling a natural WW domain. This achievement constitutes one of a small number of successful de novo protein designs through fully automated computational methods and highlights the feasibility of including backbone flexibility in the design strategy.     

Oligomerization of a symmetric β‐trefoil protein in response to folding nucleus perturbation

Gene duplication and fusion events in protein evolution are postulated to be responsible for the common protein folds exhibiting internal rotational symmetry. Such evolutionary processes can also potentially yield regions of repetitive primary structure. Repetitive primary structure offers the potential for alternative definitions of critical regions, such as the folding nucleus (FN). In principle, more than one instance of the FN potentially enables an alternative folding pathway in the face of a subsequent deleterious mutation. We describe the targeted mutation of the carboxyl‐terminal region of the (internally located) FN of the de novo designed purely‐symmetric β‐trefoil protein Symfoil‐4P. This mutation involves wholesale replacement of a repeating trefoil‐fold motif with a “blade” motif from a β‐propeller protein, and postulated to trap that region of the Symfoil‐4P FN in a nonproductive folding intermediate. The resulting protein (termed “Bladefoil”) is shown to be cooperatively folding, but as a trimeric oligomer. The results illustrate how symmetric protein architectures have potentially diverse folding alternatives available to them, including oligomerization, when preferred pathways are perturbed.     

An evolutionary strategy for all-atom folding of the 60-amino-acid bacterial ribosomal protein l20

We have investigated an evolutionary algorithm for de novo all-atom folding of the bacterial ribosomal protein L20. We report results of two simulations that converge to near-native conformations of this 60-amino-acid, four-helix protein. We observe a steady increase of "native content" in both simulated ensembles and a large number of near-native conformations in their final populations. We argue that these structures represent a significant fraction of the low-energy metastable conformations, which characterize the folding funnel of this protein. These data validate our all-atom free-energy force field PFF01 for tertiary structure prediction of a previously inaccessible structural family of proteins. We also compare folding simulations of the evolutionary algorithm with the basin-hopping technique for the Trp-cage protein. We find that the evolutionary algorithm generates a dynamic memory in the simulated population, which leads to faster overall convergence.     

Chaperone-assisted folding of newly synthesized proteins in the cytosol

The way in which a newly synthesized polypeptide chain folds into its unique three-dimensional structure remains one of the fundamental questions in molecular biology. Protein folding in the cell is a problematic process and, in many cases, requires the assistance of a network of molecular chaperones to support productive protein foldingin vivo. During protein biosynthesis, ribosome-associated chaperones guide the folding of the nascent polypeptide emerging from the ribosomal tunnel. In this review we summarize the basic principles of the protein-folding process and the involved chaperones, and focus on the role of ribosome-associated chaperones. Our discussion emphasizes the bacterial Trigger Factor, which is the best studied chaperone of this type. Recent advances have determined the atomic structure of the Trigger Factor, providing new, exciting insights into the role of ribosome-associated chaperones in co-translational protein folding.     

Extending the concept of template-assembled synthetic proteins.

The creation of native-like macromolecules in copying nature's way represents a fascinating challenge in protein chemistry today. In the absence of a detailed knowledge of the complex folding pathway the ultimate goal in protein de novo design, the construction of artificial proteins with predetermined three-dimensional structure and tailor-made functions based on a defined, generally valid set of rules, appears to be still out of reach. With progress in synthesis strategies and biostructural characterization methods, topological templates have become a versatile tool for inducing and stabilizing secondary and tertiary structures, such as protein loops, beta-turns, alpha-helices, beta-sheets and a variety of folding motifs. In this article, we extend the concept of template-assembled synthetic proteins for the construction of protein-like topologies with multiply bridged, oligocyclic chain architectures termed locked-in tertiary folds that exhibit unique physicochemical and folding properties because of the highly confined conformational space. Furthermore, we show that some fundamental questions in protein assembly can be approached applying the template concept. Using covalent template trapping of self-associated peptide assemblies in aqueous solution the structural and physical forces guiding protein folding, supramolecular assembly and molecular recognition processes can be studied on a molecular level.