What method to use for protein-protein docking?

Quantum Chemistry-Based Protein-Protein Docking without Empirical Parameters

This study developed a novel protein-protein docking approach based on quantum chemistry. To judge the appropriateness of complex structures, we introduced two criterion values, EV1 and EV2, computed using the fragment molecular orbital method without any empirical parameters. These criterion values enable us to search complex structures in which patterns of the electrostatic potential of the two proteins are optimally aligned at their interface. The performance of our method was validated using 53 complexes in a benchmark set provided for protein-protein docking. When employing bound state structures, docking success rates reached 64% for EV1 and 76% for EV2. On the other hand, when employing unbound state structures, docking success rates reached 13% for EV1 and 17% for EV2.

Simulation and Computational Study of RING Domain Mutants of BRCA1 and Ube2k in AD/PD Pathophysiology

Lysine-based post-translational modification (PTM) such as acylation, acetylation, deamination, methylation, SUMOylation, and ubiquitination has proven to be a major regulator of gene expression, chromatin structure, protein stability, protein-protein interaction, protein degradation, and cellular localization. However, besides all the PTMs, ubiquitination stands as the second most common PTM after phosphorylation that is involved in the etiology of neurodegenerative diseases (NDDs) namely, Alzheimer's disease (AD) and Parkinson's disease (PD). NDDs are characterized by the accumulation of misfolded protein aggregates in the brain that lead to disease-related gene mutation and irregular protein homeostasis. The ubiquitin-proteasome system (UPS) is in charge of degrading these misfolded proteins, which involve an interplay of E1, E2, E3, and deubiquitinase enzymes. Impaired UPS has been commonly observed in NDDs and E3 ligases are the key members of the UPS, thus, dysfunction of the same can accelerate the neurodegeneration process. Therefore, the aim of this study is firstly, to find E3 ligases that are common in both AD and PD through data mining. Secondly, to study the impact of mutation on its structure and function. The study deciphered 74 E3 ligases that were common in both AD and PD. Later, 10 hub genes were calculated of which protein-protein interaction, pathway enrichment, lysine site prediction, domain, and motif analysis were performed. The results predicted BRCA1, PML, and TRIM33 as the top three putative lysine-modified E3 ligases involved in AD and PD pathogenesis. However, based on structural characterization, BRCA1 was taken further to study RING domain mutation that inferred K32Y, K32L, K32C, K45V, K45Y, and K45G as potential mutants that alter the structural and functional ability of BRCA1 to interact with Ube2k, E2-conjugating enzyme. The most probable mutant observed after molecular dynamics simulation of 50 ns is K32L. Therefore, our study concludes BRCA1, a potential E3 ligase common in AD and PD, and RING domain mutation at sites K32 and K45 possibly disturbs its interaction with its E2, Ube2k.

Computational simulations of bispecific T cell engagers by a multiscale model

The use of bispecific antibodies as T cell engagers can bypass the normal T cell receptor-major histocompatibility class interaction, redirect the cytotoxic activity of T cells, and lead to highly efficient tumor cell killing. However, this immunotherapy also causes significant on-target off-tumor toxicologic effects, especially when it is used to treat solid tumors. To avoid these adverse events, it is necessary to understand the fundamental mechanisms involved in the physical process of T cell engagement. We developed a multiscale computational framework to reach this goal. The framework combines simulations on the intercellular and multicellular levels. On the intercellular level, we simulated the spatial-temporal dynamics of three-body interactions among bispecific antibodies, CD3 and tumor-associated antigens (TAAs). The derived number of intercellular bonds formed between CD3 and TAAs was further transferred to the multicellular simulations as the input parameter of adhesive density between cells. Through the simulations under various molecular and cellular conditions, we were able to gain new insights into how to adopt the most appropriate strategy to maximize the drug efficacy and avoid the off-target effect. For instance, we discovered that the low antibody-binding affinity resulted in the formation of large clusters at the cell-cell interface, which could be important to control the downstream signaling pathways. We also tested different molecular architectures of the bispecific antibody and suggested the existence of an optimal length in regulating the T cell engagement. Overall, the current multiscale simulations serve as a proof-of-concept study to help in the future design of new biological therapeutics.

Recent Advances and Challenges in Protein Structure Prediction

Artificial intelligence has made significant advances in the field of protein structure prediction in recent years. In particular, DeepMind's end-to-end model, AlphaFold2, has demonstrated the capability to predict three-dimensional structures of numerous unknown proteins with accuracy levels comparable to those of experimental methods. This breakthrough has opened up new possibilities for understanding protein structure and function as well as accelerating drug discovery and other applications in the field of biology and medicine. Despite the remarkable achievements of artificial intelligence in the field, there are still some challenges and limitations. In this Review, we discuss the recent progress and some of the challenges in protein structure prediction. These challenges include predicting multidomain protein structures, protein complex structures, multiple conformational states of proteins, and protein folding pathways. Furthermore, we highlight directions in which further improvements can be conducted.

Molecular Dynamics and Other HPC Simulations for Drug Discovery

High performance computing (HPC) is taking an increasingly important place in drug discovery. It makes possible the simulation of complex biochemical systems with high precision in a short time, thanks to the use of sophisticated algorithms. It promotes the advancement of knowledge in fields that are inaccessible or difficult to access through experimentation and it contributes to accelerating the discovery of drugs for unmet medical needs while reducing costs. Herein, we report how computational performance has evolved over the past years, and then we detail three domains where HPC is essential. Molecular dynamics (MD) is commonly used to explore the flexibility of proteins, thus generating a better understanding of different possible approaches to modulate their activity. Modeling and simulation of biopolymer complexes enables the study of protein-protein interactions (PPI) in healthy and disease states, thus helping the identification of targets of pharmacological interest. Virtual screening (VS) also benefits from HPC to predict in a short time, among millions or billions of virtual chemical compounds, the best potential ligands that will be tested in relevant assays to start a rational drug design process.

Accurate ligand-protein docking in CASP15 using the ClusPro LigTBM server

In the ligand prediction category of CASP15, the challenge was to predict the positions and conformations of small molecules binding to proteins that were provided as amino acid sequences or as models generated by the AlphaFold2 program. For most targets, we used our template-based ligand docking program ClusPro ligTBM, also implemented as a public server available at https://ligtbm.cluspro.org/. Since many targets had multiple chains and a number of ligands, several templates, and some manual interventions were required. In a few cases, no templates were found, and we had to use direct docking using the Glide program. Nevertheless, ligTBM was shown to be a very useful tool, and by any ranking criteria, our group was ranked among the top five best-performing teams. In fact, all the best groups used template-based docking methods. Thus, it appears that the AlphaFold2-generated models, despite the high accuracy of the predicted backbone, have local differences from the x-ray structure that make the use of direct docking methods more challenging. The results of CASP15 confirm that this limitation can be frequently overcome by homology-based docking.

TT3D: Leveraging precomputed protein 3D sequence models to predict protein-protein interactions

