Water-mediated Hydrogen Bond Network Research Articles

Graph-based algorithms to dissect long-distance water-mediated H-bond networks for conformational couplings in GPCRs.

Changes in structure and dynamics elicited by agonist ligand binding at the extracellular side of G protein coupled receptors (GPCRs) must be relayed to the cytoplasmic G protein binding side of the receptors. To decipher the role of water-mediated hydrogen-bond networks in this relay mechanism, we have developed graph-based algorithms and analysis methodologies applicable to datasets of static structures of distinct GPCRs. For a reference dataset of static structures of bovine rhodopsin solved at the same resolution, we show that graph analyses capture the internal protein-water hydrogen-bond network. The extended analyses of static structures of rhodopsins and opioid receptors suggest a relay mechanism whereby inactive receptors have in place much of the internal core hydrogen-bond network required for long-distance relay of structural change, with extensive local H-bond clusters observed in structures solved at high resolution and with internal water molecules.

Water-Mediated Hydrogen Bond Network Drives Highly Crystalline Structure Formation of Crown Ether-Based Covalent Organic Framework for Sr Adsorption.

Covalent organic frameworks (COFs) with crown ether units have drawn great attention due to their potential applications in adsorption, catalysis, and sensing. However, employing crown ethers to construct COFs is still challenging in light of the flexible nature of macrocycles. Here, a highly crystalline one-dimensional covalent organic framework (1D-18C6-COF) with crown ether units on the ribbon edge was synthesized. The water-mediated hydrogen bond network and π-π stacking hold the 1D COF ribbons together. The combination of experimental and DFT studies demonstrated that the hydrogen bond network plays a crucial role in the structure crystallinity. The 1D-18C6-COF was applied as an adsorbent for strontium, and it exhibited rapid kinetics with good selectivity. In the competitive adsorption experiment, a separation factor of 1900 was achieved, representing one of the largest values for cesium/strontium separation. This work provides new insights into the design and functional exploration of crystalline COFs with flexible units.

A Proteorhodopsin-Related Photosensor Expands the Repertoire of Structural Motifs Employed by Sensory Rhodopsins.

Microbial rhodopsins are light-activated retinal-binding membrane proteins that perform a variety of ion transport and photosensory functions. They display several cases of convergent evolution where the same function is present in unrelated or very distant protein groups. Here we report another possible case of such convergent evolution, describing the biophysical properties of a new group of sensory rhodopsins. The first representative of this group was identified in 2004 but none of the members had been expressed and characterized. The well-studied haloarchaeal sensory rhodopsins interacting with methyl-accepting Htr transducers are close relatives of the halobacterial proton pump bacteriorhodopsin. In contrast, the sensory rhodopsins we describe here are relatives of proteobacterial proton pumps, proteorhodopsins, but appear to interact with Htr-like transducers likewise, even though they do not conserve the residues important for the interaction of haloarchaeal sensory rhodopsins with their transducers. The new sensory rhodopsins display many unusual amino acid residues, including those around the retinal chromophore; most strikingly, a tyrosine in place of a carboxyl counterion of the retinal Schiff base on helix C. To characterize their unique sequence motifs, we augment the spectroscopy and biochemistry data by structural modeling of the wild-type and three mutants. Taken together, the experimental data, bioinformatics sequence analyses, and structural modeling suggest that the tyrosine/aspartate complex counterion contributes to a complex water-mediated hydrogen-bonding network that couples the protonated retinal Schiff base to an extracellular carboxylic dyad.

A water-mediated hydrogen-bond network facilitates an exceptionally strong preference for asparagine in the S2 subsite of dipeptidyl peptidase 7 from Stenotrophomonas maltophilia

A water-mediated hydrogen-bond network facilitates an exceptionally strong preference for asparagine in the S2 subsite of dipeptidyl peptidase 7 from <i>Stenotrophomonas maltophilia</i>

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes.

Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

Mechanisms by which water-mediated hydrogen-bond networks govern pH-dependent reactions at membrane interfaces.

