KPP: KEGG Pathway Painter

An integrated investigation of structural and pathway alteration caused by PIK3CA and TP53 mutations identified in cfDNA of metastatic breast cancer

In peripheral blood, cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA), which indicates molecular abnormalities in metastatic breast tumor tissue. The sequencing of cfDNA of Metastatic Breast Cancer (MBC) patients allows assessment of therapy response and noninvasive treatment. In the proposed study, clinically significant alterations in PIK3CA and TP53 genes associated with MBC resulting in a missense substitution of His1047Arg and Arg282Trp from an next-generation sequencing-based multi-gene panel were reported in a cfDNA of a patient with MBC. To investigate the impact of the reported mutation, we used molecular docking, molecular dynamics simulation, network analysis, and pathway analysis. Molecular Docking analysis determined the distinct binding pattern revealing H1047R-ATP complex has a higher number of Hydrogen bonds (H-bonds) and binding affinity with a slight difference compared to the PIK3CA-ATP complex. Following, molecular dynamics simulation for 200 ns, of which H1047R-ATP complex resulted in the instability of PIK3CA. Similarly, for TP53 mutant R282W, the zinc-free state (apo) and zinc-bounded (holo) complexes were investigated for conformational change between apo and holo complexes, of which the holo complex mutant R282W was unstable. To validate the conformational change of PIK3CA and TP53, 80% mutation of H1047R in the kinase domain of p110α expressed ubiquitously in PIK3CA protein that alters PI3K pathway, while R282W mutation in DNA binding helix (H2) region of P53 protein inhibits the transcription factor in P53 pathway causing MBC. According to our findings, the extrinsic (hypoxia, oxidative stress, and acidosis); intrinsic factors (MYC amplification) in PIK3CA and TP53 mutations will provide potential insights for developing novel therapeutic methods for MBC therapy.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

Testicular Localization and Potential Function of Vimentin Positive Cells during Spermatogonial Differentiation Stages

Simple Summary

Spermatogonia stem cells (SSCs) in the testis are responsible for transmitting genetic information to subsequent generations. During the differentiation of SSCs, various cytoskeletons are involved in giving the cell shape and internal organization. One of the essential cytoskeletons that play functions in these spermatogenic processes is vimentin. This study examined the vimentin expression in vivo and in vitro by immunocytochemistry (ICC), immunohistochemistry (IMH), Fluidigm real-time polymerase chain reaction, and bioinformatics analysis. IMH analysis demonstrated that the high vimentin expression was localized in the middle and central testicular cords cells of seminiferous tubules and low expression in the basal region under in vivo conditions. To evaluate the expression of vimentin in vitro, we first isolated SSCs and then cultured cells from the testis. Our results showed that vimentin plays important roles in the differentiation of testicular germ cells.

Abstract

Vimentin is a type of intermediate filament (IF) and one of the first filaments expressed in spermatogenesis. Vimentin plays numerous roles, consisting of the determination of cell shape, differentiation, cell motility, the maintenance of cell junctions, intracellular trafficking, and assisting in keeping normal differentiating germ cell morphology. This study investigated the vimentin expression in two populations of undifferentiated and differentiated spermatogonia. We examined vimentin expression in vivo and in vitro by immunocytochemistry (ICC), immunohistochemistry (IMH), and Fluidigm real-time polymerase chain reaction. IMH data showed that the high vimentin expression was localized in the middle of seminiferous tubules, and low expression was in the basal membrane. ICC analysis of the colonies by isolated differentiated spermatogonia indicated the positive expression for the vimentin antibody, but vimentin’s expression level in the undifferentiated population was negative under in vitro conditions. Fluidigm real-time PCR analysis showed significant vimentin expression in differentiated spermatogonia compared to undifferentiated spermatogonia (p < 0.05). Our results showed that vimentin is upregulated in the differentiation stages of spermatogenesis, proving that vimentin is an intermediate filament with crucial roles in the differentiation stages of testicular germ cells. These results support the advanced investigations of the spermatogenic process, both in vitro and in vivo.

Analysis of Microarray Data from Medulloblastoma Tissue Samples

