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 505 records · Page 28Linked to original sources

PubMed2005

The relationship between the flexibility of proteins and their conformational states on forming protein-protein complexes with an application to protein-protein docking.

We investigate the extent to which the conformational fluctuations of proteins in solution reflect the conformational changes that they undergo when they form binary protein-protein complexes. To do this, we study a set of 41 proteins that form such complexes and whose three-dimensional structures are known, both bound in the complex and unbound. We carry out molecular dynamics simulations of each protein, starting from the unbound structure, and analyze the resulting conformational fluctuations in trajectories of 5 ns in length, comparing with the structure in the complex. It is found that fluctuations take some parts of the molecules into regions of conformational space close to the bound state (or give information about it), but at no point in the simulation does each protein as whole sample the complete bound state. Subsequent use of conformations from a clustered MD ensemble in rigid-body docking is nevertheless partially successful when compared to docking the unbound conformations, as long as the unbound conformations are themselves included with the MD conformations and the whole globally rescored. For one key example where sub-domain motion is present, a ribonuclease inhibitor, principal components analysis of the MD was applied and was also able to produce conformations for docking that gave enhanced results compared to the unbound. The most significant finding is that core interface residues show a tendency to be less mobile (by size of fluctuation or entropy) than the rest of the surface even when the other binding partner is absent, and conversely the peripheral interface residues are more mobile. This surprising result, consistent across up to 40 of the 41 proteins, suggests different roles for these regions in protein recognition and binding, and suggests ways that docking algorithms could be improved by treating these regions differently in the docking process.

Cluster Analysis↗

PubMed2003

Persistent conformational heterogeneity of triosephosphate isomerase: separation and characterization of conformational isomers in solution.

Subunit dissociation of dimeric rabbit muscle triosephosphate isomerase (TIM) by hydrostatic pressure has previously been shown not to follow the expected dependence on protein concentration [Rietveld and Ferreira (1996) Biochemistry 35, 7743-7751]. This anomalous behavior was attributed to persistent conformational heterogeneity (i.e., the coexistence of long-lived conformational isomers) in the ensemble of TIM dimers. Here, we initially show that subunit dissociation/unfolding of TIM by guanidine hydrochloride (GdnHCl) also exhibits an anomalous dependence on protein concentration. Dissociation/unfolding of TIM by GdnHCl was investigated by intrinsic fluorescence and circular dichroism spectroscopies and was found to be a highly cooperative transition in which the tertiary and secondary structures of the protein were concomitantly lost. A procedure based on size-exclusion chromatography in the presence of intermediate (0.6 M) GdnHCl concentrations was developed to isolate two conformational isomers of TIM that exhibit significantly different stabilities and kinetics of unfolding by GdnHCl. Complete unfolding of the two isolated conformers at a high GdnHCl concentration (1.5 M), followed by refolding by removal of the denaturant, completely abolished the differences in their unfolding kinetics. These results indicate that such differences stem from conformational heterogeneity of TIM and are not related to any chemical modification of the protein. Furthermore, they add support to the notion that long-lived conformational isomers of TIM coexist in solution and provide a basis for the interpretation of the persistent heterogeneity of this protein.

Animals↗

PubMed2003

TOUCHSTONE II: a new approach to ab initio protein structure prediction.

We have developed a new combined approach for ab initio protein structure prediction. The protein conformation is described as a lattice chain connecting C(alpha) atoms, with attached C(beta) atoms and side-chain centers of mass. The model force field includes various short-range and long-range knowledge-based potentials derived from a statistical analysis of the regularities of protein structures. The combination of these energy terms is optimized through the maximization of correlation for 30 x 60,000 decoys between the root mean square deviation (RMSD) to native and energies, as well as the energy gap between native and the decoy ensemble. To accelerate the conformational search, a newly developed parallel hyperbolic sampling algorithm with a composite movement set is used in the Monte Carlo simulation processes. We exploit this strategy to successfully fold 41/100 small proteins (36 approximately 120 residues) with predicted structures having a RMSD from native below 6.5 A in the top five cluster centroids. To fold larger-size proteins as well as to improve the folding yield of small proteins, we incorporate into the basic force field side-chain contact predictions from our threading program PROSPECTOR where homologous proteins were excluded from the data base. With these threading-based restraints, the program can fold 83/125 test proteins (36 approximately 174 residues) with structures having a RMSD to native below 6.5 A in the top five cluster centroids. This shows the significant improvement of folding by using predicted tertiary restraints, especially when the accuracy of side-chain contact prediction is >20%. For native fold selection, we introduce quantities dependent on the cluster density and the combination of energy and free energy, which show a higher discriminative power to select the native structure than the previously used cluster energy or cluster size, and which can be used in native structure identification in blind simulations. These procedures are readily automated and are being implemented on a genomic scale.

