General Prediction of Peptide-MHC Binding Modes Using Incremental Docking: A Proof of Concept

APE-Gen2.0: Expanding Rapid Class I Peptide-Major Histocompatibility Complex Modeling to Post-Translational Modifications and Noncanonical Peptide Geometries

The recognition of peptides bound to class I major histocompatibility complex (MHC-I) receptors by T-cell receptors (TCRs) is a determinant of triggering the adaptive immune response. While the exact molecular features that drive the TCR recognition are still unknown, studies have suggested that the geometry of the joint peptide-MHC (pMHC) structure plays an important role. As such, there is a definite need for methods and tools that accurately predict the structure of the peptide bound to the MHC-I receptor. In the past few years, many pMHC structural modeling tools have emerged that provide high-quality modeled structures in the general case. However, there are numerous instances of non-canonical cases in the immunopeptidome that the majority of pMHC modeling tools do not attend to, most notably, peptides that exhibit non-standard amino acids and post-translational modifications (PTMs) or peptides that assume non-canonical geometries in the MHC binding cleft. Such chemical and structural properties have been shown to be present in neoantigens; therefore, accurate structural modeling of these instances can be vital for cancer immunotherapy. To this end, we have developed APE-Gen2.0, a tool that improves upon its predecessor and other pMHC modeling tools, both in terms of modeling accuracy and the available modeling range of non-canonical peptide cases. Some of the improvements include (i) the ability to model peptides that have different types of PTMs such as phosphorylation, nitration, and citrullination; (ii) a new and improved anchor identification routine in order to identify and model peptides that exhibit a non-canonical anchor conformation; and (iii) a web server that provides a platform for easy and accessible pMHC modeling. We further show that structures predicted by APE-Gen2.0 can be used to assess the effects that PTMs have in binding affinity in a more accurate manner than just using solely the sequence of the peptide. APE-Gen2.0 is freely available at https://apegen.kavrakilab.org.

Riding the wave of innovation: immunoinformatics in fish disease control

The spread of infectious illnesses has been a significant factor restricting aquaculture production. To maximise aquatic animal health, vaccination tactics are very successful and cost-efficient for protecting fish and aquaculture animals against many disease pathogens. However, due to the increasing number of immunological cases and their complexity, it is impossible to manage, analyse, visualise, and interpret such data without the assistance of advanced computational techniques. Hence, the use of immunoinformatics tools is crucial, as they not only facilitate the management of massive amounts of data but also greatly contribute to the creation of fresh hypotheses regarding immune responses. In recent years, advances in biotechnology and immunoinformatics have opened up new research avenues for generating novel vaccines and enhancing existing vaccinations against outbreaks of infectious illnesses, thereby reducing aquaculture losses. This review focuses on understanding in silico epitope-based vaccine design, the creation of multi-epitope vaccines, the molecular interaction of immunogenic vaccines, and the application of immunoinformatics in fish disease based on the frequency of their application and reliable results. It is believed that it can bridge the gap between experimental and computational approaches and reduce the need for experimental research, so that only wet laboratory testing integrated with in silico techniques may yield highly promising results and be useful for the development of vaccines for fish.

Computational Approaches Drive Developments in Immune-Oncology Therapies for PD-1/PD-L1 Immune Checkpoint Inhibitors

Computational approaches in immune-oncology therapies focus on using data-driven methods to identify potential immune targets and develop novel drug candidates. In particular, the search for PD-1/PD-L1 immune checkpoint inhibitors (ICIs) has enlivened the field, leveraging the use of cheminformatics and bioinformatics tools to analyze large datasets of molecules, gene expression and protein-protein interactions. Up to now, there is still an unmet clinical need for improved ICIs and reliable predictive biomarkers. In this review, we highlight the computational methodologies applied to discovering and developing PD-1/PD-L1 ICIs for improved cancer immunotherapies with a greater focus in the last five years. The use of computer-aided drug design structure- and ligand-based virtual screening processes, molecular docking, homology modeling and molecular dynamics simulations methodologies essential for successful drug discovery campaigns focusing on antibodies, peptides or small-molecule ICIs are addressed. A list of recent databases and web tools used in the context of cancer and immunotherapy has been compilated and made available, namely regarding a general scope, cancer and immunology. In summary, computational approaches have become valuable tools for discovering and developing ICIs. Despite significant progress, there is still a need for improved ICIs and biomarkers, and recent databases and web tools have been compiled to aid in this pursuit.

Peptide-binding specificity prediction using fine-tuned protein structure prediction networks

