PubMed Health⌕ Search

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 793 records · Page 44Linked to original sources

Cooperative ordering in homeodomain-DNA recognition: solution structure and dynamics of the MATa1 homeodomain.

The mating type homeodomain proteins, MATa1 and MATalpha2, combine to form a heterodimer to bind DNA in diploid yeast cells. The a1-alpha2 heterodimer tightly and specifically binds haploid-specific gene operators to repress transcription. On its own, however, the a1 homeodomain does not bind DNA in a sequence-specific manner. To help understand this interaction, we describe the solution structure and backbone dynamics of the free a1 homeodomain. Free a1 in solution is an ensemble of structures having flexible hinges at the two turns in the small protein fold. Conformational changes in the a1 homeodomain upon ternary complex formation are located in the loop between helix 1 and helix 2, where the C-terminal tail of alpha2 binds to form the heterodimer, and at the C-terminus of helix 3, the DNA recognition helix. The observed differences, comparing the free and bound a1 structures, suggest a mechanism linking van der Waals stacking changes to the ordering of a final turn in the DNA-binding helix of a1. The tail of alpha2 induces changes in loop 1 of a1 that push it toward a properly folded DNA binding conformation.

Amino Acid Sequence↗

Adsorption of polyelectrolyte molecules to a nanostructured monolayer of amphiphiles.

We investigate the conformation of atactic poly(sodium 4-styrene sulfonate) (PSS) molecules both in water-ethanol solutions and adsorbed to a template provided by a nanostructured monolayer of flat-lying octadecylamines on graphite. The conformations of the PSS molecules in these solutions undergo two conformational transitions, a re-entry transition at 40 vol % and a coil-globule conformational transition at 80 vol % ethanol. On the surface the conformations are characterized by two different length scales: on a short length scale they are defined by the template, whereas the large-scale structures reflect the conformations in solution. We conclude that the conformational transitions are discontinuous on the level of single molecules and appear continuous in the ensemble.

Adsorption↗

Molecular simulation study of the bonded-phase structure in reversed-phase liquid chromatography with neat aqueous solvent.

The dramatic loss of retention in reversed-phase liquid chromatography when switching to 100% aqueous solvent and stopping flow (depressurizing) has long intrigued separation scientists. Recent experimental evidence suggests that the observed loss of retention is due to the loss of pore wetting with subsequent loss of solvent penetration in the porous matrix. One of the prevalent explanations of this phenomenon has been that the bonded phase chains, typically octadecyl silane bound to porous silica, would undergo significant conformational changes, viz. collapse, under pure aqueous conditions. As a definitive means toward elucidating the conformation of bonded-phase chains under pure aqueous conditions, configurational-bias Monte Carlo simulations in the Gibbs ensemble were carried out for a system of dimethyl octadecyl silane of intermediate coverage bound to the (111) face of beta-cristobalite and immersed in pure water. The results of two sets of simulations, which were started from two entirely different starting configurations as a validity check toward reaching the same equilibrium distribution of states, show that chains are neither clustering together nor laying on the surface but rather have a broad distribution of orientations and of conformational states. The interface between the bonded and solvent phases is rough on a molecular level, and clusters of water molecules are sometimes found to adsorb at the silica surface. This computational study lends further evidence that the driving force for the loss of retention when switching to pure aqueous conditions and depressurizing is not the collapse of bonded-phase chains.

Alkanes↗

Trp zipper folding kinetics by molecular dynamics and temperature-jump spectroscopy.

