Design of a protein kinase-inducible domain

Protein phosphorylation is a critical regulatory strategy. New tools are necessary which may be used to interrogate and are responsive to the activities of protein kinases and phosphatases. We have used protein design to develop a protein motif, termed a protein kinase-inducible domain, whose structure is dependent on its phosphorylation state. Based on an EF hand calcium-binding loop, the key design element is the replacement of a structurally critical Glu residue, which binds metal in a bidentate manner, with a serine residue, which is expected to bind metal tightly when phosphorylated but poorly when not phosphorylated. The design comprises an EF hand consensus sequence, a tryptophan at residue 7 to sensitize lanthanide luminescence, and the recognition sequence of a serine/threonine kinase. Designed peptides, which contain minimal substrate recognition motifs of the protein kinases PKA, PKC, or the MAP kinase Erk, form complexes with Tb3+ when phosphorylated, showing strong Tb3+ luminescence emission at 544 nm, but show weak luminescence when not phosphorylated. The change in fluorescence on phosphorylation is comparable to or greater than that observed in described kinase sensors. Site-specific lanthanide binding was confirmed by NMR with diamagnetic and paramagnetic metals. The kinase-inducible domain peptides comprise an expressible sequence, potentially enabling their use as genetically encoded tags of protein kinase activity. The motif is general and potentially applicable to the majority of serine/threonine kinases.     

Modeling loop backbone flexibility in receptor‐ligand docking simulations

The relevance of receptor conformational change during ligand binding is well documented for many pharmaceutically relevant receptors, but is still not fully accounted for in in silico docking methods. While there has been significant progress in treatment of receptor side chain flexibility sampling of backbone flexibility remains challenging because the conformational space expands dramatically and the scoring function must balance protein–protein and protein–ligand contributions. Here, we investigate an efficient multistage backbone reconstruction algorithm for large loop regions in the receptor and demonstrate that treatment of backbone receptor flexibility significantly improves binding mode prediction starting from apo structures and in cross docking simulations. For three different kinase receptors in which large flexible loops reconstruct upon ligand binding, we demonstrate that treatment of backbone flexibility results in accurate models of the complexes in simulations starting from the apo structure. At the example of the DFG‐motif in the p38 kinase, we also show how loop reconstruction can be used to model allosteric binding. Our approach thus paves the way to treat the complex process of receptor reconstruction upon ligand binding in docking simulations and may help to design new ligands with high specificity by exploitation of allosteric mechanisms. © 2012 Wiley Periodicals, Inc.     

Computational fragment-based design of Wee1 kinase inhibitors with tricyclic core scaffolds

Wee1 is cell cycle protein comprising a kinase domain and is a validated cancer target. We have designed molecules with variable tricyclic core scaffolds [6-6-5] system and extended them based on the chemical space available in the active site of Wee1 kinase using de novo drug design. The core scaffolds and linking fragments were extracted from pharmacophore-based virtual screening of ZINC and PubChem databases and Ludi library. These molecules bind the hinge region of kinase active site and form hydrogen bonds as confirmed from molecular docking, molecular dynamics simulations, and MM_PBSA calculations. When compared with reference inhibitors, AZD1775 and PHA-848125, the de novo designed molecules also show good docking scores and stability, retained non-covalent interactions, and high binding free energies contributed from active site residues.

     

Modified AutoDock for accurate docking of protein kinase inhibitors

Protein kinases are an important class of enzymes controlling virtually all cellular signaling pathways. Consequently, selective inhibitors of protein kinases have attracted significant interest as potential new drugs for many diseases. Computational methods, including molecular docking, have increasingly been used in the inhibitor design process [1]. We have considered several docking packages in order to strengthen our kinase inhibitor work with computational capabilities. In our experience, AutoDock offered a reasonable combination of accuracy and speed, as opposed to methods that specialize either in fast database searches or detailed and computationally intensive calculations.However, AutoDock did not perform well in cases where extensive hydrophobic contacts were involved, such as docking of SB203580 to its target protein kinase p38. Another shortcoming was a hydrogen bonding energy function, which underestimated the attraction component and, thus, did not allow for sufficiently accurate modeling of the key hydrogen bonds in the kinase-inhibitor complexes.We have modified the parameter set used to model hydrogen bonds, which increased the accuracy of AutoDock and appeared to be generally applicable to many kinase-inhibitor pairs without customization. Binding to largely hydrophobic sites, such as the active site of p38, was significantly improved by introducing a correction factor selectively affecting only carbon and hydrogen energy grids, thus, providing an effective, although approximate, treatment of solvation.     

