Comprehensive machine learning boosts structure-based virtual screening for PARP1 inhibitors.

Poly ADP-ribose polymerase 1 (PARP1) is an attractive therapeutic target for cancer treatment. Machine-learning scoring functions constitute a promising approach to discovering novel PARP1 inhibitors. Cutting-edge PARP1-specific machine-learning scoring functions were investigated using semi-synthetic training data from docking activity-labelled molecules: known PARP1 inhibitors, hard-to-discriminate decoys property-matched to them with generative graph neural networks and confirmed inactives. We further made test sets harder by including only molecules dissimilar to those in the training set. Comprehensive analysis of these datasets using five supervised learning algorithms, and protein-ligand fingerprints extracted from docking poses and ligand only features revealed one highly predictive scoring function. This is the PARP1-specific support vector machine-based regressor, when employing PLEC fingerprints, which achieved a high Normalized Enrichment Factor at the top 1% on the hardest test set (NEF1% = 0.588, median of 10 repetitions), and was more predictive than any other investigated scoring function, especially the classical scoring function employed as baseline.

Coordinating activation of endo-lysosomal two-pore channels and TRP mucolipins.

Two-pore channels and TRP mucolipins are ubiquitous endo-lysosomal cation channels of pathophysiological relevance. Both are Ca2+-permeable and regulated by phosphoinositides, principally PI(3,5)P2. Accumulating evidence has uncovered synergistic channel activation by PI(3,5)P2 and endogenous metabolites such as the Ca2+ mobilizing messenger NAADP, synthetic agonists including approved drugs and physical cues such as voltage and osmotic pressure. Here, we provide an overview of this coordination.

Unraveling the Behavior of Intrinsically Disordered Protein c-Myc: A Study Utilizing Gaussian-Accelerated Molecular Dynamics.

The protein c-Myc is a transcription factor that remains largely intrinsically disordered and is known to be involved in various biological processes and is overexpressed in various cancers, making it an attractive drug target. However, intrinsically disordered proteins such as c-Myc do not show funnel-like basins in their free-energy landscapes; this makes their druggability a challenge. For the first time, we propose a heterodimer model of c-Myc/Max in full length in this work. We used Gaussian-accelerated molecular dynamics (GaMD) simulations to explore the behavior of c-Myc and its various regions, including the transactivation domain (TAD) and the basic helix-loop-helix-leucine-zipper (bHLH-Zipper) motif in three different conformational states: (a) monomeric c-Myc, (b) c-Myc when bound to its partner protein, Max, and (c) when Max was removed after binding. We analyzed the GaMD trajectories using root-mean-square deviation (RMSD), radius of gyration, root-mean-square fluctuation, and free-energy landscape (FEL) calculations to elaborate the behaviors of these regions. The results showed that the monomeric c-Myc structure showed a higher RMSD fluctuation as compared with the c-Myc/Max heterodimer in the bHLH-Zipper motif. This indicated that the bHLH-Zipper motif of c-Myc is more stable when it is bound to Max. The TAD region in both monomeric and Max-bound states showed similar plasticity in terms of RMSD. We also conducted residue decomposition calculations and showed that the c-Myc and Max interaction could be driven mainly by electrostatic interactions and the residues Arg299, Ile403, and Leu420 seemed to play important roles in the interaction. Our work provides insights into the behavior of c-Myc and its regions that could support the development of drugs that target c-Myc and other intrinsically disordered proteins.

Combination of QSAR Modeling and Hybrid-Based Consensus Scoring to Identify Dual-Targeting Inhibitors of PLK1 and p38γ.

