In Silico and In Vivo Evaluation of Synthesized SCP-2 Inhibiting Compounds on Life Table Parameters of Helicoverpa armigera (Hübner)

Simple Summary Different crops and vegetables are attacked by a complex of pest insects causing severe losses in their yield and quality, among them Helicoverpa armigera is one of the key pests. Its polyphagous status and extended environmental tolerance, coupled with high fertility, fecundity rate and short generation time, enables it to quickly attain a primary pest status in any suitable host. To minimize losses caused by this notorious pest, multiple applications of chemical pesticides are required. Due to such indiscriminate use of synthetic insecticides, environmental and human health hazards have become a major concern of the modern-day plant protection industry. Novel scaffolds that may regulate insect growth can offer a sustainable alternative to conventional insecticides. In this study, synthesized small molecules, with a tendency to disrupt insect molting, were evaluated against a Helicoverpa armigera. One of the tested compounds significantly reduced larval and pupal weight accumulations and prolonged stadia lengths resulting in disrupted population growth. At the same time, the emerged females had significantly reduced fertility. These findings suggest that further optimization of tested scaffold may lead to help finding new insecticide-like molecules that will reduce the dependence on traditional chemical insecticides. Abstract For environment-friendly, safe and nonpersistent chemical control of a significant polyphagous insect pest, Helicoverpa armigera, discovery of growth-regulating xenobiotics can offer a sustainable alternative to conventional insecticides. For this purpose, chemically synthesized compounds to inhibit sterol carrier protein (SCP-2) function using in silico and in vivo assays were evaluated to estimate their impact on the survivals and lifetable indices of H. armigera. From nine chemically synthesized compounds, OA-02, OA-06 and OA-09 were selected for this study based on binding poses mimicking cholesterol, a natural substrate of sterol carrier protein and molecular dynamics simulations. In vivo bioassays revealed that all compounds significantly reduced the larval and pupal weight accumulations and stadia lengths. Subsequently, the pupal periods were prolonged upon treatment with higher doses of the selected compounds. Moreover, OA-09 significantly reduced pupation and adult emergence rates as well as the fertility of female moths; however, fecundity remained unaffected, in general. The life table parameters of H. armigera were significantly reduced when treated with OA-09 at higher doses. The population treated with 450 μM of OA-09 had the least net reproductive rates (Ro) and gross reproductive rate (GRR) compared to the control population. The same compound resulted in a declining survival during the early stages of development coupled with reduced larval and pupal durations, and fertility. These results have a significant implication for developing an effective and sustainable chemical treatment against H. armigera infestation.


