Deorphanizing Peptides Using Structure Prediction

Inhibitory effect of truncated isoforms on GPCR dimerization predicted by combinatorial computational strategy

G protein-coupled receptors (GPCRs) play a pivotal role in fundamental biological processes and disease development. GPCR isoforms, derived from alternative splicing, can exhibit distinct signaling patterns. Some highly-truncated isoforms can impact functional performance of full-length receptors, suggesting their intriguing regulatory roles. However, how these truncated isoforms interact with full-length counterparts remains largely unexplored. Here, we computationally investigated the interaction patterns of three human GPCRs from three different classes, ADORA1 (Class A), mGlu2 (Class C) and SMO (Class F) with their respective truncated isoforms because their homodimer structures have been experimentally determined, and they have truncated isoforms deposited and identified at protein level in Uniprot database. Combining the neural network-based AlphaFold2 and two physics-based protein-protein docking tools, we generated multiple complex structures and assessed the binding affinity in the context of atomistic molecular dynamics simulations. Our computational results suggested all the four studied truncated isoforms showed potent binding to their counterparts and overlapping interfaces with homodimers, indicating their strong potential to block homodimerization of their counterparts. Our study offers insights into functional significance of GPCR truncated isoforms and supports the ubiquity of their regulatory roles.

Identification of type VI secretion system effector-immunity pairs using structural bioinformatics

The type VI secretion system (T6SS) is an important mediator of microbe-microbe and microbe-host interactions. Gram-negative bacteria use the T6SS to inject T6SS effectors (T6Es), which are usually proteins with toxic activity, into neighboring cells. Antibacterial effectors have cognate immunity proteins that neutralize self-intoxication. Here, we applied novel structural bioinformatic tools to perform systematic discovery and functional annotation of T6Es and their cognate immunity proteins from a dataset of 17,920 T6SS-encoding bacterial genomes. Using structural clustering, we identified 517 putative T6E families, outperforming sequence-based clustering. We developed a logistic regression model to reliably quantify protein-protein interaction of new T6E-immunity pairs, yielding candidate immunity proteins for 231 out of the 517 T6E families. We used sensitive structure-based annotation which yielded functional annotations for 51% of the T6E families, again outperforming sequence-based annotation. Next, we validated four novel T6E-immunity pairs using basic experiments in E. coli. In particular, we showed that the Pfam domain DUF3289 is a homolog of Colicin M and that DUF943 acts as its cognate immunity protein. Furthermore, we discovered a novel T6E that is a structural homolog of SleB, a lytic transglycosylase, and identified a specific glutamate that acts as its putative catalytic residue. Overall, this study applies novel structural bioinformatic tools to T6E-immunity pair discovery, and provides an extensive database of annotated T6E-immunity pairs.

PPIscreenML: Structure-based screening for protein-protein interactions using AlphaFold

Protein-protein interactions underlie nearly all cellular processes. With the advent of protein structure prediction methods such as AlphaFold2 (AF2), models of specific protein pairs can be built extremely accurately in most cases. However, determining the relevance of a given protein pair remains an open question. It is presently unclear how to use best structure-based tools to infer whether a pair of candidate proteins indeed interact with one another: ideally, one might even use such information to screen amongst candidate pairings to build up protein interaction networks. Whereas methods for evaluating quality of modeled protein complexes have been co-opted for determining which pairings interact (e.g., pDockQ and iPTM), there have been no rigorously benchmarked methods for this task. Here we introduce PPIscreenML, a classification model trained to distinguish AF2 models of interacting protein pairs from AF2 models of compelling decoy pairings. We find that PPIscreenML out-performs methods such as pDockQ and iPTM for this task, and further that PPIscreenML exhibits impressive performance when identifying which ligand/receptor pairings engage one another across the structurally conserved tumor necrosis factor superfamily (TNFSF). Analysis of benchmark results using complexes not seen in PPIscreenML development strongly suggest that the model generalizes beyond training data, making it broadly applicable for identifying new protein complexes based on structural models built with AF2.

An Effective Plant Small Secretory Peptide Recognition Model Based on Feature Correction Strategy

