WatFinder: A ProDy tool for protein-water interactions.

SUMMARY: We introduce WatFinder, a tool designed to identify and visualize protein-water interactions (water bridges, water-mediated associations, or water channels, fluxes, and clusters) relevant to protein stability, dynamics, and function. WatFinder is integrated into ProDy, a Python API broadly used for structure-based prediction of protein dynamics. WatFinder provides a suite of functions for generating raw data as well as outputs from statistical analyses. The ProDy framework facilitates comprehensive automation and efficient analysis of the ensembles of structures resolved for a given protein or the time-evolved conformations from simulations in explicit water, as illustrated in five case studies presented in the Supplementary Material. AVAILABILITY AND IMPLEMENTATION: ProDy is open-source and freely available under MIT License from https://github.com/ProDy/ProDy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online, and tutorials at http://www.bahargroup.org/prody.

Characterization of hyperactive mutations in the renal potassium channel ROMK uncovers unique effects on channel biogenesis and ion conductance.

Hypertension affects one billion people worldwide and is the most common risk factor for cardiovascular disease, yet a comprehensive picture of its underlying genetic factors is incomplete. Amongst regulators of blood pressure is the renal outer medullary potassium (ROMK) channel. While select ROMK mutants are prone to premature degradation and lead to disease, heterozygous carriers of some of these same alleles are protected from hypertension. Therefore, we hypothesized that gain-of-function (GoF) ROMK variants which increase potassium flux may predispose people to hypertension. To begin to test this hypothesis, we employed genetic screens and a candidate-based approach to identify six GoF variants in yeast. Subsequent functional assays in higher cells revealed two variant classes. The first group exhibited greater stability in the endoplasmic reticulum, enhanced channel assembly, and/or increased protein at the cell surface. The second group of variants resided in the PIP2-binding pocket, and computational modeling coupled with patch-clamp studies demonstrated lower free energy for channel opening and slowed current rundown, consistent with an acquired PIP2-activated state. Together, these findings advance our understanding of ROMK structure-function, suggest the existence of hyperactive ROMK alleles in humans, and establish a system to facilitate the development of ROMK-targeted antihypertensives.

Allosteric modulation of serotonin and dopamine transporters: New insights from computations and experiments.

Human monoamine transporters (MATs) are critical to regulating monoaminergic neurotransmission by translocating their substrates from the synaptic space back into the presynaptic neurons. As such, their primary substrate binding site S1 has been targeted by a wide range of compounds for treating neuropsychiatric and neurodegenerative disorders including depression, ADHD, neuropathic pain, and anxiety disorders. We present here a comparative study of the structural dynamics and ligand-binding properties of two MATs, dopamine transporter (DAT) and serotonin transporter (SERT), with focus on the allosteric modulation of their transport function by drugs or substrates that consistently bind a secondary site S2, proposed to serve as an allosteric site. Our systematic analysis of the conformational space and dynamics of a dataset of 50 structures resolved for DAT and SERT in the presence of one or more ligands/drugs reveals the specific residues playing a consistent role in coordinating the small molecules bound to subsites S2-I and S2-II within S2, such as R476 and Y481 in dDAT and E494, P561, and F556 in hSERT. Further analysis reveals how DAT and SERT differ in their two principal modes of structural changes, PC1 and PC2. Notably, PC1 underlies the transition between outward- and inward-facing states of the transporters as well as their gating; whereas PC2 supports the rearrangements of TM helices near the S2 site. Finally, the examination of cross-correlations between structural elements lining the respective sites S1 and S2 point to the crucial role of coupled motions between TM6a and TM10. In particular, we note the involvement of hSERT residues F335 and G338, and E493-E494-T497 belonging to these two respective helices, in establishing the allosteric communication between S1 and S2. These results help understand the molecular basis of the action of drugs that bind to the S2 site of DAT or SERT. They also provide a basis for designing allosteric modulators that may provide better control of specific interactions and cellular pathways, rather than indiscriminately inhibiting the transporter by targeting its orthosteric site.

Influence of point mutations on PR65 conformational adaptability: Insights from molecular simulations and nanoaperture optical tweezers.

PR65 is the HEAT repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem repeat protein. Its conformational mechanics plays a crucial role in PP2A function by opening/closing substrate binding/catalysis interface. Using in silico saturation mutagenesis, we identified PR65 "hinge" residues whose substitutions could alter its conformational adaptability and thereby PP2A function, and selected six mutations that were verified to be expressed and soluble. Molecular simulations and nanoaperture optical tweezers revealed consistent results on the specific effects of the mutations on the structure and dynamics of PR65. Two mutants observed in simulations to stabilize extended/open conformations exhibited higher corner frequencies and lower translational scattering in experiments, indicating a shift toward extended conformations, whereas another displayed the opposite features, confirmed by both simulations and experiments. The study highlights the power of single-molecule nanoaperture-based tweezers integrated with in silico approaches for exploring the effect of mutations on protein structure and dynamics.

