Crystal structures of Babesia microti lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis

Novel computational and drug design strategies for inhibition of monkeypox virus and Babesia microti: molecular docking, molecular dynamic simulation and drug design approach by natural compounds

Background

The alarming increase in tick-borne pathogens such as human Babesia microti is an existential threat to global public health. It is a protozoan parasitic infection transmitted by numerous species of the genus Babesia. Second, monkeypox has recently emerged as a public health crisis, and the virus has spread around the world in the post-COVID-19 period with a very rapid transmission rate. These two novel pathogens are a new concern for human health globally and have become a significant obstacle to the development of modern medicine and the economy of the whole world. Currently, there are no approved drugs for the treatment of this disease. So, this research gap encourages us to find a potential inhibitor from a natural source.

Methods and materials

In this study, a series of natural plant-based biomolecules were subjected to in-depth computational investigation to find the most potent inhibitors targeting major pathogenic proteins responsible for the diseases caused by these two pathogens.

Results

Among them, most of the selected natural compounds are predicted to bind tightly to the targeted proteins that are crucial for the replication of these novel pathogens. Moreover, all the molecules have outstanding ADMET properties such as high aqueous solubility, a higher human gastrointestinal absorption rate, and a lack of any carcinogenic or hepatotoxic effects; most of them followed Lipinski's rule. Finally, the stability of the compounds was determined by molecular dynamics simulations (MDs) for 100 ns. During MDs, we observed that the mentioned compounds have exceptional stability against selected pathogens.

Conclusion

These advanced computational strategies reported that 11 lead compounds, including dieckol and amentoflavone, exhibited high potency, excellent drug-like properties, and no toxicity. These compounds demonstrated strong binding affinities to the target enzymes, especially dieckol, which displayed superior stability during molecular dynamics simulations. The MM/PBSA method confirmed the favorable binding energies of amentoflavone and dieckol. However, further in vitro and in vivo studies are necessary to validate their efficacy. Our research highlights the role of Dieckol and Amentoflavone as promising candidates for inhibiting both monkeypox and Babesia microti, demonstrating their multifaceted roles in the control of these pathogens.

Lactate Dehydrogenase as a Potential Therapeutic Drug Target to Control Babesia bigemina

Babesia bigemina is a tick-borne apicomplexan hemoprotozoan responsible for bovine babesiosis. The current drugs used for bovine babesiosis treatment have several drawbacks, including toxicity, the lack of effectiveness to clear the parasite, and potential to develop resistance. Identifying compounds that target essential and unique parasite metabolic pathways is a rational approach toward finding alternative drug treatments. Based on the genome sequence and transcriptomics analysis, it can be inferred that anaerobic glycolysis is the dominant adenosine triphosphate (ATP) supply for Babesia, and lactate dehydrogenase (LDH) is one of the essential enzymes in this pathway. Furthermore, the Babesia LDH sequence is distinct from its bovine homologue and thus a potential chemotherapeutic target that would result in decreasing the ATP supply to the parasite but not to the host. Gossypol is a known efficient specific inhibitor of LDH in the sensu stricto B. bovis and the sensu lato B. microti, among other related parasites, but no such data are currently available in the sensu stricto B. bigemina parasites. Hereby, we show that the LDH amino acid sequence is highly conserved among sensu stricto but not in sensu lato Babesia spp. A predictive structural analysis of B. bigemina LDH showed the conservation of the key amino acids involved in the binding to gossypol compared to B. bovis. Gossypol has a significant (P < 0.0001) inhibitory effect on the in vitro growth of B. bigemina, with IC₅₀ of 43.97 mM after 72 h of treatment. The maximum IC (IC₉₈) was observed at 60 mM gossypol. However, a significant effect on the viability of cattle PBMC was observed when the cells were cultured with 60 mM (IC₉₈) gossypol compared with DMSO-exposed control cells. Interestingly, B. bigemina cultured at 3% oxygen expresses significantly higher levels of LDH and is more resistant to gossypol than the parasites maintained at ambient conditions containing ~20% oxygen. Altogether, the results suggest the potential of gossypol as an effective drug against B. bigemina infection, but the risk of host toxicity at therapeutic doses should be further evaluated in in vivo studies.

The Structural Basis of Babesia orientalis Lactate Dehydrogenase

Glycolytic enzymes play a crucial role in the anaerobic glycolysis of apicomplexan parasites for energy generation. Consequently, they are considered as potential targets for new drug development. Previous studies revealed that lactate dehydrogenase (LDH), a glycolytic enzyme, is a potential drug target in different parasites, such as Plasmodium, Toxoplasma, Cryptosporidium, and Piroplasma. Herein, in order to investigate the structural basis of LDH in Babesia spp., we determined the crystal structure of apo Babesia orientalis (Bo) LDH at 2.67-Å resolution in the space group P1. A five-peptide insertion appears in the active pocket loop of BoLDH to create a larger catalytic pocket, like other protozoa (except for Babesia microti LDH) and unlike its mammalian counterparts, and the absence of this extra insertion inactivates BoLDH. Without ligands, the apo BoLDH takes R-state (relaxed) with the active-site loop open. This feature is obviously different from that of allosteric LDHs in T-state (tense) with the active-site loop open. Compared with allosteric LDHs, the extra salt bridges and hydrogen bonds make the subunit interfaces of BoLDH more stable, and that results in the absence of T-state. Interestingly, BoLDH differs significantly from BmLDH, as it exhibits the ability to adapt quickly to the synthetic co-factor APAD⁺. In addition, the enzymatic activity of BoLDH was inhibited non-competitively by polyphenolic gossypol with a K_i value of 4.25 μM, indicating that BoLDH is sensitive to the inhibition of gossypol and possibly to its new derivative compounds. The current work provides the structural basis of BoLDH for the first time and suggests further investigation on the LDH structure of other Babesia spp. That knowledge would indeed facilitate the screening and designing of new LDH inhibitors to control the intracellular proliferation of Babesia spp.

