Mapping long-range contacts in a highly unfolded protein

Insights into the earliest events in protein folding can be obtained by analysis of the conformational propensities of unfolded or partly folded states. The structure of the acid-unfolded state of apomyoglobin has been characterized using paramagnetic spin labeling and NMR. Nitroxide side-chains, introduced by coupling to mutant cysteine residues at positions 18, 77, and 133, were used as probes of chain compaction and long-range tertiary contacts. Significant interactions are observed within and between the N and C termini, while the central region of the polypeptide chain behaves as a random polymer. Even in this highly denatured form, the protein samples transient compact states in which there are native-like contacts between the N and C-terminal regions.     

Silkworm, Bombyx mori larvae expressed the spider silk protein through a novel Bac-to-Bac/BmNPV baculovirus

Abstract:  The silkworm has become an ideal multicellular eukaryotic model system for basic research. The major advantages of expressing foreign genes in silkworm larvae are the low cost of feeding, the extremely high levels of expression achievable compared with expression in cell lines and increased safety because the baculovirus is noninfectious to vertebrates. In this study, we used a recently developed Bombyx mori Nucleopolyhedrovirus (BmNPV) bacmid to express the spider flagelliform silk gene in silkworm larvae. The recombinant bacmid baculoviruses (rBacmid/BmNPV/Flag) were introduced into the first-day larvae of the fifth instar by subcutaneous injection. The worms presented symptoms typical of NPV infection from 72 h after injection compared with control. The haemolymph was collected from the infected larvae 120 h post-infection and the recombinant 6× His-tagged Flag protein was purified by the Ni-NTA spin kit under denaturing conditions with 8  m urea. A 37.0-kDa protein was visualized both in rBacmid/BmNPV/Flag-infected haemolymph and eluting fraction. The results showed that the Bac-to-Bac/BmNPV baculovirus expression system is an efficient tool to express the target gene in silkworm larvae, which takes only 7–10 days for generating recombinant baculovirus, compared with the traditional homologous recombination method, which needs at least 40 days for multiple rounds of purification and amplification of viruses.     

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.     

pH-dependent stability and folding kinetics of a protein with an unusual alpha-beta topology: the C-terminal domain of the ribosomal protein L9

The folding kinetics and thermodynamics of the isolated C-terminal domain of the ribosomal protein L9 (CTL9) have been studied as a function of pH. CTL9 is an alpha-beta protein that contains a single beta-sheet with an unusual mixed parallel, anti-parallel topology. The folding is fully reversible and two-state over the entire pH range. Stopped-flow fluorescence and CD experiments yield the same folding rate, and the chevron plots have the characteristic V-shape expected for two-state folding. The values of DeltaG*(H2O) and the m value calculated from the kinetic experiments are in excellent agreement with the equilibrium measurements. The extrapolated initial amplitudes of both the stopped-flow fluorescence and CD measurements show that there is no detectable burst phase intermediate. The domain contains three histidine residues, two of which are largely buried in the native state. They do not participate in salt-bridges or take part in a hydrogen bonded network. NMR measurements reveal that the buried histidine residues have significantly perturbed pK(a) values in the native state. The equilibrium stability and the folding rate are found to be strongly dependent upon their ionization state. There is a linear relationship between the log of the folding rate and DeltaG* (H2O) . The protein is much more stable and folds noticeably faster at pH values above the native state pK(a) values. DeltaG*(H2O) of unfolding increases from 2.90 kcal mol(-1) at pH 5.0 to 6.40 kcal mol(-1) at pH 8.0 while the folding rate increases from 0.60 to 18.7 s(-1). Tanford linkage analysis revealed that the interactions involving the two histidine residues are largely developed in the transition state. The results are compared to other studies of the pH-dependence of folding.     

Intra- versus intermolecular interactions in monellin: contribution of surface charges to protein assembly

The relative significance of weak non-covalent interactions in biological context has been much debated. Here, we have addressed the contribution of Coulombic interactions to protein stability and assembly experimentally. The sweet protein monellin, a non-covalently linked heterodimeric protein, was chosen for this study because of its ability to spontaneously reconstitute from separated fragments. The reconstitution of monellin mutants containing large surface charge perturbations was compared to the thermostability of structurally equivalent single-chain monellin containing the same sets of mutations under varying salt concentrations. The affinity between monellin fragments is found to correlate with the thermostability of single chain monellin, indicating the involvement of the same underlying Coulombic interactions. This confirms that there are no principal differences in the interactions involved in folding and binding. Based on comparison with a previous mutational study involving hydrophobic core residues, the relative contribution of Coulombic interactions to stability and affinity is modest. However, the Coulombic perturbations only affect the association rates of reconstitution in contrast to perturbations involving hydrophobic residues, which affect primarily the dissociation rates. These results indicate that Coulombic interactions are likely to be of main importance for the association of protein assembly, relevant for functions of proteins.     

Very fast folding and association of a trimerization domain from bacteriophage T4 fibritin