Peptide-binding proteins play key roles in biology, and predicting their binding specificity is a long-standing challenge. While considerable protein structural information is available, the most successful current methods use sequence information alone, in part because it has been a challenge to model the subtle structural changes accompanying sequence substitutions. Protein structure prediction networks such as AlphaFold model sequence-structure relationships very accurately, and we reasoned that if it were possible to specifically train such networks on binding data, more generalizable models could be created. We show that placing a classifier on top of the AlphaFold network and fine-tuning the combined network parameters for both classification and structure prediction accuracy leads to a model with strong generalizable performance on a wide range of Class I and Class II peptide-MHC interactions that approaches the overall performance of the state-of-the-art NetMHCpan sequence-based method. The peptide-MHC optimized model shows excellent performance in distinguishing binding and non-binding peptides to SH3 and PDZ domains. This ability to generalize well beyond the training set far exceeds that of sequence-only models and should be particularly powerful for systems where less experimental data are available.

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

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

Binding Affinity Calculations of Gluten Peptides to HLA Risk Modifiers: DQ2.5 versus DQ7.5

Free energy perturbation (FEP) calculations can predict relative binding affinities of an antigen and its point mutants to the same human leukocyte antigen (HLA) with high accuracy (e.g., within 1.0 kcal/mol to experiment); however, a more challenging task is to compare binding affinities of wholly different antigens binding to completely different HLAs using FEP. Researchers have used a variety of different FEP schemes to compute and compare absolute binding affinities, with varied success. Here, we propose and assess a unifying scheme to compute the relative binding affinities of different antigens binding to completely different HLAs using absolute binding affinity FEP calculations. We apply our affinity calculation technique to HLA-antigen-T-cell receptor (TCR) systems relevant to celiac disease (CeD) by investigating binding affinity differences between HLA-DQ2.5 (enhanced CeD risk) and HLA-DQ7.5 (CeD protective) in the binary (HLA-gliadin) and ternary (HLA-gliadin-TCR) binding complexes for three gliadin derived epitopes: glia-α1, glia-α2, and glia-ω1. Based on FEP calculations with our carefully designed thermodynamic cycles, we demonstrate that HLA-DQ2.5 has higher binding affinity than HLA-DQ7.5 for gliadin and enhanced binding affinity with a common TCR, agreeing with known results that the HLA-DQ2.5 serotype exhibits increased risk for CeD. Our findings reveal that our proposed absolute binding affinity FEP method is appropriate for predicting HLA binding for disparate antigens with different genotypes. We also discuss atomic-level details of HLA genotypes interacting with gluten peptides and TCRs in regard to the pathogenesis of CeD.

In silico epitope-based vaccine design against influenza a neuraminidase protein: Computational analysis established on B- and T-cell epitope predictions

Objective

Influenza A virus belongs to the most studied virus and its mutant initiates epidemic and pandemics outbreaks. Inoculation is the significant foundation to diminish the risk of infection. To prevent an incidence of influenza from the transmission, various practical approaches require more advancement and progress. More efforts and research must take in front to enhance vaccine efficacy.

Methods

The present research emphasizes the development and expansion of a universal vaccine for the influenza virus. Research focuses on vaccine design with high efficacy. In this study, numerous computational approaches were used, covering a wide range of elements and ideas in bioinformatics methodology. Various B and T-cell epitopic peptides derived from the Neuraminidase protein N1 are recognized by these approaches. With the implementation of numerous obtained databases and bioinformatics tools, the different immune framework methods of the conserved sequences of N1 neuraminidase were analyzed. NCBI databases were employed to retrieve amino acid sequences. The antigenic nature of the neuraminidase sequence was achieved by the VaxiJen server and Kolaskar and Tongaonkar method. After screening of various B and T cell epitopes, one efficient peptide each from B cell epitope and T cell epitopes was assessed for their antigenic determinant vaccine efficacy. Identical two B cell epitopes were recognized from the N1 protein when analyzed using B-cell epitope prediction servers. The detailed examination of amino acid sequences for interpretation of B and T cell epitopes was achieved with the help of the ABCPred and Immune Epitope Database.

Results

Computational immunology via immunoinformatic study exhibited RPNDKTG as having its high conservancy efficiency and demonstrated as a good antigenic, accessible surface hydrophilic B-cell epitope. Among T cell epitope analysis, YVNISNTNF was selected for being a conserved epitope. T cell epitope was also analyzed for its allergenicity and cytotoxicity evaluation. YVNISNTNF epitope was found to be a non-allergen and not toxic for cells as well. This T-cell epitope with maximum world populace coverages was scrutinized for its association with the HLA-DRB1*0401 molecule. Results from docking simulation analyses showed YVNISNTNF having lower binding energy, the radius of gyration (Rg), RMSD values, and RMSE values which make the protein structure more stable and increase its ability to become an epitopic peptide for influenza virus vaccination.

