Challenges in protein docking

Highly Optimized Simulation of Atomic Resolution Cell-Like Protein Environment

Computational approaches can provide details of molecular mechanisms in a crowded environment inside cells. Protein docking predicts stable configurations of molecular complexes, which correspond to deep energy minima. Systematic docking approaches, such as those based on fast Fourier transform (FFT), also map the entire intermolecular energy landscape by determining the position and depth of the full spectrum of the energy minima. Such mapping allows speeding up simulations by precalculating the intermolecular energy values. Our earlier study combined FFT docking with the Monte Carlo protocol, enabling simulation of cell-size, crowded protein systems with seconds, and longer trajectories at atomic resolution, several orders of magnitude longer than those achievable by alternative approaches. In this study, we present a further drastic extension of the modeling capabilities by parallelized implementation of the simulation protocol. The procedure was applied to a panel of Death Fold Domains that form nucleated polymers in human innate immune signaling, recapitulating their homooligomerization tendencies and providing insights into the molecular mechanisms of polymer nucleation. The parallelized protocol allows extension of the simulation trajectories by orders of magnitude beyond the previously reported implementation, reaching into the uncharted territory of atomic resolution simulation of cell-sized systems.

Dockground: The resource expands to protein-RNA interactome

RNA is a master regulator of cellular processes and will bind to many different proteins throughout its life cycle. Dysregulation of RNA and RNA-binding proteins can lead to various diseases, including cancer. To better understand molecular mechanisms of the cellular processes, it is important to characterize protein-RNA interactions at the structural level. There is a lack of experimental structures available for protein-RNA complexes due to the RNA inherent flexibility, which complicates the experimental structure determination. The scarcity of structures can be made up for with computational modeling. Dockground is a resource for development and benchmarking of structure-based modeling of protein interactions. It contains datasets focusing on different aspects of protein recognition. The foundation of all the datasets is the database of experimentally determined protein complexes, which previously contained only protein-protein assemblies. To further expand the utility of the Dockground resource, we extended the database to protein-RNA interactions. The new functionalities are available on the Dockground website at https://dockground.compbio.ku.edu/. The database can be searched using a number of criteria, including removal of redundancies at various sequence and structure similarity thresholds. The database updates with new structures from the Protein Data Bank on a weekly basis.

AlphaFold and Docking Approaches for Antibody-Antigen and Other Targets: Insights From CAPRI Rounds 47-55

Accurate modeling of the structures of protein-protein complexes and other biomolecular interactions represents a longstanding and important challenge for computational biology. The Critical Assessment of PRedicted Interactions (CAPRI) experiment has served for over two decades as a key means to assess and compare current approaches and methods through blind predictive scenarios, highlighting useful strategies, and new developments. Here we describe the performance of our laboratory's team in recent CAPRI rounds, which included submissions for 10 modeling targets. Our team utilized a range of docking and modeling approaches, including ZDOCK, Rosetta, and ZRANK, to model, refine, and score protein-protein and protein-DNA complexes. For recent targets we utilized adaptations of AlphaFold to generate models, leading to near-native models for an antibody-peptide target, and a highly accurate (but low ranked) model for an antibody-MHC complex. These results underscore the utility of AlphaFold-based protocols for predictive protein complex modeling, including for immune recognition, and highlight considerations regarding the use of AlphaFold confidence metrics in model selection.

Steered molecular dynamics simulation as a post-process to optimize the iBRAB-designed Fab model

Therapeutic monoclonal antibodies are an effective method of treating acute infectious diseases. However, knowing which of the produced antibodies in the vast number of human antibodies can cure the disease requires a long time and advanced technology. The previously introduced iBRAB method relies on studied antibodies to design a broad-spectrum antibody capable of neutralizing antigens of many different Influenza A viral strains. To evaluate the antigen-binding fragment as an applicable drug, the therapeutic antibody profiles providing guidelines collected from clinically staged therapeutic antibodies were used to access different measurements. Although the evaluated values were within an accepted range, the modification in the amino acid sequence is required for better properties. Thus, using the steered molecular dynamics (SMD) simulation to determine the binding capacity of amino acids in the functional region, the profile of interacted amino acids of Fab with the antigen was established for modified reference. As a result, the model was modified with amino acids elimination at positions 96-97 in the heavy chain and 26-27, 91, 96-97, and 102-103 in the light chain, which has better Therapeutic Antibody Profiler evaluations than the original designation. Thus again, SMD simulation is a promising computational approach for post-modification in rational drug design.

