Peptide-based drug discovery through artificial intelligence: towards an autonomous design of therapeutic peptides.

With their diverse biological activities, peptides are promising candidates for therapeutic applications, showing antimicrobial, antitumour and hormonal signalling capabilities. Despite their advantages, therapeutic peptides face challenges such as short half-life, limited oral bioavailability and susceptibility to plasma degradation. The rise of computational tools and artificial intelligence (AI) in peptide research has spurred the development of advanced methodologies and databases that are pivotal in the exploration of these complex macromolecules. This perspective delves into integrating AI in peptide development, encompassing classifier methods, predictive systems and the avant-garde design facilitated by deep-generative models like generative adversarial networks and variational autoencoders. There are still challenges, such as the need for processing optimization and careful validation of predictive models. This work outlines traditional strategies for machine learning model construction and training techniques and proposes a comprehensive AI-assisted peptide design and validation pipeline. The evolving landscape of peptide design using AI is emphasized, showcasing the practicality of these methods in expediting the development and discovery of novel peptides within the context of peptide-based drug discovery.

Protein Language Models and Machine Learning Facilitate the Identification of Antimicrobial Peptides.

Peptides are bioactive molecules whose functional versatility in living organisms has led to successful applications in diverse fields. In recent years, the amount of data describing peptide sequences and function collected in open repositories has substantially increased, allowing the application of more complex computational models to study the relations between the peptide composition and function. This work introduces AMP-Detector, a sequence-based classification model for the detection of peptides' functional biological activity, focusing on accelerating the discovery and de novo design of potential antimicrobial peptides (AMPs). AMP-Detector introduces a novel sequence-based pipeline to train binary classification models, integrating protein language models and machine learning algorithms. This pipeline produced 21 models targeting antimicrobial, antiviral, and antibacterial activity, achieving average precision exceeding 83%. Benchmark analyses revealed that our models outperformed existing methods for AMPs and delivered comparable results for other biological activity types. Utilizing the Peptide Atlas, we applied AMP-Detector to discover over 190,000 potential AMPs and demonstrated that it is an integrative approach with generative learning to aid in de novo design, resulting in over 500 novel AMPs. The combination of our methodology, robust models, and a generative design strategy offers a significant advancement in peptide-based drug discovery and represents a pivotal tool for therapeutic applications.

BEST: a Shiny/R web-based application to easily retrieve cross-related enzyme functional parameters and information from BRENDA.

MOTIVATION: BRENDA is the largest enzyme functional database, containing information of 84 000 experimentally characterized enzyme entries. This database is an invaluable resource for researchers in the biological field, which classifies enzyme-related information in categories that are very useful to obtain specific functional and protein engineering information for enzyme families. However, the BRENDA web interface, the most used by researchers with a non-informatic background, does not allow the user to cross-reference data from different categories or sub-categories in the database. Obtaining information in an easy and fast way, in a friendly web interface, without the necessity to have a deep informatics knowledge, will facilitate and improve research in the enzymology and protein engineering field. RESULTS: We developed the Brenda Easy Search Tool (BEST), an interactive Shiny/R application that enables querying the BRENDA database for complex cross-tabulated characteristics, and retrieving enzyme-related parameters and information readily and efficiently, which can be used for the study of enzyme function or as an input for other bioinformatics tools. AVAILABILITY AND IMPLEMENTATION: BEST and its tutorial are freely available from https://pesb2.cl/best/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Total mutational load and clinical features as predictors of the metastatic status in lung adenocarcinoma and squamous cell carcinoma patients.