Conclusions

We propose that this epitope analysis may be successfully used as a measurement tool for the robustness of an antigen-antibody reaction between mutant strains in the annual design of the influenza vaccine.

Physicochemical Heuristics for Identifying High Fidelity, Near-Native Structural Models of Peptide/MHC Complexes

There is long-standing interest in accurately modeling the structural features of peptides bound and presented by class I MHC proteins. This interest has grown with the advent of rapid genome sequencing and the prospect of personalized, peptide-based cancer vaccines, as well as the development of molecular and cellular therapeutics based on T cell receptor recognition of peptide-MHC. However, while the speed and accessibility of peptide-MHC modeling has improved substantially over the years, improvements in accuracy have been modest. Accuracy is crucial in peptide-MHC modeling, as T cell receptors are highly sensitive to peptide conformation and capturing fine details is therefore necessary for useful models. Studying nonameric peptides presented by the common class I MHC protein HLA-A*02:01, here we addressed a key question common to modern modeling efforts: from a set of models (or decoys) generated through conformational sampling, which is best? We found that the common strategy of decoy selection by lowest energy can lead to substantial errors in predicted structures. We therefore adopted a data-driven approach and trained functions capable of predicting near native decoys with exceptionally high accuracy. Although our implementation is limited to nonamer/HLA-A*02:01 complexes, our results serve as an important proof of concept from which improvements can be made and, given the significance of HLA-A*02:01 and its preference for nonameric peptides, should have immediate utility in select immunotherapeutic and other efforts for which structural information would be advantageous.

Structural Prediction of Peptide-MHC Binding Modes

The immune system is constantly protecting its host from the invasion of pathogens and the development of cancer cells. The specific CD8⁺ T-cell immune response against virus-infected cells and tumor cells is based on the T-cell receptor recognition of antigenic peptides bound to class I major histocompatibility complexes (MHC) at the surface of antigen presenting cells. Consequently, the peptide binding specificities of the highly polymorphic MHC have important implications for the design of vaccines, for the treatment of autoimmune diseases, and for personalized cancer immunotherapy. Evidence-based machine-learning approaches have been successfully used for the prediction of peptide binders and are currently being developed for the prediction of peptide immunogenicity. However, understanding and modeling the structural details of peptide/MHC binding is crucial for a better understanding of the molecular mechanisms triggering the immunological processes, estimating peptide/MHC affinity using universal physics-based approaches, and driving the design of novel peptide ligands. Unfortunately, due to the large diversity of MHC allotypes and possible peptides, the growing number of 3D structures of peptide/MHC (pMHC) complexes in the Protein Data Bank only covers a small fraction of the possibilities. Consequently, there is a growing need for rapid and efficient approaches to predict 3D structures of pMHC complexes. Here, we review the key characteristics of the 3D structure of pMHC complexes before listing databases and other sources of information on pMHC structures and MHC specificities. Finally, we discuss some of the most prominent pMHC docking software.

DINC-COVID: A webserver for ensemble docking with flexible SARS-CoV-2 proteins

An unprecedented research effort has been undertaken in response to the ongoing COVID-19 pandemic. This has included the determination of hundreds of crystallographic structures of SARS-CoV-2 proteins, and numerous virtual screening projects searching large compound libraries for potential drug inhibitors. Unfortunately, these initiatives have had very limited success in producing effective inhibitors against SARS-CoV-2 proteins. A reason might be an often overlooked factor in these computational efforts: receptor flexibility. To address this issue we have implemented a computational tool for ensemble docking with SARS-CoV-2 proteins. We have extracted representative ensembles of protein conformations from the Protein Data Bank and from in silico molecular dynamics simulations. Twelve pre-computed ensembles of SARS-CoV-2 protein conformations have now been made available for ensemble docking via a user-friendly webserver called DINC-COVID (dinc-covid.kavrakilab.org). We have validated DINC-COVID using data on tested inhibitors of two SARS-CoV-2 proteins, obtaining good correlations between docking-derived binding energies and experimentally-determined binding affinities. Some of the best results have been obtained on a dataset of large ligands resolved via room temperature crystallography, and therefore capturing alternative receptor conformations. In addition, we have shown that the ensembles available in DINC-COVID capture different ranges of receptor flexibility, and that this diversity is useful in finding alternative binding modes of ligands. Overall, our work highlights the importance of accounting for receptor flexibility in docking studies, and provides a platform for the identification of new inhibitors against SARS-CoV-2 proteins.

Dynamics-Based Peptide-MHC Binding Optimization by a Convolutional Variational Autoencoder: A Use-Case Model for CASTELO

