Do electrostatic interactions make a difference in physics-based AutoDock4 scoring function?

Physics-based scoring function AutoDock4 is one of the most successfully applied tools in the area of structure-based drug design. However, current scoring functions are still far from being perfect. In a recent work highlighting the strengths and deficiencies of current scoring functions, we discovered that the residual error of ΔGbind predictions made by AutoDock4 is highly correlated to the presence of formally charged fragments in a ligand. In this work, we study how the use of the high-quality atomic charges, applied for contemporary force fields calculation, affects the quality of the experimental ΔGbind prediction by means of AutoDock4. We initially expected that the previously found discrepancy could be attributed to the Gasteiger charges used within AutoDock4. We show that AutoDock4 is, surprisingly, not sensitive to the charges used, and the use of QC-derived atomic charges does not lead to any statistical improvements. We also briefly discuss the role of the explicit empirical hydrogen bond term along with the electrostatic term.

Temperature-Induced Restructuring of Mycolic Acid Bilayers Modeling the Mycobacterium tuberculosis Outer Membrane: A Molecular Dynamics Study

The emergence of new drug-resistant strains of the tuberculosis pathogen Mycobacterium tuberculosis (Mtb) is a new challenge for modern medicine. Its resistance capacity is closely related to the properties of the outer membrane of the Mtb cell wall, which is a bilayer membrane formed by mycolic acids (MAs) and their derivatives. To date, the molecular mechanisms of the response of the Mtb outer membrane to external factors and, in particular, elevated temperatures have not been sufficiently studied. In this work, we consider the temperature-induced changes in the structure, ordering, and molecular mobility of bilayer MA membranes of various chemical and conformational compositions. Using all-atom long-term molecular dynamics simulations of various MA membranes, we report the kinetic parameters of temperature-dependent changes in the MA self-diffusion coefficients and conformational compositions, including the apparent activation energies of these processes, as well as the characteristic times of ordering changes and the features of phase transitions occurring over a wide range of elevated temperatures. Understanding these effects could be useful for the prevention of drug resistance and the development of membrane-targeting pharmaceuticals, as well as in the design of membrane-based materials.

Free Energy Barriers for Passive Drug Transport through the Mycobacterium tuberculosis Outer Membrane: A Molecular Dynamics Study

The emergence of multi-drug-resistant tuberculosis strains poses a significant challenge to modern medicine. The development of new antituberculosis drugs is hindered by the low permeability of many active compounds through the extremely strong bacterial cell wall of mycobacteria. In order to estimate the ability of potential antimycobacterial agents to diffuse through the outer mycolate membrane, the free energy profiles, the corresponding activation barriers, and possible permeability modes of passive transport for a series of known antibiotics, modern antituberculosis drugs, and prospective active drug-like molecules were determined using molecular dynamics simulations with the all-atom force field and potential of mean-force calculations. The membranes of different chemical and conformational compositions, density, thickness, and ionization states were examined. The typical activation barriers for the low-mass molecules penetrating through the most realistic membrane model were 6-13 kcal/mol for isoniazid, pyrazinamide, and etambutol, and 19 and 25 kcal/mol for bedaquilin and rifampicin. The barriers for the ionized molecules are usually in the range of 37-63 kcal/mol. The linear regression models were derived from the obtained data, allowing one to estimate the permeability barriers from simple physicochemical parameters of the diffusing molecules, notably lipophilicity and molecular polarizability.

Combining Experimental and Computational Methods to Produce Conjugates of Anticholinesterase and Antioxidant Pharmacophores with Linker Chemistries Affecting Biological Activities Related to Treatment of Alzheimer's Disease

Effective therapeutics for Alzheimer's disease (AD) are in great demand worldwide. In our previous work, we responded to this need by synthesizing novel drug candidates consisting of 4-amino-2,3-polymethylenequinolines conjugated with butylated hydroxytoluene via fixed-length alkylimine or alkylamine linkers (spacers) and studying their bioactivities pertaining to AD treatment. Here, we report significant extensions of these studies, including the use of variable-length spacers and more detailed biological characterizations. Conjugates were potent inhibitors of acetylcholinesterase (AChE, the most active was 17d IC50 15.1 ± 0.2 nM) and butyrylcholinesterase (BChE, the most active was 18d: IC50 5.96 ± 0.58 nM), with weak inhibition of off-target carboxylesterase. Conjugates with alkylamine spacers were more effective cholinesterase inhibitors than alkylimine analogs. Optimal inhibition for AChE was exhibited by cyclohexaquinoline and for BChE by cycloheptaquinoline. Increasing spacer length elevated the potency against both cholinesterases. Structure-activity relationships agreed with docking results. Mixed-type reversible AChE inhibition, dual docking to catalytic and peripheral anionic sites, and propidium iodide displacement suggested the potential of hybrids to block AChE-induced β-amyloid (Aβ) aggregation. Hybrids also exhibited the inhibition of Aβ self-aggregation in the thioflavin test; those with a hexaquinoline ring and C8 spacer were the most active. Conjugates demonstrated high antioxidant activity in ABTS and FRAP assays as well as the inhibition of luminol chemiluminescence and lipid peroxidation in mouse brain homogenates. Quantum-chemical calculations explained antioxidant results. Computed ADMET profiles indicated favorable blood-brain barrier permeability, suggesting the CNS activity potential. Thus, the conjugates could be considered promising multifunctional agents for the potential treatment of AD.

5-Nitroisoxazoles in SNAr Reactions: A Novel Chemo- and Regioselective Approach to Isoxazole-Based Bivalent Ligands of AMPA Receptors

