Hierarchical Flexible Peptide Docking by Conformer Generation and Ensemble Docking of Peptides

Investigation of the MDM2-binding potential of de novo designed peptides using enhanced sampling simulations

The tumor suppressor p53 plays a crucial role in cellular responses to various stresses, regulating key processes such as apoptosis, senescence, and DNA repair. Dysfunctional p53, prevalent in approximately 50 % of human cancers, contributes to tumor development and resistance to treatment. This study employed deep learning-based protein design and structure prediction methods to identify novel high-affinity peptide binders (Pep1 and Pep2) targeting MDM2, with the aim of disrupting its interaction with p53. Extensive all-atom molecular dynamics simulations highlighted the stability of the designed peptide in complex with the target, supported by several structural analyses, including RMSD, RMSF, Rg, SASA, PCA, and free energy landscapes. Using the steered molecular dynamics and umbrella sampling simulations, we elucidate the dissociation dynamics of p53, Pep1, and Pep2 from MDM2. Notable differences in interaction profiles were observed, emphasizing the distinct dissociation patterns of each peptide. In conclusion, the results of our umbrella sampling simulations suggest Pep1 as a higher-affinity MDM2 binder compared to p53 and Pep2, positioning it as a potential inhibitor of the MDM2-p53 interaction. Using state-of-the-art protein design tools and advanced MD simulations, this study provides a comprehensive framework for rational in silico design of peptide binders with therapeutic implications in disrupting MDM2-p53 interactions for anticancer interventions.

Antimicrobial Peptide Recognition Motif of the Substrate Binding Protein SapA from Nontypeable Haemophilus influenzae

Nontypeable Haemophilus influenzae (NTHi) is an opportunistic pathogen associated with respiratory diseases, including otitis media and exacerbations of chronic obstructive pulmonary disease. NTHi exhibits resistance to killing by host antimicrobial peptides (AMPs) mediated by SapA, the substrate binding protein of the sensitivity to antimicrobial peptides (Sap) transporter. However, the specific mechanisms by which SapA selectively binds various AMPs such as defensins and cathelicidin are unknown. In this study, we report mutational analyses of both defensin AMPs and the SapA binding pocket to define the specificity of AMP recognition. Bactericidal assays revealed that NTHi lacking SapA are more susceptible to human beta defensins and LL-37, while remaining highly resistant to a human alpha defensin. In contrast to homologues, our research underscores the distinct specificity of NTHi SapA, which selectively recognizes and binds to peptides containing the charged-hydrophobic motif PKE and RRY. These findings provide valuable insight into the divergence of SapA among bacterial species and NTHi SapA's ability to selectively interact with specific AMPs to mediate resistance.

Shotgun proteomics and molecular simulations on multifunctional bioactive peptides derived from the whey of unexplored "Gaddi" goat of Himalayas

The very first time, whey protein from the Himalayan goat breed "Gaddi" was hydrolyzed with alcalase, flavourzyme, and a combination of both in this study. The degree of hydrolysis (DH) ranged from 28 to 53%, with sequential hydrolysis by combination achieving the highest DH. The sequential hydrolysis demonstrated antimicrobial activity against all pathogens used with 3 kDa permeate showed significantly higher (p < 0.05) activity against S. aureus, E. coli, B. cereus and C. sakazakii. The antioxidant activity was in the range of IC50 = 0.49 to 2.00 mg protein/mL, flavourzyme and sequential hydrolysates showed significant ABTS radical and FRAP inhibition. The α-amylase inhibitory activity was highest in 3 kDa permeate of flavourzyme with IC50 values of 0.34 mg protein/mL. Bioactive peptides DDSPDLPK, EMPFPK and TPEVDKEALEK were identified most significant in the hydrolysates. In molecular docking, the DDSPDLPK interacted most efficiently with enzymes involved in microbial growth, oxidative stress, and hyperglycemia.

Assessing the fast prediction of peptide conformers and the impact of non-natural modifications

We present an assessment of different approaches to predict peptide structures using modeling tools. Several small molecule, protein, and peptide-focused methodologies were used for the fast prediction of conformers for peptides shorter than 30 amino acids. We assessed the effect of including restraints based on annotated or predicted secondary structure motifs. A number of peptides in bound conformations and in solution were collected to compare the tools. In addition, we studied the impact of changing single amino acids to non-natural residues using molecular dynamics simulations. Deep learning methods such as AlphaFold2, or the combination of physics-based approaches with secondary structure information, produce the most accurate results for natural sequences. In the case of peptides with non-natural modifications, modeling the peptide containing natural amino acids first and then modifying and simulating the peptide using benchmarked force fields is a recommended pipeline. The results can guide the modeling of oligopeptides for drug discovery projects.

Biocomputational characterisation of MBO_200107 protein of Mycobacterium tuberculosis variant caprae: a molecular docking and simulation study