Introduction
Insecticides have always been a main component of conventional and integrated pest management practices [1]. In the high-income countries, the application of pesticides has helped to increase crop development and yield significantly by reducing losses caused by noxious insect pests. However, intensification of agriculture with excessive use of pesticides has resulted in serious health and environment issues and among them environmental pollution and food contaminations are at the top [2]. During the 1950s, use of pesticides in agricultural landscapes increased exponentially due to plant diseases and cultivation of new land along with the need to boost yields to cater for the world's increasing dietary requirements [3,4]. The United States alone imports more than 350 MT of pesticides, whereas the import of chemical pesticides has touched the figure of 250 MT per annum in recent years [5,6]. The world population approaching 11 billion demands an increase in quality and quantity of available agricultural commodities, and if such pattern of indiscriminate use of insecticides persists, environmental impact will increase many folds in no time [7,8]. A multitude of different formulations are currently available and consumption of more than 65 MT of active pesticide ingredients in the world [9] strongly mobilize the agricultural communities and environmental activists to struggle for the use of safer alternatives.
Cotton is treated as one of the main cash crops [10] in subtropical and tropical parts of Asia, Africa and the Americas [11][12][13]. One of the major factors in crop yield decline can be attributed to insect losses [14,15]. During the whole cropping season, cotton crop is attacked by a complex of pest insects causing severe losses in yield and quality and Helicoverpa armigera (Hübner) (Lepidoptera; Noctuidae) is one of the key pests that require multiple applications of chemical pesticides [16]. Its polyphagous status [17] and extended environmental tolerance, coupled with high fertility and fecundity rates and short generation times, enables it to quickly attain a primary pest status in any suitable host [18,19]. Being a polyphagous pest, H. armigera is reported to cause severe losses to cotton, corn and soybean in US agroecosystems [20,21], tomato and tobacco in Brazil [22,23] and also soybean, sorghum, maize, vegetables and citrus in south American regions [17,24,25]. It has more than 180 wild and cultivated host species worldwide with at least forty cultivated families including Solanaceae, Poaceae, Malvaceae, Asteraceae and Fabaceae [26,27]. Since the larvae primarily feed on flowers and fruit, a direct yield loss is observed [28][29][30][31][32].
Farmers, in general, invest more to protect cash crops and the same is the case here in Pakistan, where according to an estimate, more than 80% of the pesticides that are imported to the country are applied in cotton agroecosystem [33]. Although this usage has increased cotton yields three-fold, it also poses a serious threat of environmental contamination and health hazards [34]. Due to excessive and indiscriminate use of chemical pesticides, H. armigera has developed strong resistance against organophosphates, organochlorines, carbamates and pyrethroids and even to Bt-gene Cry-1 in transgenic cultivars of cotton [35][36][37][38][39][40]. The farmers, therefore, have now started using parallel management practices for its control by managing crop rotation and cultivation adjustments [41], using resistant varieties, installing pheromone traps, shifting to genetically modified varieties, removal of alternate hosts and using biological control agents [24,27]. However, discovery of a novel mode of action synthetic molecules that may target insect physiological pathways is often neglected. Use of synthetic and target specific compounds such as the insect growth regulators (IGR) could be one of the alternative and safe approaches to develop new pesticides. Such compounds provide effective and environmentally safe insecticides that are selective and nonpersistent with the least health risks for humans and animals [42][43][44]. Already available IGRs in the market belong to the family of chitin synthesis inhibitors such as hexaflumuron and chlorfluazuron [42], and insects' juvenile hormone mimic, e.g., pyriproxyfen [45]. Further expansion in the field of identification of insect metabolic pathways will expand the options for pesticide discoveries. One such newly identified pathway, is the sterol production pathway [46]. Insects and other organisms use sterols as a structural and functional component of membrane rigidity and fluidity [16,47]. In this pathway, the sterols work as a precursor of molting hormones including 20-hydroxyecdysone, ecdysone and makisterone A. The most common form of a sterol is cholesterol that is often found in the mid-guts of the insects [28,48]. Like plants and mammals, insects are not able to produce sterol in the body, they fulfill the requirement of sterol from their food [49,50] and then use a sterol carrier protein-2 (SCP-2) to extract cholesterol. It is a nonspecific intracellular lipids carrier that is reported to play a vital role in the transport of cholesterol in humans, rats and insects [51]. The existence of sterol carrier proteins (SCP-2) has already been reported from major insect orders, including Diptera, Hymenoptera, Lepidoptera, Hemiptera and Coleoptera [52,53]. Blocking SCP-2 can result in delayed larval development, interrupted molting, deformities and reduced fecundity [54,55]. Moreover, some studies have reported the potential of SCP-2 inhibitor molecules to reduce the activities of insects' chemical detoxifying enzymes, i.e., GST and P450, etc. [53,56].

Toxicity Predictions
The toxicity risks of the OA series compounds were predicted by DataWarrior [60]. The predicted properties contained tumorigenicity, reproductive effect, irritant effect and mutagenicity.

Molecular Docking
The NMR structure of sterol carrier protein (PDB ID: 4UEI) was prepared by the protein preparation wizard in Schrodinger Maestro [61]. The protein preparation involved assigning bond orders, addition of hydrogens, creating disulfide bonds and zero-order bonds to metals, removing water molecules from hetero groups beyond 5 Å and leaving the hetero state in its default pH (7.0). The hydrogen bonds were optimized by using PROPKA at pH 7.0. Finally, restrained energy minimization was performed using the OPLS_2005 forcefield [62]. After preparation of the receptor, a 3D grid was generated at particular sites to conduct site-specific docking. The values of internal coordinates for X, Y and Z were 7.67, −8.81 and −11.16, respectively. Similarly, the compounds were prepared using the LigPrep tool [63]. Various conformations of rings and stereoisomers of ligands were generated. The OPLS_2005 forcefield was used to optimize and minimize the 3D conformers of ligands. The prepared ligands were docked to the receptor grid by using the Glide docking module in SP (Standard Precision) mode [64]. The docked ligands were analyzed based on their binding modes and glide score. The compounds with plausible binding modes were selected for MD simulation to estimate their binding stability.