An efficient regioselective approach to novel functionalized bis(isoxazoles) with a variety of aromatic and aliphatic linkers was elaborated, based on the heterocyclization reaction of electrophilic alkenes under the treatment with tetranitromethane-triethylamine complex affording 3-EWG-5-nitroisoxazoles. The subsequent SNAr reactions of 5-nitroisoxazoles with various O,O-, N,N- and S,S-bis(nucleophiles) provide a wide range of bis(isoxazole) derivatives in good isolated yields. Employing an elaborated method, a series of novel bis(3-EWG-isoxazoles) as the promising allosteric modulators of AMPA receptors were designed and synthesized. The effect of the compounds on the kainate-induced currents was studied in the patch clamp experiments, revealing modulator properties for several of them. The best positive modulator potency was found for dimethyl 5,5'-(ethane-1,2-diylbis(sulfanediyl))bis(isoxazole-3-carboxylate), which potentiated the kainate-induced currents in a wide concentration range (10-12-10-6 M) with maximum potentiation of 77% at 10-10 M. The results were rationalized using molecular docking and molecular dynamics simulations of modulator complexes with the dimeric ligand-binding domain of the GluA2 AMPA receptor. The predicted physicochemical, ADMET, and PAINS properties confirmed that the AMPA receptor modulators based on the bis(isoxazole) scaffold may serve as potential lead compounds for the development of neuroprotective drugs.

Derivatives of 9-phosphorylated acridine as butyrylcholinesterase inhibitors with antioxidant activity and the ability to inhibit β-amyloid self-aggregation: potential therapeutic agents for Alzheimer's disease

We investigated the inhibitory activities of novel 9-phosphoryl-9,10-dihydroacridines and 9-phosphorylacridines against acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and carboxylesterase (CES). We also studied the abilities of the new compounds to interfere with the self-aggregation of β-amyloid (Aβ42) in the thioflavin test as well as their antioxidant activities in the ABTS and FRAP assays. We used molecular docking, molecular dynamics simulations, and quantum-chemical calculations to explain experimental results. All new compounds weakly inhibited AChE and off-target CES. Dihydroacridines with aryl substituents in the phosphoryl moiety inhibited BChE; the most active were the dibenzyloxy derivative 1d and its diphenethyl bioisostere 1e (IC50 = 2.90 ± 0.23 µM and 3.22 ± 0.25 µM, respectively). Only one acridine, 2d, an analog of dihydroacridine, 1d, was an effective BChE inhibitor (IC50 = 6.90 ± 0.55 μM), consistent with docking results. Dihydroacridines inhibited Aβ42 self-aggregation; 1d and 1e were the most active (58.9% ± 4.7% and 46.9% ± 4.2%, respectively). All dihydroacridines 1 demonstrated high ABTS•+-scavenging and iron-reducing activities comparable to Trolox, but acridines 2 were almost inactive. Observed features were well explained by quantum-chemical calculations. ADMET parameters calculated for all compounds predicted favorable intestinal absorption, good blood-brain barrier permeability, and low cardiac toxicity. Overall, the best results were obtained for two dihydroacridine derivatives 1d and 1e with dibenzyloxy and diphenethyl substituents in the phosphoryl moiety. These compounds displayed high inhibition of BChE activity and Aβ42 self-aggregation, high antioxidant activity, and favorable predicted ADMET profiles. Therefore, we consider 1d and 1e as lead compounds for further in-depth studies as potential anti-AD preparations.

New Allosteric Modulators of AMPA Receptors: Synthesis and Study of Their Functional Activity by Radioligand-Receptor Binding Analysis

The synthetic approaches to three new AMPA receptor modulators-derivatives of 1,11-dimethyl-3,6,9-triazatricyclo[7.3.1.1^3,11]tetradecane-4,8,12-trione-had been developed and all steps of synthesis were optimized. The structures of the compounds contain tricyclic cage and indane fragments necessary for binding with the target receptor. Their physiological activity was studied by radioligand-receptor binding analysis using [³H]PAM-43 as a reference ligand, which is a highly potent positive allosteric modulator of AMPA receptors. The results of radioligand-binding studies indicated the high potency of two synthesized compounds to bind with the same targets as positive allosteric modulator PAM-43 (at least on AMPA receptors). We suggest that the Glu-dependent specific binding site of [³H]PAM-43 or the receptor containing this site may be one of the targets of the new compounds. We also suggest that enhanced radioligand binding may indicate the existence of synergistic effects of compounds 11b and 11c with respect to PAM-43 binding to the targets. At the same time, these compounds may not compete directly with PAM-43 for its specific binding sites but bind to other specific sites of this biotarget, changing its conformation and thereby causing a synergistic effect of cooperative interaction. It can be expected that the newly synthesized compounds will also have pronounced effects on the glutamatergic system of the mammalian brain.

On hidden anisotropy of formally charged fragments

The proper and precise reproduction of the molecular electrostatic potential (MEP) is crucial to describe correctly electrostatic interactions in molecular modeling. Most of the classical molecular mechanics force fields for biomolecules and drug-like molecules use the atom-centered point charges to describe MEP. However, it has been systematically pointed out in literature that such an approximation is not always enough, and some groups, like amino group or heavy halogens, require the use of anisotropic model for better description of their MEP. At the same time, the formally charged groups have not been as extensively and systematically studied as their neutral counterparts. In this report, we demonstrate that the anisotropic models for formally charged groups do bring improvements in the reference MEP reproduction, that are comparable in magnitude to those for neutral groups.

5-(Perylen-3-ylethynyl)uracil as an antiviral scaffold: Potent suppression of enveloped virus reproduction by 3-methyl derivatives in vitro

Amphipathic nucleoside and non-nucleoside derivatives of pentacyclic aromatic hydrocarbon perylene are known as potent non-cytotoxic broad-spectrum antivirals. Here we report 3-methyl-5-(perylen-3-ylethynyl)-uracil-1-acetic acid and its amides, a new series of compounds based on a 5-(perylen-3-ylethynyl)-uracil scaffold. The compounds demonstrate pronounced in vitro activity against arthropod-borne viruses, namely tick-borne encephalitis virus (TBEV) and yellow fever virus (YFV), in plaque reduction assays with EC₅₀ values below 1.9 and 1.3 nM, respectively, and Chikungunya virus (CHIKV) in cytopathic effect inhibition test with EC₅₀ values below 3.2 μM. The compounds are active against respiratory viruses as well: severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) in cytopathic effect inhibition test and influenza A virus (IAV) in virus titer reduction experiments are inhibited - EC₅₀ values below 51 nM and 2.2 μM, respectively. The activity stems from the presence of a hydrophobic perylene core, and all of the synthesized compounds exhibit comparable ¹O₂ generation rates. Nonetheless, activity can vary by orders of magnitude depending on the hydrophilic part of the molecule, suggesting a complex mode of action. A time-of-addition experiment and fluorescent imaging indicate that the compounds inhibit viral fusion in a dose-dependent manner. The localization of the compound in the lipid bilayers and visible damage to the viral envelope suggest the membrane as the primary target. Dramatic reduction of antiviral activity with limited irradiation or under treatment with antioxidants further cements the idea of photoinduced ROS-mediated viral envelope damage being the mode of antiviral action.

