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 577 records · Page 32Linked to original sources

PubMed2004

Protein conformational transitions coupled to binding in molecular recognition of unstructured proteins: hierarchy of structural loss from all-atom Monte Carlo simulations of p27Kip1 unfolding-unbinding and structural determinants of the binding mechanism.

Conformational transitions coupled to binding are studied for the p27(Kip1) protein which undergoes a functional disorder-to-order folding transition during tertiary complex formation with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) binary complex. Temperature-induced Monte Carlo simulations of p27(Kip1) unfolding-unbinding carried out from the crystal structure of the tertiary complex have revealed a systematic trend in the hierarchy of structural loss for p27(Kip1) and a considerable difference in mobility of p27(Kip1) secondary structure elements. The most persistent interactions of p27(Kip1) at the intermolecular interface during unfolding-unbinding simulations are formed by beta-hairpin and beta-strand that on average maintain their structural integrity considerably longer than other p27(Kip1) elements. We have found that the ensemble of unfolded p27(Kip1) conformations is characterized by transitions between mostly unbound, collapsed conformations and entropically favorable p27(Kip1) conformations, which are weakly bound to the cyclin A side of the binary complex. The results of this study are consistent with the experimental evidence pointing to this region of the intermolecular interface as a potential initiation docking site during binding reaction and may reconcile conflicting experimental hypotheses on the recognition of substrate recruitment motifs.

Cell Cycle Proteins↗

PubMed2000

Correlation between changes in nuclear magnetic resonance order parameters and conformational entropy: molecular dynamics simulations of native and denatured staphylococcal nuclease.

Recent work has suggested that changes in NMR order parameters may quantitatively reflect changes in the conformational entropy of a protein ensemble. The extent of the mathematical relationship between local entropy changes as seen by NMR order parameters and the full protein entropy change is a complex issue. As a step towards a fuller understanding of this problem, molecular dynamics calculations of both native and denatured staphylococcal nuclease were performed. The N-H bond vector motion, in both explicit and implicit solvent, was analyzed to estimate local and global entropy changes. The calculated N-H bond vector order parameters from simulation agreed on average with experimental values for both native and denatured structures. However, the inverted-U profile of order parameters versus residue number observed experimentally for denatured nuclease was only partially reproduced by simulation of compact denatured structures. Comparisons made across the full set of simulations revealed a correlation between the N-H order parameter-based conformational entropy change and the total quasiharmonic-based conformational entropy change between the native and denatured structures. The calculations showed that about 25% of the total entropy change was reflected by changes in simulated S2 values. This result suggests that NMR-derived order parameters may be used to provide a reasonable estimate of the total conformational entropy change on protein folding.

Magnetic Resonance Spectroscopy↗

PubMed2005

Conformation coupled enzyme catalysis: single-molecule and transient kinetics investigation of dihydrofolate reductase.

Ensemble kinetics and single-molecule fluorescence microscopy were used to study conformational transitions associated with enzyme catalysis by dihydrofolate reductase (DHFR). The active site loop of DHFR was labeled with a fluorescence quencher, QSY35, at amino acid position 17, and the fluorescent probe, Alexa555, at amino acid 37, by introducing cysteines at these sites with site-specific mutagenesis. The distance between the probes was such that approximately 50% fluorescence resonance energy transfer (FRET) occurred. The double-labeled enzyme retained essentially full catalytic activity, and stopped-flow studies of both the forward and reverse reactions revealed that the distance between probes increased prior to hydride transfer. A fluctuation in fluorescence intensity of single molecules of DHFR was observed in an equilibrium mixture of substrates but not in their absence. Ensemble rate constants were derived from the distributions of lifetimes observed and attributed to a reversible conformational change. Studies were carried out with both NADPH and NADPD as substrates, with no measurable isotope effect. Similar studies with a G121V mutant DHFR resulted in smaller rate constants. This mutant DHFR has reduced catalytic activity, so that the collective data for the conformational change suggest that the conformational change being observed is associated with catalysis and probably represents a conformational change prior to hydride transfer. If the change in fluorescence is attributed to a change in FRET, the distance change associated with the conformational change is approximately 1-2 A. These results are correlated with other measurements related to conformation coupled catalysis.

