Molecular dynamics simulations to provide insights into epitopes coupled to the soluble and membrane-bound MHC-II complexes

Reverse vaccinology approach for identification of epitopes from E1 protein as peptide vaccine against HCV: A proof of concept

The development of effective vaccines against Hepatitis C Virus (HCV) remains a global health priority and challenge. In this study, we employed an integrative approach combining computational epitope prediction with experimental validation to identify immunogenic peptides targeting the E1 glycoprotein of HCV. In the present report, computational data from various epitope prediction algorithms such as IEDB and SYFPEITHI, followed by molecular dynamics (MD) simulations and immuno-informatics analysis is presented. Through computational screening, we identified potential epitope candidates, with QVRNSSGLY (P3) and QLFTFSPRH (P7) emerging as promising candidates. MD simulations revealed stable interactions between these epitopes and MHC molecule, further validated by free energy estimations using MMPBSA method. Immuno-informatics analysis supported these findings, showing high binding potential and immunogenicity scores for the selected peptides. Subsequent synthesis and characterization of epitope peptides confirmed their structural integrity and purity required for conducting immune activation assays. Experimental immunological assays carried out in this study involved epitope peptide induced activation of CD8 + and CD4 + T cells from healthy human subjects and HCV- recovered patients. Data from experimental validation revealed significant cytokine release upon exposure to epitope peptides, particularly TNF-a, IL-6, and GM-CSF, indicative of robust immune responses. Notably, peptides P3 and P7 exhibited the most pronounced cytokine induction profiles, underscoring their potential as vaccine candidates. Further investigations addressing the mechanism of action of these epitope peptides under preclinical and clinical settings may help in developing effective vaccine against HCV.

A review on medicinally important heterocyclic compounds and importance of biophysical approach of underlying the insight mechanism in biological environment

Heterocyclic derivatives have more interesting biological properties which hold a remarkable place in pharmaceutical industries due to their unique physiochemical properties and ease of adaption in various biological environments. Of many, the above-said derivatives have been recently examined for their promising action against a few malignancies. Specifically, anti-cancer research has benefited from these derivatives' natural flexibility and dynamic core scaffold. In any case, concerning some other promising anti-cancer drugs, heterocyclic derivative doesn't come without deficiencies. To be a successful drug candidate it should poses Absorption, Distribution, Metabolism and Eliminations (ADME) parameter, and must also have good binding interaction towards carrier protein as well as DNA and less in toxic nature, economically feasible. In this review, we described the overview of biologically important heterocyclic derivatives and their main application in medicine. Further, we focus types of biophysical techniques to understand the binding interaction mechanism.Communicated by Ramaswamy H. Sarma.

Investigations on the binding properties of hydroxylated polybrominated diphenyl ethers with lysozyme using the multispectral techniques and molecular modeling

As a kind of phenolic chemical with endocrine disrupting potency, hydroxylated polybrominated diphenyl ethers (OH-PBDEs) cause a latent threat to human health from their residue in the environment. Their binding efficiency with lysozyme (LYSO) was studied by molecular simulation combined with fluorescence, UV-vis absorption and circular dichroism (CD), so as to assess their toxicity at the molecular level. Molecular docking data indicate that van der Waals force is the principal interaction force between OH-PBDEs and LYSO. The binding site for 5'-OH-BDE-25 in LYSO is ascertained as the active site, which interaction with the TRP63 and TRP108 residues of LYSO to take shape a strong face-to-face stacked rank (F-shaped). Both 4'-OH-BDE-99 and 3'-OH-BDE-154 display a certain degree of deviation from the active site. Nevertheless, their F-shaped interaction with TRP63 conduce to bind LYSO and stabilize the docking conformation. Combined with dynamics simulation and spectral analysis, the secondary structure of LYSO can be induced by the three kinds of OH-PBDEs. CD spectrum shows that the combination of LYSO and OH-PBDEs will make α- Helix content increased. The combination of OH-PBDEs and LYSO touch upon a static quenching mechanism as a result of steady state fluorescence. The energy decomposition analysis exhibited that LYSO-OH-PBDEs binding site was stable by van der Waals and hydrophobic interaction. As enzyme activity experiments demonstrate that OH-PBDEs can inhibit the activity of LYSO, which is helpful to clarify the molecular toxicity mechanism of OH-PBDEs.

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.

