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.

Comprehensive mapping of mutations in TDP-43 and α-Synuclein that affect stability and binding

Abnormal aggregation and amyloid inclusions of TAR DNA-binding protein 43 (TDP-43) and α-Synuclein (α-Syn) are frequently co-observed in amyotrophic lateral sclerosis, Parkinson's disease, and Alzheimer's disease. Several reports showed TDP-43 C-terminal domain (CTD) and α-Syn interact with each other and the aggregates of these two proteins colocalized together in different cellular and animal models. Molecular dynamics simulation was conducted to elucidate the stability of the TDP-43 and Syn complex structure. The interfacial mutations in protein complexes changes the stability and binding affinity of the protein that may cause diseases. Here, we have utilized the computational saturation mutagenesis approach including structure-based stability and binding energy calculations to compute the systemic effects of missense mutations of TDP-43 CTD and α-Syn on protein stability and binding affinity. Most of the interfacial mutations of CTD and α-Syn were found to destabilize the protein and reduced the protein binding affinity. The results thus shed light on the functional consequences of missense mutations observed in TDP-43 associated proteinopathies and may provide the mechanisms of co-morbidities involving these two proteins.Communicated by Ramaswamy H. Sarma.

ComplexQA: a deep graph learning approach for protein complex structure assessment

MOTIVATION

In recent years, the end-to-end deep learning method for single-chain protein structure prediction has achieved high accuracy. For example, the state-of-the-art method AlphaFold, developed by Google, has largely increased the accuracy of protein structure predictions to near experimental accuracy in some of the cases. At the same time, there are few methods that can evaluate the quality of protein complexes at the residue level. In particular, evaluating the quality of residues at the interface of protein complexes can lead to a wide range of applications, such as protein function analysis and drug design. In this paper, we introduce a new deep graph neural network-based method ComplexQA, to evaluate the local quality of interfaces for protein complexes by utilizing the residue-level structural information in 3D space and the sequence-level constraints.

RESULTS

We benchmark our method to other state-of-the-art quality assessment approaches on the HAF2 and DBM55-AF2 datasets (high-quality structural models predicted by AlphaFold-Multimer), and the BM5 docking dataset. The experimental results show that our proposed method achieves better or similar performance compared with other state-of-the-art methods, especially on difficult targets which only contain a few acceptable models. Our method is able to suggest a score for each interfac e residue, which demonstrates a powerful assessment tool for the ever-increasing number of protein complexes.

AVAILABILITY

https://github.com/Cao-Labs/ComplexQA.git. Contact: caora@plu.edu.

Interface engineering of cellobiose dehydrogenase improves interdomain electron transfer

Cellobiose dehydrogenase (CDH) is a bioelectrocatalyst that enables direct electron transfer (DET) in biosensors and biofuel cells. The application of this bidomain hemoflavoenzyme for physiological glucose measurements is limited by its acidic pH optimum and slow interdomain electron transfer (IET) at pH 7.5. The reason for this rate-limiting electron transfer step is electrostatic repulsion at the interface between the catalytic dehydrogenase domain and the electron mediating cytochrome domain (CYT). We applied rational interface engineering to accelerate the IET for the pH prevailing in blood or interstitial fluid. Phylogenetic and structural analyses guided the design of 17 variants in which acidic amino acids were mutated at the CYT domain. Five mutations (G71K, D160K, Q174K, D177K, M180K) increased the pH optimum and IET rate. Structure-based analysis of the variants suggested two mechanisms explaining the improvements: electrostatic steering and stabilization of the closed state by hydrogen bonding. Combining the mutations into six combinatorial variants with up to five mutations shifted the pH optimum from 4.5 to 7.0 and increased the IET at pH 7.5 over 12-fold from 0.1 to 1.24 s^-1 . While the mutants sustained a high enzymatic activity and even surpassed the IET of the wild-type enzyme, the accumulated positive charges on the CYT domain decreased DET, highlighting the importance of CYT for IET and DET. This study shows that interface engineering is an effective strategy to shift the pH optimum and improve the IET of CDH, but future work needs to maintain the DET of the CYT domain for bioelectronic applications. This article is protected by copyright. All rights reserved.

Evolution of homo-oligomerization of methionine S-adenosyltransferases is replete with structure-function constrains

