Dissecting subunit interfaces in homodimeric proteins

Assessing the functional impact of protein binding site definition

Many biomedical applications, such as classification of binding specificities or bioengineering, depend on the accurate definition of protein binding interfaces. Depending on the choice of method used, substantially different sets of residues can be classified as belonging to the interface of a protein. A typical approach used to verify these definitions is to mutate residues and measure the impact of these changes on binding. Besides the lack of exhaustive data, this approach also suffers from the fundamental problem that a mutation introduces an unknown amount of alteration into an interface, which potentially alters the binding characteristics of the interface. In this study we explore the impact of alternative binding site definitions on the ability of a protein to recognize its cognate ligand using a pharmacophore approach, which does not affect the interface. The study also shows that methods for protein binding interface predictions should perform above approximately F-score = 0.7 accuracy level to capture the biological function of a protein.

Interfacial residues in protein-protein complexes are in the eyes of the beholder

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner.

Lateral interactions govern self-assembly of the bacterial biofilm matrix protein BslA

The soil bacterium Bacillus subtilis is a model organism to investigate the formation of biofilms, the predominant form of microbial life. The secreted protein BslA self-assembles at the surface of the biofilm to give the B. subtilis biofilm its characteristic hydrophobicity. To understand the mechanism of BslA self-assembly at interfaces, here we built a molecular model based on the previous BslA crystal structure and the crystal structure of the BslA paralogue YweA that we determined. Our analysis revealed two conserved protein-protein interaction interfaces supporting BslA self-assembly into an infinite 2-dimensional lattice that fits previously determined transmission microscopy images. Molecular dynamics simulations and in vitro protein assays further support our model of BslA elastic film formation, while mutagenesis experiments highlight the importance of the identified interactions for biofilm structure. Based on this knowledge, YweA was engineered to form more stable elastic films and rescue biofilm structure in bslA deficient strains. These findings shed light on protein film assembly and will inform the development of BslA technologies which range from surface coatings to emulsions in fast-moving consumer goods.

Assessing the functional impact of protein binding site definition

Many biomedical applications, such as classification of binding specificities or bioengineering, depend on the accurate definition of protein binding interfaces. Depending on the choice of method used, substantially different sets of residues can be classified as belonging to the interface of a protein. A typical approach used to verify these definitions is to mutate residues and measure the impact of these changes on binding. Besides the lack of exhaustive data this approach generates, it also suffers from the fundamental problem that a mutation introduces an unknown amount of alteration into an interface, which potentially alters the binding characteristics of the interface. In this study we explore the impact of alternative binding site definitions on the ability of a protein to recognize its cognate ligand using a pharmacophore approach, which does not affect the interface. The study also provides guidance on the minimum expected accuracy of interface definition that is required to capture the biological function of a protein.

Cold spots are universal in protein–protein interactions

Proteins interact with each other through binding interfaces that differ greatly in size and physico‐chemical properties. Within the binding interface, a few residues called hot spots contribute the majority of the binding free energy and are hence irreplaceable. In contrast, cold spots are occupied by suboptimal amino acids, providing possibility for affinity enhancement through mutations. In this study, we identify cold spots due to cavities and unfavorable charge interactions in multiple protein–protein interactions (PPIs). For our cold spot analysis, we first use a small affinity database of PPIs with known structures and affinities and then expand our search to nearly 4000 homo‐ and heterodimers in the Protein Data Bank (PDB). We observe that cold spots due to cavities are present in nearly all PPIs unrelated to their binding affinity, while unfavorable charge interactions are relatively rare. We also find that most cold spots are located in the periphery of the binding interface, with high‐affinity complexes showing fewer centrally located colds spots than low‐affinity complexes. A larger number of cold spots is also found in non‐cognate interactions compared to their cognate counterparts. Furthermore, our analysis reveals that cold spots are more frequent in homo‐dimeric complexes compared to hetero‐complexes, likely due to symmetry constraints imposed on sequences of homodimers. Finally, we find that glycines, glutamates, and arginines are the most frequent amino acids appearing at cold spot positions. Our analysis emphasizes the importance of cold spot positions to protein evolution and facilitates protein engineering studies directed at enhancing binding affinity and specificity in a wide range of applications.

The Hidden Intricacies of Aquaporins: Remarkable Details in a Common Structural Scaffold

Evolution turned aquaporins (AQPs) into the most efficient facilitators of passive water flow through cell membranes at no expense of solute discrimination. In spite of a plethora of solved AQP structures, many structural details remain hidden. Here, by combining extensive sequence- and structural-based analysis of a unique set of 20 non-redundant high-resolution structures and molecular dynamics simulations of four representatives, key aspects of AQP stability, gating, selectivity, pore geometry, and oligomerization, with a potential impact on channel functionality, are identified. The general view of AQPs possessing a continuous open water pore is challenged and it is depicted that AQPs' selectivity is not exclusively shaped by pore-lining residues but also by the relative arrangement of transmembrane helices. Moreover, this analysis reveals that hydrophobic interactions constitute the main determinant of protein thermal stability. Finally, a numbering scheme of the conserved AQP scaffold is established, facilitating direct comparison of, for example, disease-causing mutations and prediction of potential structural consequences. Additionally, the results pave the way for the design of optimized AQP water channels to be utilized in biotechnological applications.

Over the rainbow: structural characterization of the chromoproteins gfasPurple, amilCP, spisPink and eforRed

Anthozoan chromoproteins are highly pigmented, diversely coloured and readily produced in recombinant expression systems. While they are a versatile and powerful building block in synthetic biology for applications such as biosensor development, they are not widely used in comparison to the related fluorescent proteins, partly due to a lack of structural characterization to aid protein engineering. Here, high-resolution X-ray crystal structures of four open-source chromoproteins, gfasPurple, amilCP, spisPink and eforRed, are presented. These proteins are dimers in solution, and mutation at the conserved dimer interface leads to loss of visible colour development in gfasPurple. The chromophores are trans and noncoplanar in gfasPurple, amilCP and spisPink, while that in eforRed is cis and noncoplanar, and also emits fluorescence. Like other characterized chromoproteins, gfasPurple, amilCP and eforRed contain an sp 2-hybridized N-acylimine in the peptide bond preceding the chromophore, while spisPink is unusual and demonstrates a true sp 3-hybridized trans-peptide bond at this position. It was found that point mutations at the chromophore-binding site in gfasPurple that substitute similar amino acids to those in amilCP and spisPink generate similar colours. These features and observations have implications for the utility of these chromoproteins in protein engineering and synthetic biology applications.

Characterization of Hydrophilic α-Helical Hot Spots on the Protein-Protein Interaction Interfaces for the Design of α-Helix Mimetics

