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 613 records · Page 34Linked to original sources

Structural distortion of p53 by the mutation R249S and its rescue by a designed peptide: implications for "mutant conformation".

Missense mutations in the DNA-binding core domain of the tumour suppressor protein p53 are frequent in cancer. Many of them result in loss of native structure. The mutation R249S is one of the six most common cancer-associated p53 mutations ("hot-spots"). As it is highly frequent in hepatocellular carcinoma, its rescue is an important therapeutic target. We have used NMR techniques to study the structural effects of the R249S mutation. The overall fold of the core domain is retained in R249S, and it does not take up a denatured "mutant conformation". However, the beta-sandwich had increased flexibility and, according to changes in chemical shift, there was local distortion throughout the DNA-binding interface. It is likely that the R249S mutation resulted in an ensemble of native and native-like conformations in a dynamic equilibrium. The peptide FL-CDB3 that was designed to rescue mutants of p53 by binding specifically to its native structure was found to revert the chemical shifts of R249S back towards the wild-type values and so reverse the structural effects of mutation. We discuss the implications for a rescue strategy and also for the analysis of antibody-binding data.

Binding Sites↗

Long-lived conformational isomerism of protein dimers: the role of the free energy of subunit association.

The association of protein subunits to form N-mers (N >or= 3) does not follow the dependence on the law of mass action predicted by the classical thermodynamic description used for the equilibrium of association of small molecules. For those anomalous cases, a so-called deterministic model has been previously proposed. The latter model was based on the empirical observation that the dynamics of subunit exchange between protein oligomers can be very slow, leading to the existence of long-lived conformational isomers and to a persistently heterogeneous ensemble of oligomers in solution. Contrary to the expectation for a protein dimer, we have recently shown that the subunit association of triosephosphate isomerase (TIM) could also be described as a deterministic process and that long-lived conformational isomers of TIM could be isolated in solution. Here we show that a), observation of hysteresis in pressure dissociation curves is an additional indicator of deterministic behavior; b), the extent of deviation from the classical thermodynamic behavior correlates with the free-energy change of subunit association; and c), experimental manipulation of the free energy of subunit association through the addition of a subdenaturing concentration of a chaotropic agent restores the concentration dependence of subunit association of TIM. A model that explains these features and its biological relevance is discussed.

Dimerization↗

Transition-state ensemble in enzyme catalysis: possibility, reality, or necessity?

Proteins are not rigid structures; they are dynamic entities, with numerous conformational isomers (substates). The dynamic nature of protein structures amplifies the structural variation of the transition state for chemical reactions performed by proteins. This suggests that utilizing a transition state ensemble to describe chemical reactions involving proteins may be a useful representation. Here we re-examine the nature of the transition state of protein chemical reactions (enzyme catalysis), considering both recent developments in chemical reaction theory (Marcus theory for SN2 reactions), and protein dynamics effects. The classical theory of chemical reactions relies on the assumption that a reaction must pass through an obligatory transition-state structure. The widely accepted view of enzymatic catalysis holds that there is tight binding of the substrate to the transition-state structure, lowering the activation energy. This picture, may, however, be oversimplified. The real meaning of a transition state is a surface, not a single saddle point on the potential energy surface. In a reaction with a "loose" transition-state structure, the entire transition-state region, rather than a single saddle point, contributes to reaction kinetics. Consequently, here we explore the validity of such a model, namely, the enzymatic modulation of the transition-state surface. We examine its utility in explaining enzyme catalysis. We analyse the possibility that instead of optimizing binding to a well-defined transition-state structure, enzymes are optimized by evolution to bind efficiently with a transition-state ensemble, with a broad range of activated conformations. For enzyme catalysis, the key issue is still transition state (ensemble) stabilization. The source of the catalytic power is the modulation of the transition state. However, our definition of the transition state is the entire transition-state surface rather just than a single well-defined structure. This view of the transition-state ensemble is consistent with the nature of the protein molecule, as embodied and depicted in the protein energy landscape of folding, and binding, funnels.

Animals↗

Improved recognition of native-like protein structures using a combination of sequence-dependent and sequence-independent features of proteins.