We studied the microsecond folding dynamics of three beta hairpins (Trp zippers 1-3, TZ1-TZ3) by using temperature-jump fluorescence and atomistic molecular dynamics in implicit solvent. In addition, we studied TZ2 by using time-resolved IR spectroscopy. By using distributed computing, we obtained an aggregate simulation time of 22 ms. The simulations included 150, 212, and 48 folding events at room temperature for TZ1, TZ2, and TZ3, respectively. The all-atom optimized potentials for liquid simulations (OPLS(aa)) potential set predicted TZ1 and TZ2 properties well; the estimated folding rates agreed with the experimentally determined folding rates and native conformations were the global potential-energy minimum. The simulations also predicted reasonable unfolding activation enthalpies. This work, directly comparing large simulated folding ensembles with multiple spectroscopic probes, revealed both the surprising predictive ability of current models as well as their shortcomings. Specifically, for TZ1-TZ3, OPLS for united atom models had a nonnative free-energy minimum, and the folding rate for OPLS(aa) TZ3 was sensitive to the initial conformation. Finally, we characterized the transition state; all TZs fold by means of similar, native-like transition-state conformations.

Computer Simulation↗

Single-molecule Michaelis-Menten equations.

This paper summarizes our present theoretical understanding of single-molecule kinetics associated with the Michaelis-Menten mechanism of enzymatic reactions. Single-molecule enzymatic turnover experiments typically measure the probability density f(t) of the stochastic waiting time t for individual turnovers. While f(t) can be reconciled with ensemble kinetics, it contains more information than the ensemble data; in particular, it provides crucial information on dynamic disorder, the apparent fluctuation of the catalytic rates due to the interconversion among the enzyme's conformers with different catalytic rate constants. In the presence of dynamic disorder, f(t) exhibits a highly stretched multiexponential decay at high substrate concentrations and a monoexponential decay at low substrate concentrations. We derive a single-molecule Michaelis-Menten equation for the reciprocal of the first moment of f(t), 1/ , which shows a hyperbolic dependence on the substrate concentration [S], similar to the ensemble enzymatic velocity. We prove that this single-molecule Michaelis-Menten equation holds under many conditions, in particular when the intercoversion rates among different enzyme conformers are slower than the catalytic rate. However, unlike the conventional interpretation, the apparent catalytic rate constant and the apparent Michaelis constant in this single-molecule Michaelis-Menten equation are complicated functions of the catalytic rate constants of individual conformers. We also suggest that the randomness parameter r, defined as <(t - )2> / t2, can serve as an indicator for dynamic disorder in the catalytic step of the enzymatic reaction, as it becomes larger than unity at high substrate concentrations in the presence of dynamic disorder.

Algorithms↗

Foldamer simulations: novel computational methods and applications to poly-phenylacetylene oligomers.

We apply several methods to probe the ensemble kinetic and structural properties of a model system of poly-phenylacetylene (pPA) oligomer folding trajectories. The kinetic methods employed included a brute force accounting of conformations, a Markovian state matrix method, and a nonlinear least squares fit to a minimalist kinetic model used to extract the folding time. Each method gave similar measures for the folding time of the 12-mer chain, calculated to be on the order of 7 ns for the complete folding of the chain from an extended conformation. Utilizing both a linear and a nonlinear scaling relationship between the viscosity and the folding time to correct for a low simulation viscosity, we obtain an upper and a lower bound for the approximate folding time within the range 70 ns<tau<350 ns. This is in agreement with the experimentally measured folding time on the order of 160 ns. The kinetic model used to fit the kinetic behavior of the ensemble of trajectories provides a framework to describe the bulk folding mechanism. We were able to identify two unique clusters of conformations that provide a structural basis to account for the appearance of a kinetic intermediate in the mechanism. We discuss the implications of these findings in the context of helix-coil theory.

Journal Article↗

Ca2+ binding and conformational changes in a calmodulin domain.

