SEARCH · PubMed Health

Results for “conformational ensemble”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 379 records · Page 21Linked to original sources

PubMed2005

Stitchprofiles.uio.no: analysis of partly melted DNA conformations using stitch profiles.

In this study, we describe a web server that performs computations on DNA melting, thus predicting the localized separation of the two strands for sequences provided by the users. The output types are stitch profiles, melting curves, probability profiles, etc. Stitch profile diagrams visualize the ensemble of alternative conformations that DNA can adopt with different probabilities. For example, a stitch profile shows the possible loop openings in terms of their locations, sizes, probabilities and fluctuations at a given temperature. Sequences with lengths up to several tens or hundreds of kilobase pairs can be analysed. The tools are freely available at http://stitchprofiles.uio.no.

Algorithms↗

PubMed1995

De novo prediction of polypeptide conformations using dihedral probability grid Monte Carlo methodology.

We tested the dihedral probability grid Monte Carlo (DPG-MC) methodology to determine optimal conformations of polypeptides by applying it to predict the low energy ensemble for two peptides whose solution NMR structures are known: integrin receptor peptide (YGRGDSP, Type II beta-turn) and S3 alpha-helical peptide (YMSEDEL KAAEAAFKRHGPT). DPG-MC involves importance sampling, local random stepping in the vicinity of a current local minima, and Metropolis sampling criteria for acceptance or rejection of new structures. Internal coordinate values are based on side-chain-specific dihedral angle probability distributions (from analysis of high-resolution protein crystal structures). Important features of DPG-MC are: (1) Each DPG-MC step selects the torsion angles (phi, psi, chi) from a discrete grid that are then applied directly to the structure. The torsion angle increments can be taken as S = 60, 30, 15, 10, or 5 degrees, depending on the application. (2) DPG-MC utilizes a temperature-dependent probability function (P) in conjunction with Metropolis sampling to accept or reject new structures. For each peptide, we found close agreement with the known structure for the low energy conformational ensemble located with DPG-MC. This suggests that DPG-MC will be useful for predicting conformations of other polypeptides.

Amino Acid Sequence↗

PubMed2025

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations.

MOTIVATION: Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data. RESULTS: We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization. AVAILABILITY AND IMPLEMENTATION: It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Quantitative Trait Loci↗

PubMed2004

Sparsely populated folding intermediates of the Fyn SH3 domain: matching native-centric essential dynamics and experiment.

A complete description of how a protein folds requires the characterization of intermediate conformations traversed during the folding transition. We have calculated dynamics trajectories of a simplified model of the Fyn SH3 domain with a native-centric potential energy function. Analysis of the resulting site-resolved energy trajectory identifies an ensemble of intermediate conformations for folding and another for unfolding. The model's folding intermediate is structured in the three beta-strands that make up the protein's core and is strikingly similar to intermediates detected in a recent NMR study of Fyn SH3 folding and to folding transition states elucidated in mutagenesis studies of SH3 domains. The unfolding intermediate is formed by dissociation of the folded protein's two terminal beta-strands from its core. The presence of such an intermediate is consistent with the results of a protein-engineering study on the src SH3 domain showing that these strands separate before the rate-limiting step of unfolding. Despite the presence of these conformations intermediate between the native and fully unfolded states, the computed heat capacity vs. temperature profile of the model protein indicates that its thermodynamics satisfies the usual calorimetric criterion for two-state folding. This observation highlights the fact that, if not properly interpreted, methods such as calorimetry that do not probe multiple sites in a molecule can lead to an oversimplified view of folding. The close agreement between results from this simplified model and experimental work underscores the important contributions that computational methods can make in providing insights into protein folding.

Adaptor Proteins, Signal Transducing↗

PubMed1999

The folding transition state between SH3 domains is conformationally restricted and evolutionarily conserved.

The protein engineering analysis of the alpha-spectrin SH3 domain at three different stability conditions (pH 7.0, 3.5 and 2.5) reveals a folding transition state structured around the distal loop beta-hairpin and the 310-helix. This region is impervious to overall changes in protein stability, suggesting a transition state ensemble with little conformational variability. Comparison with the Src SH3 domain (36% sequence homology) indicates that the transition state in this protein family may be conserved. Discrepancies at some positions can be rationalized in terms of the different interactions made by the different side chains in both domains. Brønsted plot analysis confirms the straight phi(doubledagger-U) results and shows two folding subdomains for this small protein. These results, together with previous data on circular permutants of the alpha-spectrin SH3 domain, indicate that polypeptide topology and chain connectivity play a major role in the folding reaction of this protein family.