The cooperativity index, K_c, was developed to examine the binding synergy between hot spots of the ligand-protein. For the first time, the convergence of the side-chain spatial arrangements of hydrophilic α-helical hot spots Thr, Tyr, Asp, Asn, Ser, Cys, and His in protein-protein interaction (PPI) complex structures was disclosed and quantified by developing novel clustering models. In-depth analyses revealed the driving force for the protein-protein binding conformation convergence of hydrophilic α-helical hot spots. This observation allows deriving pharmacophore models to design new mimetics for hydrophilic α-helical hot spots. A computational protocol was developed to search amino acid analogues and small-molecule mimetics for each hydrophilic α-helical hot spot. As a pilot study, diverse building blocks of commercially available nonstandard L-type α-amino acids and the phenyl ring-containing small-molecule fragments were obtained, which serve as a fragment collection to mimic hydrophilic α-helical hot spots for the improvement of binding affinity, selectivity, physicochemical properties, and synthesis accessibility of α-helix mimetics.

How sticky are our proteins? Quantifying hydrophobicity of the human proteome

Summary

Proteins tend to bury hydrophobic residues inside their core during the folding process to provide stability to the protein structure and to prevent aggregation. Nevertheless, proteins do expose some 'sticky' hydrophobic residues to the solvent. These residues can play an important functional role, e.g. in protein-protein and membrane interactions. Here, we first investigate how hydrophobic protein surfaces are by providing three measures for surface hydrophobicity: the total hydrophobic surface area, the relative hydrophobic surface area and-using our MolPatch method-the largest hydrophobic patch. Secondly, we analyze how difficult it is to predict these measures from sequence: by adapting solvent accessibility predictions from NetSurfP2.0, we obtain well-performing prediction methods for the THSA and RHSA, while predicting LHP is more challenging. Finally, we analyze implications of exposed hydrophobic surfaces: we show that hydrophobic proteins typically have low expression, suggesting cells avoid an overabundance of sticky proteins.

Availability and implementation

The data underlying this article are available in GitHub at https://github.com/ibivu/hydrophobic_patches.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

PDB-wide identification of physiological hetero-oligomeric assemblies based on conserved quaternary structure geometry

An accurate understanding of biomolecular mechanisms and diseases requires information on protein quaternary structure (QS). A critical challenge in inferring QS information from crystallography data is distinguishing biological interfaces from fortuitous crystal-packing contacts. Here, we employ QS conservation across homologs to infer the biological relevance of hetero-oligomers. We compare the structures and compositions of hetero-oligomers, which allow us to annotate 7,810 complexes as physiologically relevant, 1,060 as likely errors, and 1,432 with comparative information on subunit stoichiometry and composition. Excluding immunoglobulins, these annotations encompass over 51% of hetero-oligomers in the PDB. We curate a dataset of 577 hetero-oligomeric complexes to benchmark these annotations, which reveals an accuracy >94%. When homology information is not available, we compare QS across repositories (PDB, PISA, and EPPIC) to derive confidence estimates. This work provides high-quality annotations along with a large benchmark dataset of hetero-assemblies.

Marker residue types at the structural regions of transmembrane alpha-helical and beta-barrel interfaces

Membrane proteins play a variety of biological functions to the survival of organisms and functionalities of these proteins are often due to their homo- or hetero-complexation. Encoded by ~30% of the genome in most organisms, they represent the target of over half of nowadays drugs. Spanning the entirety of the cell membrane, transmembrane proteins are the most common type of membrane proteins and can be classified by secondary structures: alpha-helical and beta-barrel structures. Protein-protein interaction (PPI) have been widely studied for globular proteins and many computational tools are available for predicting PPI sites and construct models of complexes. Here, the structural regions of a non-redundant set of 232 alpha-helical and 37 beta-barrel transmembrane complexes and their interfaces are analyzed. Using the residue composition, frequency and propensity, this study brings the light on the marker residue types located at the structural regions of alpha-helical and beta-barrel transmembrane homomeric protein complexes and of their interfaces. This study also shows the necessity to relate the frequency to the composition into a ratio for immediately figuring out residue types presenting high frequencies at the interface and/or at one of its structural regions despite being a minor contributor compared to other residue types to that location's residue composition.

Membrane transporter dimerization driven by differential lipid solvation energetics of dissociated and associated states

Over two-thirds of integral membrane proteins of known structure assemble into oligomers. Yet, the forces that drive the association of these proteins remain to be delineated, as the lipid bilayer is a solvent environment that is both structurally and chemically complex. In this study, we reveal how the lipid solvent defines the dimerization equilibrium of the CLC-ec1 Cl^-/H⁺ antiporter. Integrating experimental and computational approaches, we show that monomers associate to avoid a thinned-membrane defect formed by hydrophobic mismatch at their exposed dimerization interfaces. In this defect, lipids are strongly tilted and less densely packed than in the bulk, with a larger degree of entanglement between opposing leaflets and greater water penetration into the bilayer interior. Dimerization restores the membrane to a near-native state and therefore, appears to be driven by the larger free-energy cost of lipid solvation of the dissociated protomers. Supporting this theory, we demonstrate that addition of short-chain lipids strongly shifts the dimerization equilibrium toward the monomeric state, and show that the cause of this effect is that these lipids preferentially solvate the defect. Importantly, we show that this shift requires only minimal quantities of short-chain lipids, with no measurable impact on either the macroscopic physical state of the membrane or the protein's biological function. Based on these observations, we posit that free-energy differentials for local lipid solvation define membrane-protein association equilibria. With this, we argue that preferential lipid solvation is a plausible cellular mechanism for lipid regulation of oligomerization processes, as it can occur at low concentrations and does not require global changes in membrane properties.

A cell’s outer membrane is made of molecules called lipids, which band together to form a flexible thin film, just two molecules thick. This membrane is dotted with proteins that transport materials in to and out of cells. Most of these membrane proteins join with other proteins to form structures known as oligomers. Except, how membrane-bound proteins assemble into oligomers – the physical forces driving these molecules to take shape – remains unclear.

This is partly because the structural, physical and chemical properties of fat-like lipid membranes are radically different to the cell’s watery interior. Consequently, the conditions under which membrane oligomers form are distinct from those surrounding proteins inside cells. Membrane proteins are also more difficult to study and characterize than water-soluble proteins inside the cell, and yet many therapeutic drugs such as antibiotics specifically target membrane proteins. Overall, our understanding of how the unique properties of lipid membranes affect the formation of protein structures embedded within, is lacking and warrants further investigation.

Now, Chadda, Bernhardt et al. focused on one membrane protein, known as CLC, which tends to exist in pairs – or dimers. To understand why these proteins form dimers (a process called dimerization) Chadda, Bernhardt et al. first used computer simulations, and then validated the findings in experimental tests.

These complementary approaches demonstrated that the main reason CLC proteins ‘dimerize’ lies in their interaction with the lipid membrane, and not the attraction of one protein to its partner. When CLC proteins are on their own, they deform the surrounding membrane and create structural defects that put the membrane under strain. But when two CLC proteins join as a dimer, this membrane strain disappears – making dimerization the more stable and energetically favorable option.

Chadda, Bernhardt et al. also showed that with the addition of a few certain lipids, specifically smaller lipids, cell membranes become more tolerant of protein-induced structural changes. This might explain how cells could use various lipids to fine-tune the activity of membrane proteins by controlling how oligomers form. However, the theory needs to be examined further. Altogether, this work has provided fundamental insights into the physical forces shaping membrane-bound proteins, relevant to researchers studying cell biology and pharmacology alike.