Computational simulation of ligand docking to L-type pyruvate kinase subunit

Computational blind docking approach was used for mapping of possible binding sites in L-type pyruvate kinase subunit for peptides, RRASVA and the phosphorylated derivative RRAS(P_i)VA, which model the phosphorylatable N-terminal regulatory domain of the enzyme. In parallel, the same docking analysis was done for both substrates of this enzyme, phosphoenolpyruvate (PEP) and adenosine diphosphate (ADP), and for docking of fructose 1,6-bisphosphate (FBP), which is the allosteric activator of the enzyme. The binding properties of the entire surface of the protein were scanned and several possible binding sites were identified in domains A and C of the protein, while domain B revealed no docking sites for peptides or for substrates or the allosteric regulator. It was found that the docking sites of different ligands were partially overlapping, pointing to the possibility that some regulatory effects, observed in the case of L-type pyruvate kinase, may be caused by the competition of different ligands for the same binding sites.     

Versatile fluorescence probes of protein kinase activity

We introduce a versatile fluorescent peptide reporter of protein kinase activity. The probe can be modified to target a desired kinase by changing the kinase recognition motif in the peptide sequence. The reporter motif contains the Sox amino acid, which generates a fluorescence signal when bound to Mg2+ present in the reaction mixture. The phosphorylated peptide exhibits a much greater affinity for Mg2+ than its unphosphorylated analogue and, thus, a greater fluorescence intensity. Product formation during phosphorylation by the kinase is easily followed by the increase in fluorescence intensity over time. These probes exhibit a 3-5-fold increase in fluorescence intensity upon phosphorylation, the magnitude of which depends on the substrate. Peptides containing the reporter functionality are phosphorylated on serine by Protein Kinase C and cAMP-dependent protein kinase and are shown to be good substrates for these enzymes. The principle of this design extends to peptides phosphorylated on threonine and tyrosine.     

Studying enzyme binding specificity in acetylcholinesterase using a combined molecular dynamics and multiple docking approach

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the binding specificity of acetylcholinesterase (AChE) with its natural substrate acetylcholine (ACh), a family of substrate analogues, and choline. Calculated docking energies are well correlated to experimental k(cat)/K(M) values, as well as to experimental binding affinities of a related series of TMTFA inhibitors. The "esteratic" and "anionic" subsites are found to act together to achieve substrate binding specificity. We find that the presence of ACh in the active site of AChE not only stabilizes the setup of the catalytic triad but also tightens both subsites to achieve better binding. The docking energy gained from this induced fit is 0.7 kcal/mol for ACh. For the binding of the substrate tailgroup to the anionic subsite, both the size and the positive charge of the tailgroup are important. The removal of the positive charge leads to a weaker binding of 1 kcal/mol loss in docking energy. Substituting each tail methyl group with hydrogen results in both an incremental loss in docking energy and also a decrease in the percentage of structures docked in the active site correctly set up for catalysis.     

重组嗜热β-葡萄糖苷酶转化稀有人参皂苷Rd和CK

利用高效液相色谱(HPLC)法,对重组嗜热β-葡萄糖苷酶(Fpglu1)转化稀有人参皂苷(Rd和CK)进行研究,并表征了其催化动力学参数.利用同源模建和分子动力学模拟等生物信息学技术,探究了Fpglu1转化人参皂苷的结构基础及其相互作用.结果表明,Fpglu1能够水解人参总皂苷生成稀有皂苷Rd和CK,其催化人参皂苷Rb_1,Rb_2和Rc的K_m值分别为0.318,1.840和5.269 mmol/L;酶的转换数(k_(cat))值分别为144.191,0.572和0.011 s~(-1).当转化时间分别为6和102 h时,Rd和CK的产率达到最大,分别为60%和93%.通过对该酶的结构预测及皂苷分子的对接研究发现,底物位于由疏水性氨基酸构成的底物口袋中,氨基酸残基Glu194和Glu367是参与催化作用的关键,且实验测得的酶促反应动力学参数(K_m)与对接的相互作用能量值存在线性关系.     

RNA aptamers as pathway-specific MAP kinase inhibitors

