A White-Box Machine Learning Approach for Revealing Antibiotic Mechanisms of Action.

Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated "white-box" biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.

Integrative omics indicate FMRP sequesters mRNA from translation and deadenylation in human neuronal cells.

How fragile X syndrome protein (FMRP) binds mRNAs and regulates mRNA metabolism remains unclear. Our previous work using human neuronal cells focused on mRNAs targeted for nonsense-mediated mRNA decay (NMD), which we showed are generally bound by FMRP and destabilized upon FMRP loss. Here, we identify >400 high-confidence FMRP-bound mRNAs, only ∼35% of which are NMD targets. Integrative transcriptomics together with SILAC-LC-MS/MS reveal that FMRP loss generally results in mRNA destabilization and more protein produced per FMRP target. We use our established RIP-seq technology to show that FMRP footprints are independent of protein-coding potential, target GC-rich and structured sequences, and are densest in 5' UTRs. Regardless of where within an mRNA FMRP binds, we find that FMRP protects mRNAs from deadenylation and directly binds the cytoplasmic poly(A)-binding protein. Our results reveal how FMRP sequesters polyadenylated mRNAs into stabilized and translationally repressed complexes, whose regulation is critical for neurogenesis and synaptic plasticity.

Multi-omics analysis reveals the potential pathogenesis and therapeutic targets of diabetic kidney disease.

Clinicians have long been interested in understanding the molecular basis of diabetic kidney disease (DKD)and its potential treatment targets. Its pathophysiology involves protein phosphorylation, one of the most recognizable post-transcriptional modifications, that can take part in many cellular functions and control different metabolic processes. In order to recognize the molecular and protein changes of DKD kidney, this study applied Tandem liquid chromatography-mass spectrometry (LC-MS/MS) and Next-Generation Sequencing, along with Tandem Mass Tags (TMT) labeling techniques to evaluate the mRNA, protein and modified phosphorylation sites between DKD mice and model ones. Based on Gene Ontology (GO) and KEGG pathway analyses of transcriptome and proteome, The molecular changes of DKD include accumulation of extracellular matrix, abnormally activated inflammatory microenvironment, oxidative stress and lipid metabolism disorders, leading to glomerulosclerosis and tubulointerstitial fibrosis. Oxidative stress has been emphasized as an important factor in DKD and progression to ESKD, which is directly related to podocyte injury, albuminuria and renal tubulointerstitial fibrosis. A histological study of phosphorylation further revealed that kinases were crucial. Three groups of studies have found that RAS signaling pathway, RAP1 signaling pathway, AMPK signaling pathway, PPAR signaling pathway and HIF-1 signaling pathway were crucial for the pathogenesis of DKD. Through this approach, it was discovered that targeting specific molecules, proteins, kinases and critical pathways could be a promising approach for treating DKD.

Systematic Proteome Profiling of Maternal Plasma for Development of Preeclampsia Biomarkers.

Preeclampsia (PE) is a hypertensive disorder of pregnancy with various clinical symptoms. However, traditional markers for the disease including high blood pressure and proteinuria are poor indicators of the related adverse outcomes. Here, we performed systematic proteome profiling of plasma samples obtained from pregnant women with PE to identify clinically effective diagnostic biomarkers. Proteome profiling was performed using TMT-based liquid chromatography-mass spectrometry (LC-MS/MS) followed by subsequent verification by multiple reaction monitoring (MRM) analysis on normal and PE maternal plasma samples. Functional annotations of differentially expressed proteins (DEPs) in PE were predicted using bioinformatic tools. The diagnostic accuracies of the biomarkers for PE were estimated according to the area under the receiver-operating characteristics curve (AUC). A total of 1307 proteins were identified, and 870 proteins of them were quantified from plasma samples. Significant differences were evident in 138 DEPs, including 71 upregulated DEPs and 67 downregulated DEPs in the PE group, compared with those in the control group. Upregulated proteins were significantly associated with biological processes including platelet degranulation, proteolysis, lipoprotein metabolism, and cholesterol efflux. Biological processes including blood coagulation and acute-phase response were enriched for down-regulated proteins. Of these, 40 proteins were subsequently validated in an independent cohort of 26 PE patients and 29 healthy controls. APOM, LCN2, and QSOX1 showed high diagnostic accuracies for PE detection (AUC >0.9 and p < 0.001, for all) as validated by MRM and ELISA. Our data demonstrate that three plasma biomarkers, identified by systematic proteomic profiling, present a possibility for the assessment of PE, independent of the clinical characteristics of pregnant women.

