Slow conformational dynamics in the hamster prion protein

Although the mechanism of the conformational conversion from the cellular (PrP(C)) to the scrapie (PrP(Sc)) form of animal prion proteins has yet to be elucidated, evidence is accumulating that may provide insight into the conversion process at atomic resolution. Here we show critical aspects of the slow fluctuation dynamics of the recombinant hamster prion protein, rPrP(90-231), based on NMR relaxation analysis using Carr-Purcell-Meiboom-Gill (CPMG) experiments, and compare them in detail with results from high-pressure NMR. Residues exhibiting slow fluctuations on the time scale of microseconds to milliseconds are mainly localized on helices B and C (172-193 and 200-227), which include locally disordered regions in an intermediate conformer, PrP*, identified previously by high-pressure NMR [Kuwata, K., et al., (2002) Biochemistry 41, 12277-12283]. Moreover, chemical shift differences between two putative exchanging conformers obtained by the CPMG relaxation analysis and the linear component of the pressure-induced chemical shift changes are reasonably correlated at individual residue sites. These observations suggest that both the CMPG relaxation and the pressure shifts reflect slow conformational fluctuations and that these slow motions in PrP(C) are related to the trajectories leading to the transition to PrP*.     

CLEMAPS: multiple alignment of protein structures based on conformational letters

CLEMAPS is a tool for multiple alignment of protein structures. It distinguishes itself from other existing algorithms for multiple structure alignment by the use of conformational letters, which are discretized states of 3D segmental structural states. A letter corresponds to a cluster of combinations of three angles formed by C(alpha) pseudobonds of four contiguous residues. A substitution matrix called CLESUM is available to measure the similarity between any two such letters. The input 3D structures are first converted to sequences of conformational letters. Each string of a fixed length is then taken as the center seed to search other sequences for neighbors of the seed, which are strings similar to the seed. A seed and its neighbors form a center-star, which corresponds to a fragment set of local structural similarity shared by many proteins. The detection of center-stars using CLESUM is extremely efficient. Local similarity is a necessary, but insufficient, condition for structural alignment. Once center-stars are found, the spatial consistency between any two stars are examined to find consistent star duads using atomic coordinates. Consistent duads are later joined to create a core for multiple alignment, which is further polished to produce the final alignment. The utility of CLEMAPS is tested on various protein structure ensembles.     

Protein under pressure: molecular dynamics simulation of the arc repressor

Experimental nuclear magnetic resonance results for the Arc Repressor have shown that this dimeric protein dissociates into a molten globule at high pressure. This structural change is accompanied by a modification of the hydrogen-bonding pattern of the intermolecular beta-sheet: it changes its character from intermolecular to intramolecular with respect to the two monomers. Molecular dynamics simulations of the Arc Repressor, as a monomer and a dimer, at elevated pressure have been performed with the aim to study this hypothesis and to identify the major structural and dynamical changes of the protein under such conditions. The monomer appears less stable than the dimer. However, the complete dissociation has not been seen because of the long timescale needed to observe this phenomenon. In fact, the protein structure altered very little when increasing the pressure. It became slightly compressed and the dynamics of the side-chains and the unfolding process slowed down. Increasing both, temperature and pressure, a tendency of conversion of intermolecular into intramolecular hydrogen bonds in the beta-sheet region has been detected, supporting the mentioned hypothesis. Also, the onset of denaturation of the separated chains was observed.     

Slow conformational dynamics of the guanine nucleotide-binding protein Ras complexed with the GTP analogue GTPgammaS

