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 955 records · Page 53Linked to original sources

Direct observation of amyloid growth monitored by total internal reflection fluorescence microscopy.

Most morphological investigations of amyloid fibrils have been performed with various microscopic methods. Among them, direct observation of fibril growth is possible using atomic force microscopy and fluorescence microscopy. Direct observation provides information about the rate and direction of growth at the single fibril level, which cannot be obtained from averaged ensemble measurements. In this chapter, we describe a new technique for the direct observation of amyloid fibril growth using total internal reflection fluorescence microscopy (TIRFM) combined with amyloid-specific thioflavin T (ThT) fluorescence. TIRFM has been developed to monitor single molecules by effectively reducing the background fluorescence in an evanescent field. One of the advantages of TIRFM is that one can selectively monitor fibrils lying along a glass slide, so that one can obtain the exact length of fibrils. This method was used to follow the kinetics of seed-dependent fibril growth of amyloid beta (1-40). The fibril growth was a highly cooperative process, with the fibril ends extending at a constant rate. Because ThT binding is common to all amyloid fibrils, the present method will have general applicability to the real-time analysis of amyloid fibrils.

Amyloid↗

A new method for optimum dose distribution determination taking tumour mobility into account.

A method for determining the optimum dose distribution in the planning target volume is proposed when target volumes are deliberately enlarged to account for tumour mobility in external beam radiotherapy. The optimum dose distribution is a dose distribution that will result in an acceptable level of tumour control probability (TCP) in most of the arising cases of tumour dislocation. An assumption is made that the possible shifts of the tumour are subject to a Gaussian distribution with mean zero and known variance. The idea of a reduced (mean in ensemble) tumour cell density is introduced. On this basis, the target volume and dose distribution in it are determined. The tumour control probability as a function of the shift of the tumour has been calculated. The Monte Carlo method has been used to simulate TCP distributions corresponding to tumour mobility characterized by different variances. The obtained TCP distributions are independent of the variance of the mobility because the dose distribution in the planning target volume is prescribed so that the mobility variance is taken into account. For simplicity a one-dimensional model is used but three-dimensional generalization can be done.

Cell Survival↗

Ortho-Babinet polarization-interrogating filter: an interferometric approach to polarization measurement.

A novel, interferometric, polarization-interrogating filter assembly and method for the simultaneous measurement of all four Stokes parameters across a partially polarized irradiance image in a no-moving-parts, instantaneous, highly sensitive manner is described. In the reported embodiment of the filter, two spatially varying linear retarders and a linear polarizer comprise an ortho-Babinet, polarization-interrogating (OBPI) filter. The OBPI filter uniquely encodes the incident ensemble of electromagnetic wave fronts comprising a partially polarized irradiance image in a controlled, deterministic, spatially varying manner to map the complete state of polarization across the image to local variations in a superposed interference pattern. Experimental interferograms are reported along with a numerical simulation of the method.

Journal Article↗

Improved prediction of MHC class I and class II epitopes using a novel Gibbs sampling approach.

MOTIVATION: Prediction of which peptides will bind a specific major histocompatibility complex (MHC) constitutes an important step in identifying potential T-cell epitopes suitable as vaccine candidates. MHC class II binding peptides have a broad length distribution complicating such predictions. Thus, identifying the correct alignment is a crucial part of identifying the core of an MHC class II binding motif. In this context, we wish to describe a novel Gibbs motif sampler method ideally suited for recognizing such weak sequence motifs. The method is based on the Gibbs sampling method, and it incorporates novel features optimized for the task of recognizing the binding motif of MHC classes I and II. The method locates the binding motif in a set of sequences and characterizes the motif in terms of a weight-matrix. Subsequently, the weight-matrix can be applied to identifying effectively potential MHC binding peptides and to guiding the process of rational vaccine design. RESULTS: We apply the motif sampler method to the complex problem of MHC class II binding. The input to the method is amino acid peptide sequences extracted from the public databases of SYFPEITHI and MHCPEP and known to bind to the MHC class II complex HLA-DR4(B1*0401). Prior identification of information-rich (anchor) positions in the binding motif is shown to improve the predictive performance of the Gibbs sampler. Similarly, a consensus solution obtained from an ensemble average over suboptimal solutions is shown to outperform the use of a single optimal solution. In a large-scale benchmark calculation, the performance is quantified using relative operating characteristics curve (ROC) plots and we make a detailed comparison of the performance with that of both the TEPITOPE method and a weight-matrix derived using the conventional alignment algorithm of ClustalW. The calculation demonstrates that the predictive performance of the Gibbs sampler is higher than that of ClustalW and in most cases also higher than that of the TEPITOPE method.

