Exploring details about structure requirements based on antioxidant tripeptide derived from ß-Lactoglobulin by in silico approaches.

ß-Lactoglobulin is one of the proteins in milk possessing antioxidant activity. The peptides derived from ß-Lactoglobulin exhibit higher antioxidant activities than the most commonly used antioxidant. Furthermore, the detailed structure-activity relationship of these antioxidant peptides has not been elucidated. Therefore, in the present work, two-dimensional quantitative structure-activity relationship (2D-QSAR) and three-dimensional quantitative structure-activity relationship (3D-QSAR) models were constructed to investigate the structural factors affecting activities and gave information for the rational design of novel antioxidant peptides. After calculation and screening of molecular descriptors, linear and nonlinear models were developed by multiple linear regression (MLR), partial least squares regression (PLSR) and support vector machines (SVM) approaches. The statistical parameters are as follows: R2 = 0.643, Q2 = 0.553/MLR, R2 = 0.612, Q2 = 0.5278/PLSR, R2 = 0.7085, Q2 = 0.6887/SVM, indicating that the SVM model is superior to the MLR and PLSR models. In addition, in the 3D-QSAR models, the Dragon-CoMFA (R2cv = 0.537, R2pred = 0.5201) and Dragon-CoMSIA (R2cv = 0.665, R2pred = 0.6489) methods were conducted to predict the antioxidant activities. Comparison of statistical parameters illustrates that the suitability of Dragon-CoMSIA is superior to the Dragon-CoMFA model. The results show the robustness and excellent prediction of the proposed models, and would be applied for modifying and designing novel and potent antioxidant peptides.

Identification of potential Akt activators: a ligand and structure-based computational approach.

The Akt pathway plays a significant role in various diseases like Alzheimer's, Parkinson's, and Diabetes. Akt is the central protein whose phosphorylation controls many downstream pathways. Binding of small molecules to the PH domain of Akt facilitates its phosphorylation in the cytoplasm and upregulates the Akt pathway. In the current study, to identify Akt activators, ligand-based approaches like 2D QSAR, shape, and pharmacophore-based screening were used, followed by structure-based approaches such as docking, MM-GBSA, ADME prediction, and MD simulation. The top twenty-five molecules from the Asinex gold platinum database found to be active in most 2D QSAR models were used for shape and pharmacophore-based screening. Later docking was performed using the PH domain of Akt1 (PDB: 1UNQ), and 197105, 261126, 253878, 256085, and 123435 were selected based on docking score and interaction with key residues, which were druggable and formed a stable protein-ligand complex. MD simulations of 261126 and 123435 showed better stability and interactions with key residues. To further investigate the SAR of 261126 and 123435, derivatives were downloaded from PubChem, and structure-based approaches were employed. MD simulation of derivatives 12289533, 12785801, 83824832, 102479045, and 6972939 was performed, in which 83824832 and 12289533 showed interaction with key residues for a longer duration of time, proving that they may act as Akt activators.

Discovery of novel BRD4-BD2 inhibitors via in silico approaches: QSAR techniques, molecular docking, and molecular dynamics simulations.

Bromodomain-containing protein 4(BRD4) plays an important role in the occurrence and development of various malignant tumors, which has attracted the attention of scientific research institutions and pharmaceutical companies. The structural modification of most currently available BRD4 inhibitors is relatively simple, but the drug effectiveness is limited. Research has found that the inhibition of BD1 may promote the differentiation of oligodendrocyte progenitor cell; however, the inhibition of BD2 will not cause this outcome. Therefore, newly potential drugs which target BRD4-BD2 need further research. Herein, we initially built QSAR models out of 49 compounds using HQSAR, CoMFA, CoMSIA, and Topomer CoMFA technology. All of the models have shown suitable reliabilities (q2 = 0.778, 0.533, 0.640, 0.702, respectively) and predictive abilities (r2pred = 0.716, 0.6289, 0.6153, 0.7968, respectively) for BRD4-BD2 inhibitors. On the basis of QSAR results and the search of the R-group in the topomer search module, we designed 20 new compounds with high activity that showed appropriate docking score and suitable ADMET. Docking studies and MD simulation were carried out to reveal the amino acid residues (Asn351, Cys347, Tyr350, Pro293, and Asp299) at the active site of BRD4-BD2. Free energy calculations and free energy landscapes verified the stable binding results and indicated stable conformations of the complexes. These theoretical studies provide guidance and theoretical basis for designing and developing novel BRD4-BD2 inhibitors.

