Identification of oscillatory genes in somitogenesis from functional genomic analysis of a human mesenchymal stem cell model

During somitogenesis, oscillatory expression of genes in the notch and wnt signaling pathways plays a key role in regulating segmentation. These oscillations in expression levels are elements of a species-specific developmental mechanism. To date, the periodicity and components of the human clock remain unstudied. Here we show that a human mesenchymal stem/stromal cell (MSC) model can be induced to display oscillatory gene expression. We observed that the known cycling gene HES1 oscillated with a 5 h period consistent with available data on the rate of somitogenesis in humans. We also observed cycling of Hes1 expression in mouse C2C12 myoblasts with a period of 2 h, consistent with previous in vitro and embryonic studies. Furthermore, we used microarray and quantitative PCR (Q-PCR) analysis to identify additional genes that display oscillatory expression both in vitro and in mouse embryos. We confirmed oscillatory expression of the notch pathway gene Maml3 and the wnt pathway gene Nkd2 by whole mount in situ hybridization analysis and Q-PCR. Expression patterns of these genes were disrupted in Wnt3a(tm1Amc) mutants but not in Dll3(pu) mutants. Our results demonstrate that human and mouse in vitro models can recapitulate oscillatory expression observed in embryo and that a number of genes in multiple developmental pathways display dynamic expression in vitro.     

The oscillatory network of glycolysis: a model for multiple phase output derived from a single master oscillator

In oscillatory glycolysing yeast cells the reaction rates of 11 enzymic steps are almost completely in phase whereas the concentrations of intermediates fluctuate with highly different phases. Metabolite concentrations depend on the amplitudes of reaction rates even more than on their phases. In this way almost, synchronized enzymes can given rise to a variety of phases, resembling a "master oscillatory" which controls a series of peripheric event with different phases.     

Cellular-based modeling of oscillatory dynamics in brain networks

Oscillatory, population activities have long been known to occur in our brains during different behavioral states. We know that many different cell types exist and that they contribute in distinct ways to the generation of these activities. I review recent papers that involve cellular-based models of brain networks, most of which include theta, gamma and sharp wave-ripple activities. To help organize the modeling work, I present it from a perspective of three different types of cellular-based modeling: 'Generic', 'Biophysical' and 'Linking'. Cellular-based modeling is taken to encompass the four features of experiment, model development, theory/analyses, and model usage/computation. The three modeling types are shown to include these features and interactions in different ways.     

Co-operative components, spatial localization and oscillatory cellular dynamics

The consequences of spatial localization on the dynamic behavior of catalytic components whose outputs control one another's activities but are not the substrates of those activities are explored. For some simple model systems we naturally find limit cycle oscillations, multiple asymptotic states, dependence of the dynamic states upon spatial configuration of the components, and complex switching between dynamic states as a function of external perturbation to the system.     

Spontaneous calcium oscillatory patterns in mammotropes display non-random dynamics

We previously showed that primary rat mammotropes exhibited four distinct patterns of 'spontaneous' free intracellular calcium ([Ca2+]i) oscillatory behavior: a quiescent state A and three oscillatory states B,C&D, which differed in frequency/amplitude characteristics. When [Ca2+]i was monitored in 10 min windows separated by several hours, these phenotypes were frequently found to interconvert, raising the question about whether these transitions were random or ordered events. We reasoned that if such activity were random, then neither episode duration nor transitional probabilities should differ among phenotypes. We tested this logic in the current study by making long-term, continuous measurements of [Ca2+]i in mammotropes microinjected with Fura-2-dextran and identified by their ability to express a prolactin promoter-driven reporter plasmid. We found that transitions occurred in ~25% of cells (n = 36 from 9 independent experiments) once every 1-5 h and demarcated phenotype episodes of different duration (A, 1.04 +/- 0.2 h; B, 1.64 +/- 0.3 h; C, 2.45 +/- 0.62 h; D, 0.90 +/- 0.2 h, mean +/- SEM). Moreover, some transitions occurred more frequently than others and linked specific phenotypes into a common pattern: C to B to A. Our results demonstrate that the seemingly spontaneous nature of [Ca2+]i phenotype transitions are, in fact, ordered and support the view that they comprise a structured 'code' like that proposed to underlie calcium-dependent regulation of exocytosis and gene expression.     

