AutoDock CrankPep: combining folding and docking to predict protein-peptide complexes

Identification of bioactive compounds and potential mechanisms of scutellariae radix-coptidis rhizoma in the treatment of atherosclerosis by integrating network pharmacology and experimental validation

OBJECTIVE

This study aims at investigating the potential targets and functional mechanisms of Scutellariae Radix-Coptidis Rhizoma (QLYD) against atherosclerosis (AS) through network pharmacology, molecular docking, bioinformatic analysis and experimental validation.

METHODS

The compositions of QLYD were collected from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) and literature, where the main active components of QLYD and corresponding targets were identified. The potential therapeutic targets of AS were excavated using the OMIM database, DrugBank database, DisGeNET database, CTD database and GEO datasets. The protein-protein interaction (PPI) network of common targets was constructed and visualized by Cytoscape 3.7.2 software. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis were performed to analyze the function of core targets in the PPI network. Molecular docking was carried out using AutoDockTools, AutoDock Vina, and PyMOL software to verify the correlation between the main components of QLYD and the core targets. Mouse AS model was established and the results of network pharmacology were verified by in vivo experiments.

RESULTS

Totally 49 active components and 225 corresponding targets of QLYD were obtained, where 68 common targets were identified by intersecting with AS-related targets. Five hub genes including IL6, VEGFA, AKT1, TNF, and IL1B were screened from the PPI network. GO functional analysis reported that these targets had associations mainly with cellular response to oxidative stress, regulation of inflammatory response, epithelial cell apoptotic process, and blood coagulation. KEGG pathway analysis demonstrated that these targets were correlated to AGE-RAGE signaling pathway in diabetic complications, TNF signaling pathway, IL-17 signaling pathway, MAPK signaling pathway, and NF-kappa B signaling pathway. Results of molecular docking indicated good binding affinity of QLYD to FOS, AKT1, and TNF. Animal experiments showed that QLYD could inhibit inflammation, improve blood lipid levels and reduce plaque area in AS mice to prevent and treat AS.

CONCLUSION

QLYD may exert anti-inflammatory and anti-oxidative stress effects through multi-component, multi-target and multi-pathway to treat AS.

A study related to the treatment of gastric cancer with Xiang-Sha-Liu-Jun-Zi-Tang based on network analysis

Purpose

Xiang-Sha-Liu-Jun-Zi-Tang(XSLJZT) is a common formula for the treatment of Gastric Cancer(GC) and is widely used in clinical practice, however, there is a lack of investigation into its mechanism.

Methods

We collected and organized drug and disease targets, constructed the "XSLJZT-Active Ingredient-Target" visualization network, and performed GO and KEGG functional enrichment analysis of crossover genes, followed by molecular docking of active ingredients and core targets. The best docked monomers were combined with weighted gene co-expression network analysis(WGCNA) and macroscopically analyzed by GO and KEGG enrichment techniques. The results of cluster gene difference analysis, ROC evaluation, and CIBERSORT immune infiltration analysis were evaluated and finally supported by cellular experiments.

Results

The main components of XSLJZT are quercetin, stigmasterol, and naringenin, effectively treat GC by targeting STAT3, TP53 and MAPK3, which are involved in IL-17, TNF and HIF-1 signaling pathways. The results of molecular docking showed that quercetin bound better to the core targets. We performed an in-depth analysis of this monomer and found that quercetin acts on the core targets of TP53, MMP9, TIMP1 and MYC, and is involved in two key signaling pathways, TNF and IL-17, thus effectively treating GC. The experimental results are consistent with our analysis that quercetin inhibits the proliferation of GC cells and promotes apoptosis, and TP53, MYC and TIMP1 are the quercetin targets for the treatment of GC.

Conclusion

The present study tentatively suggests that quercetin, the main active ingredient in XSLJZT, can exert a therapeutic effect on GC by targeting TIMP1.

Design of Cyclic Peptides Targeting Protein-Protein Interactions Using AlphaFold

More than 930,000 protein-protein interactions (PPIs) have been identified in recent years, but their physicochemical properties differ from conventional drug targets, complicating the use of conventional small molecules as modalities. Cyclic peptides are a promising modality for targeting PPIs, but it is difficult to predict the structure of a target protein-cyclic peptide complex or to design a cyclic peptide sequence that binds to the target protein using computational methods. Recently, AlphaFold with a cyclic offset has enabled predicting the structure of cyclic peptides, thereby enabling de novo cyclic peptide designs. We developed a cyclic peptide complex offset to enable the structural prediction of target proteins and cyclic peptide complexes and found AlphaFold2 with a cyclic peptide complex offset can predict structures with high accuracy. We also applied the cyclic peptide complex offset to the binder hallucination protocol of AfDesign, a de novo protein design method using AlphaFold, and we could design a high predicted local-distance difference test and lower separated binding energy per unit interface area than the native MDM2/p53 structure. Furthermore, the method was applied to 12 other protein-peptide complexes and one protein-protein complex. Our approach shows that it is possible to design putative cyclic peptide sequences targeting PPI.

Design of Tetra-Peptide Ligands of Antibody Fc Regions Using In Silico Combinatorial Library Screening

Peptides, or short chains of amino-acid residues, are becoming increasingly important as active ingredients of drugs and as crucial probes and/or tools in medical, biotechnological, and pharmaceutical research. Situated at the interface between small molecules and larger macromolecular systems, they pose a difficult challenge for computational methods. We report an in silico peptide library generation and prioritization workflow using CmDock for identifying tetrapeptide ligands that bind to Fc regions of antibodies that is analogous to known in vitro recombinant peptide libraries' display and expression systems. The results of our in silico study are in accordance with existing scientific literature on in vitro peptides that bind to antibody Fc regions. In addition, we postulate an evolving in silico library design workflow that will help circumvent the combinatorial problem of in vitro comprehensive peptide libraries by focusing on peptide subunits that exhibit favorable interaction profiles in initial in silico peptide generation and testing.

A Non-G-Quadruplex DNA Aptamer Targeting NCL for Diagnosis and Therapy in Bladder Cancer

Bladder cancer (BC) is a highly aggressive malignant tumor affecting the urinary system, characterized by metastasis and a poor prognosis that often leads to limited therapeutic success. This study aims to develop a novel DNA aptamer for the diagnosis and treatment of BC using a tissue-based systematic evolution of ligands by an exponential enrichment (SELEX) process. By using SELEX, this work successfully generates a new aptamer named TB-5, which demonstrates a remarkable and specific affinity for nucleolin (NCL) in BC tissues and displays marked biocompatibility both in vitro and in vivo. Additionally, this work shows that NCL is a reliable tissue-specific biomarker in BC. Moreover, according to circular dichroism spectroscopy, TB-5 forms a non-G-quadruplex structure, distinguishing it from the current NCL-targeting aptamer AS1411, and exhibits a distinct binding region on NCL compared to AS1411. Notably, this study further reveals that TB-5 activates NCL function by promoting autophagy and suppressing the migration and invasion of BC cells, which occurs by disrupting mRNA transcription processes. These findings highlight the critical role of NCL in the pathological examination of BC and warrant more comprehensive investigations on anti-NCL aptamers in BC imaging and treatment.

Computational insight into stability-enhanced systems of anthocyanin with protein/peptide