Calcium activation of the C-terminal domain of calmodulin was studied using 1H and 15N NMR spectroscopy. The important role played by the conserved bidentate glutamate Ca2+ ligand in the binding loops is emphasized by the striking effects resulting from a mutation of this glutamic acid to a glutamine, i.e. E104Q in loop III and E140Q in loop IV. The study involves determination of Ca2+ binding constants, assignments, and structural characterizations of the apo, (Ca2+)1, and (Ca2+)2 states of the E104Q mutant and comparisons to the wild-type protein and the E140Q mutant [Evenäs et al. (1997) Biochemistry 36, 3448-3457]. NMR titration data show sequential Ca2+ binding in the E104Q mutant. The first Ca2+ binds to loop IV and the second to loop III, which is the order reverse to that observed for the E140Q mutant. In both mutants, the major structural changes occur upon Ca2+ binding to loop IV, which implies a different response to Ca2+ binding in the N- and C-terminal EF-hands. Spectral characteristics show that the (Ca2+)1 and (Ca2+)2 states of the E104Q mutant undergo global exchange on a 10-100 micros time scale between conformations seemingly similar to the closed and open structures of this domain in wild-type calmodulin, paralleling earlier observations for the (Ca2+)2 state of the E140Q mutant, indicating that both glutamic acid residues, E104 and E140, are required for stabilization of the open conformation in the (Ca2+)2 state. To verify that the NOE constraints cannot be fulfilled in a single structure, solution structures of the (Ca2+)2 state of the E104Q mutant are calculated. Within the ensemble of structures the precision is good. However, the clearly dynamic nature of the state, a large number of violated distance restraints, ill-defined secondary structural elements, and comparisons to the structures of calmodulin indicate that the ensemble does not provide a good picture of the (Ca2+)2 state of the E104Q mutant but rather represents the distance-averaged structure of at least two distinct different conformations.

Amino Acid Sequence↗

Force field impact and spin-probe modeling in molecular dynamics simulations of spin-labeled T4 lysozyme.

Several attempts have been made to compute electron paramagnetic resonance (EPR) spectra of biomolecules, using motional models or simulated trajectories to describe dynamics. Ideally, the simulated trajectories should capture "fast" (picosecond) snapshots of spin-probe rotations accurately, while being lengthy enough to ensure a proper Fourier integration of the time-domain signal. It is the interplay of the two criteria that poses computational challenges to the method. In this context, an analysis of the spin-probe and protein conformational sampling and equilibration, with different force fields and with explicit solvent, may be a useful attempt. The present work reports a comparative study of the effect of the molecular dynamics (MD) force field on conformational sampling and equilibration in two spin-labeled T4 lysozyme (T4L) variants, N40C and K48C. Ensembles of 10x 3 ns-trajectories per variant and per force field (OPLS/AMBER and AMBER99) are analyzed for a reliable assessment of convergence and sampling. It is found that subtle site-dependent differences in spin-probe rotations and torsions are more readily captured in the AMBER99 trajectories than in the OPLS/AMBER simulations. On the other hand, sampling and equilibration are found to be better with the OPLS/AMBER force field at equal trajectory lengths.

Bacteriophage T4↗

Assembly of the full-length recombinant mouse prion protein I. Formation of soluble oligomers.

The conversion of a monomeric alpha-helix-rich isoform to multimeric beta-sheet-rich isoforms is a prominent feature of the conversion between PrP(C) and PrP(SC). We mimicked this process in vitro by exposing an unglycosylated recombinant form of the full-length mouse prion protein ((Mo)PrP(23-231)) to an acidic pH, at 37 degrees C, and we monitored the kinetics of conformational change and assembly. In these conditions, monomeric (Mo)PrP(23-231) converts slowly to two ensembles of soluble oligomers that are separated by size exclusion chromatography. The larger oligomers (I) are unstable, and their formation involves almost no change in secondary structure content. The smaller oligomers (II) form stable spherical or annular particles containing between 8 and 15 monomers as determined by multi-angle laser light scattering (MALLS). Their formation is concomitant with the main, thought limited, change in the secondary structure content (10%) seen by Fourier Transform Infrared (FTIR) spectroscopy. Even if these oligomers conserve a large part of the secondary structure of monomeric PrP, they exhibit amyloid features with the appearance of intermolecular beta-structure as revealed by the appearance of an IR band below 1620 cm(-1).

Animals↗

Fine-grained statistical torsion angle potentials are effective in discriminating native protein structures.

