A hierarchical whole-body modeling approach elucidates the link between in Vitro insulin signaling and in Vivo glucose homeostasis

Type 2 diabetes is a metabolic disease that profoundly affects energy homeostasis. The disease involves failure at several levels and subsystems and is characterized by insulin resistance in target cells and tissues (i.e. by impaired intracellular insulin signaling). We have previously used an iterative experimental-theoretical approach to unravel the early insulin signaling events in primary human adipocytes. That study, like most insulin signaling studies, is based on in vitro experimental examination of cells, and the in vivo relevance of such studies for human beings has not been systematically examined. Herein, we develop a hierarchical model of the adipose tissue, which links intracellular insulin control of glucose transport in human primary adipocytes with whole-body glucose homeostasis. An iterative approach between experiments and minimal modeling allowed us to conclude that it is not possible to scale up the experimentally determined glucose uptake by the isolated adipocytes to match the glucose uptake profile of the adipose tissue in vivo. However, a model that additionally includes insulin effects on blood flow in the adipose tissue and GLUT4 translocation due to cell handling can explain all data, but neither of these additions is sufficient independently. We also extend the minimal model to include hierarchical dynamic links to more detailed models (both to our own models and to those by others), which act as submodules that can be turned on or off. The resulting multilevel hierarchical model can merge detailed results on different subsystems into a coherent understanding of whole-body glucose homeostasis. This hierarchical modeling can potentially create bridges between other experimental model systems and the in vivo human situation and offers a framework for systematic evaluation of the physiological relevance of in vitro obtained molecular/cellular experimental data.     

The structure of the carboxyl terminus of striated alpha-tropomyosin in solution reveals an unusual parallel arrangement of interacting alpha-helices

Coiled coils are well-known as oligomerization domains, but they are also important sites of protein-protein interactions. We determined the NMR solution structure and backbone (15)N relaxation rates of a disulfide cross-linked, two-chain, 37-residue polypeptide containing the 34 C-terminal residues of striated muscle alpha-tropomyosin, TM9a(251-284). The peptide binds to the N-terminal region of TM and to the tropomyosin-binding domain of the regulatory protein, troponin T. Comparison of the NMR solution structure of TM9a(251-284) with the X-ray structure of a related peptide [Li, Y., Mui, S., Brown, J. H., Strand, J., Reshetnikova, L., Tobacman, L. S., and Cohen, C. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 7378-7383] reveals significant differences. In solution, residues 253-269 (like most of the tropomyosin molecule) form a canonical coiled coil. Residues 270-279, however, are parallel, linear helices, novel for tropomyosin. The packing between the parallel helices results from unusual interface residues that are atypical for coiled coils. Y267 has poor packing at the coiled-coil interface and a lower R(2) relaxation rate than neighboring residues, suggesting there is conformational flexibility around this residue. The last five residues are nonhelical and flexible. The exposed surface presented by the parallel helices, and the flexibility around Y267 and the ends, may facilitate binding to troponin T and formation of complexes with the N-terminus of tropomyosin and actin. We propose that unusual packing and flexibility are general features of coiled-coil domains in proteins that are involved in intermolecular interactions.     

Estimating genetic drift and effective population size from temporal shifts in dominant gene marker frequencies

Measurement of temporal change in allele frequencies represents an indirect method for estimating the genetically effective size of populations. When allele frequencies are estimated for gene markers that display dominant gene expression, such as, e.g. random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) markers, the estimates can be seriously biased. We quantify bias for previous allele frequency estimators and present a new expression that is generally less biased and provides a more precise assessment of temporal allele frequency change. We further develop an estimator for effective population size that is appropriate when dealing with dominant gene markers. Comparison with estimates based on codominantly expressed genes, such as allozymes or microsatellites, indicates that about twice as many loci or sampled individuals are required when using dominant markers to achieve the same precision.     

Phylogenetic relationships of Betula species (Betulaceae) based on nuclear ADH and chloroplast matK sequences

The phylogenetic relationships within the genus Betula (Betulaceae) were investigated using a part of the nuclear ADH gene and DNA sequences of the chloroplast matK gene with parts of its flanking regions. Two well-supported phylogenetic groups could be identified in the chloroplast DNA sequence: one containing the three American species B. lenta, B. alleghaniensis, and B. papyrifera and the other including all the other species studied. The ADH gene displayed more variation, and three main groups could be identified. In disagreement with the classical division of the genus Betula, B. schmidtii and B. nana grouped with the species in subgenus Betula, and B. ermanii grouped with species in subgenus Chamaebetula, including B. humilis and B. fruticosa. The ADH phylogeny suggests that several independent polyploidizations within the genus Betula could have taken place. The ADH and chloroplast phylogenies were in part incongruent due to the placement of B. papyrifera. The most likely reason for this seems to be cytoplasmic introgression.     