Algorithms↗

PubMed2006

De novo ligand design to an ensemble of protein structures.

We describe a combinatorial method for de novo ligand design to an ensemble of receptor structures. Receptor conformations, protonation states, and structural water molecules are considered consistently within the framework of de novo ligand design. The method relies on Monte Carlo optimization to search the space of ligand structures, conformations, and rigid-body movements as well as receptor models. The method is applied to an ensemble of HIV protease and human collagenase receptor models. Ligand structures generated de novo exhibit the correct hydrogen-bonding pattern in the core of the active site, with hydrophobic groups extending into the receptor S1 and S1' pocket space. Furthermore, it is shown that known ligands are recovered in the correct binding mode and in the native, most tightly binding receptor model.

Binding Sites↗

PubMed1996

Protein folding funnels: the nature of the transition state ensemble.

BACKGROUND: Energy landscape theory predicts that the folding funnel for a small fast-folding alpha-helical protein will have a transition state half-way to the native state. Estimates of the position of the transition state along an appropriate reaction coordinate can be obtained from linear free energy relationships observed for folding and unfolding rate constants as a function of denaturant concentration. The experimental results of Huang and Oas for lambda repressor, Fersht and collaborators for C12, and Gray and collaborators for cytochrome c indicate a free energy barrier midway between the folded and unfolded regions. This barrier arises from an entropic bottleneck for the folding process. RESULTS: In keeping with the experimental results, lattice simulations based on the folding funnel description show that the transition state is not just a single conformation, but rather an ensemble of a relatively large number of configurations that can be described by specific values of one or a few order parameters (e.g. the fraction of native contacts). Analysis of this transition state or bottleneck region from our lattice simulations and from atomistic models for small alpha-helical proteins by Boczko and Brooks indicates a broad distribution for native contact participation in the transition state ensemble centered around 50%. Importantly, however, the lattice-simulated transition state ensemble does include some particularly hot contacts, as seen in the experiments, which have been termed by others a folding nucleus. CONCLUSIONS: Linear free energy relations provide a crude spectroscopy of the transition state, allowing us to infer the values of a reaction coordinate based on the fraction of native contacts. This bottleneck may be thought of as a collection of delocalized nuclei where different native contacts will have different degrees of participation. The agreement between the experimental results and the theoretical predictions provides strong support for the landscape analysis.

Animals↗

PubMed2002

TROSY experiment for refinement of backbone psi and phi by simultaneous measurements of cross-correlated relaxation rates and 3,4J(H alpha HN) coupling constants.

The TROSY principle has been introduced into a HNCA experiment, which is designed for measurements of the intraresidual and sequential H(alpha)-C(alpha)/H(N)-N dipole/dipole and H(alpha)-C(alpha)/N dipole/CSA cross-correlated relaxation rates. In addition, the new experiment provides values of the (3,4)J(H alpha HN) coupling constants measured in an E.COSY manner. The conformational restraints for the psi and phi angles are obtained through the use of the cross-correlated relaxation rates together with the Karplus-type dependencies of the coupling constants. Improved signal-to-noise is achieved through preservation of all coherence transfer pathways and application of the TROSY principle. The application of the [(15)N,(13)C]-DQ/ZQ-[(15)N,(1)H]-TROSY-E.COSY experiment to the 16 kDa apo-form of the E. coli Heme Chaperon protein CcmE is described. Overall good agreement is achieved between psi and phi angles measured with the new experiment and the average values determined from an ensemble of 20 NMR conformers.