pH-dependent reactions at membrane interfaces are used by cells to maintain pH homeostasis and to ensure normal cell function. Membrane-bound proteins and peptide systems that sense pH typically bind protons using a few titratable sidechains that communicate with the bulk. The identity of these pH-sensing sidechains is central to hypotheses about protein reaction mechanisms and for the rational design of proton detectors. But, as protein groups hypothesized to sense pH in proteins and peptide systems are part of dynamic water-mediated networks, evaluating their pH-sensing functionality may prove difficult. We developed efficient graph-based algorithms and methodologies to evaluate dynamic hydrogen-bond networks of protein and peptide systems known to have a pH-sensing functionality, and to identify sites of hydrogen-bond networks where protons are likely to bind. We find that the assembly of dynamic water-mediated hydrogen-bond networks that include potential proton-sensing groups is a key event along the reaction coordinates of pH-dependent membrane proteins and membrane-bound peptides, and implement graph-based algorithms to characterize the dynamics of these water networks. We further find that polar and basic protein sidechains largely govern the dynamics of the water-mediated hydrogen-bond networks of proton-binding groups, and help maintain acidic sidechains within clusters that collectively can bind and store protons. Taken together, the ensemble of hydrogen-bond network computations we performed for pH-coupled proteins and peptide systems provide a framework to decipher the roles of dynamic water-mediated hydrogen-bond networks in proton binding at membrane interfaces. Research was supported in part by the European Union's Horizon 2020 Research and Innovation Program under the Marie Sklodowska-Curie grant agreement No 860592, Innovative Training Network ‘Proton and proton-coupled transport’, and by computing time from the Physics Department of the Freie Universität Berlin and from the Forschungszentrum Jülich (JURECA-DC, allocation PHDPORES).

Mechanistic Insights of Polyphenolic Compounds from Rosemary Bound to Their Protein Targets Obtained by Molecular Dynamics Simulations and Free-Energy Calculations.

Rosemary represents an important medicinal plant that has been attributed with various health-promoting properties, especially antioxidative, anti-inflammatory, and anticarcinogenic activities. Carnosic acid, carnosol, and rosmanol, as well as the phenolic acid ester rosmarinic acid, are the main compounds responsible for these actions. In our earlier research, we carried out an inverse molecular docking at the proteome scale to determine possible protein targets of the mentioned compounds. Here, we subjected the previously identified ligand-protein complexes with HIV-1 protease, K-RAS, and factor X to molecular dynamics simulations coupled with free-energy calculations. We observed that carnosic acid and rosmanol act as viable binders of the HIV-1 protease. In addition, carnosol represents a potential binder of the oncogene protein K-RAS. On the other hand, rosmarinic acid was characterized as a weak binder of factor X. We also emphasized the importance of water-mediated hydrogen-bond networks in stabilizing the binding conformation of the studied polyphenols, as well as in mechanistically explaining their promiscuous nature.

Melamine-cored glucosides for membrane protein solubilization and stabilization: importance of water-mediated intermolecular hydrogen bonding in detergent performance.

Membrane proteins play essential roles in a number of biological processes, and their structures are important in elucidating such processes at the molecular level and also for rational drug design and development. Membrane protein structure determination is notoriously challenging compared to that of soluble proteins, due largely to the inherent instability of their structures in non-lipid environments. Micelles formed by conventional detergents have been widely used for membrane protein manipulation, but they are suboptimal for long-term stability of membrane proteins, making downstream characterization difficult. Hence, there is an unmet need for the development of new amphipathic agents with enhanced efficacy for membrane protein stabilization. In this study, we designed and synthesized a set of glucoside amphiphiles with a melamine core, denoted melamine-cored glucosides (MGs). When evaluated with four membrane proteins (two transporters and two G protein-coupled receptors), MG-C11 conferred notably enhanced stability compared to the commonly used detergents, DDM and LMNG. These promising findings are mainly attributed to a unique feature of the MGs, i.e., the ability to form dynamic water-mediated hydrogen-bond networks between detergent molecules, as supported by molecular dynamics simulations. Thus, MG-C11 is the first example of a non-peptide amphiphile capable of forming intermolecular hydrogen bonds within a protein-detergent complex environment. Detergent micelles formed via a hydrogen-bond network could represent the next generation of highly effective membrane-mimetic systems useful for membrane protein structural studies.

