Residue Contacts In Proteins Research Articles

Predicting hotspots for disease-causing single nucleotide variants using sequences-based coevolution, network analysis, and machine learning.

To enable personalized medicine, it is important yet highly challenging to accurately predict disease-causing mutations in target proteins at high throughput. Previous computational methods have been developed using evolutionary information in combination with various biochemical and structural features of protein residues to discriminate neutral vs. deleterious mutations. However, the power of these methods is often limited because they either assume known protein structures or treat residues independently without fully considering their interactions. To address the above limitations, we build upon recent progress in machine learning, network analysis, and protein language models, and develop a sequences-based variant site prediction workflow based on the protein residue contact networks: 1. We employ and integrate various methods of building protein residue networks using state-of-the-art coevolution analysis tools (RaptorX, DeepMetaPSICOV, and SPOT-Contact) powered by deep learning. 2. We use machine learning algorithms (Random Forest, Gradient Boosting, and Extreme Gradient Boosting) to optimally combine 20 network centrality scores to jointly predict key residues as hot spots for disease mutations. 3. Using a dataset of 107 proteins rich in disease mutations, we rigorously evaluate the network scores individually and collectively (via machine learning). This work supports a promising strategy of combining an ensemble of network scores based on different coevolution analysis methods (and optionally predictive scores from other methods) via machine learning to predict hotspot sites of disease mutations, which will inform downstream applications of disease diagnosis and targeted drug design.

Unexpected dynamics in femtomolar complexes of binding proteins with peptides

Ultra-tight binding is usually observed for proteins associating with rigidified molecules. Previously, we demonstrated that femtomolar binders derived from the Armadillo repeat proteins (ArmRPs) can be designed to interact very tightly with fully flexible peptides. Here we show for ArmRPs with four and seven sequence-identical internal repeats that the peptide-ArmRP complexes display conformational dynamics. These dynamics stem from transient breakages of individual protein-residue contacts that are unrelated to overall unbinding. The labile contacts involve electrostatic interactions. We speculate that these dynamics allow attaining very high binding affinities, since they reduce entropic losses. Importantly, only NMR techniques can pick up these local events by directly detecting conformational exchange processes without complications from changes in solvent entropy. Furthermore, we demonstrate that the interaction surface of the repeat protein regularizes upon peptide binding to become more compatible with the peptide geometry. These results provide novel design principles for ultra-tight binders.

Protein Residue Contact Prediction Based on Deep Learning and Massive Statistical Features from Multi-Sequence Alignment

Sequence-based protein tertiary structure prediction is of fundamental importance because the function of a protein ultimately depends on its 3D structure. An accurate residue-residue contact map is one of the essential elements for current ab initio prediction protocols of 3D structure prediction. Recently, with the combination of deep learning and direct coupling techniques, the performance of residue contact prediction has achieved significant progress. However, a considerable number of current Deep-Learning (DL)-based prediction methods are usually time-consuming, mainly because they rely on different categories of data types and third-party programs. In this research, we transformed the complex biological problem into a pure computational problem through statistics and artificial intelligence. We have accordingly proposed a feature extraction method to obtain various categories of statistical information from only the multi-sequence alignment, followed by training a DL model for residue-residue contact prediction based on the massive statistical information. The proposed method is robust in terms of different test sets, showed high reliability on model confidence score, could obtain high computational efficiency and achieve comparable prediction precisions with DL methods that relying on multi-source inputs.

Dynamic Spatiotemporal Determinants Modulate GPCR:G protein Coupling Promiscuity and Biased Signaling

