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 19 recordsLinked to original sources
PubMed2000

Computational analysis of two similar neuropeptides yields distinct conformational ensembles.

Conformational states and thermodynamic properties for two similar neuropeptides, GDPFLRF-NH(2) and GYPFLRF-NH(2), have been computed by Monte Carlo simulated annealing (MCSA) conformational searches and Metropolis Monte Carlo (MMC) calculations. These peptides were recently shown to have dramatically different conformations in solution by NMR [Edison et al., J Neuroscience 1999;19:6318-6326]. Final conformations of multiple independent MCSA runs were the starting points for MMC calculations, and conformations saved at intervals during MMC runs were characterized in terms of total energy, configuration entropy, side-chain fraction population, and ensemble average inter-nuclear distances. Without the use of any NMR data-generated pseudo-potentials, the present calculations were in excellent qualitative agreement with all previous NMR experimental data and provided a foundation by which to more quantitatively interpret the experimental NMR results. Proteins 2000;40:367-377.

Computer Simulation↗
PubMed1999

The native state conformational ensemble of the SH3 domain from alpha-spectrin.

The folding/unfolding equilibrium of the alpha-spectrin SH3 domain has been measured by NMR-detected hydrogen/deuterium exchange and by differential scanning calorimetry. Protection factors against exchange have been obtained under native conditions for more than half of the residues in the domain. Most protected residues are located at the beta-strands, the short 3(10) helix, and part of the long RT loop, whereas the loops connecting secondary structure elements show no measurable protection. Apparent stability constants per residue and their corresponding Gibbs energies have been calculated from the exchange experiments. The most stable region of the SH3 domain is defined by the central portions of the beta-strands. The peptide binding region, on the other hand, is composed of a highly stable region (residues 53-57) and a highly unstable region, the loop between residues 34-41 (n-Src loop). All residues in the domain have apparent Gibbs energies lower than the global unfolding Gibbs energy measured by differential scanning calorimetry, indicating that under our experimental conditions the amide exchange of all residues in the SH3 domain occurs primarily via local unfolding reactions. A structure-based thermodynamic analysis has allowed us to predict correctly the thermodynamics of the global unfolding of the domain and to define the ensemble of conformational states that quantitatively accounts for the observed pattern of hydrogen exchange protection. These results demonstrate that under native conditions the SH3 domain needs to be considered as an ensemble of conformations and that the hydrogen exchange data obtained under those conditions cannot be interpreted by a two-state equilibrium. The observation that specific regions of a protein are able to undergo independent local folding/unfolding reactions indicates that under native conditions the scale of cooperative interactions is regional rather than global.

Calorimetry, Differential Scanning↗
PubMed1999

Conformational ensembles: the role of neuropeptide structures in receptor binding.

Conformational properties of several similar FMRFamide-like neuropeptides from mollusks were investigated by nuclear magnetic resonance (NMR) spectroscopy. It was found that amino acid substitutions in the N-terminal variable regions of the peptides had dramatic effects on the populations of reverse turns in solution. The populations of turns, as measured by two independent NMR parameters, were found to be highly correlated (r(2) = 0.93 and 0. 82) with IC(50) values using receptor membrane preparations from Helix aspersa (Payza, 1987; Payza et al., 1989). These results suggest that the amount of turn in the free peptide can influence the receptor binding affinities of that peptide. On the basis of these observations, a model was developed in which only a single species from a conformational ensemble of an unbound peptide will bind to a particular receptor. Thus, the conformational ensemble reduces the effective concentration of a particular peptide with respect to a particular receptor.

Amino Acid Sequence↗
PubMed2025

Pioneer in Molecular Biology: Conformational Ensembles in Molecular Recognition, Allostery, and Cell Function.

