Exploring the structural dynamics of the E.coli chaperonin GroEL using translation-libration-screw crystallographic refinement of intermediate states

Mushroom Tyrosinase: Six Isoenzymes Catalyzing Distinct Reactions

"Mushroom tyrosinase" from the common button mushroom is the most frequently used source of tyrosinase activity, both for basic and applied research. Here, the complete tyrosinase family from Agaricus bisporus var. bisporus (abPPO1-6) was cloned from mRNA and expressed heterologously using a single protocol. All six isoenzymes accept a wide range of phenolic and catecholic substrates, but display pronounced differences in their specificity and enzymatic reaction rate. AbPPO3 ignores γ-l-glutaminyl-4-hydroxybenzene (GHB), a natural phenol present in mM concentrations in A. bisporus, while AbPPO4 processes 100 μM GHB at 4-times the rate of the catechol l-DOPA. All six AbPPOs are biochemically distinct enzymes fit for different roles in the fungal life cycle, which challenges the traditional concept of isoenzymes as catalyzing the same physiological reaction and varying only in secondary properties. Transferring this approach to other enzymes and organisms will greatly stimulate both the study of the in vivo function(s) of enzymes and the application of these highly efficient catalysts.

MDFF_NM: Improved Molecular Dynamics Flexible Fitting into Cryo-EM Density Maps with a Multireplica Normal Mode-Based Search

Molecular Dynamics Flexible Fitting (MDFF) is a widely used tool to refine high-resolution structures into cryo-EM density maps. Despite many successful applications, MDFF is still limited by its high computational cost, overfitting, accuracy, and performance issues due to entrapment within wrong local minima. Modern ensemble-based MDFF tools have generated promising results in the past decade. In line with these studies, we present MDFF_NM, a stochastic hybrid flexible fitting algorithm combining Normal Mode Analysis (NMA) and simulation-based flexible fitting. Initial tests reveal that, besides accelerating the fitting process, MDFF_NM increases the diversity of fitting routes leading to the target, uncovering ensembles of conformations in closer agreement with experimental data. The potential integration of MDFF_NM with other existing methods and integrative modeling pipelines is also discussed.

Asymmetric apical domain states of mitochondrial Hsp60 coordinate substrate engagement and chaperonin assembly

The mitochondrial chaperonin, mitochondrial heat shock protein 60 (mtHsp60), promotes the folding of newly imported and transiently misfolded proteins in the mitochondrial matrix, assisted by its co-chaperone mtHsp10. Despite its essential role in mitochondrial proteostasis, structural insights into how this chaperonin progresses through its ATP-dependent client folding cycle are not clear. Here, we determined cryo-EM structures of a hyperstable disease-associated human mtHsp60 mutant, V72I. Client density is identified in three distinct states, revealing interactions with the mtHsp60 apical domains and C termini that coordinate client positioning in the folding chamber. We further identify an asymmetric arrangement of the apical domains in the ATP state, in which an alternating up/down configuration positions interaction surfaces for simultaneous recruitment of mtHsp10 and client retention. Client is then fully encapsulated in mtHsp60-10, revealing prominent contacts at two discrete sites that potentially support maturation. These results identify distinct roles for the apical domains in coordinating client capture and progression through the chaperone cycle, supporting a conserved mechanism of group I chaperonin function.

Native Capillary Electrophoresis-Mass Spectrometry of Near 1 MDa Non-Covalent GroEL/GroES/Substrate Protein Complexes

Protein complexes are essential for proteins' folding and biological function. Currently, native analysis of large multimeric protein complexes remains challenging. Structural biology techniques are time-consuming and often cannot monitor the proteins' dynamics in solution. Here, a capillary electrophoresis-mass spectrometry (CE-MS) method is reported to characterize, under near-physiological conditions, the conformational rearrangements of ∽1 MDa GroEL upon complexation with binding partners involved in a protein folding cycle. The developed CE-MS method is fast (30 min per run), highly sensitive (low-amol level), and requires ∽10 000-fold fewer samples compared to biochemical/biophysical techniques. The method successfully separates GroEL14 (∽800 kDa), GroEL7 (∽400 kDa), GroES7 (∽73 kDa), and NanA4 (∽130 kDa) oligomers. The non-covalent binding of natural substrate proteins with GroEL14 can be detected and quantified. The technique allows monitoring of GroEL14 conformational changes upon complexation with (ATPγS)4-14 and GroES7 (∽876 kDa). Native CE-pseudo-MS3 analyses of wild-type (WT) GroEL and two GroEL mutants result in up to 60% sequence coverage and highlight subtle structural differences between WT and mutated GroEL. The presented results demonstrate the superior CE-MS performance for multimeric complexes' characterization versus direct infusion ESI-MS. This study shows the CE-MS potential to provide information on binding stoichiometry and kinetics for various protein complexes.

Observation of a single protein by ultrafast X-ray diffraction

The idea of using ultrashort X-ray pulses to obtain images of single proteins frozen in time has fascinated and inspired many. It was one of the arguments for building X-ray free-electron lasers. According to theory, the extremely intense pulses provide sufficient signal to dispense with using crystals as an amplifier, and the ultrashort pulse duration permits capturing the diffraction data before the sample inevitably explodes. This was first demonstrated on biological samples a decade ago on the giant mimivirus. Since then, a large collaboration has been pushing the limit of the smallest sample that can be imaged. The ability to capture snapshots on the timescale of atomic vibrations, while keeping the sample at room temperature, may allow probing the entire conformational phase space of macromolecules. Here we show the first observation of an X-ray diffraction pattern from a single protein, that of Escherichia coli GroEL which at 14 nm in diameter is the smallest biological sample ever imaged by X-rays, and demonstrate that the concept of diffraction before destruction extends to single proteins. From the pattern, it is possible to determine the approximate orientation of the protein. Our experiment demonstrates the feasibility of ultrafast imaging of single proteins, opening the way to single-molecule time-resolved studies on the femtosecond timescale.

Structural basis of substrate progression through the bacterial chaperonin cycle

The bacterial chaperonin GroEL-GroES promotes protein folding through ATP-regulated cycles of substrate protein binding, encapsulation, and release. Here, we have used cryoEM to determine structures of GroEL, GroEL-ADP·BeF3, and GroEL-ADP·AlF3-GroES all complexed with the model substrate Rubisco. Our structures provide a series of snapshots that show how the conformation and interactions of non-native Rubisco change as it proceeds through the GroEL-GroES reaction cycle. We observe specific charged and hydrophobic GroEL residues forming strong initial contacts with non-native Rubisco. Binding of ATP or ADP·BeF3 to GroEL-Rubisco results in the formation of an intermediate GroEL complex displaying striking asymmetry in the ATP/ADP·BeF3-bound ring. In this ring, four GroEL subunits bind Rubisco and the other three are in the GroES-accepting conformation, suggesting how GroEL can recruit GroES without releasing bound substrate. Our cryoEM structures of stalled GroEL-ADP·AlF3-Rubisco-GroES complexes show Rubisco folding intermediates interacting with GroEL-GroES via different sets of residues.