Recent studies have shown that G protein coupled receptors (GPCRs) generally maintain a promiscuous coupling repertoire to different Ga protein families. This allows GPCRs to potentiate multiple cell signaling processes, but the mechanism behind how promiscuous GPCRs couple different G proteins remains unclear. Additionally, GPCRs have also been shown to signal through b-arrestins and ligands have been designed to show “bias” in signaling towards b-arrestins versus G proteins. I will demonstrate the importance of both the nature and lifetime of intermolecular contacts between GPCR and G protein or b-arrestin in mediating the coupling strength as well as selective coupling among the G proteins or between G proteins and b-arrestin. Using Molecular Dynamics simulations, we show the temporal persistence of GPCR:tranducer protein residue contacts at the interface play a critical role in selective, biased or promiscuous coupling of the transducer proteins. Our results show that promiscuous GPCRs gain the ability to promiscuously couple likely by shedding constraining contacts that confer selectivity observed among selective GPCRs. We have also used molecular dynamics simulations to identify the residues in the allosteric communication network from the extracellular region to the transducer coupling regions. This residue network shows conserved chemical properties across all class A GPCRs, demonstrating the importance of allosteric communication in modulating ligand bias.

Molecular dynamics, residue network analysis, and cross-correlation matrix to characterize the deleterious missense mutations in GALE causing galactosemia III.

Epimerase-deficiency galactosemia (EDG) is caused by mutations in the UDP-galactose 4'-epimerase enzyme, encoded by gene GALE. Catalyzing the last reaction in the Leloir pathway, UDP-galactose-4-epimerase catalyzes the interconversion of UDP-galactose and UDP-glucose. This study aimed to use in-depth computational strategies to prioritize the pathogenic missense mutations in GALE protein and investigate the systemic behavior, conformational spaces, atomic motions, and cross-correlation matrix of the GALE protein. We searched four databases (dbSNP, ClinVar, UniProt, and HGMD) and major biological literature databases (PubMed, Science Direct, and Google Scholar), for missense mutations that are associated with EDG patients, our search yielded 190 missense mutations. We applied a systematic computational prediction pipeline, including pathogenicity, stability, biochemical, conservational, protein residue contacts, and structural analysis, to predict the pathogenicity of these mutations. We found three mutations (p.K161N, p.R239W, and p.G302D) with a severe phenotype in patients with EDG that correlated with our computational prediction analysis; thus, they were selected for further structural and simulation analyses to compute the flexibility and stability of the mutant GALE proteins. The three mutants were subjected to molecular dynamics simulation (MDS) with native protein for 200 ns using GROMACS. The MDS demonstrated that these mutations affected the beta-sheets and helical region that are responsible for the catalytic activity; subsequently, affects the stability and flexibility of the mutant proteins along with a decrease and more deviations in compactness when compared to that of a native. Also, three mutations created major variations in the combined atomic motions of the catalytic and C-terminal regions. The network analysis of the residues in the native and three mutant protein structures showed disturbed residue contacts occurred owing to the missense mutations. Our findings help to understand the structural behavior of a protein owing to mutation and are intended to serve as a platform for prioritizing mutations, which could be potential targets for drug discovery and development of targeted therapeutics.

DeepHelicon: Accurate prediction of inter-helical residue contacts in transmembrane proteins by residual neural networks

Accurate prediction of amino acid residue contacts is an important prerequisite for generating high-quality 3D models of transmembrane (TM) proteins. While a large number of compositional, evolutionary, and structural properties of proteins can be used to train contact prediction methods, recent research suggests that coevolution between residues provides the strongest indication of their spatial proximity. We have developed a deep learning approach, DeepHelicon, to predict inter-helical residue contacts in TM proteins by considering only coevolutionary features. DeepHelicon comprises a two-stage supervised learning process by residual neural networks for a gradual refinement of contact maps, followed by variance reduction by an ensemble of models. We present a benchmark study of 12 contact predictors and conclude that DeepHelicon together with the two other state-of-the-art methods DeepMetaPSICOV and Membrain2 outperforms the 10 remaining algorithms on all datasets and at all settings. On a set of 44 TM proteins with an average length of 388 residues DeepHelicon achieves the best performance among all benchmarked methods in predicting the top L/5 and L/2 inter-helical contacts, with the mean precision of 87.42% and 77.84%, respectively. On a set of 57 relatively small TM proteins with an average length of 298 residues DeepHelicon ranks second best after DeepMetaPSICOV. DeepHelicon produces the most accurate predictions for large proteins with more than 10 transmembrane helices. Coevolutionary features alone allow to predict inter-helical residue contacts with an accuracy sufficient for generating acceptable 3D models for up to 30% of proteins using a fully automated modeling method such as CONFOLD2.