Bacterial Outer Membrane Proteins↗

PubMed2000

Synthesis and characterization of a beta-hairpin peptide that represents a 'core module' of bovine pancreatic trypsin inhibitor (BPTI).

A new strategy for the design and construction of peptide fragments that can achieve defined, nativelike secondary structure is presented. The strategy is based upon the hypothesis that 'core elements' of a protein, synthesized in a single polypeptide molecule, will favor nativelike structure, and that by incorporating a cross-link, nativelike core structure will dominate the ensemble as the more extended conformations are excluded. 'Core elements' are the elements of packed secondary structure that contain the slowest exchanging backbone amide protons in the native protein. The 'core elements' in bovine pancreatic trypsin inhibitor (BPTI) are the two long strands of antiparallel beta-sheet (residues 18-24 and 29-35) and the small beta-bridge (residues 43-44). To test the design strategy, we synthesized an 'oxidized core module', which contains the antiparallel strands connected by a modified reverse turn (A27 replaced by D), a natural disulfide cross-link at the open end of the hairpin, and N- and C-termini blocking groups. A peptide with identical sequence but lacking the disulfide cross-link at the open end was used as the 'reduced core module' control. The conformational behavior of both peptides was examined using (1)H NMR spectroscopy. Chemical shift dispersion, chemical shift deviation from random coil values, sequential and long-range NOEs, and H/D amide exchange rates were compared for the two peptides. We conclude that the ensemble of oxidized and reduced core module conformations samples both nativelike 4:4 and non-native 3:5 beta-hairpin structure, and that the oxidized module samples nativelike structure for a greater fraction of the time than the reduced module.

Amino Acid Sequence↗

PubMed1999

On the pH-conformational dependence of the unblocked SYPYD peptide.

Simulations were carried out for an unblocked pentapeptide with the sequence Ser-Tyr-Pro-Tyr-Asp (SYPYD) with explicit consideration of the coupling between the conformation of the molecule and the ionization equilibria at a given pH. The available NMR experimental data indicate a high preference for the cis isomeric turn-like form of Tyr-Pro at intermediate pH (approximately 6) and a destabilization of the cis form at both high (approximately 9) and low (approximately 3) pH. In order to identify the source of the stability of the conformation of this pentapeptide as a function of pH, Monte Carlo simulations were used to generate an ensemble of low-energy conformations at different pH values (viz. 3, 6 and 9). The total free energy function used in these calculations includes terms that account for the solvation free energy and free energy of ionization. These terms are evaluated by means of a fast multigrid boundary element (MBE) method. In good qualitative agreement with the experiments, our results indicate that the Boltzmann averaged population of the cis isomeric form of the pentapeptide has a maximum (45 %) at pH 6 and is significantly smaller (25 % and 23 %) for higher and lower pH values, respectively, following the trend of the experimental data. Also, the degree of charge for the lowest-energy conformations, as well as the contribution of electrostatic interactions to the stability of the preferred conformations, vary widely at the different pH values. Different kinds of packing of the aromatic side-chains of Tyr2 and Tyr4 against the proline ring are observed at different pH values, indicating that their contribution to the stability of the low-energy conformations is also pH-dependent. In summary, our results provide a basis for discussing the nature of the interactions that stabilize turn-like conformations of the peptide SYPYD as a function of pH.

Algorithms↗

PubMed1998

Protein folding in the landscape perspective: chevron plots and non-Arrhenius kinetics.

