Identification of Serum Biomarkers for Systemic Lupus Erythematosus Using a Library of Phage Displayed Random Peptides and Deep Sequencing.

Systemic lupus erythematosus (SLE) is one of the most serious autoimmune diseases, characterized by highly diverse clinical manifestations. A biomarker is still needed for accurate diagnostics. SLE serum autoantibodies were discovered and validated using serum samples from independent sample cohorts encompassing 306 participants divided into three groups, i.e. healthy, SLE patients, and other autoimmune-related diseases. To discover biomarkers for SLE, a phage displayed random peptide library (Ph.D. 12) and deep sequencing were applied to screen specific autoantibodies in a total of 100 serum samples from 50 SLE patients and 50 healthy controls. A statistical analysis protocol was set up for the identification of peptides as potential biomarkers. For validation, 10 peptides were analyzed using enzyme-linked immunosorbent assays (ELISA). As a result, four peptides (SLE2018Val001, SLE2018Val002, SLE2018Val006, and SLE2018Val008) were discovered with high diagnostic power to differentiate SLE patients from healthy controls. Among them, two peptides, i.e. SLE2018Val001 and SLE2018Val002, were confirmed between SLE with other autoimmune patients. The procedure we established could be easily adopted for the identification of autoantibodies as biomarkers for many other diseases.

Cell Lysate Microarray for Mapping the Network of Genetic Regulators for Histone Marks.

Proteins, as the major executer for cell progresses and functions, its abundance and the level of post-translational modifications, are tightly monitored by regulators. Genetic perturbation could help us to understand the relationships between genes and protein functions. Herein, to explore the impact of the genome-wide interruption on certain protein, we developed a cell lysate microarray on kilo-conditions (CLICK) with 4837 knockout (YKO) and 322 temperature-sensitive (ts) mutant strains of yeast (Saccharomyces cerevisiae). Taking histone marks as examples, a general workflow was established for the global identification of upstream regulators. Through a single CLICK array test, we obtained a series of regulators for H3K4me3, which covers most of the known regulators in S. cerevisiae We also noted that several group of proteins are involved in negatively regulation of H3K4me3. Further, we discovered that Cab4p and Cab5p, two key enzymes of CoA biosynthesis, play central roles in histone acylation. Because of its general applicability, CLICK array could be easily adopted to rapid and global identification of upstream protein/enzyme(s) that regulate/modify the level of a protein or the posttranslational modification of a non-histone protein.

Catalyst-free synthesis of novel dimeric ß-carboline derivatives via an unexpected [2 + 2 + 2] annulation.

We report here an unexpected catalyst-free [2 + 2 + 2] annulation reaction which allows access to novel complex dimeric ß-carboline derivatives in a single step. Various substituted ynones could react with dihydro-ß-carboline imines to deliver interesting [2 + 2 + 2] annulation products in moderate to good yields. Alkynoates can also be tolerated in this system.

[Preparation of SUN5-specific polyclonal antibody for detection of SUN5 expression in human germ cells].

OBJECTIVE: To prepare a specific polyclonal antibody against full-length SUN5 for detecting the expression of SUN5 in human germ cells. METHODS: Bioinformatic methods were used to compare the full-length SUN5 and its variant SUN5ß, and a short peptide was designed based on the differential region to prepare SUN5 antibody. The prepared antibody was used to detect the expression of SUN5 in Ntera-2 cells and in human germ cells by Western blotting and immunofluorescence assay. RESULTS: The short peptide was correctly synthesized and SUN5 antibody was obtained and purified. Western blotting showed that the prepared antibody was capable of recognizing full-length SUN5 in Ntera-2 cells, and SUN5 expression was localized on the nuclear membrane and in the cytoplasm as shown by immunofluorescence assay. Using this antibody, we detected SUN5 expression in the spermatocytes, round spermatids and sperms in human germ cells. CONCLUSION: We successfully prepared SUN5-specific antibody. SUN5 is expressed in the spermatocytes, round spermatids and sperms in human germ cells, suggesting its important role in spermatogenesis.

Systematic identification of arsenic-binding proteins reveals that hexokinase-2 is inhibited by arsenic.