Cryo-EM of a heterogeneous biochemical fraction elucidates multiple protein complexes from a multicellular thermophilic eukaryote

Biomolecular complexes and their interactions govern cellular structure and function. Understanding their architecture is a prerequisite for dissecting the cell's inner workings, but their higher-order assembly is often transient and challenging for structural analysis. Here, we performed cryo-EM on a single, highly heterogeneous biochemical fraction derived from Chaetomium thermophilum cell extracts to visualize the biomolecular content of the multicellular eukaryote. After cryo-EM single-particle image processing, results showed that a simultaneous three-dimensional structural characterization of multiple chemically diverse biomacromolecules is feasible. Namely, the thermophilic, eukaryotic complexes of (a) ATP citrate-lyase, (b) Hsp90, (c) 20S proteasome, (d) Hsp60 and (e) UDP-glucose pyrophosphorylase were characterized. In total, all five complexes have been structurally dissected in a thermophilic eukaryote in a total imaged sample area of 190.64 μm2, and two, in particular, 20S proteasome and Hsp60, exhibit side-chain resolution features. The C. thermophilum Hsp60 near-atomic model was resolved at 3.46 Å (FSC = 0.143) and shows a hinge-like conformational change of its equatorial domain, highly similar to the one previously shown for its bacterial orthologue, GroEL. This work demonstrates that cryo-EM of cell extracts will greatly accelerate the structural analysis of cellular complexes and provide unprecedented opportunities to annotate architectures of biomolecules in a holistic approach.

From Microstates to Macrostates in the Conformational Dynamics of GroEL: A Single-Molecule Förster Resonance Energy Transfer Study

The chaperonin GroEL is a multisubunit molecular machine that assists in protein folding in the Escherichia coli cytosol. Past studies have shown that GroEL undergoes large allosteric conformational changes during its reaction cycle. Here, we report single-molecule Förster resonance energy transfer measurements that directly probe the conformational transitions of one subunit within GroEL and its single-ring variant under equilibrium conditions. We find that four microstates span the conformational manifold of the protein and interconvert on the submillisecond time scale. A unique set of relative populations of these microstates, termed a macrostate, is obtained by varying solution conditions, e.g., adding different nucleotides or the cochaperone GroES. Strikingly, ATP titration studies demonstrate that the partition between the apo and ATP-ligated conformational macrostates traces a sigmoidal response with a Hill coefficient similar to that obtained in bulk experiments of ATP hydrolysis. These coinciding results from bulk measurements for an entire ring and single-molecule measurements for a single subunit provide new evidence for the concerted allosteric transition of all seven subunits.

Pathway and mechanism of tubulin folding mediated by TRiC/CCT along its ATPase cycle revealed using cryo-EM

The eukaryotic chaperonin TRiC/CCT assists the folding of about 10% of cytosolic proteins through an ATP-driven conformational cycle, and the essential cytoskeleton protein tubulin is the obligate substrate of TRiC. Here, we present an ensemble of cryo-EM structures of endogenous human TRiC throughout its ATPase cycle, with three of them revealing endogenously engaged tubulin in different folding stages. The open-state TRiC-tubulin-S1 and -S2 maps show extra density corresponding to tubulin in the cis-ring chamber of TRiC. Our structural and XL-MS analyses suggest a gradual upward translocation and stabilization of tubulin within the TRiC chamber accompanying TRiC ring closure. In the closed TRiC-tubulin-S3 map, we capture a near-natively folded tubulin-with the tubulin engaging through its N and C domains mainly with the A and I domains of the CCT3/6/8 subunits through electrostatic and hydrophilic interactions. Moreover, we also show the potential role of TRiC C-terminal tails in substrate stabilization and folding. Our study delineates the pathway and molecular mechanism of TRiC-mediated folding of tubulin along the ATPase cycle of TRiC, and may also inform the design of therapeutic agents targeting TRiC-tubulin interactions.

A molecular device for the redox quality control of GroEL/ES substrates

Hsp60 chaperonins and their Hsp10 cofactors assist protein folding in all living cells, constituting the paradigmatic example of molecular chaperones. Despite extensive investigations of their structure and mechanism, crucial questions regarding how these chaperonins promote folding remain unsolved. Here, we report that the bacterial Hsp60 chaperonin GroEL forms a stable, functionally relevant complex with the chaperedoxin CnoX, a protein combining a chaperone and a redox function. Binding of GroES (Hsp10 cofactor) to GroEL induces CnoX release. Cryoelectron microscopy provided crucial structural information on the GroEL-CnoX complex, showing that CnoX binds GroEL outside the substrate-binding site via a highly conserved C-terminal α-helix. Furthermore, we identified complexes in which CnoX, bound to GroEL, forms mixed disulfides with GroEL substrates, indicating that CnoX likely functions as a redox quality-control plugin for GroEL. Proteins sharing structural features with CnoX exist in eukaryotes, suggesting that Hsp60 molecular plugins have been conserved through evolution.

Cryo-EM Analyses Permit Visualization of Structural Polymorphism of Biological Macromolecules

The functions of biological macromolecules are often associated with conformational malleability of the structures. This phenomenon of chemically identical molecules with different structures is coined structural polymorphism. Conventionally, structural polymorphism is observed directly by structural determination at the density map level from X-ray crystal diffraction. Although crystallography approach can report the conformation of a macromolecule with the position of each atom accurately defined in it, the exploration of structural polymorphism and interpreting biological function in terms of crystal structures is largely constrained by the crystal packing. An alternative approach to studying the macromolecule of interest in solution is thus desirable. With the advancement of instrumentation and computational methods for image analysis and reconstruction, cryo-electron microscope (cryo-EM) has been transformed to be able to produce “in solution” structures of macromolecules routinely with resolutions comparable to crystallography but without the need of crystals. Since the sample preparation of single-particle cryo-EM allows for all forms co-existing in solution to be simultaneously frozen, the image data contain rich information as to structural polymorphism. The ensemble of structure information can be subsequently disentangled through three-dimensional (3D) classification analyses. In this review, we highlight important examples of protein structural polymorphism in relation to allostery, subunit cooperativity and function plasticity recently revealed by cryo-EM analyses, and review recent developments in 3D classification algorithms including neural network/deep learning approaches that would enable cryo-EM analyese in this regard. Finally, we brief the frontier of cryo-EM structure determination of RNA molecules where resolving the structural polymorphism is at dawn.

