The CDR3 region as the major driver of TREM-1 interaction with its ligands, an in silico characterization.

The triggering receptor expressed on myeloid cells-1 (TREM-1) is a pattern recognition receptor heavily investigated in infectious and non-infectious diseases. Because of its role in amplifying inflammation, TREM-1 has been explored as a diagnostic/prognostic biomarker. Further, as the receptor has been implicated in the pathophysiology of several diseases, therapies aiming at modulating its activity represent a promising strategy to constrain uncontrolled inflammatory or infectious diseases. Despite this, several aspects concerning its interaction with ligands and activation process, remain unclear. Although many molecules have been suggested as TREM-1 ligands, only five have been confirmed to interact with the receptor: actin, eCIRP, HMGB1, Hsp70 and PGLYRP1. However, the domains involved in the interaction between the receptor and these proteins are not clarified yet. Therefore, here we used in silico approaches to investigate the putative binding domains in the receptor, using hot spots analysis, molecular docking and molecular dynamics simulations between TREM-1 and its five known ligands. Our results indicated the complementarity-determining regions (CDRs) of the receptor as the main mediators of antigen recognition, especially the CDR3 loop. We believe that our study could be used as structural basis for the elucidation of TREM-1's recognition process, and may be useful for prospective in silico and biological investigations exploring the receptor in different contexts.

Immunoinformatics-guided design of a multi-valent vaccine against Rotavirus and Norovirus (ChRNV22).

