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 37 records · Page 2Linked to original sources

Binding sites in Escherichia coli dihydrofolate reductase communicate by modulating the conformational ensemble.

To explore how distal mutations affect binding sites and how binding sites in proteins communicate, an ensemble-based model of the native state was used to define the energetic connectivities between the different structural elements of Escherichia coli dihydrofolate reductase. Analysis of this model protein has allowed us to identify two important aspects of intramolecular communication. First, within a protein, pair-wise couplings exist that define the magnitude and extent to which mutational effects propagate from the point of origin. These pair-wise couplings can be identified from a quantity we define as the residue-specific connectivity. Second, in addition to the pair-wise energetic coupling between residues, there exists functional connectivity, which identifies energetic coupling between entire functional elements (i.e., binding sites) and the rest of the protein. Analysis of the energetic couplings provides access to the thermodynamic domain structure in dihydrofolate reductase as well as the susceptibility of the different regions of the protein to both small-scale (e.g., point mutations) and large-scale perturbations (e. g., binding ligand). The results point toward a view of allosterism and signal transduction wherein perturbations do not necessarily propagate through structure via a series of conformational distortions that extend from one active site to another. Instead, the observed behavior is a manifestation of the distribution of states in the ensemble and how the distribution is affected by the perturbation.

Binding Sites↗

Ultrarapid mixing experiments shed new light on the characteristics of the initial conformational ensemble during the folding of ribonuclease A.

The earliest folding events in single-tryptophan mutants of RNase A were investigated by fluorescence measurements by using a combination of stopped-flow and continuous-flow mixing experiments covering the time range from 70 micros to 10 s. An ultrarapid double-jump mixing protocol was used to study refolding from an unfolded ensemble containing only native proline isomers. The continuous-flow measurements revealed a series of kinetic events on the submillisecond time scale that account for the burst-phase signal observed in previous stopped-flow experiments. An initial increase in fluorescence within the 70-micros dead time of the continuous-flow experiment is consistent with a relatively nonspecific collapse of the polypeptide chain whereas a subsequent decrease in fluorescence with a time constant of approximately 80 micros is indicative of a more specific structural event. These rapid conformational changes are not observed if RNase A is allowed to equilibrate under denaturing conditions, resulting in formation of nonnative proline isomers. Thus, contrary to previous expectations, the isomerization state of proline peptide bonds can have a major impact on the structural events during early stages of folding.

Kinetics↗

NMR characterizations of an amyloidogenic conformational ensemble of the PI3K SH3 domain.

Amyloid formation is associated with structural changes of native polypeptides to monomeric intermediate states and their self-assembly into insoluble aggregates. Characterizations of the amyloidogenic intermediate state are, therefore, of great importance in understanding the early stage of amyloidogenesis. Here, we present NMR investigations of the structural and dynamic properties of the acid-unfolded amyloidogenic intermediate state of the phosphatidylinositol 3-kinase (PI3K) SH3 domain--a model peptide. The monomeric amyloidogenic state of the SH3 domain studied at pH 2.0 (35 degrees C) was shown to be substantially disordered with no secondary structural preferences. (15)N NMR relaxation experiments indicated that the unfolded polypeptide is highly flexible on a subnanosecond timescale when observed under the amyloidogenic condition (pH 2.0, 35 degrees C). However, more restricted motions were detected in residues located primarily in the beta-strands as well as in a loop in the native fold. In addition, nonnative long-range interactions were observed between the residues with the reduced flexibility by paramagnetic relaxation enhancement (PRE) experiments. These indicate that the acid-unfolded state of the SH3 domain adopts a partly folded conformation through nonnative long-range contacts between the dynamically restricted residues at the amyloid-forming condition.

Amino Acid Sequence↗

Targacept active conformation search: a new method for predicting the conformation of a ligand bound to its protein target.

