Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

RecNcMIC3-1-R is a microneme- and rhoptry-based chimeric antigen that protects against acute neosporosis and limits cerebral parasite load in the mouse model for Neospora caninum infection

In order to achieve host cell entry, the apicomplexan parasite Neospora caninum relies on the contents of distinct organelles, named micronemes, rhoptries and dense granules, which are secreted at defined timepoints during and after host cell entry. It was shown previously that a vaccine composed of a mixture of three recombinant antigens, corresponding to the two microneme antigens NcMIC1 and NcMIC3 and the rhoptry protein NcROP2, prevented disease and limited cerebral infection and transplacental transmission in mice. In this study, we selected predicted immunogenic domains of each of these proteins and created four different chimeric antigens, with the respective domains incorporated into these chimers in different orders. Following vaccination, mice were challenged intraperitoneally with 2 × 10(6)N. caninum tachzyoites and were then carefully monitored for clinical symptoms during 4 weeks post-infection. Of the four chimeric antigens, only recNcMIC3-1-R provided complete protection against disease with 100% survivors, compared to 40-80% of survivors in the other groups. Serology did not show any clear differences in total IgG, IgG1 and IgG2a levels between the different treatment groups. Vaccination with all four chimeric variants generated an IL-4 biased cytokine expression, which then shifted to an IFN-γ-dominated response following experimental infection. Sera of recNcMIC3-1-R vaccinated mice reacted with each individual recombinant antigen, as well as with three distinct bands in Neospora extracts with similar Mr as NcMIC1, NcMIC3 and NcROP2, and exhibited distinct apical labeling in tachyzoites. These results suggest that recNcMIC3-1-R is an interesting chimeric vaccine candidate and should be followed up in subsequent studies in a fetal infection model.

Posterior Simulation in the Generalized Linear Model with Semiparmetric Random Effects

Generalized linear mixed models with semiparametric random effects are useful in a wide variety of Bayesian applications. When the random effects arise from a mixture of Dirichlet process (MDP) model, normal base measures and Gibbs sampling procedures based on the Pólya urn scheme are often used to simulate posterior draws. These algorithms are applicable in the conjugate case when (for a normal base measure) the likelihood is normal. In the non-conjugate case, the algorithms proposed by MacEachern and Müller (1998) and Neal (2000) are often applied to generate posterior samples. Some common problems associated with simulation algorithms for non-conjugate MDP models include convergence and mixing difficulties. This paper proposes an algorithm based on the Pólya urn scheme that extends the Gibbs sampling algorithms to non-conjugate models with normal base measures and exponential family likelihoods. The algorithm proceeds by making Laplace approximations to the likelihood function, thereby reducing the procedure to that of conjugate normal MDP models. To ensure the validity of the stationary distribution in the non-conjugate case, the proposals are accepted or rejected by a Metropolis-Hastings step. In the special case where the data are normally distributed, the algorithm is identical to the Gibbs sampler.

Mixture Cure Survival Models with Dependent Censoring

A number of authors have studies the mixture survival model to analyze survival data with nonnegligible cure fractions. A key assumption made by these authors is the independence between the survival time and the censoring time. To our knowledge, no one has studies the mixture cure model in the presence of dependent censoring. To account for such dependence, we propose a more general cure model which allows for dependent censoring. In particular, we derive the cure models from the perspective of competing risks and model the dependence between the censoring time and the survival time using a class of Archimedean copula models. Within this framework, we consider the parameter estimation, the cure detection, and the two-sample comparison of latency distribution in the presence of dependent censoring when a proportion of patients is deemed cured. Large sample results using the martingale theory are obtained. We applied the proposed methodologies to the SEER prostate cancer data.

COMBUSTION INITIATION BY ELECTRICAL-DISCHARGE-INDUCED PLASMA IN LEAN AND DILUTE METHANE-AIR MIXTURE: EXPERIMENTAL AND MODELING INVESTIGATION

This dissertation represents experimental and numerical investigations of combustion initiation trigged by electrical-discharge-induced plasma within lean and dilute methane air mixture. This research topic is of interest due to its potential to further promote the understanding and prediction of spark ignition quality in high efficiency gasoline engines, which operate with lean and dilute fuel-air mixture. It is specified in this dissertation that the plasma to flame transition is the key process during the spark ignition event, yet it is also the most complicated and least understood procedure. Therefore the investigation is focused on the overlapped periods when plasma and flame both exists in the system. Experimental study is divided into two parts. Experiments in Part I focuses on the flame kernel resulting from the electrical discharge. A number of external factors are found to affect the growth of the flame kernel, resulting in complex correlations between discharge and flame kernel. Heat loss from the flame kernel to code ambient is found to be a dominant factor that quenches the flame kernel. Another experimental focus is on the plasma channel. Electrical discharges into gases induce intense and highly transient plasma. Detailed observation of the size and contents of the discharge-induced plasma channel is performed. Given the complex correlation and the multi-discipline physical/chemical processes involved in the plasma-flame transition, the modeling principle is taken to reproduce detailed transitions numerically with minimum analytical assumptions. Detailed measurement obtained from experimental work facilitates the more accurate description of initial reaction conditions. The novel and unique spark source considering both energy and species deposition is defined in a justified manner, which is the key feature of this Ignition by Plasma (IBP) model. The results of numerical simulation are intuitive and the potential of numerical simulation to better resolve the complex spark ignition mechanism is presented. Meanwhile, imperfections of the IBP model and numerical simulation have been specified and will address future attentions.