Molecular Dynamics Simulation
Three complexes were selected based on the binding modes. The selected complexes were explored to check binding stability by a 50 ns long MD simulation using the VMD [65] and NAMD tools [66]. The initial files were prepared through AMBER21 tools [67]. Antechamber was used to generate the ligand topology files, while the missing hydrogen were added to protein by the LeaP program [67]. In order to mimic the physiological environment, a 10 Å periodic boundary solvation box of TIP3P water molecules [68] was added to complexes. The solvated systems were neutralized by Na + and Cl − counter ions. The systems were minimized for 20 fs to avoid energy clashes. The equilibration of solvation systems was conducted at 310 K. To maintain the system stability, three additional equilibrations were conducted at 200 K, 250 K and 300 K, respectively. The stabilized systems were subjected to 50 ns simulation in the production run. The MD trajectories were stored at every 2 ps and analyzed by CPPTRAJ [69] and Bio3D package [70] of R program. Moreover, the binding free energies of the complexes were calculated by the MM/GBSA module of AMBER21 tools.

Larval Diet and Test Concentration Preparation
The selected compounds were assayed for their biological activity against the 3rd, 4th and 5th instars larvae of H. armigera using a diet incorporation method [71]. Artificial larval diet was prepared using the methods and ingredients described by Hamed and Nadeem (2008). Three concentrations of the test compounds, i.e., 50 µM, 150 µM and 450 µM were mixed at the rate of 40 µL/2 g of prepared diet.

Collection and Rearing of Cultures
Freshly laid egg masses and young larvae of H. armigera were collected from the cotton crop cultivated in the research fields of the Cotton Research Institute (CRI), Central Cotton Research Institute (CCRI) and research fields of the Bahauddin Zakariya University, Multan, to start the laboratory cultures for bioassays. The field collected larvae were kept individually to confine any infections they might have and minimize cannibalism [72] and were reared on freshly collected cotton leaves until they pupated. From the next generation, the larvae were fed on an artificial diet for two generations to exclude any effect of field diets. Glass containers (500 mL) covered with muslin clothes were used to house the larvae and the cultures were kept at 25 ± 2 • C and 75 ± 5% r.h. with a 14:10 h light:dark regime. To facilitate pupation, after the 4th molt, the larvae were supplied with sterile soil at the base of each jar as a medium for pupation. Later, the pupae were retrieved from the soil media and were transferred to clear jars (5L) until adult emergence. The moths, when eclosed, were provided with 10% honey solution [73] and strips of viscose tissue, hanging vertically as oviposition surfaces. Such rearing on artificial diet was continued until the second filial generation and then the larvae were exposed to the test compounds in bioassays. The master cultures in the laboratory population were augmented regularly by introducing newly collected field populations, when available.

