Structural modelling and comparative analysis of homologous, analogous and specific proteins from Trypanosoma cruzi versus Homo sapiens: putative drug targets for chagas' disease treatment

In silico approaches for predicting natural compounds with therapeutic potential and vaccine candidates against Streptococcus equi

Equine strangles is a prevalent disease that affects the upper respiratory in horses and is caused by the Gram-positive bacterium Streptococcus equi. In addition to strangles, other clinical conditions are caused by the two S. equi subspecies, equi and zooepidemicus, which present relevant zoonotic potential. Treatment of infections caused by S. equi has become challenging due to the worldwide spreading of infected horses and the unavailability of effective therapeutics and vaccines. Penicillin treatment is often recommended, but multidrug resistance issues arised. We explored the whole genome sequence of 18 S. equi isolates to identify candidate proteins to be targeted by natural drug-like compounds or explored as immunogens. We considered only proteins shared among the sequenced strains of subspecies equi and zooepidemicus, absent in the equine host and predicted to be essential and involved in virulence. Of these, 4 proteins with cytoplasmic subcellular location were selected for molecular docking with a library of 5008 compounds, while 6 proteins were proposed as prominent immunogens against S. equi due to their probabilities of behaving as adhesins. The molecular docking analyses revealed the best ten ligands for each of the 4 drug target candidates, and they were ranked according to their binding affinities and the number of hydrogen bonds for complex stability. Finally, the natural 5-ring compound C25H20F3N5O3 excelled in molecular dynamics simulations for the increased stability in the interaction with UDP-N-acetylenolpyruvoylglucosamine reductase (MurB). This research paves the way to developing new therapeutics to minimize the impacts caused by S. equi infections.Communicated by Ramaswamy H. Sarma.

Theoretical evaluation of the malathion and its chemical derivatives interaction with cytosolic phospholipase A2 from zebrafish

Cytosolic phospholipase A2 (cPLA2) belongs to a large family of proteins and plays a crucial role in the regulation of arachidonic acid metabolism and inflammation cascade in zebrafish (Danio rerio). This enzyme with a molecular weight of 85 kDa, has two distinct domains. One is the regulatory and calcium-dependent (Ca²⁺) domain called C2, the other is the catalytic α/β hydrolase Ca²⁺-independent domain, where serine and aspartic acid catalytic dyad residues are present. We investigated the interaction of malathion and their organophosphate metabolites in the cPLA2 using in silico tools. Molecular docking results showed hydrophobic interactions with the paraoxon and catalytic site residue (Ser 223). Malathion increases intracellular Ca²⁺ due to endoplasmic reticulum influx which in turn activities phospholipase A2 and arachidonic acid release. Molecular docking and homology modelling of proteins and ligands could be a complementary tool for ecotoxicology and environment pollution assessment.

New putative therapeutic targets against Serratia marcescens using reverse vaccinology and subtractive genomics

The Gram-negative bacillus Serratia marcescens, a member of Enterobacteriaceae family, is an opportunistic nosocomial pathogen commonly found in hospital outbreaks that can cause infections in the urinary tract, bloodstream, central nervous system and pneumonia. Because S. marcescens strains are resistant to several antibiotics, it is critical the need for effective treatments, including new drugs and vaccines. Here, we applied reverse vaccinology and subtractive genomic approaches for the in silico prediction of potential vaccine and drug targets against 59 strains of S. marcescens. We found 759 core non-host homologous proteins, of which 87 are putative surface-exposed proteins, 183 secreted proteins, and 80 membrane proteins. From these proteins, we predicted seven candidates vaccine targets: a sn-glycerol-3-phosphate-binding periplasmic protein UgpB, a vitamin B12 TonB-dependent receptor, a ferrichrome porin FhuA, a divisome-associated lipoprotein YraP, a membrane-bound lytic murein transglycosylase A, a peptidoglycan lytic exotransglycosylase, and a DUF481 domain-containing protein. We also predicted two drug targets: a N(4)-acetylcytidine amidohydrolase, and a DUF1428 family protein. Using the molecular docking approach for each drug target, we identified and selected ZINC04259491 and ZINC04235390 molecules as the most favorable interactions with the target active site residues. Our findings may contribute to the development of vaccines and new drug targets against S. marcescens. Communicated by Ramaswamy H. Sarma.

Comparative genomics of Bordetella pertussis and prediction of new vaccines and drug targets

Pertussis is a highly contagious respiratory disease caused by Bordetella pertussis, a Gram-negative bacterium described over a century ago. Despite broad vaccine coverage and treatment options, the disease is remerging as a public health problem especially in infants and older children. Recent data indicate re-emergence of the disease is related to bacterial resistance to immune defences and decreased vaccine effectiveness, which obviously suggests the need of new effective vaccines and drugs. In an attempt to contribute with solutions to this great challenge, bioinformatics tools were used to genetically comprehend the species of these bacteria and predict new vaccines and drug targets. In fact, approaches were used to analysis genomic plasticity, gene synteny and species similarities between the 20 genomes of Bordetella pertussis already available. Furthermore, it was conducted reverse vaccinology and docking analysis to identify proteins with potential to become vaccine and drug targets, respectively. The analyses showed the 20 genomes belongs to a homogeneous group that has preserved most of the genes over time. Besides that, were found genomics islands and good proteins to be candidates for vaccine and drugs. Taken together, these results suggests new possibilities that may be useful to develop new vaccines and drugs that will help the prevention and treatment strategies of pertussis disease caused by these Bordetella strains. Communicated by Ramaswamy H. Sarma.

