Time-dependent network analysis reveals molecular targets underlying the development of diet-induced obesity and non-alcoholic steatohepatitis

Prolonged high-fat diet leads to the development of obesity and multiple comorbidities including non-alcoholic steatohepatitis (NASH), but the underlying molecular basis is not fully understood. We combine molecular networks and time course gene expression profiles to reveal the dynamic changes in molecular networks underlying diet-induced obesity and NASH. We also identify hub genes associated with the development of NASH. Core diet-induced obesity networks were constructed using Ingenuity pathway analysis (IPA) based on 332 high-fat diet responsive genes identified in liver by time course microarray analysis (8 time points over 24 weeks) of high-fat diet-fed mice compared to normal diet-fed mice. IPA identified five core diet-induced obesity networks with time-dependent gene expression changes in liver. These networks were associated with cell-to-cell signaling and interaction (Network 1), lipid metabolism (Network 2), hepatic system disease (Network 3 and 5), and inflammatory response (Network 4). When we merged these core diet-induced obesity networks, Tlr2, Cd14, and Ccnd1 emerged as hub genes associated with both liver steatosis and inflammation and were altered in a time-dependent manner. Further, protein–protein interaction network analysis revealed Tlr2, Cd14, and Ccnd1 were interrelated through the ErbB/insulin signaling pathway. Dynamic changes occur in molecular networks underlying diet-induced obesity. Tlr2, Cd14, and Ccnd1 appear to be hub genes integrating molecular interactions associated with the development of NASH. Therapeutics targeting hub genes and core diet-induced obesity networks may help ameliorate diet-induced obesity and NASH.     

Roles of hydration water molecules in molecular packing of the killer toxin from Pichia farinosa in its crystalline state investigated by cryogenic X-ray crystallography

The hydration structures around the killer toxin from Pichia farinosa were investigated by cryogenic X-ray crystallography. In particular, those contributing to the molecular association and the crystal contacts were analyzed with respect to the geometry and the networks of hydrogen bonds. The hydration water molecules attached on the surface so as to make up the surface shape in the contact complementary and mediated the intermolecular interactions through the networks of hydrogen bonds. Careful inspection of the contact area led to a proposal as to the molecular association mode of the toxin to determine the biological function in cells. In addition, the water-associated protein-protein interactions were approximated well by a simple theoretical equation on the solvation force expected in confined geometry. The present analysis may provide a way to analyze the crystal contact and molecular recognition in macromolecules in aqueous solution.     

A simplified procedure for semi-targeted lipidomic analysis of oxidized phosphatidylcholines induced by UVA irradiation

Oxidized phospholipids (OxPLs) are increasingly recognized as signaling mediators that are not only markers of oxidative stress but are also "makers" of pathology relevant to disease pathogenesis. Understanding the biological role of individual molecular species of OxPLs requires the knowledge of their concentration kinetics in cells and tissues. In this work, we describe a straightforward "fingerprinting" procedure for analysis of a broad spectrum of molecular species generated by oxidation of the four most abundant species of polyunsaturated phosphatidylcholines (OxPCs). The approach is based on liquid-liquid extraction followed by reversed-phase HPLC coupled to electrospray ionization MS/MS. More than 500 peaks corresponding in retention properties to polar and oxidized PCs were detected within 8 min at 99 m/z precursor values using a single diagnostic product ion in extracts from human dermal fibroblasts. Two hundred seventeen of these peaks were fluence-dependently and statistically significantly increased upon exposure of cells to UVA irradiation, suggesting that these are genuine oxidized or oxidatively fragmented species. This method of semitargeted lipidomic analysis may serve as a simple first step for characterization of specific "signatures" of OxPCs produced by different types of oxidative stress in order to select the most informative peaks for identification of their molecular structure and biological role.     

Computational modeling of cellular signaling processes embedded into dynamic spatial contexts

