Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Ensemble methods”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27Linked to original sources
PubMed1993

Microspheres of human serum albumin with barbiturates: effect of drug partition coefficient on preparation and drug release.

Human serum albumin microspheres were prepared with a series of barbiturates by the thermal denaturation method. The barbiturates were of similar general physicochemical properties but different partition coefficients. The total drug content of microspheres was dependent on the partition coefficient, and an increase in the partition coefficient caused a decrease in the drug content of finished microspheres due to drug migration to the outer organic phase during preparation. The ensemble drug release from microspheres was followed by the paddle method and by potentiometric titration with a pH-stat. Nonlinear regression analysis showed the best fit for the spherical Higuchi equation, especially first-order kinetics (biphasic). The drug partition coefficient affected the release in such a way that drugs with higher partition coefficients were released faster and to a greater extent than those with lower coefficients.

Barbiturates↗
PubMed2006

Predicting protein subcellular locations using hierarchical ensemble of Bayesian classifiers based on Markov chains.

BACKGROUND: The subcellular location of a protein is closely related to its function. It would be worthwhile to develop a method to predict the subcellular location for a given protein when only the amino acid sequence of the protein is known. Although many efforts have been made to predict subcellular location from sequence information only, there is the need for further research to improve the accuracy of prediction. RESULTS: A novel method called HensBC is introduced to predict protein subcellular location. HensBC is a recursive algorithm which constructs a hierarchical ensemble of classifiers. The classifiers used are Bayesian classifiers based on Markov chain models. We tested our method on six various datasets; among them are Gram-negative bacteria dataset, data for discriminating outer membrane proteins and apoptosis proteins dataset. We observed that our method can predict the subcellular location with high accuracy. Another advantage of the proposed method is that it can improve the accuracy of the prediction of some classes with few sequences in training and is therefore useful for datasets with imbalanced distribution of classes. CONCLUSION: This study introduces an algorithm which uses only the primary sequence of a protein to predict its subcellular location. The proposed recursive scheme represents an interesting methodology for learning and combining classifiers. The method is computationally efficient and competitive with the previously reported approaches in terms of prediction accuracies as empirical results indicate. The code for the software is available upon request.

Amino Acid Sequence↗
PubMed2001

Quantum tunneling using entangled classical trajectories.

In this Letter, we present a new method for simulating quantum processes in the context of classical molecular dynamics simulations. The approach is based on solving numerically the quantum Liouville equation in the Wigner representation using ensembles of classical trajectories. Quantum effects arise in this formulation as a breakdown of the statistical independence of the ensemble. New interaction forces between ensemble members are derived, which require the trajectory ensemble representing the state to evolve as an entangled, unified whole. The method is applied to the simulation of quantum tunneling in a one-dimensional model system, yielding excellent agreement with exact quantum calculations.

Journal Article↗
PubMed2005

Melting temperature of ice Ih calculated from coexisting solid-liquid phases.

We carried out molecular-dynamics simulations by using the two-phase coexistence method with the constant pressure, particle number, and enthalpy ensemble to compute the melting temperature of proton-disordered hexagonal ice I(h) at 1-bar pressure. Four models of water were considered, including the widely used TIP4P [W. L. Jorgensen, J. Chandrasekha, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys.79, 926 (1983)] and TIP5P [M. W. Mahoney and W. L. Jorgensen J. Chem. Phys.112, 8910 (2000)] models, as well as recently improved TIP4P and TIP5P models for use with Ewald techniques-the TIP4P-Ew [W. Horn, W. C. Swope, J. W. Pitera, J. C. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon, J. Chem. Phys.120, 9665 (2004)] and TIP5P-Ew [S. W. Rick, J. Chem. Phys.120, 6085 (2004)] models. The calculated melting temperature at 1 bar is T(m) = 229 +/- 1 K for the TIP4P and T(m) = 272.0 +/- 0.6 K for the TIP5P ice I(h), both are consistent with previous simulations based on free-energy methods. For the TIP4P-Ew and TIP5P-Ew models, the calculated melting temperature is T(m) = 257.0 +/- 1.1 K and T(m) = 253.9 +/- 1.1 K, respectively.

Journal Article↗
PubMed2004

The canonical ensemble via symplectic integrators using Nosé and Nosé-Poincaré chains.