RE, Marín-Villa M. In Silico, In Vitro, and Pharmacokinetic Studies of UBMC-4, a Potential Novel Compound for Treating against Trypanosoma cruzi

The lack of therapeutic alternatives for the treatment of Chagas disease, a neglected disease, drives the discovery of new drugs with trypanocidal activity. Consequently, we conducted in vitro studies using UBMC-4, a potential Trypanosoma cruzi AKT-like pleckstrin homology (PH) domain inhibitory compound found using bioinformatics tools. The half effective concentration (EC50) on intracellular amastigotes was determined at 1.85 ± 1 μM showing low cytotoxicity (LC50) > 40 μM on human cell lines tested. In order to study the lethal effect caused by the compound on epimastigotes, morphological changes were assessed by scanning and transmission electron microscopy. Progressive alterations such as flagellum inactivation, cell size reduction, nuclear structure alteration, condensation of chromatin towards the nuclear periphery, vacuole formation, and mitochondrial swelling with kinetoplast integrity loss were evidenced. In addition, apoptosis-like markers in T. cruzi were assessed by flow cytometry, demonstrating that the effect of UBMC-4 on T. cruzi AKT-like kinase reduced the tolerance to nutritional stress-triggered, apoptosis-like events, including DNA fragmentation, mitochondrial damage, and loss of plasma membrane integrity. After this, UBMC-4 was formulated for oral administration and pharmacokinetics were analyzed in a mouse model. Finally, upon oral administration of 200 mg/kg in mice, we found that a UBMC-4 plasma concentration remaining in circulation beyond 24 h after administration is well described by the two-compartment model. We conclude that UBMC-4 has an effective trypanocidal activity in vitro at low concentrations and this effect is evident in T. cruzi cell structures. In mice, UBMC-4 was well absorbed and reached plasma concentrations higher than the EC50, showing features that would aid in developing a new drug to treat Chagas disease.

Computational identification of putative common genomic drug and vaccine targets in Mycoplasma genitalium

Mycoplasma genitalium is an obligate intracellular bacterium that is responsible for several sexually transmitted infections, including non-gonococcal urethritis in men and several inflammatory reproductive tract syndromes in women. Here, we applied subtractive genomics and reverse vaccinology approaches for in silico prediction of potential vaccine and drug targets against five strains of M. genitalium. We identified 403 genes shared by all five strains, from which 104 non-host homologous proteins were selected, comprising of 44 exposed/secreted/membrane proteins and 60 cytoplasmic proteins. Based on the essentiality, functionality, and structure-based binding affinity, we finally predicted 19 (14 novel) putative vaccine and 7 (2 novel) candidate drug targets. The docking analysis showed six molecules from the ZINC database as promising drug candidates against the identified targets. Altogether, both vaccine candidates and drug targets identified here may contribute to the future development of therapeutic strategies to control the spread of M. genitalium worldwide.

Pan-genomic analyses of 47 complete genomes of the Rickettsia genus and prediction of new vaccine targets and virulence factors of the species

The genus Rickettsia belongs to the Proteobacteria phylum and these bacteria infect animals and humans causing a range of diseases worldwide. The genus is divided into 4 groups and despite the public health threat and the knowledge accumulated so far, the mandatory intracellular bacteria behaviour and limitation for in vitro culture makes it difficult to create new vaccines and drug targets to these bacteria. In an attempt to overcome these limitations, pan-genomic approaches has used 47 genomes of the genus Rickettsia, in order to describe species similarities and genomics islands. Moreover, we conducted reverse vaccinology and docking analysis aiming the identification of proteins that have great potential to become vaccine and drug targets. We found out that the bacteria of the four Rickettsia groups have a high similarity with each other, with about 90 to 100% of identity. A pathogenicity island and a resistance island were predicted. In addition, 8 proteins were also predicted as strong candidates for vaccine and 9 as candidates for drug targets. The prediction of the proteins leads us to believe in a possibility of prospecting potential drugs or creating a polyvalent vaccine, which could reach most strains of this large group of bacteria.Communicated by Ramaswamy H. Sarma.

NMR structure and dynamics of Q4DY78, a conserved kinetoplasid-specific protein from Trypanosoma cruzi

