CD4(+) T cell lineage integrity is controlled by the histone deacetylases HDAC1 and HDAC2.

Molecular mechanisms that maintain lineage integrity of helper T cells are largely unknown. Here we show histone deacetylases 1 and 2 (HDAC1 and HDAC2) as crucial regulators of this process. Loss of HDAC1 and HDAC2 during late T cell development led to the appearance of major histocompatibility complex (MHC) class II-selected CD4(+) helper T cells that expressed CD8-lineage genes such as Cd8a and Cd8b1. HDAC1 and HDAC2-deficient T helper type 0 (TH0) and TH1 cells further upregulated CD8-lineage genes and acquired a CD8(+) effector T cell program in a manner dependent on Runx-CBFß complexes, whereas TH2 cells repressed features of the CD8(+) lineage independently of HDAC1 and HDAC2. These results demonstrate that HDAC1 and HDAC2 maintain integrity of the CD4 lineage by repressing Runx-CBFß complexes that otherwise induce a CD8(+) effector T cell-like program in CD4(+) T cells.

Conformational flexibility of a free and TCR-bound pMHC-I protein investigated by long-term molecular dynamics simulations.

BACKGROUND: Major histocompatibility complexes (MHCs) play a crucial role in the cell-mediated adaptive immune response as they present antigenic peptides (p) which are recognized by host T cells through a complex formation of the T cell receptor (TCR) with pMHC. In the present study, we report on changes in conformational flexibility within a pMHC molecule upon TCR binding by looking at molecular dynamics (MD) simulations of the free and the TCR-bound pMHC-I protein of the LC13-HLA-B*44:05-pEEYLQAFTY complex. RESULTS: We performed long-term MD simulations with a total simulation time of 8 µs, employing 10 independent 400 ns replicas for the free and the TCR-bound pMHC system. Upon TCR ligation, we observed a reduced dynamic flexibility in the central residues of the peptide and the MHC α1-helix, altered occurrences of hydrogen bonds between the peptide and the MHC, a reduced conformational entropy of the peptide-binding groove, as well as a decreased solvent accessible surface area. CONCLUSIONS: In summary, our results from 8 µs MD simulations indicate a restricted conformational space of the MHC peptide-binding groove upon TCR ligation and suggest a minimum simulation time of approximately 100 ns for biomolecules of comparable complexity to draw meaningful conclusions. Given the relatively long total simulation time, our results contribute to a more detailed view on conformational flexibility properties of the investigated free and TCR-bound pMHC-I system.

Galectin network in osteoarthritis: galectin-4 programs a pathogenic signature of gene and effector expression in human chondrocytes in vitro.

Galectin-4 (Gal-4) is a member of the galectin family, which have been identified as galactose-binding proteins. Gal-4 possesses two tandem repeat carbohydrate recognition domains and acts as a cross-linking bridge in sulfatide-dependent glycoprotein routing. We herein document its upregulation in osteoarthritis (OA) in correlation with the extent of cartilage degradation in vivo. Primary human OA chondrocytes in vitro respond to carbohydrate-inhibitable Gal-4 binding with the upregulation of pro-degradative/-inflammatory proteins such as interleukin-1ß (IL-1ß) and matrix metalloproteinase-13 (MMP-13), as documented by RT-qPCR-based mRNA profiling and transcriptome data processing. Activation of p65 by phosphorylation of Ser536 within the NF-κB pathway and the effect of three p65 inhibitors on Gal-4 activity support downstream involvement of such signaling. In 3D (pellet) cultures, Gal-4 presence causes morphological and biochemical signs of degradation. Taken together, our findings strongly support the concept of galectins acting as a network in OA pathogenesis and suggest that blocking their activity in disease progression may become clinically relevant in the future.

Molecular dynamics of the immune checkpoint programmed cell death protein I, PD-1: conformational changes of the BC-loop upon binding of the ligand PD-L1 and the monoclonal antibody nivolumab.