Exploring the toxic effects and mechanism of methoxylated polybrominated diphenyl ethers (MeO-PBDEs) on thyroxine-binding globulin (TBG): Synergy between spectroscopic and computations

The interaction mechanism between thyroxine-binding globulin (TBG) and three methoxylated polybrominated diphenyl ethers (MeO-PBDEs) was analyzed by steady-state fluorescence, ultraviolet-visible (UV-visible) spectroscopy, circular dichroism (CD), molecular docking and molecular dynamics simulation methods. The results of the molecular docking technique revealed that 2'-MeO-BDE-3, 5-MeO-BDE-47, and 3-MeO-BDE-100 combined with TBG at the active site. The steady-state fluorescence spectra displayed that MeO-PBDEs quenched the endogenous fluorescence of TBG through static quenching mechanism, and complex formation between MeO-PBDEs and TBG was further indicated by UV-vis spectroscopy. The thermodynamic quantities showed that the binding process is spontaneous, and the major forces responsible for the binding are hydrogen bonding and hydrophobic interactions, which are consistent with the results of molecular docking to a certain extent. The results of CD confirmed that the secondary structure of TBG was changed after combining with MeO-PBDEs. The dynamic simulation results illustrated that the protein structure is more compact and changes in the secondary structure of TBG after binding to MeO-PBDEs. Additionally, we also utilized the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) method to analyze the binding free energy of TBG and MeO-PBDEs. The results suggest that van der Waals force plays an essential role in the combination.

Mapping and IgE-binding capacity analysis of heat/digested stable epitopes of mud crab allergens

Mud crab (Scylla paramamosain) is widely consumed after thermal processing. It is necessary to comprehensively evaluate of the allergenic potential and epitopes of allergens in high temperature-pressure (HTP) treated S. paramamosain. Tropomyosin and arginine kinase presented higher prevalence (30.77% and 42.13%) than the other three important crab allergens by component-resolved diagnosis. The surface expression of basophils CD63 and CD203c were decreased in HTP treated crab, an effect that was even more evident after digestion and absorption by the intestinal Caco-2 cell model. Of the 35 stable epitope, six were for the first time identified in shellfish. Seven heat/digested stable peptides of tropomyosin retained IgE-binding capacity and were shown to interact with MHC-II. Five epitopes (amino acids 19-29, 99-109, 153-162, 170-188 and 211-221) were the first identified in crab. The study provides insight into prevention and therapy of crab allergy, as well as helps to reduce crab allergenicity during thermal processing.

Immunoinformatics and molecular dynamics studies to predict T-cell-specific epitopes of four Klebsiella pneumoniae fimbriae antigens

Klebsiella pneumoniae (K. pneumoniae) is a causative agent of severe infections in humans. There is no publically available vaccine for K. pneumoniae infections yet. Here, using comprehensive immunoinformatics methods, T-cell-specific epitopes of four type 1 fimbriae antigens of K. pneumoniae were predicted and evaluated as potential vaccine candidates. Both CD8+ (class I) and CD4+ (class II) T-cell-specific epitopes were predicted and the epitopes similar to human proteome were excluded. Subsequently, the windows of class-II epitopes containing class-I epitopes were determined. The immunogenicity, IFN-γ production and population coverage were also estimated. Using the 3D structure of HLA and epitopes, molecular docking was carried out. Two best epitopes were selected for molecular dynamics studies. Our prediction and analyses resulted in the several dominant epitopes for each antigen. The docking results showed that all selected epitopes can bind to their restricted HLA molecules with high affinity. The molecular dynamics results indicated the stability of system with minimum possible deviation, suggesting the selected epitopes can be promising candidates for stably binding to HLA molecules. Altogether, our results suggest that the selected T-cell-specific epitopes of K. pneumoniae fimbriae antigens, particularly the two epitopes confirmed by molecular dynamics, can be applied for vaccine development. However, the in vitro and in vivo studies are required to authenticate the results of the present study.Communicated by Ramaswamy H. Sarma.

