Structure of the leech protein saratin and characterization of its binding to collagen

The leech protein Saratin from Hirudo medicinalis prevents thrombocyte aggregation by interfering with the first binding step of the thrombocytes to collagen by binding to collagen. We solved the three-dimensional structure of the leech protein Saratin in solution and identified its collagen binding site by NMR titration experiments. The NMR structure of Saratin consists of one α-helix and a five-stranded β-sheet arranged in the topology ββαβββ. The C-terminal region, of about 20 amino acids in length, adopts no regular structure. NMR titration experiments with collagen peptides show that the collagen interaction of Saratin takes place in a kind of notch that is formed by the end of the α-helix and the β-sheet. NMR data-driven docking experiments to collagen model peptides were used to elucidate the putative binding mode of Saratin and collagen. Mainly, parts of the first and the end of the fifth β-strand, the loop connecting the α-helix and the third β-strand, and a short part of the loop connecting the fourth and fifth β-strand participate in binding.     

Intrinsic secondary structure propensities of the amino acids,using statistical ϕ–ψ matrices: Comparison with experimental scales

Today there are several different experimental scales for the intrinsic α-helix as well as β-strand, propensities of the 20 amino acids obtained from the thermodynamic analysis of various model systems. These scales do not compare well with those extracted from statistical analysis of three-dimensional structure databases. Possible explanations for this could be the limited size of the databases used, the definitions of intrinsic propensities, or the theoretical approach. Here we report a statistical determination of α-helix and β-strand propensities derived from the analysis of a database of 279 three-dimensional structures. Contrary to what has been generally done, we have considered a particular residue as in α-helix or β-strand conformation by looking only at its dihedral angles (?–ψ matrices). Neither the identity nor the conformation of the surrounding residues in the amino acid sequence has been taken into consideration. Pseudoenergy empirical scales have been calculated from the statistical propensities. These scales agree very well with the experimental ones in relative and absolute terms. Moreover, its correlation with the average of the experimental scales for α-helix or β-strand is as good as the correlations of the individual experimental scales with the average. These results show that by using a large enough database and a proper definition for the secondary structure propensities, it is possible to obtain a scale as good as any of experimental origin. Interestingly the ?–ψ analysis of the Ramachandran plot suggests that the amino acids could have different β-strand propensities in different subregions of the β-strand area. © 1994 Wiley-Liss, Inc.     

Secondary-structure analysis of alcohol-denatured proteins by vacuum-ultraviolet circular dichroism spectroscopy

To elucidate the structural characteristics of alcohol-denatured proteins, we measured the vacuum-ultraviolet circular dichroism (VUVCD) spectra of six proteins-myoglobin, human serum albumin, α-lactalbumin, thioredoxin, β-lactoglobulin, and α-chymotrypsinogen A-down to 170 nm in trifluoroethanol solutions (TFE: 0-50%) and down to 175 nm in methanol solutions (MeOH: 0-70%) at pH 2.0 and 25°C, using a synchrotron-radiation VUVCD spectrophotometer. The contents of α-helices, β-strands, turns, poly-L-proline type II helices (PPIIs), and unordered structures of these proteins were estimated using the SELCON3 program, including the numbers of α-helix and β-strand segments. Furthermore, the positions of α-helices and β-strands on amino acid sequences were predicted by combining these secondary-structure data with a neural-network method. All alcohol-denatured proteins showed higher α-helix contents (up to ~ 90%) compared with the native states, and they consisted of several long helical segments. The helix-forming ability was higher in TFE than in MeOH, whereas small amounts of β-strands without sheets were formed in the MeOH solution. The produced α-helices were transformed dominantly from the β-strands and unordered structures, and slightly from the turns. The content and mean length of α-helix segments decreased as the number of disulfide bonds in the proteins increased, suggesting that disulfide bonds suppress helix formation by alcohols. These results demonstrate that alcohol-denatured proteins constitute an ensemble of many long α-helices, a few β-strands and PPIIs, turns, and unordered structures, depending on the types of proteins and alcohols involved.     

The Role of Arrestin α-Helix I in Receptor Binding

