Effects of ligand binding on the stability of aldo-keto reductases: Implications for stabilizer or destabilizer chaperones

A new lysine biosynthetic enzyme from a bacterial endosymbiont shaped by genetic drift and genome reduction

The effect of population bottlenecks and genome reduction on enzyme function is poorly understood. Candidatus Liberibacter solanacearum is a bacterium with a reduced genome that is transmitted vertically to the egg of an infected psyllid-a population bottleneck that imposes genetic drift and is predicted to affect protein structure and function. Here, we define the function of Ca. L. solanacearum dihydrodipicolinate synthase (CLsoDHDPS), which catalyzes the committed branchpoint reaction in diaminopimelate and lysine biosynthesis. We demonstrate that CLsoDHDPS is expressed in Ca. L. solanacearum and expression is increased ~2-fold in the insect host compared to in planta. CLsoDHDPS has decreased thermal stability and increased aggregation propensity, implying mutations have destabilized the enzyme but are compensated for through elevated chaperone expression and a stabilized oligomeric state. CLsoDHDPS uses a ternary-complex kinetic mechanism, which is to date unique among DHDPS enzymes, has unusually low catalytic ability, but an unusually high substrate affinity. Structural studies demonstrate that the active site is more open, and the structure of CLsoDHDPS with both pyruvate and the substrate analogue succinic-semialdehyde reveals that the product is both structurally and energetically different and therefore evolution has in this case fashioned a new enzyme. Our study suggests the effects of genome reduction and genetic drift on the function of essential enzymes and provides insights on bacteria-host co-evolutionary associations. We propose that bacteria with endosymbiotic lifestyles present a rich vein of interesting enzymes useful for understanding enzyme function and/or informing protein engineering efforts.

A three-level regulatory mechanism of the aldo-keto reductase subfamily AKR12D

Modulation of protein function through allosteric regulation is central in biology, but biomacromolecular systems involving multiple subunits and ligands may exhibit complex regulatory mechanisms at different levels, which remain poorly understood. Here, we discover an aldo-keto reductase termed AKRtyl and present its three-level regulatory mechanism. Specifically, by combining steady-state and transient kinetics, X-ray crystallography and molecular dynamics simulation, we demonstrate that AKRtyl exhibits a positive synergy mediated by an unusual Monod-Wyman-Changeux (MWC) paradigm of allosteric regulation at low concentrations of the cofactor NADPH, but an inhibitory effect at high concentrations is observed. While the substrate tylosin binds at a remote allosteric site with positive cooperativity. We further reveal that these regulatory mechanisms are conserved in AKR12D subfamily, and that substrate cooperativity is common in AKRs across three kingdoms of life. This work provides an intriguing example for understanding complex allosteric regulatory networks.

Biochemical and Structural Analysis of the Bacterial Enzyme Succinyl-Diaminopimelate Desuccinylase (DapE) from Acinetobacter baumannii

There is an urgent need for new antibiotics given the rise of antibiotic resistance, and succinyl-diaminopimelate desuccinylase (DapE, E.C. 3.5.1.18) has emerged as a promising bacterial enzyme target. DapE from Haemophilus influenzae (HiDapE) has been studied and inhibitors identified, but it is essential to explore DapE from different species to assess selective versus broad-spectrum therapeutics. We have determined the structure of DapE from the ESKAPE pathogen Acinetobacter baumannii (AbDapE) and studied inhibition by known inhibitors of HiDapE. AbDapE is inhibited by captopril and sulfate comparable to HiDapE, but AbDapE was not significantly inhibited by a known indoline sulfonamide HiDapE inhibitor. Captopril and sulfate both stabilize HiDapE by increasing the thermal melting temperature (Tm) in thermal shift assays. By contrast, sulfate decreases the stability of the AbDapE enzyme, whereas captopril increases the stability. Further, we report two crystal structures of selenomethionine-substituted AbDapE in the closed conformation, one with AbDapE in complex with succinate derived from enzymatic hydrolysis of N6-methyl-l,l-SDAP substrate and acetate (PDB code 7T1Q, 2.25 Å resolution), and a crystal structure of AbDapE with bound succinate along with l-(S)-lactate, a product of degradation of citric acid from the crystallization buffer during X-ray irradiation (PDB code 8F8O, 2.10 Å resolution).