Homomers are prevalent in bacterial proteomes, particularly among core metabolic enzymes. Homomerization is often key to function and regulation, and interfaces that facilitate the formation of homomeric enzymes are subject to intense evolutionary change. However, our understanding of the molecular mechanisms that drive evolutionary variation in homomeric complexes is still lacking. How is the diversification of protein interfaces linked to variation in functional regulation and structural integrity of homomeric complexes? To address this question, we studied quaternary structure evolution of bacterial methionine S‐adenosyltransferases (MATs)—dihedral homotetramers formed along a large and conserved dimeric interface harboring two active sites, and a small, recently evolved, interdimeric interface. Here, we show that diversity in the physicochemical properties of small interfaces is directly linked to variability in the kinetic stability of MAT quaternary complexes and in modes of their functional regulation. Specifically, hydrophobic interactions within the small interface of Escherichia coli MAT render the functional homotetramer kinetically stable yet impose severe aggregation constraints on complex assembly. These constraints are alleviated by electrostatic interactions that accelerate dimer‐dimer assembly. In contrast, Neisseria gonorrhoeae MAT adopts a nonfunctional dimeric state due to the low hydrophobicity of its small interface and the high flexibility of its active site loops, which perturbs small interface integrity. Remarkably, in the presence of methionine and ATP, N. gonorrhoeae MAT undergoes substrate‐induced assembly into a functional tetrameric state. We suggest that evolution acts on the interdimeric interfaces of MATs to tailor the regulation of their activity and stability to unique organismal needs.

PDB Code(s): 7R2W;7R3B;

Boosting the analysis of protein interfaces with multiple interface string alignments: Illustration on the spikes of coronaviruses

