LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

Structure of the native chemotaxis core signaling unit from phage E-protein lysed E. coli cells

IMPORTANCE

Bacterial chemotaxis is a ubiquitous behavior that enables cell movement toward or away from specific chemicals. It serves as an important model for understanding cell sensory signal transduction and motility. Characterization of the molecular mechanisms underlying chemotaxis is of fundamental interest and requires a high-resolution structural picture of the sensing machinery, the chemosensory array. In this study, we combine cryo-electron tomography and molecular simulation to present the complete structure of the core signaling unit, the basic building block of chemosensory arrays, from Escherichia coli. Our results provide new insight into previously poorly-resolved regions of the complex and offer a structural basis for designing new experiments to test mechanistic hypotheses.

Structural basis for SMAC-mediated antagonism of caspase inhibition by the giant ubiquitin ligase BIRC6

Certain inhibitor of apoptosis (IAP) family members are sentinel proteins that prevent untimely cell death by inhibiting caspases. Antagonists, including second mitochondria-derived activator of caspases (SMAC), regulate IAPs and drive cell death. Baculoviral IAP repeat-containing protein 6 (BIRC6), a giant IAP with dual E2 and E3 ubiquitin ligase activity, regulates programmed cell death through unknown mechanisms. We show that BIRC6 directly restricts executioner caspase-3 and -7 and ubiquitinates caspase-3, -7, and -9, working exclusively with noncanonical E1, UBA6. Notably, we show that SMAC suppresses both mechanisms. Cryo-electron microscopy structures of BIRC6 alone and in complex with SMAC reveal that BIRC6 is an antiparallel dimer juxtaposing the substrate-binding module against the catalytic domain. Furthermore, we discover that SMAC multisite binding to BIRC6 results in a subnanomolar affinity interaction, enabling SMAC to competitively displace caspases, thus antagonizing BIRC6 anticaspase function.

Mechanisms of E. coli chemotaxis signaling pathways visualized using cryoET and computational approaches

Chemotaxis signaling pathways enable bacteria to sense and respond to their chemical environment and, in some species, are critical for lifestyle processes such as biofilm formation and pathogenesis. The signal transduction underlying chemotaxis behavior is mediated by large, highly ordered protein complexes known as chemosensory arrays. For nearly two decades, cryo-electron tomography (cryoET) has been used to image chemosensory arrays, providing an increasingly detailed understanding of their structure and function. In this mini-review, we provide an overview of the use of cryoET to study chemosensory arrays, including imaging strategies, key results, and outstanding questions. We further discuss the application of molecular modeling and simulation techniques to complement structure determination efforts and provide insight into signaling mechanisms. We close the review with a brief outlook, highlighting promising future directions for the field.

High-resolution reconstruction of a Jumbo-bacteriophage infecting capsulated bacteria using hyperbranched tail fibers

The Klebsiella jumbo myophage ϕKp24 displays an unusually complex arrangement of tail fibers interacting with a host cell. In this study, we combine cryo-electron microscopy methods, protein structure prediction methods, molecular simulations, microbiological and machine learning approaches to explore the capsid, tail, and tail fibers of ϕKp24. We determine the structure of the capsid and tail at 4.1 Å and 3.0 Å resolution. We observe the tail fibers are branched and rearranged dramatically upon cell surface attachment. This complex configuration involves fourteen putative tail fibers with depolymerase activity that provide ϕKp24 with the ability to infect a broad panel of capsular polysaccharide (CPS) types of Klebsiella pneumoniae. Our study provides structural and functional insight into how ϕKp24 adapts to the variable surfaces of capsulated bacterial pathogens, which is useful for the development of phage therapy approaches against pan-drug resistant K. pneumoniae strains.

Structure and activity of particulate methane monooxygenase arrays in methanotrophs

Methane-oxidizing bacteria play a central role in greenhouse gas mitigation and have potential applications in biomanufacturing. Their primary metabolic enzyme, particulate methane monooxygenase (pMMO), is housed in copper-induced intracytoplasmic membranes (ICMs), of which the function and biogenesis are not known. We show by serial cryo-focused ion beam (cryoFIB) milling/scanning electron microscope (SEM) volume imaging and lamellae-based cellular cryo-electron tomography (cryoET) that these ICMs are derived from the inner cell membrane. The pMMO trimer, resolved by cryoET and subtomogram averaging to 4.8 Å in the ICM, forms higher-order hexagonal arrays in intact cells. Array formation correlates with increased enzymatic activity, highlighting the importance of studying the enzyme in its native environment. These findings also demonstrate the power of cryoET to structurally characterize native membrane enzymes in the cellular context.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

Hexameric rings of the scaffolding protein CheW enhance response sensitivity and cooperativity in Escherichia coli chemoreceptor arrays

The Escherichia coli chemoreceptor array is a supramolecular assembly that enables cells to respond to extracellular cues dynamically and with great precision and sensitivity. In the array, transmembrane receptors organized as trimers of dimers are connected at their cytoplasmic tips by hexameric rings of alternating subunits of the kinase CheA and the scaffolding protein CheW (CheA-CheW rings). Interactions of CheW molecules with the members of receptor trimers not directly bound to CheA-CheW rings may lead to the formation of hexameric CheW rings in the chemoreceptor array. Here, we detected such CheW rings with a cellular cysteine-directed cross-linking assay and explored the requirements for their formation and their participation in array assembly. We found that CheW ring formation varied with cellular CheW abundance, depended on the presence of receptors capable of a trimer-of-dimers arrangement, and did not require CheA. Cross-linking studies of a CheA~CheW fusion protein incapable of forming homomeric CheW oligomers demonstrated that CheW rings were not essential for the assembly of CheA-containing arrays. Förster resonance energy transfer (FRET)-based kinase assays of arrays containing variable amounts of CheW rings revealed that CheW rings enhanced the cooperativity and the sensitivity of the responses to attractants. We propose that six-membered CheW rings provide the additional interconnectivity required for optimal signaling and gradient tracking performance by chemosensory arrays.

