How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

Repositioning FDA-Approved Drug Against Chagas Disease and Cutaneous Leishmaniosis by Structure-Based Virtual Screening

BACKGROUND

Chagas disease and cutaneous leishmaniasis, two parasitic diseases caused by Trypanosoma cruzi (T. cruzi) and Leishmania mexicana (L. mexicana), respectively, have a major global impact. Current pharmacological treatments for these diseases are limited and can cause severe side effects; thus, there is a need for new antiprotozoal drugs.

METHODS

Using molecular docking, this work describes a structure-based virtual screening of an FDA-approved drug library against Trypanosoma cruzi and Leishmania mexicana glycolytic enzyme triosephosphate isomerase (TIM), which is highly conserved in these parasites. The selected compounds with potential dual inhibitory activity were tested in vitro to confirm their biological activity.

RESULTS

The study showed that five compounds: nilotinib, chlorhexidine, protriptyline, cyproheptadine, and montelukast, were more active against T. cruzi, than the reference drugs, nifurtimox and benznidazole while chlorhexidine and protriptyline were the most active against L. mexicana.

CONCLUSIONS

The analysis of these compounds and their structural characteristics may provide the basis for the development of new antiprotozoal agents.

Tunnel-like region observed as a potential allosteric site in Staphylococcus aureus Glyceraldehyde-3-phosphate dehydrogenase

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been investigated for allosteric features that might be used as potential target for specific inhibition of Staphylococcus aureus (S.aureus). X-ray structure of bacterial enzyme for which a tunnel-like opening passing through the center previously proposed as an allosteric site has been subjected to six independent 500 ns long Molecular Dynamics simulations. Harmonic bond restraints were employed at key residues to underline the allosteric feature of this region. A noticeable reduction was observed in the mobility of NAD+ binding domains when restrictions were applied. Also, a substantial decrease in cross-correlations between distant Cα fluctuations was detected throughout the structure. Mutual information (MI) analysis revealed a similar decrease in the degree of correspondence in positional fluctuations in all directions everywhere in the receptor. MI between backbone and side-chain torsional variations changed its distribution profile and decreased considerably around the catalytic sites when restraints were employed. Principal component analysis clearly showed that the restrained state sampled a narrower range of conformations than apo state, especially in the first principal mode due to restriction in the conformational flexibility of NAD+ binding domain. Clustering the trajectory based on catalytic site residues displayed a smaller repertoire of conformations for restrained state compared to apo. Representative snapshots subjected to k-shortest pathway analysis revealed the impact of bond restraints on the allosteric communication which displayed distinct optimal and suboptimal pathways for two states, where observed frequencies of critical residues Gln51 and Val283 at the proposed site changed considerably.

Repurposing of rabeprazole as an anti-Trypanosoma cruzi drug that targets cellular triosephosphate isomerase

Trypanosoma cruzi is the causative agent of American trypanosomiasis, which mainly affects populations in Latin America. Benznidazole is used to control the disease, with severe effects in patients receiving this chemotherapy. Previous studies have demonstrated the inhibition of triosephosphate isomerase from T. cruzi, but cellular enzyme inhibition has yet to be established. This study demonstrates that rabeprazole inhibits both cell viability and triosephosphate isomerase activity in T. cruzi epimastigotes. Our results show that rabeprazole has an IC50 of 0.4 µM, which is 14.5 times more effective than benznidazole. Additionally, we observed increased levels of methyl-glyoxal and advanced glycation end products after the inhibition of cellular triosephosphate isomerase by rabeprazole. Finally, we demonstrate that the inactivation mechanisms of rabeprazole on triosephosphate isomerase of T. cruzi can be achieved through the derivatization of three of its four cysteine residues. These results indicate that rabeprazole is a promising candidate against American trypanosomiasis.

Mechanism-based and computational modeling of hydrogen sulfide biogenesis inhibition: interfacial inhibition

