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 91 records · Page 5Linked to original sources

PubMed1998

Conformational energy penalties of protein-bound ligands.

The conformational energies required for ligands to adopt their bioactive conformations were calculated for 33 ligand-protein complexes including 28 different ligands. In order to monitor the force field dependence of the results, two force fields, MM3* and AMBER*, were employed for the calculations. Conformational analyses were performed in vacuo and in aqueous solution by using the generalized Born/solvent accessible surface (GB/SA) solvation model. The protein-bound conformations were relaxed by using flat-bottomed Cartesian constraints. For about 70% of the ligand-protein complexes studied, the conformational energies of the bioactive conformations were calculated to be < or = 3 kcal/mol. It is demonstrated that the aqueous conformational ensemble for the unbound ligand must be used as a reference state in this type of calculations. The calculations for the ligand-protein complexes with conformational energy penalties of the ligand calculated to be larger than 3 kcal/mol suffer from uncertainties in the interpretation of the experimental data or limitations of the computational methods. For example, in the case of long-chain flexible ligands (e.g. fatty acids), it is demonstrated that several conformations may be found which are very similar to the conformation determined by X-ray crystallography and which display significantly lower conformational energy penalties for binding than obtained by using the experimental conformation. For strongly polar molecules, e.g. amino acids, the results indicate that further developments of the force fields and of the dielectric continuum solvation model are required for reliable calculations on the conformational properties of this type of compounds.

Binding Sites↗

PubMed2003

13C CPMAS NMR studies of the elastin-like polypeptide (LGGVG)n.

The elucidation of structure-function relationships in insoluble elastin is often approached using elastin-like polypeptides. In this manner, the characterization of the different regions in this extensive biopolymer may be facilitated in a "piece-wise" manner. Our solid-state NMR experiments indicate that (LGGVG)n has structural similarities to elastin and some elastin peptides, providing support for the utility of the mimetic peptides. Furthermore, previous NMR and CD studies indicated that the structure of the elastin-like polypeptide (LGGVG)n in solution is best described as a "conformational ensemble" with a mixture of type I and II beta-turns, in addition to unfolded regions. Our data indicate that the peptide does not adopt a single conformation in the solid state, lending further support to models for elastin that involve significant conformational heterogeneity.

Animals↗

PubMed2007

Simulation of beta-depsipeptides: the effect of missing hydrogen-bond donors on their folding equilibria.

beta-Depsipeptides are beta-peptides in which one or more peptide linkages are replaced by ester linkages, resulting in a loss of a hydrogen-bond donor (N--H) and weakening of the corresponding carbonyl hydrogen-bond acceptor moiety. The effects of three of such peptide by ester substitutions in a hepta-beta-peptide upon its (un)folding equilibrium in methanol solution are investigated using molecular dynamics simulations and compared to experimental data from NMR spectroscopy. The simulated conformational ensembles largely reproduce the experimentally measured NOE and 3J-coupling constant data for the three different hepta-beta-peptides, and confirm the relative stabilities of the 3(14)-helical conformation, which is most weakened by substitution of the 4th peptide linkage and least by substitution of the 6th peptide linkage. The simulations are complementary to the experimental data by providing detailed insight into the conformational distributions that are compatible with the experimentally measured average values of observables.

Computational Biology↗

PubMed2006

Ensemble-based signatures of energy propagation in proteins: a new view of an old phenomenon.

The ability of a protein to transmit the energetic effects of binding from one site to another constitutes the underlying basis for allosterism and signal transduction. Despite clear experimental evidence indicating the ability of proteins to transmit the effects of binding, the means by which this propagation is facilitated is not well understood. Using our previously developed ensemble-based description of the equilibrium, we investigated the physical basis of energy propagation and identified several fundamental and general aspects of energetic coupling between residues in a protein. First, partitioning of a conformational ensemble into four distinct sub-ensembles allows for explanation of the range of experimentally observed coupling behaviors (i.e., positive, neutral, and negative coupling between various regions of the protein structure). Second, the relative thermodynamic properties of these four sub-ensembles define the energetic coupling between residues as either positive, neutral, or negative. Third, analysis of the structural and thermodynamic features of the states within each sub-ensemble reveals significant variability. This third result suggests that a quantitative description of energy propagation in proteins requires an understanding of the structural and energetic features of more than just one or a few low-energy states, but also of many high-energy states. Such findings illuminate the difficulty in interpreting energy propagation in proteins in terms of a structural pathway that physically links coupled sites.

Binding Sites↗

PubMed2004

