Dynamic computer simulation of electrophoretic enantiomer migration order and separation in presence of a neutral cyclodextrin

One-dimensional dynamic computer simulation was employed to investigate the separation and migration order change of ketoconazole enantiomers at low pH in presence of increasing amounts of (2-hydroxypropyl)-β-cyclodextrin (OHP-β-CD). The 1:1 interaction of ketoconazole with the neutral cyclodextrin was simulated under real experimental conditions and by varying input parameters for complex mobilities and complexation constants. Simulation results obtained with experimentally determined apparent ionic mobilities, complex mobilities, and complexation constants were found to compare well with the calculated separation selectivity and experimental data. Simulation data revealed that the migration order of the ketoconazole enantiomers at low (OHP-β-CD) concentrations (i.e. below migration order inversion) is essentially determined by the difference in complexation constants and at high (OHP-β-CD) concentrations (i.e. above migration order inversion) by the difference in complex mobilities. Furthermore, simulations with complex mobilities set to zero provided data that mimic migration order and separation with the chiral selector being immobilized. For the studied CEC configuration, no migration order inversion is predicted and separations are shown to be quicker and electrophoretic transport reduced in comparison to migration in free solution. The presented data illustrate that dynamic computer simulation is a valuable tool to study electrokinetic migration and separations of enantiomers in presence of a complexing agent.

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

The genetic contribution of the major histocompatibility complex to bleomycin-induced lung injury

Pulmonary fibrosis (PF) is the result of a variety of environmental and cancer treatment related insults and is characterized by excessive deposition of collagen. Gas exchange in the alveoli is impaired as the normal lung becomes dense and collapsed leading to a loss of lung volume. It is now accepted that lung injury and fibrosis are in part genetically regulated. ^ Bleomycin is a chemotherapeutic agent used for testicular cancer and lymphomas that induces significant pulmonary toxicity. We delivered bleomycin to mice subcutaneously via a miniosmotic pump in order to elicit lung injury (LI) and quantified the %LI morphometrically using video imaging software. We previously identified a quantitative trait loci, Blmpf-1(LOD=17.4), in the Major Histocompatibility Complex (MHC), but the exact genetic components involved have remained unknown. ^ In the current studies, Blmpf-1 was narrowed to an interval spanning 31.9-32.9Mb on Chromosome 17 using MHC Congenic mice. This region includes the MHC Class II and III genes, and is flanked by the TNF-alpha super locus and MHC Class I genes. Knockout mice of MHC Class I genes (B2mko), MHC Class II genes (Cl2ko), and TNF-alpha (TNF-/-) and its receptors (p55-/-, p75-/-, and p55/p75-/-) were treated with bleomycin in order to ascertain the role of these genes in the pathogenesis of lung injury. ^ Cl2ko mice had significantly better survival and %LI when compared to treated background BL/6 (B6, P<.05). In contrast, B2mko showed no differences in survival or %LI compared to B6. This suggests that the MHC Class II locus contains susceptibility genes for bleomycin-induced lung injury. ^ TNF-alpha, a Class III gene, was examined and it was found that TNF-/- and p55-/- mice had higher %LI and lower survival when compared to B6 (P<.05). In contrast, p75-/- mice had significantly reduced %LI when compared to TNF-/-, p55-/-, and B6 mice as well as higher survival (P<.01). These data contradict the current paradigm that TNF-alpha is a profibrotic mediator of lung injury and suggest a novel and distinct role for the p55 and p75 receptors in mediating lung injury. ^

Role of {\it Mycobacterium avium-intracellulare\/} complex infection in AIDS mortality

