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 343 records · Page 19Linked to original sources

Modeling Alternative Conformational States in CASP16.

The CASP16 Ensemble Prediction experiment assessed advances in methods for modeling proteins, nucleic acids, and their complexes in multiple conformational states. Targets included systems with experimental structures determined in two or three states, evaluated by direct comparison to experimental coordinates, as well as domain-linker-domain (D-L-D) targets assessed against statistical models from NMR and SAXS data. This paper focuses on the former class of multi-state targets. Ten ensembles were released as community challenges, including ligand-induced conformational changes, protein-DNA complexes, a trimeric protein, a stem-loop RNA, and multiple oligomeric states of a single RNA. For five targets, some groups produced reasonably accurate models of both reference states (best TM-score >0.75). However, with the exception of one protein-ligand complex (T1214), where an apo structure was available as a template, predictors generally failed to capture key structural details distinguishing the states. Overall, accuracy was significantly lower than for single-state targets in other CASP experiments. The most successful approaches generated multiple AlphaFold2 models using enhanced multiple sequence alignments and sampling protocols, followed by model quality based selection. While the AlphaFold3 server performed well on several targets, individual groups outperformed it in specific cases. By contrast, predictions for one protein-DNA complex, three RNA targets, and multiple oligomeric RNA states consistently fell short (TM-score <0.75). These results highlight both progress and persistent challenges in multi-state prediction. Despite recent advances, accurate modeling of conformational ensembles, particularly RNA and large multimeric assemblies, remains a critical frontier for structural biology.

AlphaFold2↗

SESAME: a least-squares approach to the evaluation of protein structures computed from NMR data.

A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.

Magnetic Resonance Spectroscopy↗

Single molecule fluorescence studies of surface-adsorbed fibronectin.

The conformation of the extracellular matrix protein fibronectin plays a critical role in regulating cell function, including cell adhesion and migration. While average conformations of large ensembles of adhesion proteins have been previously measured, cells may sensitively respond to conformational outliers. We therefore applied both single molecule imaging and spectroscopy techniques to map a range of conformational states of individual fibronectin molecules adsorbed to glass, as well as to measure their conformational fluctuations in time. Single-step photobleaching experiments confirmed single molecule sensitivity. Single molecule spectra showed fluctuations in the peak wavelength, both over a spatial ensemble of molecules and in a single molecule over time, most likely indicating the different conformational states fibronectin can assume upon surface adsorption. Single-pair fluorescence resonance energy transfer (FRET) revealed that a fraction of fibronectin molecules existed in conformations that allowed for energy transfer between the labeled cysteine residues of the two dimeric arms folded upon each other, and that fluctuations occurred in the FRET efficiency. Fluorescence polarization experiments identified two possible sources of FRET fluctuations: changes in fluorophore orientation and conformational fluctuations of fibronectin over a time scale of seconds.

Adsorption↗

Flexibility in a tetrasaccharide fragment from the high mannose type of N-linked oligosaccharides.

An analysis has been carried out of the three-dimensional structure of a tetrasaccharide, Man(alpha 1-3)Man(alpha 1-6)Man(beta 1-4)GlcN Ac beta 1-OCD3, which is a fragment from the high mannose type of N-linked oligosaccharides. Although earlier work had suggested that this fragment might adopt a stable three-dimensional structure, both n.m.r. and conformational energy calculations support the existence of an ensemble of structures. The conformational entropy calculated from the ensemble and the distribution of distances between the terminal Man(alpha 1-3) and GlcN Ac residues, however, suggests that a significant fraction of the ensemble has the two terminal residues in close proximity.

Carbohydrate Conformation↗

Solution structure of an isolated antibody VL domain.