Locally oscillatory motion of RNA helix derived from linear relationships of backbone torsion angles

A linear relationship in each of the torsion angle pairs, α-β, β-?, ?-ζ, and α-γ, has been found by applying a statistical method based on the concept of circular variates to backbone torsion angle data of helical in yeast tTNA^Phe. A series of helical dimer models generated with these relationships have been found to be stereochemically acceptable, and the models also indicate that the backbone unit in the RNA helix is geometrically capable of an oscillatory motion with the distance of about 3.4 Å between adjacent bases. The motion of the backbone unit is analogous to that of a helical spring. The adjacent bases, because of being attached to the backbone, oscillate in a manner similar to the oscillatory dimer model proposed by Davis and Tinoco [Davis, R. C. & Tinoco, I., Jr. (1968) Biopolymers 6 , 223–242]. Here, the oscillation of the backbone unit in the RNA helix is discussed in terms of two geometrical quantities: the torsion (τ) and curvature (κ) of the helix. On these lines, a stereochemical model of RNA strand separation is proposed.     

Telomere length dynamics and chromosomal instability in cells derived from telomerase null mice

To study the effect of continued telomere shortening on chromosome stability, we have analyzed the telomere length of two individual chromosomes (chromosomes 2 and 11) in fibroblasts derived from wild-type mice and from mice lacking the mouse telomerase RNA (mTER) gene using quantitative fluorescence in situ hybridization. Telomere length at both chromosomes decreased with increasing generations of mTER-/- mice. At the 6th mouse generation, this telomere shortening resulted in significantly shorter chromosome 2 telomeres than the average telomere length of all chromosomes. Interestingly, the most frequent fusions found in mTER-/- cells were homologous fusions involving chromosome 2. Immortal cultures derived from the primary mTER-/- cells showed a dramatic accumulation of fusions and translocations, revealing that continued growth in the absence of telomerase is a potent inducer of chromosomal instability. Chromosomes 2 and 11 were frequently involved in these abnormalities suggesting that, in the absence of telomerase, chromosomal instability is determined in part by chromosome-specific telomere length. At various points during the growth of the immortal mTER-/- cells, telomere length was stabilized in a chromosome-specific man-ner. This telomere-maintenance in the absence of telomerase could provide the basis for the ability of mTER-/- cells to grow indefinitely and form tumors.     

Stochastic aspects of oscillatory Ca2+ dynamics in hepatocytes

Signal-induced Ca²⁺ oscillations have been observed in many cell types and play a primary role in cell physiology. Although it is the regular character of these oscillations that first catches the attention, a closer look at time series of Ca²⁺ increases reveals that the fluctuations on the period during individual spike trains are far from negligible. Here, we perform a statistical analysis of the regularity of Ca²⁺ oscillations in norepinephrine-stimulated hepatocytes and find that the coefficient of variation lies between 10% and 15%. Stochastic simulations based on Gillespie's algorithm and considering realistic numbers of Ca²⁺ ions and inositol trisphosphate (InsP₃) receptors account for this variability if the receptors are assumed to be grouped in clusters of a few tens of channels. Given the relatively small number of clusters (∼200), the model predicts the existence of repetitive spikes induced by fluctuations (stochastic resonance). Oscillations of this type are found in hepatocytes at subthreshold concentrations of norepinephrine. We next predict with the model that the isoforms of the InsP₃ receptor can affect the variability of the oscillations. In contrast, possible accompanying InsP₃ oscillations have no impact on the robustness of signal-induced repetitive Ca²⁺ spikes.     

Structure and dynamics of interfacial water in an Lalpha phase lipid bilayer from molecular dynamics simulations