Disseminated MAC (dMAC) is the third most prevalent opportunistic infection in AIDS patients. In order to understand the role MAC infection plays in affecting survival of AIDS patients, a cohort of 203 suspected dMAC veterans seen at the Houston Veterans Affairs Medical Center between August 14, 1987 and December 31, 1991 were analyzed. The criteria for suspected dMAC infection was HIV+ men having a CD4+ level $\le$200 cells/mm$\sp3,$ on zidovudine treatment $\ge$1 month and who had any of the following: (a) a confirmed respiratory MAC infection, (b) fever $\ge$101$\sp\circ\rm F$ for $\ge$48 hours, (c) unexplained weight loss of 10 lbs or $\ge$10% BW over 3 months or (d) Hgb $\le$7.5 g/dl or decrease in Hgb $\ge$3.0 g/dl, while on 500-600 mg/day AZT. The study was conducted before the commencement of an effective MAC anti-mycobacterial therapy, so the true course of MAC infection was seen without the confounder of a therapeutic regimen. Kaplan-Meier and Cox regression survival analysis was used to compare 45 MAC culture positive and 118 MAC culture negative veterans. The 1 year survival rate of veterans with documented dMAC infection was 0.37 compared to 0.50 for veterans not acquiring dMAC infection. Significant differences between subgroups were also seen with the variables: PCP prophylaxis, the AIDS indicator disease Candida esophagitis, CD4+ lymphocyte level, CD4 percent lymphocyte level, WBC level, Hgb and Hct levels. Using multivariate modeling, it was determined that PCP prophylaxis (RR = 6.12, CI 2.24-16.68) was a predictor of survival and both CD4% lymphocytes $\le$6.0% (RR = 0.33, CI 0.17-0.68) and WBC level $\le$3000 cells/mm$\sp3$ (RR = 0.60, CI 0.39-0.93) were predictors of mortality. CD4+ level $\le$50 cells/mm$\sp3$ was not a significant predictor of mortality. Although MAC culture status was a significant predictor of mortality in the univariate model, a positive dMAC culture was not a significant predictor of AIDS mortality in the multivariate model. A positive dMAC culture, however, did affect mortality in a stratified analysis when baseline laboratory values were: CD8+ lymphocytes $>$600 cells/mm$\sp3,$ Hgb $>$11.0 g/dl, Hct $>$31.0% and WBC level $>$3000 cells/mm$\sp3.$ ^

Geochemistry of the lower sheeted dike complex of ODP Holes 137-504B and 140-504B

Sixty-three samples representing 379 m of sheeted dikes from Deep Sea Drilling Project/Ocean Drilling Program Site 504B have been analyzed for major and selected trace elements by X-ray fluorescence. The samples range from microcrystalline aphyric basalts to moderately phyric (2%-10% phenocrysts) diabase that are typically multiply saturated with plagioclase, olivine, and clinopyroxene, in order of relative abundance. All analyzed samples are classified as Group D compositions with moderate to slightly elevated compatible elements (MgÆ-value = 0.65% ± 0.03%; Al2O3 = 15.5% ± 0.8%; CaO = 13.0% ± 0.3%; Ni = 114 ± 29 ppm), and unusually depleted levels of moderate to highly incompatible elements (Nb < 1 ppm; Zr = 44 ± 7 ppm; Rb < 0.5 ppm; Ba ~ 1 ppm; P2O5 = 0.07% ± 0.02%). These compositions are consistent with a multistage melting of a normal ocean ridge basaltic mantle source followed by extensive fractionation of olivine, plagioclase, and clinopyroxene. Leg 140 aphyric to sparsely phyric (0%-2% phenocrysts) basalts and diabases are compositionally indistinguishable from similarly phyric samples at higher levels in the hole. An examination of the entire crustal section, from the overlying volcanics through the sheeted dikes observed in Leg 140, reveals no significant trends indicating the enrichment or depletion of Costa Rica Rift Zone source magmas over time. Similarly, significant trends toward increased or decreased differentiation cannot be identified, although compositional patterns reflecting variable amounts of phenocryst addition are apparent at various depths. Below ? 1700 mbsf to the bottom of the Leg 140 section, there is a broadly systematic pattern of Zn depletion with depth, the result of high-temperature hydrothermal leaching. This zone of depletion is thought to be a significant source of Zn for the hydrothermal fluids depositing metal sulfides at ridge-crest hydrothermal vents and the sulfide-mineralization zone, located in the transition between pillow lavas and sheeted dikes. Localized zones of intense alteration (60%-95% recrystallization) are present on a centimeter to meter scale in many lithologic units. Within these zones, normally immobile elements Ti, Zr, Y, and rare-earth elements are strongly depleted compared with "fresher" samples centimeters away. The extent of compositional variability of these elements tends to obscure primary igneous trends if the highly altered samples are not identified or removed. At levels up to 40% (or possibly 60%) recrystallization, Ti, Zr, and Y retain their primary signatures. Although the mechanisms are unclear, it is possible that these intense alteration zones are a source of Y and rare-earth elements for the typically rare-earth-element-enriched hydrothermal vent fluids of mid-ocean ridges.

Geochemistry and isotopic ratios of samples from the Tiger gabbroic complex, Northern Victoria Land, Antarctica