The solution structure of the isolated VL domain of the anti-digoxin antibody 26-10 has been determined using data derived from heteronuclear multi-dimensional nuclear magnetic resonance (n.m.r.) experiments. Analytical ultracentrifugation and n.m.r. data demonstrate that the VL domain is only weakly associating (Kd = 2.5 (+/- 0.7) mM) and that it experiences a rapid monomer/dimer equilibrium under the n.m.r. experimental conditions. Therefore, the results reported here represent the first structure determination of an antibody VL domain in the absence of fixed quaternary interactions. The structure determination is based on 930 proton-proton distance constraints, 113 dihedral angle constraints, and 46 hydrogen bond constraints. Eighty initial structures were calculated with the variable target function program DIANA; of these, 31 were accepted on the basis of satisfaction of constraints (no distance constraint violations > 0.5 A; target function < 3.0 A2). Accepted DIANA structures were refined by restrained energy minimization using the X-PLOR program. The 15 best energy-minimized DIANA structures were chosen as a representative ensemble of solution conformations. The average root-mean-square differences (r.m.s.d.) between the individual structures of this ensemble and the mean coordinates is 0.85 (+/- 0.10) A for all backbone atoms and 1.29 (+/- 0.10) A for all heavy atoms. For beta-strands A, B, C, D, E and F, the average backbone atom r.m.s.d. to the mean structure is 0.46 (+/- 0.06) A. A higher-resolution ensemble, with all backbone atom and all heavy atom r.m.s.d.s. to the mean coordinates of 0.54 (+/- 0.08) A and 0.98 (+/- 0.12) A, respectively, was obtained by X-PLOR simulated annealing refinement of the 15 energy-minimized DIANA structures. A detailed analysis of the original ensemble of 15 energy-minimized DIANA structures is presented, as this ensemble retains a broader, and possibly more realistic, sampling of conformation space. The backbone atom and all heavy atom r.m.s.d.s between the mean energy-minimized DIANA structure and the X-ray derived coordinates of the VL domain within the Fab/digoxin complex are 1.05 A and 1.56 A, respectively. Subtle differences between the solution and X-ray structures occur primarily in CDR2, CDR3, beta-strands A, F and G, and localized regions of hydrophobic packing. Overall, these results demonstrate that the 26-10 VL domain conformation is determined primarily by intradomain interactions, and that quaternary VL-VH association induces relatively minor conformational adjustments.

Computer Graphics↗

Defining long-range order and local disorder in native alpha-synuclein using residual dipolar couplings.

Intrinsically unstructured proteins play key biochemical roles in a vast range of normal and pathological processes. To study these systems, it is necessary to invoke an ensemble of rapidly interconverting conformations. Residual dipolar couplings (RDCs) are particularly powerful probes of the behavior of unfolded proteins, reporting on time and ensemble-averaged conformations up to and beyond the millisecond time scale. In this study, we present a novel interpretation of RDCs in unfolded systems that simultaneously defines long-range structural order and local conformational sampling. This approach is used to describe the structure and dynamics of alpha-Synuclein (alphaS), a protein that is strongly implicated in the development of Parkinson's disease (PD), allowing unambiguous detection of strongly populated conformers containing long-range contacts between the N- and C-terminal domains. The structural model combines two features required for the description of alphaS in solution: local conformational fluctuation based on random sampling of residue-specific phi/psi distributions, and long-range contacts induced by the presence of nonbonding interactions between domains that are distant in primary sequence. Both aspects are found to be necessary for the reproduction of the nonaveraged RDCs from alphaS. Although RDCs have previously been shown to report on local conformational preferences in unstructured proteins, this study shows the additional sensitivity of these measurements to the presence of long-range order in highly flexible systems.

Protein Conformation↗

Hydrogen exchange and ligand binding: ligand-dependent and ligand-independent protection in the Src SH3 domain.

Amide hydrogen-deuterium exchange has proven to be a powerful tool for detecting and characterizing high-energy conformations in protein ensembles. Since interactions with ligands can modulate these high-energy conformations, hydrogen exchange appears to be an ideal experimental probe of the physical mechanisms underlying processes like allosteric regulation. The chemical mechanism of hydrogen exchange, however, can complicate such studies. Here, we examine hydrogen exchange rates in a simple model system, the c-Src SH3 domain interacting with a short peptide ligand. Addition of ligand slows the rates of hydrogen exchange at nearly every amide for which we can obtain data. Careful analysis, however, reveals that this slowing is due primarily to a reduction in the population of free protein in the system, and not to any specific property of the complex. We present a method to separate the contributions of free and bound protein to the exchange kinetics that has allowed us to identify the subset of amides where exchange arises directly from the complex. These results demonstrate that the slowing of hydrogen exchange induced by ligand interactions should be interpreted with caution, and more extensive experiments are required to correlate changes in hydrogen exchange with changes in structure or internal dynamics.

Animals↗

The MUMO (minimal under-restraining minimal over-restraining) method for the determination of native state ensembles of proteins.