Anthocyanins, which belong to the flavonoid group, are commonly found in the organs of plants native to South and Central America. However, these pigments are unstable under conditions of varying pH, heat, etc., which limits their potential applications. One method for preserving the stability of anthocyanins is through encapsulation using proteins or peptides. Nevertheless, the complex and diverse structure of these molecules, as well as the limitation of experimental technologies, have hindered a comprehensive understanding of the encapsulation processes and the mechanisms by which stability is enhanced. To address these challenges, computational methods, such as molecular docking and molecular dynamics simulation have been used to study the binding affinity and dynamics of interactions between proteins/peptides and anthocyanins. This review summarizes the mechanisms of interaction between these systems, based on computational approaches, and highlights the role of proteins and peptides in the stability enhancement of anthocyanins. It also discusses the current limitations of these methods and suggests possible solutions.

Accelerating the Screening of Small Peptide Ligands by Combining Peptide-Protein Docking and Machine Learning

This research introduces a novel pipeline that couples machine learning (ML), and molecular docking for accelerating the process of small peptide ligand screening through the prediction of peptide-protein docking. Eight ML algorithms were analyzed for their potential. Notably, Light Gradient Boosting Machine (LightGBM), despite having comparable F1-score and accuracy to its counterparts, showcased superior computational efficiency. LightGBM was used to classify peptide-protein docking performance of the entire tetrapeptide library of 160,000 peptide ligands against four viral envelope proteins. The library was classified into two groups, 'better performers' and 'worse performers'. By training the LightGBM algorithm on just 1% of the tetrapeptide library, we successfully classified the remaining 99%with an accuracy range of 0.81-0.85 and an F1-score between 0.58-0.67. Three different molecular docking software were used to prove that the process is not software dependent. With an adjustable probability threshold (from 0.5 to 0.95), the process could be accelerated by a factor of at least 10-fold and still get 90-95% concurrence with the method without ML. This study validates the efficiency of machine learning coupled to molecular docking in rapidly identifying top peptides without relying on high-performance computing power, making it an effective tool for screening potential bioactive compounds.

Predicting Protein-Peptide Interactions: Benchmarking Deep Learning Techniques and a Comparison with Focused Docking

The accurate prediction of protein structures achieved by deep learning (DL) methods is a significant milestone and has deeply impacted structural biology. Shortly after its release, AlphaFold2 has been evaluated for predicting protein-peptide interactions and shown to significantly outperform RoseTTAfold as well as a conventional blind docking method: PIPER-FlexPepDock. Since then, new AlphaFold2 models, trained specifically to predict multimeric assemblies, have been released and a new ab initio folding model OmegaFold has become available. Here, we assess docking success rates for these new DL folding models and compare their performance with our state-of-the-art, focused peptide-docking software AutoDock CrankPep (ADCP). The evaluation is done using the same dataset and performance metric for all methods. We show that, for a set of 99 nonredundant protein-peptide complexes, the new AlphaFold2 model outperforms other Deep Learning approaches and achieves remarkable docking success rates for peptides. While the docking success rate of ADCP is more modest when considering the top-ranking solution only, it samples correct solutions for around 62% of the complexes. Interestingly, different methods succeed on different complexes, and we describe a consensus docking approach using ADCP and AlphaFold2, which achieves a remarkable 60% for the top-ranking results and 66% for the top 5 results for this set of 99 protein-peptide complexes.

Computational design of cyclic peptides to inhibit protein-peptide interactions

Many protein-protein interactions result from the binding of one folded protein with one short peptide segment, such as complexes formed by SH3 or PDZ domains. These transient protein-peptide interactions are notably involved in cellular signaling pathways and generally have low affinities, which opens the possibility to design competitive inhibitors of these complexes. We present and assess here our computational approach, called Des3PI, to design de novo cyclic peptides with potential high affinity for protein surfaces involved in interactions with peptide segments. The results were not conclusive for two receptors, the αVβ3 integrin and the CXCR4 chemokine receptor, but were promising in the case of SH3 and PDZ domains: For the former, Des3PI was able to find at least one cyclic sequence with six hotspots that binds a SH3 domain with a better theoretical affinity to the known proline-rich RLP2 peptide. For the latter, Des3PI could identify at least four cyclic sequences with four or five hotspots that have lower binding free energies computed by the MM-PBSA method than the reference peptide GKAP.

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.

Advancing our knowledge of antigen processing with computational modelling, structural biology, and immunology

Antigen processing is an immunological mechanism by which intracellular peptides are transported to the cell surface while bound to Major Histocompatibility Complex molecules, where they can be surveyed by circulating CD8+ or CD4+ T-cells, potentially triggering an immunological response. The antigen processing pathway is a complex multistage filter that refines a huge pool of potential peptide ligands derived from protein degradation into a smaller ensemble for surface presentation. Each stage presents unique challenges due to the number of ligands, the polymorphic nature of MHC and other protein constituents of the pathway and the nature of the interactions between them. Predicting the ensemble of displayed peptide antigens, as well as their immunogenicity, is critical for improving T cell vaccines against pathogens and cancer. Our predictive abilities have always been hindered by an incomplete empirical understanding of the antigen processing pathway. In this review, we highlight the role of computational and structural approaches in improving our understanding of antigen processing, including structural biology, computer simulation, and machine learning techniques, with a particular focus on the MHC-I pathway.

Pomegranate Peel in the Amelioration of High-Altitude Disease: A Network Pharmacology and Molecular Docking Study of Underlying Mechanisms

High-altitude disease (HAD) describes the failure to adapt to the lack of oxygen found at high altitudes and therapeutic antioxidant effects have been attributed to pomegranate peel (PP) extract. Network pharmacology, molecular docking, and experimental validation were used to study mechanisms responsible for the alleviation of HAD by PP. The aim was to establish a reference for future research and aid technological development, particularly in clinical settings. Network pharmacology analysis showed that PP affected many targets in HAD via the active ingredients, luteolin 7-O-glycoside, punicalagin, and ellagic acid. HNRNPA1, HSPA1B, HSPA1A, CUL4B, CLTC, PPP1CA, PARP1, RACK1, NEDD8, and MAP3K1 were all targets, responsible for effects on ribosomes, apoptosis, cell cycle, mRNA surveillance pathway, and the MAPK signaling pathway. PP had an antiapoptosis effect on H9c2 cells damaged by hypoxia, as shown by annexinV-FITC/PI double staining. Practical Applications. HAD comprises a group of diseases caused by failure to adapt to a low-oxygen environment. PP extract has previously been shown to have antioxidant effects. PP attenuated damage to H9c2 cells and reduced the apoptosis rate. The current results lay the foundation for further experimental investigations.

A retro-inverso modified peptide alleviated ovalbumin-induced asthma model by affecting glycerophospholipid and purine metabolism of immune cells

Allergic asthma is a heterogeneous disease involving a variety of inflammatory cells. Immune imbalance or changes in the immune microenvironment are the essential causes that promote inflammation in allergic asthma. Tetraspanin CD81 can be used as a platform for receptor clustering and signal transmission owing to its special transmembrane structure and is known to participate in the physiological processes of cell proliferation, differentiation, adhesion, and migration. Previous studies have shown that CD81-targeting peptidomimetics exhibit anti-allergic lung inflammation. However, due to the low metabolic stability of peptide drugs, their druggability is limited. Here, we aimed to generate a metabolically stable anti-CD81 peptide, evaluate its anti-inflammatory action and establish its mechanism of action. Based on previous reports, we applied retro-inverse peptide modification to obtain a new compound, PD00 (NH2-D-Gly-D-Ser-D-Thr-D-Tyr-D-Thr-D-Gln-D-Gly-COOH), with high metabolic stability. Enhanced ultraperformance liquid chromatography-tandem mass spectrometry was used to investigate the in vitro and in vivo metabolic stabilities of PD00. The affinities of PD00 and CD81 were studied using molecular docking and surface plasmon resonance techniques. An ovalbumin (OVA)-induced asthma model was used to evaluate the effects of PD00 in vivo. Mice were treated with different concentrations of PD00 (175 and 350 μg/kg) for 10 days. Airway hyperresponsiveness (AHR) to acetyl-β-methacholine (Mch), inflammatory cell counts in the bronchoalveolar lavage fluid, and serum OVA-specific IgE levels were detected in the mice at the end of the experiment. Lung tissues were collected for haematoxylin and eosin staining, untargeted metabolomic analysis, and single-cell transcriptome sequencing. PD00 has a high affinity for CD81; therefore, administration of PD00 markedly ameliorated AHR and airway inflammation in mice after OVA sensitisation and exposure. Serum OVA-specific IgE levels decreased considerably. In addition, PD00 treatment increased glycerophospholipid and purine metabolism in immune cells. Collectively, PD00 may regulate the glycerophospholipid and purine metabolism pathways to ameliorate the pathophysiological features of asthma. These findings suggest that PD00 is a potential compound for the treatment of asthma.