BACKGROUND: In eukaryotic cells, many intracellular signaling pathways have closely related mitogen activated protein kinase (MAPK) paralogs as central components. Although MAPKs are therefore obvious targets to control the cellular responses resulting from the activation of these signaling pathways, the development of inhibitors which target specific cell signaling pathways involving MAPKs has proven difficult. RESULTS: We used an RNA combinatorial approach to isolate RNAs that inhibit the in vitro phosphorylation activity of extracellular regulated kinase 2 (ERK2). These inhibitors block phosphorylation by ERK1 and ERK2, but do not inhibit Jun N-terminal kinase or p38 MAPKs. Kinetic analysis indicates these inhibitors function at high picomolar concentrations through the steric exclusion of substrate and ATP binding. In one case, we identified a compact RNA structural domain responsible for inhibition. CONCLUSIONS: RNA reagents can selectively recognize and inhibit MAPKs involved in a single signal transduction pathway. The methodology described here is readily generalizable, and can be used to develop inhibitors of MAPKs involved in other signal transduction pathways. Such reagents may be valuable tools to analyze and distinguish homologous effectors which regulate distinct signaling responses.     

Proteomic identification of p38 MAP kinase substrates using in vitro phosphorylation

Protein phosphorylation is a major mechanism that regulates many basic cellular processes. Identification and characterization of substrates for a given protein kinase can lead to a better understanding of signal transduction pathways. However, it is still difficult to efficiently identify substrates for protein kinases. Here, we propose an integrated proteomic approach consisting of in vitro dephosphorylation and phosphorylation, phosphoprotein enrichment, and 2D‐DIGE. Phosphatase treatment significantly reduced the complexity of the phosphoproteome, which enabled us to efficiently identify the substrates. We employed p38 mitogen‐activated protein kinase (p38 MAP kinase) as a model kinase and identified 23 novel candidate substrates for this kinase. Seven selected candidates were phosphorylated by p38 MAP kinase in vitro and in p38 MAP kinase‐activated cells. This proteomic approach can be applied to any protein kinase, allowing global identification of novel substrates.     

Distinct patterns of activation-dependent changes in conformational mobility between ERK1 and ERK2

Hydrogen/deuterium exchange measurements by mass spectrometry (HX-MS) can be used to report localized conformational mobility within folded proteins, where exchange predominantly occurs through low energy fluctuations in structure, allowing transient solvent exposure. Changes in conformational mobility may impact protein function, even in cases where structural changes are unobservable. Previous studies of the MAP kinase, ERK2, revealed increases in HX upon activation occured at the hinge between conserved N- and C-terminal domains, which could be ascribed to enhanced backbone flexibility. This implied that kinase activation modulates interdomain closure, and was supported by evidence for two modes of nucleotide binding that were consistent with closed vs open conformations in active vs inactive forms of ERK2, respectively. Thus, phosphorylation of ERK2 releases constraints to interdomain closure, by modulating hinge flexibility. In this study, we examined ERK1, which shares 90% sequence identity with ERK2. HX-MS measurements of ERK1 showed similarities with ERK2 in overall deuteration, consistent with their similar tertiary structures. However, the patterns of HX that were altered upon activation of ERK1 differed from those in ERK2. In particular, alterations in HX at the hinge region upon activation of ERK2 did not occur in ERK1, suggesting that the two enzymes differ with respect to their regulation of hinge mobility and interdomain closure. In agreement, HX-MS measurements of nucleotide binding suggested revealed domain closure in both inactive and active forms of ERK1. We conclude that although ERK1 and ERK2 are closely related with respect to primary sequence and tertiary structure, they utilize distinct mechanisms for controlling enzyme function through interdomain interactions.     

Discovery and characterization of a cell-permeable, small-molecule c-Abl kinase activator that binds to the myristoyl binding site

c-Abl kinase activity is regulated by a unique mechanism involving the formation of an autoinhibited conformation in which the N-terminal myristoyl group binds intramolecularly to the myristoyl binding site on the kinase domain and induces the bending of the αI helix that creates a docking surface for the SH2 domain. Here, we report a small-molecule c-Abl activator, DPH, that displays potent enzymatic and cellular activity in stimulating c-Abl activation. Structural analyses indicate that DPH binds to the myristoyl binding site and prevents the formation of the bent conformation of the αI helix through steric hindrance, a mode of action distinct from the previously identified allosteric c-Abl inhibitor, GNF-2, that also binds to the myristoyl binding site. DPH represents the first cell-permeable, small-molecule tool compound for c-Abl activation.     

中药活性成分对血栓素A2受体抑制作用的分子模拟

