Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Pliability”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20Linked to original sources

Change in protein flexibility upon complex formation: analysis of Ras-Raf using molecular dynamics and a molecular framework approach.

Changes in flexibility upon protein-protein complex formation of H-Ras and the Ras-binding domain of C-Raf1 have been investigated using the molecular framework approach FIRST (Floppy Inclusion and Rigid Substructure Topology) and molecular dynamics simulations (MD) of in total approximately 35 ns length. In a computational time of about one second, FIRST identifies flexible and rigid regions in a single, static three-dimensional molecular framework, whose vertices represent protein atoms and whose edges represent covalent and non-covalent (hydrogen bond and hydrophobic) constraints and fixed bond angles within the protein. The two methods show a very good agreement with respect to the identification of changes in flexibility in both binding partners on a local scale. This implies that flexibility can be successfully predicted by identifying which bonds limit motion within a molecule and how they are coupled. In particular, as identified by MD, the beta-sheet in Raf shows considerably more pronounced orientational correlations in the bound state compared to the unbound state. Similarly, FIRST assigns the beta-sheet to the largest rigid cluster of the complex. Interestingly, FIRST allows us to identify that interactions across the interface (but not conformational changes upon complex formation) result in the observed rigidification. Since regions of the beta-sheet of Raf that do not interact directly with Ras become rigidified, this also demonstrates the long-range aspect to rigidity percolation. Possible implications of the change of flexibility of the Ras-binding domain of Raf on the activation of Raf upon complex formation are discussed. Finally, the sensitivity of FIRST results with respect to the representation of non-covalent interactions used as constraints is probed.

Algorithms↗

Predictions of protein flexibility: first-order measures.

The normal modes of a molecule are utilized, in conjunction with classical conformal vector field theory, to define a function that measures the capability of the molecule to deform at each of its residues. An efficient algorithm is presented to calculate the local chain deformability from the set of normal modes of vibration. This is done by considering each mode as an off-grid sample of a deformation vector field. Predictions of deformability are compared with experimental data in the form of dihedral angle differences between two conformations of ten kinases by using a modified correlation function. Deformability calculations correlate well with experimental results and validate the applicability of this method to protein flexibility predictions.

Chemistry, Physical↗

Modeling correlated main-chain motions in proteins for flexible molecular recognition.

We describe a new method for modeling protein and ligand main-chain flexibility, and show its ability to model flexible molecular recognition. The goal is to sample the full conformational space, including large-scale motions that typically cannot be reached in molecular dynamics simulations due to the computational intensity, as well as conformations that have not been observed yet by crystallography or NMR. A secondary goal is to assess the degree of flexibility consistent with protein-ligand recognition. Flexibility analysis of the target protein is performed using the graph-theoretic algorithm FIRST, which also identifies coupled networks of covalent and noncovalent bonds within the protein. The available conformations of the flexible regions are then explored with ROCK by random-walk sampling of the rotatable bonds. ROCK explores correlated motions by only sampling dihedral angles that preserve the coupled bond networks in the protein and generates conformers with good stereochemistry, without using a computationally expensive potential function. A representative set of the conformational ensemble generated this way can be used as targets for docking with SLIDE, which handles the flexibility of protein and ligand side-chains. The realism of this protein main-chain conformational sampling is assessed by comparison with time-resolved NMR studies of cyclophilin A motions. ROCK is also effective for modeling the flexibility of large cyclic and polycyclic ligands, as demonstrated for cyclosporin and zearalenol. The use of this combined approach to perform docking with main-chain flexibility is illustrated for the cyclophilin A-cyclosporin complex and the estrogen receptor in complex with zearalenol, while addressing the question of how much flexibility is allowed without hindering molecular recognition.

Biochemical Phenomena↗

Toward an atomistic model for predicting transcription-factor binding sites.

Identifying the specific DNA-binding sites of transcription-factor proteins is essential to understanding the regulation of gene expression in the cell. Bioinformatics approaches are fast compared to experiments, but require prior knowledge of multiple binding sites for each protein. Here, we present an atomistic force-field method to predict binding sites based only on the X-ray structure of a related bound complex. Specific flexible contacts between the protein and DNA are modeled by a library of amino acid side-chain rotamers. Using the example of the mouse transcription factor, Zif268, a well-studied zinc-finger protein, we show that the protein sequence alone, without the detailed experimental structure, gives a strong bias toward the consensus binding site.