Targacept active conformation search (TACS) is a novel variation of well-established three-dimensional quantitative structure--activity relationship methodologies that seeks to determine probable conformation(s) of ligands bound to their protein targets. A combination of affinity or activity data and energetically accessible conformational ensembles, each conformer described by three-dimensional (3-D) sensitive descriptors, forms the basis of the TACS data model. Recursive pruning is used to reduce the size of both the conformational ensemble and the descriptor space until the TACS data model contains just enough information to determine probable conformation(s) of ligands bound to their protein targets. The TACS algorithm is comprised of five components: (1) conformational ensemble generation, (2) 3-D sensitive descriptor calculation, (3) ensemble descriptor preprocessing, (4) model generation, and (5) prediction of bound conformation(s). Significantly, this method precludes the need for subjective or objective molecular alignment. We report the application of this technique to five benchmark protein-ligand couples where the conformation of a bound ligand has been previously established using X-ray crystallography: 9-cis-retinoic (1) and 9-trans-retinoic acid (2), both agonists for the retinoic acid receptor gamma, compounds KH1060 (3) and MC1288 (4), which bind to the vitamin D3 receptor, and R04 (5), an inhibitor bound to human rhinovirus 14 thermolysin. The binding conformations predicted by TACS were compared to the crystallographic structures extracted from their respective binding sites using root-mean-squared deviation (rmsd) criteria. Three of the conformations found using TACS were within crystallographic error. 9-cis-Retinoic acid, 9-trans-retinoic acid, and MC1288, when superimposed on their crystallographic structures, gave rmsd values of 0.22, 0.17, and 0.34 A, respectively. The rmsd values for KH1060 (1.54 A) and R04 (1.01 A) were larger but still reasonable.

Binding Sites↗

Quantifying polypeptide conformational space: sensitivity to conformation and ensemble definition.

Quantifying the density of conformations over phase space (the conformational distribution) is needed to model important macromolecular processes such as protein folding. In this work, we quantify the conformational distribution for a simple polypeptide (N-mer polyalanine) using the cumulative distribution function (CDF), which gives the probability that two randomly selected conformations are separated by less than a "conformational" distance and whose inverse gives conformation counts as a function of conformational radius. An important finding is that the conformation counts obtained by the CDF inverse depend critically on the assignment of a conformation's distance span and the ensemble (e.g., unfolded state model): varying ensemble and conformation definition (1 --> 2 A) varies the CDF-based conformation counts for Ala(50) from 10(11) to 10(69). In particular, relatively short molecular dynamics (MD) relaxation of Ala(50)'s random-walk ensemble reduces the number of conformers from 10(55) to 10(14) (using a 1 A root-mean-square-deviation radius conformation definition) pointing to potential disconnections in comparing the results from simplified models of unfolded proteins with those from all-atom MD simulations. Explicit waters are found to roughen the landscape considerably. Under some common conformation definitions, the results herein provide (i) an upper limit to the number of accessible conformations that compose unfolded states of proteins, (ii) the optimal clustering radius/conformation radius for counting conformations for a given energy and solvent model, (iii) a means of comparing various studies, and (iv) an assessment of the applicability of random search in protein folding.

Computer Simulation↗

COREX/BEST server: a web browser-based program that calculates regional stability variations within protein structures.

SUMMARY: Utilizing the user-supplied coordinates of a protein structure, the COREX/BEST Server generates a structural thermodynamic ensemble. This conformational ensemble can then be used to calculate the regional variations in stability of a protein structure, and the stabilities are presented in units of energy (kcal/mol). The regional stabilities, which are calculated at the resolution of individual residues, can be mapped onto the protein structure for visual representation and downloaded from the site in the form of tab delimited text. The site provides an easy to follow summary of the theoretical and algorithmic approaches and provides links to references for more detailed descriptions. AVAILABILITY: The COREX/BEST Server may be accessed through a typical web browser by visiting http://best.utmb.edu/BEST/.

Algorithms↗

Descriptors you can count on? Normalized and filtered pharmacophore descriptors for virtual screening.

