In vitro and in silico cholinesterase inhibitory and antioxidant effects of essential oils and extracts of two new Salvia fruticosa mill. cultivars (Turgut and Uysal) and GC-MS analysis of the essential oils

The EtOH extracts of the leaves of two new cultivars (Uysal-SFU and Turgut-SFT) of Salvia fruticosa Mill. was tested against acetylcholinesterase (IC50: 30.62 ± 3.27 and 32.97 ± 2.33 µg/mL for SFU and SFT, respectively) and butyrylcholinesterase (IC50: 69.91 ± 1.08 µg/mL and 86.55 ± 1.26 µg/mL), respectively, relevant to Alzheimer's disease. The essential oils showed a stumpy inhibition against AChE and no inhibition against BChE. DPPH radical scavenging activity of the extracts (86.70 ± 0.17% and 86.14 ± 1.13% for SFU and SFT, respectively) was stronger than that of quercetin (85.51 ± 0.17%): Their (1.24 ± 0.05 and 1.04 ± 0.16 for SFU and SFT, respectively) ferric-reducing antioxidant power were close to that of the reference (e.g. quercetin, 1.42 ± 0.14). Molecular docking simulations were performed on their major monoterpenes. Our findings revealed that the leaf EtOH extracts of two cultivars are promising inhibitors of both AChE and BChE.

Cholinesterase Inhibitory and In Silico Toxicity Assessment of Thirty-Four Isoquinoline Alkaloids - Berberine as the Lead Compound

BACKGROUND

Cholinesterase (ChE) inhibitors used currently in clinics for the treatment of Alzheimer's disease (AD) are the most prescribed drug class with nitrogen-containing chemical formula. Galanthamine, the latest generation anti-ChE drug, contains an isoquinoline structure.

OBJECTIVE

The aim of the current study was to investigate the inhibitory potential of thirty-four isoquinoline alkaloids, e.g. (-)-adlumidine, β-allocryptopine, berberine, (+)-bicuculline, (-)-bicuculline, (+)-bulbocapnine, (-)-canadine, (±)-chelidimerine, corydaldine, (±)-corydalidzine, (-)-corydalmine, (+)-cularicine, dehydrocavidine, (+)-fumariline, (-)-fumarophycine, (+)-α-hydrastine, (+)-isoboldine, 13-methylcolumbamine, (-)-norjuziphine, norsanguinarine, (-)-ophiocarpine, (-)-ophiocarpine-N-oxide, oxocularine, oxosarcocapnine, palmatine, (+)-parfumine, protopine, (+)-reticuline, sanguinarine, (+)-scoulerine, (±)-sibiricine, (±)-sibiricine acetate, (-)-sinactine, and (-)-stylopine isolated from several Fumaria (fumitory) and Corydalis species towards acetyl- (AChE) and butyrylcholinesterase (BChE) by microtiter plate assays. The alkaloids with strong ChE inhibition were proceeded to molecular docking simulations as well as in silico toxicity screening for their mutagenic capacity through VEGA QSAR (AMES test) consensus model and VEGA platform as statistical approaches. The inputs were evaluated in a simplified molecular input-line entry system (SMILES).

RESULTS

ChE inhibition assays indicated that the highest AChE inhibition was caused by berberine (IC50: 0.72 ± 0.04 µg/mL), palmatine (IC50: 6.29 ± 0.61 µg/mL), β-allocryptopine (IC50: 10.62 ± 0.45 µg/mL), (-)-sinactine (IC50: 11.94 ± 0.44 µg/mL), and dehydrocavidine (IC50: 15.01 ± 1.87 µg/mL) as compared to that of galanthamine (IC50: 0.74 ± 0.01 µg/mL), the reference drug with isoquinoline skeleton. Less number of the tested alkaloids exhibited notable BChE inhibition. Among them, berberine (IC50: 7.67 ± 0.36 µg/mL) and (-)-corydalmine (IC50: 7.78 ± 0.38 µg/mL) displayed a stronger inhibition than that of galanthamine (IC50: 12.02 ± 0.25 µg/mL). The mutagenic activity was shown for β-allocryptopine, (+)- and (-)-bicuculline, (±)-corydalidzine, (-)-corydalmine, (+)-cularicine, (-)-fumarophycine, (-)-norjuziphine, (-)-ophiocarpine-N-oxide, (+)-scoulerine, (-)-sinactine, and (-)-stylopine by means of in silico experiments. The results obtained by molecular docking simulations of berberine, palmatine, and (-)-corydalmine suggested that the estimated free ligand-binding energies of these compounds inside the binding domains of their targets are reasonable to make them capable of establishing strong polar and nonpolar bonds with the atoms of the active site amino acids.

CONCLUSION

Our findings revealed that berberine, palmatin, and (-)-corydalmine stand out as the most promising isoquinoline alkaloids in terms of ChE inhibition. Among them, berberine has displayed a robust dual inhibition against both ChEs and could be evaluated further as a lead compound for AD.

Evaluation of Anti-Alzheimer Activity of Synthetic Coumarins by Combination of in Vitro and in Silico Approaches

Series of synthetic coumarin derivatives (1-16) were tested against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), two enzymes linked to the pathology of Alzheimer's disease (AD). Compound 16 was the most active AChE inhibitor with IC₅₀ 32.23±2.91 μM, while the reference (galantamine) had IC₅₀ =1.85±0.12 μM. Compounds 9 (IC₅₀ 75.14±1.82 μM), 13 (IC₅₀ =16.14±0.43 μM), were determined to be stronger BChE inhibitors than the reference galantamine (IC₅₀ =93.53±2.23 μM). The IC₅₀ value of compound 16 for BChE inhibition (IC₅₀ =126.56±11.96 μM) was slightly higher than galantamine. The atomic interactions between the ligands and the key amino acids inside the binding cavities were simulated to determine their ligand-binding positions and free energies. The three inhibitory coumarins (9, 13, 16) were next tested for their effects on the genes associated with AD using human neuroblastoma (SH-SY5Y) cell lines. Our data indicate that they could be considered for further evaluation as new anti-Alzheimer drug candidates.

