Targeting the Human Influenza a Virus: The Methods, Limitations, and Pitfalls of Virtual Screening for Drug-like Candidates Including Scaffold Hopping and Compound Profiling.

In this study, we describe the input data and processing steps to find antiviral lead compounds by a virtual screen. Two-dimensional and three-dimensional filters were designed based on the X-ray crystallographic structures of viral neuraminidase co-crystallized with substrate sialic acid, substrate-like DANA, and four inhibitors (oseltamivir, zanamivir, laninamivir, and peramivir). As a result, ligand-receptor interactions were modeled, and those necessary for binding were utilized as screen filters. Prospective virtual screening (VS) was carried out in a virtual chemical library of over half a million small organic substances. Orderly filtered moieties were investigated based on 2D- and 3D-predicted binding fingerprints disregarding the "rule-of-five" for drug likeness, and followed by docking and ADMET profiling. Two-dimensional and three-dimensional screening were supervised after enriching the dataset with known reference drugs and decoys. All 2D, 3D, and 4D procedures were calibrated before execution, and were then validated. Presently, two top-ranked substances underwent successful patent filing. In addition, the study demonstrates how to work around reported VS pitfalls in detail.

Inhibition of the Human Neuronal Sodium Channel Nav1.9 by Arachidonyl-2-Chloroethylamide, An Analogue of Anandamide in a hNav1.9/rNav1.4 Chimera, An Experimental and in Silico Study.

Cannabinoids regulate analgesia, which has aroused much interest in identifying new pharmacological therapies in the management of refractory pain. Voltage-gated Na+ channels (Navs) play an important role in inflammatory and neuropathic pain. In particular, Nav1.9 is involved in nociception and the understanding of its pharmacology has lagged behind because it is difficult to express in heterologous systems. Here, we utilized the chimeric channel hNav1.9_C4, that comprises the extracellular and transmembrane domains of hNav1.9, co-expressed with the ß1 subunit on CHO-K1 cells to characterize the electrophysiological effects of ACEA, a synthetic surrogate of the endogenous cannabinoid anandamide. ACEA induced a tonic block, decelerated the fast inactivation, markedly shifted steady-state inactivation in the hyperpolarized direction, decreasing the window current and showed use-dependent block, with a high affinity for the inactivated state (ki = 0.84 µM). Thus, we argue that ACEA possess a local anaesthetic-like profile. To provide a mechanistic understanding of its mode of action at the molecular level, we combined induced fit docking with Monte Carlo simulations and electrostatic complementarity. In agreement with the experimental evidence, our computer simulations revealed that ACEA binds Tyr1599 of the local anaesthetics binding site of the hNav1.9, contacting residues that bind cannabinol (CBD) in the NavMs channel. ACEA adopted a conformation remarkably similar to the crystallographic conformation of anandamide on a non-homologous protein, obstructing the Na+ permeation pathway below the selectivity filter to occupy a highly conserved binding pocket at the intracellular side. These results describe a mechanism of action, possibly involved in cannabinoid analgesia.

Allosteric Binding of MDMA to the Human Serotonin Transporter (hSERT) via Ensemble Binding Space Analysis with ΔG Calculations, Induced Fit Docking and Monte Carlo Simulations.

Despite the recent promising results of MDMA (3,4-methylenedioxy-methamphetamine) as a psychotherapeutic agent and its history of misuse, little is known about its molecular mode of action. MDMA enhances monoaminergic neurotransmission in the brain and its valuable psychoactive effects are associated to a dual action on the 5-HT transporter (SERT). This drug inhibits the reuptake of 5-HT (serotonin) and reverses its flow, acting as a substrate for the SERT, which possesses a central binding site (S1) for antidepressants as well as an allosteric (S2) one. Previously, we characterized the spatial binding requirements for MDMA at S1. Here, we propose a structure-based mechanistic model of MDMA occupation and translocation across both binding sites, applying ensemble binding space analyses, electrostatic complementarity, and Monte Carlo energy perturbation theory. Computed results were correlated with experimental data (r = 0.93 and 0.86 for S1 and S2, respectively). Simulations on all hSERT available structures with Gibbs free energy estimations (ΔG) revealed a favourable and pervasive dual binding mode for MDMA at S2, i.e., adopting either a 5-HT or an escitalopram-like orientation. Intermediate ligand conformations were identified within the allosteric site and between the two sites, outlining an internalization pathway for MDMA. Among the strongest and more frequent interactions were salt bridges with Glu494 and Asp328, a H-bond with Thr497, a π-π with Phe556, and a cation-π with Arg104. Similitudes and differences with the allosteric binding of 5-HT and antidepressants suggest that MDMA may have a distinctive chemotype. Thus, our models may provide a framework for future virtual screening studies and pharmaceutical design and to develop hSERT allosteric compounds with a unique psychoactive MDMA-like profile.