Exploring the biotransformation of N-(2-hydroxyphenyl)-2-propylpentanamide (an aryl valproic acid derivative) by CYP2C11, using in silico predictions and in vitro studies

OBJECTIVES

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA), a derivative of valproic acid (VPA), has been proposed as a potential anticancer agent due to its improved antiproliferative effects in some cancer cell lines. Although there is evidence that VPA is metabolized by cytochrome P450 2C11 rat isoform, HO-AAVPA CYP-mediated metabolism has not yet been fully explored. Therefore, in this work, the biotransformation of HO-AAVPA by CYP2C11 was investigated.

METHODS

Kinetic parameters and spectral interaction between HO-AAVPA and CYP were evaluated using rat liver microsomes. The participation of CYP2C11 in metabolism of HO-AAVPA was confirmed by cimetidine (CIM) inhibition assay. Docking and molecular dynamics simulations coupled to MMGBSA methods were used in theoretical study.

KEY FINDINGS

HO-AAVPA is metabolized by CYP enzymes (K_M  = 38.94 µm), yielding a hydroxylated metabolite according to its HPLC retention time (5.4 min) and MS analysis (252.2 m/z). In addition, CIM inhibition in rat liver microsomes (K_i  = 59.23 µm) confirmed that CYP2C11 is mainly involved in HO-AAVPA metabolism. Furthermore, HO-AAVPA interacts with CYP2C11 as a type I ligand. HO-AAVPA is stabilized at the CYP2C11 ligand recognition site through a map of interactions similar to other typical CYP2C11 substrates.

CONCLUSION

Therefore, rat liver CYP2C11 isoform is able to metabolize HO-AAVPA.

Intramolecular Domain Movements of Free and Bound pMHC and TCR Proteins: A Molecular Dynamics Simulation Study

The interaction of antigenic peptides (p) and major histocompatibility complexes (pMHC) with T-cell receptors (TCR) is one of the most important steps during the immune response. Here we present a molecular dynamics simulation study of bound and unbound TCR and pMHC proteins of the LC13-HLA-B*44:05-pEEYLQAFTY complex to monitor differences in relative orientations and movements of domains between bound and unbound states of TCR-pMHC. We generated local coordinate systems for MHC α1- and MHC α2-helices and the variable T-cell receptor regions TCR V_α and TCR V_β and monitored changes in the distances and mutual orientations of these domains. In comparison to unbound states, we found decreased inter-domain movements in the simulations of bound states. Moreover, increased conformational flexibility was observed for the MHC α2-helix, the peptide, and for the complementary determining regions of the TCR in TCR-unbound states as compared to TCR-bound states.

Multi-spectroscopic and molecular dynamics simulations investigation of the binding mechanism of polybrominated diphenyl ethers to hen egg white lysozyme

Three PBDEs (BDE25, BDE47, and BDE154) were selected to investigate the interactions between PBDEs and hen egg white lysozyme (HEWL) by molecular modeling, fluorescence spectroscopy, and FT-IR spectra. The docking results showed that hydrogen bonds were formed between BDE25 and residue TRP63 and between BDE47 and TRP63 with bond lengths of 2.178 Å and 2.146 Å, respectively. The molecular dynamics simulations indicated that van der Waals forces played a predominant role in the binding of three PBDEs to HEWL. The observed fluorescence quenching can be attributed to the formation of complexes between HEWL and PBDEs, and the quenching mechanism is a static quenching. According to Förster's non-radiative energy transfer theory, the binding distances r were < 7 nm, indicating a high probability of energy transfer from HEWL to the three PBDEs. The synchronous fluorescence showed that the emission maximum wavelength of tryptophan (TRP) residues emerged a red-shift. FT-IR spectra indicated that BDE25, BDE47 and BDE154 induced the α-helix percentage of HEWL decreased from 32.70% ± 1.64% to 28.27% ± 1.41%, 27.50% ± 1.38% and 29.78% ± 1.49%, respectively, whereas the percentage of random coil increased from 26.67% ± 1.33% to 27.60% ± 1.38%, 29.18% ± 1.46% and 30.59% ± 1.53%, respectively.