Simulations that sample from the canonical ensemble can be generated by the addition of a single degree of freedom, provided that the system is ergodic, as described by Nosé with subsequent modifications by Hoover to allow sampling in real time. Nosé-Hoover dynamics is not ergodic for small or stiff systems and the addition of auxiliary thermostats is needed to overcome this deficiency. Nosé-Hoover dynamics, like its derivatives, does not have a Hamiltonian structure, precluding the use of symplectic integrators which are noted for their long term stability and structure preservation. As an alternative to Nosé-Hoover, the Hamiltonian Nosé-Poincaré method was proposed by Bond, Laird, and Leimkuhler [J. Comput. Phys. 151, 114 (1999)], but the straightforward addition of thermostatting chains does not sample from the canonical ensemble. In this paper a method is proposed whereby additional thermostats can be applied to a Hamiltonian system while retaining sampling from the canonical ensemble. This technique has been used to construct thermostatting chains for the Nosé and Nosé-Poincaré methods.

Journal Article↗
PubMed2005

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. II. Simulation results.

The results of a series of constant pressure and temperature molecular-dynamics (MD) simulation studies based on the rigorous shell particle formulation of the isothermal-isobaric (NpT) ensemble are presented. These MD simulations validate the newly proposed constant pressure equations of motion in which a "shell" particle is used to define uniquely the volume of the system [M. J. Uline and D. S. Corti, J. Chem. Phys. (to be published), preceding paper]. Ensemble averages obtained with the new MD NpT algorithm match the ensemble averages obtained using the previously derived shell particle Monte Carlo NpT method [D. S. Corti, Mol. Phys. 100, 1887 (2002)]. In addition, we also verify that the Hoover NpT MD algorithm [W. G. Hoover, Phys. Rev. A 31, 1695 (1985); 34, 2499 (1986)] generates the correct ensemble averages, though only when periodic boundary conditions are employed. The extension of the shell particle MD algorithm to multicomponent systems is also discussed, in which we show for equilibrium properties that the identity of the shell particle is completely arbitrary when periodic boundary conditions are applied. Self-diffusion coefficients determined with the shell particle equations of motion are also identical to those obtained in other ensembles. Finally, since the mass of the shell particle is known, the system itself, and not a piston of arbitrary mass, controls the time scales for internal pressure and volume fluctuations. We therefore consider the effects of the shell particle on the dynamics of the system. Overall, the shell particle MD algorithm is an effective simulation method for studying systems exposed to a constant external pressure and may provide an advantage over other existing constant pressure approaches when developing nonequilibrium MD methods.

Journal Article↗
PubMed2005

Single cell manipulation, analytics, and label-free protein detection in microfluidic devices for systems nanobiology.

Single cell analytics for proteomic analysis is considered a key method in the framework of systems nanobiology which allows a novel proteomics without being subjected to ensemble-averaging, cell-cycle, or cell-population effects. We are currently developing a single cell analytical method for protein fingerprinting combining a structured microfluidic device with latest optical laser technology for single cell manipulation (trapping and steering), free-solution electrophoretical protein separation, and (label-free) protein detection. In this paper we report on first results of this novel analytical device focusing on three main issues. First, single biological cells were trapped, injected, steered, and deposited by means of optical tweezers in a poly(dimethylsiloxane) microfluidic device and consecutively lysed with SDS at a predefined position. Second, separation and detection of fluorescent dyes, amino acids, and proteins were achieved with LIF detection in the visible (VIS) (488 nm) as well as in the deep UV (266 nm) spectral range for label-free, native protein detection. Minute concentrations of 100 fM injected fluorescein could be detected in the VIS and a first protein separation and label-free detection could be achieved in the UV spectral range. Third, first analytical experiments with single Sf9 insect cells (Spodoptera frugiperda) in a tailored microfluidic device exhibiting distinct electropherograms of a green fluorescent protein-construct proved the validity of the concept. Thus, the presented microfluidic concept allows novel and fascinating single cell experiments for systems nanobiology in the future.

Animals↗
PubMed2005

Rosenbluth-sampled nonequilibrium work method for calculation of free energies in molecular simulation.

