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 685 records · Page 38Linked to original sources

13C NMR relaxation studies of RNA base and ribose nuclei reveal a complex pattern of motions in the RNA binding site for human U1A protein.

The widespread importance of induced fit and order-disorder transition in RNA recognition by proteins and small molecules makes it imperative that RNA motional properties are characterized quantitatively. Until now, however, very few studies have been dedicated to the systematic characterization of RNA motion and to their changes upon protein or small-molecule binding. The U1A protein-RNA complexes provide some of the best-studied examples of the role of RNA motional changes upon protein binding. Here, we report (13)C NMR relaxation studies of base and ribose dynamics for the RNA internal loop target of human U1A protein located within the 3'-untranslated region (3'-UTR) of the mRNA coding for U1A itself. We also report the semi-quantitative analysis of both fast (nano- to picosecond) and intermediate (micro- to millisecond) motions for this paradigmatic RNA system. We measure (13)C T(1), T(1rho) and heteronuclear nuclear Overhauser effects (NOEs) for sugar and base nuclei, as well as the power dependence of T(1rho) at 500 MHz and 750 MHz, and analyze these results using the model-free formalism. The results provide a much clearer picture of the type of motions experienced by this RNA in the absence of the protein than was provided by the analysis of the structure based solely on NOEs and scalar couplings. They define a model where the RNA internal loop region "breathes" on a micro- to millisecond timescale with respect to the double-helical regions. Superimposed on this slower motion, the residues at the very tip of the loop undergo faster (nano- to picosecond) motions. We hypothesize that these motions allow the RNA to sample multiple conformations so that the protein can select a structure within the ensemble that optimizes intermolecular contacts.

Base Sequence↗

Solution structure of a peptide fragment of human alpha-lactalbumin in trifluoroethanol: a model for local structure in the molten globule.

BACKGROUND: At low pH, human alpha-lactalbumin forms a partly folded molten globule state that contains a non-native clustering of the side chains of Tyr103, Trp104 and His107. In order to understand the conformation of this region of the protein in the molten globule state, we investigated the structure of a peptide corresponding to residues 101-110 of human alpha-lactalbumin in trifluoroethanol. RESULTS: We determined the structure of the 101-110 peptide from an NMR data set of 145 nuclear Overhauser effects and nine 3JHN alpha coupling constants, using an ensemble calculation approach to take into account the possibilities of conformational averaging of the data. The backbone of residues 3-10 in the peptide adopts a series of turns, that involving residues 5-8 being the best defined, while the side chains of residues 1, 3, 4, 5, 6 and 7 form a hydrophobic cluster. CONCLUSIONS: The peptide conformation differs from that previously determined for residues 101-110 in crystal structures of native alpha-lactalbumin determined at both high and low pH, particularly in the relative orientations of the side chains. The series of turns seen in the peptide could, however, be related to the alpha-helical structure seen for residues 104-111 in crystals at high pH, and may be important in the molten globule state for bringing the peptide chain into a compact conformation where favourable interactions between the side chains can occur.

Amino Acid Sequence↗

Conservation and specialization in PAS domain dynamics.

The PAS (Per-ARNT-Sim) superfamily is presented as a well-suited study case to demonstrate how comparison of functional motions among distant homologous proteins with conserved fold characteristics may give insight into their functional specialization. Based on the importance of structural flexibility of the receptive structures in anticipating the signal-induced conformational changes of these sensory systems, the dynamics of these structures were analysed. Molecular dynamics was proved to be an effective method to obtain a reliable picture of the dynamics of the crystal structures of HERG, phy3, PYP and FixL, provided that an extensive conformational space sampling is performed. Other reliable sources of dynamic information were the ensembles of NMR structures of hPASK, HIF-2alpha and PYP. Essential dynamics analysis was successfully employed to extract the relevant information from the sampled conformational spaces. Comparison of motion patterns in the essential subspaces, based on the structural alignment, allowed identification of the specialized region in each domain. This appears to be evolved in the superfamily by following a specific trend, that also suggests the presence of a limited number of general solutions adopted by the PAS domains to sense external signals. These findings may give insight into unknown mechanisms of PAS domains and guide further experimental studies.

Adiantum↗

Statistical coil model of the unfolded state: resolving the reconciliation problem.