Comparative Binding Analysis of N-Acetylneuraminic Acid in Bovine Serum Albumin and Human α-1 Acid Glycoprotein

The present study focuses on the determination of the biologically significant N-acetylneuraminic acid (NANA) drug binding interaction mechanism between bovine serum albumin (BSA) and human α-1 acid glycoprotein (HAG) using various optical spectroscopy and computational methods. The steady state fluorescence spectroscopy result suggests that the fluorescence intensity of BSA and HAG was quenched by NANA in a static mode of quenching. Further time-resolved emission spectroscopy measurements confirm that mode of quenching mechanism of NANA in the BSA and HAG system. The FT-IR, excitation-emission matrix and circular dichroism (CD) analysis confirms the presence of NANA in the HAG, BSA system, and fluorescence resonance energy transfer analysis shows that NANA transfers energy between the HAG and BSA system. The molecular docking result shows good binding affinity in both protein complexes, and further molecular dynamics simulations and charge distribution analysis were performed to gain more insight into the binding interaction mechanism of NANA in the HAG and BSA complex.

Dynamics Determine Signaling in a Multicomponent System Associated with Rheumatoid Arthritis

Strategies that target multiple components are usually required for treatment of diseases originating from complex biological systems. The multicomponent system consisting of the DR4 major histocompatibility complex type II molecule, the glycopeptide CII259-273 from type II collagen, and a T-cell receptor is associated with development of rheumatoid arthritis (RA). We introduced non-native amino acids and amide bond isosteres into CII259-273 and investigated the effect on binding to DR4 and the subsequent T-cell response. Molecular dynamics simulations revealed that complexes between DR4 and derivatives of CII259-273 were highly dynamic. Signaling in the overall multicomponent system was found to depend on formation of an appropriate number of dynamic intramolecular hydrogen bonds between DR4 and CII259-273, together with the positioning of the galactose moiety of CII259-273 in the DR4 binding groove. Interestingly, the system tolerated modifications at several positions in CII259-273, indicating opportunities to use analogues to increase our understanding of how rheumatoid arthritis develops and for evaluation as vaccines to treat RA.

Structure-based Methods for Binding Mode and Binding Affinity Prediction for Peptide-MHC Complexes

Understanding the mechanisms involved in the activation of an immune response is essential to many fields in human health, including vaccine development and personalized cancer immunotherapy. A central step in the activation of the adaptive immune response is the recognition, by T-cell lymphocytes, of peptides displayed by a special type of receptor known as Major Histocompatibility Complex (MHC). Considering the key role of MHC receptors in T-cell activation, the computational prediction of peptide binding to MHC has been an important goal for many immunological applications. Sequence- based methods have become the gold standard for peptide-MHC binding affinity prediction, but structure-based methods are expected to provide more general predictions (i.e., predictions applicable to all types of MHC receptors). In addition, structural modeling of peptide-MHC complexes has the potential to uncover yet unknown drivers of T-cell activation, thus allowing for the development of better and safer therapies. In this review, we discuss the use of computational methods for the structural modeling of peptide-MHC complexes (i.e., binding mode prediction) and for the structure-based prediction of binding affinity.

Major Histocompatibility Complex (MHC) Class I and MHC Class II Proteins: Conformational Plasticity in Antigen Presentation