As a laboratory tool, microarray is used to detect the expression of thousands of genes at the same time. Typically, microscope slides have DNA microarrays that are printed with thousands of tiny spots in specified positions. Each spot contains a known DNA sequence or gene. These slides are commonly referred to as gene chips or DNA chips. The DNA molecules printed to each slide serve as probes to detect gene expression, which is also known as the transcriptome or the set of messenger RNA (mRNA) transcripts expressed by a group of genes. The goal of this chapter is to discuss the steps involved computational analysis of data after the completion of a typical microarray experiment.

Plasma exosomal proteomic studies of corneal epithelial injury in diabetic and non-diabetic group

OBJECTIVE

Diabetic Keratopathy (DK) is one of the significant complications of type II diabetes (T2DM) with pathogenesis not yet clarified. Since hyperglycemia is able to change the protein components contained in plasma exosomes, liquid chromatography-tandem mass spectrometry (LC-MS/MS) is considered as feasible to analyze the expression of plasma exosomal proteins in patients with T2DM and non-diabetic patients respectively, find critical biological markers, and explore the mechanism of DK as well as potential therapeutic targets.

METHOD

Blood and clinical information of corneal epithelial injury in a diabetic group (the study group) and a non-diabetic group (the control group), who were patients admitted to the Department of Ophthalmology, Yangpu Hospital, Tongji University School of Medicine from July 2020 to November 2020, were collected. The qEV size exclusion method was adopted to separate exosomes from plasma. The exosomes were then identified through transmission electron microscopy (TEM), nanoparticle tracking analyzer (NTA), and Western blot. The plasma exosomes of the study group and the control group were quantitatively analyzed by proteomics. A bioinformatics method is utilized to screen differential proteins and the expression of the differential proteins was verified by Western blot.

RESULT

TEM indicated that the exosomes had a double-concave disc-like appearance, with a size of about 100 nm, and Western blot expressed as CD63 and TSG101. The plasma exosomes of the study group and the control group were analyzed by quantitative proteomics with a total number of 952 proteins detected of which 245 proteins existed in the ExoCarta exosomal protein database. Through adoption of P-value to screen credible differential proteins, the heat map displayed 28 differential proteins, 7 upregulated proteins, and 21 downregulated proteins; the volcano map displayed 7 upregulated proteins and 22 downregulated proteins; the PPI interaction map displayed 12 upregulated proteins and 18 downregulated proteins. Through GO enrichment analysis, it was identified that the differential protein participated in the main biological processes and was involved in regulating the cell's stimulation response to insulin, the insulin receptor signaling pathway, and the activity of glycosylphosphatidylinositol phospholipase D as well as anti-oxidation. The enriched cell components include main components such as exosomes, blood particles, and cytoplasm. KEGG enrichment analysis indicated that the target protein FLOT2 was mainly concentrated in insulin-related signaling pathways. Western blot indicated that the expression of FLOT2 in the study group was lower compared with the control group while the expression of Exo70 was higher.

CONCLUSION

Proteomic analysis of the study group and the control group displayed a variety of proteins in plasma exosomes. The downregulated protein FLOT2 in the study group was closely related to the occurrence, development, and complication of DK in T2DM patients. The expression status of plasma FLOT2 protein in T2DM patients is expected to be a biomarker for diagnosing and monitoring of DK.

Dysregulation of autophagy-associated microRNAs in condyloma acuminatum

Condyloma acuminatum, which is caused by low-risk human papillomavirus (lrHPV) infection, is one of the most common sexually transmitted diseases. Autophagy is thought to be associated with condyloma acuminatum, but how the autophagy process is regulated remains unclear. MicroRNAs (miRNAs) are important regulators of gene transcription that play a central role in many biological processes, including autophagy and viral infection. This study was designed to identify autophagy-related miRNAs and their targets in condyloma acuminatum and to validate their expression. The levels of the autophagy proteins microtubule-associated protein 1 light chain 3 (LC3) and P62/SQSTM1 (P62) were abnormally increased in the local lesion tissue of condyloma acuminatum patients compared with healthy controls. MiRNAs and their target mRNAs in condyloma acuminatum patients were analyzed by bioinformatics. Eighty-one differentially expressed miRNAs were identified, of which 56 were downregulated and 25 were upregulated. Two of the differentially expressed miRNAs associated with autophagy, miRNA-30a-5p and miRNA-514a-3p, were analyzed further, and their target genes were identified as autophagy-related protein (Atg) 5 and Atg12 and Atg3 and Atg12, respectively. The expression levels of miRNA-30a-5p and miRNA-514a-3p were decreased and those of Atg5, Atg12 and Atg3 were increased in condyloma acuminatum patients compared with healthy controls. In addition, miRNA-30a-5p and miRNA-514a-3p expression correlated with the proliferation index Ki-67 in condyloma acuminatum. Taken together, our results suggest that the changes in autophagy levels in patients with condyloma acuminatum may be related to the changes in miRNA-30a-5p and miRNA-514a-3p expression. This study provides a theoretical basis for identifying new mechanisms that link miRNAs, HPV infection and host autophagy in vivo.

