Dimerisation and an increase in active site aromatic groups as adaptations to high temperatures: X-ray solution scattering and substrate-bound crystal structures of Rhodothermus marinus endoglucanase Cel12A

Cellulose, a polysaccharide consisting of beta-1,4-linked glucose, is the major component of plant cell walls and consequently one of the most abundant biopolymers on earth. Carbohydrate polymers such as cellulose are molecules with vast diversity in structure and function, and a multiplicity of hydrolases operating in concert are required for depolymerisation. The bacterium Rhodothermus marinus, isolated from shallow water marine hot springs, produces a number of carbohydrate-degrading enzymes including a family 12 cellulase Cel12A. The structure of R.marinus Cel12A in the ligand-free form (at 1.54 angstroms) and structures of RmCel12A after crystals were soaked in cellopentaose for two different lengths of time, have been determined. The shorter soaked complex revealed the conformation of unhydrolysed cellotetraose, while cellopentaose had been degraded more completely during the longer soak. Comparison of these structures with those of mesophilic family 12 cellulases in complex with inhibitors and substrate revealed that RmCel12A has a more extensive aromatic network in the active site cleft which ejects products after hydrolysis. The substrate structure confirms that during hydrolysis by family 12 cellulases glucose does not pass through a (2,5)B conformation. Small-angle X-ray scattering analysis of RmCel12A showed that the enzyme forms a loosely associated antiparallel dimer in solution, which may target the enzyme to the antiparallel polymer strands in cellulose.     

Protein rigidity and thermophilic adaptation

We probe the hypothesis of corresponding states, according to which homologues from mesophilic and thermophilic organisms are in corresponding states of similar rigidity and flexibility at their respective optimal temperatures. For this, the local distribution of flexible and rigid regions in 19 pairs of homologous proteins from meso- and thermophilic organisms is analyzed and related to activity characteristics of the enzymes by constraint network analysis (CNA). Two pairs of enzymes are considered in more detail: 3-isopropylmalate dehydrogenase and thermolysin-like protease. By comparing microscopic stability features of homologues with the help of stability maps, introduced for the first time, we show that adaptive mutations in enzymes from thermophilic organisms maintain the balance between overall rigidity, important for thermostability, and local flexibility, important for activity, at the appropriate working temperature. Thermophilic adaptation in general leads to an increase of structural rigidity but conserves the distribution of functionally important flexible regions between homologues. This finding provides direct evidence for the hypothesis of corresponding states. CNA thereby implicitly captures and unifies many different mechanisms that contribute to increased thermostability and to activity at high temperatures. This allows to qualitatively relate changes in the flexibility of active site regions, induced either by a temperature change or by the introduction of mutations, to experimentally observed losses of the enzyme function. As for applications, the results demonstrate that exploiting the principle of corresponding states not only allows for successful thermostability optimization but also for guiding experiments in order to improve enzyme activity in protein engineering.     

Cofilin increases the torsional flexibility and dynamics of actin filaments

We have measured the effects of cofilin on the conformation and dynamics of actin filaments labeled at Cys374 with erythrosin-iodoacetemide (ErIA), using time-resolved phosphorescence anisotropy (TPA). Cofilin quenches the phosphorescence intensity of actin-bound ErIA, indicating that binding changes the local environment of the probe. The cofilin concentration-dependence of the phosphorescence intensity is sigmoidal, consistent with cooperative actin filament binding. Model-independent analysis of the anisotropies indicates that cofilin increases the rates of the microsecond rotational motions of actin. In contrast to the reduction in phosphorescence intensity, the changes in the rates of rotational motions display non-nearest-neighbor cooperative interactions and saturate at substoichiometric cofilin binding densities. Detailed analysis of the TPA decays indicates that cofilin decreases the torsional rigidity (C) of actin, increasing the thermally driven root-mean-square torsional angle between adjacent filament subunits from approximately 4 degrees (C = 2.30 x 10(-27) Nm2 radian(-1)) to approximately 17 degrees (C = 0.13 x 10(-27) Nm2 radian(-1)) at 25 degrees C. We favor a mechanism in which cofilin binding shifts the equilibrium between thermal ErIA-actin filament conformers, and facilitates two distinct structural changes in actin. One is local in nature, which affects the structure of actin's C terminus and is likely to mediate nearest-neighbor cooperative binding and filament severing. The second is a change in the internal dynamics of actin, which displays non-nearest-neighbor cooperativity and increases the torsional flexibility of filaments. The long-range effects of cofilin on the torsional dynamics of actin may accelerate P(i) release from filaments and modulate interactions with other regulatory actin filament binding proteins.     

The structure of Rhodothermus marinus Cel12A,a highly thermostable family 12 endoglucanase,at 1.8 A resolution

Cellulose is one of the most abundant polysaccharides in nature and microorganisms have developed a comprehensive system for enzymatic breakdown of this ubiquitous carbon source, a subject of much interest in the biotechnology industry. Rhodothermus marinus produces a hyperthermostable cellulase, with a temperature optimum of more than 90 degrees C, the structure of which is presented here to 1.8 A resolution. The enzyme has been classified into glycoside hydrolase family 12; this is the first structure of a thermophilic member of this family to have been solved. The beta-jelly roll fold observed has identical topology to those of the two mesophilic members of the family whose structures have been elucidated previously. A Hepes buffer molecule bound in the active site may have triggered a conformational change to an active configuration as the two catalytic residues Glu124 and Glu207, together with dependent residues, are observed in a conformation similar to that seen in the structure of Streptomyces lividans CelB2 complexed with an inhibitor. The structural similarity between this cellulase and the mesophilic enzymes serves to highlight features that may be responsible for its thermostability, chiefly an increase in ion pair number and the considerable stabilisation of a mobile region seen in S. lividans CelB2. Additional aromatic residues in the active site region may also contribute to the difference in thermophilicity.     