Life Table Bioassay
The life table parameters of H. armigera were studied using the test concentrations mixed with artificial diet, while the plain diet (with no tested compound mixed) was used as a control. The experiment was replicated three times with 20 larvae of each instar. After weighting and recording the size of each larva, 2 g diet/larva were provided in each Petri dish containing a single larva. To avoid any fungal infections in the larval colonies, the diet was replaced on alternate days until pupation. The larvae health and survival were observed daily, and the stadia length was recorded. Three days after pupation, the pupae were weighed and kept in sterilized jars until adult eclosion, and their pupal periods were recorded. The newly emerged moths were sex sorted based on color of the forewings (males have greenish while the females have brownish forewings) [74]. Five pairs of adults from each treatment were transferred into oviposition jars with a supply of 10% honey solution and vertically hanging viscose tissue strips, as described earlier. The data for fecundity and fertility were recorded daily, and the pre-oviposition periods and adult longevity was noticed at the end of oviposition and until the mortality of the female moths, respectively.
Age-stage, two sex life table analyses [75] using fecundity, pre-oviposition period, total pre-oviposition period, larvae, adult and eggs were conducted following the model described by [76,77]. Survival (l x ) and fecundity (m x ) rates were calculated from proportion of larvae survived and the number of eggs laid per female followed by calculation of population growth parameters, such as net reproductive rate (R o ), intrinsic rate of population increase (r), finite rate of population increase (λ) and generation time (T), using the equations given below: The population parameter intrinsic rate "r" was calculated to follow the iterative bisection methodology of Euler-Lotka formula: Mean generation time was calculated as: where "k" = number of stages, "S xj " = survival rate at specific stage (j) and age in days (x), "ê xj " = age-stage specific female fecundity, "l x " = age specific survival rate, "m x " = age specific fecundity, "R o " = net reproductive rate and "T" = mean generation time. The age index ranged from 0 to ∞ and finite rate "λ" was considered as e r . (Goodman, 1982). At the end, variance and standard error of fecundity, developmental time, adult longevity and other life table characteristics were calculated by bootstrapping [78] and the data were conformed to the assumption of equal variance. The means were analyzed with ANOVA followed by post hoc pairwise comparisons using Tukey's HSD tests [76].

Statistical Analyses
The data regarding developmental parameters (larval, pupal, and adult longevities, pupation, and adult emergence) and reproductive potentials (fecundity and fertility) were analyzed using a 2-way factorial ANOVA using the compounds and their concentrations as independent factors, followed by Tukey's HSD post-hoc pairwise comparisons. Since the data regarding larval mortality did not conform to the assumption of homogeneity even after appropriate data transformation (arcsine square root transformation) was applied, we analyzed the effects of compound type and their concentrations on larval mortality individually, using a nonparametric Kruskal-Wallis H test. The LC 50 and LC 90 values were calculated by fitting the probit model. All the above-mentioned analyses were conducted using SPSS software.

Toxicity Predictions
The toxicity properties including tumorigenicity, reproductive effect, irritant effect and mutagenicity were predicted by DataWarrior. As shown in Table 1, OA-01 showed high mutagenicity, irritant effect, and reproductive effect. Other compounds did not show high toxicity properties.

Molecular Docking
The compounds were docked to the binding sites of SCP-2 protein to find their binding modes. All the compounds were ranked based on the glide score ( Table 2)