An unsolved challenge in the development of antigen-specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen-major histocompatibility complex (MHC) binding is paramount toward achieving this goal. Here, we apply CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach, to identify per-residue antigen binding contributions and then design novel antigens of increased MHC-II binding affinity for a type 1 diabetes-implicated system. We build upon a small-molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across four antigens and four HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and root-mean-square fluctuations (RMSF), which do not reflect binding affinity data. Our work supports the role of structure-based deep learning techniques in antigen-specific immunotherapy design.

The pockets guide to HLA class I molecules

Human leukocyte antigens (HLA) are cell-surface proteins that present peptides to T cells. These peptides are bound within the peptide binding cleft of HLA, and together as a complex, are recognised by T cells using their specialised T cell receptors. Within the cleft, the peptide residue side chains bind into distinct pockets. These pockets ultimately determine the specificity of peptide binding. As HLAs are the most polymorphic molecules in humans, amino acid variants in each binding pocket influences the peptide repertoire that can be presented on the cell surface. Here, we review each of the 6 HLA binding pockets of HLA class I (HLA-I) molecules. The binding specificity of pockets B and F are strong determinants of peptide binding and have been used to classify HLA into supertypes, a useful tool to predict peptide binding to a given HLA. Over the years, peptide binding prediction has also become more reliable by using binding affinity and mass spectrometry data. Crystal structures of peptide-bound HLA molecules provide a means to interrogate the interactions between binding pockets and peptide residue side chains. We find that most of the bound peptides from these structures conform to binding motifs determined from prediction software and examine outliers to learn how these HLAs are stabilised from a structural perspective.

Mining the Mycobacterium tuberculosis proteome for identification of potential T-cell epitope based vaccine candidates

Identification of protective antigens for designing a high-efficacy tuberculosis vaccine is the need of the hour. Till date only 7% of the Mycobacterium tuberculosis proteome has been explored for discovering antigens capable of activating T-cell responses. Therefore, it becomes crucial to screen the remaining Mycobacterium tuberculosis proteome for more immunodominant T-cell epitopes. An extensive knowledge of the epitopes recognized by our immune system can aid this process of finding potential T cell antigens for development of a better TB vaccine. In the present in-silico study, 237 proteins belonging to the 'virulence, detoxification, and adaptation' category of Mycobacterium tuberculosis proteome were targeted for T-cell epitope screening. 50825 MHC Class I and 49357 MHC Class II epitopes were generated using NetMHC3.4 and IEDB servers respectively and tested for their antigenicity and cytokine stimulation. The highest antigenic epitopes were analyzed for their world population coverage and epitope conservancy. Molecular docking and molecular dynamics simulation studies were performed to corroborate the binding affinities and structural stability of the peptide-MHC complexes. We predicted a total of 3 MHC Class I (ILLKMCWPA, FAVGMNVYV, and SLAGNSAKV) and 7 MHC Class II (DLTIGFFLHIPFPPV, RPDLTIGFFLHIPFP, LTIGFFLHIPFPPVE, VLVFALVVALVYLQF, LVFALVVALVYLQFR, PNLVAARFIQLTPVY, and LVLVFALVVALVYLQ) epitopes that can be promising vaccine candidates. These predicted epitopes belong to 6 distinct proteins: Rv0169 (mce1a), Rv3490 (ostA), Rv3496 (mce4D), Rv1085c, Rv0563 (HtpX), Rv3497c (mce4C). All these proteins are expressed at different stages in the life cycle of Mycobacterium tuberculosis and thus, the predicted epitopes could be employed as candidates for designing a multistage-multiepitopic vaccine.

Impact of Structural Observables From Simulations to Predict the Effect of Single-Point Mutations in MHC Class II Peptide Binders

The prediction of peptide binders to Major Histocompatibility Complex (MHC) class II receptors is of great interest to study autoimmune diseases and for vaccine development. Most approaches predict the affinities using sequence-based models trained on experimental data and multiple alignments from known peptide substrates. However, detecting activity differences caused by single-point mutations is a challenging task. In this work, we used interactions calculated from simulations to build scoring matrices for quickly estimating binding differences by single-point mutations. We modelled a set of 837 peptides bound to an MHC class II allele, and optimized the sampling of the conformations using the Rosetta backrub method by comparing the results to molecular dynamics simulations. From the dynamic trajectories of each complex, we averaged and compared structural observables for each amino acid at each position of the 9°mer peptide core region. With this information, we generated the scoring-matrices to predict the sign of the binding differences. We then compared the performance of the best scoring-matrix to different computational methodologies that range in computational costs. Overall, the prediction of the activity differences caused by single mutated peptides was lower than 60% for all the methods. However, the developed scoring-matrix in combination with existing methods reports an increase in the performance, up to 86% with a scoring method that uses molecular dynamics.

