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 289 records · Page 16Linked to original sources
PubMed2004

Replica-exchange extensions of simulated tempering method.

In this paper we consider combinations of two well-known generalized-ensemble algorithms, namely, simulated tempering and replica-exchange method. We discuss two examples of such combinations. One is the replica-exchange simulated tempering and the other is the simulated tempering replica-exchange method. In the former method, a short replica-exchange simulation is first performed and the simulated tempering weight factor is obtained by the multiple-histogram reweighting techniques. This process of simulated tempering weight factor determination is faster and simpler than that in the usual iterative process. A long simulated tempering production run is then performed with this weight factor. The latter method is a further extension of the former in which a simulated tempering replica-exchange simulation is performed with a small number of replicas. These algorithms are particularly useful for studying frustrated systems with rough energy landscape. We give the formulations of these two methods in detail and demonstrate their effectiveness taking the example of the system of a 17-residue helical peptide.

Journal Article↗
PubMed2004

Partial rigid-body dynamics in NPT, NPAT and NPgammaT ensembles for proteins and membranes.

A partial rigid-body method of molecular dynamics simulations for proteins and membranes is presented. In this method, the symplectic integrator for rigid bodies is combined with the equations of motion for the NPT ensemble. The standard NPT ensemble is extended to the membrane-specific ensembles, the NPAT (constant normal pressure and lateral surface area of membranes and constant temperature) and NPgammaT (constant normal pressure and lateral surface tension of membranes and constant temperature) ensembles. By more than 30-ns simulations of aqueous proteins and hydrated lipid bilayers, the results of the partial rigid-body method demonstrated excellent conservation of total energy and consistent behavior with the traditional constraint method in terms of structural distribution and fluctuation of proteins and lipids. The efficient implementation of the partial rigid-body method in parallel computation is presented, which is shown to work well in large-scale molecular dynamics simulations.

Algorithms↗
PubMed2004

Efficient fluorescence labeling of a large RNA through oligonucleotide hybridization.

We present an efficient method of introducing fluorophore labels at selected locations in a large RNA. The method is based on specific and highly efficient hybridization between a fluorophore-containing DNA oligonucleotide and a modular hairpin loop replacing a functionally unimportant hairpin loop in the RNA. We demonstrate its feasibility using a 255-nucleotide RNA derived from the catalytic domain of RNase P from Bacillus subtilis. Hybridization of the DNA oligonucleotide to the modular hairpin loop minimally perturbs the structure and function of this RNA. This labeling scheme should be applicable in studies of RNA conformational dynamics by ensemble and single molecule fluorescence methods.

Bacillus subtilis↗
PubMed2005

Creating order from random fluctuations in small spin ensembles.

We demonstrate the ability to create spin order by using a magnetic resonance force microscope to harness the naturally occurring statistical fluctuations in small ensembles of electron spins. In one method, we hyperpolarized the spin system by selectively capturing the transient spin order created by the statistical fluctuations. In a second method, we took a more active approach and rectified the spin fluctuations by applying real-time feedback to the entire spin ensemble. The created spin order can be stored in the laboratory frame for a period on the order of the longitudinal relaxation time of 30 seconds and then read out.

Journal Article↗
PubMed1993

SESAME: a least-squares approach to the evaluation of protein structures computed from NMR data.

A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.

Magnetic Resonance Spectroscopy↗
PubMed2005

A method for the automated detection of venous gas bubbles in humans using empirical mode decomposition.

Doppler ultrasound signals are widely used to grade the quantity of circulating venous bubbles in divers. Current techniques rely on trained observers, making the grading process both time-consuming and subjective. The automated detection of bubbles, however, is confounded by the presence of other signals, primarily those arising from blood motion. Empirical Mode Decomposition was used here to calculate the intrinsic mode functions (IMFs) of a number of Doppler ultrasound signals from recreational divers, post-decompression. The IMFs provide a basis set for signal decomposition, each IMF corresponding to a different timescale in the signal. Each signal was found to comprise approximately 20 IMFs: the precise number being dependent upon the nature of the signal. A method is presented to detect bubbles using the IMF; features are first identified in the individual heart cycles, these having been previously determined using a robust peak detection method, by examining deviations from the ensemble averaged IMF. Bubbles are then identified as features appearing in more than one IMF, with significant energy in the original signal. This method has been applied to a subset of the available database and appears to perform with good sensitivity even when the signal has variable signal strength.

Artificial Intelligence↗
PubMed1997

Level detection in ion channel records via idealization by statistical filtering and likelihood optimization.

A parameter-free method is presented for the level detection in ion channel records via recovery of step wise current changes. No assumptions about ion channel mechanism are made. The primary detection of the transitions is made by statistical filtering the data using the Student's t-test. The event currents are calculated as the average value of the current between two adjacent transitions. An optimal ideal trace is found by maximization of a likelihood function. The distribution of event currents recovered from the raw data is then analysed, again by using the Student's t-test, for their grouping into separate statistical ensembles, defining current levels. The method is subjected to rigorous test using simulated data, and is compared with several other methods. It produces the levels of channel current, their noise amplitudes and distributions of dwell times, the desired information for constructing the channel mechanism.