Galectin-3-Binding Protein Inhibits Extracellular Heparan 6-O-Endosulfatase Sulf-2.

Human extracellular 6-O-endosulfatases Sulf-1 and Sulf-2 are the only enzymes that post-synthetically alter the 6-O sulfation of heparan sulfate proteoglycans (HSPG), which regulates interactions of HSPG with many proteins. Oncogenicity of Sulf-2 in different cancers has been documented, and we have shown that Sulf-2 is associated with poor survival outcomes in head and neck squamous cell carcinoma (HNSCC). Despite its importance, limited information is available on direct protein-protein interactions of the Sulf-2 protein in the tumor microenvironment. In this study, we used monoclonal antibody (mAb) affinity purification and mass spectrometry to identify galectin-3-binding protein (LG3BP) as a highly specific binding partner of Sulf-2 in the conditioned media of HNSCC cell lines. We validated their direct interaction in vitro using recombinant proteins and have shown that the chondroitin sulfate (CS) covalently bound to the Sulf-2 influences the binding to LG3BP. We confirmed the importance of the CS chain for the interaction by generating a mutant Sulf-2 protein that lacks the CS. Importantly, we have shown that the LG3BP inhibits Sulf-2 activity in vitro in a concentration-dependent manner. As a consequence, the addition of LG3BP to a spheroid cell culture inhibited the invasion of the HNSCC cells into Matrigel. Thus, Sulf-2 interaction with LG3BP may regulate the physiological activity of the Sulf-2 enzyme as well as its activity in the tumor microenvironment.

Direct sequencing of total Saccharomyces cerevisiae tRNAs by LC-MS/MS.

Among RNAs, transfer RNAs (tRNAs) contain the widest variety of abundant posttranscriptional chemical modifications. These modifications are crucial for tRNAs to participate in protein synthesis, promoting proper tRNA structure and aminoacylation, facilitating anticodon:codon recognition, and ensuring the reading frame maintenance of the ribosome. While tRNA modifications were long thought to be stoichiometric, it is becoming increasingly apparent that these modifications can change dynamically in response to the cellular environment. The ability to broadly characterize the fluctuating tRNA modification landscape will be essential for establishing the molecular level contributions of individual sites of tRNA modification. The locations of modifications within individual tRNA sequences can be mapped using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). In this approach, a single tRNA species is purified, treated with ribonucleases, and the resulting single-stranded RNA products are subject to LC-MS/MS analysis. The application of LC-MS/MS to study tRNAs is limited by the necessity of analyzing one tRNA at a time, because the digestion of total tRNA mixtures by commercially available ribonucleases produces many short digestion products unable to be uniquely mapped back to a single site within a tRNA. We overcame these limitations by taking advantage of the highly structured nature of tRNAs to prevent the full digestion by single-stranded RNA-specific ribonucleases. Folding total tRNA prior to digestion allowed us to sequence Saccharomyces cerevisiae tRNAs with up to 97% sequence coverage for individual tRNA species by LC-MS/MS. This method presents a robust avenue for directly detecting the distribution of modifications in total tRNAs.

O-Glycomic and Proteomic Signatures of Spontaneous and Butyrate-Stimulated Colorectal Cancer Cell Line Differentiation.

