BactPepDB: a database of predicted peptides from a exhaustive survey of complete prokaryote genomes

Integrated sequence and -omic features reveal novel small proteome of Mycobacterium tuberculosis

Bioinformatic studies on small proteins are under-represented due to difficulties in annotation posed by their small size. However, recent discoveries emphasize the functional significance of small proteins in cellular processes including cell signaling, metabolism, and adaptation to stress. In this study, we utilized a Random Forest classifier trained on sequence features, RNA-Seq, and Ribo-Seq data to uncover small proteins (smORFs) in M. tuberculosis. Independent predictions for the exponential and starvation conditions resulted in 695 potential smORFs. We examined the functional implications of these smORFs using homology searches, LC-MS/MS, and ChIP-seq data, testing their expression in diverse growth conditions, and identifying protein domains. We provide evidence that some of these smORFs could be part of operons, or exist as upstream ORFs. This expanded data resource for the proteins of M. tuberculosis would aid in fine-tuning the existing protein and gene regulatory networks, thereby improving system-wide studies. The primary goal of this study was to uncover and characterize smORFs in M. tuberculosis through bioinformatic analysis, shedding light on their functional roles and genomic organization. Further investigation of these potential smORFs would provide valuable insights into the genome organization and functional diversity of the M. tuberculosis proteome.

Bioinformatics and bioactive peptides from foods: Do they work together?

We live in the Big Data Era which affects many aspects of science, including research on bioactive peptides derived from foods, which during the last few decades have been a focus of interest for scientists. These two issues, i.e., the development of computer technologies and progress in the discovery of novel peptides with health-beneficial properties, are closely interrelated. This Chapter presents the example applications of bioinformatics for studying biopeptides, focusing on main aspects of peptide analysis as the starting point, including: (i) the role of peptide databases; (ii) aspects of bioactivity prediction; (iii) simulation of peptide release from proteins. Bioinformatics can also be used for predicting other features of peptides, including ADMET, QSAR, structure, and taste. To answer the question asked "bioinformatics and bioactive peptides from foods: do they work together?", currently it is almost impossible to find examples of peptide research with no bioinformatics involved. However, theoretical predictions are not equivalent to experimental work and always require critical scrutiny. The aspects of compatibility of in silico and in vitro results are also summarized herein.

In Silico and In Vitro Analyses Reveal Promising Antimicrobial Peptides from Myxobacteria

Antimicrobial resistance (AMR) is a global concern, and as soon as new antibiotics are introduced, resistance to those agents emerges. Therefore, there is an increased appetite for alternative antimicrobial agents to traditional antibiotics. Here, we used in silico methods to investigate potential antimicrobial peptides (AMPs) from predatory myxobacteria. Six hundred seventy-two potential AMP sequences were extracted from eight complete myxobacterial genomes. Most putative AMPs were predicted to be active against Klebsiella pneumoniae with least activity being predicted against Staphylococcus aureus. One hundred seventeen AMPs (defined here as 'potent putative AMPs') were predicted to have very good activity against more than two bacterial pathogens, and these were characterized further in silico. All potent putative AMPs were predicted to have anti-inflammatory and antifungal properties, but none was predicted to be active against viruses. Twenty six (22%) of them were predicted to be hemolytic to human erythrocytes, five were predicted to have anticancer properties, and 56 (47%) were predicted to be biofilm active. In vitro assays using four synthesized AMPs showed high MIC values (e.g. So_ce_56_913 250 µg/ml and Coral_AMP411 125 µg/ml against E. coli). However, antibiofilm assays showed a substantial reduction in numbers (e.g. Coral_AMP411 and Myxo_mac104 showed a 69% and 73% reduction, respectively, at the lowest concentration against E. coli) compared to traditional antibiotics. Fourteen putative AMPs had high sequence similarity to proteins which were functionally associated with proteins of known function. The myxobacterial genomes also possessed a variety of biosynthetic gene clusters (BGCs) that can encode antimicrobial secondary metabolites, but their numbers did not correlate with those of the AMPs. We suggest that AMPs from myxobacteria are a promising source of novel antimicrobial agents with a plethora of biological properties.

