Structural and thermodynamic characterization of endo-1,3-β-glucanase: Insights into the substrate recognition mechanism

Endo-1,3-β-glucanase from Cellulosimicrobium cellulans is composed of a catalytic domain and a carbohydrate-binding module. We have determined the X-ray crystal structure of the catalytic domain at a high resolution of 1.66 Å. The overall fold is a sandwich-like β-jelly roll architecture like the enzymes in the glycoside hydrolase family 16. The substrate-binding cleft has a length and a width of ~ 28 and ~ 15 Å, respectively, which is thought to be capable of accommodating at least six glucopyranose units. Laminarihexaose was placed into the substrate-binding cleft, namely at the subsites + 2 to − 4 from the reducing end, and the complex structure was analyzed using molecular dynamics simulations (MD) and using a rotamer search of the pocket. During the MD simulations, the substrate fluctuated more than the enzyme, where the residues at the subsites toward the non-reducing end fluctuated more than those toward the reducing end. Little conformational change of the protein was observed for the subsites + 1 and + 2, indicating that the glucose's position could be tightly restricted inside the pocket. Substrate binding experiments using isothermal titration calorimetry showed that the binding affinity of laminaritriose was higher than that of laminaribiose and similar to those of other longer laminarioligosaccharides. Taken together, the substrates mainly bind to the subsites − 1 to − 3 with the highest affinity, while the part bound to the reducing end would be hydrolyzed.     

Structure-guided studies of the SHP-1/JAK1 interaction provide new insights into phosphatase catalytic domain substrate recognition

SHP-1 (PTPN6) is a member of the SHP sub-family of protein tyrosine phosphatases and plays a critical role in the regulation of the JAK/STAT signaling pathway. Previous studies suggested that SHP-1 contains a PTP1B-like second phosphotyrosine pocket that allows for binding of tandem phosphotyrosine residues, such as those found in the activation loop of JAK kinases. To discover the structural nature of the interaction between SHP-1 and the JAK family member, JAK1, we determined the 1.8 Å co-crystal structure of the SHP-1 catalytic domain and a JAK1-derived substrate peptide. This structure reveals electron density for only one bound phosphotyrosine residue. To investigate the role of the predicted second site pocket we determined the structures of SHP-1 in complex with phosphate and sulfate to 1.37 Å and 1.7 Å, respectively, and performed anomalous scattering experiments for a selenate-soaked crystal. These crystallographic data suggest that SHP-1 does not contain a PTP1B-like second site pocket. This conclusion is further supported by analysis of the relative dephosphorylation and binding affinities of mono- and tandem-phosphorylated peptide substrates. The crystal structures instead indicate that SHP-1 contains an extended C-terminal helix α2’ incompatible with the predicted second phosphotyrosine binding site. This study suggests that SHP-1 defines a new category of PTP1B-like protein tyrosine phosphatases with a hindered second phosphotyrosine pocket.     

Probing of CD4 binding pocket of HIV-1 gp120 glycoprotein using unnatural phenylalanine analogues

CD4-gp120 interaction is the first step for HIV-1 entry into host cells. A highly conserved pocket in gp120 protein is an attractive target for developing gp120 inhibitors or novel HIV detection tools. Here we incorporate seven phenylalanine derivatives having different sizes and steric conformations into position 43 of domain 1 of CD4 (mD1.2) to explore the architecture of the ‘Phe43 cavity’ of HIV-1 gp120. The results show that the conserved hydrophobic pocket in gp120 tolerates a hydrophobic side chain of residue 43 of CD protein, which is 12.2 Å in length and 8.0 Å in width. This result provides useful information for developing novel gp120 inhibitors or new HIV detection tools.     

Characterization of a highly efficient heterodimeric xylosidase from Humicola insolens

A heterodimeric xylosidase (E.C. 3.2.1.37) with robust activity is secreted among the plant cell wall degrading enzymes produced by the saprophytic fungus Humicola insolens. The xylosidase has been purified to homogeneity by gel filtration and cation exchange chromatography, and demonstrated to be composed of two protein subunits of 68 and 17 kDa with a molecular mass in solution of approximately 85 kDa based on a combination of SDS-PAGE, size exclusion chromatography and analytical ultracentrifugation. Peptide sequence identities from the subunits indicate the 68 kDa subunit contains a catalytic protein domain and the 17 kDa subunit a carbohydrate binding module. The xylosidase has wide biotechnological potential with maximum activity exhibited at 70 °C and kinetic constants with p-nitrophenol xylopyranoside substrate that suggest it has the highest catalytic efficiency recorded to date (Vmax 22.17 μmoles/min/mg, Km 1.74 mM and Kcat 6787/s).     