BACKGROUND: Recently, extensive cancer genomic studies have revealed mutational and clinical data of large cohorts of cancer patients. For example, the Pan-Lung Cancer 2016 dataset (part of The Cancer Genome Atlas project), summarises the mutational and clinical profiles of different subtypes of Lung Cancer (LC). Mutational and clinical signatures have been used independently for tumour typification and prediction of metastasis in LC patients. Is it then possible to achieve better typifications and predictions when combining both data streams? METHODS: In a cohort of 1144 Lung Adenocarcinoma (LUAD) and Lung Squamous Cell Carcinoma (LSCC) patients, we studied the number of missense mutations (hereafter, the Total Mutational Load TML) and distribution of clinical variables, for different classes of patients. Using the TML and different sets of clinical variables (tumour stage, age, sex, smoking status, and packs of cigarettes smoked per year), we built Random Forest classification models that calculate the likelihood of developing metastasis. RESULTS: We found that LC patients different in age, smoking status, and tumour type had significantly different mean TMLs. Although TML was an informative feature, its effect was secondary to the "tumour stage" feature. However, its contribution to the classification is not redundant with the latter; models trained using both TML and tumour stage performed better than models trained using only one of these variables. We found that models trained in the entire dataset (i.e., without using dimensionality reduction techniques) and without resampling achieved the highest performance, with an F1 score of 0.64 (95%CrI [0.62, 0.66]). CONCLUSIONS: Clinical variables and TML should be considered together when assessing the likelihood of LC patients progressing to metastatic states, as the information these encode is not redundant. Altogether, we provide new evidence of the need for comprehensive diagnostic tools for metastasis.

Relaxing restrictions at the pace of vaccination increases freedom and guards against further COVID-19 waves.

Mass vaccination offers a promising exit strategy for the COVID-19 pandemic. However, as vaccination progresses, demands to lift restrictions increase, despite most of the population remaining susceptible. Using our age-stratified SEIRD-ICU compartmental model and curated epidemiological and vaccination data, we quantified the rate (relative to vaccination progress) at which countries can lift non-pharmaceutical interventions without overwhelming their healthcare systems. We analyzed scenarios ranging from immediately lifting restrictions (accepting high mortality and morbidity) to reducing case numbers to a level where test-trace-and-isolate (TTI) programs efficiently compensate for local spreading events. In general, the age-dependent vaccination roll-out implies a transient decrease of more than ten years in the average age of ICU patients and deceased. The pace of vaccination determines the speed of lifting restrictions; Taking the European Union (EU) as an example case, all considered scenarios allow for steadily increasing contacts starting in May 2021 and relaxing most restrictions by autumn 2021. Throughout summer 2021, only mild contact restrictions will remain necessary. However, only high vaccine uptake can prevent further severe waves. Across EU countries, seroprevalence impacts the long-term success of vaccination campaigns more strongly than age demographics. In addition, we highlight the need for preventive measures to reduce contagion in school settings throughout the year 2021, where children might be drivers of contagion because of them remaining susceptible. Strategies that maintain low case numbers, instead of high ones, reduce infections and deaths by factors of eleven and five, respectively. In general, policies with low case numbers significantly benefit from vaccination, as the overall reduction in susceptibility will further diminish viral spread. Keeping case numbers low is the safest long-term strategy because it considerably reduces mortality and morbidity and offers better preparedness against emerging escape or more contagious virus variants while still allowing for higher contact numbers (freedom) with progressing vaccinations.

On the heterogeneous spread of COVID-19 in Chile.

Non-pharmaceutical interventions (NPIs) have played a crucial role in controlling the spread of COVID-19. Nevertheless, NPI efficacy varies enormously between and within countries, mainly because of population and behavioral heterogeneity. In this work, we adapted a multi-group SEIRA model to study the spreading dynamics of COVID-19 in Chile, representing geographically separated regions of the country by different groups. We use national mobilization statistics to estimate the connectivity between regions and data from governmental repositories to obtain COVID-19 spreading and death rates in each region. We then assessed the effectiveness of different NPIs by studying the temporal evolution of the reproduction number R t . Analysing data-driven and model-based estimates of R t , we found a strong coupling of different regions, highlighting the necessity of organized and coordinated actions to control the spread of SARS-CoV-2. Finally, we evaluated different scenarios to forecast the evolution of COVID-19 in the most densely populated regions, finding that the early lifting of restriction probably will lead to novel outbreaks.

