João Hermínio Martins da Silva

Computational modeling of the bHLH domain of the transcription factor TWIST1 and R118C, S144R and K145E mutants

BackgroundHuman TWIST1 is a highly conserved member of the regulatory basic helix-loop-helix (bHLH) transcription factors. TWIST1 forms homo- or heterodimers with E-box proteins, such as E2A (isoforms E12 and E47), MYOD and HAND2. Haploinsufficiency germ-line mutations of the twist1 gene in humans are the main cause of Saethre-Chotzen syndrome (SCS), which is characterized by limb abnormalities and premature fusion of cranial sutures. Because of the importance of TWIST1 in the regulation of embryonic development and its relationship with SCS, along with the lack of an experimentally solved 3D structure, we performed comparative modeling for the TWIST1 bHLH region arranged into wild-type homodimers and heterodimers with E47. In addition, three mutations that promote DNA binding failure (R118C, S144R and K145E) were studied on the TWIST1 monomer. We also explored the behavior of the mutant forms in aqueous solution using molecular dynamics (MD) simulations, focusing on the structural changes of the wild-type versus mutant dimers.ResultsThe solvent-accessible surface area of the homodimers was smaller on wild-type dimers, which indicates that the cleft between the monomers remained more open on the mutant homodimers. RMSD and RMSF analyses indicated that mutated dimers presented values that were higher than those for the wild-type dimers. For a more careful investigation, the monomer was subdivided into four regions: basic, helix I, loop and helix II. The basic domain presented a higher flexibility in all of the parameters that were analyzed, and the mutant dimer basic domains presented values that were higher than the wild-type dimers. The essential dynamic analysis also indicated a higher collective motion for the basic domain.ConclusionsOur results suggest the mutations studied turned the dimers into more unstable structures with a wider cleft, which may be a reason for the loss of DNA binding capacity observed for in vitro circumstances.

Structural and Molecular Modeling Features of P2X Receptors

Currently, adenosine 5′-triphosphate (ATP) is recognized as the extracellular messenger that acts through P2 receptors. P2 receptors are divided into two subtypes: P2Y metabotropic receptors and P2X ionotropic receptors, both of which are found in virtually all mammalian cell types studied. Due to the difficulty in studying membrane protein structures by X-ray crystallography or NMR techniques, there is little information about these structures available in the literature. Two structures of the P2X4 receptor in truncated form have been solved by crystallography. Molecular modeling has proven to be an excellent tool for studying ionotropic receptors. Recently, modeling studies carried out on P2X receptors have advanced our knowledge of the P2X receptor structure-function relationships. This review presents a brief history of ion channel structural studies and shows how modeling approaches can be used to address relevant questions about P2X receptors.

Spectroscopic analysis and X-ray diffraction of trunk fossils from the Parnaíba Basin, Northeast Brazil.

The Parnaiba Sedimentary Basin is of the Paleozoic age and is located in Northeast Brazil, covering the states of Piauí, Maranhão and Tocantins and a small part of Ceará and Pará. In this work we applied several chemical analytical techniques to characterize trunk fossils found in the Parnaíba Sedimentary Basin, collected from four different sites, and discuss their fossilization process. We performed a study of the trunk fossils through X-ray diffraction, energy dispersive spectroscopy, infrared and Raman spectroscopy. The analysis allow us to identify the different compositions which are present in the trunk fossils: kaolinite (Al2Si2O5(OH)4), hematite (Fe2O3) and quartz (SiO2). Based in these results we were able to identify that the main fossilization mechanism of the trunk fossil was silicification. Furthermore, through Raman spectroscopy, we have observed the presence of carbonaceous materials in the Permian fossils, as evidenced by the D and G Raman bands. The relative intensities and bandwidths of the D and G bands indicated that the carbon has a low crystallinity. Thus, most of trunk fossils analyzed were permineralized and not petrified, because there is the presence of carbon that characterizes the partial decomposition of the organic matter in some trunks.

Molecular dynamics simulations of Zika virus NS3 helicase: Insights into RNA binding site activity