Cyclobutanone Inhibitors of Diaminopimelate Desuccinylase (DapE) as Potential New Antibiotics

Based on our previous success in using cyclobutanone derivatives as enzyme inhibitors, we have designed and prepared a 37-member library of α-aminocyclobutanone amides and sulfonamides, screened for inhibition of the bacterial enzyme diaminopimelate desuccinylase (DapE), which is a promising antibiotic target, and identified several inhibitors with micromolar inhibitory potency. Molecular docking suggests binding of the deprotonated hydrate of the strained cyclobutanone, and thermal shift analysis with the most potent inhibitor (3y, IC50 = 23.1 µM) enabled determination of a Ki value of 10.2 +/- 0.26 µM and observed two separate Tm values for H. influenzae DapE (HiDapE).

Repurposing of US-FDA approved drugs against SARS-CoV-2 main protease (Mpro) by using STD-NMR spectroscopy, in silico studies and antiviral assays

SARS-CoV-2 Main protease (M^pro) is a well-known drug target against SARS-CoV-2 infection. Identification of M^pro inhibitors is vigorously pursued due to its crucial role in viral replication. The present study was aimed to identify M^pro inhibitors via repurposing of US-FDA approved drugs by STD-NMR spectroscopy. In this study, 156 drugs and natural compounds were evaluated against M^pro. Among them, 10 drugs were found to be interacting with M^pro, including diltiazem HCl (1), mefenamic acid (2), losartan potassium (3), mexiletine HCl (4), glaucine HBr (5), trimebutine maleate (6), flurbiprofen (7), amantadine HCl (8), dextromethorphan (9), and lobeline HCl (10) in STD-NMR spectroscopy. Their interactions were compared with three standards (Repurposed anti-viral drugs), dexamethasone, chloroquine phosphate, and remdesivir. Thermal stability of M^pro and dissociation constant (K_d) of six interacting drugs were also determined using DSF. RMSD plots in MD simulation studies showed the formation of stable protein-ligand complexes. They were further examined for their antiviral activity by plaque reduction assay against SARS-CoV-2, which showed 55-100% reduction in viral plaques. This study demonstrates the importance of drug repurposing against emerging and neglected diseases. This study also exhibits successful application of STD-NMR spectroscopy combined with plaque reduction assay in rapid identification of potential anti-viral agents.

Identification of a Bile Acid-Binding Transcription Factor in Clostridioides difficile Using Chemical Proteomics

Clostridioides difficile is a Gram-positive anaerobic bacterium that is the leading cause of hospital-acquired gastroenteritis in the US. In the gut milieu, C. difficile encounters microbiota-derived, growth-inhibiting bile acids that are thought to be a significant mechanism of colonization resistance. While the levels of certain bile acids in the gut correlate with susceptibility to C. difficile infection, their molecular targets in C. difficile remain unknown. In this study, we sought to use chemical proteomics to identify bile acid-interacting proteins in C. difficile. Using photoaffinity bile acid probes and chemical proteomics, we identified a previously uncharacterized MerR family protein, CD3583 (now BapR), as a putative bile acid-sensing transcription regulator. Our data indicate that BapR specifically binds to and is stabilized by lithocholic acid (LCA) in C. difficile. Although loss of BapR did not affect C. difficile's sensitivity to LCA, ΔbapR cells elongated more in the presence of LCA compared to wild-type cells. Transcriptomics revealed that BapR regulates several gene clusters, with the expression of the mdeA-cd3573 locus being specifically de-repressed in the presence of LCA in a BapR-dependent manner. Electrophoretic mobility shift assays revealed that BapR directly binds to the mdeA promoter region. Because mdeA is involved in amino acid-related sulfur metabolism and the mdeA-cd3573 locus encodes putative transporters, we propose that BapR senses a gastrointestinal tract-specific small molecule, LCA, as an environmental cue for metabolic adaptation.