Comparative Natural History of Visual Function From Patients With Biallelic Variants in BBS1 and BBS10

Purpose

The purpose of this study was to compare the natural history of visual function change in cohorts of patients affected with retinal degeneration due to biallelic variants in Bardet-Biedl syndrome genes: BBS1 and BBS10.

Methods

Patients were recruited from nine academic centers from six countries (Belgium, Canada, France, New Zealand, Switzerland, and the United States). Inclusion criteria were: (1) female or male patients with a clinical diagnosis of retinal dystrophy, (2) biallelic disease-causing variants in BBS1 or BBS10, and (3) measures of visual function for at least one visit. Retrospective data collected included genotypes, age, onset of symptoms, and best corrected visual acuity (VA). When possible, data on refractive error, fundus images and autofluorescence (FAF), optical coherence tomography (OCT), Goldmann kinetic perimetry (VF), electroretinography (ERG), and the systemic phenotype were collected.

Results

Sixty-seven individuals had variants in BBS1 (n = 38; 20 female patients and 18 male patients); or BBS10 (n = 29; 14 female patients and 15 male patients). Missense variants were the most common type of variants for patients with BBS1, whereas frameshift variants were most common for BBS10. When ERGs were recordable, rod-cone dystrophy (RCD) was observed in 82% (23/28) of patients with BBS1 and 73% (8/11) of patients with BBS10; cone-rod dystrophy (CORD) was seen in 18% of patients with BBS1 only, and cone dystrophy (COD) was only seen in 3 patients with BBS10 (27%). ERGs were nondetectable earlier in patients with BBS10 than in patients with BBS1. Similarly, VA and VF declined more rapidly in patients with BBS10 compared to patients with BBS1.

Conclusions

Retinal degeneration appears earlier and is more severe in BBS10 cases as compared to those with BBS1 variants. The course of change of visual function appears to relate to genetic subtypes of BBS.

Crystal structure of P. falciparum Cpn60 bound to ATP reveals an open dynamic conformation before substrate binding

Plasmodium falciparum harbors group 1 and group 2 chaperonin systems to mediate the folding of cellular proteins in different cellular locations. Two distinct group 1 chaperonins operate in the organelles of mitochondria and apicoplasts, while group 2 chaperonins function in the cytosol. No structural information has been reported for any chaperonin from plasmodium. In this study, we describe the crystal structure of a double heptameric ring Plasmodium falciparum mitochondrial chaperonin 60 (Cpn60) bound with ATP, which differs significantly from any known crystal structure of chaperonin 60. The structure likely represents a unique intermediate state during conformational conversion from the closed state to the opened state. Three of the seven apical domains are highly dynamic while the equatorial domains form a stable ring. The structure implies large movements of the apical domain in the solution play a role in nucleotide-dependent regulation of substrate binding and folding. A unique 26–27 residue insertion in the equatorial domain of Plasmodium falciparum mitochondrial chaperonin greatly increases both inter-ring and intra-ring subunit–subunit interactions. The present structure provides new insights into the mechanism of Cpn60 in chaperonin assembly and function.

Structural analysis of the Sulfolobus solfataricus TF55β chaperonin by cryo-electron microscopy

Chaperonins are biomolecular complexes that assist in protein folding. Thermophilic factor 55 (TF55) is a group II chaperonin found in the archaeal genus Sulfolobus that has α, β and γ subunits. Using cryo-electron microscopy, structures of the β-only complex of S. solfataricus TF55 (TF55β) were determined to 3.6-4.2 Å resolution. The structures of the TF55β complexes formed in the presence of ADP or ATP highlighted an open state in which nucleotide exchange can occur before progressing in the refolding cycle.

Nanoscale integration of porphyrin in GroEL protein cage: Photophysical and photochemical investigation

In this paper, we introduce a new type of functional, supramolecular porphyrin conjugate created using the bacterial GroEL protein cage based on non-specific hydrophobic interaction. The synthesis, structure and property of the porphyrin conjugate were characterized by dynamic light scattering, UV-vis spectroscopy and fluorescence spectroscopy. We observed that the model zinc-tetraphenylporphyrin (Zn-TPP) with high hydrophobicity can be well-dispersed in aqueous solutions with the aid of GroEL open chamber, which is known to be a favorable nanocompartment for aggregation-prone molecules. The maximal encapsulation efficiency of Zn-TPP in GroEL was determined to be ~98%. It is further seen that the constructed double Zn-TPP-GroEL complex exhibited good photocatalytic activity in the model reactions of the production of singlet oxygen and the reduction of methyl viologen under illumination with visible light. Moreover, we found that GroEL can significantly improve the photostability of Zn-TPP molecules as a result of nanoscale assembly within its hydrophobic chamber. Hence enhanced water solubility and photostability of Zn-TPP, which are considered as the first two hurdles for the wide usage of porphyrins, were achieved simultaneously by the development of GroEL cage as a building block. Supramolecular nanostructures formed from porphyrins (or related molecules) and GroEL for photocatalysis would greatly simplify applications of such structures.

First-Principles Collision Cross Section Measurements of Large Proteins and Protein Complexes

Rotationally averaged collision cross section (CCS) values for a series of proteins and protein complexes ranging in size from 8.6 to 810 kDa are reported. The CCSs were obtained using a native electrospray ionization drift tube ion mobility-Orbitrap mass spectrometer specifically designed to enhance sensitivity while having high-resolution ion mobility and mass capabilities. Periodic focusing (PF)-drift tube (DT)-ion mobility (IM) provides first-principles determination of the CCS of large biomolecules that can then be used as CCS calibrants. The experimental, first-principles CCS values are compared to previously reported experimentally determined and computationally calculated CCS using projected superposition approximation (PSA), the Ion Mobility Projection Approximation Calculation Tool (IMPACT), and Collidoscope. Experimental CCS values are generally in agreement with previously reported CCSs, with values falling within ∼5.5%. In addition, an ion mobility resolution (CCS centroid divided by CCS fwhm) of ∼60 is obtained for pyruvate kinase (MW ∼ 233 kDa); however, ion mobility resolution for bovine serum albumin (MW ∼ 68 kDa) is less than ∼20, which arises from sample impurities and underscores the importance of sample quality. The high resolution afforded by the ion mobility-Orbitrap mass analyzer provides new opportunities to understand the intricate details of protein complexes such as the impact of post-translational modifications (PTMs), stoichiometry, and conformational changes induced by ligand binding.

