On the application of meshfree methods to the nonlinear dynamics of multibody systems

Esta tesis propone una completa formulación termo-mecánica para la simulación no-lineal de mecanismos flexibles basada en métodos libres de malla. El enfoque se basa en tres pilares principales: la formulación de Lagrangiano total para medios continuos, la discretización de Bubnov-Galerkin, y las funciones de forma libres de malla. Los métodos sin malla se caracterizan por la definición de un conjunto de funciones de forma en dominios solapados, junto con una malla de integración de las ecuaciones discretas de balance. Dos tipos de funciones de forma se han seleccionado como representación de las familias interpolantes (Funciones de Base Radial) y aproximantes (Mínimos Cuadrados Móviles). Su formulación se ha adaptado haciendo sus parámetros compatibles, y su ausencia de conectividad predefinida se ha aprovechado para interconectar múltiples dominios de manera automática, permitiendo el uso de mallas de fondo no conformes. Se propone una formulación generalizada de restricciones, juntas y contactos, válida para sólidos rígidos y flexibles, siendo estos últimos discretizados mediante elementos finitos (MEF) o libres de malla. La mayor ventaja de este enfoque reside en que independiza completamente el dominio con respecto de las uniones y acciones externas a cada sólido, permitiendo su definición incluso fuera del contorno. Al mismo tiempo, también se minimiza el número de ecuaciones de restricción necesarias para la definición de uniones realistas. Las diversas validaciones, ejemplos y comparaciones detalladas muestran como el enfoque propuesto es genérico y extensible a un gran número de sistemas. En concreto, las comparaciones con el MEF indican una importante reducción del error para igual número de nodos, tanto en simulaciones mecánicas, como térmicas y termo-mecánicas acopladas. A igualdad de error, la eficiencia numérica de los métodos libres de malla es mayor que la del MEF cuanto más grosera es la discretización. Finalmente, la formulación se aplica a un problema de diseño real sobre el mantenimiento de estructuras masivas en el interior de un reactor de fusión, demostrando su viabilidad en análisis de problemas reales, y a su vez mostrando su potencial para su uso en simulación en tiempo real de sistemas no-lineales. A new complete formulation is proposed for the simulation of nonlinear dynamic of multibody systems with thermo-mechanical behaviour. The approach is founded in three main pillars: total Lagrangian formulation, Bubnov-Galerkin discretization, and meshfree shape functions. Meshfree methods are characterized by the definition of a set of shape functions in overlapping domains, and a background grid for integration of the Galerkin discrete equations. Two different types of shape functions have been chosen as representatives of interpolation (Radial Basis Functions), and approximation (Moving Least Squares) families. Their formulation has been adapted to use compatible parameters, and their lack of predefined connectivity is used to interconnect different domains seamlessly, allowing the use of non-conforming meshes. A generalized formulation for constraints, joints, and contacts is proposed, which is valid for rigid and flexible solids, being the later discretized using either finite elements (FEM) or meshfree methods. The greatest advantage of this approach is that makes the domain completely independent of the external links and actions, allowing to even define them outside of the boundary. At the same time, the number of constraint equations needed for defining realistic joints is minimized. Validation, examples, and benchmarks are provided for the proposed formulation, demonstrating that the approach is generic and extensible to further problems. Comparisons with FEM show a much lower error for the same number of nodes, both for mechanical and thermal analyses. The numerical efficiency is also better when coarse discretizations are used. A final demonstration to a real problem for handling massive structures inside of a fusion reactor is presented. It demonstrates that the application of meshfree methods is feasible and can provide an advantage towards the definition of nonlinear real-time simulation models.

Energy-entropy-momentum time integration methods for coupled smooth dissipative problems

