Structure-Based Modeling of SARS-CoV-2 Peptide/HLA-A02 Antigens

Conformational plasticity of RAS Q61 family of neoepitopes results in distinct features for targeted recognition

The conformational landscapes of peptide/human leucocyte antigen (pHLA) protein complexes encompassing tumor neoantigens provide a rationale for target selection towards autologous T cell, vaccine, and antibody-based therapeutic modalities. Here, using complementary biophysical and computational methods, we characterize recurrent RAS55-64 Q61 neoepitopes presented by the common HLA-A*01:01 allotype. We integrate sparse NMR restraints with Rosetta docking to determine the solution structure of NRASQ61K/HLA-A*01:01, which enables modeling of other common RAS55-64 neoepitopes. Hydrogen/deuterium exchange mass spectrometry experiments alongside molecular dynamics simulations reveal differences in solvent accessibility and conformational plasticity across a panel of common Q61 neoepitopes that are relevant for recognition by immunoreceptors. Finally, we predict binding and provide structural models of NRASQ61K antigens spanning the entire HLA allelic landscape, together with in vitro validation for HLA-A*01:191, HLA-B*15:01, and HLA-C*08:02. Our work provides a basis to delineate the solution surface features and immunogenicity of clinically relevant neoepitope/HLA targets for cancer therapy.

Targeting of intracellular oncoproteins with peptide-centric CARs

The majority of oncogenic drivers are intracellular proteins, constraining their immunotherapeutic targeting to mutated peptides (neoantigens) presented by individual human leukocyte antigen (HLA) allotypes1. However, most cancers have a modest mutational burden that is insufficient for generating responses using neoantigen-based therapies2,3. Neuroblastoma is a paediatric cancer that harbours few mutations and is instead driven by epigenetically deregulated transcriptional networks4. Here we show that the neuroblastoma immunopeptidome is enriched with peptides derived from proteins essential for tumorigenesis. We focused on targeting the unmutated peptide QYNPIRTTF discovered on HLA-A*24:02, which is derived from the neuroblastoma-dependency gene and master transcriptional regulator PHOX2B. To target QYNPIRTTF, we developed peptide-centric chimeric antigen receptors (PC-CARs) through a counter panning strategy using predicted potentially cross-reactive peptides. We further proposed that PC-CARs can recognize peptides on additional HLA allotypes when presenting a similar overall molecular surface. Informed by our computational modelling results, we show that PHOX2B PC-CARs also recognize QYNPIRTTF presented by HLA-A*23:01, the most common non-A2 allele in people with African ancestry. Finally, we demonstrate potent and specific killing of neuroblastoma cells expressing these HLAs in vitro and complete tumour regression in mice. These data suggest that PC-CARs have the potential to expand the pool of immunotherapeutic targets to include non-immunogenic intracellular oncoproteins and allow targeting through additional HLA allotypes in a clinical setting.

HLA3DB: comprehensive annotation of peptide/HLA complexes enables blind structure prediction of T cell epitopes

The class I proteins of the major histocompatibility complex (MHC-I) display epitopic peptides derived from endogenous proteins on the cell surface for immune surveillance. Accurate modeling of peptides bound to the human MHC, HLA, has been mired by conformational diversity of the central peptide residues, which are critical for recognition by T cell receptors. Here, analysis of X-ray crystal structures within our curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer pHLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work may be applied towards predicting antigen immunogenicity, and receptor cross-reactivity.

HLA3DB: comprehensive annotation of peptide/HLA complexes enables blind structure prediction of T cell epitopes

The class I proteins of the major histocompatibility complex (MHC-I) display epitopic peptides derived from endogenous proteins on the cell surface for immune surveillance. Accurate modeling of peptide/HLA (pHLA, the human MHC) structures has been mired by conformational diversity of the central peptide residues, which are critical for recognition by T cell receptors. Here, analysis of X-ray crystal structures within a curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these representative backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer peptide/HLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in terms of structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work provide a framework for linking conformational diversity with antigen immunogenicity and receptor cross-reactivity.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

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.

Cross-HLA targeting of intracellular oncoproteins with peptide-centric CARs