The 106-residue protein Q4DY78 (UniProt accession number) from Trypanosoma cruzi is highly conserved in the related kinetoplastid pathogens Trypanosoma brucei and Leishmania major. Given the essentiality of its orthologue in T. brucei, the high sequence conservation with other trypanosomatid proteins, and the low sequence similarity with mammalian proteins, Q4DY78 is an attractive protein for structural characterization. Here, we solved the structure of Q4DY78 by solution NMR and evaluated its backbone dynamics. Q4DY78 is composed of five α -helices and a small, two-stranded antiparallel β-sheet. The backbone RMSD is 0.22 ± 0.05 Å for the representative ensemble of the 20 lowest-energy structures. Q4DY78 is overall rigid, except for N-terminal residues (V⁸ to I¹⁰), residues at loop 4 (K⁵⁷ to G⁶⁵) and residues at the C-terminus (F⁸⁹ to F¹¹²). Q4DY78 has a short motif FPCAP that could potentially mediate interactions with the host cytoskeleton via interaction with EVH1 (Drosophila Enabled (Ena)/Vasodilator-stimulated phosphoprotein (VASP) homology 1) domains. Albeit Q4DY78 lacks calcium-binding motifs, its fold resembles that of eukaryotic calcium-binding proteins such as calcitracin, calmodulin, and polcacin Bet V4. We characterized this novel protein with a calcium binding fold without the capacity to bind calcium.

Evolution and Diversity of Semaphorins and Plexins in Choanoflagellates

Semaphorins and plexins are cell surface ligand/receptor proteins that affect cytoskeletal dynamics in metazoan cells. Interestingly, they are also present in Choanoflagellata, a class of unicellular heterotrophic flagellates that forms the phylogenetic sister group to Metazoa. Several members of choanoflagellates are capable of forming transient colonies, whereas others reside solitary inside exoskeletons; their molecular diversity is only beginning to emerge. Here, we surveyed genomics data from 22 choanoflagellate species and detected semaphorin/plexin pairs in 16 species. Choanoflagellate semaphorins (Sema-FN1) contain several domain features distinct from metazoan semaphorins, including an N-terminal Reeler domain that may facilitate dimer stabilization, an array of fibronectin type III domains, a variable serine/threonine-rich domain that is a potential site for O-linked glycosylation, and a SEA domain that can undergo autoproteolysis. In contrast, choanoflagellate plexins (Plexin-1) harbor a domain arrangement that is largely identical to metazoan plexins. Both Sema-FN1 and Plexin-1 also contain a short homologous motif near the C-terminus, likely associated with a shared function. Three-dimensional molecular models revealed a highly conserved structural architecture of choanoflagellate Plexin-1 as compared to metazoan plexins, including similar predicted conformational changes in a segment that is involved in the activation of the intracellular Ras-GAP domain. The absence of semaphorins and plexins in several choanoflagellate species did not appear to correlate with unicellular versus colonial lifestyle or ecological factors such as fresh versus salt water environment. Together, our findings support a conserved mechanism of semaphorin/plexin proteins in regulating cytoskeletal dynamics in unicellular and multicellular organisms.

In Silico Identification of New Targets for Diagnosis, Vaccine, and Drug Candidates against Trypanosoma cruzi

Chagas disease is a neglected tropical disease caused by the parasite Trypanosoma cruzi. Despite the efforts and distinct methodologies, the search of antigens for diagnosis, vaccine, and drug targets for the disease is still needed. The present study is aimed at identifying possible antigens that could be used for diagnosis, vaccine, and drugs targets against T. cruzi using reverse vaccinology and molecular docking. The genomes of 28 T. cruzi strains available in GenBank (NCBI) were used to obtain the genomic core. Then, subtractive genomics was carried out to identify nonhomologous genes to the host in the core. A total of 2630 conserved proteins in 28 strains of T. cruzi were predicted using OrthoFinder and Diamond software, in which 515 showed no homology to the human host. These proteins were evaluated for their subcellular localization, from which 214 are cytoplasmic and 117 are secreted or present in the plasma membrane. To identify the antigens for diagnosis and vaccine targets, we used the VaxiJen software, and 14 nonhomologous proteins were selected showing high binding efficiency with MHC I and MHC II with potential for in vitro and in vivo tests. When these 14 nonhomologous molecules were compared against other trypanosomatids, it was found that the retrotransposon hot spot (RHS) protein is specific only for T. cruzi parasite suggesting that it could be used for Chagas diagnosis. Such 14 proteins were analyzed using the IEDB software to predict their epitopes in both B and T lymphocytes. Furthermore, molecular docking analysis was performed using the software MHOLline. As a result, we identified 6 possible T. cruzi drug targets that could interact with 4 compounds already known as antiparasitic activities. These 14 protein targets, along with 6 potential drug candidates, can be further validated in future studies, in vivo, regarding Chagas disease.

Design, synthesis, and evaluation of substrate - analogue inhibitors of Trypanosoma cruzi ribose 5-phosphate isomerase type B