SSCpred: Single-Sequence-Based Protein Contact Prediction Using Deep Fully Convolutional Network.

There has been a significant improvement in protein residue contact prediction in recent years. Nevertheless, state-of-the-art methods still show deficiencies in the contact prediction of proteins with low-homology information. These top methods depend largely on statistical features that derived from homologous sequences, but previous studies, along with our analyses, show that they are insufficient for inferencing an accurate contact map for nonhomology protein targets. To compensate, we proposed a brand new single-sequence-based contact predictor (SSCpred) that performs prediction through the deep fully convolutional network (Deep FCN) with only the target sequence itself, i.e., without additional homology information. The proposed pipeline makes good use of the target sequence by utilizing the pair-wise encoding technique and Deep FCN. Experimental results demonstrated that SSCpred can produce accurate predictions based on the efficient pipeline. Compared with several most recent methods, SSCpred achieves completive performance on nonhomology targets. Overall, we explored the possibilities of single-sequence-based contact prediction and designed a novel pipeline without using a complex and redundant feature set. The proposed SSCpred can compensate for current methods' disadvantages and achieves better performance on the nonhomology targets. The web server of SSCpred is freely available at http://csbio.njust.edu.cn/bioinf/sscpred/.

Protein contact prediction using metagenome sequence data and residual neural networks.

Almost all protein residue contact prediction methods rely on the availability of deep multiple sequence alignments (MSAs). However, many proteins from the poorly populated families do not have sufficient number of homologs in the conventional UniProt database. Here we aim to solve this issue by exploring the rich sequence data from the metagenome sequencing projects. Based on the improved MSA constructed from the metagenome sequence data, we developed MapPred, a new deep learning-based contact prediction method. MapPred consists of two component methods, DeepMSA and DeepMeta, both trained with the residual neural networks. DeepMSA was inspired by the recent method DeepCov, which was trained on 441 matrices of covariance features. By considering the symmetry of contact map, we reduced the number of matrices to 231, which makes the training more efficient in DeepMSA. Experiments show that DeepMSA outperforms DeepCov by 10-13% in precision. DeepMeta works by combining predicted contacts and other sequence profile features. Experiments on three benchmark datasets suggest that the contribution from the metagenome sequence data is significant with P-values less than 4.04E-17. MapPred is shown to be complementary and comparable the state-of-the-art methods. The success of MapPred is attributed to three factors: the deeper MSA from the metagenome sequence data, improved feature design in DeepMSA and optimized training by the residual neural networks. http://yanglab.nankai.edu.cn/mappred/. Supplementary data are available at Bioinformatics online.

Identification of residue pairing in interacting \u03b2-strands from a predicted residue contact map

BackgroundDespite the rapid progress of protein residue contact prediction, predicted residue contact maps frequently contain many errors. However, information of residue pairing in β strands could be extracted from a noisy contact map, due to the presence of characteristic contact patterns in β-β interactions. This information may benefit the tertiary structure prediction of mainly β proteins. In this work, we propose a novel ridge-detection-based β-β contact predictor to identify residue pairing in β strands from any predicted residue contact map.ResultsOur algorithm RDb2C adopts ridge detection, a well-developed technique in computer image processing, to capture consecutive residue contacts, and then utilizes a novel multi-stage random forest framework to integrate the ridge information and additional features for prediction. Starting from the predicted contact map of CCMpred, RDb2C remarkably outperforms all state-of-the-art methods on two conventional test sets of β proteins (BetaSheet916 and BetaSheet1452), and achieves F1-scores of ~ 62% and ~ 76% at the residue level and strand level, respectively. Taking the prediction of the more advanced RaptorX-Contact as input, RDb2C achieves impressively higher performance, with F1-scores reaching ~ 76% and ~ 86% at the residue level and strand level, respectively. In a test of structural modeling using the top 1 L predicted contacts as constraints, for 61 mainly β proteins, the average TM-score achieves 0.442 when using the raw RaptorX-Contact prediction, but increases to 0.506 when using the improved prediction by RDb2C.ConclusionOur method can significantly improve the prediction of β-β contacts from any predicted residue contact maps. Prediction results of our algorithm could be directly applied to effectively facilitate the practical structure prediction of mainly β proteins.AvailabilityAll source data and codes are available at http://166.111.152.91/Downloads.html or the GitHub address of https://github.com/wzmao/RDb2C.