MOTIVATION

High-quality computational structural models are now precomputed and available for nearly every protein in UniProt. However, the best way to leverage these models to predict which pairs of proteins interact in a high-throughput manner is not immediately clear. The recent Foldseek method of van Kempen et al. encodes the structural information of distances and angles along the protein backbone into a linear string of the same length as the protein string, using tokens from a 21-letter discretized structural alphabet (3Di).

RESULTS

We show that using both the amino acid sequence and the 3Di sequence generated by Foldseek as inputs to our recent deep-learning method, Topsy-Turvy, substantially improves the performance of predicting protein-protein interactions cross-species. Thus TT3D (Topsy-Turvy 3D) presents a way to reuse all the computational effort going into producing high-quality structural models from sequence, while being sufficiently lightweight so that high-quality binary protein-protein interaction predictions across all protein pairs can be made genome-wide.

AVAILABILITY AND IMPLEMENTATION

TT3D is available at https://github.com/samsledje/D-SCRIPT. An archived version of the code at time of submission can be found at https://zenodo.org/records/10037674.

Predictive Modeling and Structure Analysis of Genetic Variants in Familial Hypercholesterolemia: Implications for Diagnosis and Protein Interaction Studies

PURPOSE OF REVIEW

Familial hypercholesterolemia (FH) is a hereditary condition characterized by elevated levels of low-density lipoprotein cholesterol (LDL-C), which increases the risk of cardiovascular disease if left untreated. This review aims to discuss the role of bioinformatics tools in evaluating the pathogenicity of missense variants associated with FH. Specifically, it highlights the use of predictive models based on protein sequence, structure, evolutionary conservation, and other relevant features in identifying genetic variants within LDLR, APOB, and PCSK9 genes that contribute to FH.

RECENT FINDINGS

In recent years, various bioinformatics tools have emerged as valuable resources for analyzing missense variants in FH-related genes. Tools such as REVEL, Varity, and CADD use diverse computational approaches to predict the impact of genetic variants on protein function. These tools consider factors such as sequence conservation, structural alterations, and receptor binding to aid in interpreting the pathogenicity of identified missense variants. While these predictive models offer valuable insights, the accuracy of predictions can vary, especially for proteins with unique characteristics that might not be well represented in the databases used for training. This review emphasizes the significance of utilizing bioinformatics tools for assessing the pathogenicity of FH-associated missense variants. Despite their contributions, a definitive diagnosis of a genetic variant necessitates functional validation through in vitro characterization or cascade screening. This step ensures the precise identification of FH-related variants, leading to more accurate diagnoses. Integrating genetic data with reliable bioinformatics predictions and functional validation can enhance our understanding of the genetic basis of FH, enabling improved diagnosis, risk stratification, and personalized treatment for affected individuals. The comprehensive approach outlined in this review promises to advance the management of this inherited disorder, potentially leading to better health outcomes for those affected by FH.

Computational approaches for innovative anti-viral drug discovery using Orthosiphon aristatus blume miq against dengue virus

Dengue virus has emerged as infectious mosquito borne disease involved in lowering platelets and white blood cells (WBC) count particularly. The genome structure is based on several structural and non-structural proteins essential for viral replication and progeny. One of the major proteins of replication is non-structural protein 3 (NS3) that transforms polyproteins into functional proteins with a cofactor non-structural protein (NS2B). Heat Shock Protein 70 (HSP70), is a human protein that assists in replication, viral entry and virion synthesis. Therefore, to inhibit the spread of dengue infection, there is a need of antivirals targeting replication proteins and other human proteins that help in dengue virus multiplication. By systemic approach based on molecular docking, ADMET (absorption, distribution, metabolism, excretion and toxicity) properties and molecular dynamic simulation (MD), potent inhibitors can be predicted. Inhibition of NS2B/NS3 dengue and HSP70 proteins involved in multiple steps in dengue virus progression can be prevented by using different phytochemicals. Molecular docking was performed using AutoDock Vina, PatchDock, and SwissDock. Interactions of obtained complex were observed in PyMOL and PLIP. Validation was checked by PROCHEK, simulation was performed using iMODS followed by preclinical testing by admetSAR. Ladanein, a flavonoid of Orthosiphon aristatus, was obtained as the lead compound to inhibit major replication protein of dengue virus with inhibitory potential against HSP70 protein. In summary, various in silico approaches were used to obtain the best phytochemical having anti-dengue potential.Communicated by Ramaswamy H. Sarma.

Mapping synthetic binding proteins epitopes on diverse protein targets by protein structure prediction and protein-protein docking

Synthetic binding proteins (SBPs) are a class of artificial proteins engineered from privileged protein scaffolds, which can form highly specific molecular recognition interfaces with a variety of targets. Due to the characteristics of small size, high stability, and good tissue permeability, SBPs have important applications in biomedical research, disease diagnosis and treatment. However, knowledge of SBPs epitopes on the structures of target proteins is still limited, which hinder the development of novel SBPs. In this study, based on the currently available information of SBPs and their targets, 96 pairs of interacting proteins referring to 96 representative SBPs and 80 different targets, were systemically investigated using the state-of-the-art computational modeling techniques including AlphaFold2 protein structure prediction and Rosetta protein-protein docking. As a result, 71 out of the 96 pairs were successfully docked, of which 18, 33, and 20 pairs were defined as models with high, medium, and acceptable quality, respectively. In addition, the interface information was analyzed to decipher the interaction types driven SBPs and targets recognition. Overall, this work not only provides important structural information for understanding the mechanism of action of other SBPs with same protein scaffold, but also for aiding the rational protein engineering and to design of novel SBPs with biomedical applications.

Understanding the General Principles of T Cell Engagement by Multiscale Computational Simulations

The use of bispecific antibodies as T cell engagers can bypass the normal TCR-MHC interaction, redirect the cytotoxic activity of T-cells, and lead to highly efficient tumor cell killing. However, this immunotherapy also causes significant on-target off-tumor toxicologic effects, especially when they were used to treat solid tumors. In order to avoid these adverse events, it is necessary to understand the fundamental mechanisms during the physical process of T cell engagement. We developed a multiscale computational framework to reach this goal. The framework combines simulations on the intercellular and multicellular levels. On the intercellular level, we simulated the spatial-temporal dynamics of three-body interactions among bispecific antibodies, CD3 and TAA. The derived number of intercellular bonds formed between CD3 and TAA were further transferred into the multicellular simulations as the input parameter of adhesive density between cells. Through the simulations under various molecular and cellular conditions, we were able to gain new insights of how to adopt the most appropriate strategy to maximize the drug efficacy and avoid the off-target effect. For instance, we discovered that the low antibody binding affinity resulted in the formation of large clusters at the cell-cell interface, which could be important to control the downstream signaling pathways. We also tested different molecular architectures of the bispecific antibody and suggested the existence of an optimal length in regulating the T cell engagement. Overall, the current multiscale simulations serve as a prove-of-concept study to help the future design of new biological therapeutics.

