SpotOn: High Accuracy Identification of Protein-Protein Interface Hot-Spots

Computational Screening for mAb Colloidal Stability with Coarse-Grained, Molecular-Scale Simulations

Monoclonal antibodies (mAbs) are an important modality of protein therapeutics with broad applications for numerous diseases. However, colloidal instabilities occurring at high protein concentrations can limit the ability to develop stable, high-concentration liquid dosage forms that are required for patient-centric, device-mediated products. Therefore, it is advantageous to identify colloidally stable mAbs early in the discovery process to ensure that they are selected for development. Experimental screening for colloidal stability can be time- and resource-consuming and is most feasible at the later stages of drug development due to material requirements. Alternatively, computational approaches have emerging potential to provide efficient screening and focus developmental efforts on mAbs with the greatest developability potential, while providing mechanistic relationships for colloidal instability. In this work, coarse-grained, molecular-scale models were fine-tuned to screen for colloidal stability at amino-acid resolution. This model parameterization provides a framework to screen for mAb self-interactions and extrapolate to bulk solution behavior. This approach was applied to a wide array of mAbs under multiple buffer conditions, demonstrating the utility of the presented computational approach to augment early candidate screening and later formulation strategies for protein therapeutics.

Software and Databases for Protein-Protein Docking

Protein-protein interactions (PPIs) provide valuable insights for understanding the principles of biological systems and for elucidating causes of incurable diseases. One of the techniques used for computational prediction of PPIs is protein-protein docking calculations, and a variety of software has been developed. This chapter is a summary of software and databases used for protein-protein docking.

Development and validation of a clinical prediction model for glioma grade using machine learning

BACKGROUND

Histopathological evaluation is currently the gold standard for grading gliomas; however, this technique is invasive.

OBJECTIVE

This study aimed to develop and validate a diagnostic prediction model for glioma by employing multiple machine learning algorithms to identify risk factors associated with high-grade glioma, facilitating the prediction of glioma grading.

METHODS

Data from 1114 eligible glioma patients were obtained from The Cancer Genome Atlas (TCGA) database, which was divided into a training set (n= 781) and a test set (n= 333). Fifty machine learning algorithms were employed, and the optimal algorithm was selected to construct a prediction model. The performance of the machine learning prediction model was compared to the clinical prediction model in terms of discrimination, calibration, and clinical validity to assess the performance of the prediction model.

RESULTS

The area under the curve (AUC) values of the machine learning prediction models (training set: 0.870 vs. 0.740, test set: 0.863 vs. 0.718) were significantly improved from the clinical prediction models. Furthermore, significant improvement in discrimination was observed for the Integrated Discrimination Improvement (IDI) (training set: 0.230, test set: 0.270) and Net Reclassification Index (NRI) (training set: 0.170, test set: 0.170) from the clinical prognostic model. Both models showed a high goodness of fit and an increased net benefit.

CONCLUSION

A strong prediction accuracy model can be developed using machine learning algorithms to screen for high-grade glioma risk predictors, which can serve as a non-invasive prediction tool for preoperative diagnostic grading of glioma.

Distinct or Overlapping Areas of Mitochondrial Thioredoxin 2 May Be Used for Its Covalent and Strong Non-Covalent Interactions with Protein Ligands

In silico approaches were employed to examine the characteristics of interactions between human mitochondrial thioredoxin 2 (HsTrx2) and its 38 previously identified mitochondrial protein ligands. All interactions appeared driven mainly by electrostatic forces. The statistically significant residues of HsTrx2 for interactions were characterized as "contact hot spots". Since these were identical/adjacent to putative thermodynamic hot spots, an energy network approach identified their neighbors to highlight possible contact interfaces. Three distinct areas for binding emerged: (i) one around the active site for covalent interactions, (ii) another antipodal to the active site for strong non-covalent interactions, and (iii) a third area involved in both kinds of interactions. The contact interfaces of HsTrx2 were projected as respective interfaces for Escherichia coli Trx1 (EcoTrx1), 2, and HsTrx1. Comparison of the interfaces and contact hot spots of HsTrx2 to the contact residues of EcoTx1 and HsTrx1 from existing crystal complexes with protein ligands supported the hypothesis, except for a part of the cleft/groove adjacent to Trp30 preceding the active site. The outcomes of this study raise the possibility for the rational design of selective inhibitors for the interactions of HsTrx2 with specific protein ligands without affecting the entirety of the functions of the Trx system.

Role of the redox state of the Pirin-bound cofactor on interaction with the master regulators of inflammation and other pathways

Persistent cellular stress induced perpetuation and uncontrolled amplification of inflammatory response results in a shift from tissue repair toward collateral damage, significant alterations of tissue functions, and derangements of homeostasis which in turn can lead to a large number of acute and chronic pathological conditions, such as chronic heart failure, atherosclerosis, myocardial infarction, neurodegenerative diseases, diabetes, rheumatoid arthritis, and cancer. Keeping the vital role of balanced inflammation in maintaining tissue integrity in mind, the way to combating inflammatory diseases may be through identification and characterization of mediators of inflammation that can be targeted without hampering normal body function. Pirin (PIR) is a non-heme iron containing protein having two different conformations depending on the oxidation state of the iron. Through exploration of the Pirin interactome and using molecular docking approaches, we identified that the Fe2+-bound Pirin directly interacts with BCL3, NFKBIA, NFIX and SMAD9 with more resemblance to the native binding pose and higher affinity than the Fe3+-bound form. In addition, Pirin appears to have a function in the regulation of inflammation, the transition between the canonical and non-canonical NF-κB pathways, and the remodeling of the actin cytoskeleton. Moreover, Pirin signaling appears to have a critical role in tumor invasion and metastasis, as well as metabolic and neuro-pathological complications. There are regulatory variants in PIR that can influence expression of not only PIR but also other genes, including VEGFD and ACE2. Disparity exists between South Asian and European populations in the frequencies of variant alleles at some of these regulatory loci that may lead to differential occurrence of Pirin-mediated pathogenic conditions.

FPGA Implementation of IIR Notch and Anti-Notch Filters With an Application to Localization of Protein Hot-Spots