Chaperonin-assisted protein folding: a chronologue

This chronologue seeks to document the discovery and development of an understanding of oligomeric ring protein assemblies known as chaperonins that assist protein folding in the cell. It provides detail regarding genetic, physiologic, biochemical, and biophysical studies of these ATP-utilizing machines from both in vivo and in vitro observations. The chronologue is organized into various topics of physiology and mechanism, for each of which a chronologic order is generally followed. The text is liberally illustrated to provide firsthand inspection of the key pieces of experimental data that propelled this field. Because of the length and depth of this piece, the use of the outline as a guide for selected reading is encouraged, but it should also be of help in pursuing the text in direct order.

Homology Searches Using Supersecondary Structure Code

Supersecondary structure code (SSSC), which is represented as the combination of α-helix-type (SSSC: H), β-sheet-type (SSSC: S), the other (SSSC: T), and disorder residue or C-terminal (SSSC: D) patterns, has been produced by the developed concept of Ramachandran plot, in addition, with the ω angle and with the specification of positions of torsion angles in a protein by the registration of codes for torsion angles of each amino acid peptide unit, derived from the fuzzy search of structural code homology using the template patterns 3a5c4a (SSSC: H) and 6c4a4a (SSSC: S) with conformational codes. The DSSP (Dictionary of Secondary Structure in Proteins) method assigns the secondary structure including hydrogen bond well. In contrast, supersecondary structure code is very sensitive to the supersecondary structures of proteins. In this chapter, the protocol of homology search methods, the sequence alignment using supersecondary structure code, the assignment of supersecondary structure code T, the fuzzy search using supersecondary structure code, and the exact search using supersecondary structure code are described. Supersecondary structure code is variable with the conformational change. If possible, many Protein Data Bank (PDB) data of similar main chains of proteins should be used for the homology searches. The thorough check of SSSC sequences is also useful to reveal the role of target pattern.

Inter and intra-subunit interactions at the subunit interface of chaperonin GroEL are essential for its stability and assembly

Chaperonin GroEL helps in the folding of substrate proteins under normal and stress conditions. Although it remains stable and functional during stress conditions, the quantitative estimation of stability parameters and the specific amino-acid residues playing a role in its stability are not known in sufficient detail. The reason for poor understanding is its large size, multimeric nature, and irreversible unfolding process. The X-ray crystal structure reveals that equatorial domain forms almost all intra and inter-subunit interactions for assembly of GroEL. Considering all these facts, we adopted alternate strategies to use monomeric GroEL, native GroEL and equatorial domain mutants (GroEL_K4E/GroEL_D523K/GroEL_D473C) to study the assembly and stability of GroEL. Loss of inter-subunit interaction involving K4 residue of one subunit and E59, I60, E61, I62 residues of adjacent subunit due to K4E mutation affect the oligomerization efficiency of GroEL subunits while the equilibrium unfolding studies on wild-type monomeric GroEL, native GroEL, and the selected mutants together demonstrate that intra-subunit interactions involving K4 and D523 of the same subunit play a critical role in the thermodynamic stability of both native and monomeric GroEL without affecting the oligomerization of subunits. The stability order between the GroEL_wild-type(M) and its variants is GroEL_wild-type(M) ≥ GroEL_D473C(M)˃GroEL_D523K(M)˃GroEL_K4E.

Atomic-Level Description of Protein Folding inside the GroEL Cavity

Chaperonins (ubiquitous facilitators of protein folding) sequester misfolded proteins within an internal cavity, thus preventing protein aggregation during the process of refolding. GroEL, a tetradecameric bacterial chaperonin, is one of the most studied chaperonins, but the role of the internal cavity in the refolding process is still unclear. It has been suggested that rather than simply isolating proteins while they refold, the GroEL cavity actively promotes protein folding. A detailed characterization of the folding dynamics and thermodynamics of protein substrates encapsulated within the cavity, however, has been difficult to obtain by experimental means, due to the system's complexity and the many steps in the folding cycle. Here, we examine the influence of the GroEL cavity on protein folding based on the results of unbiased, atomistic molecular dynamics simulations. We first verified that the computational setup, which uses a recently developed state-of-the-art force field that more accurately reproduces the aggregation propensity of unfolded states, could recapitulate the essential structural dynamics of GroEL. In these simulations, the GroEL tetradecamer was highly dynamic, transitioning among states corresponding to most of the structures that have been observed experimentally. We then simulated a small, unfolded protein both in the GroEL cavity and in bulk solution and compared the protein's folding process within these two environments. Inside the GroEL cavity, the unfolded protein interacted strongly with the disordered residues in GroEL's C-terminal tails. These interactions stabilized the protein's unfolded states relative to its compact states and increased the roughness of its folding free-energy surface, resulting in slower folding compared to the rate in solution. For larger proteins, which are more typical GroEL substrates, we speculate that these interactions may allow substrates to more quickly escape kinetic traps associated with compact, misfolded states, thereby actively promoting folding.

Gaussian mixture model for coarse-grained modeling from XFEL

We explore the advantage of Gaussian mixture model (GMM) for interpretation of single particle diffraction patterns from X-ray free electron laser (XFEL) experiments. GMM approximates a biomolecular shape by the superposition of Gaussian distributions. As the Fourier transformation of GMM can be quickly performed, we can efficiently simulate XFEL diffraction patterns from approximated structure models. We report that the resolution that GMM can accurately reproduce is proportional to the cubic root of the number of Gaussians used in the modeling. This behavior can be attributed to the correspondence between the number of adjustable parameters in GMM and the amount of sampling points in diffraction space. Furthermore, GMMs can successfully be used to perform angular assignment and to detect conformational variation. These results demonstrate that GMMs serve as useful coarse-grained models for hybrid approach in XFEL single particle experiments.

From deep TLS validation to ensembles of atomic models built from elemental motions. II. Analysis of TLS refinement results by explicit interpretation

TLS modelling was developed by Schomaker and Trueblood to describe atomic displacement parameters through concerted (rigid-body) harmonic motions of an atomic group [Schomaker & Trueblood (1968), Acta Cryst. B24, 63-76]. The results of a TLS refinement are T, L and S matrices that provide individual anisotropic atomic displacement parameters (ADPs) for all atoms belonging to the group. These ADPs can be calculated analytically using a formula that relates the elements of the TLS matrices to atomic parameters. Alternatively, ADPs can be obtained numerically from the parameters of concerted atomic motions corresponding to the TLS matrices. Both procedures are expected to produce the same ADP values and therefore can be used to assess the results of TLS refinement. Here, the implementation of this approach in PHENIX is described and several illustrations, including the use of all models from the PDB that have been subjected to TLS refinement, are provided.

