CryoEM-based hybrid modeling approaches for structure determination

Recent structural advances in bacterial chemotaxis signalling

Bacterial chemosensory arrays have served as a model system for in-situ structure determination, clearly cataloguing the improvement of cryo-electron tomography (cryoET) over the past decade. In recent years, this has culminated in an accurately fitted atomistic model for the full-length core signalling unit (CSU) and numerous insights into the function of the transmembrane receptors responsible for signal transduction. Here, we review the achievements of the latest structural advances in bacterial chemosensory arrays and the developments which have made such advances possible.

Fostering discoveries in the era of exascale computing: How the next generation of supercomputers empowers computational and experimental biophysics alike

Over a century ago, physicists started broadly relying on theoretical models to guide new experiments. Soon thereafter, chemists began doing the same. Now, biological research enters a new era when experiment and theory walk hand in hand. Novel software and specialized hardware became essential to understand experimental data and propose new models. In fact, current petascale computing resources already allow researchers to reach unprecedented levels of simulation throughput to connect in silico and in vitro experiments. The reduction in cost and improved access allowed a large number of research groups to adopt supercomputing resources and techniques. Here, we outline how large-scale computing has evolved to expand decades-old research, spark new research efforts, and continuously connect simulation and observation. For instance, multiple publicly and privately funded groups have dedicated extensive resources to develop artificial intelligence tools for computational biophysics, from accelerating quantum chemistry calculations to proposing protein structure models. Moreover, advances in computer hardware have accelerated data processing from single-molecule experimental observations and simulations of chemical reactions occurring throughout entire cells. The combination of software and hardware has opened the way for exascale computing and the production of the first public exascale supercomputer, Frontier, inaugurated by the Oak Ridge National Laboratory in 2022. Ultimately, the popularization and development of computational techniques and the training of researchers to use them will only accelerate the diversification of tools and learning resources for future generations.

CLoNe: automated clustering based on local density neighborhoods for application to biomolecular structural ensembles

Motivation

Proteins are intrinsically dynamic entities. Flexibility sampling methods, such as molecular dynamics or those arising from integrative modeling strategies, are now commonplace and enable the study of molecular conformational landscapes in many contexts. Resulting structural ensembles increase in size as technological and algorithmic advancements take place, making their analysis increasingly demanding. In this regard, cluster analysis remains a go-to approach for their classification. However, many state-of-the-art algorithms are restricted to specific cluster properties. Combined with tedious parameter fine-tuning, cluster analysis of protein structural ensembles suffers from the lack of a generally applicable and easy to use clustering scheme.

Results

We present CLoNe, an original Python-based clustering scheme that builds on the Density Peaks algorithm of Rodriguez and Laio. CLoNe relies on a probabilistic analysis of local density distributions derived from nearest neighbors to find relevant clusters regardless of cluster shape, size, distribution and amount. We show its capabilities on many toy datasets with properties otherwise dividing state-of-the-art approaches and improves on the original algorithm in key aspects. Applied to structural ensembles, CLoNe was able to extract meaningful conformations from membrane binding events and ligand-binding pocket opening as well as identify dominant dimerization motifs or inter-domain organization. CLoNe additionally saves clusters as individual trajectories for further analysis and provides scripts for automated use with molecular visualization software.

Availability and implementation

www.epfl.ch/labs/lbm/resources, github.com/LBM-EPFL/CLoNe.

Supplementary information

Supplementary data are available at Bioinformatics online.

Integration of Cryo-EM Model Building Software in Scipion

Advances in cryo-electron microscopy (cryo-EM) have made it possible to obtain structures of large biological macromolecules at near-atomic resolution. This "resolution revolution" has encouraged the use and development of modeling tools able to produce high-quality atomic models from cryo-EM density maps. Unfortunately, many practical problems appear when combining different packages in the same processing workflow, which make difficult the use of these tools by non-experts and, therefore, reduce their utility. We present here a major extension of the image processing framework Scipion that provides inter-package integration in the model building area and full tracking of the complete workflow, from image processing to structure validation.