Chalcones identify cTXNPx as a potential antileishmanial drug target

With current drug treatments failing due to toxicity, low efficacy and resistance; leishmaniasis is a major global health challenge that desperately needs new validated drug targets. Inspired by activity of the natural chalcone 2’,6’-dihydroxy-4’-methoxychalcone (DMC), the nitro-analogue, 3-nitro-2’,4’,6’- trimethoxychalcone (NAT22, 1c) was identified as potent broad spectrum antileishmanial drug lead. Structural modification provided an alkyne containing chemical probe that labelled a protein within the parasite that was confirmed as cytosolic tryparedoxin peroxidase (cTXNPx). Crucially, labelling is observed in both promastigote and intramacrophage amastigote life forms, with no evidence of host macrophage toxicity. Incubation of the chalcone in the parasite leads to ROS accumulation and parasite death. Deletion of cTXNPx, by CRISPR-Cas9, dramatically impacts upon the parasite phenotype and reduces the antileishmanial activity of the chalcone analogue. Molecular docking studies with a homology model of in-silico cTXNPx suggest that the chalcone is able to bind in the putative active site hindering access to the crucial cysteine residue. Collectively, this work identifies cTXNPx as an important target for antileishmanial chalcones.

Leishmaniasis is an insect vector-borne parasitic disease. With >350 million people world wide considered at risk, 12 million people currently infected and an economic cost that can be estimated in terms of >3.3 million working life years lost, leishmaniasis is a major global health challenge. The disease is of particular importance in Brazil. Current treatment of leishmaniasis is difficult requiring a long, costly course of drug treatment using old drugs with poor safety indications requiring close medical supervision. Moreover, resistance to current antileishmanials is growing, emphasising a major need for new drug targets. In earlier work we had identified a naturally inspired chalcone which had promising antileishmanial activity but with no known mode of action. In this work we use an analogue of this molecule as an activity based probe to identify a protein target of the chalcone. This protein, cTXNPx, has a major role in protecting the parasite against attack by reactive oxygen species in the host cell. By inhibiting this protein the parasite can no longer survive in the host. Collectively this work validates cTXNPx as a drug target with the chalcone as a lead structure for future drug discovery programmes.

Drug target inference by mining transcriptional data using a novel graph convolutional network framework

A fundamental challenge that arises in biomedicine is the need to characterize compounds in a relevant cellular context in order to reveal potential on-target or off-target effects. Recently, the fast accumulation of gene transcriptional profiling data provides us an unprecedented opportunity to explore the protein targets of chemical compounds from the perspective of cell transcriptomics and RNA biology. Here, we propose a novel Siamese spectral-based graph convolutional network (SSGCN) model for inferring the protein targets of chemical compounds from gene transcriptional profiles. Although the gene signature of a compound perturbation only provides indirect clues of the interacting targets, and the biological networks under different experiment conditions further complicate the situation, the SSGCN model was successfully trained to learn from known compound-target pairs by uncovering the hidden correlations between compound perturbation profiles and gene knockdown profiles. On a benchmark set and a large time-split validation dataset, the model achieved higher target inference accuracy as compared to previous methods such as Connectivity Map. Further experimental validations of prediction results highlight the practical usefulness of SSGCN in either inferring the interacting targets of compound, or reversely, in finding novel inhibitors of a given target of interest.

Screening of benzenesulfonamide in combination with chemically diverse fragments against carbonic anhydrase by differential scanning fluorimetry

The differential scanning fluorimetry (DSF) screening of 5.692 fragments in combination with benzenesulfonamide (BSA) against bovine carbonic anhydrase (bCA) delivered >100 hits that either caused, on their own, a significant thermal shift (ΔTm, °C) in the protein melting temperature or significantly influenced the thermal shift observed for BSA alone. Three hits based on 1,2,3-triazole moiety represent the periphery of the recently reported potent inhibitors of hCA II, IX and XII which were efficacious in vivo. Such a re-discovery of suitable BSA periphery essentially validates the new fragment-based approach to the discovery of future CAIs. Structures of other validated fragment hits are reported.