Computational investigation of functional water molecules in GPCRs bound to G protein or arrestin.

G protein-coupled receptors (GPCRs) are membrane proteins constituting the largest family of drug targets. The activated GPCR binds either the heterotrimeric G proteins or arrestin through its activation cycle. Water molecules have been reported to play a role in GPCR activation. Nevertheless, reported studies are focused on the hydrophobic helical bundle region. How water molecules function in GPCR bound either G protein or arrestin is rarely studied. To address this issue, we carried out computational studies on water molecules in both GPCR/G protein complexes and GPCR/arrestin complexes. Using inhomogeneous fluid theory (IFT), we locate all possible hydration sites in GPCRs binding either to G protein or arrestin. We observe that the number of water molecules on the interaction surface between GPCRs and signal proteins are correlated with the insertion depths of the α5-helix from G-protein or "finger loop" from arrestin in GPCRs. In three out of the four simulation pairs, the interfaces of Rhodopsin, M2R and NTSR1 in the G protein-associated systems show more water-mediated hydrogen-bond networks when compared to these in arrestin-associated systems. This reflects that more functionally relevant water molecules may probably be attracted in G protein-associated structures than that in arrestin-associated structures. Moreover, we find the water-mediated interaction networks throughout the NPxxY region and the orthosteric pocket, which may be a key for GPCR activation. Reported studies show that non-biased agonist, which can trigger both GPCR-G protein and GPCR-arrestin activation signal, can result in pharmacologically toxicities. Our comprehensive studies of the hydration sites in GPCR/G protein complexes and GPCR/arrestin complexes may provide important insights in the design of G-protein biased agonists.

The Solvation of the E. coli CheY Phosphorylation Site Mapped by XFMS.

The Escherichia coli CheY protein belongs to a large bacterial response regulator superfamily. X-ray hydroxy radical foot-printing with mass spectroscopy (XFMS) has shown that allosteric activation of CheY by its motor target triggers a concerted internalization of aromatic sidechains. We reanalyzed the XFMS data to compare polar versus non-polar CheY residue positions. The polar residues around and including the 57D phosphorylated site had an elevated hydroxy radical reactivity. Bioinformatic measures revealed that a water-mediated hydrogen bond network connected this ring of residues with the central 57D. These residues solvated 57D to energetically stabilize the apo-CheY fold. The abundance of these reactive residues was reduced upon activation. This result was supported by the bioinformatics and consistent with the previously reported activation-induced increase in core hydrophobicity. It further illustrated XFMS detection of structural waters. Direct contacts between the ring residues and the phosphorylation site would stabilize the aspartyl phosphate. In addition, we report that the ring residue, 18R, is a constant central node in the 57D solvation network and that 18R non-polar substitutions determine CheY diversity as assessed by its evolutionary trace in bacteria with well-studied chemotaxis. These results showcase the importance of structured water dynamics for phosphorylation-mediated signal transduction.

Conformational preferences of triantennary and tetraantennary hybrid N-glycans in aqueous solution: Insights from 20 μs long atomistic molecular dynamic simulations

In the current study, we have investigated the conformational dynamics of a triantennary (N-glycan1) and tetraantennary (N-glycan2) hybrid N-glycans found on the surface of the HIV glycoprotein using 20 μs long all-atom molecular dynamics (MD) simulations. The main objective of the present study is to elucidate the influence of adding a complex branch on the overall glycan structural dynamics. Our investigation suggests that the average RMSD value increases when a complex branch is added to N-glycan1. However, the RMSD distribution is relatively wider in the case of N-glycan1 compared to N-glycan2, which indicates that multiple complex branches restrict the conformational variability of glycans. A similar observation is obtained from the principal component analysis of both glycans. All the puckering states (4C1 to 1C4) of each monosaccharide except mannose are sampled in our simulations, although the 4C1 chair form is energetically more favorable than 1C4. In N-glycan1, the 1–6 linkage in the mannose branch [Man(9)-α(1-6)-Man(5)] stays in the gauche-gauche cluster, whereas it moves towards trans-gauche in N-glycan2. For both glycans, mannose branches are more flexible than the complex branches, and adding a complex branch does not influence the dynamics of the mannose branches. We have noticed that the end-to-end distance of the complex branch shortens by ∼ 10 Å in the presence of another complex branch. This suggests that in the presence of an additional complex branch, the other complex branch adopts a close folded structure. All these conformational changes involve the selective formation of inter-residue and water-mediated hydrogen-bond networks.