Algorithms↗

A change in the apparent m value reveals a populated intermediate under equilibrium conditions in Escherichia coli ribonuclease HI.

Experimental studies of protein stability often rely on the determination of an "m value", which describes the denaturant dependence of the free energy change between two states (DeltaG = DeltaG(H2O) - m[denaturant]). Changes in the m value accompanying site specific mutations are usually attributed to structural alterations in the native or unfolded ensemble. Here, we provide an example of significant reduction in the m value resulting from a subtle deviation in two-state behavior not detected by traditional methods. The protein that is studied is a variant of Escherchia coli RNase H in which three residues predicted to be involved in a partially buried salt bridge network were mutated to alanine (R46A, D102A, and D148A). Equilibrium denaturant profiles monitored by both fluorescence and circular dichroism appeared to be cooperative, and a two-state analysis yielded a DeltaG(UN) of approximately -3 kcal/mol with an m value of 1.4 kcal mol(-1) M(-1) (vs 2.3 for RNase H). Analysis of kinetic refolding experiments suggests that the system is actually three-state at equilibrium with an appreciable concentration of an intermediate state under low denaturant concentrations. The stability of the native state determined from a fit of these kinetic data is -6.7 kcal/mol, suggesting that the stability determined by traditional two-state equilibrium analysis is a gross underestimate. The only hint to this loss of two-state behavior was a decrease in the apparent m value, and the presence of the equilibrium intermediate was only identified by a kinetic analysis. Our work serves as a cautionary note; the possibility of a three-state system should be closely addressed before interpreting a change in the m value as a change in the native or unfolded state.

Arginine↗

Structural and thermodynamic features of spiroiminodihydantoin damaged DNA duplexes.

Oxidation of guanine or 8-oxo-7,8-dihydroguanine can produce spiroiminodihydantoin (Sp) R and S stereoisomers. Both in vitro and in vivo experiments have shown that the Sp stereoisomers are highly mutagenic, causing G --> C and G --> T transversion mutations. Therefore, they are of interest as potential endogenous cancer causing lesions. However, their structural properties in DNA duplexes remain to be elucidated. We have employed computational methods to study the Sp lesions in 11-mer DNA duplexes with A, C, G, and T partners. Molecular dynamics simulations have been carried out to obtain ensembles of structures, and the trajectories were employed to analyze the structures and compute free energies. The structural and thermodynamic analyses reveal that the Sp stereoisomers energetically favor positioning in the B-DNA major groove, with minor groove conformers also low energy in some cases, depending on the partner base. The R and S stereoisomers adopt opposite orientations with respect to the 5' to 3' direction of the modified strand. Both syn and anti glycosidic bond conformations are energetically feasible, with partner base and stereochemistry determining the preference. The lesions adversely impact base stacking and Watson-Crick hydrogen bonding interactions in the duplex, and cause groove widening. The chemical nature of the partner base determines specific hydrogen bonding and stacking properties of the damaged duplexes. The structural characteristics may relate to observed mutagenic properties of the Sp stereoisomers, including possible stereoisomer-dependent differences.

Computational Biology↗

Pharmacological activity and membrane interactions of antiarrhythmics: 4D-QSAR/QSPR analysis.

