A mesoscopic simulation method for electrolyte solutions

This thesis deals with the development of a novel simulation technique for macromolecules in electrolyte solutions, with the aim of a performance improvement over current molecular-dynamics based simulation methods. In solutions containing charged macromolecules and salt ions, it is the complex interplay of electrostatic interactions and hydrodynamics that determines the equilibrium and non-equilibrium behavior. However, the treatment of the solvent and dissolved ions makes up the major part of the computational effort. Thus an efficient modeling of both components is essential for the performance of a method. With the novel method we approach the solvent in a coarse-grained fashion and replace the explicit-ion description by a dynamic mean-field treatment. Hence we combine particle- and field-based descriptions in a hybrid method and thereby effectively solve the electrokinetic equations. The developed algorithm is tested extensively in terms of accuracy and performance, and suitable parameter sets are determined. As a first application we study charged polymer solutions (polyelectrolytes) in shear flow with focus on their viscoelastic properties. Here we also include semidilute solutions, which are computationally demanding. Secondly we study the electro-osmotic flow on superhydrophobic surfaces, where we perform a detailed comparison to theoretical predictions.

Dynamic characterisation of four nine-story large-panel R.C. buildings in Bishkek (Kyrgyzstan): A comparison between experimental ambient vibration analysis and numerical finite element modeling

With the outlook of improving seismic vulnerability assessment for the city of Bishkek (Kyrgyzstan), the global dynamic behaviour of four nine-storey r.c. large-panel buildings in elastic regime is studied. The four buildings were built during the Soviet era within a serial production system. Since they all belong to the same series, they have very similar geometries both in plan and in height. Firstly, ambient vibration measurements are performed in the four buildings. The data analysis composed of discrete Fourier transform, modal analysis (frequency domain decomposition) and deconvolution interferometry, yields the modal characteristics and an estimate of the linear impulse response function for the structures of the four buildings. Then, finite element models are set up for all four buildings and the results of the numerical modal analysis are compared with the experimental ones. The numerical models are finally calibrated considering the first three global modes and their results match the experimental ones with an error of less then 20%.

Development of an LC-MS/MS method for the quantitation of 55 compounds prescribed in combined cardiovascular therapy

This paper reports an LC-MS/MS method with positive electrospray ionization for the screening of commonly prescribed cardiovascular drugs in human plasma, including compounds with antihypertensive (57), antidiabetic (12), hypolipemiant (5), anticoagulant (2) and platelet anti-aggregation (2) effects. Sample treatment consisted of a simple protein precipitation with MeOH/0.1 M ZnSO₄ (4:1, v/v) solution after the addition of internal standard, followed by evaporation and reconstitution. Analytes separation was performed on a Polar-RP column (150 m x 2 mm, 4 μm) using a gradient elution of 15 min. The MS system was operated in MRM mode, monitoring one quantitation and one confirmation transition for each analyte. The recovery of the protein precipitation step ranged from 50 to 70% for most of the compounds, while some were considerably affected by matrix effects. Since several analytes fulfilled the linearity, accuracy and precision values required by the ICH guidelines, the method proved to be suitable for their quantitative analysis. The limits of quantitation varied from 0.38 to 9.1 μg/L and the limits of detection from 0.12 to 5.34 μg/L. The method showed to be suitable for the detection of plasma samples of patients under cardiovascular treatment with the studied drugs, and for 55 compounds reliable quantitative results could be obtained.

A Near-Infrared Spectroscopic Study of Young Field Ultracool Dwarfs

We present a near-infrared (0.9-2.4 mu m) spectroscopic study of 73 field ultracool dwarfs having spectroscopic and/or kinematic evidence of youth (approximate to 10-300 Myr). Our sample is composed of 48 low-resolution (R approximate to 100) spectra and 41 moderate-resolution spectra (R greater than or similar to 750-2000). First, we establish a method for spectral typing M5-L7 dwarfs at near-IR wavelengths that is independent of gravity. We find that both visual and index-based classification in the near-IR provides consistent spectral types with optical spectral types, though with a small systematic offset in the case of visual classification at J and K band. Second, we examine features in the spectra of similar to 10 Myr ultracool dwarfs to define a set of gravity-sensitive indices based on FeH, VO, Ki, Nai, and H-band continuum shape. We then create an index-based method for classifying the gravities of M6-L5 dwarfs that provides consistent results with gravity classifications from optical spectroscopy. Our index-based classification can distinguish between young and dusty objects. Guided by the resulting classifications, we propose a set of low-gravity spectral standards for the near-IR. Finally, we estimate the ages corresponding to our gravity classifications.