In 1978, for my PhD, I developed the efficient O(n3) dynamic programming algorithm for the-then open problem of RNA secondary structure prediction. This algorithm, now dubbed the "Nussinov algorithm", "Nussinov plots", and "Nussinov diagrams", is still taught across Europe and the U.S. As sequences started coming out in the 1980s, I started seeking genome-encoded functional signals, later becoming a bioinformatics trend. In the early 1990s I transited to proteins, co-developing a powerful computer vision-based docking algorithm. In the late 1990s, I proposed the foundational role of conformational ensembles in molecular recognition and allostery. At the time, conformational ensembles and free energy landscapes were viewed as physical properties of proteins but were not associated with function. The classical view of molecular recognition and binding was based on only two conformations captured by crystallography: open and closed. I proposed that all conformational states preexist. Proteins always have not one folded form-nor two-but many folded forms. Thus, rather than inducing fit, binding can work by shifting the ensembles between states, and this shifting, or redistributing the ensembles to maintain equilibrium, is the origin of the allosteric effect and protein, thus cell, function. This transformative paradigm impacted community views in allosteric drug design, catalysis, and regulation. Dynamic conformational ensemble shifts are now acknowledged as the origin of recognition, allostery, and signaling, underscoring that conformational ensembles-not proteins-are the workhorses of the cell, pioneering the fundamental idea that dynamic ensembles are the driving force behind cellular processes. Nussinov was recognized as pioneer in molecular biology by JMB.

Molecular Biology↗
PubMed2005

A perspective on how the Vitamin D sterol/Vitamin D receptor (VDR) conformational ensemble model can potentially be used to understand the structure-function results of A-ring modified Vitamin D sterols.

The steroid hormone 1alpha,25(OH)(2)-Vitamin D(3) (1,25D) activates both genomic and non-genomic intracellular signaling cascades. It is also well recognized that co-incubation of 1,25D with its C-1 epimer, 1beta,25D (HL), suppresses the efficiency of the non-genomic signal activated by 1,25D alone and that its C-3 epimer, 3alpha-1,25D (HJ) is nearly as potent as 1,25D in suppressing PTH secretion, believed to be propagated by 1,25D's genomic signaling. Both these sterols lack the hypercalcemic effect induced by pharmacological doses of 1,25D and have reduced VDR affinity compared to 1,25D, as measured in a steroid competition assay. Recent functional studies suggest that the VDR is required for both non-genomic and genomic signaling. Along these lines we have recently proposed a Vitamin D sterol/VDR conformational ensemble model that posits the VDR contains two distinct, yet overlapping ligand binding sites, and that the potential differential stabilities of 1,25D and HL in these two pockets can be used to explain their different non-genomic signaling properties. The overlapping region is predominantly occupied by the sterol's A-ring when it is bound to either the genomic ligand binding pocket (G-pocket), defined by X-ray crystallography, or the alternative ligand binding pocket (A-pocket), discovered using in silico techniques (directed docking). Therefore, to gain further insight into the potential application of this model we docked the other A-ring diastereomer [(1beta,3alpha)=HH] of 1,25D and its 1- and 3-deoxy forms (25D and CF, respectively) to the A- and G-pockets to assess their potential stabilities in the pockets, relative to 1,25D. The models were then used to provide putative mechanistic arguments for their known structure-function experimental results. This model may provide new insights into how Vitamin D sterols that uncouple the unwanted hypercalcemic effect from attractive growth inhibitory/differentiation properties can do so by differentially stabilizing different subpopulations of VDR conformational ensemble members.

Animals↗
PubMed2026

Quantifying the peripheral surface information entropy from conformational ensembles of globular protein-peptide complexes.

Predicting favorable protein-peptide binding events remains a central challenge in biophysics, with continued uncertainty surrounding how nonlocal effects shape the global energy landscape. Here, we introduce peripheral surface information entropy, SΨ, a quantitative measure of the statistical variability in apolar and charged non-interacting surface (NIS) proportions across conformational ensembles. Within the Gibbs free-energy relation ΔG = ΔH - TΔS, SΨ is proposed as a computationally tractable entropic proxy rather than a direct thermodynamic observable or stand-alone estimator of binding affinity. Using energy-directed molecular docking via HADDOCK3 and explicit-solvent molecular dynamics simulations, it is demonstrated that favorable binding partners exhibit emergent, low-entropy N-states (discrete macrostates in NIS state space) indicative of preferential apolar/charged surface configurations. Across dozens of peptides and multiple receptor systems (WW, PDZ, and MDM2 domains), dominant N-states persisted under varied docking parameters and initial conditions. A meta-ensemble of 657 complexes from 36 experiments over 15 years confirmed the presence of dominant NIS modes independent of in silico methodology, suggesting an evolutionary selection pressure toward specific NIS fingerprints. These findings establish SΨ as a thermoinformatic descriptor that encodes favorable binding constraints into unique statistical signatures of the NIS.