tLyp-1: A peptide suitable to target NRP-1 receptor

Targeting vascular endothelial growth factor receptor (VEFGR) and its co-receptor neuropilin-1 (NRP-1) is an interesting vascular strategy. tLyp-1 is a tumor-homing and penetrating peptide of 7 amino acids (CGNKRTR). It is a truncated form of Lyp-1 (CGNKRTRGC), which is known to target NRP-1 receptor, with high affinity and specificity. It is mediated by endocytosis via C-end rule (CendR) internalization pathway. The aim of this study is to evaluate the importance of each amino acid in the tLyp-1 sequence through alanine-scanning (Ala-scan) technique, during which each of the amino acid in the sequence was systematically replaced by alanine to produce 7 different analogues. In silico approach through molecular docking and molecular dynamics are employed to understand the interaction between the peptide and its analogues with the NRP-1 receptor, followed by in vitro ligand binding assay study. The C-terminal Arg is crucial in the interaction of tLyp-1 with NRP-1 receptor. Substituting this residue dramatically reduces the affinity of this peptide which is clearly seen in this study. Lys-4 is also important in the interaction, which is confirmed via the in vitro study and the MM-PBSA analysis. The finding in this study supports the CendR, in which the presence of R/K-XX-R/K motif is essential in the binding of a ligand with NRP-1 receptor. This presented work will serve as a guide in the future work pertaining the development of active targeting agent towards NRP-1 receptor.

Disruption of 3CLpro protease self-association by short peptides as a potential route to broad spectrum coronavirus inhibitors

Coronaviruses have posed a persistent threat to human health over the last two decades. Despite the accumulated knowledge about coronavirus-related pathogens, development of an effective treatment for its new variant COVID-19 is highly challenging. For the highly-conserved and main coronavirus protease 3CLpro, dimerization is known to be essential for its catalytic activity and thereby for virus proliferation. Here, we assess the potential of short peptide segments to disrupt dimerization of the 3CLpro protease as a route to block COVID-19 proliferation. Based on the X-ray structure of the 3CLpro dimer, we identified the SPSGVY₁₂₆QCAMRP dodecapeptide segment as overlapping the hotspot regions on the 3CLpro dimer interface. Using computational blind docking of the peptide to the 3CLpro monomer, we found that the SPSGVY₁₂₆QCAMRP peptide has favourable thermodynamic binding (ΔG= -5.93 kcal/mol) to the hotspot regions at the 3CLpro dimer interface. Importantly, the peptide was also found to preferentially bind to the hotspot regions compared to other potential binding sites lying away from the dimer interface (ΔΔG=-1.31 kcal/mol). Docking of peptides corresponding to systematic mutation of the V125 and Y126 residues led to the identification of seven peptides, SPSGHAQCAMRP, SPSGVTQCAMRP, SPSGKPQCAMRP, SPSGATQCAMRP, SPSGWLQCAMRP, SPSGAPQCAMRP and SPSGHPQCAMRP, that outperform the wild-type SPSGVY₁₂₆QCAMRP peptide in terms of preferential binding to the 3CLpro dimer interface. These peptides have the potential to disrupt 3CLpro dimerization and therefore could provide lead structures for the development of broad spectrum COVID-19 inhibitors.Communicated by Ramaswamy H. Sarma.

Molecular Mechanism of the Saposhnikovia divaricata–Angelica dahurica Herb Pair in Migraine Therapy Based on Network Pharmacology and Molecular Docking

Objective. This work studied the molecular mechanism of the Saposhnikovia divaricata–Angelica dahurica herb pair (SAHP) in migraine treatment. Methods. The active ingredients of drugs were screened, and potential targets were predicted by the Traditional Chinese Medicine Systems Pharmacology (TCMSP), TCMID, ETCM, and other databases. Migraine-related targets were obtained by harnessing the GeneCards, DrugBank, OMIM, TTD, and other databases. The protein-protein interaction (PPI) network was constructed with STRING software by performing a Venn analysis with bioinformatics. Gene Ontology (GO) functional enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed with the Metascape platform. The component-target-pathway (C-T-P) network was constructed with Cytoscape 3.7.2 software, and molecular docking was assessed with AutoDockVina software. Results. A total of 183 relevant targets and 39 active ingredients in migraine therapy were obtained from SAHP. The active ingredients and targets were screened according to topological parameters: wogonin, anomalin, imperatorin, prangenin, 2-linoleoylglycerol, and methylenetanshinquinone were identified as key active ingredients. PTGS2, PIK3CA, PIK3CB, PIK3CD, F2, and AR were identified as key targets. The molecular docking results demonstrated high binding activity between the key active ingredients and key targets. A total of 20 important signaling pathways, including neural signaling pathways, calcium signaling pathways, pathways in cancer, cAMP signaling pathways, and PI3K-Akt signaling pathways, were obtained through enrichment analysis. Conclusion. Migraine with SAHP is mainly treated through anti-inflammatory and analgesic effects. The herb pair can be used for migraine using “multicomponent, multitarget, and multipathway” approaches.

Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations

The structural characterization of protein-peptide interactions is fundamental to elucidating biological processes and designing peptide drugs. Molecular dynamics (MD) simulations are extensively used to study biomolecular systems. However, simulating the protein-peptide binding process is usually quite expensive. Based on our previous studies, herein, we propose a simple and effective method to predict the binding site and pose of the peptide simultaneously using high-temperature (high-T) MD simulations with the RSFF2C force field. Thousands of binding events (nonspecific or specific) can be sampled during microseconds of high-T MD. From density-based clustering analysis, the structures of all of the 12 complexes (nine with linear peptides and three with cyclic peptides) can be successfully predicted with root-mean-square deviation (RMSD) < 2.5 Å. By directly simulating the process of the ligand binding onto the receptor, our method approaches experimental precision for the first time, significantly surpassing previous protein-peptide docking methods in terms of accuracy.

Docking cyclic peptides formed by a disulfide bond through a hierarchical strategy

MOTIVATION

Cyclization is a common strategy to enhance the therapeutic potential of peptides. Many cyclic peptide drugs have been approved for clinical use, in which the disulfide-driven cyclic peptide is one of the most prevalent categories. Molecular docking is a powerful computational method to predict the binding modes of molecules. For protein-cyclic peptide docking, a big challenge is considering the flexibility of peptides with conformers constrained by cyclization.

RESULTS