Arrestins rapidly bind phosphorylated activated forms of their cognate G protein-coupled receptors, thereby preventing G protein coupling and often switching signaling to other pathways. Amphipathic α-helix I (residues 100-111) has been implicated in receptor binding, but the mechanism of its action has not been determined yet. Here we show that several mutations in the helix itself and in adjacent hydrophobic residues in the body of the N-domain reduce arrestin1 binding to light-activated phosphorylated rhodopsin (P-Rh^?). On the background of phosphorylation-independent mutants that bind with high affinity to both P-Rh^? and light-activated unphosphorylated rhodopsin, these mutations reduce the stability of the arrestin complex with P-Rh^?, but not with light-activated unphosphorylated rhodopsin. Using site-directed spin labeling, we found that the local structure around α-helix I changes upon binding to rhodopsin. However, the intramolecular distances between α-helix I and adjacent β-strand I (or the rest of the N-domain), measured using double electron-electron resonance, do not change, ruling out relocation of the helix due to receptor binding. Collectively, these data demonstrate that α-helix I plays an indirect role in receptor binding, likely keeping β-strand I, which carries several phosphate-binding residues, in a position favorable for its interaction with receptor-attached phosphates.     

The role of negative selection in protein evolution revealed through the energetics of the native state ensemble

Knowing the determinants of conformational specificity is essential for understanding protein structure, stability, and fold evolution. To address this issue, a novel statistical measure of energetic compatibility between sequence and structure was developed using an experimentally validated model of the energetics of the native state ensemble. This approach successfully matched sequences from a diverse subset of the human proteome to their respective folds. Unexpectedly, significant energetic compatibility between ostensibly unrelated sequences and structures was also observed. Interrogation of these matches revealed a general framework for understanding the origins of conformational specificity within a proteome: specificity is a complex function of both the ability of a sequence to adopt folds other than the native, and ability of a fold to accommodate sequences other than the native. The regional variation in energetic compatibility indicates that the compatibility is dominated by incompatibility of sequence for alternative fold segments, suggesting that evolution of protein sequences has involved substantial negative selection, with certain segments serving as “gatekeepers” that presumably prevent alternative structures. Beyond these global trends, a size dependence exists in the degree to which the energetic compatibility is determined from negative selection, with smaller proteins displaying more negative selection. This partially explains how short sequences can adopt unique folds, despite the higher probability in shorter proteins for small numbers of mutations to increase compatibility with other folds. In providing evolutionary ground rules for the thermodynamic relationship between sequence and fold, this framework imparts valuable insight for rational design of unique folds or fold switches. Proteins 2016; 84:435–447. © 2016 Wiley Periodicals, Inc.     

Studies on the rules of β-strand alignment in a protein β-sheet structure

To further disclose the underlying mechanisms of protein β-sheet formation, studies were made on the rules of β-strands alignment forming β-sheet structure using statistical and machine learning approaches. Firstly, statistical analysis was performed on the sum of β-strands between each β-strand pairs in protein sequences. The results showed a propensity of near-neighbor pairing (or called “first come first pair”) in the β-strand pairs. Secondly, based on the same dataset, the pairwise cross-combinations of real β-strand pairs and four pseudo-β-strand contained pairs were classified by support vector machine (SVM). A novel feature extracting approach was designed for classification using the average amino acid pairing encoding matrix (APEM). Analytical results of the classification indicated that a segment of β-strand had the ability to distinguish β-strands from segments of α-helix and coil. However, the result also showed that a β-strand was not strongly conserved to choose its real partner from all the alternative β-strand partners, which was corresponding with the ordination results of the statistical analysis each other. Thus, the rules of “first come first pair” propensity and the non-conservative ability to choose real partner, were possible important factors affecting the β-strands alignment forming β-sheet structures.     

A Secondary Structure Model of the Integrin a Subunit N-Terminal Domain Based on Analysis of Multiple Alignments

The integrins are α/β heterodimeric proteins which mediate cell-matrix and cell-cell inter-actions. Current data indicate that the N-terminal moiety of the a subunit is involved in ligand binding. This region of the receptor is made up of a seven-fold repeated sequence of unknown structure which contains EF-hand-like putative divalent cation-binding sites. Recent studies have shown that multiple sequence alignments can be analysed to yield secondary structure predictions. Therefore, to obtain a model structure for the integrin a subunit N-terminal domain repeat, a large alignment of the seven repeats from sixteen integrin sequences was generated. Two methods of analysis were used: First, Chou and Fasman and Garnier, Osguthorpe and Robson predictions were carried out for individual sequences and the consensus predictions derived. Consensus hydrophobicity and chain flexibility data were also used to provide additional data. Second, sites of conservation and variation were analysed by a computer program STAMA (STructure After Multiple Alignment) to yield a secondary structure prediction. The two analyses gave essentially the same predicted structure: undefined region, loop, α-helix, β-strand, divalent cation-binding loop, β-strand, putative turn, loop, β-strand. This is the first model structure to be presented for an integrin domain. Its implications for integrin function are discussed.     