Algorithms↗

PubMed2003

Preferred conformations of cyclic Ac-Cys-Pro-Xaa-Cys-NHMe peptides: a model for chain reversal and active site of disulfide oxidoreductase.

The conformational study on cyclic Ac-Cys-Pro-Xaa-Cys-NHMe (Ac-CPXC-NHMe; X=Ala, Val, Leu, Aib, Gly, His, Phe, Tyr, Asn and Ser) peptides has been carried out using the Empirical Conformational Energy Program for Peptides, version 3 (ECEPP/3) force field and the hydration shell model in the unhydrated and hydrated states. This work has been undertaken to investigate structural implications of the CPXC sequence as the chain reversal for the initiation of protein folding and as the motif for active site of disulfide oxidoreductases. The backbone conformation DAAA is commonly the most feasible for cyclic CPXC peptides in the hydrated state, which has a type I beta-turn at the Pro-Xaa sequence. The proline residue and the hydrogen bond between backbones of two cystines as well as the formation of disulfide bond appear to play a role in stabilizing this preferred conformation of cyclic CPXC peptides. However, the distributions of backbone conformations and beta-turns may indicate that the cyclic CPXC peptide seems to exist as an ensemble of beta-turns and coiled conformations in aqueous solution. The intrinsic stability of the cyclic CPXC motif itself for the active conformation seems to play a role in determining electrochemical properties of disulfide oxidoreductases.

Amino Acid Sequence↗

PubMed2000

A novel energy-based stochastic method for positioning polar protons in protein structures from X-rays.

A novel automated method for the optimal placement of polar hydrogens in a protein structure is presented. The algorithm adds initially, to a protein data bank file of the protein, nonrotatable hydrogens such as peptide backbone hydrogens according to geometric considerations. Then, water protons and polar side chain protons of lysine, serine, threonine, tyrosine, aspartic acid, glutamic acid, and the C and N termini of a protein are added according to energy considerations. A unique stochastic approach has been developed to overcome a combinatorial explosion in the search for the lowest energy structure. First, the system is divided into ensembles. Each ensemble is treated separately: N conformations are sampled at random, their energies computed, whereas common components of high-energy combinations are gathered on one hand, and low-energy combinations on the other. Components that yield only high-energy conformations and do not contribute to any low energies are excluded. This is reiterated while the total amount of combinations is decreased along the iterative process. When the total number of combinations is lower than a user defined threshold, all remaining combinations are evaluated by exhaustive search. Energy evaluations use nonbonding energy expressions alone. The program was tested on five high-resolution crystal structures: bovine pancreatic trypsin inhibitor (Brookhaven Protein Data Bank file 5PTI), RNase-A (5RSA), trypsin (1NTP), and carbon monoxymyoglobin (2MB5), for which neutron diffraction structures are available, as well as phosphate binding protein (1IXH) for which very high resolution X-ray crystallography was used. The low RMS values prove the efficiency of this algorithm as a tool for positioning protons in proteins. It may be used for other biological structures.

Algorithms↗

PubMed2001

Rapid boundary element solvation electrostatics calculations in folding simulations: successful folding of a 23-residue peptide.

