In Silico Screening of Novel α1-GABAA Receptor PAMs towards Schizophrenia Based on Combined Modeling Studies of Imidazo [1,2-a]-Pyridines

The ionotropic GABA_A receptor (GABA_AR) has been proven to be an important target of atypical antipsychotics. A novel series of imidazo [1,2-a]-pyridine derivatives, as selective positive allosteric modulators (PAMs) of α1-containing GABA_ARs with potent antipsychotic activities, have been reported recently. To better clarify the pharmacological essentiality of these PAMs and explore novel antipsychotics hits, three-dimensional quantitative structure–activity relationships (3D-QSAR), molecular docking, pharmacophore modeling, and molecular dynamics (MD) were performed on 33 imidazo [1,2-a]-pyridines. The constructed 3D-QSAR models exhibited good predictive abilities. The dockings results and MD simulations demonstrated that hydrogen bonds, π–π stackings, and hydrophobic interactions play essential roles in the binding of these novel PAMs in the GABA_AR binding pocket. Four hit compounds (DS01–04) were then screened out by the combination of the constructed models and computations, including the pharmacophore model, Topomer Search, molecular dockings, ADME/T predictions, and MD simulations. The compounds DS03 and DS04, with higher docking scores and better predicted activities, were also found to be relatively stable in the binding pocket by MD simulations. These results might provide a significant theoretical direction or information for the rational design and development of novel α1-GABA_AR PAMs with antipsychotic activities.     

Structural Determination of Three Different Series of Compounds as Hsp90 Inhibitors Using 3D-QSAR Modeling, Molecular Docking and Molecular Dynamics Methods

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q(2) > 0.60 and external r(2) (pred) > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.     

噻吩并嘧啶酮类PDE7抑制剂的特征结构及作用机制研究

基于定量构效关系和分子对接研究理论,以34个噻吩并嘧啶酮类磷酸二酯酶7(PDE7)抑制剂为研究对象,采用比较分子相似性指数法(Co MSIA)和GOLD分子对接法研究其特征结构信息以及与靶标的作用机制。结果表明,所建立Co MSIA模型的预测能力较强(Q2=0.53,R2ncv=0.96,R2pre=0.89),结合分子对接结果还发现,抑制剂与PDE7间的作用力以氢键作用力和π-π堆积作用为主;大体积的亲水性氢键供体R1取代基、亲水性R2取代基有利于化合物活性的增强,且哌啶环上有益于结合氢键的受体基团。所得模型和信息为后续新型PDE7抑制剂的设计开发提供理论指导。     

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein–protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein–protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.     

BRET Analysis of GPCR Dimers in Neurons and Non-Neuronal Cells: Evidence for Inactive,Agonist, and Constitutive Conformations

G-protein-coupled receptors (GPCRs) are dimeric proteins, but the functional consequences of the process are still debated. Active GPCR conformations are promoted either by agonists or constitutive activity. Inverse agonists decrease constitutive activity by promoting inactive conformations. The histamine H₃ receptor (H₃R) is the target of choice for the study of GPCRs because it displays high constitutive activity. Here, we study the dimerization of recombinant and brain H₃R and explore the effects of H₃R ligands of different intrinsic efficacy on dimerization. Co-immunoprecipitations and Western blots showed that H₃R dimers co-exist with monomers in transfected HEK 293 cells and in rodent brains. Bioluminescence energy transfer (BRET) analysis confirmed the existence of spontaneous H₃R dimers, not only in living HEK 293 cells but also in transfected cortical neurons. In both cells, agonists and constitutive activity of the H₃R decreased BRET signals, whereas inverse agonists and GTPγS, which promote inactive conformations, increased BRET signals. These findings show the existence of spontaneous H₃R dimers not only in heterologous systems but also in native tissues, which are able to adopt a number of allosteric conformations, from more inactive to more active states.     

Combined 3D-QSAR, Molecular Docking and Molecular Dynamics Study on Derivatives of Peptide Epoxyketone and Tyropeptin-Boronic Acid as Inhibitors Against the β5 Subunit of Human 20S Proteasome