An unfolded state ensemble is generated by using a self-avoiding statistical coil model that is based on backbone conformational frequencies in a coil library, a subset of the Protein Data Bank. The model reproduces two apparently contradicting behaviors observed in the chemically denatured state for a variety of proteins, random coil scaling of the radius of gyration and the presence of significant amounts of local backbone structure (NMR residual dipolar couplings). The most stretched members of our unfolded ensemble dominate the residual dipolar coupling signal, whereas the uniformity of the sign of the couplings follows from the preponderance of polyproline II and beta conformers in the coil library. Agreement with the NMR data substantially improves when the backbone conformational preferences include correlations arising from the chemical and conformational identity of neighboring residues. Although the unfolded ensembles match the experimental observables, they do not display evidence of native-like topology. By providing an accurate representation of the unfolded state, our statistical coil model can be used to improve thermodynamic and kinetic modeling of protein folding.

Computational Biology↗

Recognition of native structure from complete enumeration of low-resolution models with constraints.

Complete sets of low-resolution conformations are generated for eight small proteins by rotating the Calpha-Calpha virtual bonds at selected flexible regions, while the remaining structural elements are assumed to move in rigid blocks. Several filtering criteria are used to reduce the ensemble size and to ensure the sampling of well-constructed conformations. These filters, based on structure and energy constraints deduced from knowledge-based studies, include the excluded volume requirement, the radius of gyration constraint, and the occurrence of sufficiently strong attractive inter-residue potentials to stabilize compact forms. About 8,000 well-constructed decoys or "probable folds" (PFs) are constructed for each protein. A correlation between root-mean-square (rms) deviations from X-ray structure and total energies is observed, revealing a decrease in energy as the rms deviation decreases. The conformation with the lowest energy exhibits an rms deviation smaller than 3.0 A, in most of the proteins considered. The results are highly sensitive to the choice of flexible regions. A strong tendency to assume native state rotational angles is revealed for some flexible bonds from the analysis of the distributions of dihedral angles in the PFs, suggesting the formation of foldons near these locally stable regions at early folding pathway.

Computer Simulation↗

Small worlds in RNA structures.

I consider conformational spaces of tRNA(phe) defined by sets of suboptimal structures from the perspective of small-world networks. Herein, the influence of modifications on typical small-world network properties and the shape of energy landscapes is discussed. Results indicate that natural modifications influence the degree of local clustering and mean path lengths far more than random or no modifications. High frequencies in the thermodynamic ensemble coincide with high numbers of neighboring structures that one conformation can adopt by one elementary move. Conformation spaces indicate the existence of modular substructures. It can be shown that modifications leave the nature of small-world topology untouched albeit natural modifications have a reasonable enhancing and streamlining effect on the degree of clustering and therefore on the substructures of the conformational space.

Algorithms↗

Simple and versatile restraints for the accurate modeling of alpha-helical coiled-coil structures of multiple strandedness, orientation and composition.

We present a minimalist approach for the modeling of the three-dimensional structure of multistranded alpha-helical coiled coils. The approach is based on empirical principles introduced by F. H. C. Crick (F. H. C. Crick, Acta Crystallogr, 1953, Vol. 6, pp. 689-697). Crick hypothesized that keeping the distance between the residues at the interacting interface of alpha-helices constant would lead to supercoiling or the formation of a coiled coil through the knobs-into-holes mode of packing. We have implemented the latter hypothesis in a simulating annealing protocol in the simple form of interhelical distance restraints (two per heptad) between Calpha at the interfacial positions and. To demonstrate the authenticity of Crick's hypothesis and the precision and accuracy of our approach, we have modeled the crystal structures of six synthetic coiled coils in dimeric, trimeric, and tetrameric states. The mean root mean square deviations (RMSDs) between the backbone atoms of the ensemble of structures calculated and those of the corresponding geometric averages is always below 0.76 A, indicating that our protocol has an excellent degree of convergence and precision. The RMSDs between the backbone atoms of each of the six geometric average structures and the backbone of the corresponding crystal structures all range between 0.43 and 0.95 A, indicative of excellent accuracy and proving the authenticity of Crick's hypothesis. Moreover, without specifying any dihedral angles, we found that in 81% of the occurrences, the most populated conformer of the side chains at positions and in the ensembles calculated were identical to those observed in the crystal structures. This shows that our simple approach, which is the simplest reported so far, can generate accurate results for the backbone and side chains. Finally, as a test case for a wider application of our approach in the field of structural proteomics, we describe the successful modeling of the overall structure of SNARE and the organization of its interfacial ionic layer known to play an important functional role.