Polo-like kinase 1 (PLK1) and p38γ mitogen-activated protein kinase (p38γ) play important roles in cancer pathogenesis by controlling cell cycle progression and are therefore attractive cancer targets. The design of multitarget inhibitors may offer synergistic inhibition of distinct targets and reduce the risk of drug-drug interactions to improve the balance between therapeutic efficacy and safety. We combined deep-learning-based quantitative structure-activity relationship (QSAR) modeling and hybrid-based consensus scoring to screen for inhibitors with potential activity against the targeted proteins. Using this combination strategy, we identified a potent PLK1 inhibitor (compound 4) that inhibited PLK1 activity and liver cancer cell growth in the nanomolar range. Next, we deployed both our QSAR models for PLK1 and p38γ on the Enamine compound library to identify dual-targeting inhibitors against PLK1 and p38γ. Likewise, the identified hits were subsequently subjected to hybrid-based consensus scoring. Using this method, we identified a promising compound (compound 14) that could inhibit both PLK1 and p38γ activities. At nanomolar concentrations, compound 14 inhibited the growth of human hepatocellular carcinoma and hepatoblastoma cells in vitro. This study demonstrates the combined screening strategy to identify novel potential inhibitors for existing targets.

Convergent activation of Ca2+ permeability in two-pore channel 2 through distinct molecular routes.

TPC2 is a pathophysiologically relevant lysosomal ion channel that is activated directly by the phosphoinositide PI(3,5)P2 and indirectly by the calcium ion (Ca2+)-mobilizing molecule NAADP through accessory proteins that associate with the channel. TPC2 toggles between PI(3,5)P2-induced, sodium ion (Na+)-selective and NAADP-induced, Ca2+-permeable states in response to these cues. To address the molecular basis of polymodal gating and ion-selectivity switching, we investigated the mechanism by which NAADP and its synthetic functional agonist, TPC2-A1-N, induced Ca2+ release through TPC2 in human cells. Whereas NAADP required the NAADP-binding proteins JPT2 and LSM12 to evoke endogenous calcium ion signals, TPC2-A1-N did not. Residues in TPC2 that bind to PI(3,5)P2 were required for channel activation by NAADP but not for activation by TPC2-A1-N. The cryptic voltage-sensing region of TPC2 was required for the actions of TPC2-A1-N and PI(3,5)P2 but not for those of NAADP. These data mechanistically distinguish natural and synthetic agonist action at TPC2 despite convergent effects on Ca2+ permeability and delineate a route for pharmacologically correcting impaired NAADP-evoked Ca2+ signals.

Human soluble CD39 displays substrate inhibition in a substrate-specific manner.

CD39 (ectonucleoside triphosphate diphosphohydrolase-1; ENTPD1) metabolizes extracellular ATP and ADP to AMP. AMP is subsequently metabolized by CD79 to adenosine. CD39 activity is therefore a key regulator of purinergic signalling in cancer, thrombosis, and autoimmune diseases. In this study we demonstrate that soluble, recombinant CD39 shows substrate inhibition with ADP or ATP as the substrate. Although CD39 activity initially increased with increasing substrate concentration, at high concentrations of ATP or ADP, CD39 activity was markedly reduced. Although the reaction product, AMP, inhibits CD39 activity, insufficient AMP was generated under our conditions to account for the substrate inhibition seen. In contrast, inhibition was not seen with UDP or UTP as substrates. 2-methylthio-ADP also showed no substrate inhibition, indicating the nucleotide base is an important determinant of substrate inhibition. Molecular dynamics simulations revealed that ADP can undergo conformational rearrangements within the CD39 active site that were not seen with UDP or 2-methylthio-ADP. Appreciating the existence of substrate inhibition of CD39 will help the interpretation of studies of CD39 activity, including investigations into drugs that modulate CD39 activity.

The ß3-subunit modulates the effect of venom peptides ProTx-II and OD1 on NaV 1.7 gating.