The foldon domain constitutes the C-terminal 30 amino acid residues of the trimeric protein fibritin from bacteriophage T4. Its function is to promote folding and trimerization of fibritin. We investigated structure, stability and folding mechanism of the isolated foldon domain. The domain folds into the same trimeric beta-propeller structure as in fibritin and undergoes a two-state unfolding transition from folded trimer to unfolded monomers. The folding kinetics involve several consecutive reactions. Structure formation in the region of the single beta-hairpin of each monomer occurs on the submillisecond timescale. This reaction is followed by two consecutive association steps with rate constants of 1.9(+/-0.5)x10(6)M(-1)s(-1) and 5.4(+/-0.3)x10(6)M(-1)s(-1) at 0.58 M GdmCl, respectively. This is similar to the fastest reported bimolecular association reactions for folding of dimeric proteins. At low concentrations of protein, folding shows apparent third-order kinetics. At high concentrations of protein, the reaction becomes almost independent of protein concentrations with a half-time of about 3 ms, indicating that a first-order folding step from a partially folded trimer to the native protein (k=210 +/- 20 s(-1)) becomes rate-limiting. Our results suggest that all steps on the folding/trimerization pathway of the foldon domain are evolutionarily optimized for rapid and specific initiation of trimer formation during fibritin assembly. The results further show that beta-hairpins allow efficient and rapid protein-protein interactions during folding.     

拟蜘蛛拖丝蛋白基因重组蛋白的纯化及其抗血清制备

为了制备可用于拟蜘蛛拖丝蛋白体外检测的抗血清,利用前期工作中获得的蜘蛛拖丝蛋白基因重组原核表达载体2S-pET-52b( + ),采用原核表达方法获得大量重组蜘蛛拖丝蛋白2S-His,对此重组蛋白进行His标签特异性的亲和纯化,再依次经SDS-PAGE、切胶和与佐剂混合后作为抗原;皮下多点注射法将抗原注入新西兰大白兔皮...     

Mapping the essential structures of human ribosomal protein L7 for nuclear entry, ribosome assembly and function

Human large subunit protein L7 carries multiple nuclear localization signals (NLS) in its structure: there are three monobasic partite NLSs at the NH2-region of the first 54 amino acid residues and a bipartite in the middle section at position of 156-167. The C-region of the last 50 amino acid residues displays membrane binding nature, and might involve in forming a nuclear microbody for pre-nucleolar ribosome assembly. The middle section covers 144 amino acid residues which are essential for the structure and function of ribosome. This is evident from findings that truncated L7 without the NH2-region or the C-region, or missing both regions, is capable of reaching nucleolus and incorporating in ribosome, however, only ribosomes bearing truncated L7 without the NH2-region is capable of engaging in polysome formation. Combining with the phylogenic findings from homologous sequence alignment, the NH2-region of L7, besides being as a eukaryotic expansion segment, can be excluded from building a functional eukaryotic ribosome.     

Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state

Comparatively little is known about the role of non-native interactions in protein folding and their role in both folding and stability is controversial. We demonstrate that non-native electrostatic interactions involving specific residues in the denatured state can have a significant effect upon protein stability and can persist in the transition state for folding. Mutation of a single surface exposed residue, Lys12 to Met, in the N-terminal domain of the ribosomal protein L9 (NTL9), significantly increased the stability of the protein and led to faster folding. Structural and energetic studies of the wild-type and K12M mutant show that the 1.9 kcal mol(-1) increase in stability is not due to native state effects, but rather is caused by modulation of specific non-native electrostatic interactions in the denatured state. pH dependent stability measurements confirm that the increased stability of the K12M is due to the elimination of favorable non-native interactions in the denatured state. Kinetic studies show that the non-native electrostatic interactions involving K12 persist in the transition state. The analysis demonstrates that canonical Phi-values can arise from the disruption of non-native interactions as well as from the development of native interactions.     

Surprisingly high correlation between early and late stages in non-two-state protein folding

We collected quantitative kinetic data on early and late stages of folding in non-two-state proteins from the literature, and studied the relationship between the kinetics of the two stages. There was a surprisingly high correlation between the rate constants of these stages. The correlation coefficient of the logarithmic rate constants was as high as 0.97, which could not be caused by chance. We also studied relationships of the logarithmic rate constants of the two stages with native three-dimensional structures represented by the residue-residue contact map. There were again surprisingly high correlations between the logarithmic rate constants and the number of non-local contact clusters obtained from the contact maps. Because the number of non-local contact clusters represents overall arrangement of substructures in a native protein, the results strongly suggested the importance of the arrangement of the substructures for the kinetics of both early and late stages of protein folding.     

Generation of the BfiI restriction endonuclease from the fusion of a DNA recognition domain to a non-specific nuclease from the phospholipase D superfamily