An abnormal ubiquitin-proteasome is found in many human diseases, especially in cancer, and has received extensive attention as a promising therapeutic target in recent years. In this work, several in silico models have been built with two classes of proteasome inhibitors (PIs) by using 3D-QSAR, homology modeling, molecular docking and molecular dynamics (MD) simulations. The study resulted in two types of satisfactory 3D-QSAR models, i.e., the CoMFA model (Q(2) = 0.462, R(2) (pred) = 0.820) for epoxyketone inhibitors (EPK) and the CoMSIA model (Q(2) = 0.622, R(2) (pred) = 0.821) for tyropeptin-boronic acid derivatives (TBA). From the contour maps, some key structural factors responsible for the activity of these two series of PIs are revealed. For EPK inhibitors, the N-cap part should have higher electropositivity; a large substituent such as a benzene ring is favored at the C6-position. In terms of TBA inhibitors, hydrophobic substituents with a larger size anisole group are preferential at the C8-position; higher electropositive substituents like a naphthalene group at the C3-position can enhance the activity of the drug by providing hydrogen bond interaction with the protein target. Molecular docking disclosed that residues Thr60, Thr80, Gly106 and Ser189 play a pivotal role in maintaining the drug-target interactions, which are consistent with the contour maps. MD simulations further indicated that the binding modes of each conformation derived from docking is stable and in accord with the corresponding structure extracted from MD simulation overall. These results can offer useful theoretical references for designing more potent PIs.     

Effect of Ions and Sequence Variants on the Antagonist Binding Properties of the Histamine H1 Receptor

The histamine H₁ receptor (H₁R) is a G protein-coupled receptor (GPCR) and represents a main target in the treatment of allergic reactions as well as inflammatory reactions and depressions. Although the overall effect of antagonists on H₁ function has been extensively investigated, rather little is known about the potential modulatory effect of ions or sequence variants on antagonist binding. We investigated the dynamics of a phosphate ion present in the crystal structure and of a sodium ion, for which we determined the position in the allosteric pocket by metadynamics simulations. Both types of ions exhibit significant dynamics within their binding site; however, some key contacts remain stable over the simulation time, which might be exploited to develop more potent drugs targeting these sites. The dynamics of the ions is almost unaffected by the presence or absence of doxepin, as also reflected in their small effect (less than 1 kcal·mol⁻¹) on doxepin binding affinity. We also examined the effect of four H₁R sequence variants observed in the human population on doxepin binding. These variants cause a reduction in doxepin affinity of up to 2.5 kcal·mol⁻¹, indicating that personalized medical treatments that take into account individual mutation patterns could increase precision in the dosage of GPCR-targeting drugs.     

Molecular Effects of Auto-Antibodies on Angiotensin II Type 1 Receptor Signaling and Cell Proliferation

The angiotensin II (Ang II) type 1 receptor (AT₁R) is involved in the regulation of blood pressure (through vasoconstriction) and water and ion homeostasis (mediated by interaction with the endogenous agonist). AT₁R can also be activated by auto-antibodies (AT₁R-Abs), which are associated with manifold diseases, such as obliterative vasculopathy, preeclampsia and systemic sclerosis. Knowledge of the molecular mechanisms related to AT₁R-Abs binding and associated signaling cascade (dys-)regulation remains fragmentary. The goal of this study was, therefore, to investigate details of the effects of AT₁R-Abs on G-protein signaling and subsequent cell proliferation, as well as the putative contribution of the three extracellular receptor loops (ELs) to Abs-AT₁R signaling. AT₁R-Abs induced nuclear factor of activated T-cells (NFAT) signaling, which reflects G_q/11 and G_i activation. The impact on cell proliferation was tested in different cell systems, as well as activation-triggered receptor internalization. Blockwise alanine substitutions were designed to potentially investigate the role of ELs in AT₁R-Abs-mediated effects. First, we demonstrate that Ang II-mediated internalization of AT₁R is impeded by binding of AT₁R-Abs. Secondly, exclusive AT₁R-Abs-induced G_q/11 activation is most significant for NFAT stimulation and mediates cell proliferation. Interestingly, our studies also reveal that ligand-independent, baseline AT₁R activation of G_i signaling has, in turn, a negative effect on cell proliferation. Indeed, inhibition of G_i basal activity potentiates proliferation triggered by AT₁R-Abs. Finally, although AT₁R containing EL1 and EL3 blockwise alanine mutations were not expressed on the human embryonic kidney293T (HEK293T) cell surface, we at least confirmed that parts of EL2 are involved in interactions between AT₁R and Abs. This current study thus provides extended insights into the molecular action of AT₁R-Abs and associated mechanisms of interrelated pathogenesis.     