A multi-group SEIRA model for the spread of COVID-19 among heterogeneous populations.

The outbreak and propagation of COVID-19 have posed a considerable challenge to modern society. In particular, the different restrictive actions taken by governments to prevent the spread of the virus have changed the way humans interact and conceive interaction. Due to geographical, behavioral, or economic factors, different sub-groups among a population are more (or less) likely to interact, and thus to spread/acquire the virus. In this work, we present a general multi-group SEIRA model for representing the spread of COVID-19 among a heterogeneous population and test it in a numerical case of study. By highlighting its applicability and the ease with which its general formulation can be adapted to particular studies, we expect our model to lead us to a better understanding of the evolution of this pandemic and to better public-health policies to control it.

Statistically-based methodology for revealing real contagion trends and correcting delay-induced errors in the assessment of COVID-19 pandemic.

COVID-19 pandemic has reshaped our world in a timescale much shorter than what we can understand. Particularities of SARS-CoV-2, such as its persistence in surfaces and the lack of a curative treatment or vaccine against COVID-19, have pushed authorities to apply restrictive policies to control its spreading. As data drove most of the decisions made in this global contingency, their quality is a critical variable for decision-making actors, and therefore should be carefully curated. In this work, we analyze the sources of error in typically reported epidemiological variables and usual tests used for diagnosis, and their impact on our understanding of COVID-19 spreading dynamics. We address the existence of different delays in the report of new cases, induced by the incubation time of the virus and testing-diagnosis time gaps, and other error sources related to the sensitivity/specificity of the tests used to diagnose COVID-19. Using a statistically-based algorithm, we perform a temporal reclassification of cases to avoid delay-induced errors, building up new epidemiologic curves centered in the day where the contagion effectively occurred. We also statistically enhance the robustness behind the discharge/recovery clinical criteria in the absence of a direct test, which is typically the case of non-first world countries, where the limited testing capabilities are fully dedicated to the evaluation of new cases. Finally, we applied our methodology to assess the evolution of the pandemic in Chile through the Effective Reproduction Number Rt , identifying different moments in which data was misleading governmental actions. In doing so, we aim to raise public awareness of the need for proper data reporting and processing protocols for epidemiological modelling and predictions.

An integrated computational approach can classify VHL missense mutations according to risk of clear cell renal carcinoma.

Mutations in the von Hippel-Lindau (VHL) gene are pathogenic in VHL disease, congenital polycythaemia and clear cell renal carcinoma (ccRCC). pVHL forms a ternary complex with elongin C and elongin B, critical for pVHL stability and function, which interacts with Cullin-2 and RING-box protein 1 to target hypoxia-inducible factor for polyubiquitination and proteasomal degradation. We describe a comprehensive database of missense VHL mutations linked to experimental and clinical data. We use predictions from in silico tools to link the functional effects of missense VHL mutations to phenotype. The risk of ccRCC in VHL disease is linked to the degree of destabilization resulting from missense mutations. An optimized binary classification system (symphony), which integrates predictions from five in silico methods, can predict the risk of ccRCC associated with VHL missense mutations with high sensitivity and specificity. We use symphony to generate predictions for risk of ccRCC for all possible VHL missense mutations and present these predictions, in association with clinical and experimental data, in a publically available, searchable web server.

Iron-induced reactive oxygen species mediate transporter DMT1 endocytosis and iron uptake in intestinal epithelial cells.