The voltage-gated sodium channel NaV 1.7 is involved in various pain phenotypes and is physiologically regulated by the NaV -ß3-subunit. Venom toxins ProTx-II and OD1 modulate NaV 1.7 channel function and may be useful as therapeutic agents and/or research tools. Here, we use patch-clamp recordings to investigate how the ß3-subunit can influence and modulate the toxin-mediated effects on NaV 1.7 function, and we propose a putative binding mode of OD1 on NaV 1.7 to rationalise its activating effects. The inhibitor ProTx-II slowed the rate of NaV 1.7 activation, whilst the activator OD1 reduced the rate of fast inactivation and accelerated recovery from inactivation. The ß3-subunit partially abrogated these effects. OD1 induced a hyperpolarising shift in the V1/2 of steady-state activation, which was not observed in the presence of ß3. Consequently, OD1-treated NaV 1.7 exhibited an enhanced window current compared with OD1-treated NaV 1.7-ß3 complex. We identify candidate OD1 residues that are likely to prevent the upward movement of the DIV S4 helix and thus impede fast inactivation. The binding sites for each of the toxins and the predicted location of the ß3-subunit on the NaV 1.7 channel are distinct. Therefore, we infer that the ß3-subunit influences the interaction of toxins with NaV 1.7 via indirect allosteric mechanisms. The enhanced window current shown by OD1-treated NaV 1.7 compared with OD1-treated NaV 1.7-ß3 is discussed in the context of differing cellular expressions of NaV 1.7 and the ß3-subunit in dorsal root ganglion (DRG) neurons. We propose that ß3, as the native binding partner for NaV 1.7 in DRG neurons, should be included during screening of molecules against NaV 1.7 in relevant analgesic discovery campaigns.

Identification of Potential p38γ Inhibitors via In Silico Screening, In Vitro Bioassay and Molecular Dynamics Simulation Studies.

Protein kinase p38γ is an attractive target against cancer because it plays a pivotal role in cancer cell proliferation by phosphorylating the retinoblastoma tumour suppressor protein. Therefore, inhibition of p38γ with active small molecules represents an attractive alternative for developing anti-cancer drugs. In this work, we present a rigorous and systematic virtual screening framework to identify potential p38γ inhibitors against cancer. We combined the use of machine learning-based quantitative structure activity relationship modelling with conventional computer-aided drug discovery techniques, namely molecular docking and ligand-based methods, to identify potential p38γ inhibitors. The hit compounds were filtered using negative design techniques and then assessed for their binding stability with p38γ through molecular dynamics simulations. To this end, we identified a promising compound that inhibits p38γ activity at nanomolar concentrations and hepatocellular carcinoma cell growth in vitro in the low micromolar range. This hit compound could serve as a potential scaffold for further development of a potent p38γ inhibitor against cancer.

In-silico discovery of inhibitors against human papillomavirus E1 protein.

High-risk (HR) Human papillomavirus (e.g. HPV16 and HPV18) causes approximately two-thirds of all cervical cancers in women. Although the first and second-generation vaccines confer some protection against individuals, there are no approved drugs to treat HR-HPV infections to-date. The HPV E1 protein is an attractive drug target because the protein is highly conserved across all HPV types and is crucial for the regulation of viral DNA replication. Hence, we used the Random Forest algorithm to construct a Quantitative-Structure Activity Relationship (QSAR) model to predict the potential inhibitors against the HPV E1 protein. Our QSAR classification model achieved an accuracy of 87.5%, area under the receiver operating characteristic curve of 1.00, and F-measure of 0.87 when evaluated using an external test set. We conducted a drug repurposing campaign by deploying the model to screen the Drugbank database. The top three compounds, namely Cinalukast, Lobeglitazone, and Efatutazone were analyzed for their cell membrane permeability, toxicity, and carcinogenicity. Finally, these three compounds were subjected to molecular docking and 200 ns-long Molecular Dynamics (MD) simulations. The predicted binding free energies for the candidates were calculated using the MM-GBSA method. The binding free energies for Cinalukast, Lobeglitazone, and Efatutazone were -37.84 kcal/mol, -25.30 kcal/mol, and -29.89 kcal/mol respectively. Therefore, we propose their chemical scaffolds for future rational design of E1 inhibitors.Communicated by Ramaswamy H. Sarma.

Repurposing FDA-approved drugs as HIV-1 integrase inhibitors: an in silico investigation.