Plant small secretory peptides (SSPs) play an important role in the regulation of biological processes in plants. Accurately predicting SSPs enables efficient exploration of their functions. Traditional experimental verification methods are very reliable and accurate, but they require expensive equipment and a lot of time. The method of machine learning speeds up the prediction process of SSPs, but the instability of feature extraction will also lead to further limitations of this type of method. Therefore, this paper proposes a new feature-correction-based model for SSP recognition in plants, abbreviated as SE-SSP. The model mainly includes the following three advantages: First, the use of transformer encoders can better reveal implicit features. Second, design a feature correction module suitable for sequences, named 2-D SENET, to adaptively adjust the features to obtain a more robust feature representation. Third, stack multiple linear modules to further dig out the deep information on the sample. At the same time, the training based on a contrastive learning strategy can alleviate the problem of sparse samples. We construct experiments on publicly available data sets, and the results verify that our model shows an excellent performance. The proposed model can be used as a convenient and effective SSP prediction tool in the future. Our data and code are publicly available at https://github.com/wrab12/SE-SSP/.

Systematic discovery of protein interaction interfaces using AlphaFold and experimental validation

Structural resolution of protein interactions enables mechanistic and functional studies as well as interpretation of disease variants. However, structural data is still missing for most protein interactions because we lack computational and experimental tools at scale. This is particularly true for interactions mediated by short linear motifs occurring in disordered regions of proteins. We find that AlphaFold-Multimer predicts with high sensitivity but limited specificity structures of domain-motif interactions when using small protein fragments as input. Sensitivity decreased substantially when using long protein fragments or full length proteins. We delineated a protein fragmentation strategy particularly suited for the prediction of domain-motif interfaces and applied it to interactions between human proteins associated with neurodevelopmental disorders. This enabled the prediction of highly confident and likely disease-related novel interfaces, which we further experimentally corroborated for FBXO23-STX1B, STX1B-VAMP2, ESRRG-PSMC5, PEX3-PEX19, PEX3-PEX16, and SNRPB-GIGYF1 providing novel molecular insights for diverse biological processes. Our work highlights exciting perspectives, but also reveals clear limitations and the need for future developments to maximize the power of Alphafold-Multimer for interface predictions.

Novel gurmarin-like peptides from Gymnema sylvestre and their interactions with the sweet taste receptor T1R2/T1R3

Gymnema sylvestre (GS) is a traditional medicinal plant known for its hypoglycemic and hypolipidemic effects. Gurmarin (hereafter Gur-1) is the only known active peptide in GS. Gur-1 has a suppressive sweet taste effect in rodents but no or only a very weak effect in humans. Here, 8 gurmarin-like peptides (Gur-2 to Gur-9) and their isoforms are reported in the GS transcriptome. The molecular mechanism of sweet taste suppression by Gur-1 is still largely unknown. Therefore, the complete architecture of human and mouse sweet taste receptors T1R2/T1R3 and their interaction with Gur-1 to Gur-9 were predicted by AlphaFold-Multimer (AF-M) and validated. Only Gur-1 and Gur-2 interact with the T1R2/T1R3 receptor. Indeed, Gur-1 and Gur-2 bind to the region of the cysteine-rich domain (CRD) and the transmembrane domain (TMD) of the mouse T1R2 subunit. In contrast, only Gur-2 binds to the TMD of the human T1R2 subunit. This result suggests that Gur-2 may have a suppressive sweet taste effect in humans. Furthermore, AF-M predicted that Gα-gustducin, a protein involved in sweet taste transduction, interacts with the intracellular domain of the T1R2 subunit. These results highlight an unexpected diversity of gurmarin-like peptides in GS and provide the complete predicted architecture of the human and mouse sweet taste receptor with the putative binding sites of Gur-1, Gur-2, and Gα-gustducin. In addition, gurmarin-like peptides may serve as promising drug scaffolds for the development of antidiabetic molecules.

Assessing the fast prediction of peptide conformers and the impact of non-natural modifications