Structural Model for the Binding Sites of Allosterically Potentiating Ligands on Nicotinic Acetylcholine Receptors

Current treatments of Alzheimer's disease include the allosteric potentiation of nicotinic acetylcholine receptor (nAChR) response. The location of the binding site for allosteric potentiating ligands (APLs) within the receptor is not yet fully understood. Based on homology models for the ligand binding domain of human α7, human α4β2, and chicken α7 receptors, as well as blind docking experiments with galanthamine, physostigmine, codeine, and 5HT, we identified T197 as an essential element of the APL binding site at the outer surface of the ligand binding domain (LBD) of nAChR. We also found the previously known galanthamine binding site in the region of K123 at the inside of the receptor funnel, which, however, was shown to not be part of the APL site. Our results are verified by site‐directed mutagenesis and electrophysiological experiments, and suggest that APL and ACh bind to different sites on nicotinic receptors and that allosteric potentiation may arise from a direct interplay between both these sites.     

Design,Synthesis and Evaluation of 2,4-Diaminoquinazoline Derivatives as Potential Tubulin Polymerization Inhibitors

Microtubules are highly dynamic polymers composed of α- and β-tubulin proteins that have been shown to be potential therapeutic targets for the development of anticancer drugs. Currently, a wide variety of chemically diverse agents that bind to β-tubulin have been reported. Nocodazole (NZ) and colchicine (COL) are well-known tubulin-depolymerizing agents that have close binding sites in the β-tubulin. In this study, we designed and synthesized a set of nine 2,4-diaminoquinazoline derivatives that could occupy both NZ and COL binding sites. The synthesized compounds were evaluated for their antiproliferative activities against five cancer cell lines (PC-3, HCT-15, MCF-7, MDA-MB-231, and SK-LU-1), a noncancerous one (COS-7), and peripheral blood mononuclear cells (PBMC). The effect of compounds 4 e and 4 i on tubulin organization and polymerization was analyzed on the SK-LU-1 cell line by indirect immunofluorescence, western blotting, and tubulin polymerization assays. Our results demonstrated that both compounds exert their antiproliferative activity by inhibiting tubulin polymerization. Finally, a possible binding pose of 4 i in the NZ/COL binding site was determined by using molecular docking and molecular dynamics (MD) approaches. To our knowledge, this is the first report of non-N-substituted 2,4-diaminoquinazoline derivatives with the ability to inhibit tubulin polymerization.     

Structural requirements of N-substituted spiropiperidine analogues as agonists of nociceptin/orphanin FQ receptor

The nociceptin/orphanin FQ (NOP) receptor is involved in a wide range of biological functions, including pain, anxiety, depression and drug abuse. Especially, its agonists have great potential to be developed into anxiolytics. In this work, both the ligand- and receptor-based three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were carried out using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques on 103 N-substituted spiropiperidine analogues as NOP agonists. The resultant optimal ligand-based CoMSIA model exhibited Q(2) of 0.501, R(2) (ncv) of 0.912 and its predictive ability was validated by using an independent test set of 26 compounds which gave R(2) (pred) value of 0.818. In addition, docking analysis and molecular dynamics simulation (MD) were also applied to elucidate the probable binding modes of these agonists. Interpretation of the 3D contour maps, in the context of the topology of the active site of NOP, provided insight into the NOP-agonist interactions. The information obtained from this work can be used to accurately predict the binding affinity of related agonists and also facilitate the future rational design of novel agonists with improved activity.     

Modulation of CD14 and TLR4⋅MD‐2 Activities by a Synthetic Lipid A Mimetic

Monosaccharide lipid A mimetics based on a glucosamine core linked to two fatty acid chains and bearing one or two phosphate groups have been synthesized. Compounds 1 and 2 , each with one phosphate group, were practically inactive in inhibiting LPS‐induced TLR4 signaling and cytokine production in HEK‐blue cells and murine macrophages, but compound 3 , with two phosphate groups, was found to be active in efficiently inhibiting TLR4 signal in both cell types. The direct interaction between compound 3 and the MD‐2 coreceptor was investigated by NMR spectroscopy and molecular modeling/docking analysis. This compound also interacts directly with the CD14 receptor, stimulating its internalization by endocytosis. Experiments on macrophages show that the effect on CD14 reinforces the activity on MD‐2 ? TLR4 because compound 3 's activity is higher when CD14 is important for TLR4 signaling (i.e., at low LPS concentration). The dual targeting of MD‐2 and CD14, accompanied by good solubility in water and lack of toxicity, suggests the use of monosaccharide 3 as a lead compound for the development of drugs directed against TLR4related syndromes.     