SIGNIFICANCE

T-cell engagers are a class of anti-cancer drugs that can directly kill tumor cells by bringing T cells next to them. However, current treatments using T-cell engagers can cause serious side-effects. In order to reduce these effects, it is necessary to understand how T cells and tumor cells interact together through the connection of T-cell engagers. Unfortunately, this process is not well studied due to the limitations in current experimental techniques. We developed computational models on two different scales to simulate the physical process of T cell engagement. Our simulation results provide new insights into the general properties of T cell engagers. The new simulation methods can therefore serve as a useful tool to design novel antibodies for cancer immunotherapy.

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.

Molecular dynamic simulations identifying the mechanism of holoenzyme formation by O-GlcNAc transferase and active p38α

O-N-Acetylglucosamine transferase (OGT) can catalyze the O-GlcNAc modification of thousands of proteins. The holoenzyme formation of OGT and adaptor protein is the precondition for further recognition and glycosylation of the target protein, while the corresponding mechanism is still open. Here, static and dynamic schemes based on statistics can successfully screen the feasible identifying, approaching, and binding mechanism of OGT and its typical adaptor protein p38α. The most favorable interface, energy contribution of hotspots, and conformational changes of fragments were discovered. The hydrogen bond interactions were verified as the main driving force for the whole process. The distinct characteristic of active and inactive p38α is explored and demonstrates that the phosphorylated tyrosine and threonine will form strong ion-pair interactions with Lys714, playing a key role in the dynamic identification stage. Multiple method combinations from different points of view may be helpful for exploring other systems of the protein-protein interactions.

Computational approaches for the design of modulators targeting protein-protein interactions

BACKGROUND

Protein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.

AREAS COVERED

In this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.

EXPERT OPINION

Computational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

In silico analysis prediction of HepTH1-5 as a potential therapeutic agent by targeting tumour suppressor protein networks

Many studies reported that the activation of tumour suppressor protein, p53 induced the human hepcidin expression. However, its expression decreased when p53 was silenced in human hepatoma cells. Contrary to Tilapia hepcidin TH1-5, HepTH1-5 was previously reported to trigger the p53 activation through the molecular docking approach. The INhibitor of Growth (ING) family members are also shown to directly interact with p53 and promote cell cycle arrest, senescence, apoptosis and participate in DNA replication and DNA damage responses to suppress the tumour initiation and progression. However, the interrelation between INGs and HepTH1-5 remains unknown. Therefore, this study aims to identify the mechanism and their protein interactions using in silico approaches. The finding revealed that HepTH1-5 and its ligands had interacted mostly on hotspot residues of ING proteins which involved in histone modifications via acetylation, phosphorylation, and methylation. This proves that HepTH1-5 might implicate in an apoptosis signalling pathway and preserve the protein structure and function of INGs by reducing the perturbation of histone binding upon oxidative stress response. This study would provide theoretical guidance for the design and experimental studies to decipher the role of HepTH1-5 as a potential therapeutic agent for cancer therapy. Communicated by Ramaswamy H. Sarma.

Molecular Simulation Study on the Interaction between Porcine CR1-like and C3b

The molecular basis of porcine red blood cell immune adhesion function stems from the complement receptor type 1-like (CR1-like) on its cell membrane. The ligand for CR1-like is C3b, which is produced by the cleavage of complement C3; however, the molecular mechanism of the immune adhesion of porcine erythrocytes is still unclear. Here, homology modeling was used to construct three-dimensional models of C3b and two fragments of CR1-like. An interaction model of C3b-CR1-like was constructed by molecular docking, and molecular structure optimization was achieved using molecular dynamics simulation. A simulated alanine mutation scan revealed that the amino acids Tyr761, Arg763, Phe765, Thr789, and Val873 of CR1-like SCR 12-14 and the amino acid residues Tyr1210, Asn1244, Val1249, Thr1253, Tyr1267, Val1322, and Val1339 of CR1-like SCR 19-21 are key residues involved in the interaction of porcine C3b with CR1-like. This study investigated the interaction between porcine CR1-like and C3b using molecular simulation to clarify the molecular mechanism of the immune adhesion of porcine erythrocytes.

In silico formulation of a next-generation multiepitope vaccine for use as a prophylactic candidate against Crimean-Congo hemorrhagic fever

BACKGROUND

Crimean-Congo hemorrhagic fever (CCHF) is a widespread disease transmitted to humans and livestock animals through the bite of infected ticks or close contact with infected persons' blood, organs, or other bodily fluids. The virus is responsible for severe viral hemorrhagic fever outbreaks, with a case fatality rate of up to 40%. Despite having the highest fatality rate of the virus, a suitable treatment option or vaccination has not been developed yet. Therefore, this study aimed to formulate a multiepitope vaccine against CCHF through computational vaccine design approaches.

METHODS

The glycoprotein, nucleoprotein, and RNA-dependent RNA polymerase of CCHF were utilized to determine immunodominant T- and B-cell epitopes. Subsequently, an integrative computational vaccinology approach was used to formulate a multi-epitopes vaccine candidate against the virus.

RESULTS

After rigorous assessment, a multiepitope vaccine was constructed, which was antigenic, immunogenic, and non-allergenic with desired physicochemical properties. Molecular dynamics (MD) simulations of the vaccine-receptor complex show strong stability of the vaccine candidates to the targeted immune receptor. Additionally, the immune simulation of the vaccine candidates found that the vaccine could trigger real-life-like immune responses upon administration to humans.

CONCLUSIONS

Finally, we concluded that the formulated multiepitope vaccine candidates would provide excellent prophylactic properties against CCHF.

Mapping of antibody epitopes based on docking and homology modeling

Antibodies are key proteins produced by the immune system to target pathogen proteins termed antigens via specific binding to surface regions called epitopes. Given an antigen and the sequence of an antibody the knowledge of the epitope is critical for the discovery and development of antibody based therapeutics. In this work, we present a computational protocol that uses template-based modeling and docking to predict epitope residues. This protocol is implemented in three major steps. First, a template-based modeling approach is used to build the antibody structures. We tested several options, including generation of models using AlphaFold2. Second, each antibody model is docked to the antigen using the fast Fourier transform (FFT) based docking program PIPER. Attention is given to optimally selecting the docking energy parameters depending on the input data. In particular, the van der Waals energy terms are reduced for modeled antibodies relative to x-ray structures. Finally, ranking of antigen surface residues is produced. The ranking relies on the docking results, that is, how often the residue appears in the docking poses' interface, and also on the energy favorability of the docking pose in question. The method, called PIPER-Map, has been tested on a widely used antibody-antigen docking benchmark. The results show that PIPER-Map improves upon the existing epitope prediction methods. An interesting observation is that epitope prediction accuracy starting from antibody sequence alone does not significantly differ from that of starting from unbound (i.e., separately crystallized) antibody structure.