中药中黄酮类化合物和白藜芦醇等活性成分对血栓素A2受体具有抑制作用,但具体机理不详.本研究通过同源模建方法,以墨鱼视紫红质蛋白为模板,构建血栓素A2受体的蛋白质结构模型.并使用分子对接方法研究中药活性成分白藜芦醇和芹菜苷元与血栓素A2受体模型的作用方式,据此建立药效团模型,筛选其他潜在的血栓素A2受体抑制剂.结果表明:白藜芦醇等中药活性成分能与血栓素A2受体活性口袋中的残基发生氢键作用,结合方式与血栓素相似.血栓素与Ser201、Leu198、Arg295和Thr298发生氢键作用,白藜芦醇等活性成分与Ser201、Leu198和Arg295发生氢键作用.建立的药效团模型由7个药效元素以及排斥性空间元素组成,经测试对高活性的血栓素A2受体抑制剂有比较好的选择性.使用该药效团模型对中药天然产物数据库进行筛选,命中了一批可能具有血栓素A2受体抑制作用的活性化合物.其中一些已经报道有抑制血小板凝聚活性.本研究表明血栓素A2受体可能是活血化瘀类中药的一个潜在的靶点.     

Computational Analysis Reveals Monomethylated Triazolopyrimidine as a Novel Inhibitor of SARS-CoV-2 RNA-Dependent RNA Polymerase (RdRp)

The human population is still facing appalling conditions due to several outbreaks of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) virus. The absence of specific drugs, appropriate vaccines for mutants, and knowledge of potential therapeutic agents makes this situation more difficult. Several 1, 2, 4-triazolo [1, 5-a] pyrimidine (TP)-derivative compounds were comprehensively studied for antiviral activities against RNA polymerase of HIV, HCV, and influenza viruses, and showed immense pharmacological interest. Therefore, TP-derivative compounds can be repurposed against the RNA-dependent RNA polymerase (RdRp) protein of SARS-CoV-2. In this study, a meta-analysis was performed to ensure the genomic variability and stability of the SARS-CoV-2 RdRp protein. The molecular docking of natural and synthetic TP compounds to RdRp and molecular dynamic (MD) simulations were performed to analyse the dynamic behaviour of TP compounds at the active site of the RdRp protein. TP compounds were also docked against other non-structural proteins (NSP1, NSP2, NSP3, NSP5, NSP8, NSP13, and NSP15) of SARS-CoV-2. Furthermore, the inhibition potential of TP compounds was compared with Remdesivir and Favipiravir drugs as a positive control. Additionally, TP compounds were analysed for inhibitory activity against SARS-CoV RdRp protein. This study demonstrates that TP analogues (monomethylated triazolopyrimidine and essramycin) represent potential lead molecules for designing an effective inhibitor to control viral replication. Furthermore, in vitro and in vivo studies will strengthen the use of these inhibitors as suitable drug candidates against SARS-CoV-2.     

Activation of epidermal growth factor receptor is responsible for pervanadate-induced phospholipase D activation

Pervanadate, a complex of vanadate and H(2)O(2), has an insulin mimetic effect, and acts as an inhibitor of protein tyrosine phosphatase. Pervanadate-induced phospholipase D (PLD) activation is known to be dependent on the tyrosine phosphorylation of cellular proteins and protein kinase C (PKC) activation, and yet underlying molecular mechanisms are not clearly understood. Here, we investigated the signaling pathway of pervanadate-induced PLD activation in Rat2 fibroblasts. Pervanadate increased PLD activity in dose- and time- dependent manner. Protein tyrosine kinase inhibitor, genistein, blocked PLD activation. Interestingly, AG-1478, a specific inhibitor of the tyrosine kinase activity of epidermal growth factor receptor (EGFR) blocked not only the PLD activation completely but also phosphorylation of p38 mitogen-activated protein kinase (MAPK). However, AG-1295, an inhibitor specific for the tyrosine kinase activity of pletlet drived growth factor receptor (PDGFR) did not show any effect on the PLD activation by pervanadate. We further found that pervanadate increased phosphorylation levels of p38, extracellular signal-regulated kinase (ERK) and c-Jun NH(2)-terminal kinase (JNK). SB203580, a p38 MAPK inhibitor, blocked the PLD activation completely. However, the inhibitions of ERK by the treatment of PD98059 or of JNK by the overexpression of JNK interacting peptide JBD did not show any effect on pervanadate-induced PLD activation. Inhibition or down-regulation of PKC did not alter the pervanadate-induced PLD activation in Rat2 cells. Thus, these results suggest that pervanadate-induced PLD activation is coupled to the transactivation of EGFR by pervanadate resulting in the activation of p38 MAP kinase.     

Effects of Aloe vera Flower Extract and Its Active Constituent Isoorientin on Skin Moisturization via Regulating Involucrin Expression: In Vitro and Molecular Docking Studies

