Structural and functional consequences of single amino acid substitutions in the pyrimidine base binding pocket of Escherichia coli CMP kinase

Multi-Enzymatic Cascade One-Pot Biosynthesis of 3'-Sialyllactose Using Engineered Escherichia coli

Among the human milk oligosaccharides (HMOs), one of the most abundant oligosaccharides and has great benefits for human health is 3′-sialyllactose (3′-SL). Given its important physiological functions and the lack of cost-effective production processes, we constructed an in vitro multi-enzymatic cofactor recycling system for the biosynthesis of 3′-SL from a low-cost substrate. First, we constructed the biosynthetic pathway and increased the solubility of cytidine monophosphate kinase (CMK) with chaperones. We subsequently identified that β-galactosidase (lacZ) affects the yield of 3′-SL, and hence with the lacZ gene knocked out, a 3.3-fold increase in the production of 3′-SL was observed. Further, temperature, pH, polyphosphate concentration, and concentration of divalent metal ions for 3′-SL production were optimized. Finally, an efficient biotransformation system was established under the optimized conditions. The maximum production of 3′-SL reached 38.7 mM, and a molar yield of 97.1% from N-acetylneuraminic acid (NeuAc, sialic acid, SA) was obtained. The results demonstrate that the multi-enzymatic cascade biosynthetic pathway with cofactor regeneration holds promise as an industrial strategy for producing 3′-SL.

Mycobacterium Cytidylate Kinase Appears to Be an Undruggable Target

New and improved drugs against tuberculosis are urgently needed as multi-drug-resistant forms of the disease become more prevalent. Mycobacterium tuberculosis cytidylate kinase is an attractive target for screening due to its essentiality and different substrate specificity to the human orthologue. However, we selected the Mycobacterium smegmatis cytidylate kinase for screening because of the availability of high-resolution X-ray crystallographic data defining its structure and the high likelihood of active site structural similarity to the M. tuberculosis orthologue. We report the development and implementation of a high-throughput luciferase-based activity assay and screening of 19,920 compounds derived from small-molecule libraries and an in silico screen predicting likely inhibitors of the cytidylate kinase enzyme. Hit validation included a counterscreen for luciferase inhibitors that would result in false positives in the initial screen. Results of this counterscreen ruled out all of the putative cytidylate kinase inhibitors identified in the initial screening, leaving no compounds as candidates for drug development. Although a negative result, this study indicates that this important drug target may in fact be undruggable and serve as a warning for future investigations.

Mutations in the KDM5C ARID Domain and Their Plausible Association with Syndromic Claes-Jensen-Type Disease

Mutations in KDM5C gene are linked to X-linked mental retardation, the syndromic Claes-Jensen-type disease. This study focuses on non-synonymous mutations in the KDM5C ARID domain and evaluates the effects of two disease-associated missense mutations (A77T and D87G) and three not-yet-classified missense mutations (R108W, N142S, and R179H). We predict the ARID domain’s folding and binding free energy changes due to mutations, and also study the effects of mutations on protein dynamics. Our computational results indicate that A77T and D87G mutants have minimal effect on the KDM5C ARID domain stability and DNA binding. In parallel, the change in the free energy unfolding caused by the mutants A77T and D87G were experimentally measured by urea-induced unfolding experiments and were shown to be similar to the in silico predictions. The evolutionary conservation analysis shows that the disease-associated mutations are located in a highly-conserved part of the ARID structure (N-terminal domain), indicating their importance for the KDM5C function. N-terminal residues’ high conservation suggests that either the ARID domain utilizes the N-terminal to interact with other KDM5C domains or the N-terminal is involved in some yet unknown function. The analysis indicates that, among the non-classified mutations, R108W is possibly a disease-associated mutation, while N142S and R179H are probably harmless.

Doing that thing that scientists do: A discovery-driven module on protein purification and characterization for the undergraduate biochemistry laboratory classroom

In traditional introductory biochemistry laboratory classes students learn techniques for protein purification and analysis by following provided, established, step-by-step procedures. Students are exposed to a variety of biochemical techniques but are often not developing procedures or collecting new, original data. In this laboratory module, students develop research skills through work on an original research project and gain confidence in their ability to design and execute an experiment while faculty can enhance their scholarly pursuits through the acquisition of original data in the classroom laboratory. Students are prepared for a 6-8 week discovery-driven project on the purification of the Escherichia coli cytidylate kinase (CMP kinase) through in class problems and other laboratory exercises on bioinformatics and protein structure analysis. After a minimal amount of guidance on how to perform the CMP kinase in vitro enzyme assay, SDS-PAGE, and the basics of protein purification, students, working in groups of three to four, develop a protein purification protocol based on the scientific literature and investigate some aspect of CMP kinase that interests them. Through this process, students learn how to implement a new but perhaps previously worked out procedure to answer their research question. In addition, they learn the importance of keeping a clear and thorough laboratory notebook and how to interpret their data and use that data to inform the next set of experiments. Following this module, students had increased confidence in their ability to do basic biochemistry techniques and reported that the "self-directed" nature of this lab increased their engagement in the project.

