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 343 records · Page 19Linked to original sources
PubMed1996

Flexibility of DNA binding domain of trp repressor required for recognition of different operator sequences.

Trp repressor (25 kDa) is a regulatory protein that controls transcription initiation in the tryptophan biosynthetic operon and at least four other operons in Escherichia coli. An alanine to valine mutation (AV77) in the DNA binding domain is known to increase repressor activity at the trp operator in vivo, but not in vitro. We report here the amide proton exchange rates for the DNA-binding domains of both the wild-type and AV77 proteins. We find that the alanine to valine change stabilizes the flexible DNA-binding domain of the repressor. We present in vivo data showing that, although the AV77 repressor is more inhibitory at the trp operator than the wild-type repressor, it does not have increased activity at the aroH or trpR operator; repression at the aroH operator is, in fact, reduced. Our results suggest that the flexibility exhibited by the wild-type repressor allows a broader range of repressor/DNA interactions, whereas the increased rigidity resulting from the AV77 change limits the repressor's effectiveness at some operators.

Amides↗
PubMed1997

Secondary and tertiary structural changes in gamma delta resolvase: comparison of the wild-type enzyme, the I110R mutant, and the C-terminal DNA binding domain in solution.

gamma delta Resolvase is a site-specific DNA recombinase (M(r) 20.5 kDa) in Escherichia coli that shares homology with a family of bacterial resolvases and invertases. We have characterized the secondary and tertiary structural behavior of the cloned DNA binding domain (DBD) and a dimerization defective mutant in solution. Low-salt conditions were found to destabilize the tertiary structure of the DBD dramatically, with concomitant changes in the secondary structure that were localized near the hinge regions between the helices. The molten tertiary fold appears to contribute significantly to productive DNA interactions and supports a mechanism of DNA-induced folding of the tertiary structure, a process that enables the DBD to adapt in conformation for each of the three imperfect palindromic sites. At high salt concentrations, the monomeric I110R resolvase shows a minimal perturbation to the three helices of the DBD structure and changes in the linker segment in comparison to the cloned DBD containing the linker. Comparative analysis of the NMR spectra suggest that the I110R mutant contains a folded catalytic core of approximately 60 residues and that the segment from residues 100 to 149 are devoid of regular structure in the I110R resolvase. No increase in the helicity of the linker region of I110R resolvase occurs on binding DNA. These results support a subunit rotation model of strand exchange that involves the partial unfolding of the catalytic domains.

Ammonium Sulfate↗
PubMed2002

Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase.

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.

Bacterial Proteins↗
PubMed2002

Automated docking to multiple target structures: incorporation of protein mobility and structural water heterogeneity in AutoDock.

Protein motion and heterogeneity of structural waters are approximated in ligand-docking simulations, using an ensemble of protein structures. Four methods of combining multiple target structures within a single grid-based lookup table of interaction energies are tested. The method is evaluated using complexes of 21 peptidomimetic inhibitors with human immunodeficiency virus type 1 (HIV-1) protease. Several of these structures show motion of an arginine residue, which is essential for binding of large inhibitors. A structural water is also present in 20 of the structures, but it must be absent in the remaining one for proper binding. Mean and minimum methods perform poorly, but two weighted average methods permit consistent and accurate ligand docking, using a single grid representation of the target protein structures.

Computer Simulation↗
PubMed2003

Backbone dynamics of the CC-chemokine eotaxin-2 and comparison among the eotaxin group chemokines.

The eotaxin group chemokines (eotaxin, eotaxin-2, and eotaxin-3) share only 35-41% sequence identity but are all agonists for the receptor CCR3. Here we present a detailed comparison between the backbone dynamics of these three chemokines. The dynamics of eotaxin-2 were determined from 15N NMR relaxation data and compared to those obtained previously for eotaxin and eotaxin-3. For all three chemokines, the majority of residues in the first two beta-strands and the alpha-helix show highly restricted motions on the subnanosecond time scale but there is dramatically higher flexibility in the N- and C-terminal regions and also substantial mobility for residues in the N-loop region and the third beta-strand. The latter two regions form a groove on the chemokine surface that is the likely binding site for the N-terminal region of the receptor. Taken together, the available data suggest a model in which conformational rearrangements of both the chemokine and the receptor are likely to occur during binding and receptor activation.