Classification and prediction of protein-protein interaction interface using machine learning algorithm

Structural insight of the protein-protein interaction (PPI) interface can provide knowledge about the kinetics, thermodynamics and molecular functions of the complex while elucidating its role in diseases and further enabling it as a potential therapeutic target. However, owing to experimental lag in solving protein-protein complex structures, three-dimensional (3D) knowledge of the PPI interfaces can be gained via computational approaches like molecular docking and post-docking analyses. Despite development of numerous docking tools and techniques, success in identification of native like interfaces based on docking score functions is limited. Hence, we employed an in-depth investigation of the structural features of the interface that might successfully delineate native complexes from non-native ones. We identify interface properties, which show statistically significant difference between native and non-native interfaces belonging to homo and hetero, protein-protein complexes. Utilizing these properties, a support vector machine (SVM) based classification scheme has been implemented to differentiate native and non-native like complexes generated using docking decoys. Benchmarking and comparative analyses suggest very good performance of our SVM classifiers. Further, protein interactions, which are proven via experimental findings but not resolved structurally, were subjected to this approach where 3D-models of the complexes were generated and most likely interfaces were predicted. A web server called Protein Complex Prediction by Interface Properties (PCPIP) is developed to predict whether interface of a given protein-protein dimer complex resembles known protein interfaces. The server is freely available at http://www.hpppi.iicb.res.in/pcpip/ .

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.

Molecular basis for protein-protein interactions

This minireview provides an overview on the current knowledge of protein-protein interactions, common characterisation methods to characterise them, and their role in protein complex formation with some examples. A deep understanding of protein-protein interactions and their molecular interactions is important for a number of applications, including drug design. Protein-protein interactions and their discovery are thus an interesting avenue for understanding how protein complexes, which make up the majority of proteins, work.

Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other's substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.

Conditional Disorder in Small Heat-shock Proteins

Small heat-shock proteins (sHSPs) are molecular chaperones that respond to cellular stresses to combat protein aggregation. HSP27 is a critical human sHSP that forms large, dynamic oligomers whose quaternary structures and chaperone activities depend on environmental factors. Upon exposure to cellular stresses, such as heat shock or acidosis, HSP27 oligomers can dissociate into dimers and monomers, which leads to significantly enhanced chaperone activity. The structured core of the protein, the α-crystallin domain (ACD), forms dimers and can prevent the aggregation of substrate proteins to a similar degree as the full-length protein. When the ACD dimer dissociates into monomers, it partially unfolds and exhibits enhanced activity. Here, we used solution-state NMR spectroscopy to characterize the structure and dynamics of the HSP27 ACD monomer. Web show that the monomer is stabilized at low pH and that its backbone chemical shifts, 15N relaxation rates, and 1H-15N residual dipolar couplings suggest structural changes and rapid motions in the region responsible for dimerization. By analyzing the solvent accessible and buried surface areas of sHSP structures in the context of a database of dimers that are known to dissociate into disordered monomers, we predict that ACD dimers from sHSPs across all kingdoms of life may partially unfold upon dissociation. We propose a general model in which conditional disorder-the partial unfolding of ACDs upon monomerization-is a common mechanism for sHSP activity.

Analyses on clustering of the conserved residues at protein-RNA interfaces and its application in binding site identification

Background

The maintenance of protein structural stability requires the cooperativity among spatially neighboring residues. Previous studies have shown that conserved residues tend to occur clustered together within enzyme active sites and protein-protein/DNA interfaces. It is possible that conserved residues form one or more local clusters in protein tertiary structures as it can facilitate the formation of functional motifs. In this work, we systematically investigate the spatial distributions of conserved residues as well as hot spot ones within protein-RNA interfaces.

Results

The analysis of 191 polypeptide chains from 160 complexes shows the polypeptides interacting with tRNAs evolve relatively rapidly. A statistical analysis of residues in different regions shows that the interface residues are often more conserved, while the most conserved ones are those occurring at protein interiors which maintain the stability of folded polypeptide chains. Additionally, we found that 77.8% of the interfaces have the conserved residues clustered within the entire interface regions. Appling the clustering characteristics to the identification of the real interface, there are 31.1% of cases where the real interfaces are ranked in top 10% of 1000 randomly generated surface patches. In the conserved clusters, the preferred residues are the hydrophobic (Leu, Ile, Met), aromatic (Tyr, Phe, Trp) and interestingly only one positively charged Arg residues. For the hot spot residues, 51.5% of them are situated in the conserved residue clusters, and they are largely consistent with the preferred residue types in the conserved clusters.

Conclusions

The protein-RNA interface residues are often more conserved than non-interface surface ones. The conserved interface residues occur more spatially clustered relative to the entire interface residues. The high consistence of hot spot residue types and the preferred residue types in the conserved clusters has important implications for the experimental alanine scanning mutagenesis study. This work deepens the understanding of the residual organization at protein-RNA interface and is of potential applications in the identification of binding site and hot spot residues.

Biological vs. Crystallographic Protein Interfaces: An Overview of Computational Approaches for Their Classification

Complexes between proteins are at the basis of almost every process in cells. Their study, from a structural perspective, has a pivotal role in understanding biological functions and, importantly, in drug development. X-ray crystallography represents the broadest source for the experimental structural characterization of protein-protein complexes. Correctly identifying the biologically relevant interface from the crystallographic ones is, however, not trivial and can be prone to errors. Over the past two decades, computational methodologies have been developed to study the differences of those interfaces and automatically classify them as biological or crystallographic. Overall, protein-protein interfaces show differences in terms of composition, energetics and evolutionary conservation between biological and crystallographic ones. Based on those observations, a number of computational methods have been developed for this classification problem, which can be grouped into three main categories: Energy-, empirical knowledge- and machine learning-based approaches. In this review, we give a comprehensive overview of the training datasets and methods so far implemented, providing useful links and a brief description of each method.

Multiple protein-DNA interfaces unravelled by evolutionary information, physico-chemical and geometrical properties

Interactions between proteins and nucleic acids are at the heart of many essential biological processes. Despite increasing structural information about how these interactions may take place, our understanding of the usage made of protein surfaces by nucleic acids is still very limited. This is in part due to the inherent complexity associated to protein surface deformability and evolution. In this work, we present a method that contributes to decipher such complexity by predicting protein-DNA interfaces and characterizing their properties. It relies on three biologically and physically meaningful descriptors, namely evolutionary conservation, physico-chemical properties and surface geometry. We carefully assessed its performance on several hundreds of protein structures and compared it to several machine-learning state-of-the-art methods. Our approach achieves a higher sensitivity compared to the other methods, with a similar precision. Importantly, we show that it is able to unravel ‘hidden’ binding sites by applying it to unbound protein structures and to proteins binding to DNA via multiple sites and in different conformations. It is also applicable to the detection of RNA-binding sites, without significant loss of performance. This confirms that DNA and RNA-binding sites share similar properties. Our method is implemented as a fully automated tool, JETDNA2, freely accessible at: http://www.lcqb.upmc.fr/JET2DNA. We also provide a new dataset of 187 protein-DNA complex structures, along with a subset of 82 associated unbound structures. The set represents the largest body of high-resolution crystallographic structures of protein-DNA complexes, use biological protein assemblies as DNA-binding units, and covers all major types of protein-DNA interactions. It is available at: http://www.lcqb.upmc.fr/PDNAbenchmarks.

