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 469 records · Page 26Linked to original sources

PubMed2003

The Ac-RGD-NH2 peptide as a probe of slow conformational exchange of short linear peptides in DMSO.

According to general belief, the conformational information on short linear peptides in solution derived at ambient temperature from NMR spectrometry represents a population-weighted average over all members of an ensemble of rapidly interconverting conformations. Usually the search for discrete conformations is concentrated at low temperatures especially when sharp NMR resonances are detected at room temperature. Using the peptide Ac-RGD-NH(2) (Ac-Arg-Gly-Asp-NH(2), Ac: acetyl) as a model system and following a new approach, we have been able to demonstrate that short linear peptides can adopt discrete conformational states in DMSO-d(6) (DMSO: dimethylsulfoxide) which vary in a way critically dependent on the reconstitution conditions used before their dissolution in DMSO-d(6). The conformers are stabilized by intramolecular hydrogen bonds, which persist at high temperatures and undergo a very slow exchange with their extended structures in the NMR chemical shift time scale. The reported findings provide clear evidence for the occurrence of solvent-induced conformational exchange and point to DMSO as a valuable medium for folding studies of short linear peptides.

Amino Acid Sequence↗

PubMed2005

Mapping long-range interactions in alpha-synuclein using spin-label NMR and ensemble molecular dynamics simulations.

The intrinsically disordered protein alpha-synuclein plays a key role in the pathogenesis of Parkinson's disease (PD). We show here that the native state of alpha-synuclein consists of a broad distribution of conformers with an ensemble-averaged hydrodynamic radius significantly smaller than that expected for a random coil structure. This partial condensation is driven by interactions between the highly charged C-terminus and a large hydrophobic central region of the protein sequence. We suggest that this structure could inhibit the formation of alpha-synuclein aggregates, which are thought to be the cytotoxic species responsible for neurodegeneration in PD.

Computer Simulation↗

PubMed2003

Giant DNA molecules exhibit on/off switching of transcriptional activity through conformational transition.

We found that the transcriptional activity of large DNAs (40 kbp) can be completely inhibited by adding condensing agents, spermine and poly(ethylene glycol), whereas under the same conditions short fragments (140 bp) still show active transcription. Fluorescence microscopic observations of large DNAs revealed clear correlation between the higher-order structure of templates and their transcriptional activity. The steep decrease in transcriptional activity leading to complete inhibition, or on/off switching, is interpreted in terms of conformational transition of the ensemble of DNA molecules.

DNA↗

PubMed2006

Analyzing the performance of conformational search programs on compound databases.

We have studied the sampling performance of conformational search programs using geometric and energetic criteria. Ideally, a conformational search algorithm should identify the largest possible number of low-energy structures (energy criterion) covering the widest possible range of molecular shapes (geometric criterion). Geometric analysis consisted in comparing the distribution of conformations within the generated ensembles by multidimensional scaling and by analysing the eigenvalue structure of the pairwise coordinate covariance matrices. The energetic comparison was carried out by assessing the energy distribution of conformers after minimizing them all using the same semi-empirical quantum mechanics optimization protocol. The present investigation focused on five conformational search programs: DGEOM, QXP, ROTATE, LMOD and OMEGA. We have applied these methodologies to a maximally diverse 604-compound subset of the LeadQuest library. The program LMOD performs best according to the energetic criterion, whereas a wider range of geometrically diverse conformations is sampled by the other programs, at the cost of higher median conformer energies. In terms of speed, OMEGA is fastest. We recommend the use of LMOD or OMEGA for high-quality conformational search applications.

Computer Simulation↗

PubMed1997

The energy landscape of a fast-folding protein mapped by Ala-->Gly substitutions.

A moderately stable protein with typical folding kinetics unfolds and refolds many times during its cellular lifetime. In monomeric lambda repressor this process is extremely rapid, with an average folded state lifetime of only 30 milliseconds. A thermostable variant of this protein (G46A/G48A) unfolds with the wild-type rate, but it folds in approximately 20 microseconds making it the fastest-folding protein yet observed. The effects of alanine to glycine substitutions on the folding and unfolding rate constants of the G46A/G48A variant, measured by dynamic NMR spectroscopy, indicate that the transition state is an ensemble comprised of a disperse range of conformations. This structural diversity in the transition state is consistent with the idea that folding chains are directed towards the native state by a smooth funnel-like conformational energy landscape. The kinetic data for the folding of monomeric lambda repressor can be understood by merging the new energy landscape view of folding with traditional models. This hybrid model incorporates the conformational diversity of denatured and transition state ensembles, a transition state activation energy, and the importance of intrinsic helical stabilities.