We use two simple models and the energy landscape perspective to study protein folding kinetics. A major challenge has been to use the landscape perspective to interpret experimental data, which requires ensemble averaging over the microscopic trajectories usually observed in such models. Here, because of the simplicity of the model, this can be achieved. The kinetics of protein folding falls into two classes: multiple-exponential and two-state (single-exponential) kinetics. Experiments show that two-state relaxation times have "chevron plot" dependences on denaturant and non-Arrhenius dependences on temperature. We find that HP and HP+ models can account for these behaviors. The HP model often gives bumpy landscapes with many kinetic traps and multiple-exponential behavior, whereas the HP+ model gives more smooth funnels and two-state behavior. Multiple-exponential kinetics often involves fast collapse into kinetic traps and slower barrier climbing out of the traps. Two-state kinetics often involves entropic barriers where conformational searching limits the folding speed. Transition states and activation barriers need not define a single conformation; they can involve a broad ensemble of the conformations searched on the way to the native state. We find that unfolding is not always a direct reversal of the folding process.

Animals↗

PubMed2006

Flexibility and adaptability in binding of E. coli cytidine repressor to different operators suggests a role in differential gene regulation.

Interactions between DNA-bound transcription factors CytR and CRP regulate the promoters of the Escherichia coli CytR regulon. A distinctive feature of the palindromic CytR operators is highly variable length central spacers (0-9 bp). Previously we demonstrated distinct modes of CytR binding to operators that differ in spacer length. These different modes are characterized by opposite enthalpic and entropic contributions at 25 degrees C. Of particular note were radically different negative DeltaCp values suggesting variable contribution from coupled protein folding and/or DNA structural transitions. We proposed that the CytR DNA binding-domain adopts either a more rigid or flexible DNA-bound conformation in response to the different spacer lengths. More recently, similar effects were shown to contribute to discrimination between operator and non-specific DNA binding by LacR, a CytR homolog. Here we have extended the thermodynamic analysis to the remaining natural CytR operators plus a set of synthetic operators designed to isolate spacing as the single variable. The thermodynamic results show a broad and monotonic range of effects that are primarily dependent on spacer length. The magnitude of effects suggests participation by more than the DNA-binding domain. 15N HSQC NMR and CD spectral analyses were employed to characterize the structural basis for these effects. The results indicate that while CytR forms a well-ordered structure in solution, it is highly dynamic. We propose a model in which a large ensemble of native state conformations narrows upon binding, to an extent governed by operator spacing. This in turn is expected to constrain intermolecular interactions in the CytR-CRP-DNA complex, thus generating operator-specific effects on repression and induction of transcription.

Cytidine↗

PubMed2002

A comparison of methods for calculating NMR cross-relaxation rates (NOESY and ROESY intensities) in small peptides.

Three methods for calculating nuclear magnetic resonance cross-relaxation rates from molecular dynamics simulations of small flexible molecules have been compared in terms of their ability to reproduce relaxation data obtained experimentally and to produce consistent descriptions of the system. The importance of the accuracy of the simulation versus the amount of sampling of phase space has also been assessed by comparing different length simulations performed with different time step schemes. A nine-residue peptide from the protein HPr of E. coli was used as a test system. The work shows that, in this case, single conformations or a limited ensemble of configurations are insufficient to properly describe the behavior of the peptide and that different approaches to incorporate molecular motions lead to significant differences in the cross-relaxation rates calculated. The correlation between the cross-relaxation rates calculated from simulations performed with different time step schemes was high and increased with increasing simulation length indicating that the extent of sampling is more important than the details of the atomic motion.

Bacterial Proteins↗

PubMed2006

Ligation-state hydrogen exchange: coupled binding and folding equilibria in ribonuclease P protein.

Bacillus subtilis ribonuclease P protein (P protein) is predominantly unfolded (D) at physiological pH and low ionic strength; however, small molecule anionic ligands (e.g., sulfate) directly bind to and stabilize the folded state (NL2). Because the D + 2L <--> NL2 transition is experimentally two-state, high-energy states such as the singly bound, folded species (NL) and the unliganded folded species (N) are generally difficult to detect at equilibrium. To study the conformational properties of these ensembles, NMR-detected amide hydrogen exchange (HX) rates of P protein were measured at four sulfate (i.e., ligand) concentrations, a method we denote "ligation-state hydrogen exchange". The ligand concentration dependence of the HX rate of 47 residues was fit to a model with four possible HX pathways, corresponding to the local and/or global opening reactions from NL2 and NL, the local opening of N, and the global opening of N to D. Data analysis permits the calculation of the residue-specific free energy of opening from each ensemble as well as the fractional amide HX flux through each pathway. Results indicate that the predominant route of HX is through the NL and N states, which represent only 0.45% and 0.0005% of the total protein population in 20 mM sodium sulfate, respectively. Despite the low population of N, a region of protected amides was identified. Therefore, exchange through unliganded forms must be accounted for prior to the interpretation of HX-based protein-interaction studies. We offer a simple test to determine if HX occurs through the liganded or unliganded form.

