Density Functional Theory Calculations and Molecular Docking Analyses of Flavonoids for Their Possible Application against the Acetylcholinesterase and Triose-Phosphate Isomerase Proteins of Rhipicephalus microplus

Ticks and tick-borne diseases constitute a substantial hazard to the livestock industry. The rising costs and lack of availability of synthetic chemical acaricides for farmers with limited resources, tick resistance to current acaricides, and residual issues in meat and milk consumed by humans further aggravate the situation. Developing innovative, eco-friendly tick management techniques, such as natural products and commodities, is vital. Similarly, searching for effective and feasible treatments for tick-borne diseases is essential. Flavonoids are a class of natural chemicals with multiple bioactivities, including the inhibition of enzymes. We selected eighty flavonoids having enzyme inhibitory, insecticide, and pesticide properties. Flavonoids’ inhibitory effects on the acetylcholinesterase (AChE1) and triose-phosphate isomerase (TIM) proteins of Rhipicephalus microplus were examined utilizing a molecular docking approach. Our research demonstrated that flavonoids interact with the active areas of proteins. Seven flavonoids (methylenebisphloridzin, thearubigin, fortunellin, quercetagetin-7-O-(6-O-caffeoyl-β-d-glucopyranoside), quercetagetin-7-O-(6-O-p-coumaroyl-β-glucopyranoside), rutin, and kaempferol 3-neohesperidoside) were the most potent AChE1 inhibitors, while the other three flavonoids (quercetagetin-7-O-(6-O-caffeoyl-β-d-glucopyranoside), isorhamnetin, and liquiritin) were the potent inhibitors of TIM. These computationally-driven discoveries are beneficial and can be utilized in assessing drug bioavailability in both in vitro and in vivo settings. This knowledge can create new strategies for managing ticks and tick-borne diseases.