BACKGROUND: The immune checkpoint receptor programmed cell death protein I (PD-1) has been identified as a key target in immunotherapy. PD-1 reduces the risk of autoimmunity by inducing apoptosis in antigen-specific T cells upon interaction with programmed cell death protein ligand I (PD-L1). Various cancer types overexpress PD-L1 to evade the immune system by inducing apoptosis in tumor-specific CD8+ T cells. The clinically used blocking antibody nivolumab binds to PD-1 and inhibits the immunosuppressive interaction with PD-L1. Even though PD-1 is already used as a drug target, the exact mechanism of the receptor is still a matter of debate. For instance, it is hypothesized that the signal transduction is based on an active conformation of PD-1. RESULTS: Here we present the results of the first molecular dynamics simulations of PD-1 with a complete extracellular domain with a focus on the role of the BC-loop of PD-1 upon binding PD-L1 or nivolumab. We could demonstrate that the BC-loop can form three conformations. Nivolumab binds to the BC-loop according to the conformational selection model whereas PD-L1 induces allosterically a conformational change of the BC-loop. CONCLUSION: Due to the structural differences of the BC-loop, a signal transduction based on active conformation cannot be ruled out. These findings will have an impact on drug design and will help to refine immunotherapy blocking antibodies.

Transforming Growth Factor-ß and Axl Induce CXCL5 and Neutrophil Recruitment in Hepatocellular Carcinoma.

Transforming growth factor (TGF)-ß suppresses early hepatocellular carcinoma (HCC) development but triggers pro-oncogenic abilities at later stages. Recent data suggest that the receptor tyrosine kinase Axl causes a TGF-ß switch toward dedifferentiation and invasion of HCC cells. Here, we analyzed two human cellular HCC models with opposing phenotypes in response to TGF-ß. Both HCC models showed reduced proliferation and clonogenic growth behavior following TGF-ß stimulation, although they exhibited differences in chemosensitivity and migratory abilities, suggesting that HCC cells evade traits of anti-oncogenic TGF-ß. Transcriptome profiling revealed differential regulation of the chemokine CXCL5, which positively correlated with TGF-ß expression in HCC patients. The expression and secretion of CXCL5 was dependent on Axl expression, suggesting that CXCL5 is a TGF-ß target gene collaborating with Axl signaling. Loss of either TGF-ß or Axl signaling abrogated CXCL5-dependent attraction of neutrophils. In mice, tumor formation of transplanted HCC cells relied on CXCL5 expression. In HCC patients, high levels of Axl and CXCL5 correlated with advanced tumor stages, recruitment of neutrophils into HCC tissue, and reduced survival. Conclusion: The synergy of TGF-ß and Axl induces CXCL5 secretion, causing the infiltration of neutrophils into HCC tissue. Intervention with TGF-ß/Axl/CXCL5 signaling may be an effective therapeutic strategy to combat HCC progression in TGF-ß-positive patients.

Galectin-8 induces functional disease markers in human osteoarthritis and cooperates with galectins-1 and -3.

The reading of glycan-encoded signals by tissue lectins is considered a major route of the flow of biological information in many (patho)physiological processes. The arising challenge for current research is to proceed from work on a distinct protein to family-wide testing of lectin function. Having previously identified homodimeric galectin-1 and chimera-type galectin-3 as molecular switches in osteoarthritis progression, we here provide proof-of-principle evidence for an intra-network cooperation of galectins with three types of modular architecture. We show that the presence of tandem-repeat-type galectin-8 significantly correlated with cartilage degeneration and that it is secreted by osteoarthritic chondrocytes. Glycan-inhibitable surface binding of galectin-8 to these cells increased gene transcription and the secretion of functional disease markers. The natural variant galectin-8 (F19Y) was less active than the prevalent form. Genome-wide array analysis revealed induction of a pro-degradative/inflammatory gene signature, largely under control of NF-κB signaling. This signature overlapped with respective gene-expression patterns elicited by galectins-1 and -3, but also presented supplementary features. Functional assays with mixtures of galectins that mimic the pathophysiological status unveiled cooperation between the three galectins. Our findings shape the novel concept to consider individual galectins as part of a so far not realized teamwork in osteoarthritis pathogenesis, with relevance beyond this disease.