Diversity of AMPA Receptor Ligands: Chemotypes, Binding Modes, Mechanisms of Action, and Therapeutic Effects

L-Glutamic acid is the main excitatory neurotransmitter in the central nervous system (CNS). Its associated receptors localized on neuronal and non-neuronal cells mediate rapid excitatory synaptic transmission in the CNS and regulate a wide range of processes in the brain, spinal cord, retina, and peripheral nervous system. In particular, the glutamate receptors selective to α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) also play an important role in numerous neurological disorders and attract close attention as targets for the creation of new classes of drugs for the treatment or substantial correction of a number of serious neurodegenerative and neuropsychiatric diseases. For this reason, the search for various types of AMPA receptor ligands and studies of their properties are attracting considerable attention both in academic institutions and in pharmaceutical companies around the world. This review focuses mainly on the advances in this area published since 2017. Particular attention is paid to the structural diversity of new chemotypes of agonists, competitive AMPA receptor antagonists, positive and negative allosteric modulators, transmembrane AMPA regulatory protein (TARP) dependent allosteric modulators, ion channel blockers as well as their binding sites. This review also presents the studies of the mechanisms of action of AMPA receptor ligands that mediate their therapeutic effects.

Novel Nanomolar Allosteric Modulators of AMPA Receptor of Bis(pyrimidine) Series: Synthesis, Biotesting and SAR Analysis

Positive allosteric modulators (PAMs) of AMPA receptors represent attractive candidates for the development of drugs for the treatment of cognitive and neurodegenerative disorders. Dimeric molecules have been reported to have an especially potent modulating effect, due to the U-shaped form of the AMPA receptor's allosteric binding site. In the present work, novel bis(pyrimidines) were studied as AMPA receptor modulators. A convenient and flexible preparative approach to bis(pyrimidines) containing a hydroquinone linker was elaborated, and a series of derivatives with varied substituents was obtained. The compounds were examined in the patch clamp experiments for their influence on the kainate-induced currents, and 10 of them were found to have potentiating properties. The best potency was found for 2-methyl-4-(4-((2-methyl-5,6,7,8-tetrahydroquinazolin-4-yl)oxy)phenoxy)-6,7,8,9-tetrahydro-5H-cyclohepta[d]pyrimidine, which potentiated the kainate-induced currents by up to 77% in all tested concentrations (10^-12-10^-6 M). The results were rationalized via the modeling of modulator complexes with the dimeric ligand binding domain of the GluA2 AMPA receptor, using molecular docking and molecular dynamics simulation. The prediction of ADMET, physicochemical, and PAINS properties of the studied bis(pyrimidines) confirmed that PAMs of this type may act as the potential lead compounds for the development of neuroprotective drugs.

New Multifunctional Agents for Potential Alzheimer's Disease Treatment Based on Tacrine Conjugates with 2-Arylhydrazinylidene-1,3-Diketones

Alzheimer's disease (AD) is considered a modern epidemic because of its increasing prevalence worldwide and serious medico-social consequences, including the economic burden of treatment and patient care. The development of new effective therapeutic agents for AD is one of the most urgent and challenging tasks. To address this need, we used an aminoalkylene linker to combine the well-known anticholinesterase drug tacrine with antioxidant 2-tolylhydrazinylidene-1,3-diketones to create 3 groups of hybrid compounds as new multifunctional agents with the potential for AD treatment. Lead compounds of the new conjugates effectively inhibited acetylcholinesterase (AChE, IC50 0.24-0.34 µM) and butyrylcholinesterase (BChE, IC50 0.036-0.0745 µM), with weak inhibition of off-target carboxylesterase. Anti-AChE activity increased with elongation of the alkylene spacer, in agreement with molecular docking, which showed compounds binding to both the catalytic active site and peripheral anionic site (PAS) of AChE, consistent with mixed type reversible inhibition. PAS binding along with effective propidium displacement suggest the potential of the hybrids to block AChE-induced β-amyloid aggregation, a disease-modifying effect. All of the conjugates demonstrated metal chelating ability for Cu2+, Fe2+, and Zn2+, as well as high antiradical activity in the ABTS test. Non-fluorinated hybrid compounds 6 and 7 also showed Fe3+ reducing activity in the FRAP test. Predicted ADMET and physicochemical properties of conjugates indicated good CNS bioavailability and safety parameters acceptable for potential lead compounds at the early stages of anti-AD drug development.

Is an Inductive Effect Explicit Account Required for Atomic Charges Aimed at Use within the Force Fields?

Polarization and inductive effects are the concepts that have been widely used in qualitative and even quantitative descriptions of experimentally observed properties in chemistry. The polarization effect has proven to be important in cases of biomolecular modeling though still the vast majority of molecular simulations use the classical non-polarizable force fields. In the last few decades, a lot of effort has been put into promoting the polarization effect and incorporating it into modern force fields and charge calculation methods. In contrast, the inductive effect has not attracted such attention and is effectively absent in both classic and modern force fields. Thus, a question is whether this difference corresponds to the difference in the physical significance of the effects and their explicit account, or is an artifact that should be corrected in the next generation of force fields. The significance of the electronic effects is studied in this paper through the prism of performance of specific models for atomic charge calculation that take into explicit account a nested set of effects: the formal charge, the nearest neighbors, the inductive effect, and finally the model, which takes into account all effects, which are possible to account for using atomic charges. The specific choice for the methods is the following: formal charges, MMFF94 bond charge increments, Dynamic Electronegativity Relaxation (DENR), and RESP. We propose a special scheme for the separate estimation of each particular effect contribution. By pairwise comparing the residual molecular electrostatic potential (MEP) errors of those charge models (aimed at best reproducing the quantum chemical reference MEP), we sequentially revealed how the account of each effect contributes to the better-quality MEP reproduction. The following relative importance of effects was estimated; thus, the natural hierarchy of the effects was established. First, the account of formal charges is of primordial importance. Second, the nearest neighbors account is the next in significance. Third, the explicit account of inductive effect in empirical charge calculation schemes was shown to significantly─both qualitatively and quantitatively─improve the quality of MEP reproduction. Fourth, the contribution of polarization is indirectly assessed. Surprisingly, it is of the order of magnitude of the inductive effect even for the molecular systems, for which it is anticipated to be more significant. Finally, the relative importance of anisotropic effects in neutral molecules was additionally reviewed.

