Centenary Award and Sir Frederick Gowland Hopkins Memorial Lecture. Protein folding, structure prediction and design

Unfolding of protein using MoS2/SnS2heterostructure for nanopore-based sequencing

Protein sequencing is crucial for understanding the complex mechanisms driving biological functions. However, proteins are usually folded in their native state and the mechanism of fast protein conformation transitions still remains unclear, which make protein sequencing challenging. Molecular dynamics simulations with accurate force field are now able to observe the entire folding/unfolding process, providing valuable insights into protein folding mechanisms. Given that proteins can be unfolded, nanopore technology shows great potential for protein sequencing. In this study, we proposed to use MoS2/SnS2heterostructures to firstly unfold proteins and then detect them by a nanopore in the heterostructural membrane. All-atom molecular dynamics simulations performed in this work provided rich atomic-level information for a comprehensive understanding of protein unfolding process and mechanism on the MoS2/SnS2heterostructure, it was found that the strong binding of protein to SnS2nanostripe and hydrogen bond breaking were the main reasons for unfolding the protein on the heterostructure. After the protein was fully unfolded, it was restrained on the nanostripe because of the affinity of protein to the SnS2nanostripe. Thus by integrating the proposed unfolding technique with nanopore technology, detection of linear unfolded peptide was realized in this work, allowing for the identification of protein components, which is essential for sequencing proteins in the near future.

Design and Engineering of Miniproteins

The potential of miniproteins in the biological and chemical sciences is constantly increasing. Significant progress in the design methodologies has been achieved over the last 30 years. Early approaches based on propensities of individual amino acid residues to form individual secondary structures were subsequently improved by structural analyses using NMR spectroscopy and crystallography. Consequently, computational algorithms were developed, which are now highly successful in designing structures with accuracy often close to atomic range. Further perspectives include construction of miniproteins incorporating non-native secondary structures derived from sequences with units other than α-amino acids. Noteworthy, miniproteins with extended structures, which are now feasibly accessible, are excellent scaffolds for construction of functional molecules.

100th Anniversary of Macromolecular Science Viewpoint: Data-Driven Protein Design

The design of synthetic proteins with the desired function is a long-standing goal in biomolecular science, with broad applications in biochemical engineering, agriculture, medicine, and public health. Rational de novo design and experimental directed evolution have achieved remarkable successes but are challenged by the requirement to find functional "needles" in the vast "haystack" of protein sequence space. Data-driven models for fitness landscapes provide a predictive map between protein sequence and function and can prospectively identify functional candidates for experimental testing to greatly improve the efficiency of this search. This Viewpoint reviews the applications of machine learning and, in particular, deep learning as part of data-driven protein engineering platforms. We highlight recent successes, review promising computational methodologies, and provide an outlook on future challenges and opportunities. The article is written for a broad audience comprising both polymer and protein scientists and computer and data scientists interested in an up-to-date review of recent innovations and opportunities in this rapidly evolving field.

Enhancing fragment-based protein structure prediction by customising fragment cardinality according to local secondary structure

BACKGROUND

Whenever suitable template structures are not available, usage of fragment-based protein structure prediction becomes the only practical alternative as pure ab initio techniques require massive computational resources even for very small proteins. However, inaccuracy of their energy functions and their stochastic nature imposes generation of a large number of decoys to explore adequately the solution space, limiting their usage to small proteins. Taking advantage of the uneven complexity of the sequence-structure relationship of short fragments, we adjusted the fragment insertion process by customising the number of available fragment templates according to the expected complexity of the predicted local secondary structure. Whereas the number of fragments is kept to its default value for coil regions, important and dramatic reductions are proposed for beta sheet and alpha helical regions, respectively.

RESULTS

The evaluation of our fragment selection approach was conducted using an enhanced version of the popular Rosetta fragment-based protein structure prediction tool. It was modified so that the number of fragment candidates used in Rosetta could be adjusted based on the local secondary structure. Compared to Rosetta's standard predictions, our strategy delivered improved first models, + 24% and + 6% in terms of GDT, when using 2000 and 20,000 decoys, respectively, while reducing significantly the number of fragment candidates. Furthermore, our enhanced version of Rosetta is able to deliver with 2000 decoys a performance equivalent to that produced by standard Rosetta while using 20,000 decoys. We hypothesise that, as the fragment insertion process focuses on the most challenging regions, such as coils, fewer decoys are needed to explore satisfactorily conformation spaces.