Gut microbiota of the gastrointestinal tract provide health benefits to the human host via bacterial metabolites. Bacterial butyrate has beneficial effects on intestinal homeostasis and is the preferred energy source of intestinal epithelial cells, capable of inducing differentiation. It was previously observed that changes in the expression of specific proteins as well as protein glycosylation occur with differentiation. In this study, specific mucin O-glycans were identified that mark butyrate-induced epithelial differentiation of the intestinal cell line CaCo-2 (Cancer Coli-2), by applying porous graphitized carbon nano-liquid chromatography with electrospray ionization tandem mass spectrometry. Moreover, a quantitative proteomic approach was used to decipher changes in the cell proteome. It was found that the fully differentiated butyrate-stimulated cells are characterized by a higher expression of sialylated O-glycan structures, whereas fucosylation is downregulated with differentiation. By performing an integrative approach, we generated hypotheses about the origin of the observed O-glycome changes. These insights pave the way for future endeavors to study the dynamic O-glycosylation patterns in the gut, either produced via cellular biosynthesis or through the action of bacterial glycosidases as well as the functional role of these patterns in homeostasis and dysbiosis at the gut-microbiota interface.

Identification of a rare [Gγ(Aγδß)0] -thalassemia using tandem mass spectrometry.

Thalassemias are a group of inherited monogenic disorders characterized by defects in the synthesis of one or more of the globin chain subunits of the hemoglobin tetramer. Delta-beta (δß-) thalassemia has large deletions in the ß globin gene cluster involving δ- and ß-globin genes, leading to absent or reduced synthesis of both δ- and ß-globin chains. Here, we used direct globin-chain analysis using tandem mass spectrometry for the diagnosis of δß-thalassemia. Two cases from unrelated families were recruited for the study based on clinical and hematological evaluation. Peptides obtained after trypsin digestion of proteins extracted from red blood cell pellets from two affected individuals and their parents were analyzed using liquid chromatography-tandem mass spectrometry (LC-MS/MS). Mass spectrometric analysis revealed a severe reduction in δ, ß, and Aγ globin proteins with increased Gγ globin protein in the affected individuals. The diagnosis of Gγ(Aγδß)0 -thalassemia in the homozygous state in the affected individuals and in the heterozygous state in the parents was made from our results. The diagnosis was confirmed at the genetic level using multiplex ligation-dependent probe amplification (MLPA). Our findings demonstrate the utility of direct globin protein quantitation using LC-MS/MS to quantify individual globin proteins reflecting changes in globin production. This approach can be utilized for accurate and timely diagnosis of hemoglobinopathies, including rare variants, where existing diagnostic methods provide inconclusive results.

On the excessive use of coefficient of variation as a metric of quantitation quality in proteomics.

The coefficient of variation (CV) is often used in proteomics as a proxy to characterize the performance of a quantitation method and/or the related software. In this note, we question the excessive reliance on this metric in quantitative proteomics that may result in erroneous conclusions. We support this note using a ground-truth Human-Yeast-E. coli dataset demonstrating in a number of cases that erroneous data processing methods may lead to a low CV which has nothing to do with these methods' performances in quantitation.

Global profiling of protein lactylation in Caenorhabditis elegans.

Lactylation, as a novel posttranslational modification, is essential for studying the functions and regulation of proteins in physiological and pathological processes, as well as for gaining in-depth knowledge on the occurrence and development of many diseases, including tumors. However, few studies have examined the protein lactylation of one whole organism. Thus, we studied the lactylation of global proteins in Caenorhabditis elegans to obtain an in vivo lactylome. Using an MS-based platform, we identified 1836 Class I (localization probabilities > 0.75) lactylated sites in 487 proteins. Bioinformatics analysis showed that lactylated proteins were mainly located in the cytoplasm and involved in the tricarboxylic acid cycle (TCA cycle) and other metabolic pathways. Then, we evaluated the conservation of lactylation in different organisms. In total, 41 C. elegans proteins were lactylated and homologous to lactylated proteins in humans and rats. Moreover, lactylation on H4K80 was conserved in three species. An additional 238 lactylated proteins were identified in C. elegans for the first time. This study establishes the first lactylome database in C. elegans and provides a basis for studying the role of lactylation.