Ribose 5-phosphate isomerase type B (RPI-B) is a key enzyme of the pentose phosphate pathway that catalyzes the isomerization of ribose 5-phosphate (R5P) and ribulose 5-phosphate (Ru5P). Trypanosoma cruzi RPI-B (TcRPI-B) appears to be a suitable drug-target mainly due to: (i) its essentiality (as previously shown in other trypanosomatids), (ii) it does not present a homologue in mammalian genomes sequenced thus far, and (iii) it participates in the production of NADPH and nucleotide/nucleic acid synthesis that are critical for parasite cell survival. In this survey, we report on the competitive inhibition of TcRPI-B by a substrate - analogue inhibitor, Compound B (K_i = 5.5 ± 0.1 μM), by the Dixon method. This compound has an iodoacetamide moiety that is susceptible to nucleophilic attack, particularly by the cysteine thiol group. Compound B was conceived to specifically target Cys-69, an important active site residue. By incubating TcRPI-B with Compound B, a trypsin digestion LC-MS/MS analysis revealed the identification of Compound B covalently bound to Cys-69. This inhibitor also exhibited notable in vitro trypanocidal activity against T. cruzi infective life-stages co-cultured in NIH-3T3 murine host cells (IC₅₀ = 17.40 ± 1.055 μM). The study of Compound B served as a proof-of-concept so that next generation inhibitors can potentially be developed with a focus on using a prodrug group in replacement of the iodoacetamide moiety, thus representing an attractive starting point for the future treatment of Chagas' disease.

In silico Prediction of New Drug Candidates Against the Multidrug-Resistant and Potentially Zoonotic Fish Pathogen Serotype III Streptococcus agalactiae

Streptococcus agalactiae is an invasive multi-host pathogen that causes invasive diseases mainly in newborns, elderly, and individuals with underlying health complications. In fish, S. agalactiae causes streptococcosis, which is characterized by septicemia and neurological signs, and leads to great economic losses to the fish farming industry worldwide. These bacteria can be classified into different serotypes based on capsular antigens, and into different sequence types (ST) based on multilocus sequence typing (MLST). In 2015, serotype III ST283 was identified to be associated with a foodborne invasive disease in non-pregnant immunocompetent humans in Singapore, and the infection was related to raw fish consumption. In addition, a serotype III strain isolated from tilapia in Brazil has been reported to be resistant to five antibiotic classes. This specific serotype can serve as a reservoir of resistance genes and pose a serious threat to public health. Thus, new approaches for the control and treatment of S. agalactiae infections are needed. In the present study, 24 S. agalactiae serotype III complete genomes, isolated from human and fish hosts, were compared. The core genome was identified, and, using bioinformatics tools and subtractive criteria, five proteins were identified as potential drug targets. Furthermore, 5,008 drug-like natural compounds were virtually screened against the identified targets. The ligands with the best binding properties are suggested for further in vitro and in vivo analysis.

Trehalose synthesis inhibitor: A molecular in silico drug design

Infectious diseases are serious public health problems, affecting a large portion of the world's population. A molecule that plays a key role in pathogenic organisms is trehalose and recently has been an interest in the metabolism of this molecule for drug development. The trehalose-6-phosphate synthase (TPS1) is an enzyme responsible for the biosynthesis of trehalose-6-phosphate (T6P) in the TPS1/TPS2 pathway, which results in the formation of trehalose. Studies carried out by our group demonstrated the inhibitory capacity of T6P in the TPS1 enzyme from Saccharomyces cerevisiae, preventing the synthesis of trehalose. By in silico techniques, we compiled sequences and experimentally determined structures of TPS1. Sequence alignments and molecular modeling were performed. The generated structures were submitted in validation of algorithms, aligned structurally and analyzed evolutionarily. Molecular docking methodology was applied to analyze the interaction between T6P and TPS1 and ADMET properties of T6P were analyzed. The results demonstrated the models created presented sequence and structural similarities with experimentally determined structures. With the molecular docking, a cavity in the protein surface was identified and the molecule T6P was interacting with the residues TYR-40, ALA-41, MET-42, and PHE-372, indicating the possible uncompetitive inhibition mechanism provided by this ligand, which can be useful in directing the molecular design of inhibitors. In ADMET analyses, T6P had acceptable risk values compared with other compounds from World Drug Index. Therefore, these results may present a promising strategy to explore to develop a broad-spectrum antibiotic of this specific target with selectivity, potency, and reduced side effects, leading to a new way to treat infectious diseases like tuberculosis and candidiasis.

Reverse vaccinology and subtractive genomics reveal new therapeutic targets against Mycoplasma pneumoniae: a causative agent of pneumonia

Pneumonia is an infectious disease caused by bacteria, viruses or fungi that results in millions of deaths globally. Despite the existence of prophylactic methods against some of the major pathogens of the disease, there is no efficient prophylaxis against atypical agents such as Mycoplasma pneumoniae, a bacterium associated with cases of community-acquired pneumonia. Because of the morphological peculiarity of M. pneumoniae, which leads to an increased resistance to antibiotics, studies that prospectively investigate the development of vaccines and drug targets appear to be one of the best ways forward. Hence, in this paper, bioinformatics tools were used for vaccine and pharmacological prediction. We conducted comparative genomic analysis on the genomes of 88 M. pneumoniae strains, as opposed to a reverse vaccinology analysis, in relation to the capacity of M. pneumoniae proteins to bind to the major histocompatibility complex, revealing seven targets with immunogenic potential. Predictive cytoplasmic proteins were tested as potential drug targets by studying their structures in relation to other proteins, metabolic pathways and molecular anchorage, which identified five possible drug targets. These findings are a valuable addition to the development of vaccines and the selection of new in vivo drug targets that may contribute to further elucidating the molecular basis of M. pneumoniae-host interactions.

