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

PubMed2006

On using oscillating time-dependent restraints in MD simulation.

The use of time-dependent restraints in molecular simulation in order to generate a conformational ensemble for molecules that is in accordance with measured ensemble averages for particular observable quantities is investigated. Using a model system consisting of liquid butane and the cyclic peptide antamanide the reproduction of particular average (3) J-coupling constant values in a molecular dynamics simulation is analysed. It is shown that the multiple-valuedness and the sizeable gradients of the Karplus curve relating (3) J-coupling constants measured in NMR experiments to the corresponding torsional-angle values cause severe problems when trying to restrain a (3) J-coupling constant to a value close to the extrema of the Karplus curve. The introduction of a factor oscillating with time into the restraining penalty function alleviates this problem and enhances the restrained conformational sampling.

Algorithms↗

PubMed2004

Analysis of the "thermodynamic information content" of a Homo sapiens structural database reveals hierarchical thermodynamic organization.

Classification of the amounts and types of lower order structural elements in proteins is a prerequisite to effective comparisons between protein folds. In an effort to provide an additional vehicle for fold comparison, we present an alternative classification scheme whereby protein folds are represented in statistical thermodynamic terms in such a way as to illuminate the energetic building blocks within protein structures. The thermodynamic relationship is examined between amino acid sequences and the conformational ensembles for a database of 159 Homo sapiens protein structures ranging from 50 to 250 amino acids. Using hierarchical clustering, it is shown through fold-recognition experiments that (1) eight thermodynamic environmental descriptors sufficiently accounts for the energetic variation within the native state ensembles of the H. sapiens structural database, (2) an amino acid library of only six residue types is sufficient to encode >90% of the thermodynamic information required for fold specificity in the entire database, and (3) structural resolution of the statistically derived environments reveals sequential cooperative segments throughout the protein, which are independent of secondary structure. As the first level of thermodynamic organization in proteins, these segments represent the thermodynamic counterpart to secondary structure.

Computer Simulation↗

PubMed2006

The changing landscape of protein allostery.

It is becoming increasingly clear that the fundamental capacity to undergo conformational change in response to ligand binding is intrinsic to proteins. This property confers on proteins the ability to be allosterically modulated in order to shift substrate binding affinities, alter enzymatic activity or regulate protein-protein interaction. How this allosteric modulation occurs--the pathways of communication, the shifting of conformational ensembles and the altered molecular dynamics--has received considerable attention during the past two years. Recent progress has helped outline the molecular origins of allostery in proteins as diverse as Hsp70 molecular chaperones and signal integrating proteins, such as WASP. In addition, allosteric properties have been successfully engineered into proteins for drug design or the development of novel biosensors. Methodological advances have provided exciting prospects for new insights and new biological roles of allosteric systems have been uncovered.

Allosteric Regulation↗

PubMed2006

Docking of aminoglycosides to hydrated and flexible RNA.

Although much effort has been devoted to the development of programs suited for the docking of ligands to proteins, much less progress has been achieved in the nucleic acid field. We have developed a unique approach for docking aminoglycosides to RNA considering the flexibility of these macromolecules using conformational ensembles and accounting for the role of the first hydration shell. This concept, successfully implemented in AutoDock, relies on the computation of the intermolecular interaction energy that accounts for the presence of dynamically bound water molecules to the RNA. As an application, a set of 11 aminoglycosides was docked with an average root-mean-square deviation (RMSD) of 1.41 A to be compared with an average RMSD of 3.25 A when the original AutoDock protocol was used.

Aminoglycosides↗

PubMed2004

Identification of an alternative ligand-binding pocket in the nuclear vitamin D receptor and its functional importance in 1alpha,25(OH)2-vitamin D3 signaling.

Structural and molecular studies have shown that the vitamin D receptor (VDR) mediates 1alpha,25(OH)2-vitamin D3 gene transactivation. Recent evidence indicates that both VDR and the estrogen receptor are localized to plasma membrane caveolae and are required for initiation of nongenomic (NG) responses. Computer docking of the NG-specific 1alpha,25(OH)2-lumisterol to the VDR resulted in identification of an alternative ligand-binding pocket that partially overlaps the genomic pocket described in the experimentally determined x-ray structure. Data obtained from docking five different vitamin D sterols in the genomic and alternative pockets were used to generate a receptor conformational ensemble model, providing an explanation for how VDR and possibly the estrogen receptor can have genomic and NG functionality. The VDR model is compatible with the following: (i) NG chloride channel agonism and antagonism; (ii) variable ligand-stabilized trypsin digest banding patterns; and (iii) differential transcriptional activity, employing different VDR point mutants and 1alpha,25(OH)2-vitamin D3 analogs.