Exploiting the Propagation of Constrained Variables for Enhanced HDX-MS Data Optimization

Nonlinear programming has found useful applications in protein biophysics to help understand the microscopic exchange kinetics of data obtained using hydrogen-deuterium exchange mass spectrometry (HDX-MS). Finding a microscopic kinetic solution for HDX-MS data provides a window into local protein stability and energetics allowing them to be quantified and understood. Optimization of HDX-MS data is a significant challenge, however, due to the requirement to solve a large number of variables simultaneously with exceptionally large variable bounds. Modeled rates are frequently uncertain with an explicate dependency on the initial guess values. In order to enhance the search for a minimum solution in HDX-MS optimization, the ability of selected constrained variables to propagate throughout the data is considered. We reveal that locally bound constrained optimization induces a global effect on all variables. The global response to local constraints is large and surprisingly long-range, but the outcome is unpredictable, unexpectedly decreasing the overall accuracy of certain data sets depending on the stringency of the constraints. Utilizing previously described in-house validation criteria based on covariance matrices, a method is described that is able to accurately determine whether constraints benefit or impair the optimization of HDX-MS data. From this, we establish a new two-stage method for our online optimizer HDXmodeller that can effectively leverage locally bound variables to enhance HDX-MS data modeling.

Characterization and Management of Noise in HDX-MS Data Modeling

Quantification of hydrogen deuterium exchange (HDX) kinetics can provide information on the stability of individual amino acids in proteins by finding the degree to which the local backbone environment corresponds to that of a random coil. When characterized by mass spectrometry, extraction of HDX kinetics is not possible because different residue exchange rates become merged depending on the peptides that are formed during proteolytic digestion. We have recently developed an advanced programming tool called HDXmodeller, which enables the exchange rates of individual amino acids to be understood by optimization of low-resolution HDX-mass spectrometry (MS) data. HDXmodeller is also uniquely able to appraise each optimization and quantify the accuracy of modeled exchange rates ab initio using a novel autovalidation method based on a covariance matrix. Here, we address the noise-handling capabilities of HDXmodeller and demonstrate the effectiveness of the algorithm on self-inconsistent datasets. Reference intervals for experimental HDX-MS data are also derived, and this information is presented in an updated online workflow for HDXmodeller, allowing users to evaluate the consistency of their data. The development of a modified version of HDXmodeller is also discussed with enhanced noise-handling capability brought about through loss function optimization. Changes in optimizer accuracy with different loss functions are also demonstrated along with the effectiveness of HDXmodeller to select the most effective optimizer for different data using currently embedded autovalidation criteria.

Synthesis, bioactivity and binding energy calculations of novel 3-ethoxysalicylaldehyde based thiosemicarbazone derivatives

In recent decade, the entrance of α-N-heterocyclic thiosemicarbazones derivates (Triapne, COTI-2 and DpC) in clinical trials for cancer and HIV-1 has vastly increased the interests of medicinal chemists towards this class of organic compounds. In the given study, a series of eighteen new (3a-r) 3-ethoxy salicylaldehyde-based thiosemicarbazones (TSC), bearing aryl and cycloalkyl substituents, were synthesized and assayed for their pharmacological potential against carbonic anhydrases (hCA I and hCA II), cholinesterases (AChE and BChE) and α-glycosidase. The hCA I isoform was inhibited by these novel 3-ethoxysalicylaldehyde thiosemicarbazone derivatives (3a-r) in low nanomolar levels, the Ki of which differed between 144.18 ± 26.74 and 454.92 ± 48.32 nM. Against the physiologically dominant isoform hCA II, the novel compounds demonstrated K_is varying from 110.54 ± 14.05 to 444.12 ± 36.08 nM. Also, these novel derivatives (3a-r) effectively inhibited AChE, with Ki values in the range of 385.38 ± 45.03 to 983.04 ± 104.64 nM. For BChE was obtained with Ki values in the range of 400.21 ± 35.68 to 1003.02 ± 154.27 nM. For α-glycosidase the most effective Ki values of 3l, 3n, and 3q were with Ki values of 12.85 ± 1.05, 16.03 ± 2.84, and 19.16 ± 2.66 nM, respectively. Moreover, the synthesized TCSs were simulated using force field methods whereas the binding energies of the selected compounds were estimated using MM-GBSA method. The findings indicate the present novel 3-ethoxy salicylaldehyde-based thiosemicarbazones to be excellent hits for pharmaceutical applications.

Synthesis, anticholinesterase activity and molecular modeling studies of novel carvacrol-substituted amide derivatives

In the present study, 23 novel carvacrol derivatives involving the amide moiety as a linker between the alkyl chains and/or the heterocycle nucleus were synthesized and tested in vitro as acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) inhibitors. 2-(5-Isopropyl-2-methylphenoxy)-N-(quinolin-8-yl)acetamide (5v) revealed the highest inhibition properties against AChE and BuChE with the IC₅₀ values of 1.93 and 0.05 µM, respectively. The blood-brain barrier (BBB) permeability of the potent inhibitor (5v) was also assessed by the widely used parallel artificial membrane permeability assay (PAMPA-BBB). The results showed that 5v is capable of crossing the BBB. Pharmacokinetic and toxicity profiles of the studied molecule predictions were investigated by MetaCore/MetaDrug comprehensive systems biology analysis suite. Bioactive conformations of the synthesized molecules, their predicted binding energies as well as structural and dynamical profiles of molecules at the binding pockets of AChE and BuChE targets were also investigated using different docking algorithms and molecular dynamics (MD) simulations.Communicated by Ramaswamy H. Sarma.