MemBrain-contact 2.0: a new two-stage machine learning model for the prediction enhancement of transmembrane protein residue contacts in the full chain.

Inter-residue contacts in proteins have been widely acknowledged to be valuable for protein 3 D structure prediction. Accurate prediction of long-range transmembrane inter-helix residue contacts can significantly improve the quality of simulated membrane protein models. In this paper, we present an updated MemBrain predictor, which aims to predict transmembrane protein residue contacts. Our new model benefits from an efficient learning algorithm that can mine latent structural features, which exist in original feature space. The new MemBrain is a two-stage inter-helix contact predictor. The first stage takes sequence-based features as inputs and outputs coarse contact probabilities for each residue pair, which will be further fed into convolutional neural network together with predictions from three direct-coupling analysis approaches in the second stage. Experimental results on the training dataset show that our method achieves an average accuracy of 81.6% for the top L/5 predictions using a strict sequence-based jackknife cross-validation. Evaluated on the test dataset, MemBrain can achieve 79.4% prediction accuracy. Moreover, for the top L/5 predicted long-range loop contacts, the prediction performance can reach an accuracy of 56.4%. These results demonstrate that the new MemBrain is promising for transmembrane protein's contact map prediction. http://www.csbio.sjtu.edu.cn/bioinf/MemBrain/. Supplementary data are available at Bioinformatics online.

ConKit: a python interface to contact predictions.

SummaryRecent advances in protein residue contact prediction algorithms have led to the emergence of many new methods and a variety of file formats. We present ConKit, an open source, modular and extensible Python interface which allows facile conversion between formats and provides an interface to analyses of sequence alignments and sets of contact predictions.Availability and Implementation ConKit is available via the Python Package Index. The documentation can be found at http://www.conkit.org. ConKit is licensed under the BSD 3-Clause.Supplementary information Supplementary data are available at Bioinformatics online.

KScons: a Bayesian approach for protein residue contact prediction using the knob-socket model of protein tertiary structure.

By simplifying the many-bodied complexity of residue packing into patterns of simple pairwise secondary structure interactions between a single knob residue with a three-residue socket, the knob-socket construct allows a more direct incorporation of structural information into the prediction of residue contacts. By modeling the preferences between the amino acid composition of a socket and knob, we undertake an investigation of the knob-socket construct's ability to improve the prediction of residue contacts. The statistical model considers three priors and two posterior estimations to better understand how the input data affects predictions. This produces six implementations of KScons that are tested on three sets: PSICOV, CASP10 and CASP11. We compare against the current leading contact prediction methods. The results demonstrate the usefulness as well as the limits of knob-socket based structural modeling of protein contacts. The construct is able to extract good predictions from known structural homologs, while its performance degrades when no homologs exist. Among our six implementations, KScons MST-MP (which uses the multiple structure alignment prior and marginal posterior incorporating structural homolog information) performs the best in all three prediction sets. An analysis of recall and precision finds that KScons MST-MP improves accuracy not only by improving identification of true positives, but also by decreasing the number of false positives. Over the CASP10 and CASP11 sets, KScons MST-MP performs better than the leading methods using only evolutionary coupling data, but not quite as well as the supervised learning methods of MetaPSICOV and CoinDCA-NN that incorporate a large set of structural features. qiwei.li@rice.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Accurate prediction of helix interactions and residue contacts in membrane proteins.