Binding Sites↗
PubMed2003

Interaction of E. coli outer-membrane protein A with sugars on the receptors of the brain microvascular endothelial cells.

Esherichia coli, the most common gram-negative bacteria, can penetrate the brain microvascular endothelial cells (BMECs) during the neonatal period to cause meningitis with significant morbidity and mortality. Experimental studies have shown that outer-membrane protein A (OmpA) of E. coli plays a key role in the initial steps of the invasion process by binding to specific sugar moieties present on the glycoproteins of BMEC. These experiments also show that polymers of chitobiose (GlcNAcbeta1-4GlcNAc) block the invasion, while epitopes substituted with the L-fucosyl group do not. We used HierDock computational technique that consists of a hierarchy of coarse grain docking method with molecular dynamics (MD) to predict the binding sites and energies of interactions of GlcNAcbeta1-4GlcNAc and other sugars with OmpA. The results suggest two important binding sites for the interaction of carbohydrate epitopes of BMEC glycoproteins to OmpA. We identify one site as the binding pocket for chitobiose (GlcNAcbeta1-4GlcNAc) in OmpA, while the second region (including loops 1 and 2) may be important for recognition of specific sugars. We find that the site involving loops 1 and 2 has relative binding energies that correlate well with experimental observations. This theoretical study elucidates the interaction sites of chitobiose with OmpA and the binding site predictions made in this article are testable either by mutation studies or invasion assays. These results can be further extended in suggesting possible peptide antagonists and drug design for therapeutic strategies.

Bacterial Outer Membrane Proteins↗
PubMed2003

Dynamic effects of mutations within two loops of cytochrome c551 from Pseudomonas aeruginosa.

In this work, we investigated the structural and dynamic consequences of two substitutions, P58A and G36P, located in two different solvent-exposed loops of cytochrome c551. The results show that both mutations affect regions that are distant from the site of mutation. Here, the two loops appear to be dynamically coupled to each other, because the substitution at one site affects the structure and the dynamics of the other site. However, the substitutions at Gly-36 and Pro-58 presented substantial differences, which were related to the mechanical (rigidity and deformability) properties of the site surrounding the mutation. Although the P58A mutant conserved a significant dynamic similarity to the wild-type protein as the immediate surroundings of position 58 became more rigid, the G36P mutant, which had deformed its flexible surroundings, presented a dynamic behavior that was markedly different from that of the wild-type protein. These results suggest that perturbation of sterically isolated and flexible regions, such as solvent-exposed loops, can have strong dynamic consequences on the protein as a whole, raising the possibility that these effects could in turn affect the stability or the function of the protein.

Amino Acid Substitution↗
PubMed2003

Protein side-chain rearrangement in regions of point mutations.

A major problem in predicting amino acid side-chain rearrangements following point mutations is the potentially large search space. We analyzed a nonredundant data set of 393 Protein Data Bank protein pairs, each consisting of structures differing in one amino acid, to determine the number of residues changing conformation in the region of mutation. In 91-95% of cases, two or fewer residues underwent side-chain conformational change. If mutation sites with backbone displacements were excluded, the number increased to 97%. The majority of rearrangements (over 60%) were due to the inherent flexibility of side-chains, as derived from analysis of a control set of protein subunits whose crystal structures were determined more than once. Different amino acids demonstrated different degrees of flexibility near mutation sites. Large polar or charged residues, and serine, are more flexible, while the aromatic amino acids, and cysteine, are less so. This pattern is common to the inherent side-chain flexibility, as well as the increased flexibility at ligand binding sites and mutation sites. The probability for conformational change was correlated with B-factor, frequency of the side-chain conformation in proteins and solvent accessibility. The last trend was stronger for aromatic and hydrophilic residues than for hydrophobic ones. We conclude that the search space for predicting side-chain conformations in the region of mutation can be effectively restricted. However, the overall ability to predict a particular side-chain conformation, or to check predictions according to individual existing structures, is limited. These findings may be useful in deriving empirical rules for modeling side-chain conformations.

