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.13,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 [3H]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 [3H]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.

Charge-changing point mutations in the E protein of tick-borne encephalitis virus.

Introduction of point mutations is one of the forces enabling arboviruses to rapidly adapt in a changing environment. The influence of these mutations on the properties of the virus is not always obvious. In this study, we attempted to clarify this influence using an in silico approach. Using molecular dynamics (MD) simulations, we investigated how the position of charge-changing point mutations influences the structure and conformational stability of the E protein for a set of variants of a single TBEV strain. The computational findings were supported by experimental evaluation of relevant properties of virions, such as binding to heparan sulfate, thermostability, and susceptibility of the viral hemagglutinating activity to detergents. Our results also point to relationships between E protein dynamics and viral neuroinvasiveness.

Conformational Dynamics and Stability of Bilayers Formed by Mycolic Acids from the Mycobacterium tuberculosis Outer Membrane.

Bilayers of mycolic acids (MAs) form the outer membrane of Mycobacterium tuberculosis that has high strength and extremely low permeability for external molecules (including antibiotics). For the first time, we were able to study them using the all-atom long-term molecular dynamic simulations (from 300 ns up to 1.2 µs) in order to investigate the conformational changes and most favorable structures of the mycobacterial membranes. The structure and properties of the membranes are crucially dependent on the initial packing of the α-mycolic acid (AMA) molecules, as well as on the presence of the secondary membrane components, keto- and methoxy mycolic acids (KMAs and MMAs). In the case of AMA-based membranes, the most labile conformation is W while other types of conformations (sU as well as sZ, eU, and eZ) are much more stable. In the multicomponent membranes, the presence of the KMA and MMA components (in the W conformation) additionally stabilizes both the W and eU conformations of AMA. The membrane in which AMA prevails in the eU conformation is much thicker and, at the same time, much denser. Such a packing of the MA molecules promotes the formation of a significantly stronger outer mycobacterial membrane that should be much more resistant to the threatening external factors.

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.

Conjugates of Tacrine and Salicylic Acid Derivatives as New Promising Multitarget Agents for Alzheimer's Disease.

A series of previously synthesized conjugates of tacrine and salicylamide was extended by varying the structure of the salicylamide fragment and using salicylic aldehyde to synthesize salicylimine derivatives. The hybrids exhibited broad-spectrum biological activity. All new conjugates were potent inhibitors of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) with selectivity toward BChE. The structure of the salicylamide moiety exerted little effect on anticholinesterase activity, but AChE inhibition increased with spacer elongation. The most active conjugates were salicylimine derivatives: IC50 values of the lead compound 10c were 0.0826 µM (AChE) and 0.0156 µM (BChE), with weak inhibition of the off-target carboxylesterase. The hybrids were mixed-type reversible inhibitors of both cholinesterases and displayed dual binding to the catalytic and peripheral anionic sites of AChE in molecular docking, which, along with experimental results on propidium iodide displacement, suggested their potential to block AChE-induced ß-amyloid aggregation. All conjugates inhibited Aß42 self-aggregation in the thioflavin test, and inhibition increased with spacer elongation. Salicylimine 10c and salicylamide 5c with (CH2)8 spacers were the lead compounds for inhibiting Aß42 self-aggregation, which was corroborated by molecular docking to Aß42. ABTS•+-scavenging activity was highest for salicylamides 5a-c, intermediate for salicylimines 10a-c, low for F-containing salicylamides 7, and non-existent for methoxybenzoylamides 6 and difluoromethoxybenzoylamides 8. In the FRAP antioxidant (AO) assay, the test compounds displayed little or no activity. Quantum chemical analysis and molecular dynamics (MD) simulations with QM/MM potentials explained the AO structure-activity relationships. All conjugates were effective chelators of Cu2+, Fe2+, and Zn2+, with molar compound/metal (Cu2+) ratios of 2:1 (5b) and ~1:1 (10b). Conjugates exerted comparable or lower cytotoxicity than tacrine on mouse hepatocytes and had favorable predicted intestinal absorption and blood-brain barrier permeability. The overall results indicate that the synthesized conjugates are promising new multifunctional agents for the potential treatment of AD.