The Principles of protein surgery and its application to the logical drug design for the treatment of neurodegenerative diseases

Here we describe the principles of protein surgery and its application to the molecular design for regulating the protein conformation. We initially describe the Poincare duality that defines the basis of complementarity in the time-dependent geometrical space. Next we introduce the theory of protein surgery consisting of "dissection" and "suture," which correspond to differentiation and integration in a phase space, respectively. Then we introduce the surgical (pan-manifold) differential equation and solve several simple cases. In this way, we constructed "medical quantization" strategy which makes a bridge between quantum mechanics and molecular biology. As the application of this theory, here we describe a logical drug design methodology in detail, and its application to anti-prion drug design. Finally, we propose a plan for logical drug (and medical device) design center, in which all the necessary procedures could be done in one place.

A missense variant remote from the active site impairs stability of human phosphoglucomutase 1

Missense variants of human phosphoglucomutase 1 (PGM1) cause the inherited metabolic disease known as PGM1 deficiency. This condition is categorised as both a glycogen storage disease and a congenital disorder of glycosylation. Approximately 20 missense variants of PGM1 are linked to PGM1 deficiency, and biochemical studies have suggested that they fall into two general categories: those affecting the active site and catalytic efficiency, and those that appear to impair protein folding and/or stability. In this study, we characterise a novel variant of Arg422, a residue distal from the active site of PGM1 and the site of a previously identified disease-related variant (Arg422Trp). In prior studies, the R422W variant was found to produce insoluble protein in a recombinant expression system, precluding further in vitro characterisation. Here we investigate an alternative variant of this residue, Arg422Gln, which is amenable to experimental characterisation presumably due to its more conservative physicochemical substitution. Biochemical, crystallographic, and computational studies of R422Q establish that this variant causes only minor changes in catalytic efficiency and 3D structure, but is nonetheless dramatically reduced in stability. Unexpectedly, binding of a substrate analog is found to further destabilise the protein, in contrast to its stabilising effect on wild-type PGM1 and several other missense variants. This work establishes Arg422 as a lynchpin residue for the stability of PGM1 and supports the impairment of protein stability as a pathomechanism for variants that cause PGM1 deficiency. SYNOPSIS: Biochemical and structural studies of a missense variant far from the active site of human PGM1 identify a residue with a key role in enzyme stability.

Design, Synthesis, and Conformation-Activity Study of Unnatural Bridged Bicyclic Depsipeptides as Highly Potent Hypoxia Inducible Factor-1 Inhibitors and Antitumor Agents

By carrying out structural modifications based on the bicyclic peptide structure of echinomycin, we successfully synthesized various powerful antitumor derivatives. The ring conformation in the obtained compounds was restricted by cross-linking with an unnatural bond. The prepared derivatives were demonstrated to strongly suppress the hypoxia inducible factor (HIF)-1 transcriptional activation and hypoxia induction of HIF-1 protein expression. Particularly, alkene-bridged derivative 12 exhibited remarkably potent cytotoxicity (IC₅₀ = 0.22 nM on the MCF-7 cell line) and HIF-1 inhibition (IC₅₀ = 0.09 nM), which considerably exceeded those of echinomycin. Conformational analyses and molecular modeling studies revealed that the biological activities were enhanced following restriction of the conformation by cross-linking through a metabolically stable and rigid bridge bond. In addition, we proposed a new globular conformation stabilized by intramolecular π stacking that can contribute to the biological effects of bicyclic depsipeptides. The developments presented in the current study serve as a useful guide to expand the chemical space of peptides in drug discovery.

Discovery of a multipotent chaperone, 1-(2,6-Difluorobenzylamino)-3-(1,2,3,4-tetrahydrocarbazol-9-yl)-propan-2-ol with the inhibitory effects on the proliferation of prion, cancer as well as influenza virus