Complex Interaction Among Immune, Inflammatory, and Carcinogenic Mechanisms in the Head and Neck Squamous Cell Carcinoma

Inflammation is deeply involved in the development of most types of cancer. Many studies focus on the interaction between immune-inflammatory mechanisms and tumorigenesis in the head and neck squamous cell carcinoma (HNSCC). In this chapter, we emphasize the complexity of processes underlying this interaction and discuss the mechanisms of carcinogenesis in HNSCC with a special focus on metabolic changes, inflammation, and the immune landscape. Unveiling complex connections between immuno-inflammatory processes and tumor initiation, promotion, and progression will open new directions in the reliable identification of predictive factors and therapeutic targets in HNSCC.

Predicting Immunogenicity Risk in Biopharmaceuticals

The assessment of immunogenicity of biopharmaceuticals is a crucial step in the process of their development. Immunogenicity is related to the activation of adaptive immunity. The complexity of the immune system manifests through numerous different mechanisms, which allows the use of different approaches for predicting the immunogenicity of biopharmaceuticals. The direct experimental approaches are sometimes expensive and time consuming, or their results need to be confirmed. In this case, computational methods for immunogenicity prediction appear as an appropriate complement in the process of drug design. In this review, we analyze the use of various In silico methods and approaches for immunogenicity prediction of biomolecules: sequence alignment algorithms, predicting subcellular localization, searching for major histocompatibility complex (MHC) binding motifs, predicting T and B cell epitopes based on machine learning algorithms, molecular docking, and molecular dynamics simulations. Computational tools for antigenicity and allergenicity prediction also are considered.

Markov state modeling reveals alternative unbinding pathways for peptide-MHC complexes

Peptide binding to major histocompatibility complexes (MHCs) is a central component of the immune system, and understanding the mechanism behind stable peptide-MHC binding will aid the development of immunotherapies. While MHC binding is mostly influenced by the identity of the so-called anchor positions of the peptide, secondary interactions from nonanchor positions are known to play a role in complex stability. However, current MHC-binding prediction methods lack an analysis of the major conformational states and might underestimate the impact of secondary interactions. In this work, we present an atomically detailed analysis of peptide-MHC binding that can reveal the contributions of any interaction toward stability. We propose a simulation framework that uses both umbrella sampling and adaptive sampling to generate a Markov state model (MSM) for a coronavirus-derived peptide (QFKDNVILL), bound to one of the most prevalent MHC receptors in humans (HLA-A24:02). While our model reaffirms the importance of the anchor positions of the peptide in establishing stable interactions, our model also reveals the underestimated importance of position 4 (p4), a nonanchor position. We confirmed our results by simulating the impact of specific peptide mutations and validated these predictions through competitive binding assays. By comparing the MSM of the wild-type system with those of the D4A and D4P mutations, our modeling reveals stark differences in unbinding pathways. The analysis presented here can be applied to any peptide-MHC complex of interest with a structural model as input, representing an important step toward comprehensive modeling of the MHC class I pathway.

Large-Scale Structure-Based Prediction of Stable Peptide Binding to Class I HLAs Using Random Forests

Prediction of stable peptide binding to Class I HLAs is an important component for designing immunotherapies. While the best performing predictors are based on machine learning algorithms trained on peptide-HLA (pHLA) sequences, the use of structure for training predictors deserves further exploration. Given enough pHLA structures, a predictor based on the residue-residue interactions found in these structures has the potential to generalize for alleles with little or no experimental data. We have previously developed APE-Gen, a modeling approach able to produce pHLA structures in a scalable manner. In this work we use APE-Gen to model over 150,000 pHLA structures, the largest dataset of its kind, which were used to train a structure-based pan-allele model. We extract simple, homogenous features based on residue-residue distances between peptide and HLA, and build a random forest model for predicting stable pHLA binding. Our model achieves competitive AUROC values on leave-one-allele-out validation tests using significantly less data when compared to popular sequence-based methods. Additionally, our model offers an interpretation analysis that can reveal how the model composes the features to arrive at any given prediction. This interpretation analysis can be used to check if the model is in line with chemical intuition, and we showcase particular examples. Our work is a significant step toward using structure to achieve generalizable and more interpretable prediction for stable pHLA binding.

HLA-Arena: A Customizable Environment for the Structural Modeling and Analysis of Peptide-HLA Complexes for Cancer Immunotherapy

PURPOSE

HLA protein receptors play a key role in cellular immunity. They bind intracellular peptides and display them for recognition by T-cell lymphocytes. Because T-cell activation is partially driven by structural features of these peptide-HLA complexes, their structural modeling and analysis are becoming central components of cancer immunotherapy projects. Unfortunately, this kind of analysis is limited by the small number of experimentally determined structures of peptide-HLA complexes. Overcoming this limitation requires developing novel computational methods to model and analyze peptide-HLA structures.