Complete structure of the chemosensory array core signalling unit in an E. coli minicell strain

Motile bacteria sense chemical gradients with transmembrane receptors organised in supramolecular signalling arrays. Understanding stimulus detection and transmission at the molecular level requires precise structural characterisation of the array building block known as a core signalling unit. Here we introduce an Escherichia coli strain that forms small minicells possessing extended and highly ordered chemosensory arrays. We use cryo-electron tomography and subtomogram averaging to provide a three-dimensional map of a complete core signalling unit, with visible densities corresponding to the HAMP and periplasmic domains. This map, combined with previously determined high resolution structures and molecular dynamics simulations, yields a molecular model of the transmembrane core signalling unit and enables spatial localisation of its individual domains. Our work thus offers a solid structural basis for the interpretation of a wide range of existing data and the design of further experiments to elucidate signalling mechanisms within the core signalling unit and larger array.

Structure and dynamics of the E. coli chemotaxis core signaling complex by cryo-electron tomography and molecular simulations

To enable the processing of chemical gradients, chemotactic bacteria possess large arrays of transmembrane chemoreceptors, the histidine kinase CheA, and the adaptor protein CheW, organized as coupled core-signaling units (CSU). Despite decades of study, important questions surrounding the molecular mechanisms of sensory signal transduction remain unresolved, owing especially to the lack of a high-resolution CSU structure. Here, we use cryo-electron tomography and sub-tomogram averaging to determine a structure of the Escherichia coli CSU at sub-nanometer resolution. Based on our experimental data, we use molecular simulations to construct an atomistic model of the CSU, enabling a detailed characterization of CheA conformational dynamics in its native structural context. We identify multiple, distinct conformations of the critical P4 domain as well as asymmetries in the localization of the P3 bundle, offering several novel insights into the CheA signaling mechanism.

Diversity of Bacterial Chemosensory Arrays

Chemotaxis is crucial for the survival of bacteria, and the signaling systems associated with it exhibit a high level of evolutionary conservation. The architecture of the chemosensory array and the signal transduction mechanisms have been extensively studied in Escherichia coli. More recent studies have revealed a vast diversity of the chemosensory system among bacteria. Unlike E. coli, some bacteria assemble more than one chemosensory array and respond to a broader spectrum of environmental and internal stimuli. These chemosensory arrays exhibit a great variability in terms of protein composition, cellular localization, and functional variability. Here, we present recent findings that emphasize the extent of diversity in chemosensory arrays and highlight the importance of studying chemosensory arrays in bacteria other than the common model organisms.

Advances in cryo-electron tomography and subtomogram averaging and classification

•

Cryo-electron tomography (cryoET) subtomogram averaging has emerged as a structural biology method for sparse and heterogenerous sampls.

•

CryoET subtomogram averaging enables in situ structure determination.

•

CryoET subtomogram classification can delineate different conformational states of macromolecular complexes.

•

Future developments in cryoET and correlative super resolution microscopy promises to bring unprecedented integration of cell biology and structural biology.

Cryo-electron tomography (cryoET) can provide 3D reconstructions, or tomograms, of pleomorphic objects such as organelles or cells in their close-to-native states. Subtomograms that contain repetitive structures can be further extracted and subjected to averaging and classification to improve resolution, and this process has become an emerging structural biology method referred to as cryoET subtomogram averaging and classification (cryoSTAC). Recent technical advances in cryoSTAC have had a profound impact on many fields in biology. Here, I review recent exciting work on several macromolecular assemblies demonstrating the power of cryoSTAC for in situ structure analysis and discuss challenges and future directions.

In Situ Conformational Changes of the Escherichia coli Serine Chemoreceptor in Different Signaling States

