Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

Dietary Energy Partition: The Central Role of Glucose

Humans have developed effective survival mechanisms under conditions of nutrient (and energy) scarcity. Nevertheless, today, most humans face a quite different situation: excess of nutrients, especially those high in amino-nitrogen and energy (largely fat). The lack of mechanisms to prevent energy overload and the effective persistence of the mechanisms hoarding key nutrients such as amino acids has resulted in deep disorders of substrate handling. There is too often a massive untreatable accumulation of body fat in the presence of severe metabolic disorders of energy utilization and disposal, which become chronic and go much beyond the most obvious problems: diabetes, circulatory, renal and nervous disorders included loosely within the metabolic syndrome. We lack basic knowledge on diet nutrient dynamics at the tissue-cell metabolism level, and this adds to widely used medical procedures lacking sufficient scientific support, with limited or nil success. In the present longitudinal analysis of the fate of dietary nutrients, we have focused on glucose as an example of a largely unknown entity. Even most studies on hyper-energetic diets or their later consequences tend to ignore the critical role of carbohydrate (and nitrogen disposal) as (probably) the two main factors affecting the substrate partition and metabolism.

Cross-Species Analysis of Glycosaminoglycan Binding Proteins Reveals Some Animal Models Are "More Equal" than Others

Glycosaminoglycan (GAG) mimetics are synthetic or semi-synthetic analogues of heparin or heparan sulfate, which are designed to interact with GAG binding sites on proteins. The preclinical stages of drug development rely on efficacy and toxicity assessment in animals and aim to apply these findings to clinical studies. However, such data may not always reflect the human situation possibly because the GAG binding site on the protein ligand in animals and humans could differ. Possible inter-species differences in the GAG-binding sites on antithrombin III, heparanase, and chemokines of the CCL and CXCL families were examined by sequence alignments, molecular modelling and assessment of surface electrostatic potentials to determine if one species of laboratory animal is likely to result in more clinically relevant data than another. For each protein, current understanding of GAG binding is reviewed from a protein structure and function perspective. This combinatorial analysis shows chemokine dimers and oligomers can present different GAG binding surfaces for the same target protein, whereas a cleft-like GAG binding site will differently influence the types of GAG structures that bind and the species preferable for preclinical work. Such analyses will allow an informed choice of animal(s) for preclinical studies of GAG mimetic drugs.

Conservation and Role of Electrostatics in Thymidylate Synthase

Conservation of function across families of orthologous enzymes is generally accompanied by conservation of their active site electrostatic potentials. To study the electrostatic conservation in the highly conserved essential enzyme, thymidylate synthase (TS), we conducted a systematic species-based comparison of the electrostatic potential in the vicinity of its active site. Whereas the electrostatics of the active site of TS are generally well conserved, the TSs from minimal organisms do not conform to the overall trend. Since the genomes of minimal organisms have a high thymidine content compared to other organisms, the observation of non-conserved electrostatics was surprising. Analysis of the symbiotic relationship between minimal organisms and their hosts, and the genetic completeness of the thymidine synthesis pathway suggested that TS from the minimal organism Wigglesworthia glossinidia (W.g.b.) must be active. Four residues in the vicinity of the active site of Escherichia coli TS were mutated individually and simultaneously to mimic the electrostatics of W.g.b TS. The measured activities of the E. coli TS mutants imply that conservation of electrostatics in the region of the active site is important for the activity of TS, and suggest that the W.g.b. TS has the minimal activity necessary to support replication of its reduced genome.

Characterization of differentially expressed genes in liver in response to the rearing temperature of rainbow trout Oncorhynchus mykiss and their heritable differences