The guanine nucleotide-binding protein Ras occurs in solution in two different conformational states, state 1 and state 2 with an equilibrium constant K(12) of 2.0, when the GTP analogue guanosine-5'-(beta,gamma-imido)triphosphate or guanosine-5'-(beta,gamma-methyleno)triphosphate is bound to the active centre. State 2 is assumed to represent a strong binding state for effectors with a conformation similar to that found for Ras complexed to effectors. In the other state (state 1), the switch regions of Ras are most probably dynamically disordered. Ras variants that exist predominantly in state 1 show a drastically reduced affinity to effectors. In contrast, Ras(wt) bound to the GTP analogue guanosine-5'-O-(3-thiotriphosphate) (GTPgammaS) leads to (31)P NMR spectra that indicate the prevalence of only one conformational state with K(12) > 10. Titration with the Ras-binding domain of Raf-kinase (Raf-RBD) shows that this state corresponds to effector binding state 2. In the GTPgammaS complex of the effector loop mutants Ras(T35S) and Ras(T35A) two conformational states different to state 2 are detected, which interconvert over a millisecond time scale. Binding studies with Raf-RBD suggest that both mutants exist mainly in low-affinity states 1a and 1b. From line-shape analysis of the spectra measured at various temperatures an activation energy DeltaH(|) (1a1b) of 61 kJ.mol(-1) and an activation entropy DeltaS(|) (1a1b) of 65 J.K(-1).mol(-1) are derived. Isothermal titration calorimetry on Ras bound to the different GTP-analogues shows that the effective affinity K(A) for the Raf-RBD to Ras(T35S) is reduced by a factor of about 20 compared to the wild-type with the strongest reduction observed for the GTPgammaS complex.     

Prediction of protein secondary structures from conformational biases

We use LINUS (the "Local Independently Nucleated Units of Structure"), a procedure developed by Srinivasan and Rose, to provide a physical interpretation of and predict the secondary structures of proteins. The secondary structure type at a given site is identified by the largest conformational bias during short simulations. We examine the rate of successful prediction as a function of temperature and the interaction window. At high temperatures, there is a large propensity for the establishment of beta-strands whereas alpha-helices appear only when the temperature is lower than a certain threshold value. It is found that there exists an optimal temperature at which the correct secondary structures are predicted most accurately. We find that this temperature is close to the peak temperature of the specific heat. Changing the interaction window or carrying out longer simulations approaching equilibrium lead to little change in the optimal success rate. Our findings are in accord with the observation by Srinivasan and Rose that the secondary structures are mainly determined by local interactions and appear in the early stage of folding.     

Slow conformational dynamics of an endonuclease persist in its complex with its natural protein inhibitor

The bacterial toxin colicin E9 is secreted by producing Escherichia coli cells with its 9.5 kDa inhibitor protein Im9 bound tightly to its 14.5 kDa C-terminal DNase domain. Double- and triple-resonance NMR spectra of the isolated DNase domain uniformly labeled with 13C/15N bound to unlabeled Im9 contain more signals than expected for a single DNase conformer, consistent with the bound DNase being present in more than one form. The presence of chemical exchange cross peaks in 750 MHz 15N-1H-15N HSQC-NOESY-HSQC spectra for backbone NH groups of Asp20, Lys21, Trp22, Leu23, Lys69, and Asn70 showed that the bound DNase was in dynamic exchange. The rate of exchange from the major to the minor form was determined to be 1.1 +/- 0.2 s(-1) at 298 K. Previous NMR studies have shown that the free DNase interchanges between two conformers with a forward rate constant of 1.61 +/- 0.11 s(-1) at 288 K, and that the bound Im9 is fixed in one conformation. The NMR studies of the bound DNase show that Im9 binds similarly to both conformers of the DNase and that the buried Trp22 is involved in the dynamic process. For the free DNase, all NH groups within a 9 A radius of any point of the Trp22 ring exhibit heterogeneity suggesting that a rearrangement of the position of this side chain is connected with the conformational interchange. The possible functional significance of this feature of the DNase is discussed.     

Molecular dynamics simulations of bovine lactoferricin: turning a helix into a sheet

Bovine lactoferricin is a 25-residue peptide that is excised through pepsin cleavage in the stomach from the intact 80 kDa bovine milk protein lactoferrin. This basic peptide contains a single disulfide crosslink and is considerably more active as an antimicrobial peptide than the intact protein. It has been suggested that the dramatic difference in potency is related to a change in the secondary and tertiary structure of this peptide, moving from a mixed alpha-helical beta-strand region in the protein to an amphipathic twisted antiparallel beta-sheet in the peptide. Here we have used equilibrium and restrained molecular dynamics calculations to compare the stability of the solution structure of the isolated peptide with that excised from the intact protein. Simulations were performed for fully solvated peptides in the absence and presence of 250 mM salt. Our results show that the peptide as released from the protein is relatively unstable, particularly in the absence of salt. However, even though the simulations extended over 60 nsecs, no interconversion could be observed between the crystal and solution structures, unless a relatively small directional force was exerted on the peptide. A pathway for the structural transition from a helical to a sheet structure was identified in this fashion.     

