PubMed Health⌕ Search

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 397 records · Page 22Linked to original sources

Dynamic structure of a highly ordered beta-sheet molten globule: multiple conformations with a stable core.

The structure of [14-38]Abu, a variant of bovine pancreatic trypsin inhibitor (BPTI) with only the 14-38 disulfide bridge intact, has been analyzed by two-dimensional 1H and 1H-15N NMR. Except for the 18-24, 29-35 antiparallel beta-sheet, residues in all regions of the molecule give two exchange cross peaks for each 1H; for one residue, Gly 37, three exchange cross peaks are observed. The presence of exchange cross peaks indicates that the residues sample conformations that interconvert on a time scale of milliseconds or longer. Over 90% of the NMR spectra have been assigned, including backbone and side chain atoms and their exchange cross peaks. Analyses of chemical shifts, chemical exchange, hydrogen isotope exchange, and NOEs indicate that [14-38]Abu at pH 4.5 and 1 degree C is an ensemble of interconverting conformations, in all of which the 18-24, 29-35 antiparallel beta-sheet is native-like and intact. Outside the antiparallel beta-sheet, residues undergo local order/disorder transitions. The stable structure of [14-38]Abu is not in the vicinity of the 14-38 disulfide bond but rather is in the slow-exchange core. NOE analysis indicates that the main tertiary interactions involve hydrophobic contacts with the rings of Tyr 21, Tyr 23, and Tyr 35. As a model for early folding intermediates, the structure of [14-38]Abu suggests that BPTI folding is initiated by stabilization of a turn existing in the unfolded protein and involves both local and nonlocal hydrophobic interactions.

Amino Acid Sequence↗

Simulations of protein folding and unfolding.

The 'new view' of protein folding is based on a statistical analysis of the landscape for protein folding. This perspective leads to the investigation of the statistical distributions of protein conformations as the folding protein approaches the native state from unfolded states. Molecular dynamics simulations of both the thermodynamics of protein folding and the kinetics of unfolding are beginning to explore the statistical nature of these distributions. They also provide connections between the theory of protein folding landscapes and the experimental observations of the properties of key ensembles of the conformations populated as folding progresses.

Models, Chemical↗

Composites of local structure propensities: evidence for local encoding of long-range structure.

To estimate how extensively the ensemble of denatured-state conformations is constrained by local side-chain-backbone interactions, propensities of each of the 20 amino acids to occur in mono- and dipeptides mapped to discrete regions of the Ramachandran map are computed from proteins of known structure. In addition, propensities are computed for the trans, gauche-, and gauche+ rotamers, with or without consideration of the values of phi and psi. These propensities are used in scoring functions for fragment threading, which estimates the energetic favorability of fragments of protein sequence to adopt the native conformation as opposed to hundreds of thousands of incorrect conformations. As finer subdivisions of the Ramachandran plot, neighboring residue phi/psi angles, and rotamers are incorporated, scoring functions become better at ranking the native conformation as the most favorable. With the best composite propensity function, the native structure can be distinguished from 300,000 incorrect structures for 71% of the 2130 arbitrary protein segments of length 40, 48% of 2247 segments of length 30, and 20% of 2368 segments of length 20. A majority of fragments of length 30-40 are estimated to be folded into the native conformation a substantial fraction of the time. These data suggest that the variations observed in amino acid frequencies in different phi/psi/chi1 environments in folded proteins reflect energetically important local side-chain-backbone interactions, interactions that may severely restrict the ensemble of conformations populated in the denatured state to a relatively small subset with nativelike structure.

Amino Acids↗

Reassessing random-coil statistics in unfolded proteins.

