Fine-grained alignment of cryo-electron subtomograms based on MPI parallel optimization.

BACKGROUND: Cryo-electron tomography (Cryo-ET) is an imaging technique used to generate three-dimensional structures of cellular macromolecule complexes in their native environment. Due to developing cryo-electron microscopy technology, the image quality of three-dimensional reconstruction of cryo-electron tomography has greatly improved. However, cryo-ET images are characterized by low resolution, partial data loss and low signal-to-noise ratio (SNR). In order to tackle these challenges and improve resolution, a large number of subtomograms containing the same structure needs to be aligned and averaged. Existing methods for refining and aligning subtomograms are still highly time-consuming, requiring many computationally intensive processing steps (i.e. the rotations and translations of subtomograms in three-dimensional space). RESULTS: In this article, we propose a Stochastic Average Gradient (SAG) fine-grained alignment method for optimizing the sum of dissimilarity measure in real space. We introduce a Message Passing Interface (MPI) parallel programming model in order to explore further speedup. CONCLUSIONS: We compare our stochastic average gradient fine-grained alignment algorithm with two baseline methods, high-precision alignment and fast alignment. Our SAG fine-grained alignment algorithm is much faster than the two baseline methods. Results on simulated data of GroEL from the Protein Data Bank (PDB ID:1KP8) showed that our parallel SAG-based fine-grained alignment method could achieve close-to-optimal rigid transformations with higher precision than both high-precision alignment and fast alignment at a low SNR (SNR=0.003) with tilt angle range ±60∘ or ±40∘. For the experimental subtomograms data structures of GroEL and GroEL/GroES complexes, our parallel SAG-based fine-grained alignment can achieve higher precision and fewer iterations to converge than the two baseline methods.

[Ligand-Induced Reassembly of GroEL/ES Chaperone In Vitro: Visualization by Electron Microscopy].

The products of the reassembly reaction of tetradecameric two-ring quaternary structure of GroEL chaperonin under the pressure of its heptameric co-chaperonin GroES have been visualized by electron microscopy. It has been shown that one-ring heptameric oligomers of GroEL have been formed at the beginning (after ~5 min) of the reaction, while at the final stage of the reaction (after ~70 min), both one-ring heptamers in complex with one GroES and two-rings tetradecamers in complexes with one (asymmetrical complex) or two (symmetrical complex) GroES heptamers are present. The relationship between the data of light scattering, native electrophoresis, and electron microscopy obtained earlier has been discussed.

Bridging between NMA and Elastic Network Models: Preserving All-Atom Accuracy in Coarse-Grained Models.

Dynamics can provide deep insights into the functional mechanisms of proteins and protein complexes. For large protein complexes such as GroEL/GroES with more than 8,000 residues, obtaining a fine-grained all-atom description of its normal mode motions can be computationally prohibitive and is often unnecessary. For this reason, coarse-grained models have been used successfully. However, most existing coarse-grained models use extremely simple potentials to represent the interactions within the coarse-grained structures and as a result, the dynamics obtained for the coarse-grained structures may not always be fully realistic. There is a gap between the quality of the dynamics of the coarse-grained structures given by all-atom models and that by coarse-grained models. In this work, we resolve an important question in protein dynamics computations--how can we efficiently construct coarse-grained models whose description of the dynamics of the coarse-grained structures remains as accurate as that given by all-atom models? Our method takes advantage of the sparseness of the Hessian matrix and achieves a high efficiency with a novel iterative matrix projection approach. The result is highly significant since it can provide descriptions of normal mode motions at an all-atom level of accuracy even for the largest biomolecular complexes. The application of our method to GroEL/GroES offers new insights into the mechanism of this biologically important chaperonin, such as that the conformational transitions of this protein complex in its functional cycle are even more strongly connected to the first few lowest frequency modes than with other coarse-grained models.

Asp-52 in combination with Asp-398 plays a critical role in ATP hydrolysis of chaperonin GroEL.