We describe the development of a scoring function based on the decomposition P(structure/sequence) proportional to P(sequence/structure) *P(structure), which outperforms previous scoring functions in correctly identifying native-like protein structures in large ensembles of compact decoys. The first term captures sequence-dependent features of protein structures, such as the burial of hydrophobic residues in the core, the second term, universal sequence-independent features, such as the assembly of beta-strands into beta-sheets. The efficacies of a wide variety of sequence-dependent and sequence-independent features of protein structures for recognizing native-like structures were systematically evaluated using ensembles of approximately 30,000 compact conformations with fixed secondary structure for each of 17 small protein domains. The best results were obtained using a core scoring function with P(sequence/structure) parameterized similarly to our previous work (Simons et al., J Mol Biol 1997;268:209-225] and P(structure) focused on secondary structure packing preferences; while several additional features had some discriminatory power on their own, they did not provide any additional discriminatory power when combined with the core scoring function. Our results, on both the training set and the independent decoy set of Park and Levitt (J Mol Biol 1996;258:367-392), suggest that this scoring function should contribute to the prediction of tertiary structure from knowledge of sequence and secondary structure.

Models, Statistical↗

The analysis of the human high affinity IgE receptor Fc epsilon Ri alpha from multiple crystal forms.

We have solved the structure of the human high affinity IgE receptor, Fc epsilon RI alpha, in six different crystal forms, showing the structure in 15 different chemical environments. This database of structures shows no change in the overall shape of the molecule, as the angle between domains 1 and 2 (D1 and D2) varies little across the ensemble. However, the receptor has local conformational variability in the C' strand of D2 and in the BC loop of D1. In every crystal form, a residue inserts between tryptophan residues 87 and 110, mimicking the position of a proline from the IgE ligand. The different crystal forms reveal a distribution of carbohydrates lining the front and back surfaces of the structure. An analysis of crystal contacts in the different forms indicates regions where the molecule interacts with other proteins, and reveals a potential new binding site distal to the IgE binding site. The results of this study point to new directions for the design of molecules to inhibit the interaction of Fc epsilon RI alpha with its natural ligand and thus to prevent a primary step in the allergic response.

Binding Sites↗

Complementarity of structure ensembles in protein-protein binding.

Protein-protein association is often accompanied by changes in receptor and ligand structure. This interplay between protein flexibility and protein-protein recognition is currently the largest obstacle both to our understanding of and to the reliable prediction of protein complexes. We performed two sets of molecular dynamics simulations for the unbound receptor and ligand structures of 17 protein complexes and applied shape-driven rigid body docking to all combinations of representative snapshots. The crossdocking of structure ensembles increased the likelihood of finding near-native solutions. The free ensembles appeared to contain multiple complementary conformations. These were in general not related to the bound structure. We suggest that protein-protein binding follows a three-step mechanism of diffusion, free conformer selection, and refolding. This model combines previously conflicting ideas and is in better agreement with the current data on interaction forces, time scales, and kinetics.

Models, Molecular↗

Direct molecular dynamics observation of protein folding transition state ensemble.

The concept of the protein transition state ensemble (TSE), a collection of the conformations that have 50% probability to convert rapidly to the folded state and 50% chance to rapidly unfold, constitutes the basis of the modern interpretation of protein engineering experiments. It has been conjectured that conformations constituting the TSE in many proteins are the expanded and distorted forms of the native state built around a specific folding nucleus. This view has been supported by a number of on-lattice and off-lattice simulations. Here we report a direct observation and characterization of the TSE by molecular dynamic folding simulations of the C-Src SH3 domain, a small protein that has been extensively studied experimentally. Our analysis reveals a set of key interactions between residues, conserved by evolution, that must be formed to enter the kinetic basin of attraction of the native state.

Computer Simulation↗

Processing multimode binding situations in simulation-based prediction of ligand-macromolecule affinities.

The linear response (LR) approximation and similar approaches belong to practical methods for estimation of ligand-receptor binding affinities. The approaches correlate experimental binding affinities with the changes upon binding of the ligand electrostatic and van der Waals energies and of solvation characteristics. These attributes are expressed as ensemble averages that are obtained by conformational sampling of the protein-ligand complex and of the free ligand by molecular dynamics or Monte Carlo simulations. We observed that outliers in the LR correlations occasionally exhibit major conformational changes of the complex during sampling. We treated the situation as a multimode binding case, for which the observed association constant is the sum of the partial association constants of individual states/modes. The resulting nonlinear expression for the binding affinities contains all the LR variables for individual modes that are scaled by the same two to four adjustable parameters as in the one-mode LR equation. The multimode method was applied to inhibitors of a matrix metalloproteinase, where this treatment improved the explained variance in experimental activity from 75% for the unimode case to about 85%. The predictive ability scaled accordingly, as verified by extensive cross-validations.

Computer Simulation↗

Novel design of bicyclic beta-turn dipeptides on solid-phase supports and synthesis of [3.3.0]-Bicyclo([2,3])-leu-enkephalin analogues.