We previously discovered three carbazole derivatives, GJP14 (1-piperidinylmethyl-2-(1-oxo-6-methyl-1,2,3,4-tetrahydrocarbazol-9-yl)-ethan-1-ol) with anti-prion activity, GJC29 (benzylamino-3-(1,2,3,4-tetrahydrocarbazol-9-yl)-propan-2-ol) with anti-cancer activity, and THC19 (1-piperidinylmethyl-2-(1,2,3,4-tetrahydrocarnazol-9-yl)-ethan-1-ol) with anti-influenza virus activity. During optimization of GJP14 for the anti-prion activity, we discovered a compound, 1-(2,6-difluorobenzylamino)-3-(1,2,3,4-tetrahydrocarbazol-9-yl)-propan-2-ol, termed 5Y, had the most strong anti-prion activity among a series of newly synthesized derivatives. Intriguingly, we noticed that 5Y had also the most strong anti-colon cancer as well as the anti-influenza virus activities among derivatives. No significant toxicity of 5Y was observed. These results demonstrate that 5Y is a multipotent lead compound with unusually wide spectrum, and may be applicable to therapeutics targeting multiple diseases.Abbreviations: MoPrP: mouse prion protein of amino acid residues of 23-231; PrP^C: cellular form of prion protein; PrP^Sc: scrapie form of prion protein.

Deploying Fluorescent Nucleoside Analogues for High-Throughput Inhibitor Screening

High-throughput small-molecule screening in drug discovery processes commonly rely on fluorescence-based methods including fluorescent polarization and fluorescence/Förster resonance energy transfer. These techniques use highly accessible instrumentation; however, they can suffer from high false-negative rates and background signals, or might involve complex schemes for the introduction of fluorophore pairs. Herein we present the synthesis and application of fluorescent nucleoside analogues as the foundation for directed approaches for competitive binding analyses. The general approach describes selective fluorescent environment-sensitive (ES) nucleoside analogues that are adaptable to diverse enzymes that act on nucleoside-based substrates. We demonstrate screening a set of uridine analogues and development of an assay for fragment-based lead discovery with the TcdB glycosyltransferase (GT), an enzyme associated with virulence in Clostridium difficile. The uridine-based probe used for this high-throughput screen has a KD value of 7.2 μm with the TcdB GT and shows a >30-fold increase in fluorescence intensity upon binding. The ES-based probe assay is benchmarked against two other screening approaches.

Human aldose reductase unfolds through an intermediate

Background: Human aldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway.  Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR.

Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. 

Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR.

Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state.

Expression, purification and stabilization of human serotonin transporter from E. coli

The serotonin transporter belongs to the family of sodium-chloride coupled neurotransmitter transporter and is related to depression in humans. It is therefore an important drug target to support treatment of depression. Recently, structures of human serotonin transporter in complex with inhibitor molecules have been published. However, the production of large protein amounts for crystallization experiments remains a bottleneck. Here we present the possibility to obtain purified serotonin transporter from E. coli. Fos-choline 12 solubilized target protein was obtained with a purity of >95% and a yield of 1.2 mg L^-1 culture in autoinduction medium. CD spectroscopic analysis of protein stability allowed identifying CHS and POPX as stabilizing components, which increased hSERT thermostability by 7 °C. The kinetic dissociation constant K_D of 2.8 μM (±0.05) for of the inhibitor Desipramine was determined with a ka of 10,848 M - 1 s-1 (±220) and a kd of 0.03 s-1 (±4.7 × 10-5).

Structural characterization of life-extending Caenorhabditis elegans Lipid Binding Protein 8

The lysosome plays a crucial role in the regulation of longevity. Lysosomal degradation is tightly coupled with autophagy that is induced by many longevity paradigms and required for lifespan extension. The lysosome also serves as a hub for signal transduction and regulates longevity via affecting nuclear transcription. One lysosome-to-nucleus retrograde signaling pathway is mediated by a lysosome-associated fatty acid binding protein LBP-8 in Caenorhabditis elegans. LBP-8 shuttles lysosomal lipids into the nucleus to activate lipid regulated nuclear receptors NHR-49 and NHR-80 and consequently promote longevity. However, the structural basis of LBP-8 action remains unclear. Here, we determined the first 1.3 Å high-resolution structure of this life-extending protein LBP-8, which allowed us to identify a structurally conserved nuclear localization signal and amino acids involved in lipid binding. Additionally, we described the range of fatty acids LBP-8 is capable of binding and show that it binds to life-extending ligands in worms such as oleic acid and oleoylethanolamide with high affinity.