The Escherichia coli chaperonin GroEL is a double-ring chaperone that assists protein folding with the aid of GroES and ATP. Asp-398 in GroEL is known as one of the critical residues on ATP hydrolysis because GroEL(D398A) mutant is deficient in ATP hydrolysis (<2% of the wild type) but not in ATP binding. In the archaeal Group II chaperonin, another aspartate residue, Asp-52 in the corresponding E. coli GroEL, in addition to Asp-398 is also important for ATP hydrolysis. We investigated the role of Asp-52 in GroEL and found that ATPase activity of GroEL(D52A) and GroEL(D52A/D398A) mutants was ∼ 20% and <0.01% of wild-type GroEL, respectively, indicating that Asp-52 in E. coli GroEL is also involved in the ATP hydrolysis. GroEL(D52A/D398A) formed a symmetric football-shaped GroEL-GroES complex in the presence of ATP, again confirming the importance of the symmetric complex during the GroEL ATPase cycle. Notably, the symmetric complex of GroEL(D52A/D398A) was extremely stable, with a half-time of ∼ 150 h (∼ 6 days), providing a good model to characterize the football-shaped complex.

Hsp10, Hsp70, and Hsp90 immunohistochemical levels change in ulcerative colitis after therapy.

Ulcerative colitis (UC) is a form of inflammatory bowel disease (IBD) characterized by damage of large bowel mucosa and frequent extra-intestinal autoimmune comorbidities. The role played in IBD pathogenesis by molecular chaperones known to interact with components of the immune system involved in inflammation is unclear. We previously demonstrated that mucosal Hsp60 decreases in UC patients treated with conventional therapies (mesalazine, probiotics), suggesting that this chaperonin could be a reliable biomarker useful for monitoring response to treatment, and that it might play a role in pathogenesis. In the present work we investigated three other heat shock protein/molecular chaperones: Hsp10, Hsp70, and Hsp90. We found that the levels of these proteins are increased in UC patients at the time of diagnosis and decrease after therapy, supporting the notion that these proteins deserve attention in the study of the mechanisms that promote the development and maintenance of IBD, and as biomarkers of this disease (e.g., to monitor response to treatment at the histological level).

Merging molecular electron microscopy and mass spectrometry by carbon film-assisted endoproteinase digestion.

Many fundamental processes in the cell are performed by complex macromolecular assemblies that comprise a large number of proteins. Numerous macromolecular assemblies are structurally rather fragile and may suffer during purification, resulting in the partial dissociation of the complexes. These limitations can be overcome by chemical fixation of the assemblies, and recently introduced protocols such as gradient fixation during ultracentrifugation (GraFix) offer advantages for the analysis of fragile macromolecular assemblies. The irreversible fixation, however, is thought to render macromolecular samples useless for studying their protein composition. We therefore developed a novel approach that possesses the advantages of fixation for structure determination by single particle electron microscopy while still allowing a correlative compositional analysis by mass spectrometry. In this method, which we call "electron microscopy carbon film-assisted digestion", macromolecular assemblies are chemically fixed and then adsorbed onto electron microscopical carbon films. Parallel, identically prepared specimens are then subjected to structural investigation by electron microscopy and proteomics analysis by mass spectrometry of the digested sample. As identical sample preparation protocols are used for electron microscopy and mass spectrometry, the results of both methods can directly be correlated. In addition, we demonstrate improved sensitivity and reproducibility of electron microscopy carbon film-assisted digestion as compared with standard protocols. We show that sample amounts of as low as 50 fmol are sufficient to obtain a comprehensive protein composition of two model complexes. We suggest our approach to be an optimization technique for the compositional analysis of macromolecules by mass spectrometry in general.

Quantitative analysis of cryo-EM density map segmentation by watershed and scale-space filtering, and fitting of structures by alignment to regions.