We introduce multiple interface string alignment (MISA), a visualization tool to display coherently various sequence and structure based statistics at protein-protein interfaces (SSE elements, buried surface area, Δ ASA , B factor values, etc). The amino acids supporting these annotations are obtained from Voronoi interface models. The benefit of MISA is to collate annotated sequences of (homologous) chains found in different biological contexts, that is, bound with different partners or unbound. The aggregated views MISA/SSE, MISA/BSA, MISA/ΔASA, and so forth, make it trivial to identify commonalities and differences between chains, to infer key interface residues, and to understand where conformational changes occur upon binding. As such, they should prove of key relevance for knowledge-based annotations of protein databases such as the Protein Data Bank. Illustrations are provided on the receptor binding domain of coronaviruses, in complex with their cognate partner or (neutralizing) antibodies. MISA computed with a minimal number of structures complement and enrich findings previously reported. The corresponding package is available from the Structural Bioinformatics Library (http://sbl.inria.frand https://sbl.inria.fr/doc/Multiple_interface_string_alignment-user-manual.html).

Structure-Guided Computational Approaches to Unravel Druggable Proteomic Landscape of Mycobacterium leprae

Leprosy, caused by Mycobacterium leprae (M. leprae), is treated with a multidrug regimen comprising Dapsone, Rifampicin, and Clofazimine. These drugs exhibit bacteriostatic, bactericidal and anti-inflammatory properties, respectively, and control the dissemination of infection in the host. However, the current treatment is not cost-effective, does not favor patient compliance due to its long duration (12 months) and does not protect against the incumbent nerve damage, which is a severe leprosy complication. The chronic infectious peripheral neuropathy associated with the disease is primarily due to the bacterial components infiltrating the Schwann cells that protect neuronal axons, thereby inducing a demyelinating phenotype. There is a need to discover novel/repurposed drugs that can act as short duration and effective alternatives to the existing treatment regimens, preventing nerve damage and consequent disability associated with the disease. Mycobacterium leprae is an obligate pathogen resulting in experimental intractability to cultivate the bacillus in vitro and limiting drug discovery efforts to repositioning screens in mouse footpad models. The dearth of knowledge related to structural proteomics of M. leprae, coupled with emerging antimicrobial resistance to all the three drugs in the multidrug therapy, poses a need for concerted novel drug discovery efforts. A comprehensive understanding of the proteomic landscape of M. leprae is indispensable to unravel druggable targets that are essential for bacterial survival and predilection of human neuronal Schwann cells. Of the 1,614 protein-coding genes in the genome of M. leprae, only 17 protein structures are available in the Protein Data Bank. In this review, we discussed efforts made to model the proteome of M. leprae using a suite of software for protein modeling that has been developed in the Blundell laboratory. Precise template selection by employing sequence-structure homology recognition software, multi-template modeling of the monomeric models and accurate quality assessment are the hallmarks of the modeling process. Tools that map interfaces and enable building of homo-oligomers are discussed in the context of interface stability. Other software is described to determine the druggable proteome by using information related to the chokepoint analysis of the metabolic pathways, gene essentiality, homology to human proteins, functional sites, druggable pockets and fragment hotspot maps.

Molecular and Structural Basis of Cross-Reactivity in M. tuberculosis Toxin-Antitoxin Systems

Mycobacterium tuberculosis genome encodes over 80 toxin-antitoxin (TA) systems. While each toxin interacts with its cognate antitoxin, the abundance of TA systems presents an opportunity for potential non-cognate interactions. TA systems mediate manifold interactions to manage pathogenicity and stress response network of the cell and non-cognate interactions may play vital roles as well. To address if non-cognate and heterologous interactions are feasible and to understand the structural basis of their interactions, we have performed comprehensive computational analyses on the available 3D structures and generated structural models of paralogous M. tuberculosis VapBC and MazEF TA systems. For a majority of the TA systems, we show that non-cognate toxin-antitoxin interactions are structurally incompatible except for complexes like VapBC15 and VapBC11, which show similar interfaces and potential for cross-reactivity. For TA systems which have been experimentally shown earlier to disfavor non-cognate interactions, we demonstrate that they are structurally and stereo-chemically incompatible. For selected TA systems, our detailed structural analysis identifies specificity conferring residues. Thus, our work improves the current understanding of TA interfaces and generates a hypothesis based on congenial binding site, geometric complementarity, and chemical nature of interfaces. Overall, our work offers a structure-based explanation for non-cognate toxin-antitoxin interactions in M. tuberculosis.

Conserved Rotavirus NSP5 and VP2 Domains Interact and Affect Viroplasm

One step of the life cycle common to all rotaviruses (RV) studied so far is the formation of viroplasms, membrane-less cytosolic inclusions providing a microenvironment for early morphogenesis and RNA replication. Viroplasm-like structures (VLS) are simplified viroplasm models consisting of complexes of nonstructural protein 5 (NSP5) with the RV core shell VP2 or NSP2. We identified and characterized the domains required for NSP5-VP2 interaction and VLS formation. VP2 mutations L124A, V865A, and I878A impaired both NSP5 hyperphosphorylation and NSP5/VP2 VLS formation. Moreover, NSP5-VP2 interaction does not depend on NSP5 hyperphosphorylation. The NSP5 tail region is required for VP2 interaction. Notably, VP2 L124A expression acts as a dominant-negative element by disrupting the formation of either VLS or viroplasms and blocking RNA synthesis. In silico analyses revealed that VP2 L124, V865, and I878 are conserved among RV species A to H. Detailed knowledge of the protein interaction interface required for viroplasm formation may facilitate the design of broad-spectrum antivirals to block RV replication.IMPORTANCE Alternative treatments to combat rotavirus infection are a requirement for susceptible communities where vaccines cannot be applied. This demand is urgent for newborn infants, immunocompromised patients, adults traveling to high-risk regions, and even for the livestock industry. Aside from structural and physiological divergences among RV species studied before now, all replicate within cytosolic inclusions termed viroplasms. These inclusions are composed of viral and cellular proteins and viral RNA. Viroplasm-like structures (VLS), composed of RV protein NSP5 with either NSP2 or VP2, are models for investigating viroplasms. In this study, we identified a conserved amino acid in the VP2 protein, L124, necessary for its interaction with NSP5 and the formation of both VLSs and viroplasms. As RV vaccines cover a narrow range of viral strains, the identification of VP2 L124 residue lays the foundations for the design of drugs that specifically block NSP5-VP2 interaction as a broad-spectrum RV antiviral.

Soft regions of protein surface are potent for stable dimer formation

By having knowledge about the characteristics of protein interaction interfaces, we will be able to manipulate protein complexes for therapies. Dimer state is considered as the primary alphabet of the most proteins' quaternary structure. The properties of binding interface between subunits and of noninterface region define the specificity and stability of the intended protein complex. Considering some topological properties and amino acids' affinity for binding in interfaces of protein dimers, we construct the interface-specific recurrence plots. The data obtained from recurrence quantitative analysis, and accessibility-related metrics help us to classify the protein dimers into four distinct classes. Some mechanical properties of binding interfaces are computed for each predefined class of the dimers. The computed mechanical characteristics of binding patch region are compared with those of nonbinding region of proteins. Our observations indicate that the mechanical properties of protein binding sites have a decisive impact on determining the dimer stability. We introduce a new concept in analyzing protein structure by considering mechanical properties of protein structure. We conclude that the interface region between subunits of stable dimers is usually mechanically softer than the interface of unstable protein dimers. AbbreviationsAABaverage affinity for bindingANManisotropic network modelAPCaffinity propagation clusteringASAaccessible surface areaCCDinter residues distanceCSCcomplex stability codeDMdistance matrixΔG_dissPISA-computed dissociation free energyGNMGaussian normal mode analysisNMAnormal mode analysisPBPprotein binding patchPISAproteins, interfaces, structures and assembliesrASArelative accessible area in respect to unfolded state of residuesRMrecurrence matrixrPrelative protrusionRPrecurrence plotRQArecurrence quantitative analysisSEMstandard error of meanCommunicated by Ramaswamy H. Sarma.

Nervous-Like Circuits in the Ribosome Facts, Hypotheses and Perspectives

In the past few decades, studies on translation have converged towards the metaphor of a “ribosome nanomachine”; they also revealed intriguing ribosome properties challenging this view. Many studies have shown that to perform an accurate protein synthesis in a fluctuating cellular environment, ribosomes sense, transfer information and even make decisions. This complex “behaviour” that goes far beyond the skills of a simple mechanical machine has suggested that the ribosomal protein networks could play a role equivalent to nervous circuits at a molecular scale to enable information transfer and processing during translation. We analyse here the significance of this analogy and establish a preliminary link between two fields: ribosome structure-function studies and the analysis of information processing systems. This cross-disciplinary analysis opens new perspectives about the mechanisms of information transfer and processing in ribosomes and may provide new conceptual frameworks for the understanding of the behaviours of unicellular organisms.

Local Interaction Signal Analysis Predicts Protein-Protein Binding Affinity

Several models estimating the strength of the interaction between proteins in a complex have been proposed. By exploring the geometry of contact distribution at protein-protein interfaces, we provide an improved model of binding energy. Local interaction signal analysis (LISA) is a radial function based on terms describing favorable and non-favorable contacts obtained by density functional theory, the support-core-rim interface residue distribution, non-interacting charged residues and secondary structures contribution. The three-dimensional organization of the contacts and their contribution on localized hot-sites over the entire interaction surface were numerically evaluated. LISA achieves a correlation of 0.81 (and a root-mean-square error of 2.35 ± 0.38 kcal/mol) when tested on 125 complexes for which experimental measurements were realized. LISA's performance is stable for subsets defined by functional composition and extent of conformational changes upon complex formation. A large-scale comparison with 17 other functions demonstrated the power of the geometrical model in the understanding of complex binding.

The relationship between relative solvent accessible surface area (rASA) and irregular structures in protean segments (ProSs)

Intrinsically Disordered Proteins (IDPs) lack a stable, three-dimensional structure under physiological conditions, yet they exhibit numerous biological activities. Protean segments (ProSs) are the functional regions of intrinsically disordered proteins that undergo disorder-to-order transitions upon binding to their partners. Example ProSs collected from the intrinsically disordered proteins with extensive annotations and literature (IDEAL) database. The interface of protean segments (ProSs) is classified into core, rim, and support, and analyzed their secondary structure elements (SSEs) based on the relative accessible surface area (rASA). The amino acid compositions and the relative solvent accessible surface areas (rASAs) of ProS secondary structural elements (SSEs) at the interface, core and rim were compared to those of heterodimers. The average number of contacts of alpha helices and irregular residues was calculated for each ProS and heterodimer. Furthermore, the ProSs were classified into high and low efficient based on their average number of contacts at the interface. The results indicate that the irregular structures of ProSs and heterodimers are significantly different. The rASA of irregular structures in the monomeric state (rASAm) is large, leads to the formation of larger ΔrASA and many contacts in ProSs.

Polypeptide modification: an improved proglycinin design to stabilise oil-in-water emulsions

β-Conglycinin and glycinin are soybean major seed storage proteins. Previous studies have shown that adding the extension region of β-conglycinin α subunit improves the emulsifying properties of proglycinin and confers more favourable characteristics than fusing the extension region of β-conglycinin α' subunit or the hypervariable regions (A4IV) of glycinin A1aB1b subunit. To evaluate the polypeptide properties, we designed mutants of A1aB1b subunits fused with truncated versions of A4IV (A4IVcut), α (αcut) or α' (α'cut) extension regions lacking the C-terminus 25 or 31 residues (A4IVC25, αC25 or α'C31), and also A4IVcut and α'cut with αC25 residues added (A4IVcut-αC25 and α'cut-αC25). All the modified proteins displayed conformations similar to the wild type. With good solubilities, the emulsion properties of the modified proteins were much better at ionic strength μ = 0.08 than at μ = 0.5. The modified A1aB1bαcut and A1aB1bα'cut showed poorer emulsion properties than those of A1aB1bα and A1aB1bα'. Replacing the hydrophobic A4IVC25 region of A1aB1bA4IV with hydrophilic αC25 created A1aB1bA4IVcut-αC25, which had the best emulsion stability among these proglycinin mutants. We found that addition of αC25 improves the emulsifying properties of two C-terminally truncated proglycinin variants, thereby illustrating its potential general utility. Our investigation showed that in order to improve the emulsifying ability and emulsion stability of a globular protein, the introduced polypeptide should (i) be highly hydrophilic, (ii) consist of multiple hydrophobic-strong hydrophilic regions comprising at least two alpha helixes, (iii) harbour a terminal α-helix at the end of the C-terminus and (iv) have properties similar to those of αC25.

