Zhongjie Liang

Protein cysteine phosphorylation of SarA/MgrA family transcriptional regulators mediates bacterial virulence and antibiotic resistance

Protein posttranslational modifications (PTMs), particularly phosphorylation, dramatically expand the complexity of cellular regulatory networks. Although cysteine (Cys) in various proteins can be subject to multiple PTMs, its phosphorylation was previously considered a rare PTM with almost no regulatory role assigned. We report here that phosphorylation occurs to a reactive cysteine residue conserved in the staphylococcal accessary regulator A (SarA)/MarR family global transcriptional regulator A (MgrA) family of proteins, and is mediated by the eukaryotic-like kinase-phosphatase pair Stk1-Stp1 in Staphylococcus aureus. Cys-phosphorylation is crucial in regulating virulence determinant production and bacterial resistance to vancomycin. Cell wall-targeting antibiotics, such as vancomycin and ceftriaxone, inhibit the kinase activity of Stk1 and lead to decreased Cys-phosphorylation of SarA and MgrA. An in vivo mouse model of infection established that the absence of stp1, which results in elevated protein Cys-phosphorylation, significantly reduces staphylococcal virulence. Our data indicate that Cys-phosphorylation is a unique PTM that can play crucial roles in bacterial signaling and regulation.

Computational drug discovery

Computational drug discovery is an effective strategy for accelerating and economizing drug discovery and development process. Because of the dramatic increase in the availability of biological macromolecule and small molecule information, the applicability of computational drug discovery has been extended and broadly applied to nearly every stage in the drug discovery and development workflow, including target identification and validation, lead discovery and optimization and preclinical tests. Over the past decades, computational drug discovery methods such as molecular docking, pharmacophore modeling and mapping, de novo design, molecular similarity calculation and sequence-based virtual screening have been greatly improved. In this review, we present an overview of these important computational methods, platforms and successful applications in this field.

Extensive crosstalk between O-GlcNAcylation and phosphorylation regulates Akt signaling.

O-linked N-acetylglucosamine glycosylations (O-GlcNAc) and O-linked phosphorylations (O-phosphate), as two important types of post-translational modifications, often occur on the same protein and bear a reciprocal relationship. In addition to the well documented phosphorylations that control Akt activity, Akt also undergoes O-GlcNAcylation, but the interplay between these two modifications and the biological significance remain unclear, largely due to the technique challenges. Here, we applied a two-step analytic approach composed of the O-GlcNAc immunoenrichment and subsequent O-phosphate immunodetection. Such an easy method enabled us to visualize endogenous glycosylated and phosphorylated Akt subpopulations in parallel and observed the inhibitory effect of Akt O-GlcNAcylations on its phosphorylation. Further studies utilizing mass spectrometry and mutagenesis approaches showed that O-GlcNAcylations at Thr 305 and Thr 312 inhibited Akt phosphorylation at Thr 308 via disrupting the interaction between Akt and PDK1. The impaired Akt activation in turn resulted in the compromised biological functions of Akt, as evidenced by suppressed cell proliferation and migration capabilities. Together, this study revealed an extensive crosstalk between O-GlcNAcylations and phosphorylations of Akt and demonstrated O-GlcNAcylation as a new regulatory modification for Akt signaling.

Molecular Basis of NDM-1, a New Antibiotic Resistance Determinant

The New Delhi Metallo-β-lactamase (NDM-1) was first reported in 2009 in a Swedish patient. A recent study reported that Klebsiella pneumonia NDM-1 positive strain or Escherichia coli NDM-1 positive strain was highly resistant to all antibiotics tested except tigecycline and colistin. These can no longer be relied on to treat infections and therefore, NDM-1 now becomes potentially a major global health threat. In this study, we performed modeling studies to obtain its 3D structure and NDM-1/antibiotics complex. It revealed that the hydrolytic mechanisms are highly conserved. In addition, the detailed analysis indicates that the more flexible and hydrophobic loop1, together with the evolution of more positive-charged loop2 leads to NDM-1 positive strain more potent and extensive in antibiotics resistance compared with other MBLs. Furthermore, through biological experiments, we revealed the molecular basis for antibiotics catalysis of NDM-1 on the enzymatic level. We found that NDM-1 enzyme was highly potent to degrade carbapenem antibiotics, while mostly susceptible to tigecycline, which had the ability to slow down the hydrolysis velocity of meropenem by NDM-1. Meanwhile, the mutagenesis experiments, including D124A, C208A, K211A and K211E, which displayed down-regulation on meropenem catalysis, proved the accuracy of our model. At present, there are no effective antibiotics against NDM-1 positive pathogen. Our study will provide clues to investigate the molecular basis of extended antibiotics resistance of NDM-1 and then accelerate the search for new antibiotics against NDM-1 positive strain in clinical studies.