Binding Sites↗
PubMed2003

Insights into the bile acid transportation system: the human ileal lipid-binding protein-cholyltaurine complex and its comparison with homologous structures.

Bile acids are generated in vivo from cholesterol in the liver, and they undergo an enterohepatic circulation involving the small intestine, liver, and kidney. To understand the molecular mechanism of this transportation, it is essential to gain insight into the three-dimensional (3D) structures of proteins involved in the bile acid recycling in free and complexed form and to compare them with homologous members of this protein family. Here we report the solution structure of the human ileal lipid-binding protein (ILBP) in free form and in complex with cholyltaurine. Both structures are compared with a previously published structure of the porcine ILBP-cholylglycine complex and with related lipid-binding proteins. Protein structures were determined in solution by using two-dimensional (2D)- and 3D-homo and heteronuclear NMR techniques, leading to an almost complete resonance assignment and a significant number of distance constraints for distance geometry and restrained molecular dynamics simulations. The identification of several intermolecular distance constraints unambiguously determines the cholyltaurine-binding site. The bile acid is deeply buried within ILBP with its flexible side-chain situated close to the fatty acid portal as entry region into the inner ILBP core. This binding mode differs significantly from the orientation of cholylglycine in porcine ILBP. A detailed analysis using the GRID/CPCA strategy reveals differences in favorable interactions between protein-binding sites and potential ligands. This characterization will allow for the rational design of potential inhibitors for this relevant system.

Animals↗
PubMed2004

Myosin flexibility: structural domains and collective vibrations.

The movement of the myosin motor along an actin filament involves a directed conformational change within the cross-bridge formed between the protein and the filament. Despite the structural data that has been obtained on this system, little is known of the mechanics of this conformational change. We have used existing crystallographic structures of three conformations of the myosin head, containing the motor domain and the lever arm, for structural comparisons and mechanical studies with a coarse-grained elastic network model. The results enable us to define structurally conserved domains within the protein and to better understand myosin flexibility. Notably they point to the role of the light chains in rigidifying the lever arm and to changes in flexibility as a consequence of nucleotide binding.

Animals↗
PubMed2004

Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP.

Most current docking methods to identify possible ligands and putative binding sites on a receptor molecule assume a rigid receptor structure to allow virtual screening of large ligand databases. However, binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a bound ligand. An approach is presented that allows relaxation of the protein conformation in precalculated soft flexible degrees of freedom during ligand-receptor docking. For the immunosuppressant FK506-binding protein FKBP, the soft flexible modes are extracted as principal components of motion from a molecular dynamics simulation. A simple penalty function for deformations in the soft flexible mode is used to limit receptor protein deformations during docking that avoids a costly recalculation of the receptor energy by summing over all receptor atom pairs at each step. Rigid docking of the FK506 ligand binding to an unbound FKBP conformation failed to identify a geometry close to experiment as favorable binding site. In contrast, inclusion of the flexible soft modes during systematic docking runs selected a binding geometry close to experiment as lowest energy conformation. This has been achieved at a modest increase of computational cost compared to rigid docking. The approach could provide a computationally efficient way to approximately account for receptor flexibility during docking of large numbers of putative ligands and putative docking geometries.

Binding Sites↗
PubMed2001

Differential effects of isomeric incorporation of fluorophenylalanines into PvuII endonuclease.