Introduction
Ticks are hematophagous vertebrate ectoparasites that transmit bacterial, viral, and protozoan infections [1]. Cattle farmers worldwide are becoming increasingly concerned about the spread of the southern cow tick, Rhipicephalus microplus (formerly Boophilus microplus). When a cow tick bites an animal, it can transmit disease-causing microbes such 2 of 20 as Babesia bovis, Babesia bigemin, and Anaplasma marginale, causing the animal to grow slowly, produce less milk, or even die [2]. R. microplus control has emerged as a central concern and area of attention in animal parasitology research and technology development due to its global distribution and impact on cattle husbandry [3]. While chemical acaricides have effectively reduced R. microplus populations, their widespread use has resulted in acquired resistance, making population control less likely [4]. Multiple cases of R. microplus acaricide resistance have emerged in recent years, making the global spread of acaricide resistance a major concern for the cattle industry [5]. As a result, research into new chemicals with distinct mechanisms of action and acaricidal efficacy against R. microplus is ongoing.
Moraes et al. investigated the physiological mechanisms underlying tick embryo formation to develop new tick population control strategies [6]. They found that morphogenetic alterations in R. microplus embryos were linked to the predominant energy source [6]. While R. microplus embryos develop in 21 days, they are syncytia until the fifth day. On the sixth day, zygotic expression becomes active, and the embryo assumes control over its development [7]. There seem to be two separate metabolic stages during cellularization, with major changes occurring in the glycolytic and gluconeogenic pathways [8].
Triose-phosphate isomerase (TIM), an enzyme involved in glycolysis and gluconeogenesis, catalyzes the interconversion of glyceraldehyde-3-phosphate and dihydroxyacetone phosphate. The therapeutic development of TIM has shown promise against several human-pathogenic parasites [9,10]. These include Plasmodium falciparum, Trypanosoma brucei, Trypanosoma cruzi, Giardia lamblia, and Fasciola hepatica. TIM is a (β/α) 8 -barrel protein, a core of eight parallel β-strands and eight α-helices on all sides. Banner et al. found TIM's "TIM-barrel" crystal structure [11]. Moraes et al. were the first to investigate TIM's kinetic and structural properties from embryos of R. microplus [6]. Their study suggested that many cysteine residues in R. microplus TIM (RmTIM) could be used to make antiparasitic drugs. In addition, selective inhibitors can be created by targeting this enzyme's less conserved interface [12].
Acetylcholinesterase (AChE1) (E.C. 3.1.1.7.) is a synaptic enzyme and a major target of both organophosphates (OPs) and carbamates [13]. It is essential for the transmission of nerve impulses in both vertebrates and arthropods [14,15]. AChE1 inhibits nerve impulses by breaking acetylcholine into acetate and choline [16]. AChE1 is inhibited by the acaricides carbamate and OPs. Carbamates and OPs pesticides both inhibit AChE1 activity by binding to the esterasic and anionic sites of the enzyme. Both suppress the enzyme's activity in a similar mode of action [17,18]. Various flavonoids and their derivatives have been found to possess inhibitory potential on AChE in Caenorhabditis elegans and Spodoptera litura [19], as well as electric eels (Electrophorus electricus) [20].
The majority of the organic compounds produced by plants do not appear to contribute to their growth and development. The distribution of secondary metabolites, or flavonoids, varies considerably between taxonomic plant groups [21]. The interaction of flavonoids with proteins and nucleic acids has resulted in the development of numerous bioactive compounds with pharmacological, antibacterial, and insecticidal properties. Since flavonoids are used as pesticides in medicine and agriculture, they are significant. Therefore, they could be useful in a pest control program [22].
The primary objective of computational approaches such as molecular docking and in silico absorption, distribution, metabolism, excretion, and toxicity (ADMET) analyses is to screen potential flavonoids against RmTIM and R. microplus AChE1 (RmAChE1) in numerous databases and libraries. Computational screening reduces the time and resources required for experimental testing in pharmaceutical research. Utilizing computational designs and molecular information on flavonoid compounds, research into the development of effective alternatives to currently used acaricides may benefit from the wide variety of flavonoids by selecting the most promising compounds. Additionally, medicinal chemistry is frequently used to identify new bioactive molecules with acaricidal applications. For example, to identify RmTIM and RmAChE1 inhibitors, natural substances from numer-ous databases and libraries were screened using knowledge of parasiticides and enzyme inhibitors.