[structure: see text] External bicyclic beta-turn dipeptide mimetics provide an excellent design approach that can offer a rich chiral ensemble of structures with different backbone conformations. We report herein a novel design of a convergent combinatorial synthetic methodology, which is illustrated by the solid-phase synthesis of a series of [3.3.0]-bicyclo([2,3])-Leu-enkephalin analogues. The reactions were optimized and the epimeric configurations were determined by 2D NMR spectroscopy. Biological assays show that these analogues have more potent delta binding affinity and bioactivity for delta vs micro opioid receptor, which may be related to the different conformations preferred by these analogues in our modeling studies.

Dipeptides↗

Satisfying turns in folding transitions.

Protein engineering studies show that conformations in the folding transition state ensemble can be structurally polarized. In two SH3 beta-sheet domains, the formation of hydrophobic contacts goes hand in hand with the formation of the solvated distal loop beta-turn, while large parts of the molecule remain unstructured in the ensemble.

Evolution, Molecular↗

Fast protein folding kinetics.

Proteins are complex molecules, yet their folding kinetics is often fast (microseconds) and simple, involving only a single exponential function of time (called two-state kinetics). The main model for two-state kinetics has been transition-state theory, where an energy barrier defines a slow step to reach an improbable structure. But how can barriers explain fast processes, such as folding? We study a simple model with rigorous kinetics that explains the high speed instead as a result of the microscopic parallelization of folding trajectories. The single exponential results from a separation of timescales; the parallelization of routes is high at the start of folding and low thereafter. The ensemble of rate-limiting chain conformations is different from in transition-state theory; it is broad, overlaps with the denatured state, is not aligned along a single reaction coordinate, and involves well populated, rather than improbable, structures.

Kinetics↗

A transient expansion of the native state precedes aggregation of recombinant human interferon-gamma.

Aggregation of proteins, even under conditions favoring the native state, is a ubiquitous problem in biotechnology and biomedical engineering. Providing a mechanistic basis for the pathways that lead to aggregation should allow development of rational approaches for its prevention. We have chosen recombinant human interferon-gamma (rhIFN-gamma) as a model protein for a mechanistic study of aggregation. In the presence of 0.9 M guanidinium hydrochloride, rhIFN-gamma aggregates with first order kinetics, a process that is inhibited by addition of sucrose. We describe a pathway that accounts for both the observed first-order aggregation of rhIFN-gamma and the effect of sucrose. In this pathway, aggregation proceeds through a transient expansion of the native state. Sucrose shifts the equilibrium within the ensemble of rhIFN-gamma native conformations to favor the most compact native species over more expanded ones, thus stabilizing rhIFN-gamma against aggregation. This phenomenon is attributed to the preferential exclusion of sucrose from the protein surface. In addition, kinetic analysis combined with solution thermodynamics shows that only a small (9%) expansion surface area is needed to form the transient native state that precedes aggregation. The approaches used here link thermodynamics and aggregation kinetics to provide a powerful tool for understanding both the pathway of protein aggregation and the rational use of excipients to inhibit the process.

Chromatography, Gel↗

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble.

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.

Animals↗

Solution structure of the main alpha-amylase inhibitor from amaranth seeds.

The most abundant alpha-amylase inhibitor (AAI) present in the seeds of Amaranthus hypochondriacus, a variety of the Mexican crop plant amaranth, is the smallest polypeptide (32 residues) known to inhibit alpha-amylase activity of insect larvae while leaving that of mammals unaffected. In solution, 1H NMR reveals that AAI isolated from amaranth seeds adopts a major trans (70%) and minor cis (30%) conformation, resulting from slow cis-trans isomerization of the Val15-Pro16 peptide bond. Both solution structures have been determined using 2D 1H-NMR spectroscopy and XPLOR followed by restrained energy refinement in the consistent-valence force field. For the major isomer, a total of 563 distance restraints, including 55 medium-range and 173 long-range ones, were available from the NOESY spectra. This rather large number of constraints from a protein of such a small size results from a compact fold, imposed through three disulfide bridges arranged in a cysteine-knot motif. The structure of the minor cis isomer has also been determined using a smaller constraint set. It reveals a different backbone conformation in the Pro10-Pro20 segment, while preserving the overall global fold. The energy-refined ensemble of the major isomer, consisting of 20 low-energy conformers with an average backbone rmsd of 0.29 +/- 0.19 A and no violations larger than 0.4 A, represents a considerable improvement in precision over a previously reported and independently performed calculation on AAI obtained through solid-phase synthesis, which was determined with only half the number of medium-range and long-range restraints reported here, and featured the trans isomer only. The resulting differences in ensemble precision have been quantified locally and globally, indicating that, for regions of the backbone and a good fraction of the side chains, the conformation is better defined in the new solution structure. Structural comparison of the solution structure with the X-ray structure of the inhibitor when bound to its alpha-amylase target in Tenebrio molitor shows that the backbone conformation is only slightly adjusted on complexation, while that of the side chains involved in protein-protein contacts is similar to those present in solution. Therefore, the overall conformation of AAI appears to be predisposed to binding to its target alpha-amylase, confirming the view that it acts as a lid on top of the alpha-amylase active site.

