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 235 records · Page 13Linked to original sources

PubMed2005

The folding energy landscape and phosphorylation: modeling the conformational switch of the NFAT regulatory domain.

An energy landscape approach predicts the conformational changes of the configurations of the regulatory domain of the protein nuclear factor of activated T cells (NFAT) caused by phosphorylation of specific multiple sites. Structurally local effects and secondary structural changes are modeled using all-atom Brownian dynamics to investigate the changes of the backbone torsional distributions upon phosphorylation. For tertiary and global changes, we employ a coarse-grained model to sample ensembles of conformations both with and without phosphorylation. At the secondary structure level, phosphorylation moderately increases the helical propensity and gives a more rigid local backbone conformation. The tertiary effects of phosphorylation caused by the extensive charge modification are more pronounced and collectively change the conformation of the regulatory domain of NFAT from a flexible globular ensemble to a rather rigid helical bundle, blocking access to the nuclear localization sequence. These studies give computational support to one scenario conjectured from experiments.

Amino Acid Sequence↗

PubMed2002

Stability of protein-bound conformations of bioactive peptides: the folded conformation of an epidermal growth factor-like thrombomodulin fragment is similar to that recognized by thrombin.

The relationship between the free and bound conformations of bioactive peptides is explored using the epidermal growth factor (EGF)-like thrombomodulin fragment hTM409-426 as a model system. The hTM409-426 peptide has a sequence of C(409)PEGYILDDGFIC(421)TDIDE (with a disulfide bond between Cys409 and Cys421) and is a selective inhibitor of thrombin. Upon binding to thrombin, hTM409-426 adopts a well-defined conformation-namely, a beta-turn followed by an antiparallel beta-sheet, similar to those found in all other EGF-like protein repeats (Hrabal et al., Protein Science, 1996, Vol. 5, 195-203). Here we demonstrate that, at pH 6.8 and at 25 degrees C, the hTM409-426 peptide in the free state is very flexible, but still populates a type II beta-turn over residues Pro410-Glu411-Gly412-Tyr413 and the clustering of some hydrophobic side chains, both of which are present in the thrombin-bound conformation. At a lower temperature of 5 degrees C, significant conformational shifts of the C alpha H proton resonances and extensive medium- and long-range NOEs are observed, indicating the presence of folded conformations with unique backbone-backbone and side-chain interactions. A comparison of the NOE patterns in the free state with transferred NOEs shows that the free-state folded and the thrombin-bound conformations of the hTM409-426 peptide are very similar, particularly over residues Pro410-Ile424. The folded conformation of hTM409-426 appears to be stabilized by two hydrophobic clusters, one formed by the side chains of residues Pro410, Tyr413, Leu415, and Phe419 and the Cys409-Cys421 disulfide bond, the second involving residues Ile414 and Ile424. These results indicate that the overall topology of the thrombin-bound conformation of the hTM409-426 peptide is prefolded in the free state and the primary sequence (including the disulfide bond) may be selective for an ensemble of conformations similar to that recognized by thrombin.

Magnetic Resonance Spectroscopy↗

PubMed1998

Improving the accuracy of protein pKa calculations: conformational averaging versus the average structure.

Several methods for including the conformational flexibility of proteins in the calculation of titration curves are compared. The methods use the linearized Poisson-Boltzmann equation to calculate the electrostatic free energies of solvation and are applied to bovine pancreatic trypsin inhibitor (BPTI) and hen egg-white lysozyme (HEWL). An ensemble of conformations is generated by a molecular dynamics simulation of the proteins with explicit solvent. The average titration curve of the ensemble is calculated in three different ways: an average structure is used for the pKa calculation; the electrostatic interaction free energies are averaged and used for the pKa calculation; and the titration curve for each structure is calculated and the curves are averaged. The three averaging methods give very similar results and improve the pKa values to approximately the same degree. This suggests, in contrast to implications from other work, that the observed improvement of pKa values in the present studies is due not to averaging over an ensemble of structures, but rather to the generation of a single properly averaged structure for the pKa calculation.

Aprotinin↗

PubMed1996

Structure-based calculation of the equilibrium folding pathway of proteins. Correlation with hydrogen exchange protection factors.