Bringing diffuse X-ray scattering into focus

X-ray crystallography is experiencing a renaissance as a method for probing the protein conformational ensemble. The inherent limitations of Bragg analysis, however, which only reveals the mean structure, have given way to a surge in interest in diffuse scattering, which is caused by structure variations. Diffuse scattering is present in all macromolecular crystallography experiments. Recent studies are shedding light on the origins of diffuse scattering in protein crystallography, and provide clues for leveraging diffuse scattering to model protein motions with atomic detail.

The versatile mutational "repertoire" of Escherichia coli GroEL, a multidomain chaperonin nanomachine

The bacterial chaperonins are highly sophisticated molecular nanomachines, controlled by the hydrolysis of ATP to dynamically trap and remove from the environment unstable protein molecules that are susceptible to denaturation and aggregation. Chaperonins also act to assist in the refolding of these unstable proteins, providing a means by which these proteins may return in active form to the complex environment of the cell. The Escherichia coli GroE chaperonin system is one of the largest protein supramolecular complexes known, whose quaternary structure is required for segregating aggregation-prone proteins. Over the course of more than two decades of research on GroE, it has become accepted that GroE, more specifically the GroEL subunit, is a "high-tolerance" molecular system, capable of accommodating numerous mutations, while retaining its molecular integrity. In some cases, a given site of mutation was revealed to be absolutely required for GroEL function, providing hints regarding the network of signals and triggers that propel this unique system. In other instances, however, a mutation has produced a more delicate response, altering only part of, or in some cases, only a single facet of, the molecular mechanism, and these mutants have often provided invaluable hints on the extent of the complexity underlying chaperonin-assisted protein folding. In this review, we highlight some examples of the latter type of GroEL mutants which compose the unique "mutational repertoire" of GroEL and touch upon the important clues that each mutant provided to the overall effort to elucidate the details of GroE action.

Disassembly/reassembly strategy for the production of highly pure GroEL, a tetradecameric supramolecular machine, suitable for quantitative NMR, EPR and mutational studies

GroEL, a prototypical member of the chaperonin class of chaperones, is a large supramocular machine that assists protein folding and plays an important role in proteostasis. GroEL comprises two heptameric rings, each of which encloses a large cavity that provides a folding chamber for protein substrates. Many questions remain regarding the mechanistic details of GroEL facilitated protein folding. Thus, data at atomic resolution of the type provided by NMR and EPR are invaluable. Such studies often require complete deuteration of GroEL, uniform or residue specific ¹³C and ¹⁵N isotope labeling, and the introduction of selective cysteine mutations for site-specific spin labeling. In addition, high purity GroEL is essential for detailed studies of substrate-GroEL interactions as quantitative interpretation is impossible if the cavities are already occupied and blocked by other protein substrates present in the bacterial expression system. Here we present a new purification protocol designed to provide highly pure GroEL devoid of non-specific protein substrate contamination.

Chaperonin GroEL accelerates protofibril formation and decorates fibrils of the Het-s prion protein

We have studied the interaction of the prototypical chaperonin GroEL with the prion domain of the Het-s protein using solution and solid-state NMR, electron and atomic force microscopies, and EPR. While GroEL accelerates Het-s protofibril formation by several orders of magnitude, the rate of appearance of fibrils is reduced. GroEL remains bound to Het-s throughout the aggregation process and densely decorates the fibrils at a regular spacing of ∼200 Å. GroEL binds to the Het-s fibrils via its apical domain located at the top of the large open ring. Thus, apo GroEL and bullet-shaped GroEL/GroES complexes in which only a single ring is capped by GroES interact with the Het-s fibrils; no evidence is seen for any interaction with football-shaped GroEL/GroES complexes in which both rings are capped by GroES. EPR spectroscopy shows that rotational motion of a nitroxide spin label, placed at the N-terminal end of the first β-strand of Het-s fibrils, is significantly reduced in both Het-s/GroEL aggregates and Het-s fibrils, but virtually completely eliminated in Het-s/GroEL fibrils, suggesting that in the latter, GroEL may come into close proximity to the nitroxide label. Solid-state NMR measurements indicate that GroEL binds to the mobile regions of the Het-s fibril comprising the N-terminal tail and a loop connecting β-strands 4 and 5, consistent with interactions involving GroEL binding consensus sequences located therein.

Subunit conformational variation within individual GroEL oligomers resolved by Cryo-EM

Single-particle electron cryo-microscopy (cryo-EM) is an emerging tool for resolving structures of conformationally heterogeneous particles; however, each structure is derived from an average of many particles with presumed identical conformations. We used a 3.5-Å cryo-EM reconstruction with imposed D7 symmetry to further analyze structural heterogeneity among chemically identical subunits in each GroEL oligomer. Focused classification of the 14 subunits in each oligomer revealed three dominant classes of subunit conformations. Each class resembled a distinct GroEL crystal structure in the Protein Data Bank. The conformational differences stem from the orientations of the apical domain. We mapped each conformation class to its subunit locations within each GroEL oligomer in our dataset. The spatial distributions of each conformation class differed among oligomers, and most oligomers contained 10-12 subunits of the three dominant conformation classes. Adjacent subunits were found to more likely assume the same conformation class, suggesting correlation among subunits in the oligomer. This study demonstrates the utility of cryo-EM in revealing structure dynamics within a single protein oligomer.

X-ray Scattering Studies of Protein Structural Dynamics

X-ray scattering is uniquely suited to the study of disordered systems and thus has the potential to provide insight into dynamic processes where diffraction methods fail. In particular, while X-ray crystallography has been a staple of structural biology for more than half a century and will continue to remain so, a major limitation of this technique has been the lack of dynamic information. Solution X-ray scattering has become an invaluable tool in structural and mechanistic studies of biological macromolecules where large conformational changes are involved. Such systems include allosteric enzymes that play key roles in directing metabolic fluxes of biochemical pathways, as well as large, assembly-line type enzymes that synthesize secondary metabolites with pharmaceutical applications. Furthermore, crystallography has the potential to provide information on protein dynamics via the diffuse scattering patterns that are overlaid with Bragg diffraction. Historically, these patterns have been very difficult to interpret, but recent advances in X-ray detection have led to a renewed interest in diffuse scattering analysis as a way to probe correlated motions. Here, we will review X-ray scattering theory and highlight recent advances in scattering-based investigations of protein solutions and crystals, with a particular focus on complex enzymes.

Suppression of amyloid fibrils using the GroEL apical domain