Hydrogen sulfide (H₂S) is a gaseous signaling molecule that participates in various signaling functions in health and diseases. The tetrameric cystathionine γ-lyase (CSE) contributes to H₂S biogenesis and several investigations provide evidence on the pharmacological modulation of CSE as a potential target for the treatment of a multitude of conditions. D-penicillamine (D-pen) has recently been reported to selectively impede CSE-catalyzed H₂S production but the molecular bases for such inhibitory effect have not been investigated. In this study, we report that D-pen follows a mixed-inhibition mechanism to inhibit both cystathionine (CST) cleavage and H₂S biogenesis by human CSE. To decipher the molecular mechanisms underlying such a mixed inhibition, we performed docking and molecular dynamics (MD) simulations. Interestingly, MD analysis of CST binding reveals a likely active site configuration prior to gem-diamine intermediate formation, particularly H-bond formation between the amino group of the substrate and the O3' of PLP. Similar analyses realized with both CST and D-pen identified three potent interfacial ligand-binding sites for D-pen and offered a rational for D-pen effect. Thus, inhibitor binding not only induces the creation of an entirely new interacting network at the vicinity of the interface between enzyme subunits, but it also exerts long range effects by propagating to the active site. Overall, our study paves the way for the design of new allosteric interfacial inhibitory compounds that will specifically modulate H₂S biogenesis by cystathionine γ-lyase.

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC₅₀) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.

Altered Dynamics of S. aureus Phosphofructokinase via Bond Restraints at Two Distinct Allosteric Binding Sites

The effect of perturbation at the allosteric site was investigated through several replicas of molecular dynamics (MD) simulations conducted on bacterial phosphofructokinase (SaPFK). In our previous work, an alternative binding site was estimated to be allosteric in addition to the experimentally reported one. To highlight the effect of both allosteric sites on receptor's dynamics, MD runs were carried out on apo forms with and without perturbation. Perturbation was achieved via incorporating multiple bond restraints for residue pairs located at the allosteric site. Restraints applied to the predicted site caused one dimer to stiffen, whereas an increase in mobility was detected in the same dimer when the experimentally resolved site was restrained. Fluctuations in C_α-C_α distances which is used to disclose residues with high potential of communication indicated a marked increase in signal transmission within each dimer as the receptor switched to a restrained state. Cross-correlation of positional fluctuations indicated an overall decrease in the magnitude of both positive and negative correlations when restraints were employed on the predicted allosteric site whereas an exact opposite effect was observed for the reported site. Finally, mutual correspondence between positional fluctuations noticeably increased with restraints on predicted allosteric site, whereas an opposite effect was observed for restraints applied on experimentally reported one. In view of these findings, it is clear that the perturbation of either one of two allosteric sites effected the dynamics of the receptor with a distinct and contrasting character.

Precise druggability of the PTH type 1 receptor

Class B G protein-coupled receptors (GPCRs) are notoriously difficult to target by small molecules because their large orthosteric peptide-binding pocket embedded deep within the transmembrane domain limits the identification and development of nonpeptide small molecule ligands. Using the parathyroid hormone type 1 receptor (PTHR) as a prototypic class B GPCR target, and a combination of molecular dynamics simulations and elastic network model-based methods, we demonstrate that PTHR druggability can be effectively addressed. Here we found a key mechanical site that modulates the collective dynamics of the receptor and used this ensemble of PTHR conformers to identify selective small molecules with strong negative allosteric and biased properties for PTHR signaling in cell and PTH actions in vivo. This study provides a computational pipeline to detect precise druggable sites and identify allosteric modulators of PTHR signaling that could be extended to GPCRs to expedite discoveries of small molecules as novel therapeutic candidates.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Potential allosteric sites captured in glycolytic enzymes via residue-based network models: Phosphofructokinase, glyceraldehyde-3-phosphate dehydrogenase and pyruvate kinase