Combining Molecular Docking and Molecular Dynamics to Predict the Binding Modes of Flavonoid Derivatives with the Neuraminidase of the 2009 H1N1 Influenza A Virus

Control of flavonoid derivatives inhibitors release through the inhibition of neuraminidase has been identified as a potential target for the treatment of H1N1 influenza disease. We have employed molecular dynamics simulation techniques to optimize the 2009 H1N1 influenza neuraminidase X-ray crystal structure. Molecular docking of the compounds revealed the possible binding mode. Our molecular dynamics simulations combined with the solvated interaction energies technique was applied to predict the docking models of the inhibitors in the binding pocket of the H1N1 influenza neuraminidase. In the simulations, the correlation of the predicted and experimental binding free energies of all 20 flavonoid derivatives inhibitors is satisfactory, as indicated by R(2) = 0.75.     

Recognition and Conformational Properties of an Alternative Antithrombin Binding Sequence Obtained by Chemoenzymatic Synthesis

Heparin is a highly sulfated glycosaminoglycan (GAG) of natural origin used as an anticoagulant and antithrombotic drug. These properties are principally based on the binding and activation of antithrombin (AT) through the pentasaccharide sequence GlcNAc/NS,6S‐GlcA‐GlcNS,3,6S‐IdoA2S‐GlcNS,6S (AGA*IA). Literature data show that the population of the ²S₀ ring conformation of the 2‐O‐sulfo‐α‐l ‐iduronic acid (IdoA2S) motif correlates with the affinity and activation of AT. It was recently demonstrated that two synthetic AGA*IA‐containing hexasaccharides (one G unit added at the reducing end), differing in the degree of sulfation of the IdoA unit, show comparable affinity and ability to activate AT, despite a different conformation of the IdoA residue. In this paper, the binding of these two glycans to AT was studied by isothermal titration microcalorimetry (ITC), transferred (tr‐) NOESY, saturation transfer difference (STD) NMR spectroscopy and molecular dynamics (MD) simulations. Results indicated that both the IdoA2S and the IdoA units assume a ²S₀ conformation when bound with AT, and so present a common binding epitope for the two glycans, centred on the AGA*IA sequence.     

Specific Engineered G Protein Coupling to Histamine Receptors Revealed from Cellular Assay Experiments and Accelerated Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are targets of extracellular stimuli and hence occupy a key position in drug discovery. By specific and not yet fully elucidated coupling profiles with α subunits of distinct G protein families, they regulate cellular responses. The histamine H₂ and H₄ receptors (H₂R and H₄R) are prominent members of Gs- and Gi-coupled GPCRs. Nevertheless, promiscuous G protein and selective Gi signaling have been reported for the H₂R and H₄R, respectively, the molecular mechanism of which remained unclear. Using a combination of cellular experimental assays and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigated the coupling profiles of the H₂R and H₄R to engineered mini-G proteins (mG). We obtained coupling profiles of the mGs, mGsi, or mGsq proteins to the H₂R and H₄R from the mini-G protein recruitment assays using HEK293T cells. Compared to H₂R–mGs expressing cells, histamine responses were weaker (pEC₅₀, E_max) for H₂R–mGsi and –mGsq. By contrast, the H₄R selectively bound to mGsi. Similarly, in all-atom GaMD simulations, we observed a preferential binding of H₂R to mGs and H₄R to mGsi revealed by the structural flexibility and free energy landscapes of the complexes. Although the mG α5 helices were consistently located within the HR binding cavity, alternative binding orientations were detected in the complexes. Due to the specific residue interactions, all mG α5 helices of the H₂R complexes adopted the Gs-like orientation toward the receptor transmembrane (TM) 6 domain, whereas in H₄R complexes, only mGsi was in the Gi-like orientation toward TM2, which was in agreement with Gs- and Gi-coupled GPCRs structures resolved by X-ray/cryo-EM. These cellular and molecular insights support (patho)physiological profiles of the histamine receptors, especially the hitherto little studied H₂R function in the brain, as well as of the pharmacological potential of H₄R selective drugs.     

Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.     

