Re-mapping the molecular features of the human immunodeficiency virus type 1 and human T-cell lymphotropic virus type 1 Brazilian sequences using a bioinformatics unit established in Salvador, Bahia, Brazil, to give support to the viral epidemiology studies

The analysis of genetic data for human immunodeficiency virus type 1 (HIV-1) and human T-cell lymphotropic virus type 1 (HTLV-1) is essential to improve treatment and public health strategies as well as to select strains for vaccine programs. However, the analysis of large quantities of genetic data requires collaborative efforts in bioinformatics, computer biology, molecular biology, evolution, and medical science. The objective of this study was to review and improve the molecular epidemiology of HIV-1 and HTLV-1 viruses isolated in Brazil using bioinformatic tools available in the Laboratório Avançado de Sáude Pública (Lasp) bioinformatics unit. The analysis of HIV-1 isolates confirmed a heterogeneous distribution of the viral genotypes circulating in the country. The Brazilian HIV-1 epidemic is characterized by the presence of multiple subtypes (B, F1, C) and B/F1 recombinant virus while, on the other hand, most of the HTLV-1 sequences were classified as Transcontinental subgroup of the Cosmopolitan subtype. Despite the high variation among HIV-1 subtypes, protein glycosylation and phosphorylation domains were conserved in the pol, gag, and env genes of the Brazilian HIV-1 strains suggesting constraints in the HIV-1 evolution process. As expected, the functional protein sites were highly conservative in the HTLV-1 env gene sequences. Furthermore, the presence of these functional sites in HIV-1 and HTLV-1 strains could help in the development of vaccines that pre-empt the viral escape process.

Structural bioinformatics study of cyclin-dependent kinases complexed with inhibitors

The present work describes molecular models for the binary complexes CDK9, CDK5 and CDK1 complexed with Flavopiridol and Roscovitine. These structural models indicate that the inhibitors strongly bind to the ATP-binding pocket of CDKs and the structural comparison with the complexes CDK2:Flavopiridol and CDK2:Roscovitine correlates the structural differences with differences in inhibition of these CDKs by the inhibitors. These structures open the possibility of testing new inhibitor families, in addition to new substituents for the already known lead structures such as flavones and adenine derivatives.

Cloning, bioinformatics analysis, and expression of the dust mite allergen Der f 5 of Dermatophagoides farinae

Crude extracts of house dust mites are used clinically for diagnosis and immunotherapy of allergic diseases, including bronchial asthma, perennial rhinitis, and atopic dermatitis. However, crude extracts are complexes with non-allergenic antigens and lack effective concentrations of important allergens, resulting in several side effects. Dermatophagoides farinae (Hughes; Acari: Pyroglyphidae) is one of the predominant sources of dust mite allergens, which has more than 30 groups of allergen. The cDNA coding for the group 5 allergen of D. farinae from China was cloned, sequenced and expressed. According to alignment using the VECTOR NTI 9.0 software, there were eight mismatched nucleotides in five cDNA clones resulting in seven incompatible amino acid residues, suggesting that the Der f 5 allergen might have sequence polymorphism. Bioinformatics analysis revealed that the matured Der f 5 allergen has a molecular mass of 13604.03 Da, a theoretical pI of 5.43 and is probably hydrophobic and cytoplasmic. Similarities in amino acid sequences between Der f 5 and allergens of other domestic mite species, viz. Der p 5, Blo t 5, Sui m 5, and Lep d 5, were 79, 48, 53, and 37%, respectively. Phylogenetic analysis indicated that Der f 5 and Der p 5 clustered together. Blo t 5 and Ale o 5 also clustered together, although Blomia tropicalis and Aleuroglyphus ovatus belong to different mite families, viz. Echimyopodidae and Acaridae, respectively.

Bioinformatics analysis of biomarkers and transcriptional factor motifs in Down syndrome

