Eleonora Gianti

Flipping in the Pore: Discovery of Dual Inhibitors That Bind in Different Orientations to the Wild-Type versus the Amantadine-Resistant S31N Mutant of the Influenza A Virus M2 Proton Channel.

Influenza virus infections lead to numerous deaths and millions of hospitalizations each year. One challenge facing anti-influenza drug development is the heterogeneity of the circulating influenza viruses, which comprise several strains with variable susceptibility to antiviral drugs. For example, the wild-type (WT) influenza A viruses, such as the seasonal H1N1, tend to be sensitive to antiviral drugs, amantadine and rimantadine, while the S31N mutant viruses, such as the pandemic 2009 H1N1 (H1N1pdm09) and seasonal H3N2, are resistant to this class of drugs. Thus, drugs targeting both WT and the S31N mutant are highly desired. We report our design of a novel class of dual inhibitors along with their ion channel blockage and antiviral activities. The potency of the most active compound 11 in inhibiting WT and the S31N mutant influenza viruses is comparable with that of amantadine in inhibiting WT influenza virus. Solution NMR studies and molecular dynamics (MD) simulations of drug-M2 interactions supported our design hypothesis: namely, the dual inhibitor binds in the WT M2 channel with an aromatic group facing down toward the C-terminus, while the same drug binds in the S31N M2 channel with its aromatic group facing up toward the N-terminus. The flip-flop mode of drug binding correlates with the structure–activity relationship (SAR) and has paved the way for the next round of rational design of broad-spectrum antiviral drugs.

Hydrogen-Bonded Water Molecules in the M2 Channel of the Influenza A Virus Guide the Binding Preferences of Ammonium-Based Inhibitors

The tetrameric M2 proton channel of influenza A virus is an integral membrane protein responsible for the acidification of the viral interior. Drugs such as amantadine target the transmembrane region of wild type M2 by acting as pore blockers. However, a number of mutations affecting this domain confer drug resistance, prompting the need for alternative inhibitors. The availability of high-resolution structures of drug-bound M2, paired with computational investigations, revealed that inhibitors can bind at different sites, and provided useful insights in understanding the principles governing proton conduction. Here, we investigated by computation the energetic and geometric factors determining the relative stability of pore blockers at individual sites of different M2 strains. We found that local free energy minima along the translocation pathway of positively charged chemical species correspond to experimentally determined binding sites of inhibitors. Then, by examining the structure of water clusters hydrating each site, as well as of those displaced by binding of hydrophobic scaffolds, we predicted the binding preferences of M2 ligands. This information can be used to guide the identification of novel classes of inhibitors.

Does Proton Conduction in the Voltage-Gated H+ Channel hHv1 Involve Grotthuss-Like Hopping via Acidic Residues?

Hv1s are ubiquitous highly selective voltage-gated proton channels involved in male fertility, immunology, and the invasiveness of certain forms of breast cancer. The mechanism of proton extrusion in Hv1 is not yet understood, while it constitutes the first step toward the design of high-affinity drugs aimed at this important pharmacological target. In this contribution, we explore the details of the mechanism via an integrative approach, using classical and QM/MM molecular dynamics simulations of a monomeric hHv1 model. We propose that protons localize in three binding sites along the channel lumen, formed by three pairs of conserved negatively charged residues lining the pore: D174/E153, D112/D185, and E119/D123. Local rearrangements, involving notably a dihedral transition of F150, a conserved phenylalanine lining the permeation pathway, appear to allow protons to hop from one acidic residue to the next through a bridging water molecule. These results constitute a first attempt at rationalizing hHv1 selectivity for H+ and the role played by D112 in this process. They pave the way for further quantitative characterization of H+ transport in hHv1.

On the role of water density fluctuations in the inhibition of a proton channel

Significance Hv1, a voltage-gated proton channel, is an emerging pharmacological target implicated in many pathological conditions, including cancer and ischemic brain damage. We used the recently published experimental structure of Hv1 to generate structural models of relevant conformational states. Thermodynamic analyses of pore waters shed light on the molecular underpinnings of Hv1 druggability. We exploit this information to suggest possible optimizations of known inhibitors and identify a potential binding site located at the exit of the proton path. The resulting molecular picture paves the way for the discovery of novel Hv1 inhibitors and outlines a general approach for identifying druggable binding sites in ion channels. Hv1 is a transmembrane four-helix bundle that transports protons in a voltage-controlled manner. Its crucial role in many pathological conditions, including cancer and ischemic brain damage, makes Hv1 a promising drug target. Starting from the recently solved crystal structure of Hv1, we used structural modeling and molecular dynamics simulations to characterize the channel’s most relevant conformations along the activation cycle. We then performed computational docking of known Hv1 inhibitors, 2-guanidinobenzimidazole (2GBI) and analogs. Although salt-bridge patterns and electrostatic potential profiles are well-defined and distinctive features of activated versus nonactivated states, the water distribution along the channel lumen is dynamic and reflects a conformational heterogeneity inherent to each state. In fact, pore waters assemble into intermittent hydrogen-bonded clusters that are replaced by the inhibitor moieties upon ligand binding. The entropic gain resulting from releasing these conformationally restrained waters to the bulk solvent is likely a major contributor to the binding free energy. Accordingly, we mapped the water density fluctuations inside the pore of the channel and identified the regions of maximum fluctuation within putative binding sites. Two sites appear as outstanding: One is the already known binding pocket of 2GBI, which is accessible to ligands from the intracellular side; the other is a site located at the exit of the proton permeation pathway. Our analysis of the waters confined in the hydrophobic cavities of Hv1 suggests a general strategy for drug discovery that can be applied to any ion channel.