Entropy↗
PubMed2006

Characterizing the conformational ensemble of monomeric polyglutamine.

Studies of synthetic polyglutamine peptides in vitro have established that polyglutamine peptides aggregate via a classic nucleation and growth mechanism. Chen and colleagues [Proc Natl Acad Sci U S A 2002;99:11884-11889] have found that monomeric polyglutamine, which is a disordered statistical coil in solution, is the critical nucleus for aggregation. Therefore, nucleation of beta-sheet-rich aggregates requires an initial disorder to order transition, which is a highly unfavorable thermodynamic reaction. The questions of interest to us are as follows: What are the statistical fluctuations that drive beta-sheet formation in monomeric polyglutamine? How do these fluctuations vary with chain length? And why is this process thermodynamically unfavorable, that is, why is monomeric polyglutamine disordered? To answer these questions we use multiple molecular dynamics simulations to provide quantitative characterization of conformational ensembles for two short polyglutamine peptides. We find that the ensemble for polyglutamine is indeed disordered. However, the disorder is inherently different from that of denatured proteins and the average compactness and magnitude of conformational fluctuations increase with chain length. Most importantly, the effective concentration of sidechain primary amides around backbone units is inherently high and peptide units are solvated either by hydrogen bonds to sidechains or surrounding water molecules. Due to the multiplicity of backbone solvation modes the probability associated with any specific backbone conformation is small, resulting in a conformational entropy bottleneck which makes beta-sheet formation in monomeric polyglutamine thermodynamically unfavorable.

Algorithms↗
PubMed2005

Applications of the Vitamin D sterol-Vitamin D receptor (VDR) conformational ensemble model.

Over the past 20 years much has been learned about the cellular actions of the steroid hormone 1alpha,25(OH)2-Vitamin D3 (1,25D). Perhaps most importantly structure-function studies led to the discovery that different chemical and physical features of 1,25D are preferred to initiate either exonuclear, non-genomic or endonuclear, genomic cellular signaling. It is well documented that both a 1alpha-OH and 25-OH, and a 6-s-trans, bowl-shaped, sterol conformation are absolutely required for efficient gene transcription, while 6-s-cis locked analogs and 1-deoxy, 25(OH)D3 metabolites activate a variety of non-genomic, rapid responses. These results and the observation that S237 (helix-3; H3) and R274 (H5) are the most static residues in the human 1,25D-Vitamin D receptor (VDR) X-ray construct (see B-values in pdb: 1DB1) and form H-bonds with the 1alpha-OH of 1,25D in the X-ray, genomic pocket (G-pocket), provided the basis for the molecular modeling experiments that led to the discovery of a putative VDR alternative ligand binding pocket (A-pocket). The conformational ensemble model generated from the in silico results provides an explanation for how the VDR can function as a receptor propagating both genomic and non-genomic signaling events. In this report the theoretical gating properties controlling ligand access to the A- and G-pockets will be compared and the model will be used to provide a molecular explanation for the confusing structure-function results pertaining to 1,25D, its side-chain metabolite, 23S,25R-1alpha,25(OH)2-D3-26,23-lactone (BS), and its synthetic two side-chain analog, 21-(3'-hydroxy-3'-methylbutyl)-1alpha,25(OH)2-D3 (KH or Gemini). In addition, evidence that the model is consistent with the pH requirement for Vitamin D sterol-VDR crystallization will be presented.

Arginine↗
PubMed1999

NMR characterization of partially folded and unfolded conformational ensembles of proteins.

Studies of unfolded and partially folded proteins provide important insight into the initiation and process of protein folding. This review focuses on the use of nmr in characterization of ensembles of proteins that model the early stages of folding. Analysis of an ensemble includes description of the number of conformations, their structure, relative populations, interconversion rates, and dynamics of subconformations. A chemically synthesized analogue of bovine pancreatic trypsin inhibitor (BPTI), [14-38](Abu), has provided a rare system for characterization of multiple partially folded conformations in slow exchange at near physiological conditions. Multidimensional nmr techniques coupled with selective labeling were used to probe different segments of the polypeptide chain. At each labeled site, there is evidence of slow interconversion between two families of partially folded conformations that in themselves are ensembles of rapidly interconverting conformers. All these conformers display significantly more order in the core relative to the rest of the molecule. For other variants of BPTI that are unfolded at equilibrium, the most ordered structure is also favored in the hydrophobic core residues of the native protein. This is consistent with the hypothesis that the residues that are the first to fold go on to form the most stable, structure-determining part of the protein.