Modeling correlated main-chain motions in proteins for flexible molecular recognition.

We describe a new method for modeling protein and ligand main-chain flexibility, and show its ability to model flexible molecular recognition. The goal is to sample the full conformational space, including large-scale motions that typically cannot be reached in molecular dynamics simulations due to the computational intensity, as well as conformations that have not been observed yet by crystallography or NMR. A secondary goal is to assess the degree of flexibility consistent with protein-ligand recognition. Flexibility analysis of the target protein is performed using the graph-theoretic algorithm FIRST, which also identifies coupled networks of covalent and noncovalent bonds within the protein. The available conformations of the flexible regions are then explored with ROCK by random-walk sampling of the rotatable bonds. ROCK explores correlated motions by only sampling dihedral angles that preserve the coupled bond networks in the protein and generates conformers with good stereochemistry, without using a computationally expensive potential function. A representative set of the conformational ensemble generated this way can be used as targets for docking with SLIDE, which handles the flexibility of protein and ligand side-chains. The realism of this protein main-chain conformational sampling is assessed by comparison with time-resolved NMR studies of cyclophilin A motions. ROCK is also effective for modeling the flexibility of large cyclic and polycyclic ligands, as demonstrated for cyclosporin and zearalenol. The use of this combined approach to perform docking with main-chain flexibility is illustrated for the cyclophilin A-cyclosporin complex and the estrogen receptor in complex with zearalenol, while addressing the question of how much flexibility is allowed without hindering molecular recognition.

Biochemical Phenomena↗

PubMed1992

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.

Amino Acid Sequence↗

PubMed1991

Conformations of neurotensin in solution and in membrane environments studied by 2-D NMR spectroscopy.

Two-dimensional HOHAHA and ROESY nuclear magnetic resonance techniques are used to obtain complete proton resonance assignments and to perform a conformational investigation of the neuropeptide neurotensin (pGlu-Leu-Tyr-Glu-Asn-Lys-Pro-Arg-Arg-Pro-Tyr-Ile-Leu) in aqueous solution, methanol, and membrane-mimetic [deuterated sodium dodecylsulfate (SDS)] environments. Results suggest the absence of discernible elements of secondary structure in water and methanol. ROESY spectra confirm that Lys-Pro and Arg-Pro peptide bonds are all-trans, but that a significant population of cis Arg-Pro bonds arises in aqueous solution, which increases in the environment of SDS micelles. The conformational ensemble of the peptide is observed to narrow as it becomes bound through its cationic mid-region to SDS micelles, with the accompanying advent of local extended structure. The overall results indicate the inherent conformational flexibility of neurotensin, and emphasize the environmental dependence of conformation in peptides of medium length.

Amino Acid Sequence↗

PubMed2006

Spin pair geometry revealed by high-field DEER in the presence of conformational distributions.

Orientation selection on two nitroxide-labelled shape-persistent molecules is demonstrated by high-field pulsed electron-electron double resonance experiments at a frequency of 95 GHz with a commercial spectrometer. The experiments are performed with fixed observer and pump frequencies by variation of the magnetic field, so that the variation of both the dipolar frequencies and the modulation depths can be analyzed. By applying the deadtime-free four-pulse double electron-electron resonance (DEER) sequence, the lineshapes of the dipolar spectra are obtained. In the investigated linear biradical and equilateral triradical the nitroxide labels undergo restricted dynamics, so that their relative orientations are not fixed, but are correlated to some extent. In this situation, the general dependence of the dipolar spectra on the observer field can be satisfyingly modelled by simple geometrical models that involve only one rotational degree of freedom for the biradical and two rotational degrees of freedom for the triradical. A somewhat better agreement of the dipolar lineshapes for the biradical is obtained by simulations based on a molecular dynamics trajectory. For the triradical, small but significant deviations of the lineshape are observed with both models, indicating that the technique can reveal deficiencies in modelling of the conformational ensemble of a macromolecule.

Algorithms↗

PubMed2002

Principal component method for assessing structural heterogeneity across multiple alignment media.

The recent availability of residual dipolar coupling measurements in a variety of different alignment media raises the question to what extent biomolecular structure and dynamics are differentially affected by their presence. A computational method is presented that allows the sensitive assessment of such changes using dipolar couplings measured in six or more alignment media. The method is based on a principal component analysis of the covariance matrix of the dipolar couplings. It does not require a priori structural or dynamic information nor knowledge of the alignment tensors and their orientations. In the absence of experimental errors, the covariance matrix has at most five nonzero eigenvalues if the structure and dynamics of the biomolecule is the same in all media. In contrast, differential structural and dynamic changes lead to additional nonzero eigenvalues. Characteristic features of the eigenvalue distribution in the absence and presence of noise are discussed using dipolar coupling data calculated from conformational ensembles taken from a molecular dynamics trajectory of native ubiquitin.