CONCLUSIONS

Taking advantage of the high accuracy of sequence-based secondary structure predictions, we showed the value of that information to customise the number of candidates used during the fragment insertion process of fragment-based protein structure prediction. Experimentations conducted using standard Rosetta showed that, when using the recommended number of decoys, i.e. 20,000, our strategy produces better results. Alternatively, similar results can be achieved using only 2000 decoys. Consequently, we recommend the adoption of this strategy to either improve significantly model quality or reduce processing times by a factor 10.

Coarse-Grained Modeling of the Interplay between Secondary Structure Propensities and Protein Fold Assembly

We recently developed a new coarse-grained model of protein structure and dynamics [ Dawid et al. J. Chem. Theory Comput. 2017 , 13 ( 11 ), 5766 - 5779 ]. The model assumed a single bead representation of amino acid residues, where positions of such united residues were defined by centers of mass of four amino acid fragments. Replica exchange Monte Carlo sampling of the model chains provided good pictures of modeled structures and their dynamics. In its generic form the statistical knowledge-based force field of the model has been dedicated for single-domain globular proteins. Sequence-specific interactions are defined by three-letter secondary structure data. In the present work we demonstrate that different assignments and/or predictions of secondary structures are usually sufficient for enforcing cooperative formation of native-like folds of SURPASS chains for the majority of single-domain globular proteins. Simulations of native-like structure assembly for a representative set of globular proteins have shown that the accuracy of secondary structure data is usually not crucial for model performance, although some specific errors can strongly distort the obtained three-dimensional structures.

Hierarchical design of artificial proteins and complexes toward synthetic structural biology

In multiscale structural biology, synthetic approaches are important to demonstrate biophysical principles and mechanisms underlying the structure, function, and action of bio-nanomachines. A central goal of "synthetic structural biology" is the design and construction of artificial proteins and protein complexes as desired. In this paper, I review recent remarkable progress of an array of approaches for hierarchical design of artificial proteins and complexes that signpost the path forward toward synthetic structural biology as an emerging interdisciplinary field. Topics covered include combinatorial and protein-engineering approaches for directed evolution of artificial binding proteins and membrane proteins, binary code strategy for structural and functional de novo proteins, protein nanobuilding block strategy for constructing nano-architectures, protein-metal-organic frameworks for 3D protein complex crystals, and rational and computational approaches for design/creation of artificial proteins and complexes, novel protein folds, ideal/optimized protein structures, novel binding proteins for targeted therapeutics, and self-assembling nanomaterials. Protein designers and engineers look toward a bright future in synthetic structural biology for the next generation of biophysics and biotechnology.

Exploration of Protein Unfolding by Modelling Calorimetry Data from Reheating

Studies of protein unfolding mechanisms are critical for understanding protein functions inside cells, de novo protein design as well as defining the role of protein misfolding in neurodegenerative disorders. Calorimetry has proven indispensable in this regard for recording full energetic profiles of protein unfolding and permitting data fitting based on unfolding pathway models. While both kinetic and thermodynamic protein stability are analysed by varying scan rates and reheating, the latter is rarely used in curve-fitting, leading to a significant loss of information from experiments. To extract this information, we propose fitting both first and second scans simultaneously. Four most common single-peak transition models are considered: (i) fully reversible, (ii) fully irreversible, (iii) partially reversible transitions, and (iv) general three-state models. The method is validated using calorimetry data for chicken egg lysozyme, mutated Protein A, three wild-types of haloalkane dehalogenases, and a mutant stabilized by protein engineering. We show that modelling of reheating increases the precision of determination of unfolding mechanisms, free energies, temperatures, and heat capacity differences. Moreover, this modelling indicates whether alternative refolding pathways might occur upon cooling. The Matlab-based data fitting software tool and its user guide are provided as a supplement.

Functional analyses of a human vascular tumor FOS variant identify a novel degradation mechanism and a link to tumorigenesis