Increasing the structural coverage of tuberculosis drug targets

High-resolution three-dimensional structures of essential Mycobacterium tuberculosis (Mtb) proteins provide templates for TB drug design, but are available for only a small fraction of the Mtb proteome. Here we evaluate an intra-genus "homolog-rescue" strategy to increase the structural information available for TB drug discovery by using mycobacterial homologs with conserved active sites. Of 179 potential TB drug targets selected for x-ray structure determination, only 16 yielded a crystal structure. By adding 1675 homologs from nine other mycobacterial species to the pipeline, structures representing an additional 52 otherwise intractable targets were solved. To determine whether these homolog structures would be useful surrogates in TB drug design, we compared the active sites of 106 pairs of Mtb and non-TB mycobacterial (NTM) enzyme homologs with experimentally determined structures, using three metrics of active site similarity, including superposition of continuous pharmacophoric property distributions. Pair-wise structural comparisons revealed that 19/22 pairs with >55% overall sequence identity had active site Cα RMSD <1 Å, >85% side chain identity, and ≥80% PSAPF (similarity based on pharmacophoric properties) indicating highly conserved active site shape and chemistry. Applying these results to the 52 NTM structures described above, 41 shared >55% sequence identity with the Mtb target, thus increasing the effective structural coverage of the 179 Mtb targets over three-fold (from 9% to 32%). The utility of these structures in TB drug design can be tested by designing inhibitors using the homolog structure and assaying the cognate Mtb enzyme; a promising test case, Mtb cytidylate kinase, is described. The homolog-rescue strategy evaluated here for TB is also generalizable to drug targets for other diseases.

Kinetic mechanism and energetics of binding of phosphoryl group acceptors to Mycobacterium tuberculosis cytidine monophosphate kinase

Cytidine monophosphate kinase from Mycobacterium tuberculosis (MtCMK) likely plays a role in supplying precursors for nucleic acid synthesis. MtCMK catalyzes the ATP-dependent phosphoryl group transfer preferentially to CMP and dCMP. Initial velocity studies and Isothermal titration calorimetry (ITC) measurements showed that MtCMK follows a random-order mechanism of substrate (CMP and ATP) binding, and an ordered mechanism for product release, in which ADP is released first followed by CDP. The thermodynamic signatures of CMP and CDP binding to MtCMK showed favorable enthalpy and unfavorable entropy, and ATP binding was characterized by favorable changes in enthalpy and entropy. The contribution of linked protonation events to the energetics of MtCMK:phosphoryl group acceptor binary complex formation suggested a net gain of protons. Values for the pKa of a likely chemical group involved in proton exchange and for the intrinsic binding enthalpy were calculated. The Asp187 side chain of MtCMK is suggested as the likely candidate for the protonation event. Data on thermodynamics of binary complex formation were collected to evaluate the contribution of 2'-OH group to intermolecular interactions. The data are discussed in light of functional and structural comparisons between CMP/dCMP kinases and UMP/CMP ones.

Predicting folding free energy changes upon single point mutations

MOTIVATION

The folding free energy is an important characteristic of proteins stability and is directly related to protein's wild-type function. The changes of protein's stability due to naturally occurring mutations, missense mutations, are typically causing diseases. Single point mutations made in vitro are frequently used to assess the contribution of given amino acid to the stability of the protein. In both cases, it is desirable to predict the change of the folding free energy upon single point mutations in order to either provide insights of the molecular mechanism of the change or to design new experimental studies.

RESULTS

We report an approach that predicts the free energy change upon single point mutation by utilizing the 3D structure of the wild-type protein. It is based on variation of the molecular mechanics Generalized Born (MMGB) method, scaled with optimized parameters (sMMGB) and utilizing specific model of unfolded state. The corresponding mutations are built in silico and the predictions are tested against large dataset of 1109 mutations with experimentally measured changes of the folding free energy. Benchmarking resulted in root mean square deviation = 1.78 kcal/mol and slope of the linear regression fit between the experimental data and the calculations was 1.04. The sMMGB is compared with other leading methods of predicting folding free energy changes upon single mutations and results discussed with respect to various parameters.

AVAILABILITY