Cellular signaling processes depend on spatiotemporal distributions of molecular components. Multicolor, high-resolution microscopy permits detailed assessment of such distributions, providing input for fine-grained computational models that explore mechanisms governing dynamic assembly of multimolecular complexes and their role in shaping cellular behavior. However, it is challenging to incorporate into such models both complex molecular reaction cascades and the spatial localization of signaling components in dynamic cellular morphologies. Here we introduce an approach to address these challenges by automatically generating computational representations of complex reaction networks based on simple bimolecular interaction rules embedded into detailed, adaptive models of cellular morphology. Using examples of receptor-mediated cellular adhesion and signal-induced localized mitogen-activated protein kinase (MAPK) activation in yeast, we illustrate the capacity of this simulation technique to provide insights into cell biological processes. The modeling algorithms, implemented in a new version of the Simmune toolset, are accessible through intuitive graphical interfaces and programming libraries.     

Proteomic patterns for classification of ovarian cancer and CTCL serum samples utilizing peak pairs indicative of post-translational modifications

Proteomic patterns as a potential diagnostic technology has been well established for several cancer conditions and other diseases. The use of machine learning techniques such as decision trees, neural networks, genetic algorithms, and other methods has been the basis for pattern determination. Cancer is known to involve signaling pathways that are regulated through PTM of proteins. These modifications are also detectable with high confidence using high-resolution MS. We generated data using a prOTOF mass spectrometer on two sets of patient samples: ovarian cancer and cutaneous t-cell lymphoma (CTCL) with matched normal samples for each disease. Using the knowledge of mass shifts caused by common modifications, we built models using peak pairs and compared this to a conventional technique using individual peaks. The results for each disease showed that a small number of peak pairs gave classification equal to or better than the conventional technique that used multiple individual peaks. This simple peak picking technique could be used to guide identification of important peak pairs involved in the disease process.     

Analysis of MALDI-TOF mass spectrometry data for discovery of peptide and glycan biomarkers of hepatocellular carcinoma

This paper presents computational methods to analyze MALDI-TOF mass spectrometry data for quantitative comparison of peptides and glycans in serum. The methods are applied to identify candidate biomarkers in serum samples of 203 participants from Egypt; 73 hepatocellular carcinoma (HCC) cases, 52 patients with chronic liver disease (CLD) consisting of cirrhosis and fibrosis cases, and 78 population controls. Two complementary sample preparation methods were applied prior to generating mass spectra: (1) low molecular weight (LMW) enrichment of each serum sample was carried out for MALDI-TOF quantification of peptides, and (2) glycans were enzymatically released from proteins in each serum sample and permethylated for MALDI-TOF quantification of glycans. A peak selection algorithm was applied to identify the most useful peptide and glycan peaks for accurate detection of HCC cases from high-risk population of patients with CLD. In addition to global peaks selected by the whole population based approach, where identically labeled patients are treated as a single group, subgroup-specific peaks were identified by searching for peaks that are differentially abundant in a subgroup of patients only. The peak selection process was preceded by peak screening, where we eliminated peaks that have significant association with covariates such as age, gender, and viral infection based on the peptide and glycan spectra from population controls. The performance of the selected peptide and glycan peaks was evaluated in terms of their ability in detecting HCC cases from patients with CLD in a blinded validation set and through the cross-validation method. Finally, we investigated the possibility of using both peptides and glycans in a panel to enhance the diagnostic capability of these candidate markers. Further evaluation is needed to examine the potential clinical utility of the candidate peptide and glycan markers identified in this study.     

Combined frequency- and time-domain NMR spectroscopy. Application to fast protein resonance assignment

A simple and general method is presented to simplify multi-dimensional NMR spectra of isotope-labeled bio-molecules. The approach is based on band-selective Hadamard-type frequency encoding, which disperses the correlation peaks into different sub-spectra. This makes it possible to apply low-dimensionality-based NMR techniques to larger molecular systems. Here we demonstrate the use of band-selective Hadamard frequency labeling for fast protein resonance assignment, based on our recently proposed suite of 2D experiments (Bersch et al., 2003).     