In E. coli cells, rescue of non-native proteins and promotion of native state structure is assisted by the chaperonin GroEL. An important key to this activity lies in the structure of the apical domain of GroEL (GroEL-AD) (residue 191–376), which recognizes and binds non-native protein molecules through hydrophobic interactions. In this study, we investigated the effects of GroEL-AD on the aggregation of various client proteins (α-Synuclein, Aβ42, and GroES) that lead to the formation of distinct protein fibrils in vitro. We found that GroEL-AD effectively inhibited the fibril formation of these three proteins when added at concentrations above a critical threshold; the specific ratio differed for each client protein, reflecting the relative affinities. The effect of GroEL-AD in all three cases was to decrease the concentration of aggregate-forming unfolded client protein or its early intermediates in solution, thereby preventing aggregation and fibrillation. Binding affinity assays revealed some differences in the binding mechanisms of GroEL-AD toward each client. Our findings suggest a possible applicability of this minimal functioning derivative of the chaperonins (the “minichaperones”) as protein fibrillation modulators and detectors.

Use of evolutionary information in the fitting of atomic level protein models in low resolution cryo-EM map of a protein assembly improves the accuracy of the fitting

Protein-protein interface residues, especially those at the core of the interface, exhibit higher conservation than residues in solvent exposed regions. Here, we explore the ability of this differential conservation to evaluate fittings of atomic models in low-resolution cryo-EM maps and select models from the ensemble of solutions that are often proposed by different model fitting techniques. As a prelude, using a non-redundant and high-resolution structural dataset involving 125 permanent and 95 transient complexes, we confirm that core interface residues are conserved significantly better than nearby non-interface residues and this result is used in the cryo-EM map analysis. From the analysis of inter-component interfaces in a set of fitted models associated with low-resolution cryo-EM maps of ribosomes, chaperones and proteasomes we note that a few poorly conserved residues occur at interfaces. Interestingly a few conserved residues are not in the interface, though they are close to the interface. These observations raise the potential requirement of refitting the models in the cryo-EM maps. We show that sampling an ensemble of models and selection of models with high residue conservation at the interface and in good agreement with the density helps in improving the accuracy of the fit. This study indicates that evolutionary information can serve as an additional input to improve and validate fitting of atomic models in cryo-EM density maps.

Single-Particle Refinement and Variability Analysis in EMAN2.1

CryoEM single-particle reconstruction has been growing rapidly over the last 3 years largely due to the development of direct electron detectors, which have provided data with dramatic improvements in image quality. It is now possible in many cases to produce near-atomic resolution structures, and yet 2/3 of published structures remain at substantially lower resolutions. One important cause for this is compositional and conformational heterogeneity, which is both a resolution-limiting factor and presenting a unique opportunity to better relate structure to function. This manuscript discusses the canonical methods for high-resolution refinement in EMAN2.12, and then considers the wide range of available methods within this package for resolving structural variability, targeting both improved resolution and additional knowledge about particle dynamics.

Chaperonin-enhanced Escherichia coli cell-free expression of functional CXCR4

G protein-coupled receptors (GPCRs) are important therapeutic targets for a broad spectrum of diseases and disorders. Obtaining milligram quantities of functional receptors through the development of robust production methods are highly demanded to probe GPCR structure and functions. In this study, we analyzed synergies of the bacterial chaperonin GroEL-GroES and cell-free expression for the production of functionally folded C-X-C chemokine GPCR type 4 (CXCR4). The yield of soluble CXCR4 in the presence of detergent Brij-35 reached ∼1.1mg/ml. The chaperonin complex added was found to significantly enhance the productive folding of newly synthesized CXCR4, by increasing both the rate (∼30-fold) and the yield (∼1.3-fold) of folding over its spontaneous behavior. Meanwhile, the structural stability of CXCR4 was also improved with supplied GroEL-GroES, as was the soluble expression of biologically active CXCR4 with a ∼1.4-fold increase. The improved stability together with the higher ligand binding affinity suggests more efficient folding. The essential chaperonin GroEL was shown to be partially effective on its own, but for maximum efficiency both GroEL and its co-chaperonin GroES were necessary. The method reported here should prove generally useful for cell-free production of large amounts of natively folded GPCRs, and even other classes of membrane proteins.

Measuring and modeling diffuse scattering in protein X-ray crystallography

X-ray diffraction has the potential to provide rich information about the structural dynamics of macromolecules. To realize this potential, both Bragg scattering, which is currently used to derive macromolecular structures, and diffuse scattering, which reports on correlations in charge density variations, must be measured. Until now, measurement of diffuse scattering from protein crystals has been scarce because of the extra effort of collecting diffuse data. Here, we present 3D measurements of diffuse intensity collected from crystals of the enzymes cyclophilin A and trypsin. The measurements were obtained from the same X-ray diffraction images as the Bragg data, using best practices for standard data collection. To model the underlying dynamics in a practical way that could be used during structure refinement, we tested translation-libration-screw (TLS), liquid-like motions (LLM), and coarse-grained normal-modes (NM) models of protein motions. The LLM model provides a global picture of motions and was refined against the diffuse data, whereas the TLS and NM models provide more detailed and distinct descriptions of atom displacements, and only used information from the Bragg data. Whereas different TLS groupings yielded similar Bragg intensities, they yielded different diffuse intensities, none of which agreed well with the data. In contrast, both the LLM and NM models agreed substantially with the diffuse data. These results demonstrate a realistic path to increase the number of diffuse datasets available to the wider biosciences community and indicate that dynamics-inspired NM structural models can simultaneously agree with both Bragg and diffuse scattering.

Engineering a nanopore with co-chaperonin function

The emergence of an enzymatic function can reveal functional insights and allows the engineering of biological systems with enhanced properties. We engineered an alpha hemolysin nanopore to function as GroES, a protein that, in complex with GroEL, forms a two-stroke protein-folding nanomachine. The transmembrane co-chaperonin was prepared by recombination of GroES functional elements with the nanopore, suggesting that emergent functions in molecular machines can be added bottom-up by incorporating modular elements into preexisting protein scaffolds. The binding of a single-ring version of GroEL to individual GroES nanopores prompted large changes to the unitary nanopore current, most likely reflecting the allosteric transitions of the chaperonin apical domains. One of the GroEL-induced current levels showed fast fluctuations (<1 ms), a characteristic that might be instrumental for efficient substrate encapsulation or folding. In the presence of unfolded proteins, the pattern of current transitions changed, suggesting a possible mechanism in which the free energy of adenosine triphosphate binding and hydrolysis is expended only when substrate proteins are occupied.

Folding of newly translated membrane protein CCR5 is assisted by the chaperonin GroEL-GroES

