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 757 records · Page 42Linked to original sources

Unraveling the symmetry ambiguity in a hexamer: calculation of the R6 human insulin structure.

Crystallographic and NMR studies of insulin have revealed a highly flexible molecule with a range of different aggregation and structural states; the importance of these states for the function of the hormone is still unclear. To address this question, we have studied the solution structure of the insulin R6 symmetric hexamer using NMR spectroscopy. Structure determination of symmetric oligomers by NMR is complicated due to 'symmetry ambiguity' between intra- and intermonomer NOEs, and between different classes of intermonomer NOEs. Hence, to date, only two symmetric tetramers and one symmetric pentamer (VTB, B subunit of verotoxin) have been solved by NMR: there has been no other symmetric hexamer or higher-order oligomer. Recently, we reported a solution structure for R6 insulin hexamer. However, in that study, a crystal structure was used as a reference to resolve ambiguities caused by the threefold symmetry; the same method was used in solving VTB. Here, we have successfully recalculated R6 insulin using the symmetry-ADR method, a computational strategy in which ambiguities are resolved using the NMR data alone. Thus the obtained structure is a refinement of the previous R6 solution structure. Correlated motions in the final structural ensemble were analysed using a recently developed principal component method; this suggests the presence of two major conformational substates. The study demonstrates that the solution structure of higher-order symmetric oligomers can be determined unambiguously from NMR data alone, using the symmetry-ADR method. This success bodes well for future NMR studies of higher-order symmetric oligomers. The correlated motions observed in the structural ensemble suggest a new insight into the mechanism of phenol exchange and the T6 <--> R6 transition of insulin in solution.

Binding Sites↗

Conformational analysis of protein structures derived from NMR data.

A study is presented of the conformational characteristics of NMR-derived protein structures in the Protein Data Bank compared to X-ray structures. Both ensemble and energy-minimized average structures are analyzed. We have addressed the problem using the methods developed for crystal structures by examining the distribution of phi, psi, and chi angles as indicators of global conformational irregularity. All these features in NMR structures occur to varying degrees in multiple conformational states. Some measures of local geometry are very tightly constrained by the methods used to generate the structure, e.g., proline phi angles, alpha-helix phi,psi angles, omega angles, and C alpha chirality. The more lightly restrained torsion angles do show increased clustering as the number of overall experimental observations increases. phi, psi, and chi 1 angle conformational heterogeneity is strongly correlated with accessibility but shows additional differences which reflect the differing number of observations possible in NMR for the various side chains (e.g., many for Trp, few for Ser). In general, we find that the core is defined to a notional resolution of 2.0 to 2.3 A. Of real interest is the behavior of surface residues and in particular the side chains where multiple rotameric states in different structures can vary from 10% to 88%. Later generation structures show a much tighter definition which correlates with increasing use of J-coupling information, stereospecific assignments, and heteronuclear techniques. A suite of programs is being developed to address the special needs of NMR-derived structures which will take into account the existence of increased mobility in solution.

Analysis of Variance↗

Efficient light harvesting in dye-endcapped conjugated polymers probed by single molecule spectroscopy.

The development of sophisticated microscopic models of energy transfer in linear multichromophoric systems such as conjugated polymers is rarely matched by suitable experimental studies on the microscopic level. To assess the roles of structural, temporal, and energetic disorder in energy transfer, single molecule spectroscopic investigations of the elementary processes leading to energetic relaxation in conjugated polymers are desirable. We present a detailed study of energy transfer processes occurring in dye-endcapped conjugated polymer molecules on the single molecule level. These processes are mostly masked in ensemble investigations. Highly efficient intramolecular energy transfer along a single polyindenofluorene chain to a perylene endcap occurs in many instances and is resolved in real time. We further consider the spectral emission characteristics of the single molecule, the polarization anisotropy which reveals the chain conformation, the fluorescence intermittency, and the temperature dependence and conclude that the efficiency of energy transfer in the ensemble is controlled by the statistics of the individual molecules. The weak thermal activation of energy transfer indicates the involvement of vibrational modes in interchromophoric coupling. Whereas backbone-endcap coupling is strong, the rate-limiting step for intramolecular energy transfer is the migration along the backbone. The results are particularly relevant to understanding undesired exciton trapping on fluorenone defects in polyfluorenes.

Journal Article↗

Conformation and solvation structure for an isolated n-octadecane chain in water, methanol, and their mixtures.

