Cyclosporine population pharmacokinetics in pediatric renal transplant recipients

Cyclosporine is an immunosuppressant drug with a narrow therapeutic index and large variability in pharmacokinetics. To improve cyclosporine dose individualization in children, we used population pharmacokinetic modeling to study the effects of developmental, clinical, and genetic factors on cyclosporine pharmacokinetics in altogether 176 subjects (age range: 0.36–20.2 years) before and up to 16 years after renal transplantation. Pre-transplantation test doses of cyclosporine were given intravenously (3 mg/kg) and orally (10 mg/kg), on separate occasions, followed by blood sampling for 24 hours (n=175). After transplantation, in a total of 137 patients, cyclosporine concentration was quantified at trough, two hours post-dose, or with dose-interval curves. One-hundred-four of the studied patients were genotyped for 17 putatively functionally significant sequence variations in the ABCB1, SLCO1B1, ABCC2, CYP3A4, CYP3A5, and NR1I2 genes. Pharmacokinetic modeling was performed with the nonlinear mixed effects modeling computer program, NONMEM. A 3-compartment population pharmacokinetic model with first order absorption without lag-time was used to describe the data. The most important covariate affecting systemic clearance and distribution volume was allometrically scaled body weight i.e. body weight**3/4 for clearance and absolute body weight for volume of distribution. The clearance adjusted by absolute body weight declined with age and pre-pubertal children (< 8 years) had an approximately 25% higher clearance/body weight (L/h/kg) than did older children. Adjustment of clearance for allometric body weight removed its relationship to age after the first year of life. This finding is consistent with a gradual reduction in relative liver size towards adult values, and a relatively constant CYP3A content in the liver from about 6–12 months of age to adulthood. The other significant covariates affecting cyclosporine clearance and volume of distribution were hematocrit, plasma cholesterol, and serum creatinine, explaining up to 20%–30% of inter-individual differences before transplantation. After transplantation, their predictive role was smaller, as the variations in hematocrit, plasma cholesterol, and serum creatinine were also smaller. Before transplantation, no clinical or demographic covariates were found to affect oral bioavailability, and no systematic age-related changes in oral bioavailability were observed. After transplantation, older children receiving cyclosporine twice daily as the gelatine capsule microemulsion formulation had an about 1.25–1.3 times higher bioavailability than did the younger children receiving the liquid microemulsion formulation thrice daily. Moreover, cyclosporine oral bioavailability increased over 1.5-fold in the first month after transplantation, returning thereafter gradually to its initial value in 1–1.5 years. The largest cyclosporine doses were administered in the first 3–6 months after transplantation, and thereafter the single doses of cyclosporine were often smaller than 3 mg/kg. Thus, the results suggest that cyclosporine displays dose-dependent, saturable pre-systemic metabolism even at low single doses, whereas complete saturation of CYP3A4 and MDR1 (P-glycoprotein) renders cyclosporine pharmacokinetics dose-linear at higher doses. No significant associations were found between genetic polymorphisms and cyclosporine pharmacokinetics before transplantation in the whole population for which genetic data was available (n=104). However, in children older than eight years (n=22), heterozygous and homozygous carriers of the ABCB1 c.2677T or c.1236T alleles had an about 1.3 times or 1.6 times higher oral bioavailability, respectively, than did non-carriers. After transplantation, none of the ABCB1 SNPs or any other SNPs were found to be associated with cyclosporine clearance or oral bioavailability in the whole population, in the patients older than eight years, or in the patients younger than eight years. In the whole population, in those patients carrying the NR1I2 g.-25385C–g.-24381A–g.-205_-200GAGAAG–g.7635G–g.8055C haplotype, however, the bioavailability of cyclosporine was about one tenth lower, per allele, than in non-carriers. This effect was significant also in a subgroup of patients older than eight years. Furthermore, in patients carrying the NR1I2 g.-25385C–g.-24381A–g.-205_-200GAGAAG–g.7635G–g.8055T haplotype, the bioavailability was almost one fifth higher, per allele, than in non-carriers. It may be possible to improve individualization of cyclosporine dosing in children by accounting for the effects of developmental factors (body weight, liver size), time after transplantation, and cyclosporine dosing frequency/formulation. Further studies are required on the predictive value of genotyping for individualization of cyclosporine dosing in children.