Epithelioid hemangioma is a locally aggressive vascular neoplasm, found in bones and soft tissue, whose cause is currently unknown, but may involve oncogene activation. FOS is one of the earliest viral oncogenes to be characterized, and normal cellular FOS forms part of the activator protein 1 (AP-1) transcription factor complex, which plays a pivotal role in cell growth, differentiation, and survival as well as the DNA damage response. Despite this, a causal link between aberrant FOS function and naturally occurring tumors has not yet been established. Here, we describe a thorough molecular and biochemical analysis of a mutant FOS protein we identified in these vascular tumors. The mutant protein lacks a highly conserved helix consisting of the C-terminal four amino acids of FOS, which we show is indispensable for fast, ubiquitin-independent FOS degradation via the 20S proteasome. Our work reveals that FOS stimulates endothelial sprouting and that perturbation of normal FOS degradation could account for the abnormal vessel growth typical of epithelioid hemangioma. To the best of our knowledge, this is the first functional characterization of mutant FOS proteins found in tumors.

3D Computational Modeling of Proteins Using Sparse Paramagnetic NMR Data

Computational modeling of proteins using evolutionary or de novo approaches offers rapid structural characterization, but often suffers from low success rates in generating high quality models comparable to the accuracy of structures observed in X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. A computational/experimental hybrid approach incorporating sparse experimental restraints in computational modeling algorithms drastically improves reliability and accuracy of 3D models. This chapter discusses the use of structural information obtained from various paramagnetic NMR measurements and demonstrates computational algorithms implementing pseudocontact shifts as restraints to determine the structure of proteins at atomic resolution.

Side-chain dynamics analysis of KE07 series

The significant improvement of KE07 series in catalytic activities shows the great success of computational design approaches combined with directed evolution in protein design. Understanding the protein dynamics in the evolutionary optimization process of computationally designed enzyme will provide profound implication to study enzyme function and guide protein design. Here, side chain squared generalized order parameters and entropy of each protein are calculated using 50ns molecular dynamics simulation data in both apo and bound states. Our results show a correlation between the increase of side chain motion amplitude and catalytic efficiency. By analyzing the relationship between these two values, we find side chain squared generalized order parameter is linearly related to side chain entropy, which indicates the computationally designed KE07 series have similar dynamics property with natural enzymes.

Rational design of TNFα binding proteins based on the de novo designed protein DS119

De novo protein design offers templates for engineering tailor-made protein functions and orthogonal protein interaction networks for synthetic biology research. Various computational methods have been developed to introduce functional sites in known protein structures. De novo designed protein scaffolds provide further opportunities for functional protein design. Here we demonstrate the rational design of novel tumor necrosis factor alpha (TNFα) binding proteins using a home-made grafting program AutoMatch. We grafted three key residues from a virus 2L protein to a de novo designed small protein, DS119, with consideration of backbone flexibility. The designed proteins bind to TNFα with micromolar affinities. We further optimized the interface residues with RosettaDesign and significantly improved the binding capacity of one protein Tbab1-4. These designed proteins inhibit the activity of TNFα in cellular luciferase assays. Our work illustrates the potential application of the de novo designed protein DS119 in protein engineering, biomedical research, and protein sequence-structure-function studies.

Design of Self-Assembling Protein-Polymer Conjugates

Protein-polymer conjugates are of particular interest for nanobiotechnology applications because of the various and complementary roles that each component may play in composite hybrid-materials. This chapter focuses on the design principles and applications of self-assembling protein-polymer conjugate materials. We address the general design methodology, from both synthetic and genetic perspective, conjugation strategies, protein vs. polymer driven self-assembly and finally, emerging applications for conjugate materials. By marrying proteins and polymers into conjugated bio-hybrid materials, materials scientists, chemists, and biologists alike, have at their fingertips a vast toolkit for material design. These inherently hierarchical structures give rise to useful patterning, mechanical and transport properties that may help realize new, more efficient materials for energy generation, catalysis, nanorobots, etc.

Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering

Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.

Improving Protein Fold Recognition by Deep Learning Networks

For accurate recognition of protein folds, a deep learning network method (DN-Fold) was developed to predict if a given query-template protein pair belongs to the same structural fold. The input used stemmed from the protein sequence and structural features extracted from the protein pair. We evaluated the performance of DN-Fold along with 18 different methods on Lindahl’s benchmark dataset and on a large benchmark set extracted from SCOP 1.75 consisting of about one million protein pairs, at three different levels of fold recognition (i.e., protein family, superfamily, and fold) depending on the evolutionary distance between protein sequences. The correct recognition rate of ensembled DN-Fold for Top 1 predictions is 84.5%, 61.5%, and 33.6% and for Top 5 is 91.2%, 76.5%, and 60.7% at family, superfamily, and fold levels, respectively. We also evaluated the performance of single DN-Fold (DN-FoldS), which showed the comparable results at the level of family and superfamily, compared to ensemble DN-Fold. Finally, we extended the binary classification problem of fold recognition to real-value regression task, which also show a promising performance. DN-Fold is freely available through a web server at http://iris.rnet.missouri.edu/dnfold.