The majority of oncogenic drivers are intracellular proteins, thus constraining their immunotherapeutic targeting to mutated peptides (neoantigens) presented by individual human leukocyte antigen (HLA) allotypes1. However, most cancers have a modest mutational burden that is insufficient to generate responses using neoantigen-based therapies2,3. Neuroblastoma is a paediatric cancer that harbours few mutations and is instead driven by epigenetically deregulated transcriptional networks4. Here we show that the neuroblastoma immunopeptidome is enriched with peptides derived from proteins that are essential for tumourigenesis and focus on targeting the unmutated peptide QYNPIRTTF, discovered on HLA-A*24:02, which is derived from the neuroblastoma dependency gene and master transcriptional regulator PHOX2B. To target QYNPIRTTF, we developed peptide-centric chimeric antigen receptors (CARs) using a counter-panning strategy with predicted potentially cross-reactive peptides. We further hypothesized that peptide-centric CARs could recognize peptides on additional HLA allotypes when presented in a similar manner. Informed by computational modelling, we showed that PHOX2B peptide-centric CARs also recognize QYNPIRTTF presented by HLA-A*23:01 and the highly divergent HLA-B*14:02. Finally, we demonstrated potent and specific killing of neuroblastoma cells expressing these HLAs in vitro and complete tumour regression in mice. These data suggest that peptide-centric CARs have the potential to vastly expand the pool of immunotherapeutic targets to include non-immunogenic intracellular oncoproteins and widen the population of patients who would benefit from such therapy by breaking conventional HLA restriction.

[In silico specificity determination of neoantigen-reactive T-lymphocytes]

Effective personalized immunotherapies of the future will need to capture not only the peculiarities of the patient's tumor but also of his immune response to it. In this study, using results of in vitro high-throughput specificity assays, and combining comparative models of pMHCs and TCRs using molecular docking, we have constructed all-atom models for the putative complexes of all their possible pairwise TCR-pMHC combinations. For the models obtained we have calculated a dataset of physics-based scores and have trained binary classifiers that perform better compared to their solely sequence-based counterparts. These structure-based classifiers pinpoint the most prominent energetic terms and structural features characterizing the type of protein-protein interactions that underlies the immune recognition of tumors by T cells.

Myositis/Myasthenia after Pembrolizumab in a Bladder Cancer Patient with an Autoimmunity-Associated HLA: Immune-Biological Evaluation and Case Report

Pembrolizumab (mAb to PD-1) has been recently approved for the therapy of pretreated urothelial cancer. Despite the efficacy, it is often accompanied by unpredictable and sometime severe immune-related (ir) adverse events (AEs). Here, we report the clinical and immune-biological characterization of a patient with a metastatic bladder cancer who developed myositis signs (M) and a myasthenia-like syndrome (MLS) during treatment with pembrolizumab. The patient presented an autoimmunity-associated HLA haplotype (HLA-A*02/HLA-B*08/HLA-C*07/HLA-DRB1*03) and experienced an increase in activated CD8 T-cells along the treatment. The symptomatology regressed after pembrolizumab discontinuation and a pyridostigmine and steroids-based therapy. This is the first report of concurrent M and MLS appearance in cancer patients receiving pembrolizumab. More efforts are needed to define early the risk and the clinical meaning of irAEs in this setting.

Bioinformatic HLA Studies in the Context of SARS-CoV-2 Pandemic and Review on Association of HLA Alleles with Preexisting Medical Conditions

After the announcement of a new coronavirus in China in December 2019, which was then called SARS-CoV-2, this virus changed to a global concern and it was then declared as a pandemic by WHO. Human leukocyte antigen (HLA) alleles, which are one of the most polymorphic genes, play a pivotal role in both resistance and vulnerability of the body against viruses and other infections as well as chronic diseases. The association between HLA alleles and preexisting medical conditions such as cardiovascular diseases and diabetes mellitus is reported in various studies. In this review, we focused on the bioinformatic HLA studies to summarize the HLA alleles which responded to SARS-CoV-2 peptides and have been used to design vaccines. We also reviewed HLA alleles that are associated with comorbidities and might be related to the high mortality rate among COVID-19 patients. Since both genes and patients' medical conditions play a key role in both severity of the disease and the mortality rate in COVID-19 patients, a better understanding of the connection between HLA alleles and SARS-CoV-2 can provide a wider perspective on the behavior of the virus. Such understanding can help scientists, especially in terms of protecting healthcare workers and designing effective vaccines.

Immunoinformatics and molecular modeling approach to design universal multi-epitope vaccine for SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmittable and pathogenic human coronavirus that caused a pandemic situation of acute respiratory syndrome, called COVID-19, which has posed a significant threat to global health security. The aim of the present study is to computationally design an effective peptide-based multi-epitope vaccine (MEV) against SARS-CoV-2. The overall model quality of the vaccine candidate, immunogenicity, allergenicity, and physiochemical analysis have been conducted and validated. Molecular dynamics studies confirmed the stability of the candidate vaccine. The docked complexes during the simulation revealed a strong and stable binding interactions of MEV with human and mice toll-like receptors (TLR), TLR3 and TLR4. Finally, candidate vaccine codons have been optimized for their in silico cloning in E. coli expression system, to confirm increased expression. The proposed MEV can be a potential candidate against SARS-CoV-2, but experimental validation is needed to ensure its safety and immunogenicity status.