Solvation effects play a profound role in the energetics of protein folding. While a continuum dielectric model of solvation may provide a sufficiently accurate estimate of the solvation effects, until now this model was too computationally expensive and unstable for folding simulations. Here we proposed a fast yet accurate and robust implementation of the boundary element solution of the Poisson equation, the REBEL algorithm. Using our earlier double-energy scheme, we included for the first time the mathematically rigorous continuous REBEL solvation term in our Biased Probability Monte Carlo (BPMC) simulations of the peptide folding. The free energy of a 23-residue beta beta alpha-peptide was then globally optimized with and without the solvation electrostatics contribution. An ensemble of beta beta alpha conformations was found at and near the global minimum of the energy function with the REBEL electrostatic solvation term. Much poorer correspondence to the native solution structure was found in the "control" simulations with a traditional method to account for solvation via a distance-dependent dielectric constant. Each simulation took less than 40 h (21 h without electrostatic solvation calculation) on a single Alpha 677 MHz CPU and involved more than 40,000 solvation energy evaluations. This work demonstrates for the first time that such a simulation can be performed in a realistic time frame. The proposed procedure may eliminate the energy evaluation accuracy bottleneck in folding simulations.

Algorithms↗

PubMed2001

The role of backbone motions in ligand binding to the c-Src SH3 domain.

The Src homology 3 (SH3) domain of pp60(c-src) (Src) plays dual roles in signal transduction, through stabilizing the repressed form of the Src kinase and through mediating the formation of activated signaling complexes. Transition of the Src SH3 domain between a variety of binding partners during progression through the cell cycle requires adjustment of a delicate free energy balance. Although numerous structural and functional studies of SH3 have provided an in-depth understanding of structural determinants for binding, the origins of binding energy in SH3-ligand interactions are not fully understood. Considering only the protein-ligand interface, the observed favorable change in standard enthalpy (DeltaH=-9.1 kcal/mol) and unfavorable change in standard entropy (TDeltaS=-2.7 kcal/mol) upon binding the proline-rich ligand RLP2 (RALPPLPRY) are inconsistent with the predominantly hydrophobic interaction surface. To investigate possible origins of ligand binding energy, backbone dynamics of free and RLP2-bound SH3 were performed via (15)N NMR relaxation and hydrogen-deuterium (H/(2)H) exchange measurements. On the ps-ns time scale, assuming uncorrelated motions, ligand binding results in a significant reduction in backbone entropy (-1.5(+/-0.6) kcal/mol). Binding also suppresses motions on the micros-ms time scale, which may additionally contribute to an unfavorable change in entropy. A large increase in protection from H/(2)H exchange is observed upon ligand binding, providing evidence for entropy loss due to motions on longer time scales, and supporting the notion that stabilization of pre-existing conformations within a native state ensemble is a fundamental paradigm for ligand binding. Observed changes in motion on all three time scales occur at locations both near and remote from the protein-ligand interface. The propagation of ligand binding interactions across the SH3 domain has potential consequences in target selection through altering both free energy and geometry in intact Src, and suggests that looking beyond the protein-ligand interface is essential in understanding ligand binding energetics.

Animals↗

PubMed2005

Scaling behavior and structure of denatured proteins.

An ensemble of random-coil conformations with no persistent structures has long been accepted as the classical model of denatured proteins due to its consistency with the experimentally determined scaling of protein sizes. However, recent NMR spectroscopy studies on proteins at high chemical denaturant concentrations suggest the presence of significant amounts of native-like structures, in contrast to the classical random-coil picture. To reconcile these seemingly controversial observations, we examine thermally denatured states of experimentally characterized proteins by using molecular dynamics simulations. For all studied proteins, we find that denatured states indeed have strong local conformational bias toward native states while a random-coil power law scaling of protein sizes is preserved. In addition, we explain why experimentally determined size of the protein creatine kinase does not follow general scaling. In simulations, we observe that this protein exhibits a stable intermediate state, the size of which is consistent with the reported experimental observation.

Computer Simulation↗

PubMed2003

Stability and folding kinetics of a ubiquitin mutant with a strong propensity for nonnative beta-hairpin conformation in the unfolded state.

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.

Amino Acid Sequence↗

PubMed2003

Role of entropy in increased rates of intramolecular reactions.

