The development of an adequate immune response against pathogens is mediated by molecular interactions between different cell types. Among them, binding of antigenic peptides to the Major Histocompatibility Complex (MHC) molecule expressed on the membrane of antigen presenting cells (APCs), and their subsequent recognition by the T cell receptor have been demonstrated to be crucial for developing an adequate immune response. The present review compiles computational quantum chemistry studies about the electrostatic potential variations induced on the MHC binding region by peptide’s amino acids, carried out with the aim of describing MHC–peptide binding interactions. The global idea is that the electrostatic potential can be represented in terms of a series expansion (charge, dipole, quadrupole, hexadecapole, etc.) whose three first terms provide a good local approximation to the molecular electrostatic ‘landscape’ and to the variations induced on such landscape by targeted modifications on the residues of the antigenic peptide. Studies carried out in four MHC class II human allele molecules, which are the most representative alleles of their corresponding haplotypes, showed that each of these molecules have conserved as well as specific electrostatic characteristics, which can be correlated at a good extent with the peptide binding profiles reported experimentally for these molecules. The information provided by such characteristics would help increase our knowledge about antigen binding and presentation, and could ultimately contribute to developing a logical and rational methodology for designing chemically synthesized, multi-antigenic, subunit-based vaccines, through the application of quantum chemistry methods.

Iron-protein interactions are involved in electron transfer reactions. Alterations of these processes are present in a number of human pathologies; among them, in Friedreich's ataxia, in which a deficiency of functional frataxin, an iron-binding protein, leads to progressive neuromuscular degenerative disease. The putative iron-binding motif of acidic residues EExxED was selected from the first α-helical stretch of the frataxin protein family and grafted onto a foreign peptide scaffold corresponding to the C-terminal α-helix from E. coli thioredoxin. The resulting grafted peptide named GRAP was studied by applying experimental (circular dichroism, isothermal titration calorimetry, capillary zone electrophoresis, thermal denaturation, NMR) and computational approaches (docking, molecular dynamics simulations). Although isolated GRAP lacks a stable secondary structure in solution, when iron is added, the peptide acquires an α-helical structure. Here we have shown that the designed peptide is able to specifically bind Fe(3+) with a moderate affinity (KD = 1.9 ± 0.2 μM) and a 1 : 1 stoichiometry. Remarkably, the GRAP/Fe(3+) interaction is entropically driven (ΔH° = -1.53 ± 0.03 kcal mol(-1) and TΔS° = 6.26 kcal mol(-1)). Experiments and simulations indicate that Fe(3+) interacts with the peptide through three acidic side chains, inducing an α-helical conformation of the grafted motif. In addition, the acidic side chains involved undergo significant conformational rearrangements upon binding, as judged by the analysis of MDs. Altogether, these results contribute to an understanding of the iron-binding mechanisms in proteins and, in particular, in the case of human frataxin.

Thioredoxin is a ubiquitous small protein that catalyzes redox reactions of protein thiols. Additionally, thioredoxin from E. coli (EcTRX) is a widely-used model for structure-function studies. In a previous paper, we characterized several single-point mutants of the C-terminal helix (CTH) that alter global stability of EcTRX. However, spectroscopic signatures and enzymatic activity for some of these mutants were found essentially unaffected. A comprehensive structural characterization at the atomic level of these near-invariant mutants can provide detailed information about structural variability of EcTRX. We address this point through the determination of the crystal structures of four point-mutants, whose mutations occurs within or near the CTH, namely L94A, E101G, N106A and L107A. These structures are mostly unaffected compared with the wild-type variant. Notably, the E101G mutant presents a large region with two alternative traces for the backbone of the same chain. It represents a significant shift in backbone positions. Enzymatic activity measurements and conformational dynamics studies monitored by NMR and molecular dynamic simulations show that E101G mutation results in a small effect in the structural features of the protein. We hypothesize that these alternative conformations represent samples of the native-state ensemble of EcTRX, specifically the magnitude and location of conformational heterogeneity.