Alanine↗

PubMed2002

The role of heterologous receptors in McpB-mediated signalling in Bacillus subtilis chemotaxis.

Asparagine chemotaxis in Bacillus subtilis appears to involve two partially redundant adaptation mechanisms: a receptor methylation-independent process that operates at low attractant concentrations and a receptor methylation-dependent process that is required for optimal responses to high concentrations. In order to elucidate these processes, chemotactic responses were assessed for strains expressing methylation-defective mutations in the asparagine receptor, McpB, in which all 10 putative receptors (10del), five receptors (5del) or only the native copy of mcpB were deleted. This was done in both the presence and the absence of the methylesterase CheB. We found that: (i) only responses to high concentrations of asparagine were impaired; (ii) the presence of all heterologous receptors fully compensated for this defect, whereas responses progressively worsened as more receptors were taken away; (iii) methyl-group turnover occurred on heterologous receptors after the addition of asparagine, and these methylation changes were required for the restoration of normal swimming behaviour; (iv) in the absence of the methyleste-rase, the presence of heterologous receptors in some cases caused impaired chemotaxis; and (v) either a certain threshold number of receptors must be present to promote basal CheA activity, or one or more of the receptors missing in the 10del background (but present in the 5del background) is required for establishing basal CheA activity. Taken together, these findings suggest that many or all chemoreceptors work as an ensemble that constitutes a robust chemotaxis system. We propose that the ability of non-McpB receptors to compensate for the methylation-defective McpB mutations involves lateral transmission of the adapted conformational change across the ensemble.

Adaptation, Physiological↗

PubMed1988

Solution and membrane structure of enkephalins as studied by infrared spectroscopy.

The backbone conformation of the two opioid pentapeptides Leu5-enkephalin and Met5-enkephalin was studied by the technique of resolution-enhanced infrared spectroscopy. In aqueous solution, the conformation-sensitive amide I bands of the two peptides are identical. The positions of these bands are consistent with the view that in aqueous solution both enkephalins exist as an ensemble of largely unfolded conformers. Interaction of Leu5- and Met5-enkephalins with bilayer membranes of ditetradecylphosphatidylcholine results in a substantial refolding of the peptide backbones. The conformation stabilized by the membrane environment is a hydrogen-bonded turn structure. Conformational transitions in enkephalins induced by a lipid environment may play a role in the specific interactions between these hormones and their receptor sites.

Enkephalin, Leucine↗

PubMed1998

Vancomycin: conformational consequences of the sugar substituent.

High-resolution, three-dimensional structures of vancomycin and aglyco-vancomycin in DMSO were determined by nuclear magnetic resonance, metric matrix distance geometry, and molecular dynamics calculations. Conformational flexibility fast on the NMR time scale was examined by ensemble-based calculations which apply the experimentally derived restraints as an ensemble average. Two families of conformations of vancomycin, differing in the positioning of the vancosamine substituent, were observed. In contrast, the aglyco-vancomycin adopts only one conformation in solution. The conformations of vancomycin and the aglyco-vancomycin differ in the alignment of the amide protons which participate in the hydrogen-bonding network with the cell-wall precursor and orientation of the aromatic rings relative to the backbone. Therefore, the high-resolution structural characterization provides insight into a possible role of glycosylation on the activity of this important family of antibiotics.

Anti-Bacterial Agents↗

PubMed2001

Sensitivity of an empirical affinity scoring function to changes in receptor-ligand complex conformations.

A combination of empirical scoring and conformational sampling for ligand binding affinity prediction is examined. The behaviour of a scoring function with respect to the sensitivity to conformational changes is investigated using ensembles of structures generated by molecular dynamics simulation. The correlation between the calculated score and the coordinate deviation from the experimental structure is clear for the complex of arabinose with arabinose-binding protein, which is dominated by hydrogen bond interactions, while the score calculated for the hydrophobic complex between retinol and retinol binding protein is rather insensitive to ligand conformational changes. For typical ensembles of structures generated by molecular dynamics at 300 K, the variation of the calculated score is considerably smaller than that of the underlying molecular mechanics interaction energies.