Configurational-bias Monte Carlo simulations in the isobaric-isothermal ensemble (T = 323 K and p = 10 atm) were carried out to probe structural properties of an isolated n-octadecane chain solvated in water, methanol, water-rich, or methanol-rich mixtures and, for comparison, of an isolated chain in the gas phase and for neat liquid n-octadecane. The united-atom version of the TraPPE (transferable potentials for phase equilibria) force field was used to represent n-octadecane and methanol and the TIP-4P model was used for water. In all six environments, broad conformational distributions are observed and the n-octadecane chains are found to predominantly adopt extended, but not all-trans conformations. In addition, a small fraction of more collapsed conformations in which the chain ends approach each other is observed for aqueous hydration, the water-rich solvent mixture and the gas phase, but the simulation data do not support a simple two-state picture with folded and unfolded basins of attraction. For chains in these three "poor" solvent environments, the dihedral angles near the center of the chain show an enhancement of the gauche population. The ensemble of water-solvated chains with end-to-end contacts is preferentially found in a U-shaped conformation rather than a more globular state. An analysis of the local solvation structures in the water-methanol mixtures shows, as expected, an enrichment of the methyl group of methanol near the methylene and methyl segments of the n-octadecane chain. Interestingly, these local bead fractions are enhanced by factors of 2.5 and 1.5 for methyl and methylene segments reflecting the more hydrophobic nature of the former segments.

Alkanes↗

Structure and backbone dynamics of Apo-CBFbeta in solution.

Runx proteins constitute a family of mammalian transcription factors that interact with DNA through their evolutionarily conserved Runt domain. CBFbeta, alternatively denoted PEBP2beta, is the non-DNA-binding heterodimer partner and acts to enhance the DNA binding affinity of Runx proteins. Runx proteins and CBFbeta are associated with a variety of biological functions and human diseases; they are, for example, together the most frequent targets for chromosomal rearrangements in acute human leukemias. We have determined the solution structure and characterized the backbone dynamics of C-terminally truncated fragments containing residues 1-141 of CBFbeta. The present apo-CBFbeta structure is very similar to that seen in a Runt-CBFbeta complex. An evaluation of backbone (15)N NMR relaxation parameters shows that CBFbeta is a rigid molecule with high order parameters throughout the backbone; the only regions displaying significant dynamics are a long loop and the C-terminal alpha-helix. A few residues display relaxation behavior indicating conformational exchange on microsecond to millisecond time scales, but only one of these is located at the Runt binding surface. Our structure and dynamics analysis of CBFbeta therefore suggests that the protein binds to Runt without large conformational changes or induced folding ("lock-and-key" interaction). The apo-CBFbeta structure presented here exhibits several significant differences with two other published NMR ensembles of very similar protein fragments. The differences are located in four regions outside of the central beta-barrel, whereas the beta-barrel itself is almost identical in the three NMR structures. The comparison illustrates that independently determined NMR structures may display rather large differences in backbone conformation in regions that appear to be well-defined in each of the calculated NMR ensembles.

Animals↗

NMR-derived model of interconverting conformations of an ICAM-1 inhibitory cyclic nonapeptide.

We have produced by phage-display a disulfide-linked cyclic nonapeptide (inhibitory peptide-01, IP01), CLLRMRSIC, that binds to intracellular adhesion molecule-1 (ICAM-1) and blocks binding to its counter-structure, leukocyte functional antigen-1 (LFA-1). As a first step towards improving its pharmacologic properties, we have performed a structural and functional analysis of this peptide inhibitor to determine the features relevant to ICAM-1 binding. We report here the solution model of our initial product, IP01, as derived from two-dimensional nuclear magnetic resonance (NMR) restraints and molecular modeling. Distance and dihedral angle restraints, generated from nuclear Overhauser effect spectroscopy (NOESY) and one-dimensional-NMR experiments respectively, were used to generate an ensemble of structures using distance geometry and simulated annealing. Molecular dynamic simulations produced three interconverting conformational families consistent with the NMR-derived constraints. We describe these conformations and their mechanism of interconversion. Furthermore, we have measured the IC50 s of a series of inhibitors generated from IP01 through alanine substitution of each residue. These results show that the L2-L3-R4-M5-R6 segment is functionally active, conformationally flexible, and contains a beta-turn involving residues R4-S7, while the C1-C9-I8-S7 segment is less functionally-active but adopts a more defined solution conformation, consistent with a scaffolding function. This model will be useful for designing nonpeptide-based organic inhibitors with improved pharmacologic properties.

Carrier Proteins↗

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core.

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.

Amino Acid Substitution↗