Co-expressed genes enhance precision of receptor status identification in breast cancer patients.

PURPOSE: Therapeutic decisions in breast cancer patients crucially depend on the status of estrogen receptor, progesterone receptor and HER2, obtained by immunohistochemistry (IHC). These are known to be inaccurate sometimes, and we demonstrate how to use gene-expression to increase precision of receptor status. METHODS: We downloaded data from 3241 breast cancer patients out of 36 clinical studies. For each receptor, we modelled the mRNA expression of the receptor gene and a co-gene by logistic regression. For each patient, predictions from logistic regression were merged with information from IHC on a probabilistic basis to arrive at a fused prediction result. RESULTS: We introduce Sankey diagrams to visualize the step by step increase of precision as information is added from gene expression: IHC-estimates are qualified as 'confirmed', 'rejected' or 'corrected'. Additionally, we introduce the category 'inconclusive' to spot those patients in need for additional assessments so as to increase diagnostic precision and safety. CONCLUSIONS: We demonstrate a sound mathematical basis for the fusion of information, even if partly contradictive. The concept is extendable to more than three sources of information, as particularly important for OMICS data. The overall number of undecidable cases is reduced as well as those assessed falsely. We outline how decision rules may be extended to also weigh consequences, being different in severity for false-positive and false-negative assessments, respectively. The possible benefit is demonstrated by comparing the disease free survival between patients whose IHC could be confirmed versus those for which it was corrected.

Maternal embryonic leucine zipper kinase is a novel target for proliferation-associated high-risk myeloma.

Treatment of high-risk patients is a major challenge in multiple myeloma. This is especially true for patients assigned to the gene expression profiling-defined proliferation subgroup. Although recent efforts have identified some key players of proliferative myeloma, genetic interactions and players that can be targeted with clinically effective drugs have to be identified in order to overcome the poor prognosis of these patients. We therefore examined maternal embryonic leucine zipper kinase (MELK) for its implications in hyper-proliferative myeloma and analyzed the activity of the MELK inhibitor OTSSP167 both in vitro and in vivoMELK was found to be significantly overexpressed in the proliferative subgroup of myeloma. This finding translated into poor overall survival in patients with high vs low MELK expression. Enrichment analysis of upregulated genes in myeloma cells of MELKhigh patients confirmed the strong implications in myeloma cell proliferation. Targeting MELK with OTSSP167 impaired the growth and survival of myeloma cells, thereby affecting central survival factors such as MCL-1 and IRF4 This activity was also observed in the 5TGM.1 murine model of myeloma. OTSSP167 reduced bone marrow infiltration and serum paraprotein levels in a dose-dependent manner. In addition, we revealed a strong link between MELK and other proliferation-associated high-risk genes (PLK-1, EZH2, FOXM1, DEPDC1) and MELK inhibition also impaired the expression of those genes. We therefore conclude that MELK is an essential component of a proliferative gene signature and that pharmacological inhibition of MELK represents an attractive novel approach to overcome the poor prognosis of high-risk patients with a proliferative expression pattern.

Specificities of human CD4+ T cell responses to an inactivated flavivirus vaccine and infection: correlation with structure and epitope prediction.