The principal objective of this study was to delineate the potentiality of the MBO_200107 protein from the Mycobacterium tuberculosis variant caprae in cancer research. It is a cytoplasmic protein, comprised of a 354-long amino acid chain, alkaline, had a molecular weight of 39089.37 Da, an isoelectric point of 9.62 and a grand average of hydropathicity of -0.345. One of the functional domains was predicted as Gammaglutamylcyclotransferase (GGCT). Among tertiary structures, the Modeller and Phyre2 model satisfied all the quality parameters, though they are truncated; contrarily, the I-TASSER model is full length and contains the sequence for the GGCT domain, though it did not meet all the quality parameters. It also has significant sequence similarities (47.5% by EMBOSS Water and 72.4% by EMBOSS Matcher) with a human GGCT, and the conserved sequences are confined to the GGCT domain of the MBO_200107. According to molecular docking analyses, the protein has a binding affinity of -4.8 kcal/mol by Autodock Vina and -56.465 kcal/mol by HPEPDOCK to the human glutathione (GSH), an essential metabolite for GGCT metabolism. The Molecular dynamic simulation of the docked complex showed the binding efficiency of the GSH to MBO_200107 with a minimal structural alteration. The in silico findings mentioned above revealed that the protein could be used as a supplementary tool in cancer research, such as designing vaccines or drugs where the role of GGCT has been implicated. Further, we recommend fully characterising the protein and conducting essential in vitro and in vivo experiments to determine its detailed usefulness.Communicated by Ramaswamy H. Sarma.

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.

Selective events at individual sites underlie the evolution of monkeypox virus clades

In endemic regions (West Africa and the Congo Basin), the genetic diversity of monkeypox virus (MPXV) is geographically structured into two major clades (Clades I and II) that differ in virulence and host associations. Clade IIb is closely related to the B.1 lineage, which is dominating a worldwide outbreak initiated in 2022. Lineage B.1 has however accumulated mutations of unknown significance that most likely result from apolipoprotein B mRNA editing catalytic polypeptide-like 3 (APOBEC3) editing. We applied a population genetics-phylogenetics approach to investigate the evolution of MPXV during historical viral spread in Africa and to infer the distribution of fitness effects. We observed a high preponderance of codons evolving under strong purifying selection, particularly in viral genes involved in morphogenesis and replication or transcription. However, signals of positive selection were also detected and were enriched in genes involved in immunomodulation and/or virulence. In particular, several genes showing evidence of positive selection were found to hijack different steps of the cellular pathway that senses cytosolic DNA. Also, a few selected sites in genes that are not directly involved in immunomodulation are suggestive of antibody escape or other immune-mediated pressures. Because orthopoxvirus host range is primarily determined by the interaction with the host immune system, we suggest that the positive selection signals represent signatures of host adaptation and contribute to the different virulence of Clade I and II MPXVs. We also used the calculated selection coefficients to infer the effects of mutations that define the predominant human MPXV1 (hMPXV1) lineage B.1, as well as the changes that have been accumulating during the worldwide outbreak. Results indicated that a proportion of deleterious mutations were purged from the predominant outbreak lineage, whose spread was not driven by the presence of beneficial changes. Polymorphic mutations with a predicted beneficial effect on fitness are few and have a low frequency. It remains to be determined whether they have any significance for ongoing virus evolution.

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.

HCovDock: an efficient docking method for modeling covalent protein-ligand interactions

Covalent inhibitors have received extensive attentions in the past few decades because of their long residence time, high binding efficiency and strong selectivity. Therefore, it is valuable to develop computational tools like molecular docking for modeling of covalent protein-ligand interactions or screening of potential covalent drugs. Meeting the needs, we have proposed HCovDock, an efficient docking algorithm for covalent protein-ligand interactions by integrating a ligand sampling method of incremental construction and a scoring function with covalent bond-based energy. Tested on a benchmark containing 207 diverse protein-ligand complexes, HCovDock exhibits a significantly better performance than seven other state-of-the-art covalent docking programs (AutoDock, Cov_DOX, CovDock, FITTED, GOLD, ICM-Pro and MOE). With the criterion of ligand root-mean-squared distance < 2.0 Å, HCovDock obtains a high success rate of 70.5% and 93.2% in reproducing experimentally observed structures for top 1 and top 10 predictions. In addition, HCovDock is also validated in virtual screening against 10 receptors of three proteins. HCovDock is computationally efficient and the average running time for docking a ligand is only 5 min with as fast as 1 sec for ligands with one rotatable bond and about 18 min for ligands with 23 rotational bonds. HCovDock can be freely assessed at http://huanglab.phys.hust.edu.cn/hcovdock/.

Rational Discovery of Microtubule-Stabilizing Peptides

Microtubule (MT) stabilization is an attractive pharmacological strategy to hamper the progress of neurodegenerative diseases. In this regard, seeking peptides with MT-stabilizing properties has awoken great interest. This work reports the rational discovery of two structurally related MT-stabilizing octapeptides using a combination of protein-peptide docking, conventional molecular dynamics, Gaussian accelerated molecular dynamics (GaMD), and tubulin polymerization assays. FASTA sequences for ∼1000 peptides were crafted from single and double mutants of davunetide (NAP) and docked against the Taxol (TX) site on an octameric MT model representing a portion of the MT wall. Docked peptides were rescored after MM minimization and binding free energy refinement through single-point MM/GBSA calculations. The 60 best-ranked peptides were subjected to 50 ns MD simulations on peptide-MT complexes at the terminal TX site in the octameric Tau-MT model resulting in 11 complexes with occupancies greater than 99% and peptide-protein binding free energies less than -40 kcal/mol. Selected peptides were then examined through 300 ns GaMD simulations in complexes containing two identical ligands at the terminal and intermediate TX sites in the Tau-MT model to account for the differential association of MT-binding peptides to different regions of the MT structure. Six candidates showed a favorable MT-binding potential based on the analysis of interaction frequencies and relative mobilities of the complex components, suggesting a pivotal role of Arg278, Gln281, and Arg369 residues for peptides recognition. Four candidates were predicted to preserve an adequate balance of longitudinal and lateral interactions between tubulin dimers in peptide-MT complexes such that MT-stabilizing effects could be expected. MT polymerization experiments confirmed that four peptides (HAPVSIHQ, NYPVSIHQ, NWPVSIWQ, HAPVSIIQ) exhibit MT-stabilizing activity in vitro with NWPVSIWQ (P43) and HAPVSIIQ (P52) being the most active. Tryptophan quenching assays verified that P43 and P52 bind to nonpolymeric tubulin, whereas viability experiments on HEK cells confirmed their safety to pursue future pharmacological studies. The results herein presented are valuable to making progress in the rational design of MT-stabilizing peptides.

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.