Computational Biology↗

PubMed2001

Evaluation of site-directed spin labeling for characterizing protein-ligand complexes using simulated restraints.

Simulation studies have been performed to evaluate the utility of site-directed spin labeling for determining the structures of protein-ligand complexes, given a known protein structure. Two protein-ligand complexes were used as model systems for these studies: a 1.9-A-resolution x-ray structure of a dihydrofolate reductase mutant complexed with methotrexate, and a 1.5-A-resolution x-ray structure of the V-Src tyrosine kinase SH2 domain complexed with a five-residue phosphopeptide. Nitroxide spin labels were modeled at five dihydrofolate reductase residue positions and at four SH2 domain residue positions. For both systems, after energy minimization, conformational ensembles of the spin-labeled residues were generated by simulated annealing while holding the remainder of the protein-ligand complex fixed. Effective distances, simulating those that could be obtained from (1)H-NMR relaxation measurements, were calculated between ligand protons and the spin labels. These were converted to restraints with several different levels of precision. Restrained simulated annealing calculations were then performed with the aim of reproducing target ligand-binding modes. The effects of incorporating a few supplementary short-range (< or =5.0 A) distance restraints were also examined. For the dihydrofolate reductase-methotrexate complex, the ligand-binding mode was reproduced reasonably well using relatively tight spin-label restraints, but methotrexate was poorly localized using loose spin-label restraints. Short-range and spin-label restraints proved to be complementary. For the SH2 domain-phosphopeptide complex without the short-range restraints, the peptide did not localize to the correct depth in the binding groove; nevertheless, the orientation and internal conformation of the peptide was reproduced moderately well. Use of the spin-label restraints in conjunction with the short-range restraints resulted in relatively well defined structural ensembles. These results indicate that restraints derived from site-directed spin labeling can contribute significantly to defining the orientations and conformations of bound ligands. Accurate ligand localization appears to require either a few supplementary short-range distance restraints, or relatively tight spin-label restraints, with at least one spin label positioned so that some of the restraints draw the ligand into the binding pocket in the latter case.

Binding Sites↗

PubMed2002

Hydrogen exchange, core modules, and new designed proteins.

A strategy for design of new proteins that mimic folding properties of native proteins is based on peptides modeled on the slow exchange cores of natural proteins. We have synthesized peptides, called core modules, that correspond to the elements of secondary structure that carry the very slowest exchanging amides in a protein. The expectation is that, if soluble in water, core modules will form conformational ensembles that favor native-like structure. Core modules modeled on natural bovine pancreatic trypsin inhibitor have been shown by NMR studies to meet this expectation. The next step toward production of a native state mimic is to further shift the conformational bias of a core module toward more ordered structure by promoting module-module interactions that are mutually stabilizing. For this, two core modules were incorporated into a single molecule by means of a long cross-link. From a panel of several two-module peptides, one very promising lead emerged; it is called BetaCore. BetaCore is monomeric in water and forms a new fold composed of a four-stranded, antiparallel beta-sheet. The single, dominant conformation of BetaCore is characterized by various NMR experiments. Here we compare the individual core module to the two-module BetaCore and discuss the progressive stabilization of intramodule structure and the formation of new intermodule interactions.

Amino Acid Sequence↗

PubMed1999

Apparent radii of the native, stable intermediates and unfolded conformers of the alpha-subunit of tryptophan synthase from E. coli, a TIM barrel protein.