A new statistical thermodynamic formalism has been developed in order to describe the equilibrium folding pathway of proteins. The resulting formalism allows calculation of the probabilities that individual amino acid residues will be in a native or native-like conformation for any given degree of folding of the protein molecule. The residue probabilities are defined by the probability distribution of conformational states and can be used to calculate experimental quantities like native-state, hydrogen exchange protection factors. A combinatorial algorithm aimed at generating a large ensemble of conformational states (10(4) to 10(6)) using the native structure as a template has been developed. The Gibbs energy and corresponding probability of each conformational state is estimated by using a previously developed structural parametrization of the energetics. The approach has been applied to five different proteins: hen egg-white lysozyme, equine lysozyme, bovine pancreatic trypsin inhibitor, staphylococcal nuclease and turkey ovomucoid third domain. The validity of the approach has been tested by comparing predicted and experimental hydrogen exchange protection factors. It is shown that for the above proteins 76%, 73%, 74%, 78% and 81% of all observed protection factors are predicted correctly. Furthermore, on average, the magnitude of the predicted protection factors, expressed as apparent free energies per residue deviate less than 1 kcal/mol from those obtained experimentally. These results represent the first attempt at predicting both the location and magnitude of hydrogen exchange protection factors from the high-resolution structure of a protein. The good agreement between experimental and predicted values has permitted a close examination of the nature of the equilibrium folding intermediates existing under conditions of maximal stability of the native state.

Animals↗

PubMed2003

Potent side-chain to side-chain cyclized dermorphin analogues containing a carbonyl bridge.

A new family of cyclic opioid peptide analogues related to the 1-4 sequence of dermorphin/deltorphin (Tyr-D-Aaa2-Phe-Aaa4-NH2) has been synthesized. The synthesis of the linear precursor peptides was accomplished by the solid-phase method and ring formation was achieved via a ureido group incorporating the side chain amino functions of D-Aaa2 (D-Lys, D-Orn) and Aaa4 (Lys, Orn, Dab, Dap). The peptides were tested in the guinea-pig ileum (GPI) and mouse vas deferens (MVD) assays. Most showed very high agonist potency in the GPI assay. The peptide containing D-Lys in position 2 and Dab in position 4 was 210 times more active than enkephalin, and that containing Orn and Dab, respectively, was 150 times more active than enkephalin. The latter peptide was also very active in the MVD assay, and showed an IC50 MVD/GPI ratio of 0.816. NMR spectra of selected peptides were recorded, and structural parameters were determined. The conformational space of the peptides was examined using the electrostatically driven Monte Carlo method. With the help of the NMR spectra each peptide was described as an ensemble of conformations. The conformations have been interpreted with regard to the opioid activities, and comparisons have been made with a model proposed earlier for enkephalin analogues.

Analgesics, Opioid↗

PubMed2005

Simultaneous determination of protein structure and dynamics.

We present a protocol for the experimental determination of ensembles of protein conformations that represent simultaneously the native structure and its associated dynamics. The procedure combines the strengths of nuclear magnetic resonance spectroscopy--for obtaining experimental information at the atomic level about the structural and dynamical features of proteins--with the ability of molecular dynamics simulations to explore a wide range of protein conformations. We illustrate the method for human ubiquitin in solution and find that there is considerable conformational heterogeneity throughout the protein structure. The interior atoms of the protein are tightly packed in each individual conformation that contributes to the ensemble but their overall behaviour can be described as having a significant degree of liquid-like character. The protocol is completely general and should lead to significant advances in our ability to understand and utilize the structures of native proteins.

Computer Simulation↗

PubMed2000

The Ensemble/Legacy Chimera extension: standardized user and programmer interface to molecular Ensemble data and Legacy modeling programs.

Ensemble/Legacy is a toolkit extension of the Object Technology Framework (OTF) that exposes an object oriented interface for accessing and manipulating ensembles (collections of molecular conformations that share a common chemical topology) and driving Legacy programs (such as MSMS, AMBER, X-PLOR, CORMA/MARDIGRAS, Dials and Windows, and CURVES). Ensemble/Legacy provides a natural programming interface for running Legacy programs on ensembles of molecules and accessing the resulting data. Using the OTF reduces the time cost of developing a new library to store and manipulate molecular data and also allows Ensemble/Legacy to integrate into the Chimera visualization program. The extension to Chimera exposes the Legacy functionality using a graphical user interface that greatly simplifies the process of modeling and analyzing conformational ensembles. Furthermore, all the C++ functionality of the Ensemble/Legacy toolkit is "wrapped" for use in the Python programming language. More detailed documentation on using Ensemble/Legacy is available online (http:¿picasso.nmr.ucsf.edu/dek/ensemble. html).

