Exact solutions of coupled multispecies linear reaction–diffusion equations on a uniformly growing domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.

Survival probability for a diffusive process on a growing domain

We consider the motion of a diffusive population on a growing domain, 0 < x < L(t ), which is motivated by various applications in developmental biology. Individuals in the diffusing population, which could represent molecules or cells in a developmental scenario, undergo two different kinds of motion: (i) undirected movement, characterized by a diffusion coefficient, D, and (ii) directed movement, associated with the underlying domain growth. For a general class of problems with a reflecting boundary at x = 0, and an absorbing boundary at x = L(t ), we provide an exact solution to the partial differential equation describing the evolution of the population density function, C(x,t ). Using this solution, we derive an exact expression for the survival probability, S(t ), and an accurate approximation for the long-time limit, S = limt→∞ S(t ). Unlike traditional analyses on a nongrowing domain, where S ≡ 0, we show that domain growth leads to a very different situation where S can be positive. The theoretical tools developed and validated in this study allow us to distinguish between situations where the diffusive population reaches the moving boundary at x = L(t ) from other situations where the diffusive population never reaches the moving boundary at x = L(t ). Making this distinction is relevant to certain applications in developmental biology, such as the development of the enteric nervous system (ENS). All theoretical predictions are verified by implementing a discrete stochastic model.

Exact solutions of linear reaction-diffusion on a uniformly growing domain: criteria for successful colonization

Many processes during embryonic development involve transport and reaction of molecules, or transport and proliferation of cells, within growing tissues. Mathematical models of such processes usually take the form of a reaction-diffusion partial differential equation (PDE) on a growing domain. Previous analyses of such models have mainly involved solving the PDEs numerically. Here, we present a framework for calculating the exact solution of a linear reaction-diffusion PDE on a growing domain. We derive an exact solution for a general class of one-dimensional linear reaction—diffusion process on 0domain. Comparing our exact solutions with numerical approximations confirms the veracity of the method. Furthermore, our examples illustrate a delicate interplay between: (i) the rate at which the domain elongates, (ii) the diffusivity associated with the spreading density profile, (iii) the reaction rate, and (iv) the initial condition. Altering the balance between these four features leads to different outcomes in terms of whether an initial profile, located near x = 0, eventually overcomes the domain growth and colonizes the entire length of the domain by reaching the boundary where x = L(t).

Transcriptional regulation by the orphan nuclear receptor ERRgamma