Computational analysis of EBNA1 “druggability” suggests novel insights for Epstein-Barr virus inhibitor design

The Epstein-Barr Nuclear Antigen 1 (EBNA1) is a critical protein encoded by the Epstein-Barr Virus (EBV). During latent infection, EBNA1 is essential for DNA replication and transcription initiation of viral and cellular genes and is necessary to immortalize primary B-lymphocytes. Nonetheless, the concept of EBNA1 as drug target is novel. Two EBNA1 crystal structures are publicly available and the first small-molecule EBNA1 inhibitors were recently discovered. However, no systematic studies have been reported on the structural details of EBNA1 “druggable” binding sites. We conducted computational identification and structural characterization of EBNA1 binding pockets, likely to accommodate ligand molecules (i.e. “druggable” binding sites). Then, we validated our predictions by docking against a set of compounds previously tested in vitro for EBNA1 inhibition (PubChem AID-2381). Finally, we supported assessments of pocket druggability by performing induced fit docking and molecular dynamics simulations paired with binding affinity predictions by Molecular Mechanics Generalized Born Surface Area calculations for a number of hits belonging to druggable binding sites. Our results establish EBNA1 as a target for drug discovery, and provide the computational evidence that active AID-2381 hits disrupt EBNA1:DNA binding upon interacting at individual sites. Lastly, structural properties of top scoring hits are proposed to support the rational design of the next generation of EBNA1 inhibitors.

Modeling androgen receptor flexibility: a binding mode hypothesis of CYP17 inhibitors/antiandrogens for prostate cancer therapy.

Prostate Cancer (PCa), a leading cause of cancer death worldwide (www.cancer.gov), is a complex malignancy where a spectrum of targets leads to a diversity of PCa forms. A widely pursued therapeutic target is the Androgen Receptor (AR). As a Steroid Hormone Receptor, AR serves as activator of transcription upon binding to androgens and plays a central role in the development of PCa. AR is a structurally flexible protein, and conformational plasticity of residues in the binding-pocket is a key to its ability to accommodate ligands from various chemical classes. Besides direct modulation of AR activity by antagonists, inhibition of cytochrome CYP17 (17α-hydroxylase/17,20-lyase), essential in androgen biosynthesis, has widely been considered an effective strategy against PCa. Interestingly, Handratta et al. (2005) discovered new, potent inhibitors of CYP17 (C-17 steroid derivatives) with pure AR antagonistic properties. Although the antiandrogenic activity of their lead compound (VN/124-1) has been experimentally proven both in vitro and in vivo, no structural data are currently available to elucidate the molecular determinants responsible for these desirable dual inhibitory properties. We implemented a Structure-based Drug Design (SBDD) approach to generate a valuable hypothesis as to the binding modes of steroidal CYP17 inhibitors/antiandrogens against the AR. To deal with the plasticity of residues buried in the Ligand Binding Domain (LBD), we developed a flexible-receptor Docking protocol based on Induced-Fit (IFD) methodology (www.schrodinger.com/). Our results constitute an ideal starting point for the rational design of next-generation analogues of CYP17 inhibitors/antiandrogens as well as an attractive tool to suggest novel chemical classes of AR antagonists.

Sites Contributing to TRPA1 Activation by the Anesthetic Propofol Identified by Photoaffinity Labeling

In addition to inducing anesthesia, propofol activates a key component of the pain pathway, the transient receptor potential ankyrin 1 ion channel (TRPA1). Recent mutagenesis studies suggested a potential activation site within the transmembrane domain, near the A-967079 cavity. However, mutagenesis cannot distinguish between protein-based and ligand-based mechanisms, nor can this site explain the complex modulation by propofol. Thus more direct approaches are required to reveal potentially druggable binding sites. Here we apply photoaffinity labeling using a propofol derivative, meta-azipropofol, for direct identification of binding sites in mouse TRPA1. We confirm that meta-azipropofol activates TRPA1 like the parent anesthetic, and identify two photolabeled residues (V954 and E969) in the S6 helix. In combination with docking to closed and open state models of TRPA1, photoaffinity labeling suggested that the A-967079 cavity is a positive modulatory site for propofol. Further, the photoaffinity labeling of E969 indicated pore block as a likely mechanism for propofol inhibition at high concentrations. The direct identification of drug-binding sites clarifies the molecular mechanisms of important TRPA1 agonists, and will facilitate drug design efforts to modulate TRPA1.