Multiscale simulations of protein folding: application to formation of secondary structures

A multiscale simulation method of protein folding is proposed, using atomic representation of protein and solvent, combing genetic algorithms to determine the key protein structures from a global view, with molecular dynamic simulations to reveal the local folding pathways, thus providing an integrated landscape of protein folding. The method is found to be superior to previously investigated global search algorithms or dynamic simulations alone. For secondary structure formation of a selected peptide, RN24, the structures and dynamics produced by this method agree well with corresponding experimental results. Three most populated conformations are observed, including hairpin, β-sheet and α-helix. The energetic barriers separating these three structures are comparable to the kinetic energy of the atoms of the peptide, implying that the transition between these states can be easily triggered by kinetic perturbations, mainly through electrostatic interactions between charged atoms. Transitions between α-helix and β-sheet should jump over at least two energy barriers and may stay in the energetic trap of hairpin. It is proposed that the structure of proteins should be jointly governed by thermodynamic and dynamic factors; free energy is not the exclusive dominant for stability of proteins.     

The Cancer Mutation D83V Induces an α-Helix to β-Strand Conformation Switch in MEF2B

MEF2B is a major target of somatic mutations in non-Hodgkin lymphoma. Most of these mutations are non-synonymous substitutions of surface residues in the MADS-box/MEF2 domain. Among them, D83V is the most frequent mutation found in tumor cells. The link between this hotspot mutation and cancer is not well understood. Here we show that the D83V mutation induces a dramatic α-helix to β-strand switch in the MEF2 domain. Located in an α-helix region rich in β-branched residues, the D83V mutation not only removes the extensive helix stabilization interactions but also introduces an additional β-branched residue that further shifts the conformation equilibrium from α-helix to β-strand. Cross-database analyses of cancer mutations and chameleon sequences revealed a number of well-known cancer targets harboring β-strand favoring mutations in chameleon α-helices, suggesting a commonality of such conformational switch in certain cancers and a new factor to consider when stratifying the rapidly expanding cancer mutation data.     

Folding simulations of three proteins having all α-helix,all β-strand and α/β-structures

We performed folding simulations of three proteins using four force fields, AMBER parm96, AMBER parm99, CHARMM 27 and OPLS-AA/L, in order to examine the features of these force fields. We studied three proteins, protein A (all α-helix), cold-shock protein (all β-strand) and protein G (α/β-structures), for the folding simulations. For the simulation, we used the simulated annealing molecular dynamics method, which was performed 50 times for each protein using the four force fields. The results showed that the secondary-structure-forming tendencies are largely different among the four force fields. AMBER parm96 favours β-bridge structures and extended β-strand structures, and AMBER parm99 favours α-helix structures and 3₁₀-helix structures. CHARMM 27 slightly favours α-helix structures, and there are also π-helix and β-bridge structures. OPLS-AA/L favours α-helix structures and 3₁₀-helix structures.     

Secondary structure perturbations in salt-induced protein precipitates

The secondary structure implications of precipitation induced by a chaotropic salt, KSCN, and a structure stabilizing salt, Na₂SO₄, were studied for twelve different proteins. α-helix and β-sheet content of precipitate and native structures were estimated from the analysis of amide I band Raman spectra. A statistical analysis of the estimated perturbations in the secondary structure contents indicated that the most significant event is the formation of β-sheet structures with a concomitant loss of α-helix on precipitation with KSCN. The conformational changes for each protein were also analyzed with respect to elements of primary, secondary and tertiary structure existing in the native protein; primary structure was quantified by the fractions of hydrophobic and charged amino acids, secondary structure by x-ray estimates of α-helix and β-sheet contents of native proteins and tertiary structure by the dipole moment and solvent-accessible surface area. For the KSCN precipitates, factors affecting β-sheet content included the fraction of charged amino acids in the primary sequence and the surface area. Changes in α-helix content were influenced by the initial helical content and the dipole moment. The enhanced β-sheet contents of precipitates observed in this work parallel protein structural changes occurring in other aggregative phenomena.     

Amyloid-like fibrils from a domain-swapping protein feature a parallel, in-register conformation without native-like interactions