Computational Molecular Characterization of the Interaction of Acetylcholine and the NMDA Receptor to Explain the Direct Glycine-Competitive Potentiation of NMDA-Mediated Neuronal Currents.

The activation of N-methyl-d-aspartate receptor (NMDAR) is triggered by the closure of bilobed (D1 and D2) clamshell-like clefts upon binding glycine (Gly) and glutamate. There is evidence that cholinergic compounds modulate NMDAR-mediated currents via direct receptor-ligand interactions; however, molecular bases are unknown. Here, we first propose a mechanistic structure-based explanation for the observed ACh-induced submaximal potentiation of NMDA-elicited currents in striatal neurons by predicting competitive inhibition with Gly. Then, the model was validated, in principle, by confirming that the coapplication of Gly and ACh significantly reduces these neuronal currents. Finally, we delineate the interplay of ACh with the NMDAR by a combination of computational strategies. Crystallographic ACh-bound complexes were studied, revealing a similar ACh binding environment on the GluN1 subunit of the NMDAR. We illustrate how ACh can occupy X-ray monomeric open, dimeric "semiopen" cleft conformations obtained by molecular dynamics and a full-active cryo-EM NMDAR structure, explaining the suboptimal NMDAR electrophysiological activity under the "Venus Flytrap model". At an evolutionary biology level, the binding mode of ACh coincides with that of the homologous ornithine-bound periplasmic LAO binding protein complex. Our computed results indicate an analogous mechanism of action, inasmuch as ACh may stabilize the GluN1 subunit "semiclosed" conformations by inducing direct and indirect D1-to-D2 interdomain bonds. Additionally, an alternative binding site was detected, shared by the known NMDAR allosteric modulators. Experimental and computed results strongly suggest that ACh acts as a Gly-competitive, submaximal potentiating agent of the NMDAR, possibly constituting a novel chemotype for multitarget-directed drug development, e.g., to treat Alzheimer's, and it may lead to a new understanding of glutamatergic neurotransmission.

Induced fit, ensemble binding space docking and Monte Carlo simulations of MDMA 'ecstasy' and 3D pharmacophore design of MDMA derivatives on the human serotonin transporter (hSERT).

The popular recreational drug MDMA (3,4-methylenedioxy-methamphetamine) has a documented potential as a psychopharmacological clinical and research tool. This is due to its unique ability to promote reprocessing of traumatic memories, empathetic and pro-social states. Although it is established that MDMA exerts its behavioural effects via the serotonin transporter (SERT), the ligand-protein molecular interplay remains elusive. In order to shed light on the binding of MDMA and its primary congeneric entactogens (MDA, MBDB and MDAI), we first combined induced fit with Monte Carlo simulations. The computed interaction energies of the models correlated well with experimental activities (adjR2 = 0.78). Then we carried out 'ensemble binding space docking' on trajectories generated by interpolation of experimentally derived structures of the hSERT from the outward-open, and the occluded, to the inward-open states. This approach revealed low-energy alternative binding modes, suggesting high occupancy of the central site, yet considerable MDMA mobility within it, favouring the paroxetine-like orientation. Finally, we designed a pharmacophore that may be used to recognise hSERT-mediated serotonin releasers and uptake inhibitors of diverse chemical structure, identifying their active conformations and interacting residues. We conclude that the conserved amine-Asp98 ionic and edge-to-face π-π interactions are crucial to the mode of action of MDMA on the hSERT, underscoring the contributions of Tyr95 and gating residues Phe341, Tyr176 and Phe335. Amenable to experimental testing, our modelling may aid the rational design of novel entactogenic compounds and contribute to the understanding of an action mechanism, common and typical of psychotropic agents.

Mefloquine inhibits voltage dependent Nav1.4 channel by overlapping the local anaesthetic binding site.

Mefloquine constitutes a multitarget antimalaric that inhibits cation currents. However, the effect and the binding site of this compound on Na+ channels is unknown. To address the mechanism of action of mefloquine, we employed two-electrode voltage clamp recordings on Xenopus laevis oocytes, site-directed mutagenesis of the rat Na+ channel, and a combined in silico approach using Molecular Dynamics and docking protocols. We found that mefloquine: i) inhibited Nav1.4 currents (IC50 =60µM), ii) significantly delayed fast inactivation but did not affect recovery from inactivation, iii) markedly the shifted steady-state inactivation curve to more hyperpolarized potentials. The presence of the ß1 subunit significantly reduced mefloquine potency, but the drug induced a significant frequency-independent rundown upon repetitive depolarisations. Computational and experimental results indicate that mefloquine overlaps the local anaesthetic binding site by docking at a hydrophobic cavity between domains DIII and DIV that communicates the local anaesthetic binding site with the selectivity filter. This is supported by the fact that mefloquine potency significantly decreased on mutant Nav1.4 channel F1579A and significantly increased on K1237S channels. In silico this compound docked above F1579 forming stable π-π interactions with this residue. We provide structure-activity insights into how cationic amphiphilic compounds may exert inhibitory effects by docking between the local anaesthetic binding site and the selectivity filter of a mammalian Na+ channel. Our proposed synergistic cycle of experimental and computational studies may be useful for elucidating binding sites of other drugs, thereby saving in vitro and in silico resources.

