Soumya S. Ray

Triosephosphate isomerase from Plasmodium falciparum: the crystal structure provides insights into antimalarial drug design.

BACKGROUND Malaria caused by the parasite Plasmodium falciparum is a major public health concern. The parasite lacks a functional tricarboxylic acid cycle, making glycolysis its sole energy source. Although parasite enzymes have been considered as potential antimalarial drug targets, little is known about their structural biology. Here we report the crystal structure of triosephosphate isomerase (TIM) from P. falciparum at 2.2 A resolution. RESULTS The crystal structure of P. falciparum TIM (PfTIM), expressed in Escherichia coli, was determined by the molecular replacement method using the structure of trypanosomal TIM as the starting model. Comparison of the PfTIM structure with other TIM structures, particularly human TIM, revealed several differences. In most TIMs the residue at position 183 is a glutamate but in PfTIM it is a leucine. This leucine residue is completely exposed and together with the surrounding positively charged patch, may be responsible for binding TIM to the erythrocyte membrane. Another interesting feature is the occurrence of a cysteine residue at the dimer interface of PfTIM (Cys13), in contrast to human TIM where this residue is a methionine. Finally, residue 96 of human TIM (Ser96), which occurs near the active site, has been replaced by phenylalanine in PfTIM. CONCLUSIONS Although the human and Plasmodium enzymes share 42% amino acid sequence identity, several key differences suggest that PfTIM may turn out to be a potential drug target. We have identified a region which may be responsible for binding PfTIM to cytoskeletal elements or the band 3 protein of erythrocytes; attachment to the erythrocyte membrane may subsequently lead to the extracellular exposure of parts of the protein. This feature may be important in view of a recent report that patients suffering from P. falciparum malaria mount an antibody response to TIM leading to prolonged hemolysis. A second approach to drug design may be provided by the mutation of the largely conserved residue (Ser96) to phenylalanine in PfTIM. This difference may be of importance in designing specific active-site inhibitors against the enzyme. Finally, specific inhibition of PfTIM subunit assembly might be possible by targeting Cys13 at the dimer interface. The crystal structure of PfTIM provides a framework for new therapeutic leads.

Kinetic, mechanistic, and structural modeling studies of truncated wild-type leucine-rich repeat kinase 2 and the G2019S mutant.

Leucine-rich repeat kinase 2 (LRRK2), a large and complex protein that possesses two enzymatic properties, kinase and GTPase, is one of the major genetic factors in Parkinsons disease (PD). Here, we characterize the kinetic and catalytic mechanisms of truncated wild-type (t-wt) LRRK2 and its most common mutant, G2019S (t-G2019S), with a structural interpretation of the kinase domain. First, the substitution of threonine with serine in the LRRKtide peptide results in a much less efficient substrate as demonstrated by a 26-fold decrease in k(cat) and a 6-fold decrease in binding affinity. The significant decrease in k(cat) is attributed to a slow chemical transfer step as evidenced by the inverse solvent kinetic isotope effect in the proton inventory and pL (pH or pD)-dependent studies. The shape of the proton inventory and pL profile clearly signals the involvement of a general base (pK(a) = 7.5) in the catalysis with a low fractionation factor in the ground state. We report for the first time that the increased kinase activity of the G2019S mutant is substrate-dependent. Homology modeling of the kinase domain (open and closed forms) and structural analysis of the docked peptide substrates suggest that electrostatic interactions play an important role in substrate recognition, which is affected by G2019S and may directly influence the kinetic properties of the enzyme. Finally, the GTPase activity of the t-G2019S mutant was characterized, and the mutation modestly decreases GTPase activity without significantly affecting GTP binding affinity.

Inflammatory Signaling by NOD-RIPK2 Is Inhibited by Clinically Relevant Type II Kinase Inhibitors.

Summary RIPK2 mediates pro-inflammatory signaling from the bacterial sensors NOD1 and NOD2, and is an emerging therapeutic target in autoimmune and inflammatory diseases. We observed that cellular RIPK2 can be potently inhibited by type II inhibitors that displace the kinase activation segment, whereas ATP-competitive type I inhibition was only poorly effective. The most potent RIPK2 inhibitors were the US Food and Drug Administration-approved drugs ponatinib and regorafenib. Their mechanism of action was independent of NOD2 interaction and involved loss of downstream kinase activation as evidenced by lack of RIPK2 autophosphorylation. Notably, these molecules also blocked RIPK2 ubiquitination and, consequently, inflammatory nuclear factor κB signaling. In monocytes, the inhibitors selectively blocked NOD-dependent tumor necrosis factor production without affecting lipopolysaccharide-dependent pathways. We also determined the first crystal structure of RIPK2 bound to ponatinib, and identified an allosteric site for inhibitor development. These results highlight the potential for type II inhibitors to treat indications of RIPK2 activation as well as inflammation-associated cancers.