Algorithms↗

An atomic resolution structure for human fibroblast growth factor 1.

A 1.10-A atomic resolution X-ray structure of human fibroblast growth factor 1 (FGF-1), a member of the beta-trefoil superfold, has been determined. The beta-trefoil is one of 10 fundamental protein superfolds and is the only superfold to exhibit 3-fold structural symmetry (comprising 3 "trefoil" units). The quality of the diffraction data permits unambiguous assignment of Asn, Gln, and His rotamers, Pro ring pucker, as well as refinement of atomic anisotropic displacement parameters (ADPs). The FGF-1 structure exhibits numerous core-packing defects, detectable using a 1.0-A probe radius. In addition to contributing to the relatively low thermal stability of FGF-1, these defects may also permit domain motions within the structure. The availability of refined ADPs allows a translation/libration/screw (TLS) analysis of putative rigid body domains. The TLS analysis shows that beta-strands 6-12 together form a rigid body, and there is a clear demarcation in TLS motions between the adjacent carboxyl- and amino-termini. Although separate from beta-strands 6-12, the individual beta-strands 1-5 do not exhibit correlated motions; thus, this region appears to be comparatively flexible. The heparin-binding contacts of FGF-1 are located within beta-strands 6-12; conversely, a significant portion of the receptor-binding contacts are located within beta-strands 1-5. Thus, the observed rigid body motion in FGF-1 appears related to the ligand-binding functionalities.

Anisotropy↗

Crystallographic studies of V44 mutants of Clostridium pasteurianum rubredoxin: effects of side-chain size on reduction potential.

Understanding the structural origins of differences in reduction potentials is crucial to understanding how various electron transfer proteins modulate their reduction potentials and how they evolve for diverse functional roles. Here, the high-resolution structures of several Clostridium pasteurianum rubredoxin (Cp Rd) variants with changes in the vicinity of the redox site are reported in order to increase this understanding. Our crystal structures of [V44L] (at 1.8 A resolution), [V44A] (1.6 A), [V44G] (2.0 A) and [V44A, G45P] (1.5 A) Rd (all in their oxidized states) show that there is a gradual decrease in the distance between Fe and the amide nitrogen of residue 44 upon reduction in the size of the side chain of residue 44; the decrease occurs from leucine to valine, alanine or glycine and is accompanied by a gradual increase in their reduction potentials. Mutation of Cp Rd at position 44 also changes the hydrogen-bond distance between the amide nitrogen of residue 44 and the sulfur of cysteine 42 in a size-dependent manner. Our results suggest that residue 44 is an important determinant of Rd reduction potential in a manner dictated by side-chain size. Along with the electric dipole moment of the 43-44 peptide bond and the 44-42 NH--S type hydrogen bond, a modulation mechanism for solvent accessibility through residue 41 might regulate the redox reaction of the Rds.

Bacterial Proteins↗

The impact of protein flexibility on protein-protein docking.

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.

Computer Simulation↗

Flexibility of prolyl oligopeptidase: molecular dynamics and molecular framework analysis of the potential substrate pathways.

The flexibility of prolyl oligopeptidase has been investigated using molecular dynamics (MD) and molecular framework approaches to delineate the route of the substrate to the active site. The selectivity of the enzyme is mediated by a seven-bladed beta-propeller that in the crystal structure does not indicate the possible passage for the substrate to the catalytic center. Its open topology however, could allow the blades to move apart and let the substrate into the large central cavity. Flexibility analysis of prolyl oligopeptidase structure using the FIRST (Floppy Inclusion and Rigid Substructure Topology) approach and the atomic fluctuations derived from MD simulations demonstrated the rigidity of the propeller domain, which does not permit the substrate to approach the active site through this domain. Instead, a smaller tunnel at the inter-domain region comprising the highly flexible N-terminal segment of the peptidase domain and a facing hydrophilic loop from the propeller (residues 192-205) was identified by cross-correlation analysis and essential dynamics as the only potential pathway for the substrate. The functional importance of the flexible loop has been also verified by kinetic analysis of the enzyme with a split loop. Catalytic effect of engineered disulfide bridges was rationalized by characterizing the concerted motions of the two domains.

