Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs.

In population pharmacokinetic studies, the precision of parameter estimates is dependent on the population design. Methods based on the Fisher information matrix have been developed and extended to population studies to evaluate and optimize designs. In this paper we propose simple programming tools to evaluate population pharmacokinetic designs. This involved the development of an expression for the Fisher information matrix for nonlinear mixed-effects models, including estimation of the variance of the residual error. We implemented this expression as a generic function for two software applications: S-PLUS and MATLAB. The evaluation of population designs based on two pharmacokinetic examples from the literature is shown to illustrate the efficiency and the simplicity of this theoretic approach. Although no optimization method of the design is provided, these functions can be used to select and compare population designs among a large set of possible designs, avoiding a lot of simulations.

Relationship of glycosaminoglycan and matrix changes to vascular smooth cell phenotype modulation in rabbit arteries after acute injury

Purpose: The phenotype of vascular smooth muscle cells (SMCs) is altered in several arterial pathologies, including the neointima formed after acute arterial injury. This study examined the time course of this phenotypic change in relation to changes in the amount and distribution of matrix glycosaminoglycans. Methods: The immunochemical staining of heparan sulphates (HS) and chondroitin sulphates (CS) in the extracellular matrix of the arterial wall was examined at early points after balloon catheter injury of the rabbit carotid artery. SMC phenotype was assessed by means of ultrastructural morphometry of the cytoplasmic volume fraction of myofilaments. The proportions of cell and matrix components in the media were analyzed with similar morphometric techniques. Results: HS and CS were shown in close association with SMCs of the uninjured arterial media as well as being more widespread within the matrix. Within 6 hours after arterial injury, there was loss of the regular pericellular distribution of both HS and CS, which was associated with a significant expansion in the extracellular space. This preceded the change in ultrastructural phenotype of the SMCs. The glycosaminoglycan loss was most exaggerated at 4 days, after which time the HS and CS reappeared around the medial SMCs. SMCs of the recovering media were able to rapidly replace their glycosaminoglycans, whereas SMCs of the developing neointima failed to produce HS as readily as they produced CS. Conclusions: These studies indicate that changes in glycosaminoglycans of the extracellular matrix precede changes in SMC phenotype after acute arterial injury. In the recovering arterial media, SMCs replace their matrix glycosaminoglycans rapidly, whereas the newly established neointima fails to produce similar amounts of heparan sulphates.

Subspace wavepacket evolution with Newton polynomials

Time-dependent wavepacket evolution techniques demand the action of the propagator, exp(-iHt/(h)over-bar), on a suitable initial wavepacket. When a complex absorbing potential is added to the Hamiltonian for combating unwanted reflection effects, polynomial expansions of the propagator are selected on their ability to cope with non-Hermiticity. An efficient subspace implementation of the Newton polynomial expansion scheme that requires fewer dense matrix-vector multiplications than its grid-based counterpart has been devised. Performance improvements are illustrated with some benchmark one and two-dimensional examples. (C) 2001 Elsevier Science B.V. All rights reserved.

Complex natural resonances of conducting planar objects buried in a dielectric half-space

A scheme is presented to incorporate a mixed potential integral equation (MPIE) using Michalski's formulation C with the method of moments (MoM) for analyzing the scattering of a plane wave from conducting planar objects buried in a dielectric half-space. The robust complex image method with a two-level approximation is used for the calculation of the Green's functions for the half-space. To further speed up the computation, an interpolation technique for filling the matrix is employed. While the induced current distributions on the object's surface are obtained in the frequency domain, the corresponding time domain responses are calculated via the inverse fast Fourier transform (FFT), The complex natural resonances of targets are then extracted from the late time response using the generalized pencil-of-function (GPOF) method. We investigate the pole trajectories as we vary the distance between strips and the depth and orientation of single, buried strips, The variation from the pole position of a single strip in a homogeneous dielectric medium was only a few percent for most of these parameter variations.

Freeform fabrication of aluminum metal-matrix composites

A series of metal-matrix composites were formed by extrusion freeform, fabrication of a sinterable aluminum alloy in combination with silicon carbide particles and whiskers, carbon fibers, alumina particles, and hollow flyash cenospheres. Silicon carbide particles were most successful in that the composites retained high density with up to 20 vol% of reinforcement and the strength approximately doubles over the strength of the metal matrix alone. Comparison with simple models suggests that this unexpectedly high degree of reinforcement can be attributed to the concentration of small silicon carbide particles around the larger metal powder. This fabrication method also allows composites to be formed with hollow spheres that cannot be formed by other powder or melt methods.