Recent evidence shows that iron induces the endocytosis of the iron transporter dimetal transporter 1 (DMT1) during intestinal absorption. We, and others, have proposed that iron-induced DMT1 internalization underlies the mucosal block phenomena, a regulatory response that downregulates intestinal iron uptake after a large oral dose of iron. In this work, we investigated the participation of reactive oxygen species (ROS) in the establishment of this response. By means of selective surface protein biotinylation of polarized Caco-2 cells, we determined the kinetics of DMT1 internalization from the apical membrane after an iron challenge. The initial decrease in DMT1 levels in the apical membrane induced by iron was followed at later times by increased levels of DMT1. Addition of Fe(2+), but not of Cd(2+), Zn(2+), Cu(2+), or Cu(1+), induced the production of intracellular ROS, as detected by 2',7'-dichlorofluorescein (DCF) fluorescence. Preincubation with the antioxidant N-acetyl-l-cysteine (NAC) resulted in increased DMT1 at the apical membrane before and after addition of iron. Similarly, preincubation with the hydroxyl radical scavenger dimethyl sulfoxide (DMSO) resulted in the enhanced presence of DMT1 at the apical membrane. The decrease of DMT1 levels at the apical membrane induced by iron was associated with decreased iron uptake rates. A kinetic mathematical model based on operational rate constants of DMT1 endocytosis and exocytosis is proposed. The model qualitatively captures the experimental observations and accurately describes the effect of iron, NAC, and DMSO on the apical distribution of DMT1. Taken together, our data suggest that iron uptake induces the production of ROS, which modify DMT1 endocytic cycling, thus changing the iron transport activity at the apical membrane.

Early mutational signatures and transmissibility of SARS-CoV-2 Gamma and Lambda variants in Chile.

Genomic surveillance (GS) programmes were crucial in identifying and quantifying the mutating patterns of SARS-CoV-2 during the COVID-19 pandemic. In this work, we develop a Bayesian framework to quantify the relative transmissibility of different variants tailored for regions with limited GS. We use it to study the relative transmissibility of SARS-CoV-2 variants in Chile. Among the 3443 SARS-CoV-2 genomes collected between January and June 2021, where sampling was designed to be representative, the Gamma (P.1), Lambda (C.37), Alpha (B.1.1.7), B.1.1.348, and B.1.1 lineages were predominant. We found that Lambda and Gamma variants' reproduction numbers were 5% (95% CI: [1%, 14%]) and 16% (95% CI: [11%, 21%]) larger than Alpha's, respectively. Besides, we observed a systematic mutation enrichment in the Spike gene for all circulating variants, which strongly correlated with variants' transmissibility during the studied period (r = 0.93, p-value = 0.025). We also characterised the mutational signatures of local samples and their evolution over time and with the progress of vaccination, comparing them with those of samples collected in other regions worldwide. Altogether, our work provides a reliable method for quantifying variant transmissibility under subsampling and emphasises the importance of continuous genomic surveillance.

Thyroid dyshormonogenesis is mainly caused by TPO mutations in consanguineous community.

OBJECTIVE: In this study, we aimed to investigate the genetic background of thyroid dyshormonogenesis (TDH). CONTEXT: Thyroid dyshormonogenesis comprises 10-15% of all cases of congenital hypothyroidism (CH), which is the most common neonatal endocrine disorder, and might result from disruptions at any stage of thyroid hormone biosynthesis. Currently seven genes (NIS, TPO, PDS, TG, IYD, DUOX2 and DUOXA2) have been implicated in the aetiology of the disease. DESIGN: As TDH is mostly inherited in an autosomal recessive manner, we planned to conduct the study in consanguineous/multi-case families. PATIENTS: One hundred and four patients with congenital TDH all coming from consanguineous and/or multi-case families. MEASUREMENTS: Initially, we performed potential linkage analysis of cases to all seven causative-TDH loci as well as direct sequencing of the TPO gene in cases we could not exclude linkage to this locus. In addition, in silico analyses of novel missense mutations were carried out. RESULTS: TPO had the highest potential for linkage and we identified 21 TPO mutations in 28 TDH cases showing potential linkage to this locus. Four of 10 distinct TPO mutations detected in this study were novel (A5T, Y55X, E596X, D633N). CONCLUSIONS: This study underlines the importance of molecular genetic studies in diagnosis, classification and prognosis of CH and proposes a comprehensive mutation screening by new sequencing technology in all newly diagnosed primary CH cases.

Purification of transthyretin as nutritional biomarker of selenium status.

