Protein-protein docking: from interaction to interactome

Martínez-Castilla L, García-Sandoval R, Ballester P, Ferrándiz C, Sanchez MDLP, García-Ponce B, Garay-Arroyo A, R. Álvarez-Buylla E. HDACs MADS-domain protein interaction: a case study of HDA15 and XAL1 in Arabidopsis thaliana

Cellular behavior, cell differentiation and ontogenetic development in eukaryotes result from complex interactions between epigenetic and classic molecular genetic mechanisms, with many of these interactions still to be elucidated. Histone deacetylase enzymes (HDACs) promote the interaction of histones with DNA by compacting the nucleosome, thus causing transcriptional repression. MADS-domain transcription factors are highly conserved in eukaryotes and participate in controlling diverse developmental processes in animals and plants, as well as regulating stress responses in plants. In this work, we focused on finding out putative interactions of Arabidopsis thaliana HDACs and MADS-domain proteins using an evolutionary perspective combined with bioinformatics analyses and testing the more promising predicted interactions through classic molecular biology tools. Through bioinformatic analyses, we found similarities between HDACs proteins from different organisms, which allowed us to predict a putative protein-protein interaction between the Arabidopsis thaliana deacetylase HDA15 and the MADS-domain protein XAANTAL1 (XAL1). The results of two-hybrid and Bimolecular Fluorescence Complementation analysis demonstrated in vitro and in vivo HDA15-XAL1 interaction in the nucleus. Likely, this interaction might regulate developmental processes in plants as is the case for this type of interaction in animals.

SurfPro-NN: A 3D point cloud neural network for the scoring of protein-protein docking models based on surfaces features and protein language models

Protein-protein interactions (PPI) play a crucial role in numerous key biological processes, and the structure of protein complexes provides valuable clues for in-depth exploration of molecular-level biological processes. Protein-protein docking technology is widely used to simulate the spatial structure of proteins. However, there are still challenges in selecting candidate decoys that closely resemble the native structure from protein-protein docking simulations. In this study, we introduce a docking evaluation method based on three-dimensional point cloud neural networks named SurfPro-NN, which represents protein structures as point clouds and learns interaction information from protein interfaces by applying a point cloud neural network. With the continuous advancement of deep learning in the field of biology, a series of knowledge-rich pre-trained models have emerged. We incorporate protein surface representation models and language models into our approach, greatly enhancing feature representation capabilities and achieving superior performance in protein docking model scoring tasks. Through comprehensive testing on public datasets, we find that our method outperforms state-of-the-art deep learning approaches in protein-protein docking model scoring. Not only does it significantly improve performance, but it also greatly accelerates training speed. This study demonstrates the potential of our approach in addressing protein interaction assessment problems, providing strong support for future research and applications in the field of biology.

ZEPPI: Proteome-scale sequence-based evaluation of protein-protein interaction models

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence coevolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI can be implemented on a proteome-wide scale and is applied here to millions of structural models of dimeric complexes in the Escherichia coli and human interactomes found in the PrePPI database. PrePPI's scoring function is based primarily on the evaluation of protein-protein interfaces, and ZEPPI adds a new feature to this analysis through the incorporation of evolutionary information. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. As we discuss, the standard CAPRI scores used to evaluate docking models are based on model quality and not on the ability to give yes/no answers as to whether two proteins interact. ZEPPI is able to detect weak signals from PPI models that the CAPRI scores define as incorrect and, similarly, to identify potential PPIs defined as low confidence by the current PrePPI scoring function. A number of examples that illustrate how the combination of PrePPI and ZEPPI can yield functional hypotheses are provided.

Functionalized siRNA-chitosan nanoformulations promote triple-negative breast cancer cell death via blocking the miRNA-21/AKT/ERK signaling axis: in-silico and in vitro studies

Oncogenic microRNA (miRNA), especially miRNA-21 upregulation in triple-negative breast cancer (TNBC), suggests a new class of therapeutic targets. In this study, we aimed to create GE11 peptide-conjugated small interfering RNA-loaded chitosan nanoparticles (GE11-siRNA-CSNPs) for the targeting of EGFR overexpressed TNBC and selectively inhibit miRNA-21 expression. A variety of in-silico and in vitro cellular and molecular studies were conducted to investigate the binding affinities of specific targets used as well as the anticancer efficacies and mechanisms of GE11-siRNA-CSNPs in TNBC cells. An in-silico assessment reveals a distinct binding affinity of miRNA-21 with siRNA as well as between the extracellular domain of EGFR and synthesized peptides. Notably, the in vitro results showed that GE11-siRNA-CSNPs were revealed to have better cytotoxicity against TNBC cells. It significantly inhibits miRNA-21 expression, cell migration, and colony formation. The results also indicated that GE11-siRNA-CSNPs impeded cell cycle progression. It induces cell death by reducing the expression of the antiapoptotic gene Bcl-2 and increasing the expression of the proapoptotic genes Bax, Caspase 3, and Caspase 9. Additionally, the docking analysis and immunoblot investigations verified that GE1-siRNA-CSNPs, which specifically target TNBC cells and suppress miRNA-21, can prevent the effects of miRNA-21 on the proliferation of TNBC cells via controlling EGFR and subsequently inhibiting the PI3K/AKT and ERK1/2 signaling axis. The GE11-siRNA-CSNPs design, which specifically targets TNBC cells, offers a novel approach for the treatment of breast cancer with improved effectiveness. This study suggests that GE11-siRNA-CSNPs could be a promising candidate for further assessment as an additional strategy in the treatment of TNBC.

Identifying inflammation-related targets of natural lactones using network pharmacology, molecular modeling and in vitro approaches

Natural lactones have been used in traditional and folklore medicine for centuries owing to their anti-inflammatory properties. The study uses a multifaceted approach to identify lead anti-inflammatory lactones from the SISTEMATX natural products database. The study analyzed the natural lactone database, revealing 18 lactones linked to inflammation targets. The primary targets were PTGES, PTGS1, COX-2, ALOX5 and IL1B. STX 12273 was the best hit, with the lowest binding energy and potential for inhibiting the COX-2 enzyme. The study suggested natural lactone, STX 12273, from the SISTEMATX database with anti-inflammatory potential and postulated its use for inflammation treatment or prevention.Communicated by Ramaswamy H. Sarma.

Structure-based prediction of protein-protein interaction network in rice

Comprehensive protein-protein interaction (PPI) maps are critical for understanding the functional organization of the proteome, but challenging to produce experimentally. Here, we developed a computational method for predicting PPIs based on protein docking. Evaluation of performance on benchmark sets demonstrated the ability of the docking-based method to accurately identify PPIs using predicted protein structures. By employing the docking-based method, we constructed a structurally resolved PPI network consisting of 24,653 interactions between 2,131 proteins, which greatly extends the current knowledge on the rice protein-protein interactome. Moreover, we mapped the trait-associated single nucleotide polymorphisms (SNPs) to the structural interactome, and computationally identified 14 SNPs that had significant consequences on PPI network. The protein structural interactome map provided a resource to facilitate functional investigation of PPI-perturbing alleles associated with agronomically important traits in rice.

A close-up shot of protein-protein docking, from experiment to theory and reverse with the PROTAC performers