PURPOSE: This study was done to explore the relationships of both macroscopic and molecular level physicochemical properties to in-vivo antiarrhythmic activity and interactions with phospholipid membranes for a set of cationic-amphiphilic analogs. METHODS: The 4D-QSAR method, recently developed by Hopfinger and co-workers (1), was employed to establish 3D-QSAR/QSPR models. Molecular dynamics simulations provided the set of conformational ensembles which were analyzed using partial least squares regression in combination with the Genetic Function Approximation algorithm to construct QSAR and QSPR models. RESULTS: Significant QSAR models for in-vivo antiarrhythmic activity were constructed in which logP (the partition coefficient), and specific grid cell occupancy (spatial) descriptors are the main activity correlates. LogP is the most significant QSAR descriptor. 4D-QSPR models were also developed for two analog-membrane interaction properties, the change in a membrane transition temperature and the ability of the analogs to displace adsorbed Ca(2+)-ions from phosphatidylserine monolayers. CONCLUSIONS: Spatial features, represented by grid cell occupancy descriptors, supplement partition coefficient, which is the most important determinant of in-vivo antiarrhythmic activity, to provide a comprehensive model for drug action. The QSPR models are less significant in statistical measures, and limited to interpretation of possible molecular mechanisms of action.

Animals↗

Epidural electrical stimulation of posterior structures of the human lumbosacral cord: 2. quantitative analysis by computer modeling.

OBJECTIVES: Analysis of the computed recruitment order of an ensemble of ventral and dorsal root fibers should enlighten the relation between the position of a bipolar electrode and the observed order of muscle twitches. MATERIAL AND METHODS: Thresholds of selected spinal root fibers are investigated in a two step procedure. First the electric field generated by the electrodes is computed with the Finite Element Method. In the second step the calculated voltage profile along each target neuron is used as input data for a cable model. For every electrode position the electrical excitability is analyzed for 12 large diameter ventral and dorsal root fibers of the second and fourth lumbar and first sacral segment. The predictions of the neural responses of any target fiber are based on the activating function concept and on the more accurate computer simulations of the electrical behavior of all nodes and internodes in the vicinity of the electrode. RESULTS: For epidural dorsal lumbosacral spinal cord stimulation we found the following rules. (i) The recruitment order of the spinal roots is highly related to the cathode level. (ii) Dorsal root fibers have the lowest threshold values, ventral root fibers are more difficult to excite and dorsal columns are not excitable within the clinical range of 10 V. (iii) For a cathode close to the level of the spinal cord entry of a target fiber thresholds are lowest and spike initiation is expected at the border between cerebrospinal fluid and white matter; excitation of L4 roots is not possible with 210 micros/10 V pulses when cathode is more than 2.2 cm cranial to their entry level (1.5 cm for S1 roots; standard data). (iv) Cathodes positioned (essentially) below the entry level cause spike initiation close to the cathode, in a region where the fibers follow the descending course within the cerebospinal fluid. (v) At rather low stimulation voltage twitches are expected in all investigated lower limb muscles for cathodes below L5 spinal cord level. CONCLUSIONS: Our simulations demonstrate a strong relation between electrode position and the order of muscle twitches which is based on the segmental arrangement of innervation of lower limb muscles. The proposed strategy allows the identification of the position of the electrode relative to spinal cord segments.

Electric Stimulation↗

Validation of protein-unfolding transition states identified in molecular dynamics simulations.

Experimental and simulation studies can complement each other nicely in the area of protein folding. Experiment reports on the average properties of a large ensemble (approx. 10(17)-10(19) molecules), typically over time. Molecular dynamics simulations, on the other hand, provide detailed information for a single molecule, a component of the ensemble. By combining these approaches we can obtain not only a more complete picture of folding, but we can also take advantage of the strengths of different methods. For example, experiment cannot provide molecular structures. Molecular dynamics simulations can provide such information, but the simulations are meaningless without a linked experiment. Thus, the interrelated nature of simulation in assessing experimental assumptions and in providing structures to augment energetic descriptions, and experiment in judging whether the simulations are reasonable, provides more confidence in the resulting information about folding. This combination yields tested and testable molecular models of states that evade characterization by conventional methods. Therefore, we have explored the combined use of these methods to map folding/unfolding pathways at atomic resolution, in collaboration with Alan Fersht. Here we focus on chymotrypsin inhibitor 2, a small single-domain, two-state folding protein.

In Vitro Techniques↗

Structural studies of a peptide activator of human lecithin-cholesterol acyltransferase.