Binding Sites↗

PubMed2002

On the molecular interaction between lactoferrin and the dye Red HE-3b. A novel approach for docking a charged and highly flexible molecule to protein surfaces.

Lactoferrin (Lf) is a non-heme, iron binding protein present in many physiological fluids of vertebrates where its main role is the microbicidal activity. It has been isolated by different methods, including dye-affinity chromatography. Red HE-3B is one of the most common triazinic dyes applied in protein purification, but scant knowledge is available on structural details and on the energetics of its interaction with proteins. In this work we present a computational approach useful for identifying possible binding sites for Red HE-3B in apo and holo forms of Lfs from human and bovine source. A new geometrical description of Red HE-3B is introduced which greatly simplifies the conformational analysis. This approach proved to be of particular advantage for addressing conformational ensembles of highly flexible molecules. Predictions from this analysis were correlated with experimentally observed dye-binding sites, as mapped by protection from proteolysis in Red HE-3B/Lf complexes. This method could bear relevance for the screening of possible dye-binding sites in proteins whose structure is known and as a potential tool for the design of engineered protein variants which could be purified by dye-affinity chromatography.

Amino Acid Sequence↗

PubMed2003

The role of hydrophobic interactions in amyloidogenesis: example of prion-related polypeptides.

Conversion of the non-infectious, cellular form of the prion protein (PrP(C)) to the infectious form (PrP(Sc)) is thought to be driven by an alpha-helical to beta-sheet conformational transition. To reveal the sequence determinants which encourage the transition to beta-fold, we study the synthetic peptides associated with hydrophobic conserved fragments of the N-terminal region of the prion protein. The structure of peptides in solution was probed under various thermodynamic conditions employing circular dichroism and steady state fluorescence spectroscopy as well as dye binding assays. The fluorescence methods utilized pyrene residues covalently attached to the end of the model peptides. In aqueous solutions, the structure assessments indicate the formation of metastable peptide aggregates; the molecular conformations within the peptide micelles are largely coiled. This stage in molecular assembly exists without significant beta-strand formation, i.e., before the appearance of any ordered secondary structure detectable by circular dichroism. At moderate concentrations of trifluoroethanol and/or acetonitrile, the conformational ensemble shifts towards beta-strand formation, and the population of the amorphous aggregates decreases significantly. Overall, the present data indicate that hydrophobic interactions between side chains of the peptide variants prevent, in fact, the formation of the rigid beta-sheet structures. Encouragement of beta-folds requires the destabilization of local interactions in the peptide chain, which in vivo might be possible within cell membranes as well as within partly folded molecular forms.

Amino Acid Sequence↗

PubMed2003

Many residues in cytochrome c populate alternative states under equilibrium conditions.

A curved temperature dependence of an amide proton NMR chemical shift indicates that it explores discrete alternative conformations at least 1% of the time; that is, it accesses conformations that lie within 5 kcal/mol(-1) of the ground state. The simulations presented show how curvature varies with the nature of the alternative state, and are compared to experimental results. From studies in different denaturant concentrations, it is concluded that at least 25% of residues in reduced horse cytochrome c, covering most of the protein, with the exception of the center of the N- and C-terminal helices, visit alternative states under equilibrium conditions. The conformational ensemble of the protein therefore has high structural entropy. The density of alternative states is particularly high near the heme ligand Met80, which is of interest because both redox change and the first identified stage in unfolding are associated with change in Met80 ligation. By combining theoretical and experimental approaches, it is concluded that the alternative states each comprise approximately five residues, have in general less structure than the native state, and are accessed independently. They are therefore locally unfolded structures. The locations of the alternative states match what is known of the global unfolding pathway of cytochrome c, suggesting that they may determine the pathway.

Amino Acids↗

PubMed2001

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions.

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.

Amino Acid Sequence↗

PubMed1999

Detection of nascent polyproline II helices in solution by NMR in synthetic insect kinin neuropeptide mimics containing the X-Pro-Pro-X motif.