Ion Channels↗
PubMed2006

Equilibrium distributions of dipalmitoyl phosphatidylcholine and dilauroyl phosphatidylcholine in a mixed lipid bilayer: atomistic semigrand canonical ensemble simulations.

Conventional molecular dynamics (MD) simulations are seriously limited by the slow rate of diffusive mixing in their ability to predict lateral distributions of different lipid types within mixed-lipid bilayers using atomistic models. A method to overcome this limitation, using configuration-bias Monte Carlo (MC) "mutation" moves to transform lipids from one type to another in dynamic equilibrium, is demonstrated in binary fluid-phase mixtures of lipids whose tails differ in length by four carbons. The hybrid MC-MD method operates within a semigrand canonical ensemble, so that an equilibrium composition of the mixture is determined by a constant difference in chemical potential (Delta(mu)) chosen for the components. Within several nanoseconds, bilayer structures initiated as pure dipalmitoyl phosphatidylcholine (DPPC) or pure dilauroyl phosphatidylcholine (DLPC) converge to a common composition and structure in independent simulations conducted at the same Delta(mu). Trends in bilayer thickness, area per lipid, density distributions across the bilayer, and order parameters have been investigated at three mixture compositions and compared with results from the pure bilayers at 323 K. The mixtures showed a moderate increase in DPPC acyl tail sites crossing the bilayer midplane relative to pure DPPC. Correlations between lateral positions of the two lipid types within or across the bilayer were found to be weak or absent. While the lateral distribution is consistent with nearly ideal mixing, the dependence of composition on Delta(mu) indicates a positive excess free energy of mixing.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1999

Multiple structural alignment and core detection by geometric hashing.

A Multiple Structural Alignment algorithm is presented. The algorithm accepts an ensemble of protein structures and finds the largest substructure (core) of C alpha atoms whose geometric configuration appear in all the molecules of the ensemble (core). Both the detection of this core and the resulting structural alignment are done simultaneously. Other large enough multistructural superimpositions are detected as well. Our method is based on the Geometric Hashing paradigm and a superimposition clustering technique which represents superimpositions by sets of matching atoms. The algorithm proved to be efficient on real data in a series of experiments. The same method can be applied to any ensemble of molecules (not necessarily proteins) since our basic technique is sequence order independent.

Algorithms↗
PubMed1997

OLDERADO: on-line database of ensemble representatives and domains. On Line Database of Ensemble Representatives And DOmains.

In cases where the structure of a single protein is represented by an ensemble of conformations, there is often a need to determine the common features and to choose a "representative" conformation. This occurs, for example, with structures determined by NMR spectroscopy, analysis of the trajectory from a molecular dynamics simulation, or an ensemble of structures produced by comparative modeling. We reported previously automatic methods for (1) defining the atoms with low spatial variance across an ensemble (i.e., the "core" atoms) and the domains in which these atoms lie, and (2) clustering an ensemble into conformationally related subfamilies. To extend the utility of these methods, we have developed a freely available server on the World Wide Web at http:/(/)neon.chem.le.ac.uk/olderado/. This (1) contains an automatically generated database of representative structures, core atoms, and domains determined for 449 ensembles of NMR-derived protein structures in the Protein Data Bank (PDB) in May 1997, and (2) allows the user to upload a PDB-formatted file containing the coordinates of an ensemble of structures. The server returns in real time: (1) information on the residues constituting domains: (2) the structures that constitute each conformational subfamily; and (3) an interactive java-based three-dimensional viewer to visualise the domains and clusters. Such information is useful, for example, when selecting conformations to be used in comparative modeling and when choosing parts of structures to be used in molecular replacement. Here we describe the OLDERADO server.

Computer Communication Networks↗
PubMed2005

Effective desynchronization by nonlinear delayed feedback.

We show that nonlinear delayed feedback opens up novel means for the control of synchronization. In particular, we propose a demand-controlled method for powerful desynchronization, which does not require any time-consuming calibration. Our technique distinguishes itself by its robustness against variations of system parameters, even in strongly coupled ensembles of oscillators. We suggest our method for mild and effective deep brain stimulation in neurological diseases characterized by pathological cerebral synchronization.

Deep Brain Stimulation↗
PubMed2005

Spatially adaptive grand canonical ensemble Monte Carlo simulations.

A spatially adaptive Monte Carlo method is introduced directly from the underlying microscopic mechanisms, which satisfies detailed balance, gives the correct noise, and describes accurately dynamic and equilibrium states for adsorption-desorption (grand canonical ensemble) processes. It enables simulations of large scales while capturing sharp gradients with molecular resolution at significantly reduced computational cost. A posteriori estimates, in the sense used in finite-elements methods, are developed for assessing errors (information loss) in coarse-graining and guiding mesh generation.

Journal Article↗
PubMed2002

Automated docking to multiple target structures: incorporation of protein mobility and structural water heterogeneity in AutoDock.