Arsenic is highly effective for treating acute promyelocytic leukemia (APL) and has shown significant promise against many other tumors. However, although its mechanistic effects in APL are established, its broader anticancer mode of action is not understood. In this study, using a human proteome microarray, we identified 360 proteins that specifically bind arsenic. Among the most highly enriched proteins in this set are those in the glycolysis pathway, including the rate-limiting enzyme in glycolysis, hexokinase-1. Detailed biochemical and metabolomics analyses of the highly homologous hexokinase-2 (HK2), which is overexpressed in many cancers, revealed significant inhibition by arsenic. Furthermore, overexpression of HK2 rescued cells from arsenic-induced apoptosis. Our results thus strongly implicate glycolysis, and HK2 in particular, as a key target of arsenic. Moreover, the arsenic-binding proteins identified in this work are expected to serve as a valuable resource for the development of synergistic antitumor therapeutic strategies.

Elucidation of Enzymatic Mechanism of Phenazine Biosynthetic Protein PhzF Using QM/MM and MD Simulations.

The phenazine biosynthetic pathway is of considerable importance for the pharmaceutical industry. The pathway produces two products: phenazine-1,6-dicarboxylic acid and phenazine-1-carboxylic acid. PhzF is an isomerase that catalyzes trans-2,3-dihydro-3-hydroxyanthranilic acid isomerization and plays an essential role in the phenazine biosynthetic pathway. Although the PhzF crystal structure has been determined recently, an understanding of the detailed catalytic mechanism and the roles of key catalytic residues are still lacking. In this study, a computational strategy using a combination of molecular modeling, molecular dynamics simulations, and quantum mechanics/molecular mechanics simulations was used to elucidate these important issues. The Apo enzyme, enzyme-substrate complexes with negatively charged Glu45, enzyme-transition state analog inhibitor complexes with neutral Glu45, and enzyme-product complexes with negatively charged Glu45 structures were optimized and modeled using a 200 ns molecular dynamics simulation. Residues such as Gly73, His74, Asp208, Gly212, Ser213, and water, which play important roles in ligand binding and the isomerization reaction, were comprehensively investigated. Our results suggest that the Glu45 residue at the active site of PhzF acts as a general base/acid catalyst during proton transfer. This study provides new insights into the detailed catalytic mechanism of PhzF and the results have important implications for PhzF modification.

Mycobacterium tuberculosis proteome microarray for global studies of protein function and immunogenicity.

Poor understanding of the basic biology of Mycobacterium tuberculosis (MTB), the etiological agent of tuberculosis, hampers development of much-needed drugs, vaccines, and diagnostic tests. Better experimental tools are needed to expedite investigations of this pathogen at the systems level. Here, we present a functional MTB proteome microarray covering most of the proteome and an ORFome library. We demonstrate the broad applicability of the microarray by investigating global protein-protein interactions, small-molecule-protein binding, and serum biomarker discovery, identifying 59 PknG-interacting proteins, 30 bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) binding proteins, and 14 MTB proteins that together differentiate between tuberculosis (TB) patients with active disease and recovered individuals. Results suggest that the MTB rhamnose pathway is likely regulated by both the serine/threonine kinase PknG and c-di-GMP. This resource has the potential to generate a greater understanding of key biological processes in the pathogenesis of tuberculosis, possibly leading to more effective therapies for the treatment of this ancient disease.

Research/review: Insights into the mutation-induced dysfunction of arachidonic acid metabolism from modeling of human CYP2J2.