Outward spherical solidification of a superheated melt with time dependent boundary flux

Short-time analytical solutions of solid and liquid temperatures and freezing front have been obtained for the outward radially symmetric spherical solidification of a superheated melt. Although results are presented here only for time dependent boundary flux, the method of solution can be used for other kinds of boundary conditions also. Later, the analytical solution has been compared with the numerical solution obtained with the help of a finite difference numerical scheme in which the grid points change with the freezing front position. An efficient method of execution of the numerical scheme has been discussed in details. Graphs have been drawn for the total solidification times and temperature distributions in the solid.

Time Dependent Operator-split and Unsplit Schemes for One Dimensional Premixed Flames

This paper is concerned with a study of an operator split scheme and unsplit scheme for the computation of adiabatic freely propagating one-dimensional premixed flames. The study uses unsteady method for both split and unsplit schemes employing implicit chemistry and explicit diffusion, a combination which is stable and convergent. Solution scheme is not sensitive to the initial starting estimate and provides steady state even with straight line profiles (far from steady state) in small number of time steps. Two systems H2-Air and H2-NO (involving complex nitrogen chemistry) are considered in presentinvestigation. Careful comparison shows that the operator split approach is slightly superior than the unsplit when chemistry becomes complex. Comparison of computational times with those of existing steady and unsteady methods seems to suggest that the method employing implicit-explicit algorithm is very efficient and robust.

Relationship between microscopic and macroscopic orientational relaxation times in polar liquids

Microscopic relations between single-particle orientational relaxation time (T, ) , dielectric relaxation time ( T ~ )a,n d many-body orientational relaxation time ( T ~o)f a dipolar liquid are derived. We show that both T~ and T~ are influenced significantly by many-body effects. In the present theory, these many-body effects enter through the anisotropic part of the two-particle direct correlation function of the polar liquid. We use mean-spherical approximation (MSA) for dipolar hard spheres for explicit numerical evaluation of the relaxation times. We find that, although the dipolar correlation function is biexponential, the frequency-dependent dielectric constant is of simple Debye form, with T~ equal to the transverse polarization relaxation time. The microscopic T~ falls in between Debye and Onsager-Glarum expressions at large values of the static dielectric constant.

The variants of the protein toxins abrin and ricin A useful guide to understanding the processing events in the toxin transport

Kinetic data on inhibition of protein synthesis in thymocyte by three abrins and ricin have been obtained. The intrinsic efficiencies of A chains of four toxins to inactivate ribosomes, as analyzed by k1-versus-concentration plots were abrin II, III > ricin > abrin I. The lag times were 90, 66, 75 and 105 min at a 0.0744 nM concentration of each of abrin I, II, III and ricin, respectively. To account for the observed differences in the dose-dependent lag time, functional and structural variables of toxins such as binding efficiency of B chains to receptors and low-pH-induced structural alterations have been analyzed. The association constants obtained by stopped flow studies showed that abrin-I (4.13 × 105 M−1 s−1) association with putative receptor (4-methylumbelliferyl-α-D-galactoside) is nearly two times more often than abrin III (2.6 × 105 M−1 s−1) at 20°C. Equillibrium binding constants of abrin I and II to thymocyte at 37°C were 2.26 × 107 M−1 and 2.8 × 107 M−1 respectively. pH-induced structural alterations as studied by a parallel enhancement in 8-anilino-L-naphthalene sulfonate fluorescence revealed a high degree of qualitative similarity. These results taken with a nearly identical concentration-independent lag time (minimum lag of 41–42 min) indicated that the binding efficiencies and internalization efficiencies of these toxins are the same and that the observed difference in the dose-dependent lag time is causally related to the proposed processing event. The rates of reduction of inter-subunit disulfide bond, an obligatory step in the intoxication process, have been measured and compared under a variety of conditions. Intersubunit disulfide reduction of abrin I is fourfold faster than that of abrin II at pH 7.2. The rate of disulfide reduction in abrin I could be decreased 1 I-fold by adding lactose, compared to that without lactose. The observed differences in the efficiencies of A chains, the dose-dependent lag period, the modulating effect of lactose on the rates of disulfide reduction and similarity in binding properties make the variants a valuable tool to probe the processing events in toxin transport in detail.