In this paper, high-speed second-order infinite impulse response (IIR) notch filter (NF) and anti-notch filter (ANF) are designed and realized on hardware. The improvement in speed of operation for the NF is then achieved by using the re-timing concept. The ANF is designed to specify a stability margin and minimize the amplitude area. Next, an improved approach is proposed for the detection of protein hot-spot locations using the designed second-order IIR ANF. The analytical and experimental results reported in this paper show that the proposed approach provides better hot-spot prediction compared to the reported classical filtering techniques based on the IIR Chebyshev filter and S-transform. The proposed approach also yields consistency in prediction hot-spots compared to the results based on biological methodologies. Furthermore, the presented technique reveals some new "potential" hot-spots. The proposed filters are simulated and synthesized using the Xilinx Vivado 18.3 software platform with Zynq-7000 Series (ZedBoard Zynq Evaluation and Development Kit xc7z020clg484-1) FPGA family.

Prediction of hot spots towards drug discovery by protein sequence embedding with 1D convolutional neural network

Protein hotspot residues are key sites that mediate protein-protein interactions. Accurate identification of these residues is essential for understanding the mechanism from protein to function and for designing drug targets. Current research has mostly focused on using machine learning methods to predict hot spots from known interface residues, which artificially extract the corresponding features of amino acid residues from sequence, structure, evolution, energy, and other information to train and test machine learning models. The process is cumbersome, time-consuming and laborious to some extent. This paper proposes a novel idea that develops a pre-trained protein sequence embedding model combined with a one-dimensional convolutional neural network, called Embed-1dCNN, to predict protein hotspot residues. In order to obtain large data samples, this work integrates and extracts data from the datasets of ASEdb, BID, SKEMPI and dbMPIKT to generate a new dataset, and adopts the SMOTE algorithm to expand positive samples to form the training set. The experimental results show that the method achieves an F1 score of 0.82 on the test set. Compared with other hot spot prediction methods, our model achieved better prediction performance.

Molecular determinants associated with temporal succession of SARS-CoV-2 variants in Uttar Pradesh, India