Algorithms↗

PubMed1998

Clustering of low-energy conformations near the native structures of small proteins.

Recent experimental studies of the denatured state and theoretical analyses of the folding landscape suggest that there are a large multiplicity of low-energy, partially folded conformations near the native state. In this report, we describe a strategy for predicting protein structure based on the working hypothesis that there are a greater number of low-energy conformations surrounding the correct fold than there are surrounding low-energy incorrect folds. To test this idea, 12 ensembles of 500 to 1,000 low-energy structures for 10 small proteins were analyzed by calculating the rms deviation of the Calpha coordinates between each conformation and every other conformation in the ensemble. In all 12 cases, the conformation with the greatest number of conformations within 4-A rms deviation was closer to the native structure than were the majority of conformations in the ensemble, and in most cases it was among the closest 1 to 5%. These results suggest that, to fold efficiently and retain robustness to changes in amino acid sequence, proteins may have evolved a native structure situated within a broad basin of low-energy conformations, a feature which could facilitate the prediction of protein structure at low resolution.

Computer Simulation↗

PubMed2006

Perturbation theory of Phi-value analysis of two-state protein folding: relation between p fold and Phi values.

In protein folding, the transition state ensemble is defined as the set of conformations with p(fold)=12, where the p(fold) of a conformation is the probability that starting from this conformation the protein folds before it unfolds. Experimentally, this ensemble is probed by the Phi-value analysis, where Phi is the ratio of the changes in the logarithms of the folding rate and the equilibrium constant when the system is perturbed by a mutation. We show that for a two-state protein the Phi value can be expressed in terms of the perturbation and only the first two eigenfunctions of the evolution operator (e.g., a rate matrix) of the wild-type protein. The first eigenfunction is the equilibrium probability distribution while the second is proportional to p(fold), thus establishing a formal relation between p(fold) and Phi values. In addition to providing insight into the theoretical foundation of the Phi-value analysis, our results may prove practically useful in performing such analyses within the framework of models containing a large number of states.

Models, Biological↗

PubMed2002

Ab initio prediction of helical segments in polypeptides.

An ab initio method has been developed to predict helix formation for polypeptides. The approach relies on the systematic analysis of overlapping oligopeptides to determine the helical propensity for individual residues. Detailed atomistic level modeling, including entropic contributions, and solvation/ionization energies calculated through the solution of the Poisson-Boltzmann equation, is utilized. The calculation of probabilities for helix formation is based on the generation of ensembles of low energy conformers. The approach, which is easily amenable to parallelization, is shown to perform very well for several benchmark polypeptide systems, including the bovine pancreatic trypsin inhibitor, the immunoglobulin binding domain of protein G, the chymotrypsin inhibitor 2, the R69 N-terminal domain of phage 434 repressor, and the wheat germ agglutinin.

Animals↗

PubMed2001

Detection of two partially structured species in the folding process of the amyloidogenic protein beta 2-microglobulin.