Hydrogen-bond-directed self-assembly of D-(+)-dibenzoyltartaric acid and 4-aminopyridine: optical nonlinearities and stoichiometry-dependent novel structural features

Attempts to prepare hydrogen-bond-directed nonlinear optical materials from a 1:1 molar mixture Of D-(+)-dibenzoyltartaric acid (DBT, I) and 4-aminopyridine (4-AP, II) resulted in two salts of different stoichiometry. One of them crystallizes in an unusual 1.5:1 (acid:base) monohydrate salt form III while the other one crystallizes as 1:1 (acid:base) salt IV. Crystal structures of both of the salts were determined from single-crystal X-ray diffraction data. The salt III crystallizes in a monoclinic space group C2 with a = 30.339(l), b = 7.881(2), c = 14.355(1) angstrom, beta = 97.48(1)degrees, V = 3403.1(9) angstrom3, Z = 4, R(w) = 0.058, R(w)= 0.058. The salt IV also crystallizes in a monoclinic space group P2(1) with a = 7.500(1), b = 14.968(2), c = 10.370(1) angstrom, beta = 102.67(1)degrees, V = 1135.9(2) angstrom3, Z = 2, R = 0.043, R(w) = 0.043. Interestingly, two DBT molecules with distinctly different conformation are present in the same crystal lattice of salt III. Extensive hydrogen-bonding interactions are found in both of the salts, and both of them show SHG intensity 1.4-1.6 times that of urea.

Monte-Carlo estimation of time-dependent statistical characteristics of random dynamical systems

The problem of estimating the time-dependent statistical characteristics of a random dynamical system is studied under two different settings. In the first, the system dynamics is governed by a differential equation parameterized by a random parameter, while in the second, this is governed by a differential equation with an underlying parameter sequence characterized by a continuous time Markov chain. We propose, for the first time in the literature, stochastic approximation algorithms for estimating various time-dependent process characteristics of the system. In particular, we provide efficient estimators for quantities such as the mean, variance and distribution of the process at any given time as well as the joint distribution and the autocorrelation coefficient at different times. A novel aspect of our approach is that we assume that information on the parameter model (i.e., its distribution in the first case and transition probabilities of the Markov chain in the second) is not available in either case. This is unlike most other work in the literature that assumes availability of such information. Also, most of the prior work in the literature is geared towards analyzing the steady-state system behavior of the random dynamical system while our focus is on analyzing the time-dependent statistical characteristics which are in general difficult to obtain. We prove the almost sure convergence of our stochastic approximation scheme in each case to the true value of the quantity being estimated. We provide a general class of strongly consistent estimators for the aforementioned statistical quantities with regular sample average estimators being a specific instance of these. We also present an application of the proposed scheme on a widely used model in population biology. Numerical experiments in this framework show that the time-dependent process characteristics as obtained using our algorithm in each case exhibit excellent agreement with exact results. (C) 2010 Elsevier Inc. All rights reserved.

Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Monitoring of ultraviolet pulse rate dependent photomechanical actuation in carbon nanotubes using fiber Bragg gratings