Local fluctuations and global unfolding of partially folded BPTI detected by NMR.

The protein [14-38]Abu is a chemically synthesized variant of bovine pancreatic trypsin inhibitor (BPTI) with the 14-38 disulfide bond intact and cysteines 5, 30, 51, and 55 replaced by alpha-amino-n-butyric acid (Abu). At 1-6 degrees C and pH 4.5-6.5, [14-38]Abu is partially folded with a native-like core [1]. Heteronuclear NMR spectra contain two, and in a few cases three or four, exchange cross peaks for each 15N-bound 1H, reporting the presence of two or more conformations that interconvert on a time scale of > or = milliseconds. Thermodynamic analysis of 15N-1H exchange peak volumes as a function of temperature in the range 1-35 degrees C indicates that partially folded [14-38]Abu undergoes local segmental motions as well as cooperative global unfolding. The relative abundance of more folded versus disordered conformations changes throughout the molecule, indicating that various regions of the partially folded protein are disordered to different extents prior to onset of thermal denaturation. This system is unique in providing a measure of the populations of interconverting partially folded conformations, as well as a microscopic view of cooperative folding of a fluctuating ensemble. Although global thermal denaturation is cooperative, significant deviation from simple two-state behavior is reflected in several parameters, including the difference in Tm for thermal unfolding measured by NMR versus circular dichroism.

Aprotinin↗

Conformation Design of Open-Chain Compounds.

At the brink of the 21st century, chemistry is increasingly concerned with the function that molecules fulfil as drugs, receptors, or-as ensemble of molecules-as materials. The capability of compounds to fulfil such functions cannot sufficiently be described by using only the terms composition and configuration. A decisive role is played in addition by the conformation of the molecules, which serves as the link between molecular composition and molecular function. Expressions such as "active conformation" or "competent conformation" allude to this aspect. Chemists have to develop an understanding how a flexible molecule adopts the conformation (a distinct shape) which is optimal for the function in question and how this process can be controlled. On the outset of such considerations, we may ask how nature succeeded in the process of evolution to endow flexible molecules with a preference to adopt the conformation which is optimal for the function it has to serve. In this review, I report on how we have reached a crude level of understanding of conformation design in nature with reference to the class of polyketide natural products, how we developed these insights into a conformation design of open-chain compounds, and which applications are already in sight.

Journal Article↗

Temperature dependence of the folding and unfolding kinetics of the GCN4 leucine zipper via 13C(alpha)-NMR.

Studies by one-dimensional NMR are reported on the interconversion of folded and unfolded forms of the GCN4 leucine zipper in neutral saline buffer. The peptide bears 99% 13C(alpha) labels at three sites: V9, L12, and G31. Time-domain 13C(alpha)-NMR spectra are interpreted by global Bayesian lineshape analysis to extract the rate constants for both unfolding and folding as functions of temperature in the range 47-71 degrees C. The data are well fit by the assumption that the same rate constants apply at each labeled site, confirming that only two conformational states need be considered. Results show that 1) both processes require a free energy of activation; 2) unfolding is kinetically enthalpy-opposed and entropy-driven, while folding is the opposite; and 3) the transition state dimer ensemble averages approximately 40% helical. The activation parameters for unfolding, derived from NMR data at the elevated temperatures where both conformations are populated, lead to estimates of the rate constant at low temperatures (5-15 degrees C) that agree with extant values determined by stopped-flow CD via dilution from denaturing media. However, the corresponding estimated values for the folding rate constant are larger by two to three orders of magnitude than those obtained by stopped flow. We propose that this apparent disagreement is caused by the necessity, in the stopped-flow experiment, for initiation of new helices as the highly denaturant-unfolded molecule adjusts to the newly created benign solvent conditions. This must reduce the success rate of collisions in producing the folded molecule. In the NMR determinations, however, the unfolded chains always have a small, but essential, helix content that makes such initiation unnecessary. Support for this hypothesis is adduced from recent extant experiments on the helix-coil transition in single-chain helical peptides and from demonstration that the folding rate constants for coiled coils, as obtained by stopped flow, are influenced by the nature of the denaturant used.

Amino Acid Sequence↗

A novel ensemble-based scoring and search algorithm for protein redesign and its application to modify the substrate specificity of the gramicidin synthetase a phenylalanine adenylation enzyme.