iTRAQ-based proteomics analysis on insomnia rats treated with Mongolian medical warm acupuncture

Objective: To explore the proteomic changes in the hypothalamus of rats treated with Mongolian medical warm acupuncture for insomnia therapy based proteomics. Method: We used an iTRAQ-based quantitative proteomic approach to identify proteins that potential molecular mechanisms involved in the treatment of insomnia by Mongolian medical warm acupuncture. Result: In total, 7477 proteins were identified, of which 36 proteins showed increased levels and 45 proteins showed decreased levels in insomnia model group (M) compared with healthy control group (C), 72 proteins showed increased levels and 44 proteins showed decreased levels from the warm acupuncture treated insomnia group (W) compared with healthy controls (C), 28 proteins showed increased levels and 17 proteins showed decreased levels from the warm acupuncture-treated insomnia group (W) compared with insomnia model group (M). Compared with healthy control groups, warm acupuncture-treated insomnia group showed obvious recovered. Bioinformatics analysis indicated that up-regulation of neuroactive ligand–receptor interaction and oxytocin signaling was the most significantly elevated regulate process of Mongolian medical warm acupuncture treatment for insomnia. Proteins showed that increased/decreased expression in the warm acupuncture-treated insomnia group included Prolargin (PRELP), NMDA receptor synaptonuclear-signaling and neuronal migration factor (NSMF), Transmembrane protein 41B (TMEM41B) and Microtubule-associated protein 1B (MAP1B) to adjust insomnia. Conclusion: A combination of findings in the present study suggest that warm acupuncture treatment is efficacious in improving sleep by regulating the protein expression process in an experimental rat model and may be of potential benefit in treating insomnia patients with the added advantage with no adverse effects.

Common Genes Involved in Autophagy, Cellular Senescence and the Inflammatory Response in AMD and Drug Discovery Identified via Biomedical Databases

Purpose

Retinal pigment epithelial cell autophagy dysfunction, cellular senescence, and the retinal inflammatory response are key pathogenic factors in age-related macular degeneration (AMD), which has been reviewed in our previously work in 2019. This study aims to identify genes collectively involved in these three biological processes and target drugs in AMD.

Methods

The pubmed2ensembl database was used to perform text mining. The GeneCodis database was applied to analyze gene ontology biological process and the KEGG pathway. The STRING database was used to analyze protein–protein interaction analysis and hub genes were identified by the Cytoscape software. The Drug Gene Interaction Database was used to perform drug–gene interactions.

Results

We identified 62 genes collectively involved in AMD, autophagy, cellular senescence, and inflammatory response, 19 biological processes including 42 genes, 11 enriched KEGG pathways including 37 genes, and 12 hub genes step by step via the above biomedical databases. Finally, five hub genes (IL-6, VEGF-A, TP53, IL-1β, and transforming growth factor [TGF]-β1) and their specific interaction modes were identified, corresponding with 24 target drugs with therapeutic potential for AMD.

Conclusions

IL-6, VEGF-A, TP53, IL-1β, and TGF-β1 are pivotal in autophagy, cellular senescence, and the inflammatory response in AMD, corresponding with 24 drugs with therapeutic potential for AMD, providing definite molecular mechanisms for further research and new possibilities for AMD treatment in the future.

Translational Relevance

IL-6, VEGF-A, TP53, IL-1β, and TGF-β1 may be new targets for AMD gene therapy and drug development.