The three-dimensional (3D) binary pharmacophore fingerprints find wide application as descriptors in applications ranging from virtual screening through library design. While the 3D content they capture is an intuitively attractive feature of such measures, maximizing their signal to noise ratio has proven to be a tricky balancing act. This issue surfaces primarily due to the potential of such fingerprints to create an explosion of pharmacophores as molecular complexity and flexibility increases. In this article, we describe a modification to the fingerprint generation process that normalizes pharmacophore occurrence frequency by the conformational ensemble size used to derive the descriptor. By including pharmacophore frequency and conformational count, the importance of a given pharmacophore is weighted by the probability of its existence within a given conformational ensemble, rather than treating each pharmacophore equally. In addition, a number of filters have been added to permit the removal of unwanted pharmacophores from the descriptor set. These filters are based on pharmacophore composition (e.g. permutations made up primarily of lipophilic and/or aromatic centers), and size (pharmacophore perimeter length relative to the largest perimeter length found in the molecule). The highly uneven nature of pharmacophore distributions across the conformational ensemble used to generate them is highlighted, as are enrichment comparisons with their binary fingerprint peers. In addition, the limitations in descriptor comparison validation are high-lighted as an illustration of the need for more extensive validation experiments.

Drug Design↗

Applications of time-resolved resonance energy transfer measurements in studies of the molecular crowding effect.

The native structures of many globular proteins are only weakly stabilized and form in solution ensembles of multiple conformers. The energy differences between the conformers are assumed to be small. This is the case of flexible multidomain proteins where domain motions were observed. High concentrations of inert macrosolute, which create a crowded or confined environment, can cause shifts of the distribution of the conformers of such proteins towards the more compact structures. This effect may also promote compact structures in partially folded proteins. Time-resolved dynamic non-radiative excitation energy transfer (tr-RET) is suitable for detection of either subtle or major changes in distributions of intramolecular distances in protein molecules in solutions. Two experiments were performed which demonstrated the applicability of tr-RET for detection of the effect of macrosolutes on the conformational ensembles of flexible states of protein molecules. The distribution of distances between residues 203 and 169 in the CORE domain of E. coli adenylate kinase (AK) in the denatured state was determined in the presence of high concentrations of dextran 40. A significant shift of the mean of the distribution was observed without reduction of its width. This was interpreted as a shift to compact structure without change of the degree of disorder of the chain. In a second experiment the distribution of the distance between residues 55 and 169 in AK, which spans the cleft between the CORE and the AMPbind domains, was monitored. No clear effect of high concentrations of dextran 40 was found. These experiments show the strength of the application of tr-RET in investigation of changes in the sub-states of flexible conformations of globular protein. Networks of pairs of labeled sites can be prepared and tr-RET experiments can be performed in order to search for the segments of the protein molecules, which respond to the presence of inert macromolecules in their environment.

Adenylate Kinase↗

Impact of conformational flexibility on three-dimensional similarity searching using correlation vectors.

Many three-dimensional (3D) virtual screening concepts, like automated docking or pharmacophore searching, rely on the calculation of a "bioactive" or "receptor-relevant" conformation of a molecule to assess its biological activity. We investigated the dependence of the presence of conformations near the "bioactive" conformation on three-dimensional similarity searching with pharmacophore-based correlation vectors (CATS3D approach). Cocrystal structures of 11 target classes served as queries for virtual screening of a database of annotated ligands. Different numbers of conformations were calculated. Single 3D structures were obtained using the 3D structure generator CORINA and conformational ensembles by the conformation generation program ROTATE. This approach was able to reproduce conformations for high resolution cocrystal structures. For virtual screening we found that using only the CORINA-generated single conformation already resulted in a significant enrichment of isofunctional molecules having the same biological property profile. This observation was also made for ligand classes with many rotatable bonds. Although more similar conformations were considered to be more similar in the CATS3D description, the impact of using multiple conformations on the enrichment of actives was not as high as expected. CATS3D provides an alignment-free three-dimensional virtual screening approach that is less dependent on the presence of conformations which are close to the "bioactive" conformation of a molecule compared to methods that rely on an explicit three-dimensional alignment of molecules.

Algorithms↗

Solution conformation of a peptide fragment representing a proposed RNA-binding site of a viral coat protein studied by two-dimensional NMR.