Amino Acid Sequence↗

Role of protein conformation changes and transphosphorylations in the function of Na+/K(+)-transporting adenosine triphosphatase: an attempt at an integration into the Na+/K+ pump mechanism.

The particular aim of the review on some basic facets of the mechanism of Na+/K(+)-transporting ATPase (Na/K-ATPase) has been to integrate the experimental findings concerning the Na(+)- and K(+)-elicited protein conformation changes and transphosphorylations into the perspective of an allosterically regulated, phosphoryl energy transferring enzyme. This has led the authors to the following summarizing evaluations. 1. The currently dominating hypothesis on a link between protein conformation changes ('E1 in equilibrium with E2') and Na+/K+ transport (the 'Albers-Post scheme') has been constructed from a variety of partial reactions and elementary steps, which, however, do not all unequivocally support the hypothesis. 2. The Na(+)- and K(+)-elicited protein conformation changes are inducible by a variety of other ligands and modulatory factors and therefore cannot be accepted as evidence for their direct participation in effecting cation translocation. 3. There is no evidence that the 'E1 in equilibrium with E2' protein conformation changes are moving Na+ and K+ across the plasma membrane. 4. The allosterically caused ER in equilibrium with ET ('E1 in equilibrium with E2') conformer transitions and the associated cation 'occlusion' in equilibrium with 'de-occlusion' processes regulate the actual catalytic power of an enzyme ensemble. 5. A host of experimental variables determines the proportion of functionally competent ER enzyme conformers and incompetent ET conformers so that any enzyme population, even at the start of a reaction, consists of an unknown mixture of these conformers. These circumstances account for the occurrence of contradictory observations and apparent failures in their comparability. 6. The modelling of the mechanism of the Na/K-ATPase and Na+/K+ pump from the results of reductionistically designed experiments requires the careful consideration of the physiological boundary conditions. 7. Na+ and K+ ligandation of Na/K-ATPase controls the geometry and chemical reactivity of the catalytic centre in the cycle of E1 in equilibrium with E2 state conversions. This is possibly effected by hinge-bending, concerted motions of three adjacent, intracellularly exposed peptide sequences, which shape open and closed forms of the catalytic centre in lock-and-key responses. 8. The Na(+)-dependent enzyme phosphorylation with ATP and the K(+)-dependent hydrolysis of the phosphoenzyme formed are integral steps in the transport mechanism of Na/K-ATPase, but the translocations of Na+ and K+ do not occur via a phosphate-cation symport mechanism.(ABSTRACT TRUNCATED AT 400 WORDS)

Animals↗

Identification of protein-protein interaction sites from docking energy landscapes.

Protein recognition is one of the most challenging and intriguing problems in structural biology. Despite all the available structural, sequence and biophysical information about protein-protein complexes, the physico-chemical patterns, if any, that make a protein surface likely to be involved in protein-protein interactions, remain elusive. Here, we apply protein docking simulations and analysis of the interaction energy landscapes to identify protein-protein interaction sites. The new protocol for global docking based on multi-start global energy optimization of an all-atom model of the ligand, with detailed receptor potentials and atomic solvation parameters optimized in a training set of 24 complexes, explores the conformational space around the whole receptor without restrictions. The ensembles of the rigid-body docking solutions generated by the simulations were subsequently used to project the docking energy landscapes onto the protein surfaces. We found that highly populated low-energy regions consistently corresponded to actual binding sites. The procedure was validated on a test set of 21 known protein-protein complexes not used in the training set. As much as 81% of the predicted high-propensity patch residues were located correctly in the native interfaces. This approach can guide the design of mutations on the surfaces of proteins, provide geometrical details of a possible interaction, and help to annotate protein surfaces in structural proteomics.

Binding Sites↗

Slow dynamics in folded and unfolded states of an SH3 domain.