Animals↗

Protein flexibility and rigidity predicted from sequence.

Structural flexibility has been associated with various biological processes such as molecular recognition and catalytic activity. In silico studies of protein flexibility have attempted to characterize and predict flexible regions based on simple principles. B-values derived from experimental data are widely used to measure residue flexibility. Here, we present the most comprehensive large-scale analysis of B-values. We used this analysis to develop a neural network-based method that predicts flexible-rigid residues from amino acid sequence. The system uses both global and local information (i.e., features from the entire protein such as secondary structure composition, protein length, and fraction of surface residues, and features from a local window of sequence-consecutive residues). The most important local feature was the evolutionary exchange profile reflecting sequence conservation in a family of related proteins. To illustrate its potential, we applied our method to 4 different case studies, each of which related our predictions to aspects of function. The first 2 were the prediction of regions that undergo conformational switches upon environmental changes (switch II region in Ras) and the prediction of surface regions, the rigidity of which is crucial for their function (tunnel in propeller folds). Both were correctly captured by our method. The third study established that residues in active sites of enzymes are predicted by our method to have unexpectedly low B-values. The final study demonstrated how well our predictions correlated with NMR order parameters to reflect motion. Our method had not been set up to address any of the tasks in those 4 case studies. Therefore, we expect that this method will assist in many attempts at inferring aspects of function.

Amino Acid Sequence↗

Enhanced thermal stability achieved without increased conformational rigidity at physiological temperatures: spatial propagation of differential flexibility in rubredoxin hybrids.

The extreme thermal stability of proteins from hyperthermophilic organisms is widely believed to arise from an increased conformational rigidity in the native state. In apparent contrast to this paradigm, both Pyrococcus furiosus (Pf) rubredoxin, the most thermostable protein characterized to date, and its Clostridium pasteurianum (Cp) mesophile homolog undergo a transient conformational opening of their multi-turn segments, which is more favorable in hyperthermophile proteins below room temperature. Substitution of the hyperthermophile multi-turn sequence into the mesophile protein sequence yields a hybrid, (14-33(Pf)) Cp, that exhibits a 12 degrees increase in its reversible thermal unfolding transition midpoint. Nuclear magnetic resonance (NMR) magnetization transfer-based hydrogen exchange was used to monitor backbone conformational dynamics in the subsecond time regime. Despite the substantially increased thermostability, flexibility throughout the entire main chain of the more thermostable hybrid is equal to or greater than that of the wild type mesophile rubredoxin near its normal growth temperature. In comparison to the identical core residues of the (14-33(Pf)) Cp rubredoxin hybrid, six spatially clustered residues in the parental mesophile protein exhibit a substantially larger temperature dependence of exchange. The exchange behavior of these six residues closely matches that observed in the multi-turn segment, consistent with a more extensive conformational process. These six core residues exhibit a much weaker temperature dependence of exchange in the (14-33(Pf)) Cp hybrid, similar to that observed for the multi-turn segment in its parental Pf rubredoxin. These results suggest that differential temperature dependence of flexibility can underlie variations in thermostability observed for mesophile versus hyperthermophile homologs.

Amides↗

Grouping of amino acid types and extraction of amino acid properties from multiple sequence alignments using variance maximization.

Understanding of amino acid type co-occurrence in trusted multiple sequence alignments is a prerequisite for improved sequence alignment and remote homology detection algorithms. Two objective approaches were used to investigate co-occurrence, both based on variance maximization of the weighted residue frequencies in columns taken from a large alignment database. The first approach discretely grouped amino acid types, and the second approach extracted orthogonal properties of amino acids using principal components analysis. The grouping results corresponded to amino acid physical properties such as side chain hydrophobicity, size, or backbone flexibility, and an optimal arrangement of approximately eight groups was observed. However, interpretation of the orthogonal properties was more complex. Although the principal components accounting for the largest variances exhibited modest correlations with hydrophobicity and conservation of glycine, in general principal components did not correspond to physical properties of amino acids. Although not intuitive, these amino acid mathematical properties were demonstrated to be robust and to improve local pairwise alignment accuracy, relative to 20 amino acid frequencies alone, for a simple test case.