Cryo-electron microscopy produces 3D density maps of molecular machines, which consist of various molecular components such as proteins and RNA. Segmentation of individual components in such maps is a challenging task, and is mostly accomplished interactively. We present an approach based on the immersive watershed method and grouping of the resulting regions using progressively smoothed maps. The method requires only three parameters: the segmentation threshold, a smoothing step size, and the number of smoothing steps. We first apply the method to maps generated from molecular structures and use a quantitative metric to measure the segmentation accuracy. The method does not attain perfect accuracy, however it produces single or small groups of regions that roughly match individual proteins or subunits. We also present two methods for fitting of structures into density maps, based on aligning the structures with single regions or small groups of regions. The first method aligns centers and principal axes, whereas the second aligns centers and then rotates the structure to find the best fit. We describe both interactive and automated ways of using these two methods. Finally, we show segmentation and fitting results for several experimentally-obtained density maps.

Perturbation-based Markovian transmission model for probing allosteric dynamics of large macromolecular assembling: a study of GroEL-GroES.

Large macromolecular assemblies are often important for biological processes in cells. Allosteric communications between different parts of these molecular machines play critical roles in cellular signaling. Although studies of the topology and fluctuation dynamics of coarse-grained residue networks can yield important insights, they do not provide characterization of the time-dependent dynamic behavior of these macromolecular assemblies. Here we develop a novel approach called Perturbation-based Markovian Transmission (PMT) model to study globally the dynamic responses of the macromolecular assemblies. By monitoring simultaneous responses of all residues (>8,000) across many (>6) decades of time spanning from the initial perturbation until reaching equilibrium using a Krylov subspace projection method, we show that this approach can yield rich information. With criteria based on quantitative measurements of relaxation half-time, flow amplitude change, and oscillation dynamics, this approach can identify pivot residues that are important for macromolecular movement, messenger residues that are key to signal mediating, and anchor residues important for binding interactions. Based on a detailed analysis of the GroEL-GroES chaperone system, we found that our predictions have an accuracy of 71-84% judged by independent experimental studies reported in the literature. This approach is general and can be applied to other large macromolecular machineries such as the virus capsid and ribosomal complex.

GroEL-assisted protein folding: does it occur within the chaperonin inner cavity?

The folding of protein molecules in the GroEL inner cavity under the co-chaperonin GroES lid is widely accepted as a crucial event of GroEL-assisted protein folding. This review is focused on the data showing that GroEL-assisted protein folding may proceed out of the complex with the chaperonin. The models of GroEL-assisted protein folding assuming ligand-controlled dissociation of nonnative proteins from the GroEL surface and their folding in the bulk solution are also discussed.

Probing open conformation of GroEL rings by cross-linking reveals single and double open ring structures of GroEL in ADP and ATP.

Two heptamer rings of chaperonin GroEL undergo opening-closing conformational transition in the reaction cycle with the aid of GroES and ATP. We introduced Cys into the GroEL subunit at Ala-384 and Ser-509, which are very close between adjacent GroEL subunits in the open heptamer ring but far apart in the closed heptamer ring. The open ring-specific inter-subunit cross-linking between these Cys indicated that the number of rings in open conformation in GroEL was two in ATP (GroEL(OO)), one in ADP (GroEL(O)), and none in the absence of nucleotide. ADP showed an inhibitory effect on ATP-induced generation of GroEL(OO). The isolated GroEL(O) and GroEL(OO), which lost any bound nucleotide, could bind GroES to form a bullet-shaped 1:1 GroEL-GroES complex and a football-shaped 1:2 GroEL-GroES complex, respectively, even without the addition of any nucleotide. Substrate protein was unable to form a stable complex with GroEL(OO) and did not stimulate ATPase activity of GroEL. These results favor a model of the GroEL reaction cycle that includes a football complex as a critical intermediate.

Location and flexibility of the unique C-terminal tail of Aquifex aeolicus co-chaperonin protein 10 as derived by cryo-electron microscopy and biophysical techniques.