While reliable procedures for determining the conformations of proteins are available, methods for generating ensembles of structures that also reflect their flexibility are much less well established. Here we present a systematic assessment of the ability of ensemble-averaged molecular dynamics simulations with ensemble-averaged NMR restraints to simultaneously reproduce the average structure of proteins and their associated dynamics. We discuss the effects that under-restraining (overfitting) and over-restraining (underfitting) have on the structures generated in ensemble-averaged molecular simulations. We then introduce the MUMO (minimal under-restraining minimal over-restraining) method, a procedure in which different observables are averaged over a different number of molecules. As both over-restraining and under-restraining are significantly reduced in the MUMO method, it is possible to generate ensembles of conformations that accurately characterize both the structure and the dynamics of native states of proteins. The application of the MUMO method to the protein ubiquitin yields a high-resolution structural ensemble with an RDC Q-factor of 0.19.

Amino Acid Sequence↗

Solution conformations and thermodynamics of structured peptides: molecular dynamics simulation with an implicit solvation model.

Calculations of the ensemble of solution conformations and thermodynamics of an analogue of the C-terminal helix of ribonuclease A (RN24) and of a synthetic, beta-hairpin forming peptide (BH8) are presented. For efficient sampling of conformation space, molecular dynamics simulations with an implicit solvent potential and umbrella sampling of the potential energy are performed. Starting from the fully extended chains, the simulations yield several folding and unfolding transitions between disordered (coil) conformations of the peptides and the "native" state (RN24, helix; BH8, hairpin); the simulations also lead to the occurrence of "misfolded" conformations (RN24, hairpin; BH8, helix). In agreement with experiment, the calculations predict 58% helix for RN24 at 275 K and an antiparallel-beta content of 38% at 275 K for BH8; the calculated probabilities for the misfolded species are 2% or smaller at all temperatures considered (250-1100 K). Good agreement is also shown between the calculated 3JHNalpha spin-spin coupling constants of RN24 and BH8 at 275 K, and those obtained from NMR experiments at the same temperature. From the calculated probabilities of helix (h), beta-hairpin (b), and coil (c), the free energy differences between the structured substates are DeltaGch=Gc-Gh approximately 1 kcal/mol and DeltaGbh>/=1.8 kcal/mol for RN24, and DeltaGcb approximately 0.7 kcal/mol and DeltaGhb>/=2.7 kcal/mol for BH8. The free energy difference between "correctly" folded and misfolded secondary structures are of interest for understanding the alpha to beta transition that is thought to play a role in amyloid fibril formation.

Amino Acid Sequence↗

Second-generation mimics of ganglioside GM1 oligosaccharide: a three-dimensional view of their interactions with bacterial enterotoxins by NMR and computational methods.