Amaranth Dye↗

The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation.

The structural and thermodynamic properties of a 6-residue beta-peptide that was designed to form a hairpin conformation have been studied by NMR spectroscopy and MD simulation in methanol solution. The predicted hairpin would be characterized by a 10-membered hydrogen-bonded turn involving residues 3 and 4, and two extended antiparallel strands. The interproton distances and backbone torsional dihedral angles derived from the NMR experiments at room temperature are in general terms compatible with the hairpin conformation. Two trajectories of system configurations from 100-ns molecular-dynamics simulations of the peptide in solution at 298 and 340 K have been analyzed. In both simulations reversible folding to the hairpin conformation is observed. Interestingly, there is a significant conformational overlap between the unfolded state of the peptide at each of the temperatures. As already observed in previous studies of peptide folding, the unfolded state is composed of a (relatively) small number of predominant conformers and in this case lacks any type of secondary-structure element. The trajectories provide an excellent ground for the interpretation of the NMR-derived data in terms of ensemble averages and distributions as opposed to single-conformation interpretations. From this perspective, a relative population of the hairpin conformation of 20% to 30% would suffice to explain the NMR-derived data. Surprisingly, however, the ensemble of structures from the simulation at 340 K reproduces more accurately the NMR-derived data than the ensemble from the simulation at 298 K, a question that needs further investigation.

Magnetic Resonance Spectroscopy↗

A Monte Carlo method for generating structures of short single-stranded DNA sequences.

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by +/- 10 degrees in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3(7)) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150-325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CXnG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5' end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10(6) conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar-phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques.

Algorithms↗

Thermodynamic cycle between DNA and RNA constituents for conformation of the sugar ring from nuclear magnetic resonance study.

The effect of a structural change of ribose to deoxyribose, by replacement of 2'-OH by 2'-H, on the conformational equilibrium of the sugar ring is described in terms of one thermodynamic cycle. The method is based on the observation that conformational correlations of the sugar ring--side chain ensemble in DNA and RNA components show one general pattern, reflecting an intrinsic physical property of this ensemble. The pattern determines a choice of model systems to study. The systems consist of pairs of DNA and RNA components, nucleosides and nucleotides in aqueous solution, where all conformational factors are fully controlled. This approach allowed us to describe the thermodynamic cycle and measure its fundamental parameters, equilibrium constants and free energy differences, delta delta G, from a nuclear magnetic resonance study. The delta delta G values as determined for pairs of ribo- and deoxyribo-nucleosides in classes of syn-constrained and anti-preferred models, are comparable and lie in a narrow range, delta delta G = 1.7 +/- 0.1 [kJ/mol]. For pairs of ribo- and deoxyribo-nucleotides, the delta delta G values also lie in narrow ranges, delta delta G = 1.7 +/- 0.1 [kJ/mol] for 5'-phosphate nucleotides and delta delta G = 1.9 +/- 0.1 [kJ/mol] for 3'-phosphate nucleotides, i.e. similar to those observed for nucleosides. The measured quantity, delta delta G, is generally observed in a relatively narrow range, delta delta G = 1.75 +/- 0.15 [kJ/mol], irrespective of the class of the model system. This quantity represents a "pure" constant contribution, pe one sugar moiety, as a "driving force" for the N-->S shift in the sugar ring conformational equilibrium, when one compares RNA and DNA. This important thermodynamic quantity, delta delta G, has not hitherto been determined for nucleic acids. Ultimately the delta delta G quantity is revealed in the tendency to adopt S(C2'endo) sugar puckering domain by the majority of DNA structures, whereas RNA generally adopt an N(C3'endo) puckering domain. A possible biological significance of the delta delta G quantity may include evolutionary aspects of nucleic acids.

DNA↗

Specificity of the initial collapse in the folding of the cold shock protein.

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by Förster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.

Ammonium Sulfate↗