Esta tesis aborda la formulación, análisis e implementación de métodos numéricos de integración temporal para la solución de sistemas disipativos suaves de dimensión finita o infinita de manera que su estructura continua sea conservada. Se entiende por dichos sistemas aquellos que involucran acoplamiento termo-mecánico y/o efectos disipativos internos modelados por variables internas que siguen leyes continuas, de modo que su evolución es considerada suave. La dinámica de estos sistemas está gobernada por las leyes de la termodinámica y simetrías, las cuales constituyen la estructura que se pretende conservar de forma discreta. Para ello, los sistemas disipativos se describen geométricamente mediante estructuras metriplécticas que identifican claramente las partes reversible e irreversible de la evolución del sistema. Así, usando una de estas estructuras conocida por las siglas (en inglés) de GENERIC, la estructura disipativa de los sistemas es identificada del mismo modo que lo es la Hamiltoniana para sistemas conservativos. Con esto, métodos (EEM) con precisión de segundo orden que conservan la energía, producen entropía y conservan los impulsos lineal y angular son formulados mediante el uso del operador derivada discreta introducido para asegurar la conservación de la Hamiltoniana y las simetrías de sistemas conservativos. Siguiendo estas directrices, se formulan dos tipos de métodos EEM basados en el uso de la temperatura o de la entropía como variable de estado termodinámica, lo que presenta importantes implicaciones que se discuten a lo largo de esta tesis. Entre las cuales cabe destacar que las condiciones de contorno de Dirichlet son naturalmente impuestas con la formulación basada en la temperatura. Por último, se validan dichos métodos y se comprueban sus mejores prestaciones en términos de la estabilidad y robustez en comparación con métodos estándar. This dissertation is concerned with the formulation, analysis and implementation of structure-preserving time integration methods for the solution of the initial(-boundary) value problems describing the dynamics of smooth dissipative systems, either finite- or infinite-dimensional ones. Such systems are understood as those involving thermo-mechanical coupling and/or internal dissipative effects modeled by internal state variables considered to be smooth in the sense that their evolutions follow continuos laws. The dynamics of such systems are ruled by the laws of thermodynamics and symmetries which constitutes the structure meant to be preserved in the numerical setting. For that, dissipative systems are geometrically described by metriplectic structures which clearly identify the reversible and irreversible parts of their dynamical evolution. In particular, the framework known by the acronym GENERIC is used to reveal the systems' dissipative structure in the same way as the Hamiltonian is for conserving systems. Given that, energy-preserving, entropy-producing and momentum-preserving (EEM) second-order accurate methods are formulated using the discrete derivative operator that enabled the formulation of Energy-Momentum methods ensuring the preservation of the Hamiltonian and symmetries for conservative systems. Following these guidelines, two kind of EEM methods are formulated in terms of entropy and temperature as a thermodynamical state variable, involving important implications discussed throughout the dissertation. Remarkably, the formulation in temperature becomes central to accommodate Dirichlet boundary conditions. EEM methods are finally validated and proved to exhibit enhanced numerical stability and robustness properties compared to standard ones.

The product of a thyroid hormone-responsive gene interacts with thyroid hormone receptors

Thyroid hormone is a critical mediator of central nervous system (CNS) development, acting through nuclear receptors to modulate the expression of specific genes. Transcription of the rat hairless (hr) gene is highly up-regulated by thyroid hormone in the developing CNS; we show here that hr is directly induced by thyroid hormone. By identifying proteins that interact with the hr gene product (Hr), we find that Hr interacts directly and specifically with thyroid hormone receptor (TR)—the same protein that regulates its expression. Unlike previously described receptor-interacting factors, Hr associates with TR and not with retinoic acid receptors (RAR, RXR). Hr can act as a transcriptional repressor, suggesting that its interaction with TR is part of a novel autoregulatory mechanism.

A calcium responsive element that regulates expression of two calcium binding proteins in Purkinje cells

Calbindin D28 encodes a calcium binding protein that is expressed in the cerebellum exclusively in Purkinje cells. We have used biolistic transfection of organotypic slices of P12 cerebellum to identify a 40-bp element from the calbindin promoter that is necessary and sufficient for Purkinje cell specific expression in this transient in situ assay. This element (PCE1) is also present in the calmodulin II promoter, which regulates expression of a second Purkinje cell Ca2+ binding protein. Expression of high levels of exogenous calbindin or calretinin decreased transcription mediated by PCE1 in Purkinje cells 2.5- to 3-fold, whereas the presence of 1 μM ionomycin in the extracellular medium increased expression. These results demonstrate that PCE1 is a component of a cell-specific and Ca2+-sensitive transcriptional regulatory mechanism that may play a key role in setting the Ca2+ buffering capacity of Purkinje cells.

Regulation of chloroplast protein import through a protochlorophyllide-responsive transit peptide