Subduction related mafic/ultramafic complexes marking the suture between the Wilson Terrane and the Bowers Terrane in northern Victoria Land (Antarctica) are well-suited for evaluating the magmatic and structural evolu- tion at the Palaeo-Pacific continental margin of Gondwana. One of these intru- sions is the "Tiger Gabbro Complex" (TGC), which is located at the southern end of the island-arc type Bowers Terrane. The TGC is an early Palaeozoic island-arc related layered igneous complex characterized by extraordinarly fresh sequences of ultramafic, mafic and evolved lithologies and extensive development of high-temperature high-strain zones. The goal of the present study is to establish the kinematic, petrogenetic and temporal development of the TGC in order to evaluate the magmatic and structural evolution of the deep crustal roots of this Cambrian-aged island-arc. Fieldwork during GANOVEX X was carried out to provide insight into: (i) the spatial relations between the different igneous lithologies of the TGC, (ii) the nature of the contact between the TGC and Bowers Terrane, and (iii) the high-temperature shear zones exposed in parts of the TGC. Here, we report the results of detailed field and petrological observations combined with new geochronological data. Based on these new data, we tentatively propose a petrogenetic-kinematic model for the TGC, which involves a two-phase evolution during the Ross orogeny. These phases can be summarized as: (i) an early phase (maximum age c. 530 Ma) involving tectono-magmatic processes that were active at the deep crustal level represented by the TGC within the Bowers island arc and within a general NE-SW directed contractional regime and (ii) a late phase (maximum age c. 490 Ma) attributed to the late Ross orogenic intrusion of the TGC into the higher-crustal metasedimentary country rocks of the Bowers Terrane under NE-SW directed horizontal maximum stress and subsequent cooling.

Seawater carbonate chemistry and concentration boundary layers around complex assemblages of macroalgae in a laboratory experiment

Metabolic processes have the potential to modulate the effects of ocean acidification (OA) in nearshore macroalgal beds. We investigated whether natural mixed assemblages of the articulate coralline macroalgae Arthrocardia corymbosa and understory crustose coralline algae (CCA) altered pH and O2 concentrations within and immediately above their canopies. In a unidirectional flume, we tested the effect of water velocity (0-0.1 m/s), bulk seawater pH (ambient pH 8.05, and pH 7.65), and irradiance (photosynthetically saturating light and darkness) on pH and O2 concentration gradients, and the derived concentration boundary layer (CBL) thickness. At bulk seawater pH 7.65 and slow velocities (0 and 0.015 m/s), pH at the CCA surface increased to 7.90-8.00 in the light. Although these manipulations were short term, this indicates a potential daytime buffering capacity that could alleviate the effects of OA. Photosynthetic activity also increased O2 concentrations at the surface of the CCA. However, this moderating capacity was flow dependent; the CBL thickness decreased from an average of 26.8 mm from the CCA surface at 0.015 m/s to 4.1 mm at 0.04 m/s. The reverse trends occurred in the dark, with respiration causing pH and O2 concentrations to decrease at the CCA surface. At all flow velocities the CBL thicknesses (up to 68 mm) were much greater than those previously published, indicating that the presence of canopies can alter the CBL substantially. In situ, the height of macroalgal canopies can be an order of magnitude larger than those used here, indicating that the degree of buffering to OA will be context dependent.

Coupling time-stepping numerical methods and standard aerodynamics codes for instability analysis of flows in complex geometries

The development of a global instability analysis code coupling a time-stepping approach, as applied to the solution of BiGlobal and TriGlobal instability analysis 1, 2 and finite-volume-based spatial discretization, as used in standard aerodynamics codes is presented. The key advantage of the time-stepping method over matrix-formulation approaches is that the former provides a solution to the computer-storage issues associated with the latter methodology. To-date both approaches are successfully in use to analyze instability in complex geometries, although their relative advantages have never been quantified. The ultimate goal of the present work is to address this issue in the context of spatial discretization schemes typically used in industry. The time-stepping approach of Chiba 3 has been implemented in conjunction with two direct numerical simulation algorithms, one based on the typically-used in this context high-order method and another based on low-order methods representative of those in common use in industry. The two codes have been validated with solutions of the BiGlobal EVP and it has been showed that small errors in the base flow do not have affect significantly the results. As a result, a three-dimensional compressible unsteady second-order code for global linear stability has been successfully developed based on finite-volume spatial discretization and time-stepping method with the ability to study complex geometries by means of unstructured and hybrid meshes

Order 10 4 speedup in global linear instability analysis using matrix formation

A unified solution framework is presented for one-, two- or three-dimensional complex non-symmetric eigenvalue problems, respectively governing linear modal instability of incompressible fluid flows in rectangular domains having two, one or no homogeneous spatial directions. The solution algorithm is based on subspace iteration in which the spatial discretization matrix is formed, stored and inverted serially. Results delivered by spectral collocation based on the Chebyshev-Gauss-Lobatto (CGL) points and a suite of high-order finite-difference methods comprising the previously employed for this type of work Dispersion-Relation-Preserving (DRP) and Padé finite-difference schemes, as well as the Summationby- parts (SBP) and the new high-order finite-difference scheme of order q (FD-q) have been compared from the point of view of accuracy and efficiency in standard validation cases of temporal local and BiGlobal linear instability. The FD-q method has been found to significantly outperform all other finite difference schemes in solving classic linear local, BiGlobal, and TriGlobal eigenvalue problems, as regards both memory and CPU time requirements. Results shown in the present study disprove the paradigm that spectral methods are superior to finite difference methods in terms of computational cost, at equal accuracy, FD-q spatial discretization delivering a speedup of ð (10 4). Consequently, accurate solutions of the three-dimensional (TriGlobal) eigenvalue problems may be solved on typical desktop computers with modest computational effort.