Combined molecular modeling and cholinesterase inhibition studies on some natural and semisynthetic O-alkylcoumarin derivatives

Coumarins of synthetic or natural origins are an important chemical class exerting diverse pharmacological activities. In the present study, 26 novel O-alkylcoumarin derivatives were synthesized and have been tested at 100 µM for their in vitro inhibitory potential against acetylcholinesterase (AChE) and butyrlcholinesterase (BChE) targets which are the key enzymes playing role in the pathogenesis of Alzheimer's disease. Among the tested coumarins, none of them could inhibit AChE, whereas 12 of them exerted a marked and selective inhibition against BChE as compared to the reference (galanthamine, IC₅₀ = 46.58 ± 0.91 µM). In fact, 10 of the active coumarins showed higher inhibition (IC₅₀ = 7.01 ± 0.28 µM - 43.31 ± 3.63 µM) than that of galanthamine. The most active ones were revealed to be 7-styryloxycoumarin (IC₅₀ = 7.01 ± 0.28 µM) and 7-isopentenyloxy-4-methylcoumarin (IC₅₀ = 8.18 ± 0.74 µM). In addition to the in vitro tests, MetaCore/MetaDrug binary QSAR models and docking simulations were applied to evaluate the active compounds by ligand-based and target-driven approaches. The predicted pharmacokinetic profiles of the compounds suggested that the compounds reveal lipophilic character and permeate blood brain barrier (BBB) and the ADME models predict higher human serum protein binding percentages (>50%) for the compounds. The calculated docking scores indicated that the coumarins showing remarkable BChE inhibition possessed favorable free binding energies in interacting with the ligand-binding domain of the target. Therefore, our results disclose that O-alkylcoumarins are promising selective inhibitors of cholinesterase enzymes, particularly BChE in our case, which definitely deserve further studies.

Investigating the molecular mechanism of staphylococcal DNA gyrase inhibitors: A combined ligand-based and structure-based resources pipeline

Appropriate therapeutic solutions against Staphylococcal infections are currently limited. To work out the complex task of challenging drug resistance in Staphylococcus aureus, new compounds with novel modes of action are required. In this study, we performed target-driven virtual screening to filter exhaustive phytochemical libraries that can inhibit the activity of S. aureus DNA Gyrase B (Gyr B). Three top-ranked hit molecules (Mangostenone E, Candenatenin A and 2,4,4'-trihydroxydihydrochalcone) were identified from comprehensive molecular docking studies based on their strong spatial affinity with key catalytic residues of the binding pocket of DNA GyrB, especially with the well-known crucial residue Asp81. Molecular dynamics (MD) simulations were performed for these identified hit molecules for better understanding of their dynamical and structural profiles throughout the MD simulations. These compounds can be explored as future lead optimization molecules to discover a new class of antibiotics against resistant Staphylococcus aureus strains.

Structural Investigation of the Dopamine-2 Receptor Agonist Bromocriptine Binding to Dimeric D2HighR and D2LowR States

The active (D2^HighR) and inactive (D2^LowR) states of dimeric dopamine D2 receptor (D2R) models were investigated to clarify the binding mechanisms of the dopamine agonist bromocriptine, using Molecular Dynamics (MD) simulation. The aim of this comprehensive study was to investigate the critical effects of bromocriptine binding on each distinct receptor conformation. The different binding modes of the bromocriptine ligand in the active and inactive states have a significant effect on the conformational changes of the receptor. Based on the MM/GBSA approach, the calculated binding enthalpies of bromocriptine demonstrated selectivity toward the D2^HighR active state. There is good agreement between the calculated and experimentally measured D2^HighR selectivity. In the ligand-binding site, the key amino acids identified for D2^HighR were Asp114(3.32) and Glu95(2.65), and for D2^LowR, it was Ser193(5.42). Moreover, analysis of replicate MD trajectories demonstrated that the bromocriptine structure was more rigid at the D2^HighR state and more flexible at the D2^LowR state. However, the side chains of the ligand-receptor complex of D2^HighR showed larger variations relative to the corresponding regions of D2^LowR. The present study is part of an ongoing research program to study D2R conformational changes during ligand activation and to evaluate the conformational state selectivity for ligand binding.

Oligomerization and cooperativity in GPCRs from the perspective of the angiotensin AT1 and dopamine D2 receptors

G Protein-Coupled Receptors (GPCRs) can form homo- and heterodimers or constitute higher oligomeric clusters with other heptahelical GPCRs. In this article, multiscale molecular modeling approaches as well as experimental techniques which are used to study oligomerization of GPCRs are reviewed. In particular, the effect of dimerization/oligomerization to the ligand binding affinity of individual protomers and also on the efficacy of the oligomer are discussed by including diverse examples from the literature. In addition, possible allosteric effects that may emerge upon interaction of GPCRs with membrane components, like cholesterol, is also discussed. Investigation of these above-mentioned interactions may greatly contribute to the candidate molecule screening studies and development of novel therapeutics with fewer adverse effects.

Molecular modeling and in vitro approaches towards cholinesterase inhibitory effect of some natural xanthohumol, naringenin, and acyl phloroglucinol derivatives

BACKGROUND

Many natural products, particularly phenolic compounds, have been reported to have a strong inhibition against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), the key enzymes in the pathology of Alzheimer's disease (AD).

HYPOTHESIS

Therefore, we hypothesized that some xanthahumol, naringenin, and acyl phloroglucinol derivatives (1-14) isolated from Humulus lupulus L. (hops) may have an inhibitory potential against AChE and BChE.

METHODS