In this study, biomarkers and transcriptional factor motifs were identified in order to investigate the etiology and phenotypic severity of Down syndrome. GSE 1281, GSE 1611, and GSE 5390 were downloaded from the gene expression ominibus (GEO). A robust multiarray analysis (RMA) algorithm was applied to detect differentially expressed genes (DEGs). In order to screen for biological pathways and to interrogate the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, the database for annotation, visualization, and integrated discovery (DAVID) was used to carry out a gene ontology (GO) function enrichment for DEGs. Finally, a transcriptional regulatory network was constructed, and a hypergeometric distribution test was applied to select for significantly enriched transcriptional factor motifs. CBR1, DYRK1A, HMGN1, ITSN1, RCAN1, SON, TMEM50B, and TTC3 were each up-regulated two-fold in Down syndrome samples compared to normal samples; of these, SON and TTC3 were newly reported. CBR1, DYRK1A, HMGN1, ITSN1, RCAN1, SON, TMEM50B, and TTC3 were located on human chromosome 21 (mouse chromosome 16). The DEGs were significantly enriched in macromolecular complex subunit organization and focal adhesion pathways. Eleven significantly enriched transcription factor motifs (PAX5, EGR1, XBP1, SREBP1, OLF1, MZF1, NFY, NFKAPPAB, MYCMAX, NFE2, and RP58) were identified. The DEGs and transcription factor motifs identified in our study provide biomarkers for the understanding of Down syndrome pathogenesis and progression.

DETERMINATION OF VIRAL TROPISM BY GENOTYPING AND PHENOTYPING ASSAYS IN BRAZILIAN HIV-1-INFECTED PATIENTS

The clinical application of CCR5 antagonists involves first determining the coreceptor usage by the infecting viral strain. Bioinformatics programs that predict coreceptor usage could provide an alternative method to screen candidates for treatment with CCR5 antagonists, particularly in countries with limited financial resources. Thus, the present study aims to identify the best approach using bioinformatics tools for determining HIV-1 coreceptor usage in clinical practice. Proviral DNA sequences and Trofile results from 99 HIV-1-infected subjects under clinical monitoring were analyzed in this study. Based on the Trofile results, the viral variants present were 81.1% R5, 21.4% R5X4 and 1.8% X4. Determination of tropism using a Geno2pheno[coreceptor] analysis with a false positive rate of 10% gave the most suitable performance in this sampling: the R5 and X4 strains were found at frequencies of 78.5% and 28.4%, respectively, and there was 78.6% concordance between the phenotypic and genotypic results. Further studies are needed to clarify how genetic diversity amongst virus strains affects bioinformatics-driven approaches for determining tropism. Although this strategy could be useful for screening patients in developing countries, some limitations remain that restrict the wider application of coreceptor usage tests in clinical practice.

Parasitic diseases: opportunities and challenges in the 21st century

The opportunities and challenges for the study and control of parasitic diseases in the 21st century are both exciting and daunting. Based on the contributions from this field over the last part of the 20th century, we should expect new biologic concepts will continue to come from this discipline to enrich the general area of biomedical research. The general nature of such a broad category of infections is difficult to distill, but they often depend on well-orchestrated, complex life cycles and they often involve chronic, relatively well-balanced host/parasite relationships. Such characteristics force biological systems to their limits, and this may be why studies of these diseases have made fundamental contributions to molecular biology, cell biology and immunology. However, if these findings are to continue apace, parasitologists must capitalize on the new findings being generated though genomics, bioinformatics, proteomics, and genetic manipulations of both host and parasite. Furthermore, they must do so based on sound biological insights and the use of hypothesis-driven studies of these complex systems. A major challenge over the next century will be to capitalize on these new findings and translate them into successful, sustainable strategies for control, elimination and eradication of the parasitic diseases that pose major public health threats to the physical and cognitive development and health of so many people worldwide.

Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Clustering of Schistosoma mansoni mRNA sequences and analysis of the most transcribed genes: implications in metabolism and biology of different developmental stages