PROTACs (Proteolysis Targeting Chimeras), heterobifunctional molecules, exhibit selectivity in degrading target proteins through E3 ubiquitin ligases. Designing effective PROTACs requires a deep understanding of the intricate binding interactions in the ternary complex (POI/PROTAC/E3 ligase), crucial for efficient target protein degradation. To address this challenge, we introduce a novel computational virtual screening method that considers essential amino acid interactions between the protein of interest and the chosen E3 ligase. This approach enhances accuracy and reliability, facilitating the strategic development of potent PROTACs. Utilizing a crystallized model of the VHL:PROTAC:SMARCA2BD ternary complex (PDB: 7Z6L), we assessed the effectiveness of our method. Our study reveals that increasing the number of essential restraints between the two proteins reduces the generated docking poses, leading to closer alignment with the experimental ternary complex. Specifically, utilizing three restraints showed the closest resemblance to the published complex, highlighting crucial interactions such as an H-bond between A:Gln 89 and B:Asn 67, along with two hydrophobic interactions: A:Gly 22 with B:Arg 69 and A:Glu 37 with B:Pro 99. This resulted in a significant decrease in the mean RMSD value from 31.8 and 31.0 Å to 24.4 Å, respectively. This underscores the importance of incorporating multiple essential restraints to enhance docking accuracy. Building on this progress, we introduce a systematic approach to design potential PROTACs between the Estrogen receptor and the E3 ligase, utilizing bridging intermediates with 4, 6, or 7 carbon atoms. By providing a more accurate and efficient means of identifying optimal PROTAC candidates, this approach has the potential to accelerate the development of targeted therapies and reduce the time and costs associated with drug discovery.Communicated by Ramaswamy H. Sarma.

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.

GRAMM Web Server for Protein Docking

Prediction of the structure of protein complexes by docking methods is a well-established research field. The intermolecular energy landscapes in protein-protein interactions can be used to refine docking predictions and to detect macro-characteristics, such as the binding funnel. A new GRAMM web server for protein docking predicts a spectrum of docking poses that characterize the intermolecular energy landscape in protein interaction. A user-friendly interface provides options to choose free or template-based docking, as well as other advanced features, such as clustering of the docking poses, and interactive visualization of the docked models.

Exploring the Molecular Mechanism of HongTeng Decoction against Inflammation based on Network Analysis and Experiments Validation

BACKGROUND

HongTeng Decoction (HTD) is a traditional Chinese medicine that is widely used to treat bacterial infections and chronic inflammation. However, its pharmacological mechanism is not clear. Here, network pharmacology and experimental verification were applied to investigate the drug targets and potential mechanisms of HTD in inflammation treatment.

METHODS

The active ingredients of HTD were collected from the multi-source databases and clarified by Q Exactive Orbitrap analysis in the treatment of inflammation. Then, molecular docking technology was used to explore the binding ability of key active ingredients and targets in HTD. In vitro experiments, the inflammatory factors and MAPK signaling pathways are detected to verify the anti-inflammatory effect of HTD on the RAW264.7 cells. Finally, the anti-inflammatory effect of HTD was evaluated in LPS induced mice model.

RESULTS

A total of 236 active compounds and 492 targets of HTD were obtained through database screening, and 954 potential targets of inflammation were identified. Finally, 164 possible targets of HTD acting on inflammation were obtained. The PPI analysis and KEGG enrichment analyses showed that the targets of HTD in inflammation were mostly related to the MAPK signaling pathway, the IL-17 signaling pathway, and the TNF signaling pathway. By integrating the results of the network analysis, the core targets of HTD in inflammation mainly include MAPK3, TNF, MMP9, IL6, EGFR, and NFKBIA. The molecular docking results indicated solid binding activity between MAPK3-naringenin and MAPK3-paeonol. It has been shown that HTD could inhibit the levels of inflammatory factors, IL6 and TNF-α, as well as the splenic index in the LPS-stimulated mice. Moreover, HTD could regulate protein expression levels of p-JNK1/2, and p-ERK1/2, which reflects the inhibitory effect of HTD on the MAPKS signaling pathway.

CONCLUSION

Our study is expected to provide the pharmacological mechanisms by which HTD may be a promising anti-inflammatory drug for future clinical trials.

In silico prediction of interaction between Nipah virus attachment glycoprotein and host cell receptors Ephrin-B2 and Ephrin-B3 in domestic and peridomestic mammals

Nipah virus (NiV) is a lethal bat-borne zoonotic virus that causes mild to acute respiratory distress and neurological manifestations in humans with a high mortality rate. NiV transmission to humans occurs via consumption of bat-contaminated fruit and date palm sap (DPS), or through direct contact with infected individuals and livestock. Since NiV outbreaks were first reported in pigs from Malaysia and Singapore, non-neutralizing antibodies against NiV attachment Glycoprotein (G) have also been detected in a few domestic mammals. NiV infection is initiated after NiV G binds to the host cell receptors Ephrin-B2 and Ephrin-B3. In this study, we assessed the degree of NiV host tropism in domestic and peridomestic mammals commonly found in Bangladesh that may be crucial in the transmission of NiV by serving as intermediate hosts. We carried out a protein-protein docking analysis of NiV G complexes (n = 52) with Ephrin-B2 and B3 of 13 domestic and peridomestic species using bioinformatics tools. Protein models were generated by homology modelling and the structures were validated for model quality. The different protein-protein complexes in this study were stable, and their binding affinity (ΔG) scores ranged between -8.0 to -19.1 kcal/mol. NiV Bangladesh (NiV-B) strain displayed stronger binding to Ephrin receptors, especially with Ephrin-B3 than the NiV Malaysia (NiV-M) strain, correlating with the observed higher pathogenicity of NiV-B strains. From the docking result, we found that Ephrin receptors of domestic rat (R. norvegicus) had a higher binding affinity for NiV G, suggesting greater susceptibility to NiV infections compared to other study species. Investigations for NiV exposure to domestic/peridomestic animals will help us knowing more the possible role of rats and other animals as intermediate hosts of NiV and would improve future NiV outbreak control and prevention in humans and domestic animals.

Multi-domain and complex protein structure prediction using inter-domain interactions from deep learning

Accurately capturing domain-domain interactions is key to understanding protein function and designing structure-based drugs. Although AlphaFold2 has made a breakthrough on single domain, it should be noted that the structure modeling for multi-domain protein and complex remains a challenge. In this study, we developed a multi-domain and complex structure assembly protocol, named DeepAssembly, based on domain segmentation and single domain modeling algorithms. Firstly, DeepAssembly uses a population-based evolutionary algorithm to assemble multi-domain proteins by inter-domain interactions inferred from a developed deep learning network. Secondly, protein complexes are assembled by means of domains rather than chains using DeepAssembly. Experimental results show that on 219 multi-domain proteins, the average inter-domain distance precision by DeepAssembly is 22.7% higher than that of AlphaFold2. Moreover, DeepAssembly improves accuracy by 13.1% for 164 multi-domain structures with low confidence deposited in AlphaFold database. We apply DeepAssembly for the prediction of 247 heterodimers. We find that DeepAssembly successfully predicts the interface (DockQ ≥ 0.23) for 32.4% of the dimers, suggesting a lighter way to assemble complex structures by treating domains as assembly units and using inter-domain interactions learned from monomer structures.