Amino Acid Sequence↗

PubMed2006

Improvement of statistical potentials and threading score functions using information maximization.

We show that statistical potentials and threading score functions, derived from finite data sets, are informatic functions, and that their performance depends on the manner in which data are classified and compressed. The choice of sequence and structural parameters affects estimates of the conditional probabilities P(C|S), the quantification of the effect of sequence S on conformation C, and determines the amount of information extracted from the data set, as measured by information gain. The mathematical link between information gain and mean conformational energy, established in this work using the local backbone potential as model, demonstrates that manipulation of descriptive parameters also alters the "energy" values assigned to native conformation and to decoy structures in the test pool, and consequently, the performance of such statistical potential functions in fold recognition exercises. We show that sequence and structural partitions that maximize information gain also minimize the mean energy of the ensemble of native conformations. Moreover, we establish an informatic basis for the placement of the native score within an energy spectrum given by the decoy pool in a threading exercise. We discover that, among all informatic quantities, information gain is the best predictor of threading success, even better than the standard Z-score. Consequently, the choices of sequence and structural descriptors, extent of compression, and levels of discretization that maximize information gain must also produce the best potential functions. Strategies to optimize these parameters with respect to information extraction are therefore relevant to building better statistical potentials. Last, we demonstrate that the backbone torsion potential, defined by the trimer sequence, can be an effective tool in greatly reducing the set of possible conformations from a vast decoy pool.

Amino Acid Sequence↗

PubMed1994

Modeling compact denatured states of proteins.

We propose a model for the conformations of compact denatured states of globular proteins: that they are broad ensembles of chain backbone conformations that involve common localized hydrophobic clustering and helical contacts, depending on the amino acid sequence. We construct representative ensembles for chain lengths up to 136 monomers on three-dimensional cubic lattices using the "hydrophobic zippers" method (Fiebig & Dill, 1993). We find that model conformations with radii of gyration about 20% larger than native conformations commonly have bimodal distributions of P(r), of the pairwise interatomic distances, r, and Kratky plots in agreement with recent small-angle X-ray scattering (Sosnick & Trewhella, 1992; Flanagan et al., 1992; Kataoka et al., 1993; Flanagan et al., 1993) experiments on three different proteins. We also find that the lattice model of the Shortle 1-136 fragment of staphylococcal nuclease does not appear capable of forming a single hydrophobic core by hydrophobic zippering, consistent with experiments.

Amino Acid Sequence↗

PubMed1998

Conformation of the Ras-binding domain of Raf studied by molecular dynamics and free energy simulations.

Recognition of Ras by its downstream target Raf is mediated by a Ras-recognition region in the Ras-binding domain (RBD) of Raf. Residues 78-89 in this region occupy two different conformations in the ensemble of NMR solution structures of the RBD: a fully alpha-helical one, and one where 87-90 form a type IV beta-turn. Molecular dynamics simulations of the RBD in solution were performed to explore the stability of these and other possible conformations of both the wild-type RBD and the R89K mutant, which does not bind Ras. The simulations sample a fully helical conformation for residues 78-89 similar to the NMR helical structures, a conformation where 85-89 form a 3(10)-helical turn, and a conformation where 87-90 form a type I beta-turn, whose free energies are all within 0.3 kcal/mol of each other. NOE patterns and H(alpha) chemical shifts from the simulations are in reasonable agreement with experiment. The NMR turn structure is calculated to be 3 kcal/mol higher than the three above conformations. In a simulation with the same implicit solvent model used in the NMR structure generation, the turn conformation relaxes into the fully helical conformation, illustrating possible structural artifacts introduced by the implicit solvent model. With the Raf R89K mutant, simulations sample a fully helical and a turn conformation, the turn being 0.9 kcal/mol more stable. Thus, the mutation affects the population of RBD conformations, and this is expected to affect Ras binding. For example, if the fully helical conformation of residues 78-89 is required for binding, its free energy increase in R89K will increase the binding free energy by about 0.6 kcal/mol.

Algorithms↗

PubMed2000

Solution conformational study of Scyliorhinin I analogues with conformational constraints by two-dimensional NMR and theoretical conformational analysis.

