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

PubMed2005

Interplay of charge distribution and conformation in peptides: comparison of theory and experiment.

We assessed the correlation between charge distribution and conformation of flexible peptides by comparing the theoretically calculated potentiometric-titration curves of two model peptides, Ac-Lys5-NHMe (a model of poly-L-lysine) and Ac-Lys-Ala11-Lys-Gly2-Tyr-NH2 (P1) in water and methanol, with the experimental curves. The calculation procedure consisted of three steps: (i) global conformational search of the peptide under study using the electrostatically driven Monte Carlo (EDMC) method with the empirical conformational energy program for peptides (ECEPP)/3 force field plus the surface-hydration (SRFOPT) or the generalized Born surface area (GBSA) solvation model as well as a molecular dynamics method with the assisted model building and energy refinement (AMBER)99/GBSA force field; (ii) reevaluation of the energy in the pH range considered by using the modified Poisson-Boltzmann approach and taking into account all possible protonation microstates of each conformation, and (iii) calculation of the average degree of protonation of the peptide at a given pH value by Boltzmann averaging over conformations. For Ac-Lys5-NHMe, the computed titration curve agrees qualitatively with the experimental curve of poly-L-lysine in 95% methanol. The experimental titration curves of peptide P1 in water and methanol indicate a remarkable downshift of the first pK(a) value compared to the values for reference compounds (n-butylamine and phenol, respectively), suggesting the presence of a hydrogen bond between the tyrosine hydroxyl oxygen and the H(epsilon) proton of a protonated lysine side chain. The theoretical titration curves agree well with the experimental curves, if conformations with such hydrogen bonds constitute a significant part of the ensemble; otherwise, the theory predicts too small a downward pH shift.

Electrochemistry↗

PubMed2001

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations.

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.

Circular Dichroism↗

PubMed2002

Interaction of dihydrofolate reductase with methotrexate: ensemble and single-molecule kinetics.

The thermodynamics and kinetics of the interaction of dihydrofolate reductase (DHFR) with methotrexate have been studied by using fluorescence, stopped-flow, and single-molecule methods. DHFR was modified to permit the covalent addition of a fluorescent molecule, Alexa 488, and a biotin at the N terminus of the molecule. The fluorescent molecule was placed on a protein loop that closes over methotrexate when binding occurs, thus causing a quenching of the fluorescence. The biotin was used to attach the enzyme in an active form to a glass surface for single-molecule studies. The equilibrium dissociation constant for the binding of methotrexate to the enzyme is 9.5 nM. The stopped-flow studies revealed that methotrexate binds to two different conformations of the enzyme, and the association and dissociation rate constants were determined. The single-molecule investigation revealed a conformational change in the enzyme-methotrexate complex that was not observed in the stopped-flow studies. The ensemble averaged rate constants for this conformation change in both directions is about 2-4 s(-1) and is attributed to the opening and closing of the enzyme loop over the bound methotrexate. Thus the mechanism of methotrexate binding to DHFR involves multiple steps and protein conformational changes.

Base Sequence↗

PubMed2005

Effects of conformational distributions on sigma profiles in COSMO theories.

The charge density or sigma profile of a solute molecule is an essential component in COSMO (conductor-like screen model) based solvation theories, and its generation depends on the molecular conformation used. The usual procedure is to determine the conformation of an isolated molecule, and assume that this is unchanged when the molecule is placed in solution. In this paper, the conformations of 1-hexanol and 2-methoxy-ethanol in both the liquid and vapor phases obtained from Gibbs ensemble simulation and from an isolated-molecule quantum DFT optimization are used to determine the effect of realistic conformation differences on COSMO-based properties predictions. In particular, the vapor pressure at the normal boiling temperature and the binary mixture VLE (vapor-liquid equilibrium) predictions obtained using different conformations are investigated. The results show that the sigma profile for 1-hexanol varies only slightly using the different conformations, while the sigma profile of 2-methoxy-ethanol shows a significant difference between the liquid and vapor phases. Consequently, the vapor pressure predictions for 1-hexanol are similar regardless of the manner in which the conformation population was obtained, while there is a larger difference for 2-methoxy-ethanol depending on whether the liquid or vapor conformations from simulation or the DFT-optimized structure is used. These differences in predictions are seen to be largely due to differences in the ideal solvation energy term. In mixture VLE calculations involving 1-hexanol, we again see that there is little difference in the phase equilibrium predictions among the different conformations, while for the mixture with 2-methoxy-ethanol, the differences in the sigma profiles lead to a more noticeable, though not significant, difference in the phase equilibrium predictions.