Ala258Phe substitution in Bacillus sp. YX-1 glucose dehydrogenase improves its substrate preference for xylose

Bacillus sp. YX-1 glucose dehydrogenase (BsGDH) with good solvent resistance catalyzes the oxidation of β-d-glucose to d-glucono-1,5-lactone. Xylose is a recyclable resource from hemicellulase hydrolysis. In this work, to improve the preference of BsGDH for xylose, we designed seven mutants inside or adjacent to the substrate binding pocket using site-directed mutagenesis. Among all mutants, Ala258Phe mutant displayed the highest activity of 7.59 U mg⁻¹ and nearly 8-folds higher k_cat/K_m value towards xylose than wild-type BsGDH. The kinetic constants indicated that the A258F mutation effectively altered the transition state. By analysis of modeled protein structure, Ala258Phe created a space to facilitate the reactivity towards xylose. A258F mutant retained good solvent resistance in glycol, ethyl caprylate, octane, decane, cyclohexane, nonane, etc. as with BsGDH. This work provides a protein engineering approach to modify the substrate stereo-preference of alcohol dehydrogenase and a promising enzyme for cofactor regeneration in chiral catalysis.     

Autophosphorylation Is a Mechanism of Inhibition in Twitchin Kinase

Titin-like kinases are muscle-specific kinases that regulate mechanical sensing in the sarcomere. Twitchin kinase (TwcK) is the best-characterized member of this family, both structurally and enzymatically. TwcK activity is auto-inhibited by a dual intrasteric mechanism, in which N- and C-terminal tail extensions wrap around the kinase domain, blocking the hinge region, the ATP binding pocket and the peptide substrate binding groove. Physiologically, kinase activation is thought to occur by a stretch-induced displacement of the inhibitory tails from the kinase domain. Here, we now show that TwcK inhibits its catalysis even in the absence of regulatory tails, by undergoing auto-phosphorylation at mechanistically important elements of the kinase fold. Using mass spectrometry, site-directed mutagenesis and catalytic assays on recombinant samples, we identify residues T212, T301, T316 and T401 as primary auto-phosphorylation sites in TwcK in vitro. Taken together, our results suggest that residue T316, located in the peptide substrate binding P + 1 loop, is the dominantly regulatory site in TwcK. Based on these findings, we conclude that TwcK is regulated through a triple-inhibitory mechanism consisting of phosphorylation and intrasteric blockage, which is responsive not only to mechanical cues but also to biochemical modulation. This implies that mechanically stretched conformations of TwcK do not necessarily correspond to catalytically active states, as previously postulated. This further suggests a phosphorylation-dependent desensitization of the TwcK-mediated mechanoresponse of the sarcomere in vivo.     

Computational exploration of maternal embryonic leucine zipper kinase (MELK) as a cancer drug target

Maternal embryonic leucine zipper kinase (MELK) is of vital importance due to its significant role in cancer development and its association with poor prognosis in different cancers. Here, we employed several computer aided drug design approaches to shortlist potential binding molecules of MELK. For virtual screening, asinex oncology library (containing 6334 drugs) and comprehensive marine natural products database (containing approximately 32,000 drugs) were used. The study identified two drug molecules: Top-2 and Top-3 as high affinity binding MELK molecules compared to the control co-crystalized Top-1 inhibitor. Both the shortlisted compounds and the control showed high stable binding free energy and high GOLD score. The compounds and control also reported stable dynamics with root mean square deviations (RMSD) value ～ 2 Å in 500 ns. Similarly, the MELK active site residues were observed in good stability with the compounds. Further, it was noticed the compounds/control formed multiple hydrogen bonds with the MELK active pocket residues which is the main reason of high intermolecular stability. Atomic level binding free energies determined van der Waals and electrostatic energies to play vital role in stable complex formation. From drug likeness and pharmacokinetics perspective, the compounds are ideal molecules for further investigation. Overall, the results are promising and might be tested in in vivo and in vitro studies against MELK.     