Complex and simple sequences in human repeated DNAs

Highly repeated human DNA sequences were isolated by isopycnic centrifugation, or were eluted from gels after restriction enzyme cleavage. High molecular weight DNA peaks separable from the bulk of the DNA in a variety of gradients were shown to consist of very simple sequences characteristic of simple satellite DNAs; DNA fingerprint studies indicated each of these peaks could consist of tandem repeats of a specific oligonucleotide sequence as low as 10 base pairs (bp) long. All the gradient peaks could be assigned to one of two sequence groups and several different buoyant density peaks revealed the same sequence. — Restriction fragment multimers did not share common sequences with the satellite DNAs as judged by hybridization data. They could not be separated by isopycnic centrifugation. Furthermore these highly repeated DNAs were more complex in sequence and more variable than the satellites. Even the smallest (50 bp) fragments by depurination and other direct sequencing methods were shown to be more complex than the high molecular weight satellite peaks. — The idea that subsets of repeated DNAs may be defined by sequence complexity, possibly with discrete or separable functions, is proposed.     

Quantitative analysis and molecular species fingerprinting of triacylglyceride molecular species directly from lipid extracts of biological samples by electrospray ionization tandem mass spectrometry

Herein we describe a rapid, simple, and reliable method for the quantitative analysis and molecular species fingerprinting of triacylglycerides (TAG) directly from chloroform extracts of biological samples. Previous attempts at direct TAG quantitation by positive-ion electrospray ionization mass spectrometry (ESI/MS) were confounded by the presence of overlapping peaks from choline glycerophospholipids requiring chromatographic separation of lipid extracts prior to ESI/MS analyses. By exploiting the rapid loss of phosphocholine from choline glycerophospholipids, in conjunction with neutral-loss scanning for individual fatty acids, overlapping peaks in the ESI mass spectrum were deconvoluted generating a detailed molecular species fingerprint of individual TAG molecular species directly from chloroform extracts of biological samples. This method readily detects as little as 0.1 pmol of each TAG molecular species from chloroform extracts and is linear over a 1000-fold dynamic range. The sensitivity of individual TAG molecular species to ESI/MS/MS analyses correlated with the unsaturation index and inversely correlated with total aliphatic chain length of TAG. An algorithm was developed which identifies sensitivity factors, thereby allowing the rapid quantitation and molecular species fingerprinting of TAG molecular species directly from chloroform extracts of biological samples.     

Cationic double-hydrophilic model networks: synthesis,characterization, modeling and protein adsorption studies

Group transfer polymerization (GTP) was used for the preparation of eight networks based on two hydrophilic monomers, 2-(dimethylamino)ethyl methacrylate (DMAEMA) and poly(ethylene glycol) methacrylate (PEGMA). Ethylene glycol dimethacrylate (EGDMA) served as the cross-linker, whereas 1,4-bis(methoxytrimethylsiloxymethylene)cyclohexane (MTSMC) was used as a bifunctional initiator. Seven of the networks had linear segments of accurate molecular weight between the cross-links, i.e., they were model networks, whereas the eighth was an equimolar randomly cross-linked network. Five of the seven model networks were based on ABA triblock copolymers with PEGMA midblocks and DMAEMA endblocks, in which the DMAEMA/PEGMA ratio was varied. The remaining two model networks were equimolar isomers, the one based on BAB triblocks (with a DMAEMA midblock) and the other based on the statistical copolymer. The degrees of swelling of all of the networks were measured as a function of pH and were found to increase below pH 7. The degrees of swelling at low pH values increased with the percentage of the DMAEMA monomer, which is ionized under these conditions. These swelling results were confirmed qualitatively by theoretical calculations. Finally, the pH-dependence of the adsorption of the proteins pepsin, bovine serum albumin, and lysozyme onto one of the model networks was studied.     

Automated multi-dimensional liquid chromatography: sample preparation and identification of peptides from human blood filtrate