Identification of putative non-host essential genes and novel drug targets against Acinetobacter baumannii by in silico comparative genome analysis

Acinetobacter baumannii, the gram-negative bacteria emerged as an extremely critical pathogen causing nosocomial and different kinds of infections. A. baumannii exhibit resistivity towards various classes of antibiotics that shows that there is a dire need to search more drug targets by exploiting the full genome of the bacteria. In doing so, a strategy is made with the combination of computational biology, pathogen informatics and cheminformatics. Comparative genomics analysis, modeling and docking studies have been performed for the prediction of non-host essential genes and novel drug candidates against A. baumannii. Among 37 unique and 82 common metabolic pathways, 92 genes were predicted as non-host genes. Similarly, using homology search between A. baumannii genome and essential genes of different bacteria, 293 genes were predicted as essential genes of A. baumannii. Among these predicted non-host and essential genes, 86 genes were predicted as non-host essential genes which could serve as potential novel drug and vaccine targets. Additional drug-target like physicochemical properties were estimated such as the molecular weight, subcellular localization and druggability potential. On the structural part, the crystal structures of all the non-host essential genes of A. baumannii were found except the three genes. Out of these three, a homology model of Undecaprenyl-diphosphatase was built using a PDB template by MODELLER [version 9.18]. The quality of the model was assessed by the ProSA and RAMPAGE. The built model was subjected as a receptor for the molecular docking with Adenosine diphosphate (ADP) as a ligand. The molecular docking was performed by AutoDock4 and the best conformation with lowest binding energy (-4.39 kcal/mol) was obtained. The LigPlot was used to identify the close interactions between the ligand the receptor's residues. This study will further aid for the selection of putative inhibitors against a novel drug target identified against A. baumannii and hence could lead to the better therapeutics.

Influence of static forces on the expression of selected parameters of inflammation in periodontal ligament cells and alveolar bone cells in a co-culture in vitro model

OBJECTIVE

Aim of this study was to investigate the impact of human PDL-derived fibroblasts (HPDF) and human alveolar bone-derived osteoblasts (HABO) co-culture on the expression of cytokines involved in tissue remodeling using an in vitro compressive force (CF) model.

MATERIALS AND METHODS

Static compressive force (CF) of 47.4 g/cm² was applied on mono- and co-cultured HPDFs and HABOs for 1, 2, or 4 h at 30 °C. TNFA, PTGS2, and IL6 gene expressions were determined by quantitative real-time polymerase chain reaction. TNF, PGE₂, and IL6 concentrations were measured using enzyme-linked immunosorbent assay.

RESULTS

In mono-culture, TNFA, PTGS2, and IL6 gene expressions were upregulated under CF as compared to controls for each time period in both cell types. PGE₂ increased at 1 and 2 h in both cell types, and IL6 increased only at 2 and 4 h in HPDFs. Co-culture alleviated the force-induced increase of the expression of TNFA, PTGS2, IL6, PGE2, and IL6 in HPDFs at any time point. In HABOs, co-cultivation decreased the expression of PGE₂ after 1 h and 4 h, and that of IL6 after 1 h compared to mono-culture.

CONCLUSIONS

CF application on co-cultures of HPDFs and HABOs causes significant changes of TNFA, PTGS2, and IL6 gene expressions and PGE₂ and IL6 production in comparison to mono-culture indicating intercellular communication.

CLINICAL RELEVANCE

Mechanical stimulation of HPDFs and HABOs in co-culture induces a different gene expression pattern than stimulation of individual cell types alone. Co-culture might therefore be a relevant method to elucidate periodontal regeneration during orthodontic therapy.

An exploratory study examining how nano-liquid chromatography-mass spectrometry and phosphoproteomics can differentiate patients with advanced fibrosis and higher percentage collagen in non-alcoholic fatty liver disease

BACKGROUND

Non-alcoholic steatohepatitis (NASH) is among the leading causes of liver disease worldwide. It is increasingly recognized that the phenotype of NASH may involve a number of different pathways, of which each could become important therapeutic targets. The aim of this study is to use high resolution mass spectrometry (MS) and phosphoproteomics techniques to assess the serum proteome and hepatic phosphoproteome in subjects with NASH-related fibrosis.

