In Silico Exploration and Experimental Validation of Camellia sinensis Extract against Rhipicephalus microplus and Sarcoptes scabiei: An Integrated Approach

Sarcoptes scabiei is an ectoparasite of humans and animals that causes scabies. The Rhipicephalus (Boophilus) microplus is a blood-sucking ectoparasite that transmits various pathogens. These two parasites have caused great losses to a country’s dairy and agriculture sectors. The aim of this study was to determine the in vitro and in silico efficacy of Camellia sinensis plant extracts. Different concentrations of C. sinensis ethanolic plant extracts were prepared using the maceration method and were used against mites and ticks (in adult immersion test AIT and larval packet test LPT) to evaluate their in vitro acaricidal activity. Additionally, in silico molecular docking was performed to investigate the inhibitory interactions between the phytochemicals of the plant and S. scabiei and R. microplus glutathione transferase proteins (SsGST and RmGST). This study observed that the plant extract showed high efficacy in vitro against mites and different tick stages in adult immersion and larval packet tests. Additionally, the in silico study revealed a strong binding interaction between ellagic acid and SsGST protein, with a binding energy of −7.3 kcal/mol, with respect to permethrin (−6.7 kcal/mol), whereas quercetin and RmGST resulted in a docking score of −8.6 kcal/mol compared to deltamethrin (−8.2 kcal/mol). Overall, this study explored the potential of C. sinensis as a natural alternative for controlling tick and mite infestations and provided insights into the inhibitory mechanisms of its phytochemicals.


Introduction
Sarcoptic mange, known as scabies in humans, is a highly contagious skin disease caused by Sarcoptes scabiei, a mite that burrows into the top layer of the skin (epidermis).These mites are astigmatids that belong to the Sarcoptinae family.They actively penetrate the outermost layer of the skin, known as the stratum corneum [1].The adult mites mate, and the females lay eggs on the skin.Once hatched, larvae create small burrows, referred to as molting pouches, where they undergo molting and develop into nymphs, eventually Life 2023, 13, 2040 2 of 23 becoming adults [2].These parasites are found worldwide and can infect more than 150 different host species.Interestingly, they exhibit a surprising ability to be transmitted between different hosts, showcasing their epidemiological flexibility [3].In animals, the disease can manifest as a mild infection.The common symptoms include itchy papules, redness (erythema), scales, and hair loss (alopecia).In chronic cases, hyperkeratosis (thickening of the skin) and/or the formation of crusts with discharge may occur [4,5].Specific identification of this species depends on the host it infests.For example, in humans, it is S. scabiei var.hominis, whereas, in rabbits, it is referred to as S. scabiei var.cuniculi [6].
Sarcoptic mange is a highly contagious skin disease that spreads through direct skin-toskin contact, contact with contaminated objects, and exposure to an infected environment inhabited by severely affected hosts [7,8].Nymphs and female mites have a longer off-host survival period of up to 21 days than larvae and males, indicating the need for environmental application of biocides and repellents [9,10].S. scabiei infection causes significant morbidity and mortality in wild and domestic mammals, leading to potential economic losses.Primary morbidity is associated with secondary bacterial infections caused by Streptococcus pyogenes and Staphylococcus aureus [11].Rabbits infected with S. scabiei experience weight loss, reduced productivity, and compromised wool-fiber quality [12,13].Additionally, affected rabbits may develop dermatitis, pyoderma, eczema, and urticaria [12].Ticks are of great economic significance as pests in the global livestock industry, affecting cattle and other domestic species [14].According to FAO [15], over 80% of the global cattle population is infected with ticks.Among these pests, Rhipicephalus (Boophilus) microplus is of major concern in tropical and subtropical regions.The impact of R. microplus infestation includes reduced milk production and weight gain, increased mortality rates, hide damage, disease morbidity, control costs, and transmission of tick-borne pathogens, such as Babesia bigemina, B. bovis, and Anaplasma marginale.The global financial losses caused by tick infestation are estimated to be approximately USD 14,000-18,000 million.Moreover, the livestock industry in India and Pakistan spends approximately USD 498.7 million annually on tick and TBD (tick-borne disease) control [16].
Currently, scabies and R. microplus are primarily managed by chemical acaricides.These acaricides consist of active components categorized as macrocyclic lactones, organophosphates, formamidines, synthetic pyrethroids, phenyl pyrazoles, and growth inhibitors [17].Unfortunately, excessive, and sometimes inappropriate, use of these acaricides against ticks has led to the emergence of resistant tick populations [18,19].Developing new synthetic compounds is a costly and time-consuming endeavor, underscoring the need for plant-based alternatives to effectively control tick and mite infestation [20][21][22].
Green tea, scientifically known as Camellia sinensis, holds a prominent position among beverages worldwide and has deep cultural significance in China and Japan.The numerous health benefits of C. sinensis are primarily associated with its polyphenol content.These polyphenols consist mainly of catechins and derivatives such as (−)-epigallocatechin-3gallate (EGCG), (−)-epicatechin, (−)-epigallocatechin, (−)-epicatechin gallate, and (−)gallocatechin gallate [23].Studies have shown that C. sinensis exhibits anticancer, antitrypanosomal, and anti-Plasmodium properties [24,25].Despite extensive research on the pharmacological aspects of C. sinensis, no studies have been conducted to explore its acaricidal activity.
Glutathione transferase (GST), classified under the EC 2.5.1.18enzyme superfamily, plays a pivotal role in cellular detoxification processes.These enzymes facilitate the conjugation of reduced glutathione (GSH) with a wide range of endogenous and exogenous electrophilic compounds.By doing so, they protect cells from oxidative damage, contributing to overall cellular health.
In ticks, GST enzymes are believed to support tick survival by neutralizing toxins and aiding in the detoxification of various substances.This detoxification capacity is particularly relevant as overexpression of GST has been associated with drug resistance [26,27].Several studies have observed an upregulation in the transcription of GST genes and an increase in GST enzyme activity in ticks when exposed to both endogenous and exogenous compounds [28,29].
Given these critical biological functions, GSTs emerge as promising targets for the development of novel acaricidal chemotherapeutic drugs.By selectively inhibiting GST activity, it is possible to disrupt detoxification mechanisms in ticks and mites, potentially rendering them more susceptible to acaricides.This approach holds promise in addressing challenges associated with mite and tick infestations and drug resistance, ultimately contributing to more effective mite and tick control strategies.
The pivotal role of GST enzymes in detoxification processes and their association with drug resistance make them attractive candidates for targeted acaricidal drug development, offering new avenues for parasitic Acari control and management.
The current study aimed to determine the in vitro efficacy of C. sinensis plant extract against ticks and mites.Natural products are a mix of different phytochemicals.Testing the relative contribution of individual compounds is difficult and inefficient given complex mixtures.We selected the most abundant compounds (Table 1) for the in silico evaluation of inhibitory interactions with tick and mite proteins.This provides a foundation for future bioassays to evaluate the insecticidal role of individual compounds.

