Physiologically based synthetic models of hepatic disposition

Current Physiologically based pharmacokinetic (PBPK) models are inductive. We present an additional, different approach that is based on the synthetic rather than the inductive approach to modeling and simulation. It relies on object-oriented programming A model of the referent system in its experimental context is synthesized by assembling objects that represent components such as molecules, cells, aspects of tissue architecture, catheters, etc. The single pass perfused rat liver has been well described in evaluating hepatic drug pharmacokinetics (PK) and is the system on which we focus. In silico experiments begin with administration of objects representing actual compounds. Data are collected in a manner analogous to that in the referent PK experiments. The synthetic modeling method allows for recognition and representation of discrete event and discrete time processes, as well as heterogeneity in organization, function, and spatial effects. An application is developed for sucrose and antipyrine, administered separately and together PBPK modeling has made extensive progress in characterizing abstracted PK properties but this has also been its limitation. Now, other important questions and possible extensions emerge. How are these PK properties and the observed behaviors generated? The inherent heuristic limitations of traditional models have hindered getting meaningful, detailed answers to such questions. Synthetic models of the type described here are specifically intended to help answer such questions. Analogous to wet-lab experimental models, they retain their applicability even when broken apart into sub-components. Having and applying this new class of models along with traditional PK modeling methods is expected to increase the productivity of pharmaceutical research at all levels that make use of modeling and simulation.

Theory of pseudomodes in quantum optical processes

This paper deals with non-Markovian behavior in atomic systems coupled to a structured reservoir of quantum electromagnetic field modes, with particular relevance to atoms interacting with the field in high-Q cavities or photonic band-gap materials. In cases such as the former, we show that the pseudomode theory for single-quantum reservoir excitations can be obtained by applying the Fano diagonalization method to a system in which the atomic transitions are coupled to a discrete set of (cavity) quasimodes, which in turn are coupled to a continuum set of (external) quasimodes with slowly varying coupling constants and continuum mode density. Each pseudomode can be identified with a discrete quasimode, which gives structure to the actual reservoir of true modes via the expressions for the equivalent atom-true mode coupling constants. The quasimode theory enables cases of multiple excitation of the reservoir to now be treated via Markovian master equations for the atom-discrete quasimode system. Applications of the theory to one, two, and many discrete quasimodes are made. For a simple photonic band-gap model, where the reservoir structure is associated with the true mode density rather than the coupling constants, the single quantum excitation case appears to be equivalent to a case with two discrete quasimodes.

Non-Markovian quantum optical processes and pseudomodes

Non-Markovian behaviour in atomic systems coupled to a structured reservoir of quantum EM field modes, such as in high Q cavities, is treated using a quasimode description, and the pseudo mode theory for single quantum reservoir excitations is obtained via Fano diagonalisation. The atomic transitions are coupled to a discrete set of (cavity) quasimodes, which are also coupled to a continuum set of (external) quasimodes with slowly varying coupling constants. Each pseudomode corresponds to a cavity quasimode, and the original reservoir structure is obtained in expressions for the equivalent atom-true mode coupling constants. Cases of multiple excitation of the reservoir are now treatable via Markovian master equations for the atom-discrete quasimode system.

Set-valued Markov chains and negative semitrajectories of discretized dynamical systems

Computer simulation of dynamical systems involves a phase space which is the finite set of machine arithmetic. Rounding state values of the continuous system to this grid yields a spatially discrete dynamical system, often with different dynamical behaviour. Discretization of an invertible smooth system gives a system with set-valued negative semitrajectories. As the grid is refined, asymptotic behaviour of the semitrajectories follows probabilistic laws which correspond to a set-valued Markov chain, whose transition probabilities can be explicitly calculated. The results are illustrated for two-dimensional dynamical systems obtained by discretization of fractional linear transformations of the unit disc in the complex plane.

The Bazhanov-Stroganov model from 3D approach

We apply a three-dimensional approach to describe a new parametrization of the L-operators for the two-dimensional Bazhanov-Stroganov (BS) integrable spin model related to the chiral Potts model. This parametrization is based on the solution of the associated classical discrete integrable system. Using a three-dimensional vertex satisfying a modified tetrahedron equation, we construct an operator which generalizes the BS quantum intertwining matrix S. This operator describes the isospectral deformations of the integrable BS model.

