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 1,027 records · Page 57Linked to original sources
PubMed2003

Coding of auditory space.

Behavioral, anatomical, and physiological approaches can be integrated in the study of sound localization in barn owls. Space representation in owls provides a useful example for discussion of place and ensemble coding. Selectivity for space is broad and ambiguous in low-order neurons. Parallel pathways for binaural cues and for different frequency bands converge on high-order space-specific neurons, which encode space more precisely. An ensemble of broadly tuned place-coding neurons may converge on a single high-order neuron to create an improved labeled line. Thus, the two coding schemes are not alternate methods. Owls can localize sounds by using either the isomorphic map of auditory space in the midbrain or forebrain neural networks in which space is not mapped.

Animals↗
PubMed2007

The use of high-pressure nuclear magnetic resonance to study protein folding.

Recent development of high-pressure cells for a variety of spectroscopic methods has enabled the use of pressure as one of the commonly used perturbations along with temperature and chemical perturbations to study folding/unfolding reactions of proteins. Although various high-pressure spectroscopy techniques have their own significance, high-pressure nuclear magnetic resonance (NMR) is unique in that it allows one to gain residue-specific and atom-detailed information from proteins under pressure. Furthermore, because of a peculiar volume property of a protein, high-pressure NMR allows one to obtain structural information of a protein in a wide conformational space from the bottom to the upper region of the folding funnel, giving structural reality for the "open" state of a protein proposed from hydrogen exchange. The method allows a link between equilibrium folding intermediates and the kinetic intermediates, and manifests a new view of proteins as dynamic entities amply fluctuating among the folded, intermediate, and unfolded sub ensembles. This chapter briefly summarizes the technique, the principle, and the ways to use high-pressure NMR for studying protein folding.

Models, Molecular↗
PubMed2000

Interaction of chemical defense clothing and high terrestrial altitudes on lift/carry and marksmanship performance.

BACKGROUND: The increased metabolic energy requirement imposed by a chemical defense uniform (CDU) and the lower maximal aerobic capacity associated with increased altitude should produce greater demands on the cardiopulmonary system during the performance of a given work task at increasing altitudes. We hypothesized that: a) relative to sea level, the decrements in physical work performance caused by ascending to high terrestrial altitudes would be greater in a CDU compared with a standard fatigue uniform (U.S. Army, BDU); b) the aversive subjective reactions to the CDU would be accentuated with increasing altitude; and c) that staging at moderate altitude, to induce acclimatization, would restore work performance at higher altitudes to sea level norms. METHODS: The physiological and subjective responses of 8 male soldiers to work (10-min lift-and-carry task and rifle marksmanship) were measured. Subjects wore the BDU and a CDU ensemble (U.S. Army, BDO) at sea level, intermediate (2743 m) and high (4,300 m) altitudes following rapid and staged (3 d at 1,830 m) ascents to the higher altitudes. RESULTS: Lift/carry task performance tended to be lower (p = 0.076) in the CDU vs. the BDU at altitude. The cardiopulmonary responses to the lift/carry task increased at altitude and were greater in the CDU. The subjects' perception of their ability to perform the lift/carry task at altitude was adversely impacted more in the CDU than the BDU. Rapid ascent to intermediate altitude degraded marksmanship in both uniforms. Following staged ascent, lift/carry task and marksmanship performance was restored to sea level norms. CONCLUSIONS: Personnel wearing CDU or equivalent protective clothing at intermediate to high terrestrial elevations should anticipate proportionally larger CDU-induced decrements of work performance and lower tolerance to working in a CDU than experienced near sea level. Staging at moderate altitude is an effective strategy for restoring work performance to sea level norms at higher altitudes.

Activities of Daily Living↗
PubMed1990

Secondary structure of substance P bound to liposomes in organic solvents and in solution from Raman and CD spectroscopy.