The Structural Models of TIM and RmAChE1
TIM crystallographic enzyme structures have been obtained from the Protein Data Bank (PDB) repository (PDB code: 3TH6). Using the Python molecular viewer tool of PyMOL (v1.9, http://www.pymol.org (accessed on 10 June 2022)), water molecules were removed from the protein structure, and non-polar hydrogen molecules were added. The amino acid sequence of RmAChE1, with accession number A0A0F6P2D6, was retrieved from the UniProt database (found at http://www.uniprot.org (accessed on 10 June 2022)), which is a comprehensive source for protein sequences. The Basic Local Alignment Search Tool (BLAST) analysis revealed the existence of 100 protein sequences that showed significant similarities. Out of these, 48 were AChE1s from other species, 21 were butyrylcholinesterases, and the remaining 31 were proteins that shared some degree of identity with RmAChE1. The crystal structure of mouse AChE1 (PDB: 5DTI, Chain A) was selected as the template among the 48 AChEs due to its highest resolution of 2 Å. Additionally, this structure showed a high level of similarity to RmAChE1, with 44% identity and 89% query coverage. Figure 1 demonstrates the three-dimensional (3D) structures of RmTIM and RmAChE1.
computational designs and molecular information on flavonoid compounds, research into the development of effective alternatives to currently used acaricides may benefit from the wide variety of flavonoids by selecting the most promising compounds. Additionally, medicinal chemistry is frequently used to identify new bioactive molecules with acaricidal applications. For example, to identify RmTIM and RmAChE1 inhibitors, natural substances from numerous databases and libraries were screened using knowledge of parasiticides and enzyme inhibitors.

The Structural Models of TIM and RmAChE1
TIM crystallographic enzyme structures have been obtained from the Protein Data Bank (PDB) repository (PDB code: 3TH6). Using the Python molecular viewer tool of PyMOL (v1.9, http://www.pymol.org (accessed on 10 June 2022)), water molecules were removed from the protein structure, and non-polar hydrogen molecules were added. The amino acid sequence of RmAChE1, with accession number A0A0F6P2D6, was retrieved from the UniProt database (found at http://www.uniprot.org (accessed on 10 June 2022)), which is a comprehensive source for protein sequences. The Basic Local Alignment Search Tool (BLAST) analysis revealed the existence of 100 protein sequences that showed significant similarities. Out of these, 48 were AChE1s from other species, 21 were butyrylcholinesterases, and the remaining 31 were proteins that shared some degree of identity with RmAChE1. The crystal structure of mouse AChE1 (PDB: 5DTI, Chain A) was selected as the template among the 48 AChEs due to its highest resolution of 2 Å. Additionally, this structure showed a high level of similarity to RmAChE1, with 44% identity and 89% query coverage. Figure 1 demonstrates the three-dimensional (3D) structures of RmTIM and RmAChE1.

Validation of the Modeled Structures
To optimize 3D models for structural inspections, PROCHECK created a Ramachandran map. Non-glycine and non-proline residues' contributions to the Ramachandran map's u and w distribution were summarized in Figure 2A. The protein structural models were checked against protein structures in the PDB using the ProSAweb server based on z-scores ( Figure 2C) [23]. The application provided the Z-score of the input structure and a plot of the residue energies. SWISS-MODEL's Z-score for representing AChE1 showed a −9.44, indicating excellent overall model quality. Figure

Validation of the Modeled Structures
To optimize 3D models for structural inspections, PROCHECK created a Ramachandran map. Non-glycine and non-proline residues' contributions to the Ramachandran map's u and w distribution were summarized in Figure 2A. The protein structural models were checked against protein structures in the PDB using the ProSA-web server based on z-scores ( Figure 2C) [23]. The application provided the Z-score of the input structure and a plot of the residue energies. SWISS-MODEL's Z-score for representing AChE1 showed a −9.44, indicating excellent overall model quality. Figure 2B displays the ERRAT-verified 3D structures. The Verify3D was used to evaluate the 3D model's quality. The 3D models scored 96.26% on this experiment's evaluation scale, where the minimum passing score was 80% ( Figure 2D). displays the ERRAT-verified 3D structures. The Verify3D was used to evaluate the 3D model's quality. The 3D models scored 96.26% on this experiment's evaluation scale, where the minimum passing score was 80% ( Figure 2D).

Figure 2.
The predicted 3D structures of the Rhipicephalus microplus acetylcholinesterase 1 (RmAChE1) protein were verified using (A) the Ramachandran plot produced by PROCHECK, where the red, yellow and black colors represent most favorable, favorable, and disallowed region respectively, Phi and Psi bond represent torsion angle which predict the possible conformation of the peptides. (B) the ERRAT server quality factor, * On the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. ** Expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. Good high-resolution structures generally produce values around 95% or higher. For lower resolutions (2.5 to 3 Å) the average overall quality factor is around 91%. The red and yellow bars represent the 99% and 95% confidence respectively. (C) the ProSA-web (Z-score), and (D) the Verify3D.

Active Site Prediction
In the context of docking studies, it is imperative to pinpoint the putative binding pockets within the modeled protein structures. The 3D structure of RmAChE1 has not been documented in available databases, resulting in the absence of any reported catalytic sites. As a result, it is imperative to anticipate the potential binding regions of the receptor and determine the area with the largest cavity for docking compounds. For this purpose, we utilized the site finder tool integrated into MOE to determine the probable binding sites within RmAChE1. The analysis resulted in the identification of eight high-potential binding pockets. After evaluating the potential binding pockets, we selected the one that comprises the following amino acids: GLN122, VAL123, LEU124, ASP125, THR126, LEU127, SER134, TRP137, ASN138, ALA139, TYR173, GLY174, GLY175, GLY176, TYR178, Ser179, GLY180, THR181, LEU184, TYR187, GLU255, SER256, TRP289, THR335, ASN336, SER337, GLY338, GLY339, VAL340, VAL341, ASP342, PHE343, PRO344, TRP384, and PHE385. Additionally, utilizing the procedure reported by Empereur-Mot et al. (2015), a test of activity predictability was used to confirm the active site of the modeled protein.
The aforementioned active sites were docked with active molecules [24]. Hence, as demonstrated in Figure 3, every active compound was precisely docked into the Figure 2. The predicted 3D structures of the Rhipicephalus microplus acetylcholinesterase 1 (RmA-ChE1) protein were verified using (A) the Ramachandran plot produced by PROCHECK, where the red, yellow and black colors represent most favorable, favorable, and disallowed region respectively, Phi and Psi bond represent torsion angle which predict the possible conformation of the peptides. (B) the ERRAT server quality factor, * On the error axis, two lines are drawn to indicate the confidence with which it is possible to reject regions that exceed that error value. ** Expressed as the percentage of the protein for which the calculated error value falls below the 95% rejection limit. Good highresolution structures generally produce values around 95% or higher. For lower resolutions (2.5 to 3 Å) the average overall quality factor is around 91%. The red and yellow bars represent the 99% and 95% confidence respectively. (C) the ProSA-web (Z-score), and (D) the Verify3D.

Active Site Prediction
In the context of docking studies, it is imperative to pinpoint the putative binding pockets within the modeled protein structures. The 3D structure of RmAChE1 has not been documented in available databases, resulting in the absence of any reported catalytic sites. As a result, it is imperative to anticipate the potential binding regions of the receptor and determine the area with the largest cavity for docking compounds. For this purpose, we utilized the site finder tool integrated into MOE to determine the probable binding sites within RmAChE1. The analysis resulted in the identification of eight high-potential binding pockets. After evaluating the potential binding pockets, we selected the one that comprises the following amino acids: GLN122, VAL123, LEU124, ASP125, THR126, LEU127, SER134, TRP137, ASN138, ALA139, TYR173, GLY174, GLY175, GLY176, TYR178, Ser179, GLY180, THR181, LEU184, TYR187, GLU255, SER256, TRP289, THR335, ASN336, SER337, GLY338, GLY339, VAL340, VAL341, ASP342, PHE343, PRO344, TRP384, and PHE385. Additionally, utilizing the procedure reported by Empereur-Mot et al. (2015), a test of activity predictability was used to confirm the active site of the modeled protein. The aforementioned active sites were docked with active molecules [24]. Hence, as demonstrated in Figure 3, every active compound was precisely docked into the suggested active site, while inactive ones could not be docked at the intended active site ( Figure 3). In the case of TIM, the catalytic site has been previously reported [25]. As such, we selected the corresponding residues at suggested active site, while inactive ones could not be docked at the intended a ( Figure 3). In the case of TIM, the catalytic site has been previously reported [25]. we selected the corresponding residues at the catalytic site and performed simulations with the related compounds. Figure 3 demonstrates the active sites structures of RmTIM and RmAChE1 ( Figure 4).

Interaction Pattern of Virtual Hits
The interaction modes of the RmAChE1 and RmTIM proteins with fl derivatives have been determined through PLIP https://plip-tool.biotec.tu-dr (accessed on 27 Febuary 2023), an online tool for protein-ligand interactions shows the selected flavonoid molecules having maximal interaction sco RmAChE1 and RmTIM. suggested active site, while inactive ones could not be docked at the intended ( Figure 3). In the case of TIM, the catalytic site has been previously reported [2 we selected the corresponding residues at the catalytic site and perform simulations with the related compounds. Figure 3 demonstrates the active sit structures of RmTIM and RmAChE1 (Figure 4).

Interaction Pattern of Virtual Hits
The interaction modes of the RmAChE1 and RmTIM proteins with derivatives have been determined through PLIP https://plip-tool.biotec.tu (accessed on 27 Febuary 2023), an online tool for protein-ligand interactio shows the selected flavonoid molecules having maximal interaction s RmAChE1 and RmTIM.

Interaction Pattern of Virtual Hits
The interaction modes of the RmAChE1 and RmTIM proteins with flavonoid derivatives have been determined through PLIP https://plip-tool.biotec.tu-dresden.de (accessed on 27 Febuary 2023), an online tool for protein-ligand interactions. Table 1 shows the selected flavonoid molecules having maximal interaction scores with RmAChE1 and RmTIM. Thearubigin 100945367   Figure 5 illustrates the interaction mode of (A) methylenebisphloridzin, (B) thearubigin, and (C) fortunellin with RmAChE1. The methylenebisphloridzin molecule was found to create seven hydrogen bonds with specific amino acids (Gly 180,     tionally, isorhamnetin also shows hydrophobic interactions with Ile 170, Lys 274, Ser 211, and Val 231. Liquiritin is observed to interact with specific amino acid residues in its binding site through both hydrophobic contacts and hydrogen bonds. Specifically, it forms hydrophobic contacts with Lys 174, Ile 170, Glu 239, Ser 235, and Ala 234, as well as five hydrogen bonds with Lys 13, Glu 97, Ser 211, Asn 213, and Lys 237. These interactions occur at distances of approximately 2.81, 2.26, 2.47, 2.32, and 2.65 Å, respectively. The docking interactions of other selected flavonoids with RmAChE1 is shown in Figures S2-S5.

Results of ADMET Calculation
Drug discovery and development is a complex process, and the most crucial yet challenging step is the conduct of DMPK (drug metabolism and pharmacokinetics) studies. This step accounts for the failure of approximately 60% of drugs during clinical phases. ADMET (absorption, distribution, metabolism, excretion, and toxicity) studies are a vital component of pharmacokinetics/pharmacology, describing the disposition of drug compounds in the body. The ADMET Predictor is a computer-designed program that utilizes molecular structures to estimate the pharmacokinetic properties/parameters of drug-like compounds, as outlined by [26].
The freely available Swiss ADME web tool is a software application that predicts the physicochemical properties, absorption, distribution, metabolism, elimination, and pharmacokinetic properties of molecules, all of which are critical determinants for clinical trials. This tool considers six physicochemical properties that are essential, including lipophilicity, flexibility, saturation, polarity, solubility, and size, as stated by [27].
The results of the ADMET evaluation provided insights into the physicochemical properties of the designed compounds, including compliance with the rules of five (MW, iLOGP, HBAs, and HBDs), as well as other parameters such as molecular polar surface area (TPSA), Blood-Brain Barrier (BBB) permeability, number of aromatic heavy atoms, and the presence of alerts for undesirable substructures (i.e., PAINS #alert), among others, and are presented in Table 2.

Density Functional Theory (DFT) Results
The frontier molecular orbitals HOMO and LUMO convey valuable data about any molecule's reactivity. The difference in HOMO and LUMO energies (Eg) plays a significant role in comprehending the chemical reactivity and kinetic stability. It was a proportional function between the energy band gap and chemical reactivity and stability. For example, if the band gap in a molecule has a higher value, then the molecule is less polarizable and shows high kinetic stability and low chemical reactivity (hard molecule). The HOMO and LUMO of liquiritin, isorhamnetin, and complex were made. Energy details are given in Table 3, such as EHOMO, ELUMO, and Eg. Liquiritin and isorhamnetin possess band gaps of 0.16 eV and 0.005 eV, respectively, having the constituent energies of HOMO and LUMO, −0.33641 eV, −0.32355 eV and −0.20637 eV, −0.31813, respectively. The band gap shown by the complex was lower than the obtained band gaps of liquiritin, making it more reactive towards the affected area. While in the case of isorhamnetin, the band gap of complexes is more than isorhamnetin itself, making it stable reactive. This much smaller band gap of complex predicts the need for a small quantity of energy to get excited from ground level. Furthermore, this lesser energy difference shows the extent of conjugation within the LIQ and ISO (Tables 4 and 5), whereas the Table 6 represents the HOMO-LUMO of selected flavonoids.

Discussion
Nowadays, synthetic insecticides are often employed to eliminate veterinary and pharmaceutical pests. Pest resistance, the product remains, active component withdrawal, environmental persistence, and nontarget dangers need new solutions. One approach uses plant-based chemicals' toxic, irritant, or beneficial qualities. These compounds have limited human harm, short environmental persistence, and complicated chemistry, making them ideal pesticide candidates [28]. New bioactive compounds that are more potent and selective for tick targets have been explored as potential solutions to these problems [29]. In silico methods like molecular docking, homology modeling, and molecular dynamics make it easier to find novel compounds that bind to parasites' molecular targets [30].
The neuronal enzyme AChE1, which catalyzes the breakdown of acetylcholine, is a recognized chemical target of ticks [31]. It causes neuromuscular paralysis and death by prolonging the tick's neural excitation. Moreover, TIM is an enzyme involved in both glycolysis and gluconeogenesis since it catalyzes the conversion of 3-phosphoglycerate to dihydroxyacetone phosphate. To generate selective inhibitors, one must target the enzyme's interface, which is poorly conserved between species, even if the structural similarity of the enzyme is well conserved across species [12]. However, whether a specific molecular inhibitor of RmTIM has a flavonoid structure or any other structure is still unknown.
Antiparasitic activity and enzyme inhibition are two important biological effects that piqued modern interest in the flavonoid family of natural secondary metabolites. Plant products' acaricidal efficacy has been connected to the presence of flavonoids. Ghosh et al. discovered that adult Ricinus microplus might be destroyed in vitro by an ethanolic extract of Ricinus communis leaf containing flavonoids (quercetin, flavone, and kaempferol) [32]. In addition, the active fractions of Ocotea aciphylla leaf containing a flavonoid (vitexin-2 -O-rhamnoside 9) showed in vitro anti-tick efficacy against larvae of R. microplus [33]. Flavonoids' capacity to bind to enzymes like AChE1 is another important biological characteristic. Researchers have shown that flavonoids (quercetin, rutin, narigin, and hesperidin) may inhibit AChE1 from Electrophorus electricus in vitro [34,35]. In silico and in vitro studies have shown scopoletin and diosmin to inhibit mouse AChE1 [36]. Despite this, there is a lack of research on plant flavonoids' effect on the acaricidal effectiveness of RmAChE1 and RmTIM proteins.
Plant flavonoids' docking with the RmAChE1 and RmTIM proteins was examined in this work. Docking analysis requires a protein's 3D structure. Since the structure of the RmAChE1 protein has not been identified experimentally, we constructed a model of this protein using the homology modeling approach, widely accepted as the best technique for modeling proteins. The structure of RmTIM was obtained from the PDB. Multiple webbased tools assessed the accuracy and stereochemistry of the RmAChE1 model developed by the computer. Docking analysis requires a protein structure file in (.pdb) format. Docking analysis foretells a compound's binding energy to a protein and the optimal molecular orientation for binding.
Our primary objective in this work was to determine whether natural flavonoids had any effect on the RmAChE1 and RmTIM proteins. This study retrieved 80 flavonoid structures from PubChem and PubMed. After that, these flavonoids were put through a virtual screening using MOE software, which tested them against both proteins. Ten flavonoids, including methylenebisphloridzin, thearubigin, fortunellin, quercetagetin-7-O-(6-O-caffeoylβ-D-glucopyranoside) 1, quercetagetin-7-O-(6-O-p-coumaroyl-β-glucopyranoside), rutin, kaempferol 3-neohesperidoside, quercetagetin-7-O-(6-O-caffeoyl-β-D-glucopyranoside) 2, isorhamnetin, and liquiritin were obtained as potential inhibitors (the flavonoids 1-7 against RmAChE1; the flavonoids 8-10 against RmTIM) based on best docking scores against the active sites of both proteins of R. microplus. Many different pharmacological effects, such as antidiabetic, anti-inflammatory, anti-hyperglycemic, anti-cancer, antibacterial [37], AChE1 inhibitory [19], enzymatic [38], hepatoprotective [39], and antioxidant [40], may be applicable. Furthermore, we determined the drug-likeness properties using the available online tool SwissADME comprised of pharmacodynamics and pharmacokinetics and, ultimately, provided the prediction regarding absorption, distribution, metabolism, and excretion (ADME) in the human body like a drug. We incorporated the drug-like calculation based on the Lipinski rule of five to determine oral absorption or membrane permeability ( Table 2). Following Lipinski's rule of five and the concept of QED, as outlined in Table 2, all the designed compounds violated more than one rule, except Isorhamnetin and Liquirtin, which did not violate any rule. This implies that the MW, RB, HBD, HBA, TPSA, and iLOGP of some compounds fall within the acceptable range. Additionally, only one PAINS alert was identified for compounds 2, 4, 5, 6, and 8, indicating their specificity. Thus, it can be concluded that these active anti-tick compounds (1 to 10) possess an average pharmacokinetic profile.
The DFT results also revealed that both compounds (isorhamnetin and liquiritin) have better chemical reactivity and kinetic stability with considerable intramolecular charge transfer between electron-donor and electron-acceptor groups. The molecule with a higher band gap was less polarizable and showed high kinetic stability and low chemical reactivity (hard molecule). The band gap shown by the complex was lower than the obtained band gaps of liquiritin, which makes it more reactive towards the affected area. While in the case of isorhamnetin, the band gap of complexes was more than isorhamnetin itself, which makes it stable and reactive. This much smaller band gap of complex predicts the need for a small quantity of energy to get excited from ground level [41].
The stability and biological activity of flavonoids can be influenced by temperature, and their sensitivity to heat treatment can vary based on their chemical structure. Generally, glycosylated flavonoids exhibit more excellent heat treatment resistance than aglycone flavonoids [42]. Although the mechanism is unknown, incorporating acyl groups has been found to increase the thermostability of flavonoids [43]. For instance, Ishihara and Nakajima conducted a study where quercetin-3-glucoside was monoacylated in vitro using nine different aromatic carboxylic acids. The resulting flavonoids showed improved thermostability and light resistivity [44].
Based on our findings, flavonoids may function as inhibitors of the RmAChE1 and RmTIM proteins. Researchers studying plant flavonoids may gain new insight into their abilities to inhibit the RmAchE1 and RmTIM for treating R. microplus infection. Further research on flavonoid inhibition may lead to their use as acaricides in creating novel antiparasitic medications.

Preparation of TIM and the Homology Model of RmAChE1
TIM's 3D structure [PDB: 3TH6] was retrieved from the PDB archive of the Research Collaboratory for Structural Bioinformatics (RCSB) (https://www.rcsb.org (accessed on 10 June 2022)). To prepare the TIM protein structure for molecular docking, crystal water molecules were eliminated. Some residues were missing in the structure of TIMs, so we constructed the missing residues using a loop modeler implemented in MOE software.
Since the 3D structure of the RmAChE1 protein was not reported, we needed to model the structure of RmAChE1. The amino acid sequence of RmAChE1 was retrieved from uniport https://www.uniprot.org/uniprotkb/ (accessed on 11 June 2022) under accession No. A0A0F6P2D6 RHIMP and then modeled using SWISS-MODEL https://swissmodel. expasy.org/interactive/ (accessed on 12 June 2022). To identify appropriate templates for constructing the 3D structure of RmAChE1, a sequence similarity search was carried out using the NCBI BLAST against the PDB database. Sequence alignment analysis was conducted through the Clustal Omega server [45]. In addition, the stereochemical quality of the modeled structure of the test protein was validated using the PROCHECK, Verify3D, and ERRAT services.

Prediction of the Active Site
RmAChE1 lacks a determined 3D structure, and thus its active site has not been experimentally identified. Therefore, we utilize site finder, a computational tool implemented in the molecular operating environment, to evaluate the potential active site. These site prediction tools provide a theoretical assessment of the active site in RmAChE1. In addition, TIM is a dimeric protein comprising two individual monomers, each containing a TIM barrel structure. The active site responsible for the catalytic activity is situated within this TIM barrel. Given the high similarity of the catalytic sites among all TIM structures, a set of catalytic residues was selected for further investigation through molecular docking studies.

Ligands Searching and Database Preparation
An extensive literature survey of flavonoids revealing the AChE inhibiting potential [46], easy availability, and most structures available in the PubChem database served as the rationale for the flavonoid selection for this study. Flavonoids' chemical structures were obtained from PubMed (https://www.ncbi.nlm.nih.gov/pubmed (accessed on 15 June 2022)), PubMed Central (https://www.ncbi.nlm.nih.gov/pmc (accessed on 15 June 2022)), and PubChem (https://www.ncbi.nlm.nih.gov/pccompound (accessed on 15 June 2022)) as shown in Figure S1. The structures were saved as (.mol) files. Atoms of hydrogen were incorporated into the structural optimization conducted with MOE software. We could use these characteristics to reduce the energy of specific molecules (gradient: 0.05, force field: the Merck molecular force field 94× (MMFF94×), chiral constraint, and current geometry). The revised flavonoid structures were then incorporated into the molecular docking database (.mdb) for docking analysis.

Docking Analysis
Docking analysis was performed using the MOE software. A docking study was conducted using the best available AChE1 and TIM structures via the SWISS-MODEL tool from the PDB. Energy minimization and 3D protonation were performed on the protein structure using the following settings (force field: MMFF94× solvation, chiral constraint, current geometry, gradient: 0.05). These simplified structures were then used as receptors in docking protocols. MOE's site discovery module was crucial in locating a protein's active site. The docking process was carried out using MOE's default configurations. The docking score of the MOE software was used to figure out the binding free energies of each ligand structure in a given orientation.

Ligand Interaction Calculations
The ligand interactions were calculated using the MOE Program's Ligand Interaction module (https://www.chemcomp.com/Products.htm (accessed on 25 June 2022)). It combines predicted distances between interacting atoms from 2D and 3D ligand-receptor protein models.

ADMET and Drug-Likeness Evaluation
To ensure the safety and efficacy of drug candidates, it is crucial that they exhibit desirable ADME properties and are non-toxic. Consequently, an assessment was conducted on the ADME profile of the synthesized compounds, which involved the evaluation of drug-likeness, partition coefficient, solubility, and various other parameters using the Swis-sADME module [26]. This module is available on the SIB (Swiss Institute of Bioinformatics) webserver at https://www.sib.swiss (accessed on 2 July 2022).

DFT Analysis
DFT analysis of compounds was performed using Gaussian 09 software and visualized through Gauss view 5.0. The structural coordinates of the lead compounds were optimized using B3LYP/6-31 G (d,p) level basis sets without any symmetrical constraints. Frontier molecular orbitals analysis was applied for selected flavonoid molecules.

Conclusions
Our research elucidated the structure and interactions of cattle tick proteins with plant flavonoids. In the last decade, flavonoids have been the subject of much research. Their potential health benefits have also been explained beyond their nutritional value. Newly discovered 3D structures of RmTIM and RmAChE1 have acceptable stereochemical and energetic qualities, allowing for the prediction of novel inhibitors. This research aims to evaluate the efficacy of flavonoids from medicinal plants as acaricides against proteins produced by the parasite R. microplus. The seven flavonoids (methylenebisphloridzin, thearubigin, fortunellin, quercetagetin-7-O-(6-O-caffeoyl-β-D-glucopyranoside) 1, quercetagetin-7-O-(6-O-p-coumaroyl-β-glucopyranoside), rutin, and kaempferol 3-neohesperidoside) were docked against RmAChE1, proving their strong interaction. This work will help make and evaluate flavonoids as tick acaricides to prevent tick-borne illnesses.  Data Availability Statement: The data supporting this study's findings are available from the corresponding authors upon reasonable request.