Antioxidant and Antimicrobial Peptides Derived from Food Proteins

Recently, the demand for food proteins in the market has increased due to a rise in degenerative illnesses that are associated with the excessive production of free radicals and the unwanted side effects of various drugs, for which researchers have suggested diets rich in bioactive compounds. Some of the functional compounds present in foods are antioxidant and antimicrobial peptides, which are used to produce foods that promote health and to reduce the consumption of antibiotics. These peptides have been obtained from various sources of proteins, such as foods and agri-food by-products, via enzymatic hydrolysis and microbial fermentation. Peptides with antioxidant properties exert effective metal ion (Fe2+/Cu2+) chelating activity and lipid peroxidation inhibition, which may lead to notably beneficial effects in promoting human health and food processing. Antimicrobial peptides are small oligo-peptides generally containing from 10 to 100 amino acids, with a net positive charge and an amphipathic structure; they are the most important components of the antibacterial defense of organisms at almost all levels of life-bacteria, fungi, plants, amphibians, insects, birds and mammals-and have been suggested as natural compounds that neutralize the toxicity of reactive oxygen species generated by antibiotics and the stress generated by various exogenous sources. This review discusses what antioxidant and antimicrobial peptides are, their source, production, some bioinformatics tools used for their obtainment, emerging technologies, and health benefits.

A workflow to identify novel proteins based on the direct mapping of peptide-spectrum-matches to genomic locations

Background

Small Proteins have received increasing attention in recent years. They have in particular been implicated as signals contributing to the coordination of bacterial communities. In genome annotations they are often missing or hidden among large numbers of hypothetical proteins because genome annotation pipelines often exclude short open reading frames or over-predict hypothetical proteins based on simple models. The validation of novel proteins, and in particular of small proteins (sProteins), therefore requires additional evidence. Proteogenomics is considered the gold standard for this purpose. It extends beyond established annotations and includes all possible open reading frames (ORFs) as potential sources of peptides, thus allowing the discovery of novel, unannotated proteins. Typically this results in large numbers of putative novel small proteins fraught with large fractions of false-positive predictions.

Results

We observe that number and quality of the peptide-spectrum matches (PSMs) that map to a candidate ORF can be highly informative for the purpose of distinguishing proteins from spurious ORF annotations. We report here on a workflow that aggregates PSM quality information and local context into simple descriptors and reliably separates likely proteins from the large pool of false-positive, i.e., most likely untranslated ORFs. We investigated the artificial gut microbiome model SIHUMIx, comprising eight different species, for which we validate 5114 proteins that have previously been annotated only as hypothetical ORFs. In addition, we identified 37 non-annotated protein candidates for which we found evidence at the proteomic and transcriptomic level. Half (19) of these candidates have close functional homologs in other species. Another 12 candidates have homologs designated as hypothetical proteins in other species. The remaining six candidates are short (< 100 AA) and are most likely bona fide novel proteins.

Conclusions

The aggregation of PSM quality information for predicted ORFs provides a robust and efficient method to identify novel proteins in proteomics data. The workflow is in particular capable of identifying small proteins and frameshift variants. Since PSMs are explicitly mapped to genomic locations, it furthermore facilitates the integration of transcriptomics data and other sources of genome-level information.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04159-8.

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.

From biomedicinal to in silico models and back to therapeutics: a review on the advancement of peptidic modeling

Bioactive peptides participate in numerous metabolic functions of living organisms and have emerged as potential therapeutics on a diverse range of diseases. Albeit peptide design does not go without challenges, overwhelming advancements on in silico methodologies have increased the scope of peptide-based drug design and discovery to an unprecedented amount. Within an in silico model versus an experimental validation scenario, this review aims to summarize and discuss how different in silico techniques contribute at present to the design of peptide-based molecules. Published in silico results from 2014 to 2018 were selected and discriminated in major methodological groups, allowing a transversal analysis, promoting a landscape vision and asserting its increasing value in drug design.

The pepATTRACT web server for blind, large-scale peptide-protein docking