The interaction of PRDX1 with Cofilin promotes oral squamous cell carcinoma metastasis

Peroxiredoxin 1 (PRDX1) is an important member of the peroxiredoxin family (PRDX) and is upregulated in a variety of tumors. Previous studies have found that high PRDX1 expression is closely related to the metastasis of oral squamous cell carcinoma (OSCC), but the specific molecular mechanism is elusive. To elucidate the role of PRDX1 in the metastasis process of OSCC, we evaluated the expression of PRDX1 in OSCC clinical specimens and its impact on the prognosis of OSCC patients. Then, the effect of PRDX1 on OSCC metastasis and cytoskeletal reconstruction was explored in vitro and in nude mouse tongue cancer models, and the molecular mechanisms were also investigated. PRDX1 can directly interact with the actin-binding protein Cofilin, inhibiting the phosphorylation of its Ser3 site, accelerating the depolymerization and turnover of actin, promoting OSCC cell movement, and aggravating the invasion and metastasis of OSCC. In clinical samples and mouse tongue cancer models, PRDX1 also increased lymph node metastasis of OSCC and was negatively correlated with the phosphorylation of Cofilin; PRDX1 also reduced the overall survival rate of OSCC patients. In summary, our study identified that PRDX1 may be a potential therapeutic target to inhibit OSCC metastasis.

Diffusion of proteins in crowded solutions studied by docking-based modeling

The diffusion of proteins is significantly affected by macromolecular crowding. Molecular simulations accounting for protein interactions at atomic resolution are useful for characterizing the diffusion patterns in crowded environments. We present a comprehensive analysis of protein diffusion under different crowding conditions based on our recent docking-based approach simulating an intracellular crowded environment by sampling the intermolecular energy landscape using the Markov Chain Monte Carlo protocol. The procedure was extensively benchmarked, and the results are in very good agreement with the available experimental and theoretical data. The translational and rotational diffusion rates were determined for different types of proteins under crowding conditions in a broad range of concentrations. A protein system representing most abundant protein types in the E. coli cytoplasm was simulated, as well as large systems of other proteins of varying sizes in heterogeneous and self-crowding solutions. Dynamics of individual proteins was analyzed as a function of concentration and different diffusion rates in homogeneous and heterogeneous crowding. Smaller proteins diffused faster in heterogeneous crowding of larger molecules, compared to their diffusion in the self-crowded solution. Larger proteins displayed the opposite behavior, diffusing faster in the self-crowded solution. The results show the predictive power of our structure-based simulation approach for long timescales of cell-size systems at atomic resolution.

In-vitro and in-silico analyses of the thrombolytic potential of green kiwifruit

Cardiovascular diseases (CVDs), mainly caused by thrombosis complications, are the leading cause of mortality worldwide, making the development of alternative treatments highly desirable. In this study, the thrombolytic potential of green kiwifruit (Actinidia deliciosa cultivar Hayward) was assessed using in-vitro and in-silico approaches. The crude green kiwifruit extract demonstrated the ability to reduce blood clots significantly by 73.0 ± 1.12% (P < 0.01) within 6 h, with rapid degradation of Aα and Bβ fibrin chains followed by the γ chain in fibrinolytic assays. Molecular docking revealed six favorable conformations for the kiwifruit enzyme actinidin (ADHact) and fibrin chains, supported by spontaneous binding energies and distances. Moreover, molecular dynamics simulation confirmed the binding stability of the complexes of these conformations, as indicated by the stable binding affinity, high number of hydrogen bonds, and consistent distances between the catalytic residue Cys25 of ADHact and the peptide bond. The better overall binding affinity of ADHact to fibrin chains Aα and Bβ may contribute to their faster degradation, supporting the fibrinolytic results. In conclusion, this study demonstrated the thrombolytic potential of the green kiwifruit-derived enzyme and highlighted its potential role as a natural plant-based prophylactic and therapeutic agent for CVDs.

In silico models of the macromolecular NaV1.5-KIR2.1 complex

In cardiac cells, the expression of the cardiac voltage-gated Na+ channel (NaV1.5) is reciprocally regulated with the inward rectifying K+ channel (KIR2.1). These channels can form macromolecular complexes that pre-assemble early during forward trafficking (transport to the cell membrane). In this study, we present in silico 3D models of NaV1.5-KIR2.1, generated by rigid-body protein-protein docking programs and deep learning-based AlphaFold-Multimer software. Modeling revealed that the two channels could physically interact with each other along the entire transmembrane region. Structural mapping of disease-associated mutations revealed a hotspot at this interface with several trafficking-deficient variants in close proximity. Thus, examining the role of disease-causing variants is important not only in isolated channels but also in the context of macromolecular complexes. These findings may contribute to a better understanding of the life-threatening cardiovascular diseases underlying KIR2.1 and NaV1.5 malfunctions.