Amino Acids↗

Conformational flexibility may explain multiple cellular roles of PEST motifs.

PEST sequences are one of the major motifs that serve as signal for the protein degradation and are also involved in various cellular processes such as phosphorylation and protein-protein interaction. In our earlier study, we found that these motifs contribute largely to eukaryotic protein disorder. This observation led us to evaluate their conformational variability in the nonredundant Protein Data Bank (PDB) structures. For this purpose, crystallographic temperature factors, structural alignment of multiple NMR models, and dihedral angle order parameters have been used in this study. The study has revealed the hypermobility of PEST motifs as compared to other regions of the protein. Conformational flexibility may allow them to participate in number of molecular interactions under different conditions. This analysis may explain the role of protein backbone flexibility in bringing about multiple cellular roles of PEST motifs.

Amino Acid Motifs↗

Effective handling of induced-fit motion in flexible docking.

For structure-based drug design, where various ligand structures need to be docked to a target protein structure, a docking method that can handle conformational flexibility of not only the ligand, but also the protein, is indispensable. We have developed a simple and effective approach for dealing with the local induced-fit motion of the target protein, and implemented it in our docking tool, ADAM. Our approach efficiently combines the following two strategies: a vdW-offset grid in which the protein cavity is enlarged uniformly, and structure optimization allowing the motion of ligand and protein atoms. To examine the effectiveness of our approach, we performed docking validation studies, including redocking in 18 test cases and foreign-docking, in which various ligands from foreign crystal structures of complexes are docked into a target protein structure, in 22 cases (on five target proteins). With the original ADAM, the correct docking modes (RMSD < 2.0 A) were not present among the top 20 models in one case of redocking and four cases of foreign-docking. When the handling of induced-fit motion was implemented, the correct solutions were acquired in all 40 test cases. In foreign-docking on thymidine kinase, the correct docking modes were obtained as the top-ranked solutions for all 10 test ligands by our combinatorial approach, and this appears to be the best result ever reported with any docking tool. The results of docking validation have thus confirmed the effectiveness of our approach, which can provide reliable docking models even in the case of foreign-docking, where conformational change of the target protein cannot be ignored. We expect that this approach will contribute substantially to actual drug design, including virtual screening.

Ligands↗

Conformational flexibility of RecA protein filament: transitions between compressed and stretched states.

RecA protein is a central enzyme in homologous DNA recombination, repair and other forms of DNA metabolism in bacteria. It functions as a flexible helix-shaped filament bound on stretched single-stranded or double-stranded DNA in the presence of ATP. In this work, we present an atomic level model for conformational transitions of the RecA filament. The model describes small movements of the RecA N-terminal domain due to coordinated rotation of main chain dihedral angles of two amino acid residues (Psi/Lys23 and Phi/Gly24), while maintaining unchanged the RecA intersubunit interface. The model is able to reproduce a wide range of observed helix pitches in transitions between compressed and stretched conformations of the RecA filament. Predictions of the model are in agreement with Small Angle Neutron Scattering (SANS) measurements of the filament helix pitch in RecA::ADP-AlF(4) complex at various salt concentrations.

Amino Acids↗

Kinematics and the implementation of an elephant's trunk manipulator and other continuum style robots.

Traditionally, robot manipulators have been a simple arrangement of a small number of serially connected links and actuated joints. Though these manipulators prove to be very effective for many tasks, they are not without their limitations, due mainly to their lack of maneuverability or total degrees of freedom. Continuum style (i.e., continuous "back-bone") robots, on the other hand, exhibit a wide range of maneuverability, and can have a large number of degrees of freedom. The motion of continuum style robots is generated through the bending of the robot over a given section; unlike traditional robots where the motion occurs in discrete locations, i.e., joints. The motion of continuum manipulators is often compared to that of biological manipulators such as trunks and tentacles. These continuum style robots can achieve motions that could only be obtainable by a conventionally designed robot with many more degrees of freedom. In this paper we present a detailed formulation and explanation of a novel kinematic model for continuum style robots. The design, construction, and implementation of our continuum style robot called the elephant trunk manipulator is presented. Experimental results are then provided to verify the legitimacy of our model when applied to our physical manipulator. We also provide a set of obstacle avoidance experiments that help to exhibit the practical implementation of both our manipulator and our kinematic model.