Journal Article↗

PubMed2004

Protein structure, stability and solubility in water and other solvents.

Proteins carry out the most difficult tasks in living cells. They do so by interacting specifically with other molecules. This requires that they fold to a unique, globular conformation that is only marginally more stable than the large ensemble of unfolded states. The folded state is stabilized mainly by the burial and tight packing of over 80% of the peptide groups and non-polar side chains. If life as we know it is to exist in a solvent other than water, the folded state must be stable and soluble in the new solvent. Our analysis suggests that proteins will be unstable in most polar solvents such as ethanol, extremely stable in non-polar solvents such as cyclohexane, and even more stable in a vacuum. Our solubility studies suggest that protein solubility will be markedly lower in polar solvents such as ethanol and that proteins will be essentially insoluble in non-polar solvents such as cyclohexane. For these and other reasons it seems unlikely that the life we know could exist in any solvent system other than water.

Cyclohexanes↗

PubMed2006

Amplitudes and directions of internal protein motions from a JAM analysis of 15N relaxation data.

A method has been developed for characterizing dynamic structures of proteins in solution by using nuclear magnetic resonance (NMR) restraints and 15N relaxation data. This method is based on the concept of the jumping-among-minima (JAM) model. In this model we assume that protein dynamics can be described on the basis of conformational substates, and involves intra- and inter-substate motion. A set of substates is created by picking energy-minimized conformations from the conformational space consistent with the geometric NMR restraints. Intra-substate motions, which occur on the timescale of approximately 10 ps, are simulated with molecular dynamics (MD) calculations with force-field energy terms. Statistical weights of the conformational substates are determined to reproduce the NMR relaxation parameters. The refinement procedure consists of four stages: (i) determination of the ensemble of structures that satisfy NMR restraints, (ii) determination of intra-substate fluctuation, (iii) determination of statistical weights of conformational substates to reproduce model-free relaxation parameters, and (iv) analysis of the resulting dynamic structure to determine amplitudes and directions of internal protein motions. This method was employed to investigate structure and dynamics of the adhesion domain of human CD2 (hCD2) in solution. Two major collective modes, whose contributions to atomic mean-square fluctuations are 77.1% in total, are identified by the refinement. The first mode is interpreted as a rigid-body motion of a protein segment consisting of a part of the B--C loop, a part of the F strand, and the F--G loop. Another type of smaller-amplitude mode is indicated for the C'--C'' loop. The motions affect primarily the curvature of the slightly concave counterreceptor-binding site and represent transitions between a concave (closed) and flat (open) binding face. By comparing the ensemble of structures in solution to the complex structure with counterreceptor CD58, we found that these two types of motions resemble the change upon counterreceptor binding.

Models, Molecular↗

PubMed2006

Protein design simulations suggest that side-chain conformational entropy is not a strong determinant of amino acid environmental preferences.