Purification and characterization of a flavin reductase from the biodesulfurizing bacterium Mycobacterium goodii X7B

Dibenzothiophene (DBT) in fossil fuels can be efficiently biodesulfurized by a thermophilic bacterium Mycobacterium goodii X7B. Flavin reductase DszD, which catalyzes the reduction of oxidated flavin by NAD(P)H, is indispensable for the biodesulfurization process. In this work, a flavin reductase DszD in M. goodii X7B was purified to homogeneity, and then its encoding gene dszD was amplified and expressed in Escherichia coli. DszD is a homodimer with each subunit binding one FMN as cofactor. The K_m values for FMN and NADH of the purified recombinant DszD were determined to be 6.6 ± 0.3 μM and 77.9 ± 5.4 μM, respectively. The optimal temperature for DszD activity was 55 °C. DszD can use FMN or FAD as substrate to generate FMNH₂ or FADH₂ as product. DszD was coexpressed with DBT monooxygenase DszC, the enzyme catalyzing the first step of the biodesulfurization process. It was indicated that the coexpressed DszD could effectively enhance the DszC catalyzed DBT desulfurization reaction.     

Discovery of ZAP70 inhibitors by high-throughput docking into a conformation of its kinase domain generated by molecular dynamics

Very few selective inhibitors of the zeta-chain associated protein kinase 70 kDa (ZAP70) have been reported despite its importance in autoimmune diseases. Here, to induce a fit of the so-called gatekeeper residue (Met414) and hydrophobic pocket next to it, a potent Janus kinase 2 (JAK2) inhibitor was first docked into the ATP binding site of ZAP70 by structural alignment of the kinase domains. The resulting model of the complex between ZAP70 and the JAK2 inhibitor was then relaxed by an explicit solvent molecular dynamics simulation with restraints on the backbone atoms. High-throughput docking into the induced-fit conformation of ZAP70 generated by molecular dynamics has revealed 10 low-micromolar inhibitors which correspond to six distinct chemotypes. One of these ZAP70 inhibitors has an IC₅₀ of 110 nM for JAK2.     

Genistein attenuates D-galactose-induced oxidative damage through decreased reactive oxygen species and NF-κB binding activity in neuronal PC12 cells

AimsWe investigated the mechanism of D-galactose (DG)-induced oxidative damage and the neuroprotective action of genistein in PC12 cells.Main methodsPC12 cells were treated with 40 mM DG dissolved in medium containing 85% RPMI1640, 10% HBS and 5% FBS with or without genistein. We measured the protein expression of β-amyloid (Aβ), advanced glycation end products (AGEs), IκB-α and manganese-superoxide dismutase (MnSOD) by western blotting, intracellular reactive oxygen species (ROS) by 2, 7-dichlorofluorescin-diacetate, and the binding activity of nuclear factor kappa B (NF-κB) by electrophortic mobility shift assay.Key findingsDG (40 mM) completely retarded cell growth after incubation for 72 h, and this effect was not due to osmotic changes, as 40 mM mannitol had no effect. Mechanistically, we found that DG increased intracellular ROS starting at 4 h and increased Aβ and AGEs at 24 h. DG treatment for 24 h also increased the binding activity of NF-κB but strongly decreased the expression of IκB-α protein. Furthermore, DG treatment for 48 h increased MnSOD protein expression. All these effects of DG were effectively inhibited by genistein (0.5–10 μM).SignificanceThe present study indicates that the protection of genistein against DG-induced oxidative stress in PC12 cells, and the effect is likely mediated by decreased intracellular ROS and binding activity of NF-κB.     

Molecular cloning and characterization of prohormone convertase 1 gene in abalone (Haliotis diversicolor supertexta)

