Structure-guided optimization of small molecule c-Abl activators

c-Abl kinase is maintained in its normal inactive state in the cell through an assembled, compact conformation. We describe two chemical series that bind to the myristoyl site of the c-Abl kinase domain and stimulate c-Abl activation. We hypothesize that these molecules activate c-Abl either by blocking the C-terminal helix from adopting a bent conformation that is critical for the formation of the autoinhibited conformation or by simply providing no stabilizing interactions to the bent conformation of this helix. Structure-based molecular modeling guided the optimization of binding and activation of c-Abl of these two chemical series and led to the discovery of c-Abl activators with nanomolar potency. The small molecule c-Abl activators reported herein could be used as molecular tools to investigate the biological functions of c-Abl and therapeutic implications of its activation.     

Molecular dynamics simulation study reveals polar nature of pathogenic mutations responsible for stabilizing active conformation of kinase domain in leucine-rich repeat kinase II

The kinase domain of LRRK2 is increasingly gaining attention as a promising therapeutic target due to pathogenic mutation leading to development of Parkinson’s disease. Mutation in G2019S and I2020T increases the kinase activity, while A2016T mutation causes drug resistance. Increased kinase activity of LRRK2 has been associated with deposition of tau and α-synuclein proteins. However, mechanism responsible for increase in activity due to mutation is not known. In the present study, extensive molecular dynamics study has been performed on both wild and mutant homology models of DYG-In (active) conformation of the kinase domain of LRRK2 in the absence/presence of ATP at the active site to study the behavior of DYG loop. In absence of ATP, it is observed that G2019S and I2020T mutants stabilize DYG loop by increasing formation of hydrogen bond with neighboring residues, mainly with GLU 1920 and ILE 1991, respectively. In ATP-kinase complex, DYG loop also increases hydrogen bonding with neighboring residues in mutant LRRK2. The study indicates that polar side chain of mutated residues increases the polarity of DYG loop, causing an increase in hydrogen bonding with neighboring residues to stabilize the active conformation of kinase domain in LRRK2. The binding free energy of ATP is found to be higher in mutated kinase as compared to wild, due to more stable kinase domain.     

Disruption of protein-membrane binding and identification of small-molecule inhibitors of coagulation factor VIII

Factor VIII is a critical member of the blood coagulation cascade. It binds to the membrane surfaces of activated platelets at the site of vascular injury via a highly specific interaction between factor VIII's carboxy-terminal C2 domain and their phosphatidylserine-rich lipid bilayer. We have identified small-molecule inhibitors of factor VIII's membrane binding activity that have IC50 values as low as 2.5 microM. This interaction is approximately 10(3)-fold tighter than that of free o-phospho-L-serine. These compounds also inhibit factor VIII-dependent activation of factor X, indicating that disruption of membrane lipid binding leads to inhibition of the intrinsic coagulation pathway. The tightest binding inhibitor is specific and does not prevent membrane binding by the closely related coagulation factor V. These results indicate that this and related compounds may be used as leads to develop novel antithrombotic agents.     

Allosteric inhibition abrogates dysregulated LFA-1 activation: Structural insight into mechanisms of diminished immunologic disease

Lymphocyte Function Associated antigen-1(LFA-1) has been implicated severely in the pathophysiology of inflammatory and autoimmune diseases. Its active and inactive conformations correlate with its diseased and non-diseased state respectively. This is determined by its degree of affinity for its intrinsic ligand (ICAM) at the active site and accompanying synergistic coordination at the α7 helix. This potentiates the role of inhibitors in disrupting this interaction allosterically. Herein, we present a first account of the structural dynamics which characterizes the inhibitory effect of a novel LFA-1 antagonist, Lifitegrast (SAR1118), upon binding to the I-domain allosteric site (IDAS) using molecular dynamics simulation. Findings from this study revealed that the inhibitor stabilized the closed conformation and reversed the open conformation to a low ICAM-affinity state (closed) as evidenced by the upward movement of the α7 helix and corresponding transitions at the active site. This in both cases favors the formation of the non-disease inactive form. Upon allosteric modulation, the inhibitor significantly restored protein stability, enhanced compactness and decreased residual fluctuation as crucial to its potency in the amelioration of immunological and inflammatory diseases which agrees with experimental studies. These findings could therefore serve as the basis for the exploration of the allosteric domain and its active site affinity modulation to aid the design of more specific and selective inhibitors.     