Inhibitory potential of compounds 1-14 were tested against AChE and BChE using ELISA microtiter assay. Different molecular docking simulations, including IFD and GOLD protocols, were implemented to verify the interactions between the ligands and the active site amino acids and also their binding energies inside the catalytic crevices of AChE and BChE. ADME/Tox analysis were used to determine pharmacological activities of the compounds.

RESULTS

Among them, 3‑hydroxy‑xanthohumol (IC₅₀ = 51.25 ± 0.88 µM) and xanthohumol (IC₅₀ = 71.34 ± 2.09 µM), displayed a moderate AChE inhibition in comparison to that of the reference (galanthamine, IC₅₀ = 2.52 ± 0.15 µM). In addition to 3‑hydroxy‑xanthohumol (IC₅₀ = 63.07 ± 3.76 µM) and xanthohumol (IC₅₀ = 32.67 ± 2.82 µM), 8-prenylnaringenin (IC₅₀ = 86.58 ± 3.74 µM) also showed micromolar-range inhibition against BChE (galanthamine, IC₅₀ = 46.58 ± 0.91 µM). Rest of the compounds were found to be either inactive or having inhibition below 50%. Prediction of pharmacokinetic studies suggested that all the ligands revealed acceptable drug-like profiles. Docking simulations demonstrate not only the prediction of ligand binding energies of the compounds inside the catalytic domains of the targets, but also highlight the critical amino acids contributing to stabilizations of the ligands.

CONCLUSION

Our findings revealed that xanthohumol in particular could be considered as lead molecule to explore new cholinesterase inhibitors for AD.

Discovery of selective dengue virus inhibitors using combination of molecular fingerprint-based virtual screening protocols, structure-based pharmacophore model development, molecular dynamics simulations and in vitro studies

Dengue virus is a major issue of tropical and sub-tropical regions. The proliferation of virus results in immense number of deaths each year because of unavailability of on-shelf drugs. This issue necessitates the design of novel anti-Dengue drugs. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of Dengue virus. Keeping in mind the worsening situation regarding Dengue virus, approximately eighteen million drug-like compounds from the ZINC small molecule database have been screened against Nonstructural Protein 3 (NS3) previously by our group. In this study, in order to investigate the effect of extended time of molecular dynamics (MD) simulations on structural and dynamical profiles of used complexes, simulation run time is increased from 50-ns to 100-ns for the each system. In addition, a well-known Dengue virus inhibitor (MB21) from literature is used as reference structure (positive control) to compare the proposed molecules. Post-processing MD analyses including Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculations were conducted to predict binding free energies of inhibitors from derived trajectory frames of MD simulations. Identified compounds are further directed to Quantum-Polarized Ligand Docking (QPLD), molecular fingerprint-based virtual screening of another small molecule database (Otava Drug Like small molecule database), and Structure-based Pharmacophore Modeling (E-Pharmacophore). Finally, cell proliferation and cytotoxicity tests as well as pre- and post-treatment on HUH7 cells infected with DENV2 NGC strain are applied for four identified hit molecules (ZINC36681949, ZINC44921800, ZINC95518765 and ZINC39500661) to check whether these drugs inhibit DENV2 from entry and/or exit pathways. Based on cell-based Dengue quantification assays, there is no effect seen on pre-treatment of cells with these compounds indicating that the early infection processes of virus is not affected. In contrast, the post-treatment of cells with these compounds after Dengue virus infection has resulted in a significant 1 log PFU/ml reduction of the virus infectious titre.

Structure-based design of hERG-neutral antihypertensive oxazalone and imidazolone derivatives

Angiotensin II receptor type 1 (AT1) antagonists are the most recent drug class against hypertension. Recently first crystal structure of AT1 receptor is deposited to the protein data bank (PDB ID: 4YAY). In this work, several molecular screening methods such as molecular docking and de novo design studies were performed and it is found that oxazolone and imidazolone derivatives reveal similar/better interaction energy profiles compared to the FDA approved sartan molecules at the binding site of the AT1 receptor. A database consisting of 3500-fragments were used to enumerate de novo designed imidazolone and oxazolone derivatives and hereby more than 50000 novel small molecules were generated. These derivatives were then used in high throughput virtual screening simulations (Glide/HTVS) to find potent hit molecules. In addition, virtual screening of around 18 million small drug-like compounds from ZINC database were screened at the binding pocket of the AT1 receptor via Glide/HTVS method. Filtered structures were then used in more sophisticated molecular docking simulations protocols (i.e., Glide/SP; Glide/XP; Glide/IFD; Glide/QPLD, and GOLD). However, the K⁺ ion channel/drug interactions should also be considered in studies implemented in molecular level against their cardiovascular risks. Thus, selected compounds with high docking scores via all diverse docking algorithms are also screened at the pore domain regions of human ether-a-go-go-related gene (hERG1) K⁺ channel to remove the high affinity hERG1 blocking compounds. High docking scored compounds at the AT1 with low hERG1 affinity is considered for long molecular dynamics (MD) simulations. Post-processing analysis of MD simulations assisted for better understanding of molecular mechanism of studied compounds at the binding cavity of AT1 receptor. Results of this study can be useful for designing of novel and safe AT1 inhibitors.

Synthesis, biological activity and multiscale molecular modeling studies for coumaryl-carboxamide derivatives as selective carbonic anhydrase IX inhibitors