Ligand-based design of peptide entry inhibitors targeting the endosomal receptor binding site of filoviruses

Filoviruses enter cells through micropinocytosis and trafficking into the endosomes in which they bind to the receptor Niemann-Pick C1 protein (NPC1) for membrane fusion and entry into the cytoplasm. The endosomal receptor-binding is critical step for filovirus entry. Designing inhibitors to block receptor binding will prevent viral entry. Using available binding structural information from the co-crystal structures of the viral GP with the receptor NPC1 or with monoclonal antibodies, we have conducted structure-based design of peptide inhibitors to target the receptor binding site (RBS). The designed peptides were tested for their inhibition activity against pseudo-typed or replication-competent viruses in a cell-based assay. The results indicate that these peptides exhibited strong activities against both Ebola and Marburg virus infection. It is expected that these peptides can be further developed for therapeutic use to treat filovirus infection and combat the outbreaks.

Novel engineered nanobodies specific for N‐terminal region of alpha‐synuclein recognize Lewy‐body pathology and inhibit in‐vitro seeded aggregation and toxicity

Nanobodies (Nbs), the single‐domain antigen‐binding fragments of dromedary heavy‐chain antibodies (HCAb), are excellent candidates as therapeutic and diagnostic tools in synucleinopathies because of their small size, solubility and stability. Here, we constructed an immune nanobody library specific to the monomeric form of alpha‐synuclein (α‐syn). Phage display screening of the library allowed the identification of a nanobody, Nbα‐syn01, specific for α‐syn. Unlike previously developed nanobodies, Nbα‐syn01 recognized the N‐terminal region which is critical for in vitro and in vivo aggregation and contains many point mutations involved in early PD cases. The affinity of the monovalent Nbα‐syn01 and the engineered bivalent format BivNbα‐syn01 measured by isothermal titration calorimetry revealed unexpected results where Nbα‐syn01 and its bivalent format recognized preferentially α‐syn fibrils compared to the monomeric form. Nbα‐syn01 and BivNbα‐syn01 were also able to inhibit α‐syn‐seeded aggregation in vitro and reduced α‐syn‐seeded aggregation and toxicity in cells showing their potential to reduce α‐syn pathology. Moreover, both nanobody formats were able to recognize Lewy‐body pathology in human post‐mortem brain tissue from PD and DLB cases. Additionally, we present evidence through structural docking that Nbα‐syn01 binds the N‐terminal region of the α‐syn aggregated form. Overall, these results highlight the potential of Nbα‐syn01 and BivNbα‐syn01 in developing into a diagnostic or a therapeutic tool for PD and related disorders.

Anurans against SARS-CoV-2: A review of the potential antiviral action of anurans cutaneous peptides

At the end of 2019, in China, clinical signs and symptoms of unknown etiology have been reported in several patients whose sample sequencing revealed pneumonia caused by the SARS-CoV-2 virus. COVID-19 is a disease triggered by this virus, and in 2020, the World Health Organization declared it a pandemic. Since then, efforts have been made to find effective therapeutic agents against this disease. Identifying novel natural antiviral drugs can be an alternative to treatment. For this reason, antimicrobial peptides secreted by anurans' skin have gained attention for showing a promissory antiviral effect. Hence, this review aimed to elucidate how and which peptides secreted by anurans' skin can be considered therapeutic agents to treat or prevent human viral infectious diseases. Through a literature review, we attempted to identify potential antiviral frogs' peptides to combat COVID-19. As a result, the Magainin-1 and -2 peptides, from the Magainin family, the Dermaseptin-S9, from the Dermaseptin family, and Caerin 1.6 and 1.10, from the Caerin family, are molecules that already showed antiviral effects against SARS-CoV-2 in silico. In addition to these peptides, this review suggests that future studies should use other families that already have antiviral action against other viruses, such as Brevinins, Maculatins, Esculentins, Temporins, and Urumins. To apply these peptides as therapeutic agents, experimental studies with peptides already tested in silico and new studies with other families not tested yet should be considered.

Elucidating the Therapeutic Potential of Cell-Penetrating Peptides in Human Tenon Fibroblast Cells

Cell-penetrating peptides (CPPs) have been widely used as vehicles for delivering therapeutic molecules to the site of action. Apart from their delivering potential, the biological effects of CPPs have not been explored in detail. JTS-1 is a CPP that has been reported to have gene delivery functions, although its biological role is yet to be determined. Hence, in this study, we revealed the biological mechanism such as its uptake mechanism and immunogenic potential and function using primary human tenon fibroblast (TF) cells collected from patients undergoing glaucoma trabeculectomy surgery. Our results showed that the JTS-1 peptide has an α-helical structure and is nontoxic up to 1 μM concentration. It was found to be colocalized with early endosome (Rab5), recycling endosome (Rab7), and Rab11 and interacted with major histocompatibility complex (MHC) class I and II. The peptide also affected actin polymerization, which is regulated by cofilin phosphorylation and ROCK1 localization. It also inhibited TF cell proliferation. Therefore, the JTS-1 peptide could be used as a possible therapeutic agent for modifying the fibrosis process, where TF proliferation is a key cause of surgery failure.