Approaches to Backbone Flexibility in Protein-Protein Docking

Proteins are the fundamental organic macromolecules in living systems that play a key role in a variety of biological functions including immunological detection, intracellular trafficking, and signal transduction. The docking of proteins has greatly advanced during recent decades and has become a crucial complement to experimental methods. Protein-protein docking is a helpful method for simulating protein complexes whose structures have not yet been solved experimentally. This chapter focuses on major search tactics along with various docking programs used in protein-protein docking algorithms, which include: direct search, exhaustive global search, local shape feature matching, randomized search, and broad category of post-docking approaches. As backbone flexibility predictions and interactions in high-resolution protein-protein docking remain important issues in the overall optimization context, we have put forward several methods and solutions used to handle backbone flexibility. In addition, various docking methods that are utilized for flexible backbone docking, including ATTRACT, FlexDock, FLIPDock, HADDOCK, RosettaDock, FiberDock, etc., along with their scoring functions, algorithms, advantages, and limitations are discussed. Moreover, what progress in search technology is expected, including not only the creation of new search algorithms but also the enhancement of existing ones, has been debated. As conformational flexibility is one of the most crucial factors affecting docking success, more work should be put into evaluating the conformational flexibility upon binding for a particular case in addition to developing new algorithms to replace the rigid body docking and scoring approach.

Mechanisms of Nelumbinis folium targeting PPARγ for weight management: A molecular docking and molecular dynamics simulations study

The lotus leaf, Nelumbinis folium (NF), has frequently appeared in obesity clinical trials as an intervention to promote weight loss and improve metabolic profiles. However, the molecular mechanisms by which it interacts with important obesity targets and pathways, such as the peroxisome proliferator-activated receptor gamma (PPARγ) within the PPAR signalling pathway, were not well understood. This study aims to screen for candidate compounds from NF with desirable pharmacokinetic properties and examine their binding feasibility at the PPARγ ligand-binding domain (LBD). Ligand- and structure-based screening of NF compounds were performed, and a consensus approach has been applied to identify druggable candidates. By examining the pharmacokinetic profiles, a large proportion of NF compounds exhibited favourable drug-likeness and oral bioavailability properties. Furthermore, the binding affinity scores and poses provided new insights on the distinctive binding behaviours of NF compounds at the LBD of PPARγ in its inactive form. Several NF compounds could bind strongly to PPARγ at sub-pockets where partial agonists and antagonists were found to bind and may induce conformational changes that influence co-repressor binding, trans-repression, and gene expression inhibition. Subsequent molecular dynamics simulations of a candidate compound (NF129 narcissin) bound to PPARγ revealed conformational stability, residue fluctuation, and binding behaviours comparable to that of the known inhibitor, SR1664. Therefore, it can be proposed that narcissin exhibits characteristics of a PPARγ antagonist. Further experimental validation to support the development of NF129 as a future anti-obesity agent is warranted.

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

Challenges in structural modeling of RNA-protein interactions

In the past few years, the number of RNA-binding proteins (RBP) and RNA-RBP interactions has increased significantly. Here, we review recent developments in the methodology for protein-RNA and protein-protein complex structure modeling with deep learning and co-evolution, as well as discuss the challenges and opportunities for building a reliable approach for protein-RNA complex structure modelling. Protein Data bank (PDB) and Cross-linking immunoprecipitation (CLIP) data could be combined together and used to infer 2D geometry of protein-RNA interactions by deep learning.

The CDR3 region as the major driver of TREM-1 interaction with its ligands, an in silico characterization

The triggering receptor expressed on myeloid cells-1 (TREM-1) is a pattern recognition receptor heavily investigated in infectious and non-infectious diseases. Because of its role in amplifying inflammation, TREM-1 has been explored as a diagnostic/prognostic biomarker. Further, as the receptor has been implicated in the pathophysiology of several diseases, therapies aiming at modulating its activity represent a promising strategy to constrain uncontrolled inflammatory or infectious diseases. Despite this, several aspects concerning its interaction with ligands and activation process, remain unclear. Although many molecules have been suggested as TREM-1 ligands, only five have been confirmed to interact with the receptor: actin, eCIRP, HMGB1, Hsp70 and PGLYRP1. However, the domains involved in the interaction between the receptor and these proteins are not clarified yet. Therefore, here we used in silico approaches to investigate the putative binding domains in the receptor, using hot spots analysis, molecular docking and molecular dynamics simulations between TREM-1 and its five known ligands. Our results indicated the complementarity-determining regions (CDRs) of the receptor as the main mediators of antigen recognition, especially the CDR3 loop. We believe that our study could be used as structural basis for the elucidation of TREM-1's recognition process, and may be useful for prospective in silico and biological investigations exploring the receptor in different contexts.