Algorithms↗

PubMed1998

Solution structures of DNA.RNA hybrids with purine-rich and pyrimidine-rich strands: comparison with the homologous DNA and RNA duplexes.

The structures of d(GAAGAGAAGC).d(GCTTCTCTTC), d(GAAGAGAAGC). r(GCUUCUCUUC), r(GAAGAGAAGC).d(GCTTCTCTTC), and r(GAAGAGAAGC). r(GCUUCUCUUC) have been determined in solution from NMR data. Globally, the pure DNA and RNA duplexes were in the B and A forms, respectively. The two DNA.RNA hybrids were neither A nor B, but closer globally to the A than the B form. However, the thermodynamically less stable d(GAAGAGAAGC).r(GCUUCUCUUC) duplex has a significantly different conformation from r(GAAGAGAAGC). d(GCTTCTCTTC). Structures were calculated based on the NMR data, using restrained molecular dynamics. A new approach to the treatment of conformational averaging based on a prioriprobabilities has been used. The nucleotides were treated by fitting the scalar coupling data and NOE time courses to a two-state model comprising N and S sugar puckers each with a different glycosidic torsion angle, and the mole fraction of the S state. Restraint sets for different distributions of N and S states within molecules were constructed, such that each nucleotide was weighted in the ensemble according to the mole fractions (or a prioriprobabilities). The individual nucleotide conformations were strongly restrained, whereas the internucleotide restraints were set relatively loosely. Ensembles of conformations were generated and assessed by comparison of the NOEs calculated from ensemble-averaged relaxation matrices with the experimental NOEs. The ensemble averages accounted for the experimental data much better than any individual member, or for structures calculated assuming a single unique conformation. The two hybrids populated different degrees of conformational space. There was a general trend in minor and major groove widths in the order d(GAAGAGAAGC).d(GCTTCTCTTC), d(GAAGAGAAGC).r(GCUUCUCUUC), r(GAAGAGAAGC).d(GCTTCTCTTC), r(GAAGAGAAGC).r(GCUUCUCUUC) and a similar progression in global character from B-like to A-like structures. Furthermore, r(GAAGAGAAGC).d(GCTTCTCTTC) showed a greater dispersion of conformations in the ensemble than d(GAAGAGAAGC).r(GCUUCUCUUC), reflecting the greater flexibility of the sugars. If conformational averaging of the nucleotides is ignored, incorrect virtual structures are produced that nevertheless are able to satisfy a substantial fraction of the experimental data.

Crystallography, X-Ray↗

PubMed2003

Fast protein folding on downhill energy landscape.

Proteins fold in a time range of microseconds to minutes despite the large amount of possible conformers. Molecular dynamics simulations of a three-stranded antiparallel beta-sheet peptide (for a total of 12.6 microsec and 72 folding events) show that at the melting temperature the unfolded state ensemble contains many more conformers than those sampled during a folding event.

Algorithms↗

PubMed2006

Kinetic definition of protein folding transition state ensembles and reaction coordinates.

Using distributed molecular dynamics simulations we located four distinct folding transitions for a 39-residue betabetaalphabeta protein fold. To characterize the nature of each room temperature transition, we calculated the probability of transmission for 500 points along each free energy barrier. We introduced a method for determining transition states by employing the transmission probability, Ptrans, and determined which conformations were transition state ensemble members (Ptrans approximately 0.5). The transmission probability may be used to characterize the barrier in several ways. For example, we ran simulations at 82 degrees C, determined the change in Ptrans with temperature for all 2,000 conformations, and quantified Hammond behavior directly using Ptrans correlation. Additionally, we propose that diffusion along Ptrans may provide the configurational diffusion rate at the top of the barrier. Specifically, given a transition state conformation x0 with estimated Ptrans=0.5, we selected a large set of subsequent conformations from independent trajectories, each exactly a small time deltat after x0 (250 ps). Calculating Ptrans for the new trial conformations, we generated the P(Ptrans|deltat=250 ps) distribution that reflected diffusion. This approach provides a novel perspective on the diffusive nature of a protein folding transition and provides a framework for a quantitative study of activated relaxation kinetics.