METHODS

Sixty-seven biopsy-proven NAFLD subjects with frozen sera and liver tissue were included. Reverse phase protein microarray was used to quantify the phosphorylation of key signaling proteins in liver and nano-liquid chromatography (LC)-MS was used to sequence target biomarkers in the serum. An image analysis algorithm was used to quantify the percentage of collagen (% collagen) using computer-assisted morphometry. Using multiple regression models, serum proteomes and phosphorylated hepatic proteins that were independently (p ≤ 0.05) associated with advanced fibrosis (stage ≥ 2) and higher % collagen were assessed.

RESULTS

Phosphorylated signaling pathways in the liver revealed that apoptosis signal-regulating kinase 1, mitogen-activated protein kinase (ASK1-MAPK pathway involving ASK1 S38 (p < 0.02) and p38 MAPK (p = 0.0002)) activated by the inflammatory cytokine interleukin (IL-10) (p < 0.001), were independently associated with higher % collagen. LC-MS data revealed that serum alpha-2 macroglobulin (α2M) (p = 0.0004) and coagulation factor V (p = 0.0127) were independently associated with higher % hepatic collagen.

CONCLUSIONS

Simultaneous profiling of serum proteome and hepatic phosphoproteome reveals that the activation of ASK1 S38, p38 MAPK in the liver, and serum α2M and coagulation factor V are independently associated with hepatic collagen deposition in patients with NASH. These data suggest the role of these pathways in the pathogenesis of NASH-related fibrosis as a potential therapeutic target.

E47 Governs the MYC-CDKN1B/p27KIP1-RB Network to Growth Arrest PDA Cells Independent of CDKN2A/p16INK4A and Wild-Type p53

BACKGROUND & AIMS

Oncogenic mutations in KRAS, coupled with inactivation of p53, CDKN2A/p16INK4A, and SMAD4, drive progression of pancreatic ductal adenocarcinoma (PDA). Overexpression of MYC and deregulation of retinoblastoma (RB) further promote cell proliferation and make identifying a means to therapeutically alter cell-cycle control pathways in PDA a significant challenge. We previously showed that the basic helix-loop-helix transcription factor E47 induced stable growth arrest in PDA cells in vitro and in vivo. Here, we identified molecular mechanisms that underlie E47-induced growth arrest in low-passage, patient-derived primary and established PDA cell lines.

METHODS

RNA sequencing was used to profile E47-dependent transcriptomes in 5 PDA cell lines. Gene Ontology analysis identified cell-cycle control as the most altered pathway. Small interfering RNA/short hairpin RNA knockdown, small-molecule inhibitors, and viral expression were used to examine the function of E47-dependent genes in cell-cycle arrest. Cell morphology, expression of molecular markers, and senescence-associated β-galactosidase activity assays identified cellular senescence.

RESULTS

E47 uniformly inhibited PDA cell-cycle progression by decreasing expression of MYC, increasing the level of CDKN1B/p27KIP1, and restoring RB tumor-suppressor function. The molecular mechanisms by which E47 elicited these changes included altering both RNA transcript levels and protein stability of MYC and CDKN1B/p27KIP1. At the cellular level, E47 elicited a senescence-like phenotype characterized by increased senescence-associated β-galactosidase activity and altered expression of senescence markers.

CONCLUSIONS

E47 governs a highly conserved network of cell-cycle control genes, including MYC, CDKN1B/p27KIP1, and RB, which can induce a senescence-like program in PDA cells that lack CDKN2A/p16INK4A and wild-type p53. RNA sequencing data are available at the National Center for Biotechnology Information GEO at https://www.ncbi.nlm.nih.gov/geo/; accession number: GSE100327.

Multiple Acquisitions of Pathogen-Derived Francisella Endosymbionts in Soft Ticks