Synthesis, molecular modelling and QSAR study of new N-phenylacetamide-2-oxoindole benzensulfonamide conjugates as carbonic anhydrase inhibitors with antiproliferative activity.

In continuation of our previous studies to optimise potent carbonic anhydrase inhibitors, two new series of isatin N-phenylacetamide based sulphonamides were synthesised and screened for their human (h) carbonic anhydrase (EC 4.2.1.1) inhibitory activities against four isoforms hCA I, hCA II, hCA IX and hCA XII. The indole-2,3-dione derivative 2h showed the most effective inhibition profile against hCAI and hCA II (KI = 45.10, 5.87 nM) compared to acetazolamide (AAZ) as standard inhibitor. Moreover, 2h showed appreciable inhibition activity against the tumour-associated hCA XII, similar to AAZ showing KI of 7.91 and 5.70 nM, respectively. The analogs 3c and 3d showed good cytotoxicity effects, and 3c revealed promising selectivity towards lung cell line A549. Molecular docking was carried out for 2h and 3c to predict their binding conformations and affinities towards the hCA I, II, IX and XII isoforms.

Fluorine Atoms on C6H5-Corrole Affect the Interaction with Mpro and PLpro Proteases of SARS-CoV-2: Molecular Docking and 2D-QSAR Approaches.