The first 25 amino acids of the coat protein of cowpea chlorotic mottle virus are essential for binding the encapsidated RNA. Although an alpha-helical conformation has been predicted for this highly positively charged N-terminal region [Argos, P. (1981) Virology 110, 55-62; Vriend, G., Verduin, B. J. M., & Hemminga, M. A. (1986) J. Mol. Biol. 191, 453-460], no experimental evidence for this conformation has been presented so far. In this study, two-dimensional proton NMR experiments were performed on a chemically synthesized pentacosapeptide containing the first 25 amino acids of this coat protein [Ten Kortenaar, P. B. W., Krüse, J., Hemminga, M. A., & Tesser, G. I. (1986) Int. J. Pept. Protein Res. 27, 401-413]. All resonances could be assigned by a combined use of two-dimensional correlated spectroscopy and nuclear Overhauser enhancement spectroscopy carried out at four different temperatures. Various NMR parameters indicate the presence of a conformational ensemble consisting of helical structures rapidly converting into more extended states. Differences in chemical shifts and nuclear Overhauser effects indicate that lowering the temperature induces a shift of the dynamic equilibrium toward more helical structures. At 10 degrees C, a perceptible fraction of the conformational ensemble consists of structures with an alpha-helical conformation between residues 9 and 17, likely starting with a turnlike structure around Thr9 and Arg10. Both the conformation and the position of this helical region agree well with the secondary structure predictions mentioned above.

Amino Acid Sequence↗

Modelling of third cytoplasmic loop of bovine rhodopsin by multicanonical molecular dynamics.

The third cytoplasmic loop (C3) of bovine rhodopsin (Rh) is an important site for its interaction with G-protein transducin. The tertiary structure of Rh was determined by X-ray crystallography, although the local conformation around the C3 loop (residues: 236-240) was not visible in electron density maps. We constructed a canonical conformation ensemble at 310 K for the C3 loop (residues: 227-244) using a multicanonical molecular dynamics simulation, and predicted several putative conformations. The conformation ensemble was classified by principal component analysis into several distinct structural clusters, some of which could provide the putative structural models of Rh and the activated state of Rh.

Amino Acid Sequence↗

Structural characterization of modular supramolecular architectures in solution.

Structures of modular supramolecular architectures consisting of a hexameric, diphenylethyne-linked porphyrin macrocyclic array and the corresponding host-guest complex formed by inclusion of a tripyridyl guest molecule were characterized in solution using high-angle X-ray scattering. Scattering measurements made to 6 A resolution coupled with pair distance function (PDF) analyses demonstrated that (1) the porphyrin architectures are not rigid but are distributed across a conformational ensemble with a mean diameter that is 1.5 A shorter than the diameter of a symmetric, energy-minimized model structure, (2) the conformational envelope has limits of 3 A positional dispersion and full rotational freedom for all six porphyrin groups, and (3) insertion of the tripyridyl guest molecule expands the diameter of the host conformer by 0.6 A and decreases the configurational dispersion by approximately 2-fold. These results validate the molecular design, provide a new measure of conformational ensembles in solution that cannot be obtained by other techniques, and establish a structural basis for understanding the photophysical and guest-hosting functions of the hexameric porphyrin architectures in liquids.

Journal Article↗

Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques.

Pulse field gradient NMR methods have been used to determine the effective hydrodynamic radii of a range of native and nonnative protein conformations. From these experimental data, empirical relationships between the measured hydrodynamic radius (R(h)) and the number of residues in the polypeptide chain (N) have been established; for native folded proteins R(h) = 4.75N (0.29)A and for highly denatured states R(h) = 2.21N (0.57)A. Predictions from these equations agree well with experimental data from dynamic light scattering and small-angle X-ray or neutron scattering studies reported in the literature for proteins ranging in size from 58 to 760 amino acid residues. The predicted values of the hydrodynamic radii provide a framework that can be used to analyze the conformational properties of a range of nonnative states of proteins. Several examples are given here to illustrate this approach including data for partially structured molten globule states and for proteins that are unfolded but biologically active under physiological conditions. These reveal evidence for significant coupling between local and global features of the conformational ensembles adopted in such states. In particular, the effective dimensions of the polypeptide chain are found to depend significantly on the level of persistence of regions of secondary structure or features such as hydrophobic clusters within a conformational ensemble.