UVC inactivation of pathogenic samples suitable for cryo-EM analysis

Cryo-electron microscopy has become an essential tool to understand structure and function of biological samples. Especially for pathogens, such as disease-causing bacteria and viruses, insights gained by cryo-EM can aid in developing cures. However, due to the biosafety restrictions of pathogens, samples are often treated by chemical fixation to render the pathogen inert, affecting the ultrastructure of the sample. Alternatively, researchers use in vitro or ex vivo models, which are non-pathogenic but lack the complexity of the pathogen of interest. Here we show that ultraviolet-C (UVC) radiation applied at cryogenic temperatures can be used to eliminate or dramatically reduce the infectivity of Vibrio cholerae and the bacterial virus, the ICP1 bacteriophage. We show no discernable structural impact of this treatment of either sample using two cryo-EM methods: cryo-electron tomography followed by sub-tomogram averaging, and single particle analysis (SPA). Additionally, we applied the UVC irradiation to the protein apoferritin (ApoF), which is a widely used test sample for high-resolution SPA studies. The UVC-treated ApoF sample resulted in a 2.1 Å structure indistinguishable from an untreated published map. This research demonstrates that UVC treatment is an effective and inexpensive addition to the cryo-EM sample preparation toolbox.

Depelteau et al. present a new method to inactivate cryo-EM samples from pathogenic organisms before imaging using ultraviolet-C radiation in cryogenic conditions. This method allows for the inexpensive preparation of cryo-EM samples with no discernable structural impact of the treatment.

Alternative Architecture of the E. coli Chemosensory Array

Chemotactic responses in motile bacteria are the result of sophisticated signal transduction by large, highly organized arrays of sensory proteins. Despite tremendous progress in the understanding of chemosensory array structure and function, a structural basis for the heightened sensitivity of networked chemoreceptors is not yet complete. Here, we present cryo-electron tomography visualisations of native-state chemosensory arrays in E. coli minicells. Strikingly, these arrays appear to exhibit a p2-symmetric array architecture that differs markedly from the p6-symmetric architecture previously described in E. coli. Based on this data, we propose molecular models of this alternative architecture and the canonical p6-symmetric assembly. We evaluate our observations and each model in the context of previously published data, assessing the functional implications of an alternative architecture and effects for future studies.

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.

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.

CryoEM-based hybrid modeling approaches for structure determination

Recent advances in cryo-electron microscopy (cryoEM) have dramatically improved the resolutions at which vitrified biological specimens can be studied, revealing new structural and mechanistic insights over a broad range of spatial scales. Bolstered by these advances, much effort has been directed toward the development of hybrid modeling methodologies for the construction and refinement of high-fidelity atomistic models from cryoEM data. In this brief review, we will survey the key elements of cryoEM-based hybrid modeling, providing an overview of available computational tools and strategies as well as several recent applications.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

CryoEM and computer simulations reveal a novel kinase conformational switch in bacterial chemotaxis signaling

Chemotactic responses in bacteria require large, highly ordered arrays of sensory proteins to mediate the signal transduction that ultimately controls cell motility. A mechanistic understanding of the molecular events underlying signaling, however, has been hampered by the lack of a high-resolution structural description of the extended array. Here, we report a novel reconstitution of the array, involving the receptor signaling domain, histidine kinase CheA, and adaptor protein CheW, as well as a density map of the core-signaling unit at 11.3 Å resolution, obtained by cryo-electron tomography and sub-tomogram averaging. Extracting key structural constraints from our density map, we computationally construct and refine an atomic model of the core array structure, exposing novel interfaces between the component proteins. Using all-atom molecular dynamics simulations, we further reveal a distinctive conformational change in CheA. Mutagenesis and chemical cross-linking experiments confirm the importance of the conformational dynamics of CheA for chemotactic function.

CheY's acetylation sites responsible for generating clockwise flagellar rotation in Escherichia coli

Stimulation of Escherichia coli with acetate elevates the acetylation level of the chemotaxis response regulator CheY. This elevation, in an unknown mechanism, activates CheY to generate clockwise rotation. Here, using quantitative selective reaction monitoring mass spectrometry and high-resolution targeted mass spectrometry, we identified K91 and K109 as the major sites whose acetylation level in vivo increases in response to acetate. Employing single and multiple lysine replacements in CheY, we found that K91 and K109 are also the sites mainly responsible for acetate-dependent clockwise generation. Furthermore, we showed that clockwise rotation is repressed when residue K91 is nonmodified, as evidenced by an increased ability of CheY to generate clockwise rotation when K91 was acetylated or replaced by specific amino acids. Using molecular dynamics simulations, we show that K91 repression is manifested in the conformational dynamics of the β4α4 loop, shifted toward an active state upon mutation. Removal of β4α4 loop repression may represent a general activation mechanism in CheY, pertaining also to the canonical phosphorylation activation pathway as suggested by crystal structures of active and inactive CheY from Thermotoga maritima. By way of elimination, we further suggest that K109 acetylation is actively involved in generating clockwise rotation.