The rates of intramolecular condensation of a series of monoesters of dicarboxylic acids have been shown to be highly dependent on the nature of the intervening groups. To understand the origin of this effect, we estimated DeltaS(NAC,S), the entropy difference between the ensemble of accessible ground state conformers and a single ground state conformer having transition-state-like geometry. DeltaS(NAC,S) differs from the activation entropy for the reaction by DeltaS(TS,NAC), the difference in vibrational entropy between the selected ground state conformer and the transition state. The estimated values of DeltaS(NAC,S) correlate well (R(2) = 0.96 and 0.73 using dielectric constant values of 80 and 1, respectively) with experimentally determined reaction rate constants. Normal-mode analysis performed on minimized ground state conformations of each molecule suggests that the change in vibrational entropy makes only a small contribution to the total activation entropy. These results indicate that the conformational entropy difference between the transition and the ground states contributes significantly to the free energy of activation.

Journal Article↗

PubMed2006

Dynamics of Asp23-Lys28 salt-bridge formation in Abeta10-35 monomers.

In the amyloid fibrils formed from long fragments of the amyloid beta-protein (Abeta-protein), the monomers are arranged in parallel and lie perpendicular to the fibril axis. The structure of the monomers satisfies the amyloid self-organization principle; namely, the low free energy state of the monomer maximizes the number of intra- and interpeptide contacts and salt bridges. The formation of the intramolecular salt bridge between Asp(D)23 and Lys(K)28 ensures that unpaired charges are not buried in the low-dielectric interior. We have investigated, using all-atom molecular dynamics simulations in explicit water, whether the D23-K28 interaction forms spontaneously in the isolated Abeta10-35 monomer. To validate the simulation protocol, we show, using five independent trajectories spanning a total of 100 ns, that the pKa values of the titratable groups are in good agreement with experimental measurements. The computed free energy disconnectvity graph shows that broadly the ensemble of compact random coil conformations can be clustered into four basins that are separated by free energy barriers ranging from 0.3 to 2.7 kcal/mol. There is significant residual structure in the conformation of the peptide in each of the basins. Due to the desolvation penalty, the structural motif with a stable turn involving the residues VGSN and a preformed D23-K28 contact is a minor component of the simulated structures. The extent of solvation of the peptides in the four basins varies greatly, which underscores the dynamical fluctuations in the monomer. Our results suggest that the early event in the oligomerization process must be the expulsion of discrete water molecules that facilitates the formation of interpeptide-interaction-driven stable structures with an intramolecular D23-K28 salt bridge and an intact VGSN turn.

Amyloid beta-Peptides↗

PubMed1998

Osmolyte-driven contraction of a random coil protein.

The Stokes radius characteristics of reduced and carboxamidated ribonuclease A (RCAM RNase) were determined for transfer of this "random coil" protein from water to 1 M concentrations of the naturally occurring protecting osmolytes trimethylamine N-oxide, sarcosine, sucrose, and proline and the nonprotecting osmolyte urea. The denatured ensemble of RCAM RNase expands in urea and contracts in protecting osmolytes to extents proportional to the transfer Gibbs energy of the protein from water to osmolyte. This proportionality suggests that the sum of the transfer Gibbs energies of individual parts of the protein is responsible for the dimensional changes in the denatured ensemble. The dominant term in the transfer Gibbs energy of RCAM RNase from water to protecting osmolytes is the unfavorable interaction of the osmolyte with the peptide backbone, whereas the favorable interaction of urea with the backbone dominates in RCAM RNase transfer to urea. The side chains collectively favor transfer to the osmolytes, with some protecting osmolytes solubilizing hydrophobic side chains as well as urea does, a result suggesting there is nothing special about the ability of urea to solubilize hydrophobic groups. Protecting osmolytes stabilize proteins by raising the chemical potential of the denatured ensemble, and the uniform thermodynamic force acting on the peptide backbone causes the collateral effect of contracting the denatured ensemble. The contraction decreases the conformational entropy of the denatured state while increasing the density of hydrophobic groups, two effects that also contribute to the ability of protecting osmolytes to force proteins to fold.

Models, Chemical↗

PubMed2006