Catalytic Mechanism Investigation of Lysine-Specific Demethylase 1 (LSD1): A Computational Study

Lysine-specific demethylase 1 (LSD1), the first identified histone demethylase, is a flavin-dependent amine oxidase which specifically demethylates mono- or dimethylated H3K4 and H3K9 via a redox process. It participates in a broad spectrum of biological processes and is of high importance in cell proliferation, adipogenesis, spermatogenesis, chromosome segregation and embryonic development. To date, as a potential drug target for discovering anti-tumor drugs, the medical significance of LSD1 has been greatly appreciated. However, the catalytic mechanism for the rate-limiting reductive half-reaction in demethylation remains controversial. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations, the catalytic mechanism of dimethylated H3K4 demethylation by LSD1 was characterized in details. The three-dimensional (3D) model of the complex was composed of LSD1, CoREST, and histone substrate. A 30-ns MD simulation of the model highlights the pivotal role of the conserved Tyr761 and lysine-water-flavin motif in properly orienting flavin adenine dinucleotide (FAD) with respect to substrate. The synergy of the two factors effectively stabilizes the catalytic environment and facilitated the demethylation reaction. On the basis of the reasonable consistence between simulation results and available mutagenesis data, QM/MM strategy was further employed to probe the catalytic mechanism of the reductive half-reaction in demethylation. The characteristics of the demethylation pathway determined by the potential energy surface and charge distribution analysis indicates that this reaction belongs to the direct hydride transfer mechanism. Our study provides insights into the LSD1 mechanism of reductive half-reaction in demethylation and has important implications for the discovery of regulators against LSD1 enzymes.

Identification of pentacyclic triterpenes derivatives as potent inhibitors against glycogen phosphorylase based on 3D-QSAR studies

Naturally occurring pentacyclic triterpenes (PT), a novel class of inhibitors against glycogen phosphorylase (GP), hold promise for the treatment of type-2 diabetes and other diseases with disorders in glycogen metabolism. To identify novel and more potent GP inhibitors, the receptor-based comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) approaches were performed to investigate the quantitative structure-activity relationships (QSAR) among 106 PT analogues. The validated models demonstrated that the elongated or bulky substitutions in C17 position and/or C2, C3 positions are favorable. Then based on the structural information extracted from these models, 56 derivatives were synthesized and biochemically tested in this study. The IC50 value of the most potent compound P50 was found to be 1.1 μM.

Structural Optimization and Biological Evaluation of Substituted Bisphenol A Derivatives as β-Amyloid Peptide Aggregation Inhibitors

The aggregation of Abeta is a crucial step in the etiology of Alzheimers disease. Our previous work showed that Abeta undergoes alpha-helix/beta-sheet intermediate structures during the conformational transition, and an Abeta aggregation inhibitor (1) was discovered by targeting the intermediates. Here, structure optimization toward compound 1 was performed and 34 novel derivatives were designed and synthesized. Nine compounds showed more effective inhibitory activity than the hit compound 1 in ThT fluorescence assay. Among them, compound 43 demonstrated more excellent inhibitory potency, which not only can suppress the aggregation of Abeta but also can dissolve the preformed fibrils as shown by CD spectroscopy, PICUP and AFM assays. Cellular assay indicated that 43 has no toxicity to neuronal cells, moreover, can effectively inhibit Abeta(1-42)-induced neutrotoxicity and increase the cell viability. Together, on the basis of these positive results, these novel chemical structures may provide a promising potential for therapeutic applications in AD and other types of neurodegenerative disorders.