A tale of tails: deciphering the contribution of terminal tails to the biochemical properties of two Dps proteins from Streptomyces coelicolor

Dps proteins are members of an extensive family of proteins that oxidise and deposit iron in the form of ferric oxide, and are also able to bind DNA. Ferroxidation centres are formed at the interface of anti-parallel dimers, which further assemble into dodecameric nanocages with a hollow core where ferric oxide is deposited. Streptomyces coelicolor encodes three Dps-like proteins (DpsA, B and C). Despite sharing the conserved four-helix bundle organisation observed in members of the Dps family, they display significant differences in the length of terminal extensions, or tails. DpsA possess both N- and C-terminal tails of different lengths, and their removal affects quaternary structure assembly to varying degrees. DpsC quaternary structure, on the other hand, is heavily dependent on its N-terminal tail as its removal abolishes correct protein folding. Analysis of the crystal structure of dodecamers from both proteins revealed remarkable differences in the position of tails and interface surface area; and provides insight to explain the differences in biochemical behaviour observed while comparing DpsA and DpsC.

Landscape of intertwined associations in multi-domain homo-oligomeric proteins

This study charts the landscape of multi-domain protein structures that form intertwined homodimers by exchanging structural domains between subunits. A representative dataset of such homodimers was derived from the Protein Data Bank, and their structural and topological properties were compared to those of a representative set of non-intertwined homodimers. Most of the intertwined dimers form closed assemblies with head-to-tail arrangements, where the subunit interface involves contacts between dissimilar domains. In contrast, the non-intertwined dimers form preferentially head-to-head arrangements, where the subunit interface involves contacts between identical domains. Most of these contacts engage only one structural domain from each subunit, leaving the remaining domains free to form other associations. Remarkably, we find that multi-domain proteins closely related to the intertwined homodimers are significantly more likely than relatives of the non-intertwined versions to adopt alternative intramolecular domain arrangements. In ~40% of the intertwined dimers, the plasticity in domain arrangements among relatives affords maintenance of the head-to-head or head-to-tail topology and conservation of the corresponding subunit interface. This property seems to be exploited in several systems to regulate DNA binding. In ~58%, however, intramolecular domain re-arrangements are associated with changes in oligomeric states and poorly conserved interfaces among relatives. This time, the corresponding structural plasticity appears to be exploited by evolution to modulate function by switching between active and inactive states of the protein. Surprisingly, in total, only three systems were found to undergo the classical monomer to intertwined dimer conversion associated with three-dimensional domain swapping.