Sudden Death after Chest Pain: Feasibility of Virtual Autopsy with Postmortem CT Angiography and Biopsy

Purpose:To determine the potential of minimally invasive postmortem computed tomographic (CT) angiography combined with image-guided tissue biopsy of the myocardium and lungs in decedents who were thought to have died of acute chest disease and to compare this method with conventional autopsy as the reference standard.Materials and Methods:The responsible justice department and ethics committee approved this study. Twenty corpses (four female corpses and 16 male corpses; age range, 15-80 years), all of whom were reported to have had antemortem acute chest pain, were imaged with postmortem whole-body CT angiography and underwent standardized image-guided biopsy. The standard included three biopsies of the myocardium and a single biopsy of bilateral central lung tissue. Additional biopsies of pulmonary clots for differentiation of pulmonary embolism and postmortem organized thrombus were performed after initial analysis of the cross-sectional images. Subsequent traditional autopsy with sampling of histologic specimens was performed in all cases. Thereafter, conventional histologic and autopsy reports were compared with postmortem CT angiography and CT-guided biopsy findings. A Cohen k coefficient analysis was performed to explore the effect of the clustered nature of the data.Results:In 19 of the 20 cadavers, findings at postmortem CT angiography in combination with CT-guided biopsy validated the cause of death found at traditional autopsy. In one cadaver, early myocardial infarction of the papillary muscles had been missed. The Cohen κ coefficient was 0.94. There were four instances of pulmonary embolism, three aortic dissections (Stanford type A), three myocardial infarctions, three instances of fresh coronary thrombosis, three cases of obstructive coronary artery disease, one ruptured ulcer of the ascending aorta, one ruptured aneurysm of the right subclavian artery, one case of myocarditis, and one pulmonary malignancy with pulmonary artery erosion. In seven of 20 cadavers, CT-guided biopsy provided additional histopathologic information that substantiated the final diagnosis of the cause of death.Conclusion:Postmortem CT angiography combined with image-guided biopsy, because of their minimally invasive nature, have a potential role in the detection of the cause of death after acute chest pain.© RSNA, 2012.

Determination of creatine and phosphocreatine in muscle biopsy samples by capillary electrophoresis with contactless conductivity detection

A capillary electrophoresis method with contactless conductivity detection was evaluated as a new approach for quantification of creatine and phosphocreatine in human quadriceps femoris biopsy samples. The running buffers employed consisted of 1 M acetic acid at a pH of 2.3 for the determination of creatine and 50 mM 3-(N-morpholino)propanesulfonic acid/30 mM histidine at a pH of 6.4 for the determination of phosphocreatine in the centrifuged muscle extracts. The limits of detection for creatine and phosphocreatine were found to be 2.5 and 1.0 μM, respectively. Creatine and phosphocreatine were determined in six human muscle biopsy samples and the results were found comparable to those of a standard enzymatic assay. The procedures developed for creatine and phosphocreatine also allow the determination of creatinine as well as adenosine diphosphate and adenosine triphosphate.

Fitting of Mixtures with Unspecified Number of Components Using Cross Validation Distance Estimate

Estimation of the number of mixture components (k) is an unsolved problem. Available methods for estimation of k include bootstrapping the likelihood ratio test statistics and optimizing a variety of validity functionals such as AIC, BIC/MDL, and ICOMP. We investigate the minimization of distance between fitted mixture model and the true density as a method for estimating k. The distances considered are Kullback-Leibler (KL) and “L sub 2”. We estimate these distances using cross validation. A reliable estimate of k is obtained by voting of B estimates of k corresponding to B cross validation estimates of distance. This estimation methods with KL distance is very similar to Monte Carlo cross validated likelihood methods discussed by Smyth (2000). With focus on univariate normal mixtures, we present simulation studies that compare the cross validated distance method with AIC, BIC/MDL, and ICOMP. We also apply the cross validation estimate of distance approach along with AIC, BIC/MDL and ICOMP approach, to data from an osteoporosis drug trial in order to find groups that differentially respond to treatment.