The electron transfer complex between cytochrome c552 and the CuA domain of the Thermus thermophilus ba3 oxidase. A combined NMR and computational approach.

The structural analysis of the redox complex between the soluble cytochrome c552 and the membrane-integral cytochrome ba3 oxidase of Thermus thermophilus is complicated by the transient nature of this protein-protein interaction. Using NMR-based chemical shift perturbation mapping, however, we identified the contact regions between cytochrome c552 and the CuA domain, the fully functional water-soluble fragment of subunit II of the ba3 oxidase. First we determined the complete backbone resonance assignments of both proteins for each redox state. Subsequently, two-dimensional [15N,1H]TROSY spectra recorded for each redox partner both in free and complexed state indicated those surface residues affected by complex formation between the two proteins. This chemical shift analysis performed for both redox states provided a topological description of the contact surface on each partner molecule. Remarkably, very pronounced indirect effects, which were observed on the back side of the heme cleft only in the reduced state, suggested that alterations of the electron distribution in the porphyrin ring due to formation of the protein-protein complex are apparently sensed even beyond the heme propionate groups. The contact residues of each redox partner, as derived from the chemical shift perturbation mapping, were employed for a protein-protein docking calculation that provided a structure ensemble of 10 closely related conformers representing the complex between cytochrome c552 and the CuA domain. Based on these structures, the electron transfer pathway from the heme of cytochrome c552 to the CuA center of the ba3 oxidase has been predicted.

Cytochrome b Group↗

PubMed2004

Helix formation and the unfolded state of a 52-residue helical protein.

A growing class of proteins in biological processes has been found to be unfolded on isolation under normal solution conditions. We have used NMR spectroscopy to characterize the structural and dynamic properties of the unfolded and partially folded states of a 52-residue alanine-rich protein (Ala-14) at temperatures from -5 degrees C to 40 degrees C. At 40 degrees C, alanine residues in Ala-14 adopt phi and psi angles, consistent with a significant ensemble population of polyproline II conformation. Analysis of relaxation rates in the protein reveals that a series of residues, Gln 35-Ala 36-Ala 37-Lys 38-Asp 39-Asp 40-Ala 41-Ala 42, displays slow motional dynamics at both -5 degrees C and 40 degrees C. Temperature-dependent chemical shift changes indicate that this region is the site of helix initiation. The remaining N-terminal residues become increasingly dynamic as they extend from the nucleation site. The C terminus remains dynamic and changes less with temperature, indicating it is relatively unstructured. Ala-14 provides a high-resolution portrait of the unfolded state and the process of helix nucleation and propagation in the absence of tertiary contacts, information that bears on early events in protein folding.

Alanine↗

PubMed2005

Unfolding of Green Fluorescent Protein mut2 in wet nanoporous silica gels.

Many of the effects exerted on protein structure, stability, and dynamics by molecular crowding and confinement in the cellular environment can be mimicked by encapsulation in polymeric matrices. We have compared the stability and unfolding kinetics of a highly fluorescent mutant of Green Fluorescent Protein, GFPmut2, in solution and in wet, nanoporous silica gels. In the absence of denaturant, encapsulation does not induce any observable change in the circular dichroism and fluorescence emission spectra of GFPmut2. In solution, the unfolding induced by guanidinium chloride is well described by a thermodynamic and kinetic two-state process. In the gel, biphasic unfolding kinetics reveal that at least two alternative conformations of the native protein are significantly populated. The relative rates for the unfolding of each conformer differ by almost two orders of magnitude. The slower rate, once extrapolated to native solvent conditions, superimposes to that of the single unfolding phase observed in solution. Differences in the dependence on denaturant concentration are consistent with restrictions opposed by the gel to possibly expanded transition states and to the conformational entropy of the denatured ensemble. The observed behavior highlights the significance of investigating protein function and stability in different environments to uncover structural and dynamic properties that can escape detection in dilute solution, but might be relevant for proteins in vivo.

Circular Dichroism↗

PubMed2004