All the pdb files we used in this article can be downloaded from http://compbio.clemson.edu/downloadDir/mentaldisorders/sMMGB_pdb.rar.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In silico and in vitro investigations of the mutability of disease-causing missense mutation sites in spermine synthase

Background

Spermine synthase (SMS) is a key enzyme controlling the concentration of spermidine and spermine in the cell. The importance of SMS is manifested by the fact that single missense mutations were found to cause Snyder-Robinson Syndrome (SRS). At the same time, currently there are no non-synonymous single nucleoside polymorphisms, nsSNPs (harmless mutations), found in SMS, which may imply that the SMS does not tolerate amino acid substitutions, i.e. is not mutable.

Methodology/Principal Findings

To investigate the mutability of the SMS, we carried out in silico analysis and in vitro experiments of the effects of amino acid substitutions at the missense mutation sites (G56, V132 and I150) that have been shown to cause SRS. Our investigation showed that the mutation sites have different degree of mutability depending on their structural micro-environment and involvement in the function and structural integrity of the SMS. It was found that the I150 site does not tolerate any mutation, while V132, despite its key position at the interface of SMS dimer, is quite mutable. The G56 site is in the middle of the spectra, but still quite sensitive to charge residue replacement.

Conclusions/Significance

The performed analysis showed that mutability depends on the detail of the structural and functional factors and cannot be predicted based on conservation of wild type properties alone. Also, harmless nsSNPs can be expected to occur even at sites at which missense mutations were found to cause diseases.

Computational analysis of missense mutations causing Snyder-Robinson syndrome

The Snyder-Robinson syndrome is caused by missense mutations in the spermine sythase gene that encodes a protein (SMS) of 529 amino acids. Here we investigate, in silico, the molecular effect of three missense mutations, c.267G>A (p.G56S), c.496T>G (p.V132G), and c.550T>C (p.I150T) in SMS that were clinically identified to cause the disease. Single-point energy calculations, molecular dynamics simulations, and pKa calculations revealed the effects of these mutations on SMS's stability, flexibility, and interactions. It was predicted that the catalytic residue, Asp276, should be protonated prior binding the substrates. The pKa calculations indicated the p.I150T mutation causes pKa changes with respect to the wild-type SMS, which involve titratable residues interacting with the S-methyl-5'-thioadenosine (MTA) substrate. The p.I150T missense mutation was also found to decrease the stability of the C-terminal domain and to induce structural changes in the vicinity of the MTA binding site. The other two missense mutations, p.G56S and p.V132G, are away from active site and do not perturb its wild-type properties, but affect the stability of both the monomers and the dimer. Specifically, the p.G56S mutation is predicted to greatly reduce the affinity of monomers to form a dimer, and therefore should have a dramatic effect on SMS function because dimerization is essential for SMS activity.

Modeling effects of human single nucleotide polymorphisms on protein-protein interactions

A large set of three-dimensional structures of 264 protein-protein complexes with known nonsynonymous single nucleotide polymorphisms (nsSNPs) at the interface was built using homology-based methods. The nsSNPs were mapped on the proteins' structures and their effect on the binding energy was investigated with CHARMM force field and continuum electrostatic calculations. Two sets of nsSNPs were studied: disease annotated Online Mendelian Inheritance in Man (OMIM) and nonannotated (non-OMIM). It was demonstrated that OMIM nsSNPs tend to destabilize the electrostatic component of the binding energy, in contrast with the effect of non-OMIM nsSNPs. In addition, it was shown that the change of the binding energy upon amino acid substitutions is not related to the conservation of the net charge, hydrophobicity, or hydrogen bond network at the interface. The results indicate that, generally, the effect of nsSNPs on protein-protein interactions cannot be predicted from amino acids' physico-chemical properties alone, since in many cases a substitution of a particular residue with another amino acid having completely different polarity or hydrophobicity had little effect on the binding energy. Analysis of sequence conservation showed that nsSNP at highly conserved positions resulted in a large variance of the binding energy changes. In contrast, amino acid substitutions corresponding to nsSNPs at nonconserved positions, on average, were not found to have a large effect on binding affinity. pKa calculations were performed and showed that amino acid substitutions could change the wild-type proton uptake/release and thus resulting in different pH-dependence of the binding energy.

The Rv1712 Locus from Mycobacterium tuberculosis H37Rv codes for a functional CMP kinase that preferentially phosphorylates dCMP

The Mycobacterium tuberculosis cmk gene, predicted to encode a CMP kinase (CMK), was cloned and expressed, and its product was purified to homogeneity. Steady-state kinetics confirmed that M. tuberculosis CMK is a monomer that preferentially phosphorylates CMP and dCMP by a sequential mechanism. A plausible role for CMK is discussed.