As a kind of monooxygenase with the function of catalyzing many reactions involved in drug metabolism and synthesis of cholesterol, steroids and other lipids, CYP2J2 is an important member of the cytochrome P450 superfamily. Located at the endoplasmic reticulum, CYP2J2 is responsible for epoxidation of endogenous arachidonic acid in cardiac tissue to produce cis-epoxyeicosatrienoic acids (EETs), which have anti-inflammatory and antifibrinolytic properties, and can protect endothelial cells from ischemic or hypoxic injuries. Some polymorphisms, e.g., CYP2J2 with mutation T143A, R158C, I192N or N404Y, could significantly reduce the metabolism of the arachidonic acid, causing or deteriorating the coronary artery disease. However, so far the detailed mechanism for the mutationinduced dysfunction of arachidonic metabolism is still unknown. To reveal its mechanism, a 3D (three-dimensional) structure for human CYP2J2 was developed, followed by docking the arachidonic acid ligand into the active site of the receptor. It was observed based on the binding mode thus found that Gly486 and Leu378 in the active site of the receptor played a key role in recognizing and positioning the carboxyl group of the ligand via hydrogen bonding interactions, and that any of the aforementioned five mutations might have, either directly or indirectly, impact to their role and hence causing the mutation-induced dysfunction of CYP2J2-mediated arachidonic acid metabolism. It is anticipated that the findings as reported in this review article may stimulate new strategy for finding novel therapeutic approaches to treat coronary artery disease.

Research/review: Structure and linkage disequilibrium analysis of adamantane resistant mutations in influenza virus m2 proton channel.

The M2 proton channel is translated by the M gene segment of influenza viruses, and has been adopted as an attractive target for influenza A viruses, on which a series of adamantane-based drugs act. However, recently epidemic influenza viruses have had strong resistant effects against the adamantane-based drugs. In this paper, we combined evolutionary analyses, linkage disequilibrium as well as molecular dynamics simulations to explore the drug resistance of the M2 proton channel, with an aim of providing an in-depth understanding of the resistant mechanism for adamantane-based drugs. We collected 2746 coding sequences for swine, avian, and human M2 proteins. After evolutionary and linkage disequilibrium analyses, we found that the some residues in the C-terminal were associated with the famed resistant mutation S31N. Subsequently, we constructed the 3D structures of the swine, avian as well as human M2 channel, and performed MD simulations on these channels with a typical adamantane-based drug rimantadine. From the simulation trajectories, we found that the resistance against the adamantane-based drugs for the M2 channel from 2009 A(H1N1) viruses was derived from the structural allostery in the transmembrane and C-terminal regions. The helices in the transmembrane region were irregular in formation and employed larger distances between the adjacent 2 helices, which can weaken the interactions between the adjacent 2 helices and destabilize the helix-helix assembly, resulting in a comparatively loosely structure. The helices in the C-terminal region show a disordered configuration, giving chances for solvent molecules to enter into the channel pore.

Microsatellites in the genome of the edible mushroom, Volvariella volvacea.

Using bioinformatics software and database, we have characterized the microsatellite pattern in the V. volvacea genome and compared it with microsatellite patterns found in the genomes of four other edible fungi: Coprinopsis cinerea, Schizophyllum commune, Agaricus bisporus, and Pleurotus ostreatus. A total of 1346 microsatellites have been identified, with mono-nucleotides being the most frequent motif. The relative abundance of microsatellites was lower in coding regions with 21 No./Mb. However, the microsatellites in the V. volvacea gene models showed a greater tendency to be located in the CDS regions. There was also a higher preponderance of trinucleotide repeats, especially in the kinase genes, which implied a possible role in phenotypic variation. Among the five fungal genomes, microsatellite abundance appeared to be unrelated to genome size. Furthermore, the short motifs (mono- to tri-nucleotides) outnumbered other categories although these differed in proportion. Data analysis indicated a possible relationship between the most frequent microsatellite types and the genetic distance between the five fungal genomes.

Expression sensitivity analysis of human disease related genes.

Background. Genome-wide association studies (GWAS) have shown its revolutionary power in seeking the influenced loci on complex diseases genetically. Thousands of replicated loci for common traits are helpful in diseases risk assessment. However it is still difficult to elucidate the variations in these loci that directly cause susceptibility to diseases by disrupting the expression or function of a protein currently. Results. We evaluate the expression features of disease related genes and find that different diseases related genes show different expression perturbation sensitivities in various conditions. It is worth noting that the expression of some robust disease-genes doesn't show significant change in their corresponding diseases, these genes might be easily ignored in the expression profile analysis. Conclusion. Gene ontology enrichment analysis indicates that robust disease-genes execute essential function in comparison with sensitive disease-genes. The diseases associated with robust genes seem to be relatively lethal like cancer and aging. On the other hand, the diseases associated with sensitive genes are apparently nonlethal like psych and chemical dependency diseases.