Transthyretin has been proposed as nutritional biomarker of selenium intake. Previous transthyretin purification methods used different procedures to isolate transthyretin either from plasma or from pathological urine of humans. In general, the procedure for purification of transthyretin is laborious and expensive, and extensive sample recycling is necessary for purification in appreciable amounts. This work proposes a new, promissory, and cheap two-step process to purify transthyretin from blood plasma, composed by a first aqueous two-phase system fractionation followed by affinity chromatography, using thyroxine-immobilized on epoxy-activated Sepharose CL-6B. The aqueous two-phase system fractionation was demonstrated to perform better than commercial immunoaffinity-based kits for albumin depletion in blood plasma samples and is an effective first step for transthyretin purification. Thyroxine affinity chromatography was designed to bind transthyretin with high affinity, and was demonstrated to be useful to purify transthyretin, but was unable to completely resolve transthyretin from thyroxine-binding globulin and serum albumin, although the relative amount of albumin was lowered in the eluates. This purification process could be used in nutritional diagnosis tools or as a first step in structural and functional studies.

Patient-Wise Methodology to Assess Glycemic Health Status: Applications to Quantify the Efficacy and Physiological Targets of Polyphenols on Glycemic Control.

A growing body of evidence indicates that dietary polyphenols could be used as an early intervention to treat glucose-insulin (G-I) dysregulation. However, studies report heterogeneous information, and the targets of the intervention remain largely elusive. In this work, we provide a general methodology to quantify the effects of any given polyphenol-rich food or formulae over glycemic regulation in a patient-wise manner using an Oral Glucose Tolerance Test (OGTT). We use a mathematical model to represent individual OGTT curves as the coordinated action of subsystems, each one described by a parameter with physiological interpretation. Using the parameter values calculated for a cohort of 1198 individuals, we propose a statistical model to calculate the risk of dysglycemia and the coordination among subsystems for each subject, thus providing a continuous and individual health assessment. This method allows identifying individuals at high risk of dysglycemia-which would have been missed with traditional binary diagnostic methods-enabling early nutritional intervention with a polyphenol-supplemented diet where it is most effective and desirable. Besides, the proposed methodology assesses the effectiveness of interventions over time when applied to the OGTT curves of a treated individual. We illustrate the use of this method in a case study to assess the dose-dependent effects of Delphinol® on reducing dysglycemia risk and improving the coordination between subsystems. Finally, this strategy enables, on the one hand, the use of low-cost, non-invasive methods in population-scale nutritional studies. On the other hand, it will help practitioners assess the effectiveness of an intervention based on individual vulnerabilities and adapt the treatment to manage dysglycemia and avoid its progression into disease.

Generalized Property-Based Encoders and Digital Signal Processing Facilitate Predictive Tasks in Protein Engineering.

Computational methods in protein engineering often require encoding amino acid sequences, i.e., converting them into numeric arrays. Physicochemical properties are a typical choice to define encoders, where we replace each amino acid by its value for a given property. However, what property (or group thereof) is best for a given predictive task remains an open problem. In this work, we generalize property-based encoding strategies to maximize the performance of predictive models in protein engineering. First, combining text mining and unsupervised learning, we partitioned the AAIndex database into eight semantically-consistent groups of properties. We then applied a non-linear PCA within each group to define a single encoder to represent it. Then, in several case studies, we assess the performance of predictive models for protein and peptide function, folding, and biological activity, trained using the proposed encoders and classical methods (One Hot Encoder and TAPE embeddings). Models trained on datasets encoded with our encoders and converted to signals through the Fast Fourier Transform (FFT) increased their precision and reduced their overfitting substantially, outperforming classical approaches in most cases. Finally, we propose a preliminary methodology to create de novo sequences with desired properties. All these results offer simple ways to increase the performance of general and complex predictive tasks in protein engineering without increasing their complexity.

Mutagenesis Objective Search and Selection Tool (MOSST): an algorithm to predict structure-function related mutations in proteins.