Melting Down Protein Stability: PAPS Synthase 2 in Patients and in a Cellular Environment

Within the crowded and complex environment of the cell, a protein experiences stabilizing excluded-volume effects and destabilizing quinary interactions with other proteins. Which of these prevail, needs to be determined on a case-by-case basis. PAPS synthases are dimeric and bifunctional enzymes, providing activated sulfate in the form of 3′-phosphoadenosine-5′-phosphosulfate (PAPS) for sulfation reactions. The human PAPS synthases PAPSS1 and PAPSS2 differ significantly in their protein stability as PAPSS2 is a naturally fragile protein. PAPS synthases bind a series of nucleotide ligands and some of them markedly stabilize these proteins. PAPS synthases are of biomedical relevance as destabilizing point mutations give rise to several pathologies. Genetic defects in PAPSS2 have been linked to bone and cartilage malformations as well as a steroid sulfation defect. All this makes PAPS synthases ideal to study protein unfolding, ligand binding, and the stabilizing and destabilizing factors in their cellular environment. This review provides an overview on current concepts of protein folding and stability and links this with our current understanding of the different disease mechanisms of PAPSS2-related pathologies with perspectives for future research and application.

A designer molecular chaperone against transmissible spongiform encephalopathy slows disease progression in mice and macaques