The chymotrypsin-like cysteine protease (3CLpro, also known as main protease-Mpro) and papain-like protease (PLpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been used as the main targets for screening potential synthetic inhibitors for posterior in vitro evaluation of the most promising compounds. In this sense, the present work reports for the first time the evaluation of the interaction between Mpro/PLpro with a series of 17 porphyrin analogues-corrole (C1), meso-aryl-corrole (C2), and 15 fluorinated-meso-aryl-corrole derivatives (C3-C17) via molecular docking calculations. The impact of fluorine atoms on meso-aryl-corrole structure was also evaluated in terms of binding affinity and physical-chemical properties by two-dimensional quantitative structure-activity relationship (2D-QSAR). The presence of phenyl moieties increased the binding capacity of corrole for both proteases and depending on the position of fluorine atoms might impact positively or negatively the binding capacity. For Mpro the para-fluorine atoms might decrease drastically the binding capacity, while for PLpro there was a certain increase in the binding affinity of fluorinated-corroles with the increase of fluorine atoms into meso-aryl-corrole structure mainly from tri-fluorinated insertions. The 2D-QSAR models indicated two separated regions of higher and lower affinity for Mpro:C1-C17 based on dual electronic parameters (σI and σR), as well as one model was obtained with a correlation between the docking score value of Mpro:C2-C17 and the corresponding 13C nuclear magnetic resonance (NMR) chemical shifts of the sp2 carbon atoms (δC-1 and δC-2) of C2-C17. Overall, the fluorinated-meso-aryl-corrole derivatives showed favorable in silico parameters as potential synthetic compounds for future in vitro assays on the inhibition of SARS-CoV-2 replication.

Design, synthesis, biological evaluation of 3,5-diaryl-4,5-dihydro-1H-pyrazole carbaldehydes as non-purine xanthine oxidase inhibitors: Tracing the anticancer mechanism via xanthine oxidase inhibition.

Xanthine oxidase (XO) has been primarily targeted for the development of anti-hyperuriciemic /anti-gout agents as it catalyzes the conversion of xanthine and hypoxanthine into uric acid. XO overexpression in various cancer is very well correlated due to reactive oxygen species (ROS) production and metabolic activation of carcinogenic substances during the catalysis. Herein, we report the design and synthesis of a series of 3,5-diaryl-4,5-dihydro-1H-pyrazole carbaldehyde derivatives (2a-2x) as xanthine oxidase inhibitors (XOIs). A docking model was developed for the prediction of XO inhibitory activity of our novel compounds. Furthermore, our compounds anticancer activity results in low XO expression and XO-harboring cancer cells both in 2D and 3D-culture models are presented and discussed. Among the array of synthesized compounds, 2b and 2m emerged as potent XO inhibitors having IC50 values of 9.32 ± 0.45 µM and 10.03 ± 0.43 µM, respectively. Both compounds induced apoptosis, halted the cell cycle progression at the G1 phase, elevated ROS levels, altered mitochondrial membrane potential, and inhibited antioxidant enzymes. The levels of miRNA and expression of redox sensors in cells were also altered due to increase oxidative stress induced by our compounds. Compounds 2b and 2m hold a great promise for further development of XOIs for the treatment of XO-harboring tumors.

Structural modification of 4, 5-dihydro-[1, 2, 4] triazolo [4, 3-f] pteridine derivatives as BRD4 inhibitors using 2D/3D-QSAR and molecular docking analysis.

Cancer treatment continues to be one of the most serious public health issues in the world. The overexpression of BRD4 protein has led to a series of malignant tumors, hence the development of small molecule BRD4 protease inhibitors has always been a hot spot in the field of medical research. In this study, a series of 4,5-dihydro-[1, 2, 4] triazolo [4, 3-f] pteridine derivatives were used to establish 3D/2D-QSAR models and to discuss the relationship between inhibitor structure and activity. Four ideal models were established, including the comparative molecular field analysis (CoMFA: [Formula: see text] = 0.574, [Formula: see text] = 0.947) model, comparative molecular similarity index analysis (CoMSIA: [Formula: see text]= 0.622, [Formula: see text] = 0.916) model, topomer CoMFA ([Formula: see text] = 0.691, [Formula: see text]= 0.912) model and hologram quantitative structure-activity relationship (HQSAR: [Formula: see text]= 0.759, [Formula: see text] = 0.963) model. They show quite good external predictive power for the test set, with [Formula: see text] values of 0.602, 0.624, 0.671 and 0.750, respectively. In addition, the contour and color code map given by the 2D/3D-QSAR model with the results of molecular docking analyzed to chalk up modification methods for improving inhibitory activity, which was verified by designing novel compounds. The analysis results are helpful to promote the modification of the inhibitor framework and to provide a reference for the construction of new and promising BRD4 inhibitor compounds.

Design, Synthesis, Biological Evaluation, 2D-QSAR Modeling, and Molecular Docking Studies of Novel 1H-3-Indolyl Derivatives as Significant Antioxidants.

Novel candidates of 3-(4-(thiophen-2-yl)-pyridin/pyran/pyrimidin/pyrazol-2-yl)-1H-indole derivatives (2-12) were designed by pairing the pyridine/pyrane/pyrimidine/pyrazole heterocycles with indole and thiophene to investigate their potential activities as (2,2'-azinobis (3-ethylbenzothiazoline-6-sulfonic acid) inhibitors. The purpose of these derivatives' modification is to create high-efficiency antioxidants, especially against ABTS, as a result of the efficiency of this set of key heterocycles in the inhibition of ROS. Herein, 2D QSAR modeling was performed to recommend the most promising members for further in vitro investigations. Furthermore, the pharmacological assay for antioxidant activity evaluation of the yielded indole-based heterocycles was tested against ABTS (2,2'-azinobis (3-ethylbenzothiazoline-6-sulfonic acid); by utilizing ascorbic acid as the standard. Candidate 10 showed higher antioxidant activity (IC50 = 28.23 µg/mL) than ascorbic acid itself which achieved (IC50 = 30.03 µg/mL). Moreover, molecular docking studies were performed for the newly designed and synthesized drug candidates to propose their mechanism of action as promising cytochrome c peroxidase inhibitors compared to ascorbic acid as a reference standard. Our findings could be promising in the medicinal chemistry scope for further optimization of the newly designed and synthesized compounds regarding the introduced structure-activity relationship study (SAR) in order to get a superior antioxidant lead compound in the near future.

Design, synthesis, antifungal activity, and structure-activity relationship studies of chalcones and hybrid dihydrochromane-chalcones.

A series of ten chalcones (7a-j) and five new dihydrochromane-chalcone hybrids (7k-o) were synthesized and identified using spectroscopic techniques (IR, NMR, and MS). All compounds were evaluated in vitro against the B. cinerea and M. fructicola phytopathogens that affect a wide range of crops of commercial interest. All compounds were tested against both phytopathogens using the mycelial growth inhibition test, and it was found that two and five compounds had similar activity to that of the positive control for B. cinerea (7a = 43.9, 7c = 45.5, and Captan®= 24.8 µg/mL) and M. fructicola (7a = 48.5, 7d = 78.2, 7e = 56.1, 7f = 51.8, 7n = 63.2, and Mystic®= 21.6 µg/mL), respectively. To understand the key chalcone structural features for the antifungal activity on B. cinerea and M. fructicola, we developed structure-activity models with good statistical values (r2 and q2 higher than 0.8). For B. cinerea, the hydrogen bonding donor and acceptor and the atomic charge on C5 modulate the mycelial growth inhibition activity. In contrast, dipole moment and atomic charge on C1' and the carbonyl carbon modify the inhibition activity for M. fructicola. These results allow the design of other compounds with activities superior to those of the compounds obtained in this study.

Novel thiazolidines: Synthesis, antiproliferative properties and 2D-QSAR studies.

A series of N-substituted (Z)-2-imino-(5Z)-ylidene thiazolidines/thiazolidin-4-ones were synthesized and their antiproliferative activities against colon (HCT-116) and breast (MCF7) cancer cell lines were evaluated utilizing an MTT growth assay. A 2D-QSAR investigation was conducted to probe and validate the obtained antiproliferative properties for the thiazolidine derivatives. The majority of the thiazolidines exhibit higher potency against a colon cancer cell line relative to the standard reference. The p-halophenylimino p-anisylidene derivatives exhibited the highest anti-proliferative activity against HCT116 relative to control (IC50 = 8.9-10.0 µM compared to 20.4 µM observed for 5-fluorouracil as positive control). An X-ray study confirmed the Z, Z'-configurations for two examples of the synthesized compounds.

Structural characterization of plasmodial aminopeptidase: a combined molecular docking and QSAR-based in silico approaches.

Aminopeptidase M1 (PfAM1) is one of the key enzymes involved in the development of new antimalarials. To accelerate the discovery of inhibitors with selective activity against PfAM1 and microsomal neutral aminopeptidase (pAPN), in the present work, the optimum comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were built based on PfAM1 and pAPN inhibitors. The results of the developed 3D-QSAR models were as follows: PfAM1/CoMFA: [Formula: see text] = 0.740, [Formula: see text] = 0.7781; PfAM1/CoMSIA: [Formula: see text] = 0.740, [Formula: see text] = 0.7354; pAPN/CoMFA: [Formula: see text] = 0.612, [Formula: see text] = 0.7318; pAPN/CoMSIA: [Formula: see text] = 0.609, [Formula: see text] = 0.7480, and the models derived from MLR, PLSR and SVR methods provided high R2 values of 0.6960, 0.6965, 0.7971 for PfAM1, 0.7700, 0.7697, 0.8228 for pAPN and Q2 of 0.7004, 0.7004, 0.5632 for PfAM1, 0.7551, 0.7566 and 0.8394 for pAPN, respectively, indicating that the developed 3D-QSAR and 2D-QSAR models possess good ability for prediction of the relative compound activities. Furthermore, all inhibitors were docked into the active site of the PfAM1 and pAPN receptors, the hydrogen-bond interactions between the compound 33 with Glu497, Glu463 and Arg489 of the PfAM1, and the compound 4 with Ala348, Glu384 and Phe467 of the receptor pAPN are able to help to stabilize the conformation. The above results would provide helpful clues to predicting the binding activity of novel inhibitors and the foundation for understanding the interaction mechanism between the inhibitors and the receptors.

In silico toxicity evaluation of dioxins using structure-activity relationship (SAR) and two-dimensional quantitative structure-activity relationship (2D-QSAR).

Prediction of pEC50 values of dioxins binding with the aryl hydrocarbon receptor (AhR) is of great significance for exploring how dioxins induce toxicity in human body and evaluating their environmental behaviors and risks. To reveal the factors that influence the toxicity of dioxins, provide more accurate mathematical models for predicting the pEC50 values of dioxins, and supplement the toxicity database of persistent organic pollutants, qualitative structure-activity relationship (SAR) and two-dimensional quantitative structure-activity relationship (2D-QSAR) were used in this study. The research objects in this study were 60 organic compounds with pEC50 values and 162 compounds without pEC50 values, which included polychlorinated dibenzofurans (PCDFs), polychlorinated dibenzo-p-dioxins (PCDDs), and polybrominated dibenzo-p-dioxins (PBDDs). The qualitative structure-activity relationship (SAR) was performed first and concluded that halogen substitutions at any of the 2, 3, 7, and 8 sites increased the pEC50 value of the compound. Moreover, two-dimensional quantitative structure-activity relationship (2D-QSAR) models were established by employing multiple linear regression (MLR) method and artificial neural network (ANN) algorithm to investigate the factors affecting the pEC50 values of dioxins molecules. MLR was used to establish the well-understood linear model and ANN was used to establish a more accurate non-linear model. Both models have good fitting, robustness, and predictive ability. Importantly, the ability of dioxins binding to AhR is mainly determined by molecular descriptors including E1m, SM09_AEA (dm), RDF065u, F05 [Cl-Cl], and Neoplastic-80. In addition, the pEC50 values of the 162 dioxins without toxicity data were predicted by MLR and ANN models, respectively.

Application of two-dimensional binary fingerprinting methods for the design of selective Tankyrase I inhibitors.

In the present study, five important binary fingerprinting techniques were used to model novel flavones for the selective inhibition of Tankyrase I. From the fingerprints used: the fingerprint atom pairs resulted in a statistically significant 2D QSAR model using a kernel-based partial least square regression method. This model indicates that the presence of electron-donating groups positively contributes to activity, whereas the presence of electron withdrawing groups negatively contributes to activity. This model could be used to develop more potent as well as selective analogues for the inhibition of Tankyrase I. Schematic representation of 2D QSAR work flow.

Novel coumarin-6-sulfonamides as apoptotic anti-proliferative agents: synthesis, in vitro biological evaluation, and QSAR studies.

Herein, we report the synthesis of different novel sets of coumarin-6-sulfonamide derivatives bearing different functionalities (4a, b, 8a-d, 11a-d, 13a, b, and 15a-c), and in vitro evaluation of their growth inhibitory activity towards the proliferation of three cancer cell lines; HepG2 (hepatocellular carcinoma), MCF-7 (breast cancer), and Caco-2 (colon cancer). HepG2 cells were the most sensitive cells to the influence of the target coumarins. Compounds 13a and 15a emerged as the most active members against HepG2 cells (IC50 = 3.48 ± 0.28 and 5.03 ± 0.39 µM, respectively). Compounds 13a and 15a were able to induce apoptosis in HepG2 cells, as assured by the upregulation of the Bax and downregulation of the Bcl-2, besides boosting caspase-3 levels. Besides, compound 13a induced a significant increase in the percentage of cells at Pre-G1 by 6.4-folds, with concurrent significant arrest in the G2-M phase by 5.4-folds compared to control. Also, 13a displayed significant increase in the percentage of annexin V-FITC positive apoptotic cells from 1.75-13.76%. Moreover, QSAR models were established to explore the structural requirements controlling the anti-proliferative activities.

Synthesis and 2D-QSAR Study of Active Benzofuran-Based Vasodilators.

A new series of 2-alkyloxy-pyridine-3-carbonitrile-benzofuran hybrids (4a-x) was synthesized. All the new derivatives were examined via the standard technique for their vasodilation activity. Some of the investigated compounds exhibited a remarkable activity, with compounds 4w, 4e, 4r, 4s, 4f and 4g believed to be the most active hits in this study with IC50 values 0.223, 0.253, 0.254, 0.268, 0.267 and 0.275 mM, respectively, compared with amiodarone hydrochloride, the reference standard used (IC50 = 0.300 mM). CODESSA PRO was employed to obtain a statistically significant 2-Dimensional Quantitative Structure Activity Relationship (2D-QSAR) model describing the bioactivity of the newly synthesized analogs (N = 24, n = 4, R² = 0.816, R²cvOO = 0.731, R²cvMO = 0.772, F = 21.103, s² = 6.191 × 10-8).

Quantitative insights towards the design of potent deazaxanthine antagonists of adenosine 2B receptors.

Adenosine receptors have been considered as potential targets for drug development, but one of the main obstacles to this goal is to selectively inhibit one receptor subtype over the others. This subject is particularly crucial for adenosine A2b receptor antagonists (AdoRA2B). The structure­activity relationships of xanthine derivatives which are AdoRA2B have been comprehensively investigated, but the steric and electronic requirements of deazaxanthine AdoRA2B have not been described from a quantitative standpoint of view. Herein we report our efforts to shorten this knowledge gap through 2D-QSAR (HQSAR) and 3D-QSAR (CoMFA) approaches. The good statistical quality (HQSAR--r(2) = 0.85, q(2)(LOO) = 0.77; CoMFA ­ r(2) = 0.86, q(2) = 0.70) and predictive ability (r(2) = (pred1) = 0.78, r(2)(pred2) = 0.78 and r(2) = (pred1) = 0.70, r(2) = (pred2) = 0.70,respectively) of the models, along with the information provided by contribution and contour maps hints their usefulness to the design of more potent 9-deazaxanthine derivatives.

2D-QSAR, molecular docking and MD simulation based virtual screening of the herbal molecules against Alzheimer's disorder: an approach to predict CNS activity.

Acetylcholinesterase (AChE) is one of the key enzyme targets that have been used clinically for the management of Alzheimer's Disorder (AD). Numerous reports in the literature predict and demonstrate in-vitro, and in-silico anticholinergic activity of herbal molecules, however, majority of them failed to find clinical application. To address these issues, we developed a 2D-QSAR model that could efficiently predict the AChE inhibitory activity of herbal molecules along with predicting their potential to cross the blood-brain-barrier (BBB) to exert their beneficial effects during AD. Virtual screening of the herbal molecules was performed and amentoflavone, asiaticoside, astaxanthin, bahouside, biapigenin, glycyrrhizin, hyperforin, hypericin, and tocopherol were predicted as the most promising herbal molecules for inhibiting AChE. Results were validated through molecular docking, atomistic molecular dynamics simulations and Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) studies against human AChE (PDB ID: 4EY7). To determine whether or not these molecules can cross BBB to inhibit AChE within the central nervous system (CNS) for being beneficial for the management of AD, we determined a CNS Multi-parameter Optimization (MPO) score, which was found in the range of 1 to 3.76. Overall, the best results were observed for amentoflavone and our results demonstrated a PIC50 value of 7.377 nM, molecular docking score of -11.5 kcal/mol, and CNS MPO score of 3.76. In conclusion, we successfully developed a reliable and efficient 2D-QSAR model and predicted amentoflavone to be the most promising molecule that could inhibit human AChE enzyme within the CNS and could prove beneficial for the management of AD.Communicated by Ramaswamy H. Sarma.

A combined in silico approaches of 2D-QSAR, molecular docking, molecular dynamics and ADMET prediction of anti-cancer inhibitor activity for actinonin derivatives.

Inhibition of human mitochondrial peptide deformylase (HsPDF) plays a major role in reducing growth, proliferation, and cellular cancer survival. In this work, a series of 32 actinonin derivatives for HsPDF (PDB: 3G5K) inhibitor's anticancer activity was computationally analyzed for the first time, using an in silico study considering 2D-QSAR modeling, and molecular docking studies, and validated by molecular dynamics and ADMET properties. The results of multilinear regression (MLR) and artificial neural networks (ANN) statistical analysis reveal a good correlation between pIC50 activity and the seven (7) descriptors. The developed models were highly significant with cross-validation, the Y-randomization test and their applicability range. In addition, all considered data sets show that the AC30 compound, exhibits the best binding affinity (docking score = -212.074 kcal/mol and H-bonding energy = -15.879 kcal/mol). Furthermore, molecular dynamics simulations were performed at 500 ns, confirming the stability of the studied complexes under physiological conditions and validating the molecular docking results. Five selected actinonin derivatives (AC1, AC8, AC15, AC18 and AC30), exhibiting best docking score, were rationalized as potential leads for HsPDF inhibition, in well agreement with experimental outcomes. Furthermore, based on the in silico study, new six molecules (AC32, AC33, AC34, AC35, AC36 and AC37) were suggested as HsPDF inhibition candidates, which would be combined with in-vitro and in-vivo studies to perspective validation of their anticancer activity. Indeed, the ADMET predictions indicate that these six new ligands have demonstrated a fairly good drug-likeness profile.

Risk mitigation strategy and mechanism analysis of neonicotinoid pesticides on earthworms.

Neonicotinoid insecticides (NNIs) are a new class of widely used insecticides with certain risks to non-target organisms, like earthworms. The gray correlation method was used to calculate the comprehensive risk effect value of acute toxicity (LC50) and bioaccumulation (logKow) of NNIs on earthworms. A comprehensive effects three-dimensional quantitative structure-activity relationship (3D-QSAR) model was constructed, using NNIs molecular structures and the comprehensive effect value as the independent and dependent variables, respectively. One of the representatives guadipyr (GUA) was selected as the template molecule for the molecular design and modification. A total of 63 NNIs alternatives were designed with a reduced comprehensive value higher than 10%, and as high as 42%. After screening, 15 NNIs alternatives were screened with decreased acute toxicity to earthworms, bioaccumulation effects and improved functional property. The calculated primary acute risk quotient of earthworms shows that the designed NNIs alternatives have lower earthworm risks (reduction of 70.48-99.99%). Results also found that the electronic, geometric and topological parameters of NNIs are the key descriptors that affect NNIs alternatives' toxicity. The number of hydrophobic interaction amino acid residues in NNIs molecules also contributes to the acute toxicity and the bioaccumulation of NNIs alternatives on earthworms. This study aims to design and screen functionally improved and environmentally friendly NNIs alternatives that have low risk to earthworms and provide theoretical methods and new ideas for the risk control and development of pesticides represented by NNIs.

An in silico study on human carcinogenicity mechanism of polybrominated biphenyls exposure.

Polybrominated biphenyls (PBBs) are associated with an increased risk of thyroid cancer; however, relevant mechanistic studies are lacking. In this study, we investigated the mechanisms underlying PBB-induced human thyroid cancer. Molecular docking and molecular dynamics methods were employed to investigate the metabolism of PBBs by the cytochrome P450 enzyme under aryl hydrocarbon receptor mediation into mono- and di-hydroxylated metabolites. This was taken as the molecular initiation event. Subsequently, considering the interactions of PBBs and their metabolites with the thyroxine-binding globulin protein as key events, an adverse outcome pathway for thyroid cancer caused by PBBs exposure was constructed. Based on 2D quantitative structure activity relationship (2D-QSAR) models, the contribution of amino acid residues and binding energy were analyzed to understand the mechanism underlying human carcinogenicity (adverse effect) of PBBs. Hydrogen bond and van der Waals interactions were identified as key factors influencing the carcinogenic adverse outcome pathway of PBBs. Analysis of non-bonding forces revealed that PBBs and their hydroxylation products were predominantly bound to the thyroxine-binding globulin protein through hydrophobic and hydrogen bond interactions. The key amino acids involved in hydrophobic interactions were alanine 330, arginine 381 and lysine 270, and the key amino acids involved in hydrogen bond interactions were arginine 381 and lysine 270. This study provides valuable insights into the mechanisms underlying human health risk associated with PBBs exposure.