Determination of protein structures from nuclear magnetic resonance data using a restrained molecular dynamics approach: the lac repressor DNA binding domain

A procedure is described to determine from NMR data the three-dimensional structure of biomolecules. This procedure combines model building with a restrained Molecular Dynamics algorithm, in which distance information from NOEs is incorporated in the form of pseudo potentials. The method has been applied to the N-terminal DNA-binding domain or "headpiece" (amino acids 1-51) of the lac repressor from E. coli, for which no crystal structure is available. The spatial structure of the headpiece is discussed in terms of known physical and biochemical data and of its DNA binding properties.     

A recombinant glutamine-binding protein from Escherichia coli: effect of ligand-binding on protein conformational dynamics

We have investigated the effect of the binding of glutamine on the conformational dynamics of the recombinant glutamine binding protein (GlnBP) from Escherichia coli by steady-state and time-resolved fluorescence techniques. The structural stability of the protein was also studied by far-UV circular dichroism spectroscopy in the range of temperature between 25 and 80 degrees C. The results showed that the interaction of the protein with the ligand resulted in a marked change of the structural and conformational dynamics features of the protein. In particular, the fluorescence and circular dichroism data showed that the presence of glutamine resulted in a dramatic increase of the protein thermal stability of about 10 degrees C. In addition, the fluorescence time-resolved data pointed out that both in the absence and in the presence of glutamine the protein structure was highly rigid with small amplitude of segmental motion up to 65 degrees C and a low accessibility of the protein tryptophan residues to acrylamide. The obtained results on the structural properties of the recombinant glutamine-binding protein in the absence and in the presence of glutamine can contribute to a better understanding of the transport-related functions of the protein and structurally similar periplasmic transport proteins, as well as to the design and development of new biotechnological applications of this class of proteins.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

Checking the pH-induced conformational transition of prion protein by molecular dynamics simulations: effect of protonation of histidine residues

The role of acidic pH in the conversion of human prion protein to the pathogenic isoform is investigated by means of molecular dynamics simulations, focusing the attention on the effect of protonation of histidine residues on the conformational behavior of human PrPC globular domain. Our simulations reveal a significant loss of alpha-helix content under mildly acidic conditions, due to destructuration of the C-terminal part of HB (thus suggesting a possible involvement of HB into the conformational transition leading to the pathogenic isoform) and a transient lengthening of the native beta-sheet. Protonation of His-187 and His-155 seems to be crucial for the onset of the conformational rearrangement. This finding can be related to the existence of a pathogenic mutation, H187R, which is associated with GSS syndrome. Finally, the relevance of our results for the location of a Cu2+-binding pocket in the C-terminal part of the prion is discussed.     

Tropomyosin: Double helix from the protein world

This review concerns the structure and functions of tropomyosin (TM), an actin-binding protein that plays a key role in the regulation of muscle contraction. The TM molecule is a dimer of α-helices, which form a coiled-coil. Recent views on the TM structure are analyzed, and special attention is concentrated on those structural traits of the TM molecule that distinguish it from the other coiled-coil proteins. Modern data are presented on TM functional properties, such as its interaction with actin and ability to move on the surface of actin filaments, which underlies the regulation of the actin-myosin interaction upon contraction of skeletal and cardiac muscles. Also, part of the review is devoted to analysis of the effects of mutations in TM genes associated with muscle diseases (myopathies) on the structure and functions of TM.     

Insights into the conformational flexibility of Bruton's tyrosine kinase from multiple ligand complex structures