Atomic force microscopy study of DNA deposited on poly L-ornithine-coated mica.

Analyses of individual biomolecules, like DNA, or DNA-protein complexes, via atomic force microscopy, require 'gentle' methods to immobilize DNA on surfaces, which allow the ensemble of molecules to adopt conformations dictated primarily by their physical characteristics, and which possibly permit the use of a wide selection of buffers. We show that poly-L-ornithine-coated mica is a good substrate for fast, reliable deposition of DNA for wet or dry imaging. The surface firmly secures DNA, which retains the B-form helical rise (0.34 nm bp(-1)). The conformations of DNA that result are reminiscent of three-dimensional random coils projected on to a plane. The contrast is good, especially in solution, and buffers with physiological concentrations of salt with or without divalent cations may be used. This is important for comparison of scanning probe microscopy results with those obtained by different techniques.

Aluminum Silicates↗

PubMed2007

Estimation of folding probabilities and phi values from molecular dynamics simulations of reversible Peptide folding.

Molecular dynamics simulations with an implicit model of the solvent have allowed to investigate the reversible folding of structured peptides. For a 20-residue antiparallel beta-sheet peptide, the simulation results have revealed multiple folding pathways. Moreover, the conformational heterogeneity of the denatured state has been shown to originate from high enthalpy, high entropy basins with fluctuating non-native secondary structure, as well as low enthalpy, low entropy traps. An efficient and simple approach to estimate folding probabilities from molecular dynamics simulations has allowed to isolate conformations in the transition state ensemble and to evaluate Phi values, i.e., the effects of mutations on the folding kinetics and thermodynamic stability. These molecular dynamics studies have provided evidence that, if interpreted by neglecting the non-native interactions, Phi values overestimate the amount of native-like structure in the transition state.

Computer Simulation↗

PubMed2001

Comparison of aqueous molecular dynamics with NMR relaxation and residual dipolar couplings favors internal motion in a mannose oligosaccharide.

An investigation has been performed to assess how aqueous dynamical simulations of flexible molecules can be compared against NMR data. The methodology compares state-of-the-art NMR data (residual dipolar coupling, NOESY, and (13)C relaxation) to molecular dynamics simulations in water over several nanoseconds. In contrast to many previous applications of residual dipolar coupling in structure investigations of biomolecules, the approach described here uses molecular dynamics simulations to provide a dynamic representation of the molecule. A mannose pentasaccharide, alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->2)-D-Manp, was chosen as the model compound for this study. The presence of alpha-linked mannan is common to many glycopeptides, and therefore an understanding of the structure and the dynamics of this molecule is of both chemical and biological importance. This paper sets out to address the following questions. (1) Are the single structures which have been used to interpret residual dipolar couplings a useful representation of this molecule? (2) If dynamic flexibility is included in a representation of the molecule, can relaxation and residual dipolar coupling data then be simultaneously satisfied? (3) Do aqueous molecular dynamics simulations provide a reasonable representation of the dynamics present in the molecule and its interaction with water? In summary, two aqueous molecular dynamics simulations, each of 20 ns, were computed. They were started from two distant conformations and both converged to one flexible ensemble. The measured residual dipolar couplings were in agreement with predictions made by averaging the whole ensemble and from a specific single structure selected from the ensemble. However, the inclusion of internal motion was necessary to rationalize the relaxation data. Therefore, it is proposed that although residual dipolar couplings can be interpreted as a single-structure, this may not be a correct interpretation of molecular conformation in light of other experimental data. Second, the methodology described here shows that the ensembles from aqueous molecular dynamics can be effectively tested against experimental data sets. In the simulation, significant conformational motion was observed at each of the linkages, and no evidence for intramolecular hydrogen bonds at either alpha(1-->2) or alpha(1-->3) linkages was found. This is in contrast to simulations of other linkages, such as beta(1-->4), which are often predicted to maintain intramolecular hydrogen bonds and are coincidentally predicted to have less conformational freedom in solution.

Carbohydrate Sequence↗