Plant Extract Preparation
C. sinensis aerial materials were collected from the Jalala region in Mardan, Khyber Pakhtunkhwa (coordinates: 34.3345 • N, 71.9075 • E).The collected plant materials were inspected for any physical damage and then rinsed in running tap water to remove surface debris.The leaves were then submitted to the herbarium of the Department of Botany, Abdul Wali Khan University Mardan (AWKUM) and were issued the accession number Awkum.Bot.917.They were then kept in a shady place at room temperature, away from sunlight, for air drying.
After 15 d of air drying, the plant material was ground to form a coarse powder using a plant grinder (YUEYUEHONG Model: HC-3000A Zhejiang, China).Powdered material (100 g) of each plant was then soaked in 1000 mL of 80% ethanol and agitated for 48 h using an orbital shaker (labForce Model 1165U07, Thomas Scientific, Swedesboro, NJ, USA) at 300RPM.The agitated solution was then reduced to a concentrated solution using a rotary evaporator (BUCHI Rotavapor Model: R-300 Flawil, Switzerland) at 40 • C under vacuum to remove ethanol.The solution was further concentrated in a water bath (Model WTB15, Memmert GmbH & Co. KG, Schwabach, Germany) at 45 • C until a very high concentration of the extract less than 10% of the original volume remained.The concentrated extract was used as a stock solution at different concentrations.

Mite Collection and Identification
Mites were collected from rabbits kept on rabbit farms at the Abdul Wali Khan University Mardan (AWKUM).Hay was provided as bedding to the rabbits on the farm, which was changed daily.The rabbits were first diagnosed for any signs of mange, and upon confirmation, skin scraps were taken using the procedure described by [13], and the rabbits were treated immediately.The skin scraps were then made into slides and observed under a microscope for S. scabiei.

Tick Collection and Incubation
Engorged adult female ticks were collected from the ground and cattle bodies in different locations near AWKUM in Mardan, Khyber Pakhtunkhwa, Pakistan.Standard tick identification keys were utilized to morphologically identify these ticks as R. (B.) microplus under a stereo zoom microscope [31].The engorged ticks were then transported to the parasitology lab, AWKUM, and incubated to lay the eggs.The eggs were then incubated to hatch into larvae for larval assays in an incubator (BIOBASE Model: BJPX-H50IV Shandong, China).Adult ticks and larvae were then used in AIT and LPT.

