Use of generalised additive models to categorise continuous variables in clinical prediction

13 P.

Hybrid Qm/Mm Simulation Of The Hydration Phenomena Of Dipalmitoylphosphatidylcholine Headgroup

The polar headgroup of dipalmitoylphosphatidylcholine (DPPC) molecule both in gas phase and aqueous Solution is investigated by the hybrid quantum mechanical/molecular mechanical (QM/MM) method, in which the polar head of DPPC molecule and the bound water molecules are treated with density functional theory (DFT), while the apolar hydrocarbon chain of DPPC molecule is treated with MM method. It is demonstrated that the hybrid QM/MM method is both accurate and efficient to describe the conformations of DPPC headgroup. Folded structures of headgroup are found in gas phase calculations. In this work, both monohydration and polyhydration phenomena are investigated. In monohydration, different water association sites are studied. Both the hydration energy and the quantum properties of DPPC and water molecules are calculated at the DFT level of theory after geometry optimization. The binding force of monohydration is estimated by using the scan method. In polyhydration, more extended conformations are found and hydration energies in different polyhydration styles are estimated. (C) 2008 Elsevier Inc. All rights reserved.

An Atomistic-Continuum Hybrid Simulation Of Fluid Flows Over Superhydrophobic Surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper, an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces, in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths, which cannot be obtained by molecular dynamics simulation alone.

An atomistic-continuum hybrid simulation of fluid flows over superhydrophobic surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths which cannot be obtained by molecular dynamics simulation alone.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

Sensitivity to alcohols of a planar waveguide ring resonator fabricated by a sol-gel method

A planar waveguide ring resonator was fabricated by organic-inorganic hybrid sol-gel materials; its sensitivity to ethanol vapor was experimentally investigated. It was found that dips in the transmission spectrum of the device shifted to longer wavelengths with increasing the ethanol concentration, and its sensitivity showed a linear relation with the ethanol concentration, showing a coefficient of 1.13 pm/ppm. In addition, the transmission loss of the ring resonator decreased with increasing the ethanol concentration. The measured characteristics suggest that the device may be considered as one of the candidates of alcohol vapor sensors. (c) 2006 Elsevier B.V. All rights reserved.

Mode properties of off-axis one-sided hybrid resonator for novel high average power dual-slab laser

A novel dual-slab laser with off-axis one-sided hybrid resonator is presented. The mode properties of the hybrid resonator are calculated using a fast Fourier transform method (FFT). The influence of wavefront distoration on the output beam quality is considered. Results indicate that the novel dual-slab laser is better than the normal dual-slab laser with off-axis one-sided hybrid resonator.

Optical and surface properties of hybrid TiO2/ormosil planar waveguide prepared by the sol-gel process

Two kinds of silanes, 3-glycidoxypropyltrimethoxysilane (GLYMO) and 3-trimethoxysililpropylmethacrylate (TMSPM), were used to prepare ormosil waveguide films by the sol-gel method. Thirty percent Ti(OBu)(4) and 70% silane were contained in the precursor sets. The properties of films were measured by scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR), UV/VIS/NIR spectrophotometer (UV-vis), atomic force microscopy (AFM), m-line and scattering-detection method. The films from GLYMO and TMSPM precursors exhibit similar thickness (2.58 mu m for GLYMO, 2.51 mu m for TMSPM) and refractive index (1.5438 for GLYMO, 1.5392 for TMSPM, lambda=632.8 nm), but the film from TMSPM precursor has higher propagation loss (1.024 dB/cm, lambda=632.8 nm) than the film prepared from GLYMO (0.569 dB/cm, lambda=632.8 nm). Furthermore, the film prepared from TMSPM is easy to be opaque and cracks during coating whereas the same phenomenon was not found for the film prepared with GLYMO. It is confirmed that GLYMO is a better precursor than TMSPM for waveguide film preparation. (C) 2005 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Effect of water content in sol on optical properties of hybrid sol-gel derived TiO2/SiO2/ormosil film

Sol-gel derived TiO2/SiO2/ormosil hybrid planar waveguides have been deposited on soda-lime glass slides and silicon substrates, films were heat treated at 150 degreesC for 2 h or dried at room temperature. Different amounts of water were added to sols to study their impacts on microstructures and optical properties of films. The samples were characterized by m-line spectroscopy, Fourier transform infrared spectroscopy (FT-IR), UV/VIS/NIR spectrophotometer (UV-vis), atomic force microscopy (AFM), thermal analysis instrument and scattering-detection method. The refractive index was found to have the largest value at the molar ratio H2O/OR = 1 in sol (OR means -OCH3, -OC2H5 and -OC4H9 in the sol), whereas the thickest film appears at H2O/OR = 1/2. The rms surface roughness of all the films is lower than 1.1 nm, and increases with the increase of water content in sol. Higher water content leads to higher attenuation of film. (C) 2004 Elsevier B.V. All rights reserved.

