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.

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.

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.

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.

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.

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.

Functional significance of symmetrical versus asymmetrical GroEL-GroES chaperonin complexes.

The Escherichia coli chaperonin GroEL and its regulator GroES are thought to mediate adenosine triphosphate-dependent protein folding as an asymmetrical complex, with substrate protein bound within the GroEL cylinder. In contrast, a symmetrical complex formed between one GroEL and two GroES oligomers, with substrate protein binding to the outer surface of GroEL, was recently proposed to be the functional chaperonin unit. Electron microscopic and biochemical analyses have now shown that unphysiologically high magnesium concentrations and increased pH are required to assemble symmetrical complexes, the formation of which precludes the association of unfolded polypeptide. Thus, the functional significance of GroEL:(GroES)2 particles remains to be demonstrated.

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.

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.

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.

Role of the amino terminal domain in GroES oligomerization.

Digestions of the GroES oligomer with trypsin, chymotrypsin and Glu-C protease from Staphylococcus aureus V8 (V8) have helped to locate three regions in the GroES sequence that are sensitive to limited proteolysis and have provided information of the GroES domains involved in monomer-monomer and GroEL interaction. The removal of the first 20 or 27 amino acids of the N-terminal region of each GroES monomer by trypsin or chymotrypsin respectively, abolish the oligomerization of the GroES complex and its binding to GroEL. The V8-treatment of GroES promotes the breakage of the peptide bond between Glu18 and Thr19 but not the liberation of the N-terminal fragment from the GroES oligomer, which is capable of forming with GroEL a complex active in protein folding. It is deduced from these results that the N-terminal region of the GroES monomer is involved in monomer-monomer interaction, providing experimental evidence that relates some biochemical properties of GroES with its three-dimensional structure at atomic resolution.

Structural and functional imaging with carbon nanotube AFM probes.

Atomic force microscopy (AFM) has great potential as a tool for structural biology, a field in which there is increasing demand to characterize larger and more complex biomolecular systems. However, the poorly characterized silicon and silicon nitride probe tips currently employed in AFM limit its biological applications. Carbon nanotubes represent ideal AFM tip materials due to their small diameter, high aspect ratio, large Young's modulus, mechanical robustness, well-defined structure, and unique chemical properties. Nanotube probes were first fabricated by manual assembly, but more recent methods based on chemical vapor deposition provide higher resolution probes and are geared towards mass production, including recent developments that enable quantitative preparation of individual single-walled carbon nanotube tips [J. Phys. Chem. B 105 (2001) 743]. The high-resolution imaging capabilities of these nanotube AFM probes have been demonstrated on gold nanoparticles and well-characterized biomolecules such as IgG and GroES. Using the nanotube probes, new biological structures have been investigated in the areas of amyloid-beta protein aggregation and chromatin remodeling, and new biotechnologies have been developed such as AFM-based haplotyping. In addition to measuring topography, chemically functionalized AFM probes can measure the spatial arrangement of chemical functional groups in a sample. However, standard silicon and silicon nitride tips, once functionalized, do not yield sufficient resolution to allow combined structural and functional imaging of biomolecules. The unique end-group chemistry of carbon nanotubes, which can be arbitrarily modified by established chemical methods, has been exploited for chemical force microscopy, allowing single-molecule measurements with well-defined functionalized tips.

Solubilization and delivery by GroEL of megadalton complexes of the lambda holin.

GroEL can solubilize membrane proteins by binding them in its hydrophobic cavity when detergent is removed by dialysis. The best-studied example is bacteriorhodopsin, which can bind in the GroEL chaperonin at two molecules per tetradecamer. Applying this approach to the holin and antiholin proteins of phage lambda, we find that both proteins are solubilized by GroEL, in an ATP-sensitive mode, but to vastly different extents. The antiholin product, S107, saturates the chaperonin at six molecules per tetradecameric complex, whereas the holin, S105, which is missing the two N-terminal residues of S107, forms a hyper-solubilization complex with up to 350 holin molecules per GroEL, or approximately 4 MDa of protein per 0.8 MDa tetradecamer. Gel filtration chromatography and immunoprecipitation experiments confirmed the existence of complexes of the predicted masses for both S105 and S107 solubilization. For S105, negatively stained electron microscopic images show structures consistent with protein shells of the holin assembled around the chaperonin tetradecamer. Importantly, S105 can be delivered rapidly and efficiently to artificial liposomes from these complexes. In these delivery experiments, the holin exhibits efficient membrane-permeabilizing activity. The S107 antiholin can block formation of the hypersolubilization complexes, suggesting that their formation is related to an oligomerization step intrinsic to holin function.

High resolution surface structure of E. coli GroES oligomer by atomic force microscopy.

Using atomic force microscopy (AFM) in aqueous solution, we show that the surface structure of the oligomeric GroES can be obtained up to 10 angstroms resolution. The seven subunits of the heptamer were well resolved without image averaging. The overall dimension of the GroES heptamer was 8.4 +/- 0.4 nm in diameter and 3.0 +/- 0.3 nm high. However, the AFM images further suggest that there is a central protrusion of 0.8 +/- 0.2 nm high and 4.5 +/- 0.4 nm in diameter on one side of GroES which displays a profound seven-fold symmetry. It was found that GroEL could not bind to the adsorbed GroES in the presence of AMP-PNP and Mg2+, suggesting that the side of GroES with the central protrusion faces away from the GroEL lumen, because only one side of GroES was observed under these conditions. Based on the results from both electron and atomic force microscopy, a surface model for the GroES is proposed.

Symmetric GroEL-GroES complexes can contain substrate simultaneously in both GroEL rings.

Incubation of rhodanese with hche aperonins GroEL and GroES (1:2 GroEL14:GroES7 molar ratio) under functional and steady state conditions for ATP leads to the formation of a high proportion of rhodanese-bound symmetric complexes (GroEL14(GroES7)2), as revealed by native electrophoresis. Aliquots of such samples were observed under the electron microscope, and the symmetric particles were classified using neuronal networks and multivariate statistical analysis. Three different populations of symmetric particles were obtained which contained substrate in none, one or both GroEL cavities, respectively. The presence of substrate in the symmetric complexes under functional conditions supports their role as active intermediates in the protein folding cycle. These results also suggest that symmetric GroEL-GroES complexes can use both rings simultaneously for folding, probably increasing the efficiency of the reaction.

Equatorial split of holo-chaperonin from Thermus thermophilus by ATP and K+.

Holo-chaperonin molecule from Thermus thermophilus is a bullet-shaped particle whose cylinder part and round top are composed of two stacked rings of the cpn60 heptamer and a single ring of the cpn10 heptamer, respectively. We found that it splits at the plane between two cpn60 rings into two halves under physiological conditions, that is, in the presence of ATP (but not AMP-PNP, ADP) + K+ (but not Na+) at 60 degrees C. This equatorial split could be functionally important although it has not been considered in any current mechanistic model of chaperonin functioning.