The BfiI endonuclease cleaves DNA at fixed positions downstream of an asymmetric sequence. Unlike other restriction enzymes, it functions without metal ions. The N-terminal half of BfiI is similar to Nuc, an EDTA-resistant nuclease from Salmonella typhimurium that belongs to the phosphoplipase D superfamily. Nuc is a dimer with one active site at its subunit interface, as is BfiI, but it cuts DNA non-specifically. BfiI was cleaved by thermolysin into an N-terminal domain, which forms a dimer with non-specific nuclease activity, and a C-terminal domain, which lacks catalytic activity but binds specifically to the recognition sequence as a monomer. On denaturation with guanidinium, BfiI underwent two unfolding transitions: one at a relatively low concentration of guanidinium, to a dimeric non-specific nuclease; a second at a higher concentration, to an inactive monomer. The isolated C-terminal domain unfolded at the first (relatively low) concentration, the isolated N-terminal at the second. Hence, BfiI consists of two physically separate domains, with catalytic and dimerisation functions in the N terminus and DNA recognition functions in the C terminus. It is the first example of a restriction enzyme generated by the evolutionary fusion of a DNA recognition domain to a phosphodiesterase from the phospholipase D superfamily. BfiI may consist of three structural units: a stable central core with the active site, made from two copies of the N-terminal domain, flanked by relatively unstable C-terminal domains, that each bind a copy of the recognition sequence.     

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer＇s disease, Parkinson＇s disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggre- gation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute signifi- cantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both aU-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.     

Structural and thermodynamic characterization of the adrenodoxin-like domain of the electron-transfer protein Etp1 from Schizosaccharomyces pombe

The protein Etp1 of Schizosaccharomyces pombe consists of an amino-terminal COX15-like domain and a carboxy-terminal ferredoxin-like domain, Etp1^fd, which is cleaved off after mitochondrial import. The physiological function of Etp1^fd is supposed to lie in the participation in the assembly of iron-sulfur clusters and the synthesis of heme A. In addition, the protein was shown to be the first microbial ferredoxin being able to support electron transfer in mitochondrial steroid hydroxylating cytochrome P450 systems in vivo and in vitro, replacing thereby the native redox partner, adrenodoxin. Despite a sequence similarity of 39% and the fact that fission yeast is a mesophilic organism, thermodynamic studies revealed that Etp1^fd has a melting temperature more than 20 °C higher than adrenodoxin. The three-dimensional structure of Etp1^fd has been determined by crystallography. To the best of our knowledge it represents the first three-dimensional structure of a yeast ferredoxin. The structure-based sequence alignment of Etp1^fd with adrenodoxin yields a rational explanation for their observed mutual exchangeability in the cytochrome P450 system. Analysis of the electron exchange with the S. pombe redox partner Arh1 revealed differences between Etp1^fd and adrenodoxin, which might be linked to their different physiological functions in the mitochondria of mammals and yeast.     

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.     

Fast and slow intermediate accumulation and the initial barrier mechanism in protein folding

Do stable intermediates form very early in the protein folding process? New results and a quantity of literature that bear on this issue are examined here. Results available provide little support for early intermediate accumulation before an initial search-dependent nucleation barrier.     

SnapDRAGON: a method to delineate protein structural domains from sequence data

We describe a method to identify protein domain boundaries from sequence information alone based on the assumption that hydrophobic residues cluster together in space. SnapDRAGON is a suite of programs developed to predict domain boundaries based on the consistency observed in a set of alternative ab initio three-dimensional (3D) models generated for a given protein multiple sequence alignment. This is achieved by running a distance geometry-based folding technique in conjunction with a 3D-domain assignment algorithm. The overall accuracy of our method in predicting the number of domains for a non-redundant data set of 414 multiple alignments, representing 185 single and 231 multiple-domain proteins, is 72.4 %. Using domain linker regions observed in the tertiary structures associated with each query alignment as the standard of truth, inter-domain boundary positions are delineated with an accuracy of 63.9 % for proteins comprising continuous domains only, and 35.4 % for proteins with discontinuous domains. Overall, domain boundaries are delineated with an accuracy of 51.8 %. The prediction accuracy values are independent of the pair-wise sequence similarities within each of the alignments. These results demonstrate the capability of our method to delineate domains in protein sequences associated with a wide variety of structural domain organisation.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Fusion of barnase to antiferritin antibody F11 VH domain results in a partially folded functionally active protein

A chimeric protein, VH-barnase, was obtained by fusing the VH domain of anti-human ferritin monoclonal antibody F11 to barnase, a bacterial RNase from Bacillus amyloliquefaciens. After refolding from inclusion bodies, the fusion protein formed insoluble aggregates. Off-pathway aggregation was significantly reduced by adding either purified GroEL/GroES chaperones or arginine, with 10–12-fold increase in the yield of the soluble protein. The final protein conformation was identical by calorimetric criteria and CD and fluorescence spectroscopy to that obtained without additives, thus suggesting that VH-barnase structure does not depend on folding conditions. Folding of VH-barnase resulted in a single calorimetrically revealed folding unit, the so-called “calorimetric domain”, with conformation consistent with a molten globule that possessed well-defined secondary structure and compact tertiary conformation with partial exposure of hydrophobic patches and low thermodynamic stability. The unique feature of VH-barnase is that, despite the partially unfolded conformation and coupling into a single “calorimetric domain”, this immunofusion retained both the antigen-binding and RNase activities that belong to the two heterologous domains.