Metal binding discrimination of the calmodulin Q41C/K75C mutant on Ca2+ and La3+

Calmodulin (CaM) is a multifunctional Ca²⁺-binding protein regulating the activity of many enzymes in response to fluctuation of the intracellular Ca²⁺ level. It has been shown that a CaM Q41C/K75C mutant (CaMSS) with a disulfide bond in the N-terminal domain exhibits greatly reduced affinity to Ca²⁺. In the present study, the experimental results revealed a unique metal binding pattern in CaMSS towards La³⁺ and Ca²⁺ separately: the mutant protein binds Ca²⁺ at site I, III and IV; however, it binds La³⁺ at site I, II and IV. A putative mechanism was proposed which is the conformation of site II (or site III) of CaMSS could be altered and thus loses its metal ion affinity in response to metal binding in the opposite terminal domain possibly through the long range domain interaction. The present work may offer new perspectives for understanding the mechanisms of specific metal ion affinity in CaM and for CaM-based protein design.     

Membrane-dependent effects of a cytoplasmic helix on the structure and drug binding of the influenza virus M2 protein

The influenza A M2 protein forms a proton channel for virus infection and also mediates virus assembly and budding. The minimum protein length that encodes both functions contains the transmembrane (TM) domain (roughly residues 22-46) for the amantadine-sensitive proton-channel activity and an amphipathic cytoplasmic helix (roughly residues 45-62) for curvature induction and virus budding. However, structural studies involving the TM domain with or without the amphipathic helix differed on the drug-binding site. Here we use solid-state NMR spectroscopy to determine the amantadine binding site in the cytoplasmic-helix-containing M2(21-61). (13)C-(2)H distance measurements of (13)C-labeled protein and (2)H-labeled amantadine showed that in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers, the first equivalent of drug bound S31 inside the M2(21-61) pore, similar to the behavior of M2 transmembrane peptide (M2TM) in DMPC bilayers. The nonspecific surface site of D44 observed in M2TM is disfavored in the longer peptide. Thus, the pharmacologically relevant drug-binding site in the fully functional M2(21-61) is S31 in the TM pore. Interestingly, when M2(21-61) was reconstituted into a virus-mimetic membrane containing 30% cholesterol, no chemical shift perturbation was observed for pore-lining residues, whereas M2TM in the same membrane exhibited drug-induced chemical shift changes. Reduction of the cholesterol level and the use of unsaturated phospholipids shifted the conformational equilibrium of M2TM fully to the bound state but did not rescue drug binding to M2(21-61). These results suggest that the amphipathic helix, together with cholesterol, modulates the ability of the TM helix to bind amantadine. Thus, the M2 protein interacts with the lipid membrane and small-molecule inhibitors in a complex fashion, and a careful examination of the environmental dependence of the protein conformation is required to fully understand the structure-function relation of this protein.     

Dynamic Equilibrium of the Aurora A Kinase Activation Loop Revealed by Single‐Molecule Spectroscopy

The conformation of the activation loop (T‐loop) of protein kinases underlies enzymatic activity and influences the binding of small‐molecule inhibitors. By using single‐molecule fluorescence spectroscopy, we have determined that phosphorylated Aurora A kinase is in dynamic equilibrium between a DFG‐in‐like active T‐loop conformation and a DFG‐out‐like inactive conformation, and have measured the rate constants of interconversion. Addition of the Aurora A activating protein TPX2 shifts the equilibrium towards an active T‐loop conformation whereas addition of the inhibitors MLN8054 and CD532 favors an inactive T‐loop. We show that Aurora A binds TPX2 and MLN8054 simultaneously and provide a new model for kinase conformational behavior. Our approach will enable conformation‐specific effects to be integrated into inhibitor discovery across the kinome, and we outline some immediate consequences for structure‐based drug discovery.     

The Molecular Dynamic Simulation of Zolpidem Interaction with Gamma Aminobutyric Acid Type A Receptor

Our goal was to generate the extracellular domain of gamma‐aminobutyric acid type A receptor (GABA_A receptor) by comparative modeling and to study the interaction of zolpidem with the GABA_A receptor. The modeling strategy was verified to provide reasonable 3‐dimensional coordinates. These coordinates helped to combine all the subunits well. The benzodiazepine (BZ) binding site was located in a binding pocket between the α₁ and γ₂ subunits of the GABA_A receptor. Zolpidem was selected to dock into the binding site. In our study, the residues of the binding pocket were suggested to be αHis129, αTyr187, αGly228, αThr234, αTyr237, γMet96, γPhe116, γSer130, γGly143, and γMet169. By the calculation of the docking module, the conformation of zolpidem docking in the BZ binding site was investigated. A hydrogen bond was found at γArg136 when zolpidem's conformation was in rank 2 of the docking score. The contracted binding pocket showed residues at αHis129, αTyr187, αGly228, αTyr237, γPhe116, and γMet169. Zolpidem docking in a contracted binding pocket might generate a hydrogen bond in α His 129.     