Molecular dynamics studies on the conformational transitions of adenylate kinase: a computational evidence for the conformational selection mechanism.

Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme that catalyzes reversible transfer of phosphoryl group from ATP to AMP with a large-scale domain motion. The detailed mechanism for this conformational transition remains unknown. In the current study, we performed long time-scale molecular dynamics simulations on both open and closed states of ADK. Based on the structural analyses of the simulation trajectories, we detected over 20 times conformational transitions between the open and closed states of ADK and identified two novel conformations as intermediate states in the catalytic processes. With these findings, we proposed a possible mechanism for the large-scale domain motion of Escherichia coli ADK and its catalytic process: (1) the substrate free ADK adopted an open conformation; (2) ATP bound with LID domain closure; (3) AMP bound with NMP domain closure; (4) phosphoryl transfer occurred with ATP, and AMP converted into two ADPs, and no conformational transition was detected in the enzyme; (5) LID domain opened with one ADP released; (6) another ADP released with NMP domain open. As both open and closed states sampled a wide range of conformation transitions, our simulation strongly supported the conformational selection mechanism for Escherichia coli ADK.

Structural basis for the mutation-induced dysfunction of human CYP2J2: a computational study.

Arachidonic acid is an essential fatty acid in cells, acting as a key inflammatory intermediate in inflammatory reactions. In cardiac tissues, CYP2J2 can adopt arachidonic acid as a major substrate to produce epoxyeicosatrienoic acids (EETs), which can protect endothelial cells from ischemic or hypoxic injuries and have been implicated in the pathogenesis of coronary artery disease and hypertension. However, some CYP2J2 polymorphisms, i.e., T143A and N404Y, significantly reduce the metabolism of arachidonic acid. Lacking experimental structural data for CYP2J2, the detailed mechanism for the mutation-induced dysfunction in the metabolism of arachidonic acid is still unknown. In the current study, three-dimensional structural models of the wild-type CYP2J2 and two mutants (T143A and N404Y) were constructed by a coordinate reconstruction approach and ab initio modeling using CYP2R1 as a template. The structural analysis of the computational models showed that the wild-type CYP2J2 exhibited a typical CYP fold with 12 alpha-helices and three beta-sheets on one side and with the heme group buried deeply inside the core. Due to the small and hydrophobic side-chain, T143A mutation could destabilize the C helix, further placing the water access channel in a closed state to prevent the escape of the produced water molecules during the catalytic processes. N404Y mutation could reposition the side-chain of Leu(378), making it no longer form a hydrogen bond with the carboxyl group of arachidonic acid. However, this hydrogen bond was essential for substrate recognition and positioning in a correct orientation.

Metallo-ß-lactamases: structural features, antibiotic recognition, inhibition, and inhibitor design.

Owing to their ability in destroying or slowing down the growth of bacteria, antibiotics have been widely used to treat the bacterial infections. However, because of the long-term and irresponsible use of antibiotics, resistance to antibiotics has become a serious problem directly threatening the public health worldwide. To fight against and resist ß- lactam antibiotics, bacteria usually employed ß-lactamases, especially the metallo-ß-lactamases, to hydrolyze the C-N bond of the ß-lactam ring so as to inactivate the antibiotics. In this minireview, we are to summarize the structural features of the metallo-ß-lactamases, as well as their antibiotic binding modes and resistance mechanisms, in hopes that the discussion and analysis presented in this paper can stimulate new strategies to overcome the resistance problem and find novel inhibitors against the metallo-ß-lactamases.

Molecular simulation to investigate the cofactor specificity for pichia stipitis Xylose reductase.