The Gaussian-distributed random coil has been the dominant model for denatured proteins since the 1950s, and it has long been interpreted to mean that proteins are featureless, statistical coils in 6 M guanidinium chloride. Here, we demonstrate that random-coil statistics are not a unique signature of featureless polymers. The random-coil model does predict the experimentally determined coil dimensions of denatured proteins successfully. Yet, other equally convincing experiments have shown that denatured proteins are biased toward specific conformations, in apparent conflict with the random-coil model. We seek to resolve this paradox by introducing a contrived counterexample in which largely native protein ensembles nevertheless exhibit random-coil characteristics. Specifically, proteins of known structure were used to generate disordered conformers by varying backbone torsion angles at random for approximately 8% of the residues; the remaining approximately 92% of the residues remained fixed in their native conformation. Ensembles of these disordered structures were generated for 33 proteins by using a torsion-angle Monte Carlo algorithm with hard-sphere sterics; bulk statistics were then calculated for each ensemble. Despite this extreme degree of imposed internal structure, these ensembles have end-to-end distances and mean radii of gyration that agree well with random-coil expectations in all but two cases.

Animals↗

Solution conformations of the gamma-carboxyglutamic acid domain of bovine prothrombin fragment 1, residues 1-65.

Molecular dynamics simulations have been performed (AMBER version 3.1) on solvated residues 1-65 of bovine prothrombin fragment 1 (BF1) by using the 2.8-A resolution crystallographic coordinates as the starting conformation for understanding calcium ion-induced conformational changes that precede experimentally observable phospholipid binding. Simulations were performed on the non-metal-bound crystal structure, the form resulting from addition of eight calcium ions to the 1-65 region of the crystal structure, the form resulting from removal of calcium ions after 107 ps and continuing the simulation, and an isolated hexapeptide loop (residues 18-23). In all cases, the 100-ps time scale seemed adequate to sample an ensemble of solution conformers within a particular region of conformation space. The non-metal-containing BF1 did not unfold appreciably during a 106-ps simulation starting from the crystallographic geometry. The calcium ion-containing structure (Ca-BF1) underwent an interesting conformational reorganization during its evolution from the crystal structure: during the time course of a 107-ps simulation, Ca-BF1 experienced a trans----cis isomerization of the gamma-carboxyglutamic acid-21 (Gla-21)-Pro-22 peptide bond. Removal of the calcium ions from this structure followed by 114 ps of additional molecular dynamics showed significant unfolding relative to the final 20-ps average structure of the 107-ps simulation; however, the Gla-21-Pro-22 peptide bond remained cis. A 265-ps simulation on the termini-protected hexapeptide loop (Cys-18 to Cys-23) containing two calcium ions also did not undergo a trans----cis isomerization. It is believed that the necessary activation energy for the transitional event observed in the Ca-BF1 simulation was largely supplied by global conformational events with a possible assist from relief of intermolecular crystal packing forces. The presence of a Gla preceding Pro-22, the inclusion of Pro-22 in a highly strained loop structure, and the formation of two long-lived salt bridges prior to isomerization may all contribute to this finding.

1-Carboxyglutamic Acid↗

Information content of molecular structures.

For a completely enumerated set of conformers of a macromolecule or for exhaustive lattice walks of model polymers it is straightforward to use Shannon information theory to deduce the information content of the ensemble. It is also practicable to develop numerical measures of the information content of sets of exact distance constraints applied to specific conformational ensembles. We examine the effects of experimental uncertainties by considering "noisy" constraints. The introduction of noise requires additional assumptions about noise distribution and conformational clustering protocols that make the problem of measuring information content more complex. We make use of a standard concept in communication theory, the "noise sphere," to link uncertainty in measurements to information loss. Most of our numerical results are derived from two-dimensional lattice ensembles. Expressing results in terms of information per degree of freedom removes almost all of the chain length dependence. We also explore off-lattice polyalanine chains that yield surprisingly similar results.

Amino Acid Sequence↗

Recognition properties of V3-specific antibodies to V3 loop peptides derived from HIV-1 gp120 presented in multiple conformations.

To identify structural constraints and amino acid sequences important for antibody recognition of the third variable domain (V3) of HIV-1 gp120, we have studied the solution conformation of three 35-mer circular V3 loop peptides derived from HIV-1 strains which differ in syncytium- (SI) and non-syncytium-inducing (NSI) capacity. In addition to 2D NMR and CD analyses, fluid- and solid-phase immunoassays were performed using V3-specific antibodies to V3 peptides and gp120 derived from different strains of HIV-1. NMR and CD spectroscopy indicated that circular and linear V3 loops exist in water as a dynamic ensemble of multiple conformations. Amino acid substitutions and biochemical modifications of the V3 loop were found to affect antibody binding depending on the presentation of the antigens. From NMR observations and immunological experiments, we provide evidence for a V3 loop specific monoclonal antibody interaction which is directed predominantly against linear epitopes rather than against discontinuous epitopes. The absence of a single defined solution conformation of 35-mer circular V3 peptides should be taken into account when using V3-related peptides to investigate structural elements in the V3 domain of the gp120 envelope protein of HIV-1 involved in biological processes of the virus.

