Assembly of protein tertiary structures from fragments with similar local sequences using simulated annealing and Bayesian scoring functions

Constructing effective energy functions for protein structure prediction through broadening attraction-basin and reverse Monte Carlo sampling

BACKGROUND

The ab initio approaches to protein structure prediction usually employ the Monte Carlo technique to search the structural conformation that has the lowest energy. However, the widely-used energy functions are usually ineffective for conformation search. How to construct an effective energy function remains a challenging task.

RESULTS

Here, we present a framework to construct effective energy functions for protein structure prediction. Unlike existing energy functions only requiring the native structure to be the lowest one, we attempt to maximize the attraction-basin where the native structure lies in the energy landscape. The underlying rationale is that each energy function determines a specific energy landscape together with a native attraction-basin, and the larger the attraction-basin is, the more likely for the Monte Carlo search procedure to find the native structure. Following this rationale, we constructed effective energy functions as follows: i) To explore the native attraction-basin determined by a certain energy function, we performed reverse Monte Carlo sampling starting from the native structure, identifying the structural conformations on the edge of attraction-basin. ii) To broaden the native attraction-basin, we smoothened the edge points of attraction-basin through tuning weights of energy terms, thus acquiring an improved energy function. Our framework alternates the broadening attraction-basin and reverse sampling steps (thus called BARS) until the native attraction-basin is sufficiently large. We present extensive experimental results to show that using the BARS framework, the constructed energy functions could greatly facilitate protein structure prediction in improving the quality of predicted structures and speeding up conformation search.

CONCLUSION

Using the BARS framework, we constructed effective energy functions for protein structure prediction, which could improve the quality of predicted structures and speed up conformation search as well.

Learning protein constitutive motifs from sequence data

Statistical analysis of evolutionary-related protein sequences provides information about their structure, function, and history. We show that Restricted Boltzmann Machines (RBM), designed to learn complex high-dimensional data and their statistical features, can efficiently model protein families from sequence information. We here apply RBM to 20 protein families, and present detailed results for two short protein domains (Kunitz and WW), one long chaperone protein (Hsp70), and synthetic lattice proteins for benchmarking. The features inferred by the RBM are biologically interpretable: they are related to structure (residue-residue tertiary contacts, extended secondary motifs (α-helixes and β-sheets) and intrinsically disordered regions), to function (activity and ligand specificity), or to phylogenetic identity. In addition, we use RBM to design new protein sequences with putative properties by composing and 'turning up' or 'turning down' the different modes at will. Our work therefore shows that RBM are versatile and practical tools that can be used to unveil and exploit the genotype-phenotype relationship for protein families.

Integrative Structure Modeling: Overview and Assessment

Integrative structure modeling computationally combines data from multiple sources of information with the aim of obtaining structural insights that are not revealed by any single approach alone. In the first part of this review, we survey the commonly used sources of structural information and the computational aspects of model building. Throughout the past decade, integrative modeling was applied to various biological systems, with a focus on large protein complexes. Recent progress in the field of cryo-electron microscopy (cryo-EM) has resolved many of these complexes to near-atomic resolution. In the second part of this review, we compare a range of published integrative models with their higher-resolution counterparts with the aim of critically assessing their accuracy. This comparison gives a favorable view of integrative modeling and demonstrates its ability to yield accurate and informative results. We discuss possible roles of integrative modeling in the new era of cryo-EM and highlight future challenges and directions.

Computational design of structured loops for new protein functions

The ability to engineer the precise geometries, fine-tuned energetics and subtle dynamics that are characteristic of functional proteins is a major unsolved challenge in the field of computational protein design. In natural proteins, functional sites exhibiting these properties often feature structured loops. However, unlike the elements of secondary structures that comprise idealized protein folds, structured loops have been difficult to design computationally. Addressing this shortcoming in a general way is a necessary first step towards the routine design of protein function. In this perspective, we will describe the progress that has been made on this problem and discuss how recent advances in the field of loop structure prediction can be harnessed and applied to the inverse problem of computational loop design.

Studies on a landscape of perfluoroaromatic-reactive peptides

We investigated 26 midsized peptides (∼30 amino acids in length) selected using mRNA display to perform a nucleophilic aromatic substitution reaction (SNAr). Analysis suggested a diverse set of reactive sequences with significant differences in primary sequence, secondary structure and even predicted tertiary structural features. Several of the sequences displayed rapid kinetics allowing for near complete labeling in under one hour. Rosetta ab initio structure prediction of these sequences suggested a landscape of structural features, ranging from beta-sheet-based sequences to those possessing more alpha-helical-like character. Circular dichroism spectroscopy confirmed elements of the structure predictions for the majority of peptides. This analysis additionally uncovered that several peptides underwent secondary structure alterations upon reaction. These results suggest a broad sequence and structural landscape of SNAr active peptides along with a potentially important feature of these biopolymers.

Replacement, Refinement, and Reduction in Animal Studies With Biohazardous Agents

