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

PubMed2003

Identifying residue-residue clashes in protein hybrids by using a second-order mean-field approach.

In this article, a second-order mean-field-based approach is introduced for characterizing the complete set of residue-residue couplings consistent with a given protein structure. This information is subsequently used to classify protein hybrids with respect to their potential to be functional based on the presenceabsence and severity of clashing residue-residue interactions. First, atomistic representations of both the native and denatured states are used to calculate rotamer-backbone, rotamer-intrinsic, and rotamer-rotamer conformational energies. Next, this complete conformational energy table is coupled with a second-order mean-field description to elucidate the probabilities of all possible rotamer-rotamer combinations in a minimum Helmholtz free-energy ensemble. Computational results for the dihydrofolate reductase family reveal correlation in substitution patterns between not only contacting but also distal second-order structural elements. Residue-residue clashes in hybrid proteins are quantified by contrasting the ensemble probabilities of protein hybrids against the ones of the original parental sequences. Good agreement with experimental data is demonstrated by superimposing these clashes against the functional crossover profiles of bidirectional incremental truncation libraries for Escherichia coli and human glycinamide ribonucleotide transformylases.

Protein Conformation↗

PubMed2006

Multiple pathways in conformational transitions of the alanine dipeptide: an application of dynamic importance sampling.

We present multiple dynamic transition pathways on the two-dimensional dihedral plane between conformational states of the alanine dipeptide. The method used in this study is dynamic importance sampling (DIMS). To perform DIMS, unbiased molecular dynamic simulations are used to generate equilibrium ensembles for the alanine dipeptide within different states. Free energy surfaces on the dihedral plane are calculated from the equilibrium simulations, and four energy minima defined from the surface are used as the starting and ending points for DIMS dynamics. The DIMS method represents an important step towards finding multiple transition pathways within complex biomolecular systems.

Alanine↗

PubMed2002

Evaluation of a novel shape-based computational filter for lead evolution: application to thrombin inhibitors.

A novel shape-feature-based computational method is described and used to rapidly filter compound libraries. The computational model, built using three-dimensional conformations of active and inactive molecules, consists of a collection of whole molecule shapes and chemical feature positions that are ranked according to their correlation with activity. A small ensemble of these shapes and features is used to filter virtual compound libraries. The method is applied to two thrombin data sets and is shown to be efficient in identifying novel scaffolds with enhanced hit rates.

Combinatorial Chemistry Techniques↗

PubMed2006

Explicit solvent DRF INDOs/CIS computations of charge transfer state energetics in a pyrenyldeoxyuridine nucleoside model.