Studying effects of different protonation states of His11 and His102 in ribose-5-phosphate isomerase of Trypanosoma cruzi: an example of cooperative behavior

The Trypanosoma cruzi ribose-5-phosphate isomerase B (TcRpiB) is a crucial piece in the pentose phosphate pathway and thus is a potential drug target for treatment of Chagas' disease. TcRpiB residues, such as Cys69, Asp45, Glu149 and Pro47, have confirmed their roles in substrate recognition, catalytic reaction and binding site conformation. However, the joint performance of His11 and His102, in the D-ribose-5-phosphate (R5P) in the catalysis is not well understood. In this work, we probed the influence of different protonation states of His11 and His102 on the behavior of the ligand R5P using molecular dynamics simulations, network analysis and thermodynamic integration. Simulations revealed that a protonated His11 combined with a neutral His102 (His11⁺‒His102) was able to stabilize the ligand R5P in the binding site. Moreover, calculated relative free energy differences showed that when protonated His11 was coupled to a neutral His102 an exergonic process takes place. On the other hand, neutral His11 combined with a protonated His102 (His11‒His102⁺), sampled conformations that resembled the catalyzed product D-ribulose-5-phosphate (Ru5P). Network analysis also demonstrated some peculiarities for these systems with some negatively correlated nodes in the binding site for His11‒His102⁺, and exclusive suboptimal paths for His11⁺‒His102. Therefore, the combined approach presented in this paper proposes two suitable protonation states for the TcRpiB catalytic mechanism, where an extra proton in either histidines might favor R5P binding or influence isomerization reaction to Ru5P. Our results may guide further in silico drug discovery studies. Communicated by Ramaswamy H. Sarma.

Insights about resveratrol analogs against trypanothione reductase of Leishmania braziliensis: Molecular modeling, computational docking and in vitro antileishmanial studies

In this work, we combined molecular modeling, computational docking and in vitro analysis to explore the antileishmanial effect of some resveratrol analogs (ResAn), focusing on their pro-oxidant effect. The molecular target was the trypanothione reductase of Leishmania braziliensis (LbTryR), an essential component of the antioxidant defenses in trypanosomatid parasites. Three-dimensional structures of LbTryR were modeled and molecular docking studies of ResAn1-5 compounds showed the following affinity: ResAn1 > ResAn2 > ResAn4 > ResAn5 > ResAn3. Positive correlation was observed between these compounds' affinity to the LbTryR and the IC₅₀ values against Leishmania sp (ResAn1 < ResAn2 < ResAn4), which allows for TryR being considered an important target for them. As the compound ResAn1 showed the best antileishmanial activity, and docking studies showed its high affinity for NADP binding site (NS) of TryR, plus having been able to induce ROS production in L. braziliensis promastigotes treated, ResAn1 probably occupies NS interfering in the electron transfer processes responsible for the catalytic reaction. The in silico prediction of ADMET properties suggests that ResAn1 may be a promising drug candidate with properties to cross biological membranes and high gastrointestinal absorption, not violating Lipinski's rules. Ultimately, the antileishmanial effect of ResAn can be associated with a pro-oxidant effect which, in turn, can be exploited as an antimicrobial agent. Communicated by Ramaswamy H. Sarma.

Mining of potential drug targets through the identification of essential and analogous enzymes in the genomes of pathogens of Glycine max, Zea mays and Solanum lycopersicum

Pesticides are one of the most widely used pest and disease control measures in plant crops and their indiscriminate use poses a direct risk to the health of populations and environment around the world. As a result, there is a great need for the development of new, less toxic molecules to be employed against plant pathogens. In this work, we employed an in silico approach to study the genes coding for enzymes of the genomes of three commercially important plants, soybean (Glycine max), tomato (Solanum lycopersicum) and corn (Zea mays), as well as 15 plant pathogens (4 bacteria and 11 fungi), focusing on revealing a set of essential and non-homologous isofunctional enzymes (NISEs) that could be prioritized as drug targets. By combining sequence and structural data, we obtained an initial set of 568 cases of analogy, of which 97 were validated and further refined, revealing a subset of 29 essential enzymatic activities with a total of 119 different structural forms, most belonging to central metabolic routes, including the carbohydrate metabolism, the metabolism of amino acids, among others. Further, another subset of 26 enzymatic activities possess a tertiary structure specific for the pathogen, not present in plants, men and Apis mellifera, which may be of importance for the development of specific enzymatic inhibitors against plant diseases that are less harmful to humans and the environment.

Putative vaccine candidates and drug targets identified by reverse vaccinology and subtractive genomics approaches to control Haemophilus ducreyi, the causative agent of chancroid