Double Histidine Based EPR Measurements at Physiological Temperatures Permit Site-Specific Elucidation of Hidden Dynamics in Enzymes

Protein dynamics is at the heart of all cellular processes. Here, we utilize the dHis-Cu^IINTA label to obtain site-specific information on dynamics for both an α-helix and β-sheet site of GB1, the immunoglobulin binding domain of protein G. Spectral features found in our CW-EPR measurements were consistent with the overall rigid nature of GB1 and with predictions from molecular dynamics simulations. Using this information, we show the potential of this approach to elucidate the role of dynamics in substrate binding of a functionally necessary α-helix in human glutathione transferase A1-1 (hGSTA1-1). We observe two dynamical modes for the helix. The addition of the inhibitor GS-Met and GS-Hex resulted in hGSTA1-1 to favor the more rigid active state conformation, while the faster mode potentially aids the search for substrates. Together the results illustrate the remarkable potential of the dHis-based labelling approach to measure site-specific dynamics using room temperature lineshape analysis.     

Structural rationale for the cross-resistance of tumor cells bearing the A399V variant of elongation factor eEF1A1 to the structurally unrelated didemnin B,ternatin, nannocystin A and ansatrienin B

At least four classes of structurally distinct natural products with potent antiproliferative activities target the translation elongation factor eEF1A1, which is best known as the G-protein that delivers amino acyl transfer RNAs (aa-tRNAs) to ribosomes during mRNA translation. We present molecular models in atomic detail that provide a common structural basis for the high-affinity binding of didemnin B, ternatin, ansatrienin B and nannocystin A to eEF1A1, as well as a rationale based on molecular dynamics results that accounts for the deleterious effect of replacing alanine 399 with valine. The proposed binding site, at the interface between domains I and III, is eminently hydrophobic and exists only in the GTP-bound conformation. Drug binding at this site is expected to disrupt neither loading of aa-tRNAs nor GTP hydrolysis but would give rise to stabilization of this particular conformational state, in consonance with reported experimental findings. The experimental solution of the three-dimensional structure of mammalian eEF1A1 has proved elusive so far and the highly homologous eEF1A2 from rabbit muscle has been crystallized and solved only as a homodimer in a GDP-bound conformation. Interestingly, in this dimeric structure the large interdomain cavity where the drugs studied are proposed to bind is occupied by a mostly hydrophobic α-helix from domain I of the same monomer. Since binding of this α-helix and any of these drugs to domain III of eEF1A(1/2) is, therefore, mutually exclusive and involves two distinct protein conformations, one intriguing possibility that emerges from our study is that the potent antiproliferative effect of these natural products may arise not only from inhibition of protein synthesis, which is the current dogma, but also from interference with some other non-canonical functions. From this standpoint, this type of drugs could be considered antagonists of eEF1A1/2 oligomerization, a hypothesis that opens up novel areas of research.     

Molecular dynamics simulations reveal the allosteric effect of F1174C resistance mutation to ceritinib in ALK-associated lung cancer

Anaplastic lymphoma kinase (ALK) has become as an important target for the treatment of various human cancers, especially non-small-cell lung cancer. A mutation, F1174C, suited in the C-terminal helix αC of ALK and distal from the small-molecule inhibitor ceritinib bound to the ATP-binding site, causes the emergence of drug resistance to ceritinib. However, the detailed mechanism for the allosteric effect of F1174C resistance mutation to ceritinib remains unclear. Here, molecular dynamics (MD) simulations and binding free energy calculations [Molecular Mechanics/Generalized Born Surface Area (MM/GBSA)] were carried out to explore the advent of drug resistance mutation in ALK. MD simulations observed that the exquisite aromatic-aromatic network formed by residues F1098, F1174, F1245, and F1271 in the wild-type ALK-ceritinib complex was disrupted by the F1174C mutation. The resulting mutation allosterically affected the conformational dynamic of P-loop and caused the upward movement of the P-loop from the ATP-binding site, thereby weakening the interaction between ceritinib and the P-loop. The subsequent MM/GBSA binding free energy calculations and decomposition analysis of binding free energy validated this prediction. This study provides mechanistic insight into the allosteric effect of F1174C resistance mutation to ceritinib in ALK and is expected to contribute to design the next-generation of ALK inhibitors.     