Computer Simulation↗

PubMed2006

Strategies for dealing with conformational sampling in structural calculations of flexible or kinked transmembrane peptides.

Peptides corresponding to transmembrane (TM) segments from membrane proteins provide a potential route for the determination of membrane protein structure. We have determined that 2 functionally critical TM segments from the mammalian Na+/H+ exchanger display well converged structure in regions separated by break points. The flexibility of these break points results in conformational sampling in solution. A brief review of available NMR structures of helical membrane proteins demonstrates that there are a number of published structures showing similar properties. Such flexibility is likely indicative of kinks in the full-length protein. This minireview focuses on methods and protocols for NMR structure calculation and analysis of peptide structures under conditions of conformational sampling. The methods outlined allow the identification and analysis of structured peptides containing break points owing to conformational sampling and the differentiation between oligomerization and ensemble-averaged observation of multiple peptide conformations.

Mathematics↗

PubMed1999

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives.

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5'-GGGCGNNAGCCU-3' sequences commonly fold into hairpins with four-membered loops. The combinatorial search for acceptable conformations using the MC-SYM program was restricted to loop nucleotides and yielded roughly 1500 structures being compatible with a double-stranded stem. After energy minimization by the JUMNA program (without applying any experimental constraints), these structures converged into an ensemble of 74 different conformers including 26 structures which contained the sheared G-A base pair observed in experimental studies of GNRA tetraloops. Energetic analysis shows that inclusion of solvent electrostatic effects is critically important for the selection of conformers that agree with experimentally determined structures. The continuum model accounts for solvent polarization by means of the electrostatic reaction field. In the case of GNRA loop sequences, the contributions of the reaction field shift relative stabilities towards conformations showing most of the structural features derived from NMR studies. The agreement of computed conformations with the experimental structures of GAAA, GCAA, and GAGA tetraloops suggests that the continuum treatment of the solvent represents a definitive improvement over methods using simple damping models in electrostatic energy calculations. Application of the procedure described here to the evaluation of the relative stabilities of conformers resulting from searching the conformational space of RNA structural motifs provides some progress in (non-homology based) RNA 3D-structure prediction.

Base Sequence↗

PubMed1999

Validation of NMR side-chain conformations by packing calculations.

The precision and accuracy of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy depend on the completeness of input experimental data set. Typically, rather than a single structure, an ensemble of up to 20 equally representative conformers is generated and routinely deposited in the Protein Database. There are substantially more experimentally derived restraints available to define the main-chain coordinates than those of the side chains. Consequently, the side-chain conformations among the conformers are more variable and less well defined than those of the backbone. Even when a side chain is determined with high precision and is found to adopt very similar orientations among all the conformers in the ensemble, it is possible that its orientation might still be incorrect. Thus, it would be helpful if there were a method to assess independently the side-chain orientations determined by NMR. Recently, homology modeling by side-chain packing algorithms has been shown to be successful in predicting the side-chain conformations of the buried residues for a protein when the main-chain coordinates and sequence information are given. Since the main-chain coordinates determined by NMR are consistently more reliable than those of the side-chains, we have applied the side-chain packing algorithms to predict side-chain conformations that are compatible with the NMR-derived backbone. Using four test cases where the NMR solution structures and the X-ray crystal structure of the same protein are available, we demonstrate that the side-chain packing method can provide independent validation for the side-chain conformations of NMR structures. Comparison of the side-chain conformations derived by side-chain packing prediction and by NMR spectroscopy demonstrates that when there is agreement between the NMR model and the predicted model, on average 78% of the time the X-ray structure also concurs. While the side-chain packing method can confirm the reliable residue conformations in NMR models, more importantly, it can also identify the questionable residue conformations with an accuracy of 60%. This validation method can serve to increase the confidence level for potential users of structural models determined by NMR.

Animals↗

PubMed2001

Redox-dependent conformational selection in a Cys4Fe2S2 ferredoxin.