Machine Learning Prediction of Mycobacterial Cell Wall Permeability of Drugs and Drug-like Compounds.

The cell wall of Mycobacterium tuberculosis and related organisms has a very complex and unusual organization that makes it much less permeable to nutrients and antibiotics, leading to the low activity of many potential antimycobacterial drugs against whole-cell mycobacteria compared to their isolated molecular biotargets. The ability to predict and optimize the cell wall permeability could greatly enhance the development of novel antitubercular agents. Using an extensive structure-permeability dataset for organic compounds derived from published experimental big data (5371 compounds including 2671 penetrating and 2700 non-penetrating compounds), we have created a predictive classification model based on fragmental descriptors and an artificial neural network of a novel architecture that provides better accuracy (cross-validated balanced accuracy 0.768, sensitivity 0.768, specificity 0.769, area under ROC curve 0.911) and applicability domain compared with the previously published results.

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 EC50 values below 1.9 and 1.3 nM, respectively, and Chikungunya virus (CHIKV) in cytopathic effect inhibition test with EC50 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 - EC50 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 1O2 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.

Assessing How Residual Errors of Scoring Functions Correlate to Ligand Structural Features.

Scoring functions (SFs) are ubiquitous tools for early stage drug discovery. However, their accuracy currently remains quite moderate. Despite a number of successful target-specific SFs appearing recently, up until now, no ideas on how to systematically improve the general scope of SFs have been formulated. In this work, we hypothesized that the specific features of ligands, corresponding to interactions well appreciated by medicinal chemists (e.g., hydrogen bonds, hydrophobic and aromatic interactions), might be responsible, in part, for the remaining SF errors. The latter provides direction to efforts aimed at the rational and systematic improvement of SF accuracy. In this proof-of-concept work, we took a CASF-2016 coreset of 285 ligands as a basis for comparison and calculated the values of scores for a representative panel of SFs (including AutoDock 4.2, AutoDock Vina, X-Score, NNScore2.0, ΔVina RF20, and DSX). The residual error of linear correlation of each SF value, with the experimental values of affinity and activity, was then analyzed in terms of its correlation with the presence of the fragments responsible for certain medicinal chemistry defined interactions. We showed that, despite the fact that SFs generally perform reasonably, there is room for improvement in terms of better parameterization of interactions involving certain fragments in ligands. Thus, this approach opens a potential way for the systematic improvement of SFs without their significant complication. However, the straightforward application of the proposed approach is limited by the scarcity of reliable available data for ligand-receptor complexes, which is a common problem in the field.

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.

Synthesis and Cytotoxic Activity of 1,2,4-Triazolo-Linked Bis-Indolyl Conjugates as Dual Inhibitors of Tankyrase and PI3K.

A series of new 1,2,4-triazolo-linked bis-indolyl conjugates (15a-r) were prepared by multistep synthesis and evaluated for their cytotoxic activity against various human cancer cell lines. It was observed that they were more susceptible to colon and breast cancer cells. Conjugates 15o (IC50 = 2.04 µM) and 15r (IC50 = 0.85 µM) illustrated promising cytotoxicity compared to 5-fluorouracil (5-FU, IC50 = 5.31 µM) against the HT-29 cell line. Interestingly, 15o and 15r induced cell cycle arrest at the G0/G1 phase and disrupted the mitochondrial membrane potential. Moreover, these conjugates led to apoptosis in HT-29 at 2 µM and 1 µM, respectively, and also enhanced the total ROS production as well as the mitochondrial-generated ROS. Immunofluorescence and Western blot assays revealed that these conjugates reduced the expression levels of the PI3K-P85, ß-catenin, TAB-182, ß-actin, AXIN-2, and NF-κB markers that are involved in the ß-catenin pathway of colorectal cancer. The results of the in silico docking studies of 15r and 15o further support their dual inhibitory behaviour against PI3K and tankyrase. Interestingly, the conjugates have adequate ADME-toxicity parameters based on the calculated results of the molecular dynamic simulations, as we found that these inhibitors (15r) influenced the conformational flexibility of the 4OA7 and 3L54 proteins.

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.