HLA-DP2 binding prediction by molecular dynamics simulations

High throughput virtual screening (HTVS) of peptide library: Technological advancement in ligand discovery

High-throughput virtual screening (HTVS) is a leading biopharmaceutical technology that employs computational algorithms to uncover biologically active compounds from large-scale collections of chemical compound libraries. In addition, this method often leverages the precedence of screening focused libraries for assessing their binding affinities and improving physicochemical properties. Usually, developing a drug sometimes takes ages, and lessons are learnt from FDA-approved drugs. This screening strategy saves resources and time compared to laboratory testing in certain stages of drug discovery. Yet in-silico investigations remain challenging in some cases of drug discovery. For the last few decades, peptide-based drug discoveries have received remarkable momentum for several advantages over small molecules. Therefore, developing a high-fidelity HTVS platform for chemically versatile peptide libraries is highly desired. This review summarises the modern and frequently appreciated HTVS strategies for peptide libraries from 2011 to 2021. In addition, we focus on the software used for preparing peptide libraries, their screening techniques and shortcomings. An index of various HTVS methods reported here should assist researchers in identifying tools that could be beneficial for their peptide library screening projects.

MHCII3D-Robust Structure Based Prediction of MHC II Binding Peptides

Knowledge of MHC II binding peptides is highly desired in immunological research, particularly in the context of cancer, autoimmune diseases, or allergies. The most successful prediction methods are based on machine learning methods trained on sequences of experimentally characterized binding peptides. Here, we describe a complementary approach called MHCII3D, which is based on structural scaffolds of MHC II-peptide complexes and statistical scoring functions (SSFs). The MHC II alleles reported in the Immuno Polymorphism Database are processed in a dedicated 3D-modeling pipeline providing a set of scaffold complexes for each distinct allotype sequence. Antigen protein sequences are threaded through the scaffolds and evaluated by optimized SSFs. We compared the predictive power of MHCII3D with different sequence-based machine learning methods. The Pearson correlation to experimentally determine IC₅₀ values for MHC II Automated Server Benchmarks data sets from IEDB (Immune Epitope Database) is 0.42, which is in the competitor methods range. We show that MHCII3D is quite robust in leaving one molecule out tests and is therefore not prone to overfitting. Finally, we provide evidence that MHCII3D can complement the current sequence-based methods and help to identify problematic entries in IEDB. Scaffolds and MHCII3D executables can be freely downloaded from our web pages.

Molecular mechanism of Be2+-ion binding to HLA-DP2: tetrahedral coordination, conformational changes and multi-ion binding

Be small and positive: the smaller size and higher charge of the Be²⁺-ion results in strong binding between the M2 peptide and the β-chain of HLA-DP2, which induces conformational changes at the periphery suitable for TCR binding.

Ab-initio conformational epitope structure prediction using genetic algorithm and SVM for vaccine design

BACKGROUND AND OBJECTIVE

T-cell epitope structure identification is a significant challenging immunoinformatic problem within epitope-based vaccine design. Epitopes or antigenic peptides are a set of amino acids that bind with the Major Histocompatibility Complex (MHC) molecules. The aim of this process is presented by Antigen Presenting Cells to be inspected by T-cells. MHC-molecule-binding epitopes are responsible for triggering the immune response to antigens. The epitope's three-dimensional (3D) molecular structure (i.e., tertiary structure) reflects its proper function. Therefore, the identification of MHC class-II epitopes structure is a significant step towards epitope-based vaccine design and understanding of the immune system.

METHODS

In this paper, we propose a new technique using a Genetic Algorithm for Predicting the Epitope Structure (GAPES), to predict the structure of MHC class-II epitopes based on their sequence. The proposed Elitist-based genetic algorithm for predicting the epitope's tertiary structure is based on Ab-Initio Empirical Conformational Energy Program for Peptides (ECEPP) Force Field Model. The developed secondary structure prediction technique relies on Ramachandran Plot. We used two alignment algorithms: the ROSS alignment and TM-Score alignment. We applied four different alignment approaches to calculate the similarity scores of the dataset under test. We utilized the support vector machine (SVM) classifier as an evaluation of the prediction performance.

RESULTS

The prediction accuracy and the Area Under Receiver Operating Characteristic (ROC) Curve (AUC) were calculated as measures of performance. The calculations are performed on twelve similarity-reduced datasets of the Immune Epitope Data Base (IEDB) and a large dataset of peptide-binding affinities to HLA-DRB1*0101. The results showed that GAPES was reliable and very accurate. We achieved an average prediction accuracy of 93.50% and an average AUC of 0.974 in the IEDB dataset. Also, we achieved an accuracy of 95.125% and an AUC of 0.987 on the HLA-DRB1*0101 allele of the Wang benchmark dataset.