Antiviral Activity of Two Marine Carotenoids against SARS-CoV-2 Virus Entry In Silico and In Vitro

The marine carotenoids fucoxanthin and siphonaxanthin are powerful antioxidants that are attracting focused attention to identify a variety of health benefits and industry applications. In this study, the binding energy of these carotenoids with the SARS-CoV-2 Spike-glycoprotein was predicted by molecular docking simulation, and their inhibitory activity was confirmed with SARS-CoV-2 pseudovirus on HEK293 cells overexpressing angiotensin-converting enzyme 2 (ACE2). Siphonaxanthin from Codium fragile showed significant antiviral activity with an IC50 of 87.4 μM against SARS-CoV-2 pseudovirus entry, while fucoxanthin from Undaria pinnatifida sporophyll did not. The acute toxicities were predicted to be relatively low, and pharmacokinetic predictions indicate GI absorption. Although further studies are needed to elucidate the inhibition of viral infection by siphonaxanthin, these results provide useful information in the application of these marine carotenoids for the treatment and prevention of COVID-19.     

Veratridine Can Bind to a Site at the Mouth of the Channel Pore at Human Cardiac Sodium Channel NaV1.5

The cardiac sodium ion channel (Na_V1.5) is a protein with four domains (DI-DIV), each with six transmembrane segments. Its opening and subsequent inactivation results in the brief rapid influx of Na⁺ ions resulting in the depolarization of cardiomyocytes. The neurotoxin veratridine (VTD) inhibits Na_V1.5 inactivation resulting in longer channel opening times, and potentially fatal action potential prolongation. VTD is predicted to bind at the channel pore, but alternative binding sites have not been ruled out. To determine the binding site of VTD on Na_V1.5, we perform docking calculations and high-throughput electrophysiology experiments in the present study. The docking calculations identified two distinct binding regions. The first site was in the pore, close to the binding site of Na_V1.4 and Na_V1.5 blocking drugs in experimental structures. The second site was at the “mouth” of the pore at the cytosolic side, partly solvent-exposed. Mutations at this site (L409, E417, and I1466) had large effects on VTD binding, while residues deeper in the pore had no effect, consistent with VTD binding at the mouth site. Overall, our results suggest a VTD binding site close to the cytoplasmic mouth of the channel pore. Binding at this alternative site might indicate an allosteric inactivation mechanism for VTD at Na_V1.5.     

Neuroprotective Effect of Aurantio-Obtusin,a Putative Vasopressin V1A Receptor Antagonist,on Transient Forebrain Ischemia Mice Model

Traditional Chinese medicines (TCMs) have been a rich source of novel drug discovery, and Cassia seed is one of the common TCMs with numerous biological effects. Based on the existing reports on neuroprotection by Cassia seed extract, the present study aims to search possible pharmacological targets behind the neuroprotective effects of the Cassia seeds by evaluating the functional effect of specific Cassia compounds on various G-protein-coupled receptors. Among the four test compounds (cassiaside, rubrofusarin gentiobioside, aurantio-obtusin, and 2-hydroxyemodin 1-methylether), only aurantio-obtusin demonstrated a specific V_1AR antagonist effect (71.80 ± 6.0% inhibition at 100 µM) and yielded an IC₅₀ value of 67.70 ± 2.41 μM. A molecular docking study predicted an additional interaction of the hydroxyl group at C6 and a methoxy group at C7 of aurantio-obtusin with the Ser341 residue as functional for the observed antagonist effect. In the transient brain ischemia/reperfusion injury C57BL/6 mice model, aurantio-obtusin attenuated the latency time that was reduced in the bilateral common carotid artery occlusion (BCCAO) groups. Likewise, compared to neuronal damage in the BCCAO groups, treatment with aurantio-obtusin (10 mg/kg, p.o.) significantly reduced the severity of damage in medial cornu ammonis 1 (mCA1), dorsal CA1, and cortex regions. Overall, the findings of this study highlight V_1AR as a possible target of aurantio-obtusin for neuroprotection.