Biphenyl scaffold for the design of NMDA-receptor negative modulators: molecular modeling, synthesis, and biological activity

NMDA (N-methyl-d-aspartate) receptor antagonists are promising tools for the treatment of a wide variety of central nervous system impairments including major depressive disorder. We present here the activity optimization process of a biphenyl-based NMDA negative allosteric modulator (NAM) guided by free energy calculations, which led to a 100 times activity improvement (IC50 = 50 nM) compared to a hit compound identified in virtual screening. Preliminary calculation results suggest a low affinity for the human ether-a-go-go-related gene ion channel (hERG), a high affinity for which was earlier one of the main obstacles for the development of first-generation NMDA-receptor negative allosteric modulators. The docking study and the molecular dynamics calculations suggest a completely different binding mode (ifenprodil-like) compared to another biaryl-based NMDA NAM EVT-101.

In Silico Structure-Based Approach for Group Efficiency Estimation in Fragment-Based Drug Design Using Evaluation of Fragment Contributions

The notion of a contribution of a specific group in an organic molecule’s property and/or activity is both common in our thinking and is still not strictly correct due to the inherent non-additivity of free energy with respect to molecular fragments composing a molecule. The fragment- based drug discovery (FBDD) approach has proven to be fruitful in addressing the above notions. The main difficulty of the FBDD, however, is in its reliance on the low throughput and expensive experimental means of determining the fragment-sized molecules binding. In this article we propose a way to enhance the throughput and availability of the FBDD methods by judiciously using an in silico means of assessing the contribution to ligand-receptor binding energy of fragments of a molecule under question using a previously developed in silico Reverse Fragment Based Drug Discovery (R-FBDD) approach. It has been shown that the proposed structure-based drug discovery (SBDD) type of approach fills in the vacant niche among the existing in silico approaches, which mainly stem from the ligand-based drug discovery (LBDD) counterparts. In order to illustrate the applicability of the approach, our work retrospectively repeats the findings of the use case of an FBDD hit-to-lead project devoted to the experimentally based determination of additive group efficiency (GE)—an analog of ligand efficiency (LE) for a group in the molecule—using the Free-Wilson (FW) decomposition. It is shown that in using our in silico approach to evaluate fragment contributions of a ligand and to estimate GE one can arrive at similar decisions as those made using the experimentally determined activity-based FW decomposition. It is also shown that the approach is rather robust to the choice of the scoring function, provided the latter demonstrates a decent scoring power. We argue that the proposed approach of in silico assessment of GE has a wider applicability domain and expect that it will be widely applicable to enhance the net throughput of drug discovery based on the FBDD paradigm.

Bis-Amiridines as Acetylcholinesterase and Butyrylcholinesterase Inhibitors: N-Functionalization Determines the Multitarget Anti-Alzheimer’s Activity Profile

Using two ways of functionalizing amiridine—acylation with chloroacetic acid chloride and reaction with thiophosgene—we have synthesized new homobivalent bis-amiridines joined by two different spacers—bis-N-acyl-alkylene (3) and bis-N-thiourea-alkylene (5) —as potential multifunctional agents for the treatment of Alzheimer’s disease (AD). All compounds exhibited high inhibitory activity against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) with selectivity for BChE. These new agents displayed negligible carboxylesterase inhibition, suggesting a probable lack of untoward drug–drug interactions arising from hydrolytic biotransformation. Compounds 3 with bis-N-acyl-alkylene spacers were more potent inhibitors of both cholinesterases compared to compounds 5 and the parent amiridine. The lead compounds 3a–c exhibited an IC₅₀(AChE) = 2.9–1.4 µM, IC₅₀(BChE) = 0.13–0.067 µM, and 14–18% propidium displacement at 20 μM. Kinetic studies of compounds 3a and 5d indicated mixed-type reversible inhibition. Molecular docking revealed favorable poses in both catalytic and peripheral AChE sites. Propidium displacement from the peripheral site by the hybrids suggests their potential to hinder AChE-assisted Aβ₄₂ aggregation. Conjugates 3 had no effect on Aβ₄₂ self-aggregation, whereas compounds 5c–e (m = 4, 5, 6) showed mild (13–17%) inhibition. The greatest difference between conjugates 3 and 5 was their antioxidant activity. Bis-amiridines 3 with N-acylalkylene spacers were nearly inactive in ABTS and FRAP tests, whereas compounds 5 with thiourea in the spacers demonstrated high antioxidant activity, especially in the ABTS test (TEAC = 1.2–2.1), in agreement with their significantly lower HOMO-LUMO gap values. Calculated ADMET parameters for all conjugates predicted favorable blood–brain barrier permeability and intestinal absorption, as well as a low propensity for cardiac toxicity. Thus, it was possible to obtain amiridine derivatives whose potencies against AChE and BChE equaled (5) or exceeded (3) that of the parent compound, amiridine. Overall, based on their expanded and balanced pharmacological profiles, conjugates 5c–e appear promising for future optimization and development as multitarget anti-AD agents.

Novel substituted 5-methyl-4-acylaminoisoxazoles as antimitotic agents: Evaluation of selectivity to LNCaP cancer cells