(15)N relaxation dispersion experiments were applied to the isolated N-terminal SH3 domain of the Drosophila protein drk (drkN SH3) to study microsecond to second time scale exchange processes. The drkN SH3 domain exists in equilibrium between folded (F(exch)) and unfolded (U(exch)) states under nondenaturing conditions in a ratio of 2:1 at 20 degrees C, with an average exchange rate constant, k(ex), of 2.2 s(-1) (slow exchange on the NMR chemical shift time scale). Consequently a discrete set of resonances is observed for each state in NMR spectra. Within the U(exch) ensemble there is a contiguous stretch of residues undergoing conformational exchange on a micros/ms time scale, likely due to local, non-native hydrophobic collapse. For these residues both the F(exch) <--> U(exch) conformational exchange process and the micros/ms exchange event within the U(exch) state contribute to the (15)N line width and can be analyzed using CPMG-based (15)N relaxation dispersion measurements. The contribution of both processes to the apparent relaxation rate can be deconvoluted numerically by combining the experimental (15)N relaxation dispersion data with results from an (15)N longitudinal relaxation experiment that accurately quantifies exchange rates in slow exchanging systems (Farrow, N. A.; Zhang, O.; Forman-Kay, J. D.; Kay, L. E. J. Biomol. NMR 1994, 4, 727-734). A simple, generally applicable analytical expression for the dependence of the effective transverse relaxation rate constant on the pulse spacing in CPMG experiments has been derived for a two-state exchange process in the slow exchange limit, which can be used to fit the experimental data on the global folding/unfolding transition. The results illustrate that relaxation dispersion experiments provide an extremely sensitive tool to probe conformational exchange processes in unfolded states and to obtain information on the free energy landscape of such systems.

Animals↗

Distinct early folding and aggregation properties of Alzheimer amyloid-beta peptides Abeta40 and Abeta42: stable trimer or tetramer formation by Abeta42.

The amyloid beta peptide (Abeta), composed of 40 or 42 amino acids, is a critical component in the etiology of the neurodegenerative Alzheimer disease. Abeta is prone to aggregate and forms amyloid fibrils progressively both in vitro and in vivo. To understand the process of amyloidogenesis, it is pivotal to examine the initial stages of the folding process. We examined the equilibrium folding properties, assembly states, and stabilities of the early folding stages of Abeta40 and Abeta42 prior to fibril formation. We found that Abeta40 and Abeta42 have different conformations and assembly states upon refolding from their unfolded ensembles. Abeta40 is predominantly an unstable and collapsed monomeric species, whereas Abeta42 populates a stable structured trimeric or tetrameric species at concentrations above approximately 12.5 microm. Thermodynamic analysis showed that the free energies of Abeta40 monomer and Abeta42 trimer/tetramer are approximately 1.1 and approximately 15/ approximately 22 kcal/mol, respectively. The early aggregation stages of Abeta40 and Abeta42 contain different solvent-exposed hydrophobic surfaces that are located at the sequences flanking its protease-resistant segment. The amyloidogenic folded structure of Abeta is important for the formation of spherical beta oligomeric species. However, beta oligomers are not an obligatory intermediate in the process of fibril formation because oligomerization is inhibited at concentrations of urea that have no effect on fibril formation. The distinct initial folding properties of Abeta40 and Abeta42 may play an important role in the higher aggregation potential and pathological significance of Abeta42.

Alzheimer Disease↗

Free-energy landscape of a chameleon sequence in explicit water and its inherent alpha/beta bifacial property.

A sequence in yeast MATalpha2/MCM1/DNA complex that folds into alpha-helix or beta-hairpin depending on the surroundings has been known as "chameleon" sequence. We obtained the free-energy landscape of this sequence by using a generalized-ensemble method, multicanonical molecular dynamics simulation, to sample the conformational space. The system was expressed with an all-atom model in explicit water, and the initial conformation for the simulation was a random one. The free-energy landscape demonstrated that this sequence inherently has an ability to form either alpha or beta structure: The conformational distribution in the landscape consisted of two alpha-helical clusters with different packing patterns of hydrophobic residues, and four beta-hairpin clusters with different strand-strand interaction patterns. Narrow pathways connecting the clusters were found, and analysis on the pathways showed that a compact structure formed at the N-terminal root of the chameleon sequence controls the cluster-cluster transitions. The free-energy landscape indicates that a small conditional change induces alpha-beta transitions. Additional unfolding simulations done with replacing amino acids showed that the chameleon sequence has an advantage to form an alpha-helix. Current study may be useful to understand the mechanism of diseases resulting from abnormal chain folding, such as amyloid disease.