In this Letter, we present a non-contact method of controlling and monitoring photomechanical actuation in carbon nanotubes (CNT) by exposing it to ultra-violet radiation at different pulse rates (10 to 200 Hz). This is accomplished by imparting a reversible photo induced strain (5-330 mu epsilon) on CNT coated fibre Bragg gratings; CNT undergoes an internal reversible structural change due to cyclic photon absorption that leads to the development of mechanical strain, which in turn allows reversible switching of the Bragg wavelength. The results also reveal an interesting pulse rate dependent rise and fall times of photomechanical actuation in CNT. (C) 2014 AIP Publishing LLC.

Using long time series of agricultural-derived nitrates for estimating catchment transit times

The estimation of water and solute transit times in catchments is crucial for predicting the response of hydrosystems to external forcings (climatic or anthropogenic). The hydrogeochemical signatures of tracers (either natural or anthropogenic) in streams have been widely used to estimate transit times in catchments as they integrate the various processes at stake. However, most of these tracers are well suited for catchments with mean transit times lower than about 4-5 years. Since the second half of the 20th century, the intensification of agriculture led to a general increase of the nitrogen load in rivers. As nitrate is mainly transported by groundwater in agricultural catchments, this signal can be used to estimate transit times greater than several years, even if nitrate is not a conservative tracer. Conceptual hydrological models can be used to estimate catchment transit times provided their consistency is demonstrated, based on their ability to simulate the stream chemical signatures at various time scales and catchment internal processes such as N storage in groundwater. The objective of this study was to assess if a conceptual lumped model was able to simulate the observed patterns of nitrogen concentration, at various time scales, from seasonal to pluriannual and thus if it was relevant to estimate the nitrogen transit times in headwater catchments. A conceptual lumped model, representing shallow groundwater flow as two parallel linear stores with double porosity, and riparian processes by a constant nitrogen removal function, was applied on two paired agricultural catchments which belong to the Research Observatory ORE AgrHys. The Global Likelihood Uncertainty Estimation (GLUE) approach was used to estimate parameter values and uncertainties. The model performance was assessed on (i) its ability to simulate the contrasted patterns of stream flow and stream nitrate concentrations at seasonal and inter-annual time scales, (ii) its ability to simulate the patterns observed in groundwater at the same temporal scales, and (iii) the consistency of long-term simulations using the calibrated model and the general pattern of the nitrate concentration increase in the region since the beginning of the intensification of agriculture in the 1960s. The simulated nitrate transit times were found more sensitive to climate variability than to parameter uncertainty, and average values were found to be consistent with results from others studies in the same region involving modeling and groundwater dating. This study shows that a simple model can be used to simulate the main dynamics of nitrogen in an intensively polluted catchment and then be used to estimate the transit times of these pollutants in the system which is crucial to guide mitigation plans design and assessment. (C) 2015 Elsevier B.V. All rights reserved.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

A Study of Pressure-Dependent Squeeze Film Stiffness as a Resonance Modulator Using Static and Dynamic Measurements

We report on the resonant frequency modulation of inertial microelectromechanical systems (MEMS) structures due to squeeze film stiffness over a range of working pressures. Squeeze film effects have been studied extensively, but mostly in the context of damping and Q-factor determination of dynamic MEMS structures, typically suspended over a fixed substrate with a very thin air gap. Here, we show with experimental measurements and analytical calculations how the pressure-dependent air springs (squeeze film stiffness) change the resonant frequency of an inertial MEMS structure by as much as five times. For capturing the isolated effect of the squeeze film stiffness, we first determine the static stiffness of our structure with atomic force microscope probing and then study the effect of the air spring by measuring the dynamic response of the structure, thus finding the resonant frequencies while varying the air pressure from 1 to 905 mbar. We also verify our results by analytical and Finite Element Method calculations. Our findings show that the pressure-dependent squeeze film stiffness can affect a rather huge range of frequency modulation (>400%) and, therefore, can be used as a design parameter for exploiting this effect in MEMS devices. 2014-0310]

The role of intra- and inter-site exchange correlations in the extended Falicov-Kimball model on a triangular lattice