METHODS

Here we describe a new platform for the structural modeling and analysis of peptide-HLA complexes, called HLA-Arena, which we have implemented using Jupyter Notebook and Docker. It is a customizable environment that facilitates the use of computational tools, such as APE-Gen and DINC, which we have previously applied to peptide-HLA complexes. By integrating other commonly used tools, such as MODELLER and MHCflurry, this environment includes support for diverse tasks in structural modeling, analysis, and visualization.

RESULTS

To illustrate the capabilities of HLA-Arena, we describe 3 example workflows applied to peptide-HLA complexes. Leveraging the strengths of our tools, DINC and APE-Gen, the first 2 workflows show how to perform geometry prediction for peptide-HLA complexes and structure-based binding prediction, respectively. The third workflow presents an example of large-scale virtual screening of peptides for multiple HLA alleles.

CONCLUSION

These workflows illustrate the potential benefits of HLA-Arena for the structural modeling and analysis of peptide-HLA complexes. Because HLA-Arena can easily be integrated within larger computational pipelines, we expect its potential impact to vastly increase. For instance, it could be used to conduct structural analyses for personalized cancer immunotherapy, neoantigen discovery, or vaccine development.

Uncovering the Tumor Antigen Landscape: What to Know about the Discovery Process

According to the latest available data, cancer is the second leading cause of death, highlighting the need for novel cancer therapeutic approaches. In this context, immunotherapy is emerging as a reliable first-line treatment for many cancers, particularly metastatic melanoma. Indeed, cancer immunotherapy has attracted great interest following the recent clinical approval of antibodies targeting immune checkpoint molecules, such as PD-1, PD-L1, and CTLA-4, that release the brakes of the immune system, thus reviving a field otherwise poorly explored. Cancer immunotherapy mainly relies on the generation and stimulation of cytotoxic CD8 T lymphocytes (CTLs) within the tumor microenvironment (TME), priming T cells and establishing efficient and durable anti-tumor immunity. Therefore, there is a clear need to define and identify immunogenic T cell epitopes to use in therapeutic cancer vaccines. Naturally presented antigens in the human leucocyte antigen-1 (HLA-I) complex on the tumor surface are the main protagonists in evocating a specific anti-tumor CD8+ T cell response. However, the methodologies for their identification have been a major bottleneck for their reliable characterization. Consequently, the field of antigen discovery has yet to improve. The current review is intended to define what are today known as tumor antigens, with a main focus on CTL antigenic peptides. We also review the techniques developed and employed to date for antigen discovery, exploring both the direct elution of HLA-I peptides and the in silico prediction of epitopes. Finally, the last part of the review analyses the future challenges and direction of the antigen discovery field.

pMHC Structural Comparisons as a Pivotal Element to Detect and Validate T-Cell Targets for Vaccine Development and Immunotherapy-A New Methodological Proposal

The search for epitopes that will effectively trigger an immune response remains the “El Dorado” for immunologists. The development of promising immunotherapeutic approaches requires the appropriate targets to elicit a proper immune response. Considering the high degree of HLA/TCR diversity, as well as the heterogeneity of viral and tumor proteins, this number will invariably be higher than ideal to test. It is known that the recognition of a peptide-MHC (pMHC) by the T-cell receptor is performed entirely in a structural fashion, where the atomic interactions of both structures, pMHC and TCR, dictate the fate of the process. However, epitopes with a similar composition of amino acids can produce dissimilar surfaces. Conversely, sequences with no conspicuous similarities can exhibit similar TCR interaction surfaces. In the last decade, our group developed a database and in silico structural methods to extract molecular fingerprints that trigger T-cell immune responses, mainly referring to physicochemical similarities, which could explain the immunogenic differences presented by different pMHC-I complexes. Here, we propose an immunoinformatic approach that considers a structural level of information, combined with an experimental technology that simulates the presentation of epitopes for a T cell, to improve vaccine production and immunotherapy efficacy.

Machine Learning for Cancer Immunotherapies Based on Epitope Recognition by T Cell Receptors