Effects of moist chilling and solid matrix priming on germination of loblolly pine (Pinus taeda L.) seeds

Loblolly pine ( Pinus taeda L.) seeds from sources with a mild climate under maritime influence (North Carolina) required shorter moist chilling to achieve maximum germination vigor than seeds from sources with a harsher continental climate (Oklahoma). Solid matrix priming (SMP) for 6 d achieved as much as 60 d of moist chilling to improve rapidity, synchrony and completeness of germination for three of the four families studied. SMP after moist chilling increased the rapidity, synchrony and completeness of germination. The benefit of SMP was greatest for non-stratified seeds and the benefit decreased with length of moist chilling. In general, delaying planting for one week after SMP had minor effects on germination when seeds were kept in the SMP matrix at 4 degreesC. Delayed planting after SMP can increase germination rapidity and synchrony of seeds that have received long moist chilling and reduce the benefit of SMP in non-moist-chilled seeds.

On the importance function in splitting simulation

The splitting method is a simulation technique for the estimation of very small probabilities. In this technique, the sample paths are split into multiple copies, at various stages in the simulation. Of vital importance to the efficiency of the method is the Importance Function (IF). This function governs the placement of the thresholds or surfaces at which the paths are split. We derive a characterisation of the optimal IF and show that for multi-dimensional models the natural choice for the IF is usually not optimal. We also show how nearly optimal splitting surfaces can be derived or simulated using reverse time analysis. Our numerical experiments illustrate that by using the optimal IF, one can obtain a significant improvement in simulation efficiency.

A new joint coding and modulation scheme for channels with clipping

A major limitation in any high-performance digital communication system is the linearity region of the transmitting amplifier. Nonlinearities typically lead to signal clipping. Efficient communication in such conditions requires maintaining a low peak-to-average power ratio (PAR) in the transmitted signal while achieving a high throughput of data. Excessive PAR leads either to frequent clipping or to inadequate resolution in the analog-to-digital or digital-to-analog converters. Currently proposed signaling schemes for future generation wireless communications suffer from a high PAR. This paper presents a new signaling scheme for channels with clipping which achieves a PAR as low as 3. For a given linear range in the transmitter's digital-to-analog converter, this scheme achieves a lower bit-error rate than existing multicarrier schemes, owing to increased separation between constellation points. We present the theoretical basis for this new scheme, approximations for the expected bit-error rate, and simulation results. (C) 2002 Elsevier Science (USA).

Practical scheme for quantum computation with any two-qubit entangling gate

Which gates are universal for quantum computation? Although it is well known that certain gates on two-level quantum systems (qubits), such as the controlled-NOT, are universal when assisted by arbitrary one-qubit gates, it has only recently become clear precisely what class of two-qubit gates is universal in this sense. We present an elementary proof that any entangling two-qubit gate is universal for quantum computation, when assisted by one-qubit gates. A proof of this result for systems of arbitrary finite dimension has been provided by Brylinski and Brylinski; however, their proof relies on a long argument using advanced mathematics. In contrast, our proof provides a simple constructive procedure which is close to optimal and experimentally practical.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

Thermolysin activates equine lamellar hoof matrix metalloproteinases

Cultured equine lamellar hoof explants secrete the pro-enzymes matrix metalloproteinse-2 (MMP-2, 72 kDa) and MMP-2 (92 kDa). Untreated explants remained intact tested on a calibrated force transducer, but when treated with an NIMP activator, developed in-vitro laminitis, separating at the dermal-epidermal junction. Explants treated with the bacterial protease thermolysin separated dose-dependently; this was accompanied by activation of both MMP-2 and -9. Thermolysin-mediated NIP activation did not occur in a cell-free system and was not inhibited by the addition of the MMP inhibitor and batimastat. These findings suggest that thermolysin-mediated gelatinase activation is not dependent on membrane-bound matrix metalloproteinase (MT-MMP) activation, providing further evidence that bacteria can produce potent MMP activators that probably facilitate host invasion. (C) 2002 Harcourt Publishers Ltd.