SENSAAS-Flex: a joint optimization approach for aligning 3D shapes and exploring the molecular conformation space.

MOTIVATION: Popular shape-based alignment methods handle molecular flexibility by utilizing conformational ensembles to select the most fitted conformer. However, the initial conformer library generation step is computationally intensive and limiting to the overall alignment process. In this work, we describe a method to perform flexible alignment of two molecular shapes by optimizing the 3D conformation. SENSAAS-Flex, an add-on to the SENSAAS tool, is able to proceed from a limited set of initial conformers through an iterative process where additional conformational optimizations are made at the substructure level and constrained by the target shape. RESULTS: In self- and cross-alignment experiments, SENSAAS-Flex is able to reproduce the crystal structure geometry of ligands of the AstraZeneca Molecule Overlay Test set and PDBbind refined dataset. Our study shows that the point-based representation of molecular surfaces is appropriate in terms of shape constraint to sample the conformational space and perform flexible molecular alignments. AVAILABILITY AND IMPLEMENTATION: The documentation and source code are available at https://chemoinfo.ipmc.cnrs.fr/Sensaas-flex/sensaas-flex-main.tar.gz.

Pharmacological activation of PIEZO1 in human red blood cells prevents Plasmodium falciparum invasion.

An inherited gain-of-function variant (E756del) in the mechanosensitive cationic channel PIEZO1 was shown to confer a significant protection against severe malaria. Here, we demonstrate in vitro that human red blood cell (RBC) infection by Plasmodium falciparum is prevented by the pharmacological activation of PIEZO1. Yoda1 causes an increase in intracellular calcium associated with rapid echinocytosis that inhibits RBC invasion, without affecting parasite intraerythrocytic growth, division or egress. Notably, Yoda1 treatment significantly decreases merozoite attachment and subsequent RBC deformation. Intracellular Na+/K+ imbalance is unrelated to the mechanism of protection, although delayed RBC dehydration observed in the standard parasite culture medium RPMI/albumax further enhances the resistance to malaria conferred by Yoda1. The chemically unrelated Jedi2 PIEZO1 activator similarly causes echinocytosis and RBC dehydration associated with resistance to malaria invasion. Spiky outward membrane projections are anticipated to reduce the effective surface area required for both merozoite attachment and internalization upon pharmacological activation of PIEZO1. Globally, our findings indicate that the loss of the typical biconcave discoid shape of RBCs, together with an altered optimal surface to volume ratio, induced by PIEZO1 pharmacological activation prevent efficient P. falciparum invasion.

Mambalgin-1 pain-relieving peptide locks the hinge between α4 and α5 helices to inhibit rat acid-sensing ion channel 1a.

Acid-sensing ion channels (ASICs) are proton-gated cationic channels involved in pain and other processes, underscoring the potential therapeutic value of specific inhibitors such as the three-finger toxin mambalgin-1 (Mamb-1) from snake venom. A low-resolution structure of the human-ASIC1a/Mamb-1 complex obtained by cryo-electron microscopy has been recently reported, implementing the structure of the chicken-ASIC1/Mamb-1 complex previously published. Here we combine structure-activity relationship of both the rat ASIC1a channel and the Mamb-1 toxin with a molecular dynamics simulation to obtain a detailed picture at the level of side-chain interactions of the binding of Mamb-1 on rat ASIC1a channels and of its inhibition mechanism. Fingers I and II of Mamb-1 but not the core of the toxin are required for interaction with the thumb domain of ASIC1a, and Lys-8 of finger I potentially interacts with Tyr-358 in the thumb domain. Mamb-1 does not interfere directly with the pH sensor as previously suggested, but locks by several contacts a key hinge between α4 and α5 helices in the thumb domain of ASIC1a to prevent channel opening. Our results provide an improved model of inhibition of mammalian ASIC1a channels by Mamb-1 and clues for further development of optimized ASIC blockers.

sensaas: Shape-based Alignment by Registration of Colored Point-based Surfaces.