Two analogues of Scyliorhinin I (Scyl), a tachykinin with N-MeLeu in position 8 and a 1,5-disubstituted tetrazole ring between positions 7 and 8, introduced in order to generate local conformational constraints, were synthesized using the solid-phase method. Conformational studies in water and DMSO-d6 were performed on these peptides using a combination of the two-dimensional NMR technique and theoretical conformational analysis. The algorithm of conformational search consisted of the following three stages: (i) extensive global conformational analysis in order to find all low-energy conformations; (ii) calculation of the NOE effects and vicinal coupling constants for each of the low energy conformations; (iii) determining the statistical weights of these conformations by means of a nonlinear least-squares procedure, in order to obtain the best fit of the averaged simulated spectrum to the experimental one. In both solvents the three-dimensional structure of the analogues studied can be interpreted only in terms of an ensemble of multiple conformations. For [MeLeu8]Scyl, the C-terminal 6-10 fragment adopts more rigid structure than the N-terminal one. In the case of the analogue with the tetrazole ring in DMSO-d6 the three-dimensional structure is characterized by two dominant conformers with similar geometry of their backbones. They superimpose especially well (RMSD = 0.28 A) in the 6-9 fragments. All conformers calculated in both solvents superimpose in their C-terminal fragments much better than those of the first analogue. The results obtained indicate that the introduction of the tetrazole ring into the Scyl molecule rigidifies its structure significantly more than that of MeLeu.

Algorithms↗

PubMed2005

Molecular docking of balanol to dynamics snapshots of protein kinase A.

Even if the structure of a receptor has been determined experimentally, it may not be a conformation to which a ligand would bind when induced fit effects are significant. Molecular docking using such a receptor structure may thus fail to recognize a ligand to which the receptor can bind with reasonable affinity. Here, we examine one way to alleviate this problem by using an ensemble of receptor conformations generated from a molecular dynamics simulation for molecular docking. Two molecular dynamics simulations were conducted to generate snapshots for protein kinase A: one with the ligand bound, the other without. The ligand, balanol, was then docked to conformations of the receptors presented by these trajectories. The Lamarckian genetic algorithm in Autodock [Goodsell et al. J Mol Recognit 1996;9(1):1-5; Morris et al. J Comput Chem 1998;19(14):1639-1662] was used in the docking. Three ligand models were used: rigid, flexible, and flexible with torsional potentials. When the snapshots were taken from the molecular dynamics simulation of the protein-ligand complex, the correct docking structure could be recovered easily by the docking algorithm in all cases. This was an easier case for challenging the docking algorithm because, by using the structure of the protein in a protein-ligand complex, one essentially assumed that the protein already had a pocket to which the ligand can fit well. However, when the snapshots were taken from the ligand-free protein simulation, which is more useful for a practical application when the structure of the protein-ligand complex is not known, several clusters of structures were found. Of the 10 docking runs for each snapshot, at least one structure was close to the correctly docked structure when the flexible-ligand models were used. We found that a useful way to identify the correctly docked structure was to locate the structure that appeared most frequently as the lowest energy structure in the docking experiments to different snapshots.

Azepines↗

PubMed2004

Characterization of the conformational space of a triple-stranded beta-sheet forming peptide with molecular dynamics simulations.

Molecular dynamics (MD) simulations have been performed on a series of mutants of the 20 amino acid peptide Betanova in order to critically assess the ability of MD simulations to reproduce the folding and stability of small beta-sheet-forming peptides on currently accessible timescales. Simulations were performed in both water and in 40% methanol solution, using an explicit solvent model. The simulations suggest that all mutants adopt a wide range of conformations in solution, that the structures are highly flexible, and that stabilization of compact structures is due to a delicate balance of hydrophobic and polar side-chain interactions. Simulations longer than 100 ns, although not sufficient for a complete thermodynamic and kinetic description of the system, sample an ensemble of compact conformations characterized by the loss of ordered beta-sheet secondary structure. This suggests that no significant free energy barrier separates the different conformations available.

Algorithms↗

PubMed2002

Phosphorylation-induced conformational changes in a mitogen-activated protein kinase substrate. Implications for tyrosine hydroxylase activation.