Bioinspired genotype-phenotype linkages: mimicking cellular compartmentalization for the engineering of functional proteins

The idea of compartmentalization of genotype and phenotype in cells is key for enabling Darwinian evolution. This contribution describes bioinspired systems that use in vitro compartments-water-in-oil droplets and gel-shell beads-for the directed evolution of functional proteins. Technologies based on these principles promise to provide easier access to protein-based therapeutics, reagents for processes involving enzyme catalysis, parts for synthetic biology and materials with biological components.

Self-Assembling Nano-Architectures Created from a Protein Nano-Building Block Using an Intermolecularly Folded Dimeric de Novo Protein

The design of novel proteins that self-assemble into supramolecular complexes is an important step in the development of synthetic biology and nanotechnology. Recently, we described the three-dimensional structure of WA20, a de novo protein that forms an intermolecularly folded dimeric 4-helix bundle (PDB code 3VJF ). To harness the unusual intertwined structure of WA20 for the self-assembly of supramolecular nanostructures, we created a protein nanobuilding block (PN-Block), called WA20-foldon, by fusing the dimeric structure of WA20 to the trimeric foldon domain of fibritin from bacteriophage T4. The WA20-foldon fusion protein was expressed in the soluble fraction in Escherichia coli, purified, and shown to form several homooligomeric forms. The stable oligomeric forms were further purified and characterized by a range of biophysical techniques. Size exclusion chromatography, multiangle light scattering, analytical ultracentrifugation, and small-angle X-ray scattering (SAXS) analyses indicate that the small (S form), middle (M form), and large (L form) forms of the WA20-foldon oligomers exist as hexamer (6-mer), dodecamer (12-mer), and octadecamer (18-mer), respectively. These findings suggest that the oligomers in multiples of 6-mer are stably formed by fusing the interdigitated dimer of WA20 with the trimer of foldon domain. Pair-distance distribution functions obtained from the Fourier inversion of the SAXS data suggest that the S and M forms have barrel- and tetrahedron-like shapes, respectively. These results demonstrate that the de novo WA20-foldon is an effective building block for the creation of self-assembling artificial nanoarchitectures.

Application of Enhanced Sampling Monte Carlo Methods for High-Resolution Protein-Protein Docking in Rosetta

The high-resolution refinement of docked protein-protein complexes can provide valuable structural and mechanistic insight into protein complex formation complementing experiment. Monte Carlo (MC) based approaches are frequently applied to sample putative interaction geometries of proteins including also possible conformational changes of the binding partners. In order to explore efficiency improvements of the MC sampling, several enhanced sampling techniques, including temperature or Hamiltonian replica exchange and well-tempered ensemble approaches, have been combined with the MC method and were evaluated on 20 protein complexes using unbound partner structures. The well-tempered ensemble method combined with a 2-dimensional temperature and Hamiltonian replica exchange scheme (WTE-H-REMC) was identified as the most efficient search strategy. Comparison with prolonged MC searches indicates that the WTE-H-REMC approach requires approximately 5 times fewer MC steps to identify near native docking geometries compared to conventional MC searches.

Computational design of a self-assembling symmetrical β-propeller protein

The modular structure of many protein families, such as β-propeller proteins, strongly implies that duplication played an important role in their evolution, leading to highly symmetrical intermediate forms. Previous attempts to create perfectly symmetrical propeller proteins have failed, however. We have therefore developed a new and rapid computational approach to design such proteins. As a test case, we have created a sixfold symmetrical β-propeller protein and experimentally validated the structure using X-ray crystallography. Each blade consists of 42 residues. Proteins carrying 2-10 identical blades were also expressed and purified. Two or three tandem blades assemble to recreate the highly stable sixfold symmetrical architecture, consistent with the duplication and fusion theory. The other proteins produce different monodisperse complexes, up to 42 blades (180 kDa) in size, which self-assemble according to simple symmetry rules. Our procedure is suitable for creating nano-building blocks from different protein templates of desired symmetry.