A comprehensive on-line sample clean-up with an integrated two-dimensional HPLC system was developed for the analysis of natural peptides. Samples comprised of endogenous peptides with molecular weights up to 20 kDa were generated from human hemofiltrate (HF) obtained from patients with chronic renal failure. The (poly-)peptides were separated using novel silica-based restricted access materials with strong cation-exchange functionalities (SCX-RAM). The size-selective sample fractionation step is followed by cation-exchange chromatography as the first dimension. The subsequent second dimension of separation is based on hydrophobic interaction using four parallel short reversed-phase (RP) columns implemented via a fully automated column switching technique. More than 1000 peaks were resolved within the total analysis time of 96 min. Substances of selected peaks were sampled to analyse their molecular weights by off-line MALDI-TOF mass spectrometry and to determine their amino acid sequence by Edman degradation. The potential for comprehensive peptide mapping and identification is demonstrated.     

Cationic amphiphilic model networks based on symmetrical ABCBA pentablock terpolymers: synthesis, characterization, and modeling

Eight isomeric networks based on equimolar terpolymers were synthesized using group transfer polymerization (GTP) and were characterized in terms of their swelling properties. Two hydrophilic monomers, the nonionic methoxy hexa(ethylene glycol) methacrylate (HEGMA) and the ionizable 2-(dimethylamino)ethyl methacrylate (DMAEMA), and a hydrophobic (nonionic) monomer, methyl methacrylate (MMA), were employed for the syntheses. 1,4-Bis(methoxytrimethylsiloxymethylene)cyclohexane (MTSMC) was used as the bifunctional GTP initiator, while ethylene glycol dimethacrylate (EGDMA) served as the cross-linker. Seven of the networks were model networks, six of which were based on the symmetrical pentablock terpolymers ABCBA, ACBCA, BACAB, BCACB, CBABC, and CABAC, whereas the seventh model network was based on the statistical terpolymer. The eighth network was a randomly cross-linked network based on the statistical terpolymer, prepared by the simultaneous quaterpolymerization of the three monomers and the cross-linker. The molecular weights and molecular weight distributions of the linear pentablock terpolymer precursors, as well as those of their homopolymer and ABA triblock copolymer precursors, were characterized by gel permeation chromatography (GPC) in tetrahydrofuran. The sol fraction of each network was measured and found to be relatively low. The aqueous degrees of swelling of all networks were found to increase at acidic pH due to the ionization of the DMAEMA tertiary amine units. The acidic degrees of swelling of the pentablock terpolymer networks were lower than those of their statistical counterparts due to microphase separation in the former type of networks, also confirmed by thermodynamic calculations and small-angle neutron scattering experiments.     

Determination of the amino acid sequence in a cyclic lipopeptide using MS with DHT mechanism

A TOF MS/MS method to directly determine the amino acid sequence in a cyclic lipopeptide without its hydrolysis is described. The fragments of the peptide and the hydrocarbon chains were identified through comparing the MS of two analogues of the lipopeptide; the connecting relationship of amino acid residues in the lipopeptide was determined based on the difference of mass to charge ratio between peaks in the MS spectra and the amino acid analysis; and finally, according to the mechanism of double hydrogen transfer(DHT) the C-terminal of peptide and hydroxy aliphatic acid in the lipopeptide was directly determined without the hydrolysis. The determined sequence of amino acid residues in the cyclic lipopeptide is also supported by the rest peaks in the MS spectra grounded on simple fragmenting mechanism. This method can be used to determine the amino acid sequence in any aliphatic acid loop-inlaying cyclic lipopeptides.     

B-cell ligand processing pathways detected by large-scale comparative analysis

