Flexible protein–protein docking

Protein language model-embedded geometric graphs power inter-protein contact prediction

Accurate prediction of contacting residue pairs between interacting proteins is very useful for structural characterization of protein-protein interactions. Although significant improvement has been made in inter-protein contact prediction recently, there is still a large room for improving the prediction accuracy. Here we present a new deep learning method referred to as PLMGraph-Inter for inter-protein contact prediction. Specifically, we employ rotationally and translationally invariant geometric graphs obtained from structures of interacting proteins to integrate multiple protein language models, which are successively transformed by graph encoders formed by geometric vector perceptrons and residual networks formed by dimensional hybrid residual blocks to predict inter-protein contacts. Extensive evaluation on multiple test sets illustrates that PLMGraph-Inter outperforms five top inter-protein contact prediction methods, including DeepHomo, GLINTER, CDPred, DeepHomo2, and DRN-1D2D_Inter, by large margins. In addition, we also show that the prediction of PLMGraph-Inter can complement the result of AlphaFold-Multimer. Finally, we show leveraging the contacts predicted by PLMGraph-Inter as constraints for protein-protein docking can dramatically improve its performance for protein complex structure prediction.

Identification, characterization, and CADD analysis of Plasmodium DMAP1 reveals it as a potential molecular target for new anti-malarial discovery

Developing drug resistance in the malaria parasite is a reason for apprehension compelling the scientific community to focus on identifying new molecular targets that can be exploited for developing new anti-malarial compounds. Despite the availability of the Plasmodium genome, many protein-coding genes in Plasmodium are still not characterized or very less information is available about their functions. DMAP1 protein is known to be essential for growth and plays an important role in maintaining genomic integrity and transcriptional repression in vertebrate organisms. In this study, we have identified a homolog of DMAP1 in P. falciparum. Our sequence and structural analysis showed that although PfDMAP1 possesses a conserved SANT domain, parasite protein displays significant structural dissimilarities from human homolog at full-length protein level as well as within its SANT domain. PPIN analysis of PfDMAP1 revealed it to be vital for parasite and virtual High-throughput screening of various pharmacophore libraries using BIOVIA platform-identified compounds that pass ADMET profiling and showed specific binding with PfDMAP1. Based on MD simulations and protein-ligand interaction studies two best hits were identified that could be novel potent inhibitors of PfDMAP1 protein.Communicated by Ramaswamy H. Sarma.

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.

Assessing Protein-Protein Docking Protocols: Case Studies of G-Protein-Coupled Receptor Interactions

The interactions of G-protein-coupled receptors (GPCRs) with other proteins are critical in several cellular processes but resolving their structural dynamics remains challenging. An increasing number of GPCR complexes have been experimentally resolved but others including receptor variants are yet to be characterized, necessitating computational predictions of their interactions. Although integrative approaches with multi-scale simulations would provide rigorous estimates of their conformational dynamics, protein-protein docking remains a first tool of choice of many researchers due to the availability of open-source programs and easy to use web servers with reasonable predictive power. Protein-protein docking algorithms have limited ability to consider protein flexibility, environment effects, and entropy contributions and are usually a first step towards more integrative approaches. The two critical steps of docking: the sampling and scoring algorithms have improved considerably and their performance has been validated against experimental data. In this chapter, we provide an overview and generalized protocol of a few docking protocols using GPCRs as test cases. In particular, we demonstrate the interactions of GPCRs with extracellular protein ligands and an intracellular protein effectors (G-protein) predicted from docking approaches and test their limitations. The current chapter will help researchers critically assess docking protocols and predict experimentally testable structures of GPCR complexes.

Applying molecular docking to pesticides

Pesticide creation is related to the development of sustainable agricultural and ecological safety, and molecular docking technology can effectively help in pesticide innovation. This paper introduces the basic theory behind molecular docking, pesticide databases, and docking software. It also summarizes the application of molecular docking in the pesticide field, including the virtual screening of lead compounds, detection of pesticides and their metabolites in the environment, reverse screening of pesticide targets, and the study of resistance mechanisms. Finally, problems with the use of molecular docking technology in pesticide creation are discussed, and prospects for the future use of molecular docking technology in new pesticide development are discussed. © 2023 Society of Chemical Industry.

In silico designing of multiepitope-based-peptide (MBP) vaccine against MAPK protein express for Alzheimer's disease in Zebrafish

Understanding the role of the mitogen-activated protein kinases (MAPKs) signalling pathway is essential in advancing treatments for neurodegenerative disorders like Alzheimer's. In this study, we investigate in-silico techniques involving computer-based methods to extract the MAPK1 sequence. Our applied methods enable us to analyze the protein's structure, evaluate its properties, establish its evolutionary relationships, and assess its prevalence in populations. We also predict epitopes, assess their ability to trigger immune responses, and check for allergenicity using advanced computational tools to understand their immunological properties comprehensively. We apply virtual screening, docking, and structure modelling to identify promising drug candidates, analyze their interactions, and enhance drug design processes. We identified a total of 30 cell-targeting molecules against the MAPK1 protein, where we selected top 10 CTL epitopes (PAGGGPNPG, GGGPNPGSG, SAPAGGGPN, AVSAPAGGG, AGGGPNPGS, ATAAVSAPA, TAAVSAPAG, ENIIGINDI, INDIIRTPT, and NDIIRTPTI) for further evaluation to determine their potential efficacy, safety, and suitability for vaccine design based on strong binding potential. The potential to cover a large portion of the world's population with these vaccines is substantial-88.5 % for one type and 99.99 % for another. In exploring the molecular docking analyses, we examined a library of compounds from the ZINC database. Among them, we identified twelve compounds with the lowest binding energy. Critical residues in the MAPK1 protein, such as VAL48, LYS63, CYS175, ASP176, LYS160, ALA61, LEU165, TYR45, SER162, ARG33, PRO365, PHE363, ILE40, ASN163, and GLU42, are pivotal for interactions with these compounds. Our result suggests that these compounds could influence the protein's behaviour. Moreover, our docking analyses revealed that the predicted peptides have a strong affinity for the MAPK1 protein. These peptides form stable complexes, indicating their potential as potent inhibitors. This study contributes to the identification of new drug compounds and the screening of their desired properties. These compounds could potentially help reduce the excessive activity of MAPK1, which is linked to Alzheimer's disease.