A series of novel antimitotic agents was designed using the replacement of heterocyclic cores in two tubulin-targeting lead molecules with the acylated 4-aminoisoxazole moiety. Target compounds were synthesized via heterocyclization of β-aryl-substituted vinylketones by tert-butyl nitrite in the presence of water as a key step. 4-Methyl-N-[5-methyl-3-(3,4,5-trimethoxyphenyl)isoxazol-4-yl]benzamide (1aa) was found to stimulate partial depolymerization of microtubules of human lung carcinoma A549 cells at a high concentration of 100 µM and to totally inhibit cell growth (IC₅₀  = 0.99 µM) and cell viability (IC₅₀  = 0.271 µM) in the nanomolar to submicromolar concentration range. These data provide evidence of the multitarget profile of the cytotoxic action of compound 1aa. The SAR study demonstrated that the 3,4,5-trimethoxyphenyl residue is the key structural parameter determining the efficiency both towards tubulin and other molecular targets. The cytotoxicity of 3-methyl-N-[5-methyl-3-(3,4,5-trimethoxyphenyl)isoxazol-4-yl]benzamide (1ab) to the androgen-sensitive human prostate adenocarcinoma cancer cell line LNCaP (IC₅₀  = 0.301 µM) was approximately one order of magnitude higher than that to the conditionally normal cells lines WI-26 VA4 (IC₅₀  = 2.26 µM) and human umbilical vein endothelial cells (IC₅₀  = 5.58 µM) and significantly higher than that to primary fibroblasts (IC₅₀  > 75 µM).

A Facile Approach to Bis(isoxazoles), Promising Ligands of the AMPA Receptor

A convenient synthetic approach to novel functionalized bis(isoxazoles), the promising bivalent ligands of the AMPA receptor, was elaborated. It was based on the heterocyclization reactions of readily available electrophilic alkenes with the tetranitromethane-triethylamine complex. The structural diversity of the synthesized compounds was demonstrated. In the electrophysiological experiments using the patch clamp technique on Purkinje neurons, the compound 1,4-phenylenedi(methylene)bis(5-aminoisoxazole-3-carboxylate) was shown to be highly potent positive modulator of the AMPA receptor, potentiating kainate-induced currents up to 70% at 10⁻¹¹ M.

[Ligands of the AMPA-subtype glutamate receptors: mechanisms of action and novel chemotypes]

Ionotropic glutamate receptors of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) subtype play a key role in synaptic plasticity representing one of the mechanisms for learning and memory formation. They can also serve as targets for the development of novel classes of pharmaceuticals for the treatment or substantive correction of many serious neurodegenerative and psychoneurological disorders. The search and studies of various types of AMPA receptor ligands attract considerable attention from academic organizations and pharmaceutical companies all over the world. This review mainly focuses on recent advances in this field. The architecture and operational mechanism of the receptor as well as its major binding sites and ligand types are considered. Special attention is paid to the studies of mechanisms of action and novel chemotypes of AMPA receptor agonists and competitive antagonists, positive and negative allosteric modulators, auxiliary protein-dependent allosteric modulators, and ion channel blockers.

Spectrum of antiviral activity of 4-aminopyrimidine N-oxides against a broad panel of tick-borne encephalitis virus strains

Tick-borne encephalitis is an important human arbovirus neuroinfection spread across the Northern Eurasia. Inhibitors of tick-borne encephalitis virus (TBEV) strain Absettarov, presumably targeting E protein n-octyl-β-d-glucoside (β-OG) pocket, were reported earlier. In this work, these inhibitors were tested in vitro against seven strains representing three main TBEV subtypes. The most potent compound, 2-[(2-methyl-1-oxido-5,6,7,8-tetrahydroquinazolin-4-yl)amino]-phenol, showed EC₅₀ values lower than 22 µM against all the tested strains. Nevertheless, EC₅₀ values for virus samples of certain strains demonstrated a substantial variation, which appeared to be consistent with the presence of E protein not only in infectious virions, but also in non-infectious and immature virus particles, protein aggregates, and membrane complexes.

Amiridine-piperazine hybrids as cholinesterase inhibitors and potential multitarget agents for Alzheimer's disease treatment

We synthesized eleven new amiridine-piperazine hybrids 5a-j and 7 as potential multifunctional agents for Alzheimer's disease (AD) treatment by reacting N-chloroacetylamiridine with piperazines. The compounds displayed mixed-type reversible inhibition of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Conjugates were moderate inhibitors of equine and human BChE with negligible fluctuation in anti-BChE activity, whereas anti-AChE activity was substantially dependent on N4-substitution of the piperazine ring. Compounds with para-substituted aromatic moieties (5g, 5h, and bis-amiridine 7) had the highest anti-AChE activity in the low micromolar range. Top-ranked compound 5h, N-(2,3,5,6,7,8-hexahydro-1H-cyclopenta[b]quinolin-9-yl)-2-[4-(4-nitro-phenyl)-piperazin-1-yl]-acetamide, had an IC₅₀ for AChE = 1.83 ± 0.03 μM (K_i = 1.50 ± 0.12 and αK_i = 2.58 ± 0.23 μM). The conjugates possessed low activity against carboxylesterase, indicating a likely absence of unwanted drug-drug interactions in clinical use. In agreement with analysis of inhibition kinetics and molecular modeling studies, the lead compounds were found to bind effectively to the peripheral anionic site of AChE and displace propidium, indicating their potential to block AChE-induced β-amyloid aggregation. Similar propidium displacement activity was first shown for amiridine. Two compounds, 5c (R = cyclohexyl) and 5e (R = 2-MeO-Ph), exhibited appreciable antioxidant capability with Trolox equivalent antioxidant capacity values of 0.47 ± 0.03 and 0.39 ± 0.02, respectively. Molecular docking and molecular dynamics simulations provided insights into the structure-activity relationships for AChE and BChE inhibition, including the observation that inhibitory potencies and computed pK_a values of hybrids were generally lower than those of the parent molecules. Predicted ADMET and physicochemical properties of conjugates indicated good CNS bioavailability and safety parameters comparable to those of amiridine and therefore acceptable for potential lead compounds at the early stages of anti-AD drug development.

Novel potent bifunctional carboxylesterase inhibitors based on a polyfluoroalkyl-2-imino-1,3-dione scaffold