Likelihood of new allosteric sites for glycolytic enzymes, phosphofructokinase (PFK), glyceraldehyde-3-phosphate dehydrogenase (GADPH) and pyruvate kinase (PK) was evaluated for bacterial, parasitic and human species. Allosteric effect of a ligand binding at a site was revealed on the basis of low-frequency normal modes via C_α-harmonic residue network model. In bacterial PFK, perturbation of the proposed allosteric site outperformed the known allosteric one, producing a high amount of stabilization or reduced dynamics, on all catalytic regions. Another proposed allosteric spot at the dimer interface in parasitic PFK exhibited major stabilization effect on catalytic regions. In parasitic GADPH, the most desired allosteric response was observed upon perturbation of its tunnel region which incorporated key residues for functional regulation. Proposed allosteric site in bacterial PK produced a satisfactory allosteric response on all catalytic regions, whereas in human and parasitic PKs, a partial inhibition was observed. Residue network model based solely on contact topology identified the 'hub residues' with high betweenness tracing plausible allosteric communication pathways between distant functional sites. For both bacterial PFK and PK, proposed sites accommodated hub residues twice as much as the known allosteric site. Tunnel region in parasitic GADPH with the strongest allosteric effect among species, incorporated the highest number of hub residues. These results clearly suggest a one-to-one correspondence between the degree of allosteric effect and the number of hub residues in that perturbation site, which increases the likelihood of its allosteric nature.

Recent Advances in the Development of Triose Phosphate Isomerase Inhibitors as Antiprotozoal Agents

BACKGROUND

Parasitic diseases caused by protozoa such as Chagas disease, leishmaniasis, malaria, African trypanosomiasis, amebiasis, trichomoniasis, and giardiasis are considered serious public health problems in developing countries. Drug-resistance among parasites justifies the search for new therapeutic drugs and the identification of new targets becomes a valuable approach. In this scenario, glycolysis pathway which consists of the conversion of glucose into pyruvate plays an important role in the protozoa energy supply and it is therefore considered as a promising target. In this pathway, triose phosphate isomerase (TIM) plays an essential role in efficient energy production. Furthermore, protozoa TIM show structural differences with human enzyme counterparts suggesting the possibility of obtaining selective inhibitors. Therefore, TIM is considered a valid approach to develop new antiprotozoal agents, inhibiting the glycolysis in the parasite.

OBJECTIVE

In this review, we discuss the drug design strategies, structure-activity relationship, and binding modes of outstanding TIM inhibitors against Trypanosoma cruzi, Trypanosoma brucei, Plasmodium falciparum, Giardia lamblia, Leishmania mexicana, Trichomonas vaginalis, and Entamoeba histolytica.

RESULTS

TIM inhibitors showed mainly aromatic systems and symmetrical structure, where the size and type of heteroatom are important for enzyme inhibition. This inhibition is mainly based on the interaction with i) the interfacial region of TIM inducing changes on the quaternary and tertiary structure or ii) with the TIM catalytic region were the main pathways that disabled the catalytic activity of the enzyme.

CONCLUSION

Benzothiazole, benzoxazole, benzimidazole, and sulfhydryl derivatives stand out as TIM inhibitors. In silico and in vitro studies demonstrate that the inhibitors bind mainly at the TIM dimer interface. In this review, the development of new TIM inhibitors as antiprotozoal drugs is demonstrated as an important pharmaceutical strategy that may lead to new therapies for these ancient parasitic diseases.

Phenylbenzothiazole derivatives: effects against a Trypanosoma cruzi infection and toxicological profiles