The membrane lipid phase may be an important mediator of the peptide-receptor interaction. In order to understand the mechanism of this interaction, it is important to know the peptide structure, not only in the hydrophobic lipid bilayer environment, but also at the bilayer surface and in solution. To investigate this problem we have measured the secondary structure of the 11-residue neuropeptide substance P (SP) and its fragments in aqueous solutions, in membrane mimetic solvents, and associated with lipid bilayers using Raman and CD spectroscopy. Raman and CD spectra of SP bound to liposomes indicate a less than 20% helix content. We interpret these results to indicate that SP contains virtually no helix when bound to negatively charged liposomes. These spectra are similar to spectra of peptides in type I and III beta-turns. SP forms between 10 and 30% (1-3 residues) helical structure in sodium dodecyl sulfate micelles and less than 10% helix in methanol and trifluoroethanol. The binding of SP to negatively charged liposomes significantly changes the structure of the lipid acyl chains, decreasing order in some cases and increasing it in others. Raman spectra of SP in water indicates that SP near 30 mM forms an ensemble of structures in water that is distinct from completely unfolded peptide and from the aggregated beta-sheet form observed in saline solutions. We conclude from our CD results that methods used to quantitate secondary structure from CD spectra of short peptides cannot be used to distinguish between very short helical segments and beta-turns.

Amino Acid Sequence↗
PubMed2025

In silico prediction method for plant Nucleotide-binding leucine-rich repeat- and pathogen effector interactions.

Plant Nucleotide-binding leucine-rich repeat (NLR) proteins play a crucial role in effector recognition and activation of Effector triggered immunity following pathogen infection. Genome sequencing advancements have led to the identification of a myriad of NLRs in numerous agriculturally important plant species. However, deciphering which NLRs recognize specific pathogen effectors remains challenging. Predicting NLR-effector interactions in silico will provide a more targeted approach for experimental validation, critical for elucidating function, and advancing our understanding of NLR-triggered immunity. In this study, NLR-effector protein complex structures were predicted using AlphaFold2-Multimer for all experimentally validated NLR-effector interactions reported in literature. Binding affinities- and energies were predicted using 97 machine learning models from Area-Affinity. We show that AlphaFold2-Multimer predicted structures have acceptable accuracy and can be used to investigate NLR-effector interactions in silico. Binding affinities for 58 NLR-effector complexes ranged between -8.5 and -10.6 log(K), and binding energies between -11.8 and -14.4 kcal/mol-1, depending on the Area-Affinity model used. For 2427 "forced" NLR-effector complexes, these estimates showed larger variability, enabling identification of novel NLR-effector interactions with 99% accuracy using an Ensemble machine learning model. The narrow range of binding energies- and affinities for "true" interactions suggest a specific change in Gibbs free energy, and thus conformational change, is required for NLR activation. This is the first study to provide a method for predicting NLR-effector interactions, applicable to all pathosystems. Finally, the NLR-Effector Interaction Classification (NEIC) resource can streamline research efforts by identifying NLRs important for plant-pathogen resistance, advancing our understanding of plant immunity.

Plant Proteins↗
PubMed2003

Determination of molecular weight distributions of diblock copolymers with conventional size exclusion chromatography.

A convenient method that enables determination of molecular weight distributions of diblock copolymers with conventional size exclusion chromatography (SEC) without referring to any kind of calibration curve or molar mass-sensitive detecting devices is described. It is based on the well-established procedure to calculate the number-average molecular weight of copolymers from their chemical composition and the molar mass of the first block segment (precursor), applied on an ensemble of monodisperse copolymer fractions. The chemical composition of SEC fractions is calculated from their UV and RI detector traces, and the number-average molecular weight of the precursor, which is to be known from independent measurements, is taken as the molar mass of the first block segment. The molecular weight distributions and averages obtained by this method for selected copolymer samples were found to be in reasonably good agreement with the ones determined by NMR and SEC with on-line viscosity or multiangle laser light scattering detection.

Letter↗
PubMed1990

Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins.

We present a prototype of a new approach to the folding problem of polypeptide chains. This approach is based on the analysis of known protein structures. It derives the energy potentials for the atomic interactions of all amino acid residue pairs as a function of the distance between the involved atoms. These potentials are then used to calculate the energies of all conformations that exist in the data base with respect to a given sequence. Then, by using only the most stable conformations, clusters of the most probable conformations for the given sequence are obtained. To discuss the results properly we introduce a new classification of segments based on their conformational stability. Special care is taken to allow for sparse data sets. The use of the method is demonstrated in the discussion of the identical oligopeptide sequences found in different conformations in unrelated proteins. VNTFV, for example, adopts a beta-strand in ribonuclease but it is found in an alpha-helical conformation in erythrocruorin. In the case of VNTFV the ensemble obtained consists of a single cluster of beta-strand conformations, indicating that this may be the preferred conformation for the pentapeptide. When the flanking residues are included in the calculation the hepapeptide P-VNTFV-H (ribonuclease) again yields an ensemble of beta-strands. However, in the ensemble of D-VNTFV-A (erythrocruorin) the major cluster is of alpha-helical type. In the present study we concentrate on the local aspects of protein conformations. However, the theory presented is quite general and not restricted to oligopeptides. We indicate extensions of the approach to the calculation of global conformations of proteins as well as conceivable applications to a number of molecular systems.