We present methods that introduce concepts from Rosenbluth sampling [M. N. Rosenbluth and A. W. Rosenbluth, J. Chem. Phys. 23, 356 (1955)] into the Jarzynski nonequilibrium work (NEW) free-energy calculation technique [C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997)]. The proposed hybrid modifies the way steps are taken in the NEW process. With it, each step is selected from a range of alternatives, with bias given to steps that contribute the least work. The definition of the work average is modified to account for the bias. We introduce two variants of this method, lambda-bias sampling and configuration-bias sampling, respectively; a combined lambda- and configuration-bias method is also considered. By reducing the likelihood that large nonequilibrated work values enter the ensemble average, the Rosenbluth sampling aids in remedying problems of inaccuracy of the calculation. We demonstrate the performance of the proposed methods through a model system of N independent harmonic oscillators. This model captures the difficulties involved in calculating free energies in real systems while retaining many tractable features that are helpful to the study. We examine four variants of this model that differ qualitatively in the nature of their phase-space overlap. Results indicate that the lambda-bias sampling method is most useful for systems with entropic sampling barriers, while the configuration-bias methods are best for systems with energetic sampling barriers. The Rosenbluth-sampling schemes yield much more accurate results than the unbiased nonequilibrium work method. Typically the accuracy can be improved by about an order of magnitude for a given amount of sampling; this improvement translates into two or more orders of magnitude less sampling required to obtain a given level of accuracy, owing to the generally slow convergence of the NEW calculation when the inaccuracy is large.

Algorithms↗
PubMed2002

Determination and quantification of pharmacological, physiological, or behavioral manipulations on ensembles of simultaneously recorded neurons in functionally related neural circuits.

The present report describes methods for evaluating the impact of physiological, pharmacological or behavioral manipulations on simultaneously recorded single neurons within a functional sensory network of the awake, freely moving rat. Surgical techniques were developed to implant a subcutaneous electrode at the base of a single facial whisker (mystacial vibrissae) so that uniform electrical stimuli could be routinely delivered to a discrete region of the whisker pad in the awake and freely moving animal. Multi-channel extracellular recording was used to monitor the spike train activity from ensembles of single neurons in whisker-related regions of the thalamus and neocortex. Algorithms were developed to verify the stability of individual cell recordings during extended experimental sessions. Additional analysis procedures and criteria were established for identifying and evaluating the treatment-specific changes in single neuron discharge patterns that are likely to occur under these experimental conditions. Finally, analyses for evaluating the impact of experimental manipulations on sensory representations distributed over populations of neurons are discussed. The development of these techniques has provided us with the means to investigate the influence of systemically administered drugs or broadly projecting monoamine pathways on single neurons and local circuits within primary sensory networks of the awake or anesthetized mammalian brain.

Action Potentials↗
PubMed2006

Direct calculation of solid-vapor coexistence points by thermodynamic integration: application to single component and binary systems.

We present a new thermodynamic integration method that directly connects the vapor and solid phases by a reversible path. The thermodynamic integration in the isothermal-isobaric ensemble yields the Gibbs free energy difference between the two phases, from which the sublimation temperature can be easily calculated. The method extends to the binary mixture without any modification to the integration path simply by employing the isothermal-isobaric semigrand ensemble. The thermodynamic integration, in this case, yields the chemical potential difference between the solid and vapor phases for one of the components, from which the binary sublimation temperature can be calculated. The coexistence temperatures predicted by our method agree well with those in the literature for single component and binary Lennard-Jones systems.

Journal Article↗
PubMed2005

Isobaric-isothermal monte carlo simulations from first principles: application to liquid water at ambient conditions.

A series of first-principles Monte Carlo simulations in the isobaric-isothermal ensemble were carried out for liquid water at ambient conditions (T=298 K and p=1 atm). The Becke-Lee-Yang-Parr (BLYP) exchange and correlation energy functionals and norm-conserving Goedecker-Teter-Hutter (GTH) pseudopotentials were employed with the CP2 K simulation package to examine systems consisting of 64 water molecules. The fluctuations in the system volume encountered in simulations in the isobaric-isothermal ensemble require a reconsideration of the suitability of the typical charge-density cutoff and the regular grid-generation method previously used for the computation of the electrostatic energy in first-principles simulations in the microcanonical or canonical ensembles. In particular, it is noted that a much higher cutoff is needed and that the most computationally efficient method of creating grids can result in poor simulations. Analysis of the simulation trajectories using a very large charge-density cutoff at 1200 Ry and four different grid-generation methods point to a significantly underestimated liquid density of about 0.8 g cm-3 resulting in a somewhat understructured liquid (with a value of about 2.7 for the height of the first peak in the oxygen-oxygen radial distribution function) for BLYP-GTH water at ambient conditions. In addition, a simulation using a charge-density cutoff at 280 Ry yields a higher density of 0.9 g cm-3, showing the sensitivity of the simulation outcome to this parameter.