The Human Immunodeficiency Virus (HIV) infection is a global pandemic that has claimed 33 million lives to-date. One of the most efficacious treatments for naïve or pretreated HIV patients is the HIV integrase strand transfer inhibitors (INSTIs). However, given that HIV treatment is life-long, the emergence of HIV strains resistant to INSTIs is an imminent challenge. In this work, we showed two best regression QSAR models that were constructed using a boosted Random Forest algorithm (r2 = 0.998, q210CV = 0.721, q2external_test = 0.754) and a boosted K* algorithm (r2 = 0.987, q210CV = 0.721, q2external_test = 0.758) to predict the pIC50 values of INSTIs. Subsequently, the regression QSAR models were deployed against the Drugbank database for drug repositioning. The top-ranked compounds were further evaluated for their target engagement activity using molecular docking studies and accelerated Molecular Dynamics simulation. Lastly, their potential as INSTIs were also evaluated from our literature search. Our study offers the first example of a large-scale regression QSAR modelling effort for discovering highly active INSTIs to combat HIV infection.Communicated by Ramaswamy H. Sarma.

Quantal leaps in understanding Ca2+ signaling: A "Taylored" approach.

On 2 September 2022, about 85 scientists gathered in person at Queens' College in Cambridge, UK, for a scientific meeting to celebrate the career of Colin W. Taylor of Cambridge University upon his retirement. The meeting was organized by the authors, who are all former graduate students in the Taylor laboratory, which has been at the forefront of Ca2+ signaling for more than 30 years.

Electrophysiology of Endolysosomal Two-Pore Channels: A Current Account.

Two-pore channels TPC1 and TPC2 are ubiquitously expressed pathophysiologically relevant proteins that reside on endolysosomal vesicles. Here, we review the electrophysiology of these channels. Direct macroscopic recordings of recombinant TPCs expressed in enlarged lysosomes in mammalian cells or vacuoles in plants and yeast demonstrate gating by the Ca2+-mobilizing messenger NAADP and/or the lipid PI(3,5)P2. TPC currents are regulated by H+, Ca2+, and Mg2+ (luminal and/or cytosolic), as well as protein kinases, and they are impacted by single-nucleotide polymorphisms linked to pigmentation. Bisbenzylisoquinoline alkaloids, flavonoids, and several approved drugs demonstrably block channel activity. Endogenous TPC currents have been recorded from a number of primary cell types and cell lines. Many of the properties of endolysosomal TPCs are recapitulated upon rerouting channels to the cell surface, allowing more facile recording through conventional electrophysiological means. Single-channel analyses have provided high-resolution insight into both monovalent and divalent permeability. The discovery of small-molecule activators of TPC2 that toggle the ion selectivity from a Ca2+-permeable (NAADP-like) state to a Na+-selective (PI(3,5)P2-like) state explains discrepancies in the literature relating to the permeability of TPCs. Identification of binding proteins that confer NAADP-sensitive currents confirm that indirect, remote gating likely underpins the inconsistent observations of channel activation by NAADP.

Segregated cation flux by TPC2 biases Ca2+ signaling through lysosomes.

Two-pore channels are endo-lysosomal cation channels with malleable selectivity filters that drive endocytic ion flux and membrane traffic. Here we show that TPC2 can differentially regulate its cation permeability when co-activated by its endogenous ligands, NAADP and PI(3,5)P2. Whereas NAADP rendered the channel Ca2+-permeable and PI(3,5)P2 rendered the channel Na+-selective, a combination of the two increased Ca2+ but not Na+ flux. Mechanistically, this was due to an increase in Ca2+ permeability independent of changes in ion selectivity. Functionally, we show that cell permeable NAADP and PI(3,5)P2 mimetics synergistically activate native TPC2 channels in live cells, globalizing cytosolic Ca2+ signals and regulating lysosomal pH and motility. Our data reveal that flux of different ions through the same pore can be independently controlled and identify TPC2 as a likely coincidence detector that optimizes lysosomal Ca2+ signaling.

The inositol 5-phosphatase INPP5B regulates B cell receptor clustering and signaling.