Protein-DNA interactions are essential to living organisms and their impairment is associated to many diseases. For these reasons, they have become increasingly important therapeutic targets. Experimental structure determination has revealed different binding motifs and modes, associated to different functions. Yet, the available structural data gives us only a glimpse of the multiplicity and complexity of protein surface usage by DNA. In this work, we use a three-layer model to describe and predict DNA-binding sites at protein surfaces. Given a protein, we consider the way its residues are conserved through evolution, their physico-chemical properties and geometrical shapes to decrypt its surface. We are able to detect a large portion of interacting residues with good precision, even when they are ‘hidden’ by conformational changes. We highlight cases where one protein binds DNA via distinct regions to perform different functions. We are able to uncover the alternative binding sites and relate their properties with their specific roles. Our work can help guiding mutagenesis experiments and the development of new drugs specifically targeting one site while limiting possible side effects.

Computational Analysis of Non-covalent Interactions in Phycocyanin Subunit Interfaces

Protein-protein interactions are an important phenomenon in biological processes and functions. We used the manually curated non-redundant dataset of 118 phycocyanin interfaces to gain additional insight into this phenomenon using a robust inter-atomic non-covalent interaction analyzing tool PPCheck. Our observations indicate that there is a relatively high composition of hydrophobic residues at the interfaces. Most of the interface residues are clustered at the middle of the range which we call "standard-size" interfaces. Furthermore, the multiple interaction patterns founded in the present study indicate that more than half of the residues involved in these interactions participate in multiple and water-bridged hydrogen bonds. Thus, hydrogen bonds contribute maximally towards the stability of protein-protein complexes. The analysis shows that hydrogen bond energies contribute to about 88 % to the total energy and it also increases with interface size. Van der Waals (vdW) energy contributes to 9.3 %±1.7 % on average in these complexes. Moreover, there is about 1.9 %±1.5 % contribution by electrostatic energy. Nevertheless, the role by vdW and electrostatic energy could not be ignored in interface binding. Results show that the total binding energy is more for large phycocyanin interfaces. The normalized energy per residue was less than -16 kJ mol^-1 , while most of them have energy in the range from -6 to -14 kJ mol^-1 . The non-covalent interacting residues in these proteins were found to be highly conserved. Obtained results might contribute to the understanding of structural stability of this class of evolutionary essential proteins with increased practical application and future designs of novel protein-bioactive compound interactions.

Accurate Classification of Biological and non-Biological Interfaces in Protein Crystal Structures using Subtle Covariation Signals

Proteins often work as oligomers or multimers in vivo. Therefore, elucidating their oligomeric or multimeric form (quaternary structure) is crucially important to ascertain their function. X-ray crystal structures of numerous proteins have been accumulated, providing information related to their biological units. Extracting information of biological units from protein crystal structures represents a meaningful task for modern biology. Nevertheless, although many methods have been proposed for identifying biological units appearing in protein crystal structures, it is difficult to distinguish biological protein-protein interfaces from crystallographic ones. Therefore, our simple but highly accurate classifier was developed to infer biological units in protein crystal structures using large amounts of protein sequence information and a modern contact prediction method to exploit covariation signals (CSs) in proteins. We demonstrate that our proposed method is promising even for weak signals of biological interfaces. We also discuss the relation between classification accuracy and conservation of biological units, and illustrate how the selection of sequences included in multiple sequence alignments as sources for obtaining CSs affects the results. With increased amounts of sequence data, the proposed method is expected to become increasingly useful.

Dissecting macromolecular recognition sites in ribosome: implication to its self-assembly

Interactions between macromolecules play a crucial role in ribosome assembly that follows a highly coordinated process involving RNA folding and binding of ribosomal proteins (r-proteins). Although extensive studies have been carried out to understand macromolecular interactions in ribosomes, most of them are confined to either large or small ribosomal-subunit of few species. A comparative analysis of macromolecular interactions across different domains is still missing. We have analyzed the structural and physicochemical properties of protein-protein (PP), protein-RNA (PR) and RNA-RNA (RR) interfaces in small and large subunits of ribosomes, as well as in between the two subunits. Additionally, we have also developed Random Forest (RF) classifier to catalog the r-proteins. We find significant differences as well as similarities in macromolecular recognition sites between ribosomal assemblies of prokaryotes and eukaryotes. PR interfaces are substantially larger and have more ionic interactions than PP and RR interfaces in both prokaryotes and eukaryotes. PP, PR and RR interfaces in eukaryotes are well packed compared to those in prokaryotes. However, the packing density between the large and the small subunit interfaces in the entire assembly is strikingly low in both prokaryotes and eukaryotes, indicating the periodic association and dissociation of the two subunits during the translation. The structural and physicochemical properties of PR interfaces are used to predict the r-proteins in the assembly pathway into early, intermediate and late binders using RF classifier with an accuracy of 80%. The results provide new insights into the classification of r-proteins in the assembly pathway.

Kinetics and Mechanism of Mammalian Mitochondrial Ribosome Assembly

Mammalian mtDNA encodes only 13 proteins, all essential components of respiratory complexes, synthesized by mitochondrial ribosomes. Mitoribosomes contain greatly truncated RNAs transcribed from mtDNA, including a structural tRNA in place of 5S RNA as a scaffold for binding 82 nucleus-encoded proteins, mitoribosomal proteins (MRPs). Cryoelectron microscopy (cryo-EM) studies have determined the structure of the mitoribosome, but its mechanism of assembly is unknown. Our SILAC pulse-labeling experiments determine the rates of mitochondrial import of MRPs and their assembly into intact mitoribosomes, providing a basis for distinguishing MRPs that bind at early and late stages in mitoribosome assembly to generate a working model for mitoribosome assembly. Mitoribosome assembly is a slow process initiated at the mtDNA nucleoid driven by excess synthesis of individual MRPs. MRPs that are tightly associated in the structure frequently join the complex in a coordinated manner. Clinically significant MRP mutations reported to date affect proteins that bind early on during assembly.

Principles and characteristics of biological assemblies in experimentally determined protein structures

More than half of all structures in the PDB are assemblies of two or more proteins, including both homooligomers and heterooligomers. Structural information on these assemblies comes from X-ray crystallography, NMR, and cryo-EM spectroscopy. The correct assembly in an X-ray structure is often ambiguous, and computational methods have been developed to identify the most likely biologically relevant assembly based on physical properties of assemblies and sequence conservation in interfaces. Taking advantage of the large number of structures now available, some of the most recent methods have relied on similarity of interfaces and assemblies across structures of homologous proteins.