An SH2 domain model of STAT5 in complex with phospho-peptides define "STAT5 Binding Signatures".

The signal transducer and activator of transcription 5 (STAT5) is a member of the STAT family of proteins, implicated in cell growth and differentiation. STAT activation is regulated by phosphorylation of protein monomers at conserved tyrosine residues, followed by binding to phospho-peptide pockets and subsequent dimerization. STAT5 is implicated in the development of severe pathological conditions, including many cancer forms. However, nowadays a few STAT5 inhibitors are known, and only one crystal structure of the inactive STAT5 dimer is publicly available. With a view to enabling structure-based drug design, we have: (1) analyzed phospho-peptide binding pockets on SH2 domains of STAT5, STAT1 and STAT3; (2) generated a model of STAT5 bound to phospho-peptides; (3) assessed our model by docking against a class of known STAT5 inhibitors (Müller et al. in ChemBioChem 9:723–727, 2008); (4) used molecular dynamics simulations to optimize the molecular determinants responsible for binding and (5) proposed unique “Binding Signatures” of STAT5. Our results put in place the foundations to address STAT5 as a target for rational drug design, from sequence, structural and functional perspectives.

Computational Approaches to Studying Voltage-Gated Ion Channel Modulation by General Anesthetics

Voltage-gated ion channels (VGICs) are responsible for the propagation of electrical signals in excitable cells. Small-molecule modulation of VGICs affects transmission of action potentials in neurons and thus can modulate the activity of the central nervous system. For this reason, VGICs are considered key players in the medically induced state of general anesthesia. Consistently, VGICs have been shown to respond to several general anesthetics. However, in spite of extensive electrophysiological characterizations, modulation of VGICs by anesthetics is still only partially understood. Among the challenging aspects are the presence of multiple binding sites and the observation of paradoxical effects, i.e., evidence, for the same channel, of inhibition and potentiation. In this context, molecular simulations emerged in the recent past as the tool of choice to complement electrophysiology studies with a microscopic picture of binding and allosteric regulation. In this chapter, we describe the most effective computational techniques to study VGIC modulation by general anesthetics. We start by reviewing the VGIC conduction cycle, the corresponding set of channel conformations, and the approaches used to model them. We then review the most successful strategies to identify binding sites and estimate binding affinities.

Molecular Insights into TRPV1 Polymodal Activation: Is Allosteric Coupling between the Toxin and the Vanilloid Binding Sites Mediated by Annular Lipids?

1949-Plat Molecular Insights into TRPV1 Polymodal Activation: Is Allosteric Coupling between the Toxin and the Vanilloid Binding Sites Mediated by Annular Lipids? Eleonora Gianti1, Michael Klein1, Tibor Rohács2, Vincenzo Carnevale1. Chemistry and the ICMS, Temple University, Philadelphia, PA, USA, Department of Pharmacology, Physiology & Neuroscience, New Jersey Medical School, Rutgers University, Newark, NJ, USA. The TRPV1 channel is responsible for conducting cations through the cell membrane in response to a variety of stimuli, amongst which noxious heat and chemical ligands, such as the vanilloid compounds capsaicin and its ultrapotent agonist resiniferatoxin (RTX). Predominantly expressed in sensory neurons and nociceptive fibers, TRPV1 is involved in important cellular functions, including heat-sensation and pain. One unresolved question about TRPV1 concerns activation by stimuli acting from the extracellular milieu, such as the binding of the double-knot toxin (Dk/Tx) from spider venom. How is this stimulus coupled to the binding of RTX at the vanilloid site? How do Dk/Tx and RTX synergistically act on the gate? Interestingly, both these biomolecular ligands participate in the formation of the vanilloid-lipid-channel-toxin quadripartite complex, as shown by recent cryoEM experiments. We performed molecular dynamics (MD) simulations and confirmed that annular lipids are permanently bound to and stabilize the open state. While lipid head-groups interact with Dk/Tx, the tail-groups are simultaneously in contact with RTX and with a side chain crucial for the movement of the S4-S5 linker. Importantly, mutations at this position were shown to selectively impair capsaicin ability to activate the channel without appreciably affecting ligand binding. In summary, our results shed light on the mechanism by which vanilloid binding is transduced into gating motion and rationalize the experimentally observed allosteric coupling between Dk/Tx and RTX.