Modeling of Protein Complexes

The recent advances in structural biology, combined with continuously increasing computational capabilities and development of advanced softwares, have drastically simplified the workflow for protein homology modeling. Modeling of individual proteins is nowadays quick and straightforward for a large variety of protein targets, thanks to guided pipelines relying on advanced computational tools and user-friendly interfaces, which have extended and promoted the use of modeling also to scientists not focusing on molecular structures of proteins. Nevertheless, construction of models of multi-protein complexes remains quite challenging for the non-experts, often due to the usage of specific procedures depending on the system under investigation and the need for experimental validation approaches to strengthen the generated output.In this chapter, we provide a brief overview of the approaches enabling generation of multi-protein complex models starting from homology models of individual protein components. Using real-life examples, we include two examples to guide the reader in the generation of homomeric and heteromeric protein models.

The use of predictive models to develop chromatography-based purification processes

Chromatography is the workhorse of biopharmaceutical downstream processing because it can selectively enrich a target product while removing impurities from complex feed streams. This is achieved by exploiting differences in molecular properties, such as size, charge and hydrophobicity (alone or in different combinations). Accordingly, many parameters must be tested during process development in order to maximize product purity and recovery, including resin and ligand types, conductivity, pH, gradient profiles, and the sequence of separation operations. The number of possible experimental conditions quickly becomes unmanageable. Although the range of suitable conditions can be narrowed based on experience, the time and cost of the work remain high even when using high-throughput laboratory automation. In contrast, chromatography modeling using inexpensive, parallelized computer hardware can provide expert knowledge, predicting conditions that achieve high purity and efficient recovery. The prediction of suitable conditions in silico reduces the number of empirical tests required and provides in-depth process understanding, which is recommended by regulatory authorities. In this article, we discuss the benefits and specific challenges of chromatography modeling. We describe the experimental characterization of chromatography devices and settings prior to modeling, such as the determination of column porosity. We also consider the challenges that must be overcome when models are set up and calibrated, including the cross-validation and verification of data-driven and hybrid (combined data-driven and mechanistic) models. This review will therefore support researchers intending to establish a chromatography modeling workflow in their laboratory.

Scoring of protein–protein docking models utilizing predicted interface residues

Scoring docking solutions is a difficult task, and many methods have been developed for this purpose. In docking, only a handful of the hundreds of thousands of models generated by docking algorithms are acceptable, causing difficulties when developing scoring functions. Today's best scoring functions can significantly increase the number of top‐ranked models but still fail for most targets. Here, we examine the possibility of utilizing predicted interface residues to score docking models generated during the scan stage of a docking algorithm. Many methods have been developed to infer the regions of a protein surface that interact with another protein, but most have not been benchmarked using docking algorithms. This study systematically tests different interface prediction methods for scoring >300.000 low‐resolution rigid‐body template free docking decoys. Overall we find that contact‐based interface prediction by BIPSPI is the best method to score docking solutions, with >12% of first ranked docking models being acceptable. Additional experiments indicated precision as a high‐importance metric when estimating interface prediction quality, focusing on docking constraints production. Finally, we discussed several limitations for adopting interface predictions as constraints in a docking protocol.

Protein complex structure prediction powered by multiple sequence alignments of interologs from multiple taxonomic ranks and AlphaFold2

AlphaFold2 can predict protein complex structures as long as a multiple sequence alignment (MSA) of the interologs of the target protein-protein interaction (PPI) can be provided. In this study, a simplified phylogeny-based approach was applied to generate the MSA of interologs, which was then used as the input to AlphaFold2 for protein complex structure prediction. In this extensively benchmarked protocol on nonredundant PPI dataset, including 107 bacterial PPIs and 442 eukaryotic PPIs, we show complex structures of 79.5% of the bacterial PPIs and 49.8% of the eukaryotic PPIs can be successfully predicted, which yielded significantly better performance than the application of MSA of interologs prepared by two existing approaches. Considering PPIs may not be conserved in species with long evolutionary distances, we further restricted interologs in the MSA to different taxonomic ranks of the species of the target PPI in protein complex structure prediction. We found that the success rates can be increased to 87.9% for the bacterial PPIs and 56.3% for the eukaryotic PPIs if interologs in the MSA are restricted to a specific taxonomic rank of the species of each target PPI. Finally, we show that the optimal taxonomic ranks for protein complex structure prediction can be selected with the application of the predicted template modeling (TM) scores of the output models.

Characterization and molecular docking study of cathepsin L inhibitory peptides (SnuCalCpIs) from Calotropis procera R. Br

Propeptides, released from the autocatalytic activation of its zymogen, are potential inhibitors against proteases involved in cancer cell invasion and migration. Our research team previously obtained novel propeptides (SnuCalCpIs) from transcriptome analysis of the medicinal plant Calotropis procera R. Br. and reported them as promising candidates for cancer therapeutics due to their cathepsin L inhibition activity. In the present study, inhibitory activity among SnuCalCpIs was compared with inhibition efficiency and verified by in silico molecular docking analysis. Only SnuCalCpI03 and SnuCalCpI15, expressed in Escherichia coli, showed inhibitory activity against cathepsin L as competitive inhibitors, and the half-maximal inhibitory concentrations (IC₅₀) values of 2.1 nM and 1.6 nM, respectively. They were stable below 70 °C, maintaining more than 90% inhibitory activity over a wide range of pH (2.0-10.0), except at the isoelectric point (pI). The template-based docking simulation models showed that SnuCalCpI02, SnuCalCpI12, and SnuCalCpI16 could not interact with the substrate-binding cleft of cathepsin L even though they possessed the same conserved domain. In contrast, SnuCalCpI03 and SnuCalCpI15 interacted with cathepsin L along the propeptide binding loop and substrate-binding cleft, resulting in obstruction of substrate access to the active site.

Mechanism of Zhen Wu Decoction in the Treatment of Heart Failure Based on Network Pharmacology and Molecular Docking