In this work we present calculated absorption and emission spectra in acetonitrile (MeCN) solution of N-acetyl-1-aminopyrene (PAAc, a spectroscopic model compound) and N-(1-pyrenyl)-1-methyluracil-5-carboxamide (PAU(Me), a computational model for 5-(N-carboxyl-1-aminopyrenyl)-2'-deoxyuridine (PAdU)). The computational method used--the discrete reaction field approach (DRF)--combines a quantum mechanical (QM) description of the solute (here DFT and INDOs/CIS, i.e., the INDO parametrization for spectroscopy) with a classical, molecular mechanics (MM) description of the solvent molecules. The latter are modeled with point charges representing the permanent charge distribution and polarizabilities to account for many-body interactions among the solute and other solvent molecules. Molecular dynamics is used to sample the degrees of freedom of the solution around several solute conformations each in two electronic excited states. This leads to a large number of solute/solvent configurations from which 800 are selected for each excited state and collected into a single ensemble by means of proper Boltzmann averaging. DRF INDOs/CIS applied to the selected solute/solvent configurations give simulated absorption and emission band spectra--each based on 15,200 calculated transitions--that compare well with experimental results. For example, the much broader absorption and emission bands in PAdU compared with PAAc are reproduced, and the simulated emission spectra of PAU(Me) agree well with broad (380-550 nm) charge transfer (CT) emission seen for PAdU in MeCN. The observed multiexponential fluorescence decay profiles for PAdU in different polar solvents are interpreted in terms of solute/solvent conformational heterogeneity here generated in the MD simulations for PAU(Me) in MeCN. Additionally, the simulations demonstrate the mixing of the forbidden Py*+/dU*- CT states with allowed pyrenyl 1(pi,pi*) states.

Journal Article↗

PubMed2003

Comparison of the solution structures of angiotensin I & II. Implication for structure-function relationship.

Conformational analysis of angiotensin I (AI) and II (AII) peptides has been performed through 2D 1H-NMR spectroscopy in dimethylsulfoxide and 2,2,2-trifluoroethanol/H2O. The solution structural models of AI and AII have been determined in dimethylsulfoxide using NOE distance and 3JHNHalpha coupling constants. Finally, the AI family of models resulting from restrained energy minimization (REM) refinement, exhibits pairwise rmsd values for the family ensemble 0.26 +/- 0.13 A, 1.05 +/- 0.23 A, for backbone and heavy atoms, respectively, and the distance penalty function is calculated at 0.075 +/- 0.006 A2. Comparable results have been afforded for AII ensemble (rmsd values 0.30 +/- 0.22 A, 1.38 +/- 0.48 A for backbone and heavy atoms, respectively; distance penalty function is 0.029 +/- 0.003 A2). The two peptides demonstrate similar N-terminal and different C-terminal conformation as a consequence of the presence/absence of the His9-Leu10 dipeptide, which plays an important role in the different biological function of the two peptides. Other conformational variations focused on the side-chain orientation of aromatic residues, which constitute a biologically relevant hydrophobic core and whose inter-residue contacts are strong in dimethylsulfoxide and are retained even in mixed organic-aqueous media. Detailed analysis of the peptide structural features attempts to elucidate the conformational role of the C-terminal dipeptide to the different binding affinity of AI and AII towards the AT1 receptor and sets the basis for understanding the factors that might govern free- or bound-depended AII structural differentiation.

Amino Acid Sequence↗

PubMed1999

Dynamical properties of fasciculin-2.

Fasciculin-2 (FAS2) is a potent protein inhibitor of the hydrolytic enzyme acetylcholinesterase. A 2-ns isobaric-isothermal ensemble molecular dynamics simulation of this toxin was performed to examine the dynamic structural properties which may play a role in this inhibition. Conformational fluctuations of the FAS2 protein were examined by a variety of techniques to identify flexible residues and determine their characteristic motion. The tips of the toxin "finger" loops and the turn connecting loops I and II were found to fluctuate, while the rest of the protein remained fairly rigid throughout the simulation. Finally, the structural fluctuations were compared to NMR data of fluctuations on a similar timescale in a related three-finger toxin. The molecular dynamics results were in good qualitative agreement with the experimental measurements. Proteins 1999;36:447-453.

Animals↗

PubMed1998

Structures of two novel crystal forms of Naja naja naja phospholipase A2 lacking Ca2+ reveal trimeric packing.

Three crystal forms of Naja naja naja phospholipase A2 were discovered through random crystallization screening, including two heretofore uncharacterized forms. The crystallization conditions for both of these novel crystal forms are Ca(2+)-free whereas previously reported conditions include Ca2+. One of the new crystal forms has a cubic lattice in the space group P2(1)3 (a = b = c = 69.24 A), the other has an orthorhombic lattice in the space group P2(1)2(1)2(1) (a = 67.22 A, b = 73.48 A, c = 87.52 A) and a previously characterized crystal belong to the tetragonal space group P4(3)2(1)2 (a = b = 88.6 A, c = 107.4 A). The structure from the cubic crystal form has been determined to 1.8 A and refined to an R-factor of 17% while the structure from the orthorhombic form has been determined to 2.65 A and has been refined to an R-factor of 21%. The determination of the cubic structure extends the resolution to which structures of this molecule have been determined from 2.3 A to 1.8 A. The two newly determined structures, in combination with the previously determined structure, generate an informative structural ensemble from which structural changes due to Ca2+, which is required for catalysis, and the effect of crystal contacts on side-chain conformations and oligomeric association can be inferred. Both of the newly determined structures reveal a trimeric oligomer as observed in the tetragonal structure; this appears to be a unique feature of the Naja naja naja enzyme.

Animals↗

PubMed2004

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization.

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.

Amyloid beta-Peptides↗

PubMed2005

Sequence determinants of a transmembrane proton channel: an inverse relationship between stability and function.

The driving forces behind the folding processes of integral membrane proteins after insertion into the bilayer, is currently under debate. The M2 protein from the influenza A virus is an ideal system to study lateral association of transmembrane helices. Its proton selective channel is essential for virus functioning and a target of the drug amantadine. A 25 residue transmembrane fragment of M2, M2TM, forms a four-helix bundle in vivo and in various detergents and phospholipid bilayers. Presented here are the energetic consequences for mutations made to the helix/helix interfaces of the M2TM tetramer. Analytical ultracentrifugation has been used to determine the effect of ten single-site mutations, to either alanine or phenylalanine, on the oligomeric state and the free energy of M2TM in the absence and the presence of amantadine. It was expected that many of these mutations would perturb the M2TM stability and tetrameric integrity. Interestingly, none of the mutations destabilize tetramerization. This finding suggests that M2 sacrifices stability to preserve its functions, which require rapid and specific interchange between distinct conformations involved in gating and proton conduction. Mutations might therefore restrict the full range of conformations by stabilizing a given native or non-native conformational state. In order to assess one specific conformation of the tetramer, we measured the binding of amantadine to the resting state of the channel, and examined the overall free energy of assembly of the amantadine bound tetramer. All of the mutations destabilized amantadine binding or were isoenergetic. We also find that large to small residue changes destabilize the amantadine bound tetramer whereas mutations to side-chains of similar volume stabilize this conformation. A structural model of the amantadine bound state of M2TM was generated using a novel protocol that optimizes a structure for an ensemble of neutral and disruptive mutations. The model structure is consistent with the mutational data.

Amantadine↗

PubMed2000

Rapid internal dynamics of BPTI is insensitive to pressure. (15)N spin relaxation at 2 kbar.

Pressure effects on the backbone dynamics of a native basic pancreatic trypsin inhibitor (BPTI) have been measured by (15)N spin relaxation and chemical shifts at 30 and 2000 bar. The experiments utilized the on-line variable pressure cell nuclear magnetic resonance system on (15)N-uniformly labeled BPTI at a proton frequency of 750.13 MHz at 36 degrees C. Longitudinal (R(1)) and transverse (R(2)) (15)N relaxation times and ((1)H)-(15)N nuclear Overhauser effects were measured for 41 protonated backbone nitrogens at both pressures. The model free analysis of the internal dynamics gave order parameters for individual H-N vectors at both pressures. The results indicate that rapid internal dynamics in the ps-ns range for the polypeptide backbone is not significantly affected by pressure in the range between 30 bar and 2 kbar. The result is consistent with the linear pressure dependence of (1)H and (15)N chemical shifts of BPTI, which suggests that local compressibilities and amplitudes of associated conformational fluctuation are nearly invariant in the same pressure range. Overall, we conclude that at 2 kbar BPTI remains within the same native ensemble as at 1 bar, with a small shift of population from that at 1 bar.

Aprotinin↗

PubMed2005

Dense orientationally ordered states of a single semiflexible macromolecule: an expanded ensemble Monte Carlo simulation.

Using a coarse-grained model we perform a Monte Carlo simulation of the state behavior of an individual semiflexible macromolecule. Chains consisting of N = 256 and 512 monomer units have been investigated. A recently proposed enhanced sampling Monte Carlo technique for the bond fluctuation model in an expanded ensemble in four-dimensional coordinate space was applied. The algorithm allows one to accelerate the sampling of statistically independent three-dimensional conformations in a dense globular state. We found that the temperature of the intraglobular liquid-solid transition decreases with increasing chain stiffness. We have investigated the possible intraglobular orientationally ordered (i.e., liquid-crystalline) structures and obtained a diagram of states for chains consisting of N = 256 monomer units. This diagram contains regions of stability of coil, two spherical globules (liquid and solid), and rod-like globule conformations. Transitions between the globular states are rounded first-order ones since the states of liquid, solid, and cylinder-like globules do have different internal symmetry.

Journal Article↗

PubMed2005

Single-molecule Forster resonance energy transfer study of protein dynamics under denaturing conditions.

Proteins are highly complex systems, exhibiting a substantial degree of structural variability in their folded state. In the presence of denaturants, the heterogeneity is greatly enhanced, and fluctuations among vast numbers of folded and unfolded conformations occur via many different pathways. Here, we have studied the structure and dynamics of the small enzyme ribonuclease HI (RNase H) in the presence of the chemical denaturant guanidinium chloride (GdmCl) using single-molecule fluorescence microscopy, with a particular focus on the characterization of the unfolded-state ensemble. A dye pair was specifically attached to the enzyme to measure structural changes through Förster resonance energy transfer (FRET). Enzyme immobilization on star-polymer surfaces that were specially developed for negligible interaction with folded and unfolded proteins enabled us to monitor conformational changes of individual proteins for several hundred seconds. FRET efficiency histograms were calculated from confocal scan images. They showed an expansion of the unfolded proteins with increasing GdmCl concentration. Cross-correlation analysis of donor and acceptor fluorescence intensity time traces from single molecules revealed reconfiguration of the polypeptide chain on a timescale of approximately equal to 20 micros at 1.7 M GdmCl. Slow conformational dynamics gave rise to characteristic, stepwise FRET efficiency changes. Transitions between folded and unfolded enzyme molecules occurred on the 100-s timescale, in excellent agreement with bulk denaturation experiments. Transitions between unfolded conformations were more frequent, with characteristic times of approximately equal to 2 s. These data were analyzed to obtain information on the free energy landscape of RNase H in the presence of chemical denaturants.

Fluorescence↗

PubMed2002

Conformational flexibility of three cytoplasmic segments of the angiotensin II AT1A receptor: a circular dichroism and fluorescence spectroscopy study.

The conformation of three synthetic peptides encompassing the proximal and distal half of the third intracellular loop (Ni3 and Ci3) and a portion of the cytoplasmic tail (fCT) of the angiotensin II AT1A receptor has been studied using circular dischroism and fluorescence spectroscopies. The results show that the conformation of the peptides is modulated in various ways by the environmental conditions (pH, ionic strength and dielectric constant). Indeed, Ni3 and fCT fold into helical structures that possess distinct stability and polarity due to the diverse forces involved: mainly polar interactions in the first case and a combination of polar and hydrophobic interactions in the second. The presence of these various features also produce distinct intermolecular interactions. Ci3, instead, exists as an ensemble of partially folded states in equilibrium. Since the corresponding regions of the angiotensin II AT1A receptor are known to play an important role in the receptor function, due to their ability to undergo conformational changes, these data provide some new clues about their different conformational plasticity.

Amino Acid Sequence↗

PubMed1995

Interactions among polymorphic and conserved residues in MHC class II proteins affect MHC-peptide conformation and T cell recognition.

The structural basis for MHC-restricted T cell recognition of the N-terminal peptide of myelin basic protein (MBP Ac1-11) presented by two mouse class II alleles, Ak and Au, was examined, focusing on the roles of A beta chain polymorphic residues 38 beta (in the beta sheet) and 61 beta (in the alpha helix) in controlling the responses of panels of Ak- and Au-restricted T cell hybridomas. Despite the conservative nature of the substitutions at 38 beta (k = Val, u = Leu) and 61 beta (k = Trp, u = Tyr), transfectants expressing Ak or Au proteins carrying allelic substitutions at 38 beta and/or 61 beta gave dramatically reduced T cell responses. The modest reduction in peptide binding detected using a biotinylated MBP peptide analog appear insufficient to explain the reduced responses, suggesting that changes at 38,61 beta create conformational changes in the MHC-peptide complex. The impact of allelic substitutions at 38,61 beta on T cell responses is also modulated by other residues differing between Ak and Au. To explore the structural basis for these phenomena, protein models were developed of the Ak, Au and 38,61 beta mutant proteins using self-consistent ensemble optimization methodologies. Substitutions of the alternative allelic residue at 38 beta and/or 61 beta, which are in van der Waals contact, change the configuration of this region of the peptide-binding groove, and potentially might affect the conformation of the bound peptide and its hydrogen-bonding to residue 61 beta. The models predict that this region of the groove is markedly altered by allelic differences at A beta residue 9 beta (k = His, u = Val) which determine the position of the side-chain of Tyr30 beta, adjacent to residues 38 beta and 61 beta. Thus, interactions among polymorphic and conserved residues control the antigen presentations functions of MHC class II proteins.

Alleles↗

PubMed2001

Enantioselectivity in Candida antarctica lipase B: a molecular dynamics study.

A major problem in predicting the enantioselectivity of an enzyme toward substrate molecules is that even high selectivity toward one substrate enantiomer over the other corresponds to a very small difference in free energy. However, total free energies in enzyme-substrate systems are very large and fluctuate significantly because of general protein motion. Candida antarctica lipase B (CALB), a serine hydrolase, displays enantioselectivity toward secondary alcohols. Here, we present a modeling study where the aim has been to develop a molecular dynamics-based methodology for the prediction of enantioselectivity in CALB. The substrates modeled (seven in total) were 3-methyl-2-butanol with various aliphatic carboxylic acids and also 2-butanol, as well as 3,3-dimethyl-2-butanol with octanoic acid. The tetrahedral reaction intermediate was used as a model of the transition state. Investigative analyses were performed on ensembles of nonminimized structures and focused on the potential energies of a number of subsets within the modeled systems to determine which specific regions are important for the prediction of enantioselectivity. One category of subset was based on atoms that make up the core structural elements of the transition state. We considered that a more favorable energetic conformation of such a subset should relate to a greater likelihood for catalysis to occur, thus reflecting higher selectivity. The results of this study conveyed that the use of this type of subset was viable for the analysis of structural ensembles and yielded good predictions of enantioselectivity.

Butanols↗

PubMed2006

Folding pathway dependence on energetic frustration and interaction heterogeneity for a three-dimensional hydrophobic protein model.

Monte Carlo simulations of a hydrophobic protein model of 40 monomers in the cubic lattice are used to explore the effect of energetic frustration and interaction heterogeneity on its folding pathway. The folding pathway is described by the dependence of relevant conformational averages on an appropriate reaction coordinate, pfold, defined as the probability for a given conformation to reach the native structure before unfolding. We compare the energetically frustrated and heterogeneous hydrophobic potential, according to which individual monomers have a higher or lower tendency to form contacts unspecifically depending on their hydrophobicities, to an unfrustrated homogeneous Go-type potential with uniformly attractive native interactions and neutral non-native interactions (called Go1 in this study), and to an unfrustrated heterogeneous potential with neutral non-native interactions and native interactions having the same energy as the hydrophobic potential (called Go2 in this study). Folding kinetics are slowed down dramatically when energetic frustration increases, as expected and previously observed in a two-dimensional model. Contrary to our previous results in two dimensions, however, it appears that the folding pathway and transition state ensemble can be significantly dependent on the energy function used to stabilize the native structure. The sequence of events along the reaction coordinate, or the order along this coordinate in which different regions of the native conformation become structured, turns out to be similar for the hydrophobic and Go2 potentials, but with analogous events tending to occur at lower pfold values in the first case. In particular, the transition state obtained from the ensemble around pfold = 0.5 is more structured for the hydrophobic potential. For Go1, not only the transition state ensemble but the order of events itself is modified, suggesting that interaction heterogeneity, in addition to energetic frustration, can have significant effects on the folding mechanism, most likely by modifying the probability of different contacts in the unfolded state, the starting point for the folding reaction. Although based on a simple model, these results provide interesting insight into how sequence-dependent switching between folding pathways might occur in real proteins.

Kinetics↗

PubMed2004

Self-assembly of flexible supramolecular metallacyclic ensembles: structures and adsorption properties of their nanoporous crystalline frameworks.

The syntheses, structures, and N2 adsorption properties of six new supramolecular metallacycles are reported. Flexible ditopic linkers, 1-4, with systematically varied lengths and conformational degrees of freedom were synthesized utilizing ester linkages. They were used in combination with (dppp)M(OTf)2, where M = Pt(II) and Pd(II), and cis-(Me3P)2Pt(OTf)2 to form flexible supramolecular metallacycles 5-10 in 88-98% isolated yields. Their structures were characterized via multinuclear NMR and X-ray crystallography. The metallacycles stack to form porous structures in the crystalline state. The pore dimensions depend on both the phosphorus ligands attached to the metals and the flexible linkers. Adsorption studies on the porous materials show that 5a, 6, 8, and 9 held 11.7, 16.5, 5.7, and 6.8 cm3/g STP of N2 at 77 K, respectively. A guest-exchange study with nitromethane and toluene reveals that the nanopore in 5 is flexible, a property which was transferred from the linker to the supramolecular structure in the solid state.

Adsorption↗

PubMed2006

A correlation-based method for the enhancement of scoring functions on funnel-shaped energy landscapes.

A correlation-based approach is introduced for enhancing the ability of structure-scoring methods to identify and distinguish native-like conformations. The proposed method relies on a funnel-shaped scoring function that decreases steadily toward the native state. It takes advantage of the idea that the structure from a given ensemble that is closest to the native basin leads to the highest correlation coefficient between a given score and distance to that structure as an approximation of the native state for the entire ensemble. The method is applied successfully to a number of different test cases that demonstrate substantial improvements in the correlation of the score with the distance from the true native state but also result in the selection of more native-like structures compared to the original score.

Algorithms↗