The synthetic lipid-associating peptide, LAP-20 (VSSLLSSLKEYWSSLKESFS), activates lecithin-cholesterol acyltransferase (LCAT) despite its lack of sequence homology to apolipoprotein A-I, the primary in vivo activator of LCAT. Using SDS and dodecylphosphocholine (DPC) to model the lipoprotein environment, the structural features responsible for LAP-20's ability to activate LCAT were studied by optical and two-dimensional 1H NMR spectroscopy. A large blue shift in the intrinsic fluorescence of LAP-20 with the addition of detergent suggested that the peptide formed a complex with the micelles. Analysis of the CD data shows that LAP-20 lacks well defined structure in aqueous solution but adopts helical, ordered conformations upon the addition of SDS or DPC. The helical nature of the peptides in the presence of both lipids was confirmed by upfield H alpha NMR secondary shifts relative to random coil values. Average structures for both peptides in aqueous solutions containing SDS and DPC were generated using distance geometry methods from 329 (SDS) and 309 (DPC) nuclear Overhauser effect-based distance restraints. The backbone (N, Calpha, C=O) RMSD from the average structure of an ensemble of 17 out of 20 calculated structures was 0.41 +/- 0.15 Angstrom for LAP-20 in SDS and 0.41 +/- 0.12 A for an ensemble of 20 out of 20 calculated structures for LAP-20 in DPC. In the presence of SDS, the distance geometry and simulated annealing calculations show that LAP-20 adopts a well defined class A amphipathic helix with distinct hydrophobic and hydrophilic faces. A similar structure was obtained for LAP-20 in the presence of DPC, suggesting that both detergents may be used interchangeably to model the lipoprotein environment. Conformational features of the calculated structures for LAP-20 are discussed relative to models for apolipoprotein A-I activation of LCAT.

Amino Acid Sequence↗

Structure and polymorphism of the principal neutralization site of Thailand HIV-1 isolate.

Refining the geometric parameters for the ensemble of conformers, derived earlier in terms of NMR-spectroscopy data for the immunogenic tip of Thailand HIV-1 isolate, was carried out by quantum chemical methods. As a result, (i) the energy characteristics of initial structures were significantly improved, (ii) their relative locations on the scale of formation heats were determined, and (iii) the energy barriers between conformers under study were computed. On the basis of all data obtained, the high resolution 3D structure model, describing the set of stable conformers and containing the biologically active conformation, was proposed for neutralizing epitope of Thailand HIV-1 isolate. The following major conclusions were made based on the analysis of simulated conformations: i) the Gly-Pro-Gly-Gln-Val-Phe stretch forming the immunogenic crown of Thailand HIV-1 isolate exhibits the properties characteristic for metastable oligopeptide that constitutes in solution the dominant structure with other conformations admissible; (ii) three structures out of five NMR-based starting models form the cluster of conformers which adequately describes general conformational features of this functionally important site of gp120; (iii) two structures residing in this cluster are found to be well-ground for implementing the function of immunoreactive conformation of the stretch of interest; (iv) in spite of this observation, the "global" structure which gives rise to inverse gamma-turn in the central Gly-Pro-Gly crest of Thailand HIV-1 gp120 is proposed to be the most probable conformation responsible for the formation of viral antigen-antibody complex in particular case under study.

Cystine↗

HCOP: a searchable database of human orthology predictions.

The HUGO Gene Nomenclature Committee (HGNC) Comparison of Orthology Predictions (HCOP) search tool combines the human, mouse, rat and chicken orthology assertions made by PhIGs, HomoloGene, Ensembl, Inparanoid, Mouse Genome Informatics (MGI) and HGNC, enabling users to identify predicted ortholog pairs for a specified gene or genes. The HCOP resource provides a useful method to integrate, compare and access a variety of disparate sources of human orthology data. The HCOP search tool, data and documentation are available at http://www.gene.ucl.ac.uk/hcop.

Databases, Genetic↗

Integrative Network Analysis of Bioactive Compounds from Punica granatum L. Peel: Multi-Target Mechanisms in Wound Healing.

