Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Monte Carlo simulations of metasomatic enrichment in the lithosphere and implications for the source of alkaline basalts

One hypothesis for the origin of alkaline lavas erupted on oceanic islands and in intracontinental settings is that they represent the melts of amphibole-rich veins in the lithosphere (or melts of their dehydrated equivalents if metasomatized lithosphere is recycled into the convecting mantle). Amphibole-rich veins are interpreted as cumulates produced by crystallization of low-degree melts of the underlying asthenosphere as they ascend through the lithosphere. We present the results of trace-element modelling of the formation and melting of veins formed in this way with the goal of testing this hypothesis and for predicting how variability in the formation and subsequent melting of such cumulates (and adjacent cryptically and modally metasomatized lithospheric peridotite) would be manifested in magmas generated by such a process. Because the high-pressure phase equilibria of hydrous near-solidus melts of garnet lherzolite are poorly constrained and given the likely high variability of the hypothesized accumulation and remelting processes, we used Monte Carlo techniques to estimate how uncertainties in the model parameters (e.g. the compositions of the asthenospheric sources, their trace-element contents, and their degree of melting; the modal proportions of crystallizing phases, including accessory phases, as the asthenospheric partial melts ascend and crystallize in the lithosphere; the amount of metasomatism of the peridotitic country rock; the degree of melting of the cumulates and the amount of melt derived from the metasomatized country rock) propagate through the process and manifest themselves as variability in the trace-element contents and radiogenic isotopic ratios of model vein compositions and erupted alkaline magma compositions. We then compare the results of the models with amphibole observed in lithospheric veins and with oceanic and continental alkaline magmas. While the trace-element patterns of the near-solidus peridotite melts, the initial anhydrous cumulate assemblage (clinopyroxene +/- garnet +/- olivine +/- orthopyroxene), and the modelled coexisting liquids do not match the patterns observed in alkaline lavas, our calculations show that with further crystallization and the appearance of amphibole (and accessory minerals such as rutile, ilmenite, apatite, etc.) the calculated cumulate assemblages have trace-element patterns that closely match those observed in the veins and lavas. These calculated hydrous cumulate assemblages are highly enriched in incompatible trace elements and share many similarities with the trace-element patterns of alkaline basalts observed in oceanic or continental setting such as positive Nb/La, negative Ce/Pb, and similiar slopes of the rare earth elements. By varying the proportions of trapped liquid and thus simulating the cryptic and modal metasomatism observed in peridotite that surrounds these veins, we can model the variations in Ba/Nb, Ce/Pb, and Nb/U ratios that are observed in alkaline basalts. If the isotopic compositions of the initial low-degree peridotite melts are similar to the range observed in mid-ocean ridge basalt, our model calculations produce cumulates that would have isotopic compositions similar to those observed in most alkaline ocean island basalt (OIB) and continental magmas after similar to 0 center dot 15 Gyr. However, to produce alkaline basalts with HIMU isotopic compositions requires much longer residence times (i.e. 1-2 Gyr), consistent with subduction and recycling of metasomatized lithosphere through the mantle. such as a heterogeneous asthenosphere. These modelling results support the interpretation proposed by various researchers that amphibole-bearing veins represent cumulates formed during the differentiation of a volatile-bearing low-degree peridotite melt and that these cumulates are significant components of the sources of alkaline OIB and continental magmas. The results of the forward models provide the potential for detailed tests of this class of hypotheses for the origin of alkaline magmas worldwide and for interpreting major and minor aspects of the geochemical variability of these magmas.

Access of extracellular cations to their binding sites in Na,K-ATPase: role of the second extracellular loop of the alpha subunit.

Na,K-ATPase, the main active transport system for monovalent cations in animal cells, is responsible for maintaining Na(+) and K(+) gradients across the plasma membrane. During its transport cycle it binds three cytoplasmic Na(+) ions and releases them on the extracellular side of the membrane, and then binds two extracellular K(+) ions and releases them into the cytoplasm. The fourth, fifth, and sixth transmembrane helices of the alpha subunit of Na,K-ATPase are known to be involved in Na(+) and K(+) binding sites, but the gating mechanisms that control the access of these ions to their binding sites are not yet fully understood. We have focused on the second extracellular loop linking transmembrane segments 3 and 4 and attempted to determine its role in gating. We replaced 13 residues of this loop in the rat alpha1 subunit, from E314 to G326, by cysteine, and then studied the function of these mutants using electrophysiological techniques. We analyzed the results using a structural model obtained by homology with SERCA, and ab initio calculations for the second extracellular loop. Four mutants were markedly modified by the sulfhydryl reagent MTSET, and we investigated them in detail. The substituted cysteines were more readily accessible to MTSET in the E1 conformation for the Y315C, W317C, and I322C mutants. Mutations or derivatization of the substituted cysteines in the second extracellular loop resulted in major increases in the apparent affinity for extracellular K(+), and this was associated with a reduction in the maximum activity. The changes produced by the E314C mutation were reversed by MTSET treatment. In the W317C and I322C mutants, MTSET also induced a moderate shift of the E1/E2 equilibrium towards the E1(Na) conformation under Na/Na exchange conditions. These findings indicate that the second extracellular loop must be functionally linked to the gating mechanism that controls the access of K(+) to its binding site.