Based on molecular dynamics simulations, an analysis of structure and dynamics is performed on interfacial water at a liquid crystalline dipalmitoylphosphatidycholine/water system. Water properties relevant for understanding NMR relaxation are emphasized. The first and second rank orientational order parameters of the water O-H bonds were calculated, where the second rank order parameter is in agreement with experimental determined quadrupolar splittings. Also, two different interfacial water regions (bound water regions) are revealed with respect to different signs of the second rank order parameter. The water reorientation correlation function reveals a mixture of fast and slow decaying parts. The fast (ps) part of the correlation function is due to local anisotropic water reorientation whereas the much slower part is due to more complicated processes including lateral diffusion along the interface and chemical exchange between free and bound water molecules. The 100-ns-long molecular dynamics simulation at constant pressure (1 atm) and at a temperature of 50 degrees C of 64 lipid molecules and 64 x 23 water molecules lack a slow water reorientation correlation component in the ns time scale. The (2)H(2)O powder spectrum of the dipalmitoylphosphatidycholine/water system is narrow and consequently, the NMR relaxation time T(2) is too short compared to experimental results.     

Evolution of cannibalistic traits: scenarios derived from adaptive dynamics

The evolution of cannibalistic traits in consumer populations is studied in this paper with the approach of adaptive dynamics theory. The model is kept at its minimum complexity by eliminating some environmental characteristics, like heterogeneity and seasonalities, and by hiding the size-structure of the population. Evolutionary dynamics are identified through numerical bifurcation analysis, applied both to the ecological (resident-mutant) model and to the canonical equation of adaptive dynamics. The result is a rich catalog of evolutionary scenarios involving evolutionary stable strategies and branching points both in the monomorphic and dimorphic dynamics. The possibility of evolutionary extinction of highly cannibalistic populations is also ascertained. This allows one to explain why cannibalism can be a transient stage of evolution.     

Neural rate equations for bursting dynamics derived from conductance-based equations

A method of obtaining rate equations from conductance-based equations is developed and applied to fast-spiking and bursting neocortical neurons. It involves splitting systems of conductance-based equations into fast and slow subsystems, and averaging the effects of fast terms that drive the slowly varying quantities by showing that their average is closely proportional to the firing rate. The dependence of the firing rate on the injected current is then approximated in the analysis. The resulting behavior of the slow variables is then substituted back into the fast equations, with the further approximation of replacing the fast voltages in these terms by effective values. For bursting neurons the method yields two coupled limit-cycle oscillators: a self-exciting oscillator for the slow variables that commences limit-cycle oscillations at a critical current and modulates a fast spike-generating oscillator, thereby leading to slowly modulated bursts with a group of spikes in each burst. The dynamics of these coupled oscillators are then verified against those of the conductance-based equations. Finally, it is shown how to place the results in a form suitable for use in mean-field equations for neural population dynamics.     

Insights into correlated motions and long-range interactions in CheY derived from molecular dynamics simulations

CheY is a response regulator protein involved in bacterial chemotaxis. Much is known about its active and inactive conformations, but little is known about the mechanisms underlying long-range interactions or correlated motions. To investigate these events, molecular dynamics simulations were performed on the unphosphorylated, inactive structure from Salmonella typhimurium and the CheY-BeF(3)(-) active mimic structure (with BeF(3)(-) removed) from Escherichia coli. Simulations utilized both sequences in each conformation to discriminate sequence- and structure-specific behavior. The previously identified conformational differences between the inactive and active conformations of the strand-4-helix-4 loop, which are present in these simulations, arise from the structural, and not the sequence, differences. The simulations identify previously unreported structure-specific flexibility features in this loop and sequence-specific flexibility features in other regions of the protein. Both structure- and sequence-specific long-range interactions are observed in the active and inactive ensembles. In the inactive ensemble, two distinct mechanisms based on Thr-87 or Ile-95 rotameric forms, are observed for the previously identified g+ and g- rotamer sampling by Tyr-106. These molecular dynamics simulations have thus identified both sequence- and structure-specific differences in flexibility, long-range interactions, and rotameric form of key residues. Potential biological consequences of differential flexibility and long-range correlated motion are discussed.     

Conformational dynamics of full-length inducible human Hsp70 derived from microsecond molecular dynamics simulations in explicit solvent