The study of the Schistosoma mansoni genome, one of the etiologic agents of human schistosomiasis, is essential for a better understanding of the biology and development of this parasite. In order to get an overview of all S. mansoni catalogued gene sequences, we performed a clustering analysis of the parasite mRNA sequences available in public databases. This was made using softwares PHRAP and CAP3. The consensus sequences, generated after the alignment of cluster constituent sequences, allowed the identification by database homology searches of the most expressed genes in the worm. We analyzed these genes and looked for a correlation between their high expression and parasite metabolism and biology. We observed that the majority of these genes is related to the maintenance of basic cell functions, encoding genes whose products are related to the cytoskeleton, intracellular transport and energy metabolism. Evidences are presented here that genes for aerobic energy metabolism are expressed in all the developmental stages analyzed. Some of the most expressed genes could not be identified by homology searches and may have some specific functions in the parasite.

Towards an understanding of the epigenetics of schistosomes: a comparative epigenomic study

As in perhaps all eukaryotes, schistosomes use a supplementary information transmitting system, the epigenetic inheritance system, to shape genetic information and to produce different phenotypes. In contrast to other important parasites, the study of epigenetic phenomena in schistosomes is still in its infancy. Nevertheless, we are beginning to grasp what goes on behind the epigenetic scene in this parasite. We have developed techniques of native chromatin immunoprecipitation (N-ChIP) and associated the necessary bioinformatics tools that allow us to run genome-wide comparative chromatin studies on Schistosoma mansoni at different stages of its life cycle, on different strains and on different sexes. We present here an application of such an approach to study the genetic and epigenetic basis for a phenotypic trait, the compatibility of S. mansoni with its invertebrate host Biomphalaria glabrata. We have applied the ChIP procedure to two strains that are either compatible or incompatible with their intermediate host. The precipitated DNA was sequenced and aligned to a reference genome and this information was used to determine regions in which both strands differ in their genomic sequence and/or chromatin structure. This procedure allowed us to identify candidate genes that display either genetic or epigenetic difference between the two strains.

Identification of Schistosoma mansoni candidate antigens for diagnosis of schistosomiasis

The development of a more sensitive diagnostic test for schistosomiasis is needed to overcome the limitations of the use of stool examination in low endemic areas. Using parasite antigens in enzyme linked immunosorbent assay is a promising strategy, however a more rational selection of parasite antigens is necessary. In this study we performed in silico analysis of the Schistosoma mansoni genome, using SchistoDB database and bioinformatic tools for screening immunogenic antigens. Based on evidence of expression in all parasite life stage within the definitive host, extracellular or plasmatic membrane localization, low similarity to human and other helminthic proteins and presence of predicted B cell epitopes, six candidates were selected: a glycosylphosphatidylinositol-anchored 200 kDa protein, two putative cytochrome oxidase subunits, two expressed proteins and one hypothetical protein. The recognition in unidimensional and bidimensional Western blot of protein with similar molecular weight and isoelectric point to the selected antigens by sera from S. mansoni infected mice indicate a good correlation between these two approaches in selecting immunogenic proteins.

Search for a platelet-activating factor receptor in the Trypanosoma cruzi proteome: a potential target for Chagas disease chemotherapy

Chagas disease (CD) causes the highest burden of parasitic diseases in the Western Hemisphere and is therefore a priority for drug research and development. Platelet-activating factor (PAF) causes the CD parasite Trypanosoma cruzi to differentiate, which suggests that the parasite may express PAF receptors. Here, we explored the T. cruzi proteome for PAF receptor-like proteins. From a total of 23,000 protein sequences, we identified 29 hypothetical proteins that are predicted to have seven transmembrane domains (TMDs), which is the main characteristic of the G protein-coupled receptors (GPCRs), including the PAF receptor. The TMDs of these sequences were independently aligned with domains from 25 animal PAF receptors and the sequences were analysed for conserved residues. The conservation score mean values for the TMDs of the hypothetical proteins ranged from 31.7-44.1%, which suggests that if the putative T. cruzi PAF receptor is among the sequences identified, the TMDs are not highly conserved. These results suggest that T. cruzi contains several GPCR-like proteins and that one of these GPCRs may be a PAF receptor. Future studies may further validate the PAF receptor as a target for CD chemotherapy.