Revolutionizing Drug Targeting Strategies: Integrating Artificial Intelligence and Structure-Based Methods in PROTAC Development

PROteolysis TArgeting Chimera (PROTAC) is an emerging technology in chemical biology and drug discovery. This technique facilitates the complete removal of the target proteins that are "undruggable" or challenging to target through chemical molecules via the Ubiquitin-Proteasome System (UPS). PROTACs have been widely explored and outperformed not only in cancer but also in other diseases. During the past few decades, several academic institutes and pharma companies have poured more efforts into PROTAC-related technologies, setting the stage for several major degrader trial readouts in clinical phases. Despite their promising results, the formation of robust ternary orientation, off-target activity, poor permeability, and binding affinity are some of the limitations that hinder their development. Recent advancements in computational technologies have facilitated progress in the development of PROTACs. Researchers have been able to utilize these technologies to explore a wider range of E3 ligases and optimize linkers, thereby gaining a better understanding of the effectiveness and safety of PROTACs in clinical settings. In this review, we briefly explore the computational strategies reported to date for the formation of PROTAC components and discuss the key challenges and opportunities for further research in this area.

Integrated Network Pharmacology, Molecular Docking and Animal Experiment to Explore the Efficacy and Potential Mechanism of Baiyu Decoction Against Ulcerative Colitis by Enema

Background

Baiyu Decoction (BYD), a clinical prescription of traditional Chinese medicine, has been proven to be valuable for treating ulcerative colitis (UC) by enema. However, the mechanism of BYD against UC remains unclear.

Purpose

A combination of bioinformatics methods including network pharmacology and molecular docking and animal experiments were utilized to investigate the potential mechanism of BYD in the treatment of UC.

Materials and Methods

Firstly, the representative compounds of each herb in BYD were detected by liquid chromatography-mass spectrometry. Subsequently, we predicted the core targets and potential pathways of BYD for treating UC through network pharmacology. And rat colitis model was established with dextran sodium sulfate. UC rats were subjected to BYD enema administration, during which we recorded body weight changes, disease activity index, and colon length to assess the effectiveness of BYD. Besides, quantitative real-time PCR, western blotting, ELISA and immunofluorescence were used to detect intestinal inflammatory factors, intestinal barrier biomarkers and TOLL-like receptor pathway in rats. Finally, the core components and targets of BYD were subjected to molecular docking so as to further validate the results of network pharmacology.

Results

A total of 41 active compositions and 203 targets related to BYD-UC were subjected to screening. The results of bioinformatics analysis showed that quercetin and kaempferol may be the main compounds. Additionally, AKT1, IL-6, TP53, TNF and IL-1β were regarded as potential therapeutic targets. KEGG results explained that TOLL-like receptor pathway might play a pivotal role in BYD protecting against UC. In addition, animal experiments and molecular docking validated the network pharmacology results. BYD enema treatment can reduce body weight loss, lower disease activity index score, reverse colon shortening, relieve intestinal inflammation, protect intestinal barrier, and inhibit TOLL-like receptor pathway in UC rats. Besides, molecular docking suggested that quercetin and kaempferol docked well with TLR4, AKT1, IL-6, TP53.

Conclusion

Utilizing network pharmacology, animal studies, and molecular docking, enema therapy with BYD was confirmed to have anti-UC efficacy by alleviating intestinal inflammation, protecting the intestinal barrier, and inhibiting the TOLL-like receptor pathway. Researchers should focus not only on oral medications but also on the rectal administration of medications in furtherance of the cure of ulcerative colitis.

Higher-Order Structural Organization of the Mitochondrial Proteome Charted by In Situ Cross-Linking Mass Spectrometry

Mitochondria are densely packed with proteins, of which most are involved physically or more transiently in protein-protein interactions (PPIs). Mitochondria host among others all enzymes of the Krebs cycle and the oxidative phosphorylation pathway and are foremost associated with cellular bioenergetics. However, mitochondria are also important contributors to apoptotic cell death and contain their own genome indicating that they play additionally an eminent role in processes beyond bioenergetics. Despite intense efforts in identifying and characterizing mitochondrial protein complexes by structural biology and proteomics techniques, many PPIs have remained elusive. Several of these (membrane embedded) PPIs are less stable in vitro hampering their characterization by most contemporary methods in structural biology. Particularly in these cases, cross-linking mass spectrometry (XL-MS) has proven valuable for the in-depth characterization of mitochondrial protein complexes in situ. Here, we highlight experimental strategies for the analysis of proteome-wide PPIs in mitochondria using XL-MS. We showcase the ability of in situ XL-MS as a tool to map suborganelle interactions and topologies and aid in refining structural models of protein complexes. We describe some of the most recent technological advances in XL-MS that may benefit the in situ characterization of PPIs even further, especially when combined with electron microscopy and structural modeling.

Pyrrole-Based Enaminones as Building Blocks for the Synthesis of Indolizines and Pyrrolo[1,2-a]pyrazines Showing Potent Antifungal Activity

As a new approach, pyrrolo[1,2-a]pyrazines were synthesized through the cyclization of 2-formylpyrrole-based enaminones in the presence of ammonium acetate. The enaminones were prepared with a straightforward method, reacting the corresponding alkyl 2-(2-formyl-1H-pyrrol-1-yl)acetates, 2-(2-formyl-1H-pyrrol-1-yl)acetonitrile, and 2-(2-formyl-1H-pyrrol-1-yl)acetophenones with DMFDMA. Analogous enaminones elaborated from alkyl (E)-3-(1H-pyrrol-2-yl)acrylates were treated with a Lewis acid to afford indolizines. The antifungal activity of the series of substituted pyrroles, pyrrole-based enaminones, pyrrolo[1,2-a]pyrazines, and indolizines was evaluated on six Candida spp., including two multidrug-resistant ones. Compared to the reference drugs, most test compounds produced a more robust antifungal effect. Docking analysis suggests that the inhibition of yeast growth was probably mediated by the interaction of the compounds with the catalytic site of HMGR of the Candida species.

ZEPPI: proteome-scale sequence-based evaluation of protein-protein interaction models

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence co-evolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. ZEPPI can be implemented on a proteome-wide scale as evidenced by calculations on millions of structural models of dimeric complexes in the E. coli and human interactomes found in the PrePPI database. A number of examples that illustrate how these tools can yield novel functional hypotheses are provided.

Prediction of protein interactions is essential for studying biomolecular mechanisms

Structural characterization of protein interactions is essential for our ability to understand and modulate physiological processes. Computational approaches to modeling of protein complexes provide structural information that far exceeds capabilities of the existing experimental techniques. Protein structure prediction in general, and prediction of protein interactions in particular, has been revolutionized by the rapid progress in Deep Learning techniques. The work of Schweke et al. (Proteomics 2023, 23, 2200323) presents a community-wide study of an important problem of distinguishing physiological protein-protein complexes/interfaces (experimentally determined or modeled) from non-physiological ones. The authors designed and generated a large benchmark set of physiological and non-physiological homodimeric complexes, and evaluated a large set of scoring functions, as well as AlphaFold predictions, on their ability to discriminate the non-physiological interfaces. The problem of separating physiological interfaces from non-physiological ones is very difficult, largely due to the lack of a clear distinction between the two categories in a crowded environment inside a living cell. Still, the ability to identify key physiologically significant interfaces in the variety of possible configurations of a protein-protein complex is important. The study presents a major data resource and methodological development in this important direction for molecular and cellular biology.