Reciprocal regulation of brain and muscle Arnt-like protein 1 and peroxisome proliferator-activated receptor alpha defines a novel positive feedback loop in the rodent liver circadian clock.

Recent evidence has emerged that peroxisome proliferator-activated receptor alpha (PPARalpha), which is largely involved in lipid metabolism, can play an important role in connecting circadian biology and metabolism. In the present study, we investigated the mechanisms by which PPARalpha influences the pacemakers acting in the central clock located in the suprachiasmatic nucleus and in the peripheral oscillator of the liver. We demonstrate that PPARalpha plays a specific role in the peripheral circadian control because it is required to maintain the circadian rhythm of the master clock gene brain and muscle Arnt-like protein 1 (bmal1) in vivo. This regulation occurs via a direct binding of PPARalpha on a potential PPARalpha response element located in the bmal1 promoter. Reversely, BMAL1 is an upstream regulator of PPARalpha gene expression. We further demonstrate that fenofibrate induces circadian rhythm of clock gene expression in cell culture and up-regulates hepatic bmal1 in vivo. Together, these results provide evidence for an additional regulatory feedback loop involving BMAL1 and PPARalpha in peripheral clocks.

TNF-alpha-dependent loss of IKKbeta-deficient myeloid progenitors triggers a cytokine loop culminating in granulocytosis.

Loss of IκB kinase (IKK) β-dependent NF-κB signaling in hematopoietic cells is associated with increased granulopoiesis. Here we identify a regulatory cytokine loop that causes neutrophilia in Ikkβ-deficient mice. TNF-α-dependent apoptosis of myeloid progenitor cells leads to the release of IL-1β, which promotes Th17 polarization of peripheral CD4(+) T cells. Although the elevation of IL-17 and the consecutive induction of granulocyte colony-stimulating factor compensate for the loss of myeloid progenitor cells, the facilitated induction of Th17 cells renders Ikkβ-deficient animals more susceptible to the development of experimental autoimmune encephalitis. These results unravel so far unanticipated direct and indirect functions for IKKβ in myeloid progenitor survival and maintenance of innate and Th17 immunity and raise concerns about long-term IKKβ inhibition in IL-17-mediated diseases.

A positive FGFR3/FOXN1 feedback loop underlies benign skin keratosis versus squamous cell carcinoma formation in humans.

Seborrheic keratoses (SKs) are common, benign epithelial tumors of the skin that do not, or very rarely, progress into malignancy, for reasons that are not understood. We investigated this by gene expression profiling of human SKs and cutaneous squamous cell carcinomas (SCCs) and found that several genes previously connected with keratinocyte tumor development were similarly modulated in SKs and SCCs, whereas the expression of others differed by only a few fold. In contrast, the tyrosine kinase receptor FGF receptor-3 (FGFR3) and the transcription factor forkhead box N1 (FOXN1) were highly expressed in SKs, and close to undetectable in SCCs. We also showed that increased FGFR3 activity was sufficient to induce FOXN1 expression, counteract the inhibitory effect of EGFR signaling on FOXN1 expression and differentiation, and induce differentiation in a FOXN1-dependent manner. Knockdown of FOXN1 expression in primary human keratinocytes cooperated with oncogenic RAS in the induction of SCC-like tumors, whereas increased FOXN1 expression triggered the SCC cells to shift to a benign SK-like tumor phenotype, which included increased FGFR3 expression. Thus,we have uncovered a positive regulatory loop between FGFR3 and FOXN1 that underlies a benign versus malignant skin tumor phenotype.

Résistance aux diurétiques de l'anse en clinique [Resistance to loop diuretics in clinical practice]

Loop diuretics belong to the most common medications used in ambulatory and hospitalized patients, especially in situation of hypervolemia and chronic renal failure. Prolonged used of these agents in particular medical conditions can lead to attenuation of their diuretic effect, commonly known as "resistance" to diuretics. This article intends to review the main mechanisms of resistance to loop diuretic and the ways to counteract them in clinical practice.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Macrophage migration inhibitory factor is involved in a positive feedback loop increasing aromatase expression in endometriosis.