Genome structural dynamics: insights from Gaussian network analysis of Hi-C data.

Characterization of the spatiotemporal properties of the chromatin is essential to gaining insights into the physical bases of gene co-expression, transcriptional regulation and epigenetic modifications. The Gaussian network model (GNM) has proven in recent work to serve as a useful tool for modeling chromatin structural dynamics, using as input high-throughput chromosome conformation capture data. We focus here on the exploration of the collective dynamics of chromosomal structures at hierarchical levels of resolution, from single gene loci to topologically associating domains or entire chromosomes. The GNM permits us to identify long-range interactions between gene loci, shedding light on the role of cross-correlations between distal regions of the chromosomes in regulating gene expression. Notably, GNM analysis performed across diverse cell lines highlights the conservation of the global/cooperative movements of the chromatin across different types of cells. Variations driven by localized couplings between genomic loci, on the other hand, underlie cell differentiation, underscoring the significance of the four-dimensional properties of the genome in defining cellular identity. Finally, we demonstrate the close relation between the cell type-dependent mobility profiles of gene loci and their gene expression patterns, providing a clear demonstration of the role of chromosomal 4D features in defining cell-specific differential expression of genes.

Synthesis and Structure-Activity Relationships for Glutamate Transporter Allosteric Modulators.

Excitatory amino acid transporters (EAATs) are essential CNS proteins that regulate glutamate levels. Excess glutamate release and alteration in EAAT expression are associated with several CNS disorders. Previously, we identified positive allosteric modulators (PAM) of EAAT2, the main CNS transporter, and have demonstrated their neuroprotective properties in vitro. Herein, we report on the structure-activity relationships (SAR) for the analogs identified from virtual screening and from our medicinal chemistry campaign. This work identified several selective EAAT2 positive allosteric modulators (PAMs) such as compounds 4 (DA-023) and 40 (NA-014) from a library of analogs inspired by GT949, an early generation compound. This series also provides nonselective EAAT PAMs, EAAT inhibitors, and inactive compounds that may be useful for elucidating the mechanism of EAAT allosteric modulation.

Structural mechanisms for VMAT2 inhibition by tetrabenazine.

The vesicular monoamine transporter 2 (VMAT2) is a proton-dependent antiporter responsible for loading monoamine neurotransmitters into synaptic vesicles. Dysregulation of VMAT2 can lead to several neuropsychiatric disorders including Parkinson's disease and schizophrenia. Furthermore, drugs such as amphetamine and MDMA are known to act on VMAT2, exemplifying its role in the mechanisms of actions for drugs of abuse. Despite VMAT2's importance, there remains a critical lack of mechanistic understanding, largely driven by a lack of structural information. Here, we report a 3.1 Å resolution cryo-electron microscopy (cryo-EM) structure of VMAT2 complexed with tetrabenazine (TBZ), a non-competitive inhibitor used in the treatment of Huntington's chorea. We find TBZ interacts with residues in a central binding site, locking VMAT2 in an occluded conformation and providing a mechanistic basis for non-competitive inhibition. We further identify residues critical for cytosolic and lumenal gating, including a cluster of hydrophobic residues which are involved in a lumenal gating strategy. Our structure also highlights three distinct polar networks that may determine VMAT2 conformational dynamics and play a role in proton transduction. The structure elucidates mechanisms of VMAT2 inhibition and transport, providing insights into VMAT2 architecture, function, and the design of small-molecule therapeutics.

Pharmacological suppression of the OTUD4/CD73 proteolytic axis revives antitumor immunity against immune-suppressive breast cancers.

Despite widespread utilization of immunotherapy, treating immune-cold tumors remains a challenge. Multiomic analyses and experimental validation identified the OTUD4/CD73 proteolytic axis as a promising target in treating immune-suppressive triple negative breast cancer (TNBC). Mechanistically, deubiquitylation of CD73 by OTUD4 counteracted its ubiquitylation by TRIM21, resulting in CD73 stabilization inhibiting tumor immune responses. We further demonstrated the importance of TGF-ß signaling for orchestrating the OTUD4/CD73 proteolytic axis within tumor cells. Spatial transcriptomics profiling discovered spatially resolved features of interacting malignant and immune cells pertaining to expression levels of OTUD4 and CD73. In addition, ST80, a newly developed inhibitor, specifically disrupted proteolytic interaction between CD73 and OTUD4, leading to reinvigoration of cytotoxic CD8+ T cell activities. In preclinical models of TNBC, ST80 treatment sensitized refractory tumors to anti-PD-L1 therapy. Collectively, our findings uncover what we believe to be a novel strategy for targeting the immunosuppressive OTUD4/CD73 proteolytic axis in treating immune-suppressive breast cancers with the inhibitor ST80.