The urea-induced equilibrium unfolding of the alpha-subunit of tryptophan synthase (alphaTS) from Escherichia coli can be described by a four-state model, N right harpoon over left harpoon I1 right harpoon over left harpoon I2 right harpoon over left harpoon U, involving two highly populated intermediates, I1 and I2 [Gualfetti, P. J., Bilsel, O., and Matthews, C. R. (1999) Protein Sci. 8, 1623-1635]. To extend the physical characterization of these stable forms, the apparent radius was measured by several techniques. Size-exclusion chromatography (SEC), analytical ultracentrifugation (UC), and dynamic light scattering (DLS) experiments yield an apparent Stokes radius, R(s), of approximately 24 A for the native state of alphaTS. The small-angle X-ray scattering (SAXS) experiment yields a radius of gyration, R(g), of 19.1 A, consistent with the value predicted from the X-ray structure and the Stokes radius. As the equilibrium is shifted to favor I1 at approximately 3.2 M and I2 at 5.0 M urea, SEC and UC show that R(s) increases from approximately 38 to approximately 52 A. Measurements of the radius by DLS and SAXS between 2 and 4.5 M urea were complicated by the self-association of the I1 species at the relatively high concentrations required by those techniques. Above 6 M urea, SEC and UC reveal that R(s) increases linearly with increasing urea concentration to approximately 54 A at 8 M urea. The measurements of R(s) by DLS and R(g) by SAXS are sufficiently imprecise that both values appear to be identical for the I2 and U states and, considering the errors, are in good agreement with the results from SEC and UC. Thermodynamic parameters extracted from the SEC data for the N right harpoon over left harpoon I1 and I1 right harpoon over left harpoon I2 transitions agree with those from the optical data, showing that this technique accurately monitors a part of the equilibrium model. The lack of sensitivity to the I2 right harpoon over left harpoon U transition, beyond a simple swelling of both species with increasing urea concentration, implies that the Stokes radii for the I2 and U states are not distinguishable. Surprisingly, the hydrophobic core known to stabilize I2 at 5.0 M urea [Saab-Rincón, G., Gualfetti, P. J., and Matthews, C. R. (1996) Biochemistry 35, 1988-1994] develops without a significant contraction of the polypeptide, i.e., beyond that experienced by the unfolded form at decreasing urea concentrations. Kratky plots of the SAXS data, however, reveal that I2, similar to N and I1, has a globular structure while U has a more random coil-like form. By contrast, the formation of substantial secondary structure and the burial of aromatic side chains in I1 and, eventually, N are accompanied by substantial decreases in their Stokes radii and, presumably, the size of their respective conformational ensembles.

Chromatography, Gel↗

PubMed1994

Three-dimensional structure of a type VI turn in a linear peptide in water solution. Evidence for stacking of aromatic rings as a major stabilizing factor.

Structures have been calculated for the folded conformation found at high population in the cis isomeric form of the peptide NH3(+)-Ser-Tyr-Pro-Phe-Asp-Val-COO- (SYPFDV) in aqueous solution, using distance geometry and restrained molecular dynamics. A number of NMR parameters, including NOE distance restraints and phi and chi 1 dihedral angle restraints derived from coupling constants were used in the calculation. The restraints were carefully selected to exclude those that might have a significant contribution from the unfolded states of the peptide, so that the calculated structure represents only the folded form (present at about 70% population) and contains no information on the other members of the conformational ensemble. The calculations give a tight family of structures in the region Tyr2 to Phe4, all containing a type VI turn with an unusual disposition of the aromatic rings of Tyr and Phe, which stack against the proline ring. Both backbone and side-chain conformations are well defined in this region. At the ends of the molecule the polypeptide chain is disordered. The structures are consistent with all of the available NMR information, including upfield chemical shifts observed for the proline ring protons in the cis isomer, and with the independently determined proline ring pucker. There is no evidence for cross-turn hydrogen bonding. According to the calculated structures, the major source of stabilization of the turn conformation appears to be the stacking of the aromatic and proline rings.

Amino Acid Sequence↗

PubMedDate not supplied

Metric and multidimensional scaling: efficient tools for clustering molecular conformations.

The application of metric and multidimensional scaling to conformer ensembles was demonstrated in this work. An automated process was devised to cluster and assign group memberships and cluster representatives. The method allows rapid clustering, leading to intuitive results that can be visually inspected. Multidimensional scaling was found to be superior to metric scaling for clustering conformers. The performance of different hierarchical clustering algorithms was compared using multidimensional plots, and the group average method was found to perform best.

Amino Acid Sequence↗

PubMed2004

Native-like conformations are sampled by partially folded and disordered variants of bovine pancreatic trypsin inhibitor.