Tick-borne encephalitis (TBE) virus is endemic in large parts of Europe and Central and Eastern Asia and causes more than 10,000 annual cases of neurological disease in humans. It is closely related to the mosquito-borne yellow fever, dengue, Japanese encephalitis, and West Nile viruses, and vaccination with an inactivated whole-virus vaccine can effectively prevent clinical disease. Neutralizing antibodies are directed to the viral envelope protein (E) and an accepted correlate of immunity. However, data on the specificities of CD4(+) T cells that recognize epitopes in the viral structural proteins and thus can provide direct help to the B cells producing E-specific antibodies are lacking. We therefore conducted a study on the CD4(+) T cell response against the virion proteins in vaccinated people in comparison to TBE patients. The data obtained with overlapping peptides in interleukin-2 (IL-2) enzyme-linked immunosorbent spot (ELISpot) assays were analyzed in relation to the three-dimensional structures of the capsid (C) and E proteins as well as to epitope predictions based on major histocompatibility complex (MHC) class II peptide affinities. In the C protein, peptides corresponding to two out of four alpha helices dominated the response in both vaccinees and patients, whereas in the E protein concordance of immunodominance was restricted to peptides of a single domain (domain III). Epitope predictions were much better for C than for E and were especially erroneous for the transmembrane regions. Our data provide evidence for a strong impact of protein structural features that influence peptide processing, contributing to the discrepancies observed between experimentally determined and computer-predicted CD4(+) T cell epitopes. Importance: Tick-borne encephalitis virus is endemic in large parts of Europe and Asia and causes more than 10,000 annual cases of neurological disease in humans. It is closely related to yellow fever, dengue, Japanese encephalitis, and West Nile viruses, and vaccination with an inactivated vaccine can effectively prevent disease. Both vaccination and natural infection induce the formation of antibodies to a viral surface protein that neutralize the infectivity of the virus and mediate protection. B lymphocytes synthesizing these antibodies require help from other lymphocytes (helper T cells) which recognize small peptides derived from proteins contained in the viral particle. Which of these peptides dominate immune responses to vaccination and infection, however, was unknown. In our study we demonstrate which parts of the proteins contribute most strongly to the helper T cell response, highlight specific weaknesses of currently available approaches for their prediction, and demonstrate similarities and differences between vaccination and infection.

Differential geometric analysis of alterations in MH α-helices.

Antigen presenting cells present processed peptides via their major histocompatibility (MH) complex to the T cell receptors (TRs) of T cells. If a peptide is immunogenic, a signaling cascade can be triggered within the T cell. However, the binding of different peptides and/or different TRs to MH is also known to influence the spatial arrangement of the MH α-helices which could itself be an additional level of T cell regulation. In this study, we introduce a new methodology based on differential geometric parameters to describe MH deformations in a detailed and comparable way. For this purpose, we represent MH α-helices by curves. On the basis of these curves, we calculate in a first step the curvature and torsion to describe each α-helix independently. In a second step, we calculate the distribution parameter and the conical curvature of the ruled surface to describe the relative orientation of the two α-helices. On the basis of four different test sets, we show how these differential geometric parameters can be used to describe changes in the spatial arrangement of the MH α-helices for different biological challenges. In the first test set, we illustrate on the basis of all available crystal structures for (TR)/pMH complexes how the binding of TRs influences the MH helices. In the second test set, we show a cross evaluation of different MH alleles with the same peptide and the same MH allele with different peptides. In the third test set, we present the spatial effects of different TRs on the same peptide/MH complex. In the fourth test set, we illustrate how a severe conformational change in an α-helix can be described quantitatively. Taken together, we provide a novel structural methodology to numerically describe subtle and severe alterations in MH α-helices for a broad range of applications.

Flexible Risk Evidence Combination Rules in Breast Cancer Precision Therapy.

Evidence theory by Dempster-Shafer for determination of hormone receptor status in breast cancer samples was introduced in our previous paper. One major topic pointed out here is the link between pieces of evidence found from different origins. In this paper the challenge of selecting appropriate ways of fusing evidence, depending on the type and quality of data involved is addressed. A parameterized family of evidence combination rules, covering the full range of potential needs, from emphasizing discrepancies in the measurements to aspiring accordance, is covered. The consequences for real patient samples are shown by modeling different decision strategies.

Clustering molecular dynamics conformations of the CC'-loop of the PD-1 immuno-checkpoint receptor.