In the last years, immunotherapies have shown tremendous success as treatments for multiple types of cancer. However, there are still many obstacles to overcome in order to increase response rates and identify effective therapies for every individual patient. Since there are many possibilities to boost a patient's immune response against a tumor and not all can be covered, this review is focused on T cell receptor-mediated therapies. CD8+ T cells can detect and destroy malignant cells by binding to peptides presented on cell surfaces by MHC (major histocompatibility complex) class I molecules. CD4+ T cells can also mediate powerful immune responses but their peptide recognition by MHC class II molecules is more complex, which is why the attention has been focused on CD8+ T cells. Therapies based on the power of T cells can, on the one hand, enhance T cell recognition by introducing TCRs that preferentially direct T cells to tumor sites (so called TCR-T therapy) or through vaccination to induce T cells in vivo. On the other hand, T cell activity can be improved by immune checkpoint inhibition or other means that help create a microenvironment favorable for cytotoxic T cell activity. The manifold ways in which the immune system and cancer interact with each other require not only the use of large omics datasets from gene, to transcript, to protein, and to peptide but also make the application of machine learning methods inevitable. Currently, discovering and selecting suitable TCRs is a very costly and work intensive in vitro process. To facilitate this process and to additionally allow for highly personalized therapies that can simultaneously target multiple patient-specific antigens, especially neoepitopes, breakthrough computational methods for predicting antigen presentation and TCR binding are urgently required. Particularly, potential cross-reactivity is a major consideration since off-target toxicity can pose a major threat to patient safety. The current speed at which not only datasets grow and are made available to the public, but also at which new machine learning methods evolve, is assuring that computational approaches will be able to help to solve problems that immunotherapies are still facing.

Using parallelized incremental meta-docking can solve the conformational sampling issue when docking large ligands to proteins

BACKGROUND

Docking large ligands, and especially peptides, to protein receptors is still considered a challenge in computational structural biology. Besides the issue of accurately scoring the binding modes of a protein-ligand complex produced by a molecular docking tool, the conformational sampling of a large ligand is also often considered a challenge because of its underlying combinatorial complexity. In this study, we evaluate the impact of using parallelized and incremental paradigms on the accuracy and performance of conformational sampling when docking large ligands. We use five datasets of protein-ligand complexes involving ligands that could not be accurately docked by classical protein-ligand docking tools in previous similar studies.

RESULTS

Our computational evaluation shows that simply increasing the amount of conformational sampling performed by a protein-ligand docking tool, such as Vina, by running it for longer is rarely beneficial. Instead, it is more efficient and advantageous to run several short instances of this docking tool in parallel and group their results together, in a straightforward parallelized docking protocol. Even greater accuracy and efficiency are achieved by our parallelized incremental meta-docking tool, DINC, showing the additional benefits of its incremental paradigm. Using DINC, we could accurately reproduce the vast majority of the protein-ligand complexes we considered.

CONCLUSIONS

Our study suggests that, even when trying to dock large ligands to proteins, the conformational sampling of the ligand should no longer be considered an issue, as simple docking protocols using existing tools can solve it. Therefore, scoring should currently be regarded as the biggest unmet challenge in molecular docking.

Design, Synthesis and Pharmacological Evaluation of Gastro- Protective Anti-inflammatory Analgesic Agents based on Dual Oxidative Stress / Cyclooxygenase Inhibition

BACKGROUND

Non-steroidal anti-inflammatory drugs (NSAIDs) derived local generation of reactive oxygen species (ROS) plays a crucial role in the formation of gastric ulceration.

OBJECTIVE

Therefore, anti-inflammatory analgesics with potent antioxidant activity could be a potential therapeutic strategy for the treatment of pain and inflammatory disorders without gastrointestinal (GI) side effects.

METHODS

In an effort to develop gastroprotective analgesic and anti-inflammatory agents, a series of 2-methylamino-substituted-1H-benzo[d] imidazol-1-yl) (phenyl) methanone derivatives were synthesized and evaluated in vitro for cyclooxygenase (COX) inhibition as well as anti-oxidant potential by the FRAP assay. The compounds with significant in vitro COX-1/COX-2 inhibitory activity and antioxidant activity were further screened in vivo for their anti-inflammatory and analgesic activities. Moreover, the ulcerogenic potential of test compounds was also studied. To gain insight into the plausible mode of interaction of compounds within the active sites of COX-1 and COX-2, molecular docking simulations were performed.

RESULTS

Among the various synthesized molecules, most of the compounds showed good cyclooxygenase inhibitory activity and efficient antioxidant activity in FRAP assay. After preliminary and indicative in vitro assays, three compounds exhibited most significant antiinflammatory and analgesic activity with better gastric tolerability during their in vivo evaluation. Ligand interaction studies indicated highest dock score -43.05 of 1,2- disubstituted benzimidazole derivatives in comparison to the reference ligand -30.70. Overall studies provided us (2-((4-methoxyphenylamino) methyl) -1h-benzo [d] imidazol- 1-yl) (phenyl) methanone as a lead with potent gastro-protective anti-inflammatory and analgesic activities that can be used for future research.

CONCLUSION