Initial clinical experience with the 4-F self-expanding XPERT stent system for infrapopliteal treatment of patients with severe claudication and critical limb ischemia

PURPOSE: To evaluate the primary success and short-term patency associated with a new 4-F sheath-compatible self-expanding nitinol stent after failed conventional angioplasty of distal popliteal and infrapopliteal lesions in severe lifestyle-limiting claudication (LLC) and chronic critical limb ischemia (CLI). MATERIALS AND METHODS: Between May 2003 and July 2005, 35 patients with Rutherford category 3-5 disease (16 patients with CLI, 19 patients with LLC) underwent percutaneous transluminal angioplasty (PTA) and stent implantation. Indications for stent placement were residual stenosis, flow-limiting dissections, or elastic recoil after PTA. Before and after the intervention and during the 6-month follow-up, clinical investigation, color-flow and duplex Doppler ultrasonography, and digital subtraction angiography were performed. Technical success, primary patency at 6 months, clinical improvement as defined by Rutherford with clinical and hemodynamic measures, and complications were evaluated. RESULTS: A total of 22 patients underwent distal popliteal artery stent placement and 13 underwent tibioperoneal artery stent placement. Stent implantation was successfully performed in all patients. After stent placement, the primary cumulative patency rate for the study group at 6 months was 82%. The mean resting ankle-brachial index at baseline was 0.50 +/- 0.16 and significantly increased to 0.90 +/- 0.17 at 12-24 hours after intervention and 0.82 +/- 0.24 at latest follow-up (P < .001 for both). The sustained clinical improvement rate was 80% at the 6-month follow-up. The 6-month limb salvage rate regarding major amputation was 100%. The rate of major complications was 17%. CONCLUSIONS: Infrapopliteal application of the new nitinol stent is a safe, feasible, and effective method with good short-term patency rate in the treatment of severe LLC and chronic CLI.

Computational Techniques for Spatial Logistic Regression with Large Datasets

In epidemiological work, outcomes are frequently non-normal, sample sizes may be large, and effects are often small. To relate health outcomes to geographic risk factors, fast and powerful methods for fitting spatial models, particularly for non-normal data, are required. We focus on binary outcomes, with the risk surface a smooth function of space. We compare penalized likelihood models, including the penalized quasi-likelihood (PQL) approach, and Bayesian models based on fit, speed, and ease of implementation. A Bayesian model using a spectral basis representation of the spatial surface provides the best tradeoff of sensitivity and specificity in simulations, detecting real spatial features while limiting overfitting and being more efficient computationally than other Bayesian approaches. One of the contributions of this work is further development of this underused representation. The spectral basis model outperforms the penalized likelihood methods, which are prone to overfitting, but is slower to fit and not as easily implemented. Conclusions based on a real dataset of cancer cases in Taiwan are similar albeit less conclusive with respect to comparing the approaches. The success of the spectral basis with binary data and similar results with count data suggest that it may be generally useful in spatial models and more complicated hierarchical models.

Inference on Survival Data with Covariate Measurement Error - An Imputation-based Approach

We propose a new method for fitting proportional hazards models with error-prone covariates. Regression coefficients are estimated by solving an estimating equation that is the average of the partial likelihood scores based on imputed true covariates. For the purpose of imputation, a linear spline model is assumed on the baseline hazard. We discuss consistency and asymptotic normality of the resulting estimators, and propose a stochastic approximation scheme to obtain the estimates. The algorithm is easy to implement, and reduces to the ordinary Cox partial likelihood approach when the measurement error has a degenerative distribution. Simulations indicate high efficiency and robustness. We consider the special case where error-prone replicates are available on the unobserved true covariates. As expected, increasing the number of replicate for the unobserved covariates increases efficiency and reduces bias. We illustrate the practical utility of the proposed method with an Eastern Cooperative Oncology Group clinical trial where a genetic marker, c-myc expression level, is subject to measurement error.