Integrating our efficient peptide 3D conformation sampling algorithm MODPEP2.0 and knowledge-based scoring function ITScorePP, we have proposed an extended version of our hierarchical peptide docking algorithm, named HPEPDOCK2.0, to predict the binding modes of the peptide cyclized through a disulfide against a protein. Our HPEPDOCK2.0 approach was extensively evaluated on diverse test sets and compared with the state-of-the-art cyclic peptide docking program AutoDock CrankPep (ADCP). On a benchmark dataset of 18 cyclic peptide-protein complexes, HPEPDOCK2.0 obtained a native contact fraction of above 0.5 for 61% of the cases when the top prediction was considered, compared with 39% for ADCP. On a larger test set of 25 cyclic peptide-protein complexes, HPEPDOCK2.0 yielded a success rate of 44% for the top prediction, compared with 20% for ADCP. In addition, HPEPDOCK2.0 was also validated on two other test sets of 10 and 11 complexes with apo and predicted receptor structures, respectively. HPEPDOCK2.0 is computationally efficient and the average running time for docking a cyclic peptide is about 34 min on a single CPU core, compared with 496 min for ADCP. HPEPDOCK2.0 will facilitate the study of the interaction between cyclic peptides and proteins and the development of therapeutic cyclic peptide drugs.

AVAILABILITY AND IMPLEMENTATION

http://huanglab.phys.hust.edu.cn/hpepdock/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Des3PI: a fragment-based approach to design cyclic peptides targeting protein-protein interactions

Protein-protein interactions (PPIs) play crucial roles in many cellular processes and their deregulation often leads to cellular dysfunctions. One promising way to modulate PPIs is to use peptide derivatives that bind their protein target with high affinity and high specificity. Peptide modulators are often designed using secondary structure mimics. However, fragment-based design is an alternative emergent approach in the PPI field. Most of the reported computational fragment-based libraries targeting PPIs are composed of small molecules or already approved drugs, but, according to our knowledge, no amino acid based library has been reported yet. In this context, we developed a novel fragment-based approach called Des3PI (design of peptides targeting protein-protein interactions) with a library composed of natural amino acids. All the amino acids are docked into the target surface using Autodock Vina. The resulting binding modes are geometrically clustered, and, in each cluster, the most recurrent amino acids are identified and form the hotspots that will compose the designed peptide. This approach was applied on Ras and Mcl-1 proteins, as well as on A[Formula: see text] protofibril. For each target, at least five peptides generated by Des3PI were tested in silico: the peptides were first blindly docked on their target, and then, the stability of the successfully docked complexes was verified using 200 ns MD simulations. Des3PI shows very encouraging results by yielding at least 3 peptides for each protein target that succeeded in passing the two-step assessment.

Identification of new aptamer BC-3 targeting RPS7 from rapid screening for bladder carcinoma

Aptamers, short single DNA or RNA oligonucleotides, have shown immense application potential as molecular probes for the early diagnosis and therapy of cancer. However, conventional cell-SELEX technologies for aptamer discovery are time-consuming and laborious. Here we discovered a new aptamer BC-3 by using an improved rapid X-Aptamer selection process for human bladder carcinoma, for which there is no specific molecular probe yet. We show that BC-3 exhibited excellent affinity in bladder cancer cells but not normal cells. We demonstrate that BC-3 displayed high selectivity for tumor cells over their normal counterparts in vitro, in mice, and in patient tumor tissue specimens. Further endocytosis pathway analysis revealed that BC-3 internalized into bladder cancer cells via clathrin-mediated endocytosis. Importantly, we identified ribosomal protein S7 (RPS7) as the binding target of BC-3 via an integrated methodology (mass spectrometry, colocalization assay, and immunoblotting). Together, we report that a novel aptamer BC-3 is discovered for bladder cancer and its properties in the disease are unearthed. Our findings will facilitate the discovery of novel diagnostic and therapeutic strategies for bladder cancer.

Target-based in-silico screening of basil polysaccharides against different epigenetic targets responsible for breast cancer

PURPOSE

Breast cancer (BC) is one of the leading types of cancer found in women. One of the causes reported for BC is improper regulation of epigenetic modifications. Various epigenetic targets such as histone deacetylases (HDAC) and histone acetyltransferases (HAT) regulate many types of cancer, including BC. Basil is known to possess anti-cancer properties; however, the role of its polysaccharides against different epigenetic targets is still not very clear. Therefore, the molecular docking method is used to find out the binding potential of the BPSs against different epigenetic targets responsible for BC.

METHODS

All the basil polysaccharides (BPSs) were screened against the diverse epigenetic targets reported for BC (HDAC1-2, 4-8, and HAT) using molecular docking studies alongwith swissADME studies to check the drug likeliness of the BPSs.

RESULTS

It was found that glucosamine ring, glucosamine linear, glucuronic acid linear, rhamnose linear, glucuronic acid ring, galactose ring, mannose, glucose, and xylose were exhibited consistent binding potential against the epigenetic targets (HDAC1, HDAC2, HDAC4, HDAC5, HDAC6, HDAC7, HDAC8, and HAT,) responsible for BC.

CONCLUSION

This is the first report where BPSs were reported against these epigenetic targets. These studies can help to understand the underlying mechanism of BPSs used against epigenetic targets for BC. These results can be further validated experimentally to confirm their potential as a promising inhibitor against the epigenetic targets (HDAC1-2, 4-8, and HAT) having a role in BC.

Solubility-Aware Protein Binding Peptide Design Using AlphaFold

New protein-protein interactions (PPIs) are identified, but PPIs have different physicochemical properties compared with conventional targets, making it difficult to use small molecules. Peptides offer a new modality to target PPIs, but designing appropriate peptide sequences by computation is challenging. Recently, AlphaFold and RoseTTAFold have made it possible to predict protein structures from amino acid sequences with ultra-high accuracy, enabling de novo protein design. We designed peptides likely to have PPI as the target protein using the "binder hallucination" protocol of AfDesign, a de novo protein design method using AlphaFold. However, the solubility of the peptides tended to be low. Therefore, we designed a solubility loss function using solubility indices for amino acids and developed a solubility-aware AfDesign binder hallucination protocol. The peptide solubility in sequences designed using the new protocol increased with the weight of the solubility loss function; moreover, they captured the characteristics of the solubility indices. Moreover, the new protocol sequences tended to have higher affinity than random or single residue substitution sequences when evaluated by docking binding affinity. Our approach shows that it is possible to design peptide sequences that can bind to the interface of PPI while controlling solubility.

A Study of Type II ɛ-PL Degrading Enzyme (pldII) in Streptomyces albulus through the CRISPRi System

ε-Poly-L-lysine (ε-PL) is a widely used antibacterial peptide polymerized of 25–35 L-lysine residues. The antibacterial effect of ε-PL is closely related to the polymerization degree. However, the mechanism of ε-PL degradation in S. albulus remains unclear. This study utilized the integrative plasmid pSET152-based CRISPRi system to transcriptionally repress the ε-PL degrading enzyme (pldII). The expression of pldII is regulated by changing the recognition site of dCas9. Through the ε-PL bacteriostatic experiments of repression strains, it was found that the repression of pldII improves the antibacterial effect of the ε-PL product. The consecutive MALDI-TOF-MS results confirmed that the molecular weight distribution of the ε-PL was changed after repression. The repression strain S1 showed a particular peak with a polymerization degree of 44, and other repression strains also generated ε-PL with a polymerization degree of over 40. Furthermore, the homology modeling and substrate docking of pldII, a typical endo-type metallopeptidase, were performed to resolve the degradation mechanism of ε-PL in S. albulus. The hydrolysis of ε-PL within pldII, initiated from the N-terminus by two amino acid-binding residues, Thr194 and Glu281, led to varying levels of polymerization of ε-PL.

Cyclization and Docking Protocol for Cyclic Peptide-Protein Modeling Using HADDOCK2.4