Distinguishing crystallographic from biological interfaces in protein complexes: role of intermolecular contacts and energetics for classification

BACKGROUND

Study of macromolecular assemblies is fundamental to understand functions in cells. X-ray crystallography is the most common technique to solve their 3D structure at atomic resolution. In a crystal, however, both biologically-relevant interfaces and non-specific interfaces resulting from crystallographic packing are observed. Due to the complexity of the biological assemblies currently tackled, classifying those interfaces, i.e. distinguishing biological from crystal lattice interfaces, is not trivial and often prone to errors. In this context, analyzing the physico-chemical characteristics of biological/crystal interfaces can help researchers identify possible features that distinguish them and gain a better understanding of the systems.

RESULTS

In this work, we are providing new insights into the differences between biological and crystallographic complexes by focusing on "pair-properties" of interfaces that have not yet been fully investigated. We investigated properties such intermolecular residue-residue contacts (already successfully applied to the prediction of binding affinities) and interaction energies (electrostatic, Van der Waals and desolvation). By using the XtalMany and BioMany interface datasets, we show that interfacial residue contacts, classified as a function of their physico-chemical properties, can distinguish between biological and crystallographic interfaces. The energetic terms show, on average, higher values for crystal interfaces, reflecting a less stable interface due to crystal packing compared to biological interfaces. By using a variety of machine learning approaches, we trained a new interface classification predictor based on contacts and interaction energetic features. Our predictor reaches an accuracy in classifying biological vs crystal interfaces of 0.92, compared to 0.88 for EPPIC (one of the main state-of-the-art classifiers reporting same performance as PISA).

CONCLUSION

In this work we have gained insights into the nature of intermolecular contacts and energetics terms distinguishing biological from crystallographic interfaces. Our findings might have a broader applicability in structural biology, for example for the identification of near native poses in docking. We implemented our classification approach into an easy-to-use and fast software, freely available to the scientific community from http://github.com/haddocking/interface-classifier .

Integrating co-evolutionary signals and other properties of residue pairs to distinguish biological interfaces from crystal contacts

It remains challenging to accurately discriminate between biological and crystal interfaces. Most existing analyses and algorithms focused on the features derived from a single side of the interface. However, less attention has been paid to the properties of residue pairs across protein interfaces. To address this problem, we defined a novel co-evolutionary feature for homodimers through integrating direct coupling analysis and image processing techniques. The residue pairs across biological homodimeric interfaces were significantly enriched in co-evolving residues compared to those across crystal contacts, resulting in a promising classification accuracy with area under the curves (AUCs) of >0.85. Considering the availability of co-evolutionary feature, we also designed other residue pair based features that were useful for both homodimers and heterodimers. The most informative residue pairs were identified to reflect the interaction preferences across protein interfaces. Regarding the other extant properties, we designed the new descriptors at the interface residue level as well as at the pairwise contact level. Extensive validation showed that these single properties can be used to identify biological interfaces with AUCs ranging from 0.60 to 0.88. By integrating co-evolutionary feature with other residue pair based properties, our final prediction model output excellent performance with AUCs of >0.91 on different datasets. Compared to existing methods, our algorithm not only yielded better or comparable results but also provided complementary information. An easy-to-use web server is freely accessible at http://liulab.hzau.edu.cn/RPAIAnalyst.

Evolution of In Silico Strategies for Protein-Protein Interaction Drug Discovery

The advent of advanced molecular modeling software, big data analytics, and high-speed processing units has led to the exponential evolution of modern drug discovery and better insights into complex biological processes and disease networks. This has progressively steered current research interests to understanding protein-protein interaction (PPI) systems that are related to a number of relevant diseases, such as cancer, neurological illnesses, metabolic disorders, etc. However, targeting PPIs are challenging due to their "undruggable" binding interfaces. In this review, we focus on the current obstacles that impede PPI drug discovery, and how recent discoveries and advances in in silico approaches can alleviate these barriers to expedite the search for potential leads, as shown in several exemplary studies. We will also discuss about currently available information on PPI compounds and systems, along with their usefulness in molecular modeling. Finally, we conclude by presenting the limits of in silico application in drug discovery and offer a perspective in the field of computer-aided PPI drug discovery.

Investigation of protein quaternary structure via stoichiometry and symmetry information