DockSurf: A Molecular Modeling Software for the Prediction of Protein/Surface Adhesion

The elucidation of structural interfaces between proteins and inorganic surfaces is a crucial aspect of bionanotechnology development. Despite its significance, the interfacial structures between proteins and metallic surfaces are yet to be fully understood, and the lack of experimental investigation has impeded the development of many devices. To overcome this limitation, we suggest considering the generation of protein/surface structures as a molecular docking problem with a homogenous plan as the target. To this extent, we propose a new software, DockSurf, which aims to quickly propose reliable protein/surface structures. Our approach considers the conformational exploration with Euler's angles, which provide a cartography instead of a unique structure. Interaction energies were derived from quantum mechanics computations for a set of small molecules that describe protein atom types and implemented in a Derjaguin, Landau, Verwey, and Overbeek potential for the consideration of large systems such as proteins. The validation of DockSurf software was conducted with molecular dynamics for corona proteins with gold surfaces and provided enthusiastic results. This software is implemented in the RPBS platform to facilitate widespread access to the scientific community.

Critical Assessment of Methods for Predicting the 3D Structure of Proteins and Protein Complexes

Advances in a scientific discipline are often measured by small, incremental steps. In this review, we report on two intertwined disciplines in the protein structure prediction field, modeling of single chains and modeling of complexes, that have over decades emulated this pattern, as monitored by the community-wide blind prediction experiments CASP and CAPRI. However, over the past few years, dramatic advances were observed for the accurate prediction of single protein chains, driven by a surge of deep learning methodologies entering the prediction field. We review the mainscientific developments that enabled these recent breakthroughs and feature the important role of blind prediction experiments in building up and nurturing the structure prediction field. We discuss how the new wave of artificial intelligence-based methods is impacting the fields of computational and experimental structural biology and highlight areas in which deep learning methods are likely to lead to future developments, provided that major challenges are overcome.

In Silico Methodologies to Improve Antioxidants' Characterization from Marine Organisms

Marine organisms have been reported to be valuable sources of bioactive molecules that have found applications in different industrial fields. From organism sampling to the identification and bioactivity characterization of a specific compound, different steps are necessary, which are time- and cost-consuming. Thanks to the advent of the -omic era, numerous genome, metagenome, transcriptome, metatranscriptome, proteome and microbiome data have been reported and deposited in public databases. These advancements have been fundamental for the development of in silico strategies for basic and applied research. In silico studies represent a convenient and efficient approach to the bioactivity prediction of known and newly identified marine molecules, reducing the time and costs of "wet-lab" experiments. This review focuses on in silico approaches applied to bioactive molecule discoveries from marine organisms. When available, validation studies reporting a bioactivity assay to confirm the presence of an antioxidant molecule or enzyme are reported, as well. Overall, this review suggests that in silico approaches can offer a valuable alternative to most expensive approaches and proposes them as a little explored field in which to invest.

A Computational Framework for Determining the Breadth of Antibodies Against Highly Mutable Pathogens

Highly mutable pathogens pose daunting challenges for antibody design. The usual criteria of high potency and specificity are often insufficient to design antibodies that provide long-lasting protection. This is due, in part, to the ability of the pathogen to rapidly acquire mutations that permit them to evade the designed antibodies. To overcome these limitations, design of antibodies with a larger neutralizing breadth can be pursued. Such broadly neutralizing antibodies (bnAbs) should remain targeted to a specific epitope, yet show robustness against pathogen mutability, thereby neutralizing a higher number of antigens. This is particularly important for highly mutable pathogens, like the influenza virus and the human immunodeficiency virus (HIV). The protocol describes a method for computing the "breadth" of a given antibody, an essential aspect of antibody design.

In Silico and In Vitro Antiviral Activity Evaluation of Prodigiosin from Serratia marcescens Against Enterovirus 71

Prodigiosin, a red linear tripyrrole pigment found in Serratia marcescens, is one such naturally occurring compound that has gained wide attention owing to its numerous biological activities, including antibacterial, antifungal, antimalarial, anticancer, and immunosuppressive properties. This study was conducted to evaluate the possible antiviral activity of prodigiosin against Enterovirus 71, a causative agent of hand, foot, and mouth disease (HFMD). Preliminary studies were done in silico by analyzing the interaction of prodigiosin with amino acid residues of five EV71-target proteins. Interaction refinement analysis with FireDock revealed that 2C helicase (-48.01 kcal/moL) has the most negative global energy, followed by capsid (-36.52 kcal/moL), 3C protease (-34.16 kcal/moL), 3D RNA polymerase (-30.93 kcal/moL) and 2A protease (-20.61 kcal/moL). These values are indicative of the interaction strength. Prodigiosin was shown to form chemical bonds with specific amino acid residues in capsid (Gln-30, Asn-223), 2A protease (Trp-33, Trp-142), 2C helicase (Tyr-150, His-151, Gln-169, Ser-212), 3C protease (Glu-50), and 3D RNA polymerase (Ala-239, Tyr-237). To investigate further, prodigiosin was extracted from S. marcescens using a methanolic extraction method. In vitro studies revealed that prodigiosin, with an IC50 value of 0.5112 μg/mL, reduced virus titers by 0.17 log (32.39%) in 30 min and 0.19 log (35.43%) in 60 min. The findings suggest that prodigiosin has antiviral activity with an intermediate inhibitory effect against EV71. As a result of this research, new biological activities of prodigiosin have been identified.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

Integrating diverse plant bioactive ingredients with cyclodextrins to fabricate functional films for food application: a critical review

The popularity of plant bioactive ingredients has become increasingly apparent in the food industry. However, these plant bioactive ingredients have many deficiencies, including low water solubility, poor stability, and unacceptable odor. Cyclodextrins (CDs), as cyclic molecules, have been extensively studied as superb vehicles of plant bioactive ingredients. These CD inclusion compounds could be added into various film matrices to fabricate bioactive food packaging materials. Therefore, in the present review, we summarized the extraction methods of plant bioactive ingredients, the addition of these CD inclusion compounds into thin-film materials, and their applications in food packaging. Furthermore, the release model and mechanism of active film materials based on various plant bioactive ingredients with CDs were highlighted. Finally, the current challenges and new opportunities based on these film materials have been discussed.