Cerebral correlates of muscle tone fluctuations in restless legs syndrome: A pilot study with combined functional magnetic resonance imaging and anterior tibial muscle electromyography

BACKGROUND: The pathology of restless legs syndrome (RLS) is still not understood. To investigate the pathomechanism of the disorder further we recorded a surface electromyogram (EMG) of the anterior tibial muscle during functional magnetic resonance imaging (fMRI) in patients with idiopathic RLS. METHODS: Seven subjects with moderate to severe RLS were investigated in the present pilot study. Patients were lying supine in the scanner for over 50min and were instructed not to move voluntarily. Sensory leg discomfort (SLD) was evaluated on a 10-point Likert scale. For brain image analysis, an algorithm for the calculation of tonic EMG values was developed. RESULTS: We found a negative correlation of tonic EMG and SLD (p <0.01). This finding provides evidence for the clinical experience that RLS-related subjective leg discomfort increases during muscle relaxation at rest. In the fMRI analysis, the tonic EMG was associated with activation in motor and somatosensory pathways and also in some regions that are not primarily related to motor or somatosensory functions. CONCLUSIONS: By using a newly developed algorithm for the investigation of muscle tone-related changes in cerebral activity, we identified structures that are potentially involved in RLS pathology. Our method, with some modification, may also be suitable for the investigation of phasic muscle activity that occurs during periodic leg movements.

Efficient coupling of transducin to monomeric rhodopsin in a phospholipid bilayer

G protein-coupled receptors (GPCRs) are seven transmembrane domain proteins that transduce extracellular signals across the plasma membrane and couple to the heterotrimeric family of G proteins. Like most intrinsic membrane proteins, GPCRs are capable of oligomerization, the function of which has only been established for a few different receptor systems. One challenge in understanding the function of oligomers relates to the inability to separate monomeric and oligomeric receptor complexes in membrane environments. Here we report the reconstitution of bovine rhodopsin, a GPCR expressed in the retina, into an apolipoprotein A-I phospholipid particle, derived from high density lipoprotein (HDL). We demonstrate that rhodopsin, when incorporated into these 10 nm reconstituted HDL (rHDL) particles, is monomeric and functional. Rhodopsin.rHDL maintains the appropriate spectral properties with respect to photoactivation and formation of the active form, metarhodopsin II. Additionally, the kinetics of metarhodopsin II decay is similar between rhodopsin in native membranes and rhodopsin in rHDL particles. Photoactivation of monomeric rhodopsin.rHDL also results in the rapid activation of transducin, at a rate that is comparable with that found in native rod outer segments and 20-fold faster than rhodopsin in detergent micelles. These data suggest that monomeric rhodopsin is the minimal functional unit in G protein activation and that oligomerization is not absolutely required for this process.

Algorithms for spacecraft formation flying navigation based on wireless positioning system measurements