An emerging class of therapeutic molecules are cyclic peptides with over 40 cyclic peptide drugs currently in clinical use. Their mode of action is, however, not fully understood, impeding rational drug design. Computational techniques could positively impact their design, but modeling them and their interactions remains challenging due to their cyclic nature and their flexibility. This study presents a step-by-step protocol for generating cyclic peptide conformations and docking them to their protein target using HADDOCK2.4. A dataset of 30 cyclic peptide-protein complexes was used to optimize both cyclization and docking protocols. It supports peptides cyclized via an N- and C-terminus peptide bond and/or a disulfide bond. An ensemble of cyclic peptide conformations is then used in HADDOCK to dock them onto their target protein using knowledge of the binding site on the protein side to drive the modeling. The presented protocol predicts at least one acceptable model according to the critical assessment of prediction of interaction criteria for each complex of the dataset when the top 10 HADDOCK-ranked single structures are considered (100% success rate top 10) both in the bound and unbound docking scenarios. Moreover, its performance in both bound and fully unbound docking is similar to the state-of-the-art software in the field, Autodock CrankPep. The presented cyclization and docking protocol should make HADDOCK a valuable tool for rational cyclic peptide-based drug design and high-throughput screening.

Accelerating AutoDock Vina with GPUs

AutoDock Vina is one of the most popular molecular docking tools. In the latest benchmark CASF-2016 for comparative assessment of scoring functions, AutoDock Vina won the best docking power among all the docking tools. Modern drug discovery is facing a common scenario of large virtual screening of drug hits from huge compound databases. Due to the seriality characteristic of the AutoDock Vina algorithm, there is no successful report on its parallel acceleration with GPUs. Current acceleration of AutoDock Vina typically relies on the stack of computing power as well as the allocation of resource and tasks, such as the VirtualFlow platform. The vast resource expenditure and the high access threshold of users will greatly limit the popularity of AutoDock Vina and the flexibility of its usage in modern drug discovery. In this work, we proposed a new method, Vina-GPU, for accelerating AutoDock Vina with GPUs, which is greatly needed for reducing the investment for large virtual screens and also for wider application in large-scale virtual screening on personal computers, station servers or cloud computing, etc. Our proposed method is based on a modified Monte Carlo using simulating annealing AI algorithm. It greatly raises the number of initial random conformations and reduces the search depth of each thread. Moreover, a classic optimizer named BFGS is adopted to optimize the ligand conformations during the docking progress, before a heterogeneous OpenCL implementation was developed to realize its parallel acceleration leveraging thousands of GPU cores. Large benchmark tests show that Vina-GPU reaches an average of 21-fold and a maximum of 50-fold docking acceleration against the original AutoDock Vina while ensuring their comparable docking accuracy, indicating its potential for pushing the popularization of AutoDock Vina in large virtual screens.

Rules of Chinese Herbal Intervention of Radiation Pneumonia Based on Network Pharmacology and Data Mining

Objective

To explore the mechanism and principles of traditional Chinese medicine (TCM) in the management of radiation pneumonia.

Methods

The targets of radiation pneumonia were obtained by screening the GeneCards, OMIM, TTD, DrugBank, and HERB databases, analyzing ADME parameters. In addition, compounds and Chinese herbs that can act on the targets were screened from the TCMSP database. The core target compounds for TCM were used to construct the target-compound, compound-traditional Chinese medicine, and target-compound-traditional Chinese medicine networks. These networks were further used to select the core targets, compounds, and TCM. The binding strength between the core targets and compounds was determined using AutoDock Vina. The trajectory for the molecular dynamics simulation was completed by Desmond version 2020.

Results

A total of 55 active targets in radiation pneumonia were identified. Subsequently, 137 candidate compounds and 469 Chinese herbs were matched. Frequency statistics showed that the Chinese herbs that could interfere with radiation pneumonia were mainly bitter, spicy, and sweet, with both cold and warm properties. Moreover, they mainly belonged to liver and lung channels. The core targets included TNF, IL-6, TGF-β1, and TP53. The most important components were quercetin, resveratrol, and (-)-epigallocatechin-3-gallate. Moreover, the most significant traditional Chinese herbs were Perilla pueraria, ephedra, Lonicerae japonicae, and sea buckthorn. Furthermore, analysis of 222 sets of receptor-ligand docking results suggested that the compounds had good docking activity to their core targets. By combining the docking binding energy, we determined that the chemical compounds had strong binding energy to the targets.

Conclusion

Using network pharmacology, we explored the potential mechanism of TCM in the treatment of radiation pneumonia. The general rules for application of TCM in the treatment of radiation pneumonia were summarized. This study provides baseline information for future research on the development of TCM for the management of radiation pneumonia.

Integrated unsupervised-supervised modeling and prediction of protein-peptide affinities at structural level

Cell signal networks are orchestrated directly or indirectly by various peptide-mediated protein-protein interactions, which are normally weak and transient and thus ideal for biological regulation and medicinal intervention. Here, we develop a general-purpose method for modeling and predicting the binding affinities of protein-peptide interactions (PpIs) at the structural level. The method is a hybrid strategy that employs an unsupervised approach to derive a layered PpI atom-residue interaction (ulPpI[a-r]) potential between different protein atom types and peptide residue types from thousands of solved PpI complex structures and then statistically correlates the potential descriptors with experimental affinities (KD values) over hundreds of known PpI samples in a supervised manner to create an integrated unsupervised-supervised PpI affinity (usPpIA) predictor. Although both the ulPpI[a-r] potential and usPpIA predictor can be used to calculate PpI affinities from their complex structures, the latter seems to perform much better than the former, suggesting that the unsupervised potential can be improved substantially with a further correction by supervised statistical learning. We examine the robustness and fault-tolerance of usPpIA predictor when applied to treat the coarse-grained PpI complex structures modeled computationally by sophisticated peptide docking and dynamics simulation. It is revealed that, despite developed solely based on solved structures, the integrated unsupervised-supervised method is also applicable for locally docked structures to reach a quantitative prediction but can only give a qualitative prediction on globally docked structures. The dynamics refinement seems not to change (or improve) the predictive results essentially, although it is computationally expensive and time-consuming relative to peptide docking. We also perform extrapolation of usPpIA predictor to the indirect affinity quantities of HLA-A*0201 binding epitope peptides and NHERF PDZ binding scaffold peptides, consequently resulting in a good and moderate correlation of the predicted KD with experimental IC50 and BLU on the two peptide sets, with Pearson's correlation coefficients Rp = 0.635 and 0.406, respectively.

Antimicrobial Peptides from Human Microbiome Against Multidrug Efflux Pump of Pseudomonas aeruginosa: a Computational Study

The excess use of antibiotics has led to the evolution of multidrug-resistant pathogenic strains causing worldwide havoc. These multidrug-resistant strains require potent inhibitors. Pseudomonas aeruginosa is a lead cause of nosocomial infections and also feature in the critical priority list of the world health organization (WHO) for the development of new antibiotics against their antimicrobial resistance. Antimicrobial peptides (AMPs) found in almost every life form from microorganisms to humans are known to defend their hosts against various pathogens. Owing to the diversity of the human microbiome, in this study, we have identified the cell-penetrating AMPs from the human microbiome and studied their inhibitory activity against the outer membrane protein OprM of the MexAB-OprM, a constitutively expressed multidrug efflux pump of the Ps. aeruginosa. Screening of the AMPs from the human microbiome resulted in the identification of 147 cell-penetrating AMPs (CPAMPs). The virtual screening of these CPAMPs against the OprM protein showed significant inhibitory results with the top docked AMP showing binding affinity exceeding -30 kcal/mol. The molecular dynamic simulation determined the interaction stabilities between the AMPs and the OprM at the binding site. Further, the residue interaction networks (RINs) are analyses to identify the inhibitory patterns. Later, these patterns were confirmed by MM-PBSA analysis suggesting that the AMPs are majorly stabilized by electrostatic interactions at the binding site. Thus, the high binding affinity and insights from the molecular interaction signify that the identified CPAMPs from the human microbiome can be further explored as inhibitory agents against multidrug-resistant Ps. aeruginosa.