Loss of side-chain conformational entropy is an important force opposing protein folding and the relative preferences of the amino acids for being buried or solvent exposed may be partially determined by which amino acids lose more side-chain entropy when placed in the core of a protein. To investigate these preferences, we have incorporated explicit modeling of side-chain entropy into the protein design algorithm, RosettaDesign. In the standard version of the program, the energy of a particular sequence for a fixed backbone depends only on the lowest energy side-chain conformations that can be identified for that sequence. In the new model, the free energy of a single amino acid sequence is calculated by evaluating the average energy and entropy of an ensemble of structures generated by Monte Carlo sampling of amino acid side-chain conformations. To evaluate the impact of including explicit side-chain entropy, sequences were designed for 110 native protein backbones with and without the entropy model. In general, the differences between the two sets of sequences are modest, with the largest changes being observed for the longer amino acids: methionine and arginine. Overall, the identity between the designed sequences and the native sequences does not increase with the addition of entropy, unlike what is observed when other key terms are added to the model (hydrogen bonding, Lennard-Jones energies, and solvation energies). These results suggest that side-chain conformational entropy has a relatively small role in determining the preferred amino acid at each residue position in a protein.

Amino Acid Sequence↗

PubMed2006

Use of an induced fit receptor structure in virtual screening.

Structured-based drug design has traditionally relied on a single receptor structure as a target for docking and screening studies. However, it has become increasingly clear that in many cases where protein flexibility is an issue, it is critical to accurately model ligand-induced receptor movement in order to obtain high enrichment factors. We present a novel protein-ligand docking method that accounts for both ligand and receptor flexibility and accurately predicts the conformation of protein-ligand binding complexes. This method can generate viable receptor ensembles that can be used in virtual database screens.

Databases, Protein↗

PubMed2005

Local conformational fluctuations can modulate the coupling between proton binding and global structural transitions in proteins.

Local conformational fluctuations in proteins can affect the coupling between ligand binding and global structural transitions. This finding was established by monitoring quantitatively how the population distribution in the ensemble of microstates of staphylococcal nuclease was affected by proton binding. Analysis of acid unfolding and proton-binding data with an ensemble-based model suggests that local fluctuations: (i) can be effective modulators of ligand-binding affinities, (ii) are important determinants of the cooperativity of ligand-driven global structural transitions, and (iii) are well represented thermodynamically as local unfolding processes. These studies illustrate how an ensemble-based description of proteins can be used to describe quantitatively the interdependence of local conformational fluctuations, ligand-binding processes, and global structural transitions. This level of understanding of the relationship between conformation, energy, and dynamics is required for a detailed mechanistic understanding of allostery, cooperativity, and other complex functional and regulatory properties of macromolecules.

Binding Sites↗

PubMed2005

Probing single-molecule protein conformational dynamics.

Protein conformational fluctuations and dynamics, often complex and associated with inhomogeneities, play a crucial role in biomolecular functions. It is extremely difficult to characterize such inhomogeneous dynamics in an ensemble-averaged measurement, especially when the proteins are involved in multiple-step, multiple-conformation complex chemical interactions and transformations, such as in enzymatic reactions, protein-protein interactions, and ion-channel membrane protein processes. Alternatively, single-molecule spectroscopy is a powerful approach to probing and analyzing protein conformational dynamics in real time.

Bacteriophage T4↗

PubMed1999

The structure and dynamics of rat apo-cellular retinol-binding protein II in solution: comparison with the X-ray structure.

The structure and dynamics of rat apo-cellular retinol binding protein II (apo-CRBP II) in solution has been determined by multidimensional NMR analysis of uniformly enriched recombinant rat 13C, 15N-apo-CRBP II and 15N-apo-CRBP II. The final ensemble of 24 NMR structures has been calculated from 3274 conformational restraints or 24.4 restraints/residue. The average root-mean-square deviation of the backbone atoms for the final 24 structures relative to their mean structure is 1.06 A. Although the average solution structure is very similar to the crystal structure, it differs at the putative entrance to the binding cavity, which is formed by the helix-turn-helix motif, the betaC-betaD turn and the betaE-betaF turn. The mean coordinates of the main-chain atoms of amino acid residues 28-38 are displaced in the solution structure relative to the crystal structure. The side-chain of F58, located on the betaC-betaD turn, is reoriented such that it interacts with L37 and no longer blocks entry into the ligand-binding pocket. Residues 28-35, which form the second helix of the helix-turn-helix motif in the crystal structure, do not exhibit a helical conformation in the solution structure. The solution structure of apo-CRBP II exhibits discrete regions of backbone disorder which are most pronounced at residues 28-32, 37-38 and 73-76 in the betaE-betaF turn as evaluated by the consensus chemical shift index, the root-mean-square deviation, amide 1H exchange rates and 15N relaxation studies. These studies indicate that fluctuations in protein conformation occur on the microseconds to ms time-scale in these regions of the protein. Some of these exchange processes can be directly observed in the three-dimensional 15N-resolved NOESY spectrum. These results suggest that in solution, apo-CRBP II undergoes conformational changes on the microseconds to ms time-scale which result in increased access to the binding cavity.