The initiation of B-cell ligand recognition is a critical step for the generation of an immune response against foreign bodies.We sought to identify the biochemical pathways involved in the B-cell ligand recognition cascade and sets of ligands that trigger similar immunological responses.We utilized several comparative approaches to analyze the gene coexpression networks generated from a set of microarray experiments spanning 33 different ligands.First,we compared the degree distributions of the generated networks.Second,we utilized a pairwise network alignment algorithm,BiNA,to align the networks based on the hubs in the networks.Third,we aligned the networks based on a set of K_EGG pathways.We summarized our results by constructing a consensus hierarchy of pathways that are involved in B cell ligand recognition.The resulting pathways were further validated through literature for their common physiological responses.Collectively,the results based on our comparative analyses of degree distributions,alignment of hubs,and alignment based on KEGG pathways provide a basis for molecular characterization of the immune response states of B-cells and demonstrate the power of comparative approaches(e.g.,gene coexpression network alignment algorithms) in elucidating biochemical pathways involved in complex signaling events in cells.     

Identification and differentiation of Cryptosporidium species by capillary electrophoresis single-strand conformation polymorphism

Cryptosporidium species generally lack distinguishing morphological traits, and consequently, molecular methods are commonly used for parasite identification. Various methods for Cryptosporidium identification have been proposed, each with their advantages and disadvantages. In this study, we show that capillary electrophoresis coupled with single-strand conformation polymorphism (CE-SSCP) is a rapid, simple and cost-effective method for the identification of Cryptosporidium species and genotypes. Species could be readily differentiated based on the SSCP mobility of amplified 18S rRNA gene molecules. Clones that differed by single-nucleotide polymorphisms could be distinguished on CE-SSCP mobility. Profiles of species known to have heterogenic copies of 18S rRNA gene contained multiple peaks. Cloning and sequencing of Cryptosporidium parvum, Cryptosporidium hominis, Cryptosporidium fayeri and Cryptosporidium possum genotype 18S rRNA gene amplicons confirmed that these multiple peaks represented type A and type B 18S rRNA gene copies. CE-SSCP provides a reliable and sensitive analysis for epidemiological studies, environmental detection and diversity screening.     

Imaging of functional connectivity in the mouse brain

Functional neuroimaging (e.g., with fMRI) has been difficult to perform in mice, making it challenging to translate between human fMRI studies and molecular and genetic mechanisms. A method to easily perform large-scale functional neuroimaging in mice would enable the discovery of functional correlates of genetic manipulations and bridge with mouse models of disease. To satisfy this need, we combined resting-state functional connectivity mapping with optical intrinsic signal imaging (fcOIS). We demonstrate functional connectivity in mice through highly detailed fcOIS mapping of resting-state networks across most of the cerebral cortex. Synthesis of multiple network connectivity patterns through iterative parcellation and clustering provides a comprehensive map of the functional neuroarchitecture and demonstrates identification of the major functional regions of the mouse cerebral cortex. The method relies on simple and relatively inexpensive camera-based equipment, does not require exogenous contrast agents and involves only reflection of the scalp (the skull remains intact) making it minimally invasive. In principle, fcOIS allows new paradigms linking human neuroscience with the power of molecular/genetic manipulations in mouse models.     

Multiple molecular forms of murine thymocyte-stimulating factor

Murine thymocyte stimulating factor (TSF) was found to sediment in sucrose density gradients on a broad band with peaks at about 2.60 S and 2.0 S. Two main peaks of TSF activity (with buoyant densities of 1.34 and 1.28 g/ml) were found in CsCl density gradients. Gel chromatography on Sephadex G-100 columns of the material sedimented in sucrose or CsCl density gradients originated multiple peaks of TSF activity with various molecular weights. Heterogeneity of molecular forms of TSF was also found upon dilution of the factor. The lowest molecular weights found were 4000 and 4700 daltons. When Sephadex fractions containing the low molecular weight material were pooled and rerun on Sephadex columns, molecular species with a wide range of molecular weights were found. Temperature also affects the appearance of the low molecular weight forms of TSF. Most of the experiments presented in this work were carried out with Sephadex-purified TSF. Multiple molecular forms, however, and, in particular, the forms with molecular weights of 4000 and 4700 daltons were found also with TSF-Fraction IIIa, a highly purified preparation of this factor.