Modelling of drug targets requires the reliable selection of an accurate and representative structure from large ensembles of alternate models. Statistical potentials developed to discriminate native protein structures generally represent pairwise interactions between atoms, which are less sensitive to local conformational details. The discrimination of local distortions is therefore particularly difficult. Local interaction preferences, expressed through torsion angles, are rarely used, as some controversy exists in the literature regarding their discrimination power. The present study aims to benchmark the efficiency of different implementations of torsion angle propensities for selecting the native structure from ensembles of well-constructed decoys. Several statistical potentials derived from fine-grained discretisations of torsion angle space are constructed and evaluated. Results from a comparison with nine widely used statistical scoring functions show the torsion angle potentials to be more effective in recognising native structures and to improve with the number of torsion angles considered. These data suggest local structural propensities to be important for estimating the overall quality of native-like models.

Models, Chemical↗

Defining protein ensembles with native-state NH exchange: kinetics of interconversion and cooperative units from combined NMR and MS analysis.

Previous studies of native-state peptide hydrogen atom (NH) exchange in turkey ovomucoid third domain (OMTKY3) yielded the thermodynamics and kinetics of unfolding and folding for the 14 slowest-exchanging peptide hydrogen atoms (NHs). Unfolding rate constants and free energies for nine of the NHs are very similar, suggesting that these NHs exchange during a single cooperative unfolding event. Electrospray ionization mass spectrometry (ESI-MS) has been used to test this hypothesis. ESI-MS data and MS peak simulations suggest that this hypothesis is incorrect: in spite of the similarity in their unfolding rate constants, only three to five of the nine residues exchange in a cooperative manner. Thus, residues with similar thermodynamics and kinetics of exchange are probably involved in multiple conformational equilibria. Overall, combined NMR and MS analysis of NH exchange provides a rich and complex picture of the ensemble properties of native proteins.

Animals↗

Stochastic algorithm for kinase homology model construction.

A stochastic algorithm for constructing multiple loops in homology modeling of proteins is presented. The algorithm discards variable values in iterations based on a cost function and on statistical analysis of results. Values that remain are used for constructing an ensemble of best solutions. In test cases, the stochastic algorithm retains all the best solutions, compared to an exhaustive scan of the full set of conformations. Individual loops are constructed by adding dipeptide units. Dipeptide conformations are extracted from a database of proteins and their conformations include bond lengths and all angles. Single loops are constructed from both N- and C- terminals to the center, and loop closure is evaluated by a combination of penalties for the peptide closure and Miyazawa-Jernigan (MJ) [1]. residue-residue interactions with the rest of the protein. Large ensembles of each loop are clustered and re-evaluated with a refined [2]. energy term. The reduced, clustered set of each loop is then employed to construct simultaneously all the loops. The algorithm was applied to construct simultaneously six loops in c-Src kinase family proteins, incorporating a total of 37-40 residues. The best RMSD for reconstructing the loops is 1.45 A for Lck (structure 1QPE in the Protein Data Bank) and 2.54 A for human c-Src (structure 1FMK). The multiple loop conformations with lowest energy have higher RMSD values, of 2.06 A and 3.09 A, respectively. The average RMSD values for the first 1000 conformations are 3.00 A and 3.46 A, respectively. Models for the "open" structures of c-Src and of Jak-2 were constructed on the basis of 1QPE. The Jak-2 model is found to be more flexible in the loops region than its c-Src counterpart.

Algorithms↗

Heterogeneity and inaccuracy in protein structures solved by X-ray crystallography.

Proteins are dynamic molecules, exhibiting structural heterogeneity in the form of anisotropic motion and discrete conformational substates, often of functional importance. In protein structure determination by X-ray crystallography, the observed diffraction pattern results from the scattering of X-rays by an ensemble of heterogeneous molecules, ordered and oriented by packing in a crystal lattice. The majority of proteins diffract to resolutions where heterogeneity is difficult to identify and model, and are therefore approximated by a single, average conformation with isotropic variance. Here we show that disregarding structural heterogeneity introduces degeneracy into the structure determination process, as many single, isotropic models exist that explain the diffraction data equally well. The large differences among these models imply that the accuracy of crystallographic structures has been widely overestimated. Further, it suggests that analyses that depend on small differences in the relative positions of atoms may be flawed.

