The major defense system against microbial pathogens in vertebrates is the adaptive immune response and represents an effective mechanism in cancer surveillance. T cells represent an essential component of this complex system. They can recognize myriads of antigens as short peptides (p) originated from the intracellular degradation of foreign proteins presented by major histocompatibility complex (MHC) proteins. The clonotypic T-cell antigen receptor (TCR) is specialized in recognizing pMHC and triggering T cells immune response. It is still unclear how TCR engagement to pMHC is translated into the intracellular signal that initiates T-cell immune response. Some work has suggested the possibility that pMHC binding induces in the TCR conformational changes transmitted to its companion CD3 subunits that govern signaling. The conformational changes would promote phosphorylation of the CD3 complex ζ chain that initiates signal propagation intracellularly. Here, we used all-atom molecular dynamics simulations (MDs) of 500 ns to analyze the conformational behavior of three TCRs (1G4, ILA1 and ILA1α1β1) interacting with the same MHC class I (HLA-A*02:01) bound to different peptides, and modelled in the presence of a lipid bilayer. Our data suggest a correlation between the conformations explored by the β-chain constant regions and the T-cell response experimentally determined. In particular, independently by the TCR type involved in the interaction, the TCR activation seems to be linked to a specific zone of the conformational space explored by the β-chain constant region. Moreover, TCR ligation restricts the conformational space the MHC class I groove.

In this work, we describe the application of the Zernike formalism to quantitatively characterize the binding pockets of two sets of biologically relevant systems. Such an approach, when applied to molecular dynamics trajectories, is able to pinpoint the subtle differences between very similar molecular regions and their impact on the local propensity to ligand binding, allowing us to quantify such differences. The statistical robustness of our procedure suggests that it is very suitable to describe protein binding sites and protein–ligand interactions within a rigorous and well-defined framework.

The human leukocyte antigen HLA-B27 is a strong risk factor for Ankylosing Spondylitis (AS), an immune-mediated disorder affecting axial skeleton and sacroiliac joints. Additionally, evidence exists sustaining a strong protective role for HLA-B27 in viral infections. These two aspects could stem from common molecular mechanisms. Recently, we have found that the HLA-B*2705 presents an EBV epitope (pEBNA3A-RPPIFIRRL), lacking the canonical B27 binding motif but known as immunodominant in the HLA-B7 context of presentation. Notably, 69% of B*2705 carriers, mostly patients with AS, possess B*2705-restricted, pEBNA3A-specific CD8+ T cells. Contrarily, the non-AS-associated B*2709 allele, distinguished from the B*2705 by the single His116Asp polymorphism, is unable to display this peptide and, accordingly, B*2709 healthy subjects do not unleash specific T cell responses. Herein, we investigated whether the reactivity towards pEBNA3A could be a side effect of the recognition of the natural longer peptide (pKEBNA3A) having the classical B27 consensus (KRPPIFIRRL). The stimulation of PBMC from B*2705 positive patients with AS in parallel with both pEBNA3A and pKEBNA3A did not allow to reach an unambiguous conclusion since the differences in the magnitude of the response measured as percentage of IFNγ-producing CD8+ T cells were not statistically significant. Interestingly, computational analysis suggested a structural shift of pEBNA3A as well as of pKEBNA3A into the B27 grooves, leaving the A pocket partially unfilled. To our knowledge this is the first report of a viral peptide: HLA-B27 complex recognized by TCRs in spite of a partially empty groove. This implies a rethinking of the actual B27 immunopeptidome crucial for viral immune-surveillance and autoimmunity.

Graphical abstract ABSTRACTThe self-association equilibria of doxorubicin hydrochloride (DX), at high drug and NaCl concentrations, are studied by temperature scan fluorescence spectroscopy, with the support of molecular dynamics (MD) calculations. Even though all anthracyclines show dimerization equilibria, DX only can further associate into long polymeric chains according to DXmon ⇄ DXdim ⇄ DXpol. This is reflected not only in the mechanical properties of DXpol solutions (behaving as thixotropic gels) but also in their spectroscopic behaviour. Fluorescence, in particular, is the technique of election to study this complex set of equilibria. Upon increasing the temperature, DXpol melts into DXdim, which in turn is in equilibrium with DXmon. Since DXdim is non fluorescent, with a fluorescence temperature scan experiment the DXpol⇄ DXmon equilibrium is probed. However, also information on the DX dimerization equilibrium can be derived together with the relevant thermodynamic parameters ruling the dimerization process (ΔH dim°= −56 kJ mol −1 ; ΔS dim°= −97 J mol −1 K −1 ). The residence time of DX molecules in the dimer (74.7 μs), as well as the monomers mutual orientation in the dimer, are characterized by means of theoretical and computational modelling.

Human hemoglobin (HbA) is one of the prototypal systems used to investigate structure–function relationships in proteins. Indeed, HbA has been used to develop the basic concepts of protein allostery, although the atomic-level mechanism underlying the HbA functionality is still highly debated. This is due to the fact that most of the three-dimensional structural information collected over the decades refers to the endpoints of HbA functional transition with little data available for the intermediate states. Here, we report molecular dynamics (MD) simulations by focusing on the relevance of the intermediate states of the protein functional transition unraveled by the crystallographic studies carried out on vertebrate Hbs. Fully atomistic simulations of the HbA T-state indicate that the protein undergoes a spontaneous transition toward the R-state. The inspection of the trajectory structures indicates that the protein significantly populates the intermediate HL-(C) state previously unraveled by crystallography. In the structural transition, it also assumes the intermediate states crystallographically detected in Antarctic fish Hbs. This finding suggests that HbA and Antarctic fish Hbs, in addition to the endpoints of the transitions, also share a similar deoxygenation pathway despite a distace of hundreds of millions of years in the evolution scale. Finally, using the essential dynamic sampling methodology, we gained some insights into the reverse R to T transition that is not spontaneously observed in classic MD simulations.