Probabilistic Action System Trace Semantics

Action systems are a construct for reasoning about concurrent, reactive systems, in which concurrent behaviour is described by interleaving atomic actions. Sere and Troubitsyna have proposed an extension to action systems in which actions may be expressed and composed using discrete probabilistic choice as well as demonic nondeterministic choice. In this paper we develop a trace-based semantics for probabilistic action systems. This semantics provides a simple theoretical base on which practical refinement rules for probabilistic action systems may be justified.

Expression of galectin-1 in the mouse olfactory system

Primary sensory olfactory axons arise from the olfactory neuroepithelium that lines the nasal cavity and then project via the olfactory nerve into the olfactory bulb. The P-galactoside binding lectin, galectin-1,and its laminin ligand have been implicated in the growth of these axons along this pathway. In galectin-1 null mutant mice, a subpopulation of primary sensory olfactory axons fails to reach its targets in the olfactory bulb. In the present study we examined the spatiotemporal expression pattern of galectin-1 in normal mice in order to understand its role in the development of the olfactory nerve pathway. At E15.5, when olfactory axons have already contacted the olfactory bulb, galectin-1 was expressed in the cartilage and mesenchyme surrounding the nasal cavity but was absent from the olfactory neuroepithelium, nerve and bulb. Between E16.5 and birth galectin-1 began to be expressed by olfactory nerve ensheathing cells in the lamina propria of the neuroepithelium and nerve fibre layer. Galectin-1 was neither expressed by primary sensory neurons in the olfactory neuroepithelium nor by their axons in the olfactory nerve. Laminin, a galectin-1 ligand, also exhibited a similar expression pattern in the embryonic olfactory nerve pathway. Our results reveal that galectin-1 is dynamically expressed by glial elements within the nerve fibre layer during a discrete period in the developing olfactory nerve pathway. Previous studies have reported galectin-1 acts as a substrate adhesion molecule by cross-linking primary sensory olfactory neurons to laminin. Thus, the coordinate expression of galectin-1 and laminin in the embryonic nerve fibre layer suggests that these molecules support the adhesion and fasciculation of axons en route to their glomerular targets.

Satellite observation and climate system model simulation of the St. Lawrence Island polynya

The St. Lawrence Island polynya (SLIP) is a commonly occurring winter phenomenon in the Bering Sea, in which dense saline water produced during new ice formation is thought to flow northward through the Bering Strait to help maintain the Arctic Ocean halocline. Winter darkness and inclement weather conditions have made continuous in situ and remote observation of this polynya difficult. However, imagery acquired from the European Space Agency ERS-1 Synthetic Aperture Radar (SAR) has allowed observation of the St. Lawrence Island polynya using both the imagery and derived ice displacement products. With the development of ARCSyM, a high resolution regional model of the Arctic atmosphere/sea ice system, simulation of the SLIP in a climate model is now possible. Intercomparisons between remotely sensed products and simulations can lead to additional insight into the SLIP formation process. Low resolution SAR, SSM/I and AVHRR infrared imagery for the St. Lawrence Island region are compared with the results of a model simulation for the period of 24-27 February 1992. The imagery illustrates a polynya event (polynya opening). With the northerly winds strong and consistent over several days, the coupled model captures the SLIP event with moderate accuracy. However, the introduction of a stability dependent atmosphere-ice drag coefficient, which allows feedbacks between atmospheric stability, open water, and air-ice drag, produces a more accurate simulation of the SLIP in comparison to satellite imagery. Model experiments show that the polynya event is forced primarily by changes in atmospheric circulation followed by persistent favorable conditions: ocean surface currents are found to have a small but positive impact on the simulation which is enhanced when wind forcing is weak or variable.

Role of galectin-1 in the olfactory nervous system

The olfactory nervous system is responsible for the detection of odors. Primary sensory olfactory neurons are located in a neuroepithelial sheet lining the nasal cavity. The axons from these neurons converge on to discrete loci or glomeruli in the olfactory bulb. Each glomerulus consists of the termination of thousands of primary axons on the dendrites of second-order olfactory neurons. What are the molecular mechanisms which guide growing olfactory axons to select sites in the olfactory bulb? We have shown that subpopulations of these axons differentially express cell surface carbohydrates and that these different subpopulations target and terminate in particular regions of the olfactory bulb. Interestingly, the olfactory neurons and glial components in the olfactory pathway between the nose and brain express galectin-1. By using in vitro assays of neurite outgrowth we found that both galectin-1 and it's ligands were capable of specifically stimulating neurite elongation. Examination of the olfactory system in galectin-1 null mutants revealed that a subpopulation of axons failed to navigate to their target site in the olfactory bulb. This is the first phenotypic effect observed in galectin-1 null mutants and indicates that galectin-1 has a role in the growth and/or guidance of a subpopulation of axons in the olfactory system during development.