Flow and wakes in large wind farms in complex terrain and offshore

Power losses due to wind turbine wakes are of the order of 10 and 20% of total power output in large wind farms. The focus of this research carried out within the EC funded UPWIND project is wind speed and turbulence modelling for large wind farms/wind turbines in complex terrain and offshore in order to optimise wind farm layouts to reduce wake losses and loads.

Flow and wakes in complex terrain and offshore. Model development and verification in UPWIND

The paper presents research conducted in the Flow workpackage of the EU funded UPWIND project which focuses on improving models for flow within and downwind of large wind farms in complex terrain and offshore. The main activity is modelling the behaviour of wind turbine wakes in order to improve power output predictions.

Improved heuristic drift elimination with magnetically-aided dominant directions (MiHDE) for pedestrian navigation in complex buildings

The main problem of pedestrian dead-reckoning (PDR) using only a body-attached inertial measurement unit is the accumulation of heading errors. The heading provided by magnetometers in indoor buildings is in general not reliable and therefore it is commonly not used. Recently, a new method was proposed called heuristic drift elimination (HDE) that minimises the heading error when navigating in buildings. It assumes that the majority of buildings have their corridors parallel to each other, or they intersect at right angles, and consequently most of the time the person walks along a straight path with a heading constrained to one of the four possible directions. In this article we study the performance of HDE-based methods in complex buildings, i.e. with pathways also oriented at 45°, long curved corridors, and wide areas where non-oriented motion is possible. We explain how the performance of the original HDE method can be deteriorated in complex buildings, and also, how severe errors can appear in the case of false matches with the building's dominant directions. Although magnetic compassing indoors has a chaotic behaviour, in this article we analyse large data-sets in order to study the potential use that magnetic compassing has to estimate the absolute yaw angle of a walking person. Apart from these analysis, this article also proposes an improved HDE method called Magnetically-aided Improved Heuristic Drift Elimination (MiHDE), that is implemented over a PDR framework that uses foot-mounted inertial navigation with an extended Kalman filter (EKF). The EKF is fed with the MiHDE-estimated orientation error, gyro bias corrections, as well as the confidence over that corrections. We experimentally evaluated the performance of the proposed MiHDE-based PDR method, comparing it with the original HDE implementation. Results show that both methods perform very well in ideal orthogonal narrow-corridor buildings, and MiHDE outperforms HDE for non-ideal trajectories (e.g. curved paths) and also makes it robust against potential false dominant direction matchings.

Geometric diagnostics of complex patterns: Spiral defect chaos

Motivated by the observation of spiral patterns in a wide range of physical, chemical, and biological systems, we present an automated approach that aims at characterizing quantitatively spiral-like elements in complex stripelike patterns. The approach provides the location of the spiral tip and the size of the spiral arms in terms of their arc length and their winding number. In addition, it yields the number of pattern components (Betti number of order 1), as well as their size and certain aspects of their shape. We apply the method to spiral defect chaos in thermally driven Rayleigh- Bénard convection and find that the arc length of spirals decreases monotonically with decreasing Prandtl number of the fluid and increasing heating. By contrast, the winding number of the spirals is nonmonotonic in the heating. The distribution function for the number of spirals is significantly narrower than a Poisson distribution. The distribution function for the winding number shows approximately an exponential decay. It depends only weakly on the heating, but strongly on the Prandtl number. Large spirals arise only for larger Prandtl numbers. In this regime the joint distribution for the spiral length and the winding number exhibits a three-peak structure, indicating the dominance of Archimedean spirals of opposite sign and relatively straight sections. For small Prandtl numbers the distribution function reveals a large number of small compact pattern components.

A combined algorithm for analyzing structural controllability and observability of complex networks

In this paper a combined algorithm for analyzing structural controllability and observability of complex networks is presented. The algorithm addresses the two fundamental properties to guarantee structural controllability of a system: the absence of dilations and the accessibility of all nodes. The first problem is reformulated as a Maximum Matching search and it is addressed via the Hopcroft- Karp algorithm; the second problem is solved via a new wiring algorithm. Both algorithms can be combined to efficiently determine the number of required controllers and observers as well as the new required connections in order to guarantee controllability and observability in real complex networks. An application to a Twitter social network with over 100,000 nodes illustrates the proposed algorithms.