Current treatments for Chagas disease have a limited impact during the chronic stage and trigger severe side effects. Treatments target Trypanosoma cruzi, the etiological agent of the disease. The aims of this study were to evaluate the trypanocidal activity of four 2-phenylbenzothiazole derivatives (BZT1-4) in vitro by using the infectious and non-infectious forms of T. cruzi (trypomastigotes and epimastigotes, respectively) and to test the most promising compound (BZT4) in vivo in mice. Additionally, the toxicological profile and possible neuronal damage were examined. In relation to trypomastigotes, BZT4 was more selective and effective than the reference drug (benznidazole) during this infective stage, apparently due to the synergistic action of the CF₃ and COOH substituents in the molecule. During the first few hours post-administration of BZT4, parasitemia decreased by 40% in an in vivo model of short-term treatment, but parasite levels later returned to the basal state. In the long-term assessment, the compound did not produce a significant antiparasitic effect, only attaining a 30% reduction in parasitemia by day 20 with the dose of 16 mg/kg. The toxicity test was based on repeated dosing of BZT4 (administered orally) during 21 days, which did not cause liver damage. However, the compound altered the concentration of proteins and the proteinic profile of neuronal cells in vitro, perhaps leading to an effect on the central nervous system. Further research on the low trypanocidal activity in vivo compared to the better in vitro effect could possibly facilitate molecular redesign to improve trypanocidal activity.

Virtual Screening of FDA-Approved Drugs against Triose Phosphate Isomerase from Entamoeba histolytica and Giardia lamblia Identifies Inhibitors of Their Trophozoite Growth Phase

Infectious diseases caused by intestinal protozoan, such as Entamoeba histolytica (E. histolytica) and Giardia lamblia (G. lamblia) are a worldwide public health issue. They affect more than 70 million people every year. They colonize intestines causing primarily diarrhea; nevertheless, these infections can lead to more serious complications. The treatment of choice, metronidazole, is in doubt due to adverse effects and resistance. Therefore, there is a need for new compounds against these parasites. In this work, a structure-based virtual screening of FDA-approved drugs was performed to identify compounds with antiprotozoal activity. The glycolytic enzyme triosephosphate isomerase, present in both E. histolytica and G. lamblia, was used as the drug target. The compounds with the best average docking score on both structures were selected for the in vitro evaluation. Three compounds, chlorhexidine, tolcapone, and imatinib, were capable of inhibit growth on G. lamblia trophozoites (0.05–4.935 μg/mL), while folic acid showed activity against E. histolytica (0.186 μg/mL) and G. lamblia (5.342 μg/mL).

Identification of Alternative Allosteric Sites in Glycolytic Enzymes for Potential Use as Species-Specific Drug Targets

Three allosteric glycolytic enzymes, phosphofructokinase, glyceraldehyde-3 phosphate dehydrogenase and pyruvate kinase, associated with bacterial, parasitic and human species, were explored to identify potential allosteric sites that would be used as prime targets for species-specific drug design purposes using a newly developed approach which incorporates solvent mapping, elastic network modeling, sequence and structural alignments. The majority of binding sites detected by solvent mapping overlapped with the interface regions connecting the subunits, thus appeared as promising target sites for allosteric regulation. Each binding site was then evaluated by its ability to alter the global dynamics of the receptor defined by the percentage change in the frequencies of the lowest-frequency modes most significantly and as anticipated, the most effective ones were detected in the vicinity of the well-reported catalytic and allosteric sites. Furthermore, some of our proposed regions intersected with experimentally resolved sites which are known to be critical for activity regulation, which further validated our approach. Despite the high degree of structural conservation encountered between bacterial/parasitic and human glycolytic enzymes, the majority of the newly presented allosteric sites exhibited a low degree of sequence conservation which further increased their likelihood to be used as species-specific target regions for drug design studies.

Essential site scanning analysis: A new approach for detecting sites that modulate the dispersion of protein global motions