Amino Acid Sequence↗

Localization of the N-terminal domain in light-harvesting chlorophyll a/b protein by EPR measurements.

The conformational distribution of the N-terminal domain of the major light-harvesting chlorophyll a/b protein (LHCIIb) has been characterized by electron-electron double resonance yielding distances between spin labels placed in various domains of the protein. Distance distributions involving residue 3 near the N terminus turned out to be bimodal, revealing that this domain, which is involved in regulatory functions such as balancing the energy flow through photosystems (PS) I and II, exists in at least two conformational states. Models of the conformational sub-ensembles were generated on the basis of experimental distance restraints from measurements on LHCIIb monomers and then checked for consistency with the experimental distance distribution between residues 3 in trimers. Only models where residue 3 is located above the core of the protein and extends into the aqueous phase on the stromal side fit the trimer data. In the other state, which consequently is populated only in monomers, the N-terminal domain extends sideways from the protein core. The two conformational states may correspond to two functional states of LHCIIb, namely trimeric LHCIIb associated with PSII in stacked thylakoid membranes and presumably monomeric LHCIIb associated with PSI in nonstacked thylakoids. The switch between these two is known to be triggered by phosphorylation of Thr-6. A similar phosphorylation-induced conformational change of the N-terminal domain has been observed by others in bovine annexin IV which, due to the conformational switch, also loses its membrane-aggregating property.

Animals↗

Thermodynamic mechanism and consequences of the polyproline II (PII) structural bias in the denatured states of proteins.

A quantitative characterization of the structure and energy of the denatured states of proteins represents the cornerstone to a molecular-level understanding of both protein stability and fold specificity. Recent studies have revealed a significant bias in unstructured peptides toward the polyproline II (P(II)) conformation, even when no prolines are present in the sequence. This indicates that the P(II) conformation is a dominant component of the denatured states of proteins, although a quantitative description of the component enthalpy and entropy functions associated with this conformation (i.e., the thermodynamic mechanism) has thus far proven elusive. An experimental system has been designed that, when analyzed with high-precision isothermal titration calorimetry, provides direct access to the residue-specific thermodynamics of the P(II) structure formation in disordered proteins and peptides. Here, it is shown that the P(II) bias is driven by a favorable and significant enthalpy (Deltah) of -1.7 kcal mol(-1) residue(-1), which is partially offset by an unfavorable entropy (TDeltas) of -0.7 kcal mol(-1) residue(-1), relative to the ensemble of disordered conformations of the molecule. In addition to impacting dramatically the interpretation of thermal denaturation experiments, these experimental values form the framework of a quantitative energetic description of the denatured states of proteins.

Animals↗

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins.

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.

Alanine↗

Unblocked statistical-coil tetrapeptides in aqueous solution: quantum-chemical computation of the carbon-13 NMR chemical shifts.

We recently reported a theoretical characterization of representative ensembles of statistical-coil conformations for tetrapeptides with unblocked termini in aqueous solution, at pH 7. The results showed good agreement between the computed Boltzmann-averaged and experimentally-determined values for both the vicinal coupling constants (3)J(NHalpha) and the alpha-proton chemical shifts. Here, we carry out a cluster analysis of the ensembles of conformations generated in that study, and use them to compute the Boltzmann-averaged values of the quantum-chemical (13)C chemical shifts for different amino acids in the unblocked tetrapeptides GGXA (where X stands for Phe, Arg, His, Glu, Ile, Lys, Gln, Tyr, Leu, Thr, Ala, Gly and Val). The values of the (13)C chemical shifts in these thirteen amino acids (for which experimental data are available) were computed by using Density Functional Theory with a 6-311+G(2d,p) basis set. Good agreement is found in terms of both the correlation coefficient (R) and standard deviations of the difference between the computed Bolztmann-averaged and the NMR-determined values for the (13)C chemical shifts. These results suggest that it may be possible to build a reliable theoretically-derived database of chemical shifts for statistical-coil residues. The results of the current study contribute to our understanding of the relations between chemical shifts, dihedral angles and vicinal coupling constants, (3)J(NHalpha). In addition, they can shed light as to how the statistical-coil conformation is related to the conformational preference of more structured states, such as the alpha-helical conformation.