America is still suffering with the outbreak of Zika virus (ZIKV) infection. Congenital ZIKV syndrome has already caused a public health emergency of international concern. However, there are still no vaccines to prevent or drugs to treat the infection caused by ZIKV. The ZIKV NS3 helicase (NS3h) protein is a promising target for drug discovery due to its essential role in viral genome replication. NS3h unwinds the viral RNA to enable the replication of the viral genome by the NS5 protein. NS3h contains two important binding sites: the NTPase binding site and the RNA binding site. Here, we used molecular dynamics (MD) simulations to study the molecular behavior of ZIKV NS3h in the presence and absence of ssRNA and the potential implications for NS3h activity and inhibition. Although there is conformational variability and poor electron densities of the RNA binding loop in various apo flaviviruses NS3h crystallographic structures, the MD trajectories of NS3h-ssRNA demonstrated that the RNA binding loop becomes more stable when NS3h is occupied by RNA. Our results suggest that the presence of RNA generates important interactions with the RNA binding loop, and these interactions stabilize the loop sufficiently that it remains in a closed conformation. This closed conformation likely keeps the ssRNA bound to the protein for a sufficient duration to enable the unwinding/replication activities of NS3h to occur. In addition, conformational changes of this RNA binding loop can change the nature and location of the optimal ligand binding site, according to ligand binding site prediction results. These are important findings to help guide the design and discovery of new inhibitors of NS3h as promising compounds to treat the ZIKV infection.

Large-field electron imaging and X-ray elemental mapping unveil the morphology, structure, and fractal features of a Cretaceous fossil at the centimeter scale.

We used here a scanning electron microscopy approach that detected backscattered electrons (BSEs) and X-rays (from ionization processes) along a large-field (LF) scan, applied on a Cretaceous fossil of a shrimp (area ∼280 mm(2)) from the Araripe Sedimentary Basin. High-definition LF images from BSEs and X-rays were essentially generated by assembling thousands of magnified images that covered the whole area of the fossil, thus unveiling morphological and compositional aspects at length scales from micrometers to centimeters. Morphological features of the shrimp such as pleopods, pereopods, and antennae located at near-surface layers (undetected by photography techniques) were unveiled in detail by LF BSE images and in calcium and phosphorus elemental maps (mineralized as hydroxyapatite). LF elemental maps for zinc and sulfur indicated a rare fossilization event observed for the first time in fossils from the Araripe Sedimentary Basin: the mineralization of zinc sulfide interfacing to hydroxyapatite in the fossil. Finally, a dimensional analysis of the phosphorus map led to an important finding: the existence of a fractal characteristic (D = 1.63) for the hydroxyapatite-matrix interface, a result of physical-geological events occurring with spatial scale invariance on the specimen, over millions of years.

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.

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.

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 of inhibitors of ribose 5-phosphate isomerase from Trypanosoma cruzi using ligand and structure based approaches

Chagas disease, caused by the protozoan Trypanosoma cruzi, affects approximately seven million people, mainly in Latin America, and causes about 7000 deaths annually. The available treatments are unsatisfactory and search for more effective drugs against this pathogen is critical. In this context, the ribose 5-phosphate isomerase (Rpi) enzyme is a potential drug target mainly due to its function in the pentose phosphate pathway and its essentiality (previously shown in other trypanosomatids). In this study, we propose novel potential inhibitors for the Rpi of T. cruzi (TcRpi) based on a computer-aided approach, including structure-based and ligand-based pharmacophore modeling. Along with a substructural and similarity search, the selected pharmacophore hypotheses were used to screen the purchasable subset of the ZINC Database, yielding 20,183 candidate compounds. These compounds were submitted to molecular docking studies in the TcRpi and Human Rpi (HsRpi) active sites in order to identify potential selective inhibitors for the T. cruzi enzyme. After the molecular docking and ADME-T (absorption, distribution, metabolism, excretion and toxicity)/PAINS (pan-assay interference compounds) screenings, 211 molecules were selected as potential TcRpi inhibitors. Out of these, three compounds - ZINC36975961, ZINC63480117, and ZINC43763931 - were submitted to molecular dynamics simulations and two of them - ZINC36975961 and ZINC43763931- had good performance and made interactions with important active site residues over all the simulation time. These compounds could be considered potential TcRpi inhibitors candidates and also may be used as leads for developing new TcRpi inhibitors.