Incorporation of fluorine into proteins has long served as a means of probing structure and function, yet there are few studies that examine the impact of fluorine substitution, particularly at locations distant from the active sites of enzymes. The flexibility of isomeric fluorine incorporation at Phe is used to explore subtle substitution effects on enzyme activity and conformation. The unnatural amino acids o-, m-, and p-fluorophenylalanines were incorporated biosynthetically into the representative PvuII restriction endonuclease. Interestingly, m-fluoro-Phe-PvuII endonuclease exhibits very similar conformational stability to that of the native enzyme, but it exhibits a reproducible, 2-fold higher average specific activity. Given the level of incorporation and the distribution of species, the species of modified enzyme responsible for this increase in specific activity is most likely even faster. Further, moving the fluorine atom from the meta- to the para-position of Phe results in a 4-fold decrease in specific activity and a decrease in conformational stability of 1.5 kcal/mol. Since none of the Phe residues in PvuII endonuclease lies near the DNA recognition or catalytic sites, this differential behavior alludes to the impact of subtle changes in enzyme conformation on endonuclease activity and suggests novel ways to influence catalytic behavior.

Allosteric Site↗
PubMed2004

Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations.

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.

Acetates↗
PubMed2004

Locating the stabilizing residues in (alpha/beta)8 barrel proteins based on hydrophobicity, long-range interactions, and sequence conservation.

In nature, 1 out of every 10 proteins has an (alpha/beta)(8) (TIM)-barrel fold, and in most cases, pairwise comparisons show no sequence similarity between them. Hence, delineating the key residues that induce very different sequences to share a common fold is important for understanding the folding and stability of TIM-barrel domains. In this work, we propose a new consensus approach for locating these stabilizing residues based on long-range interactions, hydrophobicity, and conservation of amino acid residues. We have identified 957 stabilizing residues in 63 proteins from a nonredundant set of 71 TIM-barrel domains. Most of these residues are located in the 8-stranded beta-sheet, with nearly one half of them oriented toward the interior of the barrel and the other half oriented toward the surrounding alpha-helices. Several stabilizing residues are found in the N- and C-terminal loops, whereas very few appear in the alpha-helices that surround the internal beta-sheet. Further, these 957 residues are placed in 434 stabilizing segments of various sizes, and each domain contains 1-10 of these segments. We found that 8 segments per domain is the most abundant one, and two thirds of the proteins have 7-9 stabilizing segments. Finally, we verified the identified residues with experimental temperature factors and found that these residues are among the ones with less mobility in the considered proteins. We suggest that our new protocol serves as a powerful tool to identify the stabilizing residues in TIM-barrel domains, which can be used as potential candidates for studying protein folding and stability by means of protein engineering experiments.

Conserved Sequence↗
PubMed2004

Characterization of Lys-698-to-Met substitution in human plasminogen catalytic domain.

Streptokinase (SK) is a human plasminogen (Pg) activator secreted by streptococci. The activation mechanism of SK differs from that of physiological Pg activators in that SK is not a protease and cannot proteolytically activate Pg. Instead, it forms a tight complex with Pg that proteolytically activates other Pg molecules. The residue Lys-698 of human Pg was hypothesized to participate in triggering activation in the SK-Pg complex. Here, we report a study of the Lys-698 to Met substitution in the catalytic domain of Pg (microPg) containing the proteolytic activation-resistant background (R561A). While it remains competent in forming a complex with SK, maintaining a comparable equilibration dissociation constant (K(D)), the recombinant protein shows a nearly 60-fold reduction in amidolytic activity relative to its R561A background when mixed with native SK. A 2.3 A crystal structure of this mutant microPg confirmed the correct folding of this recombinant protein. Combined with other biochemical data, these results support the premise that Lys-698 of human Pg plays a functional role in the so-called N-terminal insertion activation mechanism by SK.

Amino Acid Substitution↗
PubMed2004

A detailed comparison of current docking and scoring methods on systems of pharmaceutical relevance.