Computer Graphics↗

PubMed2001

Calculated pH-dependent population and protonation of carbon-monoxy-myoglobin conformers.

X-ray structures of carbonmonoxymyoglobin (MbCO) are available for different pH values. We used conventional electrostatic continuum methods to calculate the titration behavior of MbCO in the pH range from 3 to 7. For our calculations, we considered five different x-ray structures determined at pH values of 4, 5, and 6. We developed a Monte Carlo method to sample protonation states and conformations at the same time so that we could calculate the population of the considered MbCO structures at different pH values and the titration behavior of MbCO for an ensemble of conformers. To increase the sampling efficiency, we introduced parallel tempering in our Monte Carlo method. The calculated population probabilities show, as expected, that the x-ray structures determined at pH 4 are most populated at low pH, whereas the x-ray structure determined at pH 6 is most populated at high pH, and the population of the x-ray structures determined at pH 5 possesses a maximum at intermediate pH. The calculated titration behavior is in better agreement with experimental results compared to calculations using only a single conformation. The most striking feature of pH-dependent conformational changes in MbCO-the rotation of His-64 out of the CO binding pocket-is reproduced by our calculations and is correlated with a protonation of His-64, as proposed earlier.

Animals↗

PubMed2004

Improved protein structure selection using decoy-dependent discriminatory functions.

BACKGROUND: A key component in protein structure prediction is a scoring or discriminatory function that can distinguish near-native conformations from misfolded ones. Various types of scoring functions have been developed to accomplish this goal, but their performance is not adequate to solve the structure selection problem. In addition, there is poor correlation between the scores and the accuracy of the generated conformations. RESULTS: We present a simple and nonparametric formula to estimate the accuracy of predicted conformations (or decoys). This scoring function, called the density score function, evaluates decoy conformations by performing an all-against-all Calpha RMSD (Root Mean Square Deviation) calculation in a given decoy set. We tested the density score function on 83 decoy sets grouped by their generation methods (4state_reduced, fisa, fisa_casp3, lmds, lattice_ssfit, semfold and Rosetta). The density scores have correlations as high as 0.9 with the Calpha RMSDs of the decoy conformations, measured relative to the experimental conformation for each decoy. We previously developed a residue-specific all-atom probability discriminatory function (RAPDF), which compiles statistics from a database of experimentally determined conformations, to aid in structure selection. Here, we present a decoy-dependent discriminatory function called self-RAPDF, where we compiled the atom-atom contact probabilities from all the conformations in a decoy set instead of using an ensemble of native conformations, with a weighting scheme based on the density scores. The self-RAPDF has a higher correlation with Calpha RMSD than RAPDF for 76/83 decoy sets, and selects better near-native conformations for 62/83 decoy sets. Self-RAPDF may be useful not only for selecting near-native conformations from decoy sets, but also for fold simulations and protein structure refinement. CONCLUSIONS: Both the density score and the self-RAPDF functions are decoy-dependent scoring functions for improved protein structure selection. Their success indicates that information from the ensemble of decoy conformations can be used to derive statistical probabilities and facilitate the identification of near-native structures.

Computational Biology↗

PubMed2000

Twist Sense Bias Induced by Chiral Side Chains in Helically Folded Oligomers.

Cooperative interactions among the side chains of the helically folded phenylene-ethynylene oligomer shown (n=2, 4, 6, 8, 10, 12, 14, 16, 18) can induce a twist sense bias. Therefore, the side chains can play more than just an ancillary role in these conformationally ordered oligomers. The onset of the twist sense bias lags significantly behind the appearance of helical conformations, possibly because a large ensemble of "collapsed" conformations is initially formed.

Journal Article↗

PubMed2006

Modeling and active site refinement for G protein-coupled receptors: application to the beta-2 adrenergic receptor.

It is well known that G protein-coupled receptors are prime targets for drug discovery. At the present time there is only one protein from this class that has an X-ray crystal structure, bovine rhodopsin. Crystal structures of rhodopsin have become invaluable templates for the modeling of class-A G protein-coupled receptors as they likely represent the overall topology of this family of proteins. However, because of low sequence homology within the class and the inherent mobility of integral membrane proteins, it is unlikely that this single structural template reflects the ensemble of conformations accessible for any given receptor. We have devised a procedure based upon comparative modeling that uses induced fit modeling coupled with binding site expansion. The modeling protocol enables an ensemble approach to binding mode prediction. The utility of models for beta-2 adrenergic receptor will be discussed.