Realization of novel molecular function requires the ability to alter molecular complex formation. Enzymatic function can be altered by changing enzyme-substrate interactions via modification of an enzyme's active site. A redesigned enzyme may either perform a novel reaction on its native substrates or its native reaction on novel substrates. A number of computational approaches have been developed to address the combinatorial nature of the protein redesign problem. These approaches typically search for the global minimum energy conformation among an exponential number of protein conformations. We present a novel algorithm for protein redesign, which combines a statistical mechanics-derived ensemble-based approach to computing the binding constant with the speed and completeness of a branch-and-bound pruning algorithm. In addition, we developed an efficient deterministic approximation algorithm, capable of approximating our scoring function to arbitrary precision. In practice, the approximation algorithm decreases the execution time of the mutation search by a factor of ten. To test our method, we examined the Phe-specific adenylation domain of the nonribosomal peptide synthetase gramicidin synthetase A (GrsA-PheA). Ensemble scoring, using a rotameric approximation to the partition functions of the bound and unbound states for GrsA-PheA, is first used to predict binding of the wildtype protein and a previously described mutant (selective for leucine), and second, to switch the enzyme specificity toward leucine, using two novel active site sequences computationally predicted by searching through the space of possible active site mutations. The top scoring in silico mutants were created in the wetlab and dissociation/binding constants were determined by fluorescence quenching. These tested mutations exhibit the desired change in specificity from Phe to Leu. Our ensemble-based algorithm, which flexibly models both protein and ligand using rotamer-based partition functions, has application in enzyme redesign, the prediction of protein-ligand binding, and computer-aided drug design.

Adenosine Triphosphate↗

Dynamic modelling of a helical peptide in solution using NMR data: multiple conformations and multi-spin effects.

Nuclear Overhauser effect (NOE) measurements on molecules in solution provide information about only the ensemble-averaged properties of these molecules. An algorithm is presented that uses a list of NOEs to produce an ensemble of molecules that on average agrees with these NOEs, taking into account the effect of surrounding spins on the buildup of each NOE ('spin diffusion'). A simplified molecular dynamics simulation on several copies of the molecule in parallel is restrained by forces that are derived directly from differences between calculated and measured NOEs. The algorithm is tested on experimental NOE data of a helical peptide derived from bovine pancreatic trypsin inhibitor.

Algorithms↗

Revolutions in RNA secondary structure prediction.

RNA structure formation is hierarchical and, therefore, secondary structure, the sum of canonical base-pairs, can generally be predicted without knowledge of the three-dimensional structure. Secondary structure prediction algorithms evolved from predicting a single, lowest free energy structure to their current state where statistics can be determined from the thermodynamic ensemble. This article reviews the free energy minimization technique and the salient revolutions in the dynamic programming algorithm methods for secondary structure prediction. Emphasis is placed on highlighting the recently developed method, which statistically samples structures from the complete Boltzmann ensemble.

Algorithms↗

Two peptide fragments G55-I72 and K97-A109 from staphylococcal nuclease exhibit different behaviors in conformational preferences for helix formation.

Two synthetic peptides, SNasealpha1 and SNasealpha2, corresponding to residues G55-I72 and K97-A109, respectively, of staphylococcal nuclease (SNase), are adopted for detecting the role of helix alpha1 (E57-A69) and helix alpha2 (M98-Q106) in the initiation of folding of SNase. The helix-forming tendencies of the two SNase peptide fragments are investigated using circular dichroism (CD) and two-dimensional (2D) nuclear magnetic resonance (NMR) methods in water and 40% trifluoroethanol (TFE) solutions. The coil-helix conformational transitions of the two peptides in the TFE-H2O mixture are different from each other. SNasealpha1 adopts a low population of localized helical conformation in water, and shows a gradual transition to helical conformation with increasing concentrations of TFE. SNasealpha2 is essentially unstructured in water, but undergoes a cooperative transition to a predominantly helical conformation at high TFE concentrations. Using the NMR data obtained in the presence of 40% TFE, an ensemble of alpha-helical structures has been calculated for both peptides in the absence of tertiary interactions. Analysis of all the experimental data available indicates that formation of ordered alpha-helical structures in the segments E57-A69 and M98-Q106 of SNase may require nonlocal interactions through transient contact with hydrophobic residues in other parts of the protein to stabilize the helical conformations in the folding. The folding of helix alpha1 is supposed to be effective in initiating protein folding. The formation of helix alpha2 depends strongly on the hydrophobic environment created in the protein folding, and is more important in the stabilization of the tertiary conformation of SNase.

Amino Acid Sequence↗

Polyproline II helix is the preferred conformation for unfolded polyalanine in water.

