iMODFIT: efficient and robust flexible fitting based on vibrational analysis in internal coordinates

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.

Flotillin-mediated stabilization of unfolded proteins in bacterial membrane microdomains

The function of many bacterial processes depends on the formation of functional membrane microdomains (FMMs), which resemble the lipid rafts of eukaryotic cells. However, the mechanism and the biological function of these membrane microdomains remain unclear. Here, we show that FMMs in the pathogen methicillin-resistant Staphylococcus aureus (MRSA) are dedicated to confining and stabilizing proteins unfolded due to cellular stress. The FMM scaffold protein flotillin forms a clamp-shaped oligomer that holds unfolded proteins, stabilizing them and favoring their correct folding. This process does not impose a direct energy cost on the cell and is crucial to survival of ATP-depleted bacteria, and thus to pathogenesis. Consequently, FMM disassembling causes the accumulation of unfolded proteins, which compromise MRSA viability during infection and cause penicillin re-sensitization due to PBP2a unfolding. Thus, our results indicate that FMMs mediate ATP-independent stabilization of unfolded proteins, which is essential for bacterial viability during infection.

Molecular basis for transposase activation by a dedicated AAA+ ATPase

Transposases drive chromosomal rearrangements and the dissemination of drug-resistance genes and toxins1-3. Although some transposases act alone, many rely on dedicated AAA+ ATPase subunits that regulate site selectivity and catalytic function through poorly understood mechanisms. Using IS21 as a model transposase system, we show how an ATPase regulator uses nucleotide-controlled assembly and DNA deformation to enable structure-based site selectivity, transposase recruitment, and activation and integration. Solution and cryogenic electron microscopy studies show that the IstB ATPase self-assembles into an autoinhibited pentamer of dimers that tightly curves target DNA into a half-coil. Two of these decamers dimerize, which stabilizes the target nucleic acid into a kinked S-shaped configuration that engages the IstA transposase at the interface between the two IstB oligomers to form an approximately 1 MDa transpososome complex. Specific interactions stimulate regulator ATPase activity and trigger a large conformational change on the transposase that positions the catalytic site to perform DNA strand transfer. These studies help explain how AAA+ ATPase regulators-which are used by classical transposition systems such as Tn7, Mu and CRISPR-associated elements-can remodel their substrate DNA and cognate transposases to promote function.

Structural dynamics of RAF1-HSP90-CDC37 and HSP90 complexes reveal asymmetric client interactions and key structural elements

RAF kinases are integral to the RAS-MAPK signaling pathway, and proper RAF1 folding relies on its interaction with the chaperone HSP90 and the cochaperone CDC37. Understanding the intricate molecular interactions governing RAF1 folding is crucial for comprehending this process. Here, we present a cryo-EM structure of the closed-state RAF1-HSP90-CDC37 complex, where the C-lobe of the RAF1 kinase domain binds to one side of the HSP90 dimer, and an unfolded N-lobe segment of the RAF1 kinase domain threads through the center of the HSP90 dimer. CDC37 binds to the kinase C-lobe, mimicking the N-lobe with its HxNI motif. We also describe structures of HSP90 dimers without RAF1 and CDC37, displaying only N-terminal and middle domains, which we term the semi-open state. Employing 1 μs atomistic simulations, energetic decomposition, and comparative structural analysis, we elucidate the dynamics and interactions within these complexes. Our quantitative analysis reveals that CDC37 bridges the HSP90-RAF1 interaction, RAF1 binds HSP90 asymmetrically, and that HSP90 structural elements engage RAF1's unfolded region. Additionally, N- and C-terminal interactions stabilize HSP90 dimers, and molecular interactions in HSP90 dimers rearrange between the closed and semi-open states. Our findings provide valuable insight into the contributions of HSP90 and CDC37 in mediating client folding.

Cryo-EM structure and B-factor refinement with ensemble representation

Cryo-EM experiments produce images of macromolecular assemblies that are combined to produce three-dimensional density maps. Typically, atomic models of the constituent molecules are fitted into these maps, followed by a density-guided refinement. We introduce TEMPy-ReFF, a method for atomic structure refinement in cryo-EM density maps. Our method represents atomic positions as components of a Gaussian mixture model, utilising their variances as B-factors, which are used to derive an ensemble description. Extensively tested on a substantial dataset of 229 cryo-EM maps from EMDB ranging in resolution from 2.1-4.9 Å with corresponding PDB and CERES atomic models, our results demonstrate that TEMPy-ReFF ensembles provide a superior representation of cryo-EM maps. On a single-model basis, it performs similarly to the CERES re-refinement protocol, although there are cases where it provides a better fit to the map. Furthermore, our method enables the creation of composite maps free of boundary artefacts. TEMPy-ReFF is useful for better interpretation of flexible structures, such as those involving RNA, DNA or ligands.

Structural Analysis of Protein Complexes by Cryo-Electron Microscopy