Investigation of the acetylation mechanism by GCN5 histone acetyltransferase.

The histone acetylation of post-translational modification can be highly dynamic and play a crucial role in regulating cellular proliferation, survival, differentiation and motility. Of the enzymes that mediate post-translation modifications, the GCN5 of the histone acetyltransferase (HAT) proteins family that add acetyl groups to target lysine residues within histones, has been most extensively studied. According to the mechanism studies of GCN5 related proteins, two key processes, deprotonation and acetylation, must be involved. However, as a fundamental issue, the structure of hGCN5/AcCoA/pH3 remains elusive. Although biological experiments have proved that GCN5 mediates the acetylation process through the sequential mechanism pathway, a dynamic view of the catalytic process and the molecular basis for hGCN5/AcCoA/pH3 are still not available and none of theoretical studies has been reported to other related enzymes in HAT family. To explore the molecular basis for the catalytic mechanism, computational approaches including molecular modeling, molecular dynamic (MD) simulation and quantum mechanics/molecular mechanics (QM/MM) simulation were carried out. The initial hGCN5/AcCoA/pH3 complex structure was modeled and a reasonable snapshot was extracted from the trajectory of a 20 ns MD simulation, with considering post-MD analysis and reported experimental results. Those residues playing crucial roles in binding affinity and acetylation reaction were comprehensively investigated. It demonstrated Glu80 acted as the general base for deprotonation of Lys171 from H3. Furthermore, the two-dimensional QM/MM potential energy surface was employed to study the sequential pathway acetylation mechanism. Energy barriers of addition-elimination reaction in acetylation obtained from QM/MM calculation indicated the point of the intermediate ternary complex. Our study may provide insights into the detailed mechanism for acetylation reaction of GCN5, and has important implications for the discovery of regulators against GCN5 enzymes and related HAT family enzymes.

Discovery of novel 2-aminopyridine-3-carboxamides as c-Met kinase inhibitors.

A series of 2-aminopyridine-3-carboxamide derivatives against c-Met were designed and synthesized by employing bioisosteric replacement of heterocyclic moieties with the amide bond. The structure-activity relationship (SAR) at various positions of the scaffold was explored. In this study, a promising compound (S)-24o with a c-Met IC(50) of 0.022 μM was identified. The compound exhibited dose-dependent inhibition of the phosphorylation of c-Met as well as downstream signaling in EBC-1 cells. Furthermore, the interactive binding model of (S)-24o with c-Met was elucidated by virtue of a molecular modeling study.

A quantum mechanics/molecular mechanics study on the hydrolysis mechanism of New Delhi metallo-β-lactamase-1

New Delhi metallo-β-lactamase-1 (NDM-1) has emerged as a major global threat to human health for its rapid rate of dissemination and ability to make pathogenic microbes resistant to almost all known β-lactam antibiotics. In addition, effective NDM-1 inhibitors have not been identified to date. In spite of the plethora of structural and kinetic data available, the accurate molecular characteristics of and details on the enzymatic reaction of NDM-1 hydrolyzing β-lactam antibiotics remain incompletely understood. In this study, a combined computational approach including molecular docking, molecular dynamics simulations and quantum mechanics/molecular mechanics calculations was performed to characterize the catalytic mechanism of meropenem catalyzed by NDM-1. The quantum mechanics/molecular mechanics results indicate that the ionized D124 is beneficial to the cleavage of the C–N bond within the β-lactam ring. Meanwhile, it is energetically favorable to form an intermediate if no water molecule coordinates to Zn2. Moreover, according to the molecular dynamics results, the conserved residue K211 plays a pivotal role in substrate binding and catalysis, which is quite consistent with previous mutagenesis data. Our study provides detailed insights into the catalytic mechanism of NDM-1 hydrolyzing meropenem β-lactam antibiotics and offers clues for the discovery of new antibiotics against NDM-1 positive strains in clinical studies.