From the above results, it can be concluded that designing of multifunctional molecules with COX-1/COX-2 inhibitory and anti-oxidant activities could hold a great promise for further development of GI-safer NSAIDs.

PeptoGrid-Rescoring Function for AutoDock Vina to Identify New Bioactive Molecules from Short Peptide Libraries

Peptides are promising drug candidates due to high specificity and standout safety. Identification of bioactive peptides de novo using molecular docking is a widely used approach. However, current scoring functions are poorly optimized for peptide ligands. In this work, we present a novel algorithm PeptoGrid that rescores poses predicted by AutoDock Vina according to frequency information of ligand atoms with particular properties appearing at different positions in the target protein’s ligand binding site. We explored the relevance of PeptoGrid ranking with a virtual screening of peptide libraries using angiotensin-converting enzyme and GABAB receptor as targets. A reasonable agreement between the computational and experimental data suggests that PeptoGrid is suitable for discovering functional leads.

Antigenic Peptide Prediction From E6 and E7 Oncoproteins of HPV Types 16 and 18 for Therapeutic Vaccine Design Using Immunoinformatics and MD Simulation Analysis

Human papillomavirus (HPV) induced cervical cancer is the second most common cause of death, after breast cancer, in females. Three prophylactic vaccines by Merck Sharp & Dohme (MSD) and GlaxoSmithKline (GSK) have been confirmed to prevent high-risk HPV strains but these vaccines have been shown to be effective only in girls who have not been exposed to HPV previously. The constitutively expressed HPV oncoproteins E6 and E7 are usually used as target antigens for HPV therapeutic vaccines. These early (E) proteins are involved, for example, in maintaining the malignant phenotype of the cells. In this study, we predicted antigenic peptides of HPV types 16 and 18, encoded by E6 and E7 genes, using an immunoinformatics approach. To further evaluate the immunogenic potential of the predicted peptides, we studied their ability to bind to class I major histocompatibility complex (MHC-I) molecules in a computational docking study that was supported by molecular dynamics (MD) simulations and estimation of the free energies of binding of the peptides at the MHC-I binding cleft. Some of the predicted peptides exhibited comparable binding free energies and/or pattern of binding to experimentally verified MHC-I-binding epitopes that we used as references in MD simulations. Such peptides with good predicted affinity may serve as candidate epitopes for the development of therapeutic HPV peptide vaccines.

Highly conserved hemagglutinin peptides of H1N1 influenza virus elicit immune response

In the current study, two highly conserved (> 90%) H1N1 hemagglutinin peptides STDTVDTVLEKNVTVTHSVNL (H1) and KVNSVIEKMNTQFTAVGKEF (H2) containing multiple T-cell epitopes have been assessed for their immunogenic potential in vitro, subjecting peripheral blood mononuclear cells from healthy volunteers to repetitive stimulation of chemically synthesised H1 and H2 peptides, and measuring their interferon (IFN)-γ level (ELISA) and proliferation (MTT assay). Further, these peptides were analysed for their binding affinity with 18 different human leukocyte antigen (HLA) class I and II by means of molecular docking. All seven samples tested for H1- and H2-induced IFN-γ secretion were found to have enhanced IFN-γ production. Six (H1) and five (H2) samples have shown proliferative response compared to unstimulated cells. Peptide-induced IFN-γ secretion and proliferation in healthy samples represent the immunogenic potential of these peptides. Further, molecular docking results reveal that the peptides have comparable binding energy to that of native bound peptide for both HLA classes which indicates that these peptides have the capability to be presented by different HLA molecules required for T-cell response. Hence, these conserved immunogenic hemagglutinin peptides are potential candidates for influenza vaccine development.

Structure-based Methods for Binding Mode and Binding Affinity Prediction for Peptide-MHC Complexes

Understanding the mechanisms involved in the activation of an immune response is essential to many fields in human health, including vaccine development and personalized cancer immunotherapy. A central step in the activation of the adaptive immune response is the recognition, by T-cell lymphocytes, of peptides displayed by a special type of receptor known as Major Histocompatibility Complex (MHC). Considering the key role of MHC receptors in T-cell activation, the computational prediction of peptide binding to MHC has been an important goal for many immunological applications. Sequence- based methods have become the gold standard for peptide-MHC binding affinity prediction, but structure-based methods are expected to provide more general predictions (i.e., predictions applicable to all types of MHC receptors). In addition, structural modeling of peptide-MHC complexes has the potential to uncover yet unknown drivers of T-cell activation, thus allowing for the development of better and safer therapies. In this review, we discuss the use of computational methods for the structural modeling of peptide-MHC complexes (i.e., binding mode prediction) and for the structure-based prediction of binding affinity.