The nuclear receptor (NR) superfamily is comprised of receptors for small lipopfilic ligands such as steroid hormones, thyroid hormone, retinoids, and vitamin D. NRs are ligand-inducible transcription factors capable of both activating and repressing their target gene expression. They control a wide range of biological functions connected to growth, development, and homeostasis. In addition to the ligand-regulated receptors, the family includes a large group of receptors whose physiological ligands are unknown. These receptors are referred to as orphan NRs. Estrogen-related receptor gamma (ERRgamma) belongs to the ERR subfamily of orphan NRs together with the related ERRalpha and ERRbeta. ERRs share amino acid sequence homology with the classical estrogen receptors (ERs) but they are unable to bind natural estrogenic ligands. ERRgamma is expressed in several embryonic and adult tissues but its biological role is still largely unknown. ERRgamma activates reporter gene expression in transfected cells independently of added hormones implying that ERRgamma harbors constitutive activity. However, the intrinsic activity of ERRgamma can be inhibited by synthetic compounds such as the selective estrogen receptor modulator 4-hydroxytamoxifen (4-OHT). Ligands of NRs can act as agonists that activate transcription, as antagonists that prevent activation of transcription, or as inverse agonists that antagonize the constitutive transcriptional activity of receptor. Most of the synthetic ERRgamma ligands act as inverse agonists but recently, a synthetic ERRgamma agonist GSK4716 was identified. This demonstrates that it is possible to design and identify agonists for ERRgamma. Prior to this thesis work, the structural and functional characteristics of ERRgamma were largely unknown. The aim of this study was to define the functional requirements for ERRgamma-mediated transcriptional regulation and to examine the cross-talk between ERRgamma and other NRs. Due to the fact that natural physiological ligands of ERRgamma are unknown, another aim of this study was to seek new natural compounds that may affect transcriptional activity of ERRgamma. Plant-derived phytoestrogens have previously been shown to act as ligands for ERs and ERRalpha, and therefore the effects of these compounds were also studied on ERRgamma-mediated transcriptional regulation. This work demonstrated that ERRgamma-mediated transcriptional regulation was dependent on DNA-binding, dimerization and activation function-2. Heterodimerization with ERRalpha inhibited the transcriptional activity of ERRgamma. In addition to 4-OHT, another anti-estrogen, 4-hydroxytoremifene (4-OHtor), was identified as an inverse agonist of ERRgamma. Interestingly, ERRgamma activated transcription in the presence of 4-OHT and 4-OHtor on activator protein-1 binding sites. ERRgamma was found to interact with another orphan NR Nurr1 by repressing the ability of Nurr1 to activate transcription of the osteopontin gene. Transcriptional activity of ERRgamma was shown to be stimulated by the phytoestrogen equol. Structural model analysis and mutational experiments indicated that equol was able to bind to the ligand binding domain of ERRgamma. The growth inhibitory effect of ERRgamma on prostate cancer cells was found to be enhanced by equol. In summary, this study demonstrates that despite the absence of an endogenous physiological ligand, the activity of ERRgamma can be modulated in other ways such as dimerization with related receptors or by cross-talk with other transcription factors as well as by binding some synthetic or natural compounds.

Androgen receptor-dependent and independent functions of ARIP4

Androgen receptor (AR) is necessary for normal male phenotype development and essential for spermatogenesis. AR is a classical steroid receptor mediating actions of male sex steroids testosterone and 5-alpha-dihydrotestosterone. Numerous coregulators interact with the receptor and regulate AR activity on target genes. This study deals with the characterization of androgen receptor-interacting protein 4 (ARIP4). ARIP4 binds DNA, interacts with AR in vitro and in cultured yeast and mammalian cells, and modulates AR-dependent transactivation. ARIP4 is an active DNA-dependent ATPase, and this enzymatic activity is essential for the ability of ARIP4 to modulate AR function. On the basis of sequence homology in its ATPase domain, ARIP4 belongs to the SNF2 family of proteins involved in chromatin remodeling, DNA repair, and homologous recombination. Similar to its closest homologs ATRX and Rad54, ARIP4 does not seem to be a classical chromatin remodeling protein in that it does not appear to form large protein complexes in vivo or remodel mononucleosomes in vitro. However, ARIP4 is able to generate superhelical torsion on linear DNA fragments. ARIP4 is covalently modified by SUMO-1, and mutation of six potential SUMO attachment sites abolishes the ability of ARIP4 to bind DNA, hydrolyze ATP, and activate AR function. ARIP4 expression starts in early embryonic development. In mouse embryo ARIP4 is present mainly in the neural tube and limb buds. In adult mouse tissues ARIP4 expression is virtually ubiquitous. In mouse testis ARIP4 is expressed in the nuclei of Sertoli cells in a stage-dependent manner. ARIP4 is also present in the nuclei of Leydig cells, spermatogonia, pachytene and diplotene spermatocytes. Testicular expression pattern of ARIP4 does not differ significantly in wild-type, FSHRKO, and LuRKO mice. In the testis of hpg mice, ARIP4 is found mainly in interstitial cells and has very low, if any, expression in Sertoli and germ cells. Heterozygous Arip4+/ mice are fertile and appear normal; however, they are haploinsufficient with regard to androgen action in Sertoli cells. In contrast, Arip4 / embryos are not viable. They have significantly reduced body size at E9.5 and die by E11.5. Compared to wild-type littermates, Arip4 / embryos possess a higher percentage of apoptotic cells at E9.5 and E10.5. Fibroblasts derived from Arip4 / embryos cease growing after 2-3 passages and exhibit a significantly increased apoptosis and decreased proliferation rate than cells from wild-type embryos. Our findings demonstrate that ARIP4 plays an essential role in mouse embryonic development. In addition, testicular expression and AR coregulatory activity of ARIP4 suggest a role of ARIP4-AR interaction in the somatic cells of the testis.