Putidaredoxin (Pdx), a Cys4Fe2S2 ferredoxin from Pseudomonas putida, exhibits redox-dependent binding to its physiological redox partner, cytochrome P450(cam) (CYP101), with the reduced form of Pdx (Pdx(r)) binding with greater affinity to oxidized camphor-bound CYP101 than the oxidized form, Pdx(o). It has been previously shown that Pdx(o) is more dynamic than Pdx(r) on all accessible time scales, and it has been proposed that Pdx(r) samples only a fraction of the conformational substates populated by Pdx(o) on a time average. It is postulated that the ensemble subset populated by Pdx(r) is the same subset that binds CYP101, providing a mechanism for coupling the Pdx oxidation state to binding affinity for CYP101. Evidence from a variety of sources, including redox-dependent shifts of 15N and 13C resonances, indicates that the metal cluster binding loop of Pdx is the primary determinant of redox-dependent conformational selection. Patterns of paramagnetic effects suggest that the metal cluster binding loop contracts around the metal cluster upon reduction, possibly due to the strengthening of hydrogen bonds between the sulfur atoms of the metal cluster and the surrounding polypeptide NH and OH groups. Effects of this perturbation are then transmitted mechanically to other affected regions of the protein. A specific mutation has been introduced into the metal binding loop of Pdx, G40N, that slows conformational exchange sufficiently that the ensemble of conformational substates in Pdx(o) are directly observable as severe broadenings or splittings in affected NMR resonances. Many of the residues most affected by the mutation also show significant exchange contributions to 15N T(2) relaxation in wild-type Pdx(o). As predicted, G40N Pdx(r) shows a collapse of many of these multiplets and broadened lines to form much sharper resonances that are essentially identical to those observed in wild-type Pdx(r), indicating that Pdx(r) occupies fewer conformational substates than does Pdx(o). This is the first direct observation of such redox-dependent ensembles at slow exchange on the chemical shift time scale. These results confirm that conformational selection within the Fe2S2 cluster binding loop is the primary source of redox-dependent changes in protein dynamics in Pdx.

Amino Acid Sequence↗

PubMed1995

The use of CVFF and CFF91 force fields in conformational analysis of carbohydrate molecules. Comparison with AMBER molecular mechanics and dynamics calculations for methyl alpha-lactoside.

The solution conformation of methyl alpha-lactoside has been studied through molecular mechanics calculations using the AMBER/Homans, CVFF and CFF91 force fields, and compared to NMR nuclear Overhauser data. Steady-state and transient nuclear Overhauser effects (NOEs) have been interpreted in terms of the ensemble average distribution of conformers. The NOEs have been analysed using the complete relaxation matrix approach for a rigid and isotropic motion model. The molecular mechanics calculations have been performed at two dielectric constants (i.e. epsilon = 1 and 80 debyes, or epsilon = r and 80 debyes) in an exhaustive way, and, in some cases, have been complemented by specific calculations at intermediate epsilon values. Relaxed energy maps and adiabatic surfaces have been generated for the different dielectric constants. The probability distribution of conformers has been estimated from these steric energy maps. Molecular dynamics simulations in vacuo have also been performed. Our results indicate that the beta-(1-->4) glycosidic linkage shows some fluctuations between three low-energy regions, although it spends about 90% of its time in the region close to the global minimum. The observed conformation of methyl alpha-lactoside seems to be closer to that predicted by CVFF, although the AMBER/Homans results are also in qualitative agreement with the experimental data.

Carbohydrate Conformation↗

PubMed2005

Electrostatic mechanism of nucleosomal array folding revealed by computer simulation.

Although numerous experiments indicate that the chromatin fiber displays salt-dependent conformations, the associated molecular mechanism remains unclear. Here, we apply an irregular Discrete Surface Charge Optimization (DiSCO) model of the nucleosome with all histone tails incorporated to describe by Monte Carlo simulations salt-dependent rearrangements of a nucleosomal array with 12 nucleosomes. The ensemble of nucleosomal array conformations display salt-dependent condensation in good agreement with hydrodynamic measurements and suggest that the array adopts highly irregular 3D zig-zag conformations at high (physiological) salt concentrations and transitions into the extended "beads-on-a-string" conformation at low salt. Energy analyses indicate that the repulsion among linker DNA leads to this extended form, whereas internucleosome attraction drives the folding at high salt. The balance between these two contributions determines the salt-dependent condensation. Importantly, the internucleosome and linker DNA-nucleosome attractions require histone tails; we find that the H3 tails, in particular, are crucial for stabilizing the moderately folded fiber at physiological monovalent salt.

Computer Simulation↗