Amino Acid Sequence↗
PubMed2001

Comparative application of multiple receptor methods to identify aerosol sources in northern Vermont.

This study applies and compares results of four receptor modeling techniques to a common set of speciated fine particle measurement data collected ata remote site in northwestern Vermont between 1988 and 1995. Two multivariate mathematical models, positive matrix factorization and UNMIX, were applied to the measurement data and identified seven "common" sources that had similar compositions and similar fine mass contributions in both models. Two ensemble backward trajectory techniques, potential source contribution function and residence-time analysis, were also applied to evaluate and interpret the mathematical model results. The trajectory techniques indicate a strong regional character to the upwind locations associated with aerosol contributions from most of the sources identified independently by the mathematical models and help in the interpretation of those results. The process of model comparison provides insights on the strengths and limitations of the individual and combined source attribution techniques. Convergent results among the multiple methods provide a degree of confidence that each of the receptor methods may represent useful tools for future air quality management. Divergent or inconsistent results among the models can help identify limitations of the individual models and of the underlying aerosol and meteorological data sets.

Aerosols↗
PubMed2003

Simple, intuitive calculations of free energy of binding for protein-ligand complexes. 2. Computational titration and pH effects in molecular models of neuraminidase-inhibitor complexes.

One factor that can strongly influence predicted free energy of binding is the ionization state of functional groups on the ligands and at the binding site at which calculations are performed. This analysis is seldom performed except in very detailed computational simulations. In this work, we address the issues of (i) modeling the complexity resulting from the different ionization states of ligand and protein residues involved in binding, (ii) if, and how, computational methods can evaluate the pH dependence of ligand inhibition constants, and (iii) how to score the protonation-dependent models. We developed a new and fairly rapid protocol called "computational titration" that enables parallel modeling of multiple ionization ensembles for each distinct protonation level. Models for possible protonation combinations for site/ligand ionizable groups are built, and the free energy of interaction for each of them is quantified by the HINT (Hydropathic INTeractions) software. We applied this procedure to the evaluation of the binding affinity of nine inhibitors (six derived from 2,3-didehydro-2-deoxy-N-acetylneuraminic acid, DANA) of influenza virus neuraminidase (NA), a surface glycoprotein essential for virus replication and thus a pharmaceutically relevant target for the design of anti-influenza drugs. The three-dimensional structures of the NA enzyme-inhibitor complexes indicate considerable complexity as the ligand-protein recognition site contains several ionizable moieties. Each computational titration experiment reveals a peak HINT score as a function of added protons. This maximum HINT score indicates the optimum pH (or the optimum protonation state of each inhibitor-protein binding site) for binding. The pH at which inhibition is measured and/or crystals were grown and analyzed can vary from this optimum. A protonation model is proposed for each ligand that reconciles the experimental complex structure with measured inhibition and the free energy of binding. Computational titration methods allow us to analyze the effect of pH in silico and may be helpful in improving ligand binding free energy prediction when protonation or deprotonation of the residues or ligand functional groups at the binding site might be significant.

Algorithms↗
PubMed2003

In silico models for the prediction of dose-dependent human hepatotoxicity.

The liver is extremely vulnerable to the effects of xenobiotics due to its critical role in metabolism. Drug-induced hepatotoxicity may involve any number of different liver injuries, some of which lead to organ failure and, ultimately, patient death. Understandably, liver toxicity is one of the most important dose-limiting considerations in the drug development cycle, yet there remains a serious shortage of methods to predict hepatotoxicity from chemical structure. We discuss our latest findings in this area and present a new, fully general in silico model which is able to predict the occurrence of dose-dependent human hepatotoxicity with greater than 80% accuracy. Utilizing an ensemble recursive partitioning approach, the model classifies compounds as toxic or non-toxic and provides a confidence level to indicate which predictions are most likely to be correct. Only 2D structural information is required and predictions can be made quite rapidly, so this approach is entirely appropriate for data mining applications and for profiling large synthetic and/or virtual libraries.