BACKGROUND: Wound-healing agents often have limited efficacy and require prolonged recovery times, prompting growing interest in developing herbal-based formulations. Among these, Punica granatum L. has attracted considerable attention because of its high polyphenolic content. Despite its widespread use, the precise pharmacological targets underlying its wound-healing effects remain poorly understood and require systematic investigation. OBJECTIVES: This study aimed to elucidate the underlying pharmacological mechanisms of the topical wound-healing properties of P. granatum L. using a network pharmacology approach. METHODS: Bioactive compounds of P. granatum L. and their potential target genes were identified using the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP), Similarity Ensemble Approach (SEA), and SwissTargetPrediction databases. Wound healing-related genes were retrieved from the GeneCards database. Genes intersecting P. granatum L. targets and wound healing-associated genes were subjected to functional enrichment analyses, including protein-protein interaction (PPI), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. The PPI network was further analyzed using Cytoscape, and the phytoconstituent-target interaction network was visualized using Gephi. These findings were validated using molecular docking. RESULTS: A total of 40 intersecting genes were identified as potential P. granatum L. targets involved in wound healing. Among these, EGFR, PTPN11, HRAS, IGF1R, and ESR1 were identified as key hub genes. Functional enrichment analysis indicated that the most significantly enriched signaling pathways included the MAPK, PI3K-Akt, EGFR tyrosine kinase inhibitor resistance, focal adhesion, and FoxO signaling pathways. Molecular docking analysis confirmed favorable binding of quercetin and ellagic acid to the hub targets EGFR, IGF1R, and ESR1. CONCLUSIONS: These findings elucidate the pharmacological pathways underlying P. granatum-mediated wound healing and suggest that P. granatum L. acts as a multi-target modulator in the wound-healing process.

Focal Adhesion↗

Premature ventricular contractions during +Gz with and without pressure breathing and extended coverage anti-G suit.

BACKGROUND: High +Gz is known to provoke cardiac dysrhythmias. Pressure breathing during G (PBG) and extended coverage anti-G suits (ECGS) are used to enhance +Gz-endurance and reduce fatigue during high +Gz flying. It is not known whether PBG in combination with ECGS increases the risk for premature ventricular contractions (PVC). HYPOTHESIS: PBG in combination with ECGS increases the risk of PVCs during high +Gz-loads. METHODS: Retrospective data were obtained from 14 subjects exposed to three different simulated aerial combat sorties each, in the centrifuge with a standard anti-G suit ensemble or with PBG in combination with ECGS. Each sortie consisted of a gradual onset G-exposure (GOR) and three simulated aerial combat maneuvers (SACM) containing four cycles of a +4.5 to +7 GzSACM; four cycles of +4 to +9 Gz tactical aerial combat maneuver (TACM) with several rapid transitions from +4 or +5 Gz to +8 or +9 Gz; and four cycles of +5 to +9 Gz SACM (5-9 SACM) with four cycles. ECG was recorded during the +Gz exposures to reveal any cardiac dysrhythmias. RESULTS AND CONCLUSIONS: No PVCs occurred during the GORs. During 4.5-7 SACMs, TACMs, and 5-9 SACMs there were 83, 50, and 24 PVCs with standard equipment, respectively, and 63, 54, and 39 PVCs with PBG and ECGS, respectively. There was no statistically significant difference between the equipment in any of the different types of +Gz exposures. No episodes of supraventricular tachycardia or relative bradycardia were found with either equipment.

Adult↗

Effect of G suit type on cognitive performance.

BACKGROUND: Sustained acceleration protection ensembles are being developed to help pilots of high performance aircraft endure high G exposures for longer periods of time. It has been assumed that better G endurance confers better pilot task performance. This premise was studied on the Armstrong Laboratory Dynamic Environment Simulator centrifuge. METHOD: Human subjects repeatedly endured prolonged high-G simulated aerial combat on a centrifuge to the point of loss of vision or physical exhaustion. Some profiles included over 20 exposures to +9 Gz. While enduring the G exposures, subjects tracked a simulated "bogey" aircraft on a visual display and performed a secondary task. Measures of cognitive function and physiologic status were taken throughout the exposures. G protection ensembles included the standard CSU 13 B/P anti-G suit, the Advanced Technology Anti-G Suit (ATAGS), COMBAT EDGE positive pressure breathing system with the CSU 13 B/P, COMBAT EDGE with ATAGS, and the Northrop Advanced Protection System (APS). RESULTS: More advanced protective systems not only allow longer G endurance, but provide adequate support for maintained cognitive performance throughout the extended exposure. Although measures were affected by the type of protective system the subject was wearing, as well as individual ability and coping strategies, consistent target tracking task performance, rapid choice reaction time, and sufficient arterial oxygen saturation were maintained throughout extended exposures to a point preceding termination by only a second or two. CONCLUSIONS: Those anti-G protection ensembles that cover and protect the body and which employ positive pressure breathing allow longer high G exposures which provide support for maintained cognitive performance.