Ground state magnetic properties of the spin-dependent Falicov-Kimball model (FKM) are studied by incorporating the intrasite exchange correlation J (between itinerant d- and localized f-electrons) and intersite (superexchange) correlation J (between localized f-electrons) on a triangular lattice for two different fillings. Numerical diagonalization and Monte-Carlo techniques are used to determine the ground state magnetic properties. Transitions from antiferromagnetic to ferromagnetic and again to re-entrant antiferromagnetic phase is observed in a wide range of parameter space. The magnetic moments of d- and f-electrons are observed to depend strongly on the value off, J and also on the total number of d-electrons (N-d). (C) 2015 Elsevier Ltd. All rights reserved.

Toughness Of Ni/Al2O3 Interfaces As Dependent On Micron-Scale Plasticity And Atomistic-Scale Separation

Ceramic/metal interfaces were studied that fail by atomistic separation accompanied by plastic dissipation in the metal. The macroscopic toughness of the specific Ni alloy/Al2O3 interface considered is typically on the order of ten times the atomistic work of separation in mode I and even higher if combinations of mode I and mode II act on the interface. Inputs to the computational model of interface toughness are: (i) strain gradient plasticity applied to the Ni alloy with a length parameter determined by an indentation test, and (ii) a potential characterizing mixed mode separation of the interface fit to atomistic results. The roles of the several length parameters in the strain gradient plasticity are determined for indentation and crack growth. One of the parameters is shown to be of dominant importance, thus establishing that indentation can be used to measure the relevant length parameter. Recent results for separation of Ni/Al2O3 interfaces computed by atomistic methods are reviewed, including a set of results computed for mixed mode separation. An approximate potential fit to these results is characterized by the work of separation, the peak separation stress for normal separation and the traction-displacement relation in pure shearing of the interface. With these inputs, the model for steady-state crack growth is used to compute the toughness of the interface under mode I and under the full range of mode mix. The effect of interface strength and the work of separation on macroscopic toughness is computed. Fundamental implications for plasticity-enhanced toughness emerge.

CMB as a probe of new physics and old times

Cosmic birefringence (CB)---a rotation of photon-polarization plane in vacuum---is a generic signature of new scalar fields that could provide dark energy. Previously, WMAP observations excluded a uniform CB-rotation angle larger than a degree.

In this thesis, we develop a minimum-variance--estimator formalism for reconstructing direction-dependent rotation from full-sky CMB maps, and forecast more than an order-of-magnitude improvement in sensitivity with incoming Planck data and future satellite missions. Next, we perform the first analysis of WMAP-7 data to look for rotation-angle anisotropies and report null detection of the rotation-angle power-spectrum multipoles below L=512, constraining quadrupole amplitude of a scale-invariant power to less than one degree. We further explore the use of a cross-correlation between CMB temperature and the rotation for detecting the CB signal, for different quintessence models. We find that it may improve sensitivity in case of marginal detection, and provide an empirical handle for distinguishing details of new physics indicated by CB.

We then consider other parity-violating physics beyond standard models---in particular, a chiral inflationary-gravitational-wave background. We show that WMAP has no constraining power, while a cosmic-variance--limited experiment would be capable of detecting only a large parity violation. In case of a strong detection of EB/TB correlations, CB can be readily distinguished from chiral gravity waves.

We next adopt our CB analysis to investigate patchy screening of the CMB, driven by inhomogeneities during the Epoch of Reionization (EoR). We constrain a toy model of reionization with WMAP-7 data, and show that data from Planck should start approaching interesting portions of the EoR parameter space and can be used to exclude reionization tomographies with large ionized bubbles.

In light of the upcoming data from low-frequency radio observations of the redshifted 21-cm line from the EoR, we examine probability-distribution functions (PDFs) and difference PDFs of the simulated 21-cm brightness temperature, and discuss the information that can be recovered using these statistics. We find that PDFs are insensitive to details of small-scale physics, but highly sensitive to the properties of the ionizing sources and the size of ionized bubbles.

Finally, we discuss prospects for related future investigations.