Estimation of pH effect on the structure and stability of kinase domain of human integrin-linked kinase

Integrin-linked kinase (ILK) is an evolutionarily conserved Ser/Thr protein kinase, involved in many physiological functions such as signal transduction, actin rearrangement, cell proliferation, migration, polarisation, angiogenesis and apoptosis. An increased expression of ILK is associated with different cancers and thus considered as an attractive target for cancer therapy. We have successfully cloned, expressed and purified the kinase domain (193–446 residues) of ILK. To see the effect of pH on the structure and conformation, we performed circular diachroism, fluorescence and absorbance measurements in a wide range of pH conditions. We observed that within the range of pH 7.5–11.0, ILK_193–446 maintains its both secondary and tertiary structures. While visible aggregates were observed under the acidic pH 2.0–5.5 conditions, in order to complement these observations, we have performed molecular dynamics simulations of this kinase domain by mimicking diverse pH conditions which enabled us to see conformational preferences of the protein under such conditions. A significant correlation between the spectroscopic and molecular dynamics simulation was observed. These findings are useful to understand the conformation of ILK protein under certain pH condition which may be further implicated in the drug design and discovery.     

Molecular dynamics simulation of the substitution of serine for the conserved glycine in the G-loop in the cdc28-srm yeast mutant using the crystal lattice of human CDK2 kinase

Two-nanosecond molecular dynamics simulations of the crystal lattice of an active complex of pT160-CDK2 kinase/cyclin A/ATP-Mg²⁺/substrate were performed. The simulations showed that the structures of the wild-type CDK2 complex and the mutant CDK2 complex involving the substitution G16S-CDK2 corresponding to the yeast substitution G20S-CDC28 differ noticeably and the differences between the structural conformations are most pronounced in the regions that play a key role in the kinase functioning. The results of the computer calculations were used to consider the structural elements that may affect the kinase activity, the regulatory phosphorylation, and the binding of protein kinase with cyclins and substrates.     

Molecular dynamics simulations provide insights into the origin of gleevec’s selectivity toward human tyrosine kinases

Protein kinases are critical drug targets against cancer. Since the discovery of Gleevec, a specific inhibitor of Abl kinase, the capability of this drug to distinguish between Abl and other tyrosine kinases, such as Src, has been intensely investigated but the origin of Gleevec’s selectivity to Abl against Src is less studied. Here, we performed molecular dynamics (MD) simulations, dynamical cross-correlation matrices (DCCM), dynamical network analysis, and binding free energy calculations to explore Gleevec’s selectivity based on the crystal structures of Abl, Src, and their common ancestors (ANC-AS) and the two constructed mutation systems (AS→Abl and AS→Src). MD simulations revealed that the conformation of the phosphate-binding loop (P-loop) was altered significantly in the AS→Abl system. DCCM results unraveled that mutations increased anticorrelated motions in the AS→Abl system. Community network analysis suggested that the P-loop established special contacts in the AS→Abl system that are devoid in the AS→Src system. The binding free energy calculations unveiled that the affinity of Gleevec to AS→Abl increased to near the Abl level, whereas its affinity to AS→Src decreased to near the Src level. Analysis of individual residue contributions showed that the differences were located mainly at the P-loop. This study is valuable for understanding the sensitivity of Gleevec to human tyrosine kinases.

Communicated by Ramaswamy H. Sarma     

酪氨酸蛋白激酶和蛋白激酶C对N-乙酰氨基葡萄糖转移酶V的双重调控