Bacterial endosymbionts of ticks are of interest due to their close evolutionary relationships with tick-vectored pathogens. For instance, whereas many ticks contain Francisella-like endosymbionts (FLEs), others transmit the mammalian pathogen Francisella tularensis. We recently sequenced the genome of an FLE present in the hard tick Amblyomma maculatum (FLE-Am) and showed that it likely evolved from a pathogenic ancestor. In order to expand our understanding of FLEs, in the current study we sequenced the genome of an FLE in the soft tick Ornithodoros moubata and compared it to the genomes of FLE-Am, Francisella persica-an FLE in the soft tick Argus (Persicargas) arboreus, Francisella sp. MA067296-a clinical isolate responsible for an opportunistic human infection, and F. tularensis, the established human pathogen. We determined that FLEs and MA067296 belonged to a sister taxon of mammalian pathogens, and contained inactivated versions of virulence genes present in F. tularensis, indicating that the most recent common ancestor shared by FLEs and F. tularensis was a potential mammalian pathogen. Our analyses also revealed that the two soft ticks (O. moubata and A. arboreus) probably acquired their FLEs separately, suggesting that the virulence attenuation observed in FLEs are not the consequence of a single acquisition event followed by speciation, but probably due to independent transitions of pathogenic francisellae into nonpathogenic FLEs within separate tick lineages. Additionally, we show that FLEs encode intact pathways for the production of several B vitamins and cofactors, denoting that they could function as nutrient-provisioning endosymbionts in ticks.

Epigenomic annotation of noncoding mutations identifies mutated pathways in primary liver cancer

Evidence that noncoding mutation can result in cancer driver events is mounting. However, it is more difficult to assign molecular biological consequences to noncoding mutations than to coding mutations, and a typical cancer genome contains many more noncoding mutations than protein-coding mutations. Accordingly, parsing functional noncoding mutation signal from noise remains an important challenge. Here we use an empirical approach to identify putatively functional noncoding somatic single nucleotide variants (SNVs) from liver cancer genomes. Annotation of candidate variants by publicly available epigenome datasets finds that 40.5% of SNVs fall in regulatory elements. When assigned to specific regulatory elements, we find that the distribution of regulatory element mutation mirrors that of nonsynonymous coding mutation, where few regulatory elements are recurrently mutated in a patient population but many are singly mutated. We find potential gain-of-binding site events among candidate SNVs, suggesting a mechanism of action for these variants. When aggregating noncoding somatic mutation in promoters, we find that genes in the ERBB signaling and MAPK signaling pathways are significantly enriched for promoter mutations. Altogether, our results suggest that functional somatic SNVs in cancer are sporadic, but occasionally occur in regulatory elements and may affect phenotype by creating binding sites for transcriptional regulators. Accordingly, we propose that noncoding mutation should be formally accounted for when determining gene- and pathway-mutation burden in cancer.

The 'Omics Revolution in CHO Biology: Roadmap to Improved CHO Productivity

Increased understanding of Chinese hamster ovary (CHO) cell physiology has been ushered in upon availability of the parental CHO-K1 cell line genome. Free and openly accessible sequence information has complemented transcriptomic and proteomic studies. The previous decade has also seen an increase in sensitivity and accuracy of proteomic methods due to technology development. In this genomic era, high-throughput screening methods, sophisticated informatic tools, and models continually drive major innovations in cell line development and process engineering. This review describes the various achievements in 'omics techniques and their application to improve recombinant protein expression from CHO cell lines.

RNA-seq transcriptome analysis of a Pseudomonas strain with diversified catalytic properties growth under different culture medium

Biocatalysis is an emerging strategy for the production of enantio-pure organic molecules. However, lacking of commercially available enzymes restricts the widespread application of biocatalysis. In this study, we report a Pseudomonas strain which exhibited versatile oxidation activity to synthesize chiral sulfoxides when growing under M9-toluene medium and reduction activity to synthesize chiral alcohols when on Luria-Bertani (LB) medium, respectively. Further comparative transcriptome analysis on samples from these two cultural conditions has identified 1038 differentially expressed genes (DEG). Gene Ontology (GO) enrichment and KEGG pathways analysis demonstrate significant changes in protein synthesis, energy metabolism, and biosynthesis of metabolites when cells cultured under different conditions. We have identified eight candidate enzymes from this bacterial which may have the potential to be used for synthesis of chiral alcohol and sulfoxide chemicals. This work provides insights into the mechanism of diversity in catalytic properties of this Pseudomonas strain growth with different cultural conditions, as well as candidate enzymes for further biocatalysis of enantiomerically pure molecules and pharmaceuticals.