NADPH:protochlorophyllide (Pchlide) oxidoreductase (POR) is the key enzyme of chlorophyll biosynthesis in angiosperms. In barley, two POR enzymes, termed PORA and PORB, exist. Both are nucleus-encoded plastid proteins that must be imported posttranslationally from the cytosol. Whereas the import of the precursor of PORA, pPORA, previously has been shown to depend on Pchlide, the import of pPORB occurred constitutively. To study this striking difference, chimeric precursor proteins were constructed in which the transit sequences of the pPORA and pPORB were exchanged and fused to either their cognate polypeptides or to a cytosolic dihydrofolate reductase (DHFR) reporter protein of mouse. As shown here, the transit peptide of the pPORA (transA) conferred the Pchlide requirement of import onto both the mature PORB and the DHFR. By contrast, the transit peptide of the pPORB directed the reporter protein into both chloroplasts that contained or lacked translocation-active Pchlide. In vitro binding studies further demonstrated that the transit peptide of the pPORA, but not of the pPORB, is able to bind Pchlide. We conclude that the import of the authentic pPORA and that of the transA-PORB and transA-DHFR fusion proteins is regulated by a direct transit peptide-Pchlide interaction, which is likely to occur in the plastid envelope, a major site of porphyrin biosynthesis.

The problem of the spreading of a liquid film along a solid surface: A new mathematical formulation

A new mathematical model is proposed for the spreading of a liquid film on a solid surface. The model is based on the standard lubrication approximation for gently sloping films (with the no-slip condition for the fluid at the solid surface) in the major part of the film where it is not too thin. In the remaining and relatively small regions near the contact lines it is assumed that the so-called autonomy principle holds—i.e., given the material components, the external conditions, and the velocity of the contact lines along the surface, the behavior of the fluid is identical for all films. The resulting mathematical model is formulated as a free boundary problem for the classical fourth-order equation for the film thickness. A class of self-similar solutions to this free boundary problem is considered.

Functional antioxidant responsive elements

Exposure of human and rodent cells to a wide variety of chemoprotective compounds confers resistance against a broad set of carcinogens. For a subset of the chemoprotective compounds, protection is generated by an increase in the abundance of protective enzymes like glutathione S-transferases (GST). Antioxidant responsive elements (AREs) mediate the transcriptional induction of a battery of genes which comprise much of this chemoprotective response system. Past studies identified a necessary ARE “core” sequence of RTGACnnnGC, but this sequence alone is insufficient to mediate induction. In this study, the additional sequences necessary to define a sufficient, functional ARE are identified through systematic mutational analysis of the murine GST Ya ARE. Introduction of the newly identified necessary nucleotides into the regions flanking a nonresponsive, ARE-like, GST-Mu promoter sequence produced an inducible element. A screen of the GenBank database with the newly identified ARE consensus identified 16 genes which contained the functional ARE consensus sequence in their promoters. Included within this group was an ARE sequence from the murine ferritin-L promoter that mediated induction when tested. In an electrophoretic mobility-shift assay, the ferritin-L ARE was bound by ARE–binding protein 1, a protein previously identified as the likely mediator of the chemoprotective response. A three-level ARE classification system is presented to account for the distinct induction strengths observed in our mutagenesis studies. A model of the ARE as a composite regulatory site, where multiple transcription factors interact, is presented to account for the complex characteristics of ARE-mediated chemoprotective gene expression.

The Formation of an Insulin-responsive Vesicular Cargo Compartment Is an Early Event in 3T3-L1 Adipocyte Differentiation

Differentiating 3T3-L1 cells exhibit a dramatic increase in the rate of insulin-stimulated glucose transport during their conversion from proliferating fibroblasts to nonproliferating adipocytes. On day 3 of 3T3-L1 cell differentiation, basal glucose transport and cell surface transferrin binding are markedly diminished. This occurs concomitant with the formation of a distinct insulin-responsive vesicular pool of intracellular glucose transporter 1 (GLUT1) and transferrin receptors as assessed by sucrose velocity gradients. The intracellular distribution of the insulin-responsive aminopeptidase is first readily detectable on day 3, and its gradient profile and response to insulin at this time are identical to that of GLUT1. With further time of differentiation, GLUT4 is expressed and targeted to the same insulin-responsive vesicles as the other three proteins. Our data are consistent with the notion that a distinct insulin-sensitive vesicular cargo compartment forms early during fat call differentiation and its formation precedes GLUT4 expression. The development of this compartment may result from the differentiation-dependent inhibition of constitutive GLUT1 and transferrin receptor trafficking such that there is a large increase in, or the new formation of, a population of postendosomal, insulin-responsive vesicles.

Differential Regulation of Secretory Compartments Containing the Insulin-responsive Glucose Transporter 4 in 3T3-L1 Adipocytes