Animal models are critical to the advancement of our knowledge of infectious disease pathogenesis, diagnostics, therapeutics, and prevention strategies. The use of animal models requires thoughtful consideration for their well-being, as infections can significantly impact the general health of an animal and impair their welfare. Application of the 3Rs—replacement, refinement, and reduction—to animal models using biohazardous agents can improve the scientific merit and animal welfare. Replacement of animal models can use in vitro techniques such as cell culture systems, mathematical models, and engineered tissues or invertebrate animal hosts such as amoeba, worms, fruit flies, and cockroaches. Refinements can use a variety of techniques to more closely monitor the course of disease. These include the use of biomarkers, body temperature, behavioral observations, and clinical scoring systems. Reduction is possible using advanced technologies such as in vivo telemetry and imaging, allowing longitudinal assessment of animals during the course of disease. While there is no single method to universally replace, refine, or reduce animal models, the alternatives and techniques discussed are broadly applicable and they should be considered when infectious disease animal models are developed.

Rapid determination of quaternary protein structures in complex biological samples

The understanding of complex biological systems is still hampered by limited knowledge of biologically relevant quaternary protein structures. Here, we demonstrate quaternary structure determination in biological samples using a combination of chemical cross-linking, high-resolution mass spectrometry and high-accuracy protein structure modeling. This approach, termed targeted cross-linking mass spectrometry (TX-MS), relies on computational structural models to score sets of targeted cross-linked peptide signals acquired using a combination of mass spectrometry acquisition techniques. We demonstrate the utility of TX-MS by creating a high-resolution quaternary model of a 1.8 MDa protein complex composed of a pathogen surface protein and ten human plasma proteins. The model is based on a dense network of cross-link distance constraints obtained directly in a mixture of human plasma and live bacteria. These results demonstrate that TX-MS can increase the applicability of flexible backbone docking algorithms to large protein complexes by providing rich cross-link distance information from complex biological samples.

Protein structure determination in complex biological samples is still challenging. Here, the authors develop a computational modeling-guided cross-linking mass spectrometry method, obtaining a high-resolution model of a 1.8 MDa protein assembly from cross-links detected in a mixture of human plasma and bacteria.

Statistical force-field for structural modeling using chemical cross-linking/mass spectrometry distance constraints

Motivation

Chemical cross-linking/mass spectrometry (XLMS) is an experimental method to obtain distance constraints between amino acid residues which can be applied to structural modeling of tertiary and quaternary biomolecular structures. These constraints provide, in principle, only upper limits to the distance between amino acid residues along the surface of the biomolecule. In practice, attempts to use of XLMS constraints for tertiary protein structure determination have not been widely successful. This indicates the need of specifically designed strategies for the representation of these constraints within modeling algorithms.

Results

A force-field designed to represent XLMS-derived constraints is proposed. The potential energy functions are obtained by computing, in the database of known protein structures, the probability of satisfaction of a topological cross-linking distance as a function of the Euclidean distance between amino acid residues. First, the strategy suggests that XL constraints should be set to shorter distances than usually assumed. Second, the complete statistical force-field improves the models obtained and can be easily incorporated into current modeling methods and software. The force-field was implemented and is distributed to be used within the Rosetta ab initio relax protocol.

Availability and implementation