Does aqueous solvent discriminate among peptide conformers? To address this question, we computed the solvation free energy of a blocked, 12-residue polyalanyl-peptide in explicit water and analyzed its solvent structure. The peptide was modeled in each of 4 conformers: alpha-helix, antiparallel beta-strand, parallel beta-strand, and polyproline II helix (P(II)). Monte Carlo simulations in the canonical ensemble were performed at 300 K using the CHARMM 22 forcefield with TIP3P water. The simulations indicate that the solvation free energy of P(II) is favored over that of other conformers for reasons that defy conventional explanation. Specifically, in these 4 conformers, an almost perfect correlation is found between a residue's solvent-accessible surface area and the volume of its first solvent shell, but neither quantity is correlated with the observed differences in solvation free energy. Instead, solvation free energy tracks with the interaction energy between the peptide and its first-shell water. An additional, previously unrecognized contribution involves the conformation-dependent perturbation of first-shell solvent organization. Unlike P(II), beta-strands induce formation of entropically disfavored peptide:water bridges that order vicinal water in a manner reminiscent of the hydrophobic effect. The use of explicit water allows us to capture and characterize these dynamic water bridges that form and dissolve during our simulations.

Computer Simulation↗

Structure refinement with molecular dynamics and a Boltzmann-weighted ensemble.

Time-averaging restraints in molecular dynamics simulations were introduced to account for the averaging implicit in spectroscopic data. Space- or molecule-averaging restraints have been used to overcome the fact that not all molecular conformations can be visited during the finite time of a simulation of a single molecule. In this work we address the issue of using the correct Boltzmann weighting for each member of an ensemble, both in time and in space. It is shown that the molecular- or space-averaging method is simple in theory, but requires a priori knowledge of the behaviour of a system. This is illustrated using a five-atom model system and the small cycle peptide analogue somatostatin. When different molecular conformers that are separated by energy barriers insurmountable on the time scale of a simulation contribute significantly to a measured NOE intensity, the use of space- or molecule-averaged distance restraints yields a more appropriate description of the measured data than conventional single-molecule refinement with or without application of time averaging.

Algorithms↗

Limits of NMR structure determination using variable target function calculations: ribonuclease T1, a case study.

Limits of NMR structure determination using multidimensional NMR spectroscopy, variable target function calculations and relaxation matrix analysis were explored using the model protein ribonuclease T1 (RNase T1). The enzyme consists of 104 amino acid residues and has a molecular mass of approximately 11 kDa. Primary experimental data comprise 1856 assigned NOE intensities, 493 3J coupling constants and 62 values of amid proton exchange rates. From these data, 2580 distance bounds, 168 allowed ranges for torsional angles and stereospecific assignments for 75% of beta-methylene protons as well as for 80% of diastereotopic methyl groups were derived. Whenever possible, the distance restraints were refined in a relaxation matrix analysis including amid proton exchange data for improvement of lower distance limits. Description of side-chain conformations were based on various models of motional averaging of 3J coupling constants. The final structure ensemble was selected from the starting ensemble comparing the global precision of structures with order parameters derived from 15N relaxation time measurements. Significant differences between the structure of RNase T1 in solution and in the crystal became apparent from a comparison of the two highly resolved structures.

Binding Sites↗

Toward absolute density of states calculations for proteins.

The density of states (DOS), which gives the number of conformations with a particular energy E, is a prerequisite in computing most thermodynamic quantities and in elucidating important biological processes such as the mechanism of protein folding. However, current methods for computing DOS of large systems such as proteins generally yield only the ratios of microstate counts for different energies, which could yield absolute conformation counts if the total number of conformations in phase space is known, thus motivating this work. Here, the total number of energy minima of 50-mer polyalanine, whose size corresponds to naturally occurring small proteins, was estimated under an all-atom potential energy function based on the cumulative distribution function (CDF) of conformational differences to be approximately 10(38). This estimate can place any DOS function, such as the Gaussian DOS distribution in the random energy model, on an absolute scale. Comparing the absolute conformational counts from a Gaussian DOS function with those from the CDF derived from quenched molecular dynamics ensembles shows that the former are far greater than the latter, indicating far fewer low-energy minima actually exist. In addition to showing how CDF and relative DOS calculations can yield absolute DOS for a discrete system, we also show how they can yield absolute DOS for continuous variable systems to a specified atomic variance. In the context of protein folding, knowing this phase-space "volume" of conformations in a DOS function, as well as characteristic transition times, constrains the set of possible folding mechanisms.

Protein Conformation↗