Adhesins, Bacterial↗

Effects of C alpha-methyl substitution on the conformation of linear GnRH antagonist analogs.

We have examined the effect of C alpha-methyl groups on the conformational ensemble of GnRH analog peptides by comparing 1H 2D NMR data from two analogs, Ac-D-Nal1-D-4-Cl-C alpha-Me-Phe2-D-Pal3-Ser4-Tyr5-D-Arg6-Leu7-Arg8-Pro9-D-Al a10- NH2 (1) and Ac-D-Nal1-D-4-Cl-C alpha-Me-Phe2-D-Pal3-Ser4-C alpha-Me-Tyr5-D-Arg6- Leu7-C alpha-Me-Arg8-Pro9-D-Ala10-NH2 (2). The two additional C alpha-methyl groups in residues 5 and 8 of 2 do not influence significantly the pattern of the observable main chain NOE intensities, or of the backbone HN proton chemical shifts, which indicates that they do not produce global changes in the conformational ensemble of the peptide. A local change induced by the substitution was observed in the conformation at D-Arg8-Pro9.

Amino Acid Sequence↗

Lattice protein folding with two and four-body statistical potentials.

The cooperative folding of proteins implies a description by multibody potentials. Such multibody potentials can be generalized from common two-body statistical potentials through a relation to probability distributions of residue clusters via the Boltzmann condition. In this exploratory study, we compare a four-body statistical potential, defined by the Delaunay tessellation of protein structures, to the Miyazawa-Jernigan (MJ) potential for protein structure prediction, using a lattice chain growth algorithm. We use the four-body potential as a discriminatory function for conformational ensembles generated with the MJ potential and examine performance on a set of 22 proteins of 30-76 residues in length. We find that the four-body potential yields comparable results to the two-body MJ potential, namely, an average coordinate root-mean-square deviation (cRMSD) value of 8 A for the lowest energy configurations of all-alpha proteins, and somewhat poorer cRMSD values for other protein classes. For both two and four-body potentials, superpositions of some predicted and native structures show a rough overall agreement. Formulating the four-body potential using larger data sets and direct, but costly, generation of conformational ensembles with multibody potentials may offer further improvements. Proteins 2001;43:161-174.

Algorithms↗

Considerations in compound database preparation--"hidden" impact on virtual screening results.

Structure-based virtual screening (SBVS) utilizing docking algorithms has become an essential tool in the drug discovery process, and significant progress has been made in successfully applying the technique to a wide range of receptor targets. In silico validation of virtual screening protocols before application to a receptor target using a corporate or commercially available compound collection is key to establishing a successful process. Ultimately, retrieval of a set of active compounds from a database of inactives is required, and the metric of enrichment (E) is habitually used to discern the quality of separation of the two. Numerous reports have addressed the performance of docking algorithms with regard to the quality of binding mode prediction and the issue of postprocessing "hit lists" of docked ligands. However, the impact of ligand database preprocessing has yet to be examined in the context of virtual screening and prioritization of compounds for biological evaluation. We provide an insight into the implications of cheminformatic preprocessing of a validation database of compounds where multiple protonated, tautomeric, stereochemical, and conformational states have been enumerated. Several commonly used methods for the generation of ligand conformations and conformational ensembles are examined, paired with an exhaustive rigid-body algorithm for the docking of different "multimeric" compound representations to the ligand binding site of the human estrogen receptor alpha. Chemgauss, a shapegaussian scoring function with intrinsic chemical knowledge, was combined with PLP as a consensus-scoring scheme to rank output from the docking protocol and enrichment rates calculated for each screen. The overheads of CPU consumption and the effect on relative database size (disk requirement) for each of the protocols employed are considered. Assessment of these parameters indicates that SBVS enrichments are highly dependent on the initial cheminformatic treatment(s) used in database construction. The interplay of SMILES representations, stereochemical information, protonation state enumeration, and ligand conformation ensembles are critical in achieving optimum enrichment rates in such screening.

Computer Simulation↗