Theoretical analysis of the neuraminidase epitope of the Mexican A H1N1 influenza strain, and experimental studies on its interaction with rabbit and human hosts

Identification of peptide epitopes of the gp120 protein of HIV-1 capable of inducing cellular and humoral immunity

The Human Immunodeficiency Virus (HIV-1) causes Acquired Immunodeficiency Syndrome (AIDS) and a high percentage of deaths. Therefore, it is necessary to design vaccines against HIV-1 for the prevention of AIDS. Bioinformatic tools and theoretical algorisms allow us to understand the structural proteins of viruses to develop vaccines based on immunogenic peptides (epitopes). In this work, we identified the epitopes: P1, P2, P10, P27 and P30 from the gp120 protein of HIV-1. These peptides were administered intranasally alone or with cholera toxin (CT) to BALB/c mice. The population of CD4+, CD8+ T lymphocytes and B cells (CD19/CD138+, IgA+ and IgG+) from nasal-associated lymphoid tissue, nasal passages, cervical and inguinal nodes was determined by flow cytometry. In addition, anti-peptides IgG and IgA from serum, nasal and vaginal washings were measured by ELISA. The results show that peptides administered by i.n. can modulate the immune response of T and B lymphocyte populations, as well as IgA and IgG antibodies secretion in the different sites analyzed. In conclusion, bioinformatics tools help us to select peptides with physicochemical properties that allow the induction of the humoral and cellular responses that depend on the peptide sequence.

The Advantage of Using Immunoinformatic Tools on Vaccine Design and Development for Coronavirus

After the outbreak of SARS-CoV-2 by the end of 2019, the vaccine development strategies became a worldwide priority. Furthermore, the appearances of novel SARS-CoV-2 variants challenge researchers to develop new pharmacological or preventive strategies. However, vaccines still represent an efficient way to control the SARS-CoV-2 pandemic worldwide. This review describes the importance of bioinformatic and immunoinformatic tools (in silico) for guide vaccine design. In silico strategies permit the identification of epitopes (immunogenic peptides) which could be used as potential vaccines, as well as nonacarriers such as: vector viral based vaccines, RNA-based vaccines and dendrimers through immunoinformatics. Currently, nucleic acid and protein sequential as well structural analyses through bioinformatic tools allow us to get immunogenic epitopes which can induce immune response alone or in complex with nanocarriers. One of the advantages of in silico techniques is that they facilitate the identification of epitopes, while accelerating the process and helping to economize some stages of the development of safe vaccines.

Role of cathepsin B of Naegleria fowleri during primary amebic meningoencephalitis

Naegleria fowleri causes primary amoebic meningoencephalitis in humans and experimental animals. It has been suggested that cysteine proteases of parasites play key roles in metabolism, nutrient uptake, host tissue invasion, and immune evasion. The aim of this work was to evaluate the presence, expression, and role of cathepsin B from N. fowleri in vitro and during PAM. Rabbit-specific polyclonal antibodies against cathepsin B were obtained from rabbit immunization with a synthetic peptide obtained by bioinformatic design. In addition, a probe was designed from mRNA for N. fowleri cathepsin B. Both protein and messenger were detected in fixed trophozoites, trophozoites interacted with polymorphonuclear and histological sections of infected mice. The main cathepsin B distribution was observed in cytoplasm or membrane mainly pseudopods and food-cups while messenger was in nucleus and cytoplasm. Surprisingly, both the messenger and enzyme were observed in extracellular medium. To determine cathepsin B release, we used trophozoites supernatant recovered from nasal passages or brain of infected mice. We observed the highest release in supernatant from recovered brain amoebae, and when we analyzed molecular weight of secreted proteins by immunoblot, we found 30 and 37 kDa bands which were highly immunogenic. Finally, role of cathepsin B during N. fowleri infection was determined; we preincubated trophozoites with E-64, pHMB or antibodies with which we obtained 60%, 100%, and 60% of survival, respectively, in infected mice. These results suggest that cathepsin B plays a role during pathogenesis caused by N. fowleri mainly in adhesion and contributes to nervous tissue damage.

Structural and molecular analyses of functional epitopes and escape mutants in Japanese encephalitis virus envelope protein domain III