Genetic regulation and targeted reversal of lysosomal dysfunction and inflammatory sterol metabolism in pulmonary arterial hypertension.

Vascular inflammation critically regulates endothelial cell (EC) pathophenotypes, particularly in pulmonary arterial hypertension (PAH). Dysregulation of lysosomal activity and cholesterol metabolism have known inflammatory roles in disease, but their relevance to PAH is unclear. In human pulmonary arterial ECs and in PAH, we found that inflammatory cytokine induction of the nuclear receptor coactivator 7 (NCOA7) both preserved lysosomal acidification and served as a homeostatic brake to constrain EC immunoactivation. Conversely, NCOA7 deficiency promoted lysosomal dysfunction and proinflammatory oxysterol/bile acid generation that, in turn, contributed to EC pathophenotypes. In vivo, mice deficient for Ncoa7 or exposed to the inflammatory bile acid 7α-hydroxy-3-oxo-4-cholestenoic acid (7HOCA) displayed worsened PAH. Emphasizing this mechanism in human PAH, an unbiased, metabolome-wide association study (N=2,756) identified a plasma signature of the same NCOA7-dependent oxysterols/bile acids associated with PAH mortality (P<1.1x10-6). Supporting a genetic predisposition to NCOA7 deficiency, in genome-edited, stem cell-derived ECs, the common variant intronic SNP rs11154337 in NCOA7 regulated NCOA7 expression, lysosomal activity, oxysterol/bile acid production, and EC immunoactivation. Correspondingly, SNP rs11154337 was associated with PAH severity via six-minute walk distance and mortality in discovery (N=93, P=0.0250; HR=0.44, 95% CI [0.21-0.90]) and validation (N=630, P=2x10-4; HR=0.49, 95% CI [0.34-0.71]) cohorts. Finally, utilizing computational modeling of small molecule binding to NCOA7, we predicted and synthesized a novel activator of NCOA7 that prevented EC immunoactivation and reversed indices of rodent PAH. In summary, we have established a genetic and metabolic paradigm and a novel therapeutic agent that links lysosomal biology as well as oxysterol and bile acid processes to EC inflammation and PAH pathobiology. This paradigm carries broad implications for diagnostic and therapeutic development in PAH and in other conditions dependent upon acquired and innate immune regulation of vascular disease.

Strikingly High Activity of 15-Lipoxygenase Towards Di-Polyunsaturated Arachidonoyl/Adrenoyl-Phosphatidylethanolamines Generates Peroxidation Signals of Ferroptotic Cell Death.

The vast majority of membrane phospholipids (PLs) include two asymmetrically positioned fatty acyls: oxidizable polyunsaturated fatty acids (PUFA) attached predominantly at the sn2 position, and non-oxidizable saturated/monounsaturated acids (SFA/MUFA) localized at the sn1 position. The peroxidation of PUFA-PLs, particularly sn2-arachidonoyl(AA)- and sn2-adrenoyl(AdA)-containing phosphatidylethanolamines (PE), has been associated with the execution of ferroptosis, a program of regulated cell death. There is a minor subpopulation (≈1-2 mol %) of doubly PUFA-acylated phospholipids (di-PUFA-PLs) whose role in ferroptosis remains enigmatic. Here we report that 15-lipoxygenase (15LOX) exhibits unexpectedly high pro-ferroptotic peroxidation activity towards di-PUFA-PEs. We revealed that peroxidation of several molecular species of di-PUFA-PEs occurred early in ferroptosis. Ferrostatin-1, a typical ferroptosis inhibitor, effectively prevented peroxidation of di-PUFA-PEs. Furthermore, co-incubation of cells with di-AA-PE and 15LOX produced PUFA-PE peroxidation and induced ferroptotic death. The decreased contents of di-PUFA-PEs in ACSL4 KO A375 cells was associated with lower levels of di-PUFA-PE peroxidation and enhanced resistance to ferroptosis. Thus, di-PUFA-PE species are newly identified phospholipid peroxidation substrates and regulators of ferroptosis, representing a promising therapeutic target for many diseases related to ferroptotic death.