Characterization of specific allosteric effects of the Na+ channel ß1 subunit on the Nav1.4 isoform.

The mechanism of inactivation of mammalian voltage-gated Na+ channels involves transient interactions between intracellular domains resulting in direct pore occlusion by the IFM motif and concomitant extracellular interactions with the ß1 subunit. Navß1 subunits constitute single-pass transmembrane proteins that form protein-protein associations with pore-forming α subunits to allosterically modulate the Na+ influx into the cell during the action potential of every excitable cell in vertebrates. Here, we explored the role of the intracellular IFM motif of rNav1.4 (skeletal muscle isoform of the rat Na+ channel) on the α-ß1 functional interaction and showed for the first time that the modulation of ß1 is independent of the IFM motif. We found that: (1) Nav1.4 channels that lack the IFM inactivation particle can undergo a "C-type-like inactivation" albeit in an ultraslow gating mode; (2) ß1 can significantly accelerate the inactivation of Nav1.4 channels in the absence of the IFM motif. Previously, we identified two residues (T109 and N110) on the ß1 subunit that disrupt the α-ß1 allosteric modulation. We further characterized the electrophysiological effects of the double alanine substitution of these residues demonstrating that it decelerates inactivation and recovery from inactivation, abolishes the modulation of steady-state inactivation and induces a current rundown upon repetitive stimulation, thus causing a general loss of function. Our results contribute to delineating the process of the mammalian Na+ channel inactivation. These findings may be relevant to the design of pharmacological strategies, targeting ß subunits to treat pathologies associated to Na+ current dysfunction.

Predicting a double mutant in the twilight zone of low homology modeling for the skeletal muscle voltage-gated sodium channel subunit beta-1 (Nav1.4 ß1).

The molecular structure modeling of the ß1 subunit of the skeletal muscle voltage-gated sodium channel (Nav1.4) was carried out in the twilight zone of very low homology. Structural significance can per se be confounded with random sequence similarities. Hence, we combined (i) not automated computational modeling of weakly homologous 3D templates, some with interfaces to analogous structures to the pore-bearing Nav1.4 α subunit with (ii) site-directed mutagenesis (SDM), as well as (iii) electrophysiological experiments to study the structure and function of the ß1 subunit. Despite the distant phylogenic relationships, we found a 3D-template to identify two adjacent amino acids leading to the long-awaited loss of function (inactivation) of Nav1.4 channels. This mutant type (T109A, N110A, herein called TANA) was expressed and tested on cells of hamster ovary (CHO). The present electrophysiological results showed that the double alanine substitution TANA disrupted channel inactivation as if the ß1 subunit would not be in complex with the α subunit. Exhaustive and unbiased sampling of "all ß proteins" (Ig-like, Ig) resulted in a plethora of 3D templates which were compared to the target secondary structure prediction. The location of TANA was made possible thanks to another "all ß protein" structure in complex with an irreversible bound protein as well as a reversible protein-protein interface (our "Rosetta Stone" effect). This finding coincides with our electrophysiological data (disrupted ß1-like voltage dependence) and it is safe to utter that the Nav1.4 α/ß1 interface is likely to be of reversible nature.

Identification of Navß1 residues involved in the modulation of the sodium channel Nav1.4.

Voltage-gated sodium channels (VGSCs) are heteromeric protein complexes that initiate action potentials in excitable cells. The voltage-gated sodium channel accessory subunit, Navß1, allosterically modulates the α subunit pore structure upon binding. To date, the molecular determinants of the interface remain unknown. We made use of sequence, knowledge and structure-based methods to identify residues critical to the association of the α and ß1 Nav1.4 subunits. The Navß1 point mutant C43A disrupted the modulation of voltage dependence of activation and inactivation and delayed the peak current decay, the recovery from inactivation, and induced a use-dependent decay upon depolarisation at 1 Hz. The Navß1 mutant R89A selectively delayed channel inactivation and recovery from inactivation and had no effect on voltage dependence or repetitive depolarisations. Navß1 mutants Y32A and G33M selectively modified the half voltage of inactivation without altering the kinetics. Despite low sequence identity, highly conserved structural elements were identified. Our models were consistent with published data and may help relate pathologies associated with VGSCs to the Navß1 subunit.