Spacecraft formation flying navigation continues to receive a great deal of interest. The research presented in this dissertation focuses on developing methods for estimating spacecraft absolute and relative positions, assuming measurements of only relative positions using wireless sensors. The implementation of the extended Kalman filter to the spacecraft formation navigation problem results in high estimation errors and instabilities in state estimation at times. This is due tp the high nonlinearities in the system dynamic model. Several approaches are attempted in this dissertation aiming at increasing the estimation stability and improving the estimation accuracy. A differential geometric filter is implemented for spacecraft positions estimation. The differential geometric filter avoids the linearization step (which is always carried out in the extended Kalman filter) through a mathematical transformation that converts the nonlinear system into a linear system. A linear estimator is designed in the linear domain, and then transformed back to the physical domain. This approach demonstrated better estimation stability for spacecraft formation positions estimation, as detailed in this dissertation. The constrained Kalman filter is also implemented for spacecraft formation flying absolute positions estimation. The orbital motion of a spacecraft is characterized by two range extrema (perigee and apogee). At the extremum, the rate of change of a spacecraft’s range vanishes. This motion constraint can be used to improve the position estimation accuracy. The application of the constrained Kalman filter at only two points in the orbit causes filter instability. Two variables are introduced into the constrained Kalman filter to maintain the stability and improve the estimation accuracy. An extended Kalman filter is implemented as a benchmark for comparison with the constrained Kalman filter. Simulation results show that the constrained Kalman filter provides better estimation accuracy as compared with the extended Kalman filter. A Weighted Measurement Fusion Kalman Filter (WMFKF) is proposed in this dissertation. In wireless localizing sensors, a measurement error is proportional to the distance of the signal travels and sensor noise. In this proposed Weighted Measurement Fusion Kalman Filter, the signal traveling time delay is not modeled; however, each measurement is weighted based on the measured signal travel distance. The obtained estimation performance is compared to the standard Kalman filter in two scenarios. The first scenario assumes using a wireless local positioning system in a GPS denied environment. The second scenario assumes the availability of both the wireless local positioning system and GPS measurements. The simulation results show that the WMFKF has similar accuracy performance as the standard Kalman Filter (KF) in the GPS denied environment. However, the WMFKF maintains the position estimation error within its expected error boundary when the WLPS detection range limit is above 30km. In addition, the WMFKF has a better accuracy and stability performance when GPS is available. Also, the computational cost analysis shows that the WMFKF has less computational cost than the standard KF, and the WMFKF has higher ellipsoid error probable percentage than the standard Measurement Fusion method. A method to determine the relative attitudes between three spacecraft is developed. The method requires four direction measurements between the three spacecraft. The simulation results and covariance analysis show that the method’s error falls within a three sigma boundary without exhibiting any singularity issues. A study of the accuracy of the proposed method with respect to the shape of the spacecraft formation is also presented.

Microscopic wire guide-based orotracheal mouse intubation: description, evaluation and comparison with transillumination

Airway access is needed for a number of experimental animal models, and the majority of animal research is based on mouse models. Anatomical conditions in mice are small, and the narrow glottic opening allows intubation only with a subtle technique. We therefore developed a microscopic endotracheal intubation method with a wire guide technique in mice anaesthetized with halothane in oxygen. The mouse is hung perpendicularly with its incisors on a thread fixed on a vertical plate. The tongue is placed with a pair of forceps between the left hand's thumb and forefinger and slightly pulled, while the neck and thorax are positioned using the third and fourth fingers. By doing so, the neck can be slightly stretched, which allows optimal visualization of the larynx and the vocal cords. To ensure a safe intubation, a fine wire guide is placed under vision between the vocal cords and advanced about 5 mm into the trachea. An intravenous 22G x 1 in. plastic or Teflon catheter is guided over this wire. In a series of 41 mice, between 21 and 38 g, the success rate for the first intubation attempt was >95%. Certainty of the judgement procedure was 100% and success rate was higher using the described method when compared with a transillumination method in a further series. The technique is safe, less invasive than tracheostomy and suitable for controlled ventilation and pulmonary substance application.

Hard thermal loop benchmark for the extraction of the nonperturbative QQ[over ¯] potential

The extraction of the finite temperature heavy quark potential from lattice QCD relies on a spectral analysis of the Wilson loop. General arguments tell us that the lowest lying spectral peak encodes, through its position and shape, the real and imaginary parts of this complex potential. Here we benchmark this extraction strategy using leading order hard-thermal loop (HTL) calculations. In other words, we analytically calculate the Wilson loop and determine the corresponding spectrum. By fitting its lowest lying peak we obtain the real and imaginary parts and confirm that the knowledge of the lowest peak alone is sufficient for obtaining the potential. Access to the full spectrum allows an investigation of spectral features that do not contribute to the potential but can pose a challenge to numerical attempts of an analytic continuation from imaginary time data. Differences in these contributions between the Wilson loop and gauge fixed Wilson line correlators are discussed. To better understand the difficulties in a numerical extraction we deploy the maximum entropy method with extended search space to HTL correlators in Euclidean time and observe how well the known spectral function and values for the real and imaginary parts are reproduced. Possible venues for improvement of the extraction strategy are discussed.