In Silico and Experimental ADAM17 Kinetic Modeling as Basis for Future Screening System for Modulators

Understanding the mechanisms of modulators’ action on enzymes is crucial for optimizing and designing pharmaceutical substances. The acute inflammatory response, in particular, is regulated mainly by a disintegrin and metalloproteinase (ADAM) 17. ADAM17 processes several disease mediators such as TNFα and APP, releasing their soluble ectodomains (shedding). A malfunction of this process leads to a disturbed inflammatory response. Chemical protease inhibitors such as TAPI-1 were used in the past to inhibit ADAM17 proteolytic activity. However, due to ADAM17′s broad expression and activity profile, the development of active-site-directed ADAM17 inhibitor was discontinued. New ‘exosite’ (secondary substrate binding site) inhibitors with substrate selectivity raised the hope of a substrate-selective modulation as a promising approach for inflammatory disease therapy. This work aimed to develop a high-throughput screen for potential ADAM17 modulators as therapeutic drugs. By combining experimental and in silico methods (structural modeling and docking), we modeled the kinetics of ADAM17 inhibitor. The results explain ADAM17 inhibition mechanisms and give a methodology for studying selective inhibition towards the design of pharmaceutical substances with higher selectivity.

Designing a Recombinant Vaccine against Providencia rettgeri Using Immunoinformatics Approach

Antibiotic resistance (AR) is the resistance mechanism pattern in bacteria that evolves over some time, thus protecting the bacteria against antibiotics. AR is due to bacterial evolution to make itself fit to changing environmental conditions in a quest for survival of the fittest. AR has emerged due to the misuse and overuse of antimicrobial drugs, and few antibiotics are now left to deal with these superbug infections. To combat AR, vaccination is an effective method, used either therapeutically or prophylactically. In the current study, an in silico approach was applied for the design of multi-epitope-based vaccines against Providencia rettgeri, a major cause of traveler’s diarrhea. A total of six proteins: fimbrial protein, flagellar hook protein (FlgE), flagellar basal body L-ring protein (FlgH), flagellar hook-basal body complex protein (FliE), flagellar basal body P-ring formation protein (FlgA), and Gram-negative pili assembly chaperone domain proteins, were considered as vaccine targets and were utilized for B- and T-cell epitope prediction. The predicted epitopes were assessed for allergenicity, antigenicity, virulence, toxicity, and solubility. Moreover, filtered epitopes were utilized in multi-epitope vaccine construction. The predicted epitopes were joined with each other through specific GPGPG linkers and were joined with cholera toxin B subunit adjuvant via another EAAAK linker in order to enhance the efficacy of the designed vaccine. Docking studies of the designed vaccine construct were performed with MHC-I (PDB ID: 1I1Y), MHC-II (1KG0), and TLR-4 (4G8A). Findings of the docking study were validated through molecular dynamic simulations, which confirmed that the designed vaccine showed strong interactions with the immune receptors, and that the epitopes were exposed to the host immune system for proper recognition and processing. Additionally, binding free energies were estimated, which highlighted both electrostatic energy and van der Waals forces to make the complexes stable. Briefly, findings of the current study are promising and may help experimental vaccinologists to formulate a novel multi-epitope vaccine against P. rettgeri.

Application of molecular docking in elaborating molecular mechanisms and interactions of supramolecular cyclodextrin

The cyclodextrin (CD)-based supramolecular nanomedicines have attracted growing interest because of their superior characteristics, including desirable biocompatibility, low toxicity, unique molecular structure and easy functionalization. The smart structures of CD impart host-guest interaction for meeting the multifunctional needs of disease therapy. However, it faces challenges in formulation design and inclusion mechanism clarification of the functional supramolecular assemblies owing to the complicated structures and mechanisms. Fortunately, molecular docking helps the researchers to comprehend the interaction between the drug and the target molecule for achieving high-through screening from the database. In this review, we summarized the category and characteristics of molecular docking along with the properties and applications of CD. Significantly, we highlighted the application of molecular docking in elaborating molecular mechanisms and simulating complex structures at molecular levels. The issues and development of CD and molecular docking were also presented to provide beneficial reference and new insights for supramolecular nano-systems.

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.

Molecular docking-based computational platform for high-throughput virtual screening

Structure-based virtual screening is a key, routine computational method in computer-aided drug design. Such screening can be used to identify potentially highly active compounds, to speed up the progress of novel drug design. Molecular docking-based virtual screening can help find active compounds from large ligand databases by identifying the binding affinities between receptors and ligands. In this study, we analyzed the challenges of virtual screening, with the aim of identifying highly active compounds faster and more easily than is generally possible. We discuss the accuracy and speed of molecular docking software and the strategy of high-throughput molecular docking calculation, and we focus on current challenges and our solutions to these challenges of ultra-large-scale virtual screening. The development of Web services helps lower the barrier to drug virtual screening. We introduced some related web sites for docking and virtual screening, focusing on the development of pre- and post-processing interactive visualization and large-scale computing.

CGRAP: A Web Server for Coarse-Grained Rigidity Analysis of Proteins

Elucidating protein rigidity offers insights about protein conformational changes. An understanding of protein motion can help speed drug development, and provide general insights into the dynamic behaviors of biomolecules. Existing rigidity analysis techniques employ fine-grained, all-atom modeling, which has a costly run-time, particularly for proteins made up of more than 500 residues. In this work, we introduce coarse-grained rigidity analysis, and showcase that it provides flexibility information about a protein that is similar in accuracy to an all-atom modeling approach. We assess the accuracy of the coarse-grained method relative to an all-atom approach via a comparison metric that reasons about the largest rigid clusters of the two methods. The apparent symmetry between the all-atom and coarse-grained methods yields very similar results, but the coarse-grained method routinely exhibits 40% reduced run-times. The CGRAP web server outputs rigid cluster information, and provides data visualization capabilities, including a interactive protein visualizer.

PIK3CA hotspot mutations p. H1047R and p. H1047L sensitize breast cancer cells to thymoquinone treatment by regulating the PI3K/Akt1 pathway

BACKGROUND

