State-dependent sequential allostery exhibited by chaperonin TRiC/CCT revealed by network analysis of Cryo-EM maps

Scalable computation of anisotropic vibrations for large macromolecular assemblies

The Normal Mode Analysis (NMA) is a standard approach to elucidate the anisotropic vibrations of macromolecules at their folded states, where low-frequency collective motions can reveal rearrangements of domains and changes in the exposed surface of macromolecules. Recent advances in structural biology have enabled the resolution of megascale macromolecules with millions of atoms. However, the calculation of their vibrational modes remains elusive due to the prohibitive cost associated with constructing and diagonalizing the underlying eigenproblem and the current approaches to NMA are not readily adaptable for efficient parallel computing on graphic processing unit (GPU). Here, we present eigenproblem construction and diagonalization approach that implements level-structure bandwidth-reducing algorithms to transform the sparse computation in NMA to a globally-sparse-yet-locally-dense computation, allowing batched tensor products to be most efficiently executed on GPU. We map, optimize, and compare several low-complexity Krylov-subspace eigensolvers, supplemented by techniques such as Chebyshev filtering, sum decomposition, external explicit deflation and shift-and-inverse, to allow fast GPU-resident calculations. The method allows accurate calculation of the first 1000 vibrational modes of some largest structures in PDB ( > 2.4 million atoms) at least 250 times faster than existing methods.

Molecular Dynamics Mappings of the CCT/TRiC Complex-Mediated Protein Folding Cycle Using Diffracted X-ray Tracking

The CCT/TRiC complex is a type II chaperonin that undergoes ATP-driven conformational changes during its functional cycle. Structural studies have provided valuable insights into the mechanism of this process, but real-time dynamics analyses of mammalian type II chaperonins are still scarce. We used diffracted X-ray tracking (DXT) to investigate the intramolecular dynamics of the CCT complex. We focused on three surface-exposed loop regions of the CCT1 subunit: the loop regions of the equatorial domain (E domain), the E and intermediate domain (I domain) juncture near the ATP-binding region, and the apical domain (A domain). Our results showed that the CCT1 subunit predominantly displayed rotational motion, with larger mean square displacement (MSD) values for twist (χ) angles compared with tilt (θ) angles. Nucleotide binding had a significant impact on the dynamics. In the absence of nucleotides, the region between the E and I domain juncture could act as a pivotal axis, allowing for greater motion of the E domain and A domain. In the presence of nucleotides, the nucleotides could wedge into the ATP-binding region, weakening the role of the region between the E and I domain juncture as the rotational axis and causing the CCT complex to adopt a more compact structure. This led to less expanded MSD curves for the E domain and A domain compared with nucleotide-absent conditions. This change may help to stabilize the functional conformation during substrate binding. This study is the first to use DXT to probe the real-time molecular dynamics of mammalian type II chaperonins at the millisecond level. Our findings provide new insights into the complex dynamics of chaperonins and their role in the functional folding cycle.

Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing

Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.

Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0

New computational biophysics pipelines for analysing the global dynamics of structural ensembles and large, dynamic complexes resolved by cryoEM are reviewed.

Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.

Chaperonin Mechanisms: Multiple and (Mis)Understood?

The chaperonins are ubiquitous and essential nanomachines that assist in protein folding in an ATP-driven manner. They consist of two back-to-back stacked oligomeric rings with cavities in which protein (un)folding can take place in a shielding environment. This review focuses on GroEL from Escherichia coli and the eukaryotic chaperonin-containing t-complex polypeptide 1, which differ considerably in their reaction mechanisms despite sharing a similar overall architecture. Although chaperonins feature in many current biochemistry textbooks after being studied intensively for more than three decades, key aspects of their reaction mechanisms remain under debate and are discussed in this review. In particular, it is unclear whether a universal reaction mechanism operates for all substrates and whether it is passive, i.e., aggregation is prevented but the folding pathway is unaltered, or active. It is also unclear how chaperonin clients are distinguished from nonclients and what are the precise roles of the cofactors with which chaperonins interact.

Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Approximating deformation fields for the analysis of continuous heterogeneity of biological macromolecules by 3D Zernike polynomials

Structural biology has evolved greatly due to the advances introduced in fields like electron microscopy. This image-capturing technique, combined with improved algorithms and current data processing software, allows the recovery of different conformational states of a macromolecule, opening new possibilities for the study of its flexibility and dynamic events. However, the ensemble analysis of these different conformations, and in particular their placement into a common variable space in which the differences and similarities can be easily recognized, is not an easy matter. To simplify the analysis of continuous heterogeneity data, this work proposes a new automatic algorithm that relies on a mathematical basis defined over the sphere to estimate the deformation fields describing conformational transitions among different structures. Thanks to the approximation of these deformation fields, it is possible to describe the forces acting on the molecules due to the presence of different motions. It is also possible to represent and compare several structures in a low-dimensional mapping, which summarizes the structural characteristics of different states. All these analyses are integrated into a common framework, providing the user with the ability to combine them seamlessly. In addition, this new approach is a significant step forward compared with principal component analysis and normal mode analysis of cryo-electron microscopy maps, avoiding the need to select components or modes and producing localized analysis.