Amides↗

PubMed2005

Protein folding and the organization of the protein topology universe.

The mechanism by which proteins fold to their native states has been the focus of intense research in recent years. The rate-limiting event in the folding reaction is the formation of a conformation in a set known as the transition-state ensemble. The structural features present within such ensembles have now been analysed for a series of proteins using data from a combination of biochemical and biophysical experiments together with computer-simulation methods. These studies show that the topology of the transition state is determined by a set of interactions involving a small number of key residues and, in addition, that the topology of the transition state is closer to that of the native state than to that of any other fold in the protein universe. Here, we review the evidence for these conclusions and suggest a molecular mechanism that rationalizes these findings by presenting a view of protein folds that is based on the topological features of the polypeptide backbone, rather than the conventional view that depends on the arrangement of different types of secondary-structure elements. By linking the folding process to the organization of the protein structure universe, we propose an explanation for the overwhelming importance of topology in the transition states for protein folding.

Animals↗

PubMed2001

Dissecting glycoprotein quality control in the secretory pathway.

In the early secretory pathway, asparagine-linked glycosylation facilitates the conformational maturation of diverse polypeptides by promoting their physical engagement with the glycoprotein-folding machinery. Misfolded glycoproteins are selectively eliminated from the endoplasmic reticulum by a stringent process of conformation-based quality control. Recent studies indicate that a small ensemble of oligosaccharide-processing enzymes and lectins use the asparagine-linked appendage to orchestrate the selective disposal of numerous transport-defective glycoproteins from the early secretory pathway. The glycan-based disposal system functions as an evolutionarily conserved terminal checkpoint in eukaryote genome expression. That the mechanisms by which glycoprotein substrates are recruited for degradation diverge at the level of signal recognition reflects a previously unappreciated component of cellular differentiation in higher eukaryotes.

Animals↗

PubMed2004

Biomolecular free energy profiles by a shooting/umbrella sampling protocol, "BOLAS".