Xylose is one of the most abundant carbohydrates in nature, and widely used to produce bioethanol via fermentation in industry. Xylulose can produce two key enzymes: xylose reductase and xylitol dehydrogenase. Owing to the disparate cofactor specificities of xylose reductase and xylitol dehydrogenase, intracellular redox imbalance is detected during the xylose fermentation, resulting in low ethanol yields. To overcome this barrier, a common strategy is applied to artificially modify the cofactor specificity of xylose reductase. In this study, we utilized molecular simulation approaches to construct a 3D (three-dimensional) structural model for the NADP-dependent Pichia stipitis xylose reductase (PsXR). Based on the 3D model, the favourable binding modes for both cofactors NAD and NADP were obtained using the flexible docking procedure and molecular dynamics simulation. Structural analysis of the favourable binding modes showed that the cofactor binding site of PsXR was composed of 3 major components: a hydrophilic pocket, a hydrophobic pocket as well as a linker channel between the aforementioned two pockets. The hydrophilic pocket could recognize the nicotinamide moiety of the cofactors by hydrogen bonding networks, while the hydrophobic pocket functioned to position the adenine moiety of the cofactors by hydrophobic and Π-Π stacking interactions. The linker channel contained some key residues for ligand-binding; their mutation could have impact to the specificity of PsXR. Finally, it was found that any of the two single mutations, K21A and K270N, might reverse the cofactor specificity of PsXR from major NADP- to NADdependent, which was further confirmed by the additional experiments. Our findings may provide useful insights into the cofactor specificity of PsXR, stimulating new strategies for better designing xylose reductase and improving ethanol production in industry.

π-π Stacking mediated drug-drug interactions in human CYP2E1.

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU-accelerated molecular dynamics simulations to study the multiple-binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple-substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple-substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π-π stacking interactions. In the multiple-binding mode, the aforementioned π-π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple-substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs.

Computational design of glutamate dehydrogenase in Bacillus subtilis natto.

Bacillus subtilis natto is widely used in industry to produce natto, a traditional and popular Japanese soybean food. However, during its secondary fermentation, high amounts of ammonia are released to give a negative influence on the flavor of natto. Glutamate dehydrogenase (GDH) is a key enzyme for the ammonia produced and released, because it catalyzes the oxidative deamination of glutamate to alpha-ketoglutarate using NAD(+) or NADP(+) as co-factor during carbon and nitrogen metabolism processes. To solve this problem, we employed multiple computational methods model and re-design GDH from Bacillus subtilis natto. Firstly, a structure model of GDH with cofactor NADP(+) was constructed by threading and ab initio modeling. Then the substrate glutamate were flexibly docked into the structure model to form the substrate-binding mode. According to the structural analysis of the substrate-binding mode, Lys80, Lys116, Arg196, Thr200, and Ser351 in the active site were found could form a significant hydrogen bonding network with the substrate, which was thought to play a crucial role in the substrate recognition and position. Thus, these residues were then mutated into other amino acids, and the substrate binding affinities for each mutant were calculated. Finally, three single mutants (K80A, K116Q, and S351A) were found to have significant decrease in the substrate binding affinities, which was further supported by our biochemical experiments.

Scaffold-based pan-agonist design for the PPARα, PPARß and PPARγ receptors.

As important members of nuclear receptor superfamily, Peroxisome proliferator-activated receptors (PPAR) play essential roles in regulating cellular differentiation, development, metabolism, and tumorigenesis of higher organisms. The PPAR receptors have 3 identified subtypes: PPARα, PPARß and PPARγ, all of which have been treated as attractive targets for developing drugs to treat type 2 diabetes. Due to the undesirable side-effects, many PPAR agonists including PPARα/γ and PPARß/γ dual agonists are stopped by US FDA in the clinical trials. An alternative strategy is to design novel pan-agonist that can simultaneously activate PPARα, PPARß and PPARγ. Under such an idea, in the current study we adopted the core hopping algorithm and glide docking procedure to generate 7 novel compounds based on a typical PPAR pan-agonist LY465608. It was observed by the docking procedures and molecular dynamics simulations that the compounds generated by the core hopping and glide docking not only possessed the similar functions as the original LY465608 compound to activate PPARα, PPARß and PPARγ receptors, but also had more favorable conformation for binding to the PPAR receptors. The additional absorption, distribution, metabolism and excretion (ADME) predictions showed that the 7 compounds (especially Cpd#1) hold high potential to be novel lead compounds for the PPAR pan-agonist. Our findings can provide a new strategy or useful insights for designing the effective pan-agonists against the type 2 diabetes.