beta 2-Microglobulin is a small, major histocompatibility complex class I-associated protein that undergoes aggregation and accumulates as amyloid deposits in human tissues as a consequence of long-term haemodialysis. The folding process of this amyloidogenic protein has been studied in vitro by diluting the guanidine hydrochloride-denatured protein in refolding buffer at pH 7.4 and monitoring the folding process by means of a number of spectroscopic probes that allow the native structure of the protein to be detected as it develops. These techniques include fluorescence spectroscopy, far and near-UV circular dichroism, 8-anilino-1-naphthalenesulfonic acid binding and double jump assays. All spectroscopic probes indicate that a significant amount of structure forms within the dead-time of stopped-flow measurements (<5 ms). The folding reaction goes to completion through a fast phase followed by a slow phase, whose rate constants are ca 5.1 and 0.0030 s(-1) in water, respectively. Unfolding-folding double jump experiments, together with the use of peptidyl prolyl isomerase, reveal that the slow phase of folding of beta 2-microglobulin is not fundamentally determined by cis/trans isomerisation of X-Pro peptide bonds. Other folding-unfolding double jump experiments also suggest that the fast and slow phases of folding are not related to independent folding of different populations of protein molecules. Rather, we provide evidence for a sequential mechanism of folding where denatured beta 2-microglobulin collapses to an ensemble of partially folded conformations (I(1)) which fold subsequently to a more highly structured species (I(2)) and, finally, attain the native state. The partially folded species I(2) appears to be closely similar to previously studied amyloidogenic forms of beta 2-microglobulin, such as those adopted by the protein at mildly acid pH values and by a variant with six residues deleted at the N terminus. Since amyloid formation in vivo originates from partial denaturation of beta 2-microglobulin under conditions favouring the folding process, the long-lived, partially structured species detected here might be significantly populated under some physiological conditions and hence might play an important role in the process of amyloid formation.

Amyloidosis↗

PubMed2002

Dynamic docking and electron transfer between myoglobin and cytochrome b(5).

The interaction of trypsin-digested bovine cytochrome b(5) (cyt b(5)) with horse heart myoglobin (Mb) and the interprotein electron transfer (ET) between these redox partners have been studied to gain better understanding of ET processes between weakly bound protein partners. The bimolecular rate constant ( k(2)) for photo-induced ET between zinc-substituted Mb (ZnMb) and cyt b(5) decreases with increasing ionic strength, consistent with the predominantly electrostatic character of this complex. The formation of a protein-protein complex has been confirmed and the binding affinities of metMb and ZnMb for cyt b(5) have been measured by two techniques: (1)H NMR titrations at pH 6.0 give binding constants of K(a) approximately (1.0+/-0.1)x10(3) M(-1) for metMb and K(a) approximately (0.75+/-0.1)x10(3) M(-1) for ZnMb; isothermal calorimetry gives K(a) approximately (0.35+/-0.1)x10(3) M(-1) for ZnMb. Brownian dynamic (BD) simulations show that cyt b(5) binds over a broad surface of Mb that includes its heme edge. The experimental results are described in terms of a dynamic docking model which proposes that Mb binds cyt b(5) in a large ensemble of protein binding conformations, not one or a few dominant ones, but that only a small subset are ET reactive. Aided by the BD simulations, this model explains why k(2) decreases with increasing pH: increasing pH not only weakens the binding affinity but also reduces the number of binding conformations with high ET reactivity.

Animals↗

PubMed2005

Cytocompatibility of self-assembled beta-hairpin peptide hydrogel surfaces.

MAX1 is a 20 amino acid peptide that undergoes triggered self-assembly to form a rigid hydrogel. When dissolved in aqueous solutions, this peptide exists in an ensemble of random coil conformations rendering it fully soluble. The addition of an exogenous stimulus results in peptide folding into beta-hairpin conformation. This folded structure undergoes rapid assembly into a highly crosslinked hydrogel network. DMEM cell culture media is one stimulus able to initiate folding and consequent self-assembly of MAX1. The cytocompatibility of this gel towards NIH 3T3 murine fibroblasts is demonstrated. Gels were shown to be non-toxic to the fibroblast cells. MAX1 hydrogels also foster the ability of the cells to attach to the hydrogel scaffold in the absence or presence of serum proteins. Additionally MAX1 hydrogels were able to support fibroblast proliferation to confluency with little effect on the rheological properties of the scaffold. MAX1 hydrogels meet the preliminary mechanical and cytocompatibiltiy requirements of a tissue engineering scaffold.

Animals↗