Amino Acid Sequence↗

Test for cooperativity in the early kinetic intermediate in lysozyme folding.

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.

Animals↗

An empirical energy potential with a reference state for protein fold and sequence recognition.

We consider modifications of an empirical energy potential for fold and sequence recognition to represent approximately the stabilities of proteins in various environments. A potential used here includes a secondary structure potential representing short-range interactions for secondary structures of proteins, and a tertiary structure potential consisting of a long-range, pairwise contact potential and a repulsive packing potential. This potential is devised to evaluate together the total conformational energy of a protein at the coarse grained residue level. It was previously estimated from the observed frequencies of secondary structures, from contact frequencies between residues, and from the distributions of the number of residues in contact in known protein structures by regarding those distributions as the equilibrium distributions with the Boltzmann factor of these interaction energies. The stability of native structures is assumed as a primary requirement for proteins to fold into their native structures. A collapse energy is subtracted from the contact energies to remove the protein size dependence and to represent protein stabilities for monomeric and multimeric states. The free energy of the whole ensemble of protein conformations that is subtracted from the conformational energy to represent protein stability is approximated as the average energy expected for a typical native structure with the same amino acid composition. This term may be constant in fold recognition but essentially varies in sequence recognition. A simple test of threading sequences into structures without gaps is employed to demonstrate the importance of the present modifications that permit the same potential to be utilized for both fold and sequence recognition. Proteins 1999;36:357-369. Published 1999 Wiley-Liss, Inc.

Models, Chemical↗

Are there temperature-dependent structural transitions in the "intrinsically unstructured" protein prothymosin alpha?

Prothymosin alpha, a typical member of the class of the so-called "intrinsically unstructured" proteins, adopts a random-chain conformation under physiological environmental conditions. An apparent formation of ordered secondary structure and a moderate compaction are observed upon the change from neutral to acid pH at room temperature. We have addressed the question of whether there are temperature-dependent changes of the conformational state of prothymosin alpha at low pH using circular dichroism spectroscopy and static and dynamic light scattering. In contrast to previous investigations, we did not observe a heat-induced conformational transition. For comparison, we have also carried out the same experimental procedures with acid-unfolded phosphoglycerate kinase from yeast. In this case we observed a weak compaction and a slight apparent increase in ordered secondary structure with increasing temperature, probably caused by the higher average hydrophobicity as compared to prothymosin alpha. In the absence of a clear structural transition, we deduce the observed effects result mainly from a progressive redistribution in the population of phi-psi angles of the polypeptide backbone when the temperature is increased. Furthermore, the paper should demonstrate the difficulties in distinguishing between such a progressive change amongst a continuum of states within the ensemble of unfolded conformations from the formation of authentic stable secondary structures in highly unfolded proteins. This problem is not solved presently and convincing evidence can only be supplied by the combination of various experimental techniques.

Circular Dichroism↗

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.

Dimethyl Sulfoxide↗

Coupling between conformation and proton binding in proteins.