Accurate prediction of intra-molecular interactions from amino acid sequence is an important pre-requisite for obtaining high-quality protein models. Over the recent years, remarkable progress in this area has been achieved through the application of novel co-variation algorithms, which eliminate transitive evolutionary connections between residues. In this work we present a new contact prediction method for α-helical transmembrane proteins, MemConP, in which evolutionary couplings are combined with a machine learning approach. MemConP achieves a substantially improved accuracy (precision: 56.0%, recall: 17.5%, MCC: 0.288) compared to the use of either machine learning or co-evolution methods alone. The method also achieves 91.4% precision, 42.1% recall and a MCC of 0.490 in predicting helix-helix interactions based on predicted contacts. The approach was trained and rigorously benchmarked by cross-validation and independent testing on up-to-date non-redundant datasets of 90 and 30 experimental three dimensional structures, respectively. MemConP is a standalone tool that can be downloaded together with the associated training data from http://webclu.bio.wzw.tum.de/MemConP.

Amino acid network for prediction of catalytic residues in enzymes: a comparison survey

Catalytic residues play a significant role in enzyme functions. With the recent accumulation of experimentally determined enzyme 3D structures and network theory on protein structures, the prediction of catalytic residues by amino acid network (AAN, where nodes are residues and links are residue interactions) has gained much interest. Computational methods of identifying catalytic residues are traditionally divided into two groups: sequence-based and structure-based methods. Two new structure- based methods are proposed in current advances: AAN and Elastic Network Model (ENM) of enzyme structures. By concentrating on AAN-based approach, we herein summarized network properties for predictions of catalytic residues. AAN attributes were showed responsible for performance improvement, and therefore the combination of AAN with previous sequence and structural information will be a promising direction for further improvement. Advantages and limitations of AAN-based methods, future perspectives on the application of AAN to the study of protein structure-function relationships are discussed.

Modeling proteins as residue interaction networks.

Proteins are the ultimate carriers of information encoded within genes. The salient feature of proteins is their ability to perform a vast array of function with extraordinary specificity. This specificity of proteins is related directly to their structural conformations, which arises as a result of complex interactions among its constituent amino acid residues. Thus, analyzing proteins through these interactions collectively is an extremely important, interesting and effective way to understand their dynamical and structural properties. This review summarizes recent graph theoretical approaches in the field of amino acid residue interaction networks (RIN), [also known as residue interaction graphs (RIG), protein contact networks (PCN) or residue contact networks (RCN)] and the outcomes thereof. It explains the strengths of network formalism in unraveling various aspects of the complex phenomena, upon which the structure and function of proteins have a direct bearing.

Improved contact predictions using the recognition of protein like contact patterns.

Given sufficient large protein families, and using a global statistical inference approach, it is possible to obtain sufficient accuracy in protein residue contact predictions to predict the structure of many proteins. However, these approaches do not consider the fact that the contacts in a protein are neither randomly, nor independently distributed, but actually follow precise rules governed by the structure of the protein and thus are interdependent. Here, we present PconsC2, a novel method that uses a deep learning approach to identify protein-like contact patterns to improve contact predictions. A substantial enhancement can be seen for all contacts independently on the number of aligned sequences, residue separation or secondary structure type, but is largest for β-sheet containing proteins. In addition to being superior to earlier methods based on statistical inferences, in comparison to state of the art methods using machine learning, PconsC2 is superior for families with more than 100 effective sequence homologs. The improved contact prediction enables improved structure prediction.

The D-ring, Not the A-ring, Rotates in Synechococcus OS-B′ Phytochrome

Phytochrome photoreceptors in plants and microorganisms switch photochromically between two states, controlling numerous important biological processes. Although this phototransformation is generally considered to involve rotation of ring D of the tetrapyrrole chromophore, Ulijasz et al. (Ulijasz, A. T., Cornilescu, G., Cornilescu, C. C., Zhang, J., Rivera, M., Markley, J. L., and Vierstra, R. D. (2010) Nature 463, 250-254) proposed that the A-ring rotates instead. Here, we apply magic angle spinning NMR to the two parent states following studies of the 23-kDa GAF (cGMP phosphodiesterase/adenylyl cyclase/FhlA) domain fragment of phytochrome from Synechococcus OS-B'. Major changes occur at the A-ring covalent linkage to the protein as well as at the protein residue contact of ring D. Conserved contacts associated with the A-ring nitrogen rule out an A-ring photoflip, whereas loss of contact of the D-ring nitrogen to the protein implies movement of ring D. Although none of the methine bridges showed a chemical shift change comparable with those characteristic of the D-ring photoflip in canonical phytochromes, denaturation experiments showed conclusively that the same occurs in Synechococcus OS-B' phytochrome upon photoconversion. The results are consistent with the D-ring being strongly tilted in both states and the C15=C16 double bond undergoing a Z/E isomerization upon light absorption. More subtle changes are associated with the A-ring linkage to the protein. Our findings thus disprove A-ring rotation and are discussed in relation to the position of the D-ring, photoisomerization, and photochromicity in the phytochrome family.