sensaas is a tool developed for aligning and comparing molecular shapes and sub-shapes. Alignment is obtained by registration of 3D point-based representations of the van der Waals surface. The method uses local properties of the shape to identify the correspondence relationships between two point clouds containing up to several thousand colored (labeled) points. Our rigid-body superimposition method follows a two-stage approach. An initial alignment is obtained by matching pose-invariant local 3D descriptors, called FPFH, of the input point clouds. This stage provides a global superimposition of the molecular surfaces, without any knowledge of their initial pose in 3D space. This alignment is then refined by optimizing the matching of colored points. In our study, each point is colored according to its closest atom, which itself belongs to a user defined physico-chemical class. Finally, sensaas provides an alignment and evaluates the molecular similarity by using Tversky coefficients. To assess the efficiency of this approach, we tested its ability to reproduce the superimposition of X-ray structures of the benchmarking AstraZeneca (AZ) data set and, compared its results with those generated by the two shape-alignment approaches shaep and shafts. We also illustrated submatching properties of our method with respect to few substructures and bioisosteric fragments. The code is available upon request from the authors (demo version at https://chemoinfo.ipmc.cnrs.fr/SENSAAS).

Inhibition of Patched Drug Efflux Increases Vemurafenib Effectiveness against Resistant BrafV600E Melanoma.

Melanoma patients harboring the BRAFV600E mutation are treated with vemurafenib. Almost all of them ultimately acquire resistance, leading to disease progression. Here, we find that a small molecule from a marine sponge, panicein A hydroquinone (PAH), overcomes resistance of BRAFV600E melanoma cells to vemurafenib, leading to tumor elimination in corresponding human xenograft models in mice. We report the synthesis of PAH and demonstrate that this compound inhibits the drug efflux activity of the Hedgehog receptor, Patched. Our SAR study allowed identifying a key pharmacophore responsible for this activity. We showed that Patched is strongly expressed in metastatic samples from a cohort of melanoma patients and is correlated with decreased overall survival. Patched is a multidrug transporter that uses the proton motive force to efflux drugs. This makes its function specific to cancer cells, thereby avoiding toxicity issues that are commonly observed with inhibitors of ABC multidrug transporters. Our data provide strong evidence that PAH is a highly promising lead for the treatment of vemurafenib resistant BRAFV600E melanoma.

Piezo Ion Channels in Cardiovascular Mechanobiology.

Mechanotransduction has a key role in vascular development, physiology, and disease states. Piezo1 is a mechanosensitive (MS) nonselective cationic channel that occurs in endothelial and vascular smooth muscle cells. It is activated by shear stress associated with increases in local blood flow, as well as by cell membrane stretch upon elevation of blood pressure. Here, we briefly review the pharmacological modulators of Piezo and discuss current understanding of the role of Piezo1 in vascular mechanobiology and associated clinical disorders, such as atherosclerosis and hypertension. Ultimately, we believe that this research will help identify novel therapeutic strategies for the treatment of vascular diseases.

Mammalian Mechanoelectrical Transduction: Structure and Function of Force-Gated Ion Channels.

The conversion of force into an electrical cellular signal is mediated by the opening of different types of mechanosensitive ion channels (MSCs), including TREK/TRAAK K2P channels, Piezo1/2, TMEM63/OSCA, and TMC1/2. Mechanoelectrical transduction plays a key role in hearing, balance, touch, and proprioception and is also implicated in the autonomic regulation of blood pressure and breathing. Thus, dysfunction of MSCs is associated with a variety of inherited and acquired disease states. Significant progress has recently been made in identifying these channels, solving their structure, and understanding the gating of both hyperpolarizing and depolarizing MSCs. Besides prototypical activation by membrane tension, additional gating mechanisms involving channel curvature and/or tethered elements are at play.

Structure and function of polycystins: insights into polycystic kidney disease.

Mutations in the polycystins PC1 or PC2 cause autosomal dominant polycystic kidney disease (ADPKD), which is characterized by the formation of fluid-filled renal cysts that disrupt renal architecture and function, ultimately leading to kidney failure in the majority of patients. Although the genetic basis of ADPKD is now well established, the physiological function of polycystins remains obscure and a matter of intense debate. The structural determination of both the homomeric PC2 and heteromeric PC1-PC2 complexes, as well as the electrophysiological characterization of PC2 in the primary cilium of renal epithelial cells, provided new valuable insights into the mechanisms of ADPKD pathogenesis. Current findings indicate that PC2 can function independently of PC1 in the primary cilium of renal collecting duct epithelial cells to form a channel that is mainly permeant to monovalent cations and is activated by both membrane depolarization and an increase in intraciliary calcium. In addition, PC2 functions as a calcium-activated calcium release channel at the endoplasmic reticulum membrane. Structural studies indicate that the heteromeric PC1-PC2 complex comprises one PC1 and three PC2 channel subunits. Surprisingly, several positively charged residues from PC1 occlude the ionic pore of the PC1-PC2 complex, suggesting that pathogenic polycystin mutations might cause ADPKD independently of an effect on channel permeation. Emerging reports of novel structural and functional findings on polycystins will continue to elucidate the molecular basis of ADPKD.

Data Sets Representative of the Structures and Experimental Properties of FDA-Approved Drugs.

Presented here are several data sets that gather information collected from the labels of the FDA approved drugs: their molecular structures and those of the described active metabolites, their associated pharmacokinetics and pharmacodynamics data, and the history of their marketing authorization by the FDA. To date, 1852 chemical structures have been identified with a molecular weight less than 2000 of which 492 are or have active metabolites. To promote the sharing of data, the original web server was upgraded for browsing the database and downloading the data sets (http://chemoinfo.ipmc.cnrs.fr/edrug3d). It is believed that the multidimensional chemistry-oriented collections are an essential resource for a thorough analysis of the current drug chemical space. The data sets are envisioned as being used in a wide range of endeavors that include drug repurposing, drug design, privileged structures analyses, structure-activity relationship studies, and improving of absorption, distribution, metabolism, and elimination predictive models.

Computational analysis of calculated physicochemical and ADMET properties of protein-protein interaction inhibitors.

The modulation of PPIs by low molecular weight chemical compounds, particularly by orally bioavailable molecules, would be very valuable in numerous disease indications. However, it is known that PPI inhibitors (iPPIs) tend to have properties that are linked to poor Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) and in some cases to poor clinical outcomes. Previously reported in silico analyses of iPPIs have essentially focused on physicochemical properties but several other ADMET parameters would be important to assess. In order to gain new insights into the ADMET properties of iPPIs, computations were carried out on eight datasets collected from several databases. These datasets involve compounds targeting enzymes, GPCRs, ion channels, nuclear receptors, allosteric modulators, oral marketed drugs, oral natural product-derived marketed drugs and iPPIs. Several trends are reported that should assist the design and optimization of future PPI inhibitors, either for drug discovery endeavors or for chemical biology projects.

Fragment-based discovery of a new family of non-peptidic small-molecule cyclophilin inhibitors with potent antiviral activities.

Cyclophilins are peptidyl-prolyl cis/trans isomerases (PPIase) that catalyse the interconversion of the peptide bond at proline residues. Several cyclophilins play a pivotal role in the life cycle of a number of viruses. The existing cyclophilin inhibitors, all derived from cyclosporine A or sanglifehrin A, have disadvantages, including their size, potential for side effects unrelated to cyclophilin inhibition and drug-drug interactions, unclear antiviral spectrum and manufacturing issues. Here we use a fragment-based drug discovery approach using nucleic magnetic resonance, X-ray crystallography and structure-based compound optimization to generate a new family of non-peptidic, small-molecule cyclophilin inhibitors with potent in vitro PPIase inhibitory activity and antiviral activity against hepatitis C virus, human immunodeficiency virus and coronaviruses. This family of compounds has the potential for broad-spectrum, high-barrier-to-resistance treatment of viral infections.

Comparison of the ligand binding site of CYP2C8 with CYP26A1 and CYP26B1: a structural basis for the identification of new inhibitors of the retinoic acid hydroxylases.

The CYP26s are responsible for metabolizing retinoic acid and play an important role in maintaining homeostatic levels of retinoic acid. Given the ability of CYP2C8 to metabolize retinoic acid, we evaluated the potential for CYP2C8 inhibitors to also inhibit CYP26. In vitro assays were used to evaluate the inhibition potencies of CYP2C8 inhibitors against CYP26A1 and CYP26B1. Using tazarotenic acid as a substrate for CYP26, IC50 values for 17 inhibitors of CYP2C8 were determined for CYP26A1 and CYP26B1, ranging from ∼20 nM to 100 µM, with a positive correlation observed between IC50s for CYP2C8 and CYP26A1. An evaluation of IC50's versus in vivo Cmax values suggests that inhibitors such as clotrimazole or fluconazole may interact with CYP26 at clinically relevant concentrations and may alter levels of retinoic acid. These findings provide insight into drug interactions resulting in elevated retinoic acid concentrations and expand upon the pharmacophore of CYP26 inhibition.

Identification of Tazarotenic Acid as the First Xenobiotic Substrate of Human Retinoic Acid Hydroxylase CYP26A1 and CYP26B1.

Cytochrome P450 (CYP) 26A1 and 26B1 are heme-containing enzymes responsible for metabolizing all-trans retinoic acid (at-RA). No crystal structures have been solved, and therefore homology models that provide structural information are extremely valuable for the development of inhibitors of cytochrome P450 family 26 (CYP26). The objectives of this study were to use homology models of CYP26A1 and CYP26B1 to characterize substrate binding characteristics, to compare structural aspects of their active sites, and to support the role of CYP26 in the metabolism of xenobiotics. Each model was verified by dockingat-RA in the active site and comparing the results to known metabolic profiles ofat-RA. The models were then used to predict the metabolic sites of tazarotenic acid with results verified by in vitro metabolite identification experiments. The CYP26A1 and CYP26B1 homology models predicted that the benzothiopyranyl moiety of tazarotenic acid would be oriented toward the heme of each enzyme and suggested that tazarotenic acid would be a substrate of CYP26A1 and CYP26B1. Metabolite identification experiments indicated that CYP26A1 and CYP26B1 oxidatively metabolized tazarotenic acid on the predicted moiety, with in vitro rates of metabolite formation by CYP26A1 and CYP26B1 being the highest across a panel of enzymes. Molecular analysis of the active sites estimated the active-site volumes of CYP26A1 and CYP26B1 to be 918 Å(3)and 977 Å(3), respectively. Overall, the homology models presented herein describe the enzyme characteristics leading to the metabolism of tazarotenic acid by CYP26A1 and CYP26B1 and support a potential role for the CYP26 enzymes in the metabolism of xenobiotics.

Stability of the Plasmodium falciparum AMA1-RON2 Complex Is Governed by the Domain II (DII) Loop.

Plasmodium falciparum is an obligate intracellular protozoan parasite that employs a highly sophisticated mechanism to access the protective environment of the host cells. Key to this mechanism is the formation of an electron dense ring at the parasite-host cell interface called the Moving Junction (MJ) through which the parasite invades. The MJ incorporates two key parasite components: the surface protein Apical Membrane Antigen 1 (AMA1) and its receptor, the Rhoptry Neck Protein (RON) complex, the latter one being targeted to the host cell membrane during invasion. Crystal structures of AMA1 have shown that a partially mobile loop, termed the DII loop, forms part of a deep groove in domain I and overlaps with the RON2 binding site. To investigate the mechanism by which the DII loop influences RON2 binding, we measured the kinetics of association and dissociation and binding equilibria of a PfRON2sp1 peptide with both PfAMA1 and an engineered form of PfAMA1 where the flexible region of the DII loop was replaced by a short Gly-Ser linker (ΔDII-PfAMA1). The reactions were tracked by fluorescence anisotropy as a function of temperature and concentration and globally fitted to acquire the rate constants and corresponding thermodynamic profiles. Our results indicate that both PfAMA1 constructs bound to the PfRON2sp1 peptide with the formation of one intermediate in a sequential reversible reaction: A↔B↔C. Consistent with Isothermal Titration Calorimetry measurements, final complex formation was enthalpically driven and slightly entropically unfavorable. Importantly, our experimental data shows that the DII loop lengthened the complex half-life time by 18-fold (900 s and 48 s at 25°C for Pf and ΔDII-Pf complex, respectively). The longer half-life of the Pf complex appeared to be driven by a slower dissociation process. These data highlight a new influential role for the DII loop in kinetically locking the functional binary complex to enable host cell invasion.

Mambalgin-1 Pain-relieving Peptide, Stepwise Solid-phase Synthesis, Crystal Structure, and Functional Domain for Acid-sensing Ion Channel 1a Inhibition.

Mambalgins are peptides isolated from mamba venom that specifically inhibit a set of acid-sensing ion channels (ASICs) to relieve pain. We show here the first full stepwise solid phase peptide synthesis of mambalgin-1 and confirm the biological activity of the synthetic toxin both in vitro and in vivo. We also report the determination of its three-dimensional crystal structure showing differences with previously described NMR structures. Finally, the functional domain by which the toxin inhibits ASIC1a channels was identified in its loop II and more precisely in the face containing Phe-27, Leu-32, and Leu-34 residues. Moreover, proximity between Leu-32 in mambalgin-1 and Phe-350 in rASIC1a was proposed from double mutant cycle analysis. These data provide information on the structure and on the pharmacophore for ASIC channel inhibition by mambalgins that could have therapeutic value against pain and probably other neurological disorders.

Computational and biophysical approaches to protein-protein interaction inhibition of Plasmodium falciparum AMA1/RON2 complex.

Invasion of the red blood cell by Plasmodium falciparum parasites requires formation of an electron dense circumferential ring called the Moving Junction (MJ). The MJ is anchored by a high affinity complex of two parasite proteins: Apical Membrane Antigen 1 (PfAMA1) displayed on the surface of the parasite and Rhoptry Neck Protein 2 that is discharged from the parasite and imbedded in the membrane of the host cell. Structural studies of PfAMA1 revealed a conserved hydrophobic groove localized to the apical surface that coordinates RON2 and invasion inhibitory peptides. In the present work, we employed computational and biophysical methods to identify competitive P. falciparum AMA1-RON2 inhibitors with the goal of exploring the 'druggability' of this attractive antimalarial target. A virtual screen followed by molecular docking with the PfAMA1 crystal structure was performed using an eight million compound collection that included commercial molecules, the ChEMBL malaria library and approved drugs. The consensus approach resulted in the selection of inhibitor candidates. We also developed a fluorescence anisotropy assay using a modified inhibitory peptide to experimentally validate the ability of the selected compounds to inhibit the AMA1-RON2 interaction. Among those, we identified one compound that displayed significant inhibition. This study offers interesting clues to improve the throughput and reliability of screening for new drug leads.

Plasmodium falciparum CTP:phosphocholine cytidylyltransferase possesses two functional catalytic domains and is inhibited by a CDP-choline analog selected from a virtual screening.

Phosphatidylcholine is the major lipid component of the malaria parasite membranes and is required for parasite multiplication in human erythrocytes. Plasmodium falciparum CTP:phosphocholine cytidylyltransferase (PfCCT) is the rate-limiting enzyme of the phosphatidylcholine biosynthesis pathway and thus considered as a potential antimalarial target. In contrast to its mammalian orthologs, PfCCT contains a duplicated catalytic domain. Here, we show that both domains are catalytically active with similar kinetic parameters. A virtual screening strategy allowed the identification of a drug-size molecule competitively inhibiting the enzyme. This compound also prevented phosphatidylcholine biosynthesis in parasites and exerted an antimalarial effect. This study constitutes the first step towards a rationalized design of future new antimalarial agents targeting PfCCT.

Binding site and inhibitory mechanism of the mambalgin-2 pain-relieving peptide on acid-sensing ion channel 1a.

Acid-sensing ion channels (ASICs) are neuronal proton-gated cation channels associated with nociception, fear, depression, seizure, and neuronal degeneration, suggesting roles in pain and neurological and psychiatric disorders. We have recently discovered black mamba venom peptides called mambalgin-1 and mambalgin-2, which are new three-finger toxins that specifically inhibit with the same pharmacological profile ASIC channels to exert strong analgesic effects in vivo. We now combined bioinformatics and functional approaches to uncover the molecular mechanism of channel inhibition by the mambalgin-2 pain-relieving peptide. Mambalgin-2 binds mainly in a region of ASIC1a involving the upper part of the thumb domain (residues Asp-349 and Phe-350), the palm domain of an adjacent subunit, and the ß-ball domain (residues Arg-190, Asp-258, and Gln-259). This region overlaps with the acidic pocket (pH sensor) of the channel. The peptide exerts both stimulatory and inhibitory effects on ASIC1a, and we propose a model where mambalgin-2 traps the channel in a closed conformation by precluding the conformational change of the palm and ß-ball domains that follows proton activation. These data help to understand inhibition by mambalgins and provide clues for the development of new optimized blockers of ASIC channels.

Silencing of the tandem pore domain halothane-inhibited K+ channel 2 (THIK2) relies on combined intracellular retention and low intrinsic activity at the plasma membrane.

The tandem pore domain halothane-inhibited K(+) channel 1 (THIK1) produces background K(+) currents. Despite 62% amino acid identity with THIK1, THIK2 is not active upon heterologous expression. Here, we show that this apparent lack of activity is due to a unique combination of retention in the endoplasmic reticulum and low intrinsic channel activity at the plasma membrane. A THIK2 mutant containing a proline residue (THIK2-A155P) in its second inner helix (M2) produces K(+)-selective currents with properties similar to THIK1, including inhibition by halothane and insensitivity to extracellular pH variations. Another mutation in the M2 helix (I158D) further increases channel activity and affects current kinetics. We also show that the cytoplasmic amino-terminal region of THIK2 (Nt-THIK2) contains an arginine-rich motif (RRSRRR) that acts as a retention/retrieval signal. Mutation of this motif in THIK2 induces a relocation of the channel to the plasma membrane, resulting in measurable currents, even in the absence of mutations in the M2 helix. Cell surface delivery of a Nt-THIK2-CD161 chimera is increased by mutating the arginines of the retention motif but also by converting the serine embedded in this motif to aspartate, suggesting a phosphorylation-dependent regulation of THIK2 trafficking.

Biochemical characterization of Plasmodium falciparum CTP:phosphoethanolamine cytidylyltransferase shows that only one of the two cytidylyltransferase domains is active.

The intra-erythrocytic proliferation of the human malaria parasite Plasmodium falciparum requires massive synthesis of PE (phosphatidylethanolamine) that together with phosphatidylcholine constitute the bulk of the malaria membrane lipids. PE is mainly synthesized de novo by the CDP:ethanolamine-dependent Kennedy pathway. We previously showed that inhibition of PE biosynthesis led to parasite death. In the present study we characterized PfECT [P. falciparum CTP:phosphoethanolamine CT (cytidylyltransferase)], which we identified as the rate-limiting step of the PE metabolic pathway in the parasite. The cellular localization and expression of PfECT along the parasite life cycle were studied using polyclonal antibodies. Biochemical analyses showed that the enzyme activity follows Michaelis-Menten kinetics. PfECT is composed of two CT domains separated by a linker region. Activity assays on recombinant enzymes upon site-directed mutagenesis revealed that the N-terminal CT domain was the only catalytically active domain of PfECT. Concordantly, three-dimensional homology modelling of PfECT showed critical amino acid differences between the substrate-binding sites of the two CT domains. PfECT was predicted to fold as an intramolecular dimer suggesting that the inactive C-terminal domain is important for dimer stabilization. Given the absence of PE synthesis in red blood cells, PfECT represents a potential antimalarial target opening the way for a rational conception of bioactive compounds.