Heart failure (HF) is a serious manifestation or advanced stage of various cardiovascular diseases, and its mortality and rehospitalization rate are still on the rise in China. Based on the network pharmacology method, 59 components of Zhen Wu decoction (ZWD) and 83 target genes related to HF were obtained. Through the PPI network, four potential therapeutic targets were identified: AKT1, IL6, JUN, and MAPK8. The beneficial components of ZWD might intervene HF through the AGE-RAGE signalling pathway in the diabetes component, fluid shear stress and atherosclerosis, the TNF signalling pathway, TB, and Kaposi sarcoma related herpesvirus infection, according to a KEGG enrichment study. The protein interaction network of candidate targets was constructed by the STRING database, and the protein interaction network was clustered by MEODE software. GO and KEGG enrichment analyses were performed on the core modules obtained by clustering. Finally, AutoDock Vina software was used for molecular docking verification of key targets and active ingredients. The result was that 75 active ingredients and 109 genes were screened as potential active ingredients and potential targets of Shengjie Tongyu decoction for CHF treatment. The main active components were quercetin, luteolin, kaempferol, dehydrated icariin, isorhamnetin, formononetin, and other flavonoids. Il-6, MAPK1, MAPK8, AKT1, VEGFA, and JUN were selected as the core targets. Molecular docking showed that the key components were well connected with the target. GO enrichment analysis showed that Shengjie Tongyu decoction could play a role through multiple biological pathways including angiogenesis, regulation of endothelial cell proliferation, binding of cytokine receptors, negative regulation of apoptotic signalling pathways, regulation of nitric oxide synthase activity, and reactive oxygen metabolism. Key pathways mainly focus on the toll-like receptor signalling pathway, nod-like receptor signalling pathway, MAPK signalling pathway, mTOR signalling pathway, JAK-STAT signalling pathway, VEGF signalling pathway, and other pathways. Through molecular docking technology, it was found that a variety of effective components in ZWD, such as kaempferol. Molecular docking technology has preliminatively verified the network pharmacology and laid a foundation for the follow-up pharmacological research.

Cross-linking and mass spectrometry as a tool for studying the structural biology of ribonucleoproteins

Cross-linking and mass spectrometry (XL-MS) workflows represent an increasingly popular technique for low-resolution structural studies of macromolecular complexes. Cross-linking reactions take place in the solution state, capturing contact sites between components of a complex that represent the native, functionally relevant structure. Protein-protein XL-MS protocols are widely adopted, providing precise localization of cross-linking sites to single amino acid positions within a pair of cross-linked peptides. In contrast, protein-RNA XL-MS workflows are evolving rapidly and differ in their ability to localize interaction regions within the RNA sequence. Here, we review protein-protein and protein-RNA XL-MS workflows, and discuss their applications in studies of protein-RNA complexes. The examples highlight the complementary value of XL-MS in structural studies of protein-RNA complexes, where more established high-resolution techniques might be unable to produce conclusive data.

Improved prediction of protein-protein interactions using AlphaFold2

Predicting the structure of interacting protein chains is a fundamental step towards understanding protein function. Unfortunately, no computational method can produce accurate structures of protein complexes. AlphaFold2, has shown unprecedented levels of accuracy in modelling single chain protein structures. Here, we apply AlphaFold2 for the prediction of heterodimeric protein complexes. We find that the AlphaFold2 protocol together with optimised multiple sequence alignments, generate models with acceptable quality (DockQ ≥ 0.23) for 63% of the dimers. From the predicted interfaces we create a simple function to predict the DockQ score which distinguishes acceptable from incorrect models as well as interacting from non-interacting proteins with state-of-art accuracy. We find that, using the predicted DockQ scores, we can identify 51% of all interacting pairs at 1% FPR.

Predicting the structure of protein complexes is extremely difficult. Here, authors apply AlphaFold2 with optimized multiple sequence alignments to model complexes of interacting proteins, enabling prediction of both if and how proteins interact with state-of-art accuracy.

Challenges and Advantages of Accounting for Backbone Flexibility in Prediction of Protein-Protein Complexes

Predicting protein binding is a core problem of computational biophysics. That this objective can be partly achieved with some amount of success using docking algorithms based on rigid protein models is remarkable, although going further requires allowing for protein flexibility. However, accurately capturing the conformational changes upon binding remains an enduring challenge for docking algorithms. Here, we adapt our Upside folding model, where side chains are represented as multi-position beads, to explore how flexibility may impact predictions of protein-protein complexes. Specifically, the Upside model is used to investigate where backbone flexibility helps, which types of interactions are important, and what is the impact of coarse graining. These efforts also shed light on the relative challenges posed by folding and docking. After training the Upside energy function for docking, the model is competitive with the established all-atom methods. However, allowing for backbone flexibility during docking is generally detrimental, as the presence of comparatively minor (3-5 Å) deviations relative to the docked structure has a large negative effect on performance. While this issue appears to be inherent to current forcefield-guided flexible docking methods, systems involving the co-folding of flexible loops such as antibody-antigen complexes represent an interesting exception. In this case, binding is improved when backbone flexibility is allowed using the Upside model.

TRScore: a three-dimensional RepVGG-based scoring method for ranking protein docking models

MOTIVATION

Protein-protein interactions (PPI) play important roles in cellular activities. Due to the technical difficulty and high cost of experimental methods, there are considerable interests towards the development of computational approaches, such as protein docking, to decipher PPI patterns. One of the important and difficult aspects in protein docking is recognizing near-native conformations from a set of decoys, but unfortunately traditional scoring functions still suffer from limited accuracy. Therefore, new scoring methods are pressingly needed in methodological and/or practical implications.

RESULTS

We present a new deep learning-based scoring method for ranking protein-protein docking models based on a three-dimensional (3D) RepVGG network, named TRScore. To recognize near-native conformations from a set of decoys, TRScore voxelizes the protein-protein interface into a 3D grid labeled by the number of atoms in different physicochemical classes. Benefiting from the deep convolutional RepVGG architecture, TRScore can effectively capture the subtle differences between energetically favorable near-native models and unfavorable non-native decoys without needing extra information. TRScore was extensively evaluated on diverse test sets including protein-protein docking benchmark 5.0 update set, DockGround decoy set, as well as realistic CAPRI decoy set, and overall obtained a significant improvement over existing methods in cross validation and independent evaluations.

AVAILABILITY

Codes available at: https://github.com/BioinformaticsCSU/TRScore.

A strategy based on gene sequencing and molecular docking for analysis and prediction of bioactive peptides in Shuxuetong injection

Peptides are a class of protein fragments with relatively high biological activity and intense specificity, which play crucial role in the treatment of Shuxuetong injection (SXT). However, the extraordinary complexity of Chinese medicinal formulates and the lack of systematic identification methods are primary challenges for study of pharmacodynamic peptides. In addition, infinitesimal peptides contents further hinder the identification and structural characterization of polypeptide by traditional means. In this paper, we described a strategy that LC-MS combined with molecular docking to systematically illustrate the peptide components of SXT. The key to this research was used of gene sequencing to establish a SXT protein database to further achieve the separation and enrichment of chemical methods. Moreover, the ADRA2A, PAR4 and DRD3 were precisely docked with the identified peptides. The result indicated that 12 compounds had stable binding ability and were speculated to be the latent bioactive monomers for the treatment of stroke. Additionally, 12 peptides were verified by cell-based experiment. The results showed that only YLKTT could indeed protect astrocytes from oxygen glucose deprivation/reoxygenation (OGD/R). The YLKTT showed higher activity than the others in vitro. It might be a completely new compound that has never been reported before, providing the basis for further research and a new paradigm for stroke.

Immune Mechanism, Gene Module, and Molecular Subtype Identification of Astragalus Membranaceus in the Treatment of Dilated Cardiomyopathy: An Integrated Bioinformatics Study

