Structures Of Protein Complexes Research Articles

Structure and Stabilities of Solution and Gas Phase Protein Complexes.

Collision-induced unfolding (CIU) has provided new levels of understanding of the stabilities and structure(s) for gas phase protein and protein complex ions formed by electrospray ionization (ESI). Variable-temperature (vT-ESI) data provide complementary information about temperature-induced folding/unfolding (TIU) reactions of solution phase ions. Results obtained by using CIU and TIU provide complementary information about stabilities of gas phase versus solution phase ions. Such comparisons may provide the most direct experimental approach to answer a long-standing question from Fred McLafferty: "For how long, under what conditions, and to what extent, can solution structure be retained without solvent?" Answers to this question require greater understanding of the (i) structure(s), stabilities, and dynamics of proteins/protein complexes in solution prior to ESI; (ii) effects of water removal by droplet fission and "freeze-drying" by evaporation of water from the nanodroplet; and (iii) potential reactions and structural changes that may occur as the ions traverse the heated capillary, the final stage in the transition to solvent-free gas phase ions. Here, we employ vT-ESI coupled with ion mobility-mass spectrometry (IM-MS) as a means to provide more detailed answers to the above-mentioned questions. Apo- and metalated-metallothionein-2A (MT), a cysteine-rich metal binding protein, and various proteoforms of transthyretin (TTR), a homotetrameric (56 kDa) retinol and thyroxine transporter protein complex were studied to examine distinct features of CIU and TIU across two different types of protein complexes. The results in this work shed light on the capabilities of CIU, TIU, and average charge state (Zavg) for probing the rugged energy landscape of native proteins and highlights the effects of water and cofactors (metal ions) on the structure and stabilities of proteins and protein complexes.

ProAffinity-GNN: A Novel Approach to Structure-Based Protein-Protein Binding Affinity Prediction via a Curated Data Set and Graph Neural Networks.

Protein-protein interactions (PPIs) are crucial for understanding biological processes and disease mechanisms, contributing significantly to advances in protein engineering and drug discovery. The accurate determination of binding affinities, essential for decoding PPIs, faces challenges due to the substantial time and financial costs involved in experimental and theoretical methods. This situation underscores the urgent need for more effective and precise methodologies for predicting binding affinity. Despite the abundance of research on PPI modeling, the field of quantitative binding affinity prediction remains underexplored, mainly due to a lack of comprehensive data. This study seeks to address these needs by manually curating pairwise interaction labels on available 3D structures of protein complexes, with experimentally determined binding affinities, creating the largest data set for structure-based pairwise protein interaction with binding affinity to date. Subsequently, we introduce ProAffinity-GNN, a novel deep learning framework using protein language model and graph neural network (GNN) to improve the accuracy of prediction of structure-based protein-protein binding affinities. The evaluation results across several benchmark test sets and an additional case study demonstrate that ProAffinity-GNN not only outperforms existing models in terms of accuracy but also shows strong generalization capabilities.

Pharmacophore-based 3D-QSAR modeling, virtual screening, docking, molecular dynamics and biological evaluation studies for identification of potential inhibitors of alpha-glucosidase.

Alpha-glucosidase enzyme is considered an important therapeutic target for controlling hyperglycemia associated with type 2 diabetes. Novel scaffolds identified as potential alpha-glucosidase inhibitors from the Maybridge library utilizing pharmacophore modeling, molecular docking and biological evaluation are reported in this manuscript. A total of 51 xanthone series scaffolds previously reported as alpha-glucosidase inhibitors were collected and used as training and test sets. These sets were employed to develop and validate a pharmacophore-based 3D-QSAR model with statistically meaningful results using Schrodinger software. The model showed a high F value (F, 80.1) at five component partial least square factors, a high cross-validation coefficient (Q2, 0.66) and a good correlation coefficient (R2, 0.95). Pearson correlation coefficient (r) of 0.8400 indicated a greater degree of confidence in the model. Subsequently, virtual screening was performed with PHASE module of Schrodinger software using the above model to identify novel alpha-glucosidase inhibitors, and mapped compounds were evaluated for their interactions with the protein. The X-ray co-crystallised structure of the alpha-glucosidase protein in complex with acarbose (PDB Code: 5NN8) was used for molecular docking analysis using GLIDE module and a total of eight compounds were further selected for biological evaluation. Molecular dynamics analysis using GROMACS software was performed in the active site of alpha-glucosidase protein to gain insights into binding mechanism of the four active compounds which were finally found to exhibit inhibitory activity in the biological assay.