Contacts Bioassay for Mites
Scabies-infested rabbits were used for isolating mites in the current study.Skin scraps were taken from the rabbit adhering to the established animal-care guidelines [32].The infested skin was briefly cleaned and then scraped into a micro-Petri plate using a sterile surgical blade until the skin appeared red.The Petri plates containing skin scraps were incubated at 37 • C for 30 min to allow mites to emerge from the skin scraps.Ten mites were introduced into a Petri plate using a fine needle, and then 0.5 mL of samples of the plant extract were directly added onto the mites in the Petri plate.The procedure was performed in triplicate for each extract concentration.

Adult Immersion Test (AIT)
Completely engorged adult female ticks were used in the present study.The adult ticks were first rinsed in distilled water to remove any debris from their skin, and their weights were recorded.Ten (10) ticks were then dipped for 2-3 min at a concentration of the extract and then incubated for two days.This procedure was repeated for each extract concentration.Three replicates were performed for the study on separate days using a fresh new concentration of the extract for each replicate.Data were recorded for the total weight of eggs laid and tick mortality.The calculation of the percentage inhibition of oviposition (% IO) was performed using the following formula:

Percent inhibition of oviposition (%IO) =
Egg laying Index(control) − Egg laying Index(treated) Egg laying Index control × 100% The egg-laying index is calculated by dividing the mean weight of eggs laid by the mean weight of engorged females.

Larval Packet Test (LPT)
Newly hatched larvae were used in this experiment.Whatman no. 1 filter paper was immersed in a 0.6 mL extract concentration and then dried in an incubator at 37 • C. The filter paper was then made into a pouch by bending it and taping its sides with adhesive tape, with the upper side left open.Then, 100 larvae were carefully placed in the pouch, making a packet of larvae.The packet was then incubated at 28 ± 1 • C in a bio-oxygen demand incubator (BOD incubator) (BIOBASE Model: BJPX-B100 Shandong, China).Larvae were inspected for mortality at 24 and 48 h.Larvae that did not respond to light or had no appendage movement when teased with a needle were considered dead.This assay was also performed in triplicates.