As a step to delineate a strategy of ligand design for cholera toxin (CT), NMR studies were performed on several mimics of the GM1 ganglioside oligosaccharide. The conformation of these analogues was investigated first in solution and then upon binding to cholera toxin by transferred nuclear Overhauser effect (TR-NOE) measurements. It was demonstrated that CT selects a conformation similar to the global minima of the free saccharides from the ensemble of presented conformations. No evidence of major conformational distortions was obtained, but one or two of the available conformers of the hydroxyacid side chain appear to be selected in the bound state. The NMR data were interpreted with the aid of computer models, generated and analyzed by using a combination of different approaches (MacroModels' MC/EM and MC/SD, Autodock, and GRID). Analysis of the NMR data supported by computational studies allowed us to interpret the experimental observations and to derive workable models of the ligand:toxin complexes. These models suggest that the higher affinity of the (R)-lactic acid derivative 3 may stem from lipophilic interactions with a hydrophobic area in the toxin binding site located in the vicinity of the sialic acid side chain binding region of the CT:GM1 complex, and formed by the side chain of Ile-58 and Lys-34. Thus, the models obtained have allowed us to make useful design suggestions for the improvement of ligand affinity.

Journal Article↗

Second-generation mimics of ganglioside GM1 oligosaccharide: a three-dimensional view of their interactions with bacterial enterotoxins by NMR and computational methods.

As a step to delineate a strategy of ligand design for cholera toxin (CT), NMR studies were performed on several mimics of the GM1 ganglioside oligosaccharide. The conformation of these analogues was investigated first in solution and then upon binding to cholera toxin by transferred nuclear Overhauser effect (TR-NOE) measurements. It was demonstrated that CT selects a conformation similar to the global minima of the free saccharides from the ensemble of presented conformations. No evidence of major conformational distortions was obtained, but one or two of the available conformers of the hydroxyacid side chain appear to be selected in the bound state. The NMR data were interpreted with the aid of computer models, generated and analyzed by using a combination of different approaches (MacroModels' MC/EM and MC/SD, Autodock, and GRID). Analysis of the NMR data supported by computational studies allowed us to interpret the experimental observations and to derive workable models of the ligand:toxin complexes. These models suggest that the higher affinity of the (R)-lactic acid derivative 3 may stem from lipophilic interactions with a hydrophobic area in the toxin binding site located in the vicinity of the sialic acid side chain binding region of the CT:GM1 complex, and formed by the side chain of Ile-58 and Lys-34. Thus, the models obtained have allowed us to make useful design suggestions for the improvement of ligand affinity.

Binding Sites↗

Molecular docking to ensembles of protein structures.

Until recently, applications of molecular docking assumed that the macromolecular receptor exists in a single, rigid conformation. However, structural studies involving different ligands bound to the same target biomolecule frequently reveal modest but significant conformational changes in the target. In this paper, two related methods for molecular docking are described that utilize information on conformational variability from ensembles of experimental receptor structures. One method combines the information into an "energy-weighted average" of the interaction energy between a ligand and each receptor structure. The other method performs the averaging on a structural level, producing a "geometry-weighted average" of the inter-molecular force field score used in DOCK 3.5. Both methods have been applied in docking small molecules to ensembles of crystal and solution structures, and we show that experimentally determined binding orientations and computed energies of known ligands can be reproduced accurately. The use of composite grids, when conformationally different protein structures are available, yields an improvement in computational speed for database searches in proportion to the number of structures.

Computer Simulation↗

Generating ensemble averages for small proteins from extended conformations by Monte Carlo simulations.

Since it is not feasible to determine the structure of every protein by experiment, algorithms delivering the folded conformation of a protein solely from its amino acid sequence are desirable. Here the diffusion-process controlled-Monte Carlo approach has been applied to generating ensemble averages for three small proteins with 31, 36, and 46 residues. Starting from extended conformations and using an energy model that was developed on other protein models, the simulations find nativelike structures deviating by 3 A rms from the experimental structures for the main chain atoms. The balance between long-range and short-range interactions is discussed briefly in the context of stability and folding.

Computer Simulation↗

Factors modulating conformational equilibria in large modular proteins: a case study with cobalamin-dependent methionine synthase.

In the course of catalysis or signaling, large multimodular proteins often undergo conformational changes that reposition the modules with respect to one another. The mechanisms that direct the reorganization of modules in these proteins are of considerable importance, but distinguishing alternate conformations is a challenge. Cobalamin-dependent methionine synthase (MetH) is a 136-kDa multimodular enzyme with a cobalamin chromophore; the color of the cobalamin reflects the conformation of the protein. The enzyme contains four modules and catalyzes three different methyl transfer reactions that require different arrangements of these modules. Two of these methyl transfer reactions occur during turnover, when homocysteine is converted to methionine by using a methyl group derived from methyltetrahydrofolate. The third reaction is occasionally required for reactivation of the enzyme and uses S-adenosyl-L-methionine as the methyl donor. The absorbance properties of the cobalamin cofactor have been exploited to assign conformations of the protein and to probe the effect of ligands and mutations on the distribution of conformers. The results imply that the methylcobalamin form of MetH exists as an ensemble of interconverting conformational states. Differential binding of substrates or products alters the distribution of conformers. Furthermore, steric conflicts disfavor conformers that juxtapose a methyl group on substrate with one on methylcobalamin. These results suggest that the methylation state of the cobalamin will influence the distribution of conformers during turnover.

5-Methyltetrahydrofolate-Homocysteine S-Methyltran↗

Dynamics in the unfolded state of beta2-microglobulin studied by NMR.

Many proteins form amyloid-like fibrils in vitro under conditions that favour the population of partially folded conformations or denatured state ensembles. Characterising the structural and dynamic properties of these states is crucial towards understanding the mechanisms of self-assembly in amyloidosis. The aggregation of beta2-microglobulin (beta2m) into amyloid fibrils in vivo occurs in the condition known as dialysis-related amyloidosis (DRA) and the protein has been shown to form amyloid-like fibrils under acidic conditions in vitro. We have used a number of 1H-15N nuclear magnetic resonance (NMR) experiments in conjunction with site-directed mutagenesis to study the acid-unfolded state of beta2m. 15N NMR transverse relaxation experiments reveal that the acid-denatured ensemble, although predominantly unfolded at the N and C termini, contains substantial non-native structure in the central region of the polypeptide chain, stabilised by long-range interactions between aromatic residues and by the single disulphide bond. Relaxation dispersion studies indicate that the acid-unfolded ensemble involves two or more distinct species in conformational equilibrium on the micro- to millisecond time-scale. One of these species appears to be hydrophobically collapsed, as mutations in an aromatic-rich region of the protein, including residues that are solvent-exposed in the native protein, disrupt this structure and cause a consequent decrease in the population of this conformer. Thus, acid-unfolded beta2m consists of a heterogeneous ensemble of rapidly fluctuating species, some of which contain stable, non-native hydrophobic clusters. Given that amyloid assembly of beta2m proceeds with lag kinetics under the conditions of this study, a rarely populated species such as a conformer with non-native aromatic clustering could be key to the initiation of amyloidosis.

Acids↗

Modeling protein loops using a phi i + 1, psi i dimer database.

We present an automated method for modeling backbones of protein loops. The method samples a database of phi i + 1 and psi i angles constructed from a nonredundant version of the Protein Data Bank (PDB). The dihedral angles phi i + 1 and psi i completely define the backbone conformation of a dimer when standard bond lengths, bond angles, and a trans planar peptide configuration are used. For the 400 possible dimers resulting from 20 natural amino acids, a list of allowed phi i + 1, psi i pairs for each dimer is created by pooling all such pairs from the loop segments of each protein in the nonredundant version of the PDB. Starting from the N-terminus of the loop sequence, conformations are generated by assigning randomly selected pairs of phi i + 1, psi i for each dimer from the respective pool using standard bond lengths, bond angles, and a trans peptide configuration. We use this database to simulate protein loops of lengths varying from 5 to 11 amino acids in five proteins of known three-dimensional structures. Typically, 10,000-50,000 models are simulated for each protein loop and are evaluated for stereochemical consistency. Depending on the length and sequence of a given loop, 50-80% of the models generated have no stereochemical strain in the backbone atoms. We demonstrate that, when simulated loops are extended to include flanking residues from homologous segments, only very few loops from an ensemble of sterically allowed conformations orient the flanking segments consistent with the protein topology. The presence of near-native backbone conformations for loops from five different proteins suggests the completeness of the dimeric database for use in modeling loops of homologous proteins. Here, we take advantage of this observation to design a method that filters near-native loop conformations from an ensemble of sterically allowed conformations. We demonstrate that our method eliminates the need for a loop-closure algorithm and hence allows for the use of topological constraints of the homologous proteins or disulfide constraints to filter near-native loop conformations.

Amino Acid Sequence↗

Long- and short-range interactions in native protein structures are consistent/minimally frustrated in sequence space.

We show that long- and short-range interactions in almost all protein native structures are actually consistent with each other for coarse-grained energy scales; specifically we mean the long-range inter-residue contact energies and the short-range secondary structure energies based on peptide dihedral angles, which are potentials of mean force evaluated from residue distributions observed in protein native structures. This consistency is observed at equilibrium in sequence space rather than in conformational space. Statistical ensembles of sequences are generated by exchanging residues for each of 797 protein native structures with the Metropolis method. It is shown that adding the other category of interaction to either the short- or long-range interactions decreases the means and variances of those energies for essentially all protein native structures, indicating that both interactions consistently work by more-or-less restricting sequence spaces available to one of the interactions. In addition to this consistency, independence by these interaction classes is also indicated by the fact that there are almost no correlations between them when equilibrated using both interactions and significant but small, positive correlations at equilibrium using only one of the interactions. Evidence is provided that protein native sequences can be regarded approximately as samples from the statistical ensembles of sequences with these energy scales and that all proteins have the same effective conformational temperature. Designing protein structures and sequences to be consistent and minimally frustrated among the various interactions is a most effective way to increase protein stability and foldability.

Amino Acid Sequence↗

Concerted motions in HIV-1 TAR RNA may allow access to bound state conformations: RNA dynamics from NMR residual dipolar couplings.

Ground-state dynamics in RNA is a critical precursor for structural adaptation observed ubiquitously in protein-RNA recognition. A tertiary conformational analysis of the stem-loop structural element in the transactivation response element (TAR) from human immunodeficiency virus type 1 (HIV-I) RNA is presented using recently introduced NMR methods that rely on the measurement of residual dipolar couplings (RDC) in partially oriented systems. Order matrix analysis of RDC data provides evidence for inter-helical motions that are of amplitude 46(+/-4) degrees, of random directional character, and that are executed about an average conformation with an inter-helical angle between 44 degrees and 54 degrees. The generated ensemble of TAR conformations have different organizations of functional groups responsible for interaction with the trans-activator protein Tat, including conformations similar to the previously characterized bound-state conformation. These results demonstrate the utility of RDC-NMR for simultaneously characterizing RNA tertiary dynamics and average conformation, and indicate an avenue for TAR complex formation involving tertiary structure capture.

Gene Products, tat↗