The conformations of three synthetic peptide analogs containing the dPro-dPro-dXaa motif (dXaa = dThr, dGlu, dAsn) in aqueous solution were studied by a combination of NMR and molecular modeling simulations. The three compounds were identified from a random D-amino acid tripeptide library on the basis of their ability to either mimic or block the diuretic activity of neuropeptides of the insect kinin family. TOCSY and ROESY correlations, as well as abnormal secondary chemical shifts for protons on the D-proline residues were employed to obtain conformational ensembles consistent with the experimental NMR data for the three analogs using an in vacuo simulated annealing protocol. Similar secondary structures were found for the three molecules after refinement, in agreement with the similarities observed between their NMR spectra. Unrestrained molecular dynamics simulations with explicit water representation indicate that the structural motifs found in vacuo are stable in aqueous solution. The three analogs can be considered initiators of right-handed poly D-proline II helices, mirror images of the poly L-proline II left-handed helical motifs normally found in proline-rich proteins. The role of these secondary folds on binding of the analogs to the kinin receptors is discussed.

Insect Proteins↗

PubMed2001

Reorientational eigenmode dynamics: a combined MD/NMR relaxation analysis method for flexible parts in globular proteins.

An approach is presented for the interpretation of heteronuclear NMR spin relaxation data in mobile protein parts in terms of reorientational eigenmode dynamics. The method is based on the covariance matrix of the spatial functions of the nuclear spin interactions that cause relaxation expressed as spherical harmonics of rank 2. The approach was applied to characterize the dynamics of a loop region of ubiquitin. The covariance matrix was determined from a conformational ensemble generated by a 5 ns molecular dynamics simulation. It was found that the time correlation functions of the dominant eigenmodes decay in good approximation with a single correlation time. From the reorientational eigenmodes, their eigenvalues, and correlation times, NMR relaxation data were calculated in accordance with Bloch-Wangsness-Redfield relaxation theory and directly compared with experimental (15)N relaxation parameters. Using a fitting procedure, agreement between calculated and experimental data was improved significantly by adjusting eigenvalues and correlation times of the dominant modes. The presented procedure provides detailed information on correlated reorientational dynamics of flexible parts in globular proteins. The covariance matrix was linked to the covariance matrix of backbone dihedral angle fluctuations, allowing one to study the motional behavior of these degrees of freedom on nano- and subnanosecond time scales.

Models, Molecular↗

PubMed2006

Receptor-based pharmacophores for serotonin 5-HT7R antagonists-implications to selectivity.

A set of 31 diversified 5-HT7 receptor antagonists was automatically docked to a conformational ensemble of rhodopsin-based 5-HT7R models (flexible docking). It was found that ergolines, aporphines, and tricyclic psychotropic agents were always docked in a pocket formed by TMHs 4-6, and besides the main ionic bond with Asp3.32, they had specific interactions with Phe6.52, Phe6.51, Trp6.48, and Ser5.42. The arylpiperidine, arylpiperazine, or beta-carboline fragment of the complex ligands occupied the same pocket, whereas the terminal amide/imide part of those compounds reached the second cavity formed by TMHs 7-3 and interacted with Phe3.28, Arg7.36, and Tyr7.43. A similar orientation of selective antagonists of the group of arylsulfonamidoalkylamines was observed, that is, the sulfonamide part was located in the second pocket. Coherent docking results allowed the generation of two receptor-based pharmacophores: first containing features necessary for high 5-HT7R affinity and the other defining selectivity for this receptor subtype. The latter model indicated the importance of specific interactions with the residues of the TMHs 7-3 pocket (especially nonconserved Arg7.36) for selectivity over other monoamine GPCRs.

Animals↗

PubMed2004

One-state downhill versus conventional protein folding.

Classical protein folding invokes a cooperative transition between distinct thermodynamic states that are individually populated at equilibrium and separated by an energy barrier. It has been proposed, however, that the small protein, BBL, undergoes one-step downhill folding whereby it folds non-cooperatively to its native state without encountering an appreciable energy barrier. Only a single conformational ensemble is populated under given conditions, and so the denatured state ensemble progressively changes into the native structure. A wide dispersion of thermal denaturation midpoints that was observed for an extrinsically labelled fragment of BBL is proposed to be evidence for its one-state, downhill folding, a phenomenon that is also suggested to be functionally important for BBL and its homologues. We found, however, that thermal denaturation of unlabelled wild-type BBL was highly cooperative, with very similar transition midpoints for the melting of secondary and tertiary interactions, as well as for individual residues when monitored by NMR. Similar results were also observed for two other homologues, E3BD and POB. Further, the extrinsic fluorophores perturbed the unfolding energetics of labelled BBL, and complicated its equilibrium behaviour. One-step downhill folding may well occur for some proteins that do not have distinct folded states but not for BBL and its well-folded homologues.