Bruton's tyrosine kinase (BTK) plays a key role in B cell receptor signaling and is considered a promising drug target for lymphoma and inflammatory diseases. We have determined the X-ray crystal structures of BTK kinase domain in complex with six inhibitors from distinct chemical classes. Five different BTK protein conformations are stabilized by the bound inhibitors, providing insights into the structural flexibility of the Gly-rich loop, helix C, the DFG sequence, and activation loop. The conformational changes occur independent of activation loop phosphorylation and do not correlate with the structurally unchanged WEI motif in the Src homology 2-kinase domain linker. Two novel activation loop conformations and an atypical DFG conformation are observed representing unique inactive states of BTK. Two regions within the activation loop are shown to structurally transform between 3(10)- and α-helices, one of which collapses into the adenosine-5'-triphosphate binding pocket. The first crystal structure of a Tec kinase family member in the pharmacologically important DFG-out conformation and bound to a type II kinase inhibitor is described. The different protein conformations observed provide insights into the structural flexibility of BTK, the molecular basis of its regulation, and the structure-based design of specific inhibitors.     

Pressure effect on the stability and the conformational dynamics of the D-Galactose/D-Glucose-binding protein from Escherichia coli

The effect of the pressure on the structure and stability of the D-Galactose/D-Glucose binding protein (GGBP) from Escherichia coli was studied by steady-state and time-resolved fluorescence spectroscopy, and the ability of glucose ligand to stabilize the GGBP structure was also investigated. Steady-state fluorescence experiments showed a marked quenching of fluorescence emission of GGBP in the absence of glucose. Instead, the presence of glucose seems to stabilize the structure of GGBP at low and moderate pressure values. Time-resolved fluorescence measurements showed that the GGBP taumean in the absence of glucose varies significantly up to 600 bar, while in the presence of the ligand it is almost unaffected by pressure increase up to 600 bar. The effect of the pressure on GGBP was also studied by molecular dynamics simulations. The simulation data support the spectroscopic results and confirm that the presence of glucose is able to contrast the negative effects of pressure on the protein structure. Taken together, the spectroscopic and computer simulation studies suggest that at pressure values up to 2000 bar the structure of GGBP in the absence of glucose remains folded, but a significant perturbation of the protein secondary structures can be detected. The binding of glucose reduces the negative effect of pressure on protein structure and confers protection from perturbation especially at moderate pressure values.     

Prediction of transition metal-binding sites from apo protein structures

Metal ions are crucial for protein function. They participate in enzyme catalysis, play regulatory roles, and help maintain protein structure. Current tools for predicting metal-protein interactions are based on proteins crystallized with their metal ions present (holo forms). However, a majority of resolved structures are free of metal ions (apo forms). Moreover, metal binding is a dynamic process, often involving conformational rearrangement of the binding pocket. Thus, effective predictions need to be based on the structure of the apo state. Here, we report an approach that identifies transition metal-binding sites in apo forms with a resulting selectivity >95%. Applying the approach to apo forms in the Protein Data Bank and structural genomics initiative identifies a large number of previously unknown, putative metal-binding sites, and their amino acid residues, in some cases providing a first clue to the function of the protein.     

Slow conformational dynamics of the human A2A adenosine receptor are temporally ordered

1. Download : Download high-res image (226KB)
2. Download : Download full-size image

     

Increased sequence hydrophobicity reduces conformational specificity: A mutational case study of the Arc repressor protein

The amino-acid sequences of soluble, globular proteins must have hydrophobic residues to form a stable core, but excess sequence hydrophobicity can lead to loss of native state conformational specificity and aggregation. Previous studies of polar-to-hydrophobic mutations in the β-sheet of the Arc repressor dimer showed that a single substitution at position 11 (N11L) leads to population of an alternate dimeric fold in which the β-sheet is replaced by helix. Two additional hydrophobic mutations at positions 9 and 13 (Q9V and R13V) lead to population of a differently folded octamer along with both dimeric folds. Here we conduct a comprehensive study of the sequence determinants of this progressive loss of fold specificity. We find that the alternate dimer-fold specifically results from the N11L substitution and is not promoted by other hydrophobic substitutions in the β-sheet. We also find that three highly hydrophobic substitutions at positions 9, 11, and 13 are necessary and sufficient for oligomer formation, but the oligomer size depends on the identity of the hydrophobic residue in question. The hydrophobic substitutions increase thermal stability, illustrating how increased hydrophobicity can increase folding stability even as it degrades conformational specificity. The oligomeric variants are predicted to be aggregation-prone but may be hindered from doing so by proline residues that flank the β-sheet region. Loss of conformational specificity due to increased hydrophobicity can manifest itself at any level of structure, depending upon the specific mutations and the context in which they occur.