Interest centers here on whether the use of a fixed charge distribution of a protein solute, or a treatment that considers proton-binding equilibria by solving the Poisson equation, is a better approach to discriminate native from non-native conformations of proteins. In this analysis of the charge distribution of 7 proteins, we estimate the solvation free energy contribution to the total free energy by exploring the 2(zeta) possible ionization states of the whole molecule, with zeta being the number of ionizable groups in the amino acid sequence, for every conformation in the ensembles of 7 proteins. As an additional consideration of the role of electrostatic interactions in determining the charge distribution of native folds, we carried out a comparison of alternative charge assignment models for the ionizable residues in a set of 21 native-like proteins. The results of this work indicate that (1) for 6 out of 7 proteins, estimation of solvent polarization based on the Generalized Born model with a fixed charge distribution provides the optimal trade-off between accuracy, with respect to the Poisson equation, and speed when compared to the accessible surface area model; for the seventh protein, consideration of all possible ionization states of the whole molecule appears to be crucial to discriminate the native from non-native conformations; (2) significant differences in the degree of ionization and hence the charge distribution for native folds are found between the different charge models examined; (3) the stability of the native state is determined by a delicate balance of all the energy components, and (4) conformational entropy, and hence the dynamics of folding, may play a crucial role for a successful ab initio protein folding prediction.

Ions↗

A kinetic and conformational study on the interaction of tetrahydropteridines with tyrosine hydroxylase.

Tetrahydropterins are obligatory cofactors for tyrosine hydroxylase (TH), the rate-limiting enzyme of catecholamine biosynthesis. A series of synthetic analogues of 6(R)-L-erythro-5,6,7, 8-tetrahydrobiopterin (BH(4)) with different substituents in positions C2, N3, C4, N5, C6, C7, and N8 on the ring were used as active site probes for recombinant human TH. The enzyme tolerates rather bulky substituents at C6, as seen by the catalytic efficiency (V(max)/K(m)) and the coupling efficiency (mol of L-DOPA produced/mol of tetrahydropterin oxidized) of the cofactors. Substitutions at C2, C4, N5, and N8 abolish the cofactor activity of the pterin analogues. Molecular docking of BH(4) into the crystal structure of the catalytic domain of ligand-free rat TH results in complexes in which the pteridine ring pi-stacks with Phe300 and the N3 and the amino group at C2 hydrogen bonds with Glu332. The pteridine ring also establishes interactions with Leu294 and Gln310. The distance between C4a in the pteridines and the active site iron was 4.2 +/- 0.5 A for the ensemble of docked conformers. Docking of BH(4) analogues into TH also shows that the most bulky substituents at C6 can be well-accommodated within the large hydrophobic pocket surrounded by Ala297, Ser368, Tyr371, and Trp372, without altering the positioning of the ring. The pterin ring of 7-BH(4) shows proper stacking with Phe300, but the distance between the C4a and the active site iron is 0.6 A longer than for bound BH(4), a finding that may be related to the high degree of uncoupling observed for 7-BH(4).

Animals↗

A decoy set for the thermostable subdomain from chicken villin headpiece, comparison of different free energy estimators.

BACKGROUND: Estimators of free energies are routinely used to judge the quality of protein structural models. As these estimators still present inaccuracies, they are frequently evaluated by discriminating native or native-like conformations from large ensembles of so-called decoy structures. RESULTS: A decoy set is obtained from snapshots taken from 5 long (100 ns) molecular dynamics (MD) simulations of the thermostable subdomain from chicken villin headpiece. An evaluation of the energy of the decoys is given using: i) a residue based contact potential supplemented by a term for the quality of dihedral angles; ii) a recently introduced combination of four statistical scoring functions for model quality estimation (FRST); iii) molecular mechanics with solvation energy estimated either according to the generalized Born surface area (GBSA) or iv) the Poisson-Boltzmann surface area (PBSA) method. CONCLUSION: The decoy set presented here has the following features which make it attractive for testing energy scoring functions:1) it covers a broad range of RMSD values (from less than 2.0 A to more than 12 A);2) it has been obtained from molecular dynamics trajectories, starting from different non-native-like conformations which have diverse behaviour, with secondary structure elements correctly or incorrectly formed, and in one case folding to a native-like structure. This allows not only for scoring of static structures, but also for studying, using free energy estimators, the kinetics of folding;3) all structures have been obtained from accurate MD simulations in explicit solvent and after molecular mechanics (MM) energy minimization using an implicit solvent method. The quality of the covalent structure therefore does not suffer from steric or covalent problems. The statistical and physical effective energy functions tested on the set behave differently when native simulation snapshots are included or not in the set and when averaging over the trajectory is performed.

Animals↗