Molecular mechanisms within the checkpoint receptor PD-1 are essential for its activation by PD-L1 as well as for blocking such an activation via checkpoint inhibitors. We use molecular dynamics to scrutinize patterns of atomic motion in PD-1 without a ligand. Molecular dynamics is performed for the whole extracellular domain of PD-1, and the analysis focuses on its CC'-loop and some adjacent Cα-atoms. We extend previous work by applying common nearest neighbor clustering (Cnn) and compare the performance of this method with Daura clustering as well as UMAP dimension reduction and subsequent agglomerative linkage clustering. As compared to Daura clustering, we found Cnn less sensitive to cutoff selection and better able to return representative clusters for sets of different 3D atomic conformations. Interestingly, Cnn yields results quite similar to UMAP plus linkage clustering.

Association of HLA-DR1 with the allergic response to the major mugwort pollen allergen: molecular background.

BACKGROUND: Mugwort pollen allergens represent the main cause of pollinosis in late summer. The major allergen, Art v 1, contains only one single immunodominant, solely HLA-DR-restricted T cell epitope (Art v 125-36). The frequency of HLA-DRB1*01 is highly increased in mugwort-allergic individuals and HLA-DR1 serves as restriction element for Art v 125-36. However, Art v 125-36 also binds to HLA-DR4 with high affinity and DR1-restricted Art v 125-36 -specific T cell receptors can be activated by HLA-DR4 molecules. To understand the predominance of HLA-DR1 in mugwort allergy in spite of the degeneracy in HLA/peptide-binding and TCR-recognition, we investigated the molecular background of Art v 125-36 /MHC/TCR interactions in the context of HLA-DR1 compared to -DR4. RESULTS: The majority of Art v 125-36 -specific T cell lines and clones from HLA-DR1 carrying, mugwort pollen-allergic donors reacted to synthetic and naturally processed Art v 1-peptides when presented by HLA-DR1 or HLA-DR4 expressing antigen presenting cells. However, at limiting peptide concentrations DR1 was more effective in T cell stimulation. In addition, the minimal epitope for 50% of Art v 125-36 -specific T cells was shorter for DR1 than for DR4. In vitro binding assays of Art v 125-36 mutant peptides to isolated DR1- and DR4-molecules indicated similar binding capacities and use of the same register. In silico simulation of Art v 125-36 binding to HLA-DR1 and -DR4 suggested similar binding of the central part of the peptide to either molecule, but a higher flexibility of the N- and C-terminal amino acids and detachment at the C-terminus in HLA-DR1. CONCLUSIONS: The predominance of HLA-DR1 in the response to Art v 125-36 may be explained by subtle conformation changes of the peptide bound to DR1 compared to DR4. Computer simulation supported our experimental data by demonstrating differences in peptide mobility within the HLA-DR complex that may influence TCR-binding. We suggest that the minor differences observed in vitro may be more relevant in the microenvironment in vivo, so that only presentation by HLA-DR1, but not -DR4 permits successful T cell activation.

Decision Theory versus Conventional Statistics for Personalized Therapy of Breast Cancer.

Estrogen and progesterone receptors being present or not represents one of the most important biomarkers for therapy selection in breast cancer patients. Conventional measurement by immunohistochemistry (IHC) involves errors, and numerous attempts have been made to increase precision by additional information from gene expression. This raises the question of how to fuse information, in particular, if there is disagreement. It is the primary domain of Dempster-Shafer decision theory (DST) to deal with contradicting evidence on the same item (here: receptor status), obtained through different techniques. DST is widely used in technical settings, such as self-driving cars and aviation, and is also promising to deliver significant advantages in medicine. Using data from breast cancer patients already presented in previous work, we focus on comparing DST with classical statistics in this work, to pave the way for its application in medicine. First, we explain how DST not only considers probabilities (a single number per sample), but also incorporates uncertainty in a concept of 'evidence' (two numbers per sample). This allows for very powerful displays of patient data in so-called ternary plots, a novel and crucial advantage for medical interpretation. Results are obtained according to conventional statistics (ODDS) and, in parallel, according to DST. Agreement and differences are evaluated, and the particular merits of DST discussed. The presented application demonstrates how decision theory introduces new levels of confidence in diagnoses derived from medical data.