Bayesian change-point analysis in linear regression model with scale mixtures of normal distributions

In this thesis, we consider Bayesian inference on the detection of variance change-point models with scale mixtures of normal (for short SMN) distributions. This class of distributions is symmetric and thick-tailed and includes as special cases: Gaussian, Student-t, contaminated normal, and slash distributions. The proposed models provide greater flexibility to analyze a lot of practical data, which often show heavy-tail and may not satisfy the normal assumption. As to the Bayesian analysis, we specify some prior distributions for the unknown parameters in the variance change-point models with the SMN distributions. Due to the complexity of the joint posterior distribution, we propose an efficient Gibbs-type with Metropolis- Hastings sampling algorithm for posterior Bayesian inference. Thereafter, following the idea of [1], we consider the problems of the single and multiple change-point detections. The performance of the proposed procedures is illustrated and analyzed by simulation studies. A real application to the closing price data of U.S. stock market has been analyzed for illustrative purposes.

On trend estimation under monotone Gaussian subordination with long-memory: application to fossil pollen series

Fossil pollen data from stratigraphic cores are irregularly spaced in time due to non-linear age-depth relations. Moreover, their marginal distributions may vary over time. We address these features in a nonparametric regression model with errors that are monotone transformations of a latent continuous-time Gaussian process Z(T). Although Z(T) is unobserved, due to monotonicity, under suitable regularity conditions, it can be recovered facilitating further computations such as estimation of the long-memory parameter and the Hermite coefficients. The estimation of Z(T) itself involves estimation of the marginal distribution function of the regression errors. These issues are considered in proposing a plug-in algorithm for optimal bandwidth selection and construction of confidence bands for the trend function. Some high-resolution time series of pollen records from Lago di Origlio in Switzerland, which go back ca. 20,000 years are used to illustrate the methods.

Estimating and quantifying uncertainties on level sets using the Vorob'ev expectation and deviation with Gaussian process models

Several methods based on Kriging have recently been proposed for calculating a probability of failure involving costly-to-evaluate functions. A closely related problem is to estimate the set of inputs leading to a response exceeding a given threshold. Now, estimating such a level set—and not solely its volume—and quantifying uncertainties on it are not straightforward. Here we use notions from random set theory to obtain an estimate of the level set, together with a quantification of estimation uncertainty. We give explicit formulae in the Gaussian process set-up and provide a consistency result. We then illustrate how space-filling versus adaptive design strategies may sequentially reduce level set estimation uncertainty.

Statistical inference for the optimal approximating model

In the setting of high-dimensional linear models with Gaussian noise, we investigate the possibility of confidence statements connected to model selection. Although there exist numerous procedures for adaptive (point) estimation, the construction of adaptive confidence regions is severely limited (cf. Li in Ann Stat 17:1001–1008, 1989). The present paper sheds new light on this gap. We develop exact and adaptive confidence regions for the best approximating model in terms of risk. One of our constructions is based on a multiscale procedure and a particular coupling argument. Utilizing exponential inequalities for noncentral χ2-distributions, we show that the risk and quadratic loss of all models within our confidence region are uniformly bounded by the minimal risk times a factor close to one.

Transcutaneous treatment with vetdrop(®) sustains the adjacent cartilage in a microfracturing joint defect model in sheep