Nigella sativa (N. sativa) exhibits anti-inflammatory, antioxidant, antidiabetic, antimetastatic and antinociceptive effects and has been used to treat dozens of diseases. Thymoquinone (TQ) is an important and active component isolated from N. sativa seeds. Inhibition of cancer-associated activating PIK3CA mutations is a new prospective targeted therapy in personalized metastatic breast cancer (MBC). TQ is reported to be an effective inhibitor of the PI3K/Akt1 pathway in MBC. This study aimed to evaluate the in vitro antitumor effect of TQ in the context of two PIK3CA hotspot mutations, p. H1047R and p. H1047L.

METHODS AND RESULTS

Molecular dynamics, free energy landscapes and principal component analyses were also used to survey the mechanistic effects of the p. H1047R and p. H1047L mutations on the PI3K/Akt1 pathway. Our findings clearly confirmed that the p. H1047R and p. H1047L mutants could reduce the inhibitory effect of ΔNp63α on the kinase domain of PIK3CA, resulting in increased activity of PI3K downstream signals. Structurally, the partial disruption of the interaction between the ΔNp63α DNA binding domain and the PIK3CA kinase domain at residues 114-359 and 797-1068 destabilizes the conformation of the activation loop and modifies the PIK3CA/ΔNp63α complex. Alongside these structural changes, we found that TQ treatment resulted in high PI3K/Akt1 pathway inhibition in p. H1047R and p. H1047L-expressing cells versus wild-type cells.

CONCLUSIONS

These two PIK3CA hotspot mutations therefore not only contribute to tumor progression in patients with MBC but may also serve as targets for the development of novel small molecule therapeutic strategies.

Cross-Linking Mass Spectrometry for Investigating Protein Conformations and Protein-Protein Interactions─A Method for All Seasons

Mass spectrometry (MS) has become one of the key technologies of structural biology. In this review, the contributions of chemical cross-linking combined with mass spectrometry (XL-MS) for studying three-dimensional structures of proteins and for investigating protein-protein interactions are outlined. We summarize the most important cross-linking reagents, software tools, and XL-MS workflows and highlight prominent examples for characterizing proteins, their assemblies, and interaction networks in vitro and in vivo. Computational modeling plays a crucial role in deriving 3D-structural information from XL-MS data. Integrating XL-MS with other techniques of structural biology, such as cryo-electron microscopy, has been successful in addressing biological questions that to date could not be answered. XL-MS is therefore expected to play an increasingly important role in structural biology in the future.

Native or Non-Native Protein-Protein Docking Models? Molecular Dynamics to the Rescue

Molecular docking excels at creating a plethora of potential models of protein-protein complexes. To correctly distinguish the favorable, native-like models from the remaining ones remains, however, a challenge. We assessed here if a protocol based on molecular dynamics (MD) simulations would allow distinguishing native from non-native models to complement scoring functions used in docking. To this end, the first models for 25 protein-protein complexes were generated using HADDOCK. Next, MD simulations complemented with machine learning were used to discriminate between native and non-native complexes based on a combination of metrics reporting on the stability of the initial models. Native models showed higher stability in almost all measured properties, including the key ones used for scoring in the Critical Assessment of PRedicted Interaction (CAPRI) competition, namely the positional root mean square deviations and fraction of native contacts from the initial docked model. A random forest classifier was trained, reaching a 0.85 accuracy in correctly distinguishing native from non-native complexes. Reasonably modest simulation lengths of the order of 50-100 ns are sufficient to reach this accuracy, which makes this approach applicable in practice.

Reconstruction of the Fas-Based Death-Inducing Signaling Complex (DISC) Using a Protein-Protein Docking Meta-Approach

The death-inducing signaling complex (DISC) is a fundamental multiprotein complex, which triggers the extrinsic apoptosis pathway through stimulation by death ligands. DISC consists of different death domain (DD) and death effector domain (DED) containing proteins such as the death receptor Fas (CD95) in complex with FADD, procaspase-8, and cFLIP. Despite many experimental and theoretical studies in this area, there is no global agreement neither on the DISC architecture nor on the mechanism of action of the involved species. In the current work, we have tried to reconstruct the DISC structure by identifying key protein interactions using a new protein-protein docking meta-approach. We combined the benefits of five of the most employed protein-protein docking engines, HADDOCK, ClusPro, HDOCK, GRAMM-X, and ZDOCK, in order to improve the accuracy of the predicted docking complexes. Free energy of binding and hot spot interacting residues were calculated and determined for each protein-protein interaction using molecular mechanics generalized Born surface area and alanine scanning techniques, respectively. In addition, a series of in-cellulo protein-fragment complementation assays were conducted to validate the protein-protein docking procedure. The results show that the DISC formation initiates by dimerization of adjacent Fas_DD trimers followed by recruitment of FADD through homotypic DD interactions with the oligomerized death receptor. Furthermore, the in-silico outcomes indicate that cFLIP cannot bind directly to FADD; instead, cFLIP recruitment to the DISC is a hierarchical and cooperative process where FADD initially recruits procaspase-8, which in turn recruits and heterodimerizes with cFLIP. Finally, a possible structure of the entire DISC is proposed based on the docking results.

IsAb: a computational protocol for antibody design

The design of therapeutic antibodies has attracted a large amount of attention over the years. Antibodies are widely used to treat many diseases due to their high efficiency and low risk of adverse events. However, the experimental methods of antibody design are time-consuming and expensive. Although computational antibody design techniques have had significant advances in the past years, there are still some challenges that need to be solved, such as the flexibility of antigen structure, the lack of antibody structural data and the absence of standard antibody design protocol. In the present work, we elaborated on an in silico antibody design protocol for users to easily perform computer-aided antibody design. First, the Rosetta web server will be applied to generate the 3D structure of query antibodies if there is no structural information available. Then, two-step docking will be used to identify the binding pose of an antibody-antigen complex when the binding information is unknown. ClusPro is the first method to be used to conduct the global docking, and SnugDock is applied for the local docking. Sequentially, based on the predicted binding poses, in silico alanine scanning will be used to predict the potential hotspots (or key residues). Finally, computational affinity maturation protocol will be used to modify the structure of antibodies to theoretically increase their affinity and stability, which will be further validated by the bioassays in the future. As a proof of concept, we redesigned antibody D44.1 and compared it with previously reported data in order to validate IsAb protocol. To further illustrate our proposed protocol, we used cemiplimab antibody, a PD-1 checkpoint inhibitor, as an example to showcase a step-by-step tutorial.