Inverse Hydrogen-Bonding Change Between the Protonated Retinal Schiff Base and Water Molecules upon Photoisomerization in Heliorhodopsin 48C12.

Heliorhodopsin (HeR) is a new class of the rhodopsin family discovered in 2018 through functional metagenomic analysis (named 48C12). Similar to typical microbial rhodopsins, HeR possesses seven transmembrane (TM) α-helices and an all-trans-retinal covalently bonded to the lysine residue on TM7 via a protonated Schiff base. Remarkably, the HeR membrane topology is inverted compared with that of typical microbial rhodopsins. The X-ray crystal structure of HeR 48C12 was elucidated after the first report on a HeR variant from Thermoplasmatales archaeon SG8-52-1, which revealed the water-mediated hydrogen-bonding network connected to the Schiff base region in the cytoplasmic side. Herein, low-temperature light-induced FTIR spectroscopic analyses of HeR 48C12 and 15N isotopically labeled proteins were used to elucidate the structural changes during retinal photoisomerization. N-D stretching vibrations of the protonated retinal Schiff base (PRSB) at 2286 and 2302 cm-1 in the dark state, and 2239 and 2252 cm-1 in the K intermediate were observed. The frequency changes indicated that the hydrogen bond of PRSB strengthens upon photoisomerization in HeR. Moreover, O-D stretching vibration frequencies of the internal water molecules indicate that the hydrogen-bonding strength decreases concomitantly. Therefore, the PRSB hydrogen bond responds to photoisomerization in an opposite way to the hydrogen-bonding network involving water molecules. No frequency changes of the indole N-H or N-D stretching vibrations of tryptophan residues were observed upon photoisomerization, suggesting that all tryptophan residues in the HeR 48C12 maintained the hydrogen-bonding strengths in the K intermediate. These results provide insights into the molecular mechanism of the energy storage and propagation upon retinal photoisomerization in HeR.

Differential GLP-1R Binding and Activation by Peptide and Non-peptide Agonists

Peptide drugs targeting class B1 G-protein-coupled receptors (GPCRs) can treat multiple diseases; however, there remains substantial interest in the development of orally delivered non-peptide drugs. Here, we reveal unexpected overlap between signaling and regulation of the glucagon-like peptide-1 (GLP-1) receptor by the non-peptide agonist PF 06882961 and GLP-1 that was not observed for another compound, CHU-128. Compounds from these patent series, including PF 06882961, are currently in clinical trials for treatment of type 2 diabetes. High-resolution cryoelectron microscopy (cryo-EM) structures reveal that the binding sites for PF 06882961 and GLP-1 substantially overlap, whereas CHU-128 adopts a unique binding mode with a more open receptor conformation at the extracellular face. Structural differences involving extensive water-mediated hydrogen bond networks could be correlated to functional data to understand how PF 06882961, but not CHU-128, can closely mimic the pharmacological properties of GLP-1. These findings will facilitate rational structure-based discovery of non-peptide agonists targeting class B GPCRs.

Structures of Hsp90α and Hsp90β bound to a purine-scaffold inhibitor reveal an exploitable residue for drug selectivity.

Hsp90α and Hsp90β are implicated in a number of cancers and neurodegenerative disorders but the lack of selective pharmacological probes confounds efforts to identify their individual roles. Here, we analyzed the binding of an Hsp90α-selective PU compound, PU-11-trans, to the two cytosolic paralogs. We determined the co-crystal structures of Hsp90α and Hsp90β bound to PU-11-trans, as well as the structure of the apo Hsp90β NTD. The two inhibitor-bound structures reveal that Ser52, a nonconserved residue in the ATP binding pocket in Hsp90α, provides additional stability to PU-11-trans through a water-mediated hydrogen-bonding network. Mutation of Ser52 to alanine, as found in Hsp90β, alters the dissociation constant of Hsp90α for PU-11-trans to match that of Hsp90β. Our results provide a structural explanation for the binding preference of PU inhibitors for Hsp90α and demonstrate that the single nonconserved residue in the ATP-binding pocket may be exploited for α/β selectivity.