Animals↗

Polymer surface science.

Molecular level studies of the structure and mechanical properties of polymer surfaces have been carried out by sum frequency generation (SFG) surface vibrational spectroscopy and atomic force microscopy (AFM). The surfaces of different grades of polyethylene and polypropylene have been characterized-including during the glass transition and when mechanically stretched. Copolymers that have hard and soft segments with different glass transition temperatures show phase separation, an effect of hydrogen bonding between the hard and soft segments, that influences their adhesive and friction properties. AFM and SFG show that low surface energy additives migrate to the surface and alter the surface mechanical properties. Polymers, where the chemical nature of the end groups is different from the backbone, show surface segregation of the hydrophobic part of the chain in air and the hydrophilic part in water. Likewise, in miscible polymer blends, surface segregation of the more hydrophobic component in air and the more hydrophilic component in water is observed. This area of surface science requires increased attention because of the predominance of polymers as structural materials and as biomaterials.

Biocompatible Materials↗

Molecular dynamics simulation and nuclear magnetic resonance studies of the terminal glucotriose unit found in the oligosaccharide of glycoprotein precursors.

The trisaccharide alpha-d-Glcp-(1 --> 2)-alpha-d-Glcp-(1 --> 3)-alpha-d-Glcp-OMe, a model for the terminal glucotriose in Glc(3)Man(9)GlcNAc(2) in glycoprotein precursors, has been investigated by computer simulations and NMR spectroscopy. Molecular dynamics simulations were performed for 1 ns in aqueous solution and 20 ns in vacuo using the CHARMM-based force fields PARM22 and CHEAT95. An additional Monte Carlo simulation with the HSEA force field was also carried out. Experimental NMR data in water solution was obtained from measurement of long-range (1)H,(13)C heteronuclear trans-glycosidic coupling constants, (3)J(H,C), using one-dimensional Hadamard spectroscopy. Calculation of the (3)J(H,C) values from the simulations showed a varying degree of agreement to experimental data. It could be shown from simulation that the &phi; torsion angles differed, which was corroborated by the NMR measurements. Analyses were done of radial distribution functions and of hydrogen bonds. It was suggested that intermolecular hydrogen bonds were present, but in contrast to simulation the results from NMR spectroscopy did not support any major contribution. Hence, their influence on the conformation of the trisaccharide is rather small. Comparison of (1)H NMR chemical shifts for the trisaccharide and the glucotriose in Glc(3)Man(8)GlcNAc revealed high similarity. However, the derived conformation of the model substance in this work differed at one glycosidic torsion angle compared to the glucotriose on a large oligosaccharide.

Carbohydrate Conformation↗

Acquisition of ordered conformation by the N-terminal domain of the human small proline rich 2 protein.

The cornified cell envelope (CE) is a crucial structure for barrier function in terminally differentiated dead stratified squamous epithelia. It is assembled by transglutaminase enzymes (TGases) that cross-link several proteins such as loricrin and the small proline rich (SPR) proteins. Human SPR2 protein is cross-linked with widely differing efficiencies by TGases 1, 2, and 3 using exclusively residues in the N- and C-terminal domains. In order to understand if the absence of the cross-linking catalyzed by TGases in the central domain is due to the conformation adopted, we have investigated the structural properties in solution of three peptides that correspond to the N-terminal domain, to three repeats of the central domain, and to the C-terminal domain. Together, the NMR and CD data strongly indicate the presence of a highly flexible non alpha-helix, non beta-sheet structure in SPR2. Thus, SPR2 appears to function as a flexible cross-bridging protein to provide tensile strength or rigidity to the CE of the stratified squamous epithelia in which it is expressed.

Amino Acid Sequence↗