The Japanese encephalitis virus (JEV) is one of the vector borne causes of encephalitis found in southeastern Asia. This positive single-stranded RNA virus is a member of the Flaviviridae family, which notably includes dengue, tick-borne, West Nile, Zika as well as yellow fever, and transmits to humans by infected mosquitos. The main site of interactions for antibodies against this virus is the envelope protein domain III (ED3). The present report investigates the time-dependent structural and conformational changes of JEV ED3 functional epitopes and escape mutants by computer simulations. The results indicate the presence of significant structural differences between the functional epitopes and the escape mutants. Mutation-induced structural/conformational instabilities of this type can decrease the antibody neutralization activity. Among the different escape mutants studied here, Ser40Lys/Asp41Arg appear to be most unstable, while Ser40Glu/Asp41Leu exhibit the lowest structural variations. The highest level of escape mutation observed in Ser40Lys is linked to the relatively higher values of root mean square deviation/fluctuation found in the molecular dynamics simulation of this protein. Secondary-structure deviations and depletion of H bonding are other contributing factors to the protein’s increased instability. Overall, the proteins with residue 41 mutations are found to be structurally more ordered than those with residue 40 mutations. The detailed time-based structural assessment of the mutant epitopes described here may contribute to the development of novel vaccines and antiviral drugs necessary to defend against future outbreaks of JEV escape mutants.

Modeling and simulation study to identify threonine synthase as possible drug target in Leishmania major

Leishmaniasis is one of the most neglected tropical diseases that demand immediate attention to the identification of new drug targets and effective drug candidates. The present study demonstrates the possibility of using threonine synthase (TS) as a putative drug target in leishmaniasis disease management. We report the construction of an effective homology model of the enzyme that appears to be structurally as well as functionally well conserved. The 200 nanosecond molecular dynamics data on TS with and without pyridoxal phosphate (PLP) shed light on mechanistic details of PLP-induced conformational changes. Moreover, we address some important structural and dynamic interactions in the PLP binding region of TS that are in good agreement with previously speculated crystallographic estimations. Additionally, after screening more than 44,000 compounds, we propose 10 putative inhibitor candidates for TS based on virtual screening data and refined Molecular Mechanics Generalized Born Surface Area calculations. We expect that structural and functional dynamics data disclosed in this study will help initiate experimental endeavors toward establishing TS as an effective antileishmanial drug target.

Bioinformatics design and experimental validation of influenza A virus multi-epitopes that induce neutralizing antibodies

Pandemics caused by influenza A virus (IAV) are responsible for the deaths of millions of humans around the world. One of these pandemics occurred in Mexico in 2009. Despite the impact of IAV on human health, there is no effective vaccine. Gene mutations and translocation of genome segments of different IAV subtypes infecting a single host cell make the development of a universal vaccine difficult. The design of immunogenic peptides using bioinformatics tools could be an interesting strategy to increase the success of vaccines. In this work, we used the predicted amino acid sequences of the neuraminidase (NA) and hemagglutinin (HA) proteins of different IAV subtypes to perform multiple alignments, epitope predictions, molecular dynamics simulations, and experimental validation. Peptide selection was based on the following criteria: promiscuity, protein surface exposure, and the degree of conservation among different medically relevant IAV strains. These peptides were tested using immunological assays to test their ability to induce production of antibodies against IAV. We immunized rabbits and mice and measured the levels of IgG and IgA antibodies in serum samples and nasal washes. Rabbit antibodies against the peptides P11 and P14 (both of which are hybrids of NA and HA) recognized HA from both group 1 (H1, H2, and H5) and group 2 (H3 and H7) IAV and also recognized the purified NA protein from the viral stock (influenza A Puerto Rico/916/34). IgG antibodies from rabbits immunized with P11 and P14 were capable of recognizing viral particles and inhibited virus hemagglutination. Additionally, intranasal immunization of mice with P11 and P14 induced specific IgG and IgA antibodies in serum and nasal mucosa, respectively. Interestingly, the IgG antibodies were found to have neutralizing capability. In conclusion, the peptides designed through in silico studies were validated in experimental assays.

Immunity to influenza: Impact of obesity

Obesity is a health concern that is recognized as a critical factor for vulnerability to influenza A/pdmH1N1 virus infection, with epidemiological and clinical impacts. In humans, obesity induces disturbances in inflammatory and immune responses to the influenza virus and in some cases, this leads to severe complications, with fatal outcomes. Obesity impairs immunity by altering the response of cytokines, resulting in a decrease in the cytotoxic cell response of immunocompetent cells which have a key anti-viral role. Additionally, obesity seems to disturb the balance of endocrine hormones, such as leptin, that affect the interplay between metabolic and immune systems. This contribution focuses on reviewing the current epidemiologic data for the immune response to immunity in obese humans and animal models. In doing so, we aim to provide potential mechanisms to enhance immunity to influenza A/pdmH1N1 virus infection and protective factors in obese people.

Structural and functional characterization of the divergent Entamoeba Src using Src inhibitor-1

Background