Electrostatics, proton sensor, and networks governing the gating transition in GLIC, a proton-gated pentameric ion channel

The pentameric ligand-gated ion channel (pLGIC) from Gloeobacter violaceus (GLIC) has provided insightful structure-function views on the permeation process and the allosteric regulation of the pLGICs family. However, GLIC is activated by pH instead of a neurotransmitter and a clear picture for the gating transition driven by protons is still lacking. We used an electrostatics-based (finite difference Poisson-Boltzmann/Debye-Hückel) method to predict the acidities of all aspartic and glutamic residues in GLIC, both in its active and closed-channel states. Those residues with a predicted pKa close to the experimental pH50 were individually replaced by alanine and the resulting variant receptors were titrated by ATR/FTIR spectroscopy. E35, located in front of loop F far away from the orthosteric site, appears as the key proton sensor with a measured individual pKa at 5.8. In the GLIC open conformation, E35 is connected through a water-mediated hydrogen-bond network first to the highly conserved electrostatic triad R192-D122-D32 and then to Y197-Y119-K248, both located at the extracellular domain-transmembrane domain interface. The second triad controls a cluster of hydrophobic side chains from the M2-M3 loop that is remodeled during the gating transition. We solved 12 crystal structures of GLIC mutants, 6 of them being trapped in an agonist-bound but nonconductive conformation. Combined with previous data, this reveals two branches of a continuous network originating from E35 that reach, independently, the middle transmembrane region of two adjacent subunits. We conclude that GLIC's gating proceeds by making use of loop F, already known as an allosteric site in other pLGICs, instead of the classic orthosteric site.

A structural mechanism for directing corepressor-selective inverse agonism of PPAR\u03b3

Small chemical modifications can have significant effects on ligand efficacy and receptor activity, but the underlying structural mechanisms can be difficult to predict from static crystal structures alone. Here we show how a simple phenyl-to-pyridyl substitution between two common covalent orthosteric ligands targeting peroxisome proliferator-activated receptor (PPAR) gamma converts a transcriptionally neutral antagonist (GW9662) into a repressive inverse agonist (T0070907) relative to basal cellular activity. X-ray crystallography, molecular dynamics simulations, and mutagenesis coupled to activity assays reveal a water-mediated hydrogen bond network linking the T0070907 pyridyl group to Arg288 that is essential for corepressor-selective inverse agonism. NMR spectroscopy reveals that PPARγ exchanges between two long-lived conformations when bound to T0070907 but not GW9662, including a conformation that prepopulates a corepressor-bound state, priming PPARγ for high affinity corepressor binding. Our findings demonstrate that ligand engagement of Arg288 may provide routes for developing corepressor-selective repressive PPARγ ligands.

Exploiting a water network to achieve enthalpy-driven, bromodomain-selective BET inhibitors

Within the last decade, the Bromodomain and Extra-Terminal domain family (BET) of proteins have emerged as promising drug targets in diverse clinical indications including oncology, auto-immune disease, heart failure, and male contraception. The BET family consists of four isoforms (BRD2, BRD3, BRD4, and BRDT/BRDT6) which are distinguished by the presence of two tandem bromodomains (BD1 and BD2) that independently recognize acetylated-lysine (KAc) residues and appear to have distinct biological roles. BET BD1 and BD2 bromodomains differ at five positions near the substrate binding pocket: the variation in the ZA channel induces different water networks nearby. We designed a set of congeneric 2- and 3-heteroaryl substituted tetrahydroquinolines (THQ) to differentially engage bound waters in the ZA channel with the goal of achieving bromodomain selectivity. SJ830599 (9) showed modest, but consistent, selectivity for BRD2-BD2. Using isothermal titration calorimetry, we showed that the binding of all THQ analogs in our study to either of the two bromodomains was enthalpy driven. Remarkably, the binding of 9 to BRD2-BD2 was marked by negative entropy and was entirely driven by enthalpy, consistent with significant restriction of conformational flexibility and/or engagement with bound waters. Co-crystallography studies confirmed that 9 did indeed stabilize a water-mediated hydrogen bond network. Finally, we report that 9 retained cytotoxicity against several pediatric cancer cell lines with EC50 values comparable to BET inhibitor (BETi) clinical candidates.