To characterize thermal-responsive genes in fish, firstly, juvenile rainbow trout were reared in four different temperature conditions (average temperatures were 10, 14, 18, and 22 °C, respectively) and differentially expressed genes were identified. Gene expression in the liver was analyzed by the differential display method, followed by validation using real-time PCR. Subsequently, to examine whether the identified genes show heritable differences, the gene expression levels were compared among juveniles of three genetically distinct lines of rainbow trout (a strain and two closed colonies) by rearing at two different temperature conditions (average 14 and 22 °C). By rearing at 22 °C, growth retardation was observed compared with fish reared at 14 and 18 °C, and six genes were identified as differentially expressed genes in response to the rearing temperature in the gene expression analyses. With the increase in rearing temperature, gene expressions of a complement C1q and two ribosomal proteins were significantly up-regulated. On the other hand, three metabolic genes (betaine homocysteine methyltransferase, triosephosphate isomerase, and glucose-6-phosphatase) were down-regulated, indicating a metabolic depression due to high temperature. In the subsequent analyses, in response to the rearing temperature (14 and 22 °C), there was a trend that the complement C1q and glucose-6-phosphatase genes showed different expression patterns among the three rainbow trout lines, suggesting heritable differences in these genes. Our study provides information on thermal-responsive genes in fish, and we anticipate it will facilitate further investigation in the thermal biology of fish.

Knockdown of CDKN1C (p57(kip2)) and PHLDA2 results in developmental changes in bovine pre-implantation embryos

Imprinted genes have been implicated in early embryonic, placental, and neonatal development and alterations in expression levels of these genes can lead to growth abnormalities and embryonic lethality. However, little is known about the functions of bovine imprinted genes during the pre-implantation period. Therefore, the objective of this study was to assess the influence of altered expression of imprinted genes on developmental progress of embryos using small interfering RNA (siRNA). Expression levels of 18 imprinted genes (MAGEL2, UBE3A, IGF2R, NAP1L5, TSSC4, PEG3, NDN, CDKN1C, PHLDA2, MKRN3, USP29, NNAT, PEG10, RTL1, IGF2, H19, MIM1, and XIST) were compared between embryos reaching the blastocyst stage and growth-arrested embryos (degenerates) using quantitative real-time PCR (qRT-PCR). Ten genes were found to be differentially expressed between blastocysts and degenerates. The CDKN1C gene showed the highest upregulation in blastocysts whereas PHLDA2 was highly expressed in degenerates. To assess whether the observed differential gene expression was causative or resultant of embryo degeneration, these genes were selected for functional analysis using siRNA. Injection of siRNA specific to PHLDA2 into one-cell zygotes resulted in a substantial increase in blastocyst development, whereas injection of CDKN1C-specific siRNA resulted in a 45% reduction (P = 0.0006) in blastocyst development. RNA-Seq analysis of CDKN1C-siRNA-injected vs. non-injected embryos revealed 51 differentially expressed genes with functions in apoptosis, lipid metabolism, differentiation, and cell cycle regulation. Gene ontology analysis revealed nine pathways related to cell signaling, metabolism, and nucleic acid processing. Overall, our results show that proper expression levels of the imprinted genes CDKN1C and PHLDA2 are critical for embryo development, which suggests that these genes can be used as markers for normal blastocyst formation.

The interaction properties of the human Rab GTPase family--comparative analysis reveals determinants of molecular binding selectivity

Background

Rab GTPases constitute the largest subfamily of the Ras protein superfamily. Rab proteins regulate organelle biogenesis and transport, and display distinct binding preferences for effector and activator proteins, many of which have not been elucidated yet. The underlying molecular recognition motifs, binding partner preferences and selectivities are not well understood.

Methodology/Principal Findings

Comparative analysis of the amino acid sequences and the three-dimensional electrostatic and hydrophobic molecular interaction fields of 62 human Rab proteins revealed a wide range of binding properties with large differences between some Rab proteins. This analysis assists the functional annotation of Rab proteins 12, 14, 26, 37 and 41 and provided an explanation for the shared function of Rab3 and 27. Rab7a and 7b have very different electrostatic potentials, indicating that they may bind to different effector proteins and thus, exert different functions. The subfamily V Rab GTPases which are associated with endosome differ subtly in the interaction properties of their switch regions, and this may explain exchange factor specificity and exchange kinetics.

Conclusions/Significance

We have analysed conservation of sequence and of molecular interaction fields to cluster and annotate the human Rab proteins. The analysis of three dimensional molecular interaction fields provides detailed insight that is not available from a sequence-based approach alone. Based on our results, we predict novel functions for some Rab proteins and provide insights into their divergent functions and the determinants of their binding partner selectivity.

Clustering of HIV-1 Subtypes Based on gp120 V3 Loop electrostatic properties