Upon antigen binding, the B cell receptor (BCR) undergoes clustering to form a signalosome that propagates downstream signaling required for normal B cell development and physiology. BCR clustering is dependent on remodeling of the cortical actin network, but the mechanisms that regulate actin remodeling in this context remain poorly defined. In this study, we identify the inositol 5-phosphatase INPP5B as a key regulator of actin remodeling, BCR clustering, and downstream signaling in antigen-stimulated B cells. INPP5B acts via dephosphorylation of the inositol lipid PI(4,5)P2 that in turn is necessary for actin disassembly, BCR mobilization, and cell spreading on immobilized surface antigen. These effects can be explained by increased actin severing by cofilin and loss of actin linking to the plasma membrane by ezrin, both of which are sensitive to INPP5B-dependent PI(4,5)P2 hydrolysis. INPP5B is therefore a new player in BCR signaling and may represent an attractive target for treatment of B cell malignancies caused by aberrant BCR signaling.

Structure-Based Discovery of Lipoteichoic Acid Synthase Inhibitors.

Lipoteichoic acid synthase (LtaS) is a key enzyme for the cell wall biosynthesis of Gram-positive bacteria. Gram-positive bacteria that lack lipoteichoic acid (LTA) exhibit impaired cell division and growth defects. Thus, LtaS appears to be an attractive antimicrobial target. The pharmacology around LtaS remains largely unexplored with only two small-molecule LtaS inhibitors reported, namely "compound 1771" and the Congo red dye. Structure-based drug discovery efforts against LtaS remain unattempted due to the lack of an inhibitor-bound structure of LtaS. To address this, we combined the use of a molecular docking technique with molecular dynamics (MD) simulations to model a plausible binding mode of compound 1771 to the extracellular catalytic domain of LtaS (eLtaS). The model was validated using alanine mutagenesis studies combined with isothermal titration calorimetry. Additionally, lead optimization driven by our computational model resulted in an improved version of compound 1771, namely, compound 4 which showed greater affinity for binding to eLtaS than compound 1771 in biophysical assays. Compound 4 reduced LTA production in S. aureus dose-dependently, induced aberrant morphology as seen for LTA-deficient bacteria, and significantly reduced bacteria titers in the lung of mice infected with S. aureus. Analysis of our MD simulation trajectories revealed the possible formation of a transient cryptic pocket in eLtaS. Virtual screening (VS) against the cryptic pocket led to the identification of a new class of inhibitors that could potentiate ß-lactams against methicillin-resistant S. aureus. Our overall workflow and data should encourage further drug design campaign against LtaS. Finally, our work reinforces the importance of considering protein conformational flexibility to a successful VS endeavor.

Rendezvous with PI(3,5)P2 - A rapalog gets caught opening TRPML1.

TRPML1 is an endolysosomally-expressed cation channel, activated physiologically by PI(3,5)P2 and by several synthetic agonists including rapamycin. New high resolution cryo-EM- structures of TRPML1 bound to both PI(3,5)P2 and temsirolimus - a rapamycin analog provides molecular insight into how the channel integrates two agonists that bind to distal sites but act cooperatively.

Pharmacological blockade of angiotensin II receptor restores diabetes-associated reduction of store operated Ca2+ entry in adult cardiomyocytes.

The store-operated Ca2+ entry (SOCE) represents an important route for generating cellular Ca2+ signals that are implicated in physiological and various pathological scenarios that include diabetic cardiomyopathy (DM-CMP) which is well known to have Ca2+ dysregulation among other salient features. In this study, we investigated the role of SOCE in Ca2+ handling of cardiomyocytes obtained from adult male Wistar rats that were made diabetic by intraperitoneal administration of streptozotocin (STZ 50 mg/kg). We also included another group of rats with diabetes induced by STZ administration but received an angiotensin II receptor blocker - losartan. In whole cell recordings with isolated cardiomyocytes, the SOCE-representative whole-cell current ICRAC was found to be significantly reduced for the diabetic group compared to the control group and chronic losartan treatment could restore ICRAC to a level comparable to the control group. However, in contrast to the observed reduction in ICRAC, Orai1 and Orai3 proteins were found to be significantly upregulated in diabetic condition whereas no significant change in the expression levels of Stim1, Stim2 and Orai2 was observed. Also, losartan treatment did not affect the expression pattern of these key proteins for SOCE in diabetic group. The observed imbalance between the functional read out of SOCE (peak ICRAC size) and expression levels of the underlying proteins was puzzling but could be, among other possibilities, due to impairment of interaction between Stim and Orai proteins. We argue that the observed changes in SOCE with diabetes could be a contributing factor for the Ca2+ dyshomeostasis associated with diabetic cardiomyopathies and blockade of angiotensin II receptor may potentially restore normal SOCE in diabetic cardiomyocytes.