X-ray structures of endothelin ETB receptor bound to clinical antagonist bosentan and its analog.

Endothelin receptors (ETRs) have crucial roles in vascular control and are targets for drugs designed to treat circulatory-system diseases and cancer progression. The nonpeptide dual-ETR antagonist bosentan is the first oral drug approved to treat pulmonary arterial hypertension. Here we report crystal structures of human endothelin ETB receptor bound to bosentan and to the ETB-selective analog K-8794, at 3.6-Å and 2.2-Å resolution, respectively. The K-8794-bound structure reveals the detailed water-mediated hydrogen-bonding network at the transmembrane core, which could account for the weak negative allosteric modulation of ETB by Na+ ions. The bosentan-bound structure reveals detailed interactions with ETB, which are probably conserved in the ETA receptor. A comparison of the two structures shows unexpected similarity between antagonist and agonist binding. Despite this similarity, bosentan sterically prevents the inward movement of transmembrane helix 6 (TM6), and thus exerts its antagonistic activity. These structural insights will facilitate the rational design of new ETR-targeting drugs.

A water-mediated allosteric network governs activation of Aurora kinase A.

The catalytic activity of many protein kinases is controlled by conformational changes of a conserved Asp-Phe-Gly (DFG) motif. We used an infrared probe to track the DFG motif of the mitotic kinase Aurora A (AurA) and found that allosteric activation by the spindle-associated protein Tpx2 involves an equilibrium shift toward the active DFG-in state. Förster resonance energy transfer experiments show that the activation loop undergoes a nanometer-scale movement that is tightly coupled to the DFG equilibrium. Tpx2 further activates AurA by stabilizing a water-mediated allosteric network that links the C-helix to the active site through an unusual polar residue in the regulatory spine. The polar spine residue and water network of AurA are essential for phosphorylation-driven activation, but an alternative form of the water network found in related kinases can support Tpx2-driven activation, suggesting that variations in the water-mediated hydrogen bond network mediate regulatory diversification in protein kinases.

Computational and Experimental Characterization of Patient Derived Mutations Reveal an Unusual Mode of Regulatory Spine Assembly and Drug Sensitivity in EGFR Kinase.

The catalytic activation of protein kinases requires precise positioning of key conserved catalytic and regulatory motifs in the kinase core. The Regulatory Spine (RS) is one such structural motif that is dynamically assembled upon kinase activation. The RS is also a mutational hotspot in cancers; however, the mechanisms by which cancer mutations impact RS assembly and kinase activity are not fully understood. In this study, through mutational analysis of patient derived mutations in the RS of EGFR kinase, we identify an activating mutation, M766T, at the RS3 position. RS3 is located in the regulatory αC-helix, and a series of mutations at the RS3 position suggest a strong correlation between the amino acid type present at the RS3 position and ligand (EGF) independent EGFR activation. Small polar amino acids increase ligand independent activity, while large aromatic amino acids decrease kinase activity. M766T relies on the canonical asymmetric dimer for full activation. Molecular modeling and molecular dynamics simulations of WT and mutant EGFR suggest a model in which M766T activates the kinase domain by disrupting conserved autoinhibitory interactions between M766 and hydrophobic residues in the activation segment. In addition, a water mediated hydrogen bond network between T766, the conserved K745-E762 salt bridge, and the backbone amide of the DFG motif is identified as a key determinant of M766T-mediated activation. M766T is resistant to FDA approved EGFR inhibitors such as gefitinib and erlotinib, and computational estimation of ligand binding free energy identifies key residues associated with drug sensitivity. In sum, our studies suggest an unusual mode of RS assembly and oncogenic EGFR activation, and provide new clues for the design of allosteric protein kinase inhibitors.