New coumaryl-carboxamide derivatives with the thiourea moiety as a linker between the alkyl chains and/or the heterocycle nucleus were synthesized and their inhibitory activity against the human carbonic anhydrase (hCA) isoforms hCA I, II, VII and IX were evaluated. While the hCA I, II and VII isoforms were not inhibited by the investigated compounds, the tumour-associated isoform hCA IX was inhibited in the high nanomolar range. 2-Oxo-N-((2-(pyrrolidin-1-yl)ethyl)carbamothioyl)-2H-chromene-3-carboxamide (e11) exhibited a selective inhibitory action against hCA IX with the K_i of 107.9 nM. In order to better understand the inhibitory profiles of studied molecules, multiscale molecular modeling approaches were used. Different molecular docking algorithms were used to investigate binding poses and predicted binding energies of studied compounds at the active sites of the CA I, II, VII and IX isoforms.

Effects of propolis, caffeic acid phenethyl ester, and pollen on renal injury in hypertensive rat: An experimental and theoretical approach

The objective of this study was to evaluate the antioxidant effects of propolis, caffeic acid phenethyl ester (CAPE; active compound in propolis), and pollen on biochemical oxidative stress biomarkers in rat kidney tissue inhibited by N^ω -nitro-L-arginine methyl ester (L-NAME). The biomarkers evaluated were paraoxonase (PON1), oxidative stress index (OSI), total antioxidant status (TAS), total oxidant status (TOS), asymmetric dimethylarginine (ADMA), and nuclear factor kappa B (NF-κB). TAS levels and PON1 activity were significantly decreased in kidney tissue samples in the L-NAME-treated group (P < 0.05). The levels of TAS and PONI were higher in the L-NAME plus propolis, CAPE, and pollen groups compared with the L-NAME-treated group. TOS, ADMA, and NF-κB levels were significantly increased in the kidney tissue samples of the L-NAME-treated group (P < 0.05). However, these parameters were significantly lower in the L-NAME plus propolis, CAPE, and pollen groups (P < 0.05) compared with rats administered L-NAME alone (P < 0.05). Furthermore, the binding energy of CAPE within catalytic domain of glutathione reductase (GR) enzyme as well as its inhibitory mechanism was determined using molecular modeling approaches. In conclusion, experimental and theoretical data suggested that oxidative alterations occurring in the kidney tissue of chronic hypertensive rats may be prevented via active compound of propolis, CAPE administration.

Discovery of selective dengue virus inhibitors using combination of molecular fingerprint-based virtual screening protocols, structure-based pharmacophore model development, molecular dynamics simulations and in vitro studies

Dengue virus is a major issue of tropical and sub-tropical regions. The proliferation of virus results in immense number of deaths each year because of unavailability of on-shelf drugs. This issue necessitates the design of novel anti-Dengue drugs. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of Dengue virus. Keeping in mind the worsening situation regarding Dengue virus, approximately eighteen million drug-like compounds from the ZINC small molecule database have been screened against Nonstructural Protein 3 (NS3) previously by our group. In this study, in order to investigate the effect of extended time of molecular dynamics (MD) simulations on structural and dynamical profiles of used complexes, simulation run time is increased from 50-ns to 100-ns for the each system. In addition, a well-known Dengue virus inhibitor (MB21) from literature is used as reference structure (positive control) to compare the proposed molecules. Post-processing MD analyses including Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculations were conducted to predict binding free energies of inhibitors from derived trajectory frames of MD simulations. Identified compounds are further directed to Quantum-Polarized Ligand Docking (QPLD), molecular fingerprint-based virtual screening of another small molecule database (Otava Drug Like small molecule database), and Structure-based Pharmacophore Modeling (E-Pharmacophore). Finally, cell proliferation and cytotoxicity tests as well as pre- and post-treatment on HUH7 cells infected with DENV2 NGC strain are applied for four identified hit molecules (ZINC36681949, ZINC44921800, ZINC95518765 and ZINC39500661) to check whether these drugs inhibit DENV2 from entry and/or exit pathways. Based on cell-based Dengue quantification assays, there is no effect seen on pre-treatment of cells with these compounds indicating that the early infection processes of virus is not affected. In contrast, the post-treatment of cells with these compounds after Dengue virus infection has resulted in a significant 1logPFU/ml reduction of the virus infectious titre.

Analysis of the Glutamate Agonist LY404,039 Binding to Nonstatic Dopamine Receptor D2 Dimer Structures and Consensus Docking

Dopamine receptor D2 (D2R) plays an important role in the human central nervous system and is a focal target of antipsychotic agents. The D2^HighR and D2^LowR dimeric models previously developed by our group are used to investigate the prediction of binding affinity of the LY404,039 ligand and its binding mechanism within the catalytic domain. The computational data obtained using molecular dynamics simulations fit well with the experimental results. The calculated binding affinities of LY404,039 using MM/PBSA for the D2^HighR and D2^LowR targets were -12.04 and -9.11 kcal/mol, respectively. The experimental results suggest that LY404,039 binds to D2^HighR and D2^LowR with binding affinities (K_i) of 8.2 and 1640 nM, respectively. The high binding affinity of LY404,039 in terms of binding to [³H]domperidone was inhibited by the presence of a guanine nucleotide, indicating an agonist action of the drug at D2^HighR. The interaction analysis demonstrated that while Asp114 was among the most critical amino acids for D2^HighR binding, residues Ser193 and Ser197 were significantly more important within the binding cavity of D2^LowR. Molecular modeling analyses are extended to ensemble docking as well as structure-based pharmacophore model (E-pharmacophore) development using the bioactive conformation of LY404,039 at the binding pocket as a template and screening of small-molecule databases with derived pharmacophore models.

Biological Insights of the Dopaminergic Stabilizer ACR16 at the Binding Pocket of Dopamine D2 Receptor