Journal Article↗
PubMed2004

Dynamic analyses of information encoding in neural ensembles.

Neural spike train decoding algorithms and techniques to compute Shannon mutual information are important methods for analyzing how neural systems represent biological signals. Decoding algorithms are also one of several strategies being used to design controls for brain-machine interfaces. Developing optimal strategies to design decoding algorithms and compute mutual information are therefore important problems in computational neuroscience. We present a general recursive filter decoding algorithm based on a point process model of individual neuron spiking activity and a linear stochastic state-space model of the biological signal. We derive from the algorithm new instantaneous estimates of the entropy, entropy rate, and the mutual information between the signal and the ensemble spiking activity. We assess the accuracy of the algorithm by computing, along with the decoding error, the true coverage probability of the approximate 0.95 confidence regions for the individual signal estimates. We illustrate the new algorithm by reanalyzing the position and ensemble neural spiking activity of CA1 hippocampal neurons from two rats foraging in an open circular environment. We compare the performance of this algorithm with a linear filter constructed by the widely used reverse correlation method. The median decoding error for Animal 1 (2) during 10 minutes of open foraging was 5.9 (5.5) cm, the median entropy was 6.9 (7.0) bits, the median information was 9.4 (9.4) bits, and the true coverage probability for 0.95 confidence regions was 0.67 (0.75) using 34 (32) neurons. These findings improve significantly on our previous results and suggest an integrated approach to dynamically reading neural codes, measuring their properties, and quantifying the accuracy with which encoded information is extracted.

Action Potentials↗
PubMed2005

G tolerance and pulmonary effects of removing chest counterpressure during pressure breathing.

BACKGROUND: In agile fighter aircraft positive pressure breathing is commonly used as part of the anti-G ensemble. To optimize G protection and prevent over-distention of the lungs, increased airway pressure is balanced by applying a counterpressure to the chest. The aim was to investigate the efficacy of chest counterpressure. METHODS: Three series of experiments were performed using the anti-G ensemble of the 39 Gripen aircraft (AGE-39) and exposing the subjects to 20-s G time profiles; in the first (n = 12) up to +8.0 Gz, in the second (n = 9) up to + 9.0 Gz, and also to simulated aerial combat maneuvers (SACM). Central and peripheral vision, arterial and airway pressures, pressure in the lower portion of the esophagus, and chest wall distension were measured. In the third series, six subjects were exposed to up to +7.0 Gz and esophageal pressure was measured in the upper thorax. In all series, two conditions were compared: with and without pressurized chest bladder. RESULTS: During the 20-s profiles arterial and esophageal pressures, chest wall distension, and visual impairment were similar with and without pressurized chest bladder. Upper esophageal pressure was slightly higher by 10-24% with than without chest bladder (p = 0.03). During SACM, time to exhaustion and the level of perceived exertion were similar with and without pressurized chest bladder. SUMMARY: The results suggest that the chest counterpressure can be removed from the AGE-39 without diminishing G tolerance or G endurance or significantly increasing the risk of lung parenchyma disruption.

Adult↗
PubMed2005

The chemistry of organic nanomaterials.

The development of nanotechnology using organic materials is one of the most intellectually and commercially exciting stories of our times. Advances in synthetic chemistry and in methods for the investigation and manipulation of individual molecules and small ensembles of molecules have produced major advances in the field of organic nanomaterials. The new insights into the optical and electronic properties of molecules obtained by means of single-molecule spectroscopy and scanning probe microscopy have spurred chemists to conceive and make novel molecular and supramolecular designs. Methods have also been sought to exploit the properties of these materials in optoelectronic devices, and prototypes and models for new nanoscale devices have been demonstrated. This Review aims to show how the interaction between synthetic chemistry and spectroscopy has driven the field of organic nanomaterials forward towards the ultimate goal of new technology.