Computer Simulation↗
PubMed2006

Conformer generation under restraints.

Conformational sampling by direct optimization of an all-atom energy function is ineffective and inefficient because of the ruggedness of the energy landscape. Discrete sampling schemes represent an attractive alternative for generating ensembles of conformers consistent with spatial restraints derived from empirical data. Conformational sampling is becoming increasingly important for structure prediction as the bottleneck in accurate prediction shifts from energy functions to the methods used to find low-energy conformers. Experimental structure determination remains a perennial challenge as investigators tackle larger macromolecular systems, and begin to incorporate more complete descriptions of uncertainty, heterogeneity and dynamics into their models. Computational approaches that combine dense, discrete sampling with all-atom energy evaluation and refinement may help to overcome the remaining barriers to solving these problems.

Algorithms↗
PubMed1997

Pattern generation.

Central pattern generators are neuronal ensembles capable of producing the basic spatiotemporal patterns underlying 'automatic' movements (e.g. locomotion, respiration, swallowing and defense reactions), in the absence of peripheral feedback. Different experimental approaches, from classical electrophysiological and pharmacological methods to molecular and genetic ones, have been used to understand the cellular and synaptic bases of central pattern generator organization and reconfiguration of generator operation in behaviorally relevant contexts. Recently, it has been shown that the high reliability and flexibility of central pattern generators is determined by their redundant organization. Everything that is crucial for generator operation is determined by a number of complementary mechanisms acting in concert; however, various mechanisms are weighted differently in determining different aspects of central pattern generator operation.

Animals↗
PubMedDate not supplied

Identification of a minimal subset of receptor conformations for improved multiple conformation docking and two-step scoring.

Docking and scoring are critical issues in virtual drug screening methods. Fast and reliable methods are required for the prediction of binding affinity especially when applied to a large library of compounds. The implementation of receptor flexibility and refinement of scoring functions for this purpose are extremely challenging in terms of computational speed. Here we propose a knowledge-based multiple-conformation docking method that efficiently accommodates receptor flexibility thus permitting reliable virtual screening of large compound libraries. Starting with a small number of active compounds, a preliminary docking operation is conducted on a large ensemble of receptor conformations to select the minimal subset of receptor conformations that provides a strong correlation between the experimental binding affinity (e.g., Ki, IC50) and the docking score. Only this subset is used for subsequent multiple-conformation docking of the entire data set of library (test) compounds. In conjunction with the multiple-conformation docking procedure, a two-step scoring scheme is employed by which the optimal scoring geometries obtained from the multiple-conformation docking are re-scored by a molecular mechanics energy function including desolvation terms. To demonstrate the feasibility of this approach, we applied this integrated approach to the estrogen receptor alpha (ERalpha) system for which published binding affinity data were available for a series of structurally diverse chemicals. The statistical correlation between docking scores and experimental values was significantly improved from those of single-conformation dockings. This approach led to substantial enrichment of the virtual screening conducted on mixtures of active and inactive ERalpha compounds.

Protein Binding↗
PubMed1997

The dynamic NMR structure of the T psi C-loop: implications for the specificity of tRNA methylation.

tRNA (m5U54)-methyltransferase (RUMT) catalyzes the S-adenosylmethionine-dependent methylation of uridine-54 in the T psi C-loop of all transfer RNAs in E. coli to form the 54-ribosylthymine residue. However, in all tRNA structures, residue 54 is completely buried and the question arises as to how RUMT gains access to the methylation site. A 17-mer RNA hairpin consisting of nucleotides 49-65 of the T psi-loop is a substrate for RUMT. Homonuclear NMR methods in conjunction with restrained molecular dynamics (MD) methods were used to determine the solution structure of the 17-mer T-arm fragment. The loop of the hairpin exhibits enhanced flexibility which renders the conventional NMR average structure less useful compared to the more commonly found situation where a molecule exists in predominantly one major conformation. However, when resorting to softer refinement methods such as MD with time-averaged restraints, the conflicting restraints in the loop can be satisfied much better. The dynamic structure of the T-arm is represented as an ensemble of 10 time-clusters. In all of these, U54 is completely exposed. The flexibility of the T psi-loop in solution in conjunction with extensive binding studies of RUMT with the T psi C-loop and tRNA suggest that the specificity of the RUMT/ tRNA recognition is associated with tRNA tertiary structure elements. For the methylation, RUMT would simply have to break the tertiary interactions between the D- and T-loops, leading to a melting of the T-arm structure and making U54 available for methylation.