Antigen presentation by major histocompatibility complex (MHC) proteins is essential for adaptive immunity. Prior to presentation, peptides need to be generated from proteins that are either produced by the cell's own translational machinery or that are funneled into the endo-lysosomal vesicular system. The prolonged interaction between a T cell receptor and specific pMHC complexes, after an extensive search process in secondary lymphatic organs, eventually triggers T cells to proliferate and to mount a specific cellular immune response. Once processed, the peptide repertoire presented by MHC proteins largely depends on structural features of the binding groove of each particular MHC allelic variant. Additionally, two peptide editors-tapasin for class I and HLA-DM for class II-contribute to the shaping of the presented peptidome by favoring the binding of high-affinity antigens. Although there is a vast amount of biochemical and structural information, the mechanism of the catalyzed peptide exchange for MHC class I and class II proteins still remains controversial, and it is not well understood why certain MHC allelic variants are more susceptible to peptide editing than others. Recent studies predict a high impact of protein intermediate states on MHC allele-specific peptide presentation, which implies a profound influence of MHC dynamics on the phenomenon of immunodominance and the development of autoimmune diseases. Here, we review the recent literature that describe MHC class I and II dynamics from a theoretical and experimental point of view and we highlight the similarities between MHC class I and class II dynamics despite the distinct functions they fulfill in adaptive immunity.

Major Histocompatibility Complex (MHC) Class I and MHC Class II Proteins: Conformational Plasticity in Antigen Presentation

Antigen presentation by major histocompatibility complex (MHC) proteins is essential for adaptive immunity. Prior to presentation, peptides need to be generated from proteins that are either produced by the cell’s own translational machinery or that are funneled into the endo-lysosomal vesicular system. The prolonged interaction between a T cell receptor and specific pMHC complexes, after an extensive search process in secondary lymphatic organs, eventually triggers T cells to proliferate and to mount a specific cellular immune response. Once processed, the peptide repertoire presented by MHC proteins largely depends on structural features of the binding groove of each particular MHC allelic variant. Additionally, two peptide editors—tapasin for class I and HLA-DM for class II—contribute to the shaping of the presented peptidome by favoring the binding of high-affinity antigens. Although there is a vast amount of biochemical and structural information, the mechanism of the catalyzed peptide exchange for MHC class I and class II proteins still remains controversial, and it is not well understood why certain MHC allelic variants are more susceptible to peptide editing than others. Recent studies predict a high impact of protein intermediate states on MHC allele-specific peptide presentation, which implies a profound influence of MHC dynamics on the phenomenon of immunodominance and the development of autoimmune diseases. Here, we review the recent literature that describe MHC class I and II dynamics from a theoretical and experimental point of view and we highlight the similarities between MHC class I and class II dynamics despite the distinct functions they fulfill in adaptive immunity.

MHC class II complexes sample intermediate states along the peptide exchange pathway

The presentation of peptide-MHCII complexes (pMHCIIs) for surveillance by T cells is a well-known immunological concept in vertebrates, yet the conformational dynamics of antigen exchange remain elusive. By combining NMR-detected H/D exchange with Markov modelling analysis of an aggregate of 275 microseconds molecular dynamics simulations, we reveal that a stable pMHCII spontaneously samples intermediate conformations relevant for peptide exchange. More specifically, we observe two major peptide exchange pathways: the kinetic stability of a pMHCII's ground state defines its propensity for intrinsic peptide exchange, while the population of a rare, intermediate conformation correlates with the propensity of the HLA-DM-catalysed pathway. Helix-destabilizing mutants designed based on our model shift the exchange behaviour towards the HLA-DM-catalysed pathway and further allow us to conceptualize how allelic variation can shape an individual's MHC restricted immune response.

MHCII proteins bind and present both foreign and self-antigens to potentially activate CD4+ T cells via cognate T cell receptors (TCRs) during the adaptive immune response. Here, the authors combine NMR-detected H/D exchange with Markov modelling analysis to shed light on the dynamics of MHCII peptide exchange.

Conformational changes associated with L16P and T118M mutations in the membrane-embedded PMP22 protein, consequential in CMT-1A