Animals↗
PubMed2004

Co-populated conformational ensembles of beta2-microglobulin uncovered quantitatively by electrospray ionization mass spectrometry.

Ordered assembly of monomeric human beta(2)-microglobulin (beta(2)m) into amyloid fibrils is associated with the disorder hemodialysis-related amyloidosis. Previously, we have shown that under acidic conditions (pH <5.0 at 37 degrees C), wild-type beta(2)m assembles spontaneously into fibrils with different morphologies. Under these conditions, beta(2)m populates a number of different conformational states in vitro. However, this equilibrium mixture of conformationally different species is difficult to resolve using ensemble techniques such as nuclear magnetic resonance or circular dichroism. Here we use electrospray ionization mass spectrometry to resolve different species of beta(2)m populated between pH 6.0 and 2.0. We show that by linear deconvolution of the charge state distributions, the extent to which each conformational ensemble is populated throughout the pH range can be determined and quantified. Thus, at pH 3.6, conditions under which short fibrils are produced, the conformational ensemble is dominated by a charge state distribution centered on the 9+ ions. By contrast, under more acidic conditions (pH 2.6), where long straight fibrils are formed, the charge state distribution is dominated by the 10+ and 11+ ions. The data are reinforced by investigations on two variants of beta(2)m (V9A and F30A) that have reduced stability to pH denaturation and show changes in the pH dependence of the charge state distribution that correlate with the decrease in stability measured by tryptophan fluorescence. The data highlight the potential of electrospray ionization mass spectrometry to resolve and quantify complex mixtures of different conformational species, one or more of which may be important in the formation of amyloid.

Amino Acid Substitution↗
PubMed1995

Probability assessment of conformational ensembles: sugar repuckering in a DNA duplex in solution.

Conformational flexibility of molecules in solution implies that different conformers contribute to the NMR signal. This may lead to internal inconsistencies in the 2D NOE-derived interproton distance restraints and to conflict with scalar coupling-based torsion angle restraints. Such inconsistencies have been revealed and analyzed for the DNA octamer GTATAATG.CATATTAC, containing the Pribnow box consensus sequence. A number of subsets of distance restraints were constructed and used in the restrained Monte Carlo refinement of different double-helical conformers. The probabilities of conformers were then calculated by a quadratic programming algorithm, minimizing a relaxation rate-base residual index. The calculated distribution of conformers agrees with the experimental NOE data as an ensemble better than any single structure. A comparison with the results of this procedure, which we term PARSE (Probability Assessment via Relaxation rates of a Structural Ensemble), to an alternative method to generate solution ensembles showed, however, that the detailed multi-conformational description of solution DNA structure remains ambiguous at this stage. Nevertheless, some ensemble properties can be deduced with confidence, the most prominent being a distribution of sugar puckers with minor populations in the N-region and major populations in the S-region. Importantly, such a distribution is in accord with the analysis of independent experimental data--deoxyribose proton-proton scalar coupling constants.

Algorithms↗
PubMed2006

Modeling protein conformational ensembles: from missing loops to equilibrium fluctuations.

Characterizing protein flexibility is an important goal for understanding the physical-chemical principles governing biological function. This paper presents a Fragment Ensemble Method to capture the mobility of a protein fragment such as a missing loop and its extension into a Protein Ensemble Method to characterize the mobility of an entire protein at equilibrium. The underlying approach in both methods is to combine a geometric exploration of conformational space with a statistical mechanics formulation to generate an ensemble of physical conformations on which thermodynamic quantities can be measured as ensemble averages. The Fragment Ensemble Method is validated by applying it to characterize loop mobility in both instances of strongly stable and disordered loop fragments. In each instance, fluctuations measured over generated ensembles are consistent with data from experiment and simulation. The Protein Ensemble Method captures the mobility of an entire protein by generating and combining ensembles of conformations for consecutive overlapping fragments defined over the protein sequence. This method is validated by applying it to characterize flexibility in ubiquitin and protein G. Thermodynamic quantities measured over the ensembles generated for both proteins are fully consistent with available experimental data. On these proteins, the method recovers nontrivial data such as order parameters, residual dipolar couplings, and scalar couplings. Results presented in this work suggest that the proposed methods can provide insight into the interplay between protein flexibility and function.

