Exploring Novel Inhibitors for Babesia bigemina Lactate Dehydrogenase: A Computational Structural Biology Perspective

Exploring Novel Inhibitors for Babesia bigemina Lactate Dehydrogenase: A Computational Structural Biology Perspective | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Exploring Novel Inhibitors for Babesia bigemina Lactate Dehydrogenase: A Computational Structural Biology Perspective Safiye Merve Bostancioglu, Ozal Mutlu This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4951492/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 07 Jan, 2025 Read the published version in Parasitology Research → Version 1 posted 4 You are reading this latest preprint version Abstract Babesia bigemina is an apicomplexan parasite and causes “Texas fever” in bovines. Lactate dehydrogenase (LDH) is an essential enzyme in regulating anaerobic metabolism. The presence of five amino acid insertion in the active site of many apicomplexan LDH makes a significant difference between the parasitic LDH and the host counterpart. Therefore, apicomplexan LDH is an attractive drug target. In this study, a structure-based drug discovery approach was performed to find novel inhibitor candidates. In the first round, possible candidates were identified by following the virtual screening workflow. Then, the compounds with favourable docking scores were filtered using the QM-polarized ligand docking and induced fit docking methods. As a result, 20 novel compounds that bind to the active site of the BbigLDH and have low affinity with the host LDHs have been identified. Molecular dynamics simulations of the complexes (in total 8.8 µs) were performed and binding free energies were calculated. In conclusion, compounds named C09, C16 and C18 deserve further investigation to better understand their potential therapeutic effects on babesiosis. The findings of this study, elucidating the structural properties of BbigLDH enzyme and discovering its potential inhibitors, might pave the way for further research in developing LDH-targeted therapeutic interventions. Lactate dehydrogenase Babesia bigemina drug discovery Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Introduction Smith and Kilborne (1893) described the 'Texas fever' for the first time, known as babesiosis, which is caused by the apicomplexan parasite Babesia bigemina and transmitted through tick vectors [ 1 ]. The disease manifests in a spectrum of states ranging from asymptomatic carriers to severe cases marked by hemolytic anemia, fever, hemoglobinuria, and, in some instances, fatalities. Its impact extends across a diverse range of mammals, including domestic pets, farm animals, wildlife, and even humans [ 2 ], [ 3 ], [ 4 ]. B. bigemina currently prevail in tropical and subtropical regions worldwide, spanning continents such as Africa, Asia, the Americas, Southern Europe, and Australia [ 3 ]. Their presence directly influences meat and milk production, thereby influencing the competitive edge of livestock industries [ 5 ]. The primary vectors for these pathogens are Rhipicephalus (Boophilus) microplus , R. annulatus, R. australis, R. decoloratus, R. bursa, R. geigyi , and R. evertsi [ 2 ], [ 3 ]. Currently, diminazene aceturate (DA) and imidocarb dipropionate stand as the most commonly employed treatments [ 6 ], [ 7 ], [ 8 ]. However, the continued or inappropriate use of these drugs poses a significant risk of engendering drug resistance [ 6 ], [ 9 ], [ 10 ], and they also carry the potential for chemical residues in meat and milk, which can have adverse effects on human health [ 11 ], [ 12 ], [ 13 ]. This pressing concern necessitates an urgent quest for the discovery or development of alternative, effective, and economically viable drugs with minimal toxicity [ 14 ], [ 15 ]. Within apicomplexans, lactate dehydrogenase (LDH) serves as a key glycolytic enzyme, facilitating the transformation of pyruvate into lactate. [ 16 ]. It bidirectionally catalyzes the conversion of pyruvate to lactate, utilizing NADH as a co-factor and regenerating NAD + essential for glycolysis to proceed, or the conversion of lactate to pyruvate, utilizing NAD + as a co-factor and regenerating NADH [ 17 ]. The energy generated through this process is harnessed by parasites to fuel their biochemical processes and ensure their survival. The substrate specific loop, which is present in the active site, is designated for its substrate in parasitic LDH. Hosts represent canonical LDH which absence of five amino acids in the substrate specific loop differs from parasitic LDH. The five amino insertion was firstly described in LDH of Plasmodium falciparum [ 18 ]. Further studies showed that the insertion site is conserved in ancient LDH and LDH of many apicomplexan parasites including B. bigemina , but in a different amino acid sequence [ 5 ], [ 17 ], [ 19 ], [ 20 ], [ 21 ], [ 22 ], [ 23 ], [ 24 ]. The vital importance of the enzyme, and the fundamental difference between the host's LDH and the parasite's LDH in the active site makes it a more specific drug target. Further investigation of P. falciparum LDH demonstrates the potential for identifying a specific inhibitor for the enzyme [ 25 ]. Other studies related to LDH in apicomplexan parasites were published for T. gondii, C. parvum , and supports the hypothesis that apicomplexan LDH could be used as a drug target [ 26 ], [ 27 ]. The crystal structure of Babesia microti LDH (BmLDH) is the first LDH whose crystal structure has been solved in babesia genus. However, it is important to note that the BmLDH shows significant differences from apicomplexans LDHs and its structure is closely related to mammalian-like LDH [ 28 ], [ 29 ]. Yu et al., 2022 was identified a novel cDNA clone that encodes apicomplexan type LDH from Babesia orientalis (BoLDH) and solved its crystal structure in its apo form [ 24 ]. On the other hand, gossypol, an LDH inhibitor, has been explored for its impact on various Babesia species. A concentration with an IC50 value of 50 µM gossypol irreversibly halted the in vitro growth of B. bovis. Crucially, the addition of gossypol led to a noteworthy inhibition of enzymatic activity, acting as a competitive inhibitor for the binding of NADH [ 20 ]. Another study for gossypol, exhibits inhibitory impact on the in vitro growth of B. bigemina , with an IC50 of 43.97 µM after 72 hours of treatment. The maximum inhibitory concentration (IC98) is observed at 60 µM gossypol [ 5 ]. These findings demonstrate that LDH is a significant drug target for Babesias as an apicomplexan parasite. Modern progress in computer-assisted drug development has made it possible to design drugs specifically for target proteins, thereby accelerating the creation of new medications. Structure-Based Virtual Screening has notably transformed and improved the process of drug discovery, establishing itself as one of the most effective computational techniques for drug design [ 30 ]. The principles and applications of Virtual Screening encompass the entire process, from the initial stages of receptor and library preparation, through docking and scoring, to the post-processing of the highest-scoring hits. This study focused on understanding the structure of Babesia bigemina LDH (BbigLDH) enzyme and discovering candidate inhibitors by structure-based drug discovery methods. After identifying the hit compounds using virtual screening, molecular dynamics simulations were performed to assess protein-ligand complex stability and calculate end-point binding free energies. This research could contribute to the identification of potential LDH inhibitors for the treatment of Babesia bigemina infection and to improve drug design efforts against bovine babesiosis and diseases caused by other apicomplexans. Materials & Methods Sequence Analysis Protein sequences of Bos taurus LDH B (CAB96751.1), Bos taurus LDH A (BAA14170.1), Babesia bovis T2Bo (XP_001611047.1), Babesia bigemina (XP_012766822.1), Babesia ovata (XP_028866450.1), Theileria equi strain WA (XP_004833698.1), Theileria orientalis strain Shintoku (XP_009689964.1), Theileria parva strain Muguga (XP_766703.1), Theileria annulata (ADG45564.1), Plasmodium malariae (SBS88125.1), Plasmodium yoelii (XP_724101.1), Plasmodium falciparum 3D7 (XP_001349989.1), Toxoplasma gondii ME49 (XP_002368488.1) and Eimeria acervulina (ACM77785.1) were obtained from the NCBI ( https://www.ncbi.nlm.nih.gov/ ). ClustalW (1.2.4) was used for multiple sequence alignment and the aligned sequences were visualized via ESPript 3.0 webserver [ 24 ]. The phylogenetic tree which is assessed by using the Maximum Likelihood method was constructed via MEGA 11 [ 31 ]. Homology Modelling and Model Validations BbigLDH, whose experimentally verified 3-D structure is unknown, were generated by homology modelling using MODELLER 9.23 software [ 32 ]. Template for homology modelling was determined using NCBI BLAST tool [Basic Local Alignment Seach Tool ( http://www.ncbi.nlm.nih.gov/blast )]. Resulting PDB structures with the highest max score, the lowest E-value, and the highest % identity score were evaluated. T. gondii LDH1 apo form (1PZE) [ 17 ] and ancestral apicomplexan LDH with malate (4PLC) [ 19 ] were chosen as templates for the modelling of the open and closed states of the BbigLDH enzyme, respectively. A total of 100 models were generated, and the model with the lowest DOPE (Discrete Optimized Protein Energy) and z-DOPE (normalized DOPE) scores was selected for further in silico analyses. UCSF Chimera 1.15 was employed to energy minimize the BbigLDH models, adding any missing hydrogen atoms and charges [ 33 ]. The minimization process involved 2000 steps of "Steepest Descent" followed by 1000 steps of "Conjugate Gradient," utilizing the AmberFF14SB force field. Quality of the protein models was evaluated by using servers including ERRAT [ 34 ], ProSA [ 35 ], [ 36 ], ProCHECK / Ramanchandran Plot [ 37 ], Verify3D [ 38 ] and QMEAN [ 39 ]. To assess the reliability of specific atomic interactions, ERRAT was utilized. ProSA was employed to verify the correctness of the overall 3D homology model structure of BbigLDH and ensure proper energy optimization. Verify3D was used for 3D structure verification by computing position-specific 3D-1D scores for each residue, referencing a database of known accurate folding structures. The QMEAN tool provided a quantitative evaluation of significant geometric features in the protein structure. The Ramachandran Plot (ProCHECK) was used to visualize the φ (phi) and ψ (psi) dihedral angles of the amino acid residues, identifying conformational outliers and confirming that these angles are within the favorable regions, which are indicative of proper folding and stability. The LDH structures of the host ( Bos taurus ) was obtained from the Alphafold database (Uniprod ID: (LDHA - P19858, LDHB - Q5E9B1, and LDHAL6B - Q3T056) [ 40 ] ( https://alphafold.ebi.ac.uk/ ). Cognate docking and Enrichment Calculations Cognate docking was performed using the Glide XP (extra precision) scoring function and docking methodology to asses prediction quality of the docking program [ 41 ], [ 42 ]. The 3-D crystal structures of P. falciparum LDH (PfLDH) with inhibitors (PDB IDs: 1LDG [ 18 ], 1CET [ 43 ], 1T26, 1T24, 1T25 [ 44 ], 1U5A, 1U4O, 1XIV, 1U4S [ 45 ]) were obtained from the Protein Data Bank [ 46 ]. The co-crystallized ligands were split from the protein structure, then protein and ligand structures were prepared using the protein preparation workflow respectively and re-docked [ 47 ], [ 48 ]. Additionally, these ligands were analyzed for their interactions with BbigLDH using XP docking. Heavy-atom RMSD values were calculated between the docked poses and the crystal structure poses using the superposition tool from Maestro [ 49 ]. In addition, docking methodology was validated by decoys set of 1000 molecules [ 50 ]. Known 34 LDH inhibitors were obtained from the BRENDA Enzyme database [ 51 ]. 1000 decoy compounds and 34 active molecules were prepared with LigPrep [ 48 ] and docked by Glide XP mode with default settings. ROC, RIE, AUAC and BEDROC parameters were analysed by Enrichment Factor Calculation. Structure-Based Virtual Screening, QM-Polarized Ligand Docking (QPLD) and Induced Fit Docking (IFD) Compound library including about 570.000 compounds from the ChemDiv (Targeted diversity library, https://www.chemdiv.com/catalog/diversity-libraries/targeted-diversity-library/ ), Enamine (Hit Locator Library HLL-460, Discovery Diversity Set DDS-50, Discovery Diversity Set DDS-10, https://enamine.net/compound-libraries/diversity-libraries ) and OTOVAchemicals (PrimScreen1, PrimScreen2, PrimScreen3, PrimScreen5, PrimScreen10 https://otavachemicals.com/products/compound-libraries-for-hts/diversity-sets ) were prepared using the LigPrep [ 48 ]program. Ionization states of the compounds were predicted using the Epik at 7.5 pH [ 52 ]. Protein Preparation Workflow was performed to prepare protein structures for docking[ 47 ], [ 53 ]. Receptor grid generation was used to define the atomic positions of the binding site for each structure. Glide-Virtual Screening Workflow was performed in three steps: High throughput virtual screening (HTVS-10%), Standard precision (SP-10%), and Extra precision (XP- the best docking scores of 100 ligands). Compounds were ranked according to XP-GScore and the best compounds having interactions with the catalytic residues and lower XP-GScores (≤-7.0 kcal/mol) were selected. Afterwards, ligands are filtered based on their ADME properties [ 54 ]. The candidates identified are within the reference values specified in the Qikprop manual. For compounds to qualify as drug candidates, several criteria and appropriate value ranges must be considered. Among these, the key five principles of Lipinski are essential [ 55 ]. Compounds with acceptable scores were then used in Quantum Mechanical Polarized Ligand Docking (QPLD). The QPLD is a computational method that combines Glide docking with QSite calculations to determine the ligand binding positions [ 56 ]. For the initial docking, the Glide XP method is employed with accurate QM level chosen for QM charges using the Jaguar module in gas phase [ 57 ], [ 58 ] Redocking is performed using the XP method. In result, ligand conformations with the lowest Emodel values obtained from the docking was selected. Compounds with XP-GScore value equal to or lower than − 7.0 kcal/mol were filtered for the induced fit docking. IFD is a molecular docking method that involves ligand repositioning onto a relaxed receptor using Prime. In this approach, a receptor grid was generated by selecting active site residues for molecular docking. Prime was utilized to perform receptor rearrangement for residues within a 5 Å distance range. The relaxed receptor was then subjected to redocking using the XP method of Glide program. The resulting poses were evaluated based on the lowest IFD scores. Compounds identified post-IFD were subjected to Glide XP docking against the structures of Bos taurus LDH A, LDH B and LDH AL6B. Compounds displaying low affinity with the host LDH structures were selected for further analysis, considering them as potential candidates for selective inhibition. Molecular Dynamics Simulations and Binding Free Energy Calculations with MM/GBSA MD simulation was performed for open and closed states of homology models by using the Nanoscale Molecular Dynamics (NAMD) program [ 59 ]. During the preparation of closed states, ligands were processed using the CHARMM-GUI server [ 60 ] with the inclusion of hydrogen atoms facilitated by Open Babel, followed by conversion to the mol2 format [ 61 ]. The coenzyme (NAI) was then added into the system and parametrized. The protein and ligand PDB and PSF files were generated and combined using VMD 1.9.3 software to create the complex structure [ 62 ]. The complex structure was solvated in a cubic box and subsequently neutralized, and the ionic strength was adjusted to 150 mM NaCl. Before MD simulations, 2-tier of minimization (total 5000 steps) was conducted. In the first step, constraint with constraint scaling of 100 kcal/mol/A^2 was applied to the solids, while in the second step, the entire system was minimized without any constraints. Equilibration was carried out using the NVT ensemble for 250 ps at 300 K by applying constraints (25 kcal/mol/A^2) to the solids. This was followed by an addition equilibration and relaxation step with the NPT ensemble for 275 ps at 300 K and 1 atm using the Nose–Hoover Langevin piston. The constraints on the solid were gradually reduced from 10 kcal/mol/A^2 to zero at every 25 ps. Finally, a production phase consisting of a 100 ns MD simulation was performed without any constrains using the NPT ensemble (300 K and 1 atm). All MD simulations were conducted using the CHARMM36 forcefield [ 63 ]. The MD simulation trajectories were analysed for the RMSD and RMSF values and visualized using the VMD 1.9.3 software [ 62 ]. Molecular dynamics simulations of the protein-ligand complexes were performed using the Desmond software (OPLS4 forcefield) [ 64 ] or GROMACS 2022 (CHARMM36 forcefield) [ 63 ] to asses ligand stability and to calculate binding energies of compounds. The system consisted of TIP3P water molecules in an orthorhombic box. To ensure charge neutrality and achieve an ionic strength of 150 mM NaCl, appropriate ions were added. Prior to production runs, a minimization step for 100 picoseconds and a default relaxation protocol provided by Desmond was employed to minimize and relax the system. The production runs were conducted for 20, 300, or 1000 ns at the different stages of compound selection procedure for the next stages. All production runs were performed at a temperature of 300 K and a pressure of 1 atmosphere. The integration time step was set as 2 femtoseconds. Temperature and pressure were controlled using the Nose-Hoover thermostat and Martyna-Tobias-Klein (MTK) barostat, respectively. A cutoff distance of 9 Å was applied for nonbonded interactions to calculate the short-range forces efficiently. Protein-ligand interactions were investigated with simulation interaction diagram. MM/GBSA method was employed to calculate the binding free energies (ΔG bind) using the thermal_mmgbsa.py script and the Prime program [ 65 ]. In parallel, 100 ns of MD simulations were performed using the GROMACS by applying the same parameters as above. Minimization was stopped when the maximum force < 10.0 kJ/mol. NVT equilibration (100 ps) was employed by the Particle Mesh Ewald method applied to calculate long range electrostatic interaction and the modified Berendsen thermostat (V-rescale) used to control the temperature of the simulation system. NPT equilibration (100 ps) and production run was performed using the C-rescale pressure coupling at 1 atm and V-rescale used to control the temperature of the simulation system. For the analysis of MD trajectories, GROMACS modules (gmx rms and gmx rmsf) were used, focusing on parameters such as root mean squared deviation (RMSD) and fluctuation (RMSF). Results & Discussions Assessment of LDH sequence similarity among apicomplexan species Upon analysis of its nucleotide sequence, it was determined that the BbigLDH gene contains an open reading frame spanning 969 nucleotides, which encodes for a sequence of 332 amino acids. The alignment of sequences was conducted using ClustalW and visualized through ESPript 3 (Supplement file, Fig. 1 ). The highest sequence identity of 95.34% and 85.98% was observed between BbigLDH and Babesia ovata and Babesia bovis, respectively. Other species closely related to BbigLDH, such as Theileria orientalis exhibited a 70.66% match, while Plasmodium falciparum showed a 50.31% of identity. Notably, there is a 30% identity with Bos taurus LDH A, LDH B, which is recognized as the host for Babesia bigemina . The phylogenetic tree, constructed using the Neighbor Joining method, is provided in Supplement file, Fig. 3 [ 66 ]. The phylogenetic tree suggests that BbigLDH is positioned within the Babesia clade, with a closer evolutionary relationship to B. ovata . In contrast, LDHs originating from mammalian hosts form a distinct clade, indicating a clear separation in their evolutionary lineage from those of the apicomplexan parasites. The InterPro database analysis of BbigLDH sequence (XP_012766822.1) revealed the presence of LDH1N-NAD binding (PF0056 / 9-152) and LDH1C-alpha/beta C terminal (PF02866 / 157–315) domains. The 3D model of BbigLDH is expected to include a binding site for NADH or NAD + cofactors as well as a binding site for pyruvate or lactate substrates. The secondary structure of BbigLDH consists of sixteen alpha-helices, fourteen beta-strands, and various loops (Supplement file, Fig. 2 ). BbigLDH displays a common five-amino acid insertion within its substrate specificity loop. Crystallographic analysis PfLDH have shown that this insertion contributes to a distinctive architecture in the active site, which could be a target for drug development [ 18 ]. While the sequences in this region vary across species, for BbigLDH, the sequence is DDEWT, and it is predicted to form a coil structure (Fig. 2 A). BbigLDH shares identical catalytic residues with those conserved in all LDHs, which include R109, D168, H195, and R171 (residues R98, D158, H185, R161 in BbigLDH) (Supplement file, Fig. 1 ) [ 5 ]. Residues 93–97 (DDEWT) typically form a mobile loop in most apicomplexan LDHs (Fig. 2 A). In the ternary complex, the loop moves downward to form the closed catalytic site [ 19 ]. The loop opposite to the substrate specificity loop (G236, T237, G238 in BbigLDH) is called as “opposing loop” and depicted in Fig. 2 B and Supplement file, Fig. 5 [ 19 ]. The distinctive two-residue deletion within the opposing loop of modern Plasmodium LDHs, along with the shared loop structure with ancient LDH (AncLDH) and Toxoplasma LDHs, marks a significant evolutionary deviation (Supplement file, Fig. 5 ) [ 19 ]. Additionally, while ancestral LDH displays a slight oxaloacetate activity in common with modern Toxoplasma LDHs, this capability is absent in Plasmodium LDHs. This intriguing correlation strongly suggests that the opposing loop deletion, and potentially the presence of Ala236 and Pro246, play a pivotal role in shaping the remarkably stringent substrate specificity observed in modern Plasmodium LDHs [ 67 ]. Unlike in modern Plasmodium LDHs, a deletion in the opposing loop is not observed in BbigLDH, like Toxoplasma gondii LDH. This sets it apart from the modern Plasmodium LDHs in this aspect. In Babesia and Theileria (excluding T. equi ) genus, a G[QSTN]G pattern is present (Fig. 2 ). Unlike in Plasmodium and B. taurus , such a pattern is not found, and there is a distinction in the opposing loop from both ancLDH and Toxoplasma LDH. While residues preceding the opposing loop in B. taurus and Plasmodium LDHs are positively charged (H and K residues), in Babesia and Theileria , these residues are observed to be aromatic (Y). Homology modelling and model validation The crystal structures of monomers in both closed and open states of BbigLDH were previously unknown. To address this, homology models were generated using the MODELLER 9.23 program. These models were developed for both open and closed states, incorporating pyruvate and NADH, respectively. The suitable template proteins were selected from PDB-structured proteins using BLASTp, based on the highest identity percentage. Additionally, during the selection process, emphasis was placed on ensuring that the template LDH structures belonged to apicomplexans and contained NADH and pyruvate in closed state. The template LDH proteins used for modelling included 1PZE (1.95 Å) [ 68 ] - the T. gondii LDH1 apoenzyme structure- for the open state, and 4PLC (1.50 Å) - an ancestral LDH from apicomplexans [ 67 ]- for the closed state, along with pyruvate and NADH. The generated 100 models were selected based on the lowest DOPE value and subsequently validated using various methods. One of the validation programs employed was ERRAT, which examines the consistency of atomic interactions within the modelled structure. The overall quality factors obtained from the ERRAT server for the open and closed models were 90.1274 and 93.949, respectively [ 69 ]. These values indicate a relatively high quality of the modelled structures. Another validation program utilized was ProSA [ 70 ], which evaluates the quality of generated models by constructing an energy profile of structural residue areas. The Z-scores obtained from ProSA were compared with those of experimentally validated structures to assess consistency. According to ProSA results, the Z-scores for mutant and wild-type models were − 3.3 and − 3.19, respectively. Models falling within the acceptable range of residue and energy efficiency were considered reliable. Quantitative Model Energy Analysis (QMEAN) is a valuable tool for analyzing crucial geometric aspects of protein structures, such as determining natural curvatures and calculating torsion angle potentials of amino acids [ 71 ]. The QMEAN value is originally in a range [0,1] with one being good expected from high resolution X-ray structures. In the results, QMEAN4 values were − 2.47 and − 2.19 for open state and closed state, respectively. Verify3D is a validation program that evaluates the absolute secondary and tertiary folding behaviour of a model protein by analyzing the 3D-1D scores in the context of 3D structure-sequence relationships. A higher value for amino acids is expected compared to the threshold value, and the overall simplified score determines the structure's validity [ 72 ]. Homology models passed by 88.82% and 93.79% for open state and close state, respectively according to the assessment based on favourable 3D-1D scores. A Ramachandran plot graph was generated for the model proteins using ProCHECK [ 73 ]. The percentages of residues in the favoured region were 88.1% and 86.6% respectively for the open and closed states. Based on the assessment of validation values, all these validation results confirmed that the enzymes were modelled with high accuracy and reliability (Fig. 3 A). Additionally, the structural difference introduced by the 5-AA insertion in the open form of BbigLDH compared to host LDH A and B is shown in Fig. 3 B. Detailed results obtained from validation were provided in Supplement file, Table 1 . Table 1 Twenty compounds, exhibiting lower affinity with BtLDH A, B and A-like B and their average values (kcal/mol) Compound no Compound ID BtLDHA BtLDHB BtLDHA-likeB Average C1 Z1478445725 -4.311 -5.156 -4.875 -4.78067 C2 Z1246322905 -4.495 -5.448 -3.34 -4.42767 C3 Z826473056 -4.253 -4.534 -4.414 -4.40033 C4 Z509196536 -4.233 -5.03 -4.275 -4.51267 C5 Z1424270493 -4.559 -4.411 -5.247 -4.739 C6 Z278525320 -4.464 -4.047 -4.283 -4.26467 C7 Z1640818270 -4.863 -5.026 -4.731 -4.87333 C8 Z604033838 -3.745 -3.219 -4.344 -3.76933 C9 Z2373310039 -4.842 -5.14 -4.485 -4.82233 C10 Z973658262 -4.959 -3.858 -3.883 -4.23333 C11 Z1616793138 -4.507 -4.505 -4.672 -4.56133 C12 Z2953881032 -4.596 -3.577 -4.391 -4.188 C13 Z1253206711 -4.843 -4.798 -4.772 -4.80433 C14 Z19650139 -4.249 -4.068 -4.262 -4.193 C15 Z2400537901 -5.409 -4.512 -4.829 -4.91667 C16 Z2201943136 -4.8 -4.406 -5.435 -4.88033 C17 Z2060912694 -4.114 -4.636 -4.857 -4.53567 C18 Z2063989764 -3.427 -5.208 -5.287 -4.64067 C19 Z1531526443 -4.526 -5.325 -4.901 -4.91733 C20 Z2055833084 -4.401 -4.787 -4.808 -4.66533 Validation of Molecular Docking Procedure The validation of the molecular docking procedure to be used in computational biology methods is an essential preliminary step to ensure the validity and reliability of the method. In the validation process, inhibitors present in the crystal structures of PfLDH with known structures were redocked into the original protein structure using the cognate docking procedure in methods section. The superposition of the resulting protein structure with the original protein structure is examined, along with the binding mode of the ligand and the RMSD values. Low RMSD values (The superimposed values for 1U4O, 1T25 and 1XIV are 0.14 Å, 0.21 Å, and 0.32 Å, respectively.) indicate a high degree of similarity between the two structures. Other RMSD values for 1U4S and 1T26 are 1.26 Å and 1.06 Å, respectively which have more higher superimposed values then other PDB structures. Additionally, when evaluating the cognate docking results, a threshold value of -7.00 kcal/mol for the docking score has been determined as the preferred threshold for assessing the HTVS docking results. The enrichment calculation results provide valuable information about the performance of a computational model or algorithm in virtual screening or drug discovery experiments. The model was tested on a dataset consisting of 1000 ligands, including 34 known actives. The goal was to determine the model's ability to accurately rank the active compounds and enrich them at the top of the ranked list. The obtained BEDROC values ranged from 0.229 to 0.332, suggesting a moderate enrichment capability. The receiver operating characteristic (ROC) curve and its associated area under the curve (AUC) were also analyzed to evaluate the overall classification accuracy of the model. The ROC value of 0.69 suggests a reasonable discriminative ability, indicating that the model can distinguish between active and inactive compounds to a certain extent. The area under the accumulation curve (AUC) is 0.68, indicating that a significant portion of the active compounds is discovered early in the screening process. The obtained RIE value of 3.43 indicates that the model is performing better than random selection but still leaves room for improvement. It highlights the model's ability to detect structurally diverse active compounds (Supplement file, Fig. 4 and Table 2 ). Structure-Based Virtual Screening, QM-Polarized Ligand Docking and Induced Fit Docking (IFD) Four structures of BbigLDH were utilized, comprising two open and two closed states. A pair of open and closed states were generated after a 100 ns MD simulation using NAMD. For the prediction of binding sites, DoGSiteScorer method was employed on all models. [ 74 ]. In both the open and closed states of BbigLDH, druggable pockets were predicted, with a druggability score of 0.81 and simple score of 0.6. The pocket includes the pyruvate and NADH binding regions. The virtual library was initially docked to the active site of BbigLDH using GLIDE HTVS, SP and XP docking methods, respectively. The top 100 poses were then ranked based on XP GScores (Table 2 ). After comparison of compounds interacting with common residues, 256 ligands were selected. Subsequently, the remaining ligands underwent ADME analysis using the QikProp program. Compounds within the specified reference value range from the manual were identified (133 compounds). After ADME analysis, the Quantum Mechanics Polarized Ligand Docking method was applied, where the QM charges were calculated by using the QSite (Jaguar) and each ligand was redocked from its corresponding pose set. The poses with the lowest Glide Emodel values were chosen. Results with XP GScore values equal to or higher than − 7.00 kcal/mol were ranked to obtain the best outcomes for 70 models. Following this, induced fit docking was employed This involved relaxing the ligand into a conformationally relaxed receptor state, followed by redocking using Glide. The resulting 70 poses were ranked based on their IFD score, with lower scores being prioritized. Remaining compounds were also docked to the host LDH (Bos taurus LDH A, LDH B and LDH AL6B) to identify selective inhibitor candidates. The structural difference introduced by the 5-AA insertion in the open form of BbigLDH compared to host LDH A and B is shown in Fig. 3 B. Twenty compounds, exhibiting lower affinity with host LDHs and available for purchase were selected for molecular dynamics simulations and end-point binding free energy calculations (Table 1 ). The evaluation of potential drug candidates involves a thorough assessment of various criteria, with a close examination of their physicochemical properties intricately linked to their structural attributes. Minor modifications in their structures have resulted in enhanced properties. In our study, the hit compounds underwent analysis to determine their compliance with Lipinski's rule of five. [ 55 ]. According to Lipinski's rule of five, an established set of guidelines for drug likeness, the molecular weight of potential drug candidates should be less than 500 Da. Additionally, the hydrogen bond donor count should be fewer than 5, the hydrogen bond acceptor count should be less than 10, and the predicted octanol/water partition coefficient should ideally fall within the range of -2.0 to 4. These criteria help assess the likelihood of a compound's success as a drug candidate based on its physicochemical properties [ 75 ]. All values for the compounds fall within the specified range for Lipinski rules and PlogP 0/w . The Plog values for both hERG channel blocking and serum albumin binding fell within the reference ranges [ 76 ]. The predicted IC50 value for blocking HERG K + channels should fall within the accepted range, with values greater than − 5.0. In this study, the compound with the largest value is C18 (5.073), while the one with the smallest value is compound C9 (-2.504). On the other hand, the predicted apparent Caco-2 cell permeability should be assessed, with values indicating low permeability if less than 25 and high permeability if greater than 500 [ 77 ]. According to this information, it is observed that the compound with the highest permeability is C20 (2104.635), while the one with the least permeable is C19 (74.968). The predicted blood-brain barrier permeability coefficient should be within the range of -3.0 to 1.2 [ 78 ]. The predicted binding capacity to human serum albumin is considered appropriate if it falls between − 1.5 and 1.5 [ 79 ]. All values for the compounds fall within the specified range for PlogBB, PlogKhsa. The predicted human oral absorption, rated on a scale of 0 to 100, is considered favorable if it exceeds 70%. The compound with the highest human oral absorption rates is C18 and C20 (both 100), whereas C19 has a rate of 72.307 (Supplement file, Table 3). Table 2 List of compounds which have favourable XP GScores. Free binding energy (∆G) values were listed for the last 5 ns of 20 ns MD simulation. Compound no Compound ID Structure of BbigLDH HTVS; XP gscore QPLD; XP gscore IFD; XP gscore Prime energy IFD score dG Average C1 Z1478445725 Closed state after MDS -8.343 -8.515 -9.569 -13590.36 -689.09 -67.3425 C2 Z1246322905 Closed state after MDS -8.608 -8.140 -9.048 -13502.23 -684.16 -58.4571 C3 Z826473056 Closed state after MDS -8.236 -8.530 -11.017 -13460.85 -684.06 -64.0088 C4 Z509196536 Closed state -8.997 -8.389 -10.565 -13452.59 -683.19 -57.2994 C5 Z1424270493 Closed state -8.869 -7.598 -9.404 -13471.78 -682.99 -60.7591 C6 Z278525320 Closed state after MDS -9.497 -7.480 -10.557 -13439.77 -682.55 -54.8145 C7 Z1640818270 Closed state after MDS -8.479 -8.473 -7.652 -13462.95 -680.8 -64.4233 C8 Z604033838 Closed state -8.382 -8.250 -10.285 -13409.02 -680.74 -66.175 C9 Z2373310039 Closed state after MDS -8.173 -7.958 -7.998 -13442.56 -680.13 -75.825 C10 Z973658262 Closed state -9.176 -9.182 -9.803 -13401.97 -679.9 -46.3542 C11 Z1616793138 Closed state -8.648 -8.555 -9.624 -13403.46 -679.8 -46.0935 C12 Z2953881032 Closed state after MDS -8.202 -7.972 -8.058 -13428.81 -679.5 -71.5432 C13 Z1253206711 Closed state -8.925 -7.731 -10.617 -13374.26 -679.33 -53.7573 C14 Z19650139 Closed state -8.701 -9.111 -7.606 -13429.61 -679.09 -60.4939 C15 Z2400537901 Closed state -10.464 -8.919 -9.831 -13379.26 -678.79 -67.0062 C16 Z2201943136 Open state -8.992 -8.443 -8.72 -13399.65 -678.70 -60.9869 C17 Z2060912694 Open state -8.774 -8.274 -10.44 -13364.38 -678.66 -60.6208 C18 Z2063989764 Open state -9.105 -7.750 -9.56 -13377.59 -678.43 -72.3022 C19 Z1531526443 Open state -8.641 -8.706 -9.55 -13338.36 -676.47 -72.2696 C20 Z2055833084 Open state -8.482 -7.328 -9.19 -13311.35 -674.76 -44.1338 Compound prioritization based on molecular dynamics simulations and binding energy In this study, we conducted MD simulations and binding free energy calculations of 20 protein-ligand complexes obtained from the docking campaigns using the Desmond program for 20 ns. Furthermore, to determine ligand stability over a long period of time, MDS with GROMACS was performed for 100 ns, too. Based on the RMSD values from the GROMACS simulations and the ∆G average values from the 20 ns MD simulations using Desmond (Table 2 ), we identified C7, C9, C12, C14, C15, C16, C18 and C20 compounds as candidates for an extended MD simulation for 300 ns. Subsequently, we assessed the RMSD values obtained from a 300 ns MD simulation (Fig. 4 ). According to the Lig-fit-prot RMSD values in Fig. 4 B, all compounds except C20 were found in the binding cavity at the end of the simulation. Moreover, C16 demonstrated remarkable positional stability, maintaining a 3 Å RMSD value throughout the simulation. According to the MM/GBSA ∆G Distribution graph in Fig. 6 A, both C9 and subsequently C16 displayed favourable binding energies at the last 30 ns of the 300 ns simulation. In the subsequent phase, we selected C9, C15, C16, and C18 for an extensive 1000 ns MDS analysis according to the ∆G and RMSD values, and interaction with the active site residues. The RMSD analysis graph in Fig. 5 B reveals that C16 maintained stability within the range of 2–4 Å. In the case of C9, stability was observed after 400 ns, with a stability range of 4–6 Å. C18 exhibited stability after 150 ns, within a range of 6–8 Å. C15 moving away from the binding pocket, therefore RMSD analysis was not performed. In Fig. 5 C, the RMSF (Root Mean Square Fluctuation) graph indicates that residues between 80 and 100, which include the substrate-specific loop, exhibit dramatic changes for C18 compared to other compounds. Additionally, the residues between 235 and 240, which encompass the opposing loop, show changes for C9 and C16, distinct from those observed for C18. According to the MM/GBSA ∆G Distribution graph in Fig. 6 B, both C9, C16 and C18 displayed favourable binding energies in the last 100 ns of the 1000 ns simulation. The ∆G averages are − 52.5 kcal/mol, -62.1 kcal/mol and − 62.5 kcal/mol for C9, C16 and C18, respectively. Figure 7 illustrates the average energies contributing to the overall MM-GBSA binding energy. Calculations were based on snapshots taken from the last 100 ns of a 1 µs MD simulation. The MM-GBSA ΔG binding energies are quite similar for C16, and C18 (− 62.11 kcal/mol, and − 62.52 kcal/mol, respectively), but the main driving forces behind these binding energies different. For C9-protein binding, lipophilic energy and van der Waals (vdW) energy are the predominant types of energy. In contrast, the contributions of coulomb and hydrogen bond (H-bond) energies in the binding of C16 and C18 are markedly higher from those in C9-protein binding. In Fig. 8 , all three compounds are found in the expected regions after the 1000 ns simulation time. All compounds are positioned within the enzyme's catalytic pocket, which is of crucial importance, and interactions with catalytic residues are evident. C9 retained its significant H-bond interactions after IFD, but lost them after MD. In Supplementary Fig. 6A and 7A, according to BbigLDH-C9 interactions, there are water bridges which interaction mediated between a water molecule and residues including T82, L85, D93, I128, and T129. On the other hand, upon examination of Supplement file, Fig. 6 B and 7 B, it is evident that in the BbigLDH-C16 complex, there are hydrogen bond interactions with residues N17 (77%), L85 (97%), R87 (94%), and N130 (70%), along with water Bridges involving L231 and pi-pi stacking with Y241 (32%). In Supplement file, Fig. 6 C and 7 C, hydrogen bond interactions are observed with residues R87, I128, N130 (46%) and R161 (39%) in the BbigLDH-C18 complex, along with pi-cation interactions involving R161 (33%). In Fig. 9A, after the induced fit docking stage, interactions with R87, I128, and N130 residues were generally observed with all compounds. On the other hand in Fig. 9B, after the 1000 ns simulation times, interactions with R87, N130 conserved. It can be observed that the interaction with these residues is maintained throughout the MD simulation for both C16 and C18 (Supplement file, Fig. 6 ). It is believed that these residues play a significant role in inhibition against BbigLDH. In the cognate docking stage, inhibitors found in the open and closed PfLDH structures, whose PDB structures are known, were docked to the open and closed states of BbigLDH structures. The interaction of the original PDB structures with the inhibitor and the interaction of the BbigLDH structures with the same inhibitor are shown in Supplement file, Fig. 8 . It is observed that the lack of a stable region in the substrate-specific loop region of the open state of PfLDH structures results in differences with the open state of BbigLDH structure. Additionally, it is seen that the residues interacting with the ligands in the open state of BbigLDH structure and also known to have catalytic properties in some cases are R87, I128, N130, R161, and S239. On the other hand, for the closed state of BbigLDH, the most important residues interacting with the ligands are R98, N130, R161, H185, S239, and H240, as shown in Supplement file, Fig. 9. In conclusion, it is considered that compounds showing similarity in the critical residues obtained from the redocking of PfLDH inhibitors, whose PDB structure is known, to BbigLDH, and interacting with residues C09, C16, and C18, may have similar inhibitory properties on BbigLDH. These findings strongly indicate that compounds located in the pocket where the cofactor and substrate bind may inhibit enzyme activity through these specific residue interactions. Particularly, C16 may be a promising candidate for in vitro inhibitor studies, given its stable interaction with BbigLDH. Compounds Z2373310039 (C09), Z2201943136 (C16) and Z2063989764 (C18) may warrant further investigation to comprehend their potential therapeutic effects on babesiosis. Conclusions Recent reports have indicated the emergence of resistance against certain existing anti-babesial drugs. This highlights the urgency of identifying targets for the development of new drugs with distinct modes of action. The lactate dehydrogenase enzyme plays a pivotal role in the anaerobic pathway of the parasite, making it a promising candidate for further drug design endeavours. We employed homology modelling to predict the unknown structure of BbigLDH, followed by structure-based virtual screening. Docking scores, interaction with active site residues, affinity to host LDHs, ADME properties, and stability through simulations were considered during identifying of potential inhibitors. The stability of the three candidate compounds (C9, C16, and C18) was validated through long molecular dynamics simulations and further binding free energy calculations. The obtained results provided insights into the molecular foundation of interactions and, consequently, the mechanism of inhibition. Taken together, the comprehensive analysis of these findings underscores the utility of structure-based virtual screening in providing crucial insights into the essential structural and binding features necessary for the design of novel BbigLDH inhibitors, offering potential avenues for the treatment of bovine babesiosis. Declarations Ethical Approval Not applicable. Funding This work has been supported by Marmara University-Scientific Research Projects Coordination Unit under grant number: FDK-2022-10061. Author contributions Safiye Merve Bostancioglu: Investigation; Visualization; Writing-Original Draft. Ozal Mutlu: Conceptualization; Methodology; Investigation; Visualization; Writing-Original Draft. Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Disclosure statement No potential conflict of interest was reported by the authors. Availability of data and materials No datasets were generated or analysed during the current study. ORCID Safiye Merve Bostancioglu https://orcid.org/0000-0002-0239-3966 Ozal Mutlu https://orcid.org/0000-0003-4551-5780 References T. Smith and F. L. Kilborne, Investigations into the nature, causation, and prevention of Texas or southern cattle fever / . Washington, D.C. : U.S. Dept. of Agriculture, Bureau of Animal Industry, 1893. doi: 10.5962/bhl.title.124068. L. Schnittger, A. E. Rodriguez, M. Florin-Christensen, and D. A. Morrison, ‘Babesia: A world emerging’, Infect. Genet. Evol. , vol. 12, no. 8, pp. 1788–1809, Dec. 2012, doi: 10.1016/j.meegid.2012.07.004. R. Bock, L. Jackson, A. De Vos, and W. Jorgensen, ‘Babesiosis of cattle’, Parasitology , vol. 129, no. S1, pp. S247–S269, Oct. 2004, doi: 10.1017/S0031182004005190. M. Calvopiña et al. , ‘Case report: First report on human infection by tick-borne Babesia bigemina in the Amazon region of Ecuador’, Front. Public Health , vol. 11, p. 1079042, Aug. 2023, doi: 10.3389/fpubh.2023.1079042. L. He, R. G. Bastos, L. Yu, J. M. Laughery, and C. E. Suarez, ‘Lactate Dehydrogenase as a Potential Therapeutic Drug Target to Control Babesia bigemina’, Front. Cell. Infect. Microbiol. , vol. 12, p. 870852, Apr. 2022, doi: 10.3389/fcimb.2022.870852. J. Mosqueda, A. Olvera-Ramirez, G. Aguilar-Tipacamu, and G. J. Canto, ‘Current Advances in Detection and Treatment of Babesiosis’, Curr. Med. Chem. , vol. 19, no. 10, pp. 1504–1518, Apr. 2012, doi: 10.2174/092986712799828355. J. H. D. Silva et al. , ‘Chemoprophylaxis for babesiosis and anaplasmosis in cattle: case report’, Rev. Bras. Parasitol. Veterinária , vol. 29, no. 4, p. e010520, 2020, doi: 10.1590/s1984-29612020096. N. D. Shah, A. U. Bhikane, R. K. Jadhav, S. G. Chavhan, and A. Mohan, ‘Therapeutic management of Babesiosis alone and its mixed infection with Theileriosis in Cattle’, 2019. B. Tuvshintulga, T. Sivakumar, N. Yokoyama, and I. Igarashi, ‘Development of unstable resistance to diminazene aceturate in Babesia bovis’, Int. J. Parasitol. Drugs Drug Resist. , vol. 9, pp. 87–92, Apr. 2019, doi: 10.1016/j.ijpddr.2019.02.001. Todorovic R. A., ‘Chemourophylaxis (Imidocarb. Against Babesia bigemina and Babesia argentina Infections’, Am J Vet Res , vol. 34(9):1153-61, Sep. 1973, doi: PMID: 4747036. N. G. Coldham et al. , ‘Imidocarb residues in edible bovine tissues and in vitro assessment of imidocarb metabolism and cytotoxicity’, Drug Metab. Dispos. Biol. Fate Chem. , vol. 23, no. 4, pp. 501–505, Apr. 1995. H. J. Vial and A. Gorenflot, ‘Chemotherapy against babesiosis’, Vet. Parasitol. , vol. 138, no. 1–2, pp. 147–160, May 2006, doi: 10.1016/j.vetpar.2006.01.048. A. Zintl, G. Mulcahy, H. E. Skerrett, S. M. Taylor, and J. S. Gray, ‘ Babesia divergens , a Bovine Blood Parasite of Veterinary and Zoonotic Importance’, Clin. Microbiol. Rev. , vol. 16, no. 4, pp. 622–636, Oct. 2003, doi: 10.1128/CMR.16.4.622-636.2003. C. E. Suarez and S. Noh, ‘Emerging perspectives in the research of bovine babesiosis and anaplasmosis’, Vet. Parasitol. , vol. 180, no. 1–2, pp. 109–125, Aug. 2011, doi: 10.1016/j.vetpar.2011.05.032. N. McHARDY, R. M. Woollon, R. B. Clampitt, J. A. James, and R. J. Crawley, ‘Efficacy, toxicity and metabolism of imidocarb dipropionate in the treatment of Babesia ovis infection in sheep’, Res. Vet. Sci. , vol. 41, no. 1, pp. 14–20, Jul. 1986, doi: 10.1016/S0034-5288(18)30565-4. C. A. Homewood and K. D. Neame, ‘Conversion of glucose to lactate by intraerythrocytic Plasmodium berghei ’, Ann. Trop. Med. Parasitol. , vol. 77, no. 2, pp. 127–129, Apr. 1983, doi: 10.1080/00034983.1983.11811686. K. L. Kavanagh, R. A. Elling, and D. K. Wilson, ‘Structure of Toxoplasma gondii LDH1: Active-Site Differences from Human Lactate Dehydrogenases and the Structural Basis for Efficient APAD + Use , ’ Biochemistry , vol. 43, no. 4, pp. 879–889, Feb. 2004, doi: 10.1021/bi035108g. C. R. Dunn et al. , ‘The structure of lactate dehydrogenase from Plasmodium falciparum reveals a new target for anti-malarial design’, Nat. Struct. Biol. , vol. 3, no. 11, pp. 912–915, Nov. 1996, doi: 10.1038/nsb1196-912. J. I. Boucher, J. R. Jacobowitz, B. C. Beckett, S. Classen, and D. L. Theobald, ‘An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases’, eLife , vol. 3, p. e02304, Jun. 2014, doi: 10.7554/eLife.02304. S. Bork, M. Okamura, S. Boonchit, H. Hirata, N. Yokoyama, and I. Igarashi, ‘Identification of Babesia bovis l-lactate dehydrogenase as a potential chemotherapeutical target against bovine babesiosis’, Mol. Biochem. Parasitol. , vol. 136, no. 2, pp. 165–172, Aug. 2004, doi: 10.1016/j.molbiopara.2004.03.009. D. Turgut-Balik et al. , ‘Cloning, sequence and expression of the lactate dehydrogenase gene from the human malaria parasite, Plasmodium vivax’, Biotechnol. Lett. , vol. 26, no. 13, pp. 1051–1055, Jul. 2004, doi: 10.1023/B:BILE.0000032958.78158.10. O. Mutlu and D. Turgut Balik, ‘Kinetic Analysis of the Amino Terminal End of Active Site Loop of Lactate Deyhdrogenase from Plasmodium Vivax’, Balk. Med. J. , vol. 29, no. 4, pp. 364–369, Dec. 2012, doi: 10.5152/balkanmedj.2012.047. A. Erdemir and O. Mutlu, ‘Functional and structural characterization of the pentapeptide insertion of Theileria annulata lactate dehydrogenase by site-directed mutagenesis, comparative modeling and molecular dynamics simulations’, J. Mol. Graph. Model. , vol. 74, pp. 344–351, Jun. 2017, doi: 10.1016/j.jmgm.2017.04.010. L. Yu, Q. Liu, W. Luo, J. Zhao, H. F. Alzan, and L. He, ‘The Structural Basis of Babesia orientalis Lactate Dehydrogenase’, Front. Cell. Infect. Microbiol. , vol. 11, p. 790101, Jan. 2022, doi: 10.3389/fcimb.2021.790101. M. S. Gomez et al. , ‘Substrate and cofactor specificity and selective inhibition of lactate dehydrogenase from the malarial parasite P. falciparum’, Mol. Biochem. Parasitol. , vol. 90, no. 1, pp. 235–246, Dec. 1997, doi: 10.1016/S0166-6851(97)00140-0. F. Al-Anouti, S. Tomavo, S. Parmley, and S. Ananvoranich, ‘The Expression of Lactate Dehydrogenase Is Important for the Cell Cycle of Toxoplasma gondii’, J. Biol. Chem. , vol. 279, no. 50, pp. 52300–52311, Dec. 2004, doi: 10.1074/jbc.M409175200. D. Madern, ‘Evolution of Cryptosporidium parvum Lactate Dehydrogenase from Malate Dehydrogenase by a Very Recent Event of Gene Duplication’, Mol. Biol. Evol. , vol. 21, no. 3, pp. 489–497, Dec. 2003, doi: 10.1093/molbev/msh042. L. Yu et al. , ‘Crystal structures of Babesia microti lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis’, FASEB J. , vol. 33, no. 12, pp. 13669–13682, Dec. 2019, doi: 10.1096/fj.201901259R. S. Akash et al. , ‘Anti-parasitic drug discovery against Babesia microti by natural compounds: an extensive computational drug design approach’, Front. Cell. Infect. Microbiol. , vol. 13, p. 1222913, Aug. 2023, doi: 10.3389/fcimb.2023.1222913. E. Lionta, G. Spyrou, D. Vassilatis, and Z. Cournia, ‘Structure-Based Virtual Screening for Drug Discovery: Principles, Applications and Recent Advances’, Curr. Top. Med. Chem. , vol. 14, no. 16, pp. 1923–1938, Oct. 2014, doi: 10.2174/1568026614666140929124445. K. Tamura, G. Stecher, and S. Kumar, ‘MEGA11: Molecular Evolutionary Genetics Analysis Version 11’, Mol. Biol. Evol. , vol. 38, no. 7, pp. 3022–3027, Jun. 2021, doi: 10.1093/molbev/msab120. B. Webb and A. Sali, ‘Comparative Protein Structure Modeling Using MODELLER’, Curr. Protoc. Bioinforma. , vol. 54, no. 1, Jun. 2016, doi: 10.1002/cpbi.3. E. F. Pettersen et al. , ‘UCSF Chimera—A visualization system for exploratory research and analysis’, J. Comput. Chem. , vol. 25, no. 13, pp. 1605–1612, Oct. 2004, doi: 10.1002/jcc.20084. C. Colovos and T. O. Yeates, ‘Verification of protein structures: Patterns of nonbonded atomic interactions’, Protein Sci. , vol. 2, no. 9, pp. 1511–1519, Sep. 1993, doi: 10.1002/pro.5560020916. M. J. Sippl, ‘Recognition of errors in three-dimensional structures of proteins’, Proteins Struct. Funct. Genet. , vol. 17, no. 4, pp. 355–362, Dec. 1993, doi: 10.1002/prot.340170404. M. Wiederstein and M. J. Sippl, ‘ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins’, Nucleic Acids Res. , vol. 35, no. Web Server, pp. W407–W410, May 2007, doi: 10.1093/nar/gkm290. R. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, ‘PROCHECK: a program to check the stereochemical quality of protein structures’, J. Appl. Crystallogr. , vol. 26, no. 2, pp. 283–291, Apr. 1993, doi: 10.1107/S0021889892009944. J. U. Bowie, R. Lüthy, and D. Eisenberg, ‘A Method to Identify Protein Sequences That Fold into a Known Three-Dimensional Structure’, Science , vol. 253, no. 5016, pp. 164–170, Jul. 1991, doi: 10.1126/science.1853201. P. Benkert, S. C. E. Tosatto, and D. Schomburg, ‘QMEAN: A comprehensive scoring function for model quality assessment’, Proteins Struct. Funct. Bioinforma. , vol. 71, no. 1, pp. 261–277, Apr. 2008, doi: 10.1002/prot.21715. J. Jumper et al. , ‘Highly accurate protein structure prediction with AlphaFold’, Nature , vol. 596, no. 7873, pp. 583–589, Aug. 2021, doi: 10.1038/s41586-021-03819-2. R. A. Friesner et al. , ‘Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes’, J. Med. Chem. , vol. 49, no. 21, pp. 6177–6196, Oct. 2006, doi: 10.1021/jm051256o. Schrödinger Release 2022-4: Glide, Schrödinger, LLC, New York, NY, 2022 . J. A. Read, K. W. Wilkinson, R. Tranter, R. B. Sessions, and R. L. Brady, ‘Chloroquine Binds in the Cofactor Binding Site ofPlasmodium falciparum Lactate Dehydrogenase’, J. Biol. Chem. , vol. 274, no. 15, pp. 10213–10218, Apr. 1999, doi: 10.1074/jbc.274.15.10213. A. Cameron et al. , ‘Identification and Activity of a Series of Azole-based Compounds with Lactate Dehydrogenase-directed Anti-malarial Activity’, J. Biol. Chem. , vol. 279, no. 30, pp. 31429–31439, Jul. 2004, doi: 10.1074/jbc.M402433200. R. Conners et al. , ‘Mapping the binding site for gossypol-like inhibitors of Plasmodium falciparum lactate dehydrogenase’, Mol. Biochem. Parasitol. , vol. 142, no. 2, pp. 137–148, Aug. 2005, doi: 10.1016/j.molbiopara.2005.03.015. H. M. Berman, ‘The Protein Data Bank’, Nucleic Acids Res. , vol. 28, no. 1, pp. 235–242, Jan. 2000, doi: 10.1093/nar/28.1.235. Schrödinger Release 2022-4: Protein Preparation Wizard; Epik, Schrödinger, LLC, New York, NY, 2022; Impact, Schrödinger, LLC, New York, NY; Prime, Schrödinger, LLC, New York, NY, 2022 . Schrödinger Release 2022-4: LigPrep, Schrödinger, LLC, New York, NY, 2022 . Schrödinger Release 2022-4: Maestro, Schrödinger, LLC, New York, NY, 2022 . R. A. Friesner et al. , ‘Glide: A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy’, J. Med. Chem. , vol. 47, no. 7, pp. 1739–1749, Mar. 2004, doi: 10.1021/jm0306430. A. Chang et al. , ‘BRENDA, the ELIXIR core data resource in 2021: new developments and updates’, Nucleic Acids Res. , vol. 49, no. D1, pp. D498–D508, Jan. 2021, doi: 10.1093/nar/gkaa1025. J. C. Shelley, A. Cholleti, L. L. Frye, J. R. Greenwood, M. R. Timlin, and M. Uchimaya, ‘Epik: a software program for pK a prediction and protonation state generation for drug-like molecules’, J. Comput. Aided Mol. Des. , vol. 21, no. 12, pp. 681–691, Dec. 2007, doi: 10.1007/s10822-007-9133-z. G. Madhavi Sastry, M. Adzhigirey, T. Day, R. Annabhimoju, and W. Sherman, ‘Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments’, J. Comput. Aided Mol. Des. , vol. 27, no. 3, pp. 221–234, Mar. 2013, doi: 10.1007/s10822-013-9644-8. Schrödinger Release 2022-4: QikProp, Schrödinger, LLC, New York, NY, 2022-4 . C. A. Lipinski, ‘Drug-like properties and the causes of poor solubility and poor permeability’, J. Pharmacol. Toxicol. Methods , vol. 44, no. 1, pp. 235–249, Jul. 2000, doi: 10.1016/S1056-8719(00)00107-6. A. E. Cho, V. Guallar, B. J. Berne, and R. Friesner, ‘Importance of accurate charges in molecular docking: Quantum mechanical/molecular mechanical (QM/MM) approach’, J. Comput. Chem. , vol. 26, no. 9, pp. 915–931, Jul. 2005, doi: 10.1002/jcc.20222. A. D. Bochevarov et al. , ‘Jaguar: A high‐performance quantum chemistry software program with strengths in life and materials sciences’, Int. J. Quantum Chem. , vol. 113, no. 18, pp. 2110–2142, Sep. 2013, doi: 10.1002/qua.24481. Schrödinger Release 2022-4: Jaguar, Schrödinger, LLC, New York, NY, 2022 . J. C. Phillips et al. , ‘Scalable molecular dynamics on CPU and GPU architectures with NAMD’, J. Chem. Phys. , vol. 153, no. 4, p. 044130, Jul. 2020, doi: 10.1063/5.0014475. J. Lee et al. , ‘CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field’, J. Chem. Theory Comput. , vol. 12, no. 1, pp. 405–413, Jan. 2016, doi: 10.1021/acs.jctc.5b00935. N. M. O’Boyle, M. Banck, C. A. James, C. Morley, T. Vandermeersch, and G. R. Hutchison, ‘Open Babel: An open chemical toolbox’, J. Cheminformatics , vol. 3, no. 1, p. 33, Dec. 2011, doi: 10.1186/1758-2946-3-33. W. Humphrey, A. Dalke, and K. Schulten, ‘VMD: Visual molecular dynamics’, J. Mol. Graph. , vol. 14, no. 1, pp. 33–38, Feb. 1996, doi: 10.1016/0263-7855(96)00018-5. K. Vanommeslaeghe et al. , ‘CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields’, J. Comput. Chem. , p. NA-NA, 2009, doi: 10.1002/jcc.21367. Schrödinger Release 2022-4: Desmond Molecular Dynamics System, D. E. Shaw Research, New York, NY, 2022 . M. P. Jacobson, R. A. Friesner, Z. Xiang, and B. Honig, ‘On the Role of the Crystal Environment in Determining Protein Side-chain Conformations’, J. Mol. Biol. , vol. 320, no. 3, pp. 597–608, Jul. 2002, doi: 10.1016/S0022-2836(02)00470-9. Naruya Saitou and M. Nei, ‘The neighbor-joining method: a new method for reconstructing phylogenetic trees.’, Mol. Biol. Evol. , Jul. 1987, doi: 10.1093/oxfordjournals.molbev.a040454. J. I. Boucher, J. R. Jacobowitz, B. C. Beckett, S. Classen, and D. L. Theobald, ‘An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases’, eLife , vol. 3, p. e02304, Jun. 2014, doi: 10.7554/eLife.02304. K. L. Kavanagh, R. A. Elling, and D. K. Wilson, ‘Structure of Toxoplasma gondii LDH1: Active-Site Differences from Human Lactate Dehydrogenases and the Structural Basis for Efficient APAD + Use , ’ Biochemistry , vol. 43, no. 4, pp. 879–889, Feb. 2004, doi: 10.1021/bi035108g. C. Colovos and T. O. Yeates, ‘Verification of protein structures: patterns of nonbonded atomic interactions’, Protein Sci. Publ. Protein Soc. , vol. 2, no. 9, pp. 1511–1519, Sep. 1993, doi: 10.1002/pro.5560020916. M. Wiederstein and M. J. Sippl, ‘ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins’, Nucleic Acids Res. , vol. 35, no. Web Server issue, pp. W407-410, Jul. 2007, doi: 10.1093/nar/gkm290. P. Benkert, S. C. E. Tosatto, and D. Schomburg, ‘QMEAN: A comprehensive scoring function for model quality assessment’, Proteins , vol. 71, no. 1, pp. 261–277, Apr. 2008, doi: 10.1002/prot.21715. J. U. Bowie, R. Lüthy, and D. Eisenberg, ‘A method to identify protein sequences that fold into a known three-dimensional structure’, Science , vol. 253, no. 5016, pp. 164–170, Jul. 1991, doi: 10.1126/science.1853201. R. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, ‘PROCHECK: a program to check the stereochemical quality of protein structures’, J. Appl. Crystallogr. , vol. 26, no. 2, pp. 283–291, Apr. 1993, doi: 10.1107/S0021889892009944. A. Volkamer, D. Kuhn, F. Rippmann, and M. Rarey, ‘DoGSiteScorer: a web server for automatic binding site prediction, analysis and druggability assessment’, Bioinformatics , vol. 28, no. 15, pp. 2074–2075, Aug. 2012, doi: 10.1093/bioinformatics/bts310. W. L. Jorgensen and E. M. Duffy, ‘Prediction of drug solubility from Monte Carlo simulations’, Bioorg. Med. Chem. Lett. , vol. 10, no. 11, pp. 1155–1158, Jun. 2000, doi: 10.1016/S0960-894X(00)00172-4. A. Cavalli, E. Poluzzi, F. De Ponti, and M. Recanatini, ‘Toward a Pharmacophore for Drugs Inducing the Long QT Syndrome: Insights from a CoMFA Study of HERG K + Channel Blockers’, J. Med. Chem. , vol. 45, no. 18, pp. 3844–3853, Aug. 2002, doi: 10.1021/jm0208875. M. Yazdanian, S. L. Glynn, J. L. Wright, and A. Hawi, ‘Correlating Partitioning and Caco-2 Cell Permeability of Structurally Diverse Small Molecular Weight Compounds’, Pharm. Res. , vol. 15, no. 9, pp. 1490–1494, 1998, doi: 10.1023/A:1011930411574. J. M. Luco, ‘Prediction of the Brain−Blood Distribution of a Large Set of Drugs from Structurally Derived Descriptors Using Partial Least-Squares (PLS) Modeling’, J. Chem. Inf. Comput. Sci. , vol. 39, no. 2, pp. 396–404, Mar. 1999, doi: 10.1021/ci980411n. G. Colmenarejo, A. Alvarez-Pedraglio, and J.-L. Lavandera, ‘Cheminformatic Models To Predict Binding Affinities to Human Serum Albumin’, J. Med. Chem. , vol. 44, no. 25, pp. 4370–4378, Dec. 2001, doi: 10.1021/jm010960b. Additional Declarations No competing interests reported. Supplementary Files Supplementaryfiles.docx Cite Share Download PDF Status: Published Journal Publication published 07 Jan, 2025 Read the published version in Parasitology Research → Version 1 posted Editorial decision: Revision requested 01 Sep, 2024 Editor assigned by journal 27 Aug, 2024 Submission checks completed at journal 26 Aug, 2024 First submitted to journal 21 Aug, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4951492","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":347926943,"identity":"2db672a3-62c4-43ac-abab-320013d5a4a7","order_by":0,"name":"Safiye Merve Bostancioglu","email":"","orcid":"","institution":"Marmara University","correspondingAuthor":false,"prefix":"","firstName":"Safiye","middleName":"Merve","lastName":"Bostancioglu","suffix":""},{"id":347926944,"identity":"231a43e6-7da2-4414-8b3d-f3e2392ca4e0","order_by":1,"name":"Ozal Mutlu","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA+0lEQVRIiWNgGAWjYHACxgMgko+B+QBCjIeAHpBSCTYGtgQULRJEaOExIE4LP//hBwcYcw7XsUnkfN3wc882efP2A4wP3rYx1Jk3YNci2XDM4ADjtsMSbBK52272PLttOOdMArPh3DYGCZkD2LUYHGxAaLnBc+A24wyGBDZpXqAWXC4zOMz+Aaol59nNPwdu28/gf8D+G6+WYzwwW3LYbgNtSZwhkcDGjE+LZA9PwYHEbemSbTzPzG7LHLidPEPiYbPknHMSkjNwhtjxjQ8+brPm52dPfnbzzYHbtjP4kw9+eFNmw483YhJQuYwNDARichSMglEwCkYBAQAAKMFZIY/6f98AAAAASUVORK5CYII=","orcid":"","institution":"Marmara University","correspondingAuthor":true,"prefix":"","firstName":"Ozal","middleName":"","lastName":"Mutlu","suffix":""}],"badges":[],"createdAt":"2024-08-21 12:05:55","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4951492/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4951492/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s00436-024-08433-5","type":"published","date":"2025-01-07T15:58:03+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":65162245,"identity":"0fbf38b4-7431-4581-8b85-a07aded7caed","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":54399,"visible":true,"origin":"","legend":"\u003cp\u003eWorkflows for model generation (orange boxes), structure-based virtual screening (green boxes), molecular dynamics simulations (purple boxes).\u003c/p\u003e","description":"","filename":"image1.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/2dbf49328f4c191365b04adf.png"},{"id":65163400,"identity":"3f8793ba-8834-4e5c-b18c-8e652d7dfc4b","added_by":"auto","created_at":"2024-09-24 09:27:49","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":341721,"visible":true,"origin":"","legend":"\u003cp\u003eSequence alignment of the substrate specific loop (residues 93 to 97, BbigLDH numbering) \u003cstrong\u003e(A)\u003c/strong\u003e and the opposing loop (residues 236 to 238, BbigLDH numbering) \u003cstrong\u003e(B).\u003c/strong\u003e\u003c/p\u003e","description":"","filename":"image2.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/7439969245cdb619aa44eb51.png"},{"id":65162248,"identity":"7d404ec0-feca-4ea9-a117-d8990fd93d8c","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":954137,"visible":true,"origin":"","legend":"\u003cp\u003e3-D structure diagram of (A) BbigLDH open (light gray) and closed (light orange) forms. The 5-AA insertion region is indicated in black for the open form and in orange for the closed form. Additionally, the opposite loop is indicated in dark gray for the open form and in yellow for the closed form. (B) The open form of BbigLDH, host LDH A, and LDH B is indicated in light gray, yellow, and light blue, respectively. The 5-AA insertion region is indicated in black for the open form of BbigLDH. Important loop regions are highlighted in yellow boxes.\u003c/p\u003e","description":"","filename":"image3.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/fa304f16bfa936f819cb929d.png"},{"id":65162247,"identity":"a2ffab4c-63cd-4c71-9582-20bfa6a8743e","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":152915,"visible":true,"origin":"","legend":"\u003cp\u003eRMSD graphs following the 300 ns molecular dynamics simulation. (A) C-alpha RMSD (Å)/ns. (B) Lig-fit-prot RMSD (Å)/ns.\u003c/p\u003e","description":"","filename":"image4.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/37eb3555cc450333137efd7d.png"},{"id":65162254,"identity":"7a74adeb-1242-4b70-b334-a896e709b868","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":144610,"visible":true,"origin":"","legend":"\u003cp\u003eRMSD graphs following the 1000 ns molecular dynamics simulation. (A) C-alpha backbone RMSD (Å)/ns. (B) Lig-fit-prot RMSD (Å)/ns. (C) C-alpha backbone RMSF for protein.\u003c/p\u003e","description":"","filename":"image5.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/69a2d19312b4dcbb9e032a22.png"},{"id":65163399,"identity":"468c8801-1333-4ee3-a435-7b4a98fd00f9","added_by":"auto","created_at":"2024-09-24 09:27:49","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":48152,"visible":true,"origin":"","legend":"\u003cp\u003e(A) MM/GBSA ∆G distribution graph obtained for the last 30 ns of 300ns MD simulation. (B) MM/GBSA ∆G distribution graph obtained for the last 100 ns of 1000ns MD simulation.\u003c/p\u003e","description":"","filename":"image6.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/61e09f26324da129918fc015.png"},{"id":65162246,"identity":"74660449-6867-4707-99b5-7aa3da7953a2","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":34707,"visible":true,"origin":"","legend":"\u003cp\u003eAverage MM/GBSA ΔG values for compounds C9, C16 and C18 and contributions of the various energies (kcal/mol).\u003c/p\u003e","description":"","filename":"image7.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/1dcf976663c6b16ea206db5b.png"},{"id":65162253,"identity":"c677320c-b6b1-415b-b663-0aed0b560e51","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":808255,"visible":true,"origin":"","legend":"\u003cp\u003eBinding interactions of 3-D structures of BbigLDH with C9, C16 and C18 after (A) IFD (indicated gray) and (B) 1000 ns simulation time (indicated green). H-bond interactions are indicated in yellow, pi-pi stacking indicated in orange and compounds have interactions with residues in red. (C) 2-D structures of compounds C9, C16 an C18.\u003c/p\u003e","description":"","filename":"image8.png","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/537d92f61f43505af212eb14.png"},{"id":73694472,"identity":"416f66b2-89b1-44ea-af98-3c654ef7abf9","added_by":"auto","created_at":"2025-01-13 16:13:19","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4044279,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/bf7b6265-d20c-4b7e-a2de-35c22d39cad4.pdf"},{"id":65162250,"identity":"ef54f6b5-d787-4a2d-9d41-d8c199e7998d","added_by":"auto","created_at":"2024-09-24 09:19:49","extension":"docx","order_by":12,"title":"","display":"","copyAsset":false,"role":"supplement","size":5105698,"visible":true,"origin":"","legend":"","description":"","filename":"Supplementaryfiles.docx","url":"https://assets-eu.researchsquare.com/files/rs-4951492/v1/86e62378697b287d055f2614.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Exploring Novel Inhibitors for Babesia bigemina Lactate Dehydrogenase: A Computational Structural Biology Perspective","fulltext":[{"header":"Introduction","content":"\u003cp\u003eSmith and Kilborne (1893) described the 'Texas fever' for the first time, known as babesiosis, which is caused by the apicomplexan parasite \u003cem\u003eBabesia bigemina\u003c/em\u003e and transmitted through tick vectors [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. The disease manifests in a spectrum of states ranging from asymptomatic carriers to severe cases marked by hemolytic anemia, fever, hemoglobinuria, and, in some instances, fatalities. Its impact extends across a diverse range of mammals, including domestic pets, farm animals, wildlife, and even humans [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e], [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e], [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. \u003cem\u003eB. bigemina\u003c/em\u003e currently prevail in tropical and subtropical regions worldwide, spanning continents such as Africa, Asia, the Americas, Southern Europe, and Australia [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Their presence directly influences meat and milk production, thereby influencing the competitive edge of livestock industries [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. The primary vectors for these pathogens are \u003cem\u003eRhipicephalus (Boophilus) microplus\u003c/em\u003e, \u003cem\u003eR. annulatus, R. australis, R. decoloratus, R. bursa, R. geigyi\u003c/em\u003e, and \u003cem\u003eR. evertsi\u003c/em\u003e [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e], [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Currently, diminazene aceturate (DA) and imidocarb dipropionate stand as the most commonly employed treatments [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e], [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e], [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]. However, the continued or inappropriate use of these drugs poses a significant risk of engendering drug resistance [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e], [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e], [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e], and they also carry the potential for chemical residues in meat and milk, which can have adverse effects on human health [\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e], [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e], [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. This pressing concern necessitates an urgent quest for the discovery or development of alternative, effective, and economically viable drugs with minimal toxicity [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e], [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eWithin apicomplexans, lactate dehydrogenase (LDH) serves as a key glycolytic enzyme, facilitating the transformation of pyruvate into lactate. [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. It bidirectionally catalyzes the conversion of pyruvate to lactate, utilizing NADH as a co-factor and regenerating NAD\u0026thinsp;+\u0026thinsp;essential for glycolysis to proceed, or the conversion of lactate to pyruvate, utilizing NAD\u0026thinsp;+\u0026thinsp;as a co-factor and regenerating NADH [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]. The energy generated through this process is harnessed by parasites to fuel their biochemical processes and ensure their survival. The substrate specific loop, which is present in the active site, is designated for its substrate in parasitic LDH. Hosts represent canonical LDH which absence of five amino acids in the substrate specific loop differs from parasitic LDH. The five amino insertion was firstly described in LDH of \u003cem\u003ePlasmodium falciparum\u003c/em\u003e [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Further studies showed that the insertion site is conserved in ancient LDH and LDH of many apicomplexan parasites including \u003cem\u003eB. bigemina\u003c/em\u003e, but in a different amino acid sequence [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e], [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e], [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e], [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e], [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e], [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e], [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e], [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. The vital importance of the enzyme, and the fundamental difference between the host's LDH and the parasite's LDH in the active site makes it a more specific drug target. Further investigation of \u003cem\u003eP. falciparum\u003c/em\u003e LDH demonstrates the potential for identifying a specific inhibitor for the enzyme [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. Other studies related to LDH in apicomplexan parasites were published for \u003cem\u003eT. gondii, C. parvum\u003c/em\u003e, and supports the hypothesis that apicomplexan LDH could be used as a drug target [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e], [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. The crystal structure of \u003cem\u003eBabesia microti\u003c/em\u003e LDH (BmLDH) is the first LDH whose crystal structure has been solved in babesia genus. However, it is important to note that the BmLDH shows significant differences from apicomplexans LDHs and its structure is closely related to mammalian-like LDH [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e], [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e]. Yu et al., 2022 was identified a novel cDNA clone that encodes apicomplexan type LDH from \u003cem\u003eBabesia orientalis\u003c/em\u003e (BoLDH) and solved its crystal structure in its apo form [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. On the other hand, gossypol, an LDH inhibitor, has been explored for its impact on various Babesia species. A concentration with an IC50 value of 50 \u0026micro;M gossypol irreversibly halted the in vitro growth of \u003cem\u003eB. bovis.\u003c/em\u003e Crucially, the addition of gossypol led to a noteworthy inhibition of enzymatic activity, acting as a competitive inhibitor for the binding of NADH [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. Another study for gossypol, exhibits inhibitory impact on the \u003cem\u003ein vitro\u003c/em\u003e growth of \u003cem\u003eB. bigemina\u003c/em\u003e, with an IC50 of 43.97 \u0026micro;M after 72 hours of treatment. The maximum inhibitory concentration (IC98) is observed at 60 \u0026micro;M gossypol [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. These findings demonstrate that LDH is a significant drug target for Babesias as an apicomplexan parasite.\u003c/p\u003e \u003cp\u003eModern progress in computer-assisted drug development has made it possible to design drugs specifically for target proteins, thereby accelerating the creation of new medications. Structure-Based Virtual Screening has notably transformed and improved the process of drug discovery, establishing itself as one of the most effective computational techniques for drug design [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. The principles and applications of Virtual Screening encompass the entire process, from the initial stages of receptor and library preparation, through docking and scoring, to the post-processing of the highest-scoring hits. This study focused on understanding the structure of \u003cem\u003eBabesia bigemina\u003c/em\u003e LDH (BbigLDH) enzyme and discovering candidate inhibitors by structure-based drug discovery methods. After identifying the hit compounds using virtual screening, molecular dynamics simulations were performed to assess protein-ligand complex stability and calculate end-point binding free energies. This research could contribute to the identification of potential LDH inhibitors for the treatment of \u003cem\u003eBabesia bigemina\u003c/em\u003e infection and to improve drug design efforts against bovine babesiosis and diseases caused by other apicomplexans.\u003c/p\u003e"},{"header":"Materials \u0026 Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eSequence Analysis\u003c/h2\u003e \u003cp\u003eProtein sequences of \u003cem\u003eBos taurus\u003c/em\u003e LDH B (CAB96751.1), \u003cem\u003eBos taurus\u003c/em\u003e LDH A (BAA14170.1), \u003cem\u003eBabesia bovis\u003c/em\u003e T2Bo (XP_001611047.1), \u003cem\u003eBabesia bigemina\u003c/em\u003e (XP_012766822.1), \u003cem\u003eBabesia ovata\u003c/em\u003e (XP_028866450.1), \u003cem\u003eTheileria equi\u003c/em\u003e strain WA (XP_004833698.1), \u003cem\u003eTheileria orientalis\u003c/em\u003e strain Shintoku (XP_009689964.1), \u003cem\u003eTheileria parva\u003c/em\u003e strain Muguga (XP_766703.1), \u003cem\u003eTheileria annulata\u003c/em\u003e (ADG45564.1), \u003cem\u003ePlasmodium malariae\u003c/em\u003e (SBS88125.1), \u003cem\u003ePlasmodium yoelii\u003c/em\u003e (XP_724101.1), \u003cem\u003ePlasmodium falciparum\u003c/em\u003e 3D7 (XP_001349989.1), \u003cem\u003eToxoplasma gondii\u003c/em\u003e ME49 (XP_002368488.1) and \u003cem\u003eEimeria acervulina\u003c/em\u003e (ACM77785.1) were obtained from the NCBI (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). ClustalW (1.2.4) was used for multiple sequence alignment and the aligned sequences were visualized via ESPript 3.0 webserver [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. The phylogenetic tree which is assessed by using the Maximum Likelihood method was constructed via MEGA 11 [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e].\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eHomology Modelling and Model Validations\u003c/h2\u003e \u003cp\u003eBbigLDH, whose experimentally verified 3-D structure is unknown, were generated by homology modelling using MODELLER 9.23 software [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. Template for homology modelling was determined using NCBI BLAST tool [Basic Local Alignment Seach Tool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.ncbi.nlm.nih.gov/blast\u003c/span\u003e\u003cspan address=\"http://www.ncbi.nlm.nih.gov/blast\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)]. Resulting PDB structures with the highest max score, the lowest E-value, and the highest % identity score were evaluated. \u003cem\u003eT. gondii\u003c/em\u003e LDH1 apo form (1PZE) [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e] and ancestral apicomplexan LDH with malate (4PLC) [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e] were chosen as templates for the modelling of the open and closed states of the BbigLDH enzyme, respectively. A total of 100 models were generated, and the model with the lowest DOPE (Discrete Optimized Protein Energy) and z-DOPE (normalized DOPE) scores was selected for further in silico analyses. UCSF Chimera 1.15 was employed to energy minimize the BbigLDH models, adding any missing hydrogen atoms and charges [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. The minimization process involved 2000 steps of \"Steepest Descent\" followed by 1000 steps of \"Conjugate Gradient,\" utilizing the AmberFF14SB force field.\u003c/p\u003e \u003cp\u003eQuality of the protein models was evaluated by using servers including ERRAT [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e], ProSA [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e], [\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e], ProCHECK / Ramanchandran Plot [\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e], Verify3D [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e] and QMEAN [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. To assess the reliability of specific atomic interactions, ERRAT was utilized. ProSA was employed to verify the correctness of the overall 3D homology model structure of BbigLDH and ensure proper energy optimization. Verify3D was used for 3D structure verification by computing position-specific 3D-1D scores for each residue, referencing a database of known accurate folding structures. The QMEAN tool provided a quantitative evaluation of significant geometric features in the protein structure. The Ramachandran Plot (ProCHECK) was used to visualize the φ (phi) and ψ (psi) dihedral angles of the amino acid residues, identifying conformational outliers and confirming that these angles are within the favorable regions, which are indicative of proper folding and stability.\u003c/p\u003e \u003cp\u003eThe LDH structures of the host (\u003cem\u003eBos taurus\u003c/em\u003e) was obtained from the Alphafold database (Uniprod ID: (LDHA - P19858, LDHB - Q5E9B1, and LDHAL6B - Q3T056) [\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e] (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://alphafold.ebi.ac.uk/\u003c/span\u003e\u003cspan address=\"https://alphafold.ebi.ac.uk/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eCognate docking and Enrichment Calculations\u003c/h2\u003e \u003cp\u003eCognate docking was performed using the Glide XP (extra precision) scoring function and docking methodology to asses prediction quality of the docking program [\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e], [\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e]. The 3-D crystal structures of \u003cem\u003eP. falciparum\u003c/em\u003e LDH (PfLDH) with inhibitors (PDB IDs: 1LDG [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e], 1CET [\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e], 1T26, 1T24, 1T25 [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e], 1U5A, 1U4O, 1XIV, 1U4S [\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]) were obtained from the Protein Data Bank [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e]. The co-crystallized ligands were split from the protein structure, then protein and ligand structures were prepared using the protein preparation workflow respectively and re-docked [\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e], [\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e]. Additionally, these ligands were analyzed for their interactions with BbigLDH using XP docking. Heavy-atom RMSD values were calculated between the docked poses and the crystal structure poses using the \u003cem\u003esuperposition\u003c/em\u003e tool from Maestro [\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e]. In addition, docking methodology was validated by decoys set of 1000 molecules [\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e]. Known 34 LDH inhibitors were obtained from the BRENDA Enzyme database [\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e]. 1000 decoy compounds and 34 active molecules were prepared with LigPrep [\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e] and docked by Glide XP mode with default settings. ROC, RIE, AUAC and BEDROC parameters were analysed by Enrichment Factor Calculation.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eStructure-Based Virtual Screening, QM-Polarized Ligand Docking (QPLD) and Induced Fit Docking (IFD)\u003c/h2\u003e \u003cp\u003eCompound library including about 570.000 compounds from the ChemDiv (Targeted diversity library, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.chemdiv.com/catalog/diversity-libraries/targeted-diversity-library/\u003c/span\u003e\u003cspan address=\"https://www.chemdiv.com/catalog/diversity-libraries/targeted-diversity-library/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), Enamine (Hit Locator Library HLL-460, Discovery Diversity Set DDS-50, Discovery Diversity Set DDS-10, \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://enamine.net/compound-libraries/diversity-libraries\u003c/span\u003e\u003cspan address=\"https://enamine.net/compound-libraries/diversity-libraries\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and OTOVAchemicals (PrimScreen1, PrimScreen2, PrimScreen3, PrimScreen5, PrimScreen10 \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://otavachemicals.com/products/compound-libraries-for-hts/diversity-sets\u003c/span\u003e\u003cspan address=\"https://otavachemicals.com/products/compound-libraries-for-hts/diversity-sets\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) were prepared using the LigPrep [\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e]program. Ionization states of the compounds were predicted using the Epik at 7.5 pH [\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e]. Protein Preparation Workflow was performed to prepare protein structures for docking[\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e], [\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. Receptor grid generation was used to define the atomic positions of the binding site for each structure.\u003c/p\u003e \u003cp\u003eGlide-Virtual Screening Workflow was performed in three steps: High throughput virtual screening (HTVS-10%), Standard precision (SP-10%), and Extra precision (XP- the best docking scores of 100 ligands). Compounds were ranked according to XP-GScore and the best compounds having interactions with the catalytic residues and lower XP-GScores (\u0026le;-7.0 kcal/mol) were selected. Afterwards, ligands are filtered based on their ADME properties [\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. The candidates identified are within the reference values specified in the Qikprop manual. For compounds to qualify as drug candidates, several criteria and appropriate value ranges must be considered. Among these, the key five principles of Lipinski are essential [\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e]. Compounds with acceptable scores were then used in Quantum Mechanical Polarized Ligand Docking (QPLD). The QPLD is a computational method that combines Glide docking with QSite calculations to determine the ligand binding positions [\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. For the initial docking, the Glide XP method is employed with accurate QM level chosen for QM charges using the Jaguar module in gas phase [\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e], [\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e] Redocking is performed using the XP method. In result, ligand conformations with the lowest Emodel values obtained from the docking was selected. Compounds with XP-GScore value equal to or lower than \u0026minus;\u0026thinsp;7.0 kcal/mol were filtered for the induced fit docking. IFD is a molecular docking method that involves ligand repositioning onto a relaxed receptor using Prime. In this approach, a receptor grid was generated by selecting active site residues for molecular docking. Prime was utilized to perform receptor rearrangement for residues within a 5 \u0026Aring; distance range. The relaxed receptor was then subjected to redocking using the XP method of Glide program. The resulting poses were evaluated based on the lowest IFD scores. Compounds identified post-IFD were subjected to Glide XP docking against the structures of \u003cem\u003eBos taurus\u003c/em\u003e LDH A, LDH B and LDH AL6B. Compounds displaying low affinity with the host LDH structures were selected for further analysis, considering them as potential candidates for selective inhibition.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eMolecular Dynamics Simulations and Binding Free Energy Calculations with MM/GBSA\u003c/h2\u003e \u003cp\u003eMD simulation was performed for open and closed states of homology models by using the Nanoscale Molecular Dynamics (NAMD) program [\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e]. During the preparation of closed states, ligands were processed using the CHARMM-GUI server [\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e] with the inclusion of hydrogen atoms facilitated by Open Babel, followed by conversion to the mol2 format [\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e]. The coenzyme (NAI) was then added into the system and parametrized. The protein and ligand PDB and PSF files were generated and combined using VMD 1.9.3 software to create the complex structure [\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e]. The complex structure was solvated in a cubic box and subsequently neutralized, and the ionic strength was adjusted to 150 mM NaCl. Before MD simulations, 2-tier of minimization (total 5000 steps) was conducted. In the first step, constraint with constraint scaling of 100 kcal/mol/A^2 was applied to the solids, while in the second step, the entire system was minimized without any constraints. Equilibration was carried out using the NVT ensemble for 250 ps at 300 K by applying constraints (25 kcal/mol/A^2) to the solids. This was followed by an addition equilibration and relaxation step with the NPT ensemble for 275 ps at 300 K and 1 atm using the Nose\u0026ndash;Hoover Langevin piston. The constraints on the solid were gradually reduced from 10 kcal/mol/A^2 to zero at every 25 ps. Finally, a production phase consisting of a 100 ns MD simulation was performed without any constrains using the NPT ensemble (300 K and 1 atm). All MD simulations were conducted using the CHARMM36 forcefield [\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e]. The MD simulation trajectories were analysed for the RMSD and RMSF values and visualized using the VMD 1.9.3 software [\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eMolecular dynamics simulations of the protein-ligand complexes were performed using the Desmond software (OPLS4 forcefield) [\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e] or GROMACS 2022 (CHARMM36 forcefield) [\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e] to asses ligand stability and to calculate binding energies of compounds. The system consisted of TIP3P water molecules in an orthorhombic box. To ensure charge neutrality and achieve an ionic strength of 150 mM NaCl, appropriate ions were added. Prior to production runs, a minimization step for 100 picoseconds and a default relaxation protocol provided by Desmond was employed to minimize and relax the system. The production runs were conducted for 20, 300, or 1000 ns at the different stages of compound selection procedure for the next stages. All production runs were performed at a temperature of 300 K and a pressure of 1 atmosphere. The integration time step was set as 2 femtoseconds. Temperature and pressure were controlled using the Nose-Hoover thermostat and Martyna-Tobias-Klein (MTK) barostat, respectively. A cutoff distance of 9 \u0026Aring; was applied for nonbonded interactions to calculate the short-range forces efficiently. Protein-ligand interactions were investigated with simulation interaction diagram. MM/GBSA method was employed to calculate the binding free energies (ΔG bind) using the thermal_mmgbsa.py script and the Prime program [\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e]. In parallel, 100 ns of MD simulations were performed using the GROMACS by applying the same parameters as above. Minimization was stopped when the maximum force\u0026thinsp;\u0026lt;\u0026thinsp;10.0 kJ/mol. NVT equilibration (100 ps) was employed by the Particle Mesh Ewald method applied to calculate long range electrostatic interaction and the modified Berendsen thermostat (V-rescale) used to control the temperature of the simulation system. NPT equilibration (100 ps) and production run was performed using the C-rescale pressure coupling at 1 atm and V-rescale used to control the temperature of the simulation system. For the analysis of MD trajectories, GROMACS modules (gmx rms and gmx rmsf) were used, focusing on parameters such as root mean squared deviation (RMSD) and fluctuation (RMSF).\u003c/p\u003e \u003c/div\u003e "},{"header":"Results \u0026 Discussions","content":"\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\n \u003cdiv id=\"Sec9\" class=\"Section3\"\u003e\n \u003ch2\u003eAssessment of LDH sequence similarity among apicomplexan species\u003c/h2\u003e\n \u003cp\u003eUpon analysis of its nucleotide sequence, it was determined that the BbigLDH gene contains an open reading frame spanning 969 nucleotides, which encodes for a sequence of 332 amino acids. The alignment of sequences was conducted using ClustalW and visualized through ESPript 3 (Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e). The highest sequence identity of 95.34% and 85.98% was observed between BbigLDH and \u003cem\u003eBabesia ovata and Babesia bovis, respectively.\u003c/em\u003e Other species closely related to BbigLDH, such as \u003cem\u003eTheileria orientalis\u003c/em\u003e exhibited a 70.66% match, while \u003cem\u003ePlasmodium falciparum\u003c/em\u003e showed a 50.31% of identity. Notably, there is a 30% identity with \u003cem\u003eBos taurus\u003c/em\u003e LDH A, LDH B, which is recognized as the host for \u003cem\u003eBabesia bigemina\u003c/em\u003e. The phylogenetic tree, constructed using the Neighbor Joining method, is provided in Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003e [\u003cspan class=\"CitationRef\"\u003e66\u003c/span\u003e]. The phylogenetic tree suggests that BbigLDH is positioned within the Babesia clade, with a closer evolutionary relationship to \u003cem\u003eB. ovata\u003c/em\u003e. In contrast, LDHs originating from mammalian hosts form a distinct clade, indicating a clear separation in their evolutionary lineage from those of the apicomplexan parasites. The InterPro database analysis of BbigLDH sequence (XP_012766822.1) revealed the presence of LDH1N-NAD binding (PF0056 / 9-152) and LDH1C-alpha/beta C terminal (PF02866 / 157\u0026ndash;315) domains. The 3D model of BbigLDH is expected to include a binding site for NADH or NAD\u0026thinsp;+\u0026thinsp;cofactors as well as a binding site for pyruvate or lactate substrates. The secondary structure of BbigLDH consists of sixteen alpha-helices, fourteen beta-strands, and various loops (Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e\n \u003cp\u003eBbigLDH displays a common five-amino acid insertion within its substrate specificity loop. Crystallographic analysis PfLDH have shown that this insertion contributes to a distinctive architecture in the active site, which could be a target for drug development [\u003cspan class=\"CitationRef\"\u003e18\u003c/span\u003e]. While the sequences in this region vary across species, for BbigLDH, the sequence is DDEWT, and it is predicted to form a coil structure (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eA). BbigLDH shares identical catalytic residues with those conserved in all LDHs, which include R109, D168, H195, and R171 (residues R98, D158, H185, R161 in BbigLDH) (Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e) [\u003cspan class=\"CitationRef\"\u003e5\u003c/span\u003e]. Residues 93\u0026ndash;97 (DDEWT) typically form a mobile loop in most apicomplexan LDHs (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eA). In the ternary complex, the loop moves downward to form the closed catalytic site [\u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e]. The loop opposite to the substrate specificity loop (G236, T237, G238 in BbigLDH) is called as \u0026ldquo;opposing loop\u0026rdquo; and depicted in Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB and Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003e [\u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e]. The distinctive two-residue deletion within the opposing loop of modern Plasmodium LDHs, along with the shared loop structure with ancient LDH (AncLDH) and Toxoplasma LDHs, marks a significant evolutionary deviation (Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003e) [\u003cspan class=\"CitationRef\"\u003e19\u003c/span\u003e]. Additionally, while ancestral LDH displays a slight oxaloacetate activity in common with modern Toxoplasma LDHs, this capability is absent in Plasmodium LDHs. This intriguing correlation strongly suggests that the opposing loop deletion, and potentially the presence of Ala236 and Pro246, play a pivotal role in shaping the remarkably stringent substrate specificity observed in modern Plasmodium LDHs [\u003cspan class=\"CitationRef\"\u003e67\u003c/span\u003e]. Unlike in modern Plasmodium LDHs, a deletion in the opposing loop is not observed in BbigLDH, like \u003cem\u003eToxoplasma gondii\u003c/em\u003e LDH. This sets it apart from the modern Plasmodium LDHs in this aspect. In \u003cem\u003eBabesia\u003c/em\u003e and \u003cem\u003eTheileria\u003c/em\u003e (excluding \u003cem\u003eT. equi\u003c/em\u003e) genus, a G[QSTN]G pattern is present (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e). Unlike in \u003cem\u003ePlasmodium\u003c/em\u003e and \u003cem\u003eB. taurus\u003c/em\u003e, such a pattern is not found, and there is a distinction in the opposing loop from both ancLDH and \u003cem\u003eToxoplasma\u003c/em\u003e LDH. While residues preceding the opposing loop in \u003cem\u003eB. taurus\u003c/em\u003e and \u003cem\u003ePlasmodium\u003c/em\u003e LDHs are positively charged (H and K residues), in \u003cem\u003eBabesia\u003c/em\u003e and \u003cem\u003eTheileria\u003c/em\u003e, these residues are observed to be aromatic (Y).\u003c/p\u003e\n \u003c/div\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e\n \u003ch2\u003eHomology modelling and model validation\u003c/h2\u003e\n \u003cp\u003eThe crystal structures of monomers in both closed and open states of BbigLDH were previously unknown. To address this, homology models were generated using the MODELLER 9.23 program. These models were developed for both open and closed states, incorporating pyruvate and NADH, respectively. The suitable template proteins were selected from PDB-structured proteins using BLASTp, based on the highest identity percentage. Additionally, during the selection process, emphasis was placed on ensuring that the template LDH structures belonged to apicomplexans and contained NADH and pyruvate in closed state. The template LDH proteins used for modelling included 1PZE (1.95 \u0026Aring;) [\u003cspan class=\"CitationRef\"\u003e68\u003c/span\u003e] - the \u003cem\u003eT. gondii\u003c/em\u003e LDH1 apoenzyme structure- for the open state, and 4PLC (1.50 \u0026Aring;) - an ancestral LDH from apicomplexans [\u003cspan class=\"CitationRef\"\u003e67\u003c/span\u003e]- for the closed state, along with pyruvate and NADH. The generated 100 models were selected based on the lowest DOPE value and subsequently validated using various methods. One of the validation programs employed was ERRAT, which examines the consistency of atomic interactions within the modelled structure. The overall quality factors obtained from the ERRAT server for the open and closed models were 90.1274 and 93.949, respectively [\u003cspan class=\"CitationRef\"\u003e69\u003c/span\u003e]. These values indicate a relatively high quality of the modelled structures. Another validation program utilized was ProSA [\u003cspan class=\"CitationRef\"\u003e70\u003c/span\u003e], which evaluates the quality of generated models by constructing an energy profile of structural residue areas. The Z-scores obtained from ProSA were compared with those of experimentally validated structures to assess consistency. According to ProSA results, the Z-scores for mutant and wild-type models were \u0026minus;\u0026thinsp;3.3 and \u0026minus;\u0026thinsp;3.19, respectively. Models falling within the acceptable range of residue and energy efficiency were considered reliable. Quantitative Model Energy Analysis (QMEAN) is a valuable tool for analyzing crucial geometric aspects of protein structures, such as determining natural curvatures and calculating torsion angle potentials of amino acids [\u003cspan class=\"CitationRef\"\u003e71\u003c/span\u003e]. The QMEAN value is originally in a range [0,1] with one being good expected from high resolution X-ray structures. In the results, QMEAN4 values were \u0026minus;\u0026thinsp;2.47 and \u0026minus;\u0026thinsp;2.19 for open state and closed state, respectively. Verify3D is a validation program that evaluates the absolute secondary and tertiary folding behaviour of a model protein by analyzing the 3D-1D scores in the context of 3D structure-sequence relationships. A higher value for amino acids is expected compared to the threshold value, and the overall simplified score determines the structure\u0026apos;s validity [\u003cspan class=\"CitationRef\"\u003e72\u003c/span\u003e]. Homology models passed by 88.82% and 93.79% for open state and close state, respectively according to the assessment based on favourable 3D-1D scores. A Ramachandran plot graph was generated for the model proteins using ProCHECK [\u003cspan class=\"CitationRef\"\u003e73\u003c/span\u003e]. The percentages of residues in the favoured region were 88.1% and 86.6% respectively for the open and closed states. Based on the assessment of validation values, all these validation results confirmed that the enzymes were modelled with high accuracy and reliability (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eA). Additionally, the structural difference introduced by the 5-AA insertion in the open form of BbigLDH compared to host LDH A and B is shown in Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB. Detailed results obtained from validation were provided in Supplement file, Table \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e. \u0026nbsp;\u003c/p\u003e\n \u003ctable id=\"Tab2\" border=\"1\"\u003e\n \u003ccaption language=\"En\"\u003e\n \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\n \u003cdiv class=\"CaptionContent\"\u003e\n \u003cp\u003eTwenty compounds, exhibiting lower affinity with BtLDH A, B and A-like B and their average values (kcal/mol)\u003c/p\u003e\n \u003c/div\u003e\n \u003c/caption\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eCompound no\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eCompound ID\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eBtLDHA\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eBtLDHB\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eBtLDHA-likeB\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eAverage\u003c/p\u003e\n \u003c/th\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC1\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1478445725\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.311\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.156\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.875\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.78067\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC2\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1246322905\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.495\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.448\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.34\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.42767\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC3\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ826473056\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.253\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.534\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.414\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.40033\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC4\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ509196536\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.233\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.03\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.275\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.51267\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC5\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1424270493\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.559\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.411\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.247\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.739\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC6\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ278525320\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.464\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.047\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.283\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.26467\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC7\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1640818270\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.863\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.026\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.731\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.87333\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC8\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ604033838\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.745\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.219\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.344\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.76933\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC9\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2373310039\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.842\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.14\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.485\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.82233\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC10\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ973658262\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.959\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.858\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.883\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.23333\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC11\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1616793138\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.507\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.505\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.672\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.56133\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC12\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2953881032\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.596\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.577\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.391\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.188\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC13\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1253206711\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.843\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.798\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.772\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.80433\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC14\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ19650139\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.249\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.068\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.262\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.193\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC15\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2400537901\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.409\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.512\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.829\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.91667\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC16\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2201943136\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.406\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.435\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.88033\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC17\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2060912694\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.114\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.636\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.857\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.53567\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC18\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2063989764\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-3.427\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.208\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.287\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.64067\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC19\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1531526443\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.526\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-5.325\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.901\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.91733\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC20\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2055833084\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.401\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.787\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.808\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-4.66533\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n \u003c/table\u003e\n \u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n \u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec11\" class=\"Section2\"\u003e\n \u003ch2\u003eValidation of Molecular Docking Procedure\u003c/h2\u003e\n \u003cp\u003eThe validation of the molecular docking procedure to be used in computational biology methods is an essential preliminary step to ensure the validity and reliability of the method. In the validation process, inhibitors present in the crystal structures of PfLDH with known structures were redocked into the original protein structure using the cognate docking procedure in methods section. The superposition of the resulting protein structure with the original protein structure is examined, along with the binding mode of the ligand and the RMSD values. Low RMSD values (The superimposed values for 1U4O, 1T25 and 1XIV are 0.14 \u0026Aring;, 0.21 \u0026Aring;, and 0.32 \u0026Aring;, respectively.) indicate a high degree of similarity between the two structures. Other RMSD values for 1U4S and 1T26 are 1.26 \u0026Aring; and 1.06 \u0026Aring;, respectively which have more higher superimposed values then other PDB structures. Additionally, when evaluating the cognate docking results, a threshold value of -7.00 kcal/mol for the docking score has been determined as the preferred threshold for assessing the HTVS docking results.\u003c/p\u003e\n \u003cp\u003eThe enrichment calculation results provide valuable information about the performance of a computational model or algorithm in virtual screening or drug discovery experiments. The model was tested on a dataset consisting of 1000 ligands, including 34 known actives. The goal was to determine the model\u0026apos;s ability to accurately rank the active compounds and enrich them at the top of the ranked list. The obtained BEDROC values ranged from 0.229 to 0.332, suggesting a moderate enrichment capability. The receiver operating characteristic (ROC) curve and its associated area under the curve (AUC) were also analyzed to evaluate the overall classification accuracy of the model. The ROC value of 0.69 suggests a reasonable discriminative ability, indicating that the model can distinguish between active and inactive compounds to a certain extent. The area under the accumulation curve (AUC) is 0.68, indicating that a significant portion of the active compounds is discovered early in the screening process. The obtained RIE value of 3.43 indicates that the model is performing better than random selection but still leaves room for improvement. It highlights the model\u0026apos;s ability to detect structurally diverse active compounds (Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003e and Table \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec12\" class=\"Section2\"\u003e\n \u003ch2\u003eStructure-Based Virtual Screening, QM-Polarized Ligand Docking and Induced Fit Docking (IFD)\u003c/h2\u003e\n \u003cp\u003eFour structures of BbigLDH were utilized, comprising two open and two closed states. A pair of open and closed states were generated after a 100 ns MD simulation using NAMD. For the prediction of binding sites, DoGSiteScorer method was employed on all models. [\u003cspan class=\"CitationRef\"\u003e74\u003c/span\u003e]. In both the open and closed states of BbigLDH, druggable pockets were predicted, with a druggability score of 0.81 and simple score of 0.6. The pocket includes the pyruvate and NADH binding regions. The virtual library was initially docked to the active site of BbigLDH using GLIDE HTVS, SP and XP docking methods, respectively. The top 100 poses were then ranked based on XP GScores (Table \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e). After comparison of compounds interacting with common residues, 256 ligands were selected. Subsequently, the remaining ligands underwent ADME analysis using the QikProp program. Compounds within the specified reference value range from the manual were identified (133 compounds). After ADME analysis, the Quantum Mechanics Polarized Ligand Docking method was applied, where the QM charges were calculated by using the QSite (Jaguar) and each ligand was redocked from its corresponding pose set. The poses with the lowest Glide Emodel values were chosen. Results with XP GScore values equal to or higher than \u0026minus;\u0026thinsp;7.00 kcal/mol were ranked to obtain the best outcomes for 70 models. Following this, induced fit docking was employed This involved relaxing the ligand into a conformationally relaxed receptor state, followed by redocking using Glide. The resulting 70 poses were ranked based on their IFD score, with lower scores being prioritized. Remaining compounds were also docked to the host LDH (Bos taurus LDH A, LDH B and LDH AL6B) to identify selective inhibitor candidates. The structural difference introduced by the 5-AA insertion in the open form of BbigLDH compared to host LDH A and B is shown in Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB. Twenty compounds, exhibiting lower affinity with host LDHs and available for purchase were selected for molecular dynamics simulations and end-point binding free energy calculations (Table \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e\n \u003cp\u003eThe evaluation of potential drug candidates involves a thorough assessment of various criteria, with a close examination of their physicochemical properties intricately linked to their structural attributes. Minor modifications in their structures have resulted in enhanced properties. In our study, the hit compounds underwent analysis to determine their compliance with Lipinski\u0026apos;s rule of five. [\u003cspan class=\"CitationRef\"\u003e55\u003c/span\u003e]. According to Lipinski\u0026apos;s rule of five, an established set of guidelines for drug likeness, the molecular weight of potential drug candidates should be less than 500 Da. Additionally, the hydrogen bond donor count should be fewer than 5, the hydrogen bond acceptor count should be less than 10, and the predicted octanol/water partition coefficient should ideally fall within the range of -2.0 to 4. These criteria help assess the likelihood of a compound\u0026apos;s success as a drug candidate based on its physicochemical properties [\u003cspan class=\"CitationRef\"\u003e75\u003c/span\u003e]. All values for the compounds fall within the specified range for Lipinski rules and PlogP\u003csub\u003e0/w\u003c/sub\u003e. The Plog values for both hERG channel blocking and serum albumin binding fell within the reference ranges [\u003cspan class=\"CitationRef\"\u003e76\u003c/span\u003e]. The predicted IC50 value for blocking HERG K\u0026thinsp;+\u0026thinsp;channels should fall within the accepted range, with values greater than \u0026minus;\u0026thinsp;5.0. In this study, the compound with the largest value is C18 (5.073), while the one with the smallest value is compound C9 (-2.504). On the other hand, the predicted apparent Caco-2 cell permeability should be assessed, with values indicating low permeability if less than 25 and high permeability if greater than 500 [\u003cspan class=\"CitationRef\"\u003e77\u003c/span\u003e]. According to this information, it is observed that the compound with the highest permeability is C20 (2104.635), while the one with the least permeable is C19 (74.968). The predicted blood-brain barrier permeability coefficient should be within the range of -3.0 to 1.2 [\u003cspan class=\"CitationRef\"\u003e78\u003c/span\u003e]. The predicted binding capacity to human serum albumin is considered appropriate if it falls between \u0026minus;\u0026thinsp;1.5 and 1.5 [\u003cspan class=\"CitationRef\"\u003e79\u003c/span\u003e]. All values for the compounds fall within the specified range for PlogBB, PlogKhsa. The predicted human oral absorption, rated on a scale of 0 to 100, is considered favorable if it exceeds 70%. The compound with the highest human oral absorption rates is C18 and C20 (both 100), whereas C19 has a rate of 72.307 (Supplement file, Table 3).\u0026nbsp;\u003c/p\u003e\n \u003ctable id=\"Tab4\" border=\"1\"\u003e\n \u003ccaption language=\"En\"\u003e\n \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e\n \u003cdiv class=\"CaptionContent\"\u003e\n \u003cp\u003eList of compounds which have favourable XP GScores. Free binding energy (∆G) values were listed for the last 5 ns of 20 ns MD simulation.\u003c/p\u003e\n \u003c/div\u003e\n \u003c/caption\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eCompound no\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eCompound ID\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eStructure of BbigLDH\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eHTVS; XP gscore\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eQPLD; XP gscore\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eIFD; XP gscore\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003ePrime energy\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003eIFD score\u003c/p\u003e\n \u003c/th\u003e\n \u003cth align=\"left\"\u003e\n \u003cp\u003edG Average\u003c/p\u003e\n \u003c/th\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC1\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1478445725\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.343\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.515\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.569\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13590.36\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-689.09\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-67.3425\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC2\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1246322905\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.608\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.140\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.048\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13502.23\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-684.16\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-58.4571\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC3\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ826473056\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.236\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.530\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-11.017\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13460.85\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-684.06\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-64.0088\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC4\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ509196536\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.997\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.389\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-10.565\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13452.59\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-683.19\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-57.2994\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC5\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1424270493\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.869\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.598\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.404\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13471.78\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-682.99\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-60.7591\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC6\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ278525320\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.497\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.480\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-10.557\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13439.77\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-682.55\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-54.8145\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC7\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1640818270\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.479\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.473\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.652\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13462.95\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-680.8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-64.4233\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC8\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ604033838\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.382\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.250\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-10.285\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13409.02\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-680.74\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-66.175\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC9\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2373310039\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.173\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.958\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.998\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13442.56\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-680.13\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-75.825\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC10\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ973658262\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.176\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.182\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.803\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13401.97\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-679.9\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-46.3542\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC11\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1616793138\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.648\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.555\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.624\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13403.46\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-679.8\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-46.0935\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC12\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2953881032\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state after MDS\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.202\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.972\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.058\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13428.81\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-679.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-71.5432\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC13\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1253206711\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.925\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.731\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-10.617\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13374.26\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-679.33\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-53.7573\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC14\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ19650139\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.701\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.111\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.606\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13429.61\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-679.09\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-60.4939\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC15\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2400537901\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eClosed state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-10.464\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.919\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.831\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13379.26\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-678.79\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-67.0062\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC16\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2201943136\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eOpen state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.992\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.443\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.72\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13399.65\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-678.70\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-60.9869\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC17\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2060912694\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eOpen state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.774\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.274\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-10.44\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13364.38\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-678.66\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-60.6208\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC18\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2063989764\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eOpen state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.105\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.750\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.56\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13377.59\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-678.43\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-72.3022\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC19\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ1531526443\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eOpen state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.641\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.706\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.55\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13338.36\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-676.47\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-72.2696\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003e\u003cstrong\u003eC20\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eZ2055833084\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"left\"\u003e\n \u003cp\u003eOpen state\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-8.482\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-7.328\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-9.19\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-13311.35\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-674.76\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd align=\"char\"\u003e\n \u003cp\u003e-44.1338\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n \u003c/table\u003e\n \u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec13\" class=\"Section2\"\u003e\n \u003ch2\u003eCompound prioritization based on molecular dynamics simulations and binding energy\u003c/h2\u003e\n \u003cp\u003eIn this study, we conducted MD simulations and binding free energy calculations of 20 protein-ligand complexes obtained from the docking campaigns using the Desmond program for 20 ns. Furthermore, to determine ligand stability over a long period of time, MDS with GROMACS was performed for 100 ns, too.\u003c/p\u003e\n \u003cp\u003eBased on the RMSD values from the GROMACS simulations and the ∆G average values from the 20 ns MD simulations using Desmond (Table \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003e), we identified C7, C9, C12, C14, C15, C16, C18 and C20 compounds as candidates for an extended MD simulation for 300 ns. Subsequently, we assessed the RMSD values obtained from a 300 ns MD simulation (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003e). According to the Lig-fit-prot RMSD values in Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eB, all compounds except C20 were found in the binding cavity at the end of the simulation. Moreover, C16 demonstrated remarkable positional stability, maintaining a 3 \u0026Aring; RMSD value throughout the simulation. According to the MM/GBSA ∆G Distribution graph in Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eA, both C9 and subsequently C16 displayed favourable binding energies at the last 30 ns of the 300 ns simulation.\u003c/p\u003e\n \u003cp\u003e\u003cbr\u003e\u003c/p\u003e\n \u003cp\u003eIn the subsequent phase, we selected C9, C15, C16, and C18 for an extensive 1000 ns MDS analysis according to the ∆G and RMSD values, and interaction with the active site residues. The RMSD analysis graph in Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eB reveals that C16 maintained stability within the range of 2\u0026ndash;4 \u0026Aring;. In the case of C9, stability was observed after 400 ns, with a stability range of 4\u0026ndash;6 \u0026Aring;. C18 exhibited stability after 150 ns, within a range of 6\u0026ndash;8 \u0026Aring;. C15 moving away from the binding pocket, therefore RMSD analysis was not performed. In Fig. \u003cspan class=\"InternalRef\"\u003e5\u003c/span\u003eC, the RMSF (Root Mean Square Fluctuation) graph indicates that residues between 80 and 100, which include the substrate-specific loop, exhibit dramatic changes for C18 compared to other compounds. Additionally, the residues between 235 and 240, which encompass the opposing loop, show changes for C9 and C16, distinct from those observed for C18.\u003c/p\u003e\n \u003cp\u003eAccording to the MM/GBSA ∆G Distribution graph in Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB, both C9, C16 and C18 displayed favourable binding energies in the last 100 ns of the 1000 ns simulation. The ∆G averages are \u0026minus;\u0026thinsp;52.5 kcal/mol, -62.1 kcal/mol and \u0026minus;\u0026thinsp;62.5 kcal/mol for C9, C16 and C18, respectively.\u003c/p\u003e\n \u003cp\u003eFigure \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003e illustrates the average energies contributing to the overall MM-GBSA binding energy. Calculations were based on snapshots taken from the last 100 ns of a 1 \u0026micro;s MD simulation. The MM-GBSA \u0026Delta;G binding energies are quite similar for C16, and C18 (\u0026minus;\u0026thinsp;62.11 kcal/mol, and \u0026minus;\u0026thinsp;62.52 kcal/mol, respectively), but the main driving forces behind these binding energies different. For C9-protein binding, lipophilic energy and van der Waals (vdW) energy are the predominant types of energy. In contrast, the contributions of coulomb and hydrogen bond (H-bond) energies in the binding of C16 and C18 are markedly higher from those in C9-protein binding.\u003c/p\u003e\n \u003cp\u003eIn Fig. \u003cspan class=\"InternalRef\"\u003e8\u003c/span\u003e, all three compounds are found in the expected regions after the 1000 ns simulation time. All compounds are positioned within the enzyme\u0026apos;s catalytic pocket, which is of crucial importance, and interactions with catalytic residues are evident. C9 retained its significant H-bond interactions after IFD, but lost them after MD. In Supplementary Fig. 6A and 7A, according to BbigLDH-C9 interactions, there are water bridges which interaction mediated between a water molecule and residues including T82, L85, D93, I128, and T129. On the other hand, upon examination of Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB and \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eB, it is evident that in the BbigLDH-C16 complex, there are hydrogen bond interactions with residues N17 (77%), L85 (97%), R87 (94%), and N130 (70%), along with water Bridges involving L231 and pi-pi stacking with Y241 (32%). In Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eC and \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eC, hydrogen bond interactions are observed with residues R87, I128, N130 (46%) and R161 (39%) in the BbigLDH-C18 complex, along with pi-cation interactions involving R161 (33%). In Fig. 9A, after the induced fit docking stage, interactions with R87, I128, and N130 residues were generally observed with all compounds. On the other hand in Fig. 9B, after the 1000 ns simulation times, interactions with R87, N130 conserved. It can be observed that the interaction with these residues is maintained throughout the MD simulation for both C16 and C18 (Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003e). It is believed that these residues play a significant role in inhibition against BbigLDH. In the cognate docking stage, inhibitors found in the open and closed PfLDH structures, whose PDB structures are known, were docked to the open and closed states of BbigLDH structures. The interaction of the original PDB structures with the inhibitor and the interaction of the BbigLDH structures with the same inhibitor are shown in Supplement file, Fig. \u003cspan class=\"InternalRef\"\u003e8\u003c/span\u003e. It is observed that the lack of a stable region in the substrate-specific loop region of the open state of PfLDH structures results in differences with the open state of BbigLDH structure. Additionally, it is seen that the residues interacting with the ligands in the open state of BbigLDH structure and also known to have catalytic properties in some cases are R87, I128, N130, R161, and S239. On the other hand, for the closed state of BbigLDH, the most important residues interacting with the ligands are R98, N130, R161, H185, S239, and H240, as shown in Supplement file, Fig.\u0026nbsp;9. In conclusion, it is considered that compounds showing similarity in the critical residues obtained from the redocking of PfLDH inhibitors, whose PDB structure is known, to BbigLDH, and interacting with residues C09, C16, and C18, may have similar inhibitory properties on BbigLDH.\u003c/p\u003e\n \u003cp\u003eThese findings strongly indicate that compounds located in the pocket where the cofactor and substrate bind may inhibit enzyme activity through these specific residue interactions. Particularly, C16 may be a promising candidate for in vitro inhibitor studies, given its stable interaction with BbigLDH. Compounds Z2373310039 (C09), Z2201943136 (C16) and Z2063989764 (C18) may warrant further investigation to comprehend their potential therapeutic effects on babesiosis.\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Conclusions","content":"\u003cp\u003eRecent reports have indicated the emergence of resistance against certain existing anti-babesial drugs. This highlights the urgency of identifying targets for the development of new drugs with distinct modes of action. The lactate dehydrogenase enzyme plays a pivotal role in the anaerobic pathway of the parasite, making it a promising candidate for further drug design endeavours. We employed homology modelling to predict the unknown structure of BbigLDH, followed by structure-based virtual screening. Docking scores, interaction with active site residues, affinity to host LDHs, ADME properties, and stability through simulations were considered during identifying of potential inhibitors. The stability of the three candidate compounds (C9, C16, and C18) was validated through long molecular dynamics simulations and further binding free energy calculations. The obtained results provided insights into the molecular foundation of interactions and, consequently, the mechanism of inhibition. Taken together, the comprehensive analysis of these findings underscores the utility of structure-based virtual screening in providing crucial insights into the essential structural and binding features necessary for the design of novel BbigLDH inhibitors, offering potential avenues for the treatment of bovine babesiosis.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthical Approval\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work has been supported by Marmara University-Scientific Research Projects Coordination Unit under grant number:\u0026nbsp;FDK-2022-10061.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSafiye Merve Bostancioglu:\u003c/strong\u003e Investigation; Visualization; Writing-Original Draft.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eOzal Mutlu:\u003c/strong\u003e Conceptualization; Methodology; Investigation; Visualization; Writing-Original Draft.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDeclaration of competing interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDisclosure statement\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNo potential conflict of interest was reported by the authors.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNo datasets were generated or analysed during the current study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eORCID\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSafiye Merve Bostancioglu\u003c/strong\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003ehttps://orcid.org/0000-0002-0239-3966\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eOzal Mutlu\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ehttps://orcid.org/0000-0003-4551-5780\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eT. Smith and F. L. Kilborne, \u003cem\u003eInvestigations into the nature, causation, and prevention of Texas or southern cattle fever /\u003c/em\u003e. Washington, D.C. : U.S. Dept. of Agriculture, Bureau of Animal Industry, 1893. doi: 10.5962/bhl.title.124068.\u003c/li\u003e\n\u003cli\u003eL. Schnittger, A. E. Rodriguez, M. Florin-Christensen, and D. A. Morrison, \u0026lsquo;Babesia: A world emerging\u0026rsquo;, \u003cem\u003eInfect. Genet. Evol.\u003c/em\u003e, vol. 12, no. 8, pp. 1788\u0026ndash;1809, Dec. 2012, doi: 10.1016/j.meegid.2012.07.004.\u003c/li\u003e\n\u003cli\u003eR. Bock, L. Jackson, A. De Vos, and W. Jorgensen, \u0026lsquo;Babesiosis of cattle\u0026rsquo;, \u003cem\u003eParasitology\u003c/em\u003e, vol. 129, no. S1, pp. S247\u0026ndash;S269, Oct. 2004, doi: 10.1017/S0031182004005190.\u003c/li\u003e\n\u003cli\u003eM. Calvopi\u0026ntilde;a \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Case report: First report on human infection by tick-borne Babesia bigemina in the Amazon region of Ecuador\u0026rsquo;, \u003cem\u003eFront. Public Health\u003c/em\u003e, vol. 11, p. 1079042, Aug. 2023, doi: 10.3389/fpubh.2023.1079042.\u003c/li\u003e\n\u003cli\u003eL. He, R. G. Bastos, L. Yu, J. M. Laughery, and C. E. Suarez, \u0026lsquo;Lactate Dehydrogenase as a Potential Therapeutic Drug Target to Control Babesia bigemina\u0026rsquo;, \u003cem\u003eFront. Cell. Infect. Microbiol.\u003c/em\u003e, vol. 12, p. 870852, Apr. 2022, doi: 10.3389/fcimb.2022.870852.\u003c/li\u003e\n\u003cli\u003eJ. Mosqueda, A. Olvera-Ramirez, G. Aguilar-Tipacamu, and G. J. Canto, \u0026lsquo;Current Advances in Detection and Treatment of Babesiosis\u0026rsquo;, \u003cem\u003eCurr. Med. Chem.\u003c/em\u003e, vol. 19, no. 10, pp. 1504\u0026ndash;1518, Apr. 2012, doi: 10.2174/092986712799828355.\u003c/li\u003e\n\u003cli\u003eJ. H. D. Silva \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Chemoprophylaxis for babesiosis and anaplasmosis in cattle: case report\u0026rsquo;, \u003cem\u003eRev. Bras. Parasitol. Veterin\u0026aacute;ria\u003c/em\u003e, vol. 29, no. 4, p. e010520, 2020, doi: 10.1590/s1984-29612020096.\u003c/li\u003e\n\u003cli\u003eN. D. Shah, A. U. Bhikane, R. K. Jadhav, S. G. Chavhan, and A. Mohan, \u0026lsquo;Therapeutic management of Babesiosis alone and its mixed infection with Theileriosis in Cattle\u0026rsquo;, 2019.\u003c/li\u003e\n\u003cli\u003eB. Tuvshintulga, T. Sivakumar, N. Yokoyama, and I. Igarashi, \u0026lsquo;Development of unstable resistance to diminazene aceturate in Babesia bovis\u0026rsquo;, \u003cem\u003eInt. J. Parasitol. Drugs Drug Resist.\u003c/em\u003e, vol. 9, pp. 87\u0026ndash;92, Apr. 2019, doi: 10.1016/j.ijpddr.2019.02.001.\u003c/li\u003e\n\u003cli\u003eTodorovic R. A., \u0026lsquo;Chemourophylaxis (Imidocarb. Against Babesia bigemina and Babesia argentina Infections\u0026rsquo;, \u003cem\u003eAm J Vet Res\u003c/em\u003e, vol. 34(9):1153-61, Sep. 1973, doi: PMID: 4747036.\u003c/li\u003e\n\u003cli\u003eN. G. Coldham \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Imidocarb residues in edible bovine tissues and in vitro assessment of imidocarb metabolism and cytotoxicity\u0026rsquo;, \u003cem\u003eDrug Metab. Dispos. Biol. Fate Chem.\u003c/em\u003e, vol. 23, no. 4, pp. 501\u0026ndash;505, Apr. 1995.\u003c/li\u003e\n\u003cli\u003eH. J. Vial and A. Gorenflot, \u0026lsquo;Chemotherapy against babesiosis\u0026rsquo;, \u003cem\u003eVet. Parasitol.\u003c/em\u003e, vol. 138, no. 1\u0026ndash;2, pp. 147\u0026ndash;160, May 2006, doi: 10.1016/j.vetpar.2006.01.048.\u003c/li\u003e\n\u003cli\u003eA. Zintl, G. Mulcahy, H. E. Skerrett, S. M. Taylor, and J. S. Gray, \u0026lsquo;\u003cem\u003eBabesia divergens\u003c/em\u003e , a Bovine Blood Parasite of Veterinary and Zoonotic Importance\u0026rsquo;, \u003cem\u003eClin. Microbiol. Rev.\u003c/em\u003e, vol. 16, no. 4, pp. 622\u0026ndash;636, Oct. 2003, doi: 10.1128/CMR.16.4.622-636.2003.\u003c/li\u003e\n\u003cli\u003eC. E. Suarez and S. Noh, \u0026lsquo;Emerging perspectives in the research of bovine babesiosis and anaplasmosis\u0026rsquo;, \u003cem\u003eVet. Parasitol.\u003c/em\u003e, vol. 180, no. 1\u0026ndash;2, pp. 109\u0026ndash;125, Aug. 2011, doi: 10.1016/j.vetpar.2011.05.032.\u003c/li\u003e\n\u003cli\u003eN. McHARDY, R. M. Woollon, R. B. Clampitt, J. A. James, and R. J. Crawley, \u0026lsquo;Efficacy, toxicity and metabolism of imidocarb dipropionate in the treatment of Babesia ovis infection in sheep\u0026rsquo;, \u003cem\u003eRes. Vet. Sci.\u003c/em\u003e, vol. 41, no. 1, pp. 14\u0026ndash;20, Jul. 1986, doi: 10.1016/S0034-5288(18)30565-4.\u003c/li\u003e\n\u003cli\u003eC. A. Homewood and K. D. Neame, \u0026lsquo;Conversion of glucose to lactate by intraerythrocytic \u003cem\u003ePlasmodium berghei\u003c/em\u003e\u0026rsquo;, \u003cem\u003eAnn. Trop. Med. Parasitol.\u003c/em\u003e, vol. 77, no. 2, pp. 127\u0026ndash;129, Apr. 1983, doi: 10.1080/00034983.1983.11811686.\u003c/li\u003e\n\u003cli\u003eK. L. Kavanagh, R. A. Elling, and D. K. Wilson, \u0026lsquo;Structure of \u003cem\u003eToxoplasma gondii\u003c/em\u003e LDH1: Active-Site Differences from Human Lactate Dehydrogenases and the Structural Basis for Efficient APAD \u003csup\u003e+\u003c/sup\u003e Use \u003csup\u003e,\u003c/sup\u003e\u0026rsquo; \u003cem\u003eBiochemistry\u003c/em\u003e, vol. 43, no. 4, pp. 879\u0026ndash;889, Feb. 2004, doi: 10.1021/bi035108g.\u003c/li\u003e\n\u003cli\u003eC. R. Dunn \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;The structure of lactate dehydrogenase from Plasmodium falciparum reveals a new target for anti-malarial design\u0026rsquo;, \u003cem\u003eNat. Struct. Biol.\u003c/em\u003e, vol. 3, no. 11, pp. 912\u0026ndash;915, Nov. 1996, doi: 10.1038/nsb1196-912.\u003c/li\u003e\n\u003cli\u003eJ. I. Boucher, J. R. Jacobowitz, B. C. Beckett, S. Classen, and D. L. Theobald, \u0026lsquo;An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases\u0026rsquo;, \u003cem\u003eeLife\u003c/em\u003e, vol. 3, p. e02304, Jun. 2014, doi: 10.7554/eLife.02304.\u003c/li\u003e\n\u003cli\u003eS. Bork, M. Okamura, S. Boonchit, H. Hirata, N. Yokoyama, and I. Igarashi, \u0026lsquo;Identification of Babesia bovis l-lactate dehydrogenase as a potential chemotherapeutical target against bovine babesiosis\u0026rsquo;, \u003cem\u003eMol. Biochem. Parasitol.\u003c/em\u003e, vol. 136, no. 2, pp. 165\u0026ndash;172, Aug. 2004, doi: 10.1016/j.molbiopara.2004.03.009.\u003c/li\u003e\n\u003cli\u003eD. Turgut-Balik \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Cloning, sequence and expression of the lactate dehydrogenase gene from the human malaria parasite, Plasmodium vivax\u0026rsquo;, \u003cem\u003eBiotechnol. Lett.\u003c/em\u003e, vol. 26, no. 13, pp. 1051\u0026ndash;1055, Jul. 2004, doi: 10.1023/B:BILE.0000032958.78158.10.\u003c/li\u003e\n\u003cli\u003eO. Mutlu and D. Turgut Balik, \u0026lsquo;Kinetic Analysis of the Amino Terminal End of Active Site Loop of Lactate Deyhdrogenase from Plasmodium Vivax\u0026rsquo;, \u003cem\u003eBalk. Med. J.\u003c/em\u003e, vol. 29, no. 4, pp. 364\u0026ndash;369, Dec. 2012, doi: 10.5152/balkanmedj.2012.047.\u003c/li\u003e\n\u003cli\u003eA. Erdemir and O. Mutlu, \u0026lsquo;Functional and structural characterization of the pentapeptide insertion of Theileria annulata lactate dehydrogenase by site-directed mutagenesis, comparative modeling and molecular dynamics simulations\u0026rsquo;, \u003cem\u003eJ. Mol. Graph. Model.\u003c/em\u003e, vol. 74, pp. 344\u0026ndash;351, Jun. 2017, doi: 10.1016/j.jmgm.2017.04.010.\u003c/li\u003e\n\u003cli\u003eL. Yu, Q. Liu, W. Luo, J. Zhao, H. F. Alzan, and L. He, \u0026lsquo;The Structural Basis of Babesia orientalis Lactate Dehydrogenase\u0026rsquo;, \u003cem\u003eFront. Cell. Infect. Microbiol.\u003c/em\u003e, vol. 11, p. 790101, Jan. 2022, doi: 10.3389/fcimb.2021.790101.\u003c/li\u003e\n\u003cli\u003eM. S. Gomez \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Substrate and cofactor specificity and selective inhibition of lactate dehydrogenase from the malarial parasite P. falciparum\u0026rsquo;, \u003cem\u003eMol. Biochem. Parasitol.\u003c/em\u003e, vol. 90, no. 1, pp. 235\u0026ndash;246, Dec. 1997, doi: 10.1016/S0166-6851(97)00140-0.\u003c/li\u003e\n\u003cli\u003eF. Al-Anouti, S. Tomavo, S. Parmley, and S. Ananvoranich, \u0026lsquo;The Expression of Lactate Dehydrogenase Is Important for the Cell Cycle of Toxoplasma gondii\u0026rsquo;, \u003cem\u003eJ. Biol. Chem.\u003c/em\u003e, vol. 279, no. 50, pp. 52300\u0026ndash;52311, Dec. 2004, doi: 10.1074/jbc.M409175200.\u003c/li\u003e\n\u003cli\u003eD. Madern, \u0026lsquo;Evolution of Cryptosporidium parvum Lactate Dehydrogenase from Malate Dehydrogenase by a Very Recent Event of Gene Duplication\u0026rsquo;, \u003cem\u003eMol. Biol. Evol.\u003c/em\u003e, vol. 21, no. 3, pp. 489\u0026ndash;497, Dec. 2003, doi: 10.1093/molbev/msh042.\u003c/li\u003e\n\u003cli\u003eL. Yu \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Crystal structures of \u003cem\u003eBabesia microti\u003c/em\u003e lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis\u0026rsquo;, \u003cem\u003eFASEB J.\u003c/em\u003e, vol. 33, no. 12, pp. 13669\u0026ndash;13682, Dec. 2019, doi: 10.1096/fj.201901259R.\u003c/li\u003e\n\u003cli\u003eS. Akash \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Anti-parasitic drug discovery against Babesia microti by natural compounds: an extensive computational drug design approach\u0026rsquo;, \u003cem\u003eFront. Cell. Infect. Microbiol.\u003c/em\u003e, vol. 13, p. 1222913, Aug. 2023, doi: 10.3389/fcimb.2023.1222913.\u003c/li\u003e\n\u003cli\u003eE. Lionta, G. Spyrou, D. Vassilatis, and Z. Cournia, \u0026lsquo;Structure-Based Virtual Screening for Drug Discovery: Principles, Applications and Recent Advances\u0026rsquo;, \u003cem\u003eCurr. Top. Med. Chem.\u003c/em\u003e, vol. 14, no. 16, pp. 1923\u0026ndash;1938, Oct. 2014, doi: 10.2174/1568026614666140929124445.\u003c/li\u003e\n\u003cli\u003eK. Tamura, G. Stecher, and S. Kumar, \u0026lsquo;MEGA11: Molecular Evolutionary Genetics Analysis Version 11\u0026rsquo;, \u003cem\u003eMol. Biol. Evol.\u003c/em\u003e, vol. 38, no. 7, pp. 3022\u0026ndash;3027, Jun. 2021, doi: 10.1093/molbev/msab120.\u003c/li\u003e\n\u003cli\u003eB. Webb and A. Sali, \u0026lsquo;Comparative Protein Structure Modeling Using MODELLER\u0026rsquo;, \u003cem\u003eCurr. Protoc. Bioinforma.\u003c/em\u003e, vol. 54, no. 1, Jun. 2016, doi: 10.1002/cpbi.3.\u003c/li\u003e\n\u003cli\u003eE. F. Pettersen \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;UCSF Chimera\u0026mdash;A visualization system for exploratory research and analysis\u0026rsquo;, \u003cem\u003eJ. Comput. Chem.\u003c/em\u003e, vol. 25, no. 13, pp. 1605\u0026ndash;1612, Oct. 2004, doi: 10.1002/jcc.20084.\u003c/li\u003e\n\u003cli\u003eC. Colovos and T. O. Yeates, \u0026lsquo;Verification of protein structures: Patterns of nonbonded atomic interactions\u0026rsquo;, \u003cem\u003eProtein Sci.\u003c/em\u003e, vol. 2, no. 9, pp. 1511\u0026ndash;1519, Sep. 1993, doi: 10.1002/pro.5560020916.\u003c/li\u003e\n\u003cli\u003eM. J. Sippl, \u0026lsquo;Recognition of errors in three-dimensional structures of proteins\u0026rsquo;, \u003cem\u003eProteins Struct. Funct. Genet.\u003c/em\u003e, vol. 17, no. 4, pp. 355\u0026ndash;362, Dec. 1993, doi: 10.1002/prot.340170404.\u003c/li\u003e\n\u003cli\u003eM. Wiederstein and M. J. Sippl, \u0026lsquo;ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins\u0026rsquo;, \u003cem\u003eNucleic Acids Res.\u003c/em\u003e, vol. 35, no. Web Server, pp. W407\u0026ndash;W410, May 2007, doi: 10.1093/nar/gkm290.\u003c/li\u003e\n\u003cli\u003eR. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, \u0026lsquo;PROCHECK: a program to check the stereochemical quality of protein structures\u0026rsquo;, \u003cem\u003eJ. Appl. Crystallogr.\u003c/em\u003e, vol. 26, no. 2, pp. 283\u0026ndash;291, Apr. 1993, doi: 10.1107/S0021889892009944.\u003c/li\u003e\n\u003cli\u003eJ. U. Bowie, R. L\u0026uuml;thy, and D. Eisenberg, \u0026lsquo;A Method to Identify Protein Sequences That Fold into a Known Three-Dimensional Structure\u0026rsquo;, \u003cem\u003eScience\u003c/em\u003e, vol. 253, no. 5016, pp. 164\u0026ndash;170, Jul. 1991, doi: 10.1126/science.1853201.\u003c/li\u003e\n\u003cli\u003eP. Benkert, S. C. E. Tosatto, and D. Schomburg, \u0026lsquo;QMEAN: A comprehensive scoring function for model quality assessment\u0026rsquo;, \u003cem\u003eProteins Struct. Funct. Bioinforma.\u003c/em\u003e, vol. 71, no. 1, pp. 261\u0026ndash;277, Apr. 2008, doi: 10.1002/prot.21715.\u003c/li\u003e\n\u003cli\u003eJ. Jumper \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Highly accurate protein structure prediction with AlphaFold\u0026rsquo;, \u003cem\u003eNature\u003c/em\u003e, vol. 596, no. 7873, pp. 583\u0026ndash;589, Aug. 2021, doi: 10.1038/s41586-021-03819-2.\u003c/li\u003e\n\u003cli\u003eR. A. Friesner \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein\u0026minus;Ligand Complexes\u0026rsquo;, \u003cem\u003eJ. Med. Chem.\u003c/em\u003e, vol. 49, no. 21, pp. 6177\u0026ndash;6196, Oct. 2006, doi: 10.1021/jm051256o.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: Glide, Schr\u0026ouml;dinger, LLC, New York, NY, 2022\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003eJ. A. Read, K. W. Wilkinson, R. Tranter, R. B. Sessions, and R. L. Brady, \u0026lsquo;Chloroquine Binds in the Cofactor Binding Site ofPlasmodium falciparum Lactate Dehydrogenase\u0026rsquo;, \u003cem\u003eJ. Biol. Chem.\u003c/em\u003e, vol. 274, no. 15, pp. 10213\u0026ndash;10218, Apr. 1999, doi: 10.1074/jbc.274.15.10213.\u003c/li\u003e\n\u003cli\u003eA. Cameron \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Identification and Activity of a Series of Azole-based Compounds with Lactate Dehydrogenase-directed Anti-malarial Activity\u0026rsquo;, \u003cem\u003eJ. Biol. Chem.\u003c/em\u003e, vol. 279, no. 30, pp. 31429\u0026ndash;31439, Jul. 2004, doi: 10.1074/jbc.M402433200.\u003c/li\u003e\n\u003cli\u003eR. Conners \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Mapping the binding site for gossypol-like inhibitors of Plasmodium falciparum lactate dehydrogenase\u0026rsquo;, \u003cem\u003eMol. Biochem. Parasitol.\u003c/em\u003e, vol. 142, no. 2, pp. 137\u0026ndash;148, Aug. 2005, doi: 10.1016/j.molbiopara.2005.03.015.\u003c/li\u003e\n\u003cli\u003eH. M. Berman, \u0026lsquo;The Protein Data Bank\u0026rsquo;, \u003cem\u003eNucleic Acids Res.\u003c/em\u003e, vol. 28, no. 1, pp. 235\u0026ndash;242, Jan. 2000, doi: 10.1093/nar/28.1.235.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: Protein Preparation Wizard; Epik, Schr\u0026ouml;dinger, LLC, New York, NY, 2022; Impact, Schr\u0026ouml;dinger, LLC, New York, NY; Prime, Schr\u0026ouml;dinger, LLC, New York, NY, 2022\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: LigPrep, Schr\u0026ouml;dinger, LLC, New York, NY, 2022\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: Maestro, Schr\u0026ouml;dinger, LLC, New York, NY, 2022\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003eR. A. Friesner \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Glide:\u0026thinsp; A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy\u0026rsquo;, \u003cem\u003eJ. Med. Chem.\u003c/em\u003e, vol. 47, no. 7, pp. 1739\u0026ndash;1749, Mar. 2004, doi: 10.1021/jm0306430.\u003c/li\u003e\n\u003cli\u003eA. Chang \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;BRENDA, the ELIXIR core data resource in 2021: new developments and updates\u0026rsquo;, \u003cem\u003eNucleic Acids Res.\u003c/em\u003e, vol. 49, no. D1, pp. D498\u0026ndash;D508, Jan. 2021, doi: 10.1093/nar/gkaa1025.\u003c/li\u003e\n\u003cli\u003eJ. C. Shelley, A. Cholleti, L. L. Frye, J. R. Greenwood, M. R. Timlin, and M. Uchimaya, \u0026lsquo;Epik: a software program for pK a prediction and protonation state generation for drug-like molecules\u0026rsquo;, \u003cem\u003eJ. Comput. Aided Mol. Des.\u003c/em\u003e, vol. 21, no. 12, pp. 681\u0026ndash;691, Dec. 2007, doi: 10.1007/s10822-007-9133-z.\u003c/li\u003e\n\u003cli\u003eG. Madhavi Sastry, M. Adzhigirey, T. Day, R. Annabhimoju, and W. Sherman, \u0026lsquo;Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments\u0026rsquo;, \u003cem\u003eJ. Comput. Aided Mol. Des.\u003c/em\u003e, vol. 27, no. 3, pp. 221\u0026ndash;234, Mar. 2013, doi: 10.1007/s10822-013-9644-8.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: QikProp, Schr\u0026ouml;dinger, LLC, New York, NY, 2022-4\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003eC. A. Lipinski, \u0026lsquo;Drug-like properties and the causes of poor solubility and poor permeability\u0026rsquo;, \u003cem\u003eJ. Pharmacol. Toxicol. Methods\u003c/em\u003e, vol. 44, no. 1, pp. 235\u0026ndash;249, Jul. 2000, doi: 10.1016/S1056-8719(00)00107-6.\u003c/li\u003e\n\u003cli\u003eA. E. Cho, V. Guallar, B. J. Berne, and R. Friesner, \u0026lsquo;Importance of accurate charges in molecular docking: Quantum mechanical/molecular mechanical (QM/MM) approach\u0026rsquo;, \u003cem\u003eJ. Comput. Chem.\u003c/em\u003e, vol. 26, no. 9, pp. 915\u0026ndash;931, Jul. 2005, doi: 10.1002/jcc.20222.\u003c/li\u003e\n\u003cli\u003eA. D. Bochevarov \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Jaguar: A high‐performance quantum chemistry software program with strengths in life and materials sciences\u0026rsquo;, \u003cem\u003eInt. J. Quantum Chem.\u003c/em\u003e, vol. 113, no. 18, pp. 2110\u0026ndash;2142, Sep. 2013, doi: 10.1002/qua.24481.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: Jaguar, Schr\u0026ouml;dinger, LLC, New York, NY, 2022\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003eJ. C. Phillips \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;Scalable molecular dynamics on CPU and GPU architectures with NAMD\u0026rsquo;, \u003cem\u003eJ. Chem. Phys.\u003c/em\u003e, vol. 153, no. 4, p. 044130, Jul. 2020, doi: 10.1063/5.0014475.\u003c/li\u003e\n\u003cli\u003eJ. Lee \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field\u0026rsquo;, \u003cem\u003eJ. Chem. Theory Comput.\u003c/em\u003e, vol. 12, no. 1, pp. 405\u0026ndash;413, Jan. 2016, doi: 10.1021/acs.jctc.5b00935.\u003c/li\u003e\n\u003cli\u003eN. M. O\u0026rsquo;Boyle, M. Banck, C. A. James, C. Morley, T. Vandermeersch, and G. R. Hutchison, \u0026lsquo;Open Babel: An open chemical toolbox\u0026rsquo;, \u003cem\u003eJ. Cheminformatics\u003c/em\u003e, vol. 3, no. 1, p. 33, Dec. 2011, doi: 10.1186/1758-2946-3-33.\u003c/li\u003e\n\u003cli\u003eW. Humphrey, A. Dalke, and K. Schulten, \u0026lsquo;VMD: Visual molecular dynamics\u0026rsquo;, \u003cem\u003eJ. Mol. Graph.\u003c/em\u003e, vol. 14, no. 1, pp. 33\u0026ndash;38, Feb. 1996, doi: 10.1016/0263-7855(96)00018-5.\u003c/li\u003e\n\u003cli\u003eK. Vanommeslaeghe \u003cem\u003eet al.\u003c/em\u003e, \u0026lsquo;CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields\u0026rsquo;, \u003cem\u003eJ. Comput. Chem.\u003c/em\u003e, p. NA-NA, 2009, doi: 10.1002/jcc.21367.\u003c/li\u003e\n\u003cli\u003e\u003cem\u003eSchr\u0026ouml;dinger Release 2022-4: Desmond Molecular Dynamics System, D. E. Shaw Research, New York, NY, 2022\u003c/em\u003e.\u003c/li\u003e\n\u003cli\u003eM. P. Jacobson, R. A. Friesner, Z. Xiang, and B. Honig, \u0026lsquo;On the Role of the Crystal Environment in Determining Protein Side-chain Conformations\u0026rsquo;, \u003cem\u003eJ. Mol. Biol.\u003c/em\u003e, vol. 320, no. 3, pp. 597\u0026ndash;608, Jul. 2002, doi: 10.1016/S0022-2836(02)00470-9.\u003c/li\u003e\n\u003cli\u003eNaruya Saitou and M. Nei, \u0026lsquo;The neighbor-joining method: a new method for reconstructing phylogenetic trees.\u0026rsquo;, \u003cem\u003eMol. Biol. Evol.\u003c/em\u003e, Jul. 1987, doi: 10.1093/oxfordjournals.molbev.a040454.\u003c/li\u003e\n\u003cli\u003eJ. I. Boucher, J. R. Jacobowitz, B. C. Beckett, S. Classen, and D. L. Theobald, \u0026lsquo;An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases\u0026rsquo;, \u003cem\u003eeLife\u003c/em\u003e, vol. 3, p. e02304, Jun. 2014, doi: 10.7554/eLife.02304.\u003c/li\u003e\n\u003cli\u003eK. L. Kavanagh, R. A. Elling, and D. K. Wilson, \u0026lsquo;Structure of \u003cem\u003eToxoplasma gondii\u003c/em\u003e LDH1: Active-Site Differences from Human Lactate Dehydrogenases and the Structural Basis for Efficient APAD \u003csup\u003e+\u003c/sup\u003e Use \u003csup\u003e,\u003c/sup\u003e\u0026rsquo; \u003cem\u003eBiochemistry\u003c/em\u003e, vol. 43, no. 4, pp. 879\u0026ndash;889, Feb. 2004, doi: 10.1021/bi035108g.\u003c/li\u003e\n\u003cli\u003eC. Colovos and T. O. Yeates, \u0026lsquo;Verification of protein structures: patterns of nonbonded atomic interactions\u0026rsquo;, \u003cem\u003eProtein Sci. Publ. Protein Soc.\u003c/em\u003e, vol. 2, no. 9, pp. 1511\u0026ndash;1519, Sep. 1993, doi: 10.1002/pro.5560020916.\u003c/li\u003e\n\u003cli\u003eM. Wiederstein and M. J. Sippl, \u0026lsquo;ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins\u0026rsquo;, \u003cem\u003eNucleic Acids Res.\u003c/em\u003e, vol. 35, no. Web Server issue, pp. W407-410, Jul. 2007, doi: 10.1093/nar/gkm290.\u003c/li\u003e\n\u003cli\u003eP. Benkert, S. C. E. Tosatto, and D. Schomburg, \u0026lsquo;QMEAN: A comprehensive scoring function for model quality assessment\u0026rsquo;, \u003cem\u003eProteins\u003c/em\u003e, vol. 71, no. 1, pp. 261\u0026ndash;277, Apr. 2008, doi: 10.1002/prot.21715.\u003c/li\u003e\n\u003cli\u003eJ. U. Bowie, R. L\u0026uuml;thy, and D. Eisenberg, \u0026lsquo;A method to identify protein sequences that fold into a known three-dimensional structure\u0026rsquo;, \u003cem\u003eScience\u003c/em\u003e, vol. 253, no. 5016, pp. 164\u0026ndash;170, Jul. 1991, doi: 10.1126/science.1853201.\u003c/li\u003e\n\u003cli\u003eR. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, \u0026lsquo;PROCHECK: a program to check the stereochemical quality of protein structures\u0026rsquo;, \u003cem\u003eJ. Appl. Crystallogr.\u003c/em\u003e, vol. 26, no. 2, pp. 283\u0026ndash;291, Apr. 1993, doi: 10.1107/S0021889892009944.\u003c/li\u003e\n\u003cli\u003eA. Volkamer, D. Kuhn, F. Rippmann, and M. Rarey, \u0026lsquo;DoGSiteScorer: a web server for automatic binding site prediction, analysis and druggability assessment\u0026rsquo;, \u003cem\u003eBioinformatics\u003c/em\u003e, vol. 28, no. 15, pp. 2074\u0026ndash;2075, Aug. 2012, doi: 10.1093/bioinformatics/bts310.\u003c/li\u003e\n\u003cli\u003eW. L. Jorgensen and E. M. Duffy, \u0026lsquo;Prediction of drug solubility from Monte Carlo simulations\u0026rsquo;, \u003cem\u003eBioorg. Med. Chem. Lett.\u003c/em\u003e, vol. 10, no. 11, pp. 1155\u0026ndash;1158, Jun. 2000, doi: 10.1016/S0960-894X(00)00172-4.\u003c/li\u003e\n\u003cli\u003eA. Cavalli, E. Poluzzi, F. De Ponti, and M. Recanatini, \u0026lsquo;Toward a Pharmacophore for Drugs Inducing the Long QT Syndrome: Insights from a CoMFA Study of HERG K \u003csup\u003e+\u003c/sup\u003e Channel Blockers\u0026rsquo;, \u003cem\u003eJ. Med. Chem.\u003c/em\u003e, vol. 45, no. 18, pp. 3844\u0026ndash;3853, Aug. 2002, doi: 10.1021/jm0208875.\u003c/li\u003e\n\u003cli\u003eM. Yazdanian, S. L. Glynn, J. L. Wright, and A. Hawi, \u0026lsquo;Correlating Partitioning and Caco-2 Cell Permeability of Structurally Diverse Small Molecular Weight Compounds\u0026rsquo;, \u003cem\u003ePharm. Res.\u003c/em\u003e, vol. 15, no. 9, pp. 1490\u0026ndash;1494, 1998, doi: 10.1023/A:1011930411574.\u003c/li\u003e\n\u003cli\u003eJ. M. Luco, \u0026lsquo;Prediction of the Brain\u0026minus;Blood Distribution of a Large Set of Drugs from Structurally Derived Descriptors Using Partial Least-Squares (PLS) Modeling\u0026rsquo;, \u003cem\u003eJ. Chem. Inf. Comput. Sci.\u003c/em\u003e, vol. 39, no. 2, pp. 396\u0026ndash;404, Mar. 1999, doi: 10.1021/ci980411n.\u003c/li\u003e\n\u003cli\u003eG. Colmenarejo, A. Alvarez-Pedraglio, and J.-L. Lavandera, \u0026lsquo;Cheminformatic Models To Predict Binding Affinities to Human Serum Albumin\u0026rsquo;, \u003cem\u003eJ. Med. Chem.\u003c/em\u003e, vol. 44, no. 25, pp. 4370\u0026ndash;4378, Dec. 2001, doi: 10.1021/jm010960b.\u003cstrong\u003e\u003c/strong\u003e\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"parasitology-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"pare","sideBox":"Learn more about [Parasitology Research](http://link.springer.com/journal/436)","snPcode":"436","submissionUrl":"https://submission.nature.com/new-submission/436/3","title":"Parasitology Research","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Lactate dehydrogenase, Babesia bigemina, drug discovery","lastPublishedDoi":"10.21203/rs.3.rs-4951492/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4951492/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e \u003cem\u003eBabesia bigemina\u003c/em\u003e is an apicomplexan parasite and causes \u0026ldquo;Texas fever\u0026rdquo; in bovines. Lactate dehydrogenase (LDH) is an essential enzyme in regulating anaerobic metabolism. The presence of five amino acid insertion in the active site of many apicomplexan LDH makes a significant difference between the parasitic LDH and the host counterpart. Therefore, apicomplexan LDH is an attractive drug target. In this study, a structure-based drug discovery approach was performed to find novel inhibitor candidates. In the first round, possible candidates were identified by following the virtual screening workflow. Then, the compounds with favourable docking scores were filtered using the QM-polarized ligand docking and induced fit docking methods. As a result, 20 novel compounds that bind to the active site of the BbigLDH and have low affinity with the host LDHs have been identified. Molecular dynamics simulations of the complexes (in total 8.8 \u0026micro;s) were performed and binding free energies were calculated. In conclusion, compounds named C09, C16 and C18 deserve further investigation to better understand their potential therapeutic effects on babesiosis. The findings of this study, elucidating the structural properties of BbigLDH enzyme and discovering its potential inhibitors, might pave the way for further research in developing LDH-targeted therapeutic interventions.\u003c/p\u003e","manuscriptTitle":"Exploring Novel Inhibitors for Babesia bigemina Lactate Dehydrogenase: A Computational Structural Biology Perspective","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-24 09:19:44","doi":"10.21203/rs.3.rs-4951492/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-09-01T11:04:47+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-08-27T06:05:47+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-08-26T11:09:42+00:00","index":"","fulltext":""},{"type":"submitted","content":"Parasitology Research","date":"2024-08-21T12:04:42+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"parasitology-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"pare","sideBox":"Learn more about [Parasitology Research](http://link.springer.com/journal/436)","snPcode":"436","submissionUrl":"https://submission.nature.com/new-submission/436/3","title":"Parasitology Research","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"e4b73380-2cc2-494b-b0f1-0d5a2402a029","owner":[],"postedDate":"September 24th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2025-01-13T16:07:57+00:00","versionOfRecord":{"articleIdentity":"rs-4951492","link":"https://doi.org/10.1007/s00436-024-08433-5","journal":{"identity":"parasitology-research","isVorOnly":false,"title":"Parasitology Research"},"publishedOn":"2025-01-07 15:58:03","publishedOnDateReadable":"January 7th, 2025"},"versionCreatedAt":"2024-09-24 09:19:44","video":"","vorDoi":"10.1007/s00436-024-08433-5","vorDoiUrl":"https://doi.org/10.1007/s00436-024-08433-5","workflowStages":[]},"version":"v1","identity":"rs-4951492","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4951492","identity":"rs-4951492","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}