PointDE: Protein Docking Evaluation Using 3D Point Cloud Neural Network

Protein-protein interactions (PPIs) play essential roles in many vital movements and the determination of protein complex structure is helpful to discover the mechanism of PPI. Protein-protein docking is being developed to model the structure of the protein. However, there is still a challenge to selecting the near-native decoys generated by protein-protein docking. Here, we propose a docking evaluation method using 3D point cloud neural network named PointDE. PointDE transforms protein structure to the point cloud. Using the state-of-the-art point cloud network architecture and a novel grouping mechanism, PointDE can capture the geometries of the point cloud and learn the interaction information from the protein interface. On public datasets, PointDE surpasses the state-of-the-art method using deep learning. To further explore the ability of our method in different types of protein structures, we developed a new dataset generated by high-quality antibody-antigen complexes. The result in this antibody-antigen dataset shows the strong performance of PointDE, which will be helpful for the understanding of PPI mechanisms.

A New Tool to Study the Binding Behavior of Intrinsically Disordered Proteins

Understanding the binding behavior and conformational dynamics of intrinsically disordered proteins (IDPs) is crucial for unraveling their regulatory roles in biological processes. However, their lack of stable 3D structures poses challenges for analysis. To address this, we propose an algorithm that explores IDP binding behavior with protein complexes by extracting topological and geometric features from the protein surface model. Our algorithm identifies a geometrically favorable binding pose for the IDP and plans a feasible trajectory to evaluate its transition to the docking position. We focus on IDPs from Homo sapiens and Mus-musculus, investigating their interaction with the Plasmodium falciparum (PF) pathogen associated with malaria-related deaths. We compare our algorithm with HawkDock and HDOCK docking tools for quantitative (computation time) and qualitative (binding affinity) measures. Our results indicated that our method outperformed the compared methods in computation performance and binding affinity in experimental conformations.

Spatial Layouts of Low-Entropy Hydration Shells Guide Protein Binding

Protein-protein binding enables orderly biological self-organization and is therefore considered a miracle of nature. Protein‒protein binding is driven by electrostatic forces, hydrogen bonding, van der Waals force, and hydrophobic interactions. Among these physical forces, only hydrophobic interactions can be considered long-range intermolecular attractions between proteins due to the electrostatic shielding of surrounding water molecules. Low-entropy hydration shells around proteins drive hydrophobic attraction among them that essentially coordinate protein‒protein binding. Here, an innovative method is developed for identifying low-entropy regions of hydration shells of proteins by screening off pseudohydrophilic groups on protein surfaces and revealing that large low-entropy regions of the hydration shells typically cover the binding sites of individual proteins. According to an analysis of determined protein complex structures, shape matching between a large low-entropy hydration shell region of a protein and that of its partner at the binding sites is revealed as a universal law. Protein‒protein binding is thus found to be mainly guided by hydrophobic collapse between the shape-matched low-entropy hydration shells that is verified by bioinformatics analyses of hundreds of structures of protein complexes, which cover four test systems. A simple algorithm is proposed to accurately predict protein binding sites.

Narrow funnel-like interaction energy distribution is an indicator of specific protein interaction partner

Protein interaction networks underlie countless biological mechanisms. However, most protein interaction predictions are based on biological evidence that are biased to well-known protein interaction or physical evidence that exhibits low accuracy for weak interactions and requires high computational power. In this study, a novel method has been suggested to predict protein interaction partners by investigating narrow funnel-like interaction energy distribution. In this study, it was demonstrated that various protein interactions including kinases and E3 ubiquitin ligases have narrow funnel-like interaction energy distribution. To analyze protein interaction distribution, modified scores of iRMS and TM-score are introduced. Then, using these scores, algorithm and deep learning model for prediction of protein interaction partner and substrate of kinase and E3 ubiquitin ligase were developed. The prediction accuracy was similar to or even better than that of yeast two-hybrid screening. Ultimately, this knowledge-free protein interaction prediction method will broaden our understanding of protein interaction networks.

GDockScore: a graph-based protein-protein docking scoring function

Summary

Protein complexes play vital roles in a variety of biological processes, such as mediating biochemical reactions, the immune response and cell signalling, with 3D structure specifying function. Computational docking methods provide a means to determine the interface between two complexed polypeptide chains without using time-consuming experimental techniques. The docking process requires the optimal solution to be selected with a scoring function. Here, we propose a novel graph-based deep learning model that utilizes mathematical graph representations of proteins to learn a scoring function (GDockScore). GDockScore was pre-trained on docking outputs generated with the Protein Data Bank biounits and the RosettaDock protocol, and then fine-tuned on HADDOCK decoys generated on the ZDOCK Protein Docking Benchmark. GDockScore performs similarly to the Rosetta scoring function on docking decoys generated using the RosettaDock protocol. Furthermore, state-of-the-art is achieved on the CAPRI score set, a challenging dataset for developing docking scoring functions.

Availability and implementation

The model implementation is available at https://gitlab.com/mcfeemat/gdockscore.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

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.

Protein-mediated interactions in the dynamic regulation of acute inflammation

Protein-mediated interactions are the fundamental mechanism through which cells regulate health and disease. These interactions require physical contact between proteins and their respective targets of interest. These targets include not only other proteins but also nucleic acids and other important molecules as well. These proteins are often involved in multibody complexes that work dynamically to regulate cellular health and function. Various techniques have been adapted to study these important interactions, such as affinity-based assays, mass spectrometry, and fluorescent detection. The application of these techniques has led to a greater understanding of how protein interactions are responsible for both the instigation and resolution of acute inflammatory diseases. These pursuits aim to provide opportunities to target specific protein interactions to alleviate acute inflammation.

Directed targeting of B-cell maturation antigen-specific CAR T cells by bioinformatic approaches: From in-silico to in-vitro

AIMS

Chimeric Antigen Receptor (CAR) T-cell is a breakthrough in cancer immunotherapy. The primary step of successful CAR T cell therapy is designing a specific single-chain fragment variable (scFv). This study aims to verify the designed anti-BCMA (B cell maturation antigen) CAR using bioinformatic techniques with the following experimental evaluations.

MAIN METHODS