Despite the wealth of methods developed for exploring the molecular basis of allostery in biomolecular systems, there is still a need for structure-based predictive tools that can efficiently detect susceptible sites for triggering allosteric responses. Toward this goal, we introduce here an elastic network model (ENM)-based method, Essential Site Scanning Analysis (ESSA). Essential sites are here defined as residues that would significantly alter the protein's global dynamics if bound to a ligand. To mimic the crowding induced upon substrate binding, the heavy atoms of each residue are incorporated as additional network nodes into the α-carbon-based ENM, and the resulting shifts in soft mode frequencies are used as a metric for evaluating the essentiality of each residue. Results on a dataset of monomeric proteins indicate the enrichment of allosteric and orthosteric binding sites, as well as global hinge regions among essential residues, highlighting the significant role of these sites in controlling the overall structural dynamics. Further integration of ESSA with information on predicted pockets and their local hydrophobicity density enables successful predictions of allosteric pockets for both ligand-bound and -unbound structures. ESSA can be efficiently applied to large multimeric systems. Three case studies, namely (i) G-protein binding to a GPCR, (ii) heterotrimeric assembly of the Ser/Thr protein phosphatase PP2A, and (iii) allo-targeting of AMPA receptor, demonstrate the utility of ESSA for identifying essential sites and narrowing down target allosteric sites identified by druggability simulations.

Intrinsic dynamics is evolutionarily optimized to enable allosteric behavior

Allosteric behavior is central to the function of many proteins, enabling molecular machinery, metabolism, signaling and regulation. Recent years have shown that the intrinsic dynamics of allosteric proteins defined by their 3-dimensional architecture or by the topology of inter-residue contacts favors cooperative motions that bear close similarity to structural changes they undergo during their allosteric actions. These conformational motions are usually driven by energetically favorable or soft modes at the low frequency end of the mode spectrum, and they are evolutionarily conserved among orthologs. These observations brought into light evolutionary adaptation mechanisms that help maintain, optimize or regulate allosteric behavior as the evolution from bacterial to higher organisms introduces sequential heterogeneities and structural complexities.

Protein-Ligand Complexes as Constrained Dynamical Systems

This study focuses on how the low-frequency end of the vibrational spectrum related to the functional motions changes as a protein binds to a small ligand(s). Our recently proposed residue-specific (RESPEC) elastic network model provides a natural laboratory for this aim due to its systematic mixed coarse-graining approach and parametrization. Current analysis on a large data set of protein-ligand complexes reveals a universal curve enclosing the frequency distributions, which bears the features of previous computational and experimental studies. We mostly observe positive frequency shifts in the collective modes of the protein upon ligand binding. This observation, conforming to the Rayleigh-Courant-Fisher theorem, points to a constraining effect imposed by ligands on protein dynamics, which may be accompanied by a negative vibrational entropy difference. Positive frequency shifts in the global modes can thus be linked to the harmonic well getting steeper, because of interactions with the ligand(s).

Phenylpyrrole fungicides act on triosephosphate isomerase to induce methylglyoxal stress and alter hybrid histidine kinase activity

Fludioxonil, a natural product of pyrrolnitrin, is a potent fungicide used on crops worldwide. Drug action requires the presence of a group III hybrid histidine kinase (HHK) and the high osmolarity glycerol (HOG) pathway. We have reported that the drug does not act directly on HHK, but triggers the conversion of the kinase to a phosphatase, which dephosphorylates Ypd1 to constitutively activate HOG signaling. Still, the direct drug target remains unknown and mode of action ill defined. Here, we heterologously expressed a group III HHK, dimorphism-regulating kinase 1 (Drk1) in Saccharomyces cerevisae to delineate fludioxonil’s target and action. We show that the drug interferes with triosephosphate isomerase (TPI) causing release of methylglyoxal (MG). MG activates the group III HHK and thus the HOG pathway. Drug action involved Drk1 cysteine 392, as a C392S substitution increased drug resistance in vivo. Drug sensitivity was reversed by dimedone treatment, indicating Drk1 responds in vivo to an aldehydic stress. Fludioxonil treatment triggered elevated cytosolic methylglyoxal. Likewise, methylglyoxal treatment of Drk1-expressing yeast phenocopied treatment with fludioxonil. Fludioxonil directly inhibited TPI and also caused it to release methylglyoxal in vitro. Thus, TPI is a drug target of the phenylpyrrole class of fungicides, inducing elevated MG which alters HHK activity, likely converting the kinase to a phosphatase that acts on Ypd1 to trigger HOG pathway activation and fungal cell death.