Skin moisturization is very crucial for maintaining the flexibility, viscoelasticity, and differentiation of the epidermis and its deprivation causes several diseases from dry skin to dermatitis. Aloe vera, a miracle plant having diverse medicinal properties including skin moisturization effects. This study investigated for the first time the molecular mechanism targeting skin moisturization effects of the Aloe vera flower and its major active constituent. By treating human epidermal keratinocytes (HaCaT cells) with Aloe vera flower water extract (AFWE), we found that AFWE upregulated epidermal involucrin by activating the expression of protein kinase C, p38, and ERK 1/2. Additionally, it modulated filaggrin, increased aquaporin expression, and hyaluronan synthesis via a balanced regulation of HAS1 and HYAL1 protein. Similarly, it was able to protect UVB-induced photodamage. Western blot analysis, ELISA, and qRT- PCR were performed to evaluate various epidermal differentiation markers and moisturization-related factors on human epidermal keratinocytes (HaCaT cells). TLC and HPLC were used to detect and analyze the chemical constituents. Among them, we found that an active component of Aloe vera flower, isoorientin (IO) has a high binding affinity to all of its targeted proteins such as involucrin, PKC, P38, etc. through molecular docking assay. This study indicated that the Aloe vera flower and its active constituent, IO can be used as a prominent ingredient to enhance skin barrier function and improve its related pathologies.     

Phosphorylation-driven protein-protein interactions: a protein kinase sensing system

A highly flexible protein kinase sensing system is described that furnishes severalfold changes in fluorescence in response to phosphorylation. A library of Src kinase peptide substrates was prepared that contained different environmentally sensitive fluorophores positioned at various sites on the active site directed sequence. Robust changes in fluorescent intensity were observed in the presence of a phosphotyrosine binding domain protein (Lck SH2 domain), which furnishes a hydrophobic environment for the fluorophore. This protein kinase sensing system has the advantages that the fluorescent indicator can be unobtrusively positioned on the peptide substrate, and that different environmentally sensitive fluorophores with distinct photophysical properties can be employed.     

Design,Synthesis, and Biological Evaluation of Tetra-Substituted Thiophenes as Inhibitors of p38α MAPK

p38α mitogen-activated protein kinase (MAPK) plays a role in several cellular processes and consequently has been a therapeutic target in inflammatory diseases, cancer, and cardiovascular disease. A number of known p38α MAPK inhibitors contain vicinal 4-fluorophenyl/4-pyridyl rings connected to either a 5- or 6-membered heterocycle. In this study, a small library of substituted thiophene-based compounds bearing the vicinal 4-fluorophenyl/4-pyridyl rings was designed using computational docking as a visualisation tool. Compounds were synthesised and evaluated in a fluorescence polarisation binding assay. The synthesised analogues had a higher binding affinity to the active phosphorylated form of p38α MAPK than the inactive nonphosphorylated form of the protein. 4-(2-(4-fluorophenyl)thiophen-3-yl)pyridine had a K_i value of 0.6 μm to active p38α MAPK highlighting that substitution of the core ring to a thiophene retains affinity to the enzyme and can be utilised in p38α MAPK inhibitors. This compound was further elaborated using a substituted phenyl ring in order to probe the second hydrophobic pocket. Many of these analogues exhibited low micromolar affinity to active p38α MAPK. The suppression of neonatal rat fibroblast collagen synthesis was also observed suggesting that further development of these compounds may lead to potential therapeutics having cardioprotective properties.     

Molecular modeling on pyrimidine-urea inhibitors of TNF-α production: an integrated approach using a combination of molecular docking, classification techniques, and 3D-QSAR CoMSIA

Molecular docking, classification techniques, and 3D-QSAR CoMSIA were combined in a multistep framework with the ultimate goal of identifying potent pyrimidine-urea inhibitors of TNF-α production. Using the crystal structure of p38α, all the compounds were docked into the enzyme active site. The docking pose of each compound was subsequently used in a receptor-based alignment for the generation of the CoMSIA fields. "Active" and "inactive" compounds were used to build a Random Tree classification model using the docking score and the CoMSIA fields as input parameters. Domain of applicability indicated the compounds for which activity estimations can be accepted with confidence. For the active compounds, a 3D-QSAR CoMSIA model was subsequently built to accurately estimate the IC(50) values. This novel multistep framework gives insight into the structural characteristics that affect the binding and the inhibitory activity of these analogues on p38α MAP kinase, and it can be extended to other classes of small-molecule inhibitors. In addition, the simplicity of the proposed approach provides expansion to its applicability such as in virtual screening procedures.