Base Sequence↗
PubMed2004

Optimal laser control of ultrafast photodissociation of I2- in water: mixed quantum/classical molecular dynamics simulation.

A linearized optimal control method in combination with mixed quantum/classical molecular dynamics simulation is used for numerically investigating the possibility of controlling photodissociation wave packets of I(2)(-) in water. Optimal pulses are designed using an ensemble of photodissociation samples, aiming at the creation of localized dissociation wave packets. Numerical results clearly show the effectiveness of the control although the control achievement is reduced with an increase in the internuclear distance associated with a target region. We introduce effective optimal pulses that are designed using a statistically averaged effective dissociation potential, and show that they semiquantitatively reproduce the control achievements calculated by using optimal pulses. The control mechanisms are interpreted from the time- and frequency-resolved spectra of the effective optimal pulses.

Journal Article↗
PubMed2005

Constant-pressure and constant-surface tension simulations in dissipative particle dynamics.

We present a method for constant-pressure and constant-surface tension simulations in dissipative particle dynamics using a Langevin piston approach. We demonstrate that the corresponding equations of motion lead to the relevant ensembles and propose an appropriate scheme of integration. After having identified a suitable set of parameters for the approach, we demonstrate the feasibility of the approach by applying it to two different systems, a simple isotropic fluid and an anisotropic fluid lipid-bilayer membrane in water. Results are presented for, respectively, isothermal bulk compressibility, tracer diffusion coefficient, lipid head-group area, and isothermal area compressibility. We find that our Langevin piston approach leads to improvements over other approaches in terms of faster equilibration and shorter correlation times of various system variables.

Computer Simulation↗
PubMed2005

The reorganization energy of electron transfer in nonpolar solvents: molecular level treatment of the solvent.

The intermolecular electron transfer in a solute pair consisting of pyrene and dimethylaniline is investigated in a nonpolar solvent, n-hexane. The earlier elaborated approach [M. Tachiya, J. Phys Chem. 97, 5911 (1993)] is used; this method provides a physically relevant background for separating inertial and inertialess polarization responses for both nonpolarizable and polarizable molecular level simulations. The molecular-dynamics technique was implemented for obtaining the equilibrium ensemble of solvent configurations. The nonpolar solvent, n-hexane, was treated in terms of OPLS-AA parametrization. Solute Lennard-Jones parameters were taken from the same parametrization. Solute charge distributions of the initial and final states were determined using ab initio level [HF/6-31G(d,p)] quantum-chemical calculations. Configuration analysis was performed explicitly taking into account the anisotropic polarizability of n-hexane. It is shown that the Gaussian law well describes calculated distribution functions of the solvent coordinate, therefore, the rate constant of the ET reaction can be characterized by the reorganization energy. Evaluated values of the reorganization energies are in a range of 0.03-0.11 eV and significant contribution (more then 40% of magnitude) comes from anisotropic polarizability. Investigation of the reorganization energy lambda dependence on the solute pair separation distance d revealed unexpected behavior. The dependence has a very sharp peak at the distance d=7 A where solvent molecules are able to penetrate into the intermediate space between the solute pair. The reason for such behavior is clarified. This new effect has a purely molecular origin and cannot be described within conventional continuum solvent models.

Journal Article↗
PubMed1992

Markov modelling of ensemble current relaxations: bovine adrenal nicotinic receptor currents analysed.

1. A general approach to the analysis of ensemble currents of ligand-gated channels is presented, with a variety of examples that include single and multiple agonist binding steps, desensitization and several blocking pathways. 2. The use of matrix methods to describe model reaction schemes leads to a simplification if the reaction scheme is irreversible: the product of the exponential relaxation rate constants is exactly equal to the product of the forward reaction steps. 3. This method of analysis applied to the bovine adrenal nicotinic receptor suggests that in the range of acetylcholine concentrations from 1 microM to 2 mM, a model with a single kinetically relevant agonist binding step is appropriate. 4. Complex models, to explain the presence of two desensitizing components in currents recorded from whole cells, may be discounted in favour of two distinct receptor types. 5. A simple model of open channel block is discounted, and desensitization of the blocked state proposed.

Acetylcholine↗