Tsr, the serine chemoreceptor in Escherichia coli, transduces signals from a periplasmic ligand-binding site to its cytoplasmic tip, where it controls the activity of the CheA kinase. To function, Tsr forms trimers of homodimers (TODs), which associate in vivo with the CheA kinase and CheW coupling protein. Together, these proteins assemble into extended hexagonal arrays. Here, we use cryo-electron tomography and molecular dynamics simulation to study Tsr in the context of a near-native array, characterizing its signaling-related conformational changes at both the individual dimer and the trimer level. In particular, we show that individual Tsr dimers within a trimer exhibit asymmetric flexibilities that are a function of the signaling state, highlighting the effect of their different protein interactions at the receptor tips. We further reveal that the dimer compactness of the Tsr trimer changes between signaling states, transitioning at the glycine hinge from a compact conformation in the kinase-OFF state to an expanded conformation in the kinase-ON state. Hence, our results support a crucial role for the glycine hinge: to allow the receptor flexibility necessary to achieve different signaling states while also maintaining structural constraints imposed by the membrane and extended array architecture.IMPORTANCE In Escherichia coli, membrane-bound chemoreceptors, the histidine kinase CheA, and coupling protein CheW form highly ordered chemosensory arrays. In core signaling complexes, chemoreceptor trimers of dimers undergo conformational changes, induced by ligand binding and sensory adaptation, which regulate kinase activation. Here, we characterize by cryo-electron tomography the kinase-ON and kinase-OFF conformations of the E. coli serine receptor in its native array context. We found distinctive structural differences between the members of a receptor trimer, which contact different partners in the signaling unit, and structural differences between the ON and OFF signaling complexes. Our results provide new insights into the signaling mechanism of chemoreceptor arrays and suggest an important functional role for a previously postulated flexible region and glycine hinge in the receptor molecule.

FragFit: a web-application for interactive modeling of protein segments into cryo-EM density maps

Cryo-electron microscopy (cryo-EM) is a standard method to determine the three-dimensional structures of molecular complexes. However, easy to use tools for modeling of protein segments into cryo-EM maps are sparse. Here, we present the FragFit web-application, a web server for interactive modeling of segments of up to 35 amino acids length into cryo-EM density maps. The fragments are provided by a regularly updated database containing at the moment about 1 billion entries extracted from PDB structures and can be readily integrated into a protein structure. Fragments are selected based on geometric criteria, sequence similarity and fit into a given cryo-EM density map. Web-based molecular visualization with the NGL Viewer allows interactive selection of fragments. The FragFit web-application, accessible at http://proteinformatics.de/FragFit, is free and open to all users, without any login requirements.

Monomeric Intermediates Formed by Vesiculovirus Glycoprotein during Its Low-pH-induced Structural Transition

•

Vesiculovirus G is the prototype of class III viral fusion glycoproteins.

•

The structures of both G pre- and post-fusion conformation have been determined.

•

The structure of monomeric intermediates reveals the pathway of the transition.

•

A fusion-loop-exposing antiparallel dimer may initiate the fusion process.

•

Those data challenge the current model proposed for viral membrane fusion.

Rigid-Body Fitting of Atomic Models on 3D Density Maps of Electron Microscopy

Cryo electron microscopy has revolutionarily evolved for the determination of the 3D structure of macromolecular complexes. The modeling procedures on the 3D density maps of electron microscopy are roughly classified into three categories: fitting, de novo modeling and refinement. The registered atomic models from the maps have mostly been hand-built and auto-refined. Several programs aiming at automatic modeling have also been developed using various kinds of molecular representations. Among these three classes of the modeling procedures, the rigid body fitting is reviewed here, because it is the most basic modeling process applied before the other steps. The fitting problems are classified as the fittings of single subunit or multiple subunits, and the fittings on global or local parts of maps. A higher resolution map enables more local fitting. Various molecular representations have been employed in the fitting programs. A point and digital image models are generally used to represent molecules, but new representations, such as the Gaussian mixture model, have been applied recently.