A thorough evaluation of some of the most advanced docking and scoring methods currently available is described, and guidelines for the choice of an appropriate protocol for docking and virtual screening are defined. The generation of a large and highly curated test set of pharmaceutically relevant protein-ligand complexes with known binding affinities is described, and three highly regarded docking programs (Glide, GOLD, and ICM) are evaluated on the same set with respect to their ability to reproduce crystallographic binding orientations. Glide correctly identified the crystallographic pose within 2.0 A in 61% of the cases, versus 48% for GOLD and 45% for ICM. In general Glide appears to perform most consistently with respect to diversity of binding sites and ligand flexibility, while the performance of ICM and GOLD is more binding site-dependent and it is significantly poorer when binding is predominantly driven by hydrophobic interactions. The results also show that energy minimization and reranking of the top N poses can be an effective means to overcome some of the limitations of a given docking function. The same docking programs are evaluated in conjunction with three different scoring functions for their ability to discriminate actives from inactives in virtual screening. The evaluation, performed on three different systems (HIV-1 protease, IMPDH, and p38 MAP kinase), confirms that the relative performance of different docking and scoring methods is to some extent binding site-dependent. GlideScore appears to be an effective scoring function for database screening, with consistent performance across several types of binding sites, while ChemScore appears to be most useful in sterically demanding sites since it is more forgiving of repulsive interactions. Energy minimization of docked poses can significantly improve the enrichments in systems with sterically demanding binding sites. Overall Glide appears to be a safe general choice for docking, while the choice of the best scoring tool remains to a larger extent system-dependent and should be evaluated on a case-by-case basis.

Algorithms↗
PubMed2004

Prediction of transmembrane regions of beta-barrel proteins using ANN- and SVM-based methods.

This article describes a method developed for predicting transmembrane beta-barrel regions in membrane proteins using machine learning techniques: artificial neural network (ANN) and support vector machine (SVM). The ANN used in this study is a feed-forward neural network with a standard back-propagation training algorithm. The accuracy of the ANN-based method improved significantly, from 70.4% to 80.5%, when evolutionary information was added to a single sequence as a multiple sequence alignment obtained from PSI-BLAST. We have also developed an SVM-based method using a primary sequence as input and achieved an accuracy of 77.4%. The SVM model was modified by adding 36 physicochemical parameters to the amino acid sequence information. Finally, ANN- and SVM-based methods were combined to utilize the full potential of both techniques. The accuracy and Matthews correlation coefficient (MCC) value of SVM, ANN, and combined method are 78.5%, 80.5%, and 81.8%, and 0.55, 0.63, and 0.64, respectively. These methods were trained and tested on a nonredundant data set of 16 proteins, and performance was evaluated using "leave one out cross-validation" (LOOCV). Based on this study, we have developed a Web server, TBBPred, for predicting transmembrane beta-barrel regions in proteins (available at http://www.imtech.res.in/raghava/tbbpred).

Algorithms↗
PubMed2004

The structure of murine neuroglobin: Novel pathways for ligand migration and binding.

Neuroglobin, a recently discovered globin predominantly expressed in neuronal tissue of vertebrates, binds small, gaseous ligands at the sixth coordination position of the heme iron. In the absence of an exogenous ligand, the distal histidine (His64) binds to the heme iron in the ferrous and ferric states. The crystal structure of murine ferric (met) neuroglobin at 1.5 A reveals interesting features relevant to the ligand binding mechanism. Only weak selectivity is observed for the two possible heme orientations, the occupancy ratio being 70:30. Two small internal cavities are present on the heme distal side, which enable the His64(E7) side chain to move out of the way upon exogenous ligand binding. Moreover, a third, huge cavity (volume approximately 290 A3) connecting both sides of the heme, is open towards the exterior and provides a potential passageway for ligands. The CD and EF corners exhibit substantial flexibility, which may assist ligands in entering the protein and accessing the active site. Based on this high-resolution structure, further structure-function studies can be planned to elucidate the role of neuroglobin in physiological responses to hypoxia.

Animals↗