Advancing Targeted Protein Degradation via Multiomics Profiling and Artificial Intelligence

Only around 20% of the human proteome is considered to be druggable with small-molecule antagonists. This leaves some of the most compelling therapeutic targets outside the reach of ligand discovery. The concept of targeted protein degradation (TPD) promises to overcome some of these limitations. In brief, TPD is dependent on small molecules that induce the proximity between a protein of interest (POI) and an E3 ubiquitin ligase, causing ubiquitination and degradation of the POI. In this perspective, we want to reflect on current challenges in the field, and discuss how advances in multiomics profiling, artificial intelligence, and machine learning (AI/ML) will be vital in overcoming them. The presented roadmap is discussed in the context of small-molecule degraders but is equally applicable for other emerging proximity-inducing modalities.

Compounds in Indonesian Ginger Rhizome Extracts and Their Potential for Anti-Skin Aging Based on Molecular Docking

Skin aging is a condition caused by reactive oxygen species (ROS) and advanced glycation end products (AGEs). Indonesian gingers (Zingiber officinale), which consists of Gajah (GG), Red (MM), and Emprit (EE) ginger, are thought to produce anti-skin aging compounds through enzyme inhibition. The enzymes used in the molecular docking study were collagenase, hyaluronidase, elastase, and tyrosinase. This study aimed to determine the compounds contained in Indonesian ginger rhizome ethanolic extracts using liquid chromatography–mass spectrometry/mass spectrometry to differentiate metabolites contained in the different Indonesian ginger rhizome extracts. A principal component analysis (PCA) and a heat map analysis were used in order to determine which compounds and extracts contained potential anti-skin aging properties based on a molecular docking study. Ascorbic acid was used as a control ligand in the molecular docking study. Ninety-eight compounds were identified in three different ginger rhizomes extracts and were grouped into three separate quadrants. The most potent compound for anti-skin aging in the Indonesian ginger rhizome extracts was octinoxate. Octinoxate showed a high abundance in the EE ginger rhizome extract. Therefore, the EE ginger extract was the Indonesian ginger rhizome extract with the greatest potential for anti-skin aging.

Dockground resource for protein recognition studies