Amino Acid Motifs↗

PubMed1995

Conformation and hydrogen ion titration of proteins: a continuum electrostatic model with conformational flexibility.

A new method for including local conformational flexibility in calculations of the hydrogen ion titration of proteins using macroscopic electrostatic models is presented. Intrinsic pKa values and electrostatic interactions between titrating sites are calculated from an ensemble of conformers in which the positions of titrating side chains are systematically varied. The method is applied to the Asp, Glu, and Tyr residues of hen lysozyme. The effects of different minimization and/or sampling protocols for both single-conformer and multi-conformer calculations are studied. For single-conformer calculations it is found that the results are sensitive to the choice of all-hydrogen versus polar-hydrogen-only atomic models and to the minimization protocol chosen. The best overall agreement of single-conformer calculations with experiment is obtained with an all-hydrogen model and either a two-step minimization process or minimization using a high dielectric constant. Multi-conformational calculations give significantly improved agreement with experiment, slightly smaller shifts between model compound pKa values and calculated intrinsic pKa values, and reduced sensitivity of the intrinsic pKa calculations to the initial details of the structure compared to single-conformer calculations. The extent of these improvements depends on the type of minimization used during the generation of conformers, with more extensive minimization giving greater improvements. The ordering of the titrations of the active-site residues, Glu-35 and Asp-52, is particularly sensitive to the minimization and sampling protocols used. The balance of strong site-site interactions in the active site suggests a need for including site-site conformational correlations.

Animals↗

PubMed1997

Structural and dynamical properties of a denatured protein. Heteronuclear 3D NMR experiments and theoretical simulations of lysozyme in 8 M urea.

Oxidized and reduced hen lysozyme denatured in 8 M urea at low pH have been studied in detail by NMR methods. 15N correlated NOESY and TOCSY experiments have provided near complete sequential assignment for both 1H and 15N resonances. Over 900 NOEs, including 130 (i, i + 2) and 23 (i, i + 3) NOEs, could be identified by analysis of the NOESY spectra of the denatured states, and 3J(HN, Halpha) coupling constants and 15N relaxation rates have been measured. The coupling constant and NOE data were analyzed by comparisons with theoretical predictions from a random coil polypeptide model based on amino acid specific phi,psi distributions extracted from the protein data bank. There is significant agreement between predicted and experimental NMR parameters suggesting that local conformations of the denatured states are largely determined by short-range interactions within the polypeptide chain. This result is supported by the observation that the chemical shift, coupling constant, and NOE data are little affected by whether or not the four disulfide bridge cross-links are formed in the denatured protein. The relaxation data, however, show significant differences between the oxidized and reduced protein. Analysis of the relaxation data in terms of simple dynamics models provides evidence for weak clustering of hydrophobic groups near tryptophan residues and increased barriers to motion in the more compact conformers formed when the polypeptide chain is cross-linked by the disulfide bridges. Using this information, a structural description of these denatured states is given in terms of an ensemble of conformers, which have a complex relationship between their local and global characteristics.

Amino Acid Sequence↗

PubMed1979

Circular dichroism study of the solution conformation of luteinizing hormone releasing hormone.

A systematic investigation has been made into the circular dichroic behavior of luteinizing hormone releasing hormone and its peptide fragments and deletion analogues. The results are interpreted to mean that the hormone exists in solution as an ensemble of conformers with different sensitivities to temperature and solvent composition. The far-ultraviolet circular dichroic spectra exhibited by the hormone under different experimental conditions can be simulated satisfactorily by the weighted addition of the spectra of its aliphatic- and aromatic-containing halves. However, the structure of the hormone is not simply the sum of its halves, since some conformational feature of the intact molecule perturbs the near-ultraviolet circular dichroism of its aromatic residues.

Amino Acid Sequence↗

PubMed2002

Characterization of the unfolding of ribonuclease a by a pulsed hydrogen exchange study: evidence for competing pathways for unfolding.