Chancroid is a sexually transmitted infection (STI) caused by the Gram-negative bacterium Haemophilus ducreyi The control of chancroid is difficult and the only current available treatment is antibiotic therapy; however, antibiotic resistance has been reported in endemic areas. Owing to recent outbreaks of STIs worldwide, it is important to keep searching for new treatment strategies and preventive measures. Here, we applied reverse vaccinology and subtractive genomic approaches for the in silico prediction of potential vaccine and drug targets against 28 strains of H. ducreyi We identified 847 non-host homologous proteins, being 332 exposed/secreted/membrane and 515 cytoplasmic proteins. We also checked their essentiality, functionality and virulence. Altogether, we predicted 13 candidate vaccine targets and three drug targets, where two vaccines (A01_1275, ABC transporter substrate-binding protein; and A01_0690, Probable transmembrane protein) and three drug targets (A01_0698, Purine nucleoside phosphorylase; A01_0702, Transcription termination factor; and A01_0677, Fructose-bisphosphate aldolase class II) are harboured by pathogenicity islands. Finally, we applied a molecular docking approach to analyse each drug target and selected ZINC77257029, ZINC43552589 and ZINC67912117 as promising molecules with favourable interactions with the target active site residues. Altogether, the targets identified here may be used in future strategies to control chancroid worldwide.

Building protein-protein interaction networks for Leishmania species through protein structural information

Background

Systematic analysis of a parasite interactome is a key approach to understand different biological processes. It makes possible to elucidate disease mechanisms, to predict protein functions and to select promising targets for drug development. Currently, several approaches for protein interaction prediction for non-model species incorporate only small fractions of the entire proteomes and their interactions. Based on this perspective, this study presents an integration of computational methodologies, protein network predictions and comparative analysis of the protozoan species Leishmania braziliensis and Leishmania infantum. These parasites cause Leishmaniasis, a worldwide distributed and neglected disease, with limited treatment options using currently available drugs.

Results

The predicted interactions were obtained from a meta-approach, applying rigid body docking tests and template-based docking on protein structures predicted by different comparative modeling techniques. In addition, we trained a machine-learning algorithm (Gradient Boosting) using docking information performed on a curated set of positive and negative protein interaction data. Our final model obtained an AUC = 0.88, with recall = 0.69, specificity = 0.88 and precision = 0.83. Using this approach, it was possible to confidently predict 681 protein structures and 6198 protein interactions for L. braziliensis, and 708 protein structures and 7391 protein interactions for L. infantum. The predicted networks were integrated to protein interaction data already available, analyzed using several topological features and used to classify proteins as essential for network stability.

Conclusions

The present study allowed to demonstrate the importance of integrating different methodologies of interaction prediction to increase the coverage of the protein interaction of the studied protocols, besides it made available protein structures and interactions not previously reported.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2105-6) contains supplementary material, which is available to authorized users.

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.

Functional Analogy in Human Metabolism: Enzymes with Different Biological Roles or Functional Redundancy?

Since enzymes catalyze almost all chemical reactions that occur in living organisms, it is crucial that genes encoding such activities are correctly identified and functionally characterized. Several studies suggest that the fraction of enzymatic activities in which multiple events of independent origin have taken place during evolution is substantial. However, this topic is still poorly explored, and a comprehensive investigation of the occurrence, distribution, and implications of these events has not been done so far. Fundamental questions, such as how analogous enzymes originate, why so many events of independent origin have apparently occurred during evolution, and what are the reasons for the coexistence in the same organism of distinct enzymatic forms catalyzing the same reaction, remain unanswered. Also, several isofunctional enzymes are still not recognized as nonhomologous, even with substantial evidence indicating different evolutionary histories. In this work, we begin to investigate the biological significance of the cooccurrence of nonhomologous isofunctional enzymes in human metabolism, characterizing functional analogous enzymes identified in metabolic pathways annotated in the human genome. Our hypothesis is that the coexistence of multiple enzymatic forms might not be interpreted as functional redundancy. Instead, these enzymatic forms may be implicated in distinct (and probably relevant) biological roles.

An In Silico Identification of Common Putative Vaccine Candidates against Treponema pallidum: A Reverse Vaccinology and Subtractive Genomics Based Approach

Sexually transmitted infections (STIs) are caused by a wide variety of bacteria, viruses, and parasites that are transmitted from one person to another primarily by vaginal, anal, or oral sexual contact. Syphilis is a serious disease caused by a sexually transmitted infection. Syphilis is caused by the bacterium Treponema pallidum subspecies pallidum. Treponema pallidum (T. pallidum) is a motile, gram-negative spirochete, which can be transmitted both sexually and from mother to child, and can invade virtually any organ or structure in the human body. The current worldwide prevalence of syphilis emphasizes the need for continued preventive measures and strategies. Unfortunately, effective measures are limited. In this study, we focus on the identification of vaccine targets and putative drugs against syphilis disease using reverse vaccinology and subtractive genomics. We compared 13 strains of T. pallidum using T. pallidum Nichols as the reference genome. Using an in silicoapproach, four pathogenic islands were detected in the genome of T. pallidum Nichols. We identified 15 putative antigenic proteins and sixdrug targets through reverse vaccinology and subtractive genomics, respectively, which can be used as candidate therapeutic targets in the future.