We develop an efficient technique for computing free energies corresponding to conformational transitions in complex systems by combining a Monte Carlo ensemble of trajectories generated by the shooting algorithm with umbrella sampling. Motivated by the transition path sampling method, our scheme "BOLAS" (named after a cowboy's lasso) preserves microscopic reversibility and leads to the correct equilibrium distribution. This makes possible computation of free energy profiles along complex reaction coordinates for biomolecular systems with a lower systematic error compared to traditional, force-biased umbrella sampling protocols. We demonstrate the validity of BOLAS for a bistable potential, and illustrate the method's scope with an application to the sugar repuckering transition in a solvated deoxyadenosine molecule.

Algorithms↗

PubMed2003

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations.

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.

Computer Simulation↗

PubMed1995

Comparison of rotation models for describing DNA conformations: application to static and polymorphic forms.

A new method, based on a space-fixed rotation axis, or local helix axis, is proposed for the calculation of the relative orientation variables for a sequence of base pairs. With this method, orientation variables are determined through the rotation of a base pair about this axis. These variables uniquely determine a set of helical variables, similar to the roll, tilt, and twist, commonly used for a description of spatial orientations of internally rigid base pairs. The proposed identification of roll and tilt with the direction cosines of the space-fixed rotation axis agrees well with their customary definitions as the openings of the angles between adjoining base pairs toward the minor groove and toward the ascending (5' to 3') backbone strand, respectively. These new variables permit a more direct physical comprehension of DNA conformations and also the behavior of self-complementary sequences. These direction cosines, together with the rotation angle about the space-fixed axis, form a set of three independent orientation variables of the bases that afford some advantages over the variously defined twist, roll, and tilt angles, either for static or average forms. An example for the static form of these variables is shown through their use to interpret crystal coordinates. An example for the average of orientation variables is based on statistical calculations. In this example, the orientation variables, together with the translational variables that describe the relative displacements of a pair of adjacent base pairs, form a canonically distributed ensemble in phase space spanned by these variables. Two sets of conformational variables are generated by using two different methods for performing rotation operations on the sequences of base pairs. The first method is based on the new single rotation about a space-fixed axis of rotation. This space-fixed axis of rotation is, in fact, the local helical axis as constructed previously by others. The second method is based on three consecutive rotations by Euler angles. Because of large flexibilities and anisotropies along various conformational variables of DNA base pairs, the two sets of generated conformational variables, based on these two different methods of performing rotation operations, lead to slightly different sets of structurally different, but energetically equivalent, spatial arrangements of the base pairs.

Base Composition↗

PubMed2003

The role of alpha-, 3(10)-, and pi-helix in helix-->coil transitions.

The conformational equilibrium between 3(10)- and alpha-helical structure has been studied via high-resolution NMR spectroscopy by Millhauser and coworkers using the MW peptide Ac-AMAAKAWAAKA AAARA-NH2. Their 750-MHz nuclear Overhauser effect spectroscopy (NOESY) spectra were interpreted to reflect appreciable populations of 3(10)-helix throughout the peptide, with the greatest contribution at the N and C termini. The presence of simultaneous alphaN(i,i + 2) and alphaN(i,i + 4) NOE cross-peaks was proposed to represent conformational averaging between 3(10)- and alpha-helical structures. In this study, we describe 25-nsec molecular dynamics simulations of the MW peptide at 298 K, using both an 8 A and a 10 A force-shifted nonbonded cutoff. The ensemble averages of both simulations are in reasonable agreement with the experimental helical content from circular dichroism (CD), the (3)J(HNalpha) coupling constants, and the 57 observed NOEs. Analysis of the structures from both simulations revealed very little formation of contiguous i --> i + 3 hydrogen bonds (3(10)-helix); however, there was a large population of bifurcated i --> i + 3 and i --> i + 4 alpha-helical hydrogen bonds. In addition, both simulations contained considerable populations of pi-helix (i --> i + 5 hydrogen bonds). Individual turns formed over residues 1-9, which we predict contribute to the intensities of the experimentally observed alphaN(i,i + 2) NOEs. Here we show how sampling of both folded and unfolded structures can provide a structural framework for deconvolution of the conformational contributions to experimental ensemble averages.

Amino Acid Sequence↗

PubMed2006

Sub-angstrom conformational changes of a single molecule captured by AFM variance analysis.

A system's equilibrium variance can be analyzed to probe its underlying dynamics at higher resolution. Here, using single-molecule atomic-force microscope techniques, we show how the variance in the length of a single dextran molecule can be used to establish thermodynamic equilibrium and to detect conformational changes not directly observable with other methods. Dextran is comprised of a chain of pyranose rings that each undergoes an Angstrom-scale transition from a chair to boat conformation under a stretching force. Our analysis of the variance of the molecule's fluctuations verifies equilibrium throughout the force-extension curve, consistent with the expected thermodynamic ensemble. This validates further analysis of the variance in the transition region, which reveals an intermediate conformation between the chair and the boat on the sub-Angstrom scale. Our test of thermal equilibrium as well as our variance analysis can be readily extended to a wide variety of molecules, including proteins.

Analysis of Variance↗