Human 70?kDa heat shock protein (hHsp70) is an ATP-dependent chaperone and is currently an important target for developing new drugs in cancer therapy. Knowledge of the conformations of hHsp70 is central to understand the interactions between its nucleotide-binding domain (NBD) and substrate-binding domain (SBD) and is a prerequisite to design inhibitors. The conformations of ADP-bound (or nucleotide-free) hHsp70 and ATP-bound hHsp70 was investigated by using unbiased all-atom molecular dynamics (MD) simulations of homology models of hHsp70 in explicit solvent on a timescale of .5 and 2.7 μs, respectively. The conformational heterogeneity of hHsp70 was analyzed by computing effective free-energy landscapes (FELs) and distance distribution between selected pair of residues. These theoretical data were compared with those extracted from single-molecule Förster resonance energy transfer (FRET) experiments and to small-angle X-ray scattering (SAXS) data of Hsp70 homologs. The distance between a pair of residues in FRET is systematically larger than the distance computed in MD which is interpreted as an effect of the size and of the dynamics of the fluorescent probes. The origin of the conformational heterogeneity of hHsp70 in the ATP-bound state is due to different binding modes of the helix B of the SBD onto the NBD. In the ADP-bound (or nucleotide-free) state, it arises from the different closed conformations of the SBD and from the different positions of the SBD relative to the NBD. In each nucleotide-binding state, Hsp70 is better represented by an ensemble of conformations on a μs timescale corresponding to different local minima of the FEL.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:30     

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.     

Uptake of alpha-aminoisobutyric acid and phosphate by membrane vesicles derived from growing and quiescent fibroblasts.

Membrane vesicles derived principally from the plasma membrane and endoplasmic reticulum of mouse 3T3 cells transformed by Simian virus 40 take up alpha-aminoisobutyric acid (AIB) and phosphate (Pi). When NaCl is added simultaneously with AIB or Pi, uptake rises two- to three-times above the equilibrium to accumulate AIB or Pi over the control value, in the presence of a Na+ gradient, is almost lost in membrane vesicles derived from benzpyrene-transformed 3T3 cells (BP3T3) arrested in the G1 phase of the cell cycle by serum starvation. When added to the membranes with NaCl and the uptake substrate, a combination of fibroblast growth factor (FGF) and epidermal growth factor EGF restores the ability of the membranes to accumulate AIB and Pi over the control value.     

The onset of oscillatory dynamics in models of multiple disease strains

 We examine a generalised SIR model for the infection dynamics of four competing disease strains. This model contains four previously-studied models as special cases. The different strains interact indirectly by the mechanism of cross-immunity; individuals in the host population may become immune to infection by a particular strain even if they have only been infected with different but closely related strains. Several different models of cross-immunity are compared in the limit where the death rate is much smaller than the rate of recovery from infection. In this limit an asymptotic analysis of the dynamics of the models is possible, and we are able to compute the location and nature of the Takens–Bogdanov bifurcation associated with the presence of oscillatory dynamics observed by previous authors. Received: 5 December 2001 / Revised version: 5 May 2002 / Published online: 17 October 2002 Keywords or phrases: Infection – Pathogen – Epidemiology – Multiple strains – Cross-immunity – Oscillations – Dynamics – Bifurcations     

Growing media constituents determine the microbial nitrogen conversions in organic growing media for horticulture

Vegetables and fruits are an important part of a healthy food diet, however, the eco‐sustainability of the production of these can still be significantly improved. European farmers and consumers spend an estimated €15.5 billion per year on inorganic fertilizers and the production of N‐fertilizers results in a high carbon footprint. We investigated if fertilizer type and medium constituents determine microbial nitrogen conversions in organic growing media and can be used as a next step towards a more sustainable horticulture. We demonstrated that growing media constituents showed differences in urea hydrolysis, ammonia and nitrite oxidation and in carbon dioxide respiration rate. Interestingly, mixing of the growing media constituents resulted in a stimulation of the function of the microorganisms. The use of organic fertilizer resulted in an increase in amoA gene copy number by factor 100 compared to inorganic fertilizers. Our results support our hypothesis that the activity of the functional microbial community with respect to nitrogen turnover in an organic growing medium can be improved by selecting and mixing the appropriate growing media components with each other. These findings contribute to the understanding of the functional microbial community in growing media and its potential role towards a more responsible horticulture.