An expanded series of alkyl 2-arylhydrazinylidene-3-oxo-3-polyfluoroalkylpropionates (HOPs) 3 was obtained via Cu(OAc)₂-catalyzed azo coupling. All were nanomolar inhibitors of carboxylesterase (CES), while moderate or weak inhibitors of acetylcholinesterase and butyrylcholinesterase. Steady-state kinetics studies showed that HOPs 3 are mixed type inhibitors of the three esterases. Molecular docking studies demonstrated that two functional groups in the structure of HOPs, trifluoromethyl ketone (TFK) and ester groups, bind to the CES active site suggesting subsequent reactions: formation of a tetrahedral adduct, and a slow hydrolysis reaction. The results of molecular modeling allowed us to explain some structure-activity relationships of CES inhibition by HOPs 3: their selectivity toward CES in comparison with cholinesterases and the high selectivity of pentafluoroethyl-substituted HOP 3p to hCES1 compared to hCES2. All compounds were predicted to have good intestinal absorption and blood-brain barrier permeability, low cardiac toxicity, good lipophilicity and aqueous solubility, and reasonable overall drug-likeness. HOPs with a TFK group and electron-donor substituents in the arylhydrazone moiety were potent antioxidants. All compounds possessed low cytotoxicity and low acute toxicity. Overall, a new promising type of bifunctional CES inhibitors has been found that are able to interact with the active site of the enzyme with the participation of two functional groups. The results indicate that HOPs have the potential to be good candidates as human CES inhibitors for biomedicinal applications.

Towards Deep Neural Network Models for the Prediction of the Blood-Brain Barrier Permeability for Diverse Organic Compounds

Permeation through the blood–brain barrier (BBB) is among the most important processes controlling the pharmacokinetic properties of drugs and other bioactive compounds. Using the fragmental (substructural) descriptors representing the occurrence number of various substructures, as well as the artificial neural network approach and the double cross-validation procedure, we have developed a predictive in silico LogBB model based on an extensive and verified dataset (529 compounds), which is applicable to diverse drugs and drug-like compounds. The model has good predictivity parameters (Q2=0.815, RMSEcv=0.318) that are similar to or better than those of the most reliable models available in the literature. Larger datasets, and perhaps more sophisticated network architectures, are required to realize the full potential of deep neural networks. The analysis of fragment contributions reveals patterns of influence consistent with the known concepts of structural characteristics that affect the BBB permeability of organic compounds. The external validation of the model confirms good agreement between the predicted and experimental LogBB values for most of the compounds. The model enables the evaluation and optimization of the BBB permeability of potential neuroactive agents and other drug compounds.

New Multifunctional Agents Based on Conjugates of 4-Amino-2,3-polymethylenequinoline and Butylated Hydroxytoluene for Alzheimer's Disease Treatment

New hybrids of 4-amino-2,3-polymethylenequinoline with different sizes of the aliphatic ring linked to butylated hydroxytoluene (BHT) by enaminoalkyl (7) or aminoalkyl (8) spacers were synthesized as potential multifunctional agents for Alzheimer's disease (AD) treatment. All compounds were potent inhibitors of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) with selectivity toward BChE. Lead compound 8c, 2,6-di-tert-butyl-4-{[2-(7,8,9,10- tetrahydro-6H-cyclohepta[b]quinolin-11-ylamino)-ethylimino]-methyl}-phenol exhibited an IC₅₀(AChE) = 1.90 ± 0.16 µM, IC₅₀(BChE) = 0.084 ± 0.008 µM, and 13.6 ± 1.2% propidium displacement at 20 μM. Compounds possessed low activity against carboxylesterase, indicating likely absence of clinically unwanted drug-drug interactions. Kinetics were consistent with mixed-type reversible inhibition of both cholinesterases. Docking indicated binding to catalytic and peripheral AChE sites; peripheral site binding along with propidium displacement suggest the potential of the hybrids to block AChE-induced β-amyloid aggregation, a disease-modifying effect. Compounds demonstrated high antioxidant activity in ABTS and FRAP assays as well as inhibition of luminol chemiluminescence and lipid peroxidation in mouse brain homogenates. Conjugates 8 with amine-containing spacers were better antioxidants than those with enamine spacers 7. Computational ADMET profiles for all compounds predicted good blood-brain barrier distribution (permeability), good intestinal absorption, and medium cardiac toxicity risk. Overall, based on their favorable pharmacological and ADMET profiles, conjugates 8 appear promising as candidates for AD therapeutics.

Discovery of Novel Tankyrase Inhibitors through Molecular Docking-Based Virtual Screening and Molecular Dynamics Simulation Studies

Tankyrase enzymes (TNKS), a core part of the canonical Wnt pathway, are a promising target in the search for potential anti-cancer agents. Although several hundreds of the TNKS inhibitors are currently known, identification of their novel chemotypes attracts considerable interest. In this study, the molecular docking and machine learning-based virtual screening techniques combined with the physico-chemical and ADMET (absorption, distribution, metabolism, excretion, toxicity) profile prediction and molecular dynamics simulations were applied to a subset of the ZINC database containing about 1.7 M commercially available compounds. Out of seven candidate compounds biologically evaluated in vitro for their inhibition of the TNKS2 enzyme using immunochemical assay, two compounds have shown a decent level of inhibitory activity with the IC₅₀ values of less than 10 nM and 10 μM. Relatively simple scores based on molecular docking or MM-PBSA (molecular mechanics, Poisson-Boltzmann, surface area) methods proved unsuitable for predicting the effect of structural modification or for accurate ranking of the compounds based on their binding energies. On the other hand, the molecular dynamics simulations and Free Energy Perturbation (FEP) calculations allowed us to further decipher the structure-activity relationships and retrospectively analyze the docking-based virtual screening performance. This approach can be applied at the subsequent lead optimization stages.

Characterization of a New Positive Allosteric Modulator of AMPA Receptors - PAM-43: Specific Binding of the Ligand and its Ability to Potentiate AMPAR Currents

BACKGROUND

Currently, the most dynamic areas in the glutamate receptor system neurobiology are the identification and development of positive allosteric modulators (PAMs) of glutamate ionotropic receptors. PAM-based drugs are of great interest as promising candidates for the treatment of neurological diseases, such as epilepsy, Alzheimer's disease, schizophrenia, etc. Understanding the molecular mechanisms underlying the biological action of natural and synthetic PAMs is a key point for modifying the original chemical compounds as well as for new drug design.

OBJECTIVE

We are trying to elaborate a system of molecular functional screening of ionotropic glutamate receptor probable PAMs.

METHODS