Amino Acid Sequence↗

PubMed2006

Evaluation of elastic properties of atomistic DNA models.

A number of intriguing aspects in dynamics of double-helical DNA is related to the coupling between its macroscopic and microscopic states. A link between the elastic properties of long DNA chains and their atom-level dynamics can be established by comparing the worm-like chain model of polymer DNA with the conformational ensembles produced by molecular dynamics simulations. This problem is complicated by the complexity of the DNA structure, the small size of DNA fragments, and relatively short trajectory durations accessible in computer simulations of microscopic DNA dynamics. A careful study of all these aspects has been performed by using longer DNA fragments and increased durations of MD trajectories as compared to earlier such investigations. Special attention is paid to the necessary conditions and criteria of time convergence, and the possibility to increase the sampling by using constrained DNA models and simplified simulation conditions. It is found that dynamics of 25-mer duplexes with regular sequences agrees well with the worm-like chain theory and that accurate evaluation of DNA elastic parameters requires at least two turns of the double helix and approximately 20-ns duration of trajectories. Bond length and bond-angle constraints affect the estimates within numerical errors. In contrast, simplified treatment of solvation can strongly change the observed elastic parameters of DNA. The elastic parameters evaluated for AT- and GC-alternating duplexes reasonably agree with experimental data and suggest that, in different basepair sequences, the torsional and stretching elasticities vary stronger than the bending stiffness.

Base Pairing↗

PubMed2005

Design and structure of peptide and peptidomimetic antagonists of protein-protein interaction.

Peptides based on the amino acid sequences found at protein-protein interaction sites make excellent leads for antagonist development. A statistical picture of amino acids involved in protein-protein interactions indicates that proteins recognize and interact with one another through the restricted set of specialized interface amino acid residues, Pro, Ile, Tyr, Trp, Asp and Arg. These amino acids represent residues from each of the three classes of amino acids, hydrophobic, aromatic and charged, with one anionic and one cationic residue at neutral pH. The use of peptides as drug leads has been successfully used to search for antagonists of cell-surface receptors. Peptide, peptidomimetic, and non-peptide organic inhibitors of a class of cell surface receptors, the integrins, currently serve as therapeutic and diagnostic imaging agents. In this review, we discuss the structural features of protein-protein interactions as well as the design of peptides, peptidomimetics, and small organic molecules for the inhibition of protein-protein interactions. Information gained from studying inhibitors of integrin functions is now being applied to the design and testing of inhibitors of other protein-protein interactions. Most drug development progress in the past several decades has been made using the enzyme binding-pocket model of drug targets. Small molecules are designed to fit into the substrate-binding pockets of proteins based on a lock-and-key, induced-fit, or conformational ensemble model of the protein binding site. Traditionally, enzymes have been used as therapeutic drug targets because it was easier to develop rapid, sensitive screening assays, and to find low molecular weight inhibitors that blocked the active site. However, for proteins which interact with other proteins, rather than with small substrate molecules, the lack of binding pockets means that this approach will not generally succeed. There exist many diseases in which the inhibition of protein-protein interactions would provide therapeutic benefit, but there are no general methods available to address such problems. The focus of the first part of this review is to discuss the features of protein-protein interactions which may serve as general guidelines for the development and design of inhibitors for protein-protein interactions. In the second part we focus on the design of peptides (lead compounds) and their conversion into peptidomimetics or small organic molecules for the inhibition of protein-protein interactions. We draw examples from the important and emerging area of integrin-based cell adhesion and show how the principles of protein-protein interactions are followed in the discovery, optimization and usage of specific protein interface peptides as drug leads.

Amino Acids↗

PubMed1999

Engineering out motion: a surface disulfide bond alters the mobility of tryptophan 22 in cytochrome b5 as probed by time-resolved fluorescence and 1H NMR experiments.