The dopamine D2 receptor (D2R) plays an important part in the human central nervous system and it is considered to be a focal target of antipsychotic agents. It is structurally modeled in active and inactive states, in which homodimerization reaction of the D2R monomers is also applied. The ASP2314 (also known as ACR16) ligand, a D2R stabilizer, is used in tests to evaluate how dimerization and conformational changes may alter the ligand binding space and to provide information on alterations in inhibitory mechanisms upon activation. The administration of the D2R agonist ligand ACR16 [³H](+)-4-propyl-3,4,4a,5,6,10b-hexahydro-2H-naphtho[1,2-b][1,4]oxazin-9-ol ((+)PHNO) revealed K_i values of 32 nM for the D2^highR and 52 μM for the D2^lowR. The calculated binding affinities of ACR16 with post processing molecular dynamics (MD) simulations analyses using MM/PBSA for the monomeric and homodimeric forms of the D2^highR were -9.46 and -8.39 kcal/mol, respectively. The data suggests that the dimerization of the D2R leads negative cooperativity for ACR16 binding. The dimerization reaction of the D2^highR is energetically favorable by -22.95 kcal/mol. The dimerization reaction structurally and thermodynamically stabilizes the D2^highR conformation, which may be due to the intermolecular forces formed between the TM4 of each monomer, and the result strongly demonstrates dimerization essential for activation of the D2R.

Biochemical changes induced by grape seed extract and low level laser therapy administration during intraoral wound healing in rat liver: an experimental and in silico study

In the present study, the changes that occur in rat liver tissue as a result of the use of grape seed extract (GSE) and low level laser therapy (LLLT) in intraoral wound (IW) healing are analyzed using biochemical parameters. Diode laser application groups received 8 J/cm² dose LLLT once a day for 4 days (810 nm wavelength, continuous mode, 0.25 W, 9 s). As a result of the biological parameter analysis, it was determined that the oxidative damage caused by the IWs and recovery period on 7th and 14th days could be substantially removed with GSE applications that have antioxidant capacity especially in rat liver tissue. In addition, the active compound of grape seed, catechin is studied in the active site of glycogen synthase kinase 3 (GSK3) target using molecular modeling approaches. Post-processing molecular dynamics (MD) results for catechin is compared with a standard GSK3 inhibitor. MD simulations assisted for better understanding of inhibition mechanism and the crucial amino acids contributing in the ligand binding. These results along with a through free energy analysis of ligands using sophisticated simulations methods are quite striking and it suggests a greater future role for simulation in deciphering complex patterns of molecular mechanism in combination with methods for understanding drug-receptor interactions.

First universal pharmacophore model for hERG1 K+ channel activators: acthER

The intra-cavitary drug blockade of hERG1 channel has been extensively studied, both experimentally and theoretically. Structurally diverse ligands inadvertently block the hERG1 K⁺ channel currents lead to drug induced Long QT Syndrome (LQTS). Accordingly, designing either hERG1 channel openers or current activators, with the potential to target other binding pockets of the channel, has been introduced as a viable approach in modern anti-arrhythmia drug development. However, reports and investigations on the molecular mechanisms underlying activators binding to the hERG1 channel remain sparse and the overall molecular design principles are largely unknown. Most of the hERG1 activators were discovered during mandatory screening for hERG1 blockade. To fill this apparent deficit, the first universal pharmacophore model for hERG1 K⁺ channel activators was developed using PHASE. 3D structures of 18 hERG1 K⁺ channel activators and their corresponding measured binding affinity values were used in the development of pharmacophore models. These compounds spanned a range of structurally different chemotypes with moderate variation in binding affinity. A five sites AAHRR (A, hydrogen-bond accepting, H, hydrophobic, R, aromatic) pharmacophore model has shown reasonable high statistical results compared to the other developed more than 1000 hypotheses. This model was used to construct steric and electrostatic contour maps. The predictive power of the model was tested with 3 external test set compounds as true unknowns. Finally, the pharmacophore model was combined with the previously developed receptor-based model of hERG1 K⁺ channel to develop and screen novel activators. The results are quite striking and it suggests a greater future role for pharmacophore modeling and virtual drug screening simulations in deciphering complex patterns of molecular mechanisms of hERG1 channel openers at the target sites. The developed model is available upon request and it may serve as basis for the synthesis of novel therapeutic hERG1 activators.

Evaluation of the potential toxicity of unmodified and modified cyclodextrins on murine blood-brain barrier endothelial cells

In this study, we investigated the cytotoxic effects of unmodified α-cyclodextrin (α-CD) and modified cyclodextrins, including trimethyl-β-cyclodextrin (TRIMEB) and hydroxypropyl-β-cyclodextrin (HPβCD), on immortalized murine microvascular endothelial (cEND) cells of the blood-brain barrier (BBB). A CellTiter-Glo viability test, performed on the cEND cells showed significant differences among the different cyclodextrins. After 24 hr of incubation, TRIMEB was the most cytotoxic, and HPβCD was non-toxic. α-CD and TRIMEB exhibited greater cytotoxicity in the Dulbecco's modified Eagle's medium than in heat-inactivated human serum indicating protective properties of the human serum. The predicted dynamic toxicity profiles (Td) for α-CD and TRIMEB indicated higher cytotoxicity for these cyclodextrins compared to the reference compound (dimethylsulfoxide). Molecular dynamics simulation of cholesterol binding to the CDs suggested that not just cholesterol but phospholipids extraction might be involved in the cytotoxicity. Overall, the results demonstrate that HPβCD has the potential to be used as a candidate for drug delivery vector development and signify a correlation between the in vitro cytotoxic effect and cholesterol binding of cyclodextrins.

Discovery of Klotho peptide antagonists against Wnt3 and Wnt3a target proteins using combination of protein engineering, protein-protein docking, peptide docking and molecular dynamics simulations