BACKGROUND

The V3 loop of the glycoprotein gp120 of HIV-1 plays an important role in viral entry into cells by utilizing as coreceptor CCR5 or CXCR4, and is implicated in the phenotypic tropisms of HIV viruses. It has been hypothesized that the interaction between the V3 loop and CCR5 or CXCR4 is mediated by electrostatics. We have performed hierarchical clustering analysis of the spatial distributions of electrostatic potentials and charges of V3 loop structures containing consensus sequences of HIV-1 subtypes.

RESULTS

Although the majority of consensus sequences have a net charge of +3, the spatial distribution of their electrostatic potentials and charges may be a discriminating factor for binding and infectivity. This is demonstrated by the formation of several small subclusters, within major clusters, which indicates common origin but distinct spatial details of electrostatic properties. Some of this information may be present, in a coarse manner, in clustering of sequences, but the spatial details are largely lost. We show the effect of ionic strength on clustering of electrostatic potentials, information that is not present in clustering of charges or sequences. We also make correlations between clustering of electrostatic potentials and net charge, coreceptor selectivity, global prevalence, and geographic distribution. Finally, we interpret coreceptor selectivity based on the N6X7T8|S8X9 sequence glycosylation motif, the specific positive charge location according to the 11/24/25 rule, and the overall charge and electrostatic potential distribution.

CONCLUSIONS

We propose that in addition to the sequence and the net charge of the V3 loop of each subtype, the spatial distributions of electrostatic potentials and charges may also be important factors for receptor recognition and binding and subsequent viral entry into cells. This implies that the overall electrostatic potential is responsible for long-range recognition of the V3 loop with coreceptors CCR5/CXCR4, whereas the charge distribution contributes to the specific short-range interactions responsible for the formation of the bound complex. We also propose a scheme for coreceptor selectivity based on the sequence glycosylation motif, the 11/24/25 rule, and net charge.

The effect of electrostatics on factor H function and related pathologies

Factor H (FH) contributes to the regulation of the complement system by binding to polyanionic surfaces and the proteins C3b/C3c/C3d. This implicates charge and electrostatic interactions in recognition and binding of FH. Despite the large amount of experimental and pathology data the exact mechanism at molecular level is not yet known. We have implemented a computational framework for comparative analysis of the charge and electrostatic diversity of FH modules and C3b domains to identify electrostatic hotspots and predict potential binding sites. Our electrostatic potential clustering analysis shows that charge distributions and electrostatic potential distributions are more useful in understanding C3b-FH interactions than net charges alone. We present a model of non-specific electrostatic interactions of FH with polyanion-rich surfaces and specific interactions with C3b, using our computational data and existing experimental data. We discuss the electrostatic contributions to the formation of the C3b-FH complex and the competition between FH and Factor Bb (Bb) for binding to C3b. We also discuss the significance of mutations of charged amino acids in the pathobiology of FH-mediated disease, such as age-related macular degeneration, atypical hemolytic uremic syndrome, and dense deposit disease. Our data can be used to guide future experimental studies.

Proteomics analysis of flax grown in Chernobyl area suggests limited effect of contaminated environment on seed proteome

The accident at the Chernobyl Nuclear Power Plant (CNPP) on April 26, 1986 is the most serious nuclear disaster in human history. Surprisingly, while the area proximal to the CNPP remains substantially contaminated with long-lived radioisotopes including (90)Sr and (137)Cs, the local ecosystem has been able to adapt. To evaluate plant adaptation, seeds of a local flax (Linum usitatissimum) variety Kyivskyi were sown in radio-contaminated and control fields of the Chernobyl region. A total protein fraction was isolated from mature seeds, and analyzed using 2-dimensional electrophoresis combined with tandem-mass spectrometry. Interestingly, growth of the plants in the radio-contaminated environment had little effect on proteome and only 35 protein spots differed in abundance (p-value of ≤0.05) out of 720 protein spots that were quantified for seeds harvested from both radio-contaminated and control fields. Of the 35 differentially abundant spots, 28 proteins were identified using state-of-the-art MS(E) method. Based on the observed changes, the proteome of seeds from plants grown in radio-contaminated soil display minor adjustments to multiple signaling pathways.