Amino Acid Sequence↗

Unusual compactness of a polyproline type II structure.

Polyproline type II (PPII) helix has emerged recently as the dominant paradigm for describing the conformation of unfolded polypeptides. However, most experimental observables used to characterize unfolded proteins typically provide only short-range, sequence-local structural information that is both time- and ensemble-averaged, giving limited detail about the long-range structure of the chain. Here, we report a study of a long-range property: the radius of gyration of an alanine-based peptide, Ace-(diaminobutyric acid)2-(Ala)7-(ornithine)2-NH2. This molecule has previously been studied as a model for the unfolded state of proteins under folding conditions and is believed to adopt a PPII fold based on short-range techniques such as NMR and CD. By using synchrotron radiation and small-angle x-ray scattering, we have determined the radius of gyration of this peptide to be 7.4 +/- 0.5 angstroms, which is significantly less than the value expected from an ideal PPII helix in solution (13.1 angstroms). To further study this contradiction, we have used molecular dynamics simulations using six variants of the AMBER force field and the GROMOS 53A6 force field. However, in all cases, the simulated ensembles underestimate the PPII content while overestimating the experimental radius of gyration. The conformational model that we propose, based on our small angle x-ray scattering results and what is known about this molecule from before, is that of a very flexible, fluctuating structure that on the level of individual residues explores a wide basin around the ideal PPII geometry but is never, or only rarely, in the ideal extended PPII helical conformation.

Circular Dichroism↗

Automated prediction of CASP-5 structures using the Robetta server.

Robetta is a fully automated protein structure prediction server that uses the Rosetta fragment-insertion method. It combines template-based and de novo structure prediction methods in an attempt to produce high quality models that cover every residue of a submitted sequence. The first step in the procedure is the automatic detection of the locations of domains and selection of the appropriate modeling protocol for each domain. For domains matched to a homolog with an experimentally characterized structure by PSI-BLAST or Pcons2, Robetta uses a new alignment method, called K*Sync, to align the query sequence onto the parent structure. It then models the variable regions by allowing them to explore conformational space with fragments in fashion similar to the de novo protocol, but in the context of the template. When no structural homolog is available, domains are modeled with the Rosetta de novo protocol, which allows the full length of the domain to explore conformational space via fragment-insertion, producing a large decoy ensemble from which the final models are selected. The Robetta server produced quite reasonable predictions for targets in the recent CASP-5 and CAFASP-3 experiments, some of which were at the level of the best human predictions.

Algorithms↗

A coarse graining method for the identification of transition rates between molecular conformations.

The coarse graining method to be advocated in this paper consists of two main steps. First, the propagation of an ensemble of molecular states is described as a Markov chain by a transition probability matrix in a finite state space. Second, we obtain metastable conformations by an aggregation of variables via Robust Perron Cluster Analysis (PCCA+). Up to now, it has been an open question as to how this coarse graining in space can be transformed to a coarse graining of the Markov chain while preserving the essential dynamic information. In this article, we construct a coarse matrix that is the correct propagator in the space of conformations. This coarse graining procedure carries over to rate matrices and allows to extract transition rates between molecular conformations. This approach is based on the fact that PCCA+ computes molecular conformations as linear combinations of the dominant eigenvectors of the transition matrix.

Journal Article↗

The native energy landscape for interleukin-1beta. Modulation of the population ensemble through native-state topology.