BCov: a method for predicting β-sheet topology using sparse inverse covariance estimation and integer programming

Prediction of protein residue contacts, even at the coarse-grain level, can help in finding solutions to the protein structure prediction problem. Unlike α-helices that are locally stabilized, β-sheets result from pairwise hydrogen bonding of two or more disjoint regions of the protein backbone. The problem of predicting contacts among β-strands in proteins has been addressed by several supervised computational approaches. Recently, prediction of residue contacts based on correlated mutations has been greatly improved and finally allows the prediction of 3D structures of the proteins. In this article, we describe BCov, which is the first unsupervised method to predict the β-sheet topology starting from the protein sequence and its secondary structure. BCov takes advantage of the sparse inverse covariance estimation to define β-strand partner scores. Then an optimization based on integer programming is carried out to predict the β-sheet connectivity. When tested on the prediction of β-strand pairing, BCov scores with average values of Matthews Correlation Coefficient (MCC) and F1 equal to 0.56 and 0.61, respectively, on a non-redundant dataset of 916 protein chains known with atomic resolution. Our approach well compares with the state-of-the-art methods trained so far for this specific task. The method is freely available under General Public License at http://biocomp.unibo.it/savojard/bcov/bcov-1.0.tar.gz. The new dataset BetaSheet1452 can be downloaded at http://biocomp.unibo.it/savojard/bcov/BetaSheet1452.dat.

Non-local residue–residue contacts in proteins are more conserved than local ones

Non-covalent residue-residue contacts drive the folding of proteins and stabilize them. They may be local-i.e. involve residues that are close in sequence, or non-local. It has been suggested that, in most proteins, local contacts drive protein folding by providing crucial constraints of the conformational space, thus allowing proteins to fold. We compared residues that are involved in local contacts to residues that are involved in non-local contacts and found that, in most proteins, residues in non-local contacts are significantly more conserved evolutionarily than residues in local contacts. Moreover, non-local contacts are more structurally conserved: a contact between positions that are distant in sequence is more likely to exist in many structural homologues compared with a contact between positions that are close in sequence. These results provide new insights into the mechanisms of protein folding and may allow for better prediction of critical intra-chain contacts.

Prediction of Protein Folding Rates from Structural Topology and Complex Network Properties

As a fundamental biological problem, revealing the protein folding mechanism remains to be one of the most challenging problems in structural bioinformatics. Prediction of protein folding rate is an important step towards our further understanding of the protein folding mechanism and the complex sequence-structure-function relationship. In this article, we develop a novel approach to predict protein folding rates for two-state and multi-state protein folding kinetics, which combines a variety of structural topology and complex network properties that are calculated from protein three-dimensional structures. To take into account the specific correlations between network properties and protein folding rates, we define two different protein residue contact networks, based on two different scales Protein Contact Network (PCN) and Long-range Interaction Network (LIN) to characterize the corresponding network features. The leave-one-out cross-validation (LOOCV) tests indicate that this integrative strategy is more powerful in predicting the folding rates from 3D structures, with the Pearson's Correlation Coefficient (CC) of 0.88, 0.90 and 0.90 for two-state, multi-state and combined protein folding kinetics, which provides an improved performance compared with other prediction work. This study provides useful insights which shed light on the network organization of interacting residues underlying protein folding process for both two-state and multi-state folding kinetics. Moreover, our method also provides a complementary approach to the current folding rate prediction algorithms and can be used as a powerful tool for the characterization of the foldomics protein data. The implemented webserver (termed PRORATE) is freely accessible at http://sunflower.kuicr.kyoto-u.ac.jp/~sjn/folding/.