Field-Template, QSAR, Ensemble Molecular Docking, and 3D-RISM Solvation Studies Expose Potential of FDA-Approved Marine Drugs as SARS-CoVID-2 Main Protease Inhibitors

Currently, SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has infected people among all countries and is a pandemic as declared by the World Health Organization (WHO). SARS-CoVID-2 main protease is one of the therapeutic drug targets that has been shown to reduce virus replication, and its high-resolution 3D structures in complex with inhibitors have been solved. Previously, we had demonstrated the potential of natural compounds such as serine protease inhibitors eventually leading us to hypothesize that FDA-approved marine drugs have the potential to inhibit the biological activity of SARS-CoV-2 main protease. Initially, field-template and structure–activity atlas models were constructed to understand and explain the molecular features responsible for SARS-CoVID-2 main protease inhibitors, which revealed that Eribulin Mesylate, Plitidepsin, and Trabectedin possess similar characteristics related to SARS-CoVID-2 main protease inhibitors. Later, protein–ligand interactions are studied using ensemble molecular-docking simulations that revealed that marine drugs bind at the active site of the main protease. The three-dimensional reference interaction site model (3D-RISM) studies show that marine drugs displace water molecules at the active site, and interactions observed are favorable. These computational studies eventually paved an interest in further in vitro studies. Finally, these findings are new and indeed provide insights into the role of FDA-approved marine drugs, which are already in clinical use for cancer treatment as a potential alternative to prevent and treat infected people with SARS-CoV-2.

Efficient Refinement and Free Energy Scoring of Predicted Protein-Protein Complexes Using Replica Exchange with Repulsive Scaling

Accurate prediction and evaluation of protein-protein complex structures are of major importance to understand the cellular interactome. Typically, putative complexes are predicted based on docking methods, and simple force field or knowledge-based scoring functions are applied to evaluate single complex structures. We have extended a replica-exchange-based scheme employing different levels of a repulsive biasing between partners in each replica simulation (RS-REMD) to simultaneously refine and score protein-protein complexes. The bias acts specifically on the intermolecular interactions based on an increase in effective pairwise van der Waals radii (repulsive scaling (RS)-REMD) without affecting interactions within each protein or with the solvent. The method provides a free energy score that correlates quite well with experimental binding free energies on a set of 36 complexes with correlation coefficients of 0.77 and 0.55 in explicit and implicit solvent simulations, respectively. For a large set of docked decoy complexes, significant improvement of docked complexes was found in many cases with the starting structure in the vicinity (within 20 Å) of the native complex. In the majority of cases (14 out of 20 in explicit solvent), near native docking solutions were identified as the best scoring complexes. The approach is computational demanding but may offer a route for refinement and realistic ranking of predicted protein-protein docking geometries.

PRosettaC: Rosetta Based Modeling of PROTAC Mediated Ternary Complexes