BACKGROUND: Functionally relevant artificial or natural mutations are difficult to assess or predict if no structure-function information is available for a protein. This is especially important to correctly identify functionally significant non-synonymous single nucleotide polymorphisms (nsSNPs) or to design a site-directed mutagenesis strategy for a target protein. A new and powerful methodology is proposed to guide these two decision strategies, based only on conservation rules of physicochemical properties of amino acids extracted from a multiple alignment of a protein family where the target protein belongs, with no need of explicit structure-function relationships. RESULTS: A statistical analysis is performed over each amino acid position in the multiple protein alignment, based on different amino acid physical or chemical characteristics, including hydrophobicity, side-chain volume, charge and protein conformational parameters. The variances of each of these properties at each position are combined to obtain a global statistical indicator of the conservation degree of each property. Different types of physicochemical conservation are defined to characterize relevant and irrelevant positions. The differences between statistical variances are taken together as the basis of hypothesis tests at each position to search for functionally significant mutable sites and to identify specific mutagenesis targets. The outcome is used to statistically predict physicochemical consensus sequences based on different properties and to calculate the amino acid propensities at each position in a given protein. Hence, amino acid positions are identified that are putatively responsible for function, specificity, stability or binding interactions in a family of proteins. Once these key functional positions are identified, position-specific statistical distributions are applied to divide the 20 common protein amino acids in each position of the protein's primary sequence into a group of functionally non-disruptive amino acids and a second group of functionally deleterious amino acids. CONCLUSIONS: With this approach, not only conserved amino acid positions in a protein family can be labeled as functionally relevant, but also non-conserved amino acid positions can be identified to have a physicochemically meaningful functional effect. These results become a discriminative tool in the selection and elaboration of rational mutagenesis strategies for the protein. They can also be used to predict if a given nsSNP, identified, for instance, in a genomic-scale analysis, can have a functional implication for a particular protein and which nsSNPs are most likely to be functionally silent for a protein. This analytical tool could be used to rapidly and automatically discard any irrelevant nsSNP and guide the research focus toward functionally significant mutations. Based on preliminary results and applications, this technique shows promising performance as a valuable bioinformatics tool to aid in the development of new protein variants and in the understanding of function-structure relationships in proteins.

Enhanced Cellular Uptake of H-Chain Human Ferritin Containing Gold Nanoparticles.

Gold nanoparticles (AuNP) capped with biocompatible layers have functional optical, chemical, and biological properties as theranostic agents in biomedicine. The ferritin protein containing in situ synthesized AuNPs has been successfully used as an effective and completely biocompatible nanocarrier for AuNPs in human cell lines and animal experiments in vivo. Ferritin can be uptaken by different cell types through receptor-mediated endocytosis. Despite these advantages, few efforts have been made to evaluate the toxicity and cellular internalization of AuNP-containing ferritin nanocages. In this work, we study the potential of human heavy-chain (H) and light-chain (L) ferritin homopolymers as nanoreactors to synthesize AuNPs and their cytotoxicity and cellular uptake in different cell lines. The results show very low toxicity of ferritin-encapsulated AuNPs on different human cell lines and demonstrate that efficient cellular ferritin uptake depends on the specific H or L protein chains forming the ferritin protein cage and the presence or absence of metallic cargo. Cargo-devoid apoferritin is poorly internalized in all cell lines, and the highest ferritin uptake was achieved with AuNP-loaded H-ferritin homopolymers in transferrin-receptor-rich cell lines, showing more than seven times more uptake than apoferritin.

Novel Symbiotic Genome-Scale Model Reveals Wolbachia's Arboviral Pathogen Blocking Mechanism in Aedes aegypti.