Peptide–protein interactions are ubiquitous in the cell and form an important part of the interactome. Computational docking methods can complement experimental characterization of these complexes, but current protocols are not applicable on the proteome scale. pepATTRACT is a novel docking protocol that is fully blind, i.e. it does not require any information about the binding site. In various stages of its development, pepATTRACT has participated in CAPRI, making successful predictions for five out of seven protein–peptide targets. Its performance is similar or better than state-of-the-art local docking protocols that do require binding site information. Here we present a novel web server that carries out the rigid-body stage of pepATTRACT. On the peptiDB benchmark, the web server generates a correct model in the top 50 in 34% of the cases. Compared to the full pepATTRACT protocol, this leads to some loss of performance, but the computation time is reduced from ∼18 h to ∼10 min. Combined with the fact that it is fully blind, this makes the web server well-suited for large-scale in silico protein–peptide docking experiments. The rigid-body pepATTRACT server is freely available at http://bioserv.rpbs.univ-paris-diderot.fr/services/pepATTRACT.

StraPep: a structure database of bioactive peptides

Bioactive peptides, with a variety of biological activities and wide distribution in nature, have attracted great research interest in biological and medical fields, especially in pharmaceutical industry. The structural information of bioactive peptide is important for the development of peptide-based drugs. Many databases have been developed cataloguing bioactive peptides. However, to our knowledge, database dedicated to collect all the bioactive peptides with known structure is not available yet. Thus, we developed StraPep, a structure database of bioactive peptides. StraPep holds 3791 bioactive peptide structures, which belong to 1312 unique bioactive peptide sequences. About 905 out of 1312 (68%) bioactive peptides in StraPep contain disulfide bonds, which is significantly higher than that (21%) of PDB. Interestingly, 150 out of 616 (24%) bioactive peptides with three or more disulfide bonds form a structural motif known as cystine knot, which confers considerable structural stability on proteins and is an attractive scaffold for drug design. Detailed information of each peptide, including the experimental structure, the location of disulfide bonds, secondary structure, classification, post-translational modification and so on, has been provided. A wide range of user-friendly tools, such as browsing, sequence and structure-based searching and so on, has been incorporated into StraPep. We hope that this database will be helpful for the research community.

Database URL: http://isyslab.info/StraPep

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

Given the importance of peptide-mediated protein interactions in cellular processes, protein-peptide docking has received increasing attention. Here, we have developed a Hierarchical flexible Peptide Docking approach through fast generation and ensemble docking of peptide conformations, which is referred to as HPepDock. Tested on the LEADS-PEP benchmark data set of 53 diverse complexes with peptides of 3-12 residues, HPepDock performed significantly better than the 11 docking protocols of five small-molecule docking programs (DOCK, AutoDock, AutoDock Vina, Surflex, and GOLD) in predicting near-native binding conformations. HPepDock was also evaluated on the 19 bound/unbound and 10 unbound/unbound protein-peptide complexes of the Glide SP-PEP benchmark and showed an overall better performance than Glide SP-PEP+MM-GBSA and FlexPepDock in both bound and unbound docking. HPepDock is computationally efficient, and the average running time for docking a peptide is ∼15 min with the range from about 1 min for short peptides to around 40 min for long peptides.

Efficient conformational ensemble generation of protein-bound peptides

Conformation generation of protein-bound peptides is critical for the determination of protein–peptide complex structures. Despite significant progress in conformer generation of small molecules, few methods have been developed for modeling protein-bound peptide conformations. Here, we have developed a fast de novo peptide modeling algorithm, referred to as MODPEP, for conformational sampling of protein-bound peptides. Given a sequence, MODPEP builds the peptide 3D structure from scratch by assembling amino acids or helix fragments based on constructed rotamer and helix libraries. The MODPEP algorithm was tested on a diverse set of 910 experimentally determined protein-bound peptides with 3–30 amino acids from the PDB and obtained an average accuracy of 1.90 Å when 200 conformations were sampled for each peptide. On average, MODPEP obtained a success rate of 74.3% for all the 910 peptides and ≥ 90% for short peptides with 3–10 amino acids in reproducing experimental protein-bound structures. Comparative evaluations of MODPEP with three other conformer generation methods, PEP-FOLD3, RDKit, and Balloon, have also been performed in both accuracy and success rate. MODPEP is fast and can generate 100 conformations for less than one second. The fast MODPEP will be beneficial for large-scale de novo modeling and docking of peptides. The MODPEP program and libraries are available for download at http://huanglab.phys.hust.edu.cn/.