Molecular dynamics identifies semi-rigid domains in the PD-1 checkpoint receptor bound to its natural ligand PD-L1.

Cells in danger of being erroneously attacked by leucocytes express PD-L1 on their surface. These cells activate PD-1 on attacking leucocytes and send them to death, thus curbing erroneous, autoimmune attack. Unfortunately, cancer cells exploit this mechanism: By expressing PD-L1, they guard themselves against leucocyte attack and thereby evade immune clearance. Checkpoint inhibitors are drugs which re-enable immune clearance of cancer cells by blocking the binding of PD-L1 to PD-1 receptors. It is therefore of utmost interest to investigate these binding mechanisms. We use three 600 ns all-atom molecular dynamics simulations to scrutinize molecular motions of PD-1 with its binding partner, the natural ligand PD-L1. Usually, atomic motion patterns are evaluated against whole molecules as a reference, disregarding that such a reference is a dynamic entity by itself, thus degrading stability of the reference. As a remedy, we identify semi-rigid domains, lending themselves as more stable and reliable reference frames against which even minute differences in molecular motion can be quantified precisely. We propose an unsupervised three-step procedure. In previous work of our group and others, minute differences in motion patterns proved decisive for differences in function. Here, several highly reliable frames of reference are established for future investigations based on molecular motion.

SLCO4A1 expression is associated with activated inflammatory pathways in high-grade serous ovarian cancer.

Patients with high-grade serous ovarian cancer (HGSOC) have a very poor overall survival. Current therapeutic approaches do not bring benefit to all patients. Although genetic alterations and molecular mechanisms are well characterized, the molecular pathological conditions are poorly investigated. Solute carrier organic anion transporter family member 4A1 (SLCO4A1) encodes OATP4A1, which is an uptake membrane transporter of metabolic products. Its expression may influence various signaling pathways associated with the molecular pathophysiological conditions of HGSOC and consequently tumor progression. RNA sequencing of 33 patient-derived HGSOC cell lines showed that SLCO4A1 expression was diverse by individual tumors, which was further confirmed by RT-qPCR, Western blotting and immunohistochemistry. Gene Set Enrichment Analysis revealed that higher SLCO4A1 level was associated with inflammation-associated pathways including NOD-like receptor, adipocytokine, TALL1, CD40, NF-κB, and TNF-receptor 2 signaling cascades, while low SLCO4A1 expression was associated with the mitochondrial electron transport chain pathway. The overall gene expression pattern in all cell lines was specific to each patient and remained largely unchanged during tumor progression. In addition, genes encoding ABCC3 along with SLCO4A1-antisense RNA 1, were associated with higher expression of the SLCO4A1, indicating their possible involvement in inflammation-associated pathways that are downstream to the prostaglandin E2/cAMP axis. Taken together, increased SLCO4A1/OATP4A1 expression is associated with the upregulation of specific inflammatory pathways, while the decreased level is associated with mitochondrial dysfunction. These molecular pathophysiological conditions are tumor specific and should be taken into consideration by the development of therapies against HGSOC.

PeptX: using genetic algorithms to optimize peptides for MHC binding.

BACKGROUND: The binding between the major histocompatibility complex and the presented peptide is an indispensable prerequisite for the adaptive immune response. There is a plethora of different in silico techniques for the prediction of the peptide binding affinity to major histocompatibility complexes. Most studies screen a set of peptides for promising candidates to predict possible T cell epitopes. In this study we ask the question vice versa: Which peptides do have highest binding affinities to a given major histocompatibility complex according to certain in silico scoring functions? RESULTS: Since a full screening of all possible peptides is not feasible in reasonable runtime, we introduce a heuristic approach. We developed a framework for Genetic Algorithms to optimize peptides for the binding to major histocompatibility complexes. In an extensive benchmark we tested various operator combinations. We found that (1) selection operators have a strong influence on the convergence of the population while recombination operators have minor influence and (2) that five different binding prediction methods lead to five different sets of "optimal" peptides for the same major histocompatibility complex. The consensus peptides were experimentally verified as high affinity binders. CONCLUSION: We provide a generalized framework to calculate sets of high affinity binders based on different previously published scoring functions in reasonable runtime. Furthermore we give insight into the different behaviours of operators and scoring functions of the Genetic Algorithm.