Structure-activity relationship study of beta-carboline derivatives as haspin kinase inhibitors

Haspin is a serine/threonine kinase that phosphorylates Thr-3 of histone H3 in mitosis that has emerged as a possible cancer therapeutic target. High throughput screening of approximately 140,000 compounds identified the beta-carbolines harmine and harmol as moderately potent haspin kinase inhibitors. Based on information obtained from a structure-activity relationship study previously conducted for an acridine series of haspin inhibitors in conjunction with in silico docking using a recently disclosed crystal structure of the kinase, harmine analogs were designed that resulted in significantly increased haspin kinase inhibitory potency. The harmine derivatives also demonstrated less activity towards DYRK2 compared to the acridine series. In vitro mouse liver microsome stability and kinase profiling of a representative member of the harmine series (42, LDN-211898) are also presented.

An Electrospray Ionization Mass Spectrometry Investigation of 1-Anilino-8- Naphthalene-Sulfonate (ANS) Binding to Proteins

The binding of 1-anilino-8-naphthalene-sulfonic acid (ANS) to various globular proteins at acidic pH has been investigated by electrospray ionization mass spectrometry (ESI-MS). Maximal ANS binding is observed in the pH range 3–5. As many as seven species of dye-bound complexes are detected for myoglobin. Similar studies were carried out with cytochrome c, carbonic anhydrase, triosephosphate isomerase, lysozyme, α-lactalbumin, and bovine pancreatic trypsin inhibitor (BPTI). Strong ANS binding was observed wherever molten globule states were postulated in solution. ANS binding is not observed for lysozyme and BPTI, which have tightly folded structures in the native form. α-Lactalbumin, which is structurally related to lysozyme but forms a molten globule at acidic pH, exhibited ANS binding. Reduction of disulfide bonds in these proteins leads to the detection of ANS binding even at neutral pH. Binding was suppressed at very low pH (<2.5), presumably due to neutralization of the charge on the sulfonate moiety. The distribution of the relative intensities of the protein bound ANS species varies with the charge state, suggesting heterogeneity of gas phase conformations. The binding strength of these complexes was qualitatively estimated by dissociating them using enhanced nozzle skimmer potentials. The skimmer voltages also affected the lower and higher charge states of these complexes in a different manner.

The Parkinson Disease-linked LRRK2 Protein Mutation I2020T Stabilizes an Active State Conformation Leading to Increased Kinase Activity

Background: LRRK2 has emerged as one of the most relevant players in Parkinson disease pathogenesis. Results: Enzyme kinetic and modeling studies elucidate the effect of the LRRK2 I2020T mutant on kinase activity. Conclusion: The I2020T mutant stabilizes the active conformation and leads to increased kinase activity. As a result, it binds to type II inhibitors competitively. Significance: This study may contribute to the development of new classes of inhibitors of LRRK2. The effect of leucine-rich repeat kinase 2 (LRRK2) mutation I2020T on its kinase activity has been controversial, with both increased and decreased effects being reported. We conducted steady-state and pre-steady-state kinetic studies on LRRKtide and its analog LRRKtideS. Their phosphorylation differs by the rate-limiting steps: product release is rate-limiting for LRRKtide and phosphoryl transfer is rate-limiting for LRRKtideS. As a result, we observed that the I2020T mutant is more active than wild type (WT) LRRK2 for LRRKtideS phosphorylation, whereas it is less active than WT for LRRKtide phosphorylation. Our pre-steady-state kinetic data suggest that (i) the I2020T mutant accelerates the rates of phosphoryl transfer of both reactions by 3–7-fold; (ii) this increase is masked by a rate-limiting product release step for LRRKtide phosphorylation; and (iii) the observed lower activity of the mutant for LRRKtide phosphorylation is a consequence of its instability: the concentration of the active form of the mutant is 3-fold lower than WT. The I2020T mutant has a dramatically low KATP and therefore leads to resistance to ATP competitive inhibitors. Two well known DFG-out or type II inhibitors are also weaker toward the mutant because they inhibit the mutant in an unexpected ATP competitive mechanism. The I2020 residue lies next to the DYG motif of the activation loop of the LRRK2 kinase domain. Our modeling and metadynamic simulations suggest that the I2020T mutant stabilizes the DYG-in active conformation and creates an unusual allosteric pocket that can bind type II inhibitors but in an ATP competitive fashion.