Amyotrophic Lateral Sclerosis Type 20 - In Silico Analysis and Molecular Dynamics Simulation of hnRNPA1

Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disease that affects the upper and lower motor neurons. 5-10% of cases are genetically inherited, including ALS type 20, which is caused by mutations in the hnRNPA1 gene. The goals of this work are to analyze the effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on hnRNPA1 protein function, to model the complete tridimensional structure of the protein using computational methods and to assess structural and functional differences between the wild type and its variants through Molecular Dynamics simulations. nsSNP, PhD-SNP, Polyphen2, SIFT, SNAP, SNPs&GO, SNPeffect and PROVEAN were used to predict the functional effects of nsSNPs. Ab initio modeling of hnRNPA1 was made using Rosetta and refined using KoBaMIN. The structure was validated by PROCHECK, Rampage, ERRAT, Verify3D, ProSA and Qmean. TM-align was used for the structural alignment. FoldIndex, DICHOT, ELM, D2P2, Disopred and DisEMBL were used to predict disordered regions within the protein. Amino acid conservation analysis was assessed by Consurf, and the molecular dynamics simulations were performed using GROMACS. Mutations D314V and D314N were predicted to increase amyloid propensity, and predicted as deleterious by at least three algorithms, while mutation N73S was predicted as neutral by all the algorithms. D314N and D314V occur in a highly conserved amino acid. The Molecular Dynamics results indicate that all mutations increase protein stability when compared to the wild type. Mutants D314N and N319S showed higher overall dimensions and accessible surface when compared to the wild type. The flexibility level of the C-terminal residues of hnRNPA1 is affected by all mutations, which may affect protein function, especially regarding the protein ability to interact with other proteins.

First co-expression of a lipase and its specific foldase obtained by metagenomics

Background

Metagenomics is a useful tool in the search for new lipases that might have characteristics that make them suitable for application in biocatalysis. This paper reports the cloning, co-expression, purification and characterization of a new lipase, denominated LipG9, and its specific foldase, LifG9, from a metagenomic library derived from a fat-contaminated soil.

Results

Within the metagenomic library, the gene lipg9 was cloned jointly with the gene of the foldase, lifg9. LipG9 and LifG9 have 96% and 84% identity, respectively, with the corresponding proteins of Aeromonas veronii B565. LipG9 and LifG9 were co-expressed, both in N-truncated form, in Escherichia coli BL21(DE3), using the vectors pET28a(+) and pT7-7, respectively, and then purified by affinity chromatography using a Ni²⁺ column (HiTrap Chelating HP). The purified enzyme eluted from the column complexed with its foldase. The molecular masses of the N-truncated proteins were 32 kDa for LipG9, including the N-terminal His-tag with 6 residues, and 23 kDa for LifG9, which did not have a His-tag. The biochemical and kinetic characteristics of the purified lipase-foldase preparation were investigated. This preparation was active and stable over a wide range of pH values (6.5-9.5) and temperatures (10-40°C), with the highest specific activity, of 1500 U mg⁻¹, being obtained at pH 7.5 at 30°C. It also had high specific activities against tributyrin, tricaprylin and triolein, with values of 1852, 1566 and 817 U mg⁻¹, respectively. A phylogenetic analysis placed LipG9 in the lipase subfamily I.1. A comparison of the sequence of LipG9 with those of other bacterial lipases in the Protein Data Bank showed that LipG9 contains not only the classic catalytic triad (Ser¹⁰³, Asp²⁵⁰, His²⁷²), with the catalytic Ser occurring within a conserved pentapeptide, Gly-His-Ser-His-Gly, but also a conserved disulfide bridge and a conserved calcium binding site. The homology-modeled structure presents a canonical α/β hydrolase folding type I.

Conclusions

This paper is the first to report the successful co-expression of a lipase and its associated foldase from a metagenomic library. The high activity and stability of Lip-LifG9 suggest that it has a good potential for use in biocatalysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-014-0171-7) contains supplementary material, which is available to authorized users.

Proteome scale comparative modeling for conserved drug and vaccine targets identification in Corynebacterium pseudotuberculosis

Corynebacterium pseudotuberculosis (Cp) is a pathogenic bacterium that causes caseous lymphadenitis (CLA), ulcerative lymphangitis, mastitis, and edematous to a broad spectrum of hosts, including ruminants, thereby threatening economic and dairy industries worldwide. Currently there is no effective drug or vaccine available against Cp. To identify new targets, we adopted a novel integrative strategy, which began with the prediction of the modelome (tridimensional protein structures for the proteome of an organism, generated through comparative modeling) for 15 previously sequenced C. pseudotuberculosis strains. This pan-modelomics approach identified a set of 331 conserved proteins having 95-100% intra-species sequence similarity. Next, we combined subtractive proteomics and modelomics to reveal a set of 10 Cp proteins, which may be essential for the bacteria. Of these, 4 proteins (tcsR, mtrA, nrdI, and ispH) were essential and non-host homologs (considering man, horse, cow and sheep as hosts) and satisfied all criteria of being putative targets. Additionally, we subjected these 4 proteins to virtual screening of a drug-like compound library. In all cases, molecules predicted to form favorable interactions and which showed high complementarity to the target were found among the top ranking compounds. The remaining 6 essential proteins (adk, gapA, glyA, fumC, gnd, and aspA) have homologs in the host proteomes. Their active site cavities were compared to the respective cavities in host proteins. We propose that some of these proteins can be selectively targeted using structure-based drug design approaches (SBDD). Our results facilitate the selection of C. pseudotuberculosis putative proteins for developing broad-spectrum novel drugs and vaccines. A few of the targets identified here have been validated in other microorganisms, suggesting that our modelome strategy is effective and can also be applicable to other pathogens.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the “Mexicana complex”, reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.