Force-field parameters and usage instructions are freely available online (http://m3g.iqm.unicamp.br/topolink/xlff).

Supplementary information

Supplementary data are available at Bioinformatics online.

Advances in Protein Super-Secondary Structure Prediction and Application to Protein Structure Prediction

Due to the advancement in various sequencing technologies, the gap between the number of protein sequences and the number of experimental protein structures is ever increasing. Community-wide initiatives like CASP have resulted in considerable efforts in the development of computational methods to accurately model protein structures from sequences. Sequence-based prediction of super-secondary structure has direct application in protein structure prediction, and there have been significant efforts in the prediction of super-secondary structure in the last decade. In this chapter, we first introduce the protein structure prediction problem and highlight some of the important progress in the field of protein structure prediction. Next, we discuss recent methods for the prediction of super-secondary structures. Finally, we discuss applications of super-secondary structure prediction in structure prediction/analysis of proteins. We also discuss prediction of protein structures that are composed of simple super-secondary structure repeats and protein structures that are composed of complex super-secondary structure repeats. Finally, we also discuss the recent trends in the field.

Modeling of Protein Tertiary and Quaternary Structures Based on Evolutionary Information

Proteins are subject to evolutionary forces that shape their three-dimensional structure to meet specific functional demands. The knowledge of the structure of a protein is therefore instrumental to gain information about the molecular basis of its function. However, experimental structure determination is inherently time consuming and expensive, making it impossible to follow the explosion of sequence data deriving from genome-scale projects. As a consequence, computational structural modeling techniques have received much attention and established themselves as a valuable complement to experimental structural biology efforts. Among these, comparative modeling remains the method of choice to model the three-dimensional structure of a protein when homology to a protein of known structure can be detected.The general strategy consists of using experimentally determined structures of proteins as templates for the generation of three-dimensional models of related family members (targets) of which the structure is unknown. This chapter provides a description of the individual steps needed to obtain a comparative model using SWISS-MODEL, one of the most widely used automated servers for protein structure homology modeling.

Contribution of Cotranslational Folding Defects to Membrane Protein Homeostasis

Membrane proteins are prone to misfolding and degradation within the cell, yet the nature of the conformational defects involved in this process remain poorly understood. The earliest stages of membrane protein folding are mediated by the Sec61 translocon, a molecular machine that facilitates the lateral partitioning of the polypeptide into the membrane. Proper membrane integration is an essential prerequisite for folding of the nascent chain. However, the marginal energetic drivers of this reaction suggest the translocon may operate with modest fidelity. In this work, we employed biophysical modeling in conjunction with quantitative biochemical measurements in order to evaluate the extent to which cotranslational folding defects influence membrane protein homeostasis. Protein engineering was employed to selectively perturb the topological energetics of human rhodopsin, and the expression and cellular trafficking of engineered variants were quantitatively compared. Our results reveal clear relationships between topological energetics and the efficiency of rhodopsin biogenesis, which appears to be limited by the propensity of a polar transmembrane domain to achieve its correct topological orientation. Though the polarity of this segment is functionally constrained, we find that its topology can be stabilized in a manner that enhances biogenesis without compromising the functional properties of rhodopsin. Furthermore, sequence alignments reveal this topological instability has been conserved throughout the course of evolution. These results suggest that topological defects significantly contribute to the inefficiency of membrane protein folding in the cell. Additionally, our findings suggest that the marginal stability of rhodopsin may represent an evolved trait.

Sampling Native-like Structures of RNA-Protein Complexes through Rosetta Folding and Docking

RNA-protein complexes underlie numerous cellular processes including translation, splicing, and posttranscriptional regulation of gene expression. The structures of these complexes are crucial to their functions but often elude high-resolution structure determination. Computational methods are needed that can integrate low-resolution data for RNA-protein complexes while modeling de novo the large conformational changes of RNA components upon complex formation. To address this challenge, we describe RNP-denovo, a Rosetta method to simultaneously fold-and-dock RNA to a protein surface. On a benchmark set of diverse RNA-protein complexes not solvable with prior strategies, RNP-denovo consistently sampled native-like structures with better than nucleotide resolution. We revisited three past blind modeling challenges involving the spliceosome, telomerase, and a methyltransferase-ribosomal RNA complex in which previous methods gave poor results. When coupled with the same sparse FRET, crosslinking, and functional data used previously, RNP-denovo gave models with significantly improved accuracy. These results open a route to modeling global folds of RNA-protein complexes from low-resolution data.

Systematic Comparison of Amber and Rosetta Energy Functions for Protein Structure Evaluation

An accurate energy function is an essential component of biomolecular structural modeling and design. The comparison of differently derived energy functions enables analysis of the strengths and weaknesses of each energy function and provides independent benchmarks for evaluating improvements within a given energy function. We compared the molecular mechanics Amber empirical energy function to two versions of the Rosetta energy function (talaris2014 and REF2015) in decoy discrimination and loop modeling tests. In decoy discrimination tests, both Rosetta and Amber (ff14SBonlySC) energy functions performed well in scoring the native state as the lowest energy conformation in many cases, but several false minima were found in with both talaris2014 and Amber ff14SBonlySC scoring functions. The current default version of the Rosetta energy function, REF2015, which is parametrized on both small molecule and macromolecular benchmark sets to improve decoy discrimination, performs significantly better than talaris2014, highlighting the improvements made to the Rosetta scoring approach. There are no cases in Rosetta REF2015, and 8/140 cases in Amber, where a false minimum is found that is absent in the alternative landscape. In loop modeling tests, Amber ff14SBonlySC and REF2015 perform equivalently, although false minima are detected in several cases for both. The balance between dihedral, electrostatic, solvation and hydrogen bonding scores contribute to the existence of false minima. To take advantage of the semi-orthogonal nature of the Rosetta and Amber energy functions, we developed a technique that combines Amber and Rosetta conformational rankings to predict the most near-native model for a given protein. This algorithm improves upon predictions from either energy function in isolation and should aid in model selection for structure evaluation and loop modeling tasks.

Generation of Pairwise Potentials Using Multidimensional Data Mining

The rapid development of molecular structural databases provides the chemistry community access to an enormous array of experimental data that can be used to build and validate computational models. Using radial distribution functions collected from experimentally available X-ray and NMR structures, a number of so-called statistical potentials have been developed over the years using the structural data mining strategy. These potentials have been developed within the context of the two-particle Kirkwood equation by extending its original use for isotropic monatomic systems to anisotropic biomolecular systems. However, the accuracy and the unclear physical meaning of statistical potentials have long formed the central arguments against such methods. In this work, we present a new approach to generate molecular energy functions using structural data mining. Instead of employing the Kirkwood equation and introducing the "reference state" approximation, we model the multidimensional probability distributions of the molecular system using graphical models and generate the target pairwise Boltzmann probabilities using the Bayesian field theory. Different from the current statistical potentials that mimic the "knowledge-based" PMF based on the 2-particle Kirkwood equation, the graphical-model-based structure-derived potential developed in this study focuses on the generation of lower-dimensional Boltzmann distributions of atoms through reduction of dimensionality. We have named this new scoring function GARF, and in this work we focus on the mathematical derivation of our novel approach followed by validation studies on its ability to predict protein-ligand interactions.

Improved fragment-based protein structure prediction by redesign of search heuristics

Difficulty in sampling large and complex conformational spaces remains a key limitation in fragment-based de novo prediction of protein structure. Our previous work has shown that even for small-to-medium-sized proteins, some current methods inadequately sample alternative structures. We have developed two new conformational sampling techniques, one employing a bilevel optimisation framework and the other employing iterated local search. We combine strategies of forced structural perturbation (where some fragment insertions are accepted regardless of their impact on scores) and greedy local optimisation, allowing greater exploration of the available conformational space. Comparisons against the Rosetta Abinitio method indicate that our protocols more frequently generate native-like predictions for many targets, even following the low-resolution phase, using a given set of fragment libraries. By contrasting results across two different fragment sets, we show that our methods are able to better take advantage of high-quality fragments. These improvements can also translate into more reliable identification of near-native structures in a simple clustering-based model selection procedure. We show that when fragment libraries are sufficiently well-constructed, improved breadth of exploration within runs improves prediction accuracy. Our results also suggest that in benchmarking scenarios, a total exclusion of fragments drawn from homologous templates can make performance differences between methods appear less pronounced.

CS-ROSETTA

Chemical Shift-Rosetta (CS-Rosetta) is an automated method that employs NMR chemical shifts to model protein structures de novo. In this chapter, we introduce the terminology and central concepts of CS-Rosetta. We describe the architecture and functionality of automatic NOESY assignment (AutoNOE) and structure determination protocols (Abrelax and RASREC) within the CS-Rosetta framework. We further demonstrate how CS-Rosetta can discriminate near-native structures against a large conformational search space using restraints obtained from NMR data, and/or sequence and structure homology. We highlight how CS-Rosetta can be combined with alternative automated approaches to (i) model oligomeric systems and (ii) create NMR-based structure determination pipeline. To show its practical applicability, we emphasize on the computational requirements and performance of CS-Rosetta for protein targets of varying molecular weight and complexity. Finally, we discuss the current Python interface, which enables easy execution of protocols for rapid and accurate high-resolution structure determination.

Recent advances in sequence-based protein structure prediction

The most accurate characterizations of the structure of proteins are provided by structural biology experiments. However, because of the high cost and labor-intensive nature of the structural experiments, the gap between the number of protein sequences and solved structures is widening rapidly. Development of computational methods to accurately model protein structures from sequences is becoming increasingly important to the biological community. In this article, we highlight some important progress in the field of protein structure prediction, especially those related to free modeling (FM) methods that generate structure models without using homologous templates. We also provide a short synopsis of some of the recent advances in FM approaches as demonstrated in the recent Computational Assessment of Structure Prediction competition as well as recent trends and outlook for FM approaches in protein structure prediction.

Contact Potential for Structure Prediction of Proteins and Protein Complexes from Potts Model

The energy function is the key component of protein modeling methodology. This work presents a semianalytical approach to the development of contact potentials for protein structure modeling. Residue-residue and atom-atom contact energies were derived by maximizing the probability of observing native sequences in a nonredundant set of protein structures. The optimization task was formulated as an inverse statistical mechanics problem applied to the Potts model. Its solution by pseudolikelihood maximization provides consistent estimates of coupling constants at atomic and residue levels. The best performance was achieved when interacting atoms were grouped according to their physicochemical properties. For individual protein structures, the performance of the contact potentials in distinguishing near-native structures from the decoys is similar to the top-performing scoring functions. The potentials also yielded significant improvement in the protein docking success rates. The potentials recapitulated experimentally determined protein stability changes upon point mutations and protein-protein binding affinities. The approach offers a different perspective on knowledge-based potentials and may serve as the basis for their further development.

Protein structure prediction

Predicting 3D structure of protein from its amino acid sequence is one of the most important unsolved problems in biophysics and computational biology. This paper attempts to give a comprehensive introduction of the most recent effort and progress on protein structure prediction. Following the general flowchart of structure prediction, related concepts and methods are presented and discussed. Moreover, brief introductions are made to several widely-used prediction methods and the community-wide critical assessment of protein structure prediction (CASP) experiments.

ClusPro PeptiDock: efficient global docking of peptide recognition motifs using FFT

Summary

We present an approach for the efficient docking of peptide motifs to their free receptor structures. Using a motif based search, we can retrieve structural fragments from the Protein Data Bank (PDB) that are very similar to the peptide’s final, bound conformation. We use a Fast Fourier Transform (FFT) based docking method to quickly perform global rigid body docking of these fragments to the receptor. According to CAPRI peptide docking criteria, an acceptable conformation can often be found among the top-ranking predictions.

Availability and Implementation

The method is available as part of the protein-protein docking server ClusPro at https://peptidock.cluspro.org/nousername.php.

Supplementary information

Supplementary data are available at Bioinformatics online.

A random forest learning assisted "divide and conquer" approach for peptide conformation search

Computational determination of peptide conformations is challenging as it is a problem of finding minima in a high-dimensional space. The "divide and conquer" approach is promising for reliably reducing the search space size. A random forest learning model is proposed here to expand the scope of applicability of the "divide and conquer" approach. A random forest classification algorithm is used to characterize the distributions of the backbone φ-ψ units ("words"). A random forest supervised learning model is developed to analyze the combinations of the φ-ψ units ("grammar"). It is found that amino acid residues may be grouped as equivalent "words", while the φ-ψ combinations in low-energy peptide conformations follow a distinct "grammar". The finding of equivalent words empowers the "divide and conquer" method with the flexibility of fragment substitution. The learnt grammar is used to improve the efficiency of the "divide and conquer" method by removing unfavorable φ-ψ combinations without the need of dedicated human effort. The machine learning assisted search method is illustrated by efficiently searching the conformations of GGG/AAA/GGGG/AAAA/GGGGG through assembling the structures of GFG/GFGG. Moreover, the computational cost of the new method is shown to increase rather slowly with the peptide length.

Loop modelling 1.0

Engineering surface loops is a sub-topic of protein engineering that is used routinely in many research fields in academia and industry alike. We provide some tools that search in the PDB for loops satisfying a wide variety of constraints. We illustrate the usefulness of these tools by applying them to a series of recently published studies that included loop engineering or loop modelling. LoopFinder finds loops that fit between two anchor stretches of typically 2, 3, or 4 amino acids each. ProDA find loops of a given length with predefined secondary structure, residue types, hydrophobicity, etc. WHAT IF has gotten a series of new options to scan the whole PDB for loops combining the LoopFinder and ProDA techniques. The open nature of these tools will allow bioinformaticians in this field to easily design their own loop modelling software around our tools.

AVAILABILITY AND IMPLEMENTATION

LoopFinder is a stand-alone Fortran program that is likely to compile and run on every computer. The LoopFinder source code, data files, and documentation are freely available from swift.cmbi.ru.nl/gv/loops/. ProDA is free to all users. There is no login requirement. It is available at: http://bioinf.modares.ac.ir/software/linda/. WHAT IF is shareware that is available from https://swift.cmbi.ru.nl/whatif/.

Rosetta Protein Structure Prediction from Hydroxyl Radical Protein Footprinting Mass Spectrometry Data

In recent years mass spectrometry-based covalent labeling techniques such as hydroxyl radical footprinting (HRF) have emerged as valuable structural biology techniques, yielding information on protein tertiary structure. These data, however, are not sufficient to predict protein structure unambiguously, as they provide information only on the relative solvent exposure of certain residues. Despite some recent advances, no software currently exists that can utilize covalent labeling mass spectrometry data to predict protein tertiary structure. We have developed the first such tool, which incorporates mass spectrometry derived protection factors from HRF labeling as a new centroid score term for the Rosetta scoring function to improve the prediction of protein tertiary structures. We tested our method on a set of four soluble benchmark proteins with known crystal structures and either published HRF experimental results or internally acquired data. Using the HRF labeling data, we rescored large decoy sets of structures predicted with Rosetta for each of the four benchmark proteins. As a result, the model quality improved for all benchmark proteins as compared to when scored with Rosetta alone. For two of the four proteins we were even able to identify atomic resolution models with the addition of HRF data.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

MyPMFs: a simple tool for creating statistical potentials to assess protein structural models

Evaluating the model quality of protein structures that evolve in environments with particular physicochemical properties requires scoring functions that are adapted to their specific residue compositions and/or structural characteristics. Thus, computational methods developed for structures from the cytosol cannot work properly on membrane or secreted proteins. Here, we present MyPMFs, an easy-to-use tool that allows users to train statistical potentials of mean force (PMFs) on the protein structures of their choice, with all parameters being adjustable. We demonstrate its use by creating an accurate statistical potential for transmembrane protein domains. We also show its usefulness to study the influence of the physical environment on residue interactions within protein structures. Our open-source software is freely available for download at https://github.com/bibip-impmc/mypmfs.

Computer-aided design of amino acid-based therapeutics: a review

During the last two decades, the pharmaceutical industry has progressed from detecting small molecules to designing biologic-based therapeutics. Amino acid-based drugs are a group of biologic-based therapeutics that can effectively combat the diseases caused by drug resistance or molecular deficiency. Computational techniques play a key role to design and develop the amino acid-based therapeutics such as proteins, peptides and peptidomimetics. In this study, it was attempted to discuss the various elements for computational design of amino acid-based therapeutics. Protein design seeks to identify the properties of amino acid sequences that fold to predetermined structures with desirable structural and functional characteristics. Peptide drugs occupy a middle space between proteins and small molecules and it is hoped that they can target "undruggable" intracellular protein-protein interactions. Peptidomimetics, the compounds that mimic the biologic characteristics of peptides, present refined pharmacokinetic properties compared to the original peptides. Here, the elaborated techniques that are developed to characterize the amino acid sequences consistent with a specific structure and allow protein design are discussed. Moreover, the key principles and recent advances in currently introduced computational techniques for rational peptide design are spotlighted. The most advanced computational techniques developed to design novel peptidomimetics are also summarized.

A Memetic Algorithm for 3-D Protein Structure Prediction Problem

Memetic Algorithms are population-based metaheuristics intrinsically concerned with exploiting all available knowledge about the problem under study. The incorporation of problem domain knowledge is not an optional mechanism, but a fundamental feature of the Memetic Algorithms. In this paper, we present a Memetic Algorithm to tackle the three-dimensional protein structure prediction problem. The method uses a structured population and incorporates a Simulated Annealing algorithm as a local search strategy, as well as ad-hoc crossover and mutation operators to deal with the problem. It takes advantage of structural knowledge stored in the Protein Data Bank, by using an Angle Probability List that helps to reduce the search space and to guide the search strategy. The proposed algorithm was tested on nineteen protein sequences of amino acid residues, and the results show the ability of the algorithm to find native-like protein structures. Experimental results have revealed that the proposed algorithm can find good solutions regarding root-mean-square deviation and global distance total score test in comparison with the experimental protein structures. We also show that our results are comparable in terms of folding organization with state-of-the-art prediction methods, corroborating the effectiveness of our proposal.

A Novel Method for Drug Screen to Regulate G Protein-Coupled Receptors in the Metabolic Network of Alzheimer's Disease

Alzheimer's disease (AD) is a chronic and progressive neurodegenerative disorder and the pathogenesis of AD is poorly understood. G protein-coupled receptors (GPCRs) are involved in numerous key AD pathways and play a key role in the pathology of AD. To fully understand the pathogenesis of AD and design novel drug therapeutics, analyzing the connection between AD and GPCRs is of great importance. In this paper, we firstly build and analyze the AD-related pathway by consulting the KEGG pathway of AD and a mass of literature and collect 25 AD-related GPCRs for drug discovery. Then the ILbind and AutoDock Vina tools are integrated to find out potential drugs related to AD. According to the analysis of DUD-E dataset, we select five drugs, that is, Acarbose (ACR), Carvedilol (CVD), Digoxin (DGX), NADH (NAI), and Telmisartan (TLS), by sorting the ILbind scores (≥0.73). Then depending on their AutoDock Vina scores and pocket position information, the binding patterns of these five drugs are obtained. We analyze the regulation function of GPCRs in the metabolic network of AD based on the drug screen results, which may be helpful for the study of the off-target effect and the side effect of drugs.

Guiding exploration in conformational feature space with Lipschitz underestimation for ab-initio protein structure prediction

Computing conformations which are essential to associate structural and functional information with gene sequences, is challenging due to the high dimensionality and rugged energy surface of the protein conformational space. Consequently, the dimension of the protein conformational space should be reduced to a proper level, and an effective exploring algorithm should be proposed. In this paper, a plug-in method for guiding exploration in conformational feature space with Lipschitz underestimation (LUE) for ab-initio protein structure prediction is proposed. The conformational space is converted into ultrafast shape recognition (USR) feature space firstly. Based on the USR feature space, the conformational space can be further converted into Underestimation space according to Lipschitz estimation theory for guiding exploration. As a consequence of the use of underestimation model, the tight lower bound estimate information can be used for exploration guidance, the invalid sampling areas can be eliminated in advance, and the number of energy function evaluations can be reduced. The proposed method provides a novel technique to solve the exploring problem of protein conformational space. LUE is applied to differential evolution (DE) algorithm, and metropolis Monte Carlo(MMC) algorithm which is available in the Rosetta; When LUE is applied to DE and MMC, it will be screened by the underestimation method prior to energy calculation and selection. Further, LUE is compared with DE and MMC by testing on 15 small-to-medium structurally diverse proteins. Test results show that near-native protein structures with higher accuracy can be obtained more rapidly and efficiently with the use of LUE.

Accurately Predicting Disordered Regions of Proteins Using Rosetta ResidueDisorder Application

Although many proteins necessitate well-folded structures to properly instigate their biological functions, a large fraction of functioning proteins contain regions-known as intrinsically disordered protein regions-where stable structures are not likely to form. Notable functional roles of intrinsically disordered proteins are in transcriptional regulation, translation, and cellular signal transduction. Moreover, intrinsically disordered protein regions are highly abundant in many proteins associated with various human diseases, therefore these segments have become attractive drug targets for potential therapeutics. Over the past decades, numerous computational methods have been developed to accurately predict disordered regions of proteins. Here we introduce a user-friendly and reliable approach for the prediction of disordered protein regions using the structure prediction software Rosetta. Using 245 proteins from a benchmark data set (16 DisProt database proteins) and a test data set (229 proteins with NMR data), we use Rosetta to predict the global protein structures and then show that there is a statistically significant difference between Rosetta scores in disordered and ordered regions, with scores being less favorable in disordered regions. Furthermore, the difference in scores between ordered and disordered protein regions is sufficient to accurately identify disordered protein regions. As a result, our Rosetta ResidueDisorder method (benchmark data set prediction accuracy of 71.77% and independent test data set prediction accuracy of 65.37%) outperformed other established disorder prediction tools and did not exhibit a biased prediction toward either ordered or disordered regions. To facilitate usage, a Rosetta application has been developed for the Rosetta ResidueDisorder method.

Protein homology model refinement by large-scale energy optimization

Proteins fold to their lowest free-energy structures, and hence the most straightforward way to increase the accuracy of a partially incorrect protein structure model is to search for the lowest-energy nearby structure. This direct approach has met with little success for two reasons: first, energy function inaccuracies can lead to false energy minima, resulting in model degradation rather than improvement; and second, even with an accurate energy function, the search problem is formidable because the energy only drops considerably in the immediate vicinity of the global minimum, and there are a very large number of degrees of freedom. Here we describe a large-scale energy optimization-based refinement method that incorporates advances in both search and energy function accuracy that can substantially improve the accuracy of low-resolution homology models. The method refined low-resolution homology models into correct folds for 50 of 84 diverse protein families and generated improved models in recent blind structure prediction experiments. Analyses of the basis for these improvements reveal contributions from both the improvements in conformational sampling techniques and the energy function.

Shikimate Induced Transcriptional Activation of Protocatechuate Biosynthesis Genes by QuiR, a LysR-Type Transcriptional Regulator, in Listeria monocytogenes

Listeria monocytogenes is a common foodborne bacterial pathogen that contaminates plant and animal consumable products. The persistent nature of L. monocytogenes is associated with millions of dollars in food recalls annually. Here, we describe the role of shikimate in directly modulating the expression of genes encoding enzymes for the conversion of quinate and shikimate metabolites to protocatechuate. In L. monocytogenes, these genes are found within two operons, named qui1 and qui2. In addition, a gene named quiR, encoding a LysR-Type Transcriptional Regulator (QuiR), is located immediately upstream of the qui1 operon. Transcriptional lacZ-promoter fusion experiments show that QuiR induces gene expression of both qui1 and qui2 operons in the presence of shikimate. Furthermore, co-crystallization of the QuiR effector binding domain in complex with shikimate provides insights into the mechanism of activation of this regulator. Together these data show that upon shikimate accumulation, QuiR activates the transcription of genes encoding enzymes involved in shikimate and quinate utilization for the production of protocatechuate. Furthermore, the accumulation of protocatechuate leads to the inhibition of Listeria growth. Since protocatechuate is not known to be utilized by Listeria, its role is distinct from those established in other bacteria.

Principal component analysis in protein tertiary structure prediction

We discuss applicability of principal component analysis (PCA) for protein tertiary structure prediction from amino acid sequence. The algorithm presented in this paper belongs to the category of protein refinement models and involves establishing a low-dimensional space where the sampling (and optimization) is carried out via particle swarm optimizer (PSO). The reduced space is found via PCA performed for a set of low-energy protein models previously found using different optimization techniques. A high frequency term is added into this expansion by projecting the best decoy into the PCA basis set and calculating the residual model. This term is aimed at providing high frequency details in the energy optimization. The goal of this research is to analyze how the dimensionality reduction affects the prediction capability of the PSO procedure. For that purpose, different proteins from the Critical Assessment of Techniques for Protein Structure Prediction experiments were modeled. In all the cases, both the energy of the best decoy and the distance to the native structure have decreased. Our analysis also shows how the predicted backbone structure of native conformation and of alternative low energy states varies with respect to the PCA dimensionality. Generally speaking, the reconstruction can be successfully achieved with 10 principal components and the high frequency term. We also provide a computational analysis of protein energy landscape for the inverse problem of reconstructing structure from the reduced number of principal components, showing that the dimensionality reduction alleviates the ill-posed character of this high-dimensional energy optimization problem. The procedure explained in this paper is very fast and allows testing different PCA expansions. Our results show that PSO improves the energy of the best decoy used in the PCA when the adequate number of PCA terms is considered.

Finding the needle in the haystack: towards solving the protein-folding problem computationally

Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as "the protein folding problem," has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.

Data-assisted protein structure modeling by global optimization in CASP12

In CASP12, 2 types of data-assisted protein structure modeling were experimented. Either SAXS experimental data or cross-linking experimental data was provided for a selected number of CASP12 targets that the CASP12 predictor could utilize for better protein structure modeling. We devised 2 separate energy terms for SAXS data and cross-linking data to drive the model structures into more native-like structures that satisfied the given experimental data as much as possible. In CASP11, we successfully performed protein structure modeling using simulated sparse and ambiguously assigned NOE data and/or correct residue-residue contact information, where the only energy term that folded the protein into its native structure was the term which was originated from the given experimental data. However, the 2 types of experimental data provided in CASP12 were far from being sufficient enough to fold the target protein into its native structure because SAXS data provides only the overall shape of the molecule and the cross-linking contact information provides only very low-resolution distance information. For this reason, we combined the SAXS or cross-linking energy term with our regular modeling energy function that includes both the template energy term and the de novo energy terms. By optimizing the newly formulated energy function, we obtained protein models that fit better with provided SAXS data than the X-ray structure of the target. However, the improvement of the model relative to the 1 modeled without the SAXS data, was not significant. Consistent structural improvement was achieved by incorporating cross-linking data into the protein structure modeling.

Protein Tertiary Structure by Crosslinking/Mass Spectrometry

Observing the structures of proteins within the cell and tracking structural changes under different cellular conditions are the ultimate challenges for structural biology. This, however, requires an experimental technique that can generate sufficient data for structure determination and is applicable in the native environment of proteins. Crosslinking/mass spectrometry (CLMS) and protein structure determination have recently advanced to meet these requirements and crosslinking-driven de novo structure determination in native environments is now possible. In this opinion article, we highlight recent successes in the field of CLMS with protein structure modeling and challenges it still holds.

The earliest structural studies on proteins using crosslinking/mass spectrometry aimed to elucidate their tertiary three-dimensional structure.

Tertiary structure modeling using crosslinking fell out of favor for almost two decades because crosslink data were not informative to aid structure modeling.

Two game-changing trends emerged: using short-range crosslinkers that capture relevant modeling information and high-density crosslinking.

High-density crosslinking uses unspecific crosslinkers to dramatically increase crosslink numbers.

In addition, computational structure modeling methods made significant progress in exploiting CLMS data.

The combination of high-density crosslinking and computational structure modeling enables the elucidation of tertiary protein structure in native environments.

This sidesteps the key limitation of today’s structure determination methods, which are unable (except for a few, specialized methods) to probe the structure of proteins in cell lysates or even intact cells.

3D structure determination of a protein in living cells using paramagnetic NMR spectroscopy

Determining the three-dimensional structure of a protein in living cells remains particularly challenging. We demonstrated that the integration of site-specific tagging proteins and GPS-Rosetta calculations provides a fast and effective way of determining the structures of proteins in living cells, and in principle the interactions and dynamics of protein-ligand complexes.

Analytical Approaches to Improve Accuracy in Solving the Protein Topology Problem

To take advantage of recent advances in genomics and proteomics it is critical that the three-dimensional physical structure of biological macromolecules be determined. Cryo-Electron Microscopy (cryo-EM) is a promising and improving method for obtaining this data, however resolution is often not sufficient to directly determine the atomic scale structure. Despite this, information for secondary structure locations is detectable. De novo modeling is a computational approach to modeling these macromolecular structures based on cryo-EM derived data. During de novo modeling a mapping between detected secondary structures and the underlying amino acid sequence must be identified. DP-TOSS (Dynamic Programming for determining the Topology Of Secondary Structures) is one tool that attempts to automate the creation of this mapping. By treating the correspondence between the detected structures and the structures predicted from sequence data as a constraint graph problem DP-TOSS achieved good accuracy in its original iteration. In this paper, we propose modifications to the scoring methodology of DP-TOSS to improve its accuracy. Three scoring schemes were applied to DP-TOSS and tested: (i) a skeleton-based scoring function; (ii) a geometry-based analytical function; and (iii) a multi-well potential energy-based function. A test of 25 proteins shows that a combination of these schemes can improve the performance of DP-TOSS to solve the topology determination problem for macromolecule proteins.

THE-DB: a threading model database for comparative protein structure analysis of the E. coli K12 and human proteomes

New methodology must be developed to improve the ability to characterize the growing number of amino acid sequences, which vastly exceeds the number of experimentally determined protein structures. Homologous proteins can be used as structural templates for modeling proteins that do not have experimentally determined structures. However, in many cases, there are no homologous proteins (typically <30% sequence identity) with determined structures from which a query sequence can be reliably modeled. The aim of protein threading is to use features, such as secondary structure, solvent accessibility and torsional angles, in addition to sequence patterns to identify structural templates from the protein databank to assist for full-length atomic-level structural modeling. However, there are still numerous protein sequences for which correct templates cannot be recognized. This raises the question as to what attributes allow query sequences to be matched to the correct but distantly homologous templates. To aid the investigation into this question and to provide genome-score protein structure for the biological community, a database called THE-DB (threading hard and easy protein database) has been developed in which it becomes possible to analyze over 15 000 query sequences from the Escherichia coli (E. coli) K12 and human proteomes, as well as to find their three-dimensional templates derived from the state-of-the-art threading algorithms which is not feasible with existing protein template databases. The E. coli K12 and human data can be downloaded in bulk from the THE-DB page.

High-resolution global peptide-protein docking using fragments-based PIPER-FlexPepDock

Peptide-protein interactions contribute a significant fraction of the protein-protein interactome. Accurate modeling of these interactions is challenging due to the vast conformational space associated with interactions of highly flexible peptides with large receptor surfaces. To address this challenge we developed a fragment based high-resolution peptide-protein docking protocol. By streamlining the Rosetta fragment picker for accurate peptide fragment ensemble generation, the PIPER docking algorithm for exhaustive fragment-receptor rigid-body docking and Rosetta FlexPepDock for flexible full-atom refinement of PIPER docked models, we successfully addressed the challenge of accurate and efficient global peptide-protein docking at high-resolution with remarkable accuracy, as validated on a small but representative set of peptide-protein complex structures well resolved by X-ray crystallography. Our approach opens up the way to high-resolution modeling of many more peptide-protein interactions and to the detailed study of peptide-protein association in general. PIPER-FlexPepDock is freely available to the academic community as a server at http://piperfpd.furmanlab.cs.huji.ac.il.

Peptide-protein interactions are crucial components of various important biological processes in living cells. High-resolution structural information of such interactions provides insight about the underlying biophysical principles governing the interactions, and a starting point for their targeted manipulations. Accurate docking algorithms can help fill the gap between the vast number of these interactions and the small number of experimentally solved structures. However, the accuracies of the existing protocols have been limited, in particular for ab initio docking when no information about the peptide beyond its sequence is available. Here we introduce PIPER-FlexPepDock, a fragment-based global docking protocol for high-resolution modeling of peptide-protein interactions. Integration of accurate and efficient representation of the peptide using fragment ensembles, their fast and exhaustive rigid-body docking, and their subsequent accurate flexible refinement, enables peptide-protein docking of remarkable accuracy. The validation on a representative benchmark set of crystallographically solved high-resolution peptide-protein complexes demonstrates significantly improved performance over all existing docking protocols. This opens up the way to the modeling of many more peptide-protein interactions, and to a more detailed study of peptide-protein association in general.

Integration of molecular dynamics simulation and hotspot residues grafting for de novo scFv design against Salmonella Typhi TolC protein

With the development of de novo binders for protein targets from non-related scaffolds, many possibilities for therapeutics and diagnostics have been created. In this study, we described the use of de novo design approach to create single-chain fragment variable (scFv) for Salmonella enterica subspecies enterica serovar Typhi TolC protein. Typhoid fever is a global health concern in developing and underdeveloped countries. Rapid typhoid diagnostics will improve disease management and therapy. In this work, molecular dynamics simulation was first performed on a homology model of TolC protein in POPE membrane bilayer to obtain the central structure that was subsequently used as the target for scFv design. Potential hotspot residues capable of anchoring the binders to the target were identified by docking "disembodied" amino acid residues against TolC surface. Next, scFv scaffolds were selected from Protein Data Bank to harbor the computed hotspot residues. The hotspot residues were then incorporated into the scFv scaffold complementarity determining regions. The designs recapitulated binding energy, shape complementarity, and interface surface area of natural protein-antibody interfaces. This approach has yielded 5 designs with high binding affinity against TolC that may be beneficial for the future development of antigen-based detection agents for typhoid diagnostics.