Transmissible spongiform encephalopathies (TSEs) are fatal neurodegenerative diseases that lack therapeutic solutions. Here, we show that the molecular chaperone (N,N'-([cyclohexylmethylene]di-4,1-phenylene)bis(2-[1-pyrrolidinyl]acetamide)), designed via docking simulations, molecular dynamics simulations and quantum chemical calculations, slows down the progress of TSEs. In vitro, the designer molecular chaperone stabilizes the normal cellular prion protein, eradicates prions in infected cells, prevents the formation of drug-resistant strains and directly inhibits the interaction between prions and abnormal aggregates, as shown via real-time quaking-induced conversion and in vitro conversion NMR. Weekly intraperitoneal injection of the chaperone in prion-infected mice prolonged their survival, and weekly intravenous administration of the compound in macaques infected with bovine TSE slowed down the development of neurological and psychological symptoms and reduced the concentration of disease-associated biomarkers in the animals' cerebrospinal fluid. The de novo rational design of chaperone compounds could lead to therapeutics that can bind to different prion protein strains to ameliorate the pathology of TSEs.

Piperine, an alkaloid inhibiting the super-relaxed state of myosin, binds to the myosin regulatory light chain

Piperine, an alkaloid from black pepper, was found to inhibit the super-relaxed state (SRX) of myosin in fast-twitch skeletal muscle fibers. In this work we report that the piperine molecule binds heavy meromyosin (HMM), whereas it does not interact with the regulatory light chain (RLC)-free subfragment-1 (S1) or with control proteins from the same muscle molecular machinery, G-actin and tropomyosin. To further narrow down the location of piperine binding, we studied interactions between piperine and a fragment of skeletal myosin consisting of the full-length RLC and a fragment of the heavy chain (HCF). The sequence of HCF was designed to bind RLC and to dimerize via formation of a stable coiled coil, thus producing a well-folded isolated fragment of the myosin neck. Both chains were co-expressed in Escherichia coli, the RLC/HCF complex was purified and tested for stability, composition and binding to piperine. RLC and HCF chains formed a stable heterotetrameric complex (RLC/HCF)₂ which was found to bind piperine. The piperine molecule was also found to bind isolated RLC. Piperine binding to RLC in (RLC/HCF)₂ altered the compactness of the complex, suggesting that the mechanism of SRX inhibition by piperine is based on changing conformation of the myosin.

Hypoxia-selective allosteric destabilization of activin receptor-like kinases: A potential therapeutic avenue for prophylaxis of heterotopic ossification

Heterotopic ossification (HO), the pathological extraskeletal formation of bone, can arise from blast injuries, severe burns, orthopedic procedures and gain-of-function mutations in a component of the bone morphogenetic protein (BMP) signaling pathway, the ACVR1/ALK2 receptor serine-threonine (protein) kinase, causative of Fibrodysplasia Ossificans Progressiva (FOP). All three ALKs (-2, -3, -6) that play roles in bone morphogenesis contribute to trauma-induced HO, hence are well-validated pharmacological targets. That said, development of inhibitors, typically competitors of ATP binding, is inherently difficult due to the conserved nature of the active site of the 500+ human protein kinases. Since these enzymes are regulated via inherent plasticity, pharmacological chaperone-like drugs binding to another (allosteric) site could hypothetically modulate kinase conformation and activity. To test for such a mechanism, a surface pocket of ALK2 kinase formed largely by a key allosteric substructure was targeted by supercomputer docking of drug-like compounds from a virtual library. Subsequently, the effects of docked hits were further screened in vitro with purified recombinant kinase protein. A family of compounds with terminal hydrogen-bonding acceptor groups was identified that significantly destabilized the protein, inhibiting activity. Destabilization was pH-dependent, putatively mediated by ionization of a histidine within the allosteric substructure with decreasing pH. In vivo, nonnative proteins are degraded by proteolysis in the proteasome complex, or cellular trashcan, allowing for the emergence of therapeutics that inhibit through degradation of over-active proteins implicated in the pathology of diseases and disorders. Because HO is triggered by soft-tissue trauma and ensuing hypoxia, dependency of ALK destabilization on hypoxic pH imparts selective efficacy on the allosteric inhibitors, providing potential for safe prophylactic use.

Genetic Architecture of the Cardiovascular Risk Proteome

BACKGROUND

We recently identified 156 proteins in human plasma that were each associated with the net Framingham Cardiovascular Disease Risk Score using an aptamer-based proteomic platform in Framingham Heart Study Offspring participants. Here we hypothesized that performing genome-wide association studies and exome array analyses on the levels of each of these 156 proteins might identify genetic determinants of risk-associated circulating factors and provide insights into early cardiovascular pathophysiology.

METHODS

We studied the association of genetic variants with the plasma levels of each of the 156 Framingham Cardiovascular Disease Risk Score-associated proteins using linear mixed-effects models in 2 population-based cohorts. We performed discovery analyses on plasma samples from 759 participants of the Framingham Heart Study Offspring cohort, an observational study of the offspring of the original Framingham Heart Study and their spouses, and validated these findings in plasma samples from 1421 participants of the MDCS (Malmö Diet and Cancer Study). To evaluate the utility of this strategy in identifying new biological pathways relevant to cardiovascular disease pathophysiology, we performed studies in a cell-model system to experimentally validate the functional significance of an especially novel genetic association with circulating apolipoprotein E levels.

RESULTS

We identified 120 locus-protein associations in genome-wide analyses and 41 associations in exome array analyses, the majority of which have not been described previously. These loci explained up to 66% of interindividual plasma protein-level variation and, on average, accounted for 3 times the amount of variation explained by common clinical factors, such as age, sex, and diabetes mellitus status. We described overlap among many of these loci and cardiovascular disease genetic risk variants. Finally, we experimentally validated a novel association between circulating apolipoprotein E levels and the transcription factor phosphatase 1G. Knockdown of phosphatase 1G in a human liver cell model resulted in decreased apolipoprotein E transcription and apolipoprotein E protein levels in cultured supernatants.

CONCLUSIONS

We identified dozens of novel genetic determinants of proteins associated with the Framingham Cardiovascular Disease Risk Score and experimentally validated a new role for phosphatase 1G in lipoprotein biology. Further, genome-wide and exome array data for each protein have been made publicly available as a resource for cardiovascular disease research.