The expanded denatured state: an ensemble of conformations trapped in a locally encoded topological space.
Explore the source record for details and available documents.
Magnetic Resonance Spectroscopy↗
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.
Reproducing the conformations of protein-bound ligands: a critical evaluation of several popular conformational searching tools.
Several programs (Catalyst, Confort, Flo99, MacroModel, and Omega) that are commonly used to generate conformational ensembles have been tested for their ability to reproduce bioactive conformations. The ligands from thirty-two different ligand-protein complexes determined by high-resolution (< 2.0 A) X-ray crystallography have been analyzed. The Low-Mode Conformational Search method (with AMBER* and the GB/SA hydration model), as implemented in MacroModel, was found to perform better than the other algorithms. The rule-based method Omega, which is orders of magnitude faster than the other methods, also gave reasonable results but were found to be dependent on the input structure. The methods supporting diverse sampling (Catalyst, Confort) performed least well. For the seven ligands in the set having eight or more rotatable bonds, none of the bioactive conformations were ever found, save for one exception (Flo99). These ligands do not bind in a local minimum conformation according to AMBER*\GB/SA. Taking these last two observations together, it is clear that geometrically similar structures should be collected in order to increase the probability of finding the bioactive conformation among the generated ensembles. Factors influencing bioactive conformational retrieval have been identified and are discussed.
Algorithms↗
Conformational stability of a thrombin-binding peptide derived from the hirudin C-terminus.
The COOH-terminal region of hirudin represents an independent functional domain that binds to an anion-binding exosite of thrombin and inhibits the interaction of thrombin with fibrinogen and regulatory proteins in blood coagulation. The thrombin-bound structure of the peptide fragment, hirudin 55-65, has been determined by use of transferred NOE spectroscopy [Ni, F., Konishi, Y., & Scheraga, H. A. (1990) Biochemistry 29, 4479-4489]. The stability of the thrombin-bound conformation has been characterized further by a combined NMR and theoretical analysis of the conformational ensemble accessible by the hirudin peptide. Medium- and long-range NOE's were found for the free hirudin peptide in aqueous solution and in a mixture of dimethyl sulfoxide and water at both ambient (25 degrees C) and low (0 degrees C) temperatures, suggesting that ordered conformations are highly populated in solution. The global folding of these conformations is similar to that in the thrombin-bound state, as indicated by NOE's involving the side-chain protons of residues Phe(56), Ile(59), Pro(60), Tyr(63), and Leu(64). Residues Glu(61), Glu(62), Tyr(63), and Leu(64) all contain approximately 50% of helical conformations calculated from the ratio of the sequential dNN and d alpha N NOE's. Among the helical ensemble, active 3(10)-helical conformations were found by an analysis of the medium-range [(i,i+2) and (i,i+3)] NOE's involving the last six residues of the peptide. An analysis of the side-chain rotamers revealed that, upon binding to thrombin, there may be a rotation around the alpha CH-beta CH bond of Ile(59) such that Ile(59) adopts a gauche- (chi 1 = +60) conformation in contrast to the highly populated trans (chi 1 = -60) found for Ile(59) in the free peptide. However, the thrombin-bound conformation of the hirudin peptide is still an intrinsically stable conformer, and the preferred conformational ensemble of the peptide contains a large population of the active conformation. The apparent preference for a gauche- (chi 1 = +60) side-chain conformation of Ile(59) in the bound state may be explained by the existence of a positively charged arginine residue among the hydrophobic residues in the thrombin exosite.
Binding Sites↗
Circular dichroic investigation of the native and non-native conformational states of the growth factor receptor-binding protein 2 N-terminal src homology domain 3: effect of binding to a proline-rich peptide from guanine nucleotide exchange factor.
SH3 (src homology domain 3) domains are small protein modules that interact with proline-rich peptides. The structure of the N-terminal SH3 domain from growth factor receptor-binding protein 2 (Grb2), an adapter protein in the intracellular signaling pathway to Ras, was investigated by circular dichroic (CD) spectroscopy. The compact native beta-barrel conformation, previously elucidated by NMR spectroscopy, was largely predominant at pH = 4.8, in the absence of salt. From the structural changes induced by varying pH, ionic strength, temperature, or hydrophobicity of the environment, evidence for the existence of distinct nonnative conformations was obtained in the far- and near-UV domains. Along the free energy scale, these appear to distribute into two conformational ensembles, depending on the extent of structural and thermodynamic differences compared to the native conformation. The first ensemble consists of non-native conformations with a nativelike secondary structure, and the second is composed of partially unfolded conformations having short alpha-helical fragments or turnlike motifs in their nonnative secondary structure. Most of the observed nonnative conformations exist in mild or nondenaturing conditions. They probably have distinct compactness of their inner structure, depending on the strength of nonlocal interactions, but only the native all-beta conformation possesses a condensed protein exterior, appropriate for the binding to the VPPPVPPRRR decapeptide from Sos. Upon binding, the native conformation undergoes a local tertiary structure change in a hydrophobic pocket at the binding site. This is accompanied by the PP-II helix folding of the proline-rich peptide. Interestingly, in the near-UV domain, a significant change in the spectral contribution of an aromatic exciton was observed, thus allowing quantitative tracking of the binding process.
Adaptor Proteins, Signal Transducing↗
Fluctuations in ion pairs and their stabilities in proteins.
This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.
Amino Acid Sequence↗
Ligand-based prediction of active conformation by 3D-QSAR flexibility descriptors and their application in 3+3D-QSAR models.
A conceptionally new 3D molecular descriptor type and methodology are deduced by simple statistical thermodynamic reasoning, based on the free energy change encountered during a transformation of a conformational ensemble of the ligand to an active conformation. The performance of the descriptor was first tested on 37 endomorphin analogues with mu-opiate activity. The method resulted in predictive 3D-QSAR models, and the active conformation was also predicted. Generally, the methodology can be combined with the traditional 3D-QSAR techniques in a 3+3D-QSAR manner. This feature was tested on a series of 38 PGF2alphaprostaglandin analogues with antinidatory activity; the extent to which the molecular flexibility explains the variation in the biological activity was estimated and the active conformation was predicted. The novel descriptors in combination with the grid-based SOMFA descriptors resulted in 3+3D-QSAR models with good levels of predictivity leading to the approach of separation of the effect of the molecular interaction field of the active conformation and the effect of the conformational free energy loss.
Amino Acids↗
Efficient conformational sampling of local side-chain flexibility.
Side-chain flexibility of ligand-binding sites needs to be considered in the rational design of novel inhibitors. We have developed a method to generate conformational ensembles that efficiently sample local side-chain flexibility from a single crystal structure. The rotamer-based approach is tested here for the S1' pocket of human collagenase-1 (MMP-1), which is known to undergo conformational changes in multiple side-chains upon binding of certain inhibitors. First, a raw ensemble consisting of a large number of conformers of the S1' pocket was generated using an exhaustive search of rotamer combinations on a template crystal structure. A combination of principal component analysis and fuzzy clustering was then employed to successfully identify a core ensemble consisting of a low number of representatives from the raw ensemble. The core ensemble contained geometrically diverse conformers of stable nature, as indicated in several cases by a relative energy lower than that of the minimised template crystal structure. Through comparisons with X-ray crystallography and NMR structural data we show that the core ensemble occupied a conformational space similar to that observed under experimental conditions. The synthetic inhibitor RS-104966 is known to induce a conformational change in the side-chains of the S1' pocket of MMP-1 and could not be docked in the template crystal structure. However, the experimental binding mode was reproduced successfully using members of the core ensemble as the docking target, establishing the usefulness of the method in drug design.
Arginine↗
In silico prediction of blood-brain barrier permeation using the calculated molecular cross-sectional area as main parameter.
The cross-sectional area, AD, of a compound oriented in an amphiphilic gradient such as the air-water or lipid-water interface has previously been shown to be crucial for membrane partitioning and permeation, respectively. Here, we developed an algorithm that determines the molecular axis of amphiphilicity and the cross-sectional area, ADcalc, perpendicular to this axis. Starting from the conformational ensemble of each molecule, the three-dimensional conformation selected as the membrane-binding conformation was the one with the smallest cross-sectional area, ADcalcM, and the strongest amphiphilicity. The calculated, ADcalcM, and the measured, AD, cross-sectional areas correlated linearly (n=55, slope, m=1.04, determination coefficient, r2=0.95). The calculated cross-sectional areas, ADcalcM, were then used together with the calculated octanol-water distribution coefficients, log D7.4, of the 55 compounds (with a known ability to permeate the blood-brain barrier) to establish a calibration diagram for the prediction of blood-brain barrier permeation. It yielded a limiting cross-sectional area (ADcalcM=70 A2) and an optimal range of octanol-water distribution coefficients (-1.4<or=log D7.4<7.0). The calibration diagram was validated with an independent set of 43 compounds with the known ability to permeate the blood-brain barrier, yielding a prediction accuracy of 86%. The incorrectly predicted compounds exhibited log D7.4 values comprised between -0.6 and -1.4, suggesting that the limitation for log D7.4 is less rigorous than the limitation for AD. An accuracy of 83% has been obtained for a second validation set of 42 compounds which were previously shown to be difficult to predict. The calculated parameters, ADcalcM and log D7.4, thus allow for a fast and accurate prediction of blood-brain barrier permeation. Analogous calibration diagrams can be established for other membrane barriers.
Air↗
Folding funnels, binding funnels, and protein function.
Folding funnels have been the focus of considerable attention during the last few years. These have mostly been discussed in the general context of the theory of protein folding. Here we extend the utility of the concept of folding funnels, relating them to biological mechanisms and function. In particular, here we describe the shape of the funnels in light of protein synthesis and folding; flexibility, conformational diversity, and binding mechanisms; and the associated binding funnels, illustrating the multiple routes and the range of complexed conformers. Specifically, the walls of the folding funnels, their crevices, and bumps are related to the complexity of protein folding, and hence to sequential vs. nonsequential folding. Whereas the former is more frequently observed in eukaryotic proteins, where the rate of protein synthesis is slower, the latter is more frequent in prokaryotes, with faster translation rates. The bottoms of the funnels reflect the extent of the flexibility of the proteins. Rugged floors imply a range of conformational isomers, which may be close on the energy landscape. Rather than undergoing an induced fit binding mechanism, the conformational ensembles around the rugged bottoms argue that the conformers, which are most complementary to the ligand, will bind to it with the equilibrium shifting in their favor. Furthermore, depending on the extent of the ruggedness, or of the smoothness with only a few minima, we may infer nonspecific, broad range vs. specific binding. In particular, folding and binding are similar processes, with similar underlying principles. Hence, the shape of the folding funnel of the monomer enables making reasonable guesses regarding the shape of the corresponding binding funnel. Proteins having a broad range of binding, such as proteolytic enzymes or relatively nonspecific endonucleases, may be expected to have not only rugged floors in their folding funnels, but their binding funnels will also behave similarly, with a range of complexed conformations. Hence, knowledge of the shape of the folding funnels is biologically very useful. The converse also holds: If kinetic and thermodynamic data are available, hints regarding the role of the protein and its binding selectivity may be obtained. Thus, the utility of the concept of the funnel carries over to the origin of the protein and to its function.
Kinetics↗
Cytochrome c conformations resolved by the photon counting histogram: watching the alkaline transition with single-molecule sensitivity.
We apply the photon counting histogram (PCH) model, a fluorescence technique with single-molecule sensitivity, to study pH-induced conformational changes of cytochrome c. PCH is able to distinguish different protein conformations based on the brightness of a fluorophore sensitive to its local environment. We label cytochrome c through its single free cysteine with tetramethylrhodamine-5-maleimide (TMR), a fluorophore with specific brightnesses that we associate with specific protein conformations. Ensemble measurements demonstrate two different fluorescence responses with increasing pH: (i) a decrease in fluorescence intensity caused by the alkaline transition of cytochrome c (pH 7.0-9.5), and (ii) an increase in intensity when the protein unfolds (pH 9.5-10.8). The magnitudes of these two responses depend strongly on the molar ratio of TMR used to label cytochrome c. Using PCH we determine that this effect arises from the proportion of a nonfunctional conformation in the sample, which can be differentiated from the functional conformation. We further determine the causes of each ensemble fluorescence response: (i) during the alkaline transition, the fluorophore enters a dark state and discrete conformations are observed, and (ii) as cytochrome c unfolds, the fluorophore incrementally brightens, but discrete conformations are no longer resolved. Moreover, we also show that functional TMR-cytochrome c undergoes a response of identical magnitude regardless of the proportion of nonfunctional protein in the sample. As expected for a technique with single-molecule sensitivity, we demonstrate that PCH can directly observe the most relevant conformation, unlike ensemble fluorometry.
Alkalies↗
Conformation of desmopressin, an analogue of the peptide hormone vasopressin, in aqueous solution as determined by NMR spectroscopy.
Desmopressin (1-desamino-[DArg8]vasopressin, is a synthetic analogue of the neurohypophyseal peptide hormone vasopressin which has high antidiuretic and antibleeding potency. The structure of desmopressin has been determined in aqueous solution by two-dimensional NMR techniques and molecular dynamics simulations. Both standard and time-averaged distance restraints were used in structure calculations because of the inherent flexibility in small peptides. 21 models calculated with standard restraints were compared with structures refined with time-averaged distance restraints and were found to be good representatives of the conformational ensemble of desmopressin. The macrocyclic ring forms an inverse gamma-turn centered around Gln4. Residues 1 and 2, the disulphide bridge and the three-residue acyclic tail were found to be flexible in solution. Residues 4-6 in the ensemble of calculated structures contain essentially the same backbone conformation as in the crystal structure of pressinoic acid, the cyclic moiety of vasopressin, whereas residues 2-6 superimpose on the NMR-derived conformation of oxytocin bound to neurophysin. The results presented in this work suggest that, in addition to the differences in sequence between desmopressin and vasopressin, differences in conformational and dynamic properties between the two compounds explain their pharmacological differences.
Amides↗
Probing the folding capacity and residual structures in 1-79 residues fragment of staphylococcal nuclease by biophysical and NMR methods.
1-79 residues SNase fragment (SNase79) has chain length containing a sequence for helix alpha(1), omega-loop, beta(I)-sheet, and partial beta(II)-sheet of native SNase. The incomplete "beta-barrel" structural region of SNase79 makes this fragment to be interested in investigation of its conformation. For this study, we use CD, fluorescence, and NMR spectroscopy to probe the folding capacity and the residual structures in SNase79. The optical spectra obtained for SNase79 and its mutants reveal the presence of retained capacity for folding of the fragment. The NMR derived (13)C(alpha) secondary chemical shifts, (3)J(NH-Halpha) coupling constants, amide-proton temperature coefficients, interresidue NOEs, and (15)N relaxation data determine the intrinsic propensities for helix- and turn- or beta-sheet-like conformations of SNase79, which is not the result of stabilizing inter-molecular interactions by oligomerization effects. The residual turn- and helix-like structures may serve as potential local nucleation sites, whereas the residual beta(I)-sheet-like structure can be regarded as a potential non-local nucleation site in the folding of SNase79. The intrinsic local and non-local interactions in these potential initiation sites are insufficient to stabilize the folding of SNase79 due to the shortage of relevant long-range interactions from other part of the fragment. The conformational ensemble of SNase79 is a highly heterogeneous collection of interconverting conformations having transiently populated helix- and beta-sheet- or turn-like structures.
Amino Acid Sequence↗
Dichroic statistical model for prediction and analysis of peptide helicity.
Traditional statistical models for the prediction of peptide helicity are written in terms of the mean fractional helicity of the peptide residues. Far ultraviolet circular dichroic measurements of peptide solutions are converted to mean fractional helicity by partitioning the observed ellipticity between that of a perfect helix and a random coil. This partition does not adequately represent the ensemble of peptide molecules present in solution that populate imperfect helical conformations of quite variable lengths. A new dichroic statistical model has been written in terms of ellipticity rather than fractional helical content that recognizes (1) the source of ellipticity, peptide bond adsorption; (2) the differential ellipticity of peptide bonds in the terminal and interior helical turns; and (3) the contributions of each participant in a conformational ensemble to the observed ellipticity. Comparative analyses of host/guest peptides indicates that significant differences are obtained between residue w and n weights and ellipticity values using the traditional and dichroic statistical models.
Algorithms↗
Searching for folding initiation sites of staphylococcal nuclease: a study of N-terminal short fragments.
The N-terminal short fragments of staphylococcal nuclease (SNase), SNase20, SNase28, and SNase36, corresponding to the sequence regions, Ala1-Gly20, Ala1-Lys28, and Ala1-Leu36, respectively, as well as an 8-residue peptide (Ala17-Ile18-Asp19-Gly20-Asp21-Thr22-Val23-Lys24) have been synthesized. The conformational states of these fragments were investigated using CD and NMR spectroscopy in aqueous solution and in trifluoroethanol (TFE)-H(2)O mixture. SNase20 containing a sequence corresponding to a bent peptide in native SNase shows a transient population of bend-like conformation around Ala12-Thr13-Leu14 in TFE-H(2)O mixture. The sequence region of Ala17-Thr22 of SNase28 displays a localized propensity for turn-like conformation in both aqueous solution and TFE-H(2)O mixture. The conformational ensemble of SNase36 in aqueous solution includes populated turn-like conformations localized in sequence regions Ala17-Thr22 and Tyr27-Gln30. The analysis suggests that these sequence regions, which form the regular secondary structures in native protein, may serve as the folding nucleation sites of SNase fragments of different chain lengths starting from the N-terminal end. Thus, the formation of bend- and turn-like conformations of these sequence regions may be involved in the early folding events of the SNase polypeptide chain in vitro.
Amino Acid Sequence↗
Conformational equilibrium of cytochrome P450 BM-3 complexed with N-palmitoylglycine: a replica exchange molecular dynamics study.
UV-vis absorbance measurements and associated studies of cytochrome P450 BM-3 in complex with N-palmitoylglycine (NPG) indicate that a conformational change occurs in the active site of the complex where the terminal atoms of the ligand move from a site distant from the heme iron, as seen in the low temperature crystal structure to a site proximal to the heme iron at biological temperatures. We employ replica exchange molecular dynamics simulations to study this conformational change. The population of the proximal state is found to increase with temperature in agreement with UV-vis absorbance and NMR measurements. In addition to the conformations characterized by X-ray crystallography and computer modeling, this study shows that a new conformational state is significantly populated at room temperature. The observed increase in the population of conformations where the terminal atoms of NPG are proximal to the heme iron with increasing temperature indicates that the proximal state is stabilized by conformational entropy. A proposal for the origin of this entropic stabilization is provided on the basis of the structure of the newly identified state. We use the temperature weighted histogram (T-WHAM) method to characterize the transition state regions of the conformational ensemble and propose a mechanism of interconversion between these low free energy conformational states.
Crystallography, X-Ray↗
A search for sources of drug resistance by the 4D-QSAR analysis of a set of antimalarial dihydrofolate reductase inhibitors.
A set of 18 structurally diverse antifolates including pyrimethamine, cycloguanil, methotrexate, aminopterin and trimethoprim, and 13 pyrrolo[2,3-d pyrimidines were studied using four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis. The corresponding biological activities of these compounds include IC50 inhibition constants for both the wild type, and a specific mutant type of Plasmodium falciparum dihydrofolate reductase (DHFR). Two thousand conformations of each analog were sampled to generate a conformational ensemble profile (CEP) from a molecular dynamics simulation (MDS) of 100,000 conformer trajectory states. Each sampled conformation was placed in a 1 A cubic grid cell lattice for each of five trial alignments. The frequency of occupation of each grid cell was computed for each of six types of pharmacophore groups of atoms of each compound. These grid cell occupancy descriptors (GCODs) were then used as a descriptor pool to construct 4D-QSAR models. Models for inhibition of both the 'wild' type and the mutant enzyme were generated which provide detailed spatial pharmacophore requirements for inhibition in terms of atom types and their corresponding relative locations in space. The 4D-QSAR models indicate some structural features perhaps relevant to the mechanism of resistance of the Plasmodium falciparum DHFR to current antimalarials. One feature identified is a slightly different binding alignment of the ligands to the mutant form of the enzyme as compared to the wild type.
Animals↗
Assembly of polypeptide and protein backbone conformations from low energy ensembles of short fragments: development of strategies and construction of models for myoglobin, lysozyme, and thymosin beta 4.
Recently we developed methods for the construction of knowledge-based mean fields from a data base of known protein structures. As shown previously, this approach can be used to calculate ensembles of probable conformations for short fragments of polypeptide chains. Here we develop procedures for the assembly of short fragments to complete three-dimensional models of polypeptide chains. The amino acid sequence of a given protein is decomposed into all possible overlapping fragments of a given length, and an ensemble of probable conformations is calculated for each fragment. The fragments are assembled to a complete model by choosing appropriate conformations from the individual ensembles and by averaging over equivalent angles. Finally a consistent model is obtained by rebuilding the conformation from the average angles. From the average angles the local variability of the structure can be calculated, which is a useful criterion for the reliability of the model. The procedure is applied to the calculation of the local backbone conformations of myoglobin and lysozyme whose structures have been solved by X-ray analysis and thymosin beta 4, a polypeptide of 43 amino acid residues whose structure was recently investigated by NMR spectroscopy. We demonstrate that substantial fractions of the calculated local backbone conformations are similar to the experimentally determined structures.
Amino Acid Sequence↗
Inhibiting viral proteases: challenges and opportunities.
Inhibitor design against viral targets must take into account the peculiar characteristics of viral biology-in particular, the plasticity of their replicative machinery. This includes maturational cleavage of the polyprotein, which is mediated by virally encoded proteases. Designing against a movable target is particularly challenging, but at the same time it offers new opportunities. Here we describe our experience with the NS3/4A (NS: nonstructural) serine protease of human hepatitis C virus (HCV). By extensive use of combinatorial peptide libraries, various inhibitor types were generated, including product inhibitors, serine traps, P-P' inhibitors, and prime side inhibitors. The latter represent a first case for a serine protease. A key finding, derived from structural studies utilizing these inhibitors, was that NS3 is an induced-fit protease, requiring both the NS4A cofactor protein and the substrate to fully activate its catalytic machinery. In the absence of cofactor and/or substrate, NS3 exists in solution as a large conformational ensemble, which can be matched by a correspondingly large set of peptide inhibitors, each one stabilizing a given conformer. In the perspective of inhibiting viral proteases in general, we suggest that combinatorial ligand ensembles may be a powerful tool, to contrast the adaptive potential of the viral quasispecies.
Amino Acid Sequence↗