In the accompanying paper [Storch et al. (1999) Biochemistry 38, 5054-5064] equilibrium denaturation studies and molecular dynamics (MD) simulations were used to investigate localized dynamics on the surface of cytochrome b5 (cyt b5) that result in the formation of a cleft. In those studies, an S18C:R47C disulfide mutant was engineered to inhibit cleft mobility. Temperature- and urea-induced denaturation studies revealed significant differences in Trp 22 fluorescence between the wild-type and mutant proteins. On the basis of the results, it was proposed that wild type populates a conformational ensemble that is unavailable to the disulfide mutant and is mediated by cleft mobility. As a result, the solvent accessibility of Trp 22 is decreased in S18C:R47C, suggesting that the local environment of this residue is less mobile due to the constraining effects of the disulfide on cleft dynamics. To further probe the structural effects on the local environment of Trp 22 caused by inhibition of cleft formation, we report here the results of steady-state and time-resolved fluorescence quenching, differential phase/modulation fluorescence anisotropy, and 1H NMR studies. In Trp fluorescence experiments, the Stern-Volmer quenching constant increases in wild type versus the oxidized disulfide mutant with increasing temperature. At 50 degrees C, KSV is nearly 1.5-fold greater in wild type compared to the oxidized disulfide mutant. In the reduced disulfide mutant, KSV was the same as wild type. The bimolecular collisional quenching constant, kq, for acrylamide quenching of Trp 22 increases 2.7-fold for wild type and only 1.8-fold for S18C:R47C, upon increasing the temperature from 25 to 50 degrees C. The time-resolved anisotropy decay at 25 degrees C was fit to a double-exponential decay for both the wild type and S18C:R47C. Both proteins exhibited a minor contribution from a low-amplitude fast decay, consistent with local motion of Trp 22. This component was more prevalent in the wild type, and the fractional contribution increased significantly upon raising the temperature. The fast rotational component of the S18C:R47C mutant was less sensitive to increasing temperature. A comparison of the 1H NMR monitored temperature titration of the delta-methyl protons of Ile 76 for wild type and oxidized disulfide mutant, S18C:R47C, showed a significantly smaller downfield shift for the mutant protein, suggesting that Trp 22 in the mutant protein experiences comparatively decreased cleft dynamics in core 2 at higher temperatures. Furthermore, comparison of the delta'-methyl protons of Leu 25 in the two proteins revealed a difference in the ratio of the equilibrium heme conformers of 1.2:1 for S18C:R47C versus 1.5:1 for wild type at 40 degrees C. The difference in equilibrium heme orientations between wild type and S18C:R47C suggests that the disulfide bond affects heme binding within core 1, possibly through damped cleft fluctuations. Taken together, the NMR and fluorescence studies support the proposal that an engineered disulfide bond inhibits the formation of a dynamic cleft on the surface of cyt b5.

Amino Acid Substitution↗

PubMed1997

Evidence for glutamate self-capping within a peptide helix.

The thermal dependence of the carbonyl carbon chemical shift of each residue in a helical peptide may be analyzed in terms of a two-state helix/coil transition. Such analyses generate values for the chemical shift of each residue in the helical and in the coil conformational ensembles of the peptide. The sequence dependence of the difference in these two values, termed the difference chemical shift, provides a description of the mean distribution of helicity within the helical ensemble. In this report, we improve two aspects of the procedures used to analyze prior chemical shift measurements of the helical peptide acetylW (EAAAR)3Aamide. The new difference chemical shift values for 16 of the 18 residues describe a very uniform central helical ensemble with frayed ends. However, the difference chemical shift values for the two remaining residues, alanines 03 and 08, are significantly diminished relative to this uniform distribution. Each of these two alanine residues is located i-4 to a glutamate residue. It is suggested that the difference chemical shifts for these two alanine residues are diminished by a self-capping interaction within the i + 4 glutamate residues.

Glutamic Acid↗

PubMed2003

Inverse Monte Carlo procedure for conformation determination of macromolecules.

A novel numerical method for determining the conformational structure of macromolecules is applied to idealized biomacromolecules in solution. The method computes effective inter-residue interaction potentials solely from the corresponding radial distribution functions, such as would be obtained from experimental data. The interaction potentials generate conformational ensembles that reproduce thermodynamic properties of the macromolecule (mean energy and heat capacity) in addition to the target radial distribution functions. As an evaluation of its utility in structure determination, we apply the method to a homopolymer and a heteropolymer model of a three-helix bundle protein [Zhou, Y.; Karplus, M. Proc Natl Acad Sci USA 1997, 94, 14429; Zhou, Y. et al. J Chem Phys 1997, 107, 10691] at various thermodynamic state points, including the ordered globule, disordered globule, and random coil states.

Journal Article↗