Efficient 3D conformer generation of cyclic peptides formed by a disulfide bond

Cyclic peptides formed by disulfide bonds have been one large group of common drug candidates in drug development. Structural information of a peptide is essential to understand its interaction with its target. However, due to the high flexibility of peptides, it is difficult to sample the near-native conformations of a peptide. Here, we have developed an extended version of our MODPEP approach, named MODPEP2.0, to fast generate the conformations of cyclic peptides formed by a disulfide bond. MODPEP2.0 builds the three-dimensional (3D) structures of a cyclic peptide from scratch by assembling amino acids one by one onto the cyclic fragment based on the constructed rotamer and cyclic backbone libraries. Being tested on a data set of 193 diverse cyclic peptides, MODPEP2.0 obtained a considerable advantage in both accuracy and computational efficiency, compared with other sampling algorithms including PEP-FOLD, ETKDG, and modified ETKDG (mETKDG). MODPEP2.0 achieved a high sampling accuracy with an average C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α RMSD of 2.20 Å and 1.66 Å when 10 and 100 conformations were considered, respectively, compared with 3.41 Å and 2.62 Å for PEP-FOLD, 3.44 Å and 3.16 Å for ETKDG, 3.09 Å and 2.72 Å for mETKDG. MODPEP2.0 also reproduced experimental peptide structures for 81.35% of the test cases when an ensemble of 100 conformations were considered, compared with 54.95%, 37.50% and 50.00% for PEP-FOLD, ETKDG, and mETKDG. MODPEP2.0 is computationally efficient and can generate 100 peptide conformations in one second. MODPEP2.0 will be useful in sampling cyclic peptide structures and modeling related protein-peptide interactions, facilitating the development of cyclic peptide drugs.

Efficient Enzymatic Incorporation of Dehydroalanine Based on SAMDI-Assisted Identification of Optimized Tags for OspF/SpvC

Site-specific modification of proteins has important applications in biological research and drug development. Reactive tags such as azide, alkyne, and tetrazine have been used extensively to achieve the abovementioned goal. However, bulky side-chain "ligation scars" are often left after the labeling and may hinder the biological application of such engineered protein products. Conjugation chemistry via dehydroalanine (Dha) may provide an opportunity for "traceless" ligation because the activated alkene moiety on Dha can then serve as an electrophile to react with radicalophile, thiol/amine nucleophile, and reactive phosphine probe to introduce a minimal linker in the protein post-translational modifications. In this report, we present a mild and highly efficient enzymatic approach to incorporate Dha with phosphothreonine/serine lyases, OspF and SpvC. These lyases originally catalyze an irreversible elimination reaction that converts a doubly phosphorylated substrate with phosphothreonine (pT) or phosphoserine (pS) to dehydrobutyrine (Dhb) or Dha. To generate a simple monophosphorylated tag for these lyases, we conducted a systematic approach to profile the substrate specificity of OspF and SpvC using peptide arrays and self-assembled monolayers for matrix-assisted laser desorption/ionization mass spectrometry. The optimized tag, [F/Y/W]-pT/pS-[F/Y/W] (where [F/Y/W] indicates an aromatic residue), results in a ∼10-fold enhancement of the overall peptide labeling efficiency via Dha chemistry and enables the first demonstration of protein labeling as well as live cell labeling with a minimal ligation linker via enzyme-mediated incorporation of Dha.

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.

Screening and identification of a novel antidiabetic peptide from collagen hydrolysates of Chinese giant salamander skin: Network pharmacology, inhibition kinetics and protection of IR-HepG2 cells

In this study, a novel peptide GPPGPA was screened from the collagen hydrolysates of Chinese giant salamander (Andrias davidianus) skin, and anti-diabetes mechanism was predicted by network pharmacology, and inhibitory...

Computational Screening for the Anticancer Potential of Seed-Derived Antioxidant Peptides: A Cheminformatic Approach

Some seed-derived antioxidant peptides are known to regulate cellular modulators of ROS production, including those proposed to be promising targets of anticancer therapy. Nevertheless, research in this direction is relatively slow owing to the inevitable time-consuming nature of wet-lab experimentations. To help expedite such explorations, we performed structure-based virtual screening on seed-derived antioxidant peptides in the literature for anticancer potential. The ability of the peptides to interact with myeloperoxidase, xanthine oxidase, Keap1, and p47^phox was examined. We generated a virtual library of 677 peptides based on a database and literature search. Screening for anticancer potential, non-toxicity, non-allergenicity, non-hemolyticity narrowed down the collection to five candidates. Molecular docking found LYSPH as the most promising in targeting myeloperoxidase, xanthine oxidase, and Keap1, whereas PSYLNTPLL was the best candidate to bind stably to key residues in p47^phox. Stability of the four peptide-target complexes was supported by molecular dynamics simulation. LYSPH and PSYLNTPLL were predicted to have cell- and blood-brain barrier penetrating potential, although intolerant to gastrointestinal digestion. Computational alanine scanning found tyrosine residues in both peptides as crucial to stable binding to the targets. Overall, LYSPH and PSYLNTPLL are two potential anticancer peptides that deserve deeper exploration in future.

Exploration of antigenic determinants in spike glycoprotein of SARS-CoV2 and identification of five salient potential epitopes