Wolbachia are endosymbiont bacteria known to infect arthropods causing different effects, such as cytoplasmic incompatibility and pathogen blocking in Aedes aegypti. Although several Wolbachia strains have been studied, there is little knowledge regarding the relationship between this bacterium and their hosts, particularly on their obligate endosymbiont nature and its pathogen blocking ability. Motivated by the potential applications on disease control, we developed a genome-scale model of two Wolbachia strains: wMel and the strongest Dengue blocking strain known to date: wMelPop. The obtained metabolic reconstructions exhibit an energy metabolism relying mainly on amino acids and lipid transport to support cell growth that is consistent with altered lipid and cholesterol metabolism in Wolbachia-infected mosquitoes. The obtained metabolic reconstruction was then coupled with a reconstructed mosquito model to retrieve a symbiotic genome-scale model accounting for 1,636 genes and 6,408 reactions of the Aedes aegypti-Wolbachia interaction system. Simulation of an arboviral infection in the obtained novel symbiotic model represents a metabolic scenario characterized by pathogen blocking in higher titer Wolbachia strains, showing that pathogen blocking by Wolbachia infection is consistent with competition for lipid and amino acid resources between arbovirus and this endosymbiotic bacteria. IMPORTANCE Arboviral diseases such as Zika and Dengue have been on the rise mainly due to climate change, and the development of new treatments and strategies to limit their spreading is needed. The use of Wolbachia as an approach for disease control has motivated new research related to the characterization of the mechanisms that underlie its pathogen-blocking properties. In this work, we propose a new approach for studying the metabolic interactions between Aedes aegypti and Wolbachia using genome-scale models, finding that pathogen blocking is mainly influenced by competition for the resources required for Wolbachia and viral replication.

Peptipedia: a user-friendly web application and a comprehensive database for peptide research supported by Machine Learning approach.

Peptides have attracted attention during the last decades due to their extraordinary therapeutic properties. Different computational tools have been developed to take advantage of existing information, compiling knowledge and making available the information for common users. Nevertheless, most related tools available are not user-friendly, present redundant information, do not clearly display the data, and usually are specific for particular biological activities, not existing so far, an integrated database with consolidated information to help research peptide sequences. To solve these necessities, we developed Peptipedia, a user-friendly web application and comprehensive database to search, characterize and analyse peptide sequences. Our tool integrates the information from 30 previously reported databases with a total of 92 055 amino acid sequences, making it the biggest repository of peptides with recorded activities to date. Furthermore, we make available a variety of bioinformatics services and statistical modules to increase our tool's usability. Moreover, we incorporated a robust assembled binary classification system to predict putative biological activities for peptide sequences. Our tools' significant differences with other existing alternatives become a substantial contribution for developing biotechnological and bioengineering applications for peptides. Peptipedia is available for non-commercial use as an open-access software, licensed under the GNU General Public License, version GPL 3.0. The web platform is publicly available at peptipedia.cl. Database URL: Both the source code and sample data sets are available in the GitHub repository https://github.com/ProteinEngineering-PESB2/peptipedia.

Mutation in a SARS-CoV-2 Haplotype from Sub-Antarctic Chile Reveals New Insights into the Spike's Dynamics.

The emergence of SARS-CoV-2 variants, as observed with the D614G spike protein mutant and, more recently, with B.1.1.7 (501Y.V1), B.1.351 (501Y.V2) and B.1.1.28.1 (P.1) lineages, represent a continuous threat and might lead to strains of higher infectivity and/or virulence. We report on the occurrence of a SARS-CoV-2 haplotype with nine mutations including D614G/T307I double-mutation of the spike. This variant expanded and completely replaced previous lineages within a short period in the subantarctic Magallanes Region, southern Chile. The rapid lineage shift was accompanied by a significant increase of cases, resulting in one of the highest incidence rates worldwide. Comparative coarse-grained molecular dynamic simulations indicated that T307I and D614G belong to a previously unrecognized dynamic domain, interfering with the mobility of the receptor binding domain of the spike. The T307I mutation showed a synergistic effect with the D614G. Continuous surveillance of new mutations and molecular analyses of such variations are important tools to understand the molecular mechanisms defining infectivity and virulence of current and future SARS-CoV-2 strains.