Insulin and guanosine-5′-O-(3-thiotriphosphate) (GTPγS) both stimulate glucose transport and translocation of the insulin-responsive glucose transporter 4 (GLUT4) to the plasma membrane in adipocytes. Previous studies suggest that these effects may be mediated by different mechanisms. In this study we have tested the hypothesis that these agonists recruit GLUT4 by distinct trafficking mechanisms, possibly involving mobilization of distinct intracellular compartments. We show that ablation of the endosomal system using transferrin-HRP causes a modest inhibition (∼30%) of insulin-stimulated GLUT4 translocation. In contrast, the GTPγS response was significantly attenuated (∼85%) under the same conditions. Introduction of a GST fusion protein encompassing the cytosolic tail of the v-SNARE cellubrevin inhibited GTPγS-stimulated GLUT4 translocation by ∼40% but had no effect on the insulin response. Conversely, a fusion protein encompassing the cytosolic tail of vesicle-associated membrane protein-2 had no significant effect on GTPγS-stimulated GLUT4 translocation but inhibited the insulin response by ∼40%. GTPγS- and insulin-stimulated GLUT1 translocation were both partially inhibited by GST-cellubrevin (∼50%) but not by GST-vesicle-associated membrane protein-2. Incubation of streptolysin O-permeabilized 3T3-L1 adipocytes with GTPγS caused a marked accumulation of Rab4 and Rab5 at the cell surface, whereas other Rab proteins (Rab7 and Rab11) were unaffected. These data are consistent with the localization of GLUT4 to two distinct intracellular compartments from which it can move to the cell surface independently using distinct sets of trafficking molecules.

Tuning of activation thresholds explains flexibility in the selection and development of T cells in the thymus

Immature CD4+CD8+ thymocytes expressing T-cell antigen receptors (TCR) are selected by TCR-mediated recognition of peptides associated with major histocompatibility complex molecules on thymic stromal cells. Selection ensures reactivity of the mature cells to foreign antigens and tolerance to self. Although much has been learned about the factors that determine whether a thymocyte with a given specificity will be positively or negatively selected, selection as an aspect of the developmental process as a whole is less well-understood. Here we invoke a model in which thymocytes tune their response characteristics individually and dynamically in the course of development. Cellular development and selection are driven by receptor-mediated metabolic perturbations. Perturbation is a measure of the net intracellular change induced by external stimulation. It results from the integration of several signals and countersignals over time and therefore depends on the environment and the maturation stage of the cell. Individual cell adaptation limits the range of perturbations. Such adaptation renders thymocytes less sensitive to the level of stimulation per se, but responsive to environmental changes in that level. This formulation begins to explain the mechanisms that link developmental and selection events to each other.

Myeloblastin is a granulocyte colony-stimulating factor-responsive gene conferring factor-independent growth to hematopoietic cells

Hematopoiesis depends on a pool of quiescent hematopoietic stem/progenitor cells. When exposed to specific cytokines, a portion of these cells enters the cell cycle to generate an amplified progeny. Myeloblastin (MBN) initially was described as involved in proliferation of human leukemia cells. The granulocyte colony-stimulating factor (G-CSF), which stimulates the proliferation of granulocytic precursors, up-regulates MBN expression. Here we show that constitutive overexpression of MBN confers factor-independent growth to murine bone marrow-derived Ba/F3/G-CSFR cells. Our results point to MBN as a G-CSF responsive gene critical to factor-independent growth and indicate that expression of the G-CSF receptor is a prerequisite to this process. A 91-bp MBN promoter region containing PU.1, C/EBP, and c-Myb binding sites is responsive to G-CSF treatment. Although PU.1, C/EBP, and c-Myb transcription factors all were critical for expression of MBN, its up-regulation by G-CSF was associated mainly with PU.1. These findings suggest that MBN is an important target of PU.1 and a key protease for factor-independent growth of hematopoietic cells.

IFN consensus sequence binding protein potentiates STAT1-dependent activation of IFNγ-responsive promoters in macrophages