Prohormone convertases (PCs) are calcium-dependent serine endoproteases of the subtilisin family that play a key role in the posttranslational processing of precursors for bioactive peptides. In this study, the cDNA of PC1 from abalone (Haliotis diversicolor supertexta) was cloned and sequenced. The PC1 cDNA consisted of 2216 bp with an open reading frame of 2010 bp encoding a 670 amino acid peptide. Comparative structural analysis revealed that abalone PC1 shared high similarity and identity with most PC counterparts. The profile of deduced peptide of PC1 was composed of an N-terminal signal peptide, a prosegment domain, a catalytic domain and a P domain, which were common in many species. Sequence analysis indicated that the abalone PC1 was highly conserved in catalytic domain, including three conserved serine catalytic signatures that comprised a catalytic triad active center. Also conserved were the potential cleavage site for release of the mature peptide, a cognate integrin binding site RGD in P domain, and four cysteine residues involved in forming an intrachain disulfide bridge. To further investigate the functions of PC1 in abalone, real-time quantitative PCR was performed to determine the expression level of this gene at three different reproduction stages (i.e. pre-, during- and post-breeding). Results indicated that PC1 was expressed throughout the three stages but the expression levels varied with the timepoints and different tissues in abalone. The expression levels of PC1 in digestive gland were much higher than those of the gonad. In female abalone, the expression of PC1 was higher at pre-breeding and during-breeding stages (P < 0.05), and the expression declined at the subsequent stage. Whereas, the level of PC1 in male individual did not exhibit a significant difference in various reproduction stages. Also, the natural enzyme activity of PC1 partially exhibited a similar tendency with the mRNA expression. According to the results, it can be concluded that PC1 gene is involved in the abalone reproduction process (e.g. spawning or sperming). PC1 is a potential prohormone processing enzyme and it may play a critical role in abalone physiological processes related to reproduction.     

Identification of novel inhibitors of HCV RNA-dependent RNA polymerase by pharmacophore-based virtual screening and in vitro evaluation

Hepatitis C virus (HCV) is the major etiological agent of non-A, non-B hepatitis where no effective treatment is available. The HCV NS5B with RNA-dependent RNA polymerase (RdRp) activity is a key target for the treatment of HCV infection. Here we report novel NS5B polymerase inhibitors identified by virtual screening and in vitro evaluation of their inhibitory activities. On the basis of a newly identified binding pocket of NS5B, distinct from the nucleotide binding site but highly conserved among various HCV isolates, we performed virtual screening of compounds that fit this binding pocket from the available chemical database of 3.5 million compounds. The inhibitory activities of the in silico selected 119 compounds were estimated with in vitro RdRp assay. Three compounds with IC50 values of about 20 μM were identified, and their kinetic analyses suggest that these compounds are noncompetitive inhibitors with respect to the ribonucleotide substrate. Furthermore, the single-point mutations of the conserved residues in the binding pocket of NS5B resulted in the significant decrease of the RdRp activity, indicating that the binding pocket presented here might be important for the therapeutic intervention of HCV. These novel inhibitors would be useful for the development of effective anti-HCV agents.     

Cloning,characterization and validation of inosine 5′-monophosphate dehydrogenase of Babesia gibsoni as molecular drug target

The inosine monophosphate dehydrogenase (IMPDH) enzyme has been characterized and validated as a molecular drug target in other apicomplexans but not in the genus Babesia. Subsequently, we cloned and expressed a Babesia gibsoni IMPDH (BgIMPDH) cDNA in Escherichia coli. We also determined the inhibitory effect of mycophenolic acid (MPA) on recombinant BgIMPDH (rBgIMPDH) activity and the Babesia-growths in vitro. The translated BgIMPDH peptide contained thirteen amino acid residues responsible for substrate and cofactor binding in its catalytic domain with Gly374 in BgIMPDH being replaced by Ser388 in mammalian IMPDH. The native BgIMPDH enzyme in the parasite was approximately 54-kDa a mass similar to His-tag rBgIMPDH protein. The K_m values of the rBgIMPDH were 8.18 ± 0.878 (mean ± standard error of the mean) μM and 360.80 ± 43.41 μM for IMP and NAD⁺, respectively. MPA inhibited the rBgIMPDH activity yielding a K_i value of 20.93 ± 1.83 μM with respect to NAD⁺. For Babesia growths, the IC₅₀s were 0.95 ± 0.21 and 2.88 ± 0.49 μM for B. gibsoni and B. bovis, respectively. Therefore, our results suggest that MPA may inhibit the replication of Babesia parasites by targeting IMPDH enzyme of the purine pathway.     