The Klotho is known as lifespan enhancing protein involved in antagonizing the effect of Wnt proteins. Wnt proteins are stem cell regulators, and uninterrupted exposure of Wnt proteins to the cell can cause stem and progenitor cell senescence, which may lead to aging. Keeping in mind the importance of Klotho in Wnt signaling, in silico approaches have been applied to study the important interactions between Klotho and Wnt3 and Wnt3a (wingless-type mouse mammary tumor virus (MMTV) integration site family members 3 and 3a). The main aim of the study is to identify important residues of the Klotho that help in designing peptides which can act as Wnt antagonists. For this aim, a protein engineering study is performed for Klotho, Wnt3 and Wnt3a. During the theoretical analysis of homology models, unexpected role of number of disulfide bonds and secondary structure elements has been witnessed in case of Wnt3 and Wnt3a proteins. Different in silico experiments were carried out to observe the effect of correct number of disulfide bonds on 3D protein models. For this aim, total of 10 molecular dynamics (MD) simulations were carried out for each system. Based on the protein–protein docking simulations of selected protein models of Klotho with Wnt3 and Wnt3a, different peptides derived from Klotho have been designed. Wnt3 and Wnt3a proteins have three important domains: Index finger, N-terminal domain and a patch of ∼10 residues on the solvent exposed surface of palm domain. Protein–peptide docking of designed peptides of Klotho against three important domains of palmitoylated Wnt3 and Wnt3a yields encouraging results and leads better understanding of the Wnt protein inhibition by proposed Klotho peptides. Further in vitro studies can be carried out to verify effects of novel designed peptides as Wnt antagonists.

Characterization, in Vivo Evaluation, and Molecular Modeling of Different Propofol-Cyclodextrin Complexes To Assess Their Drug Delivery Potential at the Blood-Brain Barrier Level

In this study, we investigated the ability of the general anesthetic propofol (PR) to form inclusion complexes with modified β-cyclodextrins, including sulfobutylether-β-cyclodextrin (SBEβCD) and hydroxypropyl-β-cyclodextrin (HPβCD). The PR/SBEβCD and PR/HPβCD complexes were prepared and characterized, and the blood-brain barrier (BBB) permeation potential of the formulated PR was examined in vivo for the purpose of controlled drug delivery. The PR/SBEβCD complex was found to be more stable in solution with a minimal degradation constant of 0.25 h^-1, a t_1/2 of 2.82 h, and a K_c of 5.19 × 10³ M^-1 and revealed higher BBB permeability rates compared with the reference substance (PR-LIPURO) considering the calculated brain-to-blood concentration ratio (logBB) values. Additionally, the diminished PR binding affinity to SBEβCD was confirmed in molecular dynamics simulations by a maximal Gibbs free energy of binding (ΔG_bind = -18.44 kcal·mol^-1), indicating the more rapid PR/SBEβCD dissociation. Overall, the results demonstrated that SBEβCD has the potential to be used as a prospective candidate for drug delivery vector development to improve the pharmacokinetic and pharmacodynamic properties of general anesthetic agents at the BBB level.

The signaling pathway of dopamine D2 receptor (D2R) activation using normal mode analysis (NMA) and the construction of pharmacophore models for D2R ligands

G-protein-coupled receptors (GPCRs) are targets of more than 30% of marketed drugs. Investigation on the GPCRs may shed light on upcoming drug design studies. In the present study, we performed a combination of receptor- and ligand-based analysis targeting the dopamine D2 receptor (D2R). The signaling pathway of D2R activation and the construction of universal pharmacophore models for D2R ligands were also studied. The key amino acids, which contributed to the regular activation of the D2R, were in detail investigated by means of normal mode analysis (NMA). A derived cross-correlation matrix provided us an understanding of the degree of pair residue correlations. Although negative correlations were not observed in the case of the inactive D2R state, a high degree of correlation appeared between the residues in the active state. NMA results showed that the cytoplasmic side of the TM5 plays a significant role in promoting of residue-residue correlations in the active state of D2R. Tracing motions of the amino acids Arg219, Arg220, Val223, Asn224, Lys226, and Ser228 in the position of the TM5 are found to be critical in signal transduction. Complementing the receptor-based modeling, ligand-based modeling was also performed using known D2R ligands. The top-scored pharmacophore models were found as 5-sited (AADPR.671, AADRR.1398, AAPRR.3900, and ADHRR.2864) hypotheses from PHASE modeling from a pool consisting of more than 100 initial candidates. The constructed models using 38 D2R ligands (in the training set) were validated with 15 additional test set compounds. The resulting model correctly predicted the pIC₅₀ values of an additional test set compounds as true unknowns.

Virtual screening of eighteen million compounds against dengue virus: Combined molecular docking and molecular dynamics simulations study

Dengue virus is a major issue of tropical and sub-tropical regions. Dengue virus has been the cause behind the major alarming epidemics in the history with mass causalities from the decades. Unavailability of on-shelf drugs for the prevention of further proliferation of virus inside the human body results in immense number of deaths each year. This issue necessitates the design of novel anti-dengue drug. The protease enzyme pathway is the critical target for drug design due to its significance in the replication, survival and other cellular activities of dengue virus. Therefore, approximately eighteen million compounds from the ZINC database have been virtually screened against nonstructural protein 3 (NS3). The incremental construction algorithm of Glide docking program has been used with its features high throughput virtual screening (HTVS), standard precision (SP), extra precision (XP) and in combination of Prime module, induced fit docking (IFD) approach has also been applied. Five top-ranked compounds were then selected from the IFD results with better predicted binding energies with the catalytic triad residues (His51, Asp75, and Ser135) that may act as potential inhibitors for the underlying target protease enzyme. The top-ranked compounds ZINC95518765, ZINC44921800, ZINC71917414, ZINC39500661, ZINC36681949 have shown the predicted binding energies of -7.55, -7.36, -8.04, -8.41, -9.18kcal/mol, respectively, forming binding interactions with three catalytically important amino acids. Top-docking poses of compounds are then used in molecular dynamics (MD) simulations. In computational studies, our proposed compounds confirm promising results against all the four serotypes of dengue virus, strengthening the opportunity of these compounds to work as potential on-shelf drugs against dengue virus. Further experimentation on the proposed compounds can result in development of strong inhibitors.