Crystallography, X-Ray↗

Conformational analysis of oligosaccharides corresponding to the cell-wall polysaccharide of the Streptococcus group A by Metropolis Monte Carlo simulations.

Metropolis Monte Carlo simulations have been performed on four substructures from the cell-wall polysaccharide antigen of Streptococcus group A to explore the conformational behaviour of these compounds. The compounds examined are the trisaccharide, propyl 3-O-(2-acetamido-2-deoxy-beta-D-glucopyranosyl)-2-O-(alpha-L-rhamnopyranosyl)- alpha-L-rhamnopyranoside, 1, the tetrasaccharide, propyl 3-O-(3-O-(2-acetamido-2-deoxy-beta-D- glucopyranosyl)-2-O-(alpha-L-rhamnopyranosyl)-alpha-L-rhamnopyranosyl)-alpha-L- rhamnopyranoside, 2, the hexasaccharide, propyl 3-O-(2-O-(3-O-(3-O-(2-acetamido-2-deoxy-beta-D-glucopyranosyl)-alpha-L- rhamnopyranosyl)-alpha-L-rhamnopyranosyl)-3-O-(2-acetamido-2-deoxy-beta-D- glucopyranosyl)-alpha-L-rhamnopyranosyl)-alpha-L-rhamnopyranoside, 3, and the hexasaccharide, propyl 3-O-(2-acetamido-2-deoxy-beta-D-glucopyranosyl)-2-O-(3-O-(3-O-(2-acetamido-2- deoxy-beta-D-glucopyranosyl)-2-O-(alpha-L-rhamnopyranosyl)-alpha-L- rhamnopyranosyl)-alpha-L-rhamnopyranosyl)-alpha-L-rhamnopyranoside, 4. In general, the conformational flexibility of similar glycosidic linkages in different compounds is comparable. However, in a few cases, small differences in the conformations available to these linkages in different structural environments could be detected. Interestingly, a second conformation found for the beta-D-GlcNAc-(1-->3)-alpha-L-Rha linkage in three of the compounds was not populated in the hexasaccharide 4. Furthermore, a conformational locale of the alpha-L-Rha-(1-->3)-alpha-L-Rha linkage found to be populated in the trisaccharide 1, tetrasaccharide 2, and hexasaccharide 4 is negligibly populated in the hexasaccharide 3. Ensemble averaged proton-proton distances compare favourably with experimental average distances obtained from NMR spectroscopy. The trisaccharide branch point in the hexasaccharides is shown to be a highly defined conformational feature. The same unit has been found to be one of the crucial elements recognized by anti-Group A Streptococcus antibodies, a result that has implications for the design of improved immunodiagnostics and vaccines.

Carbohydrate Conformation↗

Crystallization of a designed peptide from a molten globule ensemble.

Backgound. The design of amino acid sequences that adopt a desired three-dimensional fold has been of keen interest over the past decade. However, the design of proteins that adopt unique conformations is still a considerable problem. Until very recently, all of the designed proteins that have been extensively characterized possess the hallmarks of the molten globular state. Molten globular intermediates have been observed in both equilibrium and kinetic protein folding/stability studies, and understanding the forces that determine compact non-native states is critical for a comprehensive understanding of proteins. This paper describes the solution and early solid state characterization of peptides that form molten globular ensembles. Results. Crystals diffracting to 3.5&angst; resolution have been grown of a 16-residue peptide (alpha1A) designed to form a tetramer of alpha-helices. In addition, a closely related peptide, alpha1, has previously been shown to yield crystals that diffract to 1.2&angst; resolution. The solution properties of these two peptides were examined to determine whether their well defined crystalline conformations were retained in solution. On the basis of an examination of their NMR spectra, sedimentation equilibria, thermal unfolding, and ANS binding, it is concluded that the peptides form alpha-helical aggregates with properties similar to those of the molten globule state. Thus, for these peptides, the process of crystallization bears many similarities to models of protein folding. Upon dissolution, the peptides rapidly assume compact molten globular states similar to the molten globule like intermediates that are formed at short times after refolding is initiated. Following a rate-determining nucleation step, the peptides crystallize into a single or a small number of conformations in a process that mimics the formation of native structure in proteins.