Loop optimization of Trichoderma reesei endoglucanases for balancing the activity–stability trade-off through cross-strategy between machine learning and the B-factor analysis

Trichoderma reesei endoglucanases (EGs) have limited industrial applications due to its low thermostability and activity. Here, we aimed to improve the thermostability of EGs from T. reesei without reducing its activity counteracting the activity–stability trade-off. A cross-strategy combination of machine learning and B-factor analysis was used to predict beneficial amino acid substitution in EG loop optimization. Experimental validation showed single-site mutated EG concomitantly improved enzymatic activity and thermal properties by 17.21%–18.06% and 49.85%–62.90%, respectively, compared with wild-type EGs. Furthermore, the mechanism explained mutant variants had lower root mean square deviation values and a more stable overall structure than the wild type. According to this study, EGs loop optimization is crucial for balancing the activity–stability trade-off, which may provide new insights into how loop region function interacts with enzymatic characteristics. Moreover, the cross-strategy between machine learning and B-factor analysis improved superior enzyme activity–stability performance, which integrated structure-dependent and sequence-dependent information.     

The Limited Role of Nonnative Contacts in the Folding Pathways of a Lattice Protein

Models of protein energetics that neglect interactions between amino acids that are not adjacent in the native state, such as the Gō model, encode or underlie many influential ideas on protein folding. Implicit in this simplification is a crucial assumption that has never been critically evaluated in a broad context: Detailed mechanisms of protein folding are not biased by nonnative contacts, typically argued to be a consequence of sequence design and/or topology. Here we present, using computer simulations of a well-studied lattice heteropolymer model, the first systematic test of this oft-assumed correspondence over the statistically significant range of hundreds of thousands of amino acid sequences that fold to the same native structure. Contrary to previous conjectures, we find a multiplicity of folding mechanisms, suggesting that Gō-like models cannot be justified by considerations of topology alone. Instead, we find that the crucial factor in discriminating among topological pathways is the heterogeneity of native contact energies: The order in which native contacts accumulate is profoundly insensitive to omission of nonnative interactions, provided that native contact heterogeneity is retained. This robustness holds over a surprisingly wide range of folding rates for our designed sequences. Mirroring predictions based on the principle of minimum frustration, fast-folding sequences match their Gō-like counterparts in both topological mechanism and transit times. Less optimized sequences dwell much longer in the unfolded state and/or off-pathway intermediates than do Gō-like models. For dynamics that bridge unfolded and unfolded states, however, even slow folders exhibit topological mechanisms and transit times nearly identical with those of their Gō-like counterparts. Our results do not imply a direct correspondence between folding trajectories of Gō-like models and those of real proteins, but they do help to clarify key topological and energetic assumptions that are commonly used to justify such caricatures.     

Properties of a recombinant human hemoglobin with aspartic acid 99(beta), an important intersubunit contact site, substituted by lysine.

Site-directed mutagenesis of an important subunit contact site, Asp-99(beta), by a Lys residue (D99K(beta)) was proven by sequencing the entire beta-globin gene and the mutant tryptic peptide. Oxygen equilibrium curves of the mutant hemoglobin (Hb) (2-15 mM in heme) indicated that it had an increased oxygen affinity and a lowered but significant amount of cooperativity compared to native HbA. However, in contrast to normal HbA, oxygen binding of the recombinant mutant Hb was only marginally affected by the allosteric regulators 2,3-diphosphoglycerate or inositol hexaphosphate and was not at all responsive to chloride. The efficiency of oxygen binding by HbA in the presence of allosteric regulators was limited by the mutant Hb. At concentrations of 0.2 mM or lower in heme, the mutant D99K(beta) Hb was predominantly a dimer as demonstrated by gel filtration, haptoglobin binding, fluorescence quenching, and light scattering. The purified dimeric recombinant Hb mutant exists in 2 forms that are separable on isoelectric focusing by about 0.1 pH unit, in contrast to tetrameric hemoglobin, which shows 1 band. These mutant forms, which were present in a ratio of 60:40, had the same masses for their heme and globin moieties as determined by mass spectrometry. The elution positions of the alpha- and beta-globin subunits on HPLC were identical. Circular dichroism studies showed that one form of the mutant Hb had a negative ellipticity at 410 nm and the other had positive ellipticity at this wavelength. The findings suggest that the 2 D99K(beta) recombinant mutant forms have differences in their heme-protein environments.     

蛋白质折叠速度与其拓扑结构关系的研究

以6种不同的方式来定义蛋白质内存在的接触,进而运用分子动力学模拟等不同方法,对10个小蛋白进行分析,研究了不同的接触定义及不同的拓扑参数计算方法下,蛋白质的折叠速度与其拓扑参数的关系.结果表明,用含主链重原子的方式定义接触,所计算的拓扑参数与蛋白质折叠速度的相关性较好;用含侧链原子的方式定义接触,得到的拓扑参数与β型蛋白质的折叠速度的相关性较好.对不同的蛋白质,其拓扑结构与相应折叠速度间的相关程度不同。     

Cooperativity in two-state protein folding kinetics

We present a solvable model that predicts the folding kinetics of two-state proteins from their native structures. The model is based on conditional chain entropies. It assumes that folding processes are dominated by small-loop closure events that can be inferred from native structures. For CI2, the src SH3 domain, TNfn3, and protein L, the model reproduces two-state kinetics, and it predicts well the average Phi-values for secondary structures. The barrier to folding is the formation of predominantly local structures such as helices and hairpins, which are needed to bring nonlocal pairs of amino acids into contact.