Predicting Protein-Peptide Complex Structures by Accounting for Peptide Flexibility and the Physicochemical Environment

Predicting protein-peptide complex structures is crucial to the understanding of a vast variety of peptide-mediated cellular processes and to peptide-based drug development. Peptide flexibility and binding mode ranking are the two major challenges for protein-peptide complex structure prediction. Peptides are highly flexible molecules, and therefore, brute-force modeling of peptide conformations of interest in protein-peptide docking is beyond current computing power. Inspired by the fact that the protein-peptide binding process is like protein folding, we developed a novel strategy, named MDockPeP2, which tries to address these challenges using physicochemical information embedded in abundant monomeric proteins with an exhaustive search strategy, in combination with an integrated global search and a local flexible minimization method. Only the peptide sequence and the protein crystal structure are required. The method was systemically assessed using a newly constructed structural database of 89 nonredundant protein-peptide complexes with the peptide sequence length ranging from 5 to 29 in which about half of the peptides are longer than 15 residues. MDockPeP2 yielded a total success rate of 58.4% (70.8, 79.8%) for the bound docking (i.e., with the bound receptor and fully flexible peptides) and 19.0% (44.8, 70.7%) for the challenging unbound docking when top 10 (100, 1000) models were considered for each prediction. MDockPeP2 achieved significantly higher success rates on two other datasets, peptiDB and LEADS-PEP, which contain only short- and medium-size peptides (≤ 15 residues). For peptiDB, our method obtained a success rate of 62.0% for the bound docking and 35.9% for the unbound docking when the top 10 models were considered. For LEADS-PEP, MDockPeP2 achieved a success rate of 69.8% when the top 10 models were considered. The program is available at https://zougrouptoolkit.missouri.edu/mdockpep2/download.html.

Molecular Docking of SP40 Peptide towards Cellular Receptors for Enterovirus 71 (EV-A71)

Enterovirus 71 (EV-A71) is one of the predominant etiological agents of hand, foot and mouth disease (HMFD), which can cause severe central nervous system infections in young children. There is no clinically approved vaccine or antiviral agent against HFMD. The SP40 peptide, derived from the VP1 capsid of EV-A71, was reported to be a promising antiviral peptide that targeted the host receptor(s) involved in viral attachment or entry. So far, the mechanism of action of SP40 peptide is unknown. In this study, interactions between ten reported cell receptors of EV-A71 and the antiviral SP40 peptide were evaluated through molecular docking simulations, followed by in vitro receptor blocking with specific antibodies. The preferable binding region of each receptor to SP40 was predicted by global docking using HPEPDOCK and the cell receptor-SP40 peptide complexes were refined using FlexPepDock. Local molecular docking using GOLD (Genetic Optimization for Ligand Docking) showed that the SP40 peptide had the highest binding score to nucleolin followed by annexin A2, SCARB2 and human tryptophanyl-tRNA synthetase. The average GoldScore for 5 top-scoring models of human cyclophilin, fibronectin, human galectin, DC-SIGN and vimentin were almost similar. Analysis of the nucleolin-SP40 peptide complex showed that SP40 peptide binds to the RNA binding domains (RBDs) of nucleolin. Furthermore, receptor blocking by specific monoclonal antibody was performed for seven cell receptors of EV-A71 and the results showed that the blocking of nucleolin by anti-nucleolin alone conferred a 93% reduction in viral infectivity. Maximum viral inhibition (99.5%) occurred when SCARB2 was concurrently blocked with anti-SCARB2 and the SP40 peptide. This is the first report to reveal the mechanism of action of SP40 peptide in silico through molecular docking analysis. This study provides information on the possible binding site of SP40 peptide to EV-A71 cellular receptors. Such information could be useful to further validate the interaction of the SP40 peptide with nucleolin by site-directed mutagenesis of the nucleolin binding site.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

Insights on peptide topology in the computational design of protein ligands: the example of lysozyme binding peptides

Herein, we compared the ability of linear and cyclic peptides generated in silico to target different protein sites: internal pockets and solvent-exposed sites. We selected human lysozyme (HuL) as a model target protein combined with the computational evolution of linear and cyclic peptides. The sequence evolution of these peptides was based on the PARCE algorithm. The generated peptides were screened based on their aqueous solubility and HuL binding affinity. The latter was evaluated by means of scoring functions and atomistic molecular dynamics (MD) trajectories in water, which allowed prediction of the structural features of the protein-peptide complexes. The computational results demonstrated that cyclic peptides constitute the optimal choice for solvent exposed sites, while both linear and cyclic peptides are capable of targeting the HuL pocket effectively. The most promising binders found in silico were investigated experimentally by surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and electrospray ionization mass spectrometry (ESI-MS) techniques. All tested peptides displayed dissociation constants in the micromolar range, as assessed by SPR; however, both NMR and ESI-MS suggested multiple binding modes, at least for the pocket binding peptides. A detailed NMR analysis confirmed that both linear and cyclic pocket peptides correctly target the binding site they were designed for.

Divide-and-link peptide docking: a fragment-based peptide docking protocol

Protein-peptide interactions are crucial for various important cellular regulations, and are also a basis for understanding protein-protein interactions, protein folding and peptide drug design. Due to the limited structural data obtained using experimental methods, it is necessary to predict protein-peptide interaction modes using computational methods. In the present work, we designed a fragment-based docking protocol, Divide-and-Link Peptide Docking (DLPepDock), to predict protein-peptide binding modes. This protocol contains the following steps: dividing the peptide into fragments and separately docking the fragments using a third-party small molecular docking tool, linking the docked fragmental poses to form the whole peptide conformations via fragmental coordinate transformation using our in-house program, removing unreasonable poses according to several geometrical filters, extracting representative conformations after clustering for further minimization using the steepest descent and conjugation gradient methods based on a full-atom molecular force field and finally scoring using the MM/PBSA binding energy calculation implemented in Amber. When tested on the LEADS-PEP benchmark data set of 26 diverse complexes with peptides of 6-12 residues, FlexPepDock ab initio and AutoDock CrankPep achieved superior results. DLPepDock performed better than the other 15 docking protocols implemented in nine docking programs (HPepDock, DockThor, rDock, Glide, LeDock, AutoDock, AutoDock Vina, Surflex, and GOLD). The Linux scripts to call the third-party tools and run all the calculations.

A novel multi-objective metaheuristic algorithm for protein-peptide docking and benchmarking on the LEADS-PEP dataset