Allo-targeting of the kinase domain: Insights from in silico studies and comparison with experiments.

The eukaryotic protein kinase domain has been a broadly explored target for drug discovery, despite limitations imposed by its high sequence conservation as a shared modular domain and the development of resistance to drugs. One way of addressing those limitations has been to target its potential allosteric sites, shortly called allo-targeting, in conjunction with, or separately from, its conserved catalytic/orthosteric site that has been widely exploited. Allosteric regulation has gained importance as an alternative to overcome the drawbacks associated with the indiscriminate effect of targeting the active site, and it turned out to be particularly useful for these highly promiscuous and broadly shared kinase domains. Yet, allo-targeting often faces challenges as the allosteric sites are not as clearly defined as its orthosteric sites, and the effect on the protein function may not be unambiguously assessed. A robust understanding of the consequence of site-specific allo-targeting on the conformational dynamics of the target protein is essential to design effective allo-targeting strategies. Recent years have seen important advances in in silico identification of druggable sites and distinguishing among them those sites expected to allosterically mediate conformational switches essential to signal transmission. The present opinion underscores the utility of such computational approaches applied to the kinase domain, with the help of comparison between computational predictions and experimental observations.

Structural mechanisms for VMAT2 inhibition by tetrabenazine.

The vesicular monoamine transporter 2 (VMAT2) is a proton-dependent antiporter responsible for loading monoamine neurotransmitters into synaptic vesicles. Dysregulation of VMAT2 can lead to several neuropsychiatric disorders including Parkinson's disease and schizophrenia. Furthermore, drugs such as amphetamine and MDMA are known to act on VMAT2, exemplifying its role in the mechanisms of actions for drugs of abuse. Despite VMAT2's importance, there remains a critical lack of mechanistic understanding, largely driven by a lack of structural information. Here we report a 3.1 Å resolution cryo-EM structure of VMAT2 complexed with tetrabenazine (TBZ), a non-competitive inhibitor used in the treatment of Huntington's chorea. We find TBZ interacts with residues in a central binding site, locking VMAT2 in an occluded conformation and providing a mechanistic basis for non-competitive inhibition. We further identify residues critical for cytosolic and lumenal gating, including a cluster of hydrophobic residues which are involved in a lumenal gating strategy. Our structure also highlights three distinct polar networks that may determine VMAT2 conformational dynamics and play a role in proton transduction. The structure elucidates mechanisms of VMAT2 inhibition and transport, providing insights into VMAT2 architecture, function, and the design of small-molecule therapeutics.

Tandem-repeat proteins conformational mechanics are optimized to facilitate functional interactions and complexations.

The architectures of tandem-repeat proteins are distinct from those of globular proteins. Individual modules, each comprising small structural motifs of 20-40 residues, are arrayed in a quasi one-dimensional fashion to form striking, elongated, horseshoe-like, and superhelical architectures, stabilized solely by short-range interaction. The spring-like shapes of repeat arrays point to elastic modes of action, and these proteins function as adapter molecules or 'hubs,' propagating signals within multi-subunit assemblies in diverse biological contexts. This flexibility is apparent in the dramatic variability observed in the structures of tandem-repeat proteins in different complexes. Here, using computational analysis, we demonstrate the striking ability of just one or a few global motions to recapitulate these structures. These findings show how the mechanics of repeat arrays are robustly enabled by their unique architecture. Thus, the repeating architecture has been optimized by evolution to favor functional modes of motions. The global motions enabling functional transitions can be fully visualized at http://bahargroup.org/tr_web.

Structural basis of G protein-Coupled receptor CMKLR1 activation and signaling induced by a chemerin-derived agonist.

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Influence of Point Mutations on PR65 Conformational Adaptability: Insights from Nanoaperture Optical Tweezer Experiments and Molecular Simulations.

PR65 is the HEAT-repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem-repeat protein, forming a spring-like architecture. PR65 conformational mechanics play a crucial role in PP2A function by opening/closing the substrate-binding/catalysis interface. Using in-silico saturation mutagenesis we identified "hinge" residues of PR65, whose substitutions are predicted to restrict its conformational adaptability and thereby disrupt PP2A function. Molecular simulations revealed that a subset of hinge mutations stabilized the extended/open conformation, whereas another had the opposite effect. By trapping in nanoaperture optical tweezer, we characterized PR65 motion and showed that the former mutants exhibited higher corner frequencies and lower translational scattering, indicating a shift towards extended conformations, whereas the latter showed the opposite behavior. Thus, experiments confirm the conformations predicted computationally. The study highlights the utility of nanoaperture-based tweezers for exploring structure and dynamics, and the power of integrating this single-molecule method with in silico approaches.