在人肝癌细胞７７２１中研究了酪氨酸蛋白激酶（ＴＰＫ）和蛋白激酶Ｃ（ＰＫＣ）的激活剂［分别为表皮生长因子（ＥＧＦ）和佛波酯（ＰＭＡ）］和各种蛋白激酶抑制剂对Ｎ－乙酰氨基葡萄糖转移酶Ｖ（ＧｎＴ－Ｖ）活力的影响,以探讨ＴＰＫ和ＰＫＣ对ＧｎＴ－Ｖ的调节。结果发现,ＥＧＦ或ＰＭＡ处理细胞４８ｈ后,ＧｎＴ－Ｖ的活力明显增高;蛋白激酶的非特异性抑制剂槲皮素和染料木黄酮（ｇｅｎｉｓｔｅｉｎ）在抑制ＴＰＫ和ＰＫＣ的同时,抑制ＧｎＴ－Ｖ的基础活力,并完全阻断ＥＧＦ或ＰＭＡ对ＧｎＴ－Ｖ的增高作用;ＴＰＫ的特异性抑制剂Ｔｙｒｐｈｏｓｔｉｎ－２５和ＰＫＣ的特异性抑制剂Ｄ－鞘氨醇分别应用时,各自只能部分地取消ＥＧＦ或ＰＭＡ对ＧｎＴ－Ｖ的诱导。但当Ｔｙｒｐｈｏｓｔｉｎ－２５和Ｄ－鞘氨醇同时加入培养基中则可完全阻断ＥＧＦ或ＰＭＡ对ＧｎＴ－Ｖ的诱导激活。蛋白质合成抑制剂环己亚胺和蛋白激酶抑制剂作用相仿,不但可抑制ＧｎＴ－Ｖ的基础活力,也可完全消除ＥＧＦ或ＰＭＡ对ＧｎＴ－Ｖ的激活。以上结果提示ＥＧＦ或ＰＭＡ通过蛋白激酶调节ＧｎＴ－Ｖ的酶蛋白合成,并且ＧｎＴ－Ｖ受到膜性ＴＰＫ和ＰＫＣ的双重调节,其中ｍ－ＴＰＫ较ｍ－ＰＫＣ更为重要。     

Dynamics of allostery in hemoglobin: roles of the penultimate tyrosine H bonds

The tyrosine residues adjacent to the C termini of the hemoglobin (Hb) subunits, alphaY140 and betaY145, are expected to play important structural roles, because the C termini are the loci of T-state quaternary salt-bridges, and because the tyrosine side-chains bridge the H and F helices via H bonds to the alphaV93 and betaV98 carbonyl groups. These roles have been investigated via measurements of oxygen binding, (1)H NMR spectra, resonance Raman (RR) spectra, and time-resolved resonance Raman (TR(3)) spectra on site mutants in which the Hcdots, three dots, centeredF H bonds are eliminated by replacing the tyrosine residues with phenylalanine. The TR(3) spectra confirm the hypothesis, based on TR(3) studies of wild-type Hb, that the Hcdots, three dots, centeredF H bonds break and then re-form during the sub-microsecond phase of the R-T quaternary transition. The TR(3) spectra support the inference from other mutational studies that the alphabeta dimers act as single dynamic units in this early phase, motions of the E and F helices being coupled tightly across the dimer interface. Formation of T quaternary contacts occurs at about the same rate in the mutants as in HbA. However, these contacts are weakened substantially by the Y/F substitutions. Equilibrium perturbations are apparent also, especially for the alpha-subunits, in which relaxation of the Fe-His bond, strengthening of the Acdots, three dots, centeredE interhelical H bond, and weakening of the "switch" quaternary contact in deoxyHb are all apparent. Structural effects are less marked for the beta-chain Y/F replacement, but the Bohr effect is reduced by 25%, indicating that the salt-bridge and H bond interactions of the adjacent C terminus are loosened. The alpha-chain replacement reduces the Bohr effect much more, consistent with the global perturbations detected by the structure probes.     

PKC site mutations reveal differential modulation by insulin of NMDA receptors containing NR2A or NR2B subunits