Q94 is not a selective modulator of proteinase-activated receptor 1 (PAR1) in platelets.

Thrombin is a potent platelet activator, acting through proteinase-activated receptors -1 and -4 (PAR1 and PAR4). Of these, PAR-1 is activated more rapidly and by lower thrombin concentrations. Consequently, PAR-1 has been extensively investigated as a target for anti-platelet drugs to prevent myocardial infarction. Q94 has been reported to act as an allosteric modulator of PAR1, potently and selectively inhibiting PAR1-Gαq coupling in multiple cell lines, but its effects on human platelet activation have not been previously studied. Platelet Ca2+ signaling, integrin αIIbß3 activation and α-granule secretion were monitored following stimulation by a PAR1-activating peptide (PAR1-AP). Although Q94 inhibited these responses, its potency was low compared to other PAR1 antagonists. In addition, αIIbß3 activation and α-granule secretion in response to other platelet activators were also inhibited with similar potency. Finally, in endothelial cells, Q94 did not inhibit PAR1-dependent Ca2+ signaling. Our data suggest that Q94 may have PAR1-independent off-target effects in platelets, precluding its use as a selective PAR1 allosteric modulator.

Activation of endo-lysosomal two-pore channels by NAADP and PI(3,5)P2. Five things to know.

Two-pore channels are ancient members of the voltage-gated ion channel superfamily that are expressed predominantly on acidic organelles such as endosomes and lysosomes. Here we review recent advances in understanding how TPCs are activated by their ligands and identify five salient features: (1) TPCs are Ca2+-permeable non-selective cation channels gated by NAADP. (2) NAADP activation is indirect through associated NAADP receptors. (3) TPCs are also Na+-selective channels gated by PI(3,5)P2. (4) PI(3,5)P2 activation is direct through a structurally-resolved binding site. (5) TPCs switch their ion selectivity in an agonist-dependent manner.

Structure, Function and Regulation of a Second Pyruvate Kinase Isozyme in Pseudomonas aeruginosa.

Pseudomonas aeruginosa (PA) depends on the Entner-Doudoroff pathway (EDP) for glycolysis. The main enzymatic regulator in the lower half of the EDP is pyruvate kinase. PA contains genes that encode two isoforms of pyruvate kinase, denoted PykAPA and PykFPA. In other well-characterized organisms containing two pyruvate kinase isoforms (such as Escherichia coli) each isozyme is differentially regulated. The structure, function and regulation of PykAPA has been previously characterized in detail, so in this work, we set out to assess the biochemical and structural properties of the PykFPA isozyme. We show that pykF PA expression is induced in the presence of the diureide, allantoin. In spite of their relatively low amino acid sequence identity, PykAPA and PykFPA display broadly comparable kinetic parameters, and are allosterically regulated by a very similar set of metabolites. However, the x-ray crystal structure of PykFPA revealed significant differences compared with PykAPA. Notably, although the main allosteric regulator binding-site of PykFPA was empty, the "ring loop" covering the site adopted a partially closed conformation. Site-directed mutation of the proline residues flanking the ring loop yielded apparent "locked on" and "locked off" allosteric activation phenotypes, depending on the residue mutated. Analysis of PykFPA inter-protomer interactions supports a model in which the conformational transition(s) accompanying allosteric activation involve re-orientation of the A and B domains of the enzyme and subsequent closure of the active site.