Cyclisation mechanisms in the biosynthesis of ribosomally synthesised and post-translationally modified peptides

Ribosomally synthesised and post-translationally modified peptides (RiPPs) are a large class of natural products that are remarkably chemically diverse given an intrinsic requirement to be assembled from proteinogenic amino acids. The vast chemical space occupied by RiPPs means that they possess a wide variety of biological activities, and the class includes antibiotics, co-factors, signalling molecules, anticancer and anti-HIV compounds, and toxins. A considerable amount of RiPP chemical diversity is generated from cyclisation reactions, and the current mechanistic understanding of these reactions will be discussed here. These cyclisations involve a diverse array of chemical reactions, including 1,4-nucleophilic additions, [4 + 2] cycloadditions, ATP-dependent heterocyclisation to form thiazolines or oxazolines, and radical-mediated reactions between unactivated carbons. Future prospects for RiPP pathway discovery and characterisation will also be highlighted.

PEP-FOLD3: faster de novo structure prediction for linear peptides in solution and in complex

Structure determination of linear peptides of 5–50 amino acids in aqueous solution and interacting with proteins is a key aspect in structural biology. PEP-FOLD3 is a novel computational framework, that allows both (i) de novo free or biased prediction for linear peptides between 5 and 50 amino acids, and (ii) the generation of native-like conformations of peptides interacting with a protein when the interaction site is known in advance. PEP-FOLD3 is fast, and usually returns solutions in a few minutes. Testing PEP-FOLD3 on 56 peptides in aqueous solution led to experimental-like conformations for 80% of the targets. Using a benchmark of 61 peptide–protein targets starting from the unbound form of the protein receptor, PEP-FOLD3 was able to generate peptide poses deviating on average by 3.3Å from the experimental conformation and return a native-like pose in the first 10 clusters for 52% of the targets. PEP-FOLD3 is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD3.

The Road to Metagenomics: From Microbiology to DNA Sequencing Technologies and Bioinformatics

The study of microorganisms that pervade each and every part of this planet has encountered many challenges through time such as the discovery of unknown organisms and the understanding of how they interact with their environment. The aim of this review is to take the reader along the timeline and major milestones that led us to modern metagenomics. This new and thriving area is likely to be an important contributor to solve different problems. The transition from classical microbiology to modern metagenomics studies has required the development of new branches of knowledge and specialization. Here, we will review how the availability of high-throughput sequencing technologies has transformed microbiology and bioinformatics and how to tackle the inherent computational challenges that arise from the DNA sequencing revolution. New computational methods are constantly developed to collect, process, and extract useful biological information from a variety of samples and complex datasets, but metagenomics needs the integration of several of these computational methods. Despite the level of specialization needed in bioinformatics, it is important that life-scientists have a good understanding of it for a correct experimental design, which allows them to reveal the information in a metagenome.

Multiple Simulated Annealing-Molecular Dynamics (MSA-MD) for Conformational Space Search of Peptide and Miniprotein

Protein and peptide structure predictions are of paramount importance for understanding their functions, as well as the interactions with other molecules. However, the use of molecular simulation techniques to directly predict the peptide structure from the primary amino acid sequence is always hindered by the rough topology of the conformational space and the limited simulation time scale. We developed here a new strategy, named Multiple Simulated Annealing-Molecular Dynamics (MSA-MD) to identify the native states of a peptide and miniprotein. A cluster of near native structures could be obtained by using the MSA-MD method, which turned out to be significantly more efficient in reaching the native structure compared to continuous MD and conventional SA-MD simulation.