Selection of Phytochemicals
Following a thorough examination of the current literature using Google Scholar (https://scholar.google.com/,accessed on 25 June 2023, Keywords; "Camellia sinensis" AND "phytochemistry" OR phytochemicals AND "Pakistan" OR "India") and PubMed (https://pubmed.ncbi.nlm.nih.gov/advanced/,accessed on 25 June 2023, Keyword; (Camellia sinensis [Title/Abstract]) AND (Pakistan [Title/Abstract]) AND (Phytochemistry) OR (phytochemicals) OR (GC-MS)), the C. sinensis phytochemicals were chosen.Articles that studied the phytochemical content of C. sinensis from Pakistan, India, or other related climates, utilizing GC-MS, TLC, and HPLC analytical methodologies were specifically chosen [33][34][35][36].Ten phytochemicals were chosen from extracts of C. sinensis (L) because of their ability to inhibit the R. microplus and S. scabiei glutathione transferase (GST) protein.These phytochemicals, which possess anti-inflammatory, antiseptic, and astringent properties, were identified through an extensive review of the existing literature.The corresponding structures of the ligands targeting GST proteins were acquired from the PubChem-NCBI database in the structure data format (SDF).Subsequently, these structures were converted into the protein databank (PDB) format using PyMOL for further analysis.

Preparation of Protein by Homology Modeling
The RCSB's (Research Collaboratory for Structural Bioinformatics) protein database (PDB) does not contain the available three-dimensional (3D) structures of the drug targets selected for this study.Therefore, homology modeling was used to obtain the 3D structures.The primary structures of S. scabiei glutathione transferase (Uniprot accession no: Q8I9R9) and R. microplus glutathione transferase (Uniprot accession no: O97117) were obtained in FASTA format from the UniProt Knowledgebase (UniProt KB database).To predict the 3D structure of SsGST and RmGST proteins, the SWISS-MODEL protein modeling server was employed [37].This server generates homology models by performing a targettemplate sequence alignment using the BLASTp and HHBlits programs and searching through template structures in the Protein Data Bank (PDB) [38] and SWISS-MODEL Template Library (SMTL) repositories.The top-ranked alignments of the templates are compared using the global model quality estimate (GMQE) and quaternary structure quality estimate (QSQE) to generate sets of descriptive three-dimensional structures, sequence dissimilarity, and quaternary protein structure information.The QMEAN values were utilized to predict 3D protein structures, taking into account modeling errors and quality estimation.These predicted protein structures are then assessed for stereochemical quality using the SAVES v6.0 server (https://saves.mbi.ucla.edu/,accessed on 27 June 2023).Additionally, the Ramachandran plot, which plots Ψ versus Φ conformational angles of the 3D macromolecule, is used to measure the torsion angles of Cα (ideal) -N-Cβ (obs).To predict the active site amino acids, the CastP calculation server (http://sts.bioe.uic.edu/castp/calculation.html,accessed on 28 June 2023) is employed, which calculates the delineating surface area and surface volume of the 3D protein structure [39].

Preparation of Modeled Proteins and Ligands for Docking
To facilitate docking, homology-modelled SsGST and RmGST proteins were prepared.Gasteiger charges and polar hydrogens were added, and nonpolar and polar hydrogen atoms were merged by introducing partial charges using AutoDock Vina 4.2 [40].To obtain structural information about bioactive compounds isolated from the leaves, fruit, and bark of C. sinensis, the relevant literature by [41] was reviewed.The 3D structures of these compounds were retrieved from the PubChem website (https://pubchem.ncbi.nlm.nih.gov/, accessed on 30 June 2023).Prior to docking, rotatable bonds were determined and nonpolar hydrogens were combined with polar hydrogen atoms for both phytochemical structures and conventional drugs for comparison purposes.
The docking process was conducted using AutoDock Vina software (version 1.1.2).Grid dimensions of 40 × 40 × 40 Å (grid size) with a grid-point spacing of 1.000 Å were employed for both proteins.The X, Y, and Z coordinates (grid centers) varied based on the specific receptor.The exhaustiveness parameter was set to a default value of eight (8) for all docking runs.The binding energy/affinity between the ligand and protein was calculated using the search algorithm within the AutoDock Vina software package.Following the completion of the docking runs, multiple binding modes representing different conformations of the ligands were obtained, along with their respective binding energy/affinity values.The most stable conformation, characterized by the lowest binding energy/affinity, was selected as the pose and utilized for postdocking analysis using BIOVIA Discovery Studio 2021.

Molecular Dynamics Simulation Analysis
Investigating the dynamic motion of atoms is crucial to understanding the stability and functionality of protein complexes.Molecular dynamics simulation plays a significant role in this regard [42].To perform dynamic simulations of the docking complex, we utilized the iMODS server, which is specifically designed for this purpose [43].The iMOD server (iMODS) (http://imods.chaconlab.org,accessed on 3 July 2023) was employed for conducting molecular dynamics (MD) simulations of the protein-ligand complexes.These simulations allowed us to assess the stability and molecular motion of docked complexes.
Utilizing the iMOD server (iMODS), we performed molecular dynamics simulations to analyze the structural dynamics of the docking complexes and determine molecular motion.Various parameters, such as deformability, B-factor, eigenvalues, variance, covariance map, and elastic network, were employed to evaluate the stability of the two protein-ligand complexes.The input files used for the simulations were docked PDB files that were uploaded to the iMODS server with default parameter settings.

Statistical Analysis
All statistical analyses were performed using R (version 4.3) running in RStudio (version 2023.06.1).The data were first arranged in Microsoft Excel (v.2302) and imported into the R working environment for further statistical analysis.Descriptive statistics of the data were calculated and are presented as mean ± standard deviation.The significance difference between the different concentrations was calculated using a one-way analysis of variance (ANOVA) followed by the Tukey honestly significance difference (HSD) test.Furthermore, 50% and 90% lethal concentrations and lethal times (LC and LT) were calculated in RStudio using the ecotox package and all the data were graphically presented using the ggplot2 and ggpubr R packages.

Mites Contact Bioassay
The efficacy of the C. sinensis leaf extract against S. scabiei var.cuniculi mites was evaluated in vitro by measuring mean mortality rates at various concentrations and time intervals.Permethrin and distilled water were included as control groups for comparison (Table 2).   3 presents the results of the lethal concentration (LC) calculations.As the time interval increased, the LC 50 and LC 90 values decreased, indicating a higher potency of the extract.At 6 h, the LC 50 and LC 90 values were 0.247 g/mL and 0.34 g/mL, respectively (Table 3, Figure 1A).Similarly, as the concentration increased, the LT 50 and LT 90 values decreased, indicating a faster onset of mortality (Table 4 and Figure 1B).LC: lethal concentration, LCL: lower confidence limit, UCL: upper confidence limit, SE: standard error.
Life 2023, 13, x FOR PEER REVIEW 7 of 25 3, Figure 1A).Similarly, as the concentration increased, the LT50 and LT90 values decreased, indicating a faster onset of mortality (Table 4 and Figure 1B).
The results demonstrated the concentration-and time-dependent effects of C. sinensis leaf extract on S. scabiei mites.Higher concentrations and longer exposure times were associated with reduced LC50, LC90, LT50, and LT90 values, indicating increased efficacy in inducing mortality.The slope values reflect the steepness of the dose-response or dosetime curves, while the intercept values represent the baseline response.
The results demonstrated the concentration-and time-dependent effects of C. sinensis leaf extract on S. scabiei mites.Higher concentrations and longer exposure times were associated with reduced LC 50 , LC 90 , LT 50 , and LT 90 values, indicating increased efficacy in inducing mortality.The slope values reflect the steepness of the dose-response or dose-time curves, while the intercept values represent the baseline response.

Adult Immersion Test
The percentage of oviposition inhibition (% IO) of the ticks was calculated to assess the efficacy of different concentrations of extracts derived from C. sinensis in controlling tick populations.The % IO values provided insight into the ability of the extract to inhibit tick egg hatchability at varying concentrations and time intervals.The highest % IO was observed with the highest concentration of C. sinensis extract (40 mg/mL), showing a remarkable inhibition rate of 46.071 ± 7.797%.As the concentration decreased, the % IO values also declined, with 20 mg/mL, 10 mg/mL, 5 mg/mL, and 2.5 mg/mL resulting in % IO values of 36.44 ± 4.99%, 25.32 ± 2.995%, 15.468 ± 6.065%, and 4.825 ± 2.828%, respectively, as shown in Table 5.These findings highlight the potential of the C. sinensis extract to inhibit egg-hatching ability.Higher concentrations of the extract resulted in an increased % IO, indicating greater efficacy in controlling tick populations (Table 5).

Larval Packet Test
The effects of different concentrations of C. sinensis extract on tick mortality were evaluated at 24 h and 48 h intervals.The control group consisted of ticks treated with deltamethrin and distilled water.The mean values and standard deviations (SD) of the measurements are listed in Table 5.The LC 50 -LC 90 values and LT 50 LT 90 , along with their respective 95% confidence limits (CL), were calculated for different time intervals (24 and 48 h), and various concentrations (2.5, 5, 10, 20, and 40 mg/mL) and are presented in Tables 6 and 7.At the 48 h mark, the LC50 value was 2.906 mg/mL (95% CL:2.300-3.505),while the LC90 value was 36.725mg/mL (95% CL:28.999-49.837), as shown in Table 6.LT50 values were calculated for different concentrations of the extract.At 2.5, 5, 10, 20, and 40 mg/mL, the LT50 values were 58.590, 35.788, 30.406, 28.078, and 25.913 h, respectively (Table 7, Figure 2A).The larval packet test results showed that with an increase in concentration or time, the LT and LC values decreased (Figure 2B).LT: lethal time, LCL: lower confidence limit, UCL: upper confidence limit, SE: standard error.

Homology Modeling
To analyze the structural arrangement of both proteins, 3D models were generated using the SWISS-Model online server.The SWISS-Model employs BLAST and HHblits to align the target sequences with previously characterized sequences and identify the most suitable template(s).The selected templates for the SsGST and RmGST proteins were 4q5q.1.A and 6gsv.1.A, respectively, with maximum sequence identities of 62.67% and 53.00% and coverages of 99% and 92%, respectively.The predicted model range was 219 and 223 amino acids for SsGST and RmGST, respectively.The root-mean-square deviation (RMSD) values for SsGST and RmGST, with their respective templates, were 0.067 and 0.084, respectively, indicating reasonable structural similarity.
The stereochemistry of both models was assessed using Ramachandran plots, which categorized amino acids into the core, additionally allowed, generously allowed, and disallowed regions.The Ramachandran plot for SsGST showed 92.2% amino acids in the core region (Figure 3A), while the plot for RmGST demonstrated 95.4% amino acids in the core region and one amino acid in the generously allowed region (Figure 4A).No amino acids appeared in the disallowed region of either protein.The quality of the overall protein structure was evaluated using z plots and ProSA.The z scores for both predicted models were −8.36 (Figure 3C) and −8.89 (Figure 4C), respectively, indicating acceptable overall quality comparable to NMR protein structures.
To further validate the 3D models, the ERRAT online server was used to validate.The overall quality factors calculated using ERRAT were 96.682 for SsGST and 95.814 for RmGST (Figures 3B and 4B, respectively), confirming the validity of the predicted models.These results confirm the reliability of the predicted models for the SsGST and RmGST proteins.

Molecular Docking
In this study, all ten phytochemicals exhibited significant inhibitory potential against the target proteins.Notably, ellagic acid, epigallocatechin gallate, kaempferol, and quercetin demonstrated the most promising results, targeting SsGST and RmGST proteins, respectively, with better binding scores than their respective controls (Table 8).Molecular docking analysis revealed that ellagic acid exhibited a strong binding affinity towards the SsGST protein, with a binding energy of −7.3 kcal/mol, surpassing permethrin (−6.7 kcal/mol), as shown in Figure 6A,C,E and Table 8.Ellagic acid formed hydrogen bonds with Thr-100 (three bonds), Trp-96 (one bond), and Asp-89 (one bond) residues of the SsGST protein.Furthermore, hydrophobic interactions were observed between ellagic acid and Trp-96 and Arg-80, indicating their potential as inhibitors of S. scabiei glutathione transferase.Among other ligands, epicatechin gallate also displayed a favorable binding energy (−6.8 kcal/mol) and amino acid interactions compared to permethrin (−6.7 kcal/mol).

Molecular Dynamics Simulation Analysis
To assess the stability and dynamics of the docked complexes, molecular dynamics (MD) simulations were conducted using the iMOD server.Normal mode analysis (NMA) was employed to examine the slow dynamics and large-scale conformational fluctuations of the docked complexes, namely ellagic acid-SsGST and quercetin-RmGST, as shown in Figures 7 and 8, respectively.
Deformability and B-factor profiles provide valuable information about the mobility and flexibility of docked proteins.The peaks observed in these profiles indicated regions of higher deformability, suggesting greater flexibility or mobility within these regions.The highest peaks in the profiles represent the regions with the most pronounced flexibility.Figures 7 and 8 provide illustrations of the deformability and B-factor profiles for the ellagic acid-SsGST and quercetin-RmGST complexes, respectively.These profiles enable a comparison between the results obtained from normal mode analysis (NMA) and structures obtained from the Protein Data Bank (PDB).By examining these profiles, the regions in proteins that exhibit significant flexibility or undergo conformational changes can be In contrast, quercetin demonstrated a strong binding affinity for the RmGST target protein (Table 8, Figure 6D,F), with a docking score of −8.

Molecular Dynamics Simulation Analysis
To assess the stability and dynamics of the docked complexes, molecular dynamics (MD) simulations were conducted using the iMOD server.Normal mode analysis (NMA) was employed to examine the slow dynamics and large-scale conformational fluctuations of the docked complexes, namely ellagic acid-SsGST and quercetin-RmGST, as shown in Figures 7 and 8, respectively.

Discussion
The use of plant extracts for pest and disease control has gained attention because of their natural degradation properties [44].Despite their advantages, biopesticides represent a small portion of the pesticide market [45].However, the biopesticide sector has experienced significant growth in recent years, with an annual growth rate predicted to surpass that of chemical pesticides [45,46].
This study focused on investigating the acaricidal efficacy of C. sinensis extract as a potential source for developing herbal acaricides and identifying bioactive compounds Deformability and B-factor profiles provide valuable information about the mobility and flexibility of docked proteins.The peaks observed in these profiles indicated regions of higher deformability, suggesting greater flexibility or mobility within these regions.The highest peaks in the profiles represent the regions with the most pronounced flexibility.Figures 7 and 8 provide illustrations of the deformability and B-factor profiles for the ellagic acid-SsGST and quercetin-RmGST complexes, respectively.These profiles enable a comparison between results obtained from normal mode analysis (NMA) and structures obtained from the Protein Data Bank (PDB).By examining these profiles, the regions in proteins that exhibit significant flexibility or undergo conformational changes can be identified.This information can be valuable for understanding the dynamic behavior of docked complexes and their potential impact on protein function and stability.
The eigenvalue and variance graphs reveal the characteristics of each normal mode with an inverse relationship between the eigenvalues and variances.The eigenvalue and variance graphs for the ellagic acid-SsGST and quercetin-RmGST complexes are presented in Figures 7E and 8E, respectively.Additionally, a covariance matrix was provided for the docked complexes, representing anticorrelated motion (depicted in blue), uncorrelated motion (depicted in white), and correlated motion (depicted in red).This matrix provides information about the atomic interactions and their dynamics.Atomic connections were visualized using elastic springs in a string-model representation (Figures 7F and 8F), along with a plot matrix where the connections are depicted as grey dots (Figures 7G and 8G).

Discussion
The use of plant extracts for pest and disease control has gained attention because of their natural degradation properties [44].Despite their advantages, biopesticides represent a small portion of the pesticide market [45].However, the biopesticide sector has experienced significant growth in recent years, with an annual growth rate predicted to surpass that of chemical pesticides [45,46].
This study focused on investigating the acaricidal efficacy of C. sinensis extract as a potential source for developing herbal acaricides and identifying bioactive compounds against ticks and mites.Ticks and mites were treated with different concentrations of the ethanol extracts.Previous studies have explored the acaricidal properties of various herbs against mites and ticks [13,30,47,48].C. sinensis has been reported to possess antioxidant [49], antibacterial [50], anti-inflammatory, and antihistaminic properties [51], as well as insecticidal properties [52], whereas no acaricidal potential of the plant has been documented.The plant extract has resulted in greater mortality in mites and comparable mortality in ticks as compared to the positive control permethrin and deltamethrin.The decrease in the chemical acaricides toxicity can be due to the acaricidal resistance by these mites and ticks.In a study by Gu, et al. [53], the Ailanthus altissima bark extract has resulted in similar results with the extract having more significant acaricidal potential against S. scabiei and Psoroptes scabiei compared to the available chemical acaricide fenvalerate.The findings of this study are consistent with previous research of Seddiek, et al. [54].Moreover, Dodonaea angustifolia, Eucalyptus globulus, Millettia ferruginea, and Euphorbia abyssinica plant extracts have been found to have acaricidal potential against these mites [52].
Molecular docking studies have been widely utilized to predict ligand-target interactions and gain insights into the biological activity of natural products.These studies also provide clues regarding the mechanisms of action and binding modes within the active sites of enzymes [55].In this study, ten representative compounds from C. sinensis were selected for docking analyses against two target proteins: S. scabiei glutathione transferase (UniProt accession no: Q8I9R9) and R. microplus glutathione transferases (UniProt accession no: O97117).
Docking analysis with S scabiei glutathione transferase revealed that, among the ten compounds, ellagic acid exhibited strong interactions with several amino acid residues through hydrogen bonds (Thr-100, Trp-96, Asp-89) and hydrophobic interactions (Trp-96 and Arg-80), with a docking score of −7.3 kcal/mol.Epicatechingallate, epigallocatechingallate, epicatechin, quercetin, caffeoylquinicacid, kaempferol, catechin, gallicacid, and theanine also displayed docking scores, suggesting their potential involvement in the antimite activity of C. sinensis through interactions with the target protein.
In the antitick docking study, the ten compounds were docked with R. microplus glutathione transferases, and they exhibited docking scores ranging from −5.2 to −8.82 kcal/mol.Quercetin displayed the highest score against RmGST protein, followed by kaempferol, catechin, epi-

Figure 1 .
Figure 1.(A,B) represents the lethal concentrations and lethal time (LC50, LC90 and LT50, LT90), whereas the error bars represent the lower and upper confidence limits as 90% confidence intervals.

Figure 1 .
Figure 1.(A,B) represents the lethal concentrations and lethal time (LC 50 , LC 90 and LT 50 , LT 90 ), whereas the error bars represent the lower and upper confidence limits as 90% confidence intervals.

Life 2023 , 25 Figure 2 .
Figure 2. (A,B) represents the lethal concentrations and lethal time (LC50, LC90, and LT50, LT90), whereas the error bars represent the lower confidence limit and upper confidence limit at 90% confidence intervals.

Figure 2 .
Figure 2. (A,B) represents the lethal concentrations and lethal time (LC 50, LC 90 , and LT 50, LT 90 ), whereas the error bars represent the lower confidence limit and upper confidence limit at 90% confidence intervals.

Figure 3 .
Figure 3. Validation plots and scores for the TrRosetta server's predicted 3D structure of Sarcoptes scabiei glutathione transferase (SsGST) showing (A) the Ramachandran plot where the red, yellow, and black colors represent the most favorable, favorable, and disallowed regions, respectively; Phi and Psi bonds represent torsion angles that predict the possible conformation of the peptides; (B) represents the ERRAT's overall quality factor values; and (C) represents the PROSA server's Z score values.

Figure 3 .
Figure 3. Validation plots and scores for the TrRosetta server's predicted 3D structure of Sarcoptes scabiei glutathione transferase (SsGST) showing (A) the Ramachandran plot where the red, yellow, and black colors represent the most favorable, favorable, and disallowed regions, respectively; Phi and Psi bonds represent torsion angles that predict the possible conformation of the peptides; (B) represents the ERRAT's overall quality factor values; and (C) represents the PROSA server's Z score values.

Figure 4 .
Figure 4. Validation plots and scores for the TrRosetta server's predicted 3D structure of Rhipicephalus microplus glutathione transferase (RmGST) showing (A) the Ramachandran plot where the red, yellow, and black colors represent the most favorable, favorable, and disallowed regions, respectively; Phi and Psi bonds represent torsion angles that predict the possible conformation of the peptides; (B) represents the ERRAT's overall quality factor value; and (C) represents the PROSA server's Z score values.

Figure 4 .
Figure 4. Validation plots and scores for the TrRosetta server's predicted 3D structure of Rhipicephalus microplus glutathione transferase (RmGST) showing (A) the Ramachandran plot where the red, yellow, and black colors represent the most favorable, favorable, and disallowed regions, respectively; Phi and Psi bonds represent torsion angles that predict the possible conformation of the peptides; (B) represents the ERRAT's overall quality factor value; and (C) represents the PROSA server's Z score values.
6 kcal/mol compared to deltamethrin (−8.2 kcal/mol).Quercetin formed a hydrogen bond with Asp-125 and established π-alkyl bonds with Arg-18, Val-28, Ala-196, and Ala-22, as shown in Figure 6B,D,F.In addition, a carbon-hydrogen bond was observed with Ser-200.These findings suggest that quercetin has the potential to inhibit RmGST and can serve as an antitick agent.

Figure 6 .
Figure 6.(A,C,E) show the complex three-dimensional structural interactions of the ellagic acid and SsGST and their corresponding 2D interactions; (B,D,F) show the highly complex 3D structural interactions of the quercitin and RmGST and their corresponding 2D interactions.

Figure 6 .
Figure 6.(A,C,E) show the complex three-dimensional structural interactions of the ellagic acid and SsGST and their corresponding 2D interactions; (B,D,F) show the highly complex 3D structural interactions of the quercitin and RmGST and their corresponding 2D interactions.

Life 2023 , 25 Figure 7 .
Figure 7.The iMOD server's normal mode analysis (NMA) outputted the following plots for the ellagic acid-SsGST complex.(A) represents the NMA mobility, (B) B-factor, (C) deformation plot, (D) variance plot, (E) eigenvalues, (F) covariance-matrix plot (the anticorrelated, uncorrelated, and correlated states of atomic motion are represented by the blue, white, and red hues, respectively), and (G) elastic network model, where grey colors represent atom connections.

Figure 7 .
Figure 7.The iMOD server's normal mode analysis (NMA) outputted the following plots for the ellagic acid-SsGST complex.(A) represents the NMA mobility, (B) B-factor, (C) deformation plot, (D) variance plot, (E) eigenvalues, (F) covariance-matrix plot (the anticorrelated, uncorrelated, and correlated states of atomic motion are represented by the blue, white, and red hues, respectively), and (G) elastic network model, where grey colors represent atom connections.

Figure 8 .
Figure 8.The iMOD server's normal mode analysis (NMA) outputted the following plots for the quercitin-RmGST complex.(A) represents the NMA mobility, (B) B-factor, (C) deformation plot, (D) variance plot, (E) eigenvalues, (F) covariance-matrix plot (the anticorrelated, uncorrelated, and correlated states of atomic motion are represented by the blue, white, and red hues, respectively), and (G) elastic network model, where grey colors represent atom connections.

Figure 8 .
Figure 8.The iMOD server's normal mode analysis (NMA) outputted the following plots for the quercitin-RmGST complex.(A) represents the NMA mobility, (B) B-factor, (C) deformation plot, (D) variance plot, (E) eigenvalues, (F) covariance-matrix plot (the anticorrelated, uncorrelated, and correlated states of atomic motion are represented by the blue, white, and red hues, respectively), and (G) elastic network model, where grey colors represent atom connections.

Table 1 .
Selected compounds from C. sinensis plant for its in silico screening against ticks' and mites' protein.

Table 2 .
The mean mortality ± standard deviation of the mortalities of S. scabiei var.cuniculi by Camellia sinensis leaf extract in vitro.
eMeans with no similar letters in the superscript in the same column are significantly different by Tukey's HSD test at the 5% level of significance (p < 0.05).

Table 3 .
The calculated lethal concentrations were responsible for 50% and 90% S. scabiei mortality at various time intervals for the Camellia sinensis leaf extract.

Table 2 .
The mean mortality ± standard deviation of the mortalities of S. scabiei var.cuniculi by Camellia sinensis leaf extract in vitro.

Table 3 .
The calculated lethal concentrations were responsible for 50% and 90% S. scabiei mortality at various time intervals for the Camellia sinensis leaf extract.

Table 4 .
The calculated lethal time was responsible for 50% and 90% S. scabiei mortality at varying concentration intervals for Camellia sinensis leaf extract.

Table 5 .
The % mean values ± standard deviation for different concentrations of Camellia sinensis leaf extract on R. (B.) microplus larval mortality at 24 and 48 h, as well as inhibition of oviposition in adult female R. (B.) microplus.
eMeans with no similar letters in superscript in the same column are significantly different by Tukey's HSD test at the 5% level of significance (p < 0.05).

Table 6 .
The calculated lethal concentrations were responsible for 50% and 90% mortality at various time intervals of Camellia sinensis exposure to R. (B.) microplus larvae.

Table 6 .
The calculated lethal concentrations were responsible for 50% and 90% mortality at various time intervals of Camellia sinensis exposure to R. (B.) microplus larvae.