Dendrimers↗
PubMed2001

Estimating spatio-temporal receptive fields of auditory and visual neurons from their responses to natural stimuli.

We present a generalized reverse correlation technique that can be used to estimate the spatio-temporal receptive fields (STRFs) of sensory neurons from their responses to arbitrary stimuli such as auditory vocalizations or natural visual scenes. The general solution for STRF estimation requires normalization of the stimulus-response cross-correlation by the stimulus autocorrelation matrix. When the second-order stimulus statistics are stationary, normalization involves only the diagonal elements of the Fourier-transformed auto-correlation matrix (the power spectrum). In the non-stationary case normalization requires the entire auto-correlation matrix. We present modelling studies that demonstrate the feasibility and accuracy of this method as well as neurophysiological data comparing STRFs estimated using natural versus synthetic stimulus ensembles. For both auditory and visual neurons, STRFs obtained with these different stimuli are similar, but exhibit systematic differences that may be functionally significant. This method should be useful for determining what aspects of natural signals are represented by sensory neurons and may reveal novel response properties of these neurons.

Acoustic Stimulation↗
PubMed2004

Dynamic imaging of cell, extracellular matrix, and tissue movements during avian vertebral axis patterning.

Vertebrate axis patterning depends on cell and extracellular matrix (ECM) repositioning and proper cell-ECM interactions. However, there are few in vivo data addressing how large-scale tissue deformations are coordinated with the motion of local cell ensembles or the displacement of ECM constituents. Combining the methods of dynamic imaging and experimental biology allows both cell and ECM fate-mapping to be correlated with ongoing tissue deformations. These fate-mapping studies suggest that the axial ECM components "move" both as a composite meshwork and as autonomous particles, depending on the length scale being examined. Cells are also part of this composite, and subject to passive displacements resulting from tissue deformations. However, in contrast to the ECM, cells are self-propelled. The net result of cell and ECM displacements, along with proper ECM-cell adhesion, is the assembly of new tissue architecture. Data herein show that disruption of normal cell-ECM interactions during axis formation results in developmental abnormalities and a disorganization of the ECM. Our goal in characterizing the global displacement patterns of axial cells and ECM is to provide critical information regarding existing strain fields in the segmental plate and paraxial mesoderm. Deducing the mechanical influences on cell behavior is critical, if we are to understand vertebral axis patterning. Supplementary material for this article is available online at http://www.mrw.interscience.wiley.com/suppmat/1542-975X/suppmat/72/v72.266.html.

Animals↗
PubMed2003

Dynamical clustering of red blood cells in capillary vessels.

We have modeled the dynamics of a 3-D system consisting of red blood cells (RBCs), plasma and capillary walls using a discrete-particle approach. The blood cells and capillary walls are composed of a mesh of particles interacting with harmonic forces between nearest neighbors. We employ classical mechanics to mimic the elastic properties of RBCs with a biconcave disk composed of a mesh of spring-like particles. The fluid particle method allows for modeling the plasma as a particle ensemble, where each particle represents a collective unit of fluid, which is defined by its mass, moment of inertia, translational and angular momenta. Realistic behavior of blood cells is modeled by considering RBCs and plasma flowing through capillaries of various shapes. Three types of vessels are employed: a pipe with a choking point, a curved vessel and bifurcating capillaries. There is a strong tendency to produce RBC clusters in capillaries. The choking points and other irregularities in geometry influence both the flow and RBC shapes, considerably increasing the clotting effect. We also discuss other clotting factors coming from the physical properties of blood, such as the viscosity of the plasma and the elasticity of the RBCs. Modeling has been carried out with adequate resolution by using 1 to 10 million particles. Discrete particle simulations open a new pathway for modeling the dynamics of complex, viscoelastic fluids at the microscale, where both liquid and solid phases are treated with discrete particles. Figure A snapshot from fluid particle simulation of RBCs flowing along a curved capillary. The red color corresponds to the highest velocity. We can observe aggregation of RBCs at places with the most stagnant plasma flow.

Animals↗
PubMed2005

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein.

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.

Allosteric Regulation↗