The emergence and rapid evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) caused a global crisis that required a detailed characterization of the dynamics of mutational pattern of the viral genome for comprehending its epidemiology, pathogenesis and containment. We investigated the molecular evolution of the SASR-CoV-2 genome during the first, second and third waves of COVID-19 in Uttar Pradesh, India. Nanopore sequencing of the SARS-CoV-2 genome was undertaken in 544 confirmed cases of COVID-19, which included vaccinated and unvaccinated individuals. In the first wave (unvaccinated population), the 20A clade (56.32%) was superior that was replaced by 21A Delta in the second wave, which was more often seen in vaccinated individuals in comparison to unvaccinated (75.84% versus 16.17%, respectively). Subsequently, 21A delta got outcompeted by Omicron (71.8%), especially the 21L variant, in the third wave. We noticed that Q677H appeared in 20A Alpha and stayed up to Delta, D614G appeared in 20A Alpha and stayed in Delta and Omicron variants (got fixed), and several other mutations appeared in Delta and stayed in Omicron. A cross-sectional analysis of the vaccinated and unvaccinated individuals during the second wave revealed signature combinations of E156G, F157Del, L452R, T478K, D614G mutations in the Spike protein that might have facilitated vaccination breach in India. Interestingly, some of these mutation combinations were carried forward from Delta to Omicron. In silico protein docking showed that Omicron had a higher binding affinity with the host ACE2 receptor, resulting in enhanced infectivity of Omicron over the Delta variant. This work has identified the combinations of key mutations causing vaccination breach in India and provided insights into the change of [virus's] binding affinity with evolution, resulting in more virulence in Delta and more infectivity in Omicron variants of SARS-CoV-2. Our findings will help in understanding the COVID-19 disease biology and guide further surveillance of the SARS-CoV-2 genome to facilitate the development of vaccines with better efficacies.

Structural modelling and dynamics of full-length of TLR10 sheds light on possible modes of dimerization, ligand binding and mechanism of action

Toll like receptors (TLRs) play a pivotal role in innate and adaptive immunity. There are 10 TLRs in the human genome, of which TLR10 is the least characterized. Genetic polymorphism of TLR10 has been shown to be associated with multiple diseases including tuberculosis and rheumatoid arthritis. TLR10 consists of an extracellular domain (ECD), a single-pass transmembrane (TM) helix and intracellular TIR (Toll/Interleukin-1 receptor) domain. ECD is employed for ligand recognition and the intracellular domain interacts with other TIR domain-containing adapter proteins for signal transduction. Experimental structure of ECD or TM domain is not available for TLR10. In this study, we have modelled multiple forms of TLR10-ECD dimers, such as closed and open forms, starting from available structures of homologues. Subsequently, multiple full-length TLR10 homodimer models were generated by utilizing homology modelling and protein-protein docking. The dynamics of these models in membrane-aqueous environment revealed the global motion of ECD and TIR domain towards membrane bilayer. The TIR domain residues exhibited high root mean square fluctuation compared to ECD. The 'closed form' model was observed to be energetically more favorable than 'open form' model. The evaluation of persistent interchain interactions, along with their conservation score, unveiled critical residues for each model. Further, the binding of dsRNA to TLR10 was modelled by defined and blind docking approaches. Differential binding of dsRNA to the protomers of TLR10 was observed upon simulation that could provide clues on ligand disassociation. Dynamic network analysis revealed that the 'open form' model can be the functional form while 'closed form' model can be the apo form of TLR10.

Modeling of Protein Complexes

The recent advances in structural biology, combined with continuously increasing computational capabilities and development of advanced softwares, have drastically simplified the workflow for protein homology modeling. Modeling of individual proteins is nowadays quick and straightforward for a large variety of protein targets, thanks to guided pipelines relying on advanced computational tools and user-friendly interfaces, which have extended and promoted the use of modeling also to scientists not focusing on molecular structures of proteins. Nevertheless, construction of models of multi-protein complexes remains quite challenging for the non-experts, often due to the usage of specific procedures depending on the system under investigation and the need for experimental validation approaches to strengthen the generated output.In this chapter, we provide a brief overview of the approaches enabling generation of multi-protein complex models starting from homology models of individual protein components. Using real-life examples, we include two examples to guide the reader in the generation of homomeric and heteromeric protein models.

Large interfacial relocation in RBD-ACE2 complex may explain fast-spreading property of Omicron

The Omicron variant of SARS-CoV-2 emerged in South African in late 2021. This variant has a large number of mutations, and regarded as fastest-spreading Covid variant. The spike RBD region of SARS-CoV-2 and its interaction with human ACE2 play fundamental role in viral infection and transmission. To explore the reason of fast-spreading properties of Omicron variant, we have modeled the interactions of Omicron RBD and human ACE2 using docking and molecular dynamics simulations. Results show that RBD-ACE2 binding site may drastically relocate with an enlarged interface. The predicted interface has large negative binding energies and shows stable conformation in molecular dynamics simulations. It was found that the interfacial area in Omicron RBD-ACE2 complex is increased up to 40% in comparison to wild-type Sars-Cov-2. Moreover, the number of hydrogen bonds significantly increased up to 80%. The key interacting residues become also very different in Omicron variant. The new binding interface can significantly accommodate R403, as a key RBD residue, near ACE2 surface which leads to two new strong salt bridges. The exploration of the new binding interface can help to understand the reasons of high transmission rate of Omicron.

Computational Docking Reveals Co-Evolution of C4 Carbon Delivery Enzymes in Diverse Plants

Proteins are modular functionalities regulating multiple cellular activities in prokaryotes and eukaryotes. As a consequence of higher plants adapting to arid and thermal conditions, C4 photosynthesis is the carbon fixation process involving multi-enzymes working in a coordinated fashion. However, how these enzymes interact with each other and whether they co-evolve in parallel to maintain interactions in different plants remain elusive to date. Here, we report our findings on the global protein co-evolution relationship and local dynamics of co-varying site shifts in key C4 photosynthetic enzymes. We found that in most of the selected key C4 photosynthetic enzymes, global pairwise co-evolution events exist to form functional couplings. Besides, protein-protein interactions between these enzymes may suggest their unknown functionalities in the carbon delivery process. For PEPC and PPCK regulation pairs, pocket formation at the interactive interface are not necessary for their function. This feature is distinct from another well-known regulation pair in C4 photosynthesis, namely, PPDK and PPDK-RP, where the pockets are necessary. Our findings facilitate the discovery of novel protein regulation types and contribute to expanding our knowledge about C4 photosynthesis.

Exploring the Surface of the Ectodomain of the PD-L1 Immune Checkpoint with Small-Molecule Fragments

Development of small molecules targeting the PD-L1/PD-1 interface is advancing both in industry and academia, but only a few have reached early-stage clinical trials. Here, we take a closer look at the general druggability of PD-L1 using in silico hot spot mapping and nuclear magnetic resonance (NMR)-based characterization. We found that the conformational elasticity of the PD-L1 surface strongly influences the formation of hot spots. We deconstructed several generations of known inhibitors into fragments and examined their binding properties using differential scanning fluorimetry (DSF) and protein-based nuclear magnetic resonance (NMR). These biophysical analyses showed that not all fragments bind to the PD-L1 ectodomain despite having the biphenyl scaffold. Although most of the binding fragments induced PD-L1 oligomerization, two compounds, TAH35 and TAH36, retain the monomeric state of proteins upon binding. Additionally, the presence of the entire ectodomain did not affect the binding of the hit compounds and dimerization of PD-L1. The data demonstrated here provide important information on the PD-L1 druggability and the structure-activity relationship of the biphenyl core moiety and therefore may aid in the design of novel inhibitors and focused fragment libraries for PD-L1.

PPI-HotspotDB: Database of Protein-Protein Interaction Hot Spots

Single-point mutations of certain residues (so-called hot spots) impair/disrupt protein-protein interactions (PPIs), leading to pathogenesis and drug resistance. Conventionally, a PPI-hot spot is identified when its replacement decreased the binding free energy significantly, generally by ≥2 kcal/mol. The relatively few mutations with such a significant binding free energy drop limited the number of distinct PPI-hot spots. By defining PPI-hot spots based on mutations that have been manually curated in UniProtKB to significantly impair/disrupt PPIs in addition to binding free energy changes, we have greatly expanded the number of distinct PPI-hot spots by an order of magnitude. These experimentally determined PPI-hot spots along with available structures have been collected in a database called PPI-Hotspot^DB. We have applied the PPI-Hotspot^DB to create a nonredundant benchmark, PPI-Hotspot+PDB^BM, for assessing methods to predict PPI-hot spots using the free structure as input. PPI-Hotspot^DB will benefit the design of mutagenesis experiments and development of PPI-hot spot prediction methods. The database and benchmark are freely available at https://ppihotspot.limlab.dnsalias.org.

Big data and artificial intelligence (AI) methodologies for computer-aided drug design (CADD)

There have been numerous advances in the development of computational and statistical methods and applications of big data and artificial intelligence (AI) techniques for computer-aided drug design (CADD). Drug design is a costly and laborious process considering the biological complexity of diseases. To effectively and efficiently design and develop a new drug, CADD can be used to apply cutting-edge techniques to various limitations in the drug design field. Data pre-processing approaches, which clean the raw data for consistent and reproducible applications of big data and AI methods are introduced. We include the current status of the applicability of big data and AI methods to drug design areas such as the identification of binding sites in target proteins, structure-based virtual screening (SBVS), and absorption, distribution, metabolism, excretion and toxicity (ADMET) property prediction. Data pre-processing and applications of big data and AI methods enable the accurate and comprehensive analysis of massive biomedical data and the development of predictive models in the field of drug design. Understanding and analyzing biological, chemical, or pharmaceutical architectures of biomedical entities related to drug design will provide beneficial information in the biomedical big data era.

Artificial intelligence based methods for hot spot prediction

Proteins interact through their interfaces to fulfill essential functions in the cell. They bind to their partners in a highly specific manner and form complexes that have a profound effect on understanding the biological pathways they are involved in. Any abnormal interactions may cause diseases. Therefore, the identification of small molecules which modulate protein interactions through their interfaces has high therapeutic potential. However, discovering such molecules is challenging. Most protein-protein binding affinity is attributed to a small set of amino acids found in protein interfaces known as hot spots. Recent studies demonstrate that drug-like small molecules specifically may bind to hot spots. Therefore, hot spot prediction is crucial. As experimental data accumulates, artificial intelligence begins to be used for computational hot spot prediction. First, we review machine learning and deep learning for computational hot spot prediction and then explain the significance of hot spots toward drug design.

Native or Non-Native Protein-Protein Docking Models? Molecular Dynamics to the Rescue

Molecular docking excels at creating a plethora of potential models of protein-protein complexes. To correctly distinguish the favorable, native-like models from the remaining ones remains, however, a challenge. We assessed here if a protocol based on molecular dynamics (MD) simulations would allow distinguishing native from non-native models to complement scoring functions used in docking. To this end, the first models for 25 protein-protein complexes were generated using HADDOCK. Next, MD simulations complemented with machine learning were used to discriminate between native and non-native complexes based on a combination of metrics reporting on the stability of the initial models. Native models showed higher stability in almost all measured properties, including the key ones used for scoring in the Critical Assessment of PRedicted Interaction (CAPRI) competition, namely the positional root mean square deviations and fraction of native contacts from the initial docked model. A random forest classifier was trained, reaching a 0.85 accuracy in correctly distinguishing native from non-native complexes. Reasonably modest simulation lengths of the order of 50-100 ns are sufficient to reach this accuracy, which makes this approach applicable in practice.

Identification of invariant chain CD74 as a functional receptor of tissue inhibitor of metalloproteinases-1 (TIMP-1)

Multifunctionality of tissue inhibitor of metalloproteinases-1 (TIMP-1) comprising antiproteolytic as well as cytokinic activity has been attributed to its N-terminal and C-terminal domains, respectively. The molecular basis of the emerging proinflammatory cytokinic activity of TIMP-1 is still not completely understood. The cytokine receptor invariant chain (CD74) is involved in many inflammation-associated diseases and is highly expressed by immune cells. CD74 triggers zeta chain–associated protein kinase-70 (ZAP-70) signaling–associated activation upon interaction with its only known ligand, the macrophage migration inhibitory factor. Here, we demonstrate TIMP-1–CD74 interaction by coimmunoprecipitation and confocal microscopy in cells engineered to overexpress CD74. In silico docking in HADDOCK predicted regions of the N-terminal domain of TIMP-1 (N-TIMP-1) to interact with CD74. This was experimentally confirmed by confocal microscopy demonstrating that recombinant N-TIMP-1 lacking the entire C-terminal domain was sufficient to bind CD74. Interaction of TIMP-1 with endogenously expressed CD74 was demonstrated in the Namalwa B lymphoma cell line by dot blot binding assays as well as confocal microscopy. Functionally, we demonstrated that TIMP-1–CD74 interaction triggered intracellular ZAP-70 activation. N-TIMP-1 was sufficient to induce ZAP-70 activation and interference with the cytokine-binding site of CD74 using a synthetic peptide–abrogated TIMP-1-mediated ZAP-70 activation. Altogether, we here identified CD74 as a receptor and mediator of cytokinic TIMP-1 activity and revealed TIMP-1 as moonlighting protein harboring both cytokinic and antiproteolytic activity within its N-terminal domain. Recognition of this functional TIMP-1–CD74 interaction may shed new light on clinical attempts to therapeutically target ligand-induced CD74 activity in cancer and other inflammatory diseases.

Protein-protein interactions at a glance: Protocols for the visualization of biomolecular interactions

Protein-protein interactions (PPIs) play a key role in many biological processes and are intriguing targets for drug discovery campaigns. Advancements in experimental and computational techniques are leading to a growth of data accessibility, and, with it, an increased need for the analysis of PPIs. In this respect, visualization tools are essential instruments to represent and analyze biomolecular interactions. In this chapter, we reviewed some of the available tools, highlighting their features, and describing their functions with practical information on their usage.

Structural Biology in the Clouds: The WeNMR-EOSC Ecosystem

Structural biology aims at characterizing the structural and dynamic properties of biological macromolecules at atomic details. Gaining insight into three dimensional structures of biomolecules and their interactions is critical for understanding the vast majority of cellular processes, with direct applications in health and food sciences. Since 2010, the WeNMR project (www.wenmr.eu) has implemented numerous web-based services to facilitate the use of advanced computational tools by researchers in the field, using the high throughput computing infrastructure provided by EGI. These services have been further developed in subsequent initiatives under H2020 projects and are now operating as Thematic Services in the European Open Science Cloud portal (www.eosc-portal.eu), sending >12 millions of jobs and using around 4,000 CPU-years per year. Here we review 10 years of successful e-infrastructure solutions serving a large worldwide community of over 23,000 users to date, providing them with user-friendly, web-based solutions that run complex workflows in structural biology. The current set of active WeNMR portals are described, together with the complex backend machinery that allows distributed computing resources to be harvested efficiently.

Integrating Resonant Recognition Model and Stockwell Transform for Localization of Hotspots in Tubulin

Tubulin is a promising target for designing anti-cancer drugs. Identification of hotspots in multifunctional Tubulin protein provides insights for new drug discovery. Although machine learning techniques have shown significant results in prediction, they fail to identify the hotspots corresponding to a particular biological function. This paper presents a signal processing technique combining resonant recognition model (RRM) and Stockwell Transform (ST) for the identification of hotspots corresponding to a particular functionality. The characteristic frequency (CF) representing a specific biological function is determined using the RRM. Then the spectrum of the protein sequence is computed using ST. The CF is filtered from the ST spectrum using a time-frequency mask. The energy peaks in the filtered sequence represent the hotspots. The hotspots predicted by the proposed method are compared with the experimentally detected binding residues of Tubulin stabilizing drug Taxol and destabilizing drug Colchicine present in the Tubulin protein. Out of the 53 experimentally identified hotspots, 60% are predicted by the proposed method whereas around 20% are predicted by existing machine learning based methods. Additionally, the proposed method predicts some new hot spots, which may be investigated.

Centrality Measures in Residue Interaction Networks to Highlight Amino Acids in Protein–Protein Binding

Residue interaction networks (RINs) describe a protein structure as a network of interacting residues. Central nodes in these networks, identified by centrality analyses, highlight those residues that play a role in the structure and function of the protein. However, little is known about the capability of such analyses to identify residues involved in the formation of macromolecular complexes. Here, we performed six different centrality measures on the RINs generated from the complexes of the SKEMPI 2 database of changes in protein–protein binding upon mutation in order to evaluate the capability of each of these measures to identify major binding residues. The analyses were performed with and without the crystallographic water molecules, in addition to the protein residues. We also investigated the use of a weight factor based on the inter-residue distances to improve the detection of these residues. We show that for the identification of major binding residues, closeness, degree, and PageRank result in good precision, whereas betweenness, eigenvector, and residue centrality analyses give a higher sensitivity. Including water in the analysis improves the sensitivity of all measures without losing precision. Applying weights only slightly raises the sensitivity of eigenvector centrality analysis. We finally show that a combination of multiple centrality analyses is the optimal approach to identify residues that play a role in protein–protein interaction.

Computational Analysis of Domains Vulnerable to HPV-16 E6 Oncoprotein and Corresponding Hot Spot Residues

BACKGROUND

Human Papilloma Virus (HPV) is the primary cause of cancers in cervix, head and neck regions. Oncoprotein E6 of HPV-16, after infecting human body, alters host protein- protein interaction networks. E6 interacts with several proteins, causing the infection to progress into cervical cancer. The molecular basis for these interactions is the presence of short linear peptide motifs on E6 identical to those on human proteins.

METHODS

Motifs of LXXLL and E/DLLL/V-G after identification on E6, were analyzed for their dynamic fluctuations by use of elastic network models. Correlation analysis of amino acid residues of E6 was also performed in specific regions of motifs.

RESULTS

Arginine, Leucine, Glutamine, Threonine and Glutamic acid have been identified as hot spot residues of E6 which can subsequently provide a platform for drug designing and understanding of pathogenesis of cervical cancer. These amino acids play a significant role in stabilizing interactions with host proteins, ultimately causing infections and cancers.

CONCLUSION

Our study validates the role of linear binding motifs of E6 of HPV in interacting with these proteins as an important event in the propagation of HPV in human cells and its transformation into cervical cancer. The study further predicts the domains of protein kinase and armadillo as part of the regions involved in the interaction of E6AP, Paxillin and TNF R1, with viral E6.

SAAMBE-SEQ: a sequence-based method for predicting mutation effect on protein-protein binding affinity

MOTIVATION

Vast majority of human genetic disorders are associated with mutations that affect protein-protein interactions by altering wild-type binding affinity. Therefore, it is extremely important to assess the effect of mutations on protein-protein binding free energy to assist the development of therapeutic solutions. Currently, the most popular approaches use structural information to deliver the predictions, which precludes them to be applicable on genome-scale investigations. Indeed, with the progress of genomic sequencing, researchers are frequently dealing with assessing effect of mutations for which there is no structure available.

RESULTS

Here, we report a Gradient Boosting Decision Tree machine learning algorithm, the SAAMBE-SEQ, which is completely sequence-based and does not require structural information at all. SAAMBE-SEQ utilizes 80 features representing evolutionary information, sequence-based features and change of physical properties upon mutation at the mutation site. The approach is shown to achieve Pearson correlation coefficient (PCC) of 0.83 in 5-fold cross validation in a benchmarking test against experimentally determined binding free energy change (ΔΔG). Further, a blind test (no-STRUC) is compiled collecting experimental ΔΔG upon mutation for protein complexes for which structure is not available and used to benchmark SAAMBE-SEQ resulting in PCC in the range of 0.37-0.46. The accuracy of SAAMBE-SEQ method is found to be either better or comparable to most advanced structure-based methods. SAAMBE-SEQ is very fast, available as webserver and stand-alone code, and indeed utilizes only sequence information, and thus it is applicable for genome-scale investigations to study the effect of mutations on protein-protein interactions.

AVAILABILITY AND IMPLEMENTATION

SAAMBE-SEQ is available at http://compbio.clemson.edu/saambe_webserver/indexSEQ.php#started.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combining Experimental Isotherms, Minimalistic Simulations, and a Model to Understand and Predict Chemical Adsorption onto Montmorillonite Clays

An attractive approach to minimize human and animal exposures to toxic environmental contaminants is the use of safe and effective sorbent materials to sequester them. Montmorillonite clays have been shown to tightly bind diverse toxic chemicals. Due to their promise as sorbents to mitigate chemical exposures, it is important to understand their function and rapidly screen and predict optimal clay-chemical combinations for further testing. We derived adsorption free-energy values for a structurally and physicochemically diverse set of toxic chemicals using experimental adsorption isotherms performed in the current and previous studies. We studied the diverse set of chemicals using minimalistic MD simulations and showed that their interaction energies with calcium montmorillonite clays calculated using simulation snapshots in combination with their net charge and their corresponding solvent's dielectric constant can be used as inputs to a minimalistic model to predict adsorption free energies in agreement with experiments. Additionally, experiments and computations were used to reveal structural and physicochemical properties associated with chemicals that can be adsorbed to calcium montmorillonite clay. These properties include positively charged groups, phosphine groups, halide-rich moieties, hydrogen bond donor/acceptors, and large, rigid structures. The combined experimental and computational approaches used in this study highlight the importance and potential applicability of analogous methods to study and design novel advanced sorbent systems in the future, broadening their applicability for environmental contaminants.

Hot spots localization in proteins by optimized short time Ramanujan Fourier transform

Specific functions in biological processes are dependent on protein-protein interactions. Hot spot residues play a key role in the determination of these interactions and have wide applications in engineering proteins and drug discovery. Experimental techniques to identify hotspots are often labor intensive and expensive. Also, most of the computational methods which have been developed are structure based and need some training. In this work, hotspots have been identified by sequence information alone using the Resonant Recognition Model (RRM). The proposed method uses characteristic period in place of traditionally used characteristic frequency by RRM-based methods. The characteristic period has been extracted from the consensus spectrum of protein families using the Ramanujan Fourier Transform (RFT). Position-period plots for proteins have been generated using Short Time RFT (ST-RFT) with a Gaussian window. Hot spots have been identified by thresholding of the signal corresponding to the protein's characteristic period in the ST-RFT. To enhance the performance of the ST-RFT, Gaussian window shape parameter has been optimized using concentration measure as a metric. Better sensitivity of this method has been observed compared to other reported RRM-based methods. Since the method is model independent it does not requires any training and can be readily used for any protein sequence provided its interface residues and protein family are known.

MENSAdb: a thorough structural analysis of membrane protein dimers

Membrane proteins (MPs) are key players in a variety of different cellular processes and constitute the target of around 60% of all Food and Drug Administration-approved drugs. Despite their importance, there is still a massive lack of relevant structural, biochemical and mechanistic information mainly due to their localization within the lipid bilayer. To help fulfil this gap, we developed the MEmbrane protein dimer Novel Structure Analyser database (MENSAdb). This interactive web application summarizes the evolutionary and physicochemical properties of dimeric MPs to expand the available knowledge on the fundamental principles underlying their formation. Currently, MENSAdb contains features of 167 unique MPs (63% homo- and 37% heterodimers) and brings insights into the conservation of residues, accessible solvent area descriptors, average B-factors, intermolecular contacts at 2.5 Å and 4.0 Å distance cut-offs, hydrophobic contacts, hydrogen bonds, salt bridges, π-π stacking, T-stacking and cation-π interactions. The regular update and organization of all these data into a unique platform will allow a broad community of researchers to collect and analyse a large number of features efficiently, thus facilitating their use in the development of prediction models associated with MPs. Database URL: http://www.moreiralab.com/resources/mensadb.

The UBA domain of conjugating enzyme Ubc1/Ube2K facilitates assembly of K48/K63-branched ubiquitin chains

The assembly of a specific polymeric ubiquitin chain on a target protein is a key event in the regulation of numerous cellular processes. Yet, the mechanisms that govern the selective synthesis of particular polyubiquitin signals remain enigmatic. The homologous ubiquitin-conjugating (E2) enzymes Ubc1 (budding yeast) and Ube2K (mammals) exclusively generate polyubiquitin linked through lysine 48 (K48). Uniquely among E2 enzymes, Ubc1 and Ube2K harbor a ubiquitin-binding UBA domain with unknown function. We found that this UBA domain preferentially interacts with ubiquitin chains linked through lysine 63 (K63). Based on structural modeling, in vitro ubiquitination experiments, and NMR studies, we propose that the UBA domain aligns Ubc1 with K63-linked polyubiquitin and facilitates the selective assembly of K48/K63-branched ubiquitin conjugates. Genetic and proteomics experiments link the activity of the UBA domain, and hence the formation of this unusual ubiquitin chain topology, to the maintenance of cellular proteostasis.

Robust principal component analysis-based prediction of protein-protein interaction hot spots

Proteins often exert their function by binding to other cellular partners. The hot spots are key residues for protein-protein binding. Their identification may shed light on the impact of disease associated mutations on protein complexes and help design protein-protein interaction inhibitors for therapy. Unfortunately, current machine learning methods to predict hot spots, suffer from limitations caused by gross errors in the data matrices. Here, we present a novel data pre-processing pipeline that overcomes this problem by recovering a low rank matrix with reduced noise using Robust Principal Component Analysis. Application to existing databases shows the predictive power of the method.

SAAFEC-SEQ: A Sequence-Based Method for Predicting the Effect of Single Point Mutations on Protein Thermodynamic Stability

Modeling the effect of mutations on protein thermodynamics stability is useful for protein engineering and understanding molecular mechanisms of disease-causing variants. Here, we report a new development of the SAAFEC method, the SAAFEC-SEQ, which is a gradient boosting decision tree machine learning method to predict the change of the folding free energy caused by amino acid substitutions. The method does not require the 3D structure of the corresponding protein, but only its sequence and, thus, can be applied on genome-scale investigations where structural information is very sparse. SAAFEC-SEQ uses physicochemical properties, sequence features, and evolutionary information features to make the predictions. It is shown to consistently outperform all existing state-of-the-art sequence-based methods in both the Pearson correlation coefficient and root-mean-squared-error parameters as benchmarked on several independent datasets. The SAAFEC-SEQ has been implemented into a web server and is available as stand-alone code that can be downloaded and embedded into other researchers’ code.

Decoding Protein-protein Interactions: An Overview

Drug discovery has focused on the paradigm "one drug, one target" for a long time. However, small molecules can act at multiple macromolecular targets, which serves as the basis for drug repurposing. In an effort to expand the target space, and given advances in X-ray crystallography, protein-protein interactions have become an emerging focus area of drug discovery enterprises. Proteins interact with other biomolecules and it is this intricate network of interactions that determines the behavior of the system and its biological processes. In this review, we briefly discuss networks in disease, followed by computational methods for protein-protein complex prediction. Computational methodologies and techniques employed towards objectives such as protein-protein docking, protein-protein interactions, and interface predictions are described extensively. Docking aims at producing a complex between proteins, while interface predictions identify a subset of residues on one protein that could interact with a partner, and protein-protein interaction sites address whether two proteins interact. In addition, approaches to predict hot spots and binding sites are presented along with a representative example of our internal project on the chemokine CXC receptor 3 B-isoform and predictive modeling with IP10 and PF4.

Predicting Hot Spots Using a Deep Neural Network Approach

Targeting protein-protein interactions is a challenge and crucial task of the drug discovery process. A good starting point for rational drug design is the identification of hot spots (HS) at protein-protein interfaces, typically conserved residues that contribute most significantly to the binding. In this chapter, we depict point-by-point an in-house pipeline used for HS prediction using only sequence-based features from the well-known SpotOn dataset of soluble proteins (Moreira et al., Sci Rep 7:8007, 2017), through the implementation of a deep neural network. The presented pipeline is divided into three steps: (1) feature extraction, (2) deep learning classification, and (3) model evaluation. We present all the available resources, including code snippets, the main dataset, and the free and open-source modules/packages necessary for full replication of the protocol. The users should be able to develop an HS prediction model with accuracy, precision, recall, and AUROC of 0.96, 0.93, 0.91, and 0.86, respectively.

Developing a machine learning model to identify protein–protein interaction hotspots to facilitate drug discovery

Throughout the history of drug discovery, an enzymatic-based approach for identifying new drug molecules has been primarily utilized. Recently, protein–protein interfaces that can be disrupted to identify small molecules that could be viable targets for certain diseases, such as cancer and the human immunodeficiency virus, have been identified. Existing studies computationally identify hotspots on these interfaces, with most models attaining accuracies of ~70%. Many studies do not effectively integrate information relating to amino acid chains and other structural information relating to the complex. Herein, (1) a machine learning model has been created and (2) its ability to integrate multiple features, such as those associated with amino-acid chains, has been evaluated to enhance the ability to predict protein–protein interface hotspots. Virtual drug screening analysis of a set of hotspots determined on the EphB2-ephrinB2 complex has also been performed. The predictive capabilities of this model offer an AUROC of 0.842, sensitivity/recall of 0.833, and specificity of 0.850. Virtual screening of a set of hotspots identified by the machine learning model developed in this study has identified potential medications to treat diseases caused by the overexpression of the EphB2-ephrinB2 complex, including prostate, gastric, colorectal and melanoma cancers which are linked to EphB2 mutations. The efficacy of this model has been demonstrated through its successful ability to predict drug-disease associations previously identified in literature, including cimetidine, idarubicin, pralatrexate for these conditions. In addition, nadolol, a beta blocker, has also been identified in this study to bind to the EphB2-ephrinB2 complex, and the possibility of this drug treating multiple cancers is still relatively unexplored.

Cyclic Peptides as Drugs for Intracellular Targets: The Next Frontier in Peptide Therapeutic Development

Developing macrocyclic peptides that can reach intracellular targets is a significant challenge. This review discusses the most recent strategies used to develop cell permeable cyclic peptides that maintain binding to their biological target inside the cell. Macrocyclic peptides are unique from small molecules because traditional calculated physical properties are unsuccessful for predicting cell membrane permeability. Peptide synthesis and experimental membrane permeability is the only strategy that effectively differentiates between cell permeable and cell impermeable molecules. Discussed are chemical strategies, including backbone N-methylation and stereochemical changes, which have produced molecular scaffolds with improved cell permeability. However, these improvements often come at the expense of biological activity as chemical modifications alter the peptide conformation, frequently impacting the compound's ability to bind to the target. Highlighted is the most promising approach, which involves side-chain alterations that improve cell permeability without impact binding events.

Molecular analysis of a fibrin-degrading enzyme from Bacillus subtilis K2 isolated from the Indonesian soybean-based fermented food moromi

The screening of proteolytic and fibrinolytic bacteria from moromi (an Indonesian soybean-based fermented food) yielded a number of isolates. Based on morphological and biochemical analyses and sequencing of the 16S rRNA gene, the isolate that exhibited the highest proteolytic and fibrinolytic activity was identified as Bacillus subtilis K2. The study was performed to analyze molecular characteristic of a fibrin-degrading enzyme from B. subtilis K2. BLASTn analysis of the nucleotide sequence encoding this fibrinolytic protein demonstrated 73.6% homology with the gene encoding the fibrin-degrading enzyme nattokinase of the B. subtilis subsp. natto, which was isolated from fermented soybean in Japan. An analysis of the putative amino-acid sequence of this protein indicated that it is a serine protease enzyme with aspartate, histidine, and serine in the catalytic triad. This enzyme was determined to be a 26-kDa molecule, as confirmed with a zymogram assay. Further bioinformatic analysis using Protparam demonstrated that the enzyme has a pI of 6.02, low instability index, high aliphatic index, and low GRAVY value. Molecular docking analysis using HADDOCK indicated that there are favorable interactions between subtilisin K2 and the fibrin substrate, as demonstrated by a high binding affinity (ΔG: - 19.4 kcal/mol) and low Kd value (6.3E-15 M). Overall, the study concluded that subtilisin K2 belong to serine protease enzyme has strong interactions with its fibrin substrate and fibrin can be rapidly degraded by this enzyme, suggesting its application as a treatment for thrombus diseases.

In silico analysis of ACE2 orthologues to predict animal host range with high susceptibility to SARS-CoV-2

SARS-CoV-2, which causes severe pneumonia epidemics, probably originated from Chinese horseshoe bats, but the intermediate and host range is still unknown. ACE2 is the entry receptor for SARS-CoV-2. The binding capacity of SARS-CoV-2 spike protein to ACE2 is the critical determinant of viral host range and cross-species infection. Here, we used an in silico approach to predict the potential animals range with high susceptibility to SARS-CoV-2 by modelling and studying the Spike-ACE2 interaction of 22 domestic and wild animals. Our results showed that all studied animals are potentially susceptible to SARS-CoV-2 infection with a slight difference in the binding affinity and stability of their ACE2-RBD complexes. Furthermore, we identified a specific substitution of tyrosine to histidine at position 41 in ACE2 that likely reduces the affinity to SARS-CoV-2 in horses and greater horseshoe bats. These results may help to provide important insights into SARS-CoV-2 host range which will make it possible to control the spread of the virus and identify animal models that could be used for screening antiviral drugs or vaccine candidates against SARS-CoV-2.

SPOTONE: Hot Spots on Protein Complexes with Extremely Randomized Trees via Sequence-Only Features

Protein Hot-Spots (HS) are experimentally determined amino acids, key to small ligand binding and tend to be structural landmarks on protein–protein interactions. As such, they were extensively approached by structure-based Machine Learning (ML) prediction methods. However, the availability of a much larger array of protein sequences in comparison to determined tree-dimensional structures indicates that a sequence-based HS predictor has the potential to be more useful for the scientific community. Herein, we present SPOTONE, a new ML predictor able to accurately classify protein HS via sequence-only features. This algorithm shows accuracy, AUROC, precision, recall and F1-score of 0.82, 0.83, 0.91, 0.82 and 0.85, respectively, on an independent testing set. The algorithm is deployed within a free-to-use webserver, only requiring the user to submit a FASTA file with one or more protein sequences.

Ab-initio binding of barnase–barstar with DelPhiForce steered Molecular Dynamics (DFMD) approach

Receptor–ligand interactions are involved in various biological processes, therefore understanding the binding mechanism and ability to predict the binding mode are essential for many biological investigations. While many computational methods exist to predict the 3D structure of the corresponding complex provided the knowledge of the monomers, here we use the newly developed DelPhiForce steered Molecular Dynamics (DFMD) approach to model the binding of barstar to barnase to demonstrate that first-principles methods are also capable of modeling the binding. Essential component of DFMD approach is enhancing the role of long-range electrostatic interactions to provide guiding force of the monomers toward their correct binding orientation and position. Thus, it is demonstrated that the DFMD can successfully dock barstar to barnase even if the initial positions and orientations of both are completely different from the correct ones. Thus, the electrostatics provides orientational guidance along with pulling force to deliver the ligand in close proximity to the receptor.

iPNHOT: a knowledge-based approach for identifying protein-nucleic acid interaction hot spots

Background

The interaction between proteins and nucleic acids plays pivotal roles in various biological processes such as transcription, translation, and gene regulation. Hot spots are a small set of residues that contribute most to the binding affinity of a protein-nucleic acid interaction. Compared to the extensive studies of the hot spots on protein-protein interfaces, the hot spot residues within protein-nucleic acids interfaces remain less well-studied, in part because mutagenesis data for protein-nucleic acids interaction are not as abundant as that for protein-protein interactions.

Results

In this study, we built a new computational model, iPNHOT, to effectively predict hot spot residues on protein-nucleic acids interfaces. One training data set and an independent test set were collected from dbAMEPNI and some recent literature, respectively. To build our model, we generated 97 different sequential and structural features and used a two-step strategy to select the relevant features. The final model was built based only on 7 features using a support vector machine (SVM). The features include two unique features such as ∆SASsa^1/2 and esp3, which are newly proposed in this study. Based on the cross validation results, our model gave F1 score and AUROC as 0.725 and 0.807 on the subset collected from ProNIT, respectively, compared to 0.407 and 0.670 of mCSM-NA, a state-of-the art model to predict the thermodynamic effects of protein-nucleic acid interaction. The iPNHOT model was further tested on the independent test set, which showed that our model outperformed other methods.

Conclusion

In this study, by collecting data from a recently published database dbAMEPNI, we proposed a new model, iPNHOT, to predict hotspots on both protein-DNA and protein-RNA interfaces. The results show that our model outperforms the existing state-of-art models. Our model is available for users through a webserver: http://zhulab.ahu.edu.cn/iPNHOT/.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Hot Spots and Their Contribution to the Self-Assembly of the Viral Capsid: In Silico Prediction and Analysis

The viral capsid is a macromolecular complex formed by a defined number of self-assembled proteins, which, in many cases, are biopolymers with an identical amino acid sequence. Specific protein-protein interactions (PPI) drive the capsid self-assembly process, leading to several distinct protein interfaces. Following the PPI hot spot hypothesis, we present a conservation-based methodology to identify those interface residues hypothesized to be crucial elements on the self-assembly and thermodynamic stability of the capsid. We validate the predictions through a rigorous physical framework which integrates molecular dynamics simulations and free energy calculations by Umbrella sampling and the potential of mean force using an all-atom molecular representation of the capsid proteins of an icosahedral virus in an explicit solvent. Our results show that a single mutation in any of the structure-conserved hot spots significantly perturbs the quaternary protein-protein interaction, decreasing the absolute value of the binding free energy, without altering the protein's secondary nor tertiary structure. Our conservation-based hot spot prediction methodology can lead to strategies to rationally modulate the capsid's thermodynamic properties.

Machine learning techniques for protein function prediction

Proteins play important roles in living organisms, and their function is directly linked with their structure. Due to the growing gap between the number of proteins being discovered and their functional characterization (in particular as a result of experimental limitations), reliable prediction of protein function through computational means has become crucial. This paper reviews the machine learning techniques used in the literature, following their evolution from simple algorithms such as logistic regression to more advanced methods like support vector machines and modern deep neural networks. Hyperparameter optimization methods adopted to boost prediction performance are presented. In parallel, the metamorphosis in the features used by these algorithms from classical physicochemical properties and amino acid composition, up to text-derived features from biomedical literature and learned feature representations using autoencoders, together with feature selection and dimensionality reduction techniques, are also reviewed. The success stories in the application of these techniques to both general and specific protein function prediction are discussed.

Discovery of novel helix binding sites at protein-protein interfaces

Protein-protein interactions (PPIs) play a key role in numerous biological processes. Many efforts have been undertaken to develop PPI modulators for therapeutic applications; however, to date, most of the peptide binders designed to target PPIs are derived from native binding helices or using the native helix binding site, which has limited the applications of protein-protein interface binding peptide design. Here, we developed a general computational algorithm, HPer (Helix Positioner), that locates single-helix binding sites at protein-protein interfaces based on the structure of protein targets. HPer performed well on known single-helix-mediated PPIs and recaptured the key interactions and hot-spot residues of native helical binders. We also screened non-helical-mediated PPIs in the PDBbind database and identified 17 PPIs that were suitable for helical peptide binding, and the helical binding sites in these PPIs were also predicted for designing novel peptide ligands. The L2 domain of EGFR, which was the top ranked, was selected as an example to show the protocol and results of designing novel helical peptide ligands on the searched binding site. The binding stability of the designed sequences were further investigated using molecular dynamics simulations.

Small protein-protein interfaces rich in electrostatic are often linked to regulatory function

Protein-protein interaction (PPI) is critical for several biological functions in living cells through the formation of an interface. Therefore, it is of interest to characterize protein-protein interfaces using an updated non-redundant structural dataset of 2557 homo (identical subunits) and 393 hetero (different subunits) dimer protein complexes determined by X-ray crystallography. We analyzed the interfaces using van der Waals (vdW), hydrogen bonding and electrostatic energies. Results show that on average homo and hetero interfaces are similar. Hence, we further grouped the 2950 interfaces based on percentage vdW to total energies into dominant (≥60%) and sub-dominant (<60%) vdW interfaces. Majority (92%) of interfaces have dominant vdW energy with large interface size (146 ± 87 (homo) and 137 ± 76 (hetero) residues) and interface area (1622 ± 1135 Ų (homo) and 1579 ± 1060 Ų (hetero)). However, a proportion (8%) of interfaces have sub-dominant vdW energy with small interface size (85 ± 46 (homo) and 88 ± 36 (hetero) residues) and interface area (823 ± 538 Ų (homo) and 881 ± 377 Ų (hetero)). It is found that large interfaces have two-fold more interface area and interface size than small interfaces with increasing hydrogen bonding energy to interface size. However, small interfaces have three-fold more electrostatics energy than large interfaces with increasing electrostatics to interface size. Thus, 8% of complexes having small interfaces with limited interface area and sub-dominant vdW energy are rich in electrostatics. It is interesting to observe that complexes having small interfaces are often associated with regulatory function. Hence, the observed structural features with known molecular function provide insights for the better understanding of PPI.Communicated by Ramaswamy H. Sarma.

The Subtle Trade-Off between Evolutionary and Energetic Constraints in Protein-Protein Interactions

Life machinery, although overwhelmingly complex, is rooted on a rather limited number of molecular processes. One of the most important is protein-protein interaction. Metabolic regulation, protein folding control, and cellular motility are examples of processes based on the fine-tuned interaction of several protein partners. The region on the protein surface devoted to the recognition of a specific partner is essential for the function of the protein and is, therefore, likely to be conserved during evolution. On the other hand, the physical chemistry of amino acids underlies the mechanism of interactions. Both evolutionary and energetic constraints can then be used to build scoring functions capable of recognizing interaction sites. Our working hypothesis is that residues within the interaction interface tend at the same time to be evolutionarily conserved (to preserve their function) and to provide little contribution to the internal stabilization of the structure of their cognate protein, to facilitate conformational adaptation to the partner. Here, we show that for some classes of protein partners (for example, those involved in signal transduction and in enzymes) evolutionary constraints play the key role in defining the interaction surface. In contrast, energetic constraints emerge as more important in protein partners involved in immune response, in inhibitor proteins, and in structural proteins. Our results indicate that a general-purpose scoring function for protein-protein interaction should not be agnostic of the biological function of the partners.