Immune-endocrine interplay may play a major role in the pathogenesis of endometriosis. In the present study, we have investigated the interaction between macrophage migration inhibitory factor (MIF), a major pro-inflammatory and growth-promoting factor markedly expressed in active endometriotic lesions, and estradiol (E(2)) in ectopic endometrial cells. Our data showed a significant increase of MIF protein secretion and mRNA expression in endometriotic cells in response to E(2). MIF production was blocked by Fulvestrant, an estrogen receptor (ER) antagonist, and induced by ERα and ERβ selective agonists propyl-pyrazole-triol (PPT) and diarylpropionrile (DPN), respectively, thus demonstrating a specific receptor-mediated effect. Cell transfection with MIF promoter construct showed that E(2) significantly stimulates MIF promoter activity. Interestingly, our data further revealed that MIF reciprocally stimulates aromatase protein and mRNA expression via a posttranscriptional mRNA stabilization mechanism, that E(2) itself can upregulate aromatase expression, and that inhibition of endogenous MIF, using MIF specific siRNA, significantly inhibits E(2)-induced aromatase. Thus, the present study revealed the existence of a local positive feedback loop by which estrogen acts directly on ectopic endometrial cells to upregulate the expression of MIF, which, in turn, displays the capability of inducing the expression of aromatase, the key and rate-limiting enzyme for estrogen synthesis. Such interplay may have a considerable impact on the development of endometriosis.

Simulations in evolution. II. Relative fitness and the propagation of mutants.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous article presented a histogram model [1] consisting in populations of individuals whose number changed under the influence of variation and/or fitness, the total population remaining constant. Individuals are classified into bins, and the content of each bin is calculated generation after generation by an Excel spreadsheet. Here, we apply the histogram model to a stable population with fitness F(1)=1.00 in which one or two fitter mutants emerge. In a first scenario, a single mutant emerged in the population whose fitness was greater than 1.00. The simulations ended when the original population was reduced to a single individual. The histogram model was validated by excellent agreement between its predictions and those of a classical continuous function (Eqn. 1) which predicts the number of generations needed for a favorable mutation to spread throughout a population. But in contrast to Eqn. 1, our histogram model is adaptable to more complex scenarios, as demonstrated here. In the second and third scenarios, the original population was present at time zero together with two mutants which differed from the original population by two higher and distinct fitness values. In the fourth scenario, the large original population was present at time zero together with one fitter mutant. After a number of generations, when the mutant offspring had multiplied, a second mutant was introduced whose fitness was even greater. The histogram model also allows Shannon entropy (SE) to be monitored continuously as the information content of the total population decreases or increases. The results of these simulations illustrate, in a graphically didactic manner, the influence of natural selection, operating through relative fitness, in the emergence and dominance of a fitter mutant.

A nucleosome-dependent static loop potentiates estrogen-regulated transcription from the Xenopus vitellogenin B1 promoter in vitro.

We describe the transcriptional potentiation in estrogen responsive transcription extracts of the Xenopus vitellogenin B1 gene promoter through the formation of a positioned nucleosome. Nuclease digestion and hydroxyl radical cleavage indicate that strong, DNA sequence-directed positioning of a nucleosome occurs between -300 and -140 relative to the start site of transcription. Deletion of this DNA sequence abolishes the potentiation of transcription due to nucleosome assembly. The wrapping of DNA around the histone core of the nucleosome positioned between -300 and -140 creates a static loop in which distal estrogen receptor binding sites are brought close to proximal promoter elements. This might facilitate interactions between the trans-acting factors themselves and/or RNA polymerase. Such a nucleosome provides an example of how chromatin structure might have a positive effect on the transcription process.

Cooperative kinking at distant sites in mechanically stressed DNA.

In cells, DNA is routinely subjected to significant levels of bending and twisting. In some cases, such as under physiological levels of supercoiling, DNA can be so highly strained, that it transitions into non-canonical structural conformations that are capable of relieving mechanical stress within the template. DNA minicircles offer a robust model system to study stress-induced DNA structures. Using DNA minicircles on the order of 100 bp in size, we have been able to control the bending and torsional stresses within a looped DNA construct. Through a combination of cryo-EM image reconstructions, Bal31 sensitivity assays and Brownian dynamics simulations, we have been able to analyze the effects of biologically relevant underwinding-induced kinks in DNA on the overall shape of DNA minicircles. Our results indicate that strongly underwound DNA minicircles, which mimic the physical behavior of small regulatory DNA loops, minimize their free energy by undergoing sequential, cooperative kinking at two sites that are located about 180° apart along the periphery of the minicircle. This novel form of structural cooperativity in DNA demonstrates that bending strain can localize hyperflexible kinks within the DNA template, which in turn reduces the energetic cost to tightly loop DNA.