Anomalous peroxidase activity of cytochrome c is the primary pathogenic target in Barth syndrome.

Barth syndrome (BTHS) is a life-threatening genetic disorder with unknown pathogenicity caused by mutations in TAFAZZIN (TAZ) that affect remodeling of mitochondrial cardiolipin (CL). TAZ deficiency leads to accumulation of mono-lyso-CL (MLCL), which forms a peroxidase complex with cytochrome c (cyt c) capable of oxidizing polyunsaturated fatty acid-containing lipids. We hypothesized that accumulation of MLCL facilitates formation of anomalous MLCL-cyt c peroxidase complexes and peroxidation of polyunsaturated fatty acid phospholipids as the primary BTHS pathogenic mechanism. Using genetic, biochemical/biophysical, redox lipidomic and computational approaches, we reveal mechanisms of peroxidase-competent MLCL-cyt c complexation and increased phospholipid peroxidation in different TAZ-deficient cells and animal models and in pre-transplant biopsies from hearts of patients with BTHS. A specific mitochondria-targeted anti-peroxidase agent inhibited MLCL-cyt c peroxidase activity, prevented phospholipid peroxidation, improved mitochondrial respiration of TAZ-deficient C2C12 myoblasts and restored exercise endurance in a BTHS Drosophila model. Targeting MLCL-cyt c peroxidase offers therapeutic approaches to BTHS treatment.

Membrane regulation of 15LOX-1/PEBP1 complex prompts the generation of ferroptotic signals, oxygenated PEs.

Ferroptosis is a regulated form of cell death, the mechanism of which is still to be understood. 15-lipoxygenase (15LOX) complex with phosphatidylethanolamine (PE)-binding protein 1 (PEBP1) catalyzes the generation of pro-ferroptotic cell death signals, hydroperoxy-polyunsaturated PE. We focused on gaining new insights into the molecular basis of these pro-ferroptotic interactions using computational modeling and liquid chromatography-mass spectrometry experiments. Simulations of 15LOX-1/PEBP1 complex dynamics and interactions with lipids revealed that association with the membrane triggers a conformational change in the complex. This conformational change facilitates the access of stearoyl/arachidonoyl-PE (SAPE) substrates to the catalytic site. Furthermore, the binding of SAPE promotes tight interactions within the complex and induces further conformational changes that facilitate the oxidation reaction. The reaction yields two hydroperoxides as products, 15-HpETE-PE and 12-HpETE-PE, at a ratio of 5:1. A significant effect of PEBP1 is observed only on the predominant product. Moreover, combined experiments and simulations consistently demonstrate the significance of PEBP1 P112E mutation in generating ferroptotic cell death signals.

Structural basis of CMKLR1 signaling induced by chemerin9.

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Discovering selective antiferroptotic inhibitors of the 15LOX/PEBP1 complex noninterfering with biosynthesis of lipid mediators.

Programmed ferroptotic death eliminates cells in all major organs and tissues with imbalanced redox metabolism due to overwhelming iron-catalyzed lipid peroxidation under insufficient control by thiols (Glutathione (GSH)). Ferroptosis has been associated with the pathogenesis of major chronic degenerative diseases and acute injuries of the brain, cardiovascular system, liver, kidneys, and other organs, and its manipulation offers a promising new strategy for anticancer therapy. This explains the high interest in designing new small-molecule-specific inhibitors against ferroptosis. Given the role of 15-lipoxygenase (15LOX) association with phosphatidylethanolamine (PE)-binding protein 1 (PEBP1) in initiating ferroptosis-specific peroxidation of polyunsaturated PE, we propose a strategy of discovering antiferroptotic agents as inhibitors of the 15LOX/PEBP1 catalytic complex rather than 15LOX alone. Here we designed, synthesized, and tested a customized library of 26 compounds using biochemical, molecular, and cell biology models along with redox lipidomic and computational analyses. We selected two lead compounds, FerroLOXIN-1 and 2, which effectively suppressed ferroptosis in vitro and in vivo without affecting the biosynthesis of pro-/anti-inflammatory lipid mediators in vivo. The effectiveness of these lead compounds is not due to radical scavenging or iron-chelation but results from their specific mechanisms of interaction with the 15LOX-2/PEBP1 complex, which either alters the binding pose of the substrate [eicosatetraenoyl-PE (ETE-PE)] in a nonproductive way or blocks the predominant oxygen channel thus preventing the catalysis of ETE-PE peroxidation. Our successful strategy may be adapted to the design of additional chemical libraries to reveal new ferroptosis-targeting therapeutic modalities.