Insights into the conformational switching mechanism of the human vascular endothelial growth factor receptor type 2 kinase domain

Human vascular endothelial growth factor receptor type 2 (h-VEFGR2) is a receptor tyrosine kinase involved in the angiogenesis process and regarded as an interesting target for the design of anticancer drugs. Its activation/inactivation mechanism is related to conformational changes in its cytoplasmatic kinase domain, involving first among all the αC-helix in N-lobe and the A-loop in C-lobe. Affinity of inhibitors for the active or inactive kinase form could dictate the open or closed conformation of the A-loop, thus making the different conformations of the kinase domain receptor (KDR) domain different drug targets in drug discovery. In this view, a detailed knowledge of the conformational landscape of KDR domain is of central relevance to rationalize the efficiency and selectivity of kinase inhibitors. Here, molecular dynamics simulations were used to gain insight into the conformational switching activity of the KDR domain and to identify intermediate conformations between the two limiting active and inactive conformations. Specific energy barriers have been selectively removed to induce, and hence highlight at the atomistic level, the regulation mechanism of the A-loop opening. The proposed strategy allowed to repeatedly observe the escape of the KDR domain from the DFG-out free energy basin and to identify rare intermediate conformations between the DFG-out and the DFG-in structures to be employed in a structure-based drug discovery process.     

Solvent- and pH-induced self-assembly of cationic meta-linked poly(phenylene ethynylene): effects of helix formation on amplified fluorescence quenching and Förster resonance energy transfer

We reported here the synthesis and characterization of a novel water-soluble, meta-linked poly(phenylene ethynylene) (m-PPE-NEt(2)Me(+)) featuring quaternized side groups. We studied the solvent-induced self-assembly of m-PPE-NEt(2)Me(+) in MeOH/H(2)O solvent mixtures by using UV-vis absorption and fluorescence spectroscopies. The results showed that the polymer folded into a helical conformation and that the extent of helical folding increased with the volume % water in the solvent. This cationic polymer also exhibited unique pH-induced helix formation, which was attributed to the partial neutralization of quaternized side groups at high pH and the meta-links in the main chain of the polymer. Studies on the fluorescence quenching of m-PPE-NEt(2)Me(+) by anthraquinone-2,6-disulfonate (AQS) and Fe(CN)(6)(4-), two small-molecule anionic quenchers with different typical structures, revealed more efficient quenching of helical conformation by AQS than by Fe(CN)(6)(4-). We proposed that the two quenchers most likely interacted with the polymer helix in two different modes; that was, AQS featuring large planar aromatic ring could intercalate within adjacent π-stacked phenylene ethynylene units in the polymer helix, whereas Fe(CN)(6)(4-) mainly bound to the periphery of polymer helix through ion-pair formation. Finally, the results of FRET from the helical polymer to the fluorescein (C*)-labeled polyanions, ssDNA-C* (ssDNA: single-stranded DNA) and dsDNA-C* (dsDNA: double-stranded DNA) also suggested two different modes of interactions. As compared with the FRET to dsDNA-C*, the FRET to ssDNA-C* was slightly more efficient, which was believed to arise from the additional binding of ssDNA-C* with the polymer via intercalation of its exposed hydrophobic bases into the π stack of adjacent phenylene ethynylene units in the polymer helix.     

Using nitrile-derivatized amino acids as infrared probes of local environment

It is well-known that the C=N stretching vibration in acetonitrile is sensitive to solvent. Therefore, we proposed in this contribution to use this vibrational mode to report local environment of a particular amino acid in proteins or local environmental changes upon binding or folding. We have studied the solvent-induced frequency shift of two nitrile-derivatized amino acids, which are, AlaCN and PheCN, in H(2)O and tetrahydrofuran (THF), respectively. Here, THF was used to approximate a protein's hydrophobic interior because of its low dielectric constant. As expected, the C=N stretching vibrations of both AlaCN and PheCN shift as much as approximately 10 cm(-1) toward higher frequency when THF was replaced with H2O, indicative of the sensitivity of this vibration to solvation. To further test the utility of nitrile-derivatized amino acids as probes of the environment within a peptide, we have studied the binding between calmodulin (CaM) and a peptide from the CaM binding domain of skeletal muscle myosin light chain kinase (MLCK(579-595)), which contains a single PheCN. MLCK(579-595) binds to CaM in a helical conformation. When the PheCN was substituted on the polar side of the helix, which was partially exposed to water, the C=N stretching vibration is similar to that of PheCN in water. In constrast, when PheCN is introduced at a site that becomes buried in the interior of the protein, the C=N stretch is similar to that of PheCN in THF. Together, these results suggest that the C=N stretching vibration of nitrile-derivatized amino acids can indeed be used as local internal environmental markers, especially for protein conformational studies.     