Target Binding to S100B Reduces Dynamic Properties and Increases Ca2 +-Binding Affinity for Wild Type and EF-Hand Mutant Proteins

Mutations in the second EF-hand (D61N, D63N, D65N, and E72A) of S100B were used to study its Ca^2 + binding and dynamic properties in the absence and presence of a bound target, TRTK-12. With ^D63NS100B as an exception (^D63NK_D = 50 ± 9 μM), Ca^2 + binding to EF2-hand mutants were reduced by more than 8-fold in the absence of TRTK-12 (^D61NK_D = 412 ± 67 μM, ^D65NK_D = 968 ± 171 μM, and ^E72AK_D = 471 ± 133 μM), when compared to wild-type protein (^WTK_D = 56 ± 9 μM). For the TRTK-12 complexes, the Ca^2 +-binding affinity to wild type (^WT + TRTKK_D = 12 ± 10 μM) and the EF2 mutants was increased by 5- to 14-fold versus in the absence of target (^{D61N + TRTK}K_D = 29 ± 1.2 μM, ^{D63N + TRTK}K_D = 10 ± 2.2 μM, ^{D65N + TRTK}K_D = 73 ± 4.4 μM, and ^{E72A + TRTK}K_D = 18 ± 3.7 μM). In addition, R_ex, as measured using relaxation dispersion for side‐chain ¹⁵N resonances of Asn63 (^D63NS100B), was reduced upon TRTK-12 binding when measured by NMR. Likewise, backbone motions on multiple timescales (picoseconds to milliseconds) throughout wild type, ^D61NS100B, ^D63NS100B, and ^D65NS100B were lowered upon binding TRTK-12. However, the X-ray structures of Ca^2 +-bound (2.0 Å) and TRTK-bound (1.2 Å) ^D63NS100B showed no change in Ca^2 + coordination; thus, these and analogous structural data for the wild-type protein could not be used to explain how target binding increased Ca^2 +-binding affinity in solution. Therefore, a model for how S100B–TRTK‐12 complex formation increases Ca^2 + binding is discussed, which considers changes in protein dynamics upon binding the target TRTK-12.     

Role of the two ATPase domains of Escherichia coli UvrA in binding non-bulky DNA lesions and interaction with UvrB

The UvrA protein is the initial DNA damage-sensing protein in bacterial nucleotide excision repair and detects a wide variety of structurally unrelated lesions. After initial recognition of DNA damage, UvrA loads the UvrB protein onto the DNA. This protein then verifies the presence of a lesion, after which UvrA is released from the DNA.UvrA contains two ATPase domains, both belonging to the ABC ATPase superfamily. We have determined the activities of two mutants, in which a single domain was deactivated. Inactivation of either one ATPase domain in Escherichia coli UvrA results in a complete loss of ATPase activity, indicating that both domains function in a cooperative way. We could show that this ATPase activity is not required for the recognition of bulky lesions by UvrA, but it does promote the specific binding to the less distorting cyclobutane–pyrimidine dimer (CPD). The two ATPase mutants also show a difference in UvrB-loading, depending on the length of the DNA substrate. The ATPase domain I mutant was capable of loading UvrB on a lesion in a 50 bp fragment, but this loading was reduced on a longer substrate. For the ATPase domain II mutant the opposite was found: UvrB could not be loaded on a 50 bp substrate, but this loading was rescued when the length of the fragment was increased. This differential loading of UvrB by the two ATPase mutants could be related to different interactions between the UvrA and UvrB subunits.     

Molecular Recognition of Canonical and Deaminated Bases by P. abyssi Family B DNA Polymerase