Kinetic Characterization and Inhibitor Screening of Pyruvate Kinase I From Babesia microti

The apicomplexan Babesia microti is a main pathogenic parasite causing human babesiosis, which is one of the most widely distributed tick-borne diseases in humans. Pyruvate kinase (PYK) plays a central metabolic regulatory role in most living organisms and catalyzes the essentially irreversible step in glycolysis that converts phosphoenolpyruvate (PEP) to pyruvate. Hence, PYK is recognized as an attractive therapeutic target in cancer and human pathogens such as apicomplexans. In this study, we cloned, expressed, and purified B. microti PYK I (BmPYKI). Western blotting illustrated that anti-rBmPYKI antibody could specifically recognize the native BmPYKI protein in the lysate of B. microti with a 54-kDa band, which is consistent with the predicted size. In addition, the enzymatic activity of the purified recombinant PYKI (rPYKI) was tested under a range of pH values. The results showed that the maximum catalytic activity could be achieved at pH 7.0. The saturation curves for substrates demonstrated that the K_m value for PEP was 0.655 ± 0.117 mM and that for ADP was 0.388 ± 0.087 mM. We further investigated the effect of 13 compounds on rBmPYKI. Kinetic analysis indicated that six inhibitors (tannic acid, shikonin, apigenin, PKM2 inhibitor, rosiglitazone, and pioglitazone) could significantly inhibit the catalytic activity of PYKI, among which tannic acid is the most efficient inhibitor with an IC₅₀ value 0.49 μM. Besides, four inhibitors (tannic acid, apigenin, shikonin, and PKM2 inhibitor) could significantly decrease the growth of in vitro-cultured B. microti with IC₅₀ values of 0.77, 2.10, 1.73, and 1.15 μM. Overall, the present study provides a theoretical basis for the design and development of new anti-Babesia drugs.

Identifying the Naphthalene-Based Compound 3,5-Dihydroxy 2-Napthoic Acid as a Novel Lead Compound for Designing Lactate Dehydrogenase-Specific Antibabesial Drug

Human babesiosis is caused by apicomplexan Babesia parasites, including Babesia microti, Babesia crassa, Babesia sp. MOI, Babesia divergens, Babesia duncani, and Babesia venatorum. Among them, B. microti is the most common cause of human and rodent babesiosis. Currently, no vaccine is available, and drugs for the treatment have high failure rates and side effects. Due to lack of a traditional tricarboxylic acid cycle (TCA cycle) and its dominant dependence on anaerobic metabolism to produce ATP, B. microti lactate dehydrogenase (BmLDH) was assumed to play a critical role in B. microti ATP supply. Our previous study demonstrated that BmLDH is a potential drug target and Arg99 is a crucial site. Herein, a molecular docking was performed based on the crystal structure of BmLDH from a series of gossypol derivatives or structural analogs to find the potent inhibitors interacting with the residue Arg99, and three naphthalene-based compounds 2,6-naphthalenedicarboxylic acid (NDCA), 1,6-dibromo-2-hydroxynapthalene 3-carboxylic acid (DBHCA), and 3,5-dihydroxy 2-napthoic acid (DHNA) were selected for further tests. Enzyme activity inhibitory experiments show that DBHCA and DHNA inhibit recombinant BmLDH (rBmLDH) catalysis with ~109-fold and ~5,000-fold selectivity over human LDH, respectively. Surface plasmon resonance (SPR) assays demonstrate that DHNA has a lower K _D value to BmLDH (3.766 x 10^-5 M), in contrast to a higher value for DBHCA (3.988 x 10^-8 M). A comparison of the kinetic parameters [association constant (k _a) and dissociation constant (k _d) values] reveals that DBHCA can bind the target faster than DHNA, while the complex of DHNA with the target dissociates slower than that of DBHCA. Both DBHCA and DHNA can inhibit the growth of B. microti in vitro with half-maximal inhibitory concentration (IC₅₀) values of 84.83 and 85.65 μM, respectively. Cytotoxicity tests in vitro further indicate that DBHCA and DHNA have selectivity indexes (SI) of 2.6 and 22.1 between B. microti and Vero cells, respectively. Although the two naphthalene-based compounds only display modest inhibitory activity against both rBmLDH and the growth of B. microti, the compound DHNA features high selectivity and could serve as a novel lead compound for designing LDH-specific antibabesial drug.