Astragalus membranaceus has complex components as a natural drug and has multilevel, multitarget, and multichannel effects on dilated cardiomyopathy (DCM). However, the immune mechanism, gene module, and molecular subtype of astragalus membranaceus in the treatment of DCM are still not revealed. Microarray information of GSE84796 was downloaded from the GEO database, including RNA sequencing data of seven normal cardiac tissues and ten DCM cardiac tissues. A total of 4029 DCM differentially expressed genes were obtained, including 1855 upregulated genes and 2174 downregulated genes. GO/KEGG/GSEA analysis suggested that the activation of T cells and B cells was the primary cause of DCM. WGCNA was used to obtain blue module genes. The blue module genes are primarily ADCY7, BANK1, CD1E, CD19, CD38, CD300LF, CLEC4E, FLT3, GPR18, HCAR3, IRF4, LAMP3, MRC1, SYK, and TLR8, which successfully divided DCM into three molecular subtypes. Based on the CIBERSORT algorithm, the immune infiltration profile of DCM was analyzed. Many immune cell subtypes, including the abovementioned immune cells, showed different levels of increased infiltration in the myocardial tissue of DCM. However, this infiltration pattern was not obviously correlated with clinical characteristics, such as age, EF, and sex. Based on network pharmacology and ClueGO, 20 active components of Astragalus membranaceus and 40 components of DMCTGS were obtained from TCMSP. Through analysis of the immune regulatory network, we found that Astragalus membranaceus effectively regulates the activation of immune cells, such as B cells and T cells, cytokine secretion, and other processes and can intervene in DCM at multiple components, targets, and levels. The above mechanisms were verified by molecular docking results, which confirmed that AKT1, VEGFA, MMP9, and RELA are promising potential targets of DCM.

Integration of Mass Spectrometry Data for Structural Biology

Mass spectrometry (MS) is increasingly being used to probe the structure and dynamics of proteins and the complexes they form with other macromolecules. There are now several specialized MS methods, each with unique sample preparation, data acquisition, and data processing protocols. Collectively, these methods are referred to as structural MS and include cross-linking, hydrogen-deuterium exchange, hydroxyl radical footprinting, native, ion mobility, and top-down MS. Each of these provides a unique type of structural information, ranging from composition and stoichiometry through to residue level proximity and solvent accessibility. Structural MS has proved particularly beneficial in studying protein classes for which analysis by classic structural biology techniques proves challenging such as glycosylated or intrinsically disordered proteins. To capture the structural details for a particular system, especially larger multiprotein complexes, more than one structural MS method with other structural and biophysical techniques is often required. Key to integrating these diverse data are computational strategies and software solutions to facilitate this process. We provide a background to the structural MS methods and briefly summarize other structural methods and how these are combined with MS. We then describe current state of the art approaches for the integration of structural MS data for structural biology. We quantify how often these methods are used together and provide examples where such combinations have been fruitful. To illustrate the power of integrative approaches, we discuss progress in solving the structures of the proteasome and the nuclear pore complex. We also discuss how information from structural MS, particularly pertaining to protein dynamics, is not currently utilized in integrative workflows and how such information can provide a more accurate picture of the systems studied. We conclude by discussing new developments in the MS and computational fields that will further enable in-cell structural studies.

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.

GalaxyHeteromer: protein heterodimer structure prediction by template-based and ab initio docking