Rotavirus (RV) and Norovirus (NV) are the main viral etiologic agents of acute gastroenteritis (AG), a serious pediatric condition associated with significant death rates and long-term complications. Anti-RV vaccination has been proved efficient in the reduction of severe AG worldwide, however, the available vaccines are all attenuated and have suboptimal efficiencies in developing countries, where AG leads to substantial disease burden. On the other hand, no NV vaccine has been licensed so far. Therefore, we used immunoinformatics tools to develop a multi-epitope vaccine (ChRNV22) to prevent severe AG by RV and NV. Epitopes were predicted against 17 prevalent genotypes of four structural proteins (NV's VP1, RV's VP4, VP6 and VP7), and then assembled in a chimeric protein, with two small adjuvant sequences (tetanus toxin P2 epitope and a conserved sequence of RV's enterotoxin, NSP4). Simulations of the immune response and interactions with immune receptors indicated the immunogenic properties of ChRNV22, including a Th1-biased response. In silico search for putative host-homologous, allergenic and toxic regions also indicated the vaccine safety. In summary, we developed a multi-epitope vaccine against different NV and RV genotypes that seems promising for the prevention of severe AG, which will be further assessed by in vivo tests.

Polymorphisms in alternative oxidase genes from ecotypes of Arabidopsis and rice revealed an environment-induced linkage to altitude and rainfall.

We investigated SNPs in alternative oxidase (AOX) genes and their connection to ecotype origins (climate, altitude, and rainfall) by using genomic data sets of Arabidopsis and rice populations from 1190 and 90 ecotypes, respectively. Parameters were defined to detect non-synonymous SNPs in the AOX ORF, which revealed amino acid (AA) changes in AOX1c, AOX1d, and AOX2 from Arabidopsis and AOX1c from rice in comparison to AOX references from Columbia-0 and Japonica ecotypes, respectively. Among these AA changes, Arabidopsis AOX1c_A161E&G165R and AOX1c_R242S revealed a link to high rainfall and high altitude, respectively, while all other changes in Arabidopsis and rice AOX was connected to high altitude and rainfall. Comparative 3D modeling showed that all mutant AOX presented structural differences in relation to the respective references. Molecular docking analysis uncovered lower binding affinity values between AOX and the substrate ubiquinol for most of the identified structures compared to their reference, indicating better enzyme-substrate binding affinities. Thus, our in silico data suggest that the majority of the AA changes found in the available ecotypes will confer better enzyme-subtract interactions and thus indicate environment-related, more efficient AOX activity.

Optimization of Resveratrol Used as a Scaffold to Design Histone Deacetylase (HDAC-1 and HDAC-2) Inhibitors.

Histone deacetylases (HDAC) are epigenetic enzymes responsible for repressing gene expression through the deacetylation of histone lysine residues. Therefore, inhibition of HDACs has become an interesting approach for the treatment of several diseases, including cancer, hematology, neurodegenerative, immune diseases, bacterial infections, and more. Resveratrol (RVT) has pleiotropic effects, including pan-inhibition of HDAC isoforms; however, its ability to interfere with membranes requires additional optimization to eliminate nonspecific and off-target effects. Thus, to explore RVT as a scaffold, we designed a series of novel HDAC-1 and -2 inhibitors containing the 2-aminobenzamide subunit. Using molecular modeling, all compounds, except unsaturated compounds (4) and (7), exhibited a similar mode of interaction at the active sites of HDAC 1 and 2. The docking score values obtained from the study ranged from -12.780 to -10.967 Kcal/mol. All compounds were synthesized, with overall yields ranging from 33% to 67.3%. In an initial screening, compounds (4), (5), (7), and (20)-(26), showed enzymatic inhibitory effects ranging from 1 to 96% and 6 to 93% against HDAC-1 and HDAC-2, respectively. Compound (5), the most promising HDAC inhibitor in this series, was selected for IC50 assays, resulting in IC50 values of 0.44 µM and 0.37 µM against HDAC-1 and HDAC-2, respectively. In a panel of selectivity against HDACs 3-11, compound (5) presented selectivity towards Class I, mainly HDAC-1, 2, and 3. All compounds exhibited suitable physicochemical and ADMET properties as determined using in silico simulations. In conclusion, the optimization of the RVT structure allows the design of selective HDAC inhibitors, mainly targeting HDAC-1 and HDAC-2 isoforms.

Discovery of New Zika Protease and Polymerase Inhibitors through the Open Science Collaboration Project OpenZika.

The Zika virus (ZIKV) is a neurotropic arbovirus considered a global threat to public health. Although there have been several efforts in drug discovery projects for ZIKV in recent years, there are still no antiviral drugs approved to date. Here, we describe the results of a global collaborative crowdsourced open science project, the OpenZika project, from IBM's World Community Grid (WCG), which integrates different computational and experimental strategies for advancing a drug candidate for ZIKV. Initially, molecular docking protocols were developed to identify potential inhibitors of ZIKV NS5 RNA-dependent RNA polymerase (NS5 RdRp), NS3 protease (NS2B-NS3pro), and NS3 helicase (NS3hel). Then, a machine learning (ML) model was built to distinguish active vs inactive compounds for the cytoprotective effect against ZIKV infection. We performed three independent target-based virtual screening campaigns (NS5 RdRp, NS2B-NS3pro, and NS3hel), followed by predictions by the ML model and other filters, and prioritized a total of 61 compounds for further testing in enzymatic and phenotypic assays. This yielded five non-nucleoside compounds which showed inhibitory activity against ZIKV NS5 RdRp in enzymatic assays (IC50 range from 0.61 to 17 µM). Two compounds thermally destabilized NS3hel and showed binding affinity in the micromolar range (Kd range from 9 to 35 µM). Moreover, the compounds LabMol-301 inhibited both NS5 RdRp and NS2B-NS3pro (IC50 of 0.8 and 7.4 µM, respectively) and LabMol-212 thermally destabilized the ZIKV NS3hel (Kd of 35 µM). Both also protected cells from death induced by ZIKV infection in in vitro cell-based assays. However, while eight compounds (including LabMol-301 and LabMol-212) showed a cytoprotective effect and prevented ZIKV-induced cell death, agreeing with our ML model for prediction of this cytoprotective effect, no compound showed a direct antiviral effect against ZIKV. Thus, the new scaffolds discovered here are promising hits for future structural optimization and for advancing the discovery of further drug candidates for ZIKV. Furthermore, this work has demonstrated the importance of the integration of computational and experimental approaches, as well as the potential of large-scale collaborative networks to advance drug discovery projects for neglected diseases and emerging viruses, despite the lack of available direct antiviral activity and cytoprotective effect data, that reflects on the assertiveness of the computational predictions. The importance of these efforts rests with the need to be prepared for future viral epidemic and pandemic outbreaks.

Development of New Potential Inhibitors of ß1 Integrins through In Silico Methods-Screening and Computational Validation.

Integrins are transmembrane receptors that play a critical role in many biological processes which can be therapeutically modulated using integrin blockers, such as peptidomimetic ligands. This work aimed to develop new potential ß1 integrin antagonists using modeled receptors based on the aligned crystallographic structures and docked with three lead compounds (BIO1211, BIO5192, and TCS2314), widely known as α4ß1 antagonists. Lead-compound complex optimization was performed by keeping intact the carboxylate moiety of the ligand, adding substituents in two other regions of the molecule to increase the affinity with the target. Additionally, pharmacokinetic predictions were performed for the ten best ligands generated, with the lowest docking interaction energy obtained for α4ß1 and BIO5192. Results revealed an essential salt bridge between the BIO5192 carboxylate group and the Mg2+ MIDAS ion of the integrin. We then generated more than 200 new BIO5192 derivatives, some with a greater predicted affinity to α4ß1. Furthermore, the significance of retaining the pyrrolidine core of the ligand and increasing the therapeutic potential of the new compounds is emphasized. Finally, one novel molecule (1592) was identified as a potential drug candidate, with appropriate pharmacokinetic profiles, similar dynamic behavior at the integrin interaction site compared with BIO5192, and a higher predicted affinity to VLA-4.

Investigation of Unprecedented Sites and Proposition of New Ligands for Programmed Cell Death Protein I through Molecular Dynamics with Probes and Virtual Screening.

Cancer immunotherapy has attracted increasing attention over the last few years. Programmed cell death protein 1 (PD-1) promotes self-tolerance and inhibits immune responses by modulating the T-cell function. The interaction between PD-1 and programmed cell death ligand-1 (PD-L1) leads to immune exhaustion, protecting cancer cells from destruction. Here, we computationally designed a novel ligand named 1508 that binds to an unprecedented PD-1 cavity identified by MixMD and defined by amino acid residues Lys78 to Val97. We showed through a set of MD simulations totaling 12.5 µs that ligand 1508 establishes frequent cation-π and hydrogen bonding interactions with amino acid residues Lys78 and Arg86, respectively, and stabilizes the PD-1 C'D loop in a conformation that does not favor PD-1-PD-L1 complex formation. This study highlights the power of MixMD in exposing new cavities prone to protein-protein complex inhibition and establishes the basis for the design of new molecules that target the PD-1 C'D cavity as an alternative for exploring the modulation of the PD-1-PD-L1 complex in cancer therapy.

Naturally acquired antibody response to a Plasmodium falciparum chimeric vaccine candidate GMZ2.6c and its components (MSP-3, GLURP, and Pfs48/45) in individuals living in Brazilian malaria-endemic areas.

BACKGROUND: The GMZ2.6c malaria vaccine candidate is a multi-stage Plasmodium falciparum chimeric protein which contains a fragment of the sexual-stage Pfs48/45-6C protein genetically fused to GMZ2, a fusion protein of GLURP and MSP-3, that has been shown to be well tolerated, safe and immunogenic in clinical trials performed in a malaria-endemic area of Africa. However, there is no data available on the antigenicity or immunogenicity of GMZ2.6c in humans. Considering that circulating parasites can be genetically distinct in different malaria-endemic areas and that host genetic factors can influence the immune response to vaccine antigens, it is important to verify the antigenicity, immunogenicity and the possibility of associated protection in individuals living in malaria-endemic areas with different epidemiological scenarios. Herein, the profile of antibody response against GMZ2.6c and its components (MSP-3, GLURP and Pfs48/45) in residents of the Brazilian Amazon naturally exposed to malaria, in areas with different levels of transmission, was evaluated. METHODS: This study was performed using serum samples from 352 individuals from Cruzeiro do Sul and Mâncio Lima, in the state of Acre, and Guajará, in the state of Amazonas. Specific IgG, IgM, IgA and IgE antibodies and IgG subclasses were detected by Enzyme-Linked Immunosorbent Assay. RESULTS: The results showed that GMZ2.6c protein was widely recognized by naturally acquired antibodies from individuals of the Brazilian endemic areas with different levels of transmission. The higher prevalence of individuals with antibodies against GMZ2.6c when compared to its individual components may suggest an additive effect of GLURP, MSP-3, and Pfs48/45 when inserted in a same construct. Furthermore, naturally malaria-exposed individuals predominantly had IgG1 and IgG3 cytophilic anti-GMZ2.6c antibodies, an important fact considering that the acquisition of anti-malaria protective immunity results from a delicate balance between cytophilic/non-cytophilic antibodies. Interestingly, anti-GMZ2.6c antibodies seem to increase with exposure to malaria infection and may contribute to parasite immunity. CONCLUSIONS: The data showed that GMZ2.6c protein is widely recognized by naturally acquired antibodies from individuals living in malaria-endemic areas in Brazil and that these may contribute to parasite immunity. These data highlight the importance of GMZ2.6c as a candidate for an anti-malarial vaccine.

Computationally-obtained structural insights into the molecular interactions between Pidilizumab and binding partners DLL1 and PD-1.

Pidilizumab is a monoclonal antibody tested against several types of malignancies, such as lymphoma and metastatic melanoma, showing promising results. In 2016, the FDA put Pidilizumab's clinical studies on partial hold due to emerging evidence pointing to the antibody target uncertainty. Although initial studies indicated an interaction with the PD-1 checkpoint receptor, recent updates assert that Pidilizumab binds primarily to Notch ligand DLL1. However, a detailed description of which interactions coordinate antibody-antigen complex formation is lacking. Therefore, this study uses computational tools to identify molecular interactions between Pidilizumab and its reported targets PD-1 and DLL1. A docking methodology was validated and applied to determine the binding modes between modeled Pidilizumab scFvs and the two antigens. We used Molecular Dynamics (MD) simulations to verify the complexes' stability and submitted the resulting trajectory files to MM/PBSA and Principal Component Analysis. A set of different prediction tools determined scFv interface hot-spots. Whereas docking and MD simulations revealed that the antibody fragments do not interact straightforwardly with PD-1, ten scFv hot-spots, including Met93 and Leu112, mediated the interaction with the DLL1 C2 domain. The interaction triggered a conformational selection-like effect on DLL1, allowing new hydrogen bonds on the ß3-ß4 interface loop. The unprecedented structural data on Pidilizumab's interactions provided novel evidence that its legitimate target is the DLL1 protein and offered structural insight on how these molecules interact, shedding light on the pathways that could be affected by the use of this essential immunobiological. Communicated by Ramaswamy H. Sarma.

VLA-4 as a Central Target for Modulating Neuroinflammatory Disorders.

The complex steps leading to the central nervous system (CNS) inflammation and the progress to neuroinflammatory and neurodegenerative disorders have opened up new research and intervention avenues. This review focuses on the therapeutic targeting of the VLA-4 integrin to discuss the clear-cut effect on immune cell trafficking into brain tissues. Besides, we explore the possibility that blocking VLA-4 may have a relevant impact on nonmigratory activities of immune cells, such as antigen presentation and T-cell differentiation, during the neuroinflammatory process. Lastly, the recent refinement of computational techniques is highlighted as a way to increase specificity and to reduce the detrimental side effects of VLA-4 immunotherapies aiming at developing better clinical interventions.

Plasmodium vivax metacaspase 1 (PvMCA1) catalytic domain is conserved in field isolates from Brazilian Amazon.

In the present study, we investigated the genetic diversity of Plasmodium vivax metacaspase 1 (PvMCA1) catalytic domain in two municipalities of the main malaria hotspot in Brazil, i.e., the Juruá Valley, and observed complete sequence identity among all P. vivax field isolates and the Sal-1 reference strain. Analysis of PvMCA1 catalytic domain in different P. vivax genomic sequences publicly available also revealed a high degree of conservation worldwide, with very few amino acid substitutions that were not related to putative histidine and cysteine catalytic residues, whose involvement with the active site of protease was herein predicted by molecular modeling. The genetic conservation presented by PvMCA1 may contribute to its eligibility as a druggable target candidate in vivax malaria.

New data on Beurlenia araripensis Martins-Neto & Mezzalira, 1991, a lacustrine shrimp from Crato Formation, and its morphological variations based on the shape and the number of rostral spines.

Fossil freshwater carideans are very rare worldwide. Here, we present new taxonomic remarks about Beurlenia araripensis from the Early Cretaceous laminated limestones of the Crato Formation, Araripe Basin, northeastern Brazil. We analyzed five fossil samples, testing the morphological variations such as, rostrum with 5 to 14 supra-rostral spines and 2 to 3 sub-rostral spines, which appears as serrate for Caridea. This variation demonstrates a morphologic plasticity also seen in extant species of the group, such as those of the genera Macrobrachium and Palaemon.

Plasmodium vivax metacaspase 1 (PvMCA1) catalytic domain is conserved in field isolates from Brazilian Amazon

In the present study, we investigated the genetic diversity of Plasmodium vivax metacaspase 1 (PvMCA1) catalytic domain in two municipalities of the main malaria hotspot in Brazil, i.e., the Juruá Valley, and observed complete sequence identity among all P. vivax field isolates and the Sal-1 reference strain. Analysis of PvMCA1 catalytic domain in different P. vivax genomic sequences publicly available also revealed a high degree of conservation worldwide, with very few amino acid substitutions that were not related to putative histidine and cysteine catalytic residues, whose involvement with the active site of protease was herein predicted by molecular modeling. The genetic conservation presented by PvMCA1 may contribute to its eligibility as a druggable target candidate in vivax malaria.

Equilibrium, kinetics and thermodynamics of lead (II) adsorption in bioadsorvent composed by Caryocar coriaceum Wittm barks.

The improper disposal of waste containing toxic metals promotes serious environmental problems. Biosorption stands out as a powerful technique in removing toxic metals from aqueous solutions, making it an efficient and low-cost technique, since it makes use of abundant materials and of natural origins. This work aimed to study the application of the bark of Caryocar coriaceum Wittm, as a possible bioabsorbent for the removal of Lead (II) present in aqueous solutions, for this, residues of this material were collected in a public market, where they were washed and dried, then crushed, sieved (32mesh) and washed again, to then be used in the kinetic, isotherm and thermodynamic tests of lead (II) ions adsorption. After the experiments, it was found that the increase of the temperature disfavored the rate of lead (II) adsorption, ranging from 15.6 to 13.5 mg.L-1, for temperatures of 5 °C and 35 °C, respectively. The pseudo-second order kinetic model best characterizes the adsorption velocity obtaining a linear regression coefficient of 0.9997, while the pseudo-first order obtained a value of 0.7117. Adsorption occurred by chemisorption, with intra-film diffusion. The isothermal model that best represented the adsorption was Langmuir's (R2 0.99). It was concluded that the bark of Caryocar coriaceum Wittm is a good alternative to be used as a bioabsorbent for solutions containing dissolved lead (II).

Antileishmanial drugs activate inflammatory signaling pathways via toll-like receptors (docking approach) from Leishmania amazonensis-infected macrophages.

The activation of proinflammatory cellular processes and signals such as those linked to NF-kB in macrophages are involved in the control of infection by Leishmania ssp. However, little is known about the influence of the drugs used in the treatment on the host cellular inflammatory signaling pathways. This study aimed to evaluate the effects of different drugs used in the treatment of leishmaniasis on inflammatory profile related to Toll-like receptors (TLRs) from L. amazonensis-infected macrophages. J774 macrophage-like cells were infected with the promastigote forms (5:1) and 24 hs incubated with Amphotericin B (AmB), Glucantime® (GLU) or Pentamidine (Pent). The following inflammatory pathways were evaluated: NF-κB p65, NF-κB p65 phosphorylated (Ser536), stress-activated protein kinase/c-Jun NH(2)-terminal kinase (SAPK/JNK) phosphorylated (Thr183/Tyr185), p38 mitogen activated protein kinase (MAPK p38) phosphorylated (Thr180/Tyr182), signal transducer and activator of transcription-3 (Stat3) phosphorylated (Tyr705) and inhibitor kappa B-α (IκB-α) phosphorylated (Ser32). In silico tests were performed to evaluate the molecular affinity between TLRs and antileishmanial drugs. Molecular docking showed that affinities varied significantly among the binders evaluated. The lowest affinity (-8.6 Kcal/Mol) was calculated for AmB in complex with TLR4. Pent showed higher values for TLR1, TLR2 and TLR3, while for TLR4 the affinity value was lower (5.5 Kcal/Mol). The values obtained for GLU were the highest for the set of binders tested. From the infected macrophages, treatments inhibited NF-kB p65 for GLU (65.44%), for Pent (46.43%) and for AmB (54.07%) compared to untreated infected macrophages. The activation of the signaling pathway of NF-kB, SAPK/JNK and IκB-α caused by AmB and Pent may potentiate the microbicidal mechanisms of the infected macrophages.

Differential expression of recently duplicated PTOX genes in Glycine max during plant development and stress conditions.

Plastid terminal oxidase (PTOX) is a chloroplast enzyme that catalyzes oxidation of plastoquinol (PQH2) and reduction of molecular oxygen to water. Its function has been associated with carotenoid biosynthesis, chlororespiration and environmental stress responses in plants. In the majority of plant species, a single gene encodes the protein and little is known about events of PTOX gene duplication and their implication to plant metabolism. Previously, two putative PTOX (PTOX1 and 2) genes were identified in Glycine max, but the evolutionary origin and the specific function of each gene was not explored. Phylogenetic analyses revealed that this gene duplication occurred apparently during speciation involving the Glycine genus ancestor, an event absent in all other available plant leguminous genomes. Gene expression evaluated by RT-qPCR and RNA-seq data revealed that both PTOX genes are ubiquitously expressed in G. max tissues, but their mRNA levels varied during development and stress conditions. In development, PTOX1 was predominant in young tissues, while PTOX2 was more expressed in aged tissues. Under stress conditions, the PTOX transcripts varied according to stress severity, i.e., PTOX1 mRNA was prevalent under mild or moderate stresses while PTOX2 was predominant in drastic stresses. Despite the high identity between proteins (97%), molecular docking revealed that PTOX1 has higher affinity to substrate plastoquinol than PTOX2. Overall, our results indicate a functional relevance of this gene duplication in G. max metabolism, whereas PTOX1 could be associated with chloroplast effectiveness and PTOX2 to senescence and/or apoptosis.

Polymorphisms in plastoquinol oxidase (PTOX) from Arabidopsis accessions indicate SNP-induced structural variants associated with altitude and rainfall.

Plant plastoquinol oxidase (PTOX) is a chloroplast oxidoreductase involved in carotenoid biosynthesis, chlororespiration, and response to environmental stresses. The present study aimed to gain insight of the potential role of nucleotide/amino acid changes linked to environmental adaptation in PTOX gene/protein from Arabidopsis thaliana accessions. SNPs in the single-copy PTOX gene were identified in 1190 accessions of Arabidopsis using the Columbia-0 PTOX as a reference. The identified SNPs were correlated with geographical distribution of the accessions according to altitude, climate, and rainfall. Among the 32 identified SNPs in the coding region of the PTOX gene, 16 of these were characterized as non-synonymous SNPs (in which an AA is altered). A higher incidence of AA changes occurred in the mature protein at positions 78 (31%), 81 (31.4%), and 323 (49.9%). Three-dimensional structure prediction indicated that the AA change at position 323 (D323N) leads to a PTOX structure with the most favorable interaction with the substrate plastoquinol, when compared with the reference PTOX structure (Columbia-0). Molecular docking analysis suggested that the most favorable D323N PTOX-plastoquinol interaction is due to a better enzyme-substrate binding affinity. The molecular dynamics revealed that plastoquinol should be more stable in complex with D323N PTOX, likely due a restraint mechanism in this structure that stabilize plastoquinol inside of the reaction center. The integrated analysis made from accession geographical distribution and PTOX SNPs indicated that AA changes in PTOX are related to altitude and rainfall, potentially due to an adaptive positive environmental selection.

Plasmodium vivax Cell-Traversal Protein for Ookinetes and Sporozoites: Naturally Acquired Humoral Immune Response and B-Cell Epitope Mapping in Brazilian Amazon Inhabitants.

The cell-traversal protein for ookinetes and sporozoites (CelTOS), a highly conserved antigen involved in sporozoite motility, plays an important role in the traversal of host cells during the preerythrocytic stage of Plasmodium species. Recently, it has been considered an alternative target when designing novel antimalarial vaccines against Plasmodium falciparum. However, the potential of Plasmodium vivax CelTOS as a vaccine target is yet to be explored. This study evaluated the naturally acquired immune response against a recombinant P. vivax CelTOS (PvCelTOS) (IgG and IgG subclass) in 528 individuals from Brazilian Amazon, as well as the screening of B-cell epitopes in silico and peptide assays to associate the breadth of antibody responses of those individuals with exposition and/or protection correlates. We show that PvCelTOS is naturally immunogenic in Amazon inhabitants with 94 individuals (17.8%) showing specific IgG antibodies against the recombinant protein. Among responders, the IgG reactivity indexes (RIs) presented a direct correlation with the number of previous malaria episodes (p = 0.003; r = 0.315) and inverse correlation with the time elapsed from the last malaria episode (p = 0.031; r = -0.258). Interestingly, high responders to PvCelTOS (RI > 2) presented higher number of previous malaria episodes, frequency of recent malaria episodes, and ratio of cytophilic/non-cytophilic antibodies than low responders (RI < 2) and non-responders (RI < 1). Moreover, a high prevalence of the cytophilic antibody IgG1 over all other IgG subclasses (p < 0.0001) was observed. B-cell epitope mapping revealed five immunogenic regions in PvCelTOS, but no associations between the specific IgG response to peptides and exposure/protection parameters were found. However, the epitope (PvCelTOSI136-E143) was validated as a main linear B-cell epitope, as 92% of IgG responders to PvCelTOS were also responders to this peptide sequence. This study describes for the first time the natural immunogenicity of PvCelTOS in Amazon individuals and identifies immunogenic regions in a full-length protein. The IgG magnitude was mainly composed of cytophilic antibodies (IgG1) and associated with recent malaria episodes. The data presented in this paper add further evidence to consider PvCelTOS as a vaccine candidate.

Plasmodium vivax Cell Traversal Protein for Ookinetes and Sporozoites (PvCelTOS) gene sequence and potential epitopes are highly conserved among isolates from different regions of Brazilian Amazon.

The Plasmodium vivax Cell-traversal protein for ookinetes and sporozoites (PvCelTOS) plays an important role in the traversal of host cells. Although essential to PvCelTOS progress as a vaccine candidate, its genetic diversity remains uncharted. Therefore, we investigated the PvCelTOS genetic polymorphism in 119 field isolates from five different regions of Brazilian Amazon (Manaus, Novo Repartimento, Porto Velho, Plácido de Castro and Oiapoque). Moreover, we also evaluated the potential impact of non-synonymous mutations found in the predicted structure and epitopes of PvCelTOS. The field isolates showed high similarity (99.3% of bp) with the reference Sal-1 strain, presenting only four Single-Nucleotide Polymorphisms (SNP) at positions 24A, 28A, 109A and 352C. The frequency of synonymous C109A (82%) was higher than all others (p<0.0001). However, the non-synonymous G28A and G352C were observed in 9.2% and 11.7% isolates. The great majority of the isolates (79.8%) revealed complete amino acid sequence homology with Sal-1, 10.9% presented complete homology with Brazil I and two undescribed PvCelTOS sequences were observed in 9.2% field isolates. Concerning the prediction analysis, the N-terminal substitution (Gly10Ser) was predicted to be within a B-cell epitope (PvCelTOS Accession Nos. AB194053.1) and exposed at the protein surface, while the Val118Leu substitution was not a predicted epitope. Therefore, our data suggest that although G28A SNP might interfere in potential B-cell epitopes at PvCelTOS N-terminal region the gene sequence is highly conserved among the isolates from different geographic regions, which is an important feature to be taken into account when evaluating its potential as a vaccine candidate.

In silico Identification and Validation of a Linear and Naturally Immunogenic B-Cell Epitope of the Plasmodium vivax Malaria Vaccine Candidate Merozoite Surface Protein-9.

Synthetic peptide vaccines provide the advantages of safety, stability and low cost. The success of this approach is highly dependent on efficient epitope identification and synthetic strategies for efficacious delivery. In malaria, the Merozoite Surface Protein-9 of Plasmodium vivax (PvMSP9) has been considered a vaccine candidate based on the evidence that specific antibodies were able to inhibit merozoite invasion and recombinant proteins were highly immunogenic in mice and humans. However the identities of linear B-cell epitopes within PvMSP9 as targets of functional antibodies remain undefined. We used several publicly-available algorithms for in silico analyses and prediction of relevant B cell epitopes within PMSP9. We show that the tandem repeat sequence EAAPENAEPVHENA (PvMSP9E795-A808) present at the C-terminal region is a promising target for antibodies, given its high combined score to be a linear epitope and located in a putative intrinsically unstructured region of the native protein. To confirm the predictive value of the computational approach, plasma samples from 545 naturally exposed individuals were screened for IgG reactivity against the recombinant PvMSP9-RIRII729-972 and a synthetic peptide representing the predicted B cell epitope PvMSP9E795-A808. 316 individuals (58%) were responders to the full repetitive region PvMSP9-RIRII, of which 177 (56%) also presented total IgG reactivity against the synthetic peptide, confirming it validity as a B cell epitope. The reactivity indexes of anti-PvMSP9-RIRII and anti-PvMSP9E795-A808 antibodies were correlated. Interestingly, a potential role in the acquisition of protective immunity was associated with the linear epitope, since the IgG1 subclass against PvMSP9E795-A808 was the prevalent subclass and this directly correlated with time elapsed since the last malaria episode; however this was not observed in the antibody responses against the full PvMSP9-RIRII. In conclusion, our findings identified and experimentally confirmed the potential of PvMSP9E795-A808 as an immunogenic linear B cell epitope within the P. vivax malaria vaccine candidate PvMSP9 and support its inclusion in future subunit vaccines.