IFNγ, once called the macrophage-activating factor, stimulates many genes in macrophages, ultimately leading to the elicitation of innate immunity. IFNγ's functions depend on the activation of STAT1, which stimulates transcription of IFNγ-inducible genes through the GAS element. The IFN consensus sequence binding protein (icsbγ or IFN regulatory factor 8), encoding a transcription factor of the IFN regulatory factor family, is one of such IFNγ-inducible genes in macrophages. We found that macrophages from ICSBP−/− mice were defective in inducing some IFNγ-responsive genes, even though they were capable of activating STAT1 in response to IFNγ. Accordingly, IFNγ activation of luciferase reporters fused to the GAS element was severely impaired in ICSBP−/− macrophages, but transfection of ICSBP resulted in marked stimulation of these reporters. Consistent with its role in activating IFNγ-responsive promoters, ICSBP stimulated reporter activity in a GAS-specific manner, even in the absence of IFNγ treatment, and in STAT1 negative cells. Indicative of a mechanism for this stimulation, DNA affinity binding assays revealed that endogenous ICSBP was recruited to a multiprotein complex that bound to GAS. These results suggest that ICSBP, when induced by IFNγ through STAT1, in turn generates a second wave of transcription from GAS-containing promoters, thereby contributing to the elicitation of IFNγ's unique activities in immune cells.

The biological clock of very premature primate infants is responsive to light

Each year more than 250,000 infants in the United States are exposed to artificial lighting in hospital nurseries with little consideration given to environmental lighting cycles. Essential in determining whether environmental lighting cycles need to be considered in hospital nurseries is identifying when the infant’s endogenous circadian clock becomes responsive to light. Using a non-human primate model of the developing human, we examined when the circadian clock, located in the hypothalamic suprachiasmatic nuclei (SCN), becomes responsive to light. Preterm infant baboons of different ages were exposed to light (5,000 lux) at night, and then changes in SCN metabolic activity and gene expression were assessed. After exposure to bright light at night, robust increases in SCN metabolic activity and gene expression were seen at ages that were equivalent to human infants at 24 weeks after conception. These data provide direct evidence that the biological clock of very premature primate infants is responsive to light.

Tonicity-responsive enhancer binding protein, a Rel-like protein that stimulates transcription in response to hypertonicity

Hypertonicity (most often present as high salinity) is stressful to the cells of virtually all organisms. Cells survive in a hypertonic environment by increasing the transcription of genes whose products catalyze cellular accumulation of compatible osmolytes. In mammals, the kidney medulla is normally hypertonic because of the urinary concentrating mechanism. Cellular accumulation of compatible osmolytes in the renal medulla is catalyzed by the sodium/myo-inositol cotransporter (SMIT), the sodium/chloride/betaine cotransporter, and aldose reductase (synthesis of sorbitol). The importance of compatible osmolytes is underscored by the necrotic injury of the renal medulla and subsequent renal failure that results from the inhibition of SMIT in vivo by administration of a specific inhibitor. Tonicity-responsive enhancers (TonE) play a key role in hypertonicity-induced transcriptional stimulation of SMIT, sodium/chloride/betaine cotransporter, and aldose reductase. We report the cDNA cloning of human TonE binding protein (TonEBP), a transcription factor that stimulates transcription through its binding to TonE sequences via a Rel-like DNA binding domain. Western blot and immunohistochemical analyses of cells cultured in hypertonic medium reveal that exposure to hypertonicity elicits slow activation of TonEBP, which is the result of an increase in TonEBP amount and translocation to the nucleus.

LXRα functions as a cAMP-responsive transcriptional regulator of gene expression

LXRα is a member of a nuclear receptor superfamily that regulates transcription. LXRα forms a heterodimer with RXRα, another member of this family, to regulate the expression of cholesterol 7α-hydroxylase by means of binding to the DR4-type cis-element. Here, we describe a function for LXRα as a cAMP-responsive regulator of renin and c-myc gene transcriptions by the interaction with a specific cis-acting DNA element, CNRE (an overlapping cAMP response element and a negative response element). Our previous studies showed that renin gene expression is regulated by cAMP, at least partly, through the CNRE sequence in its 5′-flanking region. This sequence is also found in c-myc and several other genes. Based on our cloning results using the yeast one-hybrid system, we discovered that the mouse homologue of human LXRα binds to the CNRE and demonstrated that it binds as a monomer. To define the function of LXRα on gene expression, we transfected the renin-producing renal As4.1 cells with LXRα expression plasmid. Overexpression of LXRα in As4.1 cells confers cAMP inducibility to reporter constructs containing the renin CNRE. After stable transfection of LXRα, As4.1 cells show a cAMP-inducible up-regulation of renin mRNA expression. In parallel experiments, we demonstrated that LXRα can also bind to the homologous CNRE in the c-myc promoter. cAMP promotes transcription through c-myc/CNRE:LXRα interaction in LXRα transiently transfected cells and increases c-myc mRNA expression in stably transfected cells. Identification of LXRα as a cAMP-responsive nuclear modulator of renin and c-myc expression not only has cardiovascular significance but may have generalized implication in the regulation of gene transcription.