The different flexibility of c-Src and c-Abl kinases regulates the accessibility of a druggable inactive conformation

c-Src and c-Abl are two closely related protein kinases that constitute important anticancer targets. Despite their high sequence identity, they show different sensitivities to the anticancer drug imatinib, which binds specifically to a particular inactive conformation in which the Asp of the conserved DFG motif points outward (DFG-out). We have analyzed the DFG conformational transition of the two kinases using massive molecular dynamics simulations, free energy calculations, and isothermal titration calorimetry. On the basis of the reconstruction of the free energy surfaces for the DFG-in to DFG-out conformational changes of c-Src and c-Abl, we propose that the different flexibility of the two kinases results in a different stability of the DFG-out conformation and might be the main determinant of imatinib selectivity.     

Primary processes during the light-signal transduction of phototropin

Phototropin is a blue-light photoreceptor in plants that mediates phototropism, chloroplast relocation, stomata opening and leaf expansion. Phototropin molecule has two photoreceptive domains named LOV1 (light-oxygen-voltage) and LOV2 in the N-terminus and a serine/threonine kinase domain in the C-terminus, and acts as a blue light-regulated kinase. Each LOV domain binds a flavin mononucleotide as a chromophore and undergoes unique cyclic reactions upon blue-light absorption that comprises a cysteinyl-flavin adduct formation through a triplet-excited state and a successive adduct break to revert to the initial ground state. The molecular reactions underlying the photocycle are reviewed and one of the probable molecular schemes is presented. Adduct formation alters the secondary protein structure of the LOV domains. This structural change could be transferred to the linker between the kinase domain and involved in the photoregulation of the kinase activity. The structural changes as well as the oligomeric structures seem to differ between LOV1 and LOV2, which may explain the proposed roles of each domain in the photoregulation of the kinase activity. The photoregulation mechanism of phototropin kinase is reviewed and discussed in reference to the regulation mechanism of protein kinase A, which it resembles.     

Rewiring kinase specificity with a synthetic adaptor protein

Signaling cascades are managed in time and space by interactions between and among proteins. These interactions are often aided by adaptor proteins, which guide enzyme-substrate pairs into proximity. Miniature proteins are a class of small, well-folded protein domains possessing engineered binding properties. Here we made use of two miniature proteins with complementary binding properties to create a synthetic adaptor protein that effectively redirects a ubiquitous signaling event: tyrosine phosphorylation. We report that miniature-protein-based adaptor 3 uses templated catalysis to redirect the Src family kinase Hck to phosphorylate hDM2, a negative regulator of the p53 tumor suppressor and a poor Hck substrate. Phosphorylation occurs with multiple turnover and at a single site targeted by c-Abl kinase in the cell.     

The Mechanism of p53 Rescue by SUSP4

p53 is an important tumor‐suppressor protein deactivation of which by mdm2 results in cancers. A SUMO‐specific protease 4 (SUSP4) was shown to rescue p53 from mdm2‐mediated deactivation, but the mechanism is unknown. The discovery by NMR spectroscopy of a “p53 rescue motif” in SUSP4 that disrupts p53‐mdm2 binding is presented. This 29‐residue motif is pre‐populated with two transient helices connected by a hydrophobic linker. The helix at the C‐terminus binds to the well‐known p53‐binding pocket in mdm2 whereas the N‐terminal helix serves as an affinity enhancer. The hydrophobic linker binds to a previously unidentified hydrophobic crevice in mdm2. Overall, SUSP4 appears to use two synergizing modules, the p53 rescue motif described here and a globular‐structured SUMO‐binding catalytic domain, to stabilize p53. A p53 rescue motif peptide exhibits an anti‐tumor activity in cancer cell lines expressing wild‐type p53. A pre‐structures motif in the intrinsically disordered proteins is thus important for target recognition.