Mitogen-activated protein (MAP) kinase-mediated phosphorylation of specific residues in tyrosine hydroxylase leads to an increase in enzyme activity. However, the mechanism whereby phosphorylation affects enzyme turnover is not well understood. We used a combination of fluorescence resonance energy transfer (FRET) measurements and molecular dynamics simulations to explore the conformational free energy landscape of a 10-residue MAP kinase substrate found near the N terminus of the enzyme. This region is believed to be part of an autoregulatory sequence that overlies the active site of the enzyme. FRET was used to measure the effect of phosphorylation on the ensemble of peptide conformations, and molecular dynamics simulations generated free energy profiles for both the unphosphorylated and phosphorylated peptides. We demonstrate how FRET transfer efficiencies can be calculated from molecular dynamics simulations. For both the unphosphorylated and phosphorylated peptides, the calculated FRET efficiencies are in excellent agreement with the experimentally determined values. Moreover, the FRET measurements and molecular simulations suggest that phosphorylation causes the peptide backbone to change direction and fold into a compact structure relative to the unphosphorylated state. These results are consistent with a model of enzyme activation where phosphorylation of the MAP kinase substrate causes the N-terminal region to adopt a compact structure away from the active site. The methods we employ provide a general framework for analyzing the accessible conformational states of peptides and small molecules. Therefore, they are expected to be applicable to a variety of different systems.

Animals↗

PubMed2000

Mutational analysis of the propensity for amyloid formation by a globular protein.

Acylphosphatase can be converted in vitro, by addition of trifluoroethanol (TFE), into amyloid fibrils of the type observed in a range of human diseases. The propensity to form fibrils has been investigated for a series of mutants of acylphosphatase by monitoring the range of TFE concentrations that result in aggregation. We have found that the tendency to aggregate correlates inversely with the conformational stability of the native state of the protein in the different mutants. In accord with this, the most strongly destabilized acylphosphatase variant forms amyloid fibrils in aqueous solution in the absence of TFE. These results show that the aggregation process that leads to amyloid deposition takes place from an ensemble of denatured conformations under conditions in which non-covalent interactions are still favoured. These results support the hypothesis that the stability of the native state of globular proteins is a major factor preventing the in vivo conversion of natural proteins into amyloid fibrils under non-pathological conditions. They also suggest that stabilizing the native states of amyloidogenic proteins could aid prevention of amyloidotic diseases.

Acid Anhydride Hydrolases↗

PubMed1989

Low-temperature photoreactions of halorhodopsin. 1. Detection of conformational substates of the chromoprotein.

Absorption spectra of halorhodopsin (HR), a retinal protein in the halobacterial membrane, and its photostationary states were determined at 80 K. The absorption lines appear to narrow upon cooling, thereby revealing complex spectral fine structure of the main absorption band in the visible region, characteristic of conformational substates of HR. Illumination causes (1) the redistribution of these substrates and consequent changing of the fine structure ("hole-burning") and (2) the appearance of a hypsoproduct of undefined nature, in addition to the previously described bathoproduct HR600. Bacteriorhodopsin, a related retinal pigment, gives rise only to the bathointermediate (i.e., K590) under these conditions. After warming of illuminated HR to 110 K, and recooling to 80 K, relaxation of the illumination-induced change in spectral fine structure, and decay of the hypsoproduct but not the bathoproduct, was observed. The results are explained with a model in which one ensemble of HR conformational substates at 80 K is converted to another in a photoequilibrium via the excited state, which also produces the batho- and hypsoproducts. The original ensemble can be regained through thermal pathways at a somewhat higher temperature, and only the bathoproduct will decay thermally into the next intermediate of the HR photocycle.

Bacteriorhodopsins↗

PubMed2002

Protein-protein docking with multiple residue conformations and residue substitutions.

The protein docking problem has two major aspects: sampling conformations and orientations, and scoring them for fit. To investigate the extent to which the protein docking problem may be attributed to the sampling of ligand side-chain conformations, multiple conformations of multiple residues were calculated for the uncomplexed (unbound) structures of protein ligands. These ligand conformations were docked into both the complexed (bound) and unbound conformations of the cognate receptors, and their energies were evaluated using an atomistic potential function. The following questions were considered: (1) does the ensemble of precalculated ligand conformations contain a structure similar to the bound form of the ligand? (2) Can the large number of conformations that are calculated be efficiently docked into the receptors? (3) Can near-native complexes be distinguished from non-native complexes? Results from seven test systems suggest that the precalculated ensembles do include side-chain conformations similar to those adopted in the experimental complexes. By assuming additivity among the side chains, the ensemble can be docked in less than 12 h on a desktop computer. These multiconformer dockings produce near-native complexes and also non-native complexes. When docked against the bound conformations of the receptors, the near-native complexes of the unbound ligand were always distinguishable from the non-native complexes. When docked against the unbound conformations of the receptors, the near-native dockings could usually, but not always, be distinguished from the non-native complexes. In every case, docking the unbound ligands with flexible side chains led to better energies and a better distinction between near-native and non-native fits. An extension of this algorithm allowed for docking multiple residue substitutions (mutants) in addition to multiple conformations. The rankings of the docked mutant proteins correlated with experimental binding affinities. These results suggest that sampling multiple residue conformations and residue substitutions of the unbound ligand contributes to, but does not fully provide, a solution to the protein docking problem. Conformational sampling allows a classical atomistic scoring function to be used; such a function may contribute to better selectivity between near-native and non-native complexes. Allowing for receptor flexibility may further extend these results.