Peripheral myelin protein 22 (PMP22) resides in the plasma membrane and is required for myelin formation in the peripheral nervous system. Excess PMP22 mutants accumulate in the endoplasmic reticulum (ER) resulting in the inherited neuropathies of Charcot-Marie-Tooth disease. However, there was no evidence of the structure of PMP22 or how mutations affect its folding. Therefore, in this study, we combined bioinformatics and homology modeling approaches to obtain three-dimensional native and mutated PMP22 models and its anchoring to a POPC membrane, submitted to .5-μs MD simulations, to determine how the L16P and T118M mutations affect the conformational behavior of PMP22. In addition, we investigated the ability of the native and mutated species to accumulate in the ER, via interaction with RER1, by combining protein-protein docking and MD simulations, taking the conformations that were most representative of the native and mutated PMP22 systems and RER1 conformations. Principal component analysis over MD simulations revealed that L16P and T118M mutations resulted in increased structural instability compared to the native form, which is consistent with previous experimental findings of increased structural fluctuations along a loop connecting transmembrane α-helix1 and α-helix2. Docking and MD simulations coupled with the MMGBSA approach allowed the identification that the binding interface for the PMP22-RER1 complex takes place through transmembrane α-helix1 and α-helix2, with higher effective binding free energy values between the mutated PMP22 systems and RER1 than for the native PMP22, mainly through van der Waals interactions.

Understanding the molecular basis of agonist/antagonist mechanism of GPER1/GPR30 through structural and energetic analyses

The G-protein coupled receptors (GPCRs) represent the largest superfamily of membrane proteins in charge to pass the cell signaling after binding with their cognate ligands to the cell interior. In breast cancer, a GPCR named GPER1 plays a key role in the process of growth and the proliferation of cancer cells. In a previous study, theoretical methods were applied to construct a model of GPER1, which later was submitted to molecular dynamics (MD) simulations to perform a docking calculation. Based on this preceding work, it is known that GPER1 is sensitive to structural differences in its binding site. However, due to the nature of that past study, conformational changes linked to the ligand binding were not observed. Therefore, in this study, in order to explore the conformational changes coupled to the agonist/antagonist binding, MD simulations of about 0.25μs were performed for the free and bound states, summarizing 0.75μs of MD simulation in total. For the bound states, one agonist (G-1) and antagonist (G-15) were chosen since is widely known that these two molecules cause an impact on GPER1 mobility. Based on the conformational ensemble generated through MD simulations, we found that despite G-1 and G-15 being stabilized by similar map of residues, the structural differences between both ligands impact the hydrogen bond pattern not only at the GPER1 binding site but also along the seven-helix bundle, causing significant differences in the conformational mobility along the extracellular and cytoplasmic domain, and to a lesser degree in the curvatures of helix 2, helix 3 and helix 7 between the free and bound states, which is in agreement with reported literature, and might be linked to microscopic characteristics of the activated-inactivated transition. Furthermore, binding free energy calculations using the MM/GBSA method for the bound states, followed by an alanine scanning analysis allowed us to identify some important residues for the complex stabilization.

Energetic and flexibility properties captured by long molecular dynamics simulations of a membrane-embedded pMHCII-TCR complex

Although crystallographic data have provided important molecular insight into the interactions in the pMHC-TCR complex, the inherent features of this structural approach cause it to only provide a static picture of the interactions. While unbiased molecular dynamics simulations (UMDSs) have provided important information about the dynamic structural behavior of the pMHC-TCR complex, most of them have modeled the pMHC-TCR complex as soluble, when in physiological conditions, this complex is membrane bound; therefore, following this latter UMDS protocol might hamper important dynamic results. In this contribution, we performed three independent 300 ns-long UMDSs of the pMHCII-TCR complex anchored in two opposing membranes to explore the structural and energetic properties of the recognition of pMHCII by the TCR. The conformational ensemble generated through UMDSs was subjected to clustering and Cartesian principal component analyses (cPCA) to explore the dynamical behavior of the pMHCII-TCR association. Furthermore, based on the conformational population sampled through UMDSs, the effective binding free energy, per-residue free energy decomposition, and alanine scanning mutations were explored for the native pMHCII-TCR complex, as well as for 12 mutations (p1-p12MHCII-TCR) introduced in the native peptide. Clustering analyses and cPCA provide insight into the rocking motion of the TCR onto pMHCII, together with the presence of new electrostatic interactions not observed through crystallographic methods. Energetic results provide evidence of the main contributors to the pMHC-TCR complex formation as well as the key residues involved in this molecular recognition process.

