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.

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.

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.

π-π 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.

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.

Autoinhibitory mechanism for the mutation-induced impaired FGF9 signaling.

Fibroblast growth factor 9 (FGF9), an important member of the fibroblast growth factor (FGF) family, can bind with high affinity to FGFR3 in a heparin-dependent approach. In humans, the deletions and mutations resulting in dysfunction of the FGF9 signaling can cause human skeletal dysplasia and cancers. A mutation (S99N) in this protein has been identified to be associated with significantly impaired FGF signaling considered as a potential cause of synostoses syndrome. However, the detailed mechanism for this observation still remains unknown. In this study, we used molecular dynamics simulations and free energy calculations to study the interactions of FGF9(WT/S99N), FGFR3c, and heparin, with an aim of providing atomic sights into the detailed mechanism for the impaired FGF signaling caused by the S99N mutation. We found that the S99N mutation has a well-ordered C-terminal structure, which can reduce its homodimerization ability so as to break the monomer-dimer equilibrium in the FGF signaling, which is considered as a key factor to regulate extracellular matrix affinity and tissue diffusion in the FGF signaling pathway. The FGF9(WT) monomer can preferentially form a homodimer owing to its comparatively favorable binding free energy. In contrast, the FGF9(S99N) monomer is preferred to bind with the FGFR3c receptor to form an inactive complex, leading to impair FGF signaling. To support our computational findings, we also performed biochemical experiments, which confirm the computational results mentioned above. The impaired FGF signaling is believed to be a potential cause of human synostoses syndrome, implicating an important role for FGF9 in normal joint development.

A negative cooperativity mechanism of human CYP2E1 inferred from molecular dynamics simulations and free energy calculations.

Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. This residue is a key factor in the positioning of substrates and in proton delivery to the active site. Additionally, a long-range effect of the effector substrate is identified in which F478 is proposed to play a key role. As located in the interface between the active and effector sites, this residue structurally links the active and effector sites and is found to play a significant role in affecting substrate access and ligand positioning within the active site. In the negative cooperative binding, this residue can decrease the interactions of the active substrate with the active site by π-π stacking which then lowers the hydroxylation activity for the active substrate. These findings are in agreement with previous experimental observations and thus provide detailed atomistic insight into the poorly understood mechanism of the negative cooperativity in human CYP2E1.

Molecular dynamics studies on T1 lipase: insight into a double-flap mechanism.

T1 lipase is isolated from the palm Geobacillus zalihae strain T1 in Malaysia, functioning as a secreted protein responsible for the catalyzing hydrolysis of long-chain triglycerides into fatty acids and glycerol at high temperatures. In the current study, using 30 ns molecular dynamics simulations at different temperatures, an aqueous activation was detected for T1 lipase. This aqueous activation in T1 lipase was mainly caused by a double-flap movement mechanism. The double flaps were constituted by the hydrophobic helices 6 and 9. Helix 6 employed two major components with the hydrophilic part at the surface and the hydrophobic part inside. In the aqueous solution, the hydrophobic part could provide enough power for helix 6 to move away, driving the protein into an open configuration and exposing the catalytic triad. Our findings could provide structural evidence to support the double-flap movement, revealing the catalytic mechanism for T1 lipase.

Structure of cytochrome p450s and personalized drug.

Cytochrome P450s are the most important enzymes responsible for phase I drug metabolism. The polymorphic nature of cytochrome P450s largely influences individual drug responses, drug-drug interactions and induces adverse drug reactions. By far, thirty crystal structures of eight mammalian cytochrome P450s (CYP 2C5, 2C8, 2C9, 3A4, 2D6, 2B4, 2A6 and 1A2) have been published. This review focuses on the recent studies on the structures of cytochrome P450s: some characteristic features of these enzymes and many essential, conserved amino acids in the active sites have been identified. These results are of fundamental importance for drug development and understanding the metabolism for both endogenous and xenobiotic substrates. With the help of computational methods, the structural information will provide insights into personalization of drug treatments in both proper drug therapy and appropriate dosage of a certain drug.

Insight into the molecular switch mechanism of human Rab5a from molecular dynamics simulations.

Rab5a is currently a most interesting target because it is responsible for regulating the early endosome fusion in endocytosis and possibly the budding process. We utilized longtime-scale molecular dynamics simulations to investigate the internal motion of the wild-type Rab5a and its A30P mutant. It was observed that, after binding with GTP, the global flexibility of the two proteins is increasing, while the local flexibility in their sensitive sites (P-loop, switch I and II regions) is decreasing. Also, the mutation of Ala30 to Pro30 can cause notable flexibility variations in the sensitive sites. However, this kind of variations is dramatically reduced after binding with GTP. Such a remarkable feature is mainly caused by the water network rearrangements in the sensitive sites. These findings might be of use for revealing the profound mechanism of the displacements of Rab5a switch regions, as well as the mechanism of the GDP dissociation and GTP association.

Insights from investigating the interactions of adamantane-based drugs with the M2 proton channel from the H1N1 swine virus.

The M2 proton channel is one of indispensable components for the influenza A virus that plays a vital role in its life cycle and hence is an important target for drug design against the virus. In view of this, the three-dimensional structure of the H1N1-M2 channel was developed based on the primary sequence taken from a patient recently infected by the H1N1 (swine flu) virus. With an explicit water-membrane environment, molecular docking studies were performed for amantadine and rimantadine, the two commercial drugs generally used to treat influenza A infection. It was found that their binding affinity to the H1N1-M2 channel is significantly lower than that to the H5N1-M2 channel, fully consistent with the recent report that the H1N1 swine virus was resistant to the two drugs. The findings and the relevant analysis reported here might provide useful structural insights for developing effective drugs against the new swine flu virus.

Reversal of coenzyme specificity and improvement of catalytic efficiency of Pichia stipitis xylose reductase by rational site-directed mutagenesis.

A major problem when xylose is used for ethanol production is the intercellular redox imbalance arising from different coenzyme specificities of xylose reductase (XR) and xylitol dehydrogenase. The residue Lys21 in XR from Pichia stipitis was subjected to site-directed mutagenesis to alter its coenzyme specificity. The N272D mutant exhibited improved catalytic efficiency when NADH was the coenzyme. Both K21A and K21A/N272D preferred NADH to NADPH, their catalytic efficiencies for NADPH were almost zero. The catalytic efficiency of K21A/N272D for NADH was almost 9-fold and 2-fold that of K21A and the wild-type enzyme, respectively. Complete reversal of coenzyme specificity toward NADH and improved catalytic efficiency were achieved.

Molecular modeling of two CYP2C19 SNPs and its implications for personalized drug design.

CYP2C19 is an important member of the cytochrome P-450 enzyme superfamily and plays a significant role in the drug metabolism. In order to gain insights for developing personalized drugs, the structure-activity relationships of two SNPs, W120R and I331V, with the ligands of CEC, Fluvoxamine, Lescol and Ticlopidine were investigated through the structure-activity relationship approach. By means of a series of docking studies, the binding pockets of the two SNPs for the four compounds are explicitly defined that will be very useful for conducting mutagenesis studies, providing insights into personalization of drug treatments and stimulating novel strategies for finding desired personalized drugs.

[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.

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.

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.

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.