Amino Acid Substitution↗

PubMed1997

Conformational substates in enzyme mechanism: the 120 K structure of alpha-lytic protease at 1.5 A resolution.

Insight into the dynamic properties of alpha-lytic protease (alpha LP) has been obtained through the use of low-temperature X-ray crystallography and multiple-conformation refinement. Previous studies of alpha LP have shown that the residues around the active site are able to move significantly to accommodate substrates of different sizes. Here we show a link between the ability to accommodate ligands and the dynamics of the binding pocket. Although the structure of alpha LP at 120 K has B-factors with a uniformly low value of 4.8 A2 for the main chain, four regions stand out as having significantly higher B-factors. Because thermal motion should be suppressed at cryogenic temperatures, the high B-factors are interpreted as the result of trapped conformational substates. The active site residues that are perturbed during accommodation of different substrates are precisely those showing conformational substates, implying that substrate binding selects a subset of conformations from the ensemble of accessible states. To better characterize the precise nature of these substates, a protein model consisting of 16 structures has been refined and evaluated. The model reveals a number of features that could not be well-described by conventional B-factors: for example, 40% of the main-chain residue conformations are distributed asymmetrically or in discrete clusters. Furthermore, these data demonstrate an unexpected correlation between motions on either side of the binding pocket that we suggest is a consequence of "dynamic close packing." These results provide strong evidence for the role of protein dynamics in substrate binding and are consistent with the results of dynamic studies of ligand binding in myoglobin and ribonuclease A.

Bacterial Proteins↗

PubMed2006

NMR-derived dynamic aspects of N-type inactivation of a Kv channel suggest a transient interaction with the T1 domain.

Some eukaryotic voltage-gated K+ (Kv) channels contain an N-terminal inactivation peptide (IP), which mediates a fast inactivation process that limits channel function during membrane depolarization and thus shapes the action potential. We obtained sequence-specific nuclear magnetic resonance (NMR) assignments for the polypeptide backbone of a tetrameric N-terminal fragment (amino acids 1-181) of the Aplysia Kv1.1 channel. Additional NMR measurements show that the tetramerization domain 1 (T1) has the same globular structure in solution as previously determined by crystallography and that the IP (residues 1-20) and the linker (residues 21-65) are in a flexibly disordered, predominantly extended conformation. A potential contact site between the T1 domain and the flexible tail (residues 1-65) has been identified on the basis of chemical-shift changes of individual T1 domain amino acids, which map to the T1 surface near the interface between adjacent subunits. Paramagnetic perturbation experiments further indicate that, in the ensemble of solution conformers, there is at least a small population of species with the IP localized in close proximity to the proposed interacting residues of the T1 tetramer. Electrophysiological measurements show that all three mutations in this pocket that we tested slow the rate of inactivation and speed up recovery, as predicted from the preinactivation site model. These results suggest that specific, short-lived transient interactions between the T1 domain and the IP or the linker segment may play a role in defining the regulatory kinetics of fast channel inactivation.

Amino Acid Sequence↗

PubMed2006

Flexible protein-protein docking.

Predicting the structure of protein-protein complexes using docking approaches is a difficult problem whose major challenges include identifying correct solutions, and properly dealing with molecular flexibility and conformational changes. Flexibility can be addressed at several levels: implicitly, by smoothing the protein surfaces or allowing some degree of interpenetration (soft docking) or by performing multiple docking runs from various conformations (cross or ensemble docking); or explicitly, by allowing sidechain and/or backbone flexibility. Although significant improvements have been achieved in the modeling of sidechains, methods for the explicit inclusion of backbone flexibility in docking are still being developed. A few novel approaches have emerged involving collective degrees of motion, multicopy representations and multibody docking, which should allow larger conformational changes to be modeled.

Binding Sites↗