The Protein Data Bank (PDB) is the single worldwide archive of experimentally-determined three-dimensional (3D) structures of proteins and nucleic acids. As of January 2017, the PDB housed more than 125,000 structures and was growing by more than 11,000 structures annually. Since the 3D structure of a protein is vital to understand the mechanisms of biological processes, diseases, and drug design, correct oligomeric assembly information is of critical importance. Unfortunately, the biologically relevant oligomeric form of a 3D structure is not directly obtainable by X-ray crystallography, whilst in solution methods (NMR or single particle EM) it is known from the experiment. Instead, this information may be provided by the PDB Depositor as metadata coming from additional experiments, be inferred by sequence-sequence comparisons with similar proteins of known oligomeric state, or predicted using software, such as PISA (Proteins, Interfaces, Structures and Assemblies) or EPPIC (Evolutionary Protein Protein Interface Classifier). Despite significant efforts by professional PDB Biocurators during data deposition, there remain a number of structures in the archive with incorrect quaternary structure descriptions (or annotations). Further investigation is, therefore, needed to evaluate the correctness of quaternary structure annotations. In this study, we aim to identify the most probable oligomeric states for proteins represented in the PDB. Our approach evaluated the performance of four independent prediction methods, including text mining of primary publications, inference from homologous protein structures, and two computational methods (PISA and EPPIC). Aggregating predictions to give consensus results outperformed all four of the independent prediction methods, yielding 83% correct, 9% wrong, and 8% inconclusive predictions, when tested with a well-curated benchmark dataset. We have developed a freely-available web-based tool to make this approach accessible to researchers and PDB Biocurators (http://quatstruct.rcsb.org/).

Comparing side chain packing in soluble proteins, protein-protein interfaces, and transmembrane proteins

We compare side chain prediction and packing of core and non-core regions of soluble proteins, protein-protein interfaces, and transmembrane proteins. We first identified or created comparable databases of high-resolution crystal structures of these 3 protein classes. We show that the solvent-inaccessible cores of the 3 classes of proteins are equally densely packed. As a result, the side chains of core residues at protein-protein interfaces and in the membrane-exposed regions of transmembrane proteins can be predicted by the hard-sphere plus stereochemical constraint model with the same high prediction accuracies (>90%) as core residues in soluble proteins. We also find that for all 3 classes of proteins, as one moves away from the solvent-inaccessible core, the packing fraction decreases as the solvent accessibility increases. However, the side chain predictability remains high (80% within 30°) up to a relative solvent accessibility, rSASA≲0.3, for all 3 protein classes. Our results show that ≈40% of the interface regions in protein complexes are "core", that is, densely packed with side chain conformations that can be accurately predicted using the hard-sphere model. We propose packing fraction as a metric that can be used to distinguish real protein-protein interactions from designed, non-binding, decoys. Our results also show that cores of membrane proteins are the same as cores of soluble proteins. Thus, the computational methods we are developing for the analysis of the effect of hydrophobic core mutations in soluble proteins will be equally applicable to analyses of mutations in membrane proteins.

Crystal structures of eukaryote glycosyltransferases reveal biologically relevant enzyme homooligomers

Glycosyltransferases (GTases) transfer sugar moieties to proteins, lipids or existing glycan or polysaccharide molecules. GTases form an important group of enzymes in the Golgi, where the synthesis and modification of glycoproteins and glycolipids take place. Golgi GTases are almost invariably type II integral membrane proteins, with the C-terminal globular catalytic domain residing in the Golgi lumen. The enzymes themselves are divided into 103 families based on their sequence homology. There is an abundance of published crystal structures of GTase catalytic domains deposited in the Protein Data Bank (PDB). All of these represent either of the two main characteristic structural folds, GT-A or GT-B, or present a variation thereof. Since GTases can function as homomeric or heteromeric complexes in vivo, we have summarized the structural features of the dimerization interfaces in crystal structures of GTases, as well as considered the biochemical data available for these enzymes. For this review, we have considered all 898 GTase crystal structures in the Protein Data Bank and highlight the dimer formation characteristics of various GTases based on 24 selected structures.

Exploring the potential of 3D Zernike descriptors and SVM for protein-protein interface prediction

Background

The correct determination of protein–protein interaction interfaces is important for understanding disease mechanisms and for rational drug design. To date, several computational methods for the prediction of protein interfaces have been developed, but the interface prediction problem is still not fully understood. Experimental evidence suggests that the location of binding sites is imprinted in the protein structure, but there are major differences among the interfaces of the various protein types: the characterising properties can vary a lot depending on the interaction type and function. The selection of an optimal set of features characterising the protein interface and the development of an effective method to represent and capture the complex protein recognition patterns are of paramount importance for this task.

Results

In this work we investigate the potential of a novel local surface descriptor based on 3D Zernike moments for the interface prediction task. Descriptors invariant to roto-translations are extracted from circular patches of the protein surface enriched with physico-chemical properties from the HQI8 amino acid index set, and are used as samples for a binary classification problem. Support Vector Machines are used as a classifier to distinguish interface local surface patches from non-interface ones. The proposed method was validated on 16 classes of proteins extracted from the Protein–Protein Docking Benchmark 5.0 and compared to other state-of-the-art protein interface predictors (SPPIDER, PrISE and NPS-HomPPI).

Conclusions

The 3D Zernike descriptors are able to capture the similarity among patterns of physico-chemical and biochemical properties mapped on the protein surface arising from the various spatial arrangements of the underlying residues, and their usage can be easily extended to other sets of amino acid properties. The results suggest that the choice of a proper set of features characterising the protein interface is crucial for the interface prediction task, and that optimality strongly depends on the class of proteins whose interface we want to characterise. We postulate that different protein classes should be treated separately and that it is necessary to identify an optimal set of features for each protein class.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2043-3) contains supplementary material, which is available to authorized users.

Inferring and Using Protein Quaternary Structure Information from Crystallographic Data

A precise knowledge of the quaternary structure of proteins is essential to illuminate both their function and their evolution. The major part of our knowledge on quaternary structure is inferred from X-ray crystallography data, but this inference process is hard and error-prone. The difficulty lies in discriminating fortuitous protein contacts, which make up the lattice of protein crystals, from biological protein contacts that exist in the native cellular environment. Here, we review methods devised to discriminate between both types of contacts and describe resources for downloading protein quaternary structure information and identifying high-confidence quaternary structures. The use of high-confidence datasets of quaternary structures will be critical for the analysis of structural, functional, and evolutionary properties of proteins.

In silico structure-based approaches to discover protein-protein interaction-targeting drugs

A core concept behind modern drug discovery is finding a small molecule that modulates a function of a target protein. This concept has been successfully applied since the mid-1970s. However, the efficiency of drug discovery is decreasing because the druggable target space in the human proteome is limited. Recently, protein-protein interaction (PPI) has been identified asan emerging target space for drug discovery. PPI plays a pivotal role in biological pathways including diseases. Current human interactome research suggests that the number of PPIs is between 130,000 and 650,000, and only a small number of them have been targeted as drug targets. For traditional drug targets, in silico structure-based methods have been successful in many cases. However, their performance suffers on PPI interfaces because PPI interfaces are different in five major aspects: From a geometric standpoint, they have relatively large interface regions, flat geometry, and the interface surface shape tends to fluctuate upon binding. Also, their interactions are dominated by hydrophobic atoms, which is different from traditional binding-pocket-targeted drugs. Finally, PPI targets usually lack natural molecules that bind to the target PPI interface. Here, we first summarize characteristics of PPI interfaces and their known binders. Then, we will review existing in silico structure-based approaches for discovering small molecules that bind to PPI interfaces.

Protein-protein interactions leave evolutionary footprints: High molecular coevolution at the core of interfaces

Protein-protein interactions are essential to all aspects of life. Specific interactions result from evolutionary pressure at the interacting interfaces of partner proteins. However, evolutionary pressure is not homogeneous within the interface: for instance, each residue does not contribute equally to the binding energy of the complex. To understand functional differences between residues within the interface, we analyzed their properties in the core and rim regions. Here, we characterized protein interfaces with two evolutionary measures, conservation and coevolution, using a comprehensive dataset of 896 protein complexes. These scores can detect different selection pressures at a given position in a multiple sequence alignment. We also analyzed how the number of interactions in which a residue is involved influences those evolutionary signals. We found that the coevolutionary signal is higher in the interface core than in the interface rim region. Additionally, the difference in coevolution between core and rim regions is comparable to the known difference in conservation between those regions. Considering proteins with multiple interactions, we found that conservation and coevolution increase with the number of different interfaces in which a residue is involved, suggesting that more constraints (i.e., a residue that must satisfy a greater number of interactions) allow fewer sequence changes at those positions, resulting in higher conservation and coevolution values. These findings shed light on the evolution of protein interfaces and provide information useful for identifying protein interfaces and predicting protein-protein interactions.

Dimerization of the adaptor Gads facilitates antigen receptor signaling by promoting the cooperative binding of Gads to the adaptor LAT

The accurate assembly of signalosomes centered on the adaptor protein LAT (linker of activated T cells) is required for antigen receptor signaling in T cells and mast cells. During signalosome assembly, members of the growth factor receptor-bound protein 2 (Grb2) family of cytosolic adaptor proteins bind cooperatively to LAT through interactions with its phosphorylated tyrosine (pTyr) residues. We demonstrated the Src homology 2 (SH2) domain-mediated dimerization of the Grb2 family member, Grb2-related adaptor downstream of Shc (Gads). Gads dimerization was mediated by an SH2 domain interface, which is distinct from the pTyr binding pocket and which promoted cooperative, preferential binding of paired Gads to LAT. This SH2 domain-intrinsic mechanism of cooperativity, which we quantified by mathematical modeling, enabled Gads to discriminate between dually and singly phosphorylated LAT molecules. Mutational inactivation of the dimerization interface reduced cooperativity and abrogated Gads signaling in T cells and mast cells. The dimerization-dependent, cooperative binding of Gads to LAT may increase antigen receptor sensitivity by reducing signalosome formation at incompletely phosphorylated LAT molecules, thereby prioritizing the formation of complete signalosomes.