The system will be based on the radioligand - receptor method of analysis and will allow rapid quantification of new AMPAR probable PAMs molecular activity. We plan to use a tritiumlabeled analogue of recently elaborated ionotropic GluR probable PAM ([3H]PAM-43) as the main radioligand.

RESULTS

Here, we characterized the specific binding of the ligand and its ability to potentiate ionotropic GluR currents. The existence of at least two different sites of [3H]PAM-43 specific binding has been shown. One of the above sites is glutamate-dependent and is characterized by higher affinity. "Patchclamp" technique showed the ability of PAM-43 to potentiate ionotropic GluR currents in rat cerebellar Purkinje neurons in a concentration-dependent manner.

CONCLUSION

The possibility of using PAM-43 as a model compound to study different allosteric effects of potential regulatory drugs (AMPAR allosteric regulators) was shown. [3H]PAM-43 based screening system will allow rapid selection of new AMPAR probable PAM structures and quantification of their molecular activity.

New Positive Allosteric Modulator of AMPA Receptors: in vitro and in vivo Studies

A new derivative of 3,7-diazabicyclo[3.3.1]nonane, which showed a high activity as a positive allosteric modulator of AMPA receptors of the CNS, was studied in electrophysiological experiments. At doses of 0.01 mg/kg, this compound significantly improved the memory of experimental animals disturbed by maximal electric shock. The results indicate that this compound is a promising candidate for preclinical trials and clinical studies as a drug for treatment of a number of psychoneurological diseases.

Bivalent AMPA receptor positive allosteric modulators of the bis(pyrimidine) series

The first example of a novel class of AMPA receptor positive allosteric modulators of the bis(pyrimidine) series having a hydroquinone linker has been obtained and showed a potency to increase kainate-induced currents at subnanomolar concentrations.

Examination of molecular space and feasible structures of bioactive components of humic substances by FTICR MS data mining in ChEMBL database

Humic substances (HS) are complex natural mixtures comprising a large variety of compounds produced during decomposition of decaying biomass. The molecular composition of HS is extremely diverse as it was demonstrated with the use of high resolution mass spectrometry. The building blocks of HS are mostly represented by plant-derived biomolecules (lignins, lipids, tannins, carbohydrates, etc.). As a result, HS show a wide spectrum of biological activity. Despite that, HS remain a 'biological activity black-box' due to unknown structures of constituents responsible for the interaction with molecular targets. In this study, we investigated the antiviral activity of eight HS fractions isolated from peat and coal, as well as of two synthetic humic-like materials. We determined molecular compositions of the corresponding samples using ultra-high resolution Fourier-transform ion cyclotron resonance mass-spectrometry (FTICR MS). Inhibitory activity of HS was studied with respect to reproduction of tick-borne encephalitis virus (TBEV), which is a representative of Flavivirus genus, and to a panel of enteroviruses (EVs). The samples of natural HS inhibited TBEV reproduction already at a concentration of 1 µg/mL, but they did not inhibit reproduction of EVs. We found that the total relative intensity of FTICR MS formulae within elemental composition range commonly attributed to flavonoid-like structures is correlating with the activity of the samples. In order to surmise on possible active structural components of HS, we mined formulae within FTICR MS assignments in the ChEMBL database. Out of 6502 formulae within FTICR MS assignments, 3852 were found in ChEMBL. There were more than 71 thousand compounds related to these formulae in ChEMBL. To support chemical relevance of these compounds to natural HS we applied the previously developed approach of selective isotopic exchange coupled to FTICR MS to obtain structural information on the individual components of HS. This enabled to propose compounds from ChEMBL, which corroborated the labeling data. The obtained results provide the first insight onto the possible structures, which comprise antiviral components of HS and, respectively, can be used for further disclosure of antiviral activity mechanism of HS.

Time-Domain Analysis of Molecular Dynamics Trajectories Using Deep Neural Networks: Application to Activity Ranking of Tankyrase Inhibitors

Molecular dynamics simulations provide valuable insights into the behavior of molecular systems. Extending the recent trend of using machine learning techniques to predict physicochemical properties from molecular dynamics data, we propose to consider the trajectories as multidimensional time series represented by 2D tensors containing the ligand-protein interaction descriptor values for each time step. Similar in structure to the time series encountered in modern approaches for signal, speech, and natural language processing, these time series can be directly analyzed using long short-term memory (LSTM) recurrent neural networks or convolutional neural networks (CNNs). The predictive regression models for the ligand-protein affinity were built for a subset of the PDBbind v.2017 database and applied to inhibitors of tankyrase, an enzyme of the poly(ADP-ribose)-polymerase (PARP) family that can be used in the treatment of colorectal cancer. As an additional test set, a subset of the Community Structure-Activity Resource (CSAR) data set was used. For comparison, the random forest and simple neural network models based on the crystal pose or the trajectory-averaged descriptors were used, as well as the commonly employed docking and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) scores. Convolutional neural networks based on the 2D tensors of ligand-protein interaction descriptors for short (2 ns) trajectories provide the best accuracy and predictive power, reaching the Spearman rank correlation coefficient of 0.73 and Pearson correlation coefficient of 0.70 for the tankyrase test set. Taking into account the recent increase in computational power of modern GPUs and relatively low computational complexity of the proposed approach, it can be used as an advanced virtual screening filter for compound prioritization.

Iterative Solvers for Empirical Partial Atomic Charges: Breaking the Curse of Cubic Numerical Complexity

Rational drug design involves a vast amount of computations to get thermodynamically reliable results and often relies on atomic charges as a means to model electrostatic interactions within the system. Computational inefficiency often hampers the development of new and wider dissemination of the known methods; thus, any source to speed up the calculations without a sacrifice in quality is warranted. At the heart of many empirical methods of calculating atomic charges is the solution of a system of linear algebraic equations (SLAE). The classical method of solving SLAE-the Gauss method-has in general case a cubic computational complexity. It is shown that the use of iterative methods for solving SLAE, characteristic to typical empirical atomic charge calculation methods, makes it possible to significantly reduce the amount of calculations and to obtain a computational complexity approaching a quadratic one. Despite the fact that this phenomenon is well-known in numerical methods, iterative solvers surprisingly do not seem to have been systematically applied to calculation of atomic charges via empirical schemes. Another finding is the relative values of the matrix elements, determined by the physical grounds of the interactions within the empirical system, generally lead to SLAE's with well-defined matrices, suited to use with iterative solvers to fasten computation compared to using the noniterative solvers. This finding broadens the applicability range of atomic charges obtained with empirical methods for such cases as, e.g., account of polarizability via "on-the-fly" recalculation of charges in changing surroundings within the force fields in molecular dynamics settings.