Predicting peptide vaccine candidates against H1N1 influenza virus through theoretical approaches

Identification of potential epitopes that might activate the immune system has been facilitated by the employment of algorithms that use experimental data as templates. However, in order to prove the affinity and the map of interactions between the receptor (major histocompatibility complex, MHC, or T-cell receptor) and the potential epitope, further computational studies are required. Docking and molecular dynamics (MDs) simulations have been an effective source of generating structural information at molecular level in immunology. Herein, in order to provide a detailed understanding of the origins of epitope recognition and to select the best peptide candidate to develop an epitope-based vaccine, docking and MDs simulations in combination with MMGBSA free energy calculations and per-residue free energy decomposition were performed, taking as starting complexes those formed between four designed epitopes (P1-P4) from hemagglutinin (HA) of the H1N1 influenza virus and MHC-II anchored in POPC membrane. Our results revealed that the energetic contributions of individual amino acids within the pMHC-II complexes are mainly dictated by van der Waals interactions and the nonpolar part of solvation energy, whereas the electrostatic interactions corresponding to hydrogen bonds and salt bridges determine the binding specificity, being the most favorable interactions formed between p4 and MHC-II. Then, P1-P4 epitopes were synthesized and tested experimentally to compare theoretical and experimental results. Experimental results show that P4 elicited the highest strong humoral immune response to HA of the H1N1 and may induce antibodies that are cross-reactive to other influenza subtypes, suggesting that it could be a good candidate for the development of a peptide-based vaccine.

Three-dimensional structure and molecular dynamics studies of prorrenin/renin receptor: description of the active site

The recent finding of a specific receptor for prorrenin/renin (PRR) has brought new insights into the physiology of the renin-angiotensin-aldosterone system. No undoubtable role has been described for this receptor so far. Its role seems to be important in chronic illnesses such as hypertension, possibly participating in the cardiovascular remodeling process, and diabetes where participation in inflammation development has been described. It is not possible, however, to explore the PRR function using classical pharmacological approaches due to the lack of specific agonists or antagonists. Two synthetic peptides have been described to accomplish these roles, but no conclusive data have been provided. There are no X-ray crystallography studies available to describe the structure and potential sites for drug development. So, the aim of this work was to model and theoretically describe the PRR. We describe and characterize the whole receptor protein, its spatial conformation and the potential interactions of PRR with the synthetic peptides available, describing the amino acid residues responsible for these interactions. This information provides the basis for directed development of drugs, seeking to agonize or antagonize PRR activity and study its function in health and ill stages.

Large scale characterization of the LC13 TCR and HLA-B8 structural landscape in reaction to 172 altered peptide ligands: a molecular dynamics simulation study

The interplay between T cell receptors (TCRs) and peptides bound by major histocompatibility complexes (MHCs) is one of the most important interactions in the adaptive immune system. Several previous studies have computationally investigated their structural dynamics. On the basis of these simulations several structural and dynamical properties have been proposed as effectors of the immunogenicity. Here we present the results of a large scale Molecular Dynamics simulation study consisting of 100 ns simulations of 172 different complexes. These complexes consisted of all possible point mutations of the Epstein Barr Virus peptide FLRGRAYGL bound by HLA-B*08:01 and presented to the LC13 TCR. We compare the results of these 172 structural simulations with experimental immunogenicity data. We found that simulations with more immunogenic peptides and those with less immunogenic peptides are in fact highly similar and on average only minor differences in the hydrogen binding footprints, interface distances, and the relative orientation between the TCR chains are present. Thus our large scale data analysis shows that many previously suggested dynamical and structural properties of the TCR/peptide/MHC interface are unlikely to be conserved causal factors for peptide immunogenicity.