Structural information of protein-protein interactions is essential for characterization of life processes at the molecular level. While a small fraction of known protein interactions has experimentally determined structures, computational modeling of protein complexes (protein docking) has to fill the gap. The Dockground resource (http://dockground.compbio.ku.edu) provides a collection of datasets for the development and testing of protein docking techniques. Currently, Dockground contains datasets for the bound and the unbound (experimentally determined and simulated) protein structures, model-model complexes, docking decoys of experimentally determined and modeled proteins, and templates for comparative docking. The Dockground bound proteins dataset is a core set, from which other Dockground datasets are generated. It is devised as a relational PostgreSQL database containing information on experimentally determined protein-protein complexes. This report on the Dockground resource describes current status of the datasets, new automated update procedures and further development of the core datasets. We also present a new Dockground interactive web interface, which allows search by various parameters, such as release date, multimeric state, complex type, structure resolution, and so on, visualization of the search results with a number of customizable parameters, as well as downloadable datasets with predefined levels of sequence and structure redundancy.

CoVM2: Molecular Biological Data Integration of SARS-CoV-2 Proteins in a Macro-to-Micro Method

The COVID-19 pandemic has been a major public health event since 2020. Multiple variant strains of SARS-CoV-2, the causative agent of COVID-19, were detected based on the mutation sites in their sequences. These sequence mutations may lead to changes in the protein structures and affect the binding states of SARS-CoV-2 and human proteins. Experimental research on SARS-CoV-2 has accumulated a large amount of structural data and protein-protein interactions (PPIs), but the studies on the SARS-CoV-2–human PPI networks lack integration of physical associations with possible protein docking information. In addition, the docking structures of variant viral proteins with human receptor proteins are still insufficient. This study constructed SARS-CoV-2–human protein–protein interaction network with data integration methods. Crystal structures were collected to map the interaction pairs. The pairs of direct interactions and physical associations were selected and analyzed for variant docking calculations. The study examined the structures of spike (S) glycoprotein of variants Delta B.1.617.2, Omicron BA.1, and Omicron BA.2. The calculated docking structures of S proteins and potential human receptors were obtained. The study integrated binary protein interactions with 3D docking structures to fulfill an extended view of SARS-CoV-2 proteins from a macro- to micro-scale.

DOCKGROUND membrane protein-protein set

Membrane proteins are significantly underrepresented in Protein Data Bank despite their essential role in cellular mechanisms and the major progress in experimental protein structure determination. Thus, computational approaches are especially valuable in the case of membrane proteins and their assemblies. The main focus in developing structure prediction techniques has been on soluble proteins, in part due to much greater availability of the structural data. Currently, structure prediction of protein complexes (protein docking) is a well-developed field of study. However, the generic protein docking approaches are not optimal for the membrane proteins because of the differences in physicochemical environment and the spatial constraints imposed by the membranes. Thus, docking of the membrane proteins requires specialized computational methods. Development and benchmarking of the membrane protein docking approaches has to be based on high-quality sets of membrane protein complexes. In this study we present a new dataset of 456 non-redundant alpha helical binary interfaces. The set is significantly larger and more representative than the previously developed sets. In the future, it will become the basis for the development of docking and scoring benchmarks, similar to the ones for soluble proteins in the Dockground resource http://dockground.compbio.ku.edu.

GWYRE: A resource for mapping variants onto experimental and modeled structures of human protein complexes

•

Structure of protein complexes is important for interpreting genetic variation.

•

Data on single amino acid variants is available from high-throughput sequencing.

•

Integrated modeling approach was applied to proteins and their complexes.

•

GWYRE resource incorporates predicted protein complexes with mapped mutations.

Rapid progress in structural modeling of proteins and their interactions is powered by advances in knowledge-based methodologies along with better understanding of physical principles of protein structure and function. The pool of structural data for modeling of proteins and protein–protein complexes is constantly increasing due to the rapid growth of protein interaction databases and Protein Data Bank. The GWYRE (Genome Wide PhYRE) project capitalizes on these developments by advancing and applying new powerful modeling methodologies to structural modeling of protein–protein interactions and genetic variation. The methods integrate knowledge-based tertiary structure prediction using Phyre2 and quaternary structure prediction using template-based docking by a full-structure alignment protocol to generate models for binary complexes. The predictions are incorporated in a comprehensive public resource for structural characterization of the human interactome and the location of human genetic variants. The GWYRE resource facilitates better understanding of principles of protein interaction and structure/function relationships. The resource is available at http://www.gwyre.org.

Oncological drug discovery: AI meets structure-based computational research

The integration of machine learning and structure-based methods has proven valuable in the past as a way to prioritize targets and compounds in early drug discovery. In oncological research, these methods can be highly beneficial in addressing the diversity of neoplastic diseases portrayed by the different hallmarks of cancer. Here, we review six use case scenarios for integrated computational methods, namely driver prediction, computational mutagenesis, (off)-target prediction, binding site prediction, virtual screening, and allosteric modulation analysis. We address the heterogeneity of integration approaches and individual methods, while acknowledging their current limitations and highlighting their potential to bring drugs for personalized oncological therapies to the market faster.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics

Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.

Challenges and frontiers of computational modelling of biomolecular recognition

Biomolecular recognition including binding of small molecules, peptides and proteins to their target receptors plays a key role in cellular function and has been targeted for therapeutic drug design. However, the high flexibility of biomolecules and slow binding and dissociation processes have presented challenges for computational modelling. Here, we review the challenges and computational approaches developed to characterise biomolecular binding, including molecular docking, molecular dynamics simulations (especially enhanced sampling) and machine learning. Further improvements are still needed in order to accurately and efficiently characterise binding structures, mechanisms, thermodynamics and kinetics of biomolecules in the future.

From complete cross-docking to partners identification and binding sites predictions

Proteins ensure their biological functions by interacting with each other. Hence, characterising protein interactions is fundamental for our understanding of the cellular machinery, and for improving medicine and bioengineering. Over the past years, a large body of experimental data has been accumulated on who interacts with whom and in what manner. However, these data are highly heterogeneous and sometimes contradictory, noisy, and biased. Ab initio methods provide a means to a "blind" protein-protein interaction network reconstruction. Here, we report on a molecular cross-docking-based approach for the identification of protein partners. The docking algorithm uses a coarse-grained representation of the protein structures and treats them as rigid bodies. We applied the approach to a few hundred of proteins, in the unbound conformations, and we systematically investigated the influence of several key ingredients, such as the size and quality of the interfaces, and the scoring function. We achieved some significant improvement compared to previous works, and a very high discriminative power on some specific functional classes. We provide a readout of the contributions of shape and physico-chemical complementarity, interface matching, and specificity, in the predictions. In addition, we assessed the ability of the approach to account for protein surface multiple usages, and we compared it with a sequence-based deep learning method. This work may contribute to guiding the exploitation of the large amounts of protein structural models now available toward the discovery of unexpected partners and their complex structure characterisation.

Assessment of the CASP14 assembly predictions

In CASP14, 39 research groups submitted more than 2500 3D models on 22 protein complexes. In general, the community performed well in predicting the fold of the assemblies (for 80% of the targets), although it faced significant challenges in reproducing the native contacts. This is especially the case for the complexes without whole-assembly templates. The leading predictor, BAKER-experimental, used a methodology combining classical techniques (template-based modeling, protein docking) with deep learning-based contact predictions and a fold-and-dock approach. The Venclovas team achieved the runner-up position with template-based modeling and docking. By analyzing the target interfaces, we showed that the complexes with depleted charged contacts or dominating hydrophobic interactions were the most challenging ones to predict. We also demonstrated that if AlphaFold2 predictions were at hand, the interface prediction challenge could be alleviated for most of the targets. All in all, it is evident that new approaches are needed for the accurate prediction of assemblies, which undoubtedly will expand on the significant improvements in the tertiary structure prediction field.

InterEvDock3: a combined template-based and free docking server with increased performance through explicit modeling of complex homologs and integration of covariation-based contact maps

The InterEvDock3 protein docking server exploits the constraints of evolution by multiple means to generate structural models of protein assemblies. The server takes as input either several sequences or 3D structures of proteins known to interact. It returns a set of 10 consensus candidate complexes, together with interface predictions to guide further experimental validation interactively. Three key novelties were implemented in InterEvDock3 to help obtain more reliable models: users can (i) generate template-based structural models of assemblies using close and remote homologs of known 3D structure, detected through an automated search protocol, (ii) select the assembly models most consistent with contact maps from external methods that implement covariation-based contact prediction with or without deep learning and (iii) exploit a novel coevolution-based scoring scheme at atomic level, which leads to significantly higher free docking success rates. The performance of the server was validated on two large free docking benchmark databases, containing respectively 230 unbound targets (Weng dataset) and 812 models of unbound targets (PPI4DOCK dataset). Its effectiveness has also been proven on a number of challenging examples. The InterEvDock3 web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock3/.

Moving pictures: Reassessing docking experiments with a dynamic view of protein interfaces

The modeling of protein assemblies at the atomic level remains a central issue in structural biology, as protein interactions play a key role in numerous cellular processes. This problem is traditionally addressed using docking tools, where the quality of the models is based on their similarity to a single reference experimental structure. However, using a static reference does not take into account the dynamic quality of the protein interface. Here, we used all-atom classical Molecular Dynamics simulations to investigate the stability of the reference interface for three complexes that previously served as targets in the CAPRI competition. For each one of these targets, we also ran MD simulations for ten models that are distributed over the High, Medium and Acceptable accuracy categories. To assess the quality of these models from a dynamic perspective, we set up new criteria which take into account the stability of the reference experimental protein interface. We show that, when the protein interfaces are allowed to evolve along time, the original ranking based on the static CAPRI criteria no longer holds as over 50% of the docking models undergo a category change (which can be either toward a better or a lower accuracy group) when reassessing their quality using dynamic information.

Tools for the Recognition of Sorting Signals and the Prediction of Subcellular Localization of Proteins From Their Amino Acid Sequences

At the time of translation, nascent proteins are thought to be sorted into their final subcellular localization sites, based on the part of their amino acid sequences (i.e., sorting or targeting signals). Thus, it is interesting to computationally recognize these signals from the amino acid sequences of any given proteins and to predict their final subcellular localization with such information, supplemented with additional information (e.g., k-mer frequency). This field has a long history and many prediction tools have been released. Even in this era of proteomic atlas at the single-cell level, researchers continue to develop new algorithms, aiming at accessing the impact of disease-causing mutations/cell type-specific alternative splicing, for example. In this article, we overview the entire field and discuss its future direction.