Emerging pathogens have been an eternal threat to mankind. In a series of pandemics caused by notorious coronaviruses, a newly emerged SARS-CoV2 virus is creating panic among the world population. The unavailability of reliable theranostics insists the exploration of antigenic determinants in spike glycoprotein of SARS-CoV2. The four novel inserts ('70VSGTNGT76', '150KSWM153', 247SYLTPG252 and 674QTQTNSPRR682) in SARS-CoV2 spike protein were unraveled via multiple sequence alignment of spike proteins of SARS-CoV2, SARS-CoV, and MERS-CoV. The three-dimension (3D) modeling of the spike protein of the SARS-CoV2 and their interaction with the ACE2 receptor was delineated with the help of SWISS-MODEL and 3DLigandSite web servers. The predicted 3D model of SARS-CoV2 was further verified by SAVES, RAMPAGE, and ProSA-web tools. The potential B-cell immunogenic epitopes of SARS-CoV2 were predicted out by using various software viz. IEDB B-cell epitopes prediction tool, BepiPred linear epitope prediction tool, Emini Surface Accessibility Prediction tool, and Kolaskar-Tongaonkar antigenicity web tool. The five epitopes (i.e. '71SGTNGTKRFDN81, 247SYLTPG252, 634RVYST638, 675QTQTNSPRRARSV687, and 1054QSAPH1058) were selected as potent antigenic determinants. The quantum of information generated by this study will prove beneficial for the development of effective therapeutics, diagnostics, and multi-epitopic vaccines to combat this ongoing menace.

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.

A Bioluminescent Biosensor for Quantifying the Interaction of SARS-CoV-2 and Its Receptor ACE2 in Cells and In Vitro

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is currently spreading and mutating with increasing speed worldwide. Therefore, there is an urgent need for a simple, sensitive, and high-throughput (HTP) assay to quantify virus–host interactions in order to quickly evaluate the infectious ability of mutant viruses and to develop or validate virus-inhibiting drugs. Here, we developed an ultrasensitive bioluminescent biosensor to evaluate virus–cell interactions by quantifying the interaction between the SARS-CoV-2 receptor binding domain (RBD) and its cellular receptor angiotensin-converting enzyme 2 (ACE2) both in living cells and in vitro. We have successfully used this novel biosensor to analyze SARS-CoV-2 RBD mutants and evaluated candidate small molecules (SMs), antibodies, and peptides that may block RBD:ACE2 interaction. This simple, rapid, and HTP biosensor tool will significantly expedite the detection of viral mutants and the anti-COVID-19 drug discovery process.

Rational Design of Novel Glycomimetic Peptides for E-Selectin Targeting

E-selectin is a cell-adhesion receptor with specific recognition capacity toward sialo-fucosylated Lewis carbohydrates present in leukocytes and tumor cells. E-selectin interactions mediate the progress of inflammatory processes and tumor metastasis, which aroused the interest in using this protein as a biomolecular target to design glycomimetic inhibitors for active targeting or therapeutic purposes. In this work, we report the rational discovery of two novel glycomimetic peptides targeting E-selectin based on mutations of the reference selectin-binding peptide IELLQAR. Sixteen single or double mutants at Ile1, Leu3, Leu4, and Arg7 residues were evaluated as potential candidates for E-selectin targeting using 50 ns molecular dynamics (MD) simulations. Nine peptides showing a stable association with the functional pocket were modified by adding a cysteine residue to the N-terminus to confer versatility for further chemical conjugation. Subsequent 50 ns MD simulations resulted in five cysteine-modified peptides with retained or improved E-selectin binding potential. Then, 300 ns accelerated MD (aMD) simulations were used to examine the binding properties of the best five cysteine-modified peptides. CIEELQAR and CIELFQAR exhibit the most selective association with the functional pocket of E-selectin, as revealed by potential of mean force profiles. Microscale thermophoresis experiments confirmed the E-selectin binding capacity of the selected peptides with K_D values in the low micromolar range (CIEELQAR K_D = 35.0 ± 1.4 μM; CIELFQAR K_D = 16.4 ± 0.7 μM), which are 25-fold lower than the reported value for the native ligand sLe^x (K_D = 878 μM). Our findings support the potential of CIEELQAR and CIELFQAR as novel E-selectin-targeting peptides with high recognition capacity and versatility for chemical conjugation, which are critical for enabling future applications in active targeting.

Introduction of a Glycine Linker Connecting the Heavy and Light Chains in Synthetic Cardosin B-Derived Rennet Changes the Specificity of Subpocket S3'

The development of plant-based synthetic rennets is of high commercial interest, due to the current great consumer demand for animal product alternatives. A previously developed recombinant form of the aspartic protease cardosin B with a three-glycine linker showed great potential due to its good performance in milk coagulation. This enzyme was found to be more specific and less proteolytically active than the native form for milk clotting, but the underlying structural causes for these activity changes were not completely clear. Here, we have performed molecular dynamics simulations with the recombinant enzyme with and without the linker. Our results showed that the introduction of the linker changes the subpocket S3', which is located more than 4 nm away. These results showcase how small modifications in proteins can have significant effects in distant regions in the protein structure that affect their biotechnological applications.

A Novel Machine Learning Strategy for the Prediction of Antihypertensive Peptides Derived from Food with High Efficiency

Strategies to screen antihypertensive peptides with high throughput and rapid speed will doubtlessly contribute to the treatment of hypertension. Food-derived antihypertensive peptides can reduce blood pressure without side effects. In the present study, a novel model based on the eXtreme Gradient Boosting (XGBoost) algorithm was developed and compared with the dominating machine learning models. To further reflect on the reliability of the method in a real situation, the optimized XGBoost model was utilized to predict the antihypertensive degree of the k-mer peptides cutting from six key proteins in bovine milk, and the peptide-protein docking technology was introduced to verify the findings. The results showed that the XGBoost model achieved outstanding performance, with an accuracy of 86.50% and area under the receiver operating characteristic curve of 94.11%, which were better than the other models. Using the XGBoost model, the prediction of antihypertensive peptides derived from milk protein was consistent with the peptide-protein docking results, and was more efficient. Our results indicate that using the XGBoost algorithm as a novel auxiliary tool is feasible to screen for antihypertensive peptides derived from food, with high throughput and high efficiency.

Atlasing of Assembly Landscapes using Distance Geometry and Graph Rigidity

This Article describes a novel geometric methodology for analyzing free energy and kinetics of assembly driven by short-range pair-potentials in an implicit solvent and provides a proof-of-concept illustration of its unique capabilities. An atlas is a labeled partition of the assembly landscape into a roadmap of maximal, contiguous, nearly-equipotential-energy conformational regions or macrostates, together with their neighborhood relationships. The new methodology decouples the roadmap generation from sampling and produces: (1) a queryable atlas of local potential energy minima, their basin structure, energy barriers, and neighboring basins; (2) paths between a specified pair of basins, each path being a sequence of conformational regions or macrostates below a desired energy threshold; and (3) approximations of relative path lengths, basin volumes (configurational entropy), and path probabilities. Results demonstrating the core algorithm's capabilities and high computational efficiency have been generated by a resource-light, curated open source software implementation EASAL (Efficient Atlasing and Search of Assembly Landscapes, ACM Trans. Math. Softw. 2018 44, 1-48. 10.1145/3204472; see software, Efficient Atlasing and Search of Assembly Landscapes, 2016. https://bitbucket.org/geoplexity/easal; video, Video Illustrating the opensource software EASAL, 2016. https://cise.ufl.edu/~sitharam/EASALvideo.mpeg; and user guide, EASAL software user guide, 2016. https://bitbucket.org/geoplexity/easal/src/master/CompleteUserGuide.pdf). Running on a laptop with Intel(R) Core(TM) i7-7700@3.60 GHz CPU with 16GB of RAM, EASAL atlases several hundred thousand conformational regions or macrostates in minutes using a single compute core. Subsequent path and basin computations each take seconds. A parallelized EASAL version running on the same laptop with 4 cores gives a 3× speedup for atlas generation. The core algorithm's correctness, time complexity, and efficiency-accuracy trade-offs are formally guaranteed using modern distance geometry, geometric constraint systems and combinatorial rigidity. The methodology further links the shape of the input assembling units to a type of intuitive and queryable bar-code of the output atlas, which in turn determine stable assembled structures and kinetics. This succinct input-output relationship facilitates reverse analysis and control toward design. A novel feature that is crucial to both the high sampling efficiency and decoupling of roadmap generation from sampling is a recently developed theory of convex Cayley (distance-based) custom parametrizations specific to assembly, as opposed to folding. Representing microstates with macrostate-specific Cayley parameters, to generate microstate samples, avoids gradient-descent search used by all prevailing methods. Further, these parametrizations convexify conformational regions or macrostates. This ratchets up sampling efficiency, significantly reducing number of repeated and discarded samples. These features of the new stand-alone methodology can also be used to complement the strengths of prevailing methodologies including Molecular Dynamics, Monte Carlo, and Fast Fourier Transform based methods.

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

MOTIVATION

Protein-peptide interactions mediate a wide variety of cellular and biological functions. Methods for predicting these interactions have garnered a lot of interest over the past few years, as witnessed by the rapidly growing number of peptide-based therapeutic molecules currently in clinical trials. The size and flexibility of peptides has shown to be challenging for existing automated docking software programs.

RESULTS

Here we present AutoDock CrankPep or ADCP in short, a novel approach to dock flexible peptides into rigid receptors. ADCP folds a peptide in the potential field created by the protein to predict the protein-peptide complex. We show that it outperforms leading peptide docking methods on two protein-peptide datasets commonly used for benchmarking docking methods: LEADS-PEP and peptiDB, comprised of peptides with up to 15 amino acids in length. Beyond these datasets, ADCP reliably docked a set of protein-peptide complexes containing peptides ranging in lengths from 16 to 20 amino acids. The robust performance of ADCP on these longer peptides enables accurate modeling of peptide-mediated protein-protein interactions and interactions with disordered proteins.

AVAILABILITY AND IMPLEMENTATION

ADCP is distributed under the LGPL 2.0 open source license and is available at http://adcp.scripps.edu. The source code is available at https://github.com/ccsb-scripps/ADCP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Docking: Challenges, Advances and its Use in Drug Discovery Perspective

Molecular docking is a process through which small molecules are docked into the macromolecular structures for scoring its complementary values at the binding sites. It is a vibrant research area with dynamic utility in structure-based drug-designing, lead optimization, biochemical pathway and for drug designing being the most attractive tools. Two pillars for a successful docking experiment are correct pose and affinity prediction. Each program has its own advantages and drawbacks with respect to their docking accuracy, ranking accuracy and time consumption so a general conclusion cannot be drawn. Moreover, users don't always consider sufficient diversity in their test sets which results in certain programs to outperform others. In this review, the prime focus has been laid on the challenges of docking and troubleshooters in existing programs, underlying algorithmic background of docking, preferences regarding the use of docking programs for best results illustrated with examples, comparison of performance for existing tools and algorithms, state of art in docking, recent trends of diseases and current drug industries, evidence from clinical trials and post-marketing surveillance are discussed. These aspects of the molecular drug designing paradigm are quite controversial and challenging and this review would be an asset to the bioinformatics and drug designing communities.

Improving Protein-Peptide Docking Results via Pose-Clustering and Rescoring with a Combined Knowledge-Based and MM-GBSA Scoring Function

Protein-peptide docking, which predicts the complex structure between a protein and a peptide, is a valuable computational tool in peptide therapeutics development and the mechanistic investigation of peptides involved in cellular processes. Although current peptide docking approaches are often able to sample near-native peptide binding modes, correctly identifying those near-native modes from decoys is still challenging because of the extremely high complexity of the peptide binding energy landscape. In this study, we have developed an efficient postdocking rescoring protocol using a combined scoring function of knowledge-based ITScorePP potentials and physics-based MM-GBSA energies. Tested on five benchmark/docking test sets, our postdocking strategy showed an overall significantly better performance in binding mode prediction and score-rmsd correlation than original docking approaches. Specifically, our postdocking protocol outperformed original docking approaches with success rates of 15.8 versus 10.5% for pepATTRACT on the Global_57 benchmark, 5.3 versus 5.3% for CABS-dock on the Global_57 benchmark, 17.0 versus 11.3% for FlexPepDock on the LEADS-PEP data set, 40.3 versus 33.9% for HPEPDOCK on the Local_62 benchmark, and 64.2 versus 52.8% for HPEPDOCK on the LEADS-PEP data set when the top prediction was considered. These results demonstrated the efficacy and robustness of our postdocking protocol.

Peptide Combination Generator: a Tool for Generating Peptide Combinations

Peptides are used as reagents both for basic research and diagnostic purposes. Therefore, there is a need for novel methods for the design of peptide molecules with a particular specific physiochemical profile. The properties of the peptides are governed by the nature of amino acids constituting the peptide. There is a lack of a web server or tools which could predict all the possible combinations of the peptides generated because of the combinations of amino acids based on the physiochemical properties. We have developed a peptide combination generator (PepCoGen), a web server for generating all the possible combinations of peptides by varying the amino acids having similar physiochemical properties at a particular position. It also predicts other properties of the peptides including molecular weight, charge, solubility, hydrophobic plot, and isoelectric point, and random three-dimensional models for each generated combination.

Challenges and opportunities of automated protein-protein docking: HDOCK server vs human predictions in CAPRI Rounds 38-46

Protein-protein docking plays an important role in the computational prediction of the complex structure between two proteins. For years, a variety of docking algorithms have been developed, as witnessed by the critical assessment of prediction interactions (CAPRI) experiments. However, despite their successes, many docking algorithms often require a series of manual operations like modeling structures from sequences, incorporating biological information, and selecting final models. The difficulties in these manual steps have significantly limited the applications of protein-protein docking, as most of the users in the community are nonexperts in docking. Therefore, automated docking like a web server, which can give a comparable performance to human docking protocol, is pressingly needed. As such, we have participated in the blind CAPRI experiments for Rounds 38-45 and CASP13-CAPRI challenge for Round 46 with both our HDOCK automated docking web server and human docking protocol. It was shown that our HDOCK server achieved an "acceptable" or higher CAPRI-rated model in the top 10 submitted predictions for 65.5% and 59.1% of the targets in the docking experiments of CAPRI and CASP13-CAPRI, respectively, which are comparable to 66.7% and 54.5% for human docking protocol. Similar trends can also be observed in the scoring experiments. These results validated our HDOCK server as an efficient automated docking protocol for nonexpert users. Challenges and opportunities of automated docking are also discussed.

Highly Flexible Ligand Docking: Benchmarking of the DockThor Program on the LEADS-PEP Protein-Peptide Data Set

Protein-peptide interactions play a crucial role in many cellular and biological functions, which justify the increasing interest in the development of peptide-based drugs. However, predicting experimental binding modes and affinities in protein-peptide docking remains a great challenge for most docking programs due to some particularities of this class of ligands, such as the high degree of flexibility. In this paper, we present the performance of the DockThor program on the LEADS-PEP data set, a benchmarking set composed of 53 diverse protein-peptide complexes with peptides ranging from 3 to 12 residues and with up to 51 rotatable bonds. The DockThor performance for pose prediction on redocking studies was compared with some state-of-the-art docking programs that were also evaluated on the LEADS-PEP data set, AutoDock, AutoDock Vina, Surflex, GOLD, Glide, rDock, and DINC, as well as with the task-specific docking protocol HPepDock. Our results indicate that DockThor could dock 40% of the cases with an overall backbone RMSD below 2.5 Å when the top-scored docking pose was considered, exhibiting similar results to Glide and outperforming other protein-ligand docking programs, whereas rDock and HPepDock achieved superior results. Assessing the docking poses closest to the crystal structure (i.e., best-RMSD pose), DockThor achieved a success rate of 60% in pose prediction. Due to the great overall performance of handling peptidic compounds, the DockThor program can be considered as suitable for docking highly flexible and challenging ligands, with up to 40 rotatable bonds. DockThor is freely available as a virtual screening Web server at https://www.dockthor.lncc.br/ .

PepPro: A Nonredundant Structure Data Set for Benchmarking Peptide-Protein Computational Docking

We present a nonredundant benchmark, coined PepPro, for testing peptide-protein docking algorithms. Currently, PepPro contains 89 nonredundant experimentally determined peptide-protein complex structures, with peptide sequence lengths ranging from 5 to 30 amino acids. The benchmark covers peptides with distinct secondary structures, including helix, partial helix, a mixture of helix and β-sheet, β-sheet formed through binding, β-sheet formed through self-folding, and coil. In addition, unbound proteins' structures are provided for 58 complexes and can be used for testing the ability of a docking algorithm handling the conformational changes of proteins during the binding process. PepPro should benefit the docking community for the development and improvement of peptide docking algorithms. The benchmark is available at http://zoulab.dalton.missouri.edu/PepPro_benchmark. © 2019 Wiley Periodicals, Inc.

Assessing the performance of MM/PBSA and MM/GBSA methods. 9. Prediction reliability of binding affinities and binding poses for protein-peptide complexes

A significant number of protein-protein interactions (PPIs) are mediated through the interactions between proteins and peptide segments, and therefore determination of protein-peptide interactions (PpIs) is critical to gain an in-depth understanding of the PPI network and even design peptides or small molecules capable of modulating PPIs. Computational approaches, especially molecular docking, provide an efficient way to model PpIs, and a reliable scoring function that can recognize the correct binding conformations for protein-peptide complexes is one of the most important components in protein-peptide docking. The end-point binding free energy calculation methods, such as MM/GBSA and MM/PBSA, are theoretically more rigorous than most empirical and semi-empirical scoring functions designed for protein-peptide docking, but their performance in predicting binding affinities and binding poses for protein-peptide systems has not been systematically assessed. In this study, we first evaluated the capability of MM/GBSA and MM/PBSA with different solvation models, interior dielectric constants (εin) and force fields to predict the binding affinities for 53 protein-peptide complexes. For the 19 short peptides with 5-12 residues, MM/PBSA based on the minimized structures in explicit solvent with the ff99 force field and εin = 2 yields the best correlation between the predicted binding affinities and the experimental data (rp = 0.748), while for the 34 medium-size peptides with 20-25 residues, MM/GBSA based on 1 ns of molecular dynamics (MD) simulations in implicit solvent with the ff03 force field, the GBOBC1 model and a low interior dielectric constant (εin = 1) yields the best accuracy (rp = 0.735). Then, we assessed the rescoring capability of MM/PBSA and MM/GBSA to distinguish the correct binding conformations from the decoys for 112 protein-peptide systems. The results illustrate that MM/PBSA based on the minimized structures with the ff99 or ff14SB force field and MM/GBSA based on the minimized structures with the ff03 force field show excellent capability to recognize the near-native binding poses for the short and medium-size peptides, respectively, and they outperform the predictions given by two popular protein-peptide docking algorithms (pepATTRACT and HPEPDOCK). Therefore, MM/PBSA and MM/GBSA are powerful tools to predict the binding affinities and identify the correct binding poses for protein-peptide systems.

Complex macrocycle exploration: parallel, heuristic, and constraint-based conformer generation using ForceGen

ForceGen is a template-free, non-stochastic approach for 2D to 3D structure generation and conformational elaboration for small molecules, including both non-macrocycles and macrocycles. For conformational search of non-macrocycles, ForceGen is both faster and more accurate than the best of all tested methods on a very large, independently curated benchmark of 2859 PDB ligands. In this study, the primary results are on macrocycles, including results for 431 unique examples from four separate benchmarks. These include complex peptide and peptide-like cases that can form networks of internal hydrogen bonds. By making use of new physical movements (“flips” of near-linear sub-cycles and explicit formation of hydrogen bonds), ForceGen exhibited statistically significantly better performance for overall RMS deviation from experimental coordinates than all other approaches. The algorithmic approach offers natural parallelization across multiple computing-cores. On a modest multi-core workstation, for all but the most complex macrocycles, median wall-clock times were generally under a minute in fast search mode and under 2 min using thorough search. On the most complex cases (roughly cyclic decapeptides and larger) explicit exploration of likely hydrogen bonding networks yielded marked improvements, but with calculation times increasing to several minutes and in some cases to roughly an hour for fast search. In complex cases, utilization of NMR data to constrain conformational search produces accurate conformational ensembles representative of solution state macrocycle behavior. On macrocycles of typical complexity (up to 21 rotatable macrocyclic and exocyclic bonds), design-focused macrocycle optimization can be practically supported by computational chemistry at interactive time-scales, with conformational ensemble accuracy equaling what is seen with non-macrocyclic ligands. For more complex macrocycles, inclusion of sparse biophysical data is a helpful adjunct to computation.

Predicting protein-peptide interaction sites using distant protein complexes as structural templates

Protein-peptide interactions play an important role in major cellular processes, and are associated with several human diseases. To understand and potentially regulate these cellular function and diseases it is important to know the molecular details of the interactions. However, because of peptide flexibility and the transient nature of protein-peptide interactions, peptides are difficult to study experimentally. Thus, computational methods for predicting structural information about protein-peptide interactions are needed. Here we present InterPep, a pipeline for predicting protein-peptide interaction sites. It is a novel pipeline that, given a protein structure and a peptide sequence, utilizes structural template matches, sequence information, random forest machine learning, and hierarchical clustering to predict what region of the protein structure the peptide is most likely to bind. When tested on its ability to predict binding sites, InterPep successfully pinpointed 255 of 502 (50.7%) binding sites in experimentally determined structures at rank 1 and 348 of 502 (69.3%) among the top five predictions using only structures with no significant sequence similarity as templates. InterPep is a powerful tool for identifying peptide-binding sites; with a precision of 80% at a recall of 20% it should be an excellent starting point for docking protocols or experiments investigating peptide interactions. The source code for InterPred is available at http://wallnerlab.org/InterPep/ .