The unfolding of ribonuclease A was studied in 5.2 M guanidine hydrochloride at pH 8 and 10 degrees C using multiple optical probes, native-state hydrogen exchange (HX), and pulse labeling by hydrogen exchange. First, native-state HX studies were used to demonstrate that the protein exists in two slowly interconverting forms under equilibrium native conditions: a predominant exchange-incompetent N form and an alternative ensemble of conformations, N(I), in which some amide hydrogens are fully exposed to exchange. Pulsed HX studies indicated that, during unfolding, the rates of exposure to exchange with solvent protons were similar for all backbone NH probe protons. It is shown that two parallel routes of unfolding are available to the predominant N conformation as soon as it encounters strong unfolding conditions. A fraction of molecules appears to rapidly form N(I) on one route. On the other route an exchange-incompetent intermediate state ensemble, I(U)(2), is formed. The kinetics of unfolding measured by far-UV circular dichroism (CD) were faster than those measured by near-UV CD and intrinsic tyrosine fluorescence of the protein. The logarithms of the rate constants of the unfolding reaction measured by all three optical probes also showed a nonlinear dependence on GdnHCl concentration. All of the data suggest that N(I) and I(U)(2) are nativelike in their secondary and tertiary structures. While N(I) unfolds directly to the fully exchange-competent unfolded state (U), I(U)(2) forms another intermediate I(U)(3) which then unfolds to U. I(U)(3) is devoid of all native alpha-helical secondary structure and has only 30% of the tertiary interactions still intact. Since the rates of global unfolding measured by near-UV CD and fluorescence agree well with the rates of exposure determined for all of the backbone NH probe protons, it appears that the rate-limiting step for the unfolding of RNase A is the dissolution of the entire native tertiary structure and penetration of water into the hydrophobic core.

Animals↗

PubMed2006

Conformer generation under restraints.

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

Algorithms↗

PubMed1993

Analysis of the delocalized Raman modes of conformationally disordered polypeptides.

Bands associated with delocalized vibrational modes were identified in the isotropic Raman spectra of a series of polyglycine oligomers in aqueous solution as zwitterions and as cations. The dependence of these bands on conformational disorder and chain length was determined. The observed dependence is closely mimicked in spectra calculated for a series of corresponding model polypeptides. The simulated spectra were calculated in a skeletal approximation for ensembles of conformationally disordered chains. As the chain length of the conformationally disordered polypeptides increases, the observed isotropic spectra rapidly approach the spectrum of the infinitely long disordered chain. Convergence is nearly complete at the tripeptide for both the zwitterion and the cation. The stimulated spectra behave in essentially the same way. Convergence to the spectrum of the infinitely long chain is much more rapid for the conformationally disordered polyglycines than for the ordered polyglycines because of the mode localization that results from disorder. In the low-frequency region the bands in the calculated spectra have frequencies that are systematically dependent on chain length. These bands are related to the longitudinal acoustic modes of the ordered chain.

Biophysical Phenomena↗

PubMed2005

Fluorescence analysis of a dynamic loop in the PCAF/GCN5 histone acetyltransferase.

PCAF and GCN5 are histone acetyltransferase (HAT) paralogs which play roles in the remodeling of chromatin in health and disease. Previously, a conformationally flexible loop in the catalytic domain had been observed in the X-ray structures of GCN5 in different liganded states. Here, the conformation and dynamics of this PCAF/GCN5 alpha5-beta6 loop was investigated in solution using tryptophan fluorescence. A mutant human PCAF HAT domain (PCAF(Wloop)) was created in which the natural tryptophan (Trp-514) remote from the alpha5-beta6 loop was replaced with tyrosine and a glutamate within the loop (Glu-641) was substituted with tryptophan. This PCAF(Wloop) protein exhibited catalytic parameters within 3-fold of those of the wild-type PCAF catalytic domain, suggesting that the loop mutation was not deleterious for HAT activity. While saturating CoASH induced a 30% quenching of Trp fluorescence in PCAF(Wloop), binding of the high-affinity bisubstrate analogue H3-CoA-20 led to a 2-fold fluorescence increase. These different effects correlate with the different alpha5-beta6 loop conformations seen previously in X-ray structures. On the basis of stopped-flow fluorescence studies, binding of H3-CoA-20 to PCAF(Wloop) proceeds via a rapid association step followed by a slower conformational change involving loop movement. Time-resolved fluorescence measurements support a model in which the alpha5-beta6 loop in the H3-CoA-20-PCAF(Wloop) complex exists in a narrower ensemble of conformations compared to free PCAF(Wloop). The relevance of loop dynamics to PCAF/GCN5 catalysis and substrate specificity are discussed.

Animals↗