We present an assessment of different approaches to predict peptide structures using modeling tools. Several small molecule, protein, and peptide-focused methodologies were used for the fast prediction of conformers for peptides shorter than 30 amino acids. We assessed the effect of including restraints based on annotated or predicted secondary structure motifs. A number of peptides in bound conformations and in solution were collected to compare the tools. In addition, we studied the impact of changing single amino acids to non-natural residues using molecular dynamics simulations. Deep learning methods such as AlphaFold2, or the combination of physics-based approaches with secondary structure information, produce the most accurate results for natural sequences. In the case of peptides with non-natural modifications, modeling the peptide containing natural amino acids first and then modifying and simulating the peptide using benchmarked force fields is a recommended pipeline. The results can guide the modeling of oligopeptides for drug discovery projects.

Dynamic Dimerization of Chemokine Receptors and Potential Inhibitory Role of Their Truncated Isoforms Revealed through Combinatorial Prediction

Chemokine receptors play crucial roles in fundamental biological processes. Their malfunction may result in many diseases, including cancer, autoimmune diseases, and HIV. The oligomerization of chemokine receptors holds significant functional implications that directly affect their signaling patterns and pharmacological responses. However, the oligomerization patterns of many chemokine receptors remain poorly understood. Furthermore, several chemokine receptors have highly truncated isoforms whose functional role is not yet clear. Here, we computationally show homo- and heterodimerization patterns of four human chemokine receptors, namely CXCR2, CXCR7, CCR2, and CCR7, along with their interaction patterns with their respective truncated isoforms. By combining the neural network-based AlphaFold2 and physics-based protein-protein docking tool ClusPro, we predicted 15 groups of complex structures and assessed the binding affinities in the context of atomistic molecular dynamics simulations. Our results are in agreement with previous experimental observations and support the dynamic and diverse nature of chemokine receptor dimerization, suggesting possible patterns of higher-order oligomerization. Additionally, we uncover the strong potential of truncated isoforms to block homo- and heterodimerization of chemokine receptors, also in a dynamic manner. Our study provides insights into the dimerization patterns of chemokine receptors and the functional significance of their truncated isoforms.

Ligand recognition and G protein coupling of the human itch receptor MRGPRX1

MRGPRX1, a Mas-related GPCR (MRGPR), is a key receptor for itch perception and targeting MRGPRX1 may have potential to treat both chronic itch and pain. Here we report cryo-EM structures of the MRGPRX1-Gi1 and MRGPRX1-Gq trimers in complex with two peptide ligands, BAM8-22 and CNF-Tx2. These structures reveal a shallow orthosteric pocket and its conformational plasticity for sensing multiple different peptidic itch allergens. Distinct from MRGPRX2, MRGPRX1 contains a unique pocket feature at the extracellular ends of TM3 and TM4 to accommodate the peptide C-terminal "RF/RY" motif, which could serve as key mechanisms for peptidic allergen recognition. Below the ligand binding pocket, the G6.48XP6.50F6.51G6.52X(2)F/W6.55 motif is essential for the inward tilting of the upper end of TM6 to induce receptor activation. Moreover, structural features inside the ligand pocket and on the cytoplasmic side of MRGPRX1 are identified as key elements for both Gi and Gq signaling. Collectively, our studies provide structural insights into understanding itch sensation, MRGPRX1 activation, and downstream G protein signaling.

Growing ecosystem of deep learning methods for modeling protein-protein interactions

Numerous cellular functions rely on protein-protein interactions. Efforts to comprehensively characterize them remain challenged however by the diversity of molecular recognition mechanisms employed within the proteome. Deep learning has emerged as a promising approach for tackling this problem by exploiting both experimental data and basic biophysical knowledge about protein interactions. Here, we review the growing ecosystem of deep learning methods for modeling protein interactions, highlighting the diversity of these biophysically informed models and their respective trade-offs. We discuss recent successes in using representation learning to capture complex features pertinent to predicting protein interactions and interaction sites, geometric deep learning to reason over protein structures and predict complex structures, and generative modeling to design de novo protein assemblies. We also outline some of the outstanding challenges and promising new directions. Opportunities abound to discover novel interactions, elucidate their physical mechanisms, and engineer binders to modulate their functions using deep learning and, ultimately, unravel how protein interactions orchestrate complex cellular behaviors.