Study of hybrid TiO2/ormosil films doped with laser dyes

TiO2/ormosil films doped with laser dyes have been prepared by the sol-gel method. Spectroscopic properties of the entrapped dyes are studied by the absorption and emission techniques. The results indicate that the absorption and fluorescence spectra of kiton red depend strongly on the properties of the ormosil matrices. The heat-treatment of the kiton red-doped film obviously leads to the increasing fluorescence intensity and the largest fluorescence intensity is obtained after heat-treatment of 150 degrees C for 2 h. However, the fluorescence intensity of the rhodamine 6G-doped film decreases with the increase of the heat-treatment temperature. (c) 2005 Elsevier B.V. All rights reserved.

Preparation and characterization of hybrid ZnO/ormosils films

Hybrid ZnO/ormosils Elms are prepared by the sol-gel method. A FT-IR spectrometer, 900 UV/VIS/NIR spectrophotometer, atomic force microscope, and ellipsometer are employed to investigate microstructure and optical properties of the films fired at different temperatures. The results show that the films with high transmittance and low surface roughness could be obtained at the heat-treatment temperature of 150 degrees C, the refractive index and thickness of the film are 1.413, 2.11 mu m, respectively. Higher temperatures (350 degrees C, 550 degrees C) change the Elm microstructure severely, and then decrease the transmittance of the films.

Prediction and confirmation of seasonal migration of Pacific sardine (Sardinops sagax) in the California Current Ecosystem

During the last century, the population of Pacific sardine (Sardinops sagax) in the California Current Ecosystem has exhibited large fluctuations in abundance and migration behavior. From approximately 1900 to 1940, the abundance of sardine reached 3.6 million metric tons and the “northern stock” migrated from offshore of California in the spring to the coastal areas near Oregon, Washington, and Vancouver Island in the summer. In the 1940s, the sardine stock collapsed and the few remaining sardine schools concentrated in the coastal region off southern California, year-round, for the next 50 years. The stock gradually recovered in the late 1980s and resumed its seasonal migration between regions off southern California and Canada. Recently, a model was developed which predicts the potential habitat for the northern stock of Pacific sardine and its seasonal dynamics. The habitat predictions were successfully validated using data from sardine surveys using the daily egg production method; scientific trawl surveys off the Columbia River mouth; and commercial sardine landings off Oregon, Washington, and Vancouver Island. Here, the predictions of the potential habitat and seasonal migration of the northern stock of sardine are validated using data from “acoustic–trawl” surveys of the entire west coast of the United States during the spring and summer of 2008. The estimates of sardine biomass and lengths from the two surveys are not significantly different between spring and summer, indicating that they are representative of the entire stock. The results also confirm that the model of potential sardine habitat can be used to optimally apply survey effort and thus minimize random and systematic sampling error in the biomass estimates. Furthermore, the acoustic–trawl survey data are useful to estimate concurrently the distributions and abundances of other pelagic fishes.

Gene function prediction from synthetic lethality networks via ranking on demand.

MOTIVATION: Synthetic lethal interactions represent pairs of genes whose individual mutations are not lethal, while the double mutation of both genes does incur lethality. Several studies have shown a correlation between functional similarity of genes and their distances in networks based on synthetic lethal interactions. However, there is a lack of algorithms for predicting gene function from synthetic lethality interaction networks. RESULTS: In this article, we present a novel technique called kernelROD for gene function prediction from synthetic lethal interaction networks based on kernel machines. We apply our novel algorithm to Gene Ontology functional annotation prediction in yeast. Our experiments show that our method leads to improved gene function prediction compared with state-of-the-art competitors and that combining genetic and congruence networks leads to a further improvement in prediction accuracy.

Variance prediction in sea

A method is presented for predicting the variance of the energy levels in a built-up system to accompany the mean values predicted by SEA. Closed form expressions for the variance are obtained in terms of the standard SEA parameters and an additional set of parameters αk that describe the nature of the power input to each subsystem k, and αks that describe the nature of the coupling between subsystems k and s.