Euryarchaeal polymerase B can recognize deaminated bases on the template strand, effectively stalling the replication fork 4 nt downstream the modified base. Using Pyrococcus abyssi DNA B family polymerase (PabPolB), we investigated the discrimination between deaminated and natural nucleotide(s) by primer extension assays, electrophoretic mobility shift assays, and X-ray crystallography. Structures of complexes between the protein and DNA duplexes with either a dU or a dH in position + 4 were solved at 2.3 Å and 2.9 Å resolution, respectively. The PabPolB is found in the editing mode. A new metal binding site has been uncovered below the base-checking cavity where the + 4 base is flipped out; it is fully hydrated in an octahedral fashion and helps guide the strongly kinked template strand. Four other crystal structures with each of the canonical bases were also solved in the editing mode, and the presence of three nucleotides in the exonuclease site caused a shift in the coordination state of its metal A from octahedral to tetrahedral. Surprisingly, we find that all canonical bases also enter the base-checking pocket with very small differences in the binding geometry and in the calculated binding free energy compared to deaminated ones. To explain how this can lead to stalling of the replication fork, the full catalytic pathway and its branches must be taken into account, during which the base is checked several times. Our results strongly suggest a switch from elongation to editing modes right after nucleotide insertion when the modified base is at position + 5.     

Thiazolopyridone ureas as DNA gyrase B inhibitors: Optimization of antitubercular activity and efficacy

Scaffold hopping from the thiazolopyridine ureas led to thiazolopyridone ureas with potent antitubercular activity acting through inhibition of DNA GyrB ATPase activity. Structural diversity was introduced, by extension of substituents from the thiazolopyridone N-4 position, to access hydrophobic interactions in the ribose pocket of the ATP binding region of GyrB. Further optimization of hydrogen bond interactions with arginines in site-2 of GyrB active site pocket led to potent inhibition of the enzyme (IC₅₀ 2 nM) along with potent cellular activity (MIC = 0.1 μM) against Mycobacterium tuberculosis (Mtb). Efficacy was demonstrated in an acute mouse model of tuberculosis on oral administration.     

Multiple FAS1 domains and the RGD motif of TGFBI act cooperatively to bind αvβ3 integrin,leading to anti-angiogenic and anti-tumor effects

TGFBI, a transforming growth factor β-induced extracellular matrix protein, circulates at a level of ~ 300 ng/ml in humans and modulates several integrin-mediated cellular functions. The protein contains an N-terminal EMI domain, four consecutive FAS1 domains, and the RGD motif. Each FAS1 domain and the RGD motif have been known to interact with avb3 integrin. Here, we found that the binding affinity (K_d) of TGFBI for αvβ3 integrin was approximately 3.8 × 10^? 8 M, a value ~ 2300-fold higher than that of a single FAS1 domain, and demonstrated that this greater affinity was due to the cooperative action of the four FAS1 domains and the RGD motif. Moreover, TGFBI exhibited more potent anti-angiogenic and anti-tumorigenic activities, even at a 100-fold lower molar dose than the reported effective dose of the FAS1 domain. Finally, our data showed that TGFBI specifically targeted the tumor vasculature and accumulated at the tumor site. Collectively, our results support the theory that TGFBI acts as a potent endogenous anti-tumor and anti-angiogenic molecule by targeting αvβ3 integrin, and highlights the importance of physiological circulating TGFBI levels in inhibiting tumor growth.     

Enantiomer and conformer recognition of (+) and (−)-disparlure and their analogs by the pheromone binding proteins of the gypsy moth,Lymantria dispar

Adult female gypsy moths produce a sex pheromone (+)-(7R,8S)-2-methyl-7,8-epoxyoctadecane, (+)-disparlure, to attract male gypsy moths. To better understand the recognition of (+)-disparlure by the male’s olfactory system, we synthesized racemic and enantiopure oxa and thia analogs of (+)-disparlure (ee > 98%). Ab initio calculations of the conformeric landscapes around the dihedral angles C_{5–6–7–8} and C_{7–8–9–10} of (+)-disparlure and corresponding dihedral angles of analogs revealed that introduction of the heteroatom changes the conformeric landscape around these important epitopes. The energy difference between HOMO and LUMO decreased after oxygen or sulfur was introduced into the backbone. Consistent with this, an enhancement of binding affinity between sulfur analogs and the pheromone-binding proteins (PBPs) was observed in vitro. Docking of the pheromone and analogs onto models of the two known PBPs of the gypsy moth revealed that the internal binding pocket of PBP1 showed higher selectivity than that of PBP2, consistent with in vitro binding assays. Further energy analysis revealed that enantiomers adopted different conformations with different energies when docked in the internal binding pocket of PBPs, resulting in enantiomer discrimination of PBPs towards disparlure and its analogs.