Partially folded conformational ensembles of bovine pancreatic trypsin inhibitor (BPTI) are accessed by replacing Cys 5, 30, 51, and 55 by alpha-amino-n-butyric acid (Abu) while retaining the disulfide between Cys 14 and 38; the resultant variant is termed [14-38](Abu). Two new analogues with modifications in the beta-turn, P26D27[14-38](Abu) and N26G27K28[14-38](Abu), are compared to partially folded [14-38](Abu), as well as to [R](Abu), the unfolded protein with all six Cys residues replaced by Abu. Structural features of the new analogues of [14-38](Abu) have been determined by circular dichroism (CD), one-dimensional (1)H NMR, and 8-anilino-1-naphthalenesulfonic acid (ANS) fluorescence experiments. Both analogues are more disordered than the parent [14-38](Abu), but while P26D27[14-38](Abu) has a small population of native-like conformations observed by NMR, no ordered structure is detected for N26G27K28[14-38](Abu). Trypsin inhibition assays were carried out using a modified rat trypsin, C191A/C220A, that minimizes cleavage of unfolded peptides. Both [14-38](Abu) and P26D27[14-38](Abu) significantly inhibit modified trypsin. N26G27K28[14-38](Abu) has low but measurable inhibitor activity, while [R](Abu) has no activity even when in very high molar excess relative to trypsin. ANS fluorescence is enhanced by [14-38](Abu) and by both variants but not by [R](Abu). We conclude that partially folded ensembles of BPTI, even those with little or no CD- or NMR-detectable structure, contain minor populations of native-like conformations. Partially folded [14-38](Abu) and both variants, as well as [R](Abu), have enhanced negative ellipticity in CD spectra acquired in the presence of the osmolyte trimethylamine N-oxide (TMAO). TMAO-induced structure is formed cooperatively, as indicated by thermal unfolding curves. Inhibitor activity as a function of TMAO concentration implies that the osmolyte-induced structure is native-like for [14-38](Abu) and P26D27[14-38](Abu) and is probably native-like for N26G27K28[14-38](Abu). [R](Abu) also shows increased CD-detected structure in the presence of TMAO, but such structure is likely to be collapsed and non-native.

Amino Acid Substitution↗

PubMed2006

Multiple intermediate conformations of jack bean urease at low pH: anion-induced refolding.

Structural and functional characteristics of jack bean urease (JBU), a hexameric enzyme having identical subunits, were investigated under neutral as well as acidic conditions by using CD, fluorescence, ANS binding and enzyme activity measurements. At low pH and low ionic strength, JBU exists in a partially unfolded state (U(A)-state), having predominantly beta structure and no tertiary interactions along with a strong ANS binding. Addition of salts like NaCl, KCl and Na(2)SO(4) to the U(A)-state induces refolding resulting in structural propensities similar to that of native hexamer. Moreover, at low concentrations, GuHCl behaves like an anion by inducing refolding of the U(A)-state. The anion-induced refolded state (I(A)-state) is more stable than U(A)-state and the stability is nearly equal to that of the native protein against chemical-induced and thermal denaturation. Overall, these observations support a model of protein folding for a multimeric protein where certain conformations (ensembles of substates) of low energy prevail and populated under non-native conditions with different stability.

Anions↗

PubMed2005

Probing conformational disorder in neurotensin by two-dimensional solid-state NMR and comparison to molecular dynamics simulations.

An approach is introduced to characterize conformational ensembles of intrinsically unstructured peptides on the atomic level using two-dimensional solid-state NMR data and their combination with molecular dynamics simulations. For neurotensin, a peptide that binds with high affinity to a G-protein coupled receptor, this method permits the investigation of the changes in conformational preferences of a neurotransmitter transferred from a frozen aqueous solution via a lipid model phase to the receptor-bound form. The results speak against a conformational pre-organization of the ligand in detergents in which the receptor has been shown to be functional. Further extensions to the study of protein folding are possible.

Detergents↗

PubMed1995

Refined solution structure of human profilin I.

Profilin is a ubiquitous eukaryotic protein that binds to both cytosolic actin and the phospholipid phosphatidylinositol-4,5-bisphosphate. These dual competitive binding capabilities of profilin suggest that profilin serves as a link between the phosphatidyl inositol cycle and actin polymerization, and thus profilin may be an essential component in the signaling pathway leading to cytoskeletal rearrangement. The refined three-dimensional solution structure of human profilin I has been determined using multidimensional heteronuclear NMR spectroscopy. Twenty structures were selected to represent the solution conformational ensemble. This ensemble of structures has root-mean-square distance deviations from the mean structure of 0.58 A for the backbone atoms and 0.98 A for all non-hydrogen atoms. Comparison of the solution structure of human profilin to the crystal structure of bovine profilin reveals that, although profilin adopts essentially identical conformations in both states, the solution structure is more compact than the crystal structure. Interestingly, the regions that show the most structural diversity are located at or near the actin-binding site of profilin. We suggest that structural differences are reflective of dynamical properties of profilin that facilitate favorable interactions with actin. The global folding pattern of human profilin also closely resembles that of Acanthamoeba profilin I, reflective of the 22% sequence identity and approximately 45% sequence similarity between these two proteins.

Actins↗