Quantitative structure-activity relationships and the prediction of MHC supermotifs

The underlying assumption in quantitative structure–activity relationship (QSAR) methodology is that related chemical structures exhibit related biological activities. We review here two QSAR methods in terms of their applicability for human MHC supermotif definition. Supermotifs are motifs that characterise binding to more than one allele. Supermotif definition is the initial in silico step of epitope-based vaccine design. The first QSAR method we review here—the additive method—is based on the assumption that the binding affinity of a peptide depends on contributions from both amino acids and the interactions between them. The second method is a 3D-QSAR method: comparative molecular similarity indices analysis (CoMSIA). Both methods were applied to 771 peptides binding to 9 HLA alleles. Five of the alleles (A*0201, A* 0202, A*0203, A*0206 and A*6802) belong to the HLA-A2 superfamily and the other four (A*0301, A*1101, A*3101 and A*6801) to the HLA-A3 superfamily. For each superfamily, supermotifs defined by the two QSAR methods agree closely and are supported by many experimental data.

Transporter associated with antigen processing preselection of peptides binding to the MHC:a bioinformatic evaluation

TAP is responsible for the transit of peptides from the cytosol to the lumen of the endoplasmic reticulum. In an immunological context, this event is followed by the binding of peptides to MHC molecules before export to the cell surface and recognition by T cells. Because TAP transport precedes MHC binding, TAP preferences may make a significant contribution to epitope selection. To assess the impact of this preselection, we have developed a scoring function for TAP affinity prediction using the additive method, have used it to analyze and extend the TAP binding motif, and have evaluated how well this model acts as a preselection step in predicting MHC binding peptides. To distinguish between MHC alleles that are exclusively dependent on TAP and those exhibiting only a partial dependence on TAP, two sets of MHC binding peptides were examined: HLA-A*0201 was selected as a representative of partially TAP-dependent HLA alleles, and HLA-A*0301 represented fully TAP-dependent HLA alleles. TAP preselection has a greater impact on TAP-dependent alleles than on TAP-independent alleles. The reduction in the number of nonbinders varied from 10% (TAP-independent) to 33% (TAP-dependent), suggesting that TAP preselection is an important component in the successful in silico prediction of T cell epitopes.

New horizons in mouse immunoinformatics:reliable in silico prediction of mouse class I histocompatibility major complex peptide binding affinity