Insulin modulates N-methyl-d-aspartate (NMDA) receptors in the CNS and potentiates currents of recombinant NMDA receptors in a subunit-specific manner in Xenopus oocytes. Previously we identified two sites in the NR2B C-terminus as targets for direct phosphorylation by C-type protein kinases (PKCs). Mutating these sites reduced insulin potentiation of currents by one half, reflecting the PKC-mediated portion of the NR2B insulin effect. The PKC-proline rich tyrosine kinase (Pyk2)-Src family kinase pathway may also mediate insulin potentiation. A dominant negative Pyk2 mutant significantly reduced insulin potentiation when co-expressed with NR2B-containing receptors, suggesting that Pyk2 and downstream Src-family tyrosine kinases are involved, along with PKCs, in insulin potentiation of NR2B. The NR2A C-terminus contains two residues homologous to the NR2B PKC targets. Mutating both these sites eliminated insulin potentiation of NR2A-containing receptors, while co-expression of dominant negative Pyk2 had no effect. Together, these data indicate that PKCs alone mediate the NR2A insulin effect. When tested individually for importance in insulin potentiation, the two PKC sites showed an additive effect in potentiation of NR2A-containing receptors. Insulin modulation of NR2A-containing receptors is mediated solely by PKCs, whereas insulin modulation of NR2B-containing receptors is mediated by PKCs and tyrosine kinases (PTKs).     

Multicomplex-based pharmacophore modeling in conjunction with multi-target docking and molecular dynamics simulations for the identification of PfDHFR inhibitors

Abstract

Plasmodium falciparum dihydrofolate reductase enzyme (PfDHFR) is counted as one of the attractive and validated antimalarial drug targets. However, the point mutations in the active site of wild-type PfDHFR have developed resistance against the well-known antifolates. Therefore, there is a dire need for the development of inhibitors that can inhibit both wild-type and mutant-type DHFR enzyme. In the present contribution, we have constructed the common feature pharmacophore models from the available PfDHFR. A representative hypothesis was prioritized and then employed for the screening of natural product library to search for the molecules with complementary features responsible for the inhibition. The screened candidates were processed via drug-likeness filters and molecular docking studies. The docking was carried out on the wild-type PfDHFR (3QGT); double-mutant PfDHFR (3UM5 and 1J3J) and quadruple-mutant PfDHFR (1J3K) enzymes. A total of eight common hits were obtained from the docking calculations that could be the potential inhibitors for both wild and mutant type DHFR enzymes. Eventually, the stability of these candidates with the selected proteins was evaluated via molecular dynamics simulations. Except for SPECS14, all the prioritized candidates were found to be stable throughout the simulation run. Overall, the strategy employed in the present work resulted in the retrieval of seven candidates that may show inhibitory activity against PfDHFR and could be further exploited as a scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Dynamics and structural stability effects of germline PTEN mutations associated with cancer versus autism phenotypes

Individuals with germline mutations in the tumor suppressor gene phosphatase and tensin homolog (PTEN), irrespective of clinical presentation, are diagnosed with PTEN hamartoma tumor syndrome (PHTS). PHTS confers a high risk of breast, thyroid, and other cancers or autism spectrum disorder (ASD) with macrocephaly. It remains unclear why mutations in one gene can lead to seemingly disparate phenotypes. Thus, we sought to identify differences in ASD vs. cancer-associated germline PTEN missense mutations by investigating putative structural effects induced by each mutation. We utilized a theoretical computational approach combining in silico structural analysis and molecular dynamics (MD) to interrogate 17 selected mutations from our patient population: six mutations were observed in patients with ASD (only), six mutations in patients with PHTS-associated cancer (only), four mutations shared across both phenotypes, and one mutation with both ASD and cancer. We demonstrate structural stability changes where all six cancer-associated mutations showed a global decrease in structural stability and increased dynamics across the domain interface with a proclivity to unfold, mediating a closed (inactive) active site. In contrast, five of the six ASD-associated mutations showed localized destabilization that contribute to the partial opening of the active site. Our results lend insight into distinctive structural effects of germline PTEN mutations associated with PTEN-ASD vs. those associated with PTEN-cancer, potentially aiding in identification of the shared and separate molecular features that contribute to autism or cancer, thus, providing a deeper understanding of genotype–phenotype relationships for germline PTEN mutations.