Decision theory for precision therapy of breast cancer.

Correctly estimating the hormone receptor status for estrogen (ER) and progesterone (PGR) is crucial for precision therapy of breast cancer. It is known that conventional diagnostics (immunohistochemistry, IHC) yields a significant rate of wrongly diagnosed receptor status. Here we demonstrate how Dempster Shafer decision Theory (DST) enhances diagnostic precision by adding information from gene expression. We downloaded data of 3753 breast cancer patients from Gene Expression Omnibus. Information from IHC and gene expression was fused according to DST, and the clinical criterion for receptor positivity was re-modelled along DST. Receptor status predicted according to DST was compared with conventional assessment via IHC and gene-expression, and deviations were flagged as questionable. The survival of questionable cases turned out significantly worse (Kaplan Meier p < 1%) than for patients with receptor status confirmed by DST, indicating a substantial enhancement of diagnostic precision via DST. This study is not only relevant for precision medicine but also paves the way for introducing decision theory into OMICS data science.

vmdICE: a plug-in for rapid evaluation of molecular dynamics simulations using VMD.

Molecular dynamics (MD) is a powerful in silico method to investigate the interactions between biomolecules. It solves Newton's equations of motion for atoms over a specified period of time and yields a trajectory file, containing the different spatial arrangements of atoms during the simulation. The movements and energies of each single atom are recorded. For evaluating of these simulation trajectories with regard to biomedical implications, several methods are available. Three well-known ones are the root mean square deviation (RMSD), the root mean square fluctuation (RMSF) and solvent accessible surface area (SASA). Herein, we present a novel plug-in for the software "visual molecular dynamics" (VMD) that allows an interactive 3D representation of RMSD, RMSF, and SASA, directly on the molecule. On the one hand, our plug-in is easy to handle for inexperienced users, and on the other hand, it provides a fast and flexible graphical impression of the spatial dynamics of a system for experts in the field.

Pembrolizumab Induces an Unexpected Conformational Change in the CC'-loop of PD-1.

To improve cancer immunotherapy, a clearer understanding of key targets such as the immune checkpoint receptor PD-1 is essential. The PD-1 inhibitors nivolumab and pembrolizumab were recently approved by the FDA. The CC'-loop of PD-1 has been identified as a hotspot for drug targeting. Here, we investigate the influence of nivolumab and pembrolizumab on the molecular motion of the CC'-loop of PD-1. We performed molecular dynamics simulations on the complete extracellular domain of PD-1, in complex with PD-L1, and the blocking antibodies nivolumab and pembrolizumab. Conformations of the CC'-loop were analyzed unsupervised with the Daura et al. clustering algorithm and multidimensional scaling. Surprisingly, two conformations found were seen to correspond to the 'open' and 'closed' conformation of CC'-loop in apo-PD-1, already known from literature. Unsupervised clustering also surprisingly reproduced the natural ligand, PD-L1, exclusively stabilizing the 'closed' conformation, as also known from literature. Nivolumab, like PD-L1, was found to shift the equilibrium towards the 'closed' conformation, in accordance with the conformational selection model. Pembrolizumab, on the other hand, induced a third conformation of the CC'-loop which has not been described to date: Relative to the conformation 'open' the, CC'-loop turned 180° to form a new conformation which we called 'overturned'. We show that the combination of clustering and multidimensional scaling is a fast, easy, and powerful method in analyzing structural changes in proteins. Possible refined antibodies or new small molecular compounds could utilize the flexibility of the CC'-loop to improve immunotherapy.