Genome-wide homozygosity mapping localizes a gene for autosomal recessive non-progressive infantile ataxia to 20q11-q13

Autosomal recessive ataxias represent genetic and clinical heterogeneity. Unsteady gait is often accompanied by poor coordination of limbs, speech, and eye movements. To date, seven genes have been identified. In addition, five chromosomal loci have been localized in non-related families. Here, we report homozygosity mapping of a novel locus to a 19.5-cM region on chromosome 20q11–q13 in a large inbred Norwegian family with infantile non-progressive ataxia.L. Tranebjaerg, T.M. Teslovich, and M. Jones contributed equally to this work     

Scene segmentation by spike synchronization in reciprocally connected visual areas. II. Global assemblies and synchronization on larger space and time scales

 We present further simulation results of the model of two reciprocally connected visual areas proposed in the first paper [Knoblauch and Palm (2002) Biol Cybern 87:151–167]. One area corresponds to the orientation–selective subsystem of the primary visual cortex, the other is modeled as an associative memory representing stimulus objects according to Hebbian learning. We examine the scene-segmentation capability of our model on larger time and space scales, and relate it to experimental findings. Scene segmentation is achieved by attention switching on a time-scale longer than the gamma range. We find that the time-scale can vary depending on habituation parameters in the range of tens to hundreds of milliseconds. The switching process can be related to findings concerning attention and biased competition, and we reproduce experimental poststimulus time histograms (PSTHs) of single neurons under different stimulus and attentional conditions. In a larger variant the model exhibits traveling waves of activity on both slow and fast time-scales, with properties similar to those found in experiments. An apparent weakness of our standard model is the tendency to produce anti-phase correlations for fast activity from the two areas. Increasing the inter-areal delays in our model produces alternations of in-phase and anti-phase oscillations. The experimentally observed in-phase correlations can most naturally be obtained by the involvement of both fast and slow inter-areal connections; e.g., by two axon populations corresponding to fast-conducting myelinated and slow-conducting unmyelinated axons. Received: 22 August 2001 / Accepted in revised form: 8 April 2002     

Engineering crop plants: getting a handle on phosphate

In plant seeds, most of the phosphate is in the form of phytic acid. Phytic acid is largely indigestible by monogastric animals and is the single most important factor hindering the uptake of a range of minerals. Engineering crop plants to produce a heterologous phytase improves phosphate bioavailability and reduces phytic acid excretion. This reduces the phosphate load on agricultural ecosystems and thereby alleviates eutrophication of the aquatic environment. Improved phosphate availability also reduces the need to add inorganic phosphate, a non-renewable resource. Iron and zinc uptake might be improved, which is significant for human nutrition in developing countries.     

Wiener-like system identification in physiology

Summary Applications of Wiener-like identification methods to biological systems have revealed several limitations of this technique. These practical limitations correspond to conceptual and mathematical problems intrinsic to this kind of identification of nonlinear systems.     

The Archaea Monophyly Issue: A Phylogeny of Translational Elongation Factor G(2) Sequences Inferred from an Optimized Selection of Alignment Positions

A global alignment of EF-G(2) sequences was corrected by reference to protein structure. The selection of characters eligible for construction of phylogenetic trees was optimized by searching for regions arising from the artifactual matching of sequence segments unique to different phylogenetic domains. The spurious matchings were identified by comparing all sections of the global alignment with a comprehensive inventory of significant binary alignments obtained by BLAST probing of the DNA and protein databases with representative EF-G(2) sequences. In three discrete alignment blocks (one in domain II and two in domain IV), the alignment of the bacterial sequences with those of Archaea–Eucarya was not retrieved by database probing with EF-G(2) sequences, and no EF-G homologue of the EF-2 sequence segments was detected by using partial EF-G(2) sequences as probes in BLAST/FASTA searches. The two domain IV regions (one of which comprises the ADP-ribosylatable site of EF-2) are almost certainly due to the artifactual alignment of insertion segments that are unique to Bacteria and to Archaea–Eucarya. Phylogenetic trees have been constructed from the global alignment after deselecting positions encompassing the unretrieved, spuriously aligned regions, as well as positions arising from misalignment of the G′ and G″ subdomain insertion segments flanking the ``fifth' consensus motif of the G domain (?varsson, 1995). The results show inconsistencies between trees inferred by alternative methods and alternative (DNA and protein) data sets with regard to Archaea being a monophyletic or paraphyletic grouping. Both maximum-likelihood and maximum-parsimony methods do not allow discrimination (by log-likelihood difference and difference in number of inferred substitutions) between the conflicting (monophyletic vs. paraphyletic Archaea) topologies. No specific EF-2 insertions (or terminal accretions) supporting a crenarchaeal–eucaryal clade are detectable in the new EF-G(2) sequence alignment.