MD Simulation
The binding modes of selected compounds were further validated by 50 ns MD simulation. The MD trajectories were analyzed to calculate the root mean-square deviation (RMSD), root mean-square fluctuations (RMSF) and solvent accessible surface Area (SASA).
The RMSD of backbones atoms of apo protein and its complexes was calculated to estimate the stability of the protein. It can be observed that all systems were equilibrated at 5 ns. After equilibration, the RMSD of apo protein (blue) remained in the range of~3.5-4 Å. A major deviation was observed during 30 to 40 ns where RMSD increased to~6 Å but it attained the previous conformation at 40 ns and then remained in the same range until the end of simulation (Figure 2A). The RMSD plots of SCP-2 complexes were compared to apo protein. The major deviations in the RMSD were observed in SCP-2-OA-02 (black) complex, which attained a RMSD value of~4.5-5 Å at 15 ns and remained in the same range until 45 ns. A slight drop was observed after 45 ns, but it gained the same range towards the end of simulation. The other two complexes, i.e., SCP-2-OA-06 and SCP-2-OA-09 showed more stable trajectories than the apo protein. The plot of SCP-2-OA-06 (red) showed that the RMSD remained in the range of~3.5-4 Å until 35 ns and then a deviation of~0.5 Å was observed from 35 to 45 ns, but it gained the same range towards the end of the simulation. The RMSD plot of SCP-2-OA-09 (green) showed a stable trend throughout the simulation. The stabilities of protein-ligand complexes were further validated by calculating the RMSD of ligand. It can be observed that the ligands showed stable plots with RMSD values ranging from~0.5 to 1.5 Å except for OA-02 that showed a minor deviation at 30 ns, but it again attained the previous confirmations ( Figure 2B). These plots showed the stability of protein-ligand complexes. ns. The SASA values remained in the range of ~7000 to 8000 Å2 until 30 ns and then minor deviations were observed in red and green plots, while black and blue remained stable. The overall SASA analysis showed that the protein structure remained stable in all systems and no diffusion was observed as the total accessible surface area to solvent was not increased. All these analyses demonstrated that the protein remained stable when bound to the OA-02, 06 and 09 compounds and no conformational changes were observed in protein structure. The binding free energies of these complexes were also calculated by applying the MM/GBSA to the last 300 frames of MD trajectories. The total binding energy in terms of solvation energy, gas phase energy and entropic contributions were calculated. To find out the contribution of binding site residues, the systems were decomposed. The contribution of energy components, i.e., VDWAALS, EEL, EGB, ESURF, Delta Ggas, Delta Gsolv and total energy of each complex is given in Supplementary Materials Table S1. Moreover, the contribution of binding site residues was also calculated by energy decomposition The dynamic behavior of the amino acid residues was determined by calculating the RMSF values. The higher RMSF values show the flexibility of residues or the loops of the protein, while the lower values describe the rigid part of protein, i.e., alpha helices and beta sheets. Higher RMSF values at N and C terminals can be observed as these are the loop regions ( Figure 2C). A minor fluctuation was observed at residues 23 to 26 and 45 to 50 as these reside in loops. The remaining residues maintained their low flexibility, which showed the stability of protein and its complexes.
The protein structure stability was further analyzed by calculating the solvent accessible surface area. The distortion in the protein structure increased the total exposed area of protein accessible to solvent. So, if the SASA increases from the initial conformation, it shows a diffused protein structure, while low values elaborate the structure stability. The SASA values of the apo SCP-2 and its complexes are shown in Figure 2D. The SASA values of all systems were started from~8500 Å2 and then dropped gradually to~8000 Å2 at 10 ns. The SASA values remained in the range of~7000 to 8000 Å2 until 30 ns and then minor deviations were observed in red and green plots, while black and blue remained stable. The overall SASA analysis showed that the protein structure remained stable in all systems and no diffusion was observed as the total accessible surface area to solvent was not increased. All these analyses demonstrated that the protein remained stable when bound to the OA-02, 06 and 09 compounds and no conformational changes were observed in protein structure.
The binding free energies of these complexes were also calculated by applying the MM/GBSA to the last 300 frames of MD trajectories. The total binding energy in terms of solvation energy, gas phase energy and entropic contributions were calculated. To find out the contribution of binding site residues, the systems were decomposed. The contribution of energy components, i.e., VDWAALS, EEL, EGB, ESURF, Delta G gas , Delta G solv and total energy of each complex is given in Supplementary Materials Table S1. Moreover, the contribution of binding site residues was also calculated by energy decomposition function of MM/GBSA, which showed that GLN107, PHE110 and MET121 had high contribution in total binding free energy (Supplementary Materials Table S2).

Larval Mortality
When the 3rd instar larvae were exposed to different concentrations of OA-02, OA-06 and OA-09, a significantly different larval mortality, both due to type of compounds and their concentrations was observed (Kruskal-Wallis H = 6.33; df = 2; p = 0.042 and Kruskal-Wallis H = 9.30; df = 2; p = 0.026, respectively). Highest mortality of 3rd instar larvae was recorded when treated with 50 µM concentration of OA-09 ( Figure 3A). Since the insects are capable of selective feeding [79], it can be assumed that avoiding excessive feeding on a diet with higher concentrations may have resulted in low mortalities as compared to those fed with a lower dosage. Whereas, the mortality of 4th and 5th instar larvae were not significantly affected by the type of compounds  Table 3).
Insects 2022, 13, x FOR PEER REVIEW 10 of 20 function of MM/GBSA, which showed that GLN107, PHE110 and MET121 had high contribution in total binding free energy (Supplementary Materials Table S2).