Novel and Selective Rhipicephalus microplus Triosephosphate Isomerase Inhibitors with Acaricidal Activity

The cattle tick Rhipicephalus microplus is one of the most important ectoparasites causing significant economic losses for the cattle industry. The major tool of control is reducing the number of ticks, applying acaricides in cattle. However, overuse has led to selection of resistant populations of R. microplus to most of these products, some even to more than one active principle. Thus, exploration for new molecules with acaricidal activity in R. microplus has become necessary. Triosephosphate isomerase (TIM) is an essential enzyme in R. microplus metabolism and could be an interesting target for the development of new methods for tick control. In this work, we screened 227 compounds, from our in-house chemo-library, against TIM from R. microplus. Four compounds (50, 98, 14, and 161) selectively inhibited this enzyme with IC₅₀ values between 25 and 50 μM. They were also able to diminish cellular viability of BME26 embryonic cells by more than 50% at 50 μM. A molecular docking study showed that the compounds bind in different regions of the protein; compound 14 interacts with the dimer interface. Furthermore, compound 14 affected the survival of partially engorged females, fed artificially, using the capillary technique. This molecule is simple, easy to produce, and important biological data—including toxicological information—are available for it. Our results imply a promising role for compound 14 as a prototype for development of a new acaricidal involving selective TIM inhibition.

Species-Specific Inactivation of Triosephosphate Isomerase from Trypanosoma brucei: Kinetic and Molecular Dynamics Studies

Human African Trypanosomiasis (HAT), a disease that provokes 2184 new cases a year in Sub-Saharan Africa, is caused by Trypanosoma brucei. Current treatments are limited, highly toxic, and parasite strains resistant to them are emerging. Therefore, there is an urgency to find new drugs against HAT. In this context, T. brucei depends on glycolysis as the unique source for ATP supply; therefore, the enzyme triosephosphate isomerase (TIM) is an attractive target for drug design. In the present work, three new benzimidazole derivatives were found as TbTIM inactivators (compounds 1, 2 and 3) with an I₅₀ value of 84, 82 and 73 µM, respectively. Kinetic analyses indicated that the three molecules were selective when tested against human TIM (HsTIM) activity. Additionally, to study their binding mode in TbTIM, we performed a 100 ns molecular dynamics simulation of TbTIM-inactivator complexes. Simulations showed that the binding of compounds disturbs the structure of the protein, affecting the conformations of important domains such as loop 6 and loop 8. In addition, the physicochemical and drug-like parameters showed by the three compounds suggest a good oral absorption. In conclusion, these molecules will serve as a guide to design more potent inactivators that could be used to obtain new drugs against HAT.

Ligand Docking to Intermediate and Close-To-Bound Conformers Generated by an Elastic Network Model Based Algorithm for Highly Flexible Proteins

Incorporating receptor flexibility in small ligand-protein docking still poses a challenge for proteins undergoing large conformational changes. In the absence of bound structures, sampling conformers that are accessible by apo state may facilitate docking and drug design studies. For this aim, we developed an unbiased conformational search algorithm, by integrating global modes from elastic network model, clustering and energy minimization with implicit solvation. Our dataset consists of five diverse proteins with apo to complex RMSDs 4.7-15 Å. Applying this iterative algorithm on apo structures, conformers close to the bound-state (RMSD 1.4-3.8 Å), as well as the intermediate states were generated. Dockings to a sequence of conformers consisting of a closed structure and its "parents" up to the apo were performed to compare binding poses on different states of the receptor. For two periplasmic binding proteins and biotin carboxylase that exhibit hinge-type closure of two dynamics domains, the best pose was obtained for the conformer closest to the bound structure (ligand RMSDs 1.5-2 Å). In contrast, the best pose for adenylate kinase corresponded to an intermediate state with partially closed LID domain and open NMP domain, in line with recent studies (ligand RMSD 2.9 Å). The docking of a helical peptide to calmodulin was the most challenging case due to the complexity of its 15 Å transition, for which a two-stage procedure was necessary. The technique was first applied on the extended calmodulin to generate intermediate conformers; then peptide docking and a second generation stage on the complex were performed, which in turn yielded a final peptide RMSD of 2.9 Å. Our algorithm is effective in producing conformational states based on the apo state. This study underlines the importance of such intermediate states for ligand docking to proteins undergoing large transitions.