Improving binding specificity of pharmacological chaperones that target mutant superoxide dismutase-1 linked to familial amyotrophic lateral sclerosis using computational methods.

We recently described a set of drug-like molecules obtained from an in silico screen that stabilize mutant superoxide dismutase-1 (SOD-1) linked to familial amyotrophic lateral sclerosis (ALS) against unfolding and aggregation but exhibited poor binding specificity toward SOD-1 in presence of blood plasma. A reasonable but not a conclusive model for the binding of these molecules was proposed on the basis of restricted docking calculations and site-directed mutagenesis of key residues at the dimer interface. A set of hydrogen bonding constraints obtained from these experiments were used to guide docking calculations with compound library around the dimer interface. A series of chemically unrelated hits were predicted, which were experimentally tested for their ability to block aggregation. At least six of the new molecules exhibited high specificity of binding toward SOD-1 in the presence of blood plasma. These molecules represent a new class of molecules for further development into clinical candidates.

Engineered disulfide bonds restore chaperone-like function of DJ-1 mutants linked to familial Parkinson's disease.

Loss-of-function mutations such as L166P, A104T, and M26I in the DJ-1 gene (PARK7) have been linked to autosomal-recessive early onset Parkinsons disease (PD). Cellular and structural studies of the familial mutants suggest that these mutations may destabilize the dimeric structure. To look for common dynamical signatures among the DJ-1 mutants, short MD simulations of up to 1000 ps were conducted to identify the weakest region of the protein (residues 38-70). In an attempt to stabilize the protein, we mutated residue Val 51 to cysteine (V51C) to make a symmetry-related disulfide bridge with the preexisting Cys 53 on the opposite subunit. We found that the introduction of this disulfide linkage stabilized the mutants A104T and M26I against thermal denaturation, improved their ability to scavenge reactive oxygen species (ROS), and restored a chaperone-like function of blocking alpha-synuclein aggregation. The L166P mutant was far too unstable to be rescued by introduction of the V51C mutation. The results presented here point to the possible development of pharmacological chaperones, which may eventually lead to PD therapeutics.

Enzymatic Characterization of ER Stress-Dependent Kinase, PERK, and Development of a High-Throughput Assay for Identification of PERK Inhibitors.

PERK is serine/threonine kinase localized to the endoplasmic reticulum (ER) membrane. PERK is activated and contributes to cell survival in response to a variety of physiological stresses that affect protein quality control in the ER, such as hypoxia, glucose depravation, increased lipid biosynthesis, and increased protein translation. Pro-survival functions of PERK are triggered by such stresses, suggesting that development of small-molecule inhibitors of PERK may be efficacious in a variety of disease scenarios. Hence, we have conducted a detailed enzymatic characterization of the PERK kinase to develop a high-throughput-screening assay (HTS) that will permit the identification of small-molecule PERK inhibitors. In addition to establishing the Km of PERK for both its primary substrate, eIF2α, and for adenosine triphosphate, further mechanistic studies revealed that PERK targets its substrate via either a random/steady-state ordered mechanism. For HTS, we developed a time-resolved fluorescence resonance energy transfer–based assay that yielded a robust Z′ factor and percent coefficient of variation value, enabling the successful screening of 79,552 compounds. This approach yielded one compound that exhibited good in vitro and cellular activity. These results demonstrate the validity of this screen and represent starting points for drug discovery efforts.

Discovery of LRRK2 inhibitors using sequential in silico joint pharmacophore space (JPS) and ensemble docking

Joint pharmacophore space (JPS), ensemble docking and sequential JPS-ensemble docking were used to select three panels of compounds (10 per panel) for evaluation as LRRK2 inhibitors. These computational methods identified four LRRK2 inhibitors with IC50 values <12μM. The sequential JPS-ensemble docking predicted the majority of active hits. One of the inhibitors (Z-8205) identified using this method was also found to inhibit the G2019S mutant of LRRK2 25-fold better than wild-type enzyme. This bias for the G2019S mutant is proposed to arise from an interaction with S2019 in this form of the enzyme. In addition, Z-8205 was found to only inhibit one other kinase when profiled against a panel of 97 kinases at 10μM.