Larval Mortality
When the 3rd instar larvae were exposed to different concentrations of OA-02, OA-06 and OA-09, a significantly different larval mortality, both due to type of compounds and their concentrations was observed (Kruskal-Wallis H = 6.33; df = 2; p = 0.042 and Kruskal-Wallis H = 9.30; df = 2; p = 0.026, respectively). Highest mortality of 3rd instar larvae was recorded when treated with 50μM concentration of OA-09 ( Figure 3A). Since the insects are capable of selective feeding [79], it can be assumed that avoiding excessive feeding on a diet with higher concentrations may have resulted in low mortalities as compared to those fed with a lower dosage. Whereas, the mortality of 4th and 5th instar larvae were not significantly affected by the type of compounds

Growth and Longevity of Different Stages
The data have demonstrated that the 3rd-5th instar larvae treated with different compounds have significantly different weight accumulations (F = 4.309; df = 2, 24; p = 0.025, F = 15.552; df = 2, 24; p < 0.001, F = 20.407; df = 2, 24; p < 0.001, respectively), while no significant difference due to variable concentrations was observed ( Figure 4). The least weight accumulation was observed in the 3rd instar larvae treated with OA-02, while the highest weight accumulation was in the 5th instar larvae treated with the same compound ( Figure 4). In contrast, the pupal weights appeared to be significantly affected due to variable concentrations of all compounds (F = 5.640; df = 3, 24; p = 0.005) (Figure 4) where pupae in control treatment accumulated the lower weight compared to those who were fed on treated diets during larval stages (Figure 4).
When the developmental periods of larvae, pupae and adults were observed, significant effects of compound used, and their concentrations were observed. Significantly shorter stadia lengths were observed (compound × concentration: F = 140.48; df = 6, 24; p < 0.001) when the larvae were fed on diet containing 450 µM concentration of OA-09 respect to other tested compounds (Table 4). When the pupal periods were observed, a significant difference due to interaction of compound type and their concentrations were noticed (compound × concentration: F = 11.137; df = 6, 24; p < 0.001) ( Table 4). The larvae that were fed on OA-02 treated food (450 µM) ultimately spent the highest duration as pupae while all other treatments had similar pupal periods as that in the control (Table 4). In general, OA-02 and OA-06 had a similar effect on pupal period and 50 µM and 150 µM concentrations also had similar effects on pupal periods (Table 4).
Overall, when the total periods from 3rd instar larvae to adult emergence were observed, a significant difference due to interaction effect of the compounds and their concentrations were observed (compound × concentration: F = 4.044; df = 6, 24; p = 0.006) ( Table 4). The trend was similar as observed in larval periods and stadia lengths. When emerged, the compounds and concentrations also had significant impact on longevity of both the male and female moths (F = 28.226; df = 6, 24; p < 0.001 and F = 19.399; df = 6, 24; p < 0.001, respectively). The shortest-lived male and female moths were recorded in the population that was reared on the 450 µM treated diet ( Table 4). The adults from the population fed on OA-02 and OA-06 treated diet had statistically similar effects but significantly lower than control populations (Table 4).   The data are presented as mean ± SEM and the superscript uppercase and lowercase letters represent the post-hoc pairwise comparisons between compounds and the concentrations, respectively. The data with different letters in superscript differ significantly from each other (p < 0.05, Tukey's HSD).

Developmental and Reproductive Potentials
In terms of pupation (%) and adult emergence (%) significant effects of both the concentrations and compound types were observed (F = 28.926; df = 6, 24; p < 0.001 and F = 115.577; df = 6, 24; p < 0.001, respectively) ( Table 5). All the different concentrations significantly reduced pupation and adult emergence rates as compared to the control populations. The OA-02 treated populations had significantly different and higher pupation and adult emergence rates as compared to OA-06 and OA-09 (Table 5).
As far as fecundity is concerned, we observed that all compounds at different concentrations had a significant effect on fecundity and fertility (compound × concentration: F = 7.376; df = 6, 24; p < 0.001 and F = 11.591; df = 6, 24; p < 0.001, respectively) ( Table 5). The mean numbers of eggs laid did not differ significantly among the populations treated with all three concentrations but still was significantly lower as compared to the control. However, a significant effect on the egg hatch (fertility) was obvious where 450 µM of OA-09 reduced the fertility and only 58% of the eggs laid were hatched (Table 5). The data are presented as mean ± SEM and the superscript uppercase and lowercase letters atop each bar represent the post-hoc pairwise comparisons between compounds and the concentrations, respectively. The data with different letters in superscript differ significantly from each other (p < 0.05, Tukey's HSD).

Lifetable Parameters
The intrinsic (r) and finite (λ) rates of population increase were not affected by treating with different compounds, while in general, net reproductive rates (R o ) and gross reproductive rate (GRR) were reduced as compared to the control population ( Table 6). The results clearly demonstrated significantly reduced generation times among treated and control populations ( Table 6). The population that was fed on OA-09 had the least net reproductive rate coupled with shortest generation time as compared to other compounds (Table 6). Table 6. Lifetable parameters of Helicoverpa armigera populations exposed to OA-02, OA-06 and OA-09 compounds. The data presented are mean ± SEM while the superscript lowercase letters represent the post-hoc pairwise comparisons between compounds. The data with different letters in superscript differ significantly from each other (p < 0.05, Tukey's HSD). The above data were bootstrapped (100,000 times) to calculate variances and then analyzed using a one-way ANOVA.
Age-specific survival rates (l x ) and fecundity (m x ) of the populations that were fed on diets treated with OA-02, OA-06, OA-09 and control demonstrated a variable trend ( Figure 5). The population feeding on OA-02 treated diets had survived better (l x = 0.78) in the early larval stages and thereafter a quick population decline was observed until day 32 ( Figure 5). The populations which were reared on OA-06 and OA-09 treated diet had declining survival rates in the early stages of development and then they became constant before a sharp fall at the 23rd day, respectively ( Figure 5). In the control population, the decline in survival rates started at about the 37th day to the 40th day ( Figure 5). When the fecundity (m x ) was observed among these populations, a significant difference in the onset of oviposition was observed. Those reared on OA-09 diet started oviposition on the 21st day of their lifecycle, while the OA-06 reared population started egg laying on the 22nd day ( Figure 5). The control and OA-02 reared populations had oviposited at the same time, i.e., the 24th day ( Figure 5). Oviposition peaks in the population fed on OA-09 treated dyes were observed earlier, i.e., 21st day while the control population had the peak of oviposition observed around the 38th day of their lifecycle. The highest fecundity (m x ) was observed in females that emerged from the OA-06 treated population and this might be due to a sublethal effect of the compound [80]. day ( Figure 5). The control and OA-02 reared populations had oviposited at the same time, i.e., the 24th day ( Figure 5). Oviposition peaks in the population fed on OA-09 treated dyes were observed earlier, i.e., 21st day while the control population had the peak of oviposition observed around the 38th day of their lifecycle. The highest fecundity (mx) was observed in females that emerged from the OA-06 treated population and this might be due to a sublethal effect of the compound [80].

Discussion
There are many factors that affect herbivore population establishment on a particular host, availability of requisite nutrients is one of those significant factors. The variable chemical composition of the host plant impact their developmental, growth and survival [81]. Demographic evaluation [82] with either restricting the extraction of some key nutrients or blocking their metabolism after ingestion can provide an effective and ecofriendly pest control approach. In this study, we have evaluated the activity of some SCP-2 inhibitor compounds against H. armigera with the above-mentioned approach.
In previous studies, several candidate compounds from phytochemicals [83], natural products [84] and commercial databases [85] have been identified by in silico studies. So, by applying the computer aided drug design approaches, we identified three compounds from a series of quinolone azines synthetic compounds. The binding poses of the docked

Discussion
There are many factors that affect herbivore population establishment on a particular host, availability of requisite nutrients is one of those significant factors. The variable chemical composition of the host plant impact their developmental, growth and survival [81]. Demographic evaluation [82] with either restricting the extraction of some key nutrients or blocking their metabolism after ingestion can provide an effective and ecofriendly pest control approach. In this study, we have evaluated the activity of some SCP-2 inhibitor compounds against H. armigera with the above-mentioned approach.
In previous studies, several candidate compounds from phytochemicals [83], natural products [84] and commercial databases [85] have been identified by in silico studies. So, by applying the computer aided drug design approaches, we identified three compounds from a series of quinolone azines synthetic compounds. The binding poses of the docked compounds were assessed and the compounds having similar binding mode to cholesterol, a natural substrate of Ha-SCP2, were selected for the protein-ligand stability analysis by MD Simulation. The RMSD of protein backbone atoms helps to study the stability of protein [86] when a certain ligand is bound to it. The RMSD values of our studied complexes showed stable confirmations with minor deviations (Figure 2A,B). Similarly, the residual flexibility of protein was estimated by RMSF analysis [87], which indicated the stable protein structure in all three complexes ( Figure 2C). To find the binding free energy in the selected complexes, a MM/GBSA module was applied on the last 300 snapshots of MD trajectories. MM/GBSA provides accurate estimates of binding affinity in terms of binding free energy [88]. The total binding free energy values −24.42 ± 0.33, −21.35 ± 0.45, −28.19 ± 0.40 for OA-02, OA-06 and OA-09, respectively, showed that the compounds were bound with protein with reasonably good affinity. Moreover, the contribution of binding site residues in total binding free energy was estimated (Table S2) and based on the contributions, key residues Gln107, Phe110 and Met121, were identified. Further, the three selected compounds were used for in vivo evaluation.
In general, our results have demonstrated that the growth patterns and population increase are significantly affected in the populations when treated with different doses of SCP-2 inhibiting molecules. These findings have important implications in developing a targeted and environmentally safe management strategy against H. armigera. Although the higher concentrations did not kill the larvae when fed on it, the life table parameters and developmental variables were greatly affected, resulting in decreased larval and pupal weights, decreased pupation and eclosion, as well as decreased fertility (Tables 3 and 4). These findings are in accordance with earlier studies reporting a significant decrease in pupal weights of H. armigera [55,89] and fertility [90] when reared on different insecticidal admixed artificial diets.
Our findings have demonstrated that, contrary to other studies depicting the dose dependent effect on insect mortality [91,92], lower mortality was observed when the larvae were treated with the higher concentrations of the test compounds. This might be due to the nature of the compound, as the test molecules in our study were SCP-2 inhibitors that may act as insect growth regulators (IGR) rather than as a toxin [91][92][93].
We also noticed that fecundity was not significantly affected, but the fertility was greatly reduced in the populations where the larvae were fed on higher concentrations of the SCP-2 inhibitors. This clearly suggests that the compounds may have decreased palatability and cholesterol accumulation and transportation in larval bodies ultimately affecting reproductive potentials [94,95]. The second major impact of the treated compounds was observed in the form of reduced stadia lengths as even it varies when H. armigera feed on different plant species [96,97]. The larval period was shortened when treated with test compounds; however, the pupal periods remained unaffected. We assume that disruption in ample accumulation of cholesterol and sterol that is required in chitin synthesis [90,98] resulted in the larvae molting earlier than control. Clearly, OA-09 at the rate of 450 µM had demonstrated promising results where the larval periods, weights, pupal periods, weights and ultimately the female longevity and fertility were negatively affected in the test populations. Such decrement in significant life parameters of H. armigera can be attributed to reduced feeding efficiency and disruption of sterol carrier protein functions [53].

Conclusions
In conclusion, the life table parameters of H. armigera were significantly reduced when treated with OA-09. The reduced larval and pupal durations and oviposition of females may reduce population and growth. Although the generation time is also reduced due to test compounds, coupled with significantly reduced female longevity and fertility (up to 42%), it suggests that the pest establishment will be seriously compromised in terms of population growth if these novel compounds are applied. This conclusion can further be supported by at least 70% decreased net reproductive rates. Our findings have significant implications for environmentally sustainable pest control strategies by reducing the dependence on traditional chemical insecticide usage.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/insects13121169/s1, Figure S1: Interactions of the OA series compounds with the SCP2 protein; Table S1: Binding free energy calculations of complexes by MM/GBSA; Table S2: Contribution of active site residues in the binding free energies.