Inactivation of L-type calcium channels is determined by the length of the N terminus of mutant beta(1) subunits.

Voltage-dependent calcium channel (Ca(v)) pores are modulated by cytosolic beta subunits. Four beta-subunit genes and their splice variants offer a wide structural array for tissue- or disease-specific biophysical gating phenotypes. For instance, the length of the N terminus of beta(2) subunits has major effects on activation and inactivation rates. We tested whether a similar mechanism principally operates in a beta(1) subunit. Wild-type beta(1a) subunit (N terminus length 60 aa) and its newly generated N-terminal deletion mutants (51, 27 and 18 aa) were examined within recombinant L-type calcium channel complexes (Ca(v)1.2 and alpha(2)delta2) in HEK293 cells at the whole-cell and single-channel level. Whole-cell currents were enhanced by co-transfection of the full-length beta(1a) subunit and by all truncated constructs. Voltage dependence of steady-state activation and inactivation did not depend on N terminus length, but inactivation rate was diminished by N terminus truncation. This was confirmed at the single-channel level, using ensemble average currents. Additionally, gating properties were estimated by Markov modeling. In confirmation of the descriptive analysis, inactivation rate, but none of the other transition rates, was reduced by shortening of the beta(1a) subunit N terminus. Our study shows that the length-dependent mechanism of modulating inactivation kinetics of beta(2) calcium channel subunits can be confirmed and extended to the beta(1) calcium channel subunit.

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completelyabsent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and byMartín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involvedparts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method isintroduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that thetheoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approachhas reasonable properties from a compositional point of view. In particular, it is “natural” in the sense thatit recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in thesame paper a substitution method for missing values on compositional data sets is introduced

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Parallel scheduling of multiclass M/M/m queues: Approximate and heavy-traffic optimization of achievable performance

We address the problem of scheduling a multiclass $M/M/m$ queue with Bernoulli feedback on $m$ parallel servers to minimize time-average linear holding costs. We analyze the performance of a heuristic priority-index rule, which extends Klimov's optimal solution to the single-server case: servers select preemptively customers with larger Klimov indices. We present closed-form suboptimality bounds (approximate optimality) for Klimov's rule, which imply that its suboptimality gap is uniformly bounded above with respect to (i) external arrival rates, as long as they stay within system capacity;and (ii) the number of servers. It follows that its relativesuboptimality gap vanishes in a heavy-traffic limit, as external arrival rates approach system capacity (heavy-traffic optimality). We obtain simpler expressions for the special no-feedback case, where the heuristic reduces to the classical $c \mu$ rule. Our analysis is based on comparing the expected cost of Klimov's ruleto the value of a strong linear programming (LP) relaxation of the system's region of achievable performance of mean queue lengths. In order to obtain this relaxation, we derive and exploit a new set ofwork decomposition laws for the parallel-server system. We further report on the results of a computational study on the quality of the $c \mu$ rule for parallel scheduling.

Location models for airline hubs behaving as M/D/c queues

Models are presented for the optimal location of hubs in airline networks, that take into consideration the congestion effects. Hubs, which are the most congested airports, are modeled as M/D/c queuing systems, that is, Poisson arrivals, deterministic service time, and {\em c} servers. A formula is derived for the probability of a number of customers in the system, which is later used to propose a probabilistic constraint. This constraint limits the probability of {\em b} airplanes in queue, to be lesser than a value $\alpha$. Due to the computational complexity of the formulation. The model is solved using a meta-heuristic based on tabu search. Computational experience is presented.

Intensity modulated radiation therapy or stereotactic fractionated radiotherapy for infratentorial ependymoma in children: a multicentric study.

This study was to evaluate the treatment dosimetry, efficacy and toxicity of intensity modulated radiation therapy (IMRT) and fractionated stereotactic radiotherapy (FSRT) in the management of infratentorial ependymoma. Between 1999 and 2007, seven children (median age, 3.1 years) with infratentorial ependymoma were planned with either IMRT (3 patients) or SFRT (4 patients), the latter after conventional posterior fossa irradiation. Two children underwent gross total resection. Median prescribed dose was 59.4 Gy (range, 55.8-60). The median follow-up for surviving patients was 4.8 years (range, 1.3-8). IMRT (median dose, 59.4 Gy) and FSRT (median dose, 55.8 Gy) achieved similar optimal target coverage. Percentages of maximum doses delivered to the cochleae (59.5 vs 85.0% Gy; P = 0.05) were significantly inferior with IMRT, when compared to FSRT planning. Percentages of maximum doses administered to the pituitary gland (38.2 vs 20.1%; P = 0.05) and optic chiasm (38.1 vs 14.1%; P = 0.001) were, however, significantly higher with IMRT, when compared to FSRT planning. No recurrences were observed at the last follow-up. The estimated 3-year progression-free survival and overall survival were 87.5 and 100%, respectively. No grade >1 acute toxicity was observed. Two patients presented late adverse events (grade 2 hypoacousia) during follow-up, without cognitive impairment. IMRT or FSRT for infratentorial ependymomas is effective and associated with a tolerable toxicity level. Both treatment techniques were able to capitalize their intrinsic conformal ability to deliver high-dose radiation. Larger series of patients treated with these two modalities will be necessary to more fully evaluate these delivery techniques.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Turn-on-time statistics of modulated lasers subjected to resonant weak optical feedback

We present analytical calculations of the turn-on-time probability distribution of intensity-modulated lasers under resonant weak optical feedback. Under resonant conditions, the external cavity round-trip time is taken to be equal to the modulation period. The probability distribution of the solitary laser results are modified to give reduced values of the mean turn-on-time and its variance. Numerical simulations have been carried out showing good agreement with the analytical results.

Front propagation in spatially modulated media

A two-dimensional reaction-diffusion front which propagates in a modulated medium is studied. The modulation consists of a spatial variation of the local front velocity in the transverse direction to that of the front propagation. We study analytically and numerically the final steady-state velocity and shape of the front, resulting from a nontrivial interplay between the local curvature effects and the global competition process between different maxima of the control parameter. The transient dynamics of the process is also studied numerically and analytically by means of singular perturbation techniques.

Wave competition in excitable modulated media.

The propagation of an initially planar front is studied within the framework of the photosensitive Belousov-Zhabotinsky reaction modulated by a smooth spatial variation of the local front velocity in the direction perpendicular to front propagation. Under this modulation, the wave front develops several fingers corresponding to the local maxima of the modulation function. After a transient, the wave front achieves a stationary shape that does not necessarily coincide with the one externally imposed by the modulation. Theoretical predictions for the selection criteria of fingers and steady-state velocity are experimentally validated.

Stochastic processes induced by dichotomous markov noise: Some exact dynamical results

Stochastic processes defined by a general Langevin equation of motion where the noise is the non-Gaussian dichotomous Markov noise are studied. A non-FokkerPlanck master differential equation is deduced for the probability density of these processes. Two different models are exactly solved. In the second one, a nonequilibrium bimodal distribution induced by the noise is observed for a critical value of its correlation time. Critical slowing down does not appear in this point but in another one.