Binding Interactions of Dopamine and Apomorphine in D2High and D2Low States of Human Dopamine D2 Receptor Using Computational and Experimental Techniques

We have recently reported G-protein coupled receptor (GPCR) model structures for the active and inactive states of the human dopamine D2 receptor (D2R) using adrenergic crystal structures as templates. Since the therapeutic concentrations of dopamine agonists that suppress the release of prolactin are the same as those that act at the high-affinity state of the D2 receptor (D2High), D2High in the anterior pituitary gland is considered to be the functional state of the receptor. In addition, the therapeutic concentrations of anti-Parkinson drugs are also related to the dissociation constants in the D2High form of the receptor. The discrimination between the high- and low-affinity (D2Low) components of the D2R is not obvious and requires advanced computer-assisted structural biology investigations. Therefore, in this work, the derived D2High and D2Low receptor models (GPCR monomer and dimer three-dimensional structures) are used as drug-binding targets to investigate binding interactions of dopamine and apomorphine. The study reveals a match between the experimental dissociation constants of dopamine and apomorphine at their high- and low-affinity sites of the D2 receptor in monomer and dimer and their calculated dissociation constants. The allosteric receptor-receptor interaction for dopamine D2R dimer is associated with the accessibility of adjacent residues of transmembrane region 4. The measured negative cooperativity between agonist ligand at dopamine D2 receptor is also correctly predicted using the D2R homodimerization model.

Mutated form (G52E) of inactive diphtheria toxin CRM197: molecular simulations clearly display effect of the mutation to NAD binding

Mutated form (G52E) of diphtheria toxin (DT) CRM197 is an inactive and nontoxic enzyme. Here, we provided a molecular insight using comparative molecular dynamics (MD) simulations to clarify the influence of a single point mutation on overall protein and active-site loop. Post-processing MD analysis (i.e. stability, principal component analysis, hydrogen-bond occupancy, etc.) is carried out on both wild and mutated targets to investigate and to better understand the mechanistic differences of structural and dynamical properties on an atomic scale especially at nicotinamide adenine dinucleotide (NAD) binding site when a single mutation (G52E) happens at the DT. In addition, a docking simulation is performed for wild and mutated forms. The docking scoring analysis and docking poses results revealed that mutant form is not able to properly accommodate the NAD molecule.

Investigation of Inhibition Mechanism of Chemokine Receptor CCR5 by Micro-second Molecular Dynamics Simulations

Chemokine receptor 5 (CCR5) belongs to G protein coupled receptors (GPCRs) and plays an important role in treatment of human immunodeficiency virus (HIV) infection since HIV uses CCR5 protein as a co-receptor. Recently, the crystal structure of CCR5-bound complex with an approved anti-retroviral drug (maroviroc) was resolved. During the crystallization procedure, amino acid residues (i.e., Cys224, Arg225, Asn226 and Glu227) at the third intra-cellular loop were replaced by the rubredoxin for stability reasons. In the current study, we aimed to understand the impact of the incorporated rubredoxin on the conformations of TM domains of the target protein. For this reason, rubredoxin was deleted from the crystal structure and the missing amino acids were engineered. The resultant structure was subjected to long (μs) molecular dynamics (MD) simulations to shed light into the inhibitory mechanism. The derived model structure displayed a significant deviation in the cytoplasmic domain of TM5 and IC3 in the absence of rubredoxin. The principal component analyses (PCA) and MD trajectory analyses revealed important structural and dynamical differences at apo and holo forms of the CCR5.

Discovering novel carbonic anhydrase type IX (CA IX) inhibitors from seven million compounds using virtual screening and in vitro analysis

Carbonic anhydrase type IX (CA IX) enzyme is mostly over expressed in different cancer cell lines and tumor tissues. Potent CA IX inhibitors can be effective for adjusting the pH imbalance in tumor cells. In the present work, we represented the successful application of high throughput virtual screening (HTVS) of large dataset from ZINC database included of ∼7 million compounds to discover novel inhibitors of CA IX. HTVS and molecular docking were performed using consequence Glide/standard precision (SP), extra precision (XP) and induced fit docking (IFD) molecular docking protocols. For each compound, docking code calculates a set of low-energy poses and then exhaustively scans the binding pocket of the target with small compounds. Novel CA IX inhibitor candidates were suggested based on molecular modeling studies and a few of them were tested using in vitro analysis. These compounds were determined as good inhibitors against human CA IX target with Ki in the range of 0.85-1.58 μM. In order to predict the pharmaceutical properties of the selected compounds, ADME (absorption, distribution, metabolism and excretion) analysis was also carried out.

Kinetic and in silico studies of hydroxy-based inhibitors of carbonic anhydrase isoforms I and II

A series of hydroxy and phenolic compounds have been assayed for the inhibition of two physiologically relevant carbonic anhydrase (CA, EC 4.2.1.1) isozymes, the cytosolic human isozymes I and II. The investigated molecules showed inhibition constants in the range of 1.07-4003 and 0.09-31.5 μM at the hCA I and hCA II enzymes, respectively. In order to investigate the binding mechanisms of these inhibitors, in silico studies were also applied. Molecular docking scores of the studied compounds are compared using three different scoring algorithms, namely Glide/SP, Glide/XP and Glide/IFD. In addition, different ADME (absorption, distribution, metabolism and excretion) analysis was performed. All the examined compounds were found within the acceptable range of pharmacokinetic profiles.