Protein-peptide interactions have attracted the attention of many drug discovery scientists due to their possible druggability features on most key biological activities such as regulating disease-related signaling pathways and enhancing the immune system's responses. Different studies have utilized some protein-peptide-specific docking algorithms/methods to predict protein-peptide interactions. However, the existing algorithms/methods suffer from two serious limitations which make them unsuitable for protein-peptide docking problems. First, it seems that the prevalent approaches require to be modified and remodeled for weighting the unbounded forces between a protein and a peptide. Second, they do not employ state-of-the-art search algorithms for detecting the 3D pose of a peptide relative to a protein. To address these restrictions, the present study aims to introduce a novel multi-objective algorithm, which first generates some potential 3D poses of a peptide, and then, improves them through its operators. The candidate solutions are further evaluated using Multi-Objective Pareto Front (MOPF) optimization concepts. To this end, van der Waals, electrostatic, solvation, and hydrogen bond energies between the atoms of a protein and designated peptide are computed. To evaluate the algorithm, it is first applied to the LEADS-PEP dataset containing 53 protein-peptide complexes with up to 53 rotatable branches/bonds and then compared with three popular/efficient algorithms. The obtained results indicate that the MOPF-based approaches which reduce the backbone RMSD between the original and predicted states, achieve significantly better results in terms of the success rate in predicting the near-native conditions. Besides, a comparison between the different types of search algorithms reveals that efficient ones like the multi-objective Trader/differential evolution algorithm can predict protein-peptide interactions better than the popular algorithms such as the multi-objective genetic/particle swarm optimization algorithms.

Molecular Modeling of Signal Peptide Recognition by Eukaryotic Sec Complexes

Here, we review recent molecular modelling and simulation studies of the Sec translocon, the primary component/channel of protein translocation into the endoplasmic reticulum (ER) and bacterial periplasm, respectively. Our focus is placed on the eukaryotic Sec61, but we also mention modelling studies on prokaryotic SecY since both systems operate in related ways. Cryo-EM structures are now available for different conformational states of the Sec61 complex, ranging from the idle or closed state over an inhibited state with the inhibitor mycolactone bound near the lateral gate, up to a translocating state with bound substrate peptide in the translocation pore. For all these states, computational studies have addressed the conformational dynamics of the translocon with respect to the pore ring, the plug region, and the lateral gate. Also, molecular simulations are addressing mechanistic issues of insertion into the ER membrane vs. translocation into the ER, how signal-peptides are recognised at all in the translocation pore, and how accessory proteins affect the Sec61 conformation in the co- and post-translational pathways.

Genetic modifications of metallothionein enhance the tolerance and bioaccumulation of heavy metals in Escherichia coli

Metallothioneins (MTs) are low molecular weight cysteine-rich proteins that bind to metals. Owing to their high cysteine (Cys) content, MTs are effective mediators of heavy metal detoxification. To enhance the heavy metal binding ability of MT from the freshwater crab Sinopotamon henanense (ShMT), sequence-based multiple sequence alignment (MSA) and structure-based molecular docking simulation (MDS) were conducted in order to identify amino acid residues that could be mutated to bolster such metal-binding activity. Site-directed mutagenesis was then used to modify the primary structure of ShMT, and the recombinant proteins were further enhanced using the SUMO fusion expression system to yield SUMO-ShMT1, SUMO-ShMT2, and SUMO-ShMT3 harboring one-, two-, and three- point mutations, respectively. The resultant modified proteins were primarily expressed in a soluble form and exhibited the ability to readily bind to heavy metals. Importantly, these modified proteins exhibited significantly enhanced heavy metal binding capacities, and they improved Cd²⁺, Cu²⁺ and Zn²⁺ tolerance and bioaccumulation in Escherichia coli (E. coli) in a manner dependent upon the number of introduced point mutations (SUMO-ShMT3 > SUMO-ShMT2 > SUMO-ShMT1 > SUMO-ShMT > control). Indeed, E. coli cells harboring the pET28a-SUMO-ShMT3 expression vector exhibited maximal Cd²⁺, Cu²⁺, and Zn²⁺ bioaccumulation that was increased by 1.86 ± 0.02-, 1.71 ± 0.03-, and 2.13 ± 0.02-fold relative to that in E. coli harboring the pET28a-SUMO-ShMT vector. The present study offers a basis for the preparation of genetically engineered bacteria that are better able to bioaccumulate and tolerate heavy metals, thus providing a foundation for biological heavy metal water pollution treatment.

A deep-learning framework for multi-level peptide-protein interaction prediction

Peptide-protein interactions are involved in various fundamental cellular functions and their identification is crucial for designing efficacious peptide therapeutics. Recently, a number of computational methods have been developed to predict peptide-protein interactions. However, most of the existing prediction approaches heavily depend on high-resolution structure data. Here, we present a deep learning framework for multi-level peptide-protein interaction prediction, called CAMP, including binary peptide-protein interaction prediction and corresponding peptide binding residue identification. Comprehensive evaluation demonstrated that CAMP can successfully capture the binary interactions between peptides and proteins and identify the binding residues along the peptides involved in the interactions. In addition, CAMP outperformed other state-of-the-art methods on binary peptide-protein interaction prediction. CAMP can serve as a useful tool in peptide-protein interaction prediction and identification of important binding residues in the peptides, which can thus facilitate the peptide drug discovery process.

Helical and β-Turn Conformations in the Peptide Recognition Regions of the VIM1 PHD Finger Abrogate H3K4 Peptide Recognition

The PHD finger-containing VARIANT IN METHYLATION/ORTHRUS (VIM/ORTH) family of proteins in Arabidopsis consists of functional homologues of mammalian UHRF1 and is required for the maintenance of DNA methylation. Comparison of the sequence with those of other PHD fingers implied that VIM1 and VIM3 PHD could recognize lysine 4 of histone H3 (H3K4) through interactions mediated by a conserved aspartic acid. However, our calorimetric and modified histone peptide array binding studies suggested that neither H3K4 nor other histone marks are recognized by VIM1 and VIM3 PHD fingers. Here, we report a 2.6 Å resolution crystal structure of the VIM1 PHD finger and demonstrate significant structural changes in the putative H3 recognition segments in contrast to canonical H3K4 binding PHD fingers. These changes include (i) the H3A1 binding region, (ii) strand β1 that forms an intermolecular β-sheet with the H3 peptide, and (iii) an aspartate-containing motif involved in salt bridge interaction with H3K4, which together appear to abrogate recognition of H3K4 by the VIM1 PHD finger. To understand the significance of the altered structural features in the VIM1 PHD that might prevent histone H3 recognition, we modeled a chimeric VIM1 PHD (chmVIM1 PHD) by grafting the peptide binding structural features of the BHC80 PHD onto the VIM1 PHD. Molecular dynamics simulation and metadynamics analyses revealed that the chmVIM1 PHD-H3 complex is stable and also showed a network of intermolecular interactions similar to those of the BHC80 PHD-H3 complex. Collectively, this study reveals that subtle structural changes in the peptide binding region of the VIM1 PHD abrogate histone H3 recognition.

Source and exploration of the peptides used to construct peptide-drug conjugates

Peptide-drug conjugates (PDCs) are a class of novel molecules widely designed and synthesized for delivering payload drugs. The peptide part plays a vital role in the whole molecule, because they determine the ability of the molecules to penetrate the membrane and target to the specific targets. Here, we introduce the source of different kinds of cell-penetrating peptides (CPPs) and cell-targeting peptides (CTPs) that have been used or could be used in constructing PDCs as well as their latest application in delivering drugs. What's more, the approaches of developing CPPs and CTPs and the techniques to discover novel peptides are focused on and summarized in the review. This review aims to help relevant researchers fast understand the research status of peptides in PDCs and carry forward the process of novel peptides discovery.

Improving Docking Power for Short Peptides Using Random Forest