CONCLUSIONS

The results indicate that the proposed prediction technique "GAPES" is a promising technique that will help researchers and scientists to predict the protein structure and it will assist them in the intelligent design of new epitope-based vaccines.

Elucidation, functional clustering and structural characterization of βTrCP1 substrates through a molecular dynamics study

Structural knowledge of substrate recognition by SKP1–CUL1–βTrCP1 complex for targeted cancer therapeutic strategy.

Peptides and Peptide Analogs to Inhibit Protein-Protein Interactions

Protein-protein interactions are governed by relatively few amino acid residues at the binding interface. Peptides derived from these protein regions may serve as mimics of one of the interaction partners in structural studies or as inhibitors to disrupt the respective interaction and investigate its biological consequences. Inhibitory peptides may also be lead structures for drug development if the respective protein-protein interaction is essential for a pathogen or disease mechanism. Binding peptides may be systematically derived from one of the binding partners or found in the screen of combinatorial peptide libraries. Molecular modelling based on structural data helps to refine existing peptides or even design novel binding peptides. This chapter gives an outline of the binding peptide discovery process and subsequent chemical modifications to further enhance affinity and specificity and to increase stability against degradation in vivo. Examples from the past three decades illustrate the great diversity of applications for protein binding peptides and peptide analogs.

Chronic Beryllium Disease: revealing the role of beryllium ion and small peptides binding to HLA-DP2

Chronic Beryllium (Be) Disease (CBD) is a granulomatous disorder that predominantly affects the lung. The CBD is caused by Be exposure of individuals carrying the HLA-DP2 protein of the major histocompatibility complex class II (MHCII). While the involvement of Be in the development of CBD is obvious and the binding site and the sequence of Be and peptide binding were recently experimentally revealed [1], the interplay between induced conformational changes and the changes of the peptide binding affinity in presence of Be were not investigated. Here we carry out in silico modeling and predict the Be binding to be within the acidic pocket (Glu26, Glu68 and Glu69) present on the HLA-DP2 protein in accordance with the experimental work [1]. In addition, the modeling indicates that the Be ion binds to the HLA-DP2 before the corresponding peptide is able to bind to it. Further analysis of the MD generated trajectories reveals that in the presence of the Be ion in the binding pocket of HLA-DP2, all the different types of peptides induce very similar conformational changes, but their binding affinities are quite different. Since these conformational changes are distinctly different from the changes caused by peptides normally found in the cell in the absence of Be, it can be speculated that CBD can be caused by any peptide in presence of Be ion. However, the affinities of peptides for Be loaded HLA-DP2 were found to depend of their amino acid composition and the peptides carrying acidic group at positions 4 and 7 are among the strongest binders. Thus, it is proposed that CBD is caused by the exposure of Be of an individual carrying the HLA-DP2*0201 allele and that the binding of Be to HLA-DP2 protein alters the conformational and ionization properties of HLA-DP2 such that the binding of a peptide triggers a wrong signaling cascade.

Quantitative prediction of peptide binding to HLA-DP1 protein

The exogenous proteins are processed by the host antigen-processing cells. Peptidic fragments of them are presented on the cell surface bound to the major hystocompatibility complex (MHC) molecules class II and recognized by the CD4+ T lymphocytes. The MHC binding is considered as the crucial prerequisite for T-cell recognition. Only peptides able to form stable complexes with the MHC proteins are recognized by the T-cells. These peptides are known as T-cell epitopes. All T-cell epitopes are MHC binders, but not all MHC binders are T-cell epitopes. The T-cell epitope prediction is one of the main priorities of immunoinformatics. In the present study, three chemometric techniques are combined to derive a model for in silico prediction of peptide binding to the human MHC class II protein HLA-DP1. The structures of a set of known peptide binders are described by amino acid z-descriptors. Data are processed by an iterative self-consisted algorithm using the method of partial least squares, and a quantitative matrix (QM) for peptide binding prediction to HLA-DP1 is derived. The QM is validated by two sets of proteins and showed an average accuracy of 86 percent.

Structure-based prediction of Major Histocompatibility Complex (MHC) epitopes

Because of the enormous diversity of both MHC proteins and peptide epitopes, computational epitope prediction methods are needed in order to supplement limited experimental data. These prediction methods are useful for guiding experiments and have many potential biomedical applications. Unlike popular sequence-based methods, structure-based epitope prediction methods can predict epitopes for multiple MHC types with highly distinct peptide binding propensities. In this chapter, we describe in detail our previously developed structure-based epitope prediction methods for both class I and class II MHC proteins. We also discuss the relative advantages and disadvantages of sequence-based versus structure-based methods and how to evaluate prediction performance.