The significance of the adjacent cartilage in cartilage defect healing is not yet completely understood. Furthermore, it is unknown if the adjacent cartilage can somehow be influenced into responding after cartilage damage. The present study was undertaken to investigate whether the adjacent cartilage can be better sustained after microfracturing in a cartilage defect model in the stifle joint of sheep using a transcutaneous treatment concept (Vetdrop(®)). Carprofen and chito-oligosaccharids were added either as single components or as a mixture to a vehicle suspension consisting of a herbal carrier oil in a water-in-oil phase. This mixture was administered onto the skin with the aid of a specific applicator during 6 weeks in 28 sheep, allocated into 6 different groups, that underwent microfracturing surgery either on the left or the right medial femoral condyle. Two groups served as control and were either treated intravenously or sham treated with oxygen only. Sheep were sacrificed and their medial condyle histologically evaluated qualitatively and semi-quantitatively according to 4 different scoring systems (Mankin, ICRS, Little and O'Driscoll). The adjacent cartilage of animals of group 4 treated transcutaneously with vehicle, chito-oligosaccharids and carprofen had better histological scores compared to all the other groups (Mankin 3.3±0.8, ICRS 15.7±0.7, Little 9.0±1.4). Complete defect filling was absent from the transcutaneous treatment groups. The experiment suggests that the adjacent cartilage is susceptible to treatment and that the combination of vehicle, chitooligosaccharids and carprofen may sustain the adjacent cartilage during the recovery period.

Biodiversity effects on nitrate concentrations in soil solution: a Bayesian model

Ecosystems are faced with high rates of species loss which has consequences for their functions and services. To assess the effects of plant species diversity on the nitrogen (N) cycle, we developed a model for monthly mean nitrate (NO3-N) concentrations in soil solution in 0-30 cm mineral soil depth using plant species and functional group richness and functional composition as drivers and assessing the effects of conversion of arable land to grassland, spatially heterogeneous soil properties, and climate. We used monthly mean NO3-N concentrations from 62 plots of a grassland plant diversity experiment from 2003 to 2006. Plant species richness (1-60) and functional group composition (1-4 functional groups: legumes, grasses, non-leguminous tall herbs, non-leguminous small herbs) were manipulated in a factorial design. Plant community composition, time since conversion from arable land to grassland, soil texture, and climate data (precipitation, soil moisture, air and soil temperature) were used to develop one general Bayesian multiple regression model for the 62 plots to allow an in-depth evaluation using the experimental design. The model simulated NO3-N concentrations with an overall Bayesian coefficient of determination of 0.48. The temporal course of NO3-N concentrations was simulated differently well for the individual plots with a maximum plot-specific Nash-Sutcliffe Efficiency of 0.57. The model shows that NO3-N concentrations decrease with species richness, but this relation reverses if more than approx. 25 % of legume species are included in the mixture. Presence of legumes increases and presence of grasses decreases NO3-N concentrations compared to mixtures containing only small and tall herbs. Altogether, our model shows that there is a strong influence of plant community composition on NO3-N concentrations.

Do different anesthesia regimes affect hippocampal apoptosis and neurologic deficits in a rodent cardiac arrest model?

Background Different anesthesia regimes are commonly used in experimental models of cardiac arrest, but the effects of various anesthetics on clinical outcome parameters are unknown. We conducted a study in which we subjected rats to cardiac arrest under medetomidine/ketamine or sevoflurane/fentanyl anesthesia. Methods Asystolic cardiac arrest for 8 minutes was induced in 73 rats with a mixture of potassium chloride and esmolol. Daily behavioral and neurological examination included the open field test (OFT), the tape removal test (TRT) and a neurodeficit score (NDS). Animals were randomized for sacrifice on day 2 or day 5 and brains were harvested for histology in the hippocampus cornus ammonis segment CA1. The inflammatory markers IL-6, TNF-α, MCP-1 and MIP-1α were assessed in cerebrospinal fluid (CSF). Proportions of survival were tested with the Fisher’s exact test, repeated measurements were assessed with the Friedman’s test; the baseline values were tested using Mann–Whitney U test and the difference of results of repeated measures were compared. Results In 31 animals that survived beyond 24 hours neither OFT, TRT nor NDS differed between the groups; histology was similar on day 2. On day 5, significantly more apoptosis in the CA1 segment of the hippocampus was found in the sevoflurane/fentanyl group. MCP-1 was higher on day 5 in the sevoflurane/fentanyl group (p = 0.04). All other cyto- and chemokines were below detection threshold. Conclusion In our cardiac arrest model neurological function was not influenced by different anesthetic regimes; in contrast, anesthesia with sevoflurane/fentanyl results in increased CSF inflammation and histologic damage at day 5 post cardiac arrest.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (1 age) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.

An improved simple rat model for global cerebral ischaemia by induced cardiac arrest

OBJECTIVES Cerebral hypoxic-ischaemic injury following cardiac arrest is a devastating disease affecting thousands of patients each year. There is a complex interaction between post-resuscitation injury after whole-body ischaemia-reperfusion and cerebral damage which cannot be explored in in vitro systems only; there is a need for animal models. In this study, we describe and evaluate the feasibility and efficiency of our simple rodent cardiac arrest model. METHODS Ten wistar rats were subjected to 9 and 10 minutes of cardiac arrest. Cardiac arrest was introduced with a mixture of the short-acting beta-blocking drug esmolol and potassium chloride. RESULTS All animals could be resuscitated within 1 minute, and survived until day 5.General health score and neurobehavioural testing indicated substantial impairment after cardiac arrest, without differences between groups. Histological examination of the hippocampus CA1 segment, the most vulnerable segment of the cerebrum, demonstrated extensive damage in the cresyl violet staining, as well as in the Fluoro-Jade B staining and in the Iba-1 staining, indicating recruitment of microglia after the hypoxic-ischaemic event. Again, there were no differences between the 9- and 10-minute cardiac arrest groups. DISCUSSION We were able to establish a simple and reproducible 9- and 10-minute rodent cardiac arrest models with a well-defined no-flow-time. Extensive damage can be found in the hippocampus CA1 segment. The lack of difference between 9- and 10-minute cardiac arrest time in the neuropsychological, the open field test and the histological evaluations is mainly due to the small sample size.