A minimalist Go-model, with no energetic frustration in the native conformation, has been shown to describe accurately the folding pathway of the beta-trefoil protein, interleukin-1beta (IL-1beta). While it appears that these models successfully model transition states and intermediates between the unfolded and native ensembles, it is unclear how accurately they capture smaller, yet biologically relevant, structural changes within the native ensemble after energetic perturbation. Here, we address the following questions. Can a simple Go-model of interleukin-1beta, based on native topology, describe changes in structural properties of the native ensemble as the protein stability is changed? Or is it necessary to include a more explicit representation of atoms, electrostatic, hydrogen bonding, and van der Waals forces to describe these changes? The native ensemble of IL-1beta was characterized using a variety of experimental probes under native (0 M NaCl, guanidine hydrochloride (Gdn-HCl)), moderately destabilized (0 M NaCl, 0.8 M Gdn-HCl), and in moderate salt concentration (0.8 M NaCl, 0 M Gdn-HCl). Heteronuclear (1)H-(15)N nuclear Overhauser effect spectroscopy (NOESY) and heteronuclear single quantum correlation (HSQC) NMR spectra confirmed that the beta-trefoil global fold was largely intact under these three conditions. However, 25 of the 153 residues throughout the chain did demonstrate (13)C and (1)H-(15)N chemical shifts when perturbed with 0.8 M NaCl or Gdn-HCl. Despite large differences in protection factors from solvent hydrogen-deuterium exchange for all residues between stable (0 M Gdn-HCl) and destabilized (0.8 M Gdn-HCl) IL-1beta, no difference in steady-state (15)N-(1)H NOE enhancements were measured. Thus, the chemical shifts correlate with a global but limited increase in residue flexibility in the presence of Gdn-HCl. Minimalist simulations highlight the regions of greatest position shift between native and 0.8 M Gdn-HCl, which were determined experimentally. This correlation demonstrates that structural changes within the native ensemble of IL-1beta are, at least partially, governed by the principle of minimal energetic frustration.

Amides↗

Stochastic roadmap simulation: an efficient representation and algorithm for analyzing molecular motion.

Classic molecular motion simulation techniques, such as Monte Carlo (MC) simulation, generate motion pathways one at a time and spend most of their time in the local minima of the energy landscape defined over a molecular conformation space. Their high computational cost prevents them from being used to compute ensemble properties (properties requiring the analysis of many pathways). This paper introduces stochastic roadmap simulation (SRS) as a new computational approach for exploring the kinetics of molecular motion by simultaneously examining multiple pathways. These pathways are compactly encoded in a graph, which is constructed by sampling a molecular conformation space at random. This computation, which does not trace any particular pathway explicitly, circumvents the local-minima problem. Each edge in the graph represents a potential transition of the molecule and is associated with a probability indicating the likelihood of this transition. By viewing the graph as a Markov chain, ensemble properties can be efficiently computed over the entire molecular energy landscape. Furthermore, SRS converges to the same distribution as MC simulation. SRS is applied to two biological problems: computing the probability of folding, an important order parameter that measures the "kinetic distance" of a protein's conformation from its native state; and estimating the expected time to escape from a ligand-protein binding site. Comparison with MC simulations on protein folding shows that SRS produces arguably more accurate results, while reducing computation time by several orders of magnitude. Computational studies on ligand-protein binding also demonstrate SRS as a promising approach to study ligand-protein interactions.

Algorithms↗

X-ray diffraction "fingerprinting" of DNA structure in solution for quantitative evaluation of molecular dynamics simulation.

Solution state x-ray diffraction fingerprinting is demonstrated as a method for experimentally assessing the accuracy of molecular dynamics (MD) simulations. Fourier transforms of coordinate data from MD simulations are used to produce reciprocal space "fingerprints" of atomic pair distance correlations that are characteristic of the ensemble and are the direct numerical analogues of experimental solution x-ray diffraction (SXD). SXD experiments and MD simulations were carried out to test the ability of experiment and simulation to resolve sequence-dependent modifications in helix conformation for B-form DNA. SXD experiments demonstrated that solution-state poly(AT) and poly(A)-poly(T) duplex DNA sequences exist in ensembles close to canonical B-form and B'-form structures, respectively. In contrast, MD simulations analyzed in terms of SXD fingerprints are shown to deviate from experiment, most significantly for poly(A)-poly(T) duplex DNA. Compared with experiment, MD simulation shortcomings were found to include both mismatches in simulated conformer structures and number population within the ensembles. This work demonstrates an experimental approach for quantitatively evaluating MD simulations and other coordinate models to simulate biopolymer structure in solution and suggests opportunities to use solution diffraction data as experimental benchmarks for developing supramolecular force fields optimized for a range of in situ applications.

DNA↗