Following the second generation of anti-BCMA CAR designing, the protein structure, function prediction, physicochemical complementarity at the ligand-receptor interface, and biding sites analysis of anti-BCMA CAR construct were confirmed using different modeling and docking server, including Expasy, I-TASSER, HDock, and PyMOL software. To generate CAR T-cells, isolated T cells were transduced. Then, anti-BCMA CAR mRNA and its surface expression were confirmed by real-time -PCR and flow cytometry methods, respectively. To evaluate the surface expression of anti-BCMA CAR, anti-(Fab')2 and anti-CD8 antibodies were employed. Finally, anti-BCMA CAR T cells were co-cultured with BCMA+/- cell lines to assess the expression of CD69 and CD107a as activation and cytotoxicity markers.

KEY FINDINGS

In-silico results approved the suitable protein folding, perfect orientation, and correct locating of functional domains at the receptor-ligand binding site. The in-vitro results confirmed high expression of scFv (89 ± 1.15% (and CD8α (54 ± 2.88%). The expression of CD69 (91.97 ± 1.7%) and CD107a (92.05 ± 1.29%) were significantly increased, indicating appropriate activation and cytotoxicity.

SIGNIFICANCE

In-silico studies before experimental assessments are crucial for state-of-art CAR designing. Highly activation and cytotoxicity of anti-BCMA CAR T-cell revealed that our CAR construct methodology would be applicable to define the road map of CAR T cell therapy.

Reconstitution and expression of mcy gene cluster in the model cyanobacterium Synechococcus 7942 reveals a role of MC-LR in cell division

Cyanobacterial blooms pose a serious threat to public health due to the presence of cyanotoxins. Microcystin-LR (MC-LR) produced by Microcystis aeruginosa is the most common cyanotoxins. Due to the limitation of isolation, purification, and genetic manipulation techniques, it is difficult to study and verify in situ the biosynthetic pathways and molecular mechanisms of MC-LR. We reassembled the biosynthetic gene cluster (mcy cluster) of MC-LR in vitro by synthetic biology, designed and constructed the strong bidirectional promoter biPpsbA₂ , transformed it into Synechococcus 7942, and successfully expressed MC-LR at a level of 0.006-0.018 fg cell^-1  d^-1 . We found the expression of MC-LR led to abnormal cell division and cellular filamentation, further using various methods proved that by irreversibly competing its GTP-binding site, MC-LR inhibits assembly of the cell division protein FtsZ. The study represents the first reconstitution and expression of the mcy cluster and the autotrophic production of MC-LR in model cyanobacterium, which lays the foundation for resolving the microcystins biosynthesis pathway. The discovered role of MC-LR in cell division reveals a mechanism of how blooming cyanobacteria gain a competitive edge over their nonblooming counterparts.

Recent advances in predicting and modeling protein-protein interactions

Protein-protein interactions (PPIs) drive biological processes, and disruption of PPIs can cause disease. With recent breakthroughs in structure prediction and a deluge of genomic sequence data, computational methods to predict PPIs and model spatial structures of protein complexes are now approaching the accuracy of experimental approaches for permanent interactions and show promise for elucidating transient interactions. As we describe here, the key to this success is rich evolutionary information deciphered from thousands of homologous sequences that coevolve in interacting partners. This covariation signal, revealed by sophisticated statistical and machine learning (ML) algorithms, predicts physiological interactions. Accurate artificial intelligence (AI)-based modeling of protein structures promises to provide accurate 3D models of PPIs at a proteome-wide scale.

Improved inter-protein contact prediction using dimensional hybrid residual networks and protein language models

The knowledge of contacting residue pairs between interacting proteins is very useful for the structural characterization of protein-protein interactions (PPIs). However, accurately identifying the tens of contacting ones from hundreds of thousands of inter-protein residue pairs is extremely challenging, and performances of the state-of-the-art inter-protein contact prediction methods are still quite limited. In this study, we developed a deep learning method for inter-protein contact prediction, which is referred to as DRN-1D2D_Inter. Specifically, we employed pretrained protein language models to generate structural information-enriched input features to residual networks formed by dimensional hybrid residual blocks to perform inter-protein contact prediction. Extensively bechmarking DRN-1D2D_Inter on multiple datasets, including both heteromeric PPIs and homomeric PPIs, we show DRN-1D2D_Inter consistently and significantly outperformed two state-of-the-art inter-protein contact prediction methods, including GLINTER and DeepHomo, although both the latter two methods leveraged the native structures of interacting proteins in the prediction, and DRN-1D2D_Inter made the prediction purely from sequences. We further show that applying the predicted contacts as constraints for protein-protein docking can significantly improve its performance for protein complex structure prediction.

Analysis of protein-protein interface with incorporating low-frequency molecular interactions in molecular dynamics simulation

Protein-protein interactions are vital for various biological processes such as immune reaction, signal transduction, and viral infection. Molecular Dynamics (MD) simulation is a powerful tool for analyzing non-covalent interactions between two protein molecules. In general, MD simulation studies on the protein-protein interface have focused on the analysis of major and frequent molecular interactions. In this study, we demonstrate that minor interactions with low-frequency need to be incorporated to analyze the molecular interactions in the protein-protein interface more efficiently using the complex of SARS-CoV2-RBD and ACE2 receptor as a model system. It was observed that the dominance of interactions in the MD-simulated structures didn't directly correlate with the interactions in the experimentally determined structure. The interactions from the experimentally determined structure could be reproduced better in the ensemble of MD simulated structures by including the less frequent interactions compared to the norm of choosing only highly frequent interactions. Residue Interaction Networks (RINs) analysis also showed that the critical residues in the protein-protein interface could be more efficiently identified by incorporating low-frequency interactions in MD simulation. It is expected that the approach proposed in this study can be a new way of studying protein-protein interaction through MD simulation.

Structure based Drug Designing Approaches in SARS-CoV-2 Spike Inhibitor Design

The COVID-19 outbreak and the pandemic situation have hastened the research community to design a novel drug and vaccine against its causative organism, the SARS-CoV-2. The spike glycoprotein present on the surface of this pathogenic organism plays an immense role in viral entry and antigenicity. Hence, it is considered an important drug target in COVID-19 drug design. Several three-dimensional crystal structures of this SARS-CoV-2 spike protein have been identified and deposited in the Protein DataBank during the pandemic period. This accelerated the research in computer- aided drug designing, especially in the field of structure-based drug designing. This review summarizes various structure-based drug design approaches applied to this SARS-CoV-2 spike protein and its findings. Specifically, it is focused on different structure-based approaches such as molecular docking, high-throughput virtual screening, molecular dynamics simulation, drug repurposing, and target-based pharmacophore modelling and screening. These structural approaches have been applied to different ligands and datasets such as FDA-approved drugs, small molecular chemical compounds, chemical libraries, chemical databases, structural analogs, and natural compounds, which resulted in the prediction of spike inhibitors, spike-ACE-2 interface inhibitors, and allosteric inhibitors.

Artificial Intelligence and Machine Learning Technology Driven Modern Drug Discovery and Development

The discovery and advances of medicines may be considered as the ultimate relevant translational science effort that adds to human invulnerability and happiness. But advancing a fresh medication is a quite convoluted, costly, and protracted operation, normally costing USD ~2.6 billion and consuming a mean time span of 12 years. Methods to cut back expenditure and hasten new drug discovery have prompted an arduous and compelling brainstorming exercise in the pharmaceutical industry. The engagement of Artificial Intelligence (AI), including the deep-learning (DL) component in particular, has been facilitated by the employment of classified big data, in concert with strikingly reinforced computing prowess and cloud storage, across all fields. AI has energized computer-facilitated drug discovery. An unrestricted espousing of machine learning (ML), especially DL, in many scientific specialties, and the technological refinements in computing hardware and software, in concert with various aspects of the problem, sustain this progress. ML algorithms have been extensively engaged for computer-facilitated drug discovery. DL methods, such as artificial neural networks (ANNs) comprising multiple buried processing layers, have of late seen a resurgence due to their capability to power automatic attribute elicitations from the input data, coupled with their ability to obtain nonlinear input-output pertinencies. Such features of DL methods augment classical ML techniques which bank on human-contrived molecular descriptors. A major part of the early reluctance concerning utility of AI in pharmaceutical discovery has begun to melt, thereby advancing medicinal chemistry. AI, along with modern experimental technical knowledge, is anticipated to invigorate the quest for new and improved pharmaceuticals in an expeditious, economical, and increasingly compelling manner. DL-facilitated methods have just initiated kickstarting for some integral issues in drug discovery. Many technological advances, such as "message-passing paradigms", "spatial-symmetry-preserving networks", "hybrid de novo design", and other ingenious ML exemplars, will definitely come to be pervasively widespread and help dissect many of the biggest, and most intriguing inquiries. Open data allocation and model augmentation will exert a decisive hold during the progress of drug discovery employing AI. This review will address the impending utilizations of AI to refine and bolster the drug discovery operation.

Transformer-based deep learning for predicting protein properties in the life sciences

Recent developments in deep learning, coupled with an increasing number of sequenced proteins, have led to a breakthrough in life science applications, in particular in protein property prediction. There is hope that deep learning can close the gap between the number of sequenced proteins and proteins with known properties based on lab experiments. Language models from the field of natural language processing have gained popularity for protein property predictions and have led to a new computational revolution in biology, where old prediction results are being improved regularly. Such models can learn useful multipurpose representations of proteins from large open repositories of protein sequences and can be used, for instance, to predict protein properties. The field of natural language processing is growing quickly because of developments in a class of models based on a particular model-the Transformer model. We review recent developments and the use of large-scale Transformer models in applications for predicting protein characteristics and how such models can be used to predict, for example, post-translational modifications. We review shortcomings of other deep learning models and explain how the Transformer models have quickly proven to be a very promising way to unravel information hidden in the sequences of amino acids.

MetaScore: A Novel Machine-Learning-Based Approach to Improve Traditional Scoring Functions for Scoring Protein-Protein Docking Conformations

Protein-protein interactions play a ubiquitous role in biological function. Knowledge of the three-dimensional (3D) structures of the complexes they form is essential for understanding the structural basis of those interactions and how they orchestrate key cellular processes. Computational docking has become an indispensable alternative to the expensive and time-consuming experimental approaches for determining the 3D structures of protein complexes. Despite recent progress, identifying near-native models from a large set of conformations sampled by docking-the so-called scoring problem-still has considerable room for improvement. We present MetaScore, a new machine-learning-based approach to improve the scoring of docked conformations. MetaScore utilizes a random forest (RF) classifier trained to distinguish near-native from non-native conformations using their protein-protein interfacial features. The features include physicochemical properties, energy terms, interaction-propensity-based features, geometric properties, interface topology features, evolutionary conservation, and also scores produced by traditional scoring functions (SFs). MetaScore scores docked conformations by simply averaging the score produced by the RF classifier with that produced by any traditional SF. We demonstrate that (i) MetaScore consistently outperforms each of the nine traditional SFs included in this work in terms of success rate and hit rate evaluated over conformations ranked among the top 10; (ii) an ensemble method, MetaScore-Ensemble, that combines 10 variants of MetaScore obtained by combining the RF score with each of the traditional SFs outperforms each of the MetaScore variants. We conclude that the performance of traditional SFs can be improved upon by using machine learning to judiciously leverage protein-protein interfacial features and by using ensemble methods to combine multiple scoring functions.

DeepPRObind: Modular Deep Learner that Accurately Predicts Structure and Disorder-Annotated Protein Binding Residues

Current sequence-based predictors of protein-binding residues (PBRs) belong to two distinct categories: structure-trained vs. intrinsic disorder-trained. Since disordered PBRs differ from structured PBRs in several ways, including ability to bind multiple partners by folding into different conformations and enrichment in different amino acids, the structure-trained and disorder-trained predictors were shown to provide inaccurate results for the other annotation type. A simple consensus-based solution that combines structure- and disorder-trained methods provides limited levels of predictive performance and generates relatively many cross-predictions, where residues that interact with other ligand types are predicted as PBRs. We address this unsolved problem by designing a novel and fast deep-learner, DeepPRObind, that relies on carefully designed modular convolutional architecture and uses innovative aggregate input features. Comparative empirical tests on a low-similarity test dataset reveal that DeepPRObind generates accurate predictions of structured and disordered PBRs and low amounts of cross-predictions, outperforming a comprehensive collection of 12 predictors of PBRs. Given the relatively low runtime of DeepPRObind (40 seconds per protein), we further validate its results based on an analysis of putative PBRs in the yeast proteome, confirming that interactions in disordered regions are enriched among hub proteins. We release DeepPRObind as a convenient web server at https://www.csuligroup.com/DeepPRObind/.

Pathogen-driven cancers from a structural perspective: Targeting host-pathogen protein-protein interactions

Host-pathogen interactions (HPIs) affect and involve multiple mechanisms in both the pathogen and the host. Pathogen interactions disrupt homeostasis in host cells, with their toxins interfering with host mechanisms, resulting in infections, diseases, and disorders, extending from AIDS and COVID-19, to cancer. Studies of the three-dimensional (3D) structures of host-pathogen complexes aim to understand how pathogens interact with their hosts. They also aim to contribute to the development of rational therapeutics, as well as preventive measures. However, structural studies are fraught with challenges toward these aims. This review describes the state-of-the-art in protein-protein interactions (PPIs) between the host and pathogens from the structural standpoint. It discusses computational aspects of predicting these PPIs, including machine learning (ML) and artificial intelligence (AI)-driven, and overviews available computational methods and their challenges. It concludes with examples of how theoretical computational approaches can result in a therapeutic agent with a potential of being used in the clinics, as well as future directions.

Mapping the energy landscape of PROTAC-mediated protein-protein interactions

A principal challenge in computational modeling of macromolecules is the vast conformational space that arises out of large numbers of atomic degrees of freedom. Recently, growing interest in building predictive models of complexes mediated by Proteolysis Targeting Chimeras (PROTACs) has led to the application of state-of-the-art computational techniques to tackle this problem. However, repurposing existing tools to carry out protein-protein docking and linker conformer generation independently results in extensive sampling of structures incompatible with PROTAC-mediated complex formation. Here we show that it is possible to restrict the search to the space of protein-protein conformations that can be bridged by a PROTAC molecule with a given linker composition by using a cyclic coordinate descent algorithm to position PROTACs into complex-bound configurations. We use this methodology to construct potential energy and solvation energy landscapes of PROTAC-mediated interactions. Our results suggest that desolvation of amino acids at interfaces could play a dominant role in PROTAC-mediated complex formation.

Using PlaPPISite to Predict and Analyze Plant Protein-Protein Interaction Sites

Proteome-wide characterization of protein-protein interactions (PPIs) is crucial to understand the functional roles of protein machinery within cells systematically. With the accumulation of PPI data in different plants, the interaction details of binary PPIs, such as the three-dimensional (3D) structural contexts of interaction sites/interfaces, are urgently demanded. To meet this requirement, we have developed a comprehensive and easy-to-use database called PlaPPISite ( http://zzdlab.com/plappisite/index.php ) to present interaction details for 13 plant interactomes. Here, we provide a clear guide on how to search and view protein interaction details through the PlaPPISite database. Firstly, the running environment of our database is introduced. Secondly, the input file format is briefly introduced. Moreover, we discussed which information related to interaction sites can be achieved through several examples. In addition, some notes about PlaPPISite are also provided. More importantly, we would like to emphasize the importance of interaction site information in plant systems biology through this user guide of PlaPPISite. In particular, the easily accessible 3D structures of PPIs in the coming post-AlphaFold2 era will definitely boost the application of plant interactome to decipher the molecular mechanisms of many fundamental biological issues.

Modulation of interaction of BRCA1-RAD51 and BRCA1-AURKA protein complexes by natural metabolites using as possible therapeutic intervention toward cardiotoxic effects of cancer drugs: an in-silico approach

Breast cancer type 1 susceptibility protein (BRCA1) plays an important role in maintaining genome stability and is known to interact with several proteins involved in cellular pathways, gene transcription regulation and DNA damage response. More than 40% of inherited breast cancer cases are due to BRCA1 mutation. It is also a prognostic marker in non-small cell lung cancer patients as well as a gatekeeper of cardiac function. Interaction of mutant BRCA1 with other proteins is known to disrupt the tumor suppression mechanism. Two directly interacting proteins with BRCA1 namely, DNA repair protein RAD51 (RAD51) and Aurora kinase A (AURKA), known to regulate homologous recombination (HR) and G/M cell cycle transition, respectively, form protein complex with both wild and mutant BRCA1. To analyze the interactions, protein-protein complexes were generated for each pair of proteins. In order to combat the cardiotoxic effects of cancer drugs, pharmacokinetically screened natural metabolites derived from plant, marine and bacterial sources and along with FDA-approved cancer drugs as control, were subjected to molecular docking. Piperoleine B and dihydrocircumin were the best docked natural metabolites in both RAD51 and AURKA complexes, respectively. Molecular dynamics simulation (MDS) analysis and binding free energy calculations for the best docked natural metabolite and drug for both the mutant BRCA1 complexes suggested better stability for the natural metabolites piperolein B and dihydrocurcumin as compared to drug. Thus, both natural metabolites could be further analyzed for their role against the cardiotoxic effects of cancer drugs through wet lab experiments.Communicated by Ramaswamy H. Sarma.

Exploring the potential mechanism of Fritiliariae Irrhosae Bulbus on ischemic stroke based on network pharmacology and experimental validation

Objective: To study the potential targets and molecular mechanisms of Fritiliariae Irrhosae Bulbus (FIB) in the treatment of ischemic strokes based on a network pharmacology strategy, with a combination of molecular docking and animal experiments.

Methods: The active components and targets of FIB were screened by TCMSP database and TCMIP database, and the related targets of ischemic strokes were screened by GeneCards, OMIM, CTD, and DrugBank, then the intersection targets of the two were taken. The protein interaction network was constructed by STRING, the PPI network diagram was drawn by using Cytoscape software, and the key targets of FIB treatment of ischemic strokes were analyzed by MCODE. The DAVID database was used for GO and KEGG enrichment analysis, and the potential pathway of FIB against ischemic strokes was obtained. Molecular docking was performed by using AutoDock Tools 1.5.6 software. Finally, a mouse model of ischemic stroke was established, and the results of network pharmacology were verified by in vivo experiments. Realtime Polymerase Chain Reaction was used to detect the expression levels of relevant mRNAs in the mouse brain tissue. Western blot was used to detect the expression levels of related proteins in the mouse brain tissue.

Results: 13 kinds of active components of FIB were screened, 31 targets were found in the intersection of FIB and ischemic strokes, 10 key targets were obtained by MCODE analysis, 236 biological processes were involved in GO enrichment analysis, and key targets of KEGG enrichment analysis were mainly concentrated in Neuroactive light receptor interaction, Calcium signaling pathway, Cholinergic synapse, Hepatitis B, Apoptosis—multiple specifications, Pathways in cancer and other significantly related pathways. There was good binding activity between the screened main active components and target proteins when molecular docking was performed. Animal experiments showed that the infarct volume of brain tissue in the FIB treatment group was considerably reduced. RT-qPCR and the results of Western Blot showed that FIB could inhibit the expression of active-Caspase3, HSP90AA1, phosphorylated C-JUN, and COX2.

Conclusion: Based on network pharmacology, the effect of FIB in the treatment of ischemic strokes was discussed through the multi-component-multi-target-multi-pathway. The therapeutic effect and potential mechanisms of FIB on ischemic strokes were preliminarily explored, which provided a ground work for further researches on the pharmacodynamic material basis, mechanism of action and clinical application.

Vetinformatics from functional genomics to drug discovery: Insights into decoding complex molecular mechanisms of livestock systems in veterinary science

Having played important roles in human growth and development, livestock animals are regarded as integral parts of society. However, industrialization has depleted natural resources and exacerbated climate change worldwide, spurring the emergence of various diseases that reduce livestock productivity. Meanwhile, a growing human population demands sufficient food to meet their needs, necessitating innovations in veterinary sciences that increase productivity both quantitatively and qualitatively. We have been able to address various challenges facing veterinary and farm systems with new scientific and technological advances, which might open new opportunities for research. Recent breakthroughs in multi-omics platforms have produced a wealth of genetic and genomic data for livestock that must be converted into knowledge for breeding, disease prevention and management, productivity, and sustainability. Vetinformatics is regarded as a new bioinformatics research concept or approach that is revolutionizing the field of veterinary science. It employs an interdisciplinary approach to understand the complex molecular mechanisms of animal systems in order to expedite veterinary research, ensuring food and nutritional security. This review article highlights the background, recent advances, challenges, opportunities, and application of vetinformatics for quality veterinary services.

Immunoinformatics Approach to Design Novel Subunit Vaccine against the Epstein-Barr Virus

Epstein-Barr virus (EBV) is a lymphotropic virus responsible for numerous epithelial and lymphoid cell malignancies, including gastric carcinoma, Hodgkin's lymphoma, nasopharyngeal carcinoma, and Burkitt lymphoma. Hundreds of thousands of people worldwide get infected with this virus, and in most cases, this viral infection leads to cancer. Although researchers are trying to develop potential vaccines and drug therapeutics, there is still no effective vaccine to combat this virus. In this study, the immunoinformatics approach was utilized to develop a potential multiepitope subunit vaccine against the two most common subtypes of EBV, targeting three of their virulent envelope glycoproteins. Eleven cytotoxic T lymphocyte (CTL) epitopes, 11 helper T lymphocyte (HTL) epitopes, and 10 B-cell lymphocyte (BCL) epitopes were predicted to be antigenic, nonallergenic, nontoxic, and fully conserved among the two subtypes, and nonhuman homologs were used for constructing the vaccine after much analysis. Later, further validation experiments, including molecular docking with different immune receptors (e.g., Toll-like receptors [TLRs]), molecular dynamics simulation analyses (including root means square deviation [RMSD], root mean square fluctuation [RMSF], radius of gyration [Rg], principal-component analysis [PCA], dynamic cross-correlation [DCC], definition of the secondary structure of proteins [DSSP], and Molecular Mechanics Poisson-Boltzmann Surface Area [MM-PBSA]), and immune simulation analyses generated promising results, ensuring the safe and stable response of the vaccine with specific immune receptors after potential administration within the human body. The vaccine's high binding affinity with TLRs was revealed in the docking study, and a very stable interaction throughout the simulation proved the potential high efficacy of the proposed vaccine. Further, in silico cloning was also conducted to design an efficient mass production strategy for future bulk industrial vaccine production. IMPORTANCE Epstein-Barr virus (EBV) vaccines have been developing for over 30 years, but polyphyletic and therapeutic vaccines have failed to get licensed. Our vaccine surpasses the limitations of many such vaccines and remains very promising, which is crucial because the infection rate is higher than most viral infections, affecting a whopping 90% of the adult population. One of the major identifications covers a holistic analysis of populations worldwide, giving us crucial information about its effectiveness for everyone's unique immunological system. We targeted three glycoproteins that enhance the virulence of the virus to design an epitope-based polyvalent vaccine against two different strains of EBV, type 1 and 2. Our methodology in this study is nonconventional yet swift to show effective results while designing vaccines.

Docking-based long timescale simulation of cell-size protein systems at atomic resolution

Advances in computational modeling have led to an increasing focus on larger biomolecular systems, up to the level of a cell. Protein interactions are a central component of cellular processes. Techniques for modeling protein interactions have been divided between two fields: protein docking (predicting the static structures of protein complexes) and molecular simulation (modeling the dynamics of protein association, for relatively short simulation times at atomic resolution). Our study combined the two approaches to reach very long simulation times. The study makes the model more adequate to the real cells, to explore cellular processes at atomic resolution to better understand molecular mechanisms of life, and to use this knowledge to improve our ability to treat diseases.

Computational methodologies are increasingly addressing modeling of the whole cell at the molecular level. Proteins and their interactions are the key component of cellular processes. Techniques for modeling protein interactions, thus far, have included protein docking and molecular simulation. The latter approaches account for the dynamics of the interactions but are relatively slow, if carried out at all-atom resolution, or are significantly coarse grained. Protein docking algorithms are far more efficient in sampling spatial coordinates. However, they do not account for the kinetics of the association (i.e., they do not involve the time coordinate). Our proof-of-concept study bridges the two modeling approaches, developing an approach that can reach unprecedented simulation timescales at all-atom resolution. The global intermolecular energy landscape of a large system of proteins was mapped by the pairwise fast Fourier transform docking and sampled in space and time by Monte Carlo simulations. The simulation protocol was parametrized on existing data and validated on a number of observations from experiments and molecular dynamics simulations. The simulation protocol performed consistently across very different systems of proteins at different protein concentrations. It recapitulated data on the previously observed protein diffusion rates and aggregation. The speed of calculation allows reaching second-long trajectories of protein systems that approach the size of the cells, at atomic resolution.

Mechanisms of Xiaochaihu Decoction on Treating Hepatic Fibrosis Explored by Network Pharmacology

Purpose. To explore the material basis and pharmacological mechanism of Xiaochaihu Decoction (XCHD), the classic Traditional Chinese Medicine (TCM) formula in inhibiting hepatic fibrosis (HF). Methods. The main components in XCHD were screened from the TCMSP database, ETCM database, and literature, and their potential targets were detected and predicted using the Swiss Target Prediction platform. The HF-related targets were retrieved and screened through GeneCard database and OMIM database, combined with GEO gene chips. The XCHD targets and HF targets were mapped to search common targets. The protein-protein interaction (PPI) network was acquired via the STRING11.0 database and analyzed visually using Cytoscape 3.8.0 software. The potential mechanisms of the common targets identified through GO and KEGG pathway enrichment analysis were analyzed by using Metascape database. The results were visualized through OmicShare Tools. The “XCHD compound-HF target” network was visually constructed by Cytoscape 3.8.0 software. AutoDockVina1.1.2 and PyMoL software were used to verify the molecular docking of XCHD main active compounds and HF key targets. Results. A total of 164 potential active compounds from XCHD were screened to act on 95 HF-related targets. Bioinformatics analysis revealed that quercetin, β-sitosterol, and kaempferol may be candidate agents, which acted on multiple targets like PTGS2, HSP90AA1, and PTGS1 and regulate multiple key biological pathways like IL-17 signaling pathway, TNF signaling pathway and PI3K-Akt signaling pathway to relieve HF. Moreover, molecular docking suggested that quercetin and PTGS2 could statically bind and interact with each other through amino acid residues val-349, LEU-352, PHE-381, etc. Conclusion. This work provides a systems perspective to study the relationship between Chinese medicines and diseases. The therapeutic efficacy of XCHD on HF was the sum of multitarget and multi-approach effects from the bioactive ingredients. This study could be one of the cornerstones for further research.

An in silico study of how histone tail conformation affects the binding affinity of ING family proteins

Background

Due to its intrinsically disordered nature, the histone tail is conformationally heterogenic. Therefore, it provides specific binding sites for different binding proteins or factors through reversible post-translational modifications (PTMs). For instance, experimental studies stated that the ING family binds with the histone tail that has methylation on the lysine in position 4. However, numerous complexes featuring a methylated fourth lysine residue of the histone tail can be found in the UniProt database. So the question arose if other factors like the conformation of the histone tail affect the binding affinity.

Methods

The crystal structure of the PHD finger domain from the proteins ING1, ING2, ING4, and ING5 are docked to four histone H3 tails with two different conformations using Haddock 2.4 and ClusPro. The best four models for each combination are selected and a two-sample t-test is performed to compare the binding affinities of helical conformations vs. linear conformations using Prodigy. The protein-protein interactions are examined using LigPlot.

Results

The linear histone conformations in predicted INGs-histone H3 complexes exhibit statistically significant higher binding affinity than their helical counterparts (confidence level of 99%). The outputs of predicted models generated by the molecular docking programs Haddock 2.4 and ClusPro are comparable, and the obtained protein-protein interaction patterns are consistent with experimentally confirmed binding patterns.

Conclusion

The results show that the conformation of the histone tail is significantly affecting the binding affinity of the docking protein. Herewith, this in silico study demonstrated in detail the binding preference of the ING protein family to histone H3 tail. Further research on the effect of certain PTMs on the final tail conformation and the interaction between those factors seem to be promising for a better understanding of epigenetics.

Protein-protein interaction prediction with deep learning: A comprehensive review

Most proteins perform their biological function by interacting with themselves or other molecules. Thus, one may obtain biological insights into protein functions, disease prevalence, and therapy development by identifying protein-protein interactions (PPI). However, finding the interacting and non-interacting protein pairs through experimental approaches is labour-intensive and time-consuming, owing to the variety of proteins. Hence, protein-protein interaction and protein-ligand binding problems have drawn attention in the fields of bioinformatics and computer-aided drug discovery. Deep learning methods paved the way for scientists to predict the 3-D structure of proteins from genomes, predict the functions and attributes of a protein, and modify and design new proteins to provide desired functions. This review focuses on recent deep learning methods applied to problems including predicting protein functions, protein-protein interaction and their sites, protein-ligand binding, and protein design.

Precise Assembly of Proteins and Carbohydrates for Next-Generation Biomaterials

The complexity and diversity of biomacromolecules make them a unique class of building blocks for generating precise assemblies. They are particularly available to a new generation of biomaterials integrated with living systems due to their intrinsic properties such as accurate recognition, self-organization, and adaptability. Therefore, many excellent approaches have been developed, leading to a variety of quite practical outcomes. Here, we review recent advances in the fabrication and application of artificially precise assemblies by employing proteins and carbohydrates as building blocks, followed by our perspectives on some of new challenges, goals, and opportunities for the future research directions in this field.