Protein motion and heterogeneity of structural waters are approximated in ligand-docking simulations, using an ensemble of protein structures. Four methods of combining multiple target structures within a single grid-based lookup table of interaction energies are tested. The method is evaluated using complexes of 21 peptidomimetic inhibitors with human immunodeficiency virus type 1 (HIV-1) protease. Several of these structures show motion of an arginine residue, which is essential for binding of large inhibitors. A structural water is also present in 20 of the structures, but it must be absent in the remaining one for proper binding. Mean and minimum methods perform poorly, but two weighted average methods permit consistent and accurate ligand docking, using a single grid representation of the target protein structures.

Computer Simulation↗
PubMed2005

A microarray data analysis framework for postmortem tissues.

This paper will give a complete methodological approach to the processing of oligonucleotide microarray data from postmortem tissue, particularly brain matter. Attention will be drawn to each of the important stages in the process; specifically the quality control, gene expression value calculation, multiple hypothesis testing and correlation analyses. We shall initially discuss the theoretical foundations of each individual method and subsequently apply the ensemble to a sample data set to illustrate and visualise important points.

Algorithms↗
PubMed2003

Computational simulation of the statistical properties of unfolded proteins.

A simple Monte Carlo method was used to generate ensembles of simulated polypeptide conformations that are restricted only by steric repulsion. The models used for these simulations were based on the sequences of four real proteins, ranging in size from 26 to 268 amino acid residues, and included all non-hydrogen atoms. Two sets of calculations were performed, one that included only intra-residue steric repulsion terms and those between adjacent residues, and one that included repulsion terms between all possible atom pairs, so as to explicitly account for the excluded volume effect. Excluded volume was found to increase the average radius of gyration of the chains by 20-40%, with the expansion factor increasing with chain length. Contrary to recent suggestions, however, the excluded volume effect did not greatly restrict the distribution of dihedral angles or favor native-like topologies. The average dimensions of the ensembles calculated with excluded volume were consistent with those measured experimentally for unfolded proteins of similar sizes under denaturing conditions, without introducing any adjustable scaling factor. The simulations also reproduced experimentally determined effective concentrations for the formation of disulfide bonds in reduced and unfolded proteins. The statistically generated ensembles included significant numbers of conformations that were nearly as compact as the corresponding native proteins, as well as many that were as accessible to solvent as a fully extended chain. On the other hand, conformations with as much buried surface area as the native proteins were very rare, as were highly extended conformations. These results suggest that the overall properties of unfolded proteins can be usefully described by a random coil model and that an unfolded polypeptide can undergo significant collapse while losing only a relatively small fraction of its conformational entropy.

Biophysical Phenomena↗
PubMedDate not supplied

Statistically based reduced representation of amino acid side chains.

Preferred conformations of amino acid side chains have been well established through statistically obtained rotamer libraries. Typically, these provide bond torsion angles allowing a side chain to be traced atom by atom. In cases where it is desirable to reduce the complexity of a protein representation or prediction, fixing all side-chain atoms may prove unwieldy. Therefore, we introduce a general parametrization to allow positions of representative atoms (in the present study, these are terminal atoms) to be predicted directly given backbone atom coordinates. Using a large, culled data set of amino acid residues from high-resolution protein crystal structures, anywhere from 1 to 7 preferred conformations were observed for each terminal atom of the non-glycine residues. Side-chain length from the backbone C(alpha) is one of the parameters determined for each conformation, which should itself be useful. Prediction of terminal atoms was then carried out for a second, nonredundant set of protein structures to validate the data set. Using four simple probabilistic approaches, the Monte Carlo style prediction of terminal atom locations given only backbone coordinates produced an average root mean-square deviation (RMSD) of approximately 3 A from the experimentally determined terminal atom positions. With prediction using conditional probabilities based on the side-chain chi(1) rotamer, this average RMSD was improved to 1.74 A. The observed terminal atom conformations therefore provide reasonable and potentially highly accurate representations of side-chain conformation, offering a viable alternative to existing all-atom rotamers for any case where reduction in protein model complexity, or in the amount of data to be handled, is desired. One application of this representation with strong potential is the prediction of charge density in proteins. This would likely be especially valuable on protein surfaces, where side chains are much less likely to be fixed in single rotamers. Prediction of ensembles of structures provides a method to determine the probability density of charge and atom location; such a prediction is demonstrated graphically.

Amino Acids↗
PubMed2004

Monte Carlo simulation of properties of monolayers and metal islands adsorbed on metallic (111) surfaces.

To obtain the surface stress changes deltasigma due to the adsorption of metal monolayers onto metallic surfaces, a new model derived from thermodynamic considerations is presented. Such a model is based on continuum Monte Carlo simulations with embedded atom method potentials in the canonical ensemble, and it is extended to consider the behavior on different islands adsorbed onto (111) substrate surfaces. Homoepitaxial and heteroepitaxial systems are studied. Pseudomorphic growth is not observed for small metal islands with considerable positive misfit with the substrate. Instead, the islands become compressed upon increase of their size. A simple model is proposed to interpolate between the misfits of atoms in small islands and the pseudomorphic behavior of the monolayer.

Journal Article↗