Structural and functional analysis of human SOD1 in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with familial inheritance (fALS) in 5% to 10% of cases; 25% of those are caused by mutations in the superoxide dismutase 1 (SOD1) protein. More than 100 mutations in the SOD1 gene have been associated with fALS, altering the geometry of the active site, protein folding and the interaction between monomers. We performed a functional analysis of non-synonymous single nucleotide polymorphisms (nsSNPs) in 124 fALS SOD1 mutants. Eleven different algorithms were used to estimate the functional impact of the replacement of one amino acid on protein structure: SNPs&GO, PolyPhen-2, SNAP, PMUT, Sift, PhD-SNP, nsSNPAnalyzer, TANGO, WALTZ, LIMBO and FoldX. For the structural analysis, theoretical models of 124 SNPs of SOD1 were created by comparative modeling using the MHOLline workflow, which includes Modeller and Procheck. Models were aligned with the native protein by the TM-align algorithm. A human-curated database was developed using the server side include in Java, JMOL. The results of this functional analysis indicate that the majority of the 124 natural mutants are harmful to the protein structure and thus corroborate the correlation between the reported mutations and fALS. In the structural analysis, all models showed conformational changes when compared to wild-type SOD1, and the degree of structural alignment varied between them. The SOD1 database converge structural and functional analyses of SOD1; it is a vast resource for the molecular analysis of amyotrophic lateral sclerosis, which allows the user to expand his knowledge on the molecular basis of the disease. The SOD1 database is available at http://bioinfogroup.com/database.

Computational models for neglected diseases: gaps and opportunities

Neglected diseases, such as Chagas disease, African sleeping sickness, and intestinal worms, affect millions of the world's poor. They disproportionately affect marginalized populations, lack effective treatments or vaccines, or existing products are not accessible to the populations affected. Computational approaches have been used across many of these diseases for various aspects of research or development, and yet data produced by computational approaches are not integrated and widely accessible to others. Here, we identify gaps in which computational approaches have been used for some neglected diseases and not others. We also make recommendations for the broad-spectrum integration of these techniques into a neglected disease drug discovery and development workflow.

Structural modeling and in silico analysis of human superoxide dismutase 2

Aging in the world population has increased every year. Superoxide dismutase 2 (Mn-SOD or SOD2) protects against oxidative stress, a main factor influencing cellular longevity. Polymorphisms in SOD2 have been associated with the development of neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, as well as psychiatric disorders, such as schizophrenia, depression and bipolar disorder. In this study, all of the described natural variants (S10I, A16V, E66V, G76R, I82T and R156W) of SOD2 were subjected to in silico analysis using eight different algorithms: SNPeffect, PolyPhen-2, PhD-SNP, PMUT, SIFT, SNAP, SNPs&GO and nsSNPAnalyzer. This analysis revealed disparate results for a few of the algorithms. The results showed that, from at least one algorithm, each amino acid substitution appears to harmfully affect the protein. Structural theoretical models were created for variants through comparative modelling performed using the MHOLline server (which includes MODELLER and PROCHECK) and ab initio modelling, using the I-Tasser server. The predicted models were evaluated using TM-align, and the results show that the models were constructed with high accuracy. The RMSD values of the modelled mutants indicated likely pathogenicity for all missense mutations. Structural phylogenetic analysis using ConSurf revealed that human SOD2 is highly conserved. As a result, a human-curated database was generated that enables biologists and clinicians to explore SOD2 nsSNPs, including predictions of their effects and visualisation of the alignment of both the wild-type and mutant structures. The database is freely available at http://bioinfogroup.com/database and will be regularly updated.

Specific and nonhomologous isofunctional enzymes of the genetic information processing pathways as potential therapeutical targets for tritryps

Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi (Tritryps) are unicellular protozoa that cause leishmaniasis, sleeping sickness and Chagas' disease, respectively. Most drugs against them were discovered through the screening of large numbers of compounds against whole parasites. Nonhomologous isofunctional enzymes (NISEs) may present good opportunities for the identification of new putative drug targets because, though sharing the same enzymatic activity, they possess different three-dimensional structures thus allowing the development of molecules against one or other isoform. From public data of the Tritryps' genomes, we reconstructed the Genetic Information Processing Pathways (GIPPs). We then used AnEnPi to look for the presence of these enzymes between Homo sapiens and Tritryps, as well as specific enzymes of the parasites. We identified three candidates (ECs 3.1.11.2 and 6.1.1.-) in these pathways that may be further studied as new therapeutic targets for drug development against these parasites.