In recent years, therapeutic peptides have gained a lot interest as demonstrated by the 60 peptides approved as drugs in major markets and 150+ peptides currently in clinical trials. However, while small molecule docking is routinely used in rational drug design efforts, docking peptides has proven challenging partly because docking scoring functions, developed and calibrated for small molecules, perform poorly for these molecules. Here, we present random forest classifiers trained to discriminate correctly docked peptides. We show that, for a testing set of 47 protein-peptide complexes, structurally dissimilar from the training set and previously used to benchmark AutoDock Vina's ability to dock short peptides, these random forest classifiers improve docking power from ∼25% for AutoDock scoring functions to an average of ∼70%. These results pave the way for peptide-docking success rates comparable to those of small molecule docking. To develop these classifiers, we compiled the ProptPep37_2021 data set, a curated, high-quality set of 322 crystallographic protein-peptides complexes annotated with structural similarity information. The data set also provides a collection of high-quality putative poses with a range of deviations from the crystallographic pose, providing correct and incorrect poses (i.e., decoys) of the peptide for each entry. The ProptPep37_2021 data set as well as the classifiers presented here are freely available.

Cyclic Peptides as Protein Kinase Inhibitors: Structure-Activity Relationship and Molecular Modeling

Under-expression or overexpression of protein kinases has been shown to be associated with unregulated cell signal transduction in cancer cells. Therefore, there is major interest in designing protein kinase inhibitors as anticancer agents. We have previously reported [WR]5, a peptide containing alternative arginine (R) and tryptophan (W) residues as a non-competitive c-Src tyrosine kinase inhibitor. A number of larger cyclic peptides containing alternative hydrophobic and positively charged residues [WR]x (x = 6-9) and hybrid cyclic-linear peptides, [R6K]W6 and [R5K]W7, containing R and W residues were evaluated for their protein kinase inhibitory potency. Among all the peptides, cyclic peptide [WR]9 was found to be the most potent tyrosine kinase inhibitor. [WR]9 showed higher inhibitory activity (IC50 = 0.21 μM) than [WR]5, [WR]6, [WR]7, and [WR]8 with IC50 values of 0.81, 0.57, 0.35, and 0.33 μM, respectively, against c-Src kinase as determined by a radioactive assay using [γ-33P]ATP. Consistent with the result above, [WR]9 inhibited other protein kinases such as Abl kinase activity with an IC50 value of 0.35 μM, showing 2.2-fold higher inhibition than [WR]5 (IC50 = 0.79 μM). [WR]9 also inhibited PKCa kinase activity with an IC50 value of 2.86 μM, approximately threefold higher inhibition than [WR]5 (IC50 = 8.52 μM). A similar pattern was observed against Braf, c-Src, Cdk2/cyclin A1, and Lck. [WR]9 exhibited IC50 values of <0.25 μM against Akt1, Alk, and Btk. These data suggest that [WR]9 is consistently more potent than other cyclic peptides with a smaller ring size and hybrid cyclic-linear peptides [R6K]W6 and [R5K]W7 against selected protein kinases. Thus, the presence of R and W residues in the ring, ring size, and the number of amino acids in the structure of the cyclic peptide were found to be critical in protein kinase inhibitory potency. We identified three putative binding pockets through automated blind docking of cyclic peptides [WR](5-9). The most populated pocket is located between the SH2, SH3, and N-lobe domains on the opposite side of the ATP binding site. The second putative pocket is formed by the same domains and located on the ATP binding site side of the protein. Finally, a third pocket was identified between the SH2 and SH3 domains. These results are consistent with the non-competitive nature of the inhibition displayed by these molecules. Molecular dynamics simulations of the protein-peptide complexes indicate that the presence of either [WR]5 or [WR]9 affects the plasticity of the protein and in particular the volume of the ATP binding site pocket in different ways. These results suggest that the second pocket is most likely the site where these peptides bind and offer a plausible rationale for the increased affinity of [WR]9.

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.

An affinity prediction approach for the ligand of E3 ligase Cbl-b and an insight into substrate binding pattern

Protein-protein interactions (PPIs) are essentially fundamental to all cellular processes, so that developing small molecule inhibitors of PPIs have great significance despite representing a huge challenge. Studying PPIs with the help of peptide motifs could obtain the structural information and reference significance to reduce the difficulty in the development of small molecules. Computational methods are powerful tools to characterize peptide-protein interactions, especially molecular dynamics simulation and binding free energy calculation. Here, we established an affinity prediction model suitable for Casitas B lymphoma-b (Cbl-b) and phosphorylated motif system. According to the affinity data set of multiple truncated peptides, the force field, solvent model, and internal dielectric constant of molecular mechanics/generalized Born surface area (MM/GBSA) method were optimized. Further, we predicted the affinity of the rationally designed new sequences through this model and obtained a new 6-mer motif with a 7-fold increase in affinity and the comprehensive structure-activity relationship. Moreover, we proposed an insight of unexpected activity of the truncated 5-mer peptide and revealed the possible binding mode of the new highly active 6-mer motif by extended simulation. Our results showed that the activity enhancement of the truncated peptide was caused by the acetyl-mediated conformation change. The side chain of Arg and pTyr in the 6-mer motif co-occupied the site p1 to form numerous hydrogen bond interactions and increased hydrophobic interaction formed with Tyr266, leading to the higher affinity. The present work provided a reference to investigate the PPI of Cbl-b and phosphorylated substrates and guided the development of Cbl-b inhibitors.

How does Sec63 affect the conformation of Sec61 in yeast?

The Sec complex catalyzes the translocation of proteins of the secretory pathway into the endoplasmic reticulum and the integration of membrane proteins into the endoplasmic reticulum membrane. Some substrate peptides require the presence and involvement of accessory proteins such as Sec63. Recently, a structure of the Sec complex from Saccharomyces cerevisiae, consisting of the Sec61 channel and the Sec62, Sec63, Sec71 and Sec72 proteins was determined by cryo-electron microscopy (cryo-EM). Here, we show by co-precipitation that the Sec61 channel subunit Sbh1 is not required for formation of stable Sec63-Sec61 contacts. Molecular dynamics simulations started from the cryo-EM conformation of Sec61 bound to Sec63 and of unbound Sec61 revealed how Sec63 affects the conformation of Sec61 lateral gate, plug, pore region and pore ring diameter via three intermolecular contact regions. Molecular docking of SRP-dependent vs. SRP-independent signal peptide chains into the Sec61 channel showed that the pore regions affected by presence/absence of Sec63 play a crucial role in positioning the signal anchors of SRP-dependent substrates nearby the lateral gate.

Binding Ensembles of p53-MDM2 Peptide Inhibitors by Combining Bayesian Inference and Atomistic Simulations

Designing peptide inhibitors of the p53-MDM2 interaction against cancer is of wide interest. Computational modeling and virtual screening are a well established step in the rational design of small molecules. But they face challenges for binding flexible peptide molecules that fold upon binding. We look at the ability of five different peptides, three of which are intrinsically disordered, to bind to MDM2 with a new Bayesian inference approach (MELD × MD). The method is able to capture the folding upon binding mechanism and differentiate binding preferences between the five peptides. Processing the ensembles with statistical mechanics tools depicts the most likely bound conformations and hints at differences in the binding mechanism. Finally, the study shows the importance of capturing two driving forces to binding in this system: the ability of peptides to adopt bound conformations (ΔGconformation) and the interaction between interface residues (ΔGinteraction).

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.

The AutoDock suite at 30

The AutoDock suite provides a comprehensive toolset for computational ligand docking and drug design and development. The suite builds on 30 years of methods development, including empirical free energy force fields, docking engines, methods for site prediction, and interactive tools for visualization and analysis. Specialized tools are available for challenging systems, including covalent inhibitors, peptides, compounds with macrocycles, systems where ordered hydration plays a key role, and systems with substantial receptor flexibility. All methods in the AutoDock suite are freely available for use and reuse, which has engendered the continued growth of a diverse community of primary users and third-party developers.