The in vitro folding of newly translated human CC chemokine receptor type 5 (CCR5), which belongs to the physiologically important family of G protein-coupled receptors (GPCRs), has been studied in a cell-free system supplemented with the surfactant Brij-35. The freshly synthesized CCR5 can spontaneously fold into its biologically active state but only slowly and inefficiently. However, on addition of the GroEL-GroES molecular chaperone system, the folding of the nascent CCR5 was significantly enhanced, as was the structural stability and functional expression of the soluble form of CCR5. The chaperonin GroEL was partially effective on its own, but for maximum efficiency both the GroEL and its GroES lid were necessary. These results are direct evidence for chaperone-assisted membrane protein folding and therefore demonstrate that GroEL-GroES may be implicated in the folding of membrane proteins.

Predicting X-ray diffuse scattering from translation-libration-screw structural ensembles

Identifying the intramolecular motions of proteins and nucleic acids is a major challenge in macromolecular X-ray crystallography. Because Bragg diffraction describes the average positional distribution of crystalline atoms with imperfect precision, the resulting electron density can be compatible with multiple models of motion. Diffuse X-ray scattering can reduce this degeneracy by reporting on correlated atomic displacements. Although recent technological advances are increasing the potential to accurately measure diffuse scattering, computational modeling and validation tools are still needed to quantify the agreement between experimental data and different parameterizations of crystalline disorder. A new tool, phenix.diffuse, addresses this need by employing Guinier's equation to calculate diffuse scattering from Protein Data Bank (PDB)-formatted structural ensembles. As an example case, phenix.diffuse is applied to translation-libration-screw (TLS) refinement, which models rigid-body displacement for segments of the macromolecule. To enable the calculation of diffuse scattering from TLS-refined structures, phenix.tls_as_xyz builds multi-model PDB files that sample the underlying T, L and S tensors. In the glycerophosphodiesterase GpdQ, alternative TLS-group partitioning and different motional correlations between groups yield markedly dissimilar diffuse scattering maps with distinct implications for molecular mechanism and allostery. These methods demonstrate how, in principle, X-ray diffuse scattering could extend macromolecular structural refinement, validation and analysis.

From deep TLS validation to ensembles of atomic models built from elemental motions

The translation-libration-screw model first introduced by Cruickshank, Schomaker and Trueblood describes the concerted motions of atomic groups. Using TLS models can improve the agreement between calculated and experimental diffraction data. Because the T, L and S matrices describe a combination of atomic vibrations and librations, TLS models can also potentially shed light on molecular mechanisms involving correlated motions. However, this use of TLS models in mechanistic studies is hampered by the difficulties in translating the results of refinement into molecular movement or a structural ensemble. To convert the matrices into a constituent molecular movement, the matrix elements must satisfy several conditions. Refining the T, L and S matrix elements as independent parameters without taking these conditions into account may result in matrices that do not represent concerted molecular movements. Here, a mathematical framework and the computational tools to analyze TLS matrices, resulting in either explicit decomposition into descriptions of the underlying motions or a report of broken conditions, are described. The description of valid underlying motions can then be output as a structural ensemble. All methods are implemented as part of the PHENIX project.

Dynamics, flexibility, and allostery in molecular chaperonins

The chaperonins are a family of molecular chaperones present in all three kingdoms of life. They are classified into Group I and Group II. Group I consists of the bacterial variants (GroEL) and the eukaryotic ones from mitochondria and chloroplasts (Hsp60), while Group II consists of the archaeal (thermosomes) and eukaryotic cytosolic variants (CCT or TRiC). Both groups assemble into a dual ring structure, with each ring providing a protective folding chamber for nascent and denatured proteins. Their functional cycle is powered by ATP binding and hydrolysis, which drives a series of structural rearrangements that enable encapsulation and subsequent release of the substrate protein. Chaperonins have elaborate allosteric mechanisms to regulate their functional cycle. Long-range negative cooperativity between the two rings ensures alternation of the folding chambers. Positive intra-ring cooperativity, which facilitates concerted conformational transitions within the protein subunits of one ring, has only been demonstrated for Group I chaperonins. In this review, we describe our present understanding of the underlying mechanisms and the structure-function relationships in these complex protein systems with a particular focus on the structural dynamics, allostery, and associated conformational rearrangements.

Conformational dynamics of a crystalline protein from microsecond-scale molecular dynamics simulations and diffuse X-ray scattering

X-ray diffraction from protein crystals includes both sharply peaked Bragg reflections and diffuse intensity between the peaks. The information in Bragg scattering is limited to what is available in the mean electron density. The diffuse scattering arises from correlations in the electron density variations and therefore contains information about collective motions in proteins. Previous studies using molecular-dynamics (MD) simulations to model diffuse scattering have been hindered by insufficient sampling of the conformational ensemble. To overcome this issue, we have performed a 1.1-μs MD simulation of crystalline staphylococcal nuclease, providing 100-fold more sampling than previous studies. This simulation enables reproducible calculations of the diffuse intensity and predicts functionally important motions, including transitions among at least eight metastable states with different active-site geometries. The total diffuse intensity calculated using the MD model is highly correlated with the experimental data. In particular, there is excellent agreement for the isotropic component of the diffuse intensity, and substantial but weaker agreement for the anisotropic component. Decomposition of the MD model into protein and solvent components indicates that protein-solvent interactions contribute substantially to the overall diffuse intensity. We conclude that diffuse scattering can be used to validate predictions from MD simulations and can provide information to improve MD models of protein motions.

Diffuse X-ray scattering to model protein motions

Problems in biology increasingly need models of protein flexibility to understand and control protein function. At the same time, as they improve, crystallographic methods are marching closer to the limits of what can be learned from Bragg data in isolation. It is thus inevitable that mainstream protein crystallography will turn to diffuse scattering to model protein motions and improve crystallographic models. The time is ripe to make it happen.

B-factor Analysis and Conformational Rearrangement of Aldose Reductase

The NADPH-dependent reduction of glucose reaction that is catalyzed by Aldose Reductase (AR) follows a sequential ordered kinetic mechanism in which the co-factor NADPH binds to the enzyme prior to the aldehyde substrate. The kinetic/structural experiments have found a conformational change involving a hinge-like movement of a surface loop (residues 213-224) which is anticipated to take place upon the binding of the diphosphate moiety of NADPH. The reorientation of this loop, expected to permit the release of NADP⁺, represents the rate-limiting step of the catalytic mechanism. This study reveals: 1) The Translation/Libration/Screw (TLS) analysis of absolute B-factors of apo AR crystal structures indicates that the 212-224 loop might move as a rigid group. 2) Residues that make the flexible loop slide in the AR binary and ternary complexes. 3) The normalized B-factors separate this segment into three different clusters with fewer residues.

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.