Structural studies of bio-complexes using single particle cryo-Electron Microscopy (cryo-EM) is nowadays a well-established technique in structural biology and has become competitive with X-ray crystallography. Development of digital registration systems for electron microscopy images and algorithms for the fast and efficient processing of the recorded images and their following analysis has facilitated the determination of structures at near-atomic resolution. The latest advances in EM have enabled the determination of protein complex structures at 1.4-3 Å resolution for an extremely broad range of sizes (from ~100 kDa up to hundreds of MDa (Bartesaghi et al., Science 348(6239):1147-1151, 2015; Herzik et al., Nat Commun 10:1032, 2019; Wu et al., J Struct Biol X 4:100020, 2020; Zhang et al., Nat Commun 10:5511, 2019; Zhang et al., Cell Res 30(12):1136-1139, 2020; Yip et al., Nature 587(7832):157-161, 2020; https://www.ebi.ac.uk/emdb/statistics/emdb_resolution_year )). In 2022, nearly 1200 structures deposited to the EMDB database were at a resolution of better than 3 Å ( https://www.ebi.ac.uk/emdb/statistics/emdb_resolution_year ).To date, the highest resolutions have been achieved for apoferritin, which comprises a homo-oligomer of high point group symmetry (O432) and has rigid organization together with high stability (Zhang et al., Cell Res 30(12):1136-1139, 2020; Yip et al., Nature 587(7832):157-161, 2020). It has been used as a test object for the assessments of modern cryo-microscopes and processing methods during the last 5 years. In contrast to apoferritin bacterial secretion systems are typical examples of multi protein complexes exhibiting high flexibility owing to their functions relating to the transportation of small molecules, proteins, and DNA into the extracellular space or target cells. This makes their structural characterization extremely challenging (Barlow, Methods Mol Biol 532:397-411, 2009; Costa et al., Nat Rev Microbiol 13:343-359, 2015). The most feasible approach to reveal their spatial organization and functional modification is cryo-electron microscopy (EM). During the last decade, structural cryo-EM has become broadly used for the analysis of the bio-complexes that comprise multiple components and are not amenable to crystallization (Lyumkis, J Biol Chem 294:5181-5197, 2019; Orlova and Saibil, Methods Enzymol 482:321-341, 2010; Orlova and Saibil, Chem Rev 111(12):7710-7748, 2011).In this review, we will describe the basics of sample preparation for cryo-EM, the principles of digital data collection, and the logistics of image analysis focusing on the common steps required for reconstructions of both small and large biological complexes together with refinement of their structures to nearly atomic resolution. The workflow of processing will be illustrated by examples of EM analysis of Type IV Secretion System.

Observation of Bacteriophage Ultrastructure by Cryo-Electron Microscopy

Transmission electron microscopy (TEM) is an ideal method to observe and determine the structure of bacteriophages. From early studies by negative staining to the present atomic structure models derived from cryo-TEM, bacteriophage detection, classification, and structure determination have been mostly done by electron microscopy. Although embedding in metal salts has been a routine method for virus observation for many years, the preservation of bacteriophages in a thin layer of fast frozen buffer has proven to be the most convenient preparation method for obtaining images using cryo-electron microscopy (cryo-EM). In this technique, frozen samples are observed at liquid nitrogen temperature, and the images are acquired using different recording media. The incorporation of direct electron detectors has been a fundamental step in achieving atomic resolution images of a number of viruses. These projection images can be numerically combined using different approaches to render a three-dimensional model of the virus. For those viral components exhibiting any symmetry, averaging can nowadays achieve atomic structures in most cases. Image processing methods have also evolved to improve the resolution in asymmetric viral components or regions showing different types of symmetries (symmetry mismatch).

Structure and activation mechanism of the Makes caterpillars floppy 1 toxin

The bacterial Makes caterpillars floppy 1 (Mcf1) toxin promotes apoptosis in insects, leading to loss of body turgor and death. The molecular mechanism underlying Mcf1 intoxication is poorly understood. Here, we present the cryo-EM structure of Mcf1 from Photorhabdus luminescens, revealing a seahorse-like shape with a head and tail. While the three head domains contain two effectors, as well as an activator-binding domain (ABD) and an autoprotease, the tail consists of two putative translocation and three putative receptor-binding domains. Rearrangement of the tail moves the C-terminus away from the ABD and allows binding of the host cell ADP-ribosylation factor 3, inducing conformational changes that position the cleavage site closer to the protease. This distinct activation mechanism that is based on a hook-loop interaction results in three autocleavage reactions and the release of two toxic effectors. Unexpectedly, the BH3-like domain containing ABD is not an active effector. Our findings allow us to understand key steps of Mcf1 intoxication at the molecular level.

Application of conformational space annealing to the protein structure modeling using cryo-EM maps

Conformational space annealing (CSA), a global optimization method, has been applied to various protein structure modeling tasks. In this paper, we applied CSA to the cryo-EM structure modeling task by combining the python subroutine of CSA (PyCSA) and the fast relax (FastRelax) protocol of PyRosetta. Refinement of initial structures generated from two methods, rigid fitting of predicted structures to the Cryo-EM map and de novo protein modeling by tracing the Cryo-EM map, was performed by CSA. In the refinement of the rigid-fitted structures, the final models showed that CSA can generate reliable atomic structures of proteins, even when large movements of protein domains were required. In the de novo modeling case, although the overall structural qualities of the final models were rather dependent on the initial models, the final models generated by CSA showed improved MolProbity scores and cross-correlation coefficients to the maps. These results suggest that CSA can accomplish flexible fitting and refinement together by sampling diverse conformations effectively and thus can be utilized for cryo-EM structure modeling.

Cryo-EM structure of the human cardiac myosin filament

Pumping of the heart is powered by filaments of the motor protein myosin that pull on actin filaments to generate cardiac contraction. In addition to myosin, the filaments contain cardiac myosin-binding protein C (cMyBP-C), which modulates contractility in response to physiological stimuli, and titin, which functions as a scaffold for filament assembly1. Myosin, cMyBP-C and titin are all subject to mutation, which can lead to heart failure. Despite the central importance of cardiac myosin filaments to life, their molecular structure has remained a mystery for 60 years2. Here we solve the structure of the main (cMyBP-C-containing) region of the human cardiac filament using cryo-electron microscopy. The reconstruction reveals the architecture of titin and cMyBP-C and shows how myosin's motor domains (heads) form three different types of motif (providing functional flexibility), which interact with each other and with titin and cMyBP-C to dictate filament architecture and function. The packing of myosin tails in the filament backbone is also resolved. The structure suggests how cMyBP-C helps to generate the cardiac super-relaxed state3; how titin and cMyBP-C may contribute to length-dependent activation4; and how mutations in myosin and cMyBP-C might disturb interactions, causing disease5,6. The reconstruction resolves past uncertainties and integrates previous data on cardiac muscle structure and function. It provides a new paradigm for interpreting structural, physiological and clinical observations, and for the design of potential therapeutic drugs.

Advancing cryo-electron microscopy data analysis through accelerated simulation-based flexible fitting approaches

Flexible fitting based on molecular dynamics simulation is a technique for structure modeling from cryo-EM data. It has been utilized for nearly two decades, and while cryo-EM resolution has improved significantly, it remains a powerful approach that can provide structural and dynamical insights that are not directly accessible from experimental data alone. Molecular dynamics simulations provide a means to extract atomistic details of conformational changes that are encoded in cryo-EM data and can also assist in improving the quality of structural models. Additionally, molecular dynamics simulations enable the characterization of conformational heterogeneity in cryo-EM data. We will summarize the advancements made in these techniques and highlight recent developments in this field.

Structure of the peroxisomal Pex1/Pex6 ATPase complex bound to a substrate

The double-ring AAA+ ATPase Pex1/Pex6 is required for peroxisomal receptor recycling and is essential for peroxisome formation. Pex1/Pex6 mutations cause severe peroxisome associated developmental disorders. Despite its pathophysiological importance, mechanistic details of the heterohexamer are not yet available. Here, we report cryoEM structures of Pex1/Pex6 from Saccharomyces cerevisiae, with an endogenous protein substrate trapped in the central pore of the catalytically active second ring (D2). Pairs of Pex1/Pex6(D2) subdomains engage the substrate via a staircase of pore-1 loops with distinct properties. The first ring (D1) is catalytically inactive but undergoes significant conformational changes resulting in alternate widening and narrowing of its pore. These events are fueled by ATP hydrolysis in the D2 ring and disengagement of a "twin-seam" Pex1/Pex6(D2) heterodimer from the staircase. Mechanical forces are propagated in a unique manner along Pex1/Pex6 interfaces that are not available in homo-oligomeric AAA-ATPases. Our structural analysis reveals the mechanisms of how Pex1 and Pex6 coordinate to achieve substrate translocation.

Structure and dynamics of an archetypal DNA nanoarchitecture revealed via cryo-EM and molecular dynamics simulations

DNA can be folded into rationally designed, unique, and functional materials. To fully realise the potential of these DNA materials, a fundamental understanding of their structure and dynamics is necessary, both in simple solvents as well as more complex and diverse anisotropic environments. Here we analyse an archetypal six-duplex DNA nanoarchitecture with single-particle cryo-electron microscopy and molecular dynamics simulations in solvents of tunable ionic strength and within the anisotropic environment of biological membranes. Outside lipid bilayers, the six-duplex bundle lacks the designed symmetrical barrel-type architecture. Rather, duplexes are arranged in non-hexagonal fashion and are disorted to form a wider, less elongated structure. Insertion into lipid membranes, however, restores the anticipated barrel shape due to lateral duplex compression by the bilayer. The salt concentration has a drastic impact on the stability of the inserted barrel-shaped DNA nanopore given the tunable electrostatic repulsion between the negatively charged duplexes. By synergistically combining experiments and simulations, we increase fundamental understanding into the environment-dependent structural dynamics of a widely used nanoarchitecture. This insight will pave the way for future engineering and biosensing applications.

Atomic model for core modifying region of human fatty acid synthase in complex with Denifanstat

Fatty acid synthase (FASN) catalyzes the de novo synthesis of palmitate, a 16-carbon chain fatty acid that is the primary precursor of lipid metabolism and an important intracellular signaling molecule. FASN is an attractive drug target in diabetes, cancer, fatty liver diseases, and viral infections. Here, we develop an engineered full-length human FASN (hFASN) that enables isolation of the condensing and modifying regions of the protein post-translation. The engineered protein enables electron cryo-microscopy (cryoEM) structure determination of the core modifying region of hFASN to 2.7 Å resolution. Examination of the dehydratase dimer within this region reveals that unlike its close homolog, porcine FASN, the catalytic cavity is close-ended and is accessible only through one opening in the vicinity of the active site. The core modifying region exhibits two major global conformational variabilities that describe long-range bending and twisting motions of the complex in solution. Finally, we solved the structure of this region bound to an anti-cancer drug, Denifanstat (i.e., TVB-2640), demonstrating the utility of our approach as a platform for structure guided design of future hFASN small molecule inhibitors.

Cryo-EM structure of the human cardiac myosin filament

Pumping of the heart is powered by filaments of the motor protein myosin, which pull on actin filaments to generate cardiac contraction. In addition to myosin, the filaments contain cardiac myosin-binding protein C (cMyBP-C), which modulates contractility in response to physiological stimuli, and titin, which functions as a scaffold for filament assembly 1 . Myosin, cMyBP-C and titin are all subject to mutation, which can lead to heart failure. Despite the central importance of cardiac myosin filaments to life, their molecular structure has remained a mystery for 60 years 2 . Here, we have solved the structure of the main (cMyBP-C-containing) region of the human cardiac filament to 6 Å resolution by cryo-EM. The reconstruction reveals the architecture of titin and cMyBP-C for the first time, and shows how myosin's motor domains (heads) form 3 different types of motif (providing functional flexibility), which interact with each other and with specific domains of titin and cMyBP-C to dictate filament architecture and regulate function. A novel packing of myosin tails in the filament backbone is also resolved. The structure suggests how cMyBP-C helps generate the cardiac super-relaxed state 3 , how titin and cMyBP-C may contribute to length-dependent activation 4 , and how mutations in myosin and cMyBP-C might disrupt interactions, causing disease 5, 6 . A similar structure is likely in vertebrate skeletal myosin filaments. The reconstruction resolves past uncertainties, and integrates previous data on cardiac muscle structure and function. It provides a new paradigm for interpreting structural, physiological and clinical observations, and for the design of potential therapeutic drugs.

EGOC inhibits TOROID polymerization by structurally activating TORC1

Target of rapamycin complex 1 (TORC1) is a protein kinase controlling cell homeostasis and growth in response to nutrients and stresses. In Saccharomyces cerevisiae, glucose depletion triggers a redistribution of TORC1 from a dispersed localization over the vacuole surface into a large, inactive condensate called TOROID (TORC1 organized in inhibited domains). However, the mechanisms governing this transition have been unclear. Here, we show that acute depletion and repletion of EGO complex (EGOC) activity is sufficient to control TOROID distribution, independently of other nutrient-signaling pathways. The 3.9-Å-resolution structure of TORC1 from TOROID cryo-EM data together with interrogation of key interactions in vivo provide structural insights into TORC1-TORC1' and TORC1-EGOC interaction interfaces. These data support a model in which glucose-dependent activation of EGOC triggers binding to TORC1 at an interface required for TOROID assembly, preventing TORC1 polymerization and promoting release of active TORC1.

Insights into the structure-function relationship of the NorQ/NorD chaperones from Paracoccus denitrificans reveal shared principles of interacting MoxR AAA+/VWA domain proteins

BACKGROUND

NorQ, a member of the MoxR-class of AAA+ ATPases, and NorD, a protein containing a Von Willebrand Factor Type A (VWA) domain, are essential for non-heme iron (Fe_B) cofactor insertion into cytochrome c-dependent nitric oxide reductase (cNOR). cNOR catalyzes NO reduction, a key step of bacterial denitrification. This work aimed at elucidating the specific mechanism of NorQD-catalyzed Fe_B insertion, and the general mechanism of the MoxR/VWA interacting protein families.

RESULTS

We show that NorQ-catalyzed ATP hydrolysis, an intact VWA domain in NorD, and specific surface carboxylates on cNOR are all features required for cNOR activation. Supported by BN-PAGE, low-resolution cryo-EM structures of NorQ and the NorQD complex show that NorQ forms a circular hexamer with a monomer of NorD binding both to the side and to the central pore of the NorQ ring. Guided by AlphaFold predictions, we assign the density that "plugs" the NorQ ring pore to the VWA domain of NorD with a protruding "finger" inserting through the pore and suggest this binding mode to be general for MoxR/VWA couples.

CONCLUSIONS

Based on our results, we present a tentative model for the mechanism of NorQD-catalyzed cNOR remodeling and suggest many of its features to be applicable to the whole MoxR/VWA family.

Visualization of translation and protein biogenesis at the ER membrane

The dynamic ribosome-translocon complex, which resides at the endoplasmic reticulum (ER) membrane, produces a major fraction of the human proteome^1,2. It governs the synthesis, translocation, membrane insertion, N-glycosylation, folding and disulfide-bond formation of nascent proteins. Although individual components of this machinery have been studied at high resolution in isolation^3-7, insights into their interplay in the native membrane remain limited. Here we use cryo-electron tomography, extensive classification and molecular modelling to capture snapshots of mRNA translation and protein maturation at the ER membrane at molecular resolution. We identify a highly abundant classical pre-translocation intermediate with eukaryotic elongation factor 1a (eEF1a) in an extended conformation, suggesting that eEF1a may remain associated with the ribosome after GTP hydrolysis during proofreading. At the ER membrane, distinct polysomes bind to different ER translocons specialized in the synthesis of proteins with signal peptides or multipass transmembrane proteins with the translocon-associated protein complex (TRAP) present in both. The near-complete atomic model of the most abundant ER translocon variant comprising the protein-conducting channel SEC61, TRAP and the oligosaccharyltransferase complex A (OSTA) reveals specific interactions of TRAP with other translocon components. We observe stoichiometric and sub-stoichiometric cofactors associated with OSTA, which are likely to include protein isomerases. In sum, we visualize ER-bound polysomes with their coordinated downstream machinery.

Local Normal Mode Analysis for Fast Loop Conformational Sampling

We propose and validate a novel method to efficiently explore local protein loop conformations based on a new formalism for constrained normal mode analysis (NMA) in internal coordinates. The manifold of possible loop configurations imposed by the position and orientation of the fixed loop ends is reduced to an orthogonal set of motions (or modes) encoding concerted rotations of all the backbone dihedral angles. We validate the sampling power on a set of protein loops with highly variable experimental structures and demonstrate that our approach can efficiently explore the conformational space of closed loops. We also show an acceptable resemblance of the ensembles around equilibrium conformations generated by long molecular simulations and constrained NMA on a set of exposed and diverse loops. In comparison with other methods, the main advantage is the lack of restrictions on the number of dihedrals that can be altered simultaneously. Furthermore, the method is computationally efficient since it only requires the diagonalization of a tiny matrix, and the modes of motions are energetically contextualized by the elastic network model, which includes both the loop and the neighboring residues.

MarkovFit: Structure Fitting for Protein Complexes in Electron Microscopy Maps Using Markov Random Field

An increasing number of protein complex structures are determined by cryo-electron microscopy (cryo-EM). When individual protein structures have been determined and are available, an important task in structure modeling is to fit the individual structures into the density map. Here, we designed a method that fits the atomic structures of proteins in cryo-EM maps of medium to low resolutions using Markov random fields, which allows probabilistic evaluation of fitted models. The accuracy of our method, MarkovFit, performed better than existing methods on datasets of 31 simulated cryo-EM maps of resolution 10 Å , nine experimentally determined cryo-EM maps of resolution less than 4 Å , and 28 experimentally determined cryo-EM maps of resolution 6 to 20 Å .

Protein Structural Modeling for Electron Microscopy Maps Using VESPER and MAINMAST

An increasing number of protein structures are determined by cryo-electron microscopy (cryo-EM) and stored in the Electron Microscopy Data Bank (EMDB). To interpret determined cryo-EM maps, several methods have been developed that model the tertiary structure of biomolecules, particularly proteins. Here we show how to use two such methods, VESPER and MAINMAST, which were developed in our group. VESPER is a method mainly for two purposes: fitting protein structure models into an EM map and aligning two EM maps locally or globally to capture their similarity. VESPER represents each EM map as a set of vectors pointing toward denser points. By considering matching the directions of vectors, in general, VESPER aligns maps better than conventional methods that only consider local densities of maps. MAINMAST is a de novo protein modeling tool designed for EM maps with resolution of 3-5 Å or better. MAINMAST builds a protein main chain directly from a density map by tracing dense points in an EM map and connecting them using a tree-graph structure. This article describes how to use these two tools using three illustrative modeling examples. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Protein structure model fitting using VESPER Alternate Protocol: Atomic model fitting using VESPER web server Basic Protocol 2: Protein de novo modeling using MAINMAST.

CR-I-TASSER: assemble protein structures from cryo-EM density maps using deep convolutional neural networks

Cryo-electron microscopy (cryo-EM) has become a leading approach for protein structure determination, but it remains challenging to accurately model atomic structures with cryo-EM density maps. We propose a hybrid method, CR-I-TASSER, which integrates deep neural-network learning with I-TASSER assembly simulations for automated cryo-EM structure determination. The method is benchmarked on 778 proteins with simulated and experimental density maps, where CR-I-TASSER constructs models with a correct fold (TM-score>0.5) for 643 targets that is 64% higher than the best of other de novo and refinement-based approaches on high-resolution data samples. Detailed data analyses showed that the major advantage of CR-I-TASSER lies in the deep-learning based Cα position prediction, which significantly improves the threading template quality and therefore boosts the accuracy of final models through optimized fragment assembly simulations. These results demonstrate a new avenue to determine cryo-EM protein structures with high accuracy and robustness covering various target types and density-map resolutions.

NMMD: Efficient cryo-EM flexible fitting based on simultaneous Normal Mode and Molecular Dynamics atomic displacements

Atomic models of cryo electron microscopy (cryo-EM) maps of biomolecular conformations are often obtained by flexible fitting of the maps with available atomic structures of other conformations (e.g., obtained by X-ray crystallography). This article presents a new flexible fitting method, NMMD, which combines normal mode analysis (NMA) and molecular dynamics simulation (MD). Given an atomic structure and a cryo-EM map to fit, NMMD simultaneously estimates global atomic displacements based on NMA and local displacements based on MD. NMMD was implemented by modifying EMfit, a flexible fitting method using MD only, in GENESIS 1.4. As EMfit, NMMD can be run with replica exchange umbrella sampling procedure. The new method was tested using a variety of EM maps (synthetic and experimental, with different noise levels and resolutions). The results of the tests show that adding normal modes to MD-based fitting makes the fitting faster (40% in average) and, in the majority of cases, more accurate.

Structural basis and mode of action for two broadly neutralizing antibodies against SARS-CoV-2 emerging variants of concern

Emerging variants of concern for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can transmit more efficiently and partially evade protective immune responses, thus necessitating continued refinement of antibody therapies and immunogen design. Here, we elucidate the structural basis and mode of action for two potent SARS-CoV-2 spike (S)-neutralizing monoclonal antibodies, CV3-1 and CV3-25, which remain effective against emerging variants of concern in vitro and in vivo. CV3-1 binds to the (485-GFN-487) loop within the receptor-binding domain (RBD) in the "RBD-up" position and triggers potent shedding of the S1 subunit. In contrast, CV3-25 inhibits membrane fusion by binding to an epitope in the stem helix region of the S2 subunit that is highly conserved among β-coronaviruses. Thus, vaccine immunogen designs that incorporate the conserved regions in the RBD and stem helix region are candidates to elicit pan-coronavirus protective immune responses.

Structural mechanism for tyrosine hydroxylase inhibition by dopamine and reactivation by Ser40 phosphorylation

Tyrosine hydroxylase (TH) catalyzes the rate-limiting step in the biosynthesis of dopamine (DA) and other catecholamines, and its dysfunction leads to DA deficiency and parkinsonisms. Inhibition by catecholamines and reactivation by S40 phosphorylation are key regulatory mechanisms of TH activity and conformational stability. We used Cryo-EM to determine the structures of full-length human TH without and with DA, and the structure of S40 phosphorylated TH, complemented with biophysical and biochemical characterizations and molecular dynamics simulations. TH presents a tetrameric structure with dimerized regulatory domains that are separated 15 Å from the catalytic domains. Upon DA binding, a 20-residue α-helix in the flexible N-terminal tail of the regulatory domain is fixed in the active site, blocking it, while S40-phosphorylation forces its egress. The structures reveal the molecular basis of the inhibitory and stabilizing effects of DA and its counteraction by S40-phosphorylation, key regulatory mechanisms for homeostasis of DA and TH.

Tyrosine hydroxylase (TH) catalyzes the rate-limiting step in the synthesis of the catecholamine neurotransmitters and hormones dopamine (DA), adrenaline and noradrenaline. Here, the authors present the cryo-EM structures of full-length human TH in the apo form and bound with DA, as well as the structure of Ser40 phosphorylated TH, and discuss the inhibitory and stabilizing effects of DA on TH and its counteraction by Ser40-phosphorylation.

Using cryo-EM to uncover mechanisms of bacterial transcriptional regulation

Transcription is the principal control point for bacterial gene expression, and it enables a global cellular response to an intracellular or environmental trigger. Transcriptional regulation is orchestrated by transcription factors, which activate or repress transcription of target genes by modulating the activity of RNA polymerase. Dissecting the nature and precise choreography of these interactions is essential for developing a molecular understanding of transcriptional regulation. While the contribution of X-ray crystallography has been invaluable, the 'resolution revolution' of cryo-electron microscopy has transformed our structural investigations, enabling large, dynamic and often transient transcription complexes to be resolved that in many cases had resisted crystallisation. In this review, we highlight the impact cryo-electron microscopy has had in gaining a deeper understanding of transcriptional regulation in bacteria. We also provide readers working within the field with an overview of the recent innovations available for cryo-electron microscopy sample preparation and image reconstruction of transcription complexes.

Integrative structural modeling of macromolecular complexes using Assembline

Integrative modeling enables structure determination of macromolecular complexes by combining data from multiple experimental sources such as X-ray crystallography, electron microscopy or cross-linking mass spectrometry. It is particularly useful for complexes not amenable to high-resolution electron microscopy-complexes that are flexible, heterogeneous or imaged in cells with cryo-electron tomography. We have recently developed an integrative modeling protocol that allowed us to model multi-megadalton complexes as large as the nuclear pore complex. Here, we describe the Assembline software package, which combines multiple programs and libraries with our own algorithms in a streamlined modeling pipeline. Assembline builds ensembles of models satisfying data from atomic structures or homology models, electron microscopy maps and other experimental data, and provides tools for their analysis. Compared with other methods, Assembline enables efficient sampling of conformational space through a multistep procedure, provides new modeling restraints and includes a unique configuration system for setting up the modeling project. Our protocol achieves exhaustive sampling in less than 100-1,000 CPU-hours even for complexes in the megadalton range. For larger complexes, resources available in institutional or public computer clusters are needed and sufficient to run the protocol. We also provide step-by-step instructions for preparing the input, running the core modeling steps and assessing modeling performance at any stage.

Mechanism of actin-dependent activation of nucleotidyl cyclase toxins from bacterial human pathogens

Bacterial human pathogens secrete initially inactive nucleotidyl cyclases that become potent enzymes by binding to actin inside eukaryotic host cells. The underlying molecular mechanism of this activation is, however, unclear. Here, we report structures of ExoY from Pseudomonas aeruginosa and Vibrio vulnificus bound to their corresponding activators F-actin and profilin-G-actin. The structures reveal that in contrast to the apo-state, two flexible regions become ordered and interact strongly with actin. The specific stabilization of these regions results in an allosteric stabilization of the nucleotide binding pocket and thereby to an activation of the enzyme. Differences in the sequence and conformation of the actin-binding regions are responsible for the selective binding to either F- or G-actin. Other nucleotidyl cyclase toxins that bind to calmodulin rather than actin undergo a similar disordered-to-ordered transition during activation, suggesting that the allosteric activation-by-stabilization mechanism of ExoY is conserved in these enzymes, albeit the different activator.

Flexible open conformation of the AP-3 complex explains its role in cargo recruitment at the Golgi

Vesicle formation at endomembranes requires the selective concentration of cargo by coat proteins. Conserved adapter protein complexes at the Golgi (AP-3), the endosome (AP-1), or the plasma membrane (AP-2) with their conserved core domain and flexible ear domains mediate this function. These complexes also rely on the small GTPase Arf1 and/or specific phosphoinositides for membrane binding. The structural details that influence these processes, however, are still poorly understood. Here we present cryo-EM structures of the full-length stable 300 kDa yeast AP-3 complex. The structures reveal that AP-3 adopts an open conformation in solution, comparable to the membrane-bound conformations of AP-1 or AP-2. This open conformation appears to be far more flexible than AP-1 or AP-2, resulting in compact, intermediate, and stretched subconformations. Mass spectrometrical analysis of the cross-linked AP-3 complex further indicates that the ear domains are flexibly attached to the surface of the complex. Using biochemical reconstitution assays, we also show that efficient AP-3 recruitment to the membrane depends primarily on cargo binding. Once bound to cargo, AP-3 clustered and immobilized cargo molecules, as revealed by single-molecule imaging on polymer-supported membranes. We conclude that its flexible open state may enable AP-3 to bind and collect cargo at the Golgi and could thus allow coordinated vesicle formation at the trans-Golgi upon Arf1 activation.

Biasing RNA Coarse-Grained Folding Simulations with Small-Angle X-ray Scattering Data

RNA molecules can easily adopt alternative structures in response to different environmental conditions. As a result, a molecule's energy landscape is rough and can exhibit a multitude of deep basins. In the absence of a high-resolution structure, small-angle X-ray scattering data (SAXS) can narrow down the conformational space available to the molecule and be used in conjunction with physical modeling to obtain high-resolution putative structures to be further tested by experiments. Because of the low resolution of these data, it is natural to implement the integration of SAXS data into simulations using a coarse-grained representation of the molecule, allowing for much wider searches and faster evaluation of SAXS theoretical intensity curves than with atomistic models. We present here the theoretical framework and the implementation of a simulation approach based on our coarse-grained model HiRE-RNA combined with SAXS evaluations "on-the-fly" leading the simulation toward conformations agreeing with the scattering data, starting from partially folded structures as the ones that can easily be obtained from secondary structure prediction-based tools. We show on three benchmark systems how our approach can successfully achieve high-resolution structures with remarkable similarity with the native structure recovering not only the overall shape, as imposed by SAXS data, but also the details of initially missing base pairs.

Mutational analysis to explore long-range allosteric couplings involved in a pentameric channel receptor pre-activation and activation

Pentameric ligand-gated ion channels (pLGICs) mediate chemical signaling through a succession of allosteric transitions that are yet not completely understood as intermediate states remain poorly characterized by structural approaches. In a previous study on the prototypic bacterial proton-gated channel GLIC, we generated several fluorescent sensors of the protein conformation that report a fast transition to a pre-active state, which precedes the slower process of activation with pore opening. Here, we explored the phenotype of a series of allosteric mutations, using simultaneous steady-state fluorescence and electrophysiological measurements over a broad pH range. Our data, fitted to a three-state Monod-Wyman-Changeux model, show that mutations at the subunit interface in the extracellular domain (ECD) principally alter pre-activation, while mutations in the lower ECD and in the transmembrane domain principally alter activation. We also show that propofol alters both transitions. Data are discussed in the framework of transition pathways generated by normal mode analysis (iModFit). It further supports that pre-activation involves major quaternary compaction of the ECD, and suggests that activation involves principally a reorganization of a ‘central gating region’ involving a contraction of the ECD β-sandwich and the tilt of the channel lining M2 helix.

Unwinding of a DNA replication fork by a hexameric viral helicase

Hexameric helicases are motor proteins that unwind double-stranded DNA (dsDNA) during DNA replication but how they are optimised for strand separation is unclear. Here we present the cryo-EM structure of the full-length E1 helicase from papillomavirus, revealing all arms of a bound DNA replication fork and their interactions with the helicase. The replication fork junction is located at the entrance to the helicase collar ring, that sits above the AAA + motor assembly. dsDNA is escorted to and the 5´ single-stranded DNA (ssDNA) away from the unwinding point by the E1 dsDNA origin binding domains. The 3´ ssDNA interacts with six spirally-arranged β-hairpins and their cyclical top-to-bottom movement pulls the ssDNA through the helicase. Pulling of the RF against the collar ring separates the base-pairs, while modelling of the conformational cycle suggest an accompanying movement of the collar ring has an auxiliary role, helping to make efficient use of ATP in duplex unwinding.

Replicative hexameric helicases are fundamental components of replisomes. Here the authors resolve a cryo-EM structure of the E1 helicase from papillomavirus bound to a DNA replication fork, providing insights into the mechanism of DNA unwinding by these hexameric enzymes.

Structural analysis of receptors and actin polarity in platelet protrusions

During activation the platelet cytoskeleton is reorganized, inducing adhesion to the extracellular matrix and cell spreading. These processes are critical for wound healing and clot formation. Initially, this task relies on the formation of strong cellular-extracellular matrix interactions, exposed in subendothelial lesions. Despite the medical relevance of these processes, there is a lack of high-resolution structural information on the platelet cytoskeleton controlling cell spreading and adhesion. Here, we present in situ structural analysis of membrane receptors and the underlying cytoskeleton in platelet protrusions by applying cryoelectron tomography to intact platelets. We utilized three-dimensional averaging procedures to study receptors at the plasma membrane. Analysis of substrate interaction-free receptors yielded one main structural class resolved to 26 Å, resembling the α_IIbβ₃ integrin folded conformation. Furthermore, structural analysis of the actin network in pseudopodia indicates a nonuniform polarity of filaments. This organization would allow generation of the contractile forces required for integrin-mediated cell adhesion.

Parameterizing elastic network models to capture the dynamics of proteins

Coarse-grained normal mode analyses of protein dynamics rely on the idea that the geometry of a protein structure contains enough information for computing its fluctuations around its equilibrium conformation. This geometry is captured in the form of an elastic network (EN), namely a network of edges between its residues. The normal modes of a protein are then identified with the normal modes of its EN. Different approaches have been proposed to construct ENs, focusing on the choice of the edges that they are comprised of, and on their parameterizations by the force constants associated with those edges. Here we propose new tools to guide choices on these two facets of EN. We study first different geometric models for ENs. We compare cutoff-based ENs, whose edges have lengths that are smaller than a cutoff distance, with Delaunay-based ENs and find that the latter provide better representations of the geometry of protein structures. We then derive an analytical method for the parameterization of the EN such that its dynamics leads to atomic fluctuations that agree with experimental B-factors. To limit overfitting, we attach a parameter referred to as flexibility constant to each atom instead of to each edge in the EN. The parameterization is expressed as a non-linear optimization problem whose parameters describe both rigid-body and internal motions. We show that this parameterization leads to improved ENs, whose dynamics mimic MD simulations better than ENs with uniform force constants, and reduces the number of normal modes needed to reproduce functional conformational changes.

GluN2A and GluN2B NMDA receptors use distinct allosteric routes

Allostery represents a fundamental mechanism of biological regulation that involves long-range communication between distant protein sites. It also provides a powerful framework for novel therapeutics. NMDA receptors (NMDARs), glutamate-gated ionotropic receptors that play central roles in synapse maturation and plasticity, are prototypical allosteric machines harboring large extracellular N-terminal domains (NTDs) that provide allosteric control of key receptor properties with impact on cognition and behavior. It is commonly thought that GluN2A and GluN2B receptors, the two predominant NMDAR subtypes in the adult brain, share similar allosteric transitions. Here, combining functional and structural interrogation, we reveal that GluN2A and GluN2B receptors utilize different long-distance allosteric mechanisms involving distinct subunit-subunit interfaces and molecular rearrangements. NMDARs have thus evolved multiple levels of subunit-specific allosteric control over their transmembrane ion channel pore. Our results uncover an unsuspected diversity in NMDAR molecular mechanisms with important implications for receptor physiology and precision drug development.

NMDA receptors are glutamate-gated ion channels essential for synapse maturation and plasticity. Here the authors show that GluN2A and GluN2B NMDA receptors — the two principal subtypes NMDARs in the adult CNS — operate through distinct long range allosteric mechanisms.

Structural Basis and Mode of Action for Two Broadly Neutralizing Antibodies Against SARS-CoV-2 Emerging Variants of Concern

Emerging variants of concern for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can transmit more efficiently and partially evade protective immune responses, thus necessitating continued refinement of antibody therapies and immunogen design. Here we elucidate the structural basis and mode of action for two potent SARS-CoV-2 Spike (S) neutralizing monoclonal antibodies CV3-1 and CV3-25 that remained effective against emerging variants of concern in vitro and in vivo. CV3-1 bound to the (485-GFN-487) loop within the receptor-binding domain (RBD) in the "RBD-up" position and triggered potent shedding of the S1 subunit. In contrast, CV3-25 inhibited membrane fusion by binding to an epitope in the stem helix region of the S2 subunit that is highly conserved among β-coronaviruses. Thus, vaccine immunogen designs that incorporate the conserved regions in RBD and stem helix region are candidates to elicit pan-coronavirus protective immune responses.

Study of the DnaB:DciA interplay reveals insights into the primary mode of loading of the bacterial replicative helicase

Replicative helicases are essential proteins that unwind DNA in front of replication forks. Their loading depends on accessory proteins and in bacteria, DnaC and DnaI are well characterized loaders. However, most bacteria do not express either of these two proteins. Instead, they are proposed to rely on DciA, an ancestral protein unrelated to DnaC/I. While the DciA structure from Vibrio cholerae shares no homology with DnaC, it reveals similarities with DnaA and DnaX, two proteins involved during replication initiation. As other bacterial replicative helicases, VcDnaB adopts a toroid-shaped homo-hexameric structure, but with a slightly open dynamic conformation in the free state. We show that VcDnaB can load itself on DNA in vitro and that VcDciA stimulates this function, resulting in an increased DNA unwinding. VcDciA interacts with VcDnaB with a 3/6 stoichiometry and we show that a determinant residue, which discriminates DciA- and DnaC/I-helicases, is critical in vivo. Our work is the first step toward the understanding of the ancestral mode of loading of bacterial replicative helicases on DNA. It sheds light on the strategy employed by phage helicase loaders to hijack bacterial replicative helicases and may explain the recurrent domestication of dnaC/I through evolution in bacteria.

Integrative Structural Biology in the Era of Accurate Structure Prediction

Characterizing the three-dimensional structure of macromolecules is central to understanding their function. Traditionally, structures of proteins and their complexes have been determined using experimental techniques such as X-ray crystallography, NMR, or cryo-electron microscopy-applied individually or in an integrative manner. Meanwhile, however, computational methods for protein structure prediction have been improving their accuracy, gradually, then suddenly, with the breakthrough advance by AlphaFold2, whose models of monomeric proteins are often as accurate as experimental structures. This breakthrough foreshadows a new era of computational methods that can build accurate models for most monomeric proteins. Here, we envision how such accurate modeling methods can combine with experimental structural biology techniques, enhancing integrative structural biology. We highlight the challenges that arise when considering multiple structural conformations, protein complexes, and polymorphic assemblies. These challenges will motivate further developments, both in modeling programs and in methods to solve experimental structures, towards better and quicker investigation of structure-function relationships.

Efficient Flexible Fitting Refinement with Automatic Error Fixing for De Novo Structure Modeling from Cryo-EM Density Maps

Structural modeling of proteins from cryo-electron microscopy (cryo-EM) density maps is one of the challenging issues in structural biology. De novo modeling combined with flexible fitting refinement (FFR) has been widely used to build a structure of new proteins. In de novo prediction, artificial conformations containing local structural errors such as chirality errors, cis peptide bonds, and ring penetrations are frequently generated and cannot be easily removed in the subsequent FFR. Moreover, refinement can be significantly suppressed due to the low mobility of atoms inside the protein. To overcome these problems, we propose an efficient scheme for FFR, in which the local structural errors are fixed first, followed by FFR using an iterative simulated annealing (SA) molecular dynamics protocol with the united atom (UA) model in an implicit solvent model; we call this scheme "SAUA-FFR". The best model is selected from multiple flexible fitting runs with various biasing force constants to reduce overfitting. We apply our scheme to the decoys obtained from MAINMAST and demonstrate an improvement of the best model of eight selected proteins in terms of the root-mean-square deviation, MolProbity score, and RWplus score compared to the original scheme of MAINMAST. Fixing the local structural errors can enhance the formation of secondary structures, and the UA model enables progressive refinement compared to the all-atom model owing to its high mobility in the implicit solvent. The SAUA-FFR scheme realizes efficient and accurate protein structure modeling from medium-resolution maps with less overfitting.

Induced forms of α2-macroglobulin neutralize heparin and direct oral anticoagulant effects

Alpha₂-macroglobulin (α₂M) is a physiological macromolecule that facilitates the clearance of many proteinases, cytokines and growth factors in human. Here, we explored the effect of induced forms of α₂M on anticoagulant drugs. Gla-domainless factor Xa (GDFXa) and methylamine (MA)-induced α₂M were prepared and characterized by electrophoresis, immunonephelometry, chromogenic, clot waveform and rotational thromboelastometry assays. Samples from healthy volunteers and anticoagulated patients were included. In vivo neutralization of anticoagulants was evaluated in C57Bl/6JRj mouse bleeding-model. Anticoagulant binding sites on induced α₂M were depicted by computer-aided energy minimization modeling. GDFXa-induced α₂M neutralized dabigatran and heparins in plasma and whole blood. In mice, a single IV dose of GDFXa-induced α₂M following anticoagulant administration significantly reduced blood loss and bleeding time. Being far easier to prepare, we investigated the efficacy of MA-induced α₂M. It neutralized rivaroxaban, apixaban, dabigatran and heparins in spiked samples in a concentration-dependent manner and in samples from treated patients. Molecular docking analysis evidenced the ability of MA-induced α₂M to bind non-covalently these compounds via some deeply buried binding sites. Induced forms of α₂M have the potential to neutralize direct oral anticoagulants and heparins, and might be developed as a universal antidote in case of major bleeding or urgent surgery.

Structural basis of substrate recognition and thermal protection by a small heat shock protein

Small heat shock proteins (sHsps) bind unfolding proteins, thereby playing a pivotal role in the maintenance of proteostasis in virtually all living organisms. Structural elucidation of sHsp-substrate complexes has been hampered by the transient and heterogeneous nature of their interactions, and the precise mechanisms underlying substrate recognition, promiscuity, and chaperone activity of sHsps remain unclear. Here we show the formation of a stable complex between Arabidopsis thaliana plastid sHsp, Hsp21, and its natural substrate 1-deoxy-D-xylulose 5-phosphate synthase (DXPS) under heat stress, and report cryo-electron microscopy structures of Hsp21, DXPS and Hsp21-DXPS complex at near-atomic resolution. Monomeric Hsp21 binds across the dimer interface of DXPS and engages in multivalent interactions by recognizing highly dynamic structural elements in DXPS. Hsp21 partly unfolds its central α-crystallin domain to facilitate binding of DXPS, which preserves a native-like structure. This mode of interaction suggests a mechanism of sHsps anti-aggregation activity towards a broad range of substrates.

HOPMA: Boosting Protein Functional Dynamics with Colored Contact Maps

In light of the recent very rapid progress in protein structure prediction, accessing the multitude of functional protein states is becoming more central than ever before. Indeed, proteins are flexible macromolecules, and they often perform their function by switching between different conformations. However, high-resolution experimental techniques such as X-ray crystallography and cryogenic electron microscopy can catch relatively few protein functional states. Many others are only accessible under physiological conditions in solution. Therefore, there is a pressing need to fill this gap with computational approaches. We present HOPMA, a novel method to predict protein functional states and transitions by using a modified elastic network model. The method exploits patterns in a protein contact map, taking its 3D structure as input, and excludes some disconnected patches from the elastic network. Combined with nonlinear normal mode analysis, this strategy boosts the protein conformational space exploration, especially when the input structure is highly constrained, as we demonstrate on a set of more than 400 transitions. Our results let us envision the discovery of new functional conformations, which were unreachable previously, starting from the experimentally known protein structures. The method is computationally efficient and available at https://github.com/elolaine/HOPMA and https://team.inria.fr/nano-d/software/nolb-normal-modes.

Electron microscopy as a critical tool in the determination of pore forming mechanisms in proteins

Electron microscopy has consistently played an important role in the description of pore-forming protein systems. The discovery of pore-forming proteins has depended on visualization of the structural pores formed by their oligomeric protein complexes, and as electron microscopy has advanced technologically so has the degree of insight it has been able to give. This review considers a large number of published studies of pore-forming complexes in prepore and pore states determined using single-particle cryo-electron microscopy. Sample isolation and preparation, imaging and image analysis, structure determination and optimization of results are all discussed alongside challenges which pore-forming proteins particularly present. The review also considers the use made of cryo-electron tomography to study pores within their membrane environment and which will prove an increasingly important approach for the future.

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

The eukaryotic chaperonin TRiC/CCT plays a major role in assisting the folding of many proteins through an ATP-driven allosteric cycle. Recent structures elucidated by cryo-electron microscopy provide a broad view of the conformations visited at various stages of the chaperonin cycle, including a sequential activation of its subunits in response to nucleotide binding. But we lack a thorough mechanistic understanding of the structure-based dynamics and communication properties that underlie the TRiC/CCT machinery. In this study, we present a computational methodology based on elastic network models adapted to cryo-EM density maps to gain a deeper understanding of the structure-encoded allosteric dynamics of this hexadecameric machine. We have analysed several structures of the chaperonin resolved in different states toward mapping its conformational landscape. Our study indicates that the overall architecture intrinsically favours cooperative movements that comply with the structural variabilities observed in experiments. Furthermore, the individual subunits CCT1-CCT8 exhibit state-dependent sequential events at different states of the allosteric cycle. For example, in the ATP-bound state, subunits CCT5 and CCT4 selectively initiate the lid closure motions favoured by the overall architecture; whereas in the apo form of the heteromer, the subunit CCT7 exhibits the highest predisposition to structural change. The changes then propagate through parallel fluxes of allosteric signals to neighbours on both rings. The predicted state-dependent mechanisms of sequential activation provide new insights into TRiC/CCT intra- and inter-ring signal transduction events.

JARID2 and AEBP2 regulate PRC2 in the presence of H2AK119ub1 and other histone modifications

Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) cooperate to determine cell identity by epigenetic gene expression regulation. However, the mechanism of PRC2 recruitment by means of recognition of PRC1-mediated H2AK119ub1 remains poorly understood. Our PRC2 cryo-electron microscopy structure with cofactors JARID2 and AEBP2 bound to a H2AK119ub1-containing nucleosome reveals a bridge helix in EZH2 that connects the SET domain, H3 tail, and nucleosomal DNA. JARID2 and AEBP2 each interact with one ubiquitin and the H2A-H2B surface. JARID2 stimulates PRC2 through interactions with both the polycomb protein EED and the H2AK119-ubiquitin, whereas AEBP2 has an additional scaffolding role. The presence of these cofactors partially overcomes the inhibitory effect that H3K4me3 and H3K36me3 exert on core PRC2 (in the absence of cofactors). Our results support a key role for JARID2 and AEBP2 in the cross-talk between histone modifications and PRC2 activity.

High-Resolution Cryo-EM Structure of the Cardiac Actomyosin Complex

Heart contraction depends on a complicated array of interactions between sarcomeric proteins required to convert chemical energy into mechanical force. Cyclic interactions between actin and myosin molecules, controlled by troponin and tropomyosin, generate the sliding force between the actin-based thin and myosin-based thick filaments. Alterations in this sophisticated system due to missense mutations can lead to cardiovascular diseases. Numerous structural studies proposed pathological mechanisms of missense mutations at the myosin-myosin, actin-tropomyosin, and tropomyosin-troponin interfaces. However, despite the central role of actomyosin interactions a detailed structural description of the cardiac actomyosin interface remained unknown. Here, we report a cryo-EM structure of a cardiac actomyosin complex at 3.8 Å resolution. The structure reveals the molecular basis of cardiac diseases caused by missense mutations in myosin and actin proteins.

Hybrid methods for combined experimental and computational determination of protein structure

Knowledge of protein structure is paramount to the understanding of biological function, developing new therapeutics, and making detailed mechanistic hypotheses. Therefore, methods to accurately elucidate three-dimensional structures of proteins are in high demand. While there are a few experimental techniques that can routinely provide high-resolution structures, such as x-ray crystallography, nuclear magnetic resonance (NMR), and cryo-EM, which have been developed to determine the structures of proteins, these techniques each have shortcomings and thus cannot be used in all cases. However, additionally, a large number of experimental techniques that provide some structural information, but not enough to assign atomic positions with high certainty have been developed. These methods offer sparse experimental data, which can also be noisy and inaccurate in some instances. In cases where it is not possible to determine the structure of a protein experimentally, computational structure prediction methods can be used as an alternative. Although computational methods can be performed without any experimental data in a large number of studies, inclusion of sparse experimental data into these prediction methods has yielded significant improvement. In this Perspective, we cover many of the successes of integrative modeling, computational modeling with experimental data, specifically for protein folding, protein-protein docking, and molecular dynamics simulations. We describe methods that incorporate sparse data from cryo-EM, NMR, mass spectrometry, electron paramagnetic resonance, small-angle x-ray scattering, Förster resonance energy transfer, and genetic sequence covariation. Finally, we highlight some of the major challenges in the field as well as possible future directions.

Structural basis of Focal Adhesion Kinase activation on lipid membranes

Focal adhesion kinase (FAK) is a key component of the membrane proximal signaling layer in focal adhesion complexes, regulating important cellular processes, including cell migration, proliferation, and survival. In the cytosol, FAK adopts an autoinhibited state but is activated upon recruitment into focal adhesions, yet how this occurs or what induces structural changes is unknown. Here, we employ cryo-electron microscopy to reveal how FAK associates with lipid membranes and how membrane interactions unlock FAK autoinhibition to promote activation. Intriguingly, initial binding of FAK to the membrane causes steric clashes that release the kinase domain from autoinhibition, allowing it to undergo a large conformational change and interact itself with the membrane in an orientation that places the active site toward the membrane. In this conformation, the autophosphorylation site is exposed and multiple interfaces align to promote FAK oligomerization on the membrane. We show that interfaces responsible for initial dimerization and membrane attachment are essential for FAK autophosphorylation and resulting cellular activity including cancer cell invasion, while stable FAK oligomerization appears to be needed for optimal cancer cell proliferation in an anchorage-independent manner. Together, our data provide structural details of a key membrane bound state of FAK that is primed for efficient autophosphorylation and activation, hence revealing the critical event in integrin mediated FAK activation and signaling at focal adhesions.

Continuous flexibility analysis of SARS-CoV-2 spike prefusion structures

Using a new consensus-based image-processing approach together with principal component analysis, the flexibility and conformational dynamics of the SARS-CoV-2 spike in the prefusion state have been analysed. These studies revealed concerted motions involving the receptor-binding domain (RBD), N-terminal domain, and subdomains 1 and 2 around the previously characterized 1-RBD-up state, which have been modeled as elastic deformations. It is shown that in this data set there are not well defined, stable spike conformations, but virtually a continuum of states. An ensemble map was obtained with minimum bias, from which the extremes of the change along the direction of maximal variance were modeled by flexible fitting. The results provide a warning of the potential image-processing classification instability of these complicated data sets, which has a direct impact on the interpretability of the results.