Quantitative structure–activity relationship (QSAR) analysis is a main cornerstone of modern informatic disciplines. Predictive computational models, based on QSAR technology, of peptide-major histocompatibility complex (MHC) binding affinity have now become a vital component of modern day computational immunovaccinology. Historically, such approaches have been built around semi-qualitative, classification methods, but these are now giving way to quantitative regression methods. The additive method, an established immunoinformatics technique for the quantitative prediction of peptide–protein affinity, was used here to identify the sequence dependence of peptide binding specificity for three mouse class I MHC alleles: H2–Db, H2–Kb and H2–Kk. As we show, in terms of reliability the resulting models represent a significant advance on existing methods. They can be used for the accurate prediction of T-cell epitopes and are freely available online (http://www.jenner.ac.uk/MHCPred).

Coupling in silico and in vitro analysis of peptide-MHC binding:a bioinformatic approach enabling prediction of superbinding peptides and anchorless epitopes

The ability to define and manipulate the interaction of peptides with MHC molecules has immense immunological utility, with applications in epitope identification, vaccine design, and immunomodulation. However, the methods currently available for prediction of peptide-MHC binding are far from ideal. We recently described the application of a bioinformatic prediction method based on quantitative structure-affinity relationship methods to peptide-MHC binding. In this study we demonstrate the predictivity and utility of this approach. We determined the binding affinities of a set of 90 nonamer peptides for the MHC class I allele HLA-A*0201 using an in-house, FACS-based, MHC stabilization assay, and from these data we derived an additive quantitative structure-affinity relationship model for peptide interaction with the HLA-A*0201 molecule. Using this model we then designed a series of high affinity HLA-A2-binding peptides. Experimental analysis revealed that all these peptides showed high binding affinities to the HLA-A*0201 molecule, significantly higher than the highest previously recorded. In addition, by the use of systematic substitution at principal anchor positions 2 and 9, we showed that high binding peptides are tolerant to a wide range of nonpreferred amino acids. Our results support a model in which the affinity of peptide binding to MHC is determined by the interactions of amino acids at multiple positions with the MHC molecule and may be enhanced by enthalpic cooperativity between these component interactions.

Quantitative online prediction of peptide binding to the major histocompatibility complex

With its implications for vaccine discovery, the accurate prediction of T cell epitopes is one of the key aspirations of computational vaccinology. We have developed a robust multivariate statistical method, based on partial least squares, for the quantitative prediction of peptide binding to major histocompatibility complexes (MHC), the principal checkpoint on the antigen presentation pathway. As a service to the immunobiology community, we have made a Perl implementation of the method available via a World Wide Web server. We call this server MHCPred. Access to the server is freely available from the URL: http://www.jenner.ac.uk/MHCPred. We have exemplified our method with a model for peptides binding to the common human MHC molecule HLA-B*3501.

Selection of conserved epitopes from hepatitis C virus for pan-populational stimulation of T-cell responses

The hepatitis C virus (HCV) is able to persist as a chronic infection, which can lead to cirrhosis and liver cancer. There is evidence that clearance of HCV is linked to strong responses by CD8 cytotoxic T lymphocytes (CTLs), suggesting that eliciting CTL responses against HCV through an epitope-based vaccine could prove an effective means of immunization. However, HCV genomic plasticity as well as the polymorphisms of HLA I molecules restricting CD8 T-cell responses challenges the selection of epitopes for a widely protective vaccine. Here, we devised an approach to overcome these limitations. From available databases, we first collected a set of 245 HCV-specific CD8 T-cell epitopes, all known to be targeted in the course of a natural infection in humans. After a sequence variability analysis, we next identified 17 highly invariant epitopes. Subsequently, we predicted the epitope HLA I binding profiles that determine their potential presentation and recognition. Finally, using the relevant HLA I-genetic frequencies, we identified various epitope subsets encompassing 6 conserved HCV-specific CTL epitopes each predicted to elicit an effective T-cell response in any individual regardless of their HLA I background. We implemented this epitope selection approach for free public use at the EPISOPT web server. © 2013 Magdalena Molero-Abraham et al.

MHC class II binding prediction by molecular docking

Proteins of the Major Histocompatibility Complex (MHC) bind self and nonself peptide antigens or epitopes within the cell and present them at the cell surface for recognition by T cells. All T-cell epitopes are MHC binders but not all MCH binders are T-cell epitopes. The MHC class II proteins are extremely polymorphic. Polymorphic residues cluster in the peptide-binding region and largely determine the MHC's peptide selectivity. The peptide binding site on MHC class II proteins consist of five binding pockets. Using molecular docking, we have modelled the interactions between peptide and MHC class II proteins from locus DRB1. A combinatorial peptide library was generated by mutation of residues at peptide positions which correspond to binding pockets (so called anchor positions). The binding affinities were assessed using different scoring functions. The normalized scoring functions for each amino acid at each anchor position were used to construct quantitative matrices (QM) for MHC class II binding prediction. Models were validated by external test sets comprising 4540 known binders. Eighty percent of the known binders are identified in the best predicted 15% of all overlapping peptides, originating from one protein. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Towards the knowledge-based design of universal influenza epitope ensemble vaccines

Motivation: Influenza A viral heterogeneity remains a significant threat due to unpredictable antigenic drift in seasonal influenza and antigenic shifts caused by the emergence of novel subtypes. Annual review of multivalent influenza vaccines targets strains of influenza A and B likely to be predominant in future influenza seasons. This does not induce broad, cross protective immunity against emergent subtypes. Better strategies are needed to prevent future pandemics. Cross-protection can be achieved by activating CD8+ and CD4+ T cells against highly-conserved regions of the influenza genome. We combine available experimental data with informatics-based immunological predictions to help design vaccines potentially able to induce cross-protective T-cells against multiple influenza subtypes. Results: To exemplify our approach we designed two epitope ensemble vaccines comprising highly-conserved and experimentally-verified immunogenic influenza A epitopes as putative non-seasonal influenza vaccines; one specifically targets the US population and the other is a universal vaccine. The USA-specific vaccine comprised 6 CD8+ T cell epitopes (GILGFVFTL, FMYSDFHFI, GMDPRMCSL, SVKEKDMTK, FYIQMCTEL, DTVNRTHQY) and 3 CD4+ epitopes (KGILGFVFTLTVPSE, EYIMKGVYINTALLN, ILGFVFTLTVPSERG). The universal vaccine comprised 8 CD8+ epitopes: (FMYSDFHFI, GILGFVFTL, ILRGSVAHK, FYIQMCTEL, ILKGKFQTA, YYLEKANKI, VSDGGPNLY, YSHGTGTGY) and the same 3 CD4+ epitopes. Our USA-specific vaccine has a population protection coverage (portion of the population potentially responsive to one or more component epitopes of the vaccine, PPC) of over 96% and 95% coverage of observed influenza subtypes. The universal vaccine has a PPC value of over 97% and 88% coverage of observed subtypes.

Multiparameter flow cytometric analysis of C -reactive protein binding to human-derived cell lines and peripheral blood monocytes

C-reactive protein (CRP), a normally occurring human plasma protein may become elevated as much as 1,000 fold during disease states involving acute inflammation or tissue damage. Through its binding to phosphorylcholine in the presence of calcium, CRP has been shown to potentiate the activation of complement, stimulate phagocytosis and opsonize certain microorganisms. Utilizing a flow cytometric functional ligand binding assay I have demonstrated that a monocyte population in human peripheral blood and specific human-derived myelomonocytic cell lines reproducibly bind an evolutionarily conserved conformational pentraxin epitope on human CRP through a mechanism that does not involve its ligand, phosphorylcholine. ^ A variety of cell lines at different stages of differentiation were examined. The monocytic cell line, THP-1, bound the most CRP followed by U937 and KG-1a cells. The HL-60 cell line was induced towards either the granulocyte or monocyte pathway with DMSO or PMA, respectively. Untreated HL-60 cells or DMSO-treated cells did not bind CRP while cells treated with PMA showed increased binding of CRP, similar to U-937 cells. T cell and B-cell derived lines were negative. ^ Inhibition studies with Limulin and human SAP demonstrated that the binding site is a conserved pentraxin epitope. The calcium requirement necessary for binding to occur indicated that the cells recognize a conformational form of CRP. Phosphorylcholine did not inhibit the reaction therefore the possibility that CRP had bound to damaged membranes with exposed PC sites was discounted. ^ A study of 81 normal donors using flow cytometry demonstrated that a majority of peripheral blood monocytes (67.9 ± 1.3, mean ± sem) bound CRP. The percentage of binding was normally distributed and not affected by gender, age or ethnicity. Whole blood obtained from donors representing a variety of disease states showed a significant reduction in the level of CRP bound by monocytes in those donors classified with infection, inflammation or cancer. This reduction in monocyte populations binding CRP did not correlate with the concentration of plasma CRP. ^ The ability of monocytes to specifically bind CRP combined with the binding reactivity of the protein itself to a variety of phosphorylcholine containing substances may represent an important bridge between innate and adaptive immunity. ^

A deimmunised form of the ribotoxin, α-sarcin, lacking CD4+ T cell epitopes and its use as an immunotoxin warhead

Fungal ribotoxins that block protein synthesis can be useful warheads in the context of a targeted immunotoxin. α-Sarcin is a small (17 kDa) fungal ribonuclease produced by Aspergillus giganteus that functions by catalytically cleaving a single phosphodiester bond in the sarcin–ricin loop of the large ribosomal subunit, thus making the ribosome unrecognisable to elongation factors and leading to inhibition of protein synthesis. Peptide mapping using an ex vivo human T cell assay determined that α-sarcin contained two T cell epitopes; one in the N-terminal 20 amino acids and the other in the C-terminal 20 amino acids. Various mutations were tested individually within each epitope and then in combination to isolate deimmunised α-sarcin variants that had the desired properties of silencing T cell epitopes and retention of the ability to inhibit protein synthesis (equivalent to wild-type, WT α-sarcin). A deimmunised variant (D9T/Q142T) demonstrated a complete lack of T cell activation in in vitro whole protein human T cell assays using peripheral blood mononuclear cells from donors with diverse HLA allotypes. Generation of an immunotoxin by fusion of the D9T/Q142T variant to a single-chain Fv targeting Her2 demonstrated potent cell killing equivalent to a fusion protein comprising the WT α-sarcin. These results represent the first fungal ribotoxin to be deimmunised with the potential to construct a new generation of deimmunised immunotoxin therapeutics.

Gas Solubility and Vapor Pressure of Electrolytes Based on Alkylcarbonates Containing Litfsi, Lifap, and LiPF6 Lithium Salt: Measurement and Prediction

Most liquid electrolytes used in commercial lithium-ion batteries are composed by alkylcarbonate mixture containing lithium salt. The decomposition of these solvents by oxidation or reduction during cycling of the cell, induce generation of gases (CO2, CH4, C2H4, CO …) increasing of pressure in the sealed cell, which causes a safety problem [1]. The prior understanding of parameters, such as structure and nature of salt, temperature pressure, concentration, salting effects and solvation parameters, which influence gas solubility and vapor pressure of electrolytes is required to formulate safer and suitable electrolytes especially at high temperature.

We present in this work the CO2, CH4, C2H4, CO solubility in different pure alkyl-carbonate solvents (PC, DMC, EMC, DEC) and their binary or ternary mixtures as well as the effect of temperature and lithium salt LiX (X = LiPF6, LiTFSI or LiFAP) structure and concentration on these properties. Furthermore, in order to understand parameters that influence the choice of the structure of the solvents and their ability to dissolve gas through the addition of a salt, we firstly analyzed experimentally the transport properties (Self diffusion coefficient (D), fluidity (h-1), and conductivity (s) and lithium transport number (tLi) using the Stock-Einstein, and extended Jones-Dole equations [2]. Furthermore, measured data for the of CO2, C2H4, CH4 and CO solubility in pure alkylcarbonates and their mixtures containing LiPF6; LiFAP; LiTFSI salt, are reported as a function of temperature and concentration in salt. Based on experimental solubility data, the Henry’s law constant of gases in these solvents and electrolytes was then deduced and compared with values predicted by using COSMO-RS methodology within COSMOthermX software. From these results, the molar thermodynamic functions of dissolution such as the standard Gibbs energy, the enthalpy, and the entropy, as well as the mixing enthalpy of the solvents and electrolytes with the gases in its hypothetical liquid state were calculated and discussed [3]. Finally, the analysis of the CO2 solubility variations with the salt addition was then evaluated by determining specific ion parameters Hi by using the Setchenov coefficients in solution. This study showed that the gas solubility is entropy driven and can been influenced by the shape, charge density, and size of the anions in lithium salt.

References

[1] S.A. Freunberger, Y. Chen, Z. Peng, J.M. Griffin, L.J. Hardwick, F. Bardé, P. Novák, P.G. Bruce, Journal of the American Chemical Society 133 (2011) 8040-8047.

[2] P. Porion, Y.R. Dougassa, C. Tessier, L. El Ouatani, J. Jacquemin, M. Anouti, Electrochimica Acta 114 (2013) 95-104.

[3] Y.R. Dougassa, C. Tessier, L. El Ouatani, M. Anouti, J. Jacquemin, The Journal of Chemical Thermodynamics 61 (2013) 32-44.

Common polygenic variation enhances risk prediction for Alzheimer's disease.

The identification of subjects at high risk for Alzheimer’s disease is important for prognosis and early intervention. We investigated the polygenic architecture of Alzheimer’s disease and the accuracy of Alzheimer’s disease prediction models, including and excluding the polygenic component in the model. This study used genotype data from the powerful dataset comprising 17 008 cases and 37 154 controls obtained from the International Genomics of Alzheimer’s Project (IGAP). Polygenic score analysis tested whether the alleles identified to associate with disease in one sample set were significantly enriched in the cases relative to the controls in an independent sample. The disease prediction accuracy was investigated in a subset of the IGAP data, a sample of 3049 cases and 1554 controls (for whom APOE genotype data were available) by means of sensitivity, specificity, area under the receiver operating characteristic curve (AUC) and positive and negative predictive values. We observed significant evidence for a polygenic component enriched in Alzheimer’s disease (P = 4.9 × 10−26). This enrichment remained significant after APOE and other genome-wide associated regions were excluded (P = 3.4 × 10−19). The best prediction accuracy AUC = 78.2% (95% confidence interval 77–80%) was achieved by a logistic regression model with APOE, the polygenic score, sex and age as predictors. In conclusion, Alzheimer’s disease has a significant polygenic component, which has predictive utility for Alzheimer’s disease risk and could be a valuable research tool complementing experimental designs, including preventative clinical trials, stem cell selection and high/low risk clinical studies. In modelling a range of sample disease prevalences, we found that polygenic scores almost doubles case prediction from chance with increased prediction at polygenic extremes.

From genes to radioresistance in head and neck squamous cell carcinoma

Head and Neck Cancers (HNC) are a group of tumours located in the upper aero-digestive tract. Head and Neck Squamous Cell Carcinoma (HNSCC) represent about 90% of all HNC cases. It has been considered the sixth most malignant tumour worldwide and, despite clinical and technological advances, the five-year survival rate has not improved much in the last years. Nowadays, HNSCC is well established as a heterogeneous disease and that its development is due to accumulation of genetic events. Apart from the majority of the patients being diagnosed in an advanced stage, HNSCC is also a disease with poor therapeutic outcome. One of the therapeutic approaches is radiotherapy. However, this approach has different drawbacks like the radioresistance acquired by some tumour cells, leading to a worse prognosis. A major knowledge in radiation biology is imperative to improve this type of treatment and avoid late toxicities, maintaining patient quality of life in the subsequent years after treatment. Then, identification of genetic markers associated to radiotherapy response in patients and possible alterations in cells after radiotherapy are essential steps towards an improved diagnosis, higher survival rate and a better life quality. Not much is known about the radiation effects on cells, so, the principal aim of this study was to contribute to a more extensive knowledge about radiation treatment in HNSCC. For this, two commercial cell lines, HSC-3 and BICR-10, were used and characterized resorting to karyotyping, aCGH and MS-MLPA. These cell lines were submitted to different doses of irradiation and the resulting genetic and methylation alterations were evaluated. Our results showed a great difference in radiation response between the two cell lines, allowing the conclusion that HSC-3 was much more radiosensitive than BICR-10. Bearing this in mind, analysis of cell death, cell cycle and DNA damages was performed to try to elucidate the motifs behind this difference. The characterization of both cell lines allowed the confirmation that HSC-3 was derived from a metastatic tumour and the hypothesis that BICR-10 was derived from a dysplasia. Furthermore, this pilot study enabled the suggestion of some genetic and epigenetic alterations that cells suffer after radiation treatment. Additionally, it also allowed the association of some genetic characteristics that could be related to the differences in radiation response observable in this two cell lines. Taken together all of our results contribute to a better understanding of radiation effects on HNSCC allowing one further step towards the prediction of patients’ outcome, better choice of treatment approaches and ultimately a better quality of life.

Towards the Knowledge-based Design of Universal Influenza Epitope Ensemble Vaccines

Motivation: Influenza A viral heterogeneity remains a significant threat due to unpredictable antigenic drift in seasonal influenza and antigenic shifts caused by the emergence of novel subtypes. Annual review of multivalent influenza vaccines targets strains of influenza A and B likely to be predominant in future influenza seasons. This does not induce broad, cross protective immunity against emergent subtypes. Better strategies are needed to prevent future pandemics. Cross-protection can be achieved by activating CD8+ and CD4+ T cells against highly-conserved regions of the influenza genome. We combine available experimental data with informatics-based immunological predictions to help design vaccines potentially able to induce cross-protective T-cells against multiple influenza subtypes. Results: To exemplify our approach we designed two epitope ensemble vaccines comprising highlyconserved and experimentally-verified immunogenic influenza A epitopes as putative non-seasonal influenza vaccines; one specifically targets the US population and the other is a universal vaccine. The USA-specific vaccine comprised 6 CD8+ T cell epitopes (GILGFVFTL, FMYSDFHFI, GMDPRMCSL, SVKEKDMTK, FYIQMCTEL, DTVNRTHQY) and 3 CD4+ epitopes (KGILGFVFTLTVPSE, EYIMKGVYINTALLN, ILGFVFTLTVPSERG). The universal vaccine comprised 8 CD8+ epitopes: (FMYSDFHFI, GILGFVFTL, ILRGSVAHK, FYIQMCTEL, ILKGKFQTA, YYLEKANKI, VSDGGPNLY, YSHGTGTGY) and the same 3 CD4+ epitopes. Our USA-specific vaccine has a population protection coverage (portion of the population potentially responsive to one or more component epitopes of the vaccine, PPC) of over 96% and 95% coverage of observed influenza subtypes. The universal vaccine has a PPC value of over 97% and 88% coverage of observed subtypes.

Mineral stress response of Vitis cell wall transcriptome: identification, characterization and expression of genes from key families involved in wall biosynthesis and modification

Doutoramento em Engenharia Agronómica - Instituto Superior de Agronomia - UL