Co-chaperonin protein 10 (cpn10, GroES in Escherichia coli) is a ring-shaped heptameric protein that facilitates substrate folding when in complex with cpn60 (GroEL in E. coli). The cpn10 from the hyperthermophilic, ancient bacterium Aquifex aeolicus (Aacpn10) has a 25-residue C-terminal extension in each monomer not found in any other cpn10 protein. Earlier in vitro work has shown that this tail is not needed for heptamer assembly or protein function. Without the tail, however, the heptamers (Aacpn10del-25) readily aggregate into fibrillar stacked rings. To explain this phenomenon, we performed binding experiments with a peptide construct of the tail to establish its specificity for Aacpn10del-25 and used cryo-electron microscopy to determine the three-dimensional (3D) structure of the GroEL-Aacpn10-ADP complex at an 8-A resolution. We found that the GroEL-Aacpn10 structure is similar to the GroEL-GroES structure at this resolution, suggesting that Aacpn10 has molecular interactions with cpn60 similar to other cpn10s. The cryo-electron microscopy density map does not directly reveal the density of the Aacpn10 25-residue tail. However, the 3D statistical variance coefficient map computed from multiple 3D reconstructions with randomly selected particle images suggests that the tail is located at the Aacpn10 monomer-monomer interface and extends toward the cis-ring apical domain of GroEL. The tail at this location does not block the formation of a functional co-chaperonin/chaperonin complex but limits self-aggregation into linear fibrils at high temperatures. In addition, the 3D variance coefficient map identifies several regions inside the GroEL-Aacpn10 complex that have flexible conformations. This observation is in full agreement with the structural properties of an effective chaperonin.

Do chaperonins boost protein yields by accelerating folding or preventing aggregation?

The GroEL chaperonin has the ability to behave as an unfoldase, repeatedly denaturing proteins upon binding, which in turn can free them from kinetic traps and increase their folding rates. The complex formed by GroEL+GroES+ATP can also act as an infinite dilution cage, enclosing proteins within a protective container where they can fold without danger of aggregation. Controversy remains over which of these two properties is more critical to the GroEL/ES chaperonin's function. We probe the importance of the unfoldase nature of GroEL under conditions where aggregation is the predominant protein degradation pathway. We consider the effect of a hypothetical mutation to GroEL which increases the cycle frequency of GroEL/ES by increasing the rate of hydrolysis of GroEL-bound ATP. Using a simple kinetic model, we show that this modified chaperonin would be self-defeating: any potential reduction in folding time would be negated by an increase in time spent in the bulk, causing an increase in aggregation and a net decrease in protein folding yields.

Asymmetry of the GroEL-GroES complex under physiological conditions as revealed by small-angle x-ray scattering.

Despite the well-known functional importance of GroEL-GroES complex formation during the chaperonin cycle, the stoichiometry of the complex has not been clarified. The complex can occur either as an asymmetric 1:1 GroEL-GroES complex or as a symmetric 1:2 GroEL-GroES complex, although it remains uncertain which type is predominant under physiological conditions. To resolve this question, we studied the structure of the GroEL-GroES complex under physiological conditions by small-angle x-ray scattering, which is a powerful technique to directly observe the structure of the protein complex in solution. We evaluated molecular structural parameters, the radius of gyration and the maximum dimension of the complex, from the x-ray scattering patterns under various nucleotide conditions (3 mM ADP, 3 mM ATP gamma S, and 3 mM ATP in 10 mM MgCl(2) and 100 mM KCl) at three different temperatures (10 degrees C, 25 degrees C, and 37 degrees C). We then compared the experimentally observed scattering patterns with those calculated from the known x-ray crystallographic structures of the GroEL-GroES complex. The results clearly demonstrated that the asymmetric complex must be the major species stably present in solution under physiological conditions. On the other hand, in the presence of ATP (3 mM) and beryllium fluoride (10 mM NaF and 300 microM BeCl(2)), we observed the formation of a stable symmetric complex, suggesting the existence of a transiently formed symmetric complex during the chaperonin cycle.

Folding and assembly of co-chaperonin heptamer probed by forster resonance energy transfer.