Analysis-By-Synthesis Based Switched Transform Domain Split Vq Using Gaussian Mixture Model

Using analysis-by-synthesis (AbS) approach, we develop a soft decision based switched vector quantization (VQ) method for high quality and low complexity coding of wideband speech line spectral frequency (LSF) parameters. For each switching region, a low complexity transform domain split VQ (TrSVQ) is designed. The overall rate-distortion (R/D) performance optimality of new switched quantizer is addressed in the Gaussian mixture model (GMM) based parametric framework. In the AbS approach, the reduction of quantization complexity is achieved through the use of nearest neighbor (NN) TrSVQs and splitting the transform domain vector into higher number of subvectors. Compared to the current LSF quantization methods, the new method is shown to provide competitive or better trade-off between R/D performance and complexity.

NMR Spectroscopy of Multi-domain Proteins : Immunoglobulin-like Domains of Human Filamin A

Proteins are complex biomacromolecules playing fundamental roles in the physiological processes of all living organisms. They function as structural units, enzymes, transporters, process regulators, and signal transducers. Defects in protein functions often derive from genetic mutations altering the protein structure, and impairment of essential protein functions manifests itself as pathological conditions. Proteins operate through interactions, and all protein functions depend on protein structure. In order to understand biological mechanisms at the molecular level, one has to know the structures of the proteins involved. This thesis covers structural and functional characterization of human filamins. Filamins are actin-binding and -bundling proteins that have numerous interaction partners. In addition to their actin-organizing functions, filamins are also known to have roles in cell adhesion and locomotion, and to participate in the logistics of cell membrane receptors, and in the coordination of intracellular signaling pathways. Filamin mutations in humans induce severe pathological conditions affecting the brain, bones, limbs, and the cardiovascular system. Filamins are large modular proteins composed of an N-terminal actin-binding domain and 24 consecutive immunoglobulin-like domains (IgFLNs). Nuclear magnetic resonance (NMR) spectroscopy is a versatile method of gaining insight into protein structure, dynamics and interactions. NMR spectroscopy was employed in this thesis to study the atomic structure and interaction mechanisms of C-terminal IgFLNs, which are known to house the majority of the filamin interaction sites. The structures of IgFLN single-domains 17 and 23 and IgFLN domain pairs 16-17 and 18-19 were determined using NMR spectroscopy. The structures of domain pairs 16 17 and 18 19 both revealed novel domain domain interaction modes of IgFLNs. NMR titrations were employed to characterize the interactions of filamins with glycoprotein Ibα, FilGAP, integrin β7 and dopamine receptors. Domain packing of IgFLN domain sextet 16 21 was further characterized using residual dipolar couplings and NMR relaxation analysis. This thesis demonstrates the versatility and potential of NMR spectroscopy in structural and functional studies of multi-domain proteins.

Resistivity noise in crystalline magnetic nanowires and its implications to domain formation and kinetics

We have investigated the time-dependent fluctuations in electrical resistance, or noise, in high quality crystalline magnetic nanowires within nanoporous templates. The noise increases exponentially with increasing temperature and magnetic field, and has been analyzed in terms of domain wall depinning within the Neel-Brown framework. The frequency-dependence of noise also indicates a crossover from nondiffusive kinetics to long-range diffusion at higher temperatures, as well as a strong collective depinning, which need to be considered when implementing these nanowires in magnetoelectronic devices.