Transcriptome analysis of dormant tomonts of the marine fish ectoparasitic ciliate Cryptocaryon irritans under low temperature

BACKGROUND

Cryptocaryon irritans, a species of obligatory ciliate ectoparasite, can infect various species of marine teleost fish. Cryptocaryon irritans that fall to the seabed or aquarium bottom in winter can form "dormant tomonts" and wake up when the temperature rises the next year. Abundant studies and analyses on the dormant tomonts were carried out at the transcriptome level, in order to investigate the molecular mechanism of C. irritans tomonts entering the dormant state under low-temperature conditions.

METHODS

The paired-end sequencing strategy was used to better assemble the entire transcriptome de novo. All clean sequencing reads from each of the three libraries (Group A: untreated blank control; Group B: treated for 24 h at 12 °C; and Group C: developed for 24 h at 25 °C) were respectively mapped back to the transcriptome assembly using the bioinformatics software.

RESULTS

In this study, 25,695,034, 21,944,467, and 28,722,875 paired-end clean reads were obtained respectively from the three cDNA libraries of the C. irritans tomont by Illumina paired-end sequencing technology. A total of 25,925 unique transcript fragments (unigenes) were assembled, with an average length of 839 bp. Differentially expressed genes (DEGs) were scrutinized; in Group B/A pairwise comparison, 343 genes presented differential expression, including 265 up-regulated genes and 78 down-regulated genes in Group B; in Group C/A pairwise comparison, there were 567 DEGs, including 548 up-regulated genes and 19 down-regulated genes in Group C; and in Group B/C pairwise comparison, 185 genes showed differential expression, including 145 up-regulated genes and 40 down-regulated genes in Group B.

CONCLUSIONS

This is the first transcriptomic analytical study of the C. irritans tomonts under low temperature. It can be concluded that most of the genes required for its cell survival under low temperature, or for cell entry into a deeper dormancy state were discovered, and that they might be considered as candidate genes to develop the diagnostic and control measures for cryptocaryoniasis.

Potent and Selective Inhibitors of Trypanosoma cruzi Triosephosphate Isomerase with Concomitant Inhibition of Cruzipain: Inhibition of Parasite Growth through Multitarget Activity

Triosephosphate isomerase (TIM) is an essential Trypanosoma cruzi enzyme and one of the few validated drug targets for Chagas disease. The known inhibitors of this enzyme behave poorly or have low activity in the parasite. In this work, we used symmetrical diarylideneketones derived from structures with trypanosomicidal activity. We obtained an enzymatic inhibitor with an IC50 value of 86 nm without inhibition effects on the mammalian enzyme. These molecules also affected cruzipain, another essential proteolytic enzyme of the parasite. This dual activity is important to avoid resistance problems. The compounds were studied in vitro against the epimastigote form of the parasite, and nonspecific toxicity to mammalian cells was also evaluated. As a proof of concept, three of the best derivatives were also assayed in vivo. Some of these derivatives showed higher in vitro trypanosomicidal activity than the reference drugs and were effective in protecting infected mice. In addition, these molecules could be obtained by a simple and economic green synthetic route, which is an important feature in the research and development of future drugs for neglected diseases.