Mycobacterium leprae virulence-associated peptides are indicators of exposure to M. leprae in Brazil, Ethiopia and Nepal

Silent transmission of Mycobacterium leprae, as evidenced by stable leprosy incidence rates in various countries, remains a health challenge despite the implementation of multidrug therapy worldwide. Therefore, the development of tools for the early diagnosis of M. leprae infection should be emphasised in leprosy research. As part of the continuing effort to identify antigens that have diagnostic potential, unique M. leprae peptides derived from predicted virulence-associated proteins (group IV.A) were identified using advanced genome pattern programs and bioinformatics. Based on human leukocyte antigen (HLA)-binding motifs, we selected 21 peptides that were predicted to be promiscuous HLA-class I T-cell epitopes and eight peptides that were predicted to be HLA-class II restricted T-cell epitopes for field-testing in Brazil, Ethiopia and Nepal. High levels of interferon (IFN)-γ were induced when peripheral blood mononuclear cells (PBMCs) from tuberculoid/borderline tuberculoid leprosy patients located in Brazil and Ethiopia were stimulated with the ML2055 p35 peptide. PBMCs that were isolated from healthy endemic controls living in areas with high leprosy prevalence (EChigh) in Ethiopia also responded to the ML2055 p35 peptide. The Brazilian EChigh group recognised the ML1358 p20 and ML1358 p24 peptides. None of the peptides were recognised by PBMCs from healthy controls living in non-endemic region. In Nepal, mixtures of these peptides induced the production of IFN-γ by the PBMCs of leprosy patients and EChigh. Therefore, the M. leprae virulence-associated peptides identified in this study may be useful for identifying exposure to M. leprae in population with differing HLA polymorphisms.

Identification of the nicotinamide mononucleotide adenylyltransferase of Trypanosoma cruzi

The intracellular parasite Trypanosoma cruzi is the aetiological agent of Chagas disease, a public health concern with an increasing incidence rate. This increase is due, among other reasons, to the parasite’s drug resistance mechanisms, which require nicotinamide adenine dinucleotide (NAD+). Furthermore, this molecule is involved in metabolic and intracellular signalling processes necessary for the survival of T. cruzi throughout its life cycle. NAD+ biosynthesis is performed by de novo and salvage pathways, which converge on the step that is catalysed by the enzyme nicotinamide mononucleotide adenylyltransferase (NMNAT) (enzyme commission number: 2.7.7.1). The identification of the NMNAT of T. cruzi is important for the development of future therapeutic strategies to treat Chagas disease. In this study, a hypothetical open reading frame (ORF) for NMNAT was identified in the genome of T. cruzi. The corresponding putative protein was analysed by simulating structural models. The ORF was amplified from genomic DNA by polymerase chain reaction and was further used for the construction of a corresponding recombinant expression vector. The expressed recombinant protein was partially purified and its activity was evaluated using enzymatic assays. These results comprise the first identification of an NMNAT in T. cruzi using bioinformatics and experimental tools and hence represent the first step to understanding NAD+ metabolism in these parasites.

Assessment of Epstein-Barr virus nucleic acids in gastric but not in breast cancer by next-generation sequencing of pooled Mexican samples

Gastric (GC) and breast (BrC) cancer are two of the most common and deadly tumours. Different lines of evidence suggest a possible causative role of viral infections for both GC and BrC. Wide genome sequencing (WGS) technologies allow searching for viral agents in tissues of patients with cancer. These technologies have already contributed to establish virus-cancer associations as well as to discovery new tumour viruses. The objective of this study was to document possible associations of viral infection with GC and BrC in Mexican patients. In order to gain idea about cost effective conditions of experimental sequencing, we first carried out an in silico simulation of WGS. The next-generation-platform IlluminaGallx was then used to sequence GC and BrC tumour samples. While we did not find viral sequences in tissues from BrC patients, multiple reads matching Epstein-Barr virus (EBV) sequences were found in GC tissues. An end-point polymerase chain reaction confirmed an enrichment of EBV sequences in one of the GC samples sequenced, validating the next-generation sequencing-bioinformatics pipeline.