Presence of novel N-CAM glycoforms in the rat olfactory system

The functional activity of the neural cell adhesion molecule N-CAM can be modulated by posttranslational modifications such as glycosylation. For instance, the long polysialic acid side chains of N-CAM alter the adhesion properties of the protein backbone. In the present study, we identified two novel carbohydrates present on N-CAM, NOC-3 and NOC-4. Both carbohydrates were detected on N-CAM glycoforms expressed by subpopulations of primary sensory olfactory neurons in the rat olfactory system. Based on the expression of NOC-3 and NOC-4 and the olfactory marker protein (OMP), four independent subpopulations of primary sensory olfactory neurons were characterized. These neurons expressed: both NOC-3 and NOC-4 but not OMP; both NOC-4 and OMP but not NOC-3; NOC-3, NOC-4, and OMP together; and OMP alone. The NOC-3- and NOC-4-expressing neurons were widely dispersed in the olfactory neuroepithelium lining the nasal cavity. The axons of NOC-4 expressing neurons innervated all glomeruli in the olfactory bulb, whereas the NOC-3 expressing axons terminated in a discrete subset of glomeruli scattered throughout the whole olfactory bulb. We propose that both NOC-3 and NOC-4 are part of a chemical code of olfactory neurons which is used in establishing the topography of connections between the olfactory neuroepithelium and the olfactory bulb. (C) 1997 John Wiley & Sons, Inc.

Egr transcription factors in the nervous system

The Egr proteins, Egr-1, Egr-2, Egr-3 and Egr-4, are closely related members of a subclass of immediate early gene-encoded, inducible transcription factors. They share a highly homologous DNA-binding domain which recognises an identical DNA response element. In addition, they have several less-well conserved structural features in common. As immediate early proteins, the Egr transcription factors are rapidly induced by diverse extracellular stimuli within the nervous system in a discretely controlled manner. The basal expression of the Egr proteins in the developing and adult rat brain and the induction of Egr proteins by neurotransmitter analogue stimulation, physiological mimetic and brain injury paradigms is reviewed. We review evidence indicating that Egr proteins are subject to tight differential control through diverse mechanisms at several levels of regulation. These include transcriptional, translational and posttranslational (including glycosylation, phosphorylation and redox) mechanisms and protein-protein interaction. Ultimately the differentially co-ordinated Egr response may lead to discrete effects on target gene expression. Some of the known target genes of Egr proteins and functions of the Egr proteins in different cell types are also highlighted. Future directions for research into the control and function of the different Egr proteins are also explored. (C) 1997 Elsevier Science Ltd.

Optimal quantum trajectories for discrete measurements

Using the method of quantum trajectories we show that a known pure state can be optimally monitored through time when subject to a sequence of discrete measurements. By modifying the way that we extract information from the measurement apparatus we can minimize the average algorithmic information of the measurement record, without changing the unconditional evolution of the measured system. We define an optimal measurement scheme as one which has the lowest average algorithmic information allowed. We also show how it is possible to extract information about system operator averages from the measurement records and their probabilities. The optimal measurement scheme, in the limit of weak coupling, determines the statistics of the variance of the measured variable directly. We discuss the relevance of such measurements for recent experiments in quantum optics.

Anisotropy-based performance analysis of linear discrete time invariant control systems

The anisotropic norm of a linear discrete-time-invariant system measures system output sensitivity to stationary Gaussian input disturbances of bounded mean anisotropy. Mean anisotropy characterizes the degree of predictability (or colouredness) and spatial non-roundness of the noise. The anisotropic norm falls between the H-2 and H-infinity norms and accommodates their loss of performance when the probability structure of input disturbances is not exactly known. This paper develops a method for numerical computation of the anisotropic norm which involves linked Riccati and Lyapunov equations and an associated special type equation.