Functions of tryptophan residues in EWGWS insert of Plasmodium falciparum enolase

Plasmodium falciparum enolase (Pfeno) is a dimeric enzyme with multiple moonlighting functions. This enzyme is thus a potential target for anti-malarial treatments. A unique feature of Pfeno is the presence of a pentapeptide insert ¹⁰⁴ EWGWS ¹⁰⁸. The functional role of tryptophan residues in this insert was investigated using site-directed mutagenesis. Replacement of these two Trp residues with alanines (or lysines) resulted in a near complete loss of enolase activity and dissociation of the normal dimeric form into monomers. Molecular modeling indicated that ³⁴⁰R forms π-cation bonds with the aromatic rings of ¹⁰⁵W and ⁴⁶Y. Mutation induced changes in the interactions among these three residues were presumably relayed to the inter-subunit interface via a coil formed by ⁴⁶Y : ⁵⁹Y, resulting in the disruption of a salt bridge between ¹¹R : ⁴²⁵E and a π-cation interaction between ¹¹R : ⁵⁹Y. This led to a drop of ~ 4 kcal·mole^-1 in the inter-subunit docking energy in the mutant, causing a ~ 10³ fold decrease in affinity. Partial restoration of the inter-subunit interactions led to reformation of dimers and also restored a significant fraction of the lost enzyme activity. These results suggested that the perturbations in the conformation of the surface loop containing the insert sequence were relayed to the interface region, causing dimer dissociation that, in turn, disrupted the enzyme's active site. Since Plasmodium enolase is a moonlighting protein with multiple parasite-specific functions, it is likely that these functions may map on to the highly conserved unique insert region of this protein.

ENZYMES

Enolase(EC4.2.1.11).

Improving protein-protein interaction prediction using evolutionary information from low-quality MSAs

Evolutionary information stored in multiple sequence alignments (MSAs) has been used to identify the interaction interface of protein complexes, by measuring either co-conservation or co-mutation of amino acid residues across the interface. Recently, maximum entropy related correlated mutation measures (CMMs) such as direct information, decoupling direct from indirect interactions, have been developed to identify residue pairs interacting across the protein complex interface. These studies have focussed on carefully selected protein complexes with large, good-quality MSAs. In this work, we study protein complexes with a more typical MSA consisting of fewer than 400 sequences, using a set of 79 intramolecular protein complexes. Using a maximum entropy based CMM at the residue level, we develop an interface level CMM score to be used in re-ranking docking decoys. We demonstrate that our interface level CMM score compares favourably to the complementarity trace score, an evolutionary information-based score measuring co-conservation, when combined with the number of interface residues, a knowledge-based potential and the variability score of individual amino acid sites. We also demonstrate, that, since co-mutation and co-complementarity in the MSA contain orthogonal information, the best prediction performance using evolutionary information can be achieved by combining the co-mutation information of the CMM with co-conservation information of a complementarity trace score, predicting a near-native structure as the top prediction for 41% of the dataset. The method presented is not restricted to small MSAs, and will likely improve interface prediction also for complexes with large and good-quality MSAs.

Specificity and stability of transient protein-protein interactions

Remarkable features that are achieved in a protein-protein complex to precise levels are stability and specificity. Deviation from the normal levels of specificity and stability, which is often caused by mutations, could result in disease conditions. Chemical nature, 3-D arrangement and dynamics of interface residues code for both specificity and stability. This article reviews roles of interfacial residues in transient protein-protein complexes. It is proposed that aside from hotspot residues conferring stability to the complex, a small set of 'rigid' residues at the interface that maintain conformation between complexed and uncomplexed forms, play a major role in conferring specificity. Exceptionally, 'super hotspot' residues, which confer both stability and specificity, are attractive sites for interaction with small molecule inhibitors.

A novel index of protein-protein interface propensity improves interface residue recognition

Background

Protein-protein interface holds important information of protein-protein interactions which play key roles in most biological processes. In the past few years, a lot of efforts have been made to improve interface residue recognition by characterizing protein-protein interfaces and extracting relevant features. However, most previous studies were carried out in a qualitative level, and there are also some inconsistencies between them.

Results

In the present work, to improve interface residue recognition, we built a novel quantitative residue protein-protein interface propensity index (QIPI) and gained a comprehensive picture of protein-protein interface through analyzing protein-protein interfaces on our comprehensive protein-protein interfaces dataset (Astral2.05-40-4506). Furthermore, in order to assess the effect of QIPI in improving the protein-protein interface prediction, we developed an interface residue recognition method SPR (Single domain based Patch Recognition) based on the QIPI. The evaluation results proved that our novel QIPI is able to improve the interface residue recognition.

Conclusions

Through a comprehensive quantitative analysis of protein-protein interface, we constructed a novel quantitative protein-protein interface propensity index (QIPI), which could be easily applied to improve the interface residue recognition and helpful in understanding the protein-protein interface.

Availability

QIPI and SPR are available to non-commercial users at our website: http://www.scbit.org/QIPI/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0351-7) contains supplementary material, which is available to authorized users.

Computational design of self-assembling cyclic protein homo-oligomers

Self-assembling cyclic protein homo-oligomers play important roles in biology and the ability to generate custom homo-oligomeric structures could enable new approaches to probe biological function. Here we report a general approach to design cyclic homo-oligomers that employs a new residue pair transform method for assessing the design ability of a protein-protein interface. This method is sufficiently rapid to enable systematic enumeration of cyclically docked arrangements of a monomer followed by sequence design of the newly formed interfaces. We use this method to design interfaces onto idealized repeat proteins that direct their assembly into complexes that possess cyclic symmetry. Of 96 designs that were experimentally characterized, 21 were found to form stable monodisperse homo-oligomers in solution, and 15 (4 homodimers, 6 homotrimers, 6 homotetramers and 1 homopentamer) had solution small angle X-ray scattering data consistent with the design models. X-ray crystal structures were obtained for five of the designs and each of these were shown to be very close to their design model.

ClusPro-DC: Dimer Classification by the Cluspro Server for Protein-Protein Docking