Blind prediction of interfacial water positions in CAPRI

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

Chemical characteristics of dimer interfaces in the legume lectin family

The Erythrina corallodendron lectin (EcorL) crystallizes in monoclinic and hexagonal crystal forms. Comparison of the newly determined hexagonal form (PDB code 1fyu) with the monoclinic form shows that the dimeric structure of EcorL reflects the inherent biological structure of the protein and is not an artifact of the crystal packing. To further understand the factors determining the dimerization modes of legume lectins, EcorL, concanavalin A (ConA), and Griffonia simplicifolia (GS4) were taken as representatives of the three unique dimers found in the family. Six virtual homodimers were generated. The hydropathy, amino acid composition, and solvation energy were calculated for all nine homodimers. Each of the three native dimers has a distinct chemical composition. EcorL has a dominant hydrophobic component, and ConA has a strong polar component, but in GS4 the three components contribute equally to the interface. This distribution pattern at the interface is unique to the native dimers and distinct from the partition observed in the virtual dimers. Amino acid composition of other members of the family that dimerize like EcorL or ConA maintain the same pattern of amino acids distribution observed in EcorL and ConA. However, lectins that dimerize like GS4 do not show a particularly distinct distribution. In all cases, the calculated solvation energy of the native dimer was lower than that of the virtual dimers, suggesting that the observed mode of dimerization is the most stable organization for the given sequence and tertiary structure. The dimerization type cannot be predicted by sequence analysis.