Algorithms↗
PubMed1998

Dynamics of the conformational ensemble of partially folded bovine pancreatic trypsin inhibitor.

A single-disulfide variant of bovine pancreatic trypsin inhibitor (BPTI), [14-38]Abu, is a partially folded ensemble which includes two, and in one case three, conformations that interconvert slowly enough to exhibit separate cross-peaks in the amide region of homonuclear and heteronuclear NMR spectra. Each conformation is itself composed of many subconformations in rapid equilibrium. Partially folded BPTI undergoes local motions that are slow, noncooperative, independent fluctuations of short segments within the chain. Cooperative global unfolding of the ensemble is also observed. Heteronuclear NMR has been used to measure interconversion rate constants of partially folded conformational substates; the rate constants differ for each residue and vary over an order of magnitude. For local fluctuation, the forward rate constants for amide protons of the antiparallel beta-sheet are significantly smaller than the rest of the molecule, consistent with other indications that this is the most stable part of the partially folded protein. The reverse rate constants also vary; they are the highest for Ala 27 in the turn between the strands in the sheet and for Phe 33 in the antiparallel beta-sheet. Global unfolding interconversion rate constants vary over a 3-fold range, consistent with previously observed deviations from two-state behavior. Fast backbone dynamics, from T1, T2, and NOE relaxation parameters, are obtained for the slowly interchanging conformations in the partially folded ensemble. Clear differences are observed between the two conformations; one is more flexible and less compact than the other. In the more flexible and disordered partially folded conformation, intermediate exchange is detected for some backbone amides, namely, those in the central beta-sheet and the turn. These same sheet and turn residues are more ordered in the globally denatured ensemble as well. Our results suggest that the turn initiates formation of a partially folded ensemble in which the slow-exchange core is the most stable region and in which segmental fluctuations reflect multiple nuclei for folding of the rest of the molecule.

Animals↗
PubMed1998

Flexible ligand docking using conformational ensembles.

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.

Algorithms↗
PubMed1990

Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins.

We present a prototype of a new approach to the folding problem of polypeptide chains. This approach is based on the analysis of known protein structures. It derives the energy potentials for the atomic interactions of all amino acid residue pairs as a function of the distance between the involved atoms. These potentials are then used to calculate the energies of all conformations that exist in the data base with respect to a given sequence. Then, by using only the most stable conformations, clusters of the most probable conformations for the given sequence are obtained. To discuss the results properly we introduce a new classification of segments based on their conformational stability. Special care is taken to allow for sparse data sets. The use of the method is demonstrated in the discussion of the identical oligopeptide sequences found in different conformations in unrelated proteins. VNTFV, for example, adopts a beta-strand in ribonuclease but it is found in an alpha-helical conformation in erythrocruorin. In the case of VNTFV the ensemble obtained consists of a single cluster of beta-strand conformations, indicating that this may be the preferred conformation for the pentapeptide. When the flanking residues are included in the calculation the hepapeptide P-VNTFV-H (ribonuclease) again yields an ensemble of beta-strands. However, in the ensemble of D-VNTFV-A (erythrocruorin) the major cluster is of alpha-helical type. In the present study we concentrate on the local aspects of protein conformations. However, the theory presented is quite general and not restricted to oligopeptides. We indicate extensions of the approach to the calculation of global conformations of proteins as well as conceivable applications to a number of molecular systems.

Amino Acid Sequence↗
PubMed2001

Detection of multiple protein conformational ensembles in solution via deconvolution of charge-state distributions in ESI MS.

Monitoring the changes in charge-state distributions of protein ions in electrospray ionization (ESI) mass spectra has become one of the commonly accepted tools to detect large-scale conformational changes of proteins in solution. However, these experiments produce only qualitative, low-resolution information. Our goal is to develop a procedure that would produce quantitative data on protein conformational isomers coexisting in solution at equilibrium. To that end, we have examined the evolution of positive ion charge-state distributions in the

Animals↗