ClusPro-DC (https://cluspro.bu.edu/) implements a straightforward approach to the discrimination between crystallographic and biological dimers by docking the two subunits to exhaustively sample the interaction energy landscape. If a substantial number of low energy docked poses cluster in a narrow vicinity of the native structure of the dimer, then one can assume that there is a well-defined free energy well around the native state, which makes the interaction stable. In contrast, if the interaction sites in the docked poses do not form a large enough cluster around the native structure, then it is unlikely that the subunits form a stable biological dimer. The number of near-native structures is used to estimate the probability of a dimer being biological. Currently, the server examines only the stability of a given interface rather than generating all putative quaternary structures as accomplished by PISA or EPPIC, but it complements the information provided by these methods.

Multiple Domain Associations within the Arabidopsis Immune Receptor RPP1 Regulate the Activation of Programmed Cell Death

Upon recognition of pathogen virulence effectors, plant nucleotide-binding leucine-rich repeat (NLR) proteins induce defense responses including localized host cell death. In an effort to understand the molecular mechanisms leading to this response, we examined the Arabidopsis thaliana NLR protein RECOGNITION OF PERONOSPORA PARASITICA1 (RPP1), which recognizes the Hyaloperonospora arabidopsidis effector ARABIDOPSIS THALIANA RECOGNIZED1 (ATR1). Expression of the N-terminus of RPP1, including the Toll/interleukin-1 receptor (TIR) domain (“N-TIR”), elicited an effector-independent cell death response, and we used allelic variation in TIR domain sequences to define the key residues that contribute to this phenotype. Further biochemical characterization indicated that cell death induction was correlated with N-TIR domain self-association. In addition, we demonstrated that the nucleotide-binding (NB)-ARC1 region of RPP1 self-associates and plays a critical role in cell death activation, likely by facilitating TIR:TIR interactions. Structural homology modeling of the NB subdomain allowed us to identify a putative oligomerization interface that was shown to influence NB-ARC1 self-association. Significantly, full-length RPP1 exhibited effector-dependent oligomerization and, although mutations at the NB-ARC1 oligomerization interface eliminated cell death induction, RPP1 self-association was unaffected, suggesting that additional regions contribute to oligomerization. Indeed, the leucine-rich repeat domain of RPP1 also self-associates, indicating that multiple interaction interfaces exist within activated RPP1 oligomers. Finally, we observed numerous intramolecular interactions that likely function to negatively regulate RPP1, and present a model describing the transition to an active NLR protein.

Many plant pathogens inject proteins known as effectors into the cells of their hosts in order to suppress host immune responses and promote pathogen growth. Over time, plants have evolved receptors, described as nucleotide-binding leucine-rich repeat (NLR) proteins, which recognize the activity of pathogen effectors and stimulate defense responses. Plant NLRs contain several domains that exhibit striking functional conservation with NLRs from other eukaryotes. Despite their important contribution to plant immunity, the molecular mechanisms that underlie effector recognition and subsequent immune activation by NLRs remain to be fully elucidated. Here, we focus on RPP1, an NLR from Arabidopsis that recognizes the oomycete effector ATR1. Using transient co-expression of proteins in plants, we demonstrate that recognition of ATR1 stimulates RPP1 oligomerization. This interaction involves multiple domains of RPP1 and is critical for immune activation. In the absence of ATR1, we documented interactions between domains within an individual RPP1 protein, likely occurring to prevent inappropriate immune activation. Finally, we examined differences between RPP1 alleles as well as structural data from animal NLRs to help identify specific amino acids that mediate interactions within and between RPP1 molecules. Collectively, these data allow us to propose a model for the activation of RPP1 following ATR1 recognition.

Comparative analyses of quaternary arrangements in homo-oligomeric proteins in superfamilies: Functional implications

A comprehensive analysis of the quaternary features of distantly related homo-oligomeric proteins is the focus of the current study. This study has been performed at the levels of quaternary state, symmetry, and quaternary structure. Quaternary state and quaternary structure refers to the number of subunits and spatial arrangements of subunits, respectively. Using a large dataset of available 3D structures of biologically relevant assemblies, we show that only 53% of the distantly related homo-oligomeric proteins have the same quaternary state. Considering these homologous homo-oligomers with the same quaternary state, conservation of quaternary structures is observed only in 38% of the pairs. In 36% of the pairs of distantly related homo-oligomers with different quaternary states the larger assembly in a pair shows high structural similarity with the entire quaternary structure of the related protein with lower quaternary state and it is referred as "Russian doll effect." The differences in quaternary state and structure have been suggested to contribute to the functional diversity. Detailed investigations show that even though the gross functions of many distantly related homo-oligomers are the same, finer level differences in molecular functions are manifested by differences in quaternary states and structures. Comparison of structures of biological assemblies in distantly and closely related homo-oligomeric proteins throughout the study differentiates the effects of sequence divergence on the quaternary structures and function. Knowledge inferred from this study can provide insights for improved protein structure classification and function prediction of homo-oligomers. Proteins 2016; 84:1190-1202. © 2016 Wiley Periodicals, Inc.

Exploring the structure of glutamate racemase from Mycobacterium tuberculosis as a template for anti-mycobacterial drug discovery

Glutamate racemase (MurI) is responsible for providing D-glutamate for peptidoglycan biosynthesis in bacteria and has been a favoured target in pharmaceutical drug design efforts. It has recently been proven to be essential in Mycobacterium tuberculosis, the causative organism of tuberculosis, a disease for which new medications are urgently needed. In the present study, we have determined the protein crystal structures of MurI from both M. tuberculosis and Mycobacterium smegmatis in complex with D-glutamate to 2.3 Å and 1.8 Å resolution respectively. These structures are conserved, but reveal differences in their active site architecture compared with that of other MurI structures. Furthermore, compounds designed to target other glutamate racemases have been screened but do not inhibit mycobacterial MurI, suggesting that a new drug design effort will be needed to develop inhibitors. A new type of MurI dimer arrangement has been observed in both structures, and this arrangement becomes the third biological dimer geometry for MurI found to date. The mycobacterial MurI dimer is tightly associated, with a KD in the nanomolar range. The enzyme binds D- and L-glutamate specifically, but is inactive in solution unless the dimer interface is mutated. We created triple mutants of this interface in the M. smegmatis glutamate racemase (D26R/R105A/G194R or E) that have appreciable activity (kcat=0.056-0.160 min(-1) and KM=0.26-0.51 mM) and can be utilized to screen proposed antimicrobial candidates for inhibition.

Local Geometry and Evolutionary Conservation of Protein Surfaces Reveal the Multiple Recognition Patches in Protein-Protein Interactions

Protein-protein interactions (PPIs) are essential to all biological processes and they represent increasingly important therapeutic targets. Here, we present a new method for accurately predicting protein-protein interfaces, understanding their properties, origins and binding to multiple partners. Contrary to machine learning approaches, our method combines in a rational and very straightforward way three sequence- and structure-based descriptors of protein residues: evolutionary conservation, physico-chemical properties and local geometry. The implemented strategy yields very precise predictions for a wide range of protein-protein interfaces and discriminates them from small-molecule binding sites. Beyond its predictive power, the approach permits to dissect interaction surfaces and unravel their complexity. We show how the analysis of the predicted patches can foster new strategies for PPIs modulation and interaction surface redesign. The approach is implemented in JET2, an automated tool based on the Joint Evolutionary Trees (JET) method for sequence-based protein interface prediction. JET2 is freely available at www.lcqb.upmc.fr/JET2.