Adaptation, Psychological↗

Image interpolation by two-dimensional parametric cubic convolution.

Cubic convolution is a popular method for image interpolation. Traditionally, the piecewise-cubic kernel has been derived in one dimension with one parameter and applied to two-dimensional (2-D) images in a separable fashion. However, images typically are statistically nonseparable, which motivates this investigation of nonseparable cubic convolution. This paper derives two new nonseparable, 2-D cubic-convolution kernels. The first kernel, with three parameters (designated 2D-3PCC), is the most general 2-D, piecewise-cubic interpolator defined on [-2, 2] x [-2, 2] with constraints for biaxial symmetry, diagonal (or 90 degrees rotational) symmetry, continuity, and smoothness. The second kernel, with five parameters (designated 2D-5PCC), relaxes the constraint of diagonal symmetry, based on the observation that many images have rotationally asymmetric statistical properties. This paper also develops a closed-form solution for determining the optimal parameter values for parametric cubic-convolution kernels with respect to ensembles of scenes characterized by autocorrelation (or power spectrum). This solution establishes a practical foundation for adaptive interpolation based on local autocorrelation estimates. Quantitative fidelity analyses and visual experiments indicate that these new methods can outperform several popular interpolation methods. An analysis of the error budgets for reconstruction error associated with blurring and aliasing illustrates that the methods improve interpolation fidelity for images with aliased components. For images with little or no aliasing, the methods yield results similar to other popular methods. Both 2D-3PCC and 2D-5PCC are low-order polynomials with small spatial support and so are easy to implement and efficient to apply.

Algorithms↗

Weighted density-functional theory for simple fluids: prewetting of a Lennard-Jones fluid.

The prewetting of a Lennard-Jones fluid is studied using weighted density-functional theory. The intrinsic Helmholtz free-energy functional is separated into repulsive and attractive contributions. An accurate functional for hard spheres is used for the repulsive functional and a weighted density-functional method is used for the attractive part. The results for this theory are compared against mean-field density-functional theory, the theory of Velasco and Tarazona [E. Velasco and P. Tarazona, J. Chem. Phys. 91, 7916 (1989)] and grand canonical ensemble simulation results. The results demonstrate that the weighted density functional for attractive forces may offer a significant increase in accuracy over the other theories. The density-functional and simulation results also indicate that a previous estimate of the wetting temperature for a model of the interaction of argon with solid carbon dioxide, obtained from simulations [J. E. Finn and P. A. Monson, Phys. Rev. A, 39, 6402 (1989)], is incorrect. The weighted density-functional method indicates that triple-point prewetting is observed for this model potential.

Journal Article↗

Designing generalized statistical ensembles for numerical simulations of biopolymers.

Conformational properties of polymers, such as average dihedral angles or molecular alpha-helicity, display a rather weak dependence on the detailed arrangement of the elementary constituents (atoms). We propose a computer simulation method to explore the polymer phase space using a variant of the standard multicanonical method, in which the density of states associated to suitably chosen configurational variables is considered in place of the standard energy density of states. This configurational density of states is used in the Metropolis acceptance/rejection test when configurations are generated with the help of a hybrid Monte Carlo algorithm. The resulting configurational probability distribution is then modulated by exponential factors derived from the general principle of the maximal constrained entropy by requiring that certain average configurational quantities take preassigned (possibly temperature dependent) values. Thermal averages of other configurational quantities can be computed by using the probability distributions obtained in this way. Moments of the energy distribution require an extra canonical sampling of the system phase space at the desired temperature, in order to locally thermalize the configurational degrees of freedom. As an application of these ideas we present the study of the structural properties of two simple models: a bead-and-spring model of polyethylene with independent hindered torsions and an all-atom model of alanine and glycine oligomers with 12 amino acids in vacuum.

Biopolymers↗