The ring-shaped heptameric co-chaperonin protein 10 (cpn10) is one of few oligomeric beta-sheet proteins that unfold and disassemble reversibly in vitro. Here, we labeled human mitochondrial cpn10 with donor and acceptor dyes to obtain FRET signals. Cpn10 mixed in a 1:1:5 ratio of donor:acceptor:unlabeled monomers form heptamers that are active in an in vitro functional assay. Monomer-monomer affinity, as well as thermal and chemical stability, of the labeled cpn10 is similar to the unlabeled protein, demonstrating that the labels do not perturb the system. Using changes in FRET, we then probed for the first time cpn10 heptamer-monomer assembly/disassembly kinetics. Heptamer dissociation is very slow (1/k(diss) approximately 3h; 20 degrees C, pH 7) corresponding to an activation energy of approximately 50kJ/mol. Ring-ring mixing experiments reveal that cpn10 heptamer dissociation is rate limiting; subsequent associations events are faster. Kinetic inertness explains how cpn10 cycles on and off cpn60 as an intact heptamer in vivo.

Folding and assembly pathways of co-chaperonin proteins 10: Origin of bacterial thermostability.

To compare folding/assembly processes of heptameric co-chaperonin proteins 10 (cpn10) from different species and search for the origin of thermostability in hyper-thermostable Aquifex aeolicus cpn10 (Aacpn10), we have studied two bacterial variants-Aacpn10 and Escherichia coli cpn10 (GroES)-and compared the results to data on Homo sapiens cpn10 (hmcpn10). Equilibrium denaturation of GroES by urea, guanidine hydrochloride (GuHCl) and temperature results in coupled heptamer-to-monomer transitions in all cases. This is similar to the behavior of Aacpn10 but differs from hmcpn10 denaturation in urea. Time-resolved experiments reveal that GroES unfolds before heptamer dissociation, whereas refolding/reassembly begins with folding of individual monomers; these assemble in a slower step. The sequential folding/assembly mechanism for GroES is rather similar to that observed for Aacpn10 but contradicts the parallel paths of hmcpn10. We reveal that Aacpn10's stability profile is shifted upwards, broadened, and also moved horizontally to higher temperatures, as compared to that of GroES.

Fast-scanning atomic force microscopy reveals the ATP/ADP-dependent conformational changes of GroEL.

In order to fold non-native proteins, chaperonin GroEL undergoes numerous conformational changes and GroES binding in the ATP-dependent reaction cycle. We constructed the real-time three-dimensional-observation system at high resolution using a newly developed fast-scanning atomic force microscope. Using this system, we visualized the GroES binding to and dissociation from individual GroEL with a lifetime of 6 s (k=0.17 s(-1)). We also caught ATP/ADP-induced open-closed conformational changes of individual GroEL in the absence of qGroES and substrate proteins. Namely, the ATP/ADP-bound GroEL can change its conformation 'from closed to open' without additional ATP hydrolysis. Furthermore, the lifetime of open conformation in the presence of ADP ( approximately 1.0 s) was apparently lower than those of ATP and ATP-analogs (2-3 s), meaning that ADP-bound open-form is structurally less stable than ATP-bound open-form. These results indicate that GroEL has at least two distinct open-conformations in the presence of nucleotide; ATP-bound prehydrolysis open-form and ADP-bound open-form, and the ATP hydrolysis in open-form destabilizes its open-conformation and induces the 'from open to closed' conformational change of GroEL.

An expanded protein folding cage in the GroEL-gp31 complex.

Bacteriophage T4 produces a GroES analogue, gp31, which cooperates with the Escherichia coli GroEL to fold its major coat protein gp23. We have used cryo-electron microscopy and image processing to obtain three-dimensional structures of the E.coli chaperonin GroEL complexed with gp31, in the presence of both ATP and ADP. The GroEL-gp31-ADP map has a resolution of 8.2 A, which allows accurate fitting of the GroEL and gp31 crystal structures. Comparison of this fitted structure with that of the GroEL-GroES-ADP structure previously determined by cryo-electron microscopy shows that the folding cage is expanded. The enlarged volume for folding is consistent with the size of the bacteriophage coat protein gp23, which is the major substrate of GroEL-gp31 chaperonin complex. At 56 kDa, gp23 is close to the maximum size limit of a polypeptide that is thought to fit inside the GroEL-GroES folding cage.