L2-stability of nonstationary feedback systems: Frequency-domain criteria

A frequency-domain positivity condition is derived for linear time-varying operators in2and is used to develop2stability criteria for linear and nonlinear feedback systems. These criteria permit the use of a very general class of operators in2with nonstationary kernels, as multipliers. More specific results are obtained by using a first-order differential operator with a time-varying coefficient as multiplier. Finally, by employing periodic multipliers, improved stability criteria are derived for the nonlinear damped Mathieu equation with a forcing function.

Time-domain criteria for theL_{2}-stability of nonstationary feedback systems

Criteria for the L2-stability of linear and nonlinear time-varying feedback systems are given. These are conditions in the time domain involving the solution of certain associated matrix Riccati equations and permitting the use of a very general class of L2-operators as multipliers.

Effect of fourier domain quantization on matched filtering

It is shown that the use of a coarsely quantized binary digital hologram as a matched filter on an optical computer does not degrade signal-to-noise ratio (SNR) appreciably.

Inhibition of Human Two-Pore Domain K+ Channel TREK1 by Local Anesthetic Lidocaine: Negative Cooperativity and Half-of-Sites Saturation Kinetics

TWIK-related K+ channel TREK1, a background leak K+ channel, has been strongly implicated as the target of several general and local anesthetics. Here, using the whole-cell and single-channel patch-clamp technique, we investigated the effect of lidocaine, a local anesthetic, on the human (h) TREK1 channel heterologously expressed in human embryonic kidney 293 cells by an adenoviral-mediated expression system. Lidocaine, at clinical concentrations, produced reversible, concentration-dependent inhibition of hTREK1 current, with IC50 value of 180 mu M, by reducing the single-channel open probability and stabilizing the closed state. We have identified a strategically placed unique aromatic couplet (Tyr352 and Phe355) in the vicinity of the protein kinase A phosphorylation site, Ser348, in the C-terminal domain (CTD) of hTREK1, that is critical for the action of lidocaine. Furthermore, the phosphorylation state of Ser348 was found to have a regulatory role in lidocaine-mediated inhibition of hTREK1. It is interesting that we observed strong intersubunit negative cooperativity (Hill coefficient = 0.49) and half-of-sites saturation binding stoichiometry (half-reaction order) for the binding of lidocaine to hTREK1. Studies with the heterodimer of wild-type (wt)-hTREK1 and Delta 119 C-terminal deletion mutant (hTREK1(wt)-Delta 119) revealed that single CTD of hTREK1 was capable of mediating partial inhibition by lidocaine, but complete inhibition necessitates the cooperative interaction between both the CTDs upon binding of lidocaine. Based on our observations, we propose a model that explains the unique kinetics and provides a plausible paradigm for the inhibitory action of lidocaine on hTREK1.

Time-Domain Measurement of Voltage Induction by Transient EM Fields

An efficient measurement technique for studying the effect of transient electromagnetic fields under controlled conditions has been described. Broad-band TEM fields with a rise-time of a few nanoseconds were generated using a stripline method. Theoretical results are obtained and experimental measurements which confirm these results are described. The work will form the basis for a study of the susceptibility of digital integrated circuits and their interconnections to transient electromagnetic fields.

Higher-resolution structure of the human insulin receptor ectodomain: Multi-modal inclusion of the insert domain

Insulin receptor (IR) signaling is critical to controlling nutrient uptake and metabolism. However, only a low-resolution (3.8 Å) structure currently exists for the IR ectodomain, with some segments ill-defined or unmodeled due to disorder. Here, we revise this structure using new diffraction data to 3.3 Å resolution that allow improved modeling of the N-linked glycans, the first and third fibronectin type III domains, and the insert domain. A novel haptic interactive molecular dynamics strategy was used to aid fitting to low-resolution electron density maps. The resulting model provides a foundation for investigation of structural transitions in IR upon ligand binding.