Vesicle Picker: A tool for efficient identification of membrane protein complexes in vesicles

Electron cryomicroscopy (cryo-EM) has recently allowed determination of near-atomic resolution structures of membrane proteins and protein complexes embedded in lipid vesicles. However, particle selection from electron micrographs of these vesicles can be challenging due to the strong signal contributed from the lipid bilayer. This challenge often requires iterative and laborious particle selection workflows to generate a dataset of high-quality particle images for subsequent analysis. Here we present Vesicle Picker, an open-source program built on the Segment Anything model. Vesicle Picker enables automatic identification of vesicles in cryo-EM micrographs with high recall and precision. It then exhaustively selects all potential particle locations, either at the perimeter or uniformly over the surface of the projection of the vesicle. The program is designed to interface with cryoSPARC, which performs both upstream micrograph processing and downstream single particle image analysis. We demonstrate Vesicle Picker’s utility by determining a high-resolution map of the vacuolar-type ATPase from micrographs of native synaptic vesicles (SVs) and identifying an additional protein or protein complex in the SV membrane.

A new discrete-geometry approach for integrative docking of proteins using chemical crosslinks.

The structures of protein complexes allow us to understand and modulate the biological functions of the proteins. Integrative docking is a computational method to obtain the structures of a protein complex, given the atomic structures of the constituent proteins along with other experimental data on the complex, such as chemical crosslinks or SAXS profiles. Here, we develop a new discrete geometry-based method, wall-EASAL, for integrative rigid docking of protein pairs given the structures of the constituent proteins and chemical crosslinks. The method is an adaptation of EASAL (Efficient Atlasing and Search of Assembly Landscapes), a state-of-the-art discrete geometry method for efficient and exhaustive sampling of macromolecular configurations under pairwise inter-molecular distance constraints. We provide a mathematical proof that the method finds a structure satisfying the crosslink constraints under a natural condition satisfied by energy landscapes. We compare wall-EASAL with IMP (Integrative Modeling Platform), a commonly used integrative modeling method, on a benchmark, varying the numbers, types, and sources of input crosslinks, and sources of monomer structures. The wall-EASAL method performs better than IMP in terms of the average satisfaction of the configurations to the input crosslinks and the average similarity of the configurations to their corresponding native structures. The ensembles from IMP exhibit greater variability in these two measures. Further, wall-EASAL is more efficient than IMP. Although the current study uses crosslinks, the method is general and any source of distance constraints can be used for integrative docking with wall-EASAL. However, the current implementation only supports binary rigid protein docking, i.e., assumes that the monomer structures are known and remain rigid. Additionally, the current implementation is deterministic, i.e., it does not account for uncertainties in the crosslinking data beyond using distance bounds. Neither of these appears to be a theoretical or algorithmic limitation of the EASAL methodology. Structures from wall-EASAL can be incorporated in methods for modeling large macromolecular assemblies, for example by suggesting rigid bodies or restraints for use in these methods. This will facilitate the characterization of assemblies and cellular neighborhoods at increased efficiency, accuracy, and precision. The wall-EASAL method is available at https://bitbucket.org/geoplexity/easal-dev/src/Crosslink and the benchmark is available at https://github.com/isblab/Integrative_docking_benchmark.

A new framework for SubtiWiki, the database for the model organism Bacillus subtilis.