The abundant number of kinases that Entamoeba histolytica possesses allows us to assume that the regulation of cellular functions by phosphorylation-dephosphorylation processes is very important. However, the kinases responsible for the phosphorylation in Entamoeba spp. vary in the structure of their domains and, therefore, could be responsible for the unusual biological characteristics of this parasite. In higher eukaryotes, Src kinases share conserved structural domains and are very important in the regulation of the actin cytoskeleton. In both Entamoeba histolytica and Entamoeba invadens, the major Src kinase homologue of higher eukaryotes lacks SH3 and SH2 domains, but does have KELCH domains; the latter are part of actin cross-linking proteins in higher eukaryotic cells.

Methods

The function of the EhSrc protein kinase of Entamoeba spp. was evaluated using Src inhibitor-1, microscopy assays, Src kinase activity and western blot. In addition, to define the potential inhibitory mechanism of Src-inhibitor-1 for the amoebic EhSrc protein kinase, molecular dynamic simulations using NAnoscale Molecular Dynamics (NAMD2) program and docking studies were performed with MOE software.

Results

We demonstrate that Src inhibitor-1 is able to prevent the activity of EhSrc protein kinase, most likely by binding to the catalytic domain, which affects cell morphology via the disruption of actin cytoskeleton remodeling and the formation of phagocytic structures without an effect on cell adhesion. Furthermore, in E. invadens, Src inhibitor-1 inhibited the encystment process by blocking RhoA GTPase activity, a small GTPase protein of Rho family.

Conclusions

Even though the EhSrc molecule of Entamoeba is not a typical Src, because its divergent amino acid sequence, it is a critical factor in the biology of this parasite via the regulation of actin cytoskeleton remodeling via RhoA GTPase activation. Based on this, we conclude that EhSrc could become a target molecule for the future design of drugs that can prevent the transmission of the disease.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2461-5) contains supplementary material, which is available to authorized users.

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.

Molecular modeling studies demonstrate key mutations that could affect the ligand recognition by influenza AH1N1 neuraminidase

The goal of this study was to identify neuraminidase (NA) residue mutants from human influenza AH1N1 using sequences from 1918 to 2012. Multiple alignment studies of complete NA sequences (5732) were performed. Subsequently, the crystallographic structure of the 1918 influenza (PDB ID: 3BEQ-A) was used as a wild-type structure and three-dimensional (3-D) template for homology modeling of the mutated selected NA sequences. The 3-D mutated NAs were refined using molecular dynamics (MD) simulations (50 ns). The refined 3-D models were used to perform docking studies using oseltamivir. Multiple sequence alignment studies showed seven representative mutations (A232V, K262R, V263I, T264V, S367L, S369N, and S369K). MD simulations applied to 3-D NAs showed that each NA had different active-site shapes according to structural surface visualization and docking results. Moreover, Cartesian principal component analyses (cPCA) show structural differences among these NA structures caused by mutations. These theoretical results suggest that the selected mutations that are located outside of the active site of NA could affect oseltamivir recognition and could be associated with resistance to oseltamivir.

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

Epitope recognition by major histocompatibility complex II (MHC-II) is essential for the activation of immunological responses to infectious diseases. Several studies have demonstrated that this molecular event takes place in the MHC-II peptide-binding groove constituted by the α and β light chains of the heterodimer. This MHC-II peptide-binding groove has several pockets (P1-P11) involved in peptide recognition and complex stabilization that have been probed through crystallographic experiments and in silico calculations. However, most of these theoretical calculations have been performed without taking into consideration the heavy chains, which could generate misleading information about conformational mobility both in water and in the membrane environment. Therefore, in absence of structural information about the difference in the conformational changes between the peptide-free and peptide-bound states (pMHC-II) when the system is soluble in an aqueous environment or non-covalently bound to a cell membrane, as the physiological environment for MHC-II is. In this study, we explored the mechanistic basis of these MHC-II components using molecular dynamics (MD) simulations in which MHC-II was previously co-crystallized with a small epitope (P₇) or coupled by docking procedures to a large (P₂₂) epitope. These MD simulations were performed at 310 K over 100 ns for the water-soluble (MHC-II_w, MHC-II-P_7w, and MHC-II-P_22w) and 150 ns for the membrane-bound species (MHC-II_m, MHC-II-P_7m, and MHC-II-P_22m). Our results reveal that despite the different epitope sizes and MD simulation environments, both peptides are stabilized primarily by residues lining P1, P4, and P6-7, and similar noncovalent intermolecular energies were observed for the soluble and membrane-bound complexes. However, there were remarkably differences in the conformational mobility and intramolecular energies upon complex formation, causing some differences with respect to how the two peptides are stabilized in the peptide-binding groove.