Analyses of cobalt-ligand and potassium-ligand bond lengths in metalloproteins: trends and patterns

Cobalt and potassium are biologically important metal elements that are present in a large array of proteins. Cobalt is mostly found in vivo associated with a corrin ring, which represents the core of the vitamin B12 molecule. Potassium is the most abundant metal in the cytosol, and it plays a crucial role in maintaining membrane potential as well as correct protein function. Here, we report a thorough analysis of the geometric properties of cobalt and potassium coordination spheres that was performed with high resolution on a representative set of structures from the Protein Data Bank and complemented by quantum mechanical calculations realized at the DFT level of theory (B3LYP/ SDD) on mononuclear model systems. The results allowed us to draw interesting conclusions on the structural characteristics of both Co and K centers, and to evaluate the importance of effects such as their association energies and intrinsic thermodynamic stabilities. Overall, the results obtained provide useful data for enhancing the atomic models normally applied in theoretical and computational studies of Co or K proteins performed at the quantum mechanical level, and for developing molecular mechanical parameters for treating Co or K coordination spheres in molecular mechanics or molecular dynamics studies.

iMODS: internal coordinates normal mode analysis server

Normal mode analysis (NMA) in internal (dihedral) coordinates naturally reproduces the collective functional motions of biological macromolecules. iMODS facilitates the exploration of such modes and generates feasible transition pathways between two homologous structures, even with large macromolecules. The distinctive internal coordinate formulation improves the efficiency of NMA and extends its applicability while implicitly maintaining stereochemistry. Vibrational analysis, motion animations and morphing trajectories can be easily carried out at different resolution scales almost interactively. The server is versatile; non-specialists can rapidly characterize potential conformational changes, whereas advanced users can customize the model resolution with multiple coarse-grained atomic representations and elastic network potentials. iMODS supports advanced visualization capabilities for illustrating collective motions, including an improved affine-model-based arrow representation of domain dynamics. The generated all-heavy-atoms conformations can be used to introduce flexibility for more advanced modeling or sampling strategies. The server is free and open to all users with no login requirement at http://imods.chaconlab.org.

Probing the dynamic process of encapsulation in Escherichia coli GroEL

Kinetic analyses of GroE-assisted folding provide a dynamic sequence of molecular events that underlie chaperonin function. We used stopped-flow analysis of various fluorescent GroEL mutants to obtain details regarding the sequence of events that transpire immediately after ATP binding to GroEL and GroEL with prebound unfolded proteins. Characterization of GroEL CP86, a circularly permuted GroEL with the polypeptide ends relocated to the vicinity of the ATP binding site, showed that GroES binding and protection of unfolded protein from solution is achieved surprisingly early in the functional cycle, and in spite of greatly reduced apical domain movement. Analysis of fluorescent GroEL SR-1 and GroEL D398A variants suggested that among other factors, the presence of two GroEL rings and a specific conformational rearrangement of Helix M in GroEL contribute significantly to the rapid release of unfolded protein from the GroEL apical domain.

Bayesian analysis of individual electron microscopy images: towards structures of dynamic and heterogeneous biomolecular assemblies

We develop a method to extract structural information from electron microscopy (EM) images of dynamic and heterogeneous molecular assemblies. To overcome the challenge of disorder in the imaged structures, we analyze each image individually, avoiding information loss through clustering or averaging. The Bayesian inference of EM (BioEM) method uses a likelihood-based probabilistic measure to quantify the consistency between each EM image and given structural models. The likelihood function accounts for uncertainties in the molecular position and orientation, variations in the relative intensities and noise in the experimental images. The BioEM formalism is physically intuitive and mathematically simple. We show that for experimental GroEL images, BioEM correctly identifies structures according to the functional state. The top-ranked structure is the corresponding X-ray crystal structure, followed by an EM structure generated previously from a superset of the EM images used here. To analyze EM images of highly flexible molecules, we propose an ensemble refinement procedure, and validate it with synthetic EM maps of the ESCRT-I-II supercomplex. Both the size of the ensemble and its structural members are identified correctly. BioEM offers an alternative to 3D-reconstruction methods, extracting accurate population distributions for highly flexible structures and their assemblies. We discuss limitations of the method, and possible applications beyond ensemble refinement, including the cross-validation and unbiased post-assessment of model structures, and the structural characterization of systems where traditional approaches fail. Overall, our results suggest that the BioEM framework can be used to analyze EM images of both ordered and disordered molecular systems.

TLS from fundamentals to practice

The Translation-Libration-Screw-rotation (TLS) model of rigid-body harmonic displacements introduced in crystallography by Schomaker & Trueblood (1968) is now a routine tool in macromolecular studies and is a feature of most modern crystallographic structure refinement packages. In this review we consider a number of simple examples that illustrate important features of the TLS model. Based on these examples simplified formulae are given for several special cases that may occur in structure modeling and refinement. The derivation of general TLS formulae from basic principles is also provided. This manuscript describes the principles of TLS modeling, as well as some select algorithmic details for practical application. An extensive list of applications references as examples of TLS in macromolecular crystallography refinement is provided.

Crystal structure of a GroEL-ADP complex in the relaxed allosteric state at 2.7 Å resolution

The chaperonin proteins GroEL and GroES are cellular nanomachines driven by the hydrolysis of ATP that facilitate the folding of structurally diverse substrate proteins. In response to ligand binding, the subunits of a ring cycle in a concerted manner through a series of allosteric states (T, R, and R″), enabling work to be performed on the substrate protein. Removing two salt bridges that ordinarily break during the allosteric transitions of the WT permitted the structure of GroEL-ADP in the R state to be solved to 2.7 Å resolution. Whereas the equatorial domain displays almost perfect sevenfold symmetry, the apical domains, to which substrate proteins bind, and to a lesser extent, the intermediate domains display a remarkable asymmetry. Freed of intersubunit contacts, the apical domain of each subunit adopts a different conformation, suggesting a flexibility that permits interaction with diverse substrate proteins. This result contrasts with a previous cryo-EM study of a related allosteric ATP-bound state at lower resolution. After artificially imposing sevenfold symmetry it was concluded that a GroEL ring in the R-ATP state existed in six homogeneous but slightly different states. By imposing sevenfold symmetry on each of the subunits of the crystal structure of GroEL-ADP, we showed that the synthetic rings of (X-ray) GroEL-ADP and (cryo-EM) GroEL-ATP are structurally closely related. A deterministic model, the click stop mechanism, that implied temporal transitions between these states was proposed. Here, however, these conformational states are shown to exist as a structurally heterogeneous ensemble within a single ring.