Bacillus subtilis is a Gram-positive model bacterium and one of the most-studied and best understood organisms. The complex information resulting from its investigation is compiled in the database SubtiWiki (https://subtiwiki.uni-goettingen.de/v5) in an integrated and intuitive manner. To enhance the utility of SubtiWiki, we have added novel features such as a viewer to interrogate conserved genomic organization, a widget that shows mutant fitness data for all non-essential genes, and a widget showing protein structures, structure predictions and complex structures. Moreover, we have integrated metabolites as new entities. The new framework also includes a documented API, enabling programmatic access to data for computational tasks. Here we present the recent developments of SubtiWiki and the current state of the data for this organism.

DiffModeler: large macromolecular structure modeling for cryo-EM maps using a diffusion model.

Cryogenic electron microscopy (cryo-EM) has now been widely used for determining multichain protein complexes. However, modeling a large complex structure, such as those with more than ten chains, is challenging, particularly when the map resolution decreases. Here we present DiffModeler, a fully automated method for modeling large protein complex structures. DiffModeler employs a diffusion model for backbone tracing and integrates AlphaFold2-predicted single-chain structures for structure fitting. DiffModeler showed an average template modeling score of 0.88 and 0.91 for two datasets of cryo-EM maps of 0-5 Å resolution and 0.92 for intermediate resolution maps (5-10 Å), substantially outperforming existing methodologies. Further benchmarking at low resolutions (10-20 Å) confirms its versatility, demonstrating plausible performance.

Protein complex structure modeling by cross-modal alignment between cryo-EM maps and protein sequences

Cryo-electron microscopy (cryo-EM) technique is widely used for protein structure determination. Current automatic cryo-EM protein complex modeling methods mostly rely on prior chain separation. However, chain separation without sequence guidance often suffers from errors caused by cross-chain interaction or noise densities, which would accumulate and mislead the subsequent steps. Here, we present EModelX, a fully automated cryo-EM protein complex structure modeling method, which achieves sequence-guiding modeling through cross-modal alignments between cryo-EM maps and protein sequences. EModelX first employs multi-task deep learning to predict Cα atoms, backbone atoms, and amino acid types from cryo-EM maps, which is subsequently used to sample Cα traces with amino acid profiles. The profiles are then aligned with protein sequences to obtain initial structural models, which yielded an average RMSD of 1.17 Å in our test set, approaching atomic-level precision in recovering PDB-deposited structures. After filling unmodeled gaps through sequence-guiding Cα threading, the final models achieved an average TM-score of 0.808, outperforming the state-of-the-art method. The further combination with AlphaFold can improve the average TM-score to 0.911. Analyzes conducted by comparing some EModelX-built models and PDB structures highlight its potential to improve PDB structures. EModelX is accessible at https://bio-web1.nscc-gz.cn/app/EModelX.

De novo protein design with a denoising diffusion network independent of pretrained structure prediction models.

The recent success of RFdiffusion, a method for protein structure design with a denoising diffusion probabilistic model, has relied on fine-tuning the RoseTTAFold structure prediction network for protein backbone denoising. Here, we introduce SCUBA-diffusion (SCUBA-D), a protein backbone denoising diffusion probabilistic model freshly trained by considering co-diffusion of sequence representation to enhance model regularization and adversarial losses to minimize data-out-of-distribution errors. While matching the performance of the pretrained RoseTTAFold-based RFdiffusion in generating experimentally realizable protein structures, SCUBA-D readily generates protein structures with not-yet-observed overall folds that are different from those predictable with RoseTTAFold. The accuracy of SCUBA-D was confirmed by the X-ray structures of 16 designed proteins and a protein complex, and by experiments validating designed heme-binding proteins and Ras-binding proteins. Our work shows that deep generative models of images or texts can be fruitfully extended to complex physical objects like protein structures by addressing outstanding issues such as the data-out-of-distribution errors.

BEMM-GEN: A Toolkit for Generating a Biomolecular Environment-Mimicking Model for Molecular Dynamics Simulation.

Understanding the influence of the cellular environment on protein conformations is crucial for elucidating protein functions within living cells. In studies using molecular dynamics (MD) simulation, carbon nanotubes and hydrophobic cages have been widely used to emulate the cellular environment inside specific large biomolecules such as ribosome tunnels and chaperones. However, recent studies suggest that these uniform hydrophobic models may not adequately capture the environmental effects inside each biomolecule. Based on these facts, it is necessary to generate spherical and cylindrical models with varied chemical properties corresponding to the components within target biomolecules. We developed a biomolecular environment-mimicking model generator (BEMM-GEN) that generates spherical and cylindrical models with user-specified chemical properties and allows the integration of arbitrary protein conformations into the generated models. BEMM-GEN provides model and protein complex structures, along with the corresponding parameter files for MD simulation (AMBER and GROMACS), and users immediately run their MD simulation based on the generated input files. BEMM-GEN can be freely downloaded and installed via a Python package manager (pip install BEMM-gen). The source code files and a user manual for operation are provided on GitHub (https://github.com/y4suda/BEMM-GEN).

Molecular Mechanisms of Methamphetamine-Induced Addiction via TAAR1 Activation.

Trace amine-associated receptor 1 (TAAR1), a member of the trace amine receptor family, recognizes various trace amines in the brain, including endogenous β-phenylethylamine (PEA) and methamphetamine (METH). TAAR1 is a novel target for several neurological disorders, including schizophrenia, depression, and substance abuse. Herein, we report the structure of the human TAAR1-Gs protein complex bound to METH. Using functional studies, we reveal the molecular basis of METH recognition by TAAR1, and potential mechanisms underlying the selectivity of TAAR1 for different ligands. Molecular dynamics simulations further elucidated possible mechanisms for the binding of chiral amphetamine (AMPH)-like psychoactive drugs to TAAR1. Additionally, we discovered a hydrophobic core on the transmembrane helices (TM), TM5 and TM6, explaining the unique mechanism of TAAR1 activation. These findings reveal the ligand recognition pattern and activation mechanism of TAAR1, which has important implications for the development of next-generation treatments for substance abuse and various neurological disorders.

Advances and Challenges in Scoring Functions for RNA–Protein Complex Structure Prediction

RNA–protein complexes play a crucial role in cellular functions, providing insights into cellular mechanisms and potential therapeutic targets. However, experimental determination of these complex structures is often time-consuming and resource-intensive, and it rarely yields high-resolution data. Many computational approaches have been developed to predict RNA–protein complex structures in recent years. Despite these advances, achieving accurate and high-resolution predictions remains a formidable challenge, primarily due to the limitations inherent in current RNA–protein scoring functions. These scoring functions are critical tools for evaluating and interpreting RNA–protein interactions. This review comprehensively explores the latest advancements in scoring functions for RNA–protein docking, delving into the fundamental principles underlying various approaches, including coarse-grained knowledge-based, all-atom knowledge-based, and machine-learning-based methods. We critically evaluate the strengths and limitations of existing scoring functions, providing a detailed performance assessment. Considering the significant progress demonstrated by machine learning techniques, we discuss emerging trends and propose future research directions to enhance the accuracy and efficiency of scoring functions in RNA–protein complex prediction. We aim to inspire the development of more sophisticated and reliable computational tools in this rapidly evolving field.

Overlapping Structure Detection in Protein-Protein Interaction Networks Using a Modified Version of Particle Swarm Optimization

In today's world, the science of bioinformatics is developing rapidly, especially with regard to the analysis and study of biological networks. Scientists have used various nature-inspired algorithms to find protein complexes in protein-protein interaction (PPI) networks. These networks help scientists guess the molecular function of unknown proteins and show how cells work regularly. It is very common in PPI networks for a protein to participate in multiple functions and belong to many complexes, and as a result, complexes may overlap in the PPI networks. However, developing an efficient and reliable method to address the problem of detecting overlapping protein complexes remains a challenge since it is considered a complex and hard optimization problem. One of the main difficulties in identifying overlapping protein complexes is the accuracy of the partitioning results. In order to accurately identify the overlapping structure of protein complexes, this paper has proposed an overlapping complex detection algorithm termed OCDPSO-Net, which is based on PSO-Net (a well-known modified version of the particle swarm optimization algorithm). The framework of the OCDPSO-Net method consists of three main steps, including an initialization strategy, a movement strategy for each particle, and enhancing search ability in order to expand the solution space. The proposed algorithm has employed the partition density concept for measuring the partitioning quality in PPI network complexes and tried to optimize the value of this quantity by applying the line graph concept of the original graph representing the protein interaction network. The OCDPSO-Net algorithm is applied to a Collins PPI network and the obtained results are compared with different state-of-the-art algorithms in terms of precision ( ), recall ( ), and F-measure ( ). Experimental results confirm that the proposed algorithm has good clustering performance and has outperformed most of the existing recent overlapping algorithms. .

RAB12-LRRK2 complex suppresses primary ciliogenesis and regulates centrosome homeostasis in astrocytes

The leucine-rich repeat kinase 2 (LRRK2) phosphorylates a subset of RAB GTPases, and their phosphorylation levels are elevated by Parkinson’s disease (PD)-linked mutations of LRRK2. However, the precise function of the LRRK2-regulated RAB GTPase in the brain remains to be elucidated. Here, we identify RAB12 as a robust LRRK2 substrate in the mouse brain through phosphoproteomics profiling and solve the structure of RAB12-LRRK2 protein complex through Cryo-EM analysis. Mechanistically, RAB12 cooperates with LRRK2 to inhibit primary ciliogenesis and regulate centrosome homeostasis in astrocytes through enhancing the phosphorylation of RAB10 and recruiting RILPL1, while the functions of RAB12 require a direct interaction with LRRK2 and LRRK2 activity. Furthermore, the ciliary and centrosome defects caused by the PD-linked LRRK2-G2019S mutation are prevented by Rab12 deletion in astrocytes. Thus, our study reveals a physiological function of the RAB12-LRRK2 complex in regulating ciliogenesis and centrosome homeostasis. The RAB12-LRRK2 structure offers a guidance in the therapeutic development of PD by targeting the RAB12-LRRK2 interaction.

PCP-GC-LM: single-sequence-based protein contact prediction using dual graph convolutional neural network and convolutional neural network

BackgroundRecently, the process of evolution information and the deep learning network has promoted the improvement of protein contact prediction methods. Nevertheless, still remain some bottleneck: (1) One of the bottlenecks is the prediction of orphans and other fewer evolution information proteins. (2) The other bottleneck is the method of predicting single-sequence-based proteins mainly focuses on selecting protein sequence features and tuning the neural network architecture, However, while the deeper neural networks improve prediction accuracy, there is still the problem of increasing the computational burden. Compared with other neural networks in the field of protein prediction, the graph neural network has the following advantages: due to the advantage of revealing the topology structure via graph neural network and being able to take advantage of the hierarchical structure and local connectivity of graph neural networks has certain advantages in capturing the features of different levels of abstraction in protein molecules. When using protein sequence and structure information for joint training, the dependencies between the two kinds of information can be better captured. And it can process protein molecular structures of different lengths and shapes, while traditional neural networks need to convert proteins into fixed-size vectors or matrices for processing.ResultsHere, we propose a single-sequence-based protein contact map predictor PCP-GC-LM, with dual-level graph neural networks and convolution networks. Our method performs better with other single-sequence-based predictors in different independent tests. In addition, to verify the validity of our method against complex protein structures, we will also compare it with other methods in two homodimers protein test sets (DeepHomo test dataset and CASP-CAPRI target dataset). Furthermore, we also perform ablation experiments to demonstrate the necessity of a dual graph network. In all, our framework presents new modules to accurately predict inter-chain contact maps in protein and it’s also useful to analyze interactions in other types of protein complexes.

Prognostic importance of splicing-triggeredaberrations of protein complex interfaces in cancer.

Aberrant alternative splicing (AS) is a prominent hallmark of cancer. AS can perturb protein-protein interactions (PPIs) by adding or removing interface regions encoded by individual exons. Identifying prognostic exon-exon interactions (EEIs) from PPI interfaces can help discover AS-affected cancer-driving PPIs that can serve as potential drug targets. Here, we assessed the prognostic significance of EEIs across 15 cancer types by integrating RNA-seq data with three-dimensional (3D) structures of protein complexes. By analyzing the resulting EEI network we identified patient-specific perturbed EEIs (i.e.,EEIs present in healthy samples but absent from the paired cancer samples or vice versa) that were significantly associated with survival. We provide the first evidence that EEIs can be used as prognostic biomarkers for cancer patient survival. Our findings provide mechanistic insights into AS-affected PPI interfaces. Given the ongoing expansion of available RNA-seq data and the number of 3D structurally-resolved (or confidently predicted) protein complexes, our computational framework will help accelerate the discovery of clinically important cancer-promoting AS events.

Cryo-EM Structure of the Mnx Protein Complex Reveals a Tunnel Framework for the Mechanism of Manganese Biomineralization.

The global manganese cycle relies on microbes to oxidize soluble Mn(II) to insoluble Mn(IV) oxides. Some microbes require peroxide or superoxide as oxidants, but others can use O2 directly, via multicopper oxidase (MCO) enzymes. One of these, MnxG from Bacillus sp. strain PL-12, was isolated in tight association with small accessory proteins, MnxE and MnxF. The protein complex, called Mnx, has eluded crystallization efforts, but we now report the 3D structure of a point mutant using cryo-EM single particle analysis, cross-linking mass spectrometry, and AlphaFold Multimer prediction. The β-sheet-rich complex features MnxG enzyme, capped by a heterohexameric ring of alternating MnxE and MnxF subunits, and a tunnel that runs through MnxG and its MnxE3F3 cap. The tunnel dimensions and charges can accommodate the mechanistically inferred binuclear manganese intermediates. Comparison with the Fe(II)-oxidizing MCO, ceruloplasmin, identifies likely coordinating groups for the Mn(II) substrate, at the entrance to the tunnel. Thus, the 3D structure provides a rationale for the established manganese oxidase mechanism, and a platform for further experiments to elucidate mechanistic details of manganese biomineralization.

Deciphering 3'UTR Mediated Gene Regulation Using Interpretable Deep Representation Learning.

The 3' untranslated regions (3'UTRs) of messenger RNAs contain many important cis-regulatory elements that are under functional and evolutionary constraints. It is hypothesized that these constraints are similar to grammars and syntaxes in human languages and can be modeled by advanced natural language techniques such as Transformers, which has been very effective in modeling complex protein sequence and structures. Here 3UTRBERT is described, which implements an attention-based language model, i.e., Bidirectional Encoder Representations from Transformers (BERT). 3UTRBERT is pre-trained on aggregated 3'UTR sequences of human mRNAs in a task-agnostic manner; the pre-trained model is then fine-tuned for specific downstream tasks such as identifying RBP binding sites, m6A RNA modification sites, and predicting RNA sub-cellular localizations. Benchmark results show that 3UTRBERT generally outperformed other contemporary methods in each of these tasks. More importantly, the self-attention mechanism within 3UTRBERT allows direct visualization of the semantic relationship between sequence elements and effectively identifies regions with important regulatory potential. It is expected that 3UTRBERT model can serve as the foundational tool to analyze various sequence labeling tasks within the 3'UTR fields, thus enhancing the decipherability of post-transcriptional regulatorymechanisms.

Generation of a high confidence set of domain-domain interface types to guide protein complex structure predictions by AlphaFold.

While the release of AlphaFold (AF) represented a breakthrough for the prediction of protein complex structures, its sensitivity, especially when using full length protein sequences, still remains limited. Modeling success rates might increase if AF predictions were guided by likely interacting protein fragments. This approach requires available sets of highly confident protein-protein interface types. Computational resources, such as 3did, infer interacting globular domain types from observed contacts in protein structures. Assessing the accuracy of these predicted interface types is difficult because we lack hand-curated reference sets of verified domain-domain interface (DDI) types. To improve protein complex modeling of DDIs by AF, we manually inspected 80 randomly selected DDI types from the 3did resource to generate a first reference set of DDI types. Identified cases of DDI type nonapproval (40%) primarily resulted from inaccurate Pfam domain matches, crystal contacts, and synthetic protein constructs. Using logistic regression, we predicted a subset of 2411 out of 5724 considered DDI types in 3did to be of high confidence, which we subsequently applied to 53000 human-protein interactions to predict DDIs followed by AF modeling. We obtained highly confident AF models for 604 out of 1129 predicted DDIs. Of note, for 47% of them no confident AF structural model could be obtained using full length protein sequences. Code is available at https://github.com/KatjaLuckLab/DDI_manuscript.

MAGPIE: An interactive tool for visualizing and analyzing protein-ligand interactions.

Quantitative tools to compile and analyze biomolecular interactions among chemically diverse binding partners would improve therapeutic design and aid in studying molecular evolution. Here we present Mapping Areas of Genetic Parsimony In Epitopes (MAGPIE), a publicly available software package for simultaneously visualizing and analyzing thousands of interactions between a single protein or small molecule ligand (the "target") and all of its protein binding partners ("binders"). MAGPIE generates an interactive three-dimensional visualization from a set of protein complex structures that share the target ligand, as well as sequence logo-style amino acid frequency graphs that show all the amino acids from the set of protein binders that interact with user-defined target ligand positions or chemical groups. MAGPIE highlights all the salt bridge and hydrogen bond interactions made by the target in the visualization and as separate amino acid frequency graphs. Finally, MAGPIE collates the most common target-binder interactions as a list of "hotspots," which can be used to analyze trends or guide the de novo design of protein binders. As an example of the utility of the program, we used MAGPIE to probe how different antibody fragments bind a viral antigen; how a common metabolite binds diverse protein partners; and how two ligands bind orthologs of a well-conserved glycolytic enzyme for a detailed understanding of evolutionarily conserved interactions involved in its activation and inhibition. MAGPIE is implemented in Python 3 and freely available at https://github.com/glasgowlab/MAGPIE, along with sample datasets, usage examples, and helper scripts to prepare input structures.