Translocation boost protein-folding efficiency of double-barreled chaperonins.

Incorrect folding of proteins in living cells may lead to malfunctioning of the cell machinery. To prevent such cellular disasters from happening, all cells contain molecular chaperones that assist nonnative proteins in folding into the correct native structure. One of the most studied chaperone complexes is the GroEL-GroES complex. The GroEL part has a "double-barrel" structure, which consists of two cylindrical chambers joined at the bottom in a symmetrical fashion. The hydrophobic rim of one of the GroEL chambers captures nonnative proteins. The GroES part acts as a lid that temporarily closes the filled chamber during the folding process. Several capture-folding-release cycles are required before the nonnative protein reaches its native state. Here we report molecular simulations that suggest that translocation of the nonnative protein through the equatorial plane of the complex boosts the efficiency of the chaperonin action. If the target protein is correctly folded after translocation, it is released. However, if it is still nonnative, it is likely to remain trapped in the second chamber, which then closes to start a reverse translocation process. This shuttling back and forth continues until the protein is correctly folded. Our model provides a natural explanation for the prevalence of double-barreled chaperonins. Moreover, we argue that internal folding is both more efficient and safer than a scenario where partially refolded proteins escape from the complex before being recaptured.

Allosteric signaling of ATP hydrolysis in GroEL-GroES complexes.

The double-ring chaperonin GroEL and its lid-like cochaperonin GroES form asymmetric complexes that, in the ATP-bound state, mediate productive folding in a hydrophilic, GroES-encapsulated chamber, the so-called cis cavity. Upon ATP hydrolysis within the cis ring, the asymmetric complex becomes able to accept non-native polypeptides and ATP in the open, trans ring. Here we have examined the structural basis for this allosteric switch in activity by cryo-EM and single-particle image processing. ATP hydrolysis does not change the conformation of the cis ring, but its effects are transmitted through an inter-ring contact and cause domain rotations in the mobile trans ring. These rigid-body movements in the trans ring lead to disruption of its intra-ring contacts, expansion of the entire ring and opening of both the nucleotide pocket and the substrate-binding domains, admitting ATP and new substrate protein.

Role of the unique peptide tail in hyperthermostable Aquifex aeolicus cochaperonin protein 10.

All known cochaperonin protein 10 (cpn10) molecules are heptamers of seven identical subunits noncovalently linked by beta-strand interactions. Cpn10 from the deep-branching, hyperthermophilic bacterium Aquifex aeolicus (Aacpn10) shows high homology with mesophilic and other thermophilic cpn10 sequences, except for a 25-residue C-terminal extension not found in any other cpn10. Prior to atomic structure information, we here address the role of the tail by biophysical means. A tail-lacking variant (Aacpn10-del25) also adopts a heptameric structure in solution and exhibits nativelike substrate-refolding activity. Thermal and chemical perturbations of both Aacpn10 and Aacpn10-del25, probed by far-UV circular dichroism, demonstrate that both proteins have high thermodynamic stability. Heptamer-monomer dissociation midpoints were defined by isothermal titration calorimetry; at 25 degrees C, the values for Aacpn10 and Aacpn10-del25 are within 2-fold of each other and close to reported midpoints for mesophilic cpn10 proteins. In contrast, the monomer stabilities for the A. aeolicus proteins are significantly higher than those of mesophilic homologues at 30 degrees C; thus, heptamer thermophily is a result of more stable monomers. Electron microscopy data reveals that Aacpn10-del25 heptamers are prone to stack on top of each other forming chainlike molecules; the electrostatic surface pattern of a structural model can explain this behavior. Taken together, the unique tail in Aacpn10 is not required for heptamer structure, stability, or function; instead, it appears to be an ancient strategy to avoid cochaperonin aggregation at extreme temperatures.