Protein–protein interactions play crucial roles in diverse biological processes, including various disease progressions. Atomistic structural details of protein–protein interactions may provide important information that can facilitate the design of therapeutic agents. GalaxyHeteromer is a freely available automatic web server (http://galaxy.seoklab.org/heteromer) that predicts protein heterodimer complex structures from two subunit protein sequences or structures. When subunit structures are unavailable, they are predicted by template- or distance-prediction-based modelling methods. Heterodimer complex structures can be predicted by both template-based and ab initio docking, depending on the template's availability. Structural templates are detected from the protein structure database based on both the sequence and structure similarities. The templates for heterodimers may be selected from monomer and homo-oligomer structures, as well as from hetero-oligomers, owing to the evolutionary relationships of heterodimers with domains of monomers or subunits of homo-oligomers. In addition, the server employs one of the best ab initio docking methods when heterodimer templates are unavailable. The multiple heterodimer structure models and the associated scores, which are provided by the web server, may be further examined by user to test or develop functional hypotheses or to design new functional molecules.

Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations

Most of the protein–protein docking methods treat proteins as almost rigid objects. Only the side-chains flexibility is usually taken into account. The few approaches enabling docking with a flexible backbone typically work in two steps, in which the search for protein–protein orientations and structure flexibility are simulated separately. In this work, we propose a new straightforward approach for docking sampling. It consists of a single simulation step during which a protein undergoes large-scale backbone rearrangements, rotations, and translations. Simultaneously, the other protein exhibits small backbone fluctuations. Such extensive sampling was possible using the CABS coarse-grained protein model and Replica Exchange Monte Carlo dynamics at a reasonable computational cost. In our proof-of-concept simulations of 62 protein–protein complexes, we obtained acceptable quality models for a significant number of cases.

Modeling of protein complexes in CASP14 with emphasis on the interaction interface prediction

The goal of CASP experiments is to monitor the progress in the protein structure prediction field. During the 14th CASP edition we aimed to test our capabilities of predicting structures of protein complexes. Our protocol for modeling protein assemblies included both template‐based modeling and free docking. Structural templates were identified using sensitive sequence‐based searches. If sequence‐based searches failed, we performed structure‐based template searches using selected CASP server models. In the absence of reliable templates we applied free docking starting from monomers generated by CASP servers. We evaluated and ranked models of protein complexes using an improved version of our protein structure quality assessment method, VoroMQA, taking into account both interaction interface and global structure scores. If reliable templates could be identified, generally accurate models of protein assemblies were generated with the exception of an antibody‐antigen interaction. The success of free docking mainly depended on the accuracy of initial subunit models and on the scoring of docking solutions. To put our overall results in perspective, we analyzed our performance in the context of other CASP groups. Although the subunits in our assembly models often were not of the top quality, these models had, overall, the best‐predicted intersubunit interfaces according to several accuracy measures. We attribute our relative success primarily to the emphasis on the interaction interface when modeling and scoring.

Use of molecular docking computational tools in drug discovery

Molecular docking has become an important component of the drug discovery process. Since first being developed in the 1980s, advancements in the power of computer hardware and the increasing number of and ease of access to small molecule and protein structures have contributed to the development of improved methods, making docking more popular in both industrial and academic settings. Over the years, the modalities by which docking is used to assist the different tasks of drug discovery have changed. Although initially developed and used as a standalone method, docking is now mostly employed in combination with other computational approaches within integrated workflows. Despite its invaluable contribution to the drug discovery process, molecular docking is still far from perfect. In this chapter we will provide an introduction to molecular docking and to the different docking procedures with a focus on several considerations and protocols, including protonation states, active site waters and consensus, that can greatly improve the docking results.

Hotspot Coevolution Is a Key Identifier of Near-Native Protein Complexes

Protein-protein interactions play a key role in mediating numerous biological functions, with more than half the proteins in living organisms existing as either homo- or hetero-oligomeric assemblies. Protein subunits that form oligomers minimize the free energy of the complex, but exhaustive computational search-based docking methods have not comprehensively addressed the challenge of distinguishing a natively bound complex from non-native forms. Current protein docking approaches address this problem by sampling multiple binding modes in proteins and scoring each mode, with the lowest-energy (or highest scoring) binding mode being regarded as a near-native complex. However, high-scoring modes often match poorly with the true bound form, suggesting a need for improvement of the scoring function. In this study, we propose a scoring function, KFC-E, that accounts for both conservation and coevolution of putative binding hotspot residues at protein-protein interfaces. We tested KFC-E on four benchmark sets of unbound examples and two benchmark sets of bound examples, with the results demonstrating a clear improvement over scores that examine conservation and coevolution across the entire interface.

Exploring the Binding Interaction of Raf Kinase Inhibitory Protein With the N-Terminal of C-Raf Through Molecular Docking and Molecular Dynamics Simulation

Protein-protein interactions are indispensable physiological processes regulating several biological functions. Despite the availability of structural information on protein-protein complexes, deciphering their complex topology remains an outstanding challenge. Raf kinase inhibitory protein (RKIP) has gained substantial attention as a favorable molecular target for numerous pathologies including cancer and Alzheimer’s disease. RKIP interferes with the RAF/MEK/ERK signaling cascade by endogenously binding with C-Raf (Raf-1 kinase) and preventing its activation. In the current investigation, the binding of RKIP with C-Raf was explored by knowledge-based protein-protein docking web-servers including HADDOCK and ZDOCK and a consensus binding mode of C-Raf/RKIP structural complex was obtained. Molecular dynamics (MD) simulations were further performed in an explicit solvent to sample the conformations for when RKIP binds to C-Raf. Some of the conserved interface residues were mutated to alanine, phenylalanine and leucine and the impact of mutations was estimated by additional MD simulations and MM/PBSA analysis for the wild-type (WT) and constructed mutant complexes. Substantial decrease in binding free energy was observed for the mutant complexes as compared to the binding free energy of WT C-Raf/RKIP structural complex. Furthermore, a considerable increase in average backbone root mean square deviation and fluctuation was perceived for the mutant complexes. Moreover, per-residue energy contribution analysis of the equilibrated simulation trajectory by HawkDock and ANCHOR web-servers was conducted to characterize the key residues for the complex formation. One residue each from C-Raf (Arg398) and RKIP (Lys80) were identified as the druggable “hot spots” constituting the core of the binding interface and corroborated by additional long-time scale (300 ns) MD simulation of Arg398Ala mutant complex. A notable conformational change in Arg398Ala mutant occurred near the mutation site as compared to the equilibrated C-Raf/RKIP native state conformation and an essential hydrogen bonding interaction was lost. The thirteen binding sites assimilated from the overall analysis were mapped onto the complex as surface and divided into active and allosteric binding sites, depending on their location at the interface. The acquired information on the predicted 3D structural complex and the detected sites aid as promising targets in designing novel inhibitors to block the C-Raf/RKIP interaction.

Computational Modeling as a Tool to Investigate PPI: From Drug Design to Tissue Engineering

Protein-protein interactions (PPIs) mediate a large number of important regulatory pathways. Their modulation represents an important strategy for discovering novel therapeutic agents. However, the features of PPI binding surfaces make the use of structure-based drug discovery methods very challenging. Among the diverse approaches used in the literature to tackle the problem, linear peptides have demonstrated to be a suitable methodology to discover PPI disruptors. Unfortunately, the poor pharmacokinetic properties of linear peptides prevent their direct use as drugs. However, they can be used as models to design enzyme resistant analogs including, cyclic peptides, peptide surrogates or peptidomimetics. Small molecules have a narrower set of targets they can bind to, but the screening technology based on virtual docking is robust and well tested, adding to the computational tools used to disrupt PPI. We review computational approaches used to understand and modulate PPI and highlight applications in a few case studies involved in physiological processes such as cell growth, apoptosis and intercellular communication.

In silico Approaches for the Design and Optimization of Interfering Peptides Against Protein-Protein Interactions

Large contact surfaces of protein-protein interactions (PPIs) remain to be an ongoing issue in the discovery and design of small molecule modulators. Peptides are intrinsically capable of exploring larger surfaces, stable, and bioavailable, and therefore bear a high therapeutic value in the treatment of various diseases, including cancer, infectious diseases, and neurodegenerative diseases. Given these promising properties, a long way has been covered in the field of targeting PPIs via peptide design strategies. In silico tools have recently become an inevitable approach for the design and optimization of these interfering peptides. Various algorithms have been developed to scrutinize the PPI interfaces. Moreover, different databases and software tools have been created to predict the peptide structures and their interactions with target protein complexes. High-throughput screening of large peptide libraries against PPIs; "hotspot" identification; structure-based and off-structure approaches of peptide design; 3D peptide modeling; peptide optimization strategies like cyclization; and peptide binding energy evaluation are among the capabilities of in silico tools. In the present study, the most recent advances in the field of in silico approaches for the design of interfering peptides against PPIs will be reviewed. The future perspective of the field and its advantages and limitations will also be pinpointed.

Atomic-level evolutionary information improves protein-protein interface scoring

MOTIVATION

The crucial role of protein interactions and the difficulty in characterising them experimentally strongly motivates the development of computational approaches for structural prediction. Even when protein-protein docking samples correct models, current scoring functions struggle to discriminate them from incorrect decoys. The previous incorporation of conservation and coevolution information has shown promise for improving protein-protein scoring. Here, we present a novel strategy to integrate atomic-level evolutionary information into different types of scoring functions to improve their docking discrimination.

RESULTS

: We applied this general strategy to our residue-level statistical potential from InterEvScore and to two atomic-level scores, SOAP-PP and Rosetta interface score (ISC). Including evolutionary information from as few as ten homologous sequences improves the top 10 success rates of individual atomic-level scores SOAP-PP and Rosetta ISC by respectively 6 and 13.5 percentage points, on a large benchmark of 752 docking cases. The best individual homology-enriched score reaches a top 10 success rate of 34.4%. A consensus approach based on the complementarity between different homology-enriched scores further increases the top 10 success rate to 40%.

AVAILABILITY

All data used for benchmarking and scoring results, as well as a Singularity container of the pipeline, are available at http://biodev.cea.fr/interevol/interevdata/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Conservation of binding properties in protein models

We study the models submitted to round 12 of the Critical Assessment of protein Structure Prediction (CASP) experiment to assess how well the binding properties are conserved when the X-ray structures of the target proteins are replaced by their models. To explore small molecule binding we generate distributions of molecular probes - which are fragment-sized organic molecules of varying size, shape, and polarity - around the protein, and count the number of interactions between each residue and the probes, resulting in a vector of interactions we call a binding fingerprint. The similarity between two fingerprints, one for the X-ray structure and the other for a model of the protein, is determined by calculating the correlation coefficient between the two vectors. The resulting correlation coefficients are shown to correlate with global measures of accuracy established in CASP, and the relationship yields an accuracy threshold that has to be reached for meaningful binding surface conservation. The clusters formed by the probe molecules reliably predict binding hot spots and ligand binding sites in both X-ray structures and reasonably accurate models of the target, but ensembles of models may be needed for assessing the availability of proper binding pockets. We explored ligand docking to the few targets that had bound ligands in the X-ray structure. More targets were available to assess the ability of the models to reproduce protein-protein interactions by docking both the X-ray structures and models to their interaction partners in complexes. It was shown that this application is more difficult than finding small ligand binding sites, and the success rates heavily depend on the local structure in the potential interface. In particular, predicted conformations of flexible loops are frequently incorrect in otherwise highly accurate models, and may prevent predicting correct protein-protein interactions.

In-Silico Drug Designing of Spike Receptor with Its ACE2 Receptor and Nsp10/Nsp16 MTase Complex Against SARS-CoV-2

The realm Riboviria constitutes Coronaviruses, which led to the emergence of the pandemic COVID 19 in the twenty-first century affected millions of lives. At present, the management of COVID 19 largely depends on antiviral therapeutics along with the anti-inflammatory drug. The vaccine is under the final clinical phase, and emergency use is available. We aim at ACE2 and Nsp10/Nsp16 MTase as potential drug candidate in COVID 19 management in the present work. For drug designing, various computational simulation strategies have been employed like Swiss-Model, Hawk Dock, HDOCK, py Dock, and PockDrug for homology modeling, binding energies of the molecule with a target, simulate the conformation and binding poses, statistics of protein lock with target key and drug ability, respectively. The current in-silico screening depicts that the spike protein receptor is complementary to the target when bound to each other and forms a stable complex. The MMGBSA free energy binding property of receptor and ligand is critical. The intermolecular Statistics with the target Nsp10/Nsp16 MTase complex are plausible. We have also observed a high-affinity pocket binding site with the target. Therefore, the favorable intermolecular interactions and Physico-chemical properties emanate as a drug candidate treating COVID-19. This study has approached computational tools to analyze the conformation, binding affinity, and drug ability of receptor-ligand. Thus, the spike receptor with its ACE2 receptor with Nsp10/Nsp16 MTase complex would be a potent drug against SARS CoV-2 and can cure the infection as per consensus scoring.

Testing methods and statistical models of genomic prediction for quantitative disease resistance to Phytophthora sojae in soybean [Glycine max (L.) Merr] germplasm collections

KEY MESSAGE

Genomic prediction of quantitative resistance toward Phytophthora sojae indicated that genomic selection may increase breeding efficiency. Statistical model and marker set had minimal effect on genomic prediction with > 1000 markers. Quantitative disease resistance (QDR) toward Phytophthora sojae in soybean is a complex trait controlled by many small-effect loci throughout the genome. Along with the technical and rate-limiting challenges of phenotyping resistance to a root pathogen, the trait complexity can limit breeding efficiency. However, the application of genomic prediction to traits with complex genetic architecture, such as QDR toward P. sojae, is likely to improve breeding efficiency. We provide a novel example of genomic prediction by measuring QDR to P. sojae in two diverse panels of more than 450 plant introductions (PIs) that had previously been genotyped with the SoySNP50K chip. This research was completed in a collection of diverse germplasm and contributes to both an initial assessment of genomic prediction performance and characterization of the soybean germplasm collection. We tested six statistical models used for genomic prediction including Bayesian Ridge Regression; Bayesian LASSO; Bayes A, B, C; and reproducing kernel Hilbert spaces. We also tested how the number and distribution of SNPs included in genomic prediction altered predictive ability by varying the number of markers from less than 50 to more than 34,000 SNPs, including SNPs based on sequential sampling, random sampling, or selections from association analyses. Predictive ability was relatively independent of statistical model and marker distribution, with a diminishing return when more than 1000 SNPs were included in genomic prediction. This work estimated relative efficiency per breeding cycle between 0.57 and 0.83, which may improve the genetic gain for P. sojae QDR in soybean breeding programs.

New insights on human IRE1 tetramer structures based on molecular modeling

Inositol-Requiring Enzyme 1α (IRE1α; hereafter IRE1) is a transmembrane kinase/ribonuclease protein related with the unfolded protein response (UPR) signaling. Experimental evidence suggests that IRE1 forms several three dimensional (3D) structural variants: dimers, tetramers and higher order oligomers, where each structural variant can contain different IRE1 conformers in different arrangements. For example, studies have shown that two sets of IRE1 dimers exist; a face-to-face dimer and a back-to-back dimer, with the latter considered the important unit for UPR signaling propagation. However, the structural configuration and mechanistic details of the biologically important IRE1 tetramers are limited. Here, we combine protein–protein docking with molecular dynamics simulations to derive human IRE1 tetramer models and identify a molecular mechanism of IRE1 activation. To validate the derived models of the human IRE1 tetramer, we compare the dynamic behavior of the models with the yeast IRE1 tetramer crystallographic structure. We show that IRE1 tetramer conformational changes could be linked to the initiation of the unconventional splicing of mRNA encoding X-box binding protein-1 (XBP1), which allows for the expression of the transcription factor XBP1s (XBP1 spliced). The derived IRE1 tetrameric models bring new mechanistic insights about the IRE1 molecular activation mechanism by describing the IRE1 tetramers as active protagonists accommodating the XBP1 substrate.

Computational structure modeling for diverse categories of macromolecular interactions

Computational protein-protein docking is one of the most intensively studied topics in structural bioinformatics. The field has made substantial progress through over three decades of development. The development began with methods for rigid-body docking of two proteins, which have now been extended in different directions to cover the various macromolecular interactions observed in a cell. Here, we overview the recent developments of the variations of docking methods, including multiple protein docking, peptide-protein docking, and disordered protein docking methods.

Performance and Its Limits in Rigid Body Protein-Protein Docking

The development of fast Fourier transform (FFT) algorithms enabled the sampling of billions of complex conformations and thus revolutionized protein-protein docking. FFT-based methods are now widely available and have been used in hundreds of thousands of docking calculations. Although the methods perform "soft" docking, which allows for some overlap of component proteins, the rigid body assumption clearly introduces limitations on accuracy and reliability. In addition, the method can work only with energy expressions represented by sums of correlation functions. In this paper we use a well-established protein-protein docking benchmark set to evaluate the results of these limitations by focusing on the performance of the docking server ClusPro, which implements one of the best rigid body methods. Furthermore, we explore the theoretical limits of accuracy when using established energy terms for scoring, provide comparison with flexible docking algorithms, and review the historical performance of servers in the CAPRI docking experiment.