Compounds based on 5-(perylen-3-ylethynyl)uracil scaffold: High activity against tick-borne encephalitis virus and non-specific activity against enterovirus A

Rigid amphipathic fusion inhibitors (RAFIs) are potent antivirals based on a perylene core linked with a nucleoside moiety. Sugar-free analogues of RAFIs, 5-(perylen-3-ylethynyl)uracil-1-acetic acid 1 and its amides 2, were synthesized using combined protection group strategy. Compounds 1 and 2 appeared to have low toxicity on porcine embryo kidney (PEK) or rhabdomiosarcoma (RD) cells together with remarkable activity against enveloped tick-borne encephalitis virus (TBEV): EC₅₀ values vary from 0.077 μM to subnanomolar range. Surprisingly, 3-pivaloyloxymethyl (Pom) protected precursors 7 and 8 showed even more pronounced activity. All the compounds showed no activity against several non-enveloped enteroviruses, except 4-hydroxybutylamides 2d,g, which inhibited the reproduction of enterovirus A71 with EC₅₀ 50-100 μM, with a non-specific mode of action. The results suggest that the carbohydrate moiety of RAFI nucleosides does not play a crucial role in their antiviral action, and biological activity of the 5-(perylen-3-ylethynyl)uracil scaffold can be effectively modulated by substituents in positions 1 and 3. The high antiviral activity of these new compounds, coupled with low toxicity advocate their potential role in antiviral therapy.

Enhanced taxonomy annotation of antiviral activity data from ChEMBL

The discovery of antiviral drugs is a rapidly developing area of medicinal chemistry research. The emergence of resistant variants and outbreaks of poorly studied viral diseases make this area constantly developing. The amount of antiviral activity data available in ChEMBL consistently grows, but virus taxonomy annotation of these data is not sufficient for thorough studies of antiviral chemical space. We developed a procedure for semi-automatic extraction of antiviral activity data from ChEMBL and mapped them to the virus taxonomy developed by the International Committee for Taxonomy of Viruses (ICTV). The procedure is based on the lists of virus-related values of ChEMBL annotation fields and a dictionary of virus names and acronyms mapped to ICTV taxa. Application of this data extraction procedure allows retrieving from ChEMBL 1.6 times more assays linked to 2.5 times more compounds and data points than ChEMBL web interface allows. Mapping of these data to ICTV taxa allows analyzing all the compounds tested against each viral species. Activity values and structures of the compounds were standardized, and the antiviral activity profile was created for each standard structure. Data set compiled using this algorithm was called ViralChEMBL. As case studies, we compared descriptor and scaffold distributions for the full ChEMBL and its `viral' and `non-viral' subsets, identified the most studied compounds and created a self-organizing map for ViralChEMBL. Our approach to data annotation appeared to be a very efficient tool for the study of antiviral chemical space.

Ramified derivatives of 5-(perylen-3-ylethynyl)uracil-1-acetic acid and their antiviral properties

The propargylamide of N3-Pom-protected 5-(perylen-3-ylethynyl)uracil acetic acid, a universal precursor, was used in a CuAAC click reaction for the synthesis of several derivatives, including three ramified molecules with high activities against tick-borne encephalitis virus (TBEV). Pentaerythritol-based polyazides were used for the assembly of molecules containing 2⋯4 antiviral 5-(perylen-3-ylethynyl)uracil scaffolds, the first examples of polyvalent perylene antivirals. Cluster compounds showed enhanced absorbance, however, their fluorescence was reduced due to self-quenching. Due to the solubility issues, Pom group removal succeeded only for compounds with one peryleneethynyluracil unit. Four compounds, including one ramified cluster 9f, showed remarkable 1⋯3 nM EC₅₀ values against TBEV in cell culture.

Ramified clusters of antiviral perylenylethynyl scaffold were prepared using CuAAC reaction of 5-(perylen-3-ylethynyl)-3-pivaloyloxymethyl-1-(propargylamidomethyl)uracil with azides. Compounds inhibited TBEV reproduction at nanomolar concentrations.

Synthesis and cytotoxicity of novel simplified eleutherobin analogues as potential antitumour agents

Synthesis of title compounds based on adamantane and 8-oxabicyclo[3.2.1]octane scaffolds and evaluation of their biological activity.

Molecular design, synthesis and biological evaluation of cage compound-based inhibitors of hepatitis C virus p7 ion channels

The hepatitis C caused by the hepatitis C virus (HCV) is an acute and/or chronic liver disease ranging in severity from a mild brief ailment to a serious lifelong illness that affects up to 3% of the world population and imposes significant and increasing social, economic, and humanistic burden. Over the past decade, its treatment was revolutionized by the development and introduction into clinical practice of the direct acting antiviral (DAA) agents targeting the non-structural viral proteins NS3/4A, NS5A, and NS5B. However, the current treatment options still have important limitations, thus, the development of new classes of DAAs acting on different viral targets and having better pharmacological profile is highly desirable. The hepatitis C virus p7 viroporin is a relatively small hydrophobic oligomeric viral ion channel that plays a critical role during virus assembly and maturation, making it an attractive and validated target for the development of the cage compound-based inhibitors. Using the homology modeling, molecular dynamics, and molecular docking techniques, we have built a representative set of models of the hepatitis C virus p7 ion channels (Gt1a, Gt1b, Gt1b_L20F, Gt2a, and Gt2b), analyzed the inhibitor binding sites, and identified a number of potential broad-spectrum inhibitor structures targeting them. For one promising compound, the binding to these targets was additionally confirmed and the binding modes and probable mechanisms of action were clarified by the molecular dynamics simulations. A number of compounds were synthesized, and the tests of their antiviral activity (using the BVDV model) and cytotoxicity demonstrate their potential therapeutic usefulness and encourage further more detailed studies. The proposed approach is also suitable for the design of broad-spectrum ligands interacting with other multiple labile targets including various viroporins.