Journal Article↗

Crystallization of a designed peptide from a molten globule ensemble.

BACKGROUND: The design of amino acid sequences that adopt a desired three-dimensional fold has been of keen interest over the past decade. However, the design of proteins that adopt unique conformations is still a considerable problem. Until very recently, all of the designed proteins that have been extensively characterized possess the hallmarks of the molten globular state. Molten globular intermediates have been observed in both equilibrium and kinetic protein folding/stability studies, and understanding the forces that determine compact non-native states is critical for a comprehensive understanding of proteins. This paper describes the solution and early solid state characterization of peptides that form molten globular ensembles. RESULTS & CONCLUSIONS: Crystals diffracting to 3.5 A resolution have been grown of a 16-residue peptide (alpha 1A) designed to form a tetramer of alpha-helices. In addition, a closely related peptide, alpha 1, has previously been shown to yield crystals that diffract to 1.2 A resolution. The solution properties of these two peptides were examined to determine whether their well defined crystalline conformations were retained in solution. On the basis of an examination of their NMR spectra, sedimentation equilibria, thermal unfolding, and ANS binding, it is concluded that the peptides form alpha-helical aggregates with properties similar to those of the molten globule state. Thus, for these peptides, the process of crystallization bears many similarities to models of protein folding. Upon dissolution, the peptides rapidly assume compact molten globular states similar to the molten globule like intermediates that are formed at short times after refolding is initiated. Following a rate-determining nucleation step, the peptides crystallize into a single or a small number of conformations in a process that mimics the formation of native structure in proteins.

Amino Acid Sequence↗

A newly synthesized, ribosome-bound polypeptide chain adopts conformations dissimilar from early in vitro refolding intermediates.

Little is known about the conformations of newly synthesized polypeptide chains as they emerge from the large ribosomal subunit, or how these conformations compare with those populated immediately after dilution of polypeptide chains out of denaturant in vitro. Both in vivo and in vitro, partially folded intermediates of the tailspike protein from Salmonella typhimurium phage P22 can be trapped in the cold. A subset of monoclonal antibodies raised against tailspike recognize partially folded intermediates, whereas other antibodies recognize only later intermediates and/or the native state. We have used a pair of monoclonal antibodies to probe the conformational features of full-length, newly synthesized tailspike chains recovered on ribosomes from phage-infected cells. The antibody that recognizes early intermediates in vitro also recognizes the ribosome-bound intermediates. Surprisingly, the antibody that did not recognize early in vitro intermediates did recognize ribosome-bound tailspike chains translated in vivo. Thus, the newly synthesized, ribosome-bound tailspike chains display structured epitopes not detected upon dilution of tailspike chains from denaturant. As opposed to the random ensemble first populated when polypeptide chains are diluted out of denaturant, folding in vivo from the ribosome may begin with polypeptide conformations already directed toward the productive folding and assembly pathway.

Antibodies, Monoclonal↗

Single-molecule RNA folding.

Single-molecule experiments significantly expand our capability to characterize complex dynamics of biological processes. This relatively new approach has contributed significantly to our understanding of the RNA folding problem. Recent single-molecule experiments, together with structural and biochemical characterizations of RNA at the ensemble level, show that RNA molecules typically fold across a highly rugged energy landscape. As a result, long-lived folding intermediates, multiple folding pathways, and heterogeneous conformational dynamics are commonly found for RNA enzymes. While initial results have suggested that stable secondary structures are partly responsible for the rugged energy landscape of RNA, a complete mechanistic understanding of the complex folding behavior has not yet been obtained. A combination of single-molecule experiments, which are well suited to analyze transient and heterogeneous dynamic behaviors, with ensemble characterizations that can provide structural information at a superior resolution will likely provide more answers.

Fluorescence Resonance Energy Transfer↗