Proteolysis-targeting chimeras (PROTACs), which induce degradation by recruitment of an E3 ligase to a target protein, are gaining much interest as a new pharmacological modality. However, designing PROTACs is challenging. Formation of a ternary complex between the protein target, the PROTAC, and the recruited E3 ligase is considered paramount for successful degradation. A structural model of this ternary complex could in principle inform rational PROTAC design. Unfortunately, only a handful of structures are available for such complexes, necessitating tools for their modeling. We developed a combined protocol for the modeling of a ternary complex induced by a given PROTAC. Our protocol alternates between sampling of the protein-protein interaction space and the PROTAC molecule conformational space. Application of this protocol-PRosettaC-to a benchmark of known PROTAC ternary complexes results in near-native predictions, with often atomic accuracy prediction of the protein chains, as well as the PROTAC binding moieties. It allowed the modeling of a CRBN/BTK complex that recapitulated experimental results for a series of PROTACs. PRosettaC generated models may be used to design PROTACs for new targets, as well as improve PROTACs for existing targets, potentially cutting down time and synthesis efforts. To enable wide access to this protocol, we have made it available through a web server (https://prosettac.weizmann.ac.il/).

GWOVina: A grey wolf optimization approach to rigid and flexible receptor docking

Protein–ligand docking programs are indispensable tools for predicting the binding pose of a ligand to the receptor protein. In this paper, we introduce an efficient flexible docking method, gwovina, which is a variant of the Vina implementation using the grey wolf optimizer (GWO) and random walk for the global search, and the Dunbrack rotamer library for side‐chain sampling. The new method was validated for rigid and flexible‐receptor docking using four independent datasets. In rigid docking, gwovina showed comparable docking performance to Vina in terms of ligand pose RMSD, success rate, and affinity prediction. In flexible‐receptor docking, gwovina has improved success rate compared to Vina and AutoDockFR. It ran 2 to 7 times faster than Vina and 40 to 100 times faster than AutoDockFR. Therefore, gwovina can play a role in solving the complex flexible‐receptor docking cases and is suitable for virtual screening of compound libraries. gwovina is freely available at https://cbbio.cis.um.edu.mo/software/gwovina for testing.

Developing Computational Model to Predict Protein-Protein Interaction Sites Based on the XGBoost Algorithm

The study of protein-protein interaction is of great biological significance, and the prediction of protein-protein interaction sites can promote the understanding of cell biological activity and will be helpful for drug development. However, uneven distribution between interaction and non-interaction sites is common because only a small number of protein interactions have been confirmed by experimental techniques, which greatly affects the predictive capability of computational methods. In this work, two imbalanced data processing strategies based on XGBoost algorithm were proposed to re-balance the original dataset from inherent relationship between positive and negative samples for the prediction of protein-protein interaction sites. Herein, a feature extraction method was applied to represent the protein interaction sites based on evolutionary conservatism of proteins, and the influence of overlapping regions of positive and negative samples was considered in prediction performance. Our method showed good prediction performance, such as prediction accuracy of 0.807 and MCC of 0.614, on an original dataset with 10,455 surface residues but only 2297 interface residues. Experimental results demonstrated the effectiveness of our XGBoost-based method.

Smart Collagen Hydrogels Based on 1-Ethyl-3-methylimidazolium Acetate and Microbial Transglutaminase for Potential Applications in Tissue Engineering and Cancer Therapy

For the first time, collagen-based hydrogels were fabricated in the presence of a biocompatible ionic liquid, 1-ethyl-3-methylimidazolium acetate ([EMIM] [Ac]), by a simple biopolymer cross-linking process facilitated by the strong catalytic hydrolysis of microbial transglutaminase (MTGase). Phosphate buffer solution (PBS)-encapsulated human-like collagen (HLC) or fish bone collagen (FBC) for the composite hydrogels was simply prepared by the codissolution of biopolymers in [EMIM] [Ac] or, in the absence of the ionic liquid, by the dispersion of MTGase in the biopolymer solution, leading to the formation of MTGase-aided hydrogels (Gel1 and Gel4) and [EMIM] [Ac]/MTGase-aided hydrogels (Gel2, Gel3, and Gel5). The effects of different contents of [EMIM] [Ac] and collagens of different origins (HLC and FBC) during fabrication on a range of structural and material characteristics, including the synthesis mechanism, three-dimensional structure, swelling behavior, mechanical strength, enzymatic hydrolysis rate, cytotoxicity, fibroblast cell proliferation rate, in vitro inhibition of cancer cells and cell adhesion, and in vivo histocompatibility, were investigated. Surprisingly, fabrication with [EMIM] [Ac] had significant effects on the structure and properties of the collagen/MTGase-based hydrogels. In other words, [EMIM] [Ac] changed the underlying mechanism responsible for the advantageous properties of the hydrogels by changing the three-dimensional structure of HLC or FBC, which improved their effects on fibroblast proliferation (3T3-L1 and L929 cells) and their in vitro inhibition of cancer cells (HepG2 and MKN45 cells). The use of the ionic liquid also imbued the hydrogels with degradation resistance and anti-inflammatory properties after subcutaneous injection into mice (in vivo). The catalytic hydrolysis by MTGase and the [EMIM] [Ac] content were the major factors that influenced the properties of the collagen. This result suggests the potential application of ionic liquid-enzymatic hydrolysis in the fabrication of collagen hydrogels in circumstances where the control of the properties by an ionic liquid is desirable. Therefore, [EMIM] [Ac] could be a promising solvent for the development of collagen into smart biomaterials with controlled biodegradation rates that can meet the needs of specific potential applications, such as tissue engineering and cancer therapy.

EROS-DOCK: protein–protein docking using exhaustive branch-and-bound rotational search

Motivation

Protein–protein docking algorithms aim to predict the 3D structure of a binary complex using the structures of the individual proteins. This typically involves searching and scoring in a 6D space. Many docking algorithms use FFT techniques to exhaustively cover the search space and to accelerate the scoring calculation. However, FFT docking results often depend on the initial protein orientations with respect to the Fourier sampling grid. Furthermore, Fourier-transforming a physics-base force field can involve a serious loss of precision.

Results

Here, we present EROS-DOCK, an algorithm to rigidly dock two proteins using a series of exhaustive 3D rotational searches in which non-clashing orientations are scored using the ATTRACT coarse-grained force field model. The rotational space is represented as a quaternion ‘π-ball’, which is systematically sub-divided in a ‘branch-and-bound’ manner, allowing efficient pruning of rotations that will give steric clashes. The algorithm was tested on 173 Docking Benchmark complexes, and results were compared with those of ATTRACT and ZDOCK. According to the CAPRI quality criteria, EROS-DOCK typically gives more acceptable or medium quality solutions than ATTRACT and ZDOCK.

Availability and implementation

The EROS-DOCK program is available for download at http://erosdock.loria.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.

Evaluation of Predicted Protein-Protein Complexes by Binding Free Energy Simulations

The accurate prediction of protein-protein complex geometries is of major importance to ultimately model the complete interactome of interacting proteins in a cell. A major bottleneck is the realistic free energy evaluation of predicted docked structures. Typically, simple scoring functions applied to single-complex structures are employed that neglect conformational entropy and often solvent effects completely. The binding free energy of a predicted protein-protein complex can, however, be calculated using umbrella sampling (US) along a predefined dissociation/association coordinate of a complex. We employed atomistic US-molecular dynamics simulations including appropriate conformational and axial restraints and an implicit generalized Born solvent model to calculate binding free energies of a large set of docked decoys for 20 different complexes. Free energies associated with the restraints were calculated separately. In principle, the approach includes all energetic and entropic contributions to the binding process. The evaluation of docked complexes based on binding free energy calculation was in better agreement with experiment compared to a simple scoring based on energy minimization or MD refinement using exactly the same force field description. Even calculated absolute binding free energies of structures close to the native binding geometry showed a reasonable correlation to experiment. However, still for a number of complexes docked decoys of lower free energy than near-native geometries were found indicating inaccuracies in the force field or the implicit solvent model. Although time consuming the approach may open up a new route for realistic ranking of predicted geometries based on calculated free energy of binding.

DaReUS-Loop: accurate loop modeling using fragments from remote or unrelated proteins

Despite efforts during the past decades, loop modeling remains a difficult part of protein structure modeling. Several approaches have been developed in the framework of crystal structures. However, for homology models, the modeling of loops is still far from being solved. We propose DaReUS-Loop, a data-based approach that identifies loop candidates mining the complete set of experimental structures available in the Protein Data Bank. Candidate filtering relies on local conformation profile-profile comparison, together with physico-chemical scoring. Applied to three different template-based test sets, DaReUS-Loop shows significant increase in the number of high-accuracy loops, and significant enhancement for modeling long loops. A special advantage is that our method proposes a prediction confidence score that correlates well with the expected accuracy of the loops. Strikingly, over 50% of successful loop models are derived from unrelated proteins, indicating that fragments under similar constraints tend to adopt similar structure, beyond mere homology.

ProPOSE: Direct Exhaustive Protein-Protein Docking with Side Chain Flexibility

Despite decades of development, protein-protein docking remains a largely unsolved problem. The main difficulties are the immense space spanned by the translational and rotational degrees of freedom and the prediction of the conformational changes of proteins upon binding. FFT is generally the preferred method to exhaustively explore the translation-rotation space at a fine grid resolution, albeit with the trade-off of approximating force fields with correlation functions. This work presents a direct search alternative that samples the states in Cartesian space at the same resolution and computational cost as standard FFT methods. Operating in real space allows the use of standard force field functional forms used in typical non-FFT methods as well as the implementation of strategies for focused exploration of conformational flexibility. Currently, a few misplaced side chains can cause docking programs to fail. This work specifically addresses the problem of side chain rearrangements upon complex formation. Based on the observation that most side chains retain their unbound conformation upon binding, each rigidly docked pose is initially scored ignoring up to a limited number of side chain overlaps which are resolved in subsequent repacking and minimization steps. On test systems where side chains are altered and backbones held in their bound state, this implementation provides significantly better native pose recovery and higher quality (lower RMSD) predictions when compared with five of the most popular docking programs. The method is implemented in the software program ProPOSE (Protein Pose Optimization by Systematic Enumeration).

Development of a new benchmark for assessing the scoring functions applicable to protein–protein interactions

Aim: Scoring functions are important component of protein–protein docking methods. They need to be evaluated on high-quality benchmarks to reveal their strengths and weaknesses. Evaluation results obtained on such benchmarks can provide valuable guidance for developing more advanced scoring functions. Methodology & results: In our comparative assessment of scoring functions for protein–protein interactions benchmark, the performance of a scoring function was characterized by ‘docking power’ and ‘scoring power’. A high-quality dataset of 273 protein–protein complexes was compiled and employed in both tests. Four scoring functions, including FASTCONTACT, ZRANK, dDFIRE and ATTRACT were tested as demonstration. ZRANK and ATTRACT exhibited encouraging performance in the docking power test. However, all four scoring functions failed badly in the scoring power test. Conclusion: Our comparative assessment of scoring functions for protein–protein interaction benchmark is created especially for assessing the scoring functions applicable to protein–protein interactions. It is different from other benchmarks for assessing protein–protein docking methods. Our benchmark is available to the public at www.pdbbind-cn.org/download/CASF-PPI/ .

Conformational changes in antibody Fab fragments upon binding and their consequences on the performance of docking algorithms

BACKGROUND

The existence of conformational changes in antibodies upon binding has been previously established. However, existing analyses focus on individual cases and no quantitative study provides a more global view of potential moves and repacking, especially on recent data. The present study focuses on analyzing the conformational changes in various antibodies upon binding, providing quantitative observations to be exploited for antibody-related modeling.

METHODS

Cartesian and dihedral Root-Mean-Squared Deviations were calculated for different subparts of 27 different antibodies, for which X-ray structures in the bound and unbound states are available. Elbow angle variations were also calculated. Previously reported results of four docking algorithms were condensed into one score giving overall docking success for each of 16 antibody-antigen cases.

RESULTS

Very diverse movements are observed upon binding. While many loops stay very rigid, several others display side-chain repacking or backbone rearrangements, or both, at many different levels. Large conformational changes restricted to one or more antibody hypervariable loops were found to be a better indicator of docking difficulty than overall conformational variation at the antibody-antigen interface. However, the failure of docking algorithms on some almost-rigid cases shows that scoring is still a major bottleneck in docking pose prediction.

CONCLUSIONS

This study is aimed to help scientists working on antibody analysis and design by giving insights into the nature and the extent of conformational changes at different levels upon antigen binding.

Morphing and docking visualisation of biomolecular structures using Multi-Dimensional Scaling

Protein structures are often solved at atomic resolution in two states defining a functional movement but intervening conformations are usually unknown. Morphing methods generate intervening conformations between two known structures. When viewed as an animation using molecular graphics, a smooth, direct morph enables the eye to track changes in structure that might be otherwise missed. We present a morphing method that aims to linearly interpolate interatomic distances and which uses SMACOF (Scaling by MAjorisation of COmplicated Function) and multigrid techniques with a cut-off distance based weighting that optimizes the MolProbity score of intervening structures. The all-atom morphs are smooth, move directly between the two structures, and are shown, in general, to pass closer to a set of known intermediates than those generated using other methods. The techniques are also used for docking by putting the unbound structures in a "near-approach pose" and then morphing to the bound complex. The resulting GPU-accelerated tools are available on a webserver, Morphit_Pro, at http://morphit-pro.cmp.uea.ac.uk/ and more than 5000 domains movements available at the DynDom website can now be viewed as morphs http://morphit-pro.cmp.uea.ac.uk/dyndom/.

Site-Specific Studies of Nucleosome Interactions by Solid-State NMR Spectroscopy

Chromatin function depends on a dense network of interactions between nucleosomes and a wide range of proteins. A detailed description of these protein-nucleosome interactions is required to reach a full molecular understanding of chromatin function in both genetics and epigenetics. Herein, we show that the structure, dynamics, and interactions of nucleosomes can be interrogated in a residue-specific manner by using state-of-the-art solid-state NMR spectroscopy. Using sedimented nucleosomes, high-resolution spectra were obtained for both flexible histone tails and the non-mobile histone core. Through co-sedimentation of a nucleosome-binding peptide, we demonstrate that protein-binding sites on the nucleosome surface can be determined. We believe that this approach holds great promise as it is generally applicable, extendable to include the structure and dynamics of the bound proteins, and scalable to interactions of proteins with higher-order chromatin structures, including isolated and cellular chromatin.

Prediction of spacer-α6 complex: a novel insight into binding of ADAMTS13 with A2 domain of von Willebrand factor under forces

Force-regulated cleavage of A2 domain of von Willebrand factor (vWF) by ADAMTS13 is a key event in preventing thrombotic thrombocytopenic purpura (TTP). Recognition and cleavage depend on cooperative and modular contacts between several ADAMTS13 subdomains and discrete segments of vWF A2 domain. Spacer domain of ADAMTS13 contains an important exosite interacting with α6 helix of unfold A2 domain, but it remains unclear whether stretching of α6 regulates binding to spacer. To understand the molecular mechanism underlying the interactions between spacer and α6 under stretching, we successfully predicted spacer-α6 complex by a novel computer strategy combined the steered molecular dynamics (SMD) and flexible docking techniques. This strategy included three steps: (1) constant-velocity SMD simulation of α6; (2) zero-velocity SMD simulations of α6, and (3) flexible dockings of α6 to spacer. In our spacer-α6 complex model, 13 key residues, six in α6 and seven in spacer, were identified. Our data demonstrated a biphasic extension-regulated binding of α6 to spacer. The binding strength of the complex increased with α6 extension until it reaches its optimum of 0.25 nm, and then decreased as α6 extension further increased, meaning that spacer is in favor to binding with a partially extended α6, which may contribute to the optimal contact and proteolysis. Changes of interface area and intermolecular salt bridge may serve as the molecular basis for this characteristic. These findings provide a novel insight into mechano-chemical regulation on interaction between ADAMTS13 and vWF A2 domain under forces.

SnapDock-template-based docking by Geometric Hashing

Motivation

A highly efficient template-based protein–protein docking algorithm, nicknamed SnapDock, is presented. It employs a Geometric Hashing-based structural alignment scheme to align the target proteins to the interfaces of non-redundant protein–protein interface libraries. Docking of a pair of proteins utilizing the 22 600 interface PIFACE library is performed in < 2 min on the average. A flexible version of the algorithm allowing hinge motion in one of the proteins is presented as well.

Results

To evaluate the performance of the algorithm a blind re-modelling of 3547 PDB complexes, which have been uploaded after the PIFACE publication has been performed with success ratio of about 35%. Interestingly, a similar experiment with the template free PatchDock docking algorithm yielded a success rate of about 23% with roughly 1/3 of the solutions different from those of SnapDock. Consequently, the combination of the two methods gave a 42% success ratio.

Availability and implementation

A web server of the application is under development.

Combining Stochastic Deformation/Relaxation and Intermolecular Contacts Analysis for Extracting Pharmacophores from Ligand-Receptor Complexes

We previously combined molecular dynamics (classical or simulated annealing) with ligand-receptor contacts analysis as a means to extract valid pharmacophore model(s) from single ligand-receptor complexes. However, molecular dynamics methods are computationally expensive and time-consuming. Here we describe a novel method for extracting valid pharmacophore model(s) from a single crystallographic structure within a reasonable time scale. The new method is based on ligand-receptor contacts analysis following energy relaxation of a predetermined set of randomly deformed complexes generated from the targeted crystallographic structure. Ligand-receptor contacts maintained across many deformed/relaxed structures are assumed to be critical and used to guide pharmacophore development. This methodology was implemented to develop valid pharmacophore models for PI3K-γ, RENIN, and JAK1. The resulting pharmacophore models were validated by receiver operating characteristic (ROC) analysis against inhibitors extracted from the CHEMBL database. Additionally, we implemented pharmacophores extracted from PI3K-γ to search for new inhibitors from the National Cancer Institute list of compounds. The process culminated in new PI3K-γ/mTOR inhibitory leads of low micromolar IC₅₀s.

The computational prediction of protein assemblies

The function of proteins in the cell is almost always mediated by their interaction with different partners, including other proteins, nucleic acids or small organic molecules. The ability of identifying all of them is an essential step in our quest for understanding life at the molecular level. The inference of the protein complex composition and of its molecular details can also provide relevant clues for the development and the design of drugs. In this short review, I will discuss the computational aspects of the analysis and prediction of protein-protein assemblies and discuss some of the most recent developments as seen in the last Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiment.

Steered molecular dynamics simulation of the binding of the bovine auxilin J domain to the Hsc70 nucleotide-binding domain

The Hsp70 and Hsp40 chaperone machine plays critical roles in protein folding, membrane translocation, and protein degradation by binding and releasing protein substrates in a process that utilizes ATP. The activities of the Hsp70 family of chaperones are recruited and stimulated by the J domains of Hsp40 chaperones. However, structural information on the Hsp40-Hsp70 complex is lacking, and the molecular details of this interaction are yet to be elucidated. Here we used steered molecular dynamics (SMD) simulations to investigate the molecular interactions that occur during the dissociation of the auxilin J domain from the Hsc70 nucleotide-binding domain (NBD). The changes in energy observed during the SMD simulation suggest that electrostatic interactions are the dominant type of interaction. Additionally, we found that Hsp70 mainly interacts with auxilin through the surface residues Tyr866, Arg867, and Lys868 of helix II, His874, Asp876, Lys877, Thr879, and Gln881 of the HPD loop, and Phe891, Asn895, Asp896, and Asn903 of helix III. The conservative residues Tyr866, Arg867, Lys868, His874, Asp876, Lys877, and Phe891 were also found in a previous study to be indispensable to the catalytic activity of the DnaJ J domain and the binding of it with the NBD of DnaK. The in silico identification of the importance of auxilin residues Asn895, Asp896, and Asn903 agrees with previous mutagenesis and NMR data suggesting that helix III of the J domain of the T antigen interacts with Hsp70. Furthermore, our data indicate that Thr879 and Gln881 from the HPD loop are also important as they mediate the interaction between the bovine auxilin J domain and Hsc70.

PEPSI-Dock: a detailed data-driven protein-protein interaction potential accelerated by polar Fourier correlation

MOTIVATION

Docking prediction algorithms aim to find the native conformation of a complex of proteins from knowledge of their unbound structures. They rely on a combination of sampling and scoring methods, adapted to different scales. Polynomial Expansion of Protein Structures and Interactions for Docking (PEPSI-Dock) improves the accuracy of the first stage of the docking pipeline, which will sharpen up the final predictions. Indeed, PEPSI-Dock benefits from the precision of a very detailed data-driven model of the binding free energy used with a global and exhaustive rigid-body search space. As well as being accurate, our computations are among the fastest by virtue of the sparse representation of the pre-computed potentials and FFT-accelerated sampling techniques. Overall, this is the first demonstration of a FFT-accelerated docking method coupled with an arbitrary-shaped distance-dependent interaction potential.

RESULTS

First, we present a novel learning process to compute data-driven distant-dependent pairwise potentials, adapted from our previous method used for rescoring of putative protein-protein binding poses. The potential coefficients are learned by combining machine-learning techniques with physically interpretable descriptors. Then, we describe the integration of the deduced potentials into a FFT-accelerated spherical sampling provided by the Hex library. Overall, on a training set of 163 heterodimers, PEPSI-Dock achieves a success rate of 91% mid-quality predictions in the top-10 solutions. On a subset of the protein docking benchmark v5, it achieves 44.4% mid-quality predictions in the top-10 solutions when starting from bound structures and 20.5% when starting from unbound structures. The method runs in 5-15 min on a modern laptop and can easily be extended to other types of interactions.

AVAILABILITY AND IMPLEMENTATION

https://team.inria.fr/nano-d/software/PEPSI-Dock

CONTACT

sergei.grudinin@inria.fr.