The formation of amyloid-like fibrils is characteristic of various diseases, but the underlying mechanism and the factors that determine whether, when, and how proteins form amyloid, remain uncertain. Certain mechanisms have been proposed based on the three-dimensional or runaway domain swapping, inspired by the fact that some proteins show an apparent correlation between the ability to form domain-swapped dimers and a tendency to form fibrillar aggregates. Intramolecular β-sheet contacts present in the monomeric state could constitute intermolecular β-sheets in the dimeric and fibrillar states. One example is an amyloid-forming mutant of the immunoglobulin binding domain B1 of streptococcal protein G, which in its native conformation consists of a four-stranded β-sheet and one α-helix. Under native conditions this mutant adopts a domain-swapped dimer, and it also forms amyloid-like fibrils, seemingly in correlation to its domain-swapping ability. We employ magic angle spinning solid-state NMR and other methods to examine key structural features of these fibrils. Our results reveal a highly rigid fibril structure that lacks mobile domains and indicate a parallel in-register β-sheet structure and a general loss of native conformation within the mature fibrils. This observation contrasts with predictions that native structure, and in particular intermolecular β-strand interactions seen in the dimeric state, may be preserved in "domain-swapping" fibrils. We discuss these observations in light of recent work on related amyloid-forming proteins that have been argued to follow similar mechanisms and how this may have implications for the role of domain-swapping propensities for amyloid formation.     

β2-Strand of salivary S cystatins: A “chemeleon sequence”

Secondary structure prediction of salivary cystatins S, SA, and SN carried out by several methods label the 39-58 sequence (β2-strand) as predominantly α-helical. The helical propensity of a peptide corresponding to β2-strand of salivary SA cystatin analyzed by CD display high helical propensity in aqueous solution, whereas peptides matching the β2-strand amino acid sequence of cystatins S and SN, display random coil conformation in aqueous solution but acquire α-helical conformation in the presence of trifluoroethanol (TFE). Moreover molecular dynamics simulation performed on the homology modeling of cystatin SA constructed on the basis of recently determined three-dimensional structure of salivary cystatin D, suggests that cystatin SA does not significantly deviate from the starting structure over the course of the simulation. The results obtained indicate that the β2-strand of salivary S cystatins has high helical propensity when isolated from native protein and acquire the final β structure by interaction with the rest of the polypeptide chain.     

Beta conformation of polyglutamine track revealed by a crystal structure of Huntingtin N-terminal region with insertion of three histidine residues

Huntington disease is an autosomal-dominant neurodegenerative disorder caused by a polyglutamine (polyQ) expansion (> 35Q) in the first exon (EX1) of huntingtin protein (Htt). mHtt protein is thought to adopt one or more toxic conformation(s) that are involved in pathogenic interactions in cells . However, the structure of mHtt is not known. Here, we present a near atomic resolution structure of mHtt36Q-EX1. To facilitate crystallization, three histidine residues (3H) were introduced within the Htt36Q stretch resulting in the sequence of Q₇HQHQHQ₂₇. The Htt36Q3H region adopts α-helix, loop, β-hairpin conformations. Furthermore, we observed interactions between the backbone of the Htt36Q3H β-strand with the aromatic residues mimicking putative-toxic interactions with other proteins. Our findings support previous predictions that the expanded mHtt-polyQ region adopts a β-sheet structure. Detailed structural information about mHtt improves our understanding of the pathogenic mechanisms in HD and other polyQ expansion disorders and may form the basis for rational design of small molecules that target toxic conformations of disease-causing proteins.     

Origins of structural diversity within sequentially identical hexapeptides

Efforts to predict protein secondary structure have been hampered by the apparent structural plasticity of local amino acid sequences. Kabsch and Sander (1984, Proc. Natl. Acad. Sci. USA 81, 1075–1078) articulated this problem by demonstrating that identical pentapeptide sequences can adopt distinct structures in different proteins. With the increased size of the protein structure database and the availability of new methods to characterize structural environments, we revisit this observation of structural plasticity. Within a set of proteins with less than 50% sequence identity, 59 pairs of identical hexapeptide sequences were identified. These local structures were compared and their surrounding structural environments examined. Within a protein structural class (α/α, β/β, α/β, α + β), the structural similarity of sequentially identical hexapeptides usually is preserved. This study finds eight pairs of identical hexapeptide sequences that adopt β-strand structure in one protein and α-helical structure in the other. In none of the eight cases do the members of these sequence pairs come from proteins within the same folding class. These results have implications for class dependent secondary structure prediction algorithms.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

小麦多聚半乳糖醛酸酶抑制蛋白的部分结构

为了弄清小麦多聚半乳糖醛酸酶抑制蛋白 (polygalacturonase inhibitingprotein ,PGIP)的作用机制 ,并为其在基因工程中的应用提供依据 ,对其结构进行了研究 .用Edman降解法测得小麦PGIP的N端序列为Lys Pro Leu Leu Thr Lys Ile Thr Lys Gly Ala Ala Ser Thr .用CD谱研究其二级结构 ,发现小麦PGIP天然态含有 4 3 7%的 β折叠和 13 1%的α螺旋 .酸碱和温度变性引起了二级结构改变 .不完全变性阶段 ,二级结构的变化表现为α螺旋无明显变化 ,β折叠遭到破坏 ;活性完全丧失阶段 ,β折叠变化很小 ,α螺旋含量明显减少 .用NR R(非还原 还原 )双向对角线SDS PAGE鉴定出小麦PGIP含有链内二硫键 .用去糖基化法确证了小麦PGIP的糖含量为 2 2 %.小麦PGIP与双子叶植物PGIP相比 ,一级结构差异较大 ,同源性由 36 %变为 9%;二级结构相似 ,都是高 β 折叠的蛋白 ;均具有链内二硫键 ;在糖含量上也相似 .研究结果为进一步弄清小麦PGIP作用机理打下了基础 ,同时对于植物抗赤霉病基因工程具有重要意义 .     

NMR structure of the mengovirus Leader protein zinc-finger domain

The Leader protein is a defining feature of picornaviruses from the Cardiovirus genus. This protein was recently shown to inhibit cellular nucleocytoplasmic transport through an activity mapped to its zinc-binding region. Here we report the three-dimensional solution structure determined by nuclear magnetic resonance (NMR) spectroscopy of this domain (residues 5-28) from mengovirus. The domain forms a CHCC zinc-finger with a fold comprising a β-hairpin followed by a short α-helix that can adopt two different conformations. This structure is divergent from those of other eukaryotic zinc-fingers and instead resembles motifs found in a group of DNA-binding proteins from Archaea.     

Local secondary structure in ribonuclease A denatured by guanidine · HCl near 1 °C

Recent work has shown that a short α-helix can be stable in water near 1 °C when stabilized by specific interactions between side-chains, while earlier “host-guest” results with random copolymers have shown that a short α-helix is unstable in water at all temperatures in the absence of stabilizing side-chain interactions. As regards the mechanism of protein folding, it is now reasonable on energetic grounds to consider isolated α-helices and β-sheets as the first intermediates on the pathway of protein folding. Proton nuclear magnetic resonance is used here to detect isolated secondary structures in ribonuclease A denatured by guanidine · HCl (GuHCl). Temperatures near 1 °C are used because the low-temperature stability of the C-peptide helix may be a general property of isolated secondary structures in water.Our procedure is to titrate with GuHCl the C2H resonance lines of the four histidine residues of denatured ribonuclease A. Studies of model peptides (C-peptide (lactone) and C-peptide carboxylate, residues 1 to 13 of ribonuclease A; S-peptide, residues 1 to 20) show linear titration curves for the C2H resonance of His12 above 0.5 M-GuHCl, once helix unfolding is complete. Deviations from this line are used to monitor helix formation. The GuHCl titration curves of the other three histidine residues are also linear, once unfolding is complete. The results show that the helix found in C-peptide and S-peptide is also found in denatured ribonuclease A, where it behaves as an isolated helix not stabilized significantly by interactions with other chain segments. Studies of denatured S-protein show that the remaining three His residues, His48, His105 and His119, are involved in structure only below 1 m-GuHCl at 9 °C, pH 1.9. The nature of this structure is not known. The main conclusion from this work is that the His12 helix can be observed as a stable, isolated helix in denatured ribonuclease A near 1 °C, and that none of the other three His residues is involved in a comparably stable local structure. In native ribonuclease A, His12 is within an α-helix and the other three His residues are involved in a 3-stranded β-sheet structure.The helix-coil transition of C-peptide has also been studied for other side-chain resonances by GuHCl titration. Typically, but not always, the titration curves are linear after helix unfolding takes place and resonance lines from different residues of the same amino acid type can be resolved in GuHCl solutions. This is true of the four histidine residues of ribonuclease A although their pK values in 5 m-GuHCl are nearly the same. In C-peptide, the βCH₃ resonance of Ala6 is affected strongly by GuHCl while the lines of Ala4 and Ala5 are shifted only weakly by GuHCl. Evidently the interactions between GuHCl and side-chains in an unfolded peptide depend upon neighboring groups.