Quantitative analysis of cryo-EM density map segmentation by watershed and scale-space filtering, and fitting of structures by alignment to regions

The giant mimivirus 1.2 Mb genome is elegantly organized into a 30-nm diameter helical protein shield

Mimivirus is the prototype of the Mimiviridae family of giant dsDNA viruses. Little is known about the organization of the 1.2 Mb genome inside the membrane-limited nucleoid filling the ~0.5 µm icosahedral capsids. Cryo-electron microscopy, cryo-electron tomography, and proteomics revealed that it is encased into a ~30-nm diameter helical protein shell surprisingly composed of two GMC-type oxidoreductases, which also form the glycosylated fibrils decorating the capsid. The genome is arranged in 5- or 6-start left-handed super-helices, with each DNA-strand lining the central channel. This luminal channel of the nucleoprotein fiber is wide enough to accommodate oxidative stress proteins and RNA polymerase subunits identified by proteomics. Such elegant supramolecular organization would represent a remarkable evolutionary strategy for packaging and protecting the genome, in a state ready for immediate transcription upon unwinding in the host cytoplasm. The parsimonious use of the same protein in two unrelated substructures of the virion is unexpected for a giant virus with thousand genes at its disposal.

Structural analysis of the basal state of the Artemis:DNA-PKcs complex

Artemis nuclease and DNA-dependent protein kinase catalytic subunit (DNA-PKcs) are key components in nonhomologous DNA end joining (NHEJ), the major repair mechanism for double-strand DNA breaks. Artemis activation by DNA-PKcs resolves hairpin DNA ends formed during V(D)J recombination. Artemis deficiency disrupts development of adaptive immunity and leads to radiosensitive T- B- severe combined immunodeficiency (RS-SCID). An activated state of Artemis in complex with DNA-PK was solved by cryo-EM recently, which showed Artemis bound to the DNA. Here, we report that the pre-activated form (basal state) of the Artemis:DNA-PKcs complex is stable on an agarose-acrylamide gel system, and suitable for cryo-EM structural analysis. Structures show that the Artemis catalytic domain is dynamically positioned externally to DNA-PKcs prior to ABCDE autophosphorylation and show how both the catalytic and regulatory domains of Artemis interact with the N-HEAT and FAT domains of DNA-PKcs. We define a mutually exclusive binding site for Artemis and XRCC4 on DNA-PKcs and show that an XRCC4 peptide disrupts the Artemis:DNA-PKcs complex. All of the findings are useful in explaining how a hypomorphic L3062R missense mutation of DNA-PKcs could lead to insufficient Artemis activation, hence RS-SCID. Our results provide various target site candidates to design disruptors for Artemis:DNA-PKcs complex formation.

Bacteriophage protein PEIP is a potent Bacillus subtilis enolase inhibitor

Enolase is a highly conserved enzyme that presents in all organisms capable of glycolysis or fermentation. Its immediate product phosphoenolpyruvate is essential for other important processes like peptidoglycan synthesis and the phosphotransferase system in bacteria. Therefore, enolase inhibitors are of great interest. Here, we report that Gp60, a phage-encoded enolase inhibitor protein (PEIP) of bacteriophage SPO1 for Bacillus subtilis, is an enolase inhibitor. PEIP-expressing bacteria exhibit growth attenuation, thinner cell walls, and safranin color in Gram staining owing to impaired peptidoglycan synthesis. We solve the structure of PEIP-enolase tetramer and show that PEIP disassembles enolase by disrupting the basic dimer unit. The structure reveals that PEIP does not compete for substrate binding but induces a cascade of conformational changes that limit accessibility to the enolase catalytic site. This phage-inspired disassembly of enolase represents an alternative strategy for the development of anti-microbial drugs.

Membrane-anchored HDCR nanowires drive hydrogen-powered CO2 fixation

Filamentous enzymes have been found in all domains of life, but the advantage of filamentation is often elusive¹. Some anaerobic, autotrophic bacteria have an unusual filamentous enzyme for CO₂ fixation-hydrogen-dependent CO₂ reductase (HDCR)^2,3-which directly converts H₂ and CO₂ into formic acid. HDCR reduces CO₂ with a higher activity than any other known biological or chemical catalyst^4,5, and it has therefore gained considerable interest in two areas of global relevance: hydrogen storage and combating climate change by capturing atmospheric CO₂. However, the mechanistic basis of the high catalytic turnover rate of HDCR has remained unknown. Here we use cryo-electron microscopy to reveal the structure of a short HDCR filament from the acetogenic bacterium Thermoanaerobacter kivui. The minimum repeating unit is a hexamer that consists of a formate dehydrogenase (FdhF) and two hydrogenases (HydA2) bound around a central core of hydrogenase Fe-S subunits, one HycB3 and two HycB4. These small bacterial polyferredoxin-like proteins oligomerize through their C-terminal helices to form the backbone of the filament. By combining structure-directed mutagenesis with enzymatic analysis, we show that filamentation and rapid electron transfer through the filament enhance the activity of HDCR. To investigate the structure of HDCR in situ, we imaged T. kivui cells with cryo-electron tomography and found that HDCR filaments bundle into large ring-shaped superstructures attached to the plasma membrane. This supramolecular organization may further enhance the stability and connectivity of HDCR to form a specialized metabolic subcompartment within the cell.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Allosteric interactions prime androgen receptor dimerization and activation

The androgen receptor (AR) is a nuclear receptor that governs gene expression programs required for prostate development and male phenotype maintenance. Advanced prostate cancers display AR hyperactivation and transcriptome expansion, in part, through AR amplification and interaction with oncoprotein cofactors. Despite its biological importance, how AR domains and cofactors cooperate to bind DNA has remained elusive. Using single-particle cryo-electron microscopy, we isolated three conformations of AR bound to DNA, showing that AR forms a non-obligate dimer, with the buried dimer interface utilized by ancestral steroid receptors repurposed to facilitate cooperative DNA binding. We identify novel allosteric surfaces which are compromised in androgen insensitivity syndrome and reinforced by AR's oncoprotein cofactor, ERG, and by DNA-binding motifs. Finally, we present evidence that this plastic dimer interface may have been adopted for transactivation at the expense of DNA binding. Our work highlights how fine-tuning AR's cooperative interactions translate to consequences in development and disease.

Structural basis of SARM1 activation, substrate recognition, and inhibition by small molecules

The NADase SARM1 (sterile alpha and TIR motif containing 1) is a key executioner of axon degeneration and a therapeutic target for several neurodegenerative conditions. We show that a potent SARM1 inhibitor undergoes base exchange with the nicotinamide moiety of nicotinamide adenine dinucleotide (NAD+) to produce the bona fide inhibitor 1AD. We report structures of SARM1 in complex with 1AD, NAD+ mimetics and the allosteric activator nicotinamide mononucleotide (NMN). NMN binding triggers reorientation of the armadillo repeat (ARM) domains, which disrupts ARM:TIR interactions and leads to formation of a two-stranded TIR domain assembly. The active site spans two molecules in these assemblies, explaining the requirement of TIR domain self-association for NADase activity and axon degeneration. Our results reveal the mechanisms of SARM1 activation and substrate binding, providing rational avenues for the design of new therapeutics targeting SARM1.

Structural definition of a pan-sarbecovirus neutralizing epitope on the spike S2 subunit

Three betacoronaviruses have crossed the species barrier and established human-to-human transmission causing significant morbidity and mortality in the past 20 years. The most current and widespread of these is SARS-CoV-2. The identification of CoVs with zoonotic potential in animal reservoirs suggests that additional outbreaks could occur. Monoclonal antibodies targeting conserved neutralizing epitopes on diverse CoVs can form the basis for prophylaxis and therapeutic treatments and enable the design of vaccines aimed at providing pan-CoV protection. We previously identified a neutralizing monoclonal antibody, CV3-25 that binds to the SARS-CoV-2 spike, neutralizes the SARS-CoV-2 Beta variant comparably to the ancestral Wuhan Hu-1 strain, cross neutralizes SARS-CoV-1 and binds to recombinant proteins derived from the spike-ectodomains of HCoV-OC43 and HCoV-HKU1. Here, we show that the neutralizing activity of CV3-25 is maintained against the Alpha, Delta, Gamma and Omicron variants of concern as well as a SARS-CoV-like bat coronavirus with zoonotic potential by binding to a conserved linear peptide in the stem-helix region. Negative stain electron microscopy and a 1.74 Å crystal structure of a CV3-25/peptide complex demonstrates that CV3-25 binds to the base of the stem helix at the HR2 boundary to an epitope that is distinct from other stem-helix directed neutralizing mAbs.

Cryo-EM analysis of Ebola virus nucleocapsid-like assembly

This protocol describes the reconstitution of the filamentous Ebola virus nucleocapsid-like assembly in vitro. This is followed by solving the cryo-EM structure using helical reconstruction, and flexible fitting of the existing model into the 5.8 Å cryo-EM map. The protocol can be applied to other filamentous viral protein assemblies, particularly those with high flexibility and moderate resolution maps, which present technical challenges to model building.

For complete details on the use and execution of this profile, please refer to Su et al. (2018).

•

Preparation of Ebola nucleocapsid-like assembly for cryo-EM

•

Cryo-EM helical reconstruction of flexible filamentous protein assembly

•

Flexible fitting of protein model into cryo-EM density at moderate resolution

Application of Homology Modeling by Enhanced Profile-Profile Alignment and Flexible-Fitting Simulation to Cryo-EM Based Structure Determination

Application of cryo-electron microscopy (cryo-EM) is crucially important for ascertaining the atomic structure of large biomolecules such as ribosomes and protein complexes in membranes. Advances in cryo-EM technology and software have made it possible to obtain data with near-atomic resolution, but the method is still often capable of producing only a density map with up to medium resolution, either partially or entirely. Therefore, bridging the gap separating the density map and the atomic model is necessary. Herein, we propose a methodology for constructing atomic structure models based on cryo-EM maps with low-to-medium resolution. The method is a combination of sensitive and accurate homology modeling using our profile-profile alignment method with a flexible-fitting method using molecular dynamics simulation. As described herein, this study used benchmark applications to evaluate the model constructions of human two-pore channel 2 (one target protein in CASP13 with its structure determined using cryo-EM data) and the overall structure of Enterococcus hirae V-ATPase complex.

Cryo-ET of Toxoplasma parasites gives subnanometer insight into tubulin-based structures

Tubulin is a conserved protein that polymerizes into different forms of filamentous structures in Toxoplasma gondii, an obligate intracellular parasite in the phylum Apicomplexa. Two key tubulin-containing cytoskeletal components are subpellicular microtubules (SPMTs) and conoid fibrils (CFs). The SPMTs help maintain shape and gliding motility, while the CFs are implicated in invasion. Here, we use cryogenic electron tomography to determine the molecular structures of the SPMTs and CFs in vitrified intact and detergent-extracted parasites. Subvolume densities from detergent-extracted parasites yielded averaged density maps at subnanometer resolutions, and these were related back to their architecture in situ. An intralumenal spiral lines the interior of the 13-protofilament SPMTs, revealing a preferred orientation of these microtubules relative to the parasite's long axis. Each CF is composed of nine tubulin protofilaments that display a comma-shaped cross-section, plus additional associated components. Conoid protrusion, a crucial step in invasion, is associated with an altered pitch of each CF. The use of basic building blocks of protofilaments and different accessory proteins in one organism illustrates the versatility of tubulin to form two distinct types of assemblies, SPMTs and CFs.

Cryo-electron Tomography Reveals the Roles of FliY in Helicobacter pylori Flagellar Motor Assembly

Helicobacter pylori plays a causative role in gastric diseases. The pathogenicity of H. pylori depends on its ability to colonize the stomach guided by motility. FliY is a unique flagellar motor switch component coexisting with the classical FliG, FliM, and FliN switch proteins in some bacteria and has been shown to be essential for flagellation. However, the functional importance of FliY in H. pylori flagellar motor assembly is not well understood. Here, we applied cryo-electron tomography and subtomogram averaging to analyze the in situ structures of flagellar motors from wild-type strain, fliY-null mutant and complementation mutants expressing the N-terminal or C-terminal domain of FliY. Loss of full-length FliY or its C-terminal domain interrupted the formation of an intact C ring and soluble export apparatus, as well as the hook and flagellar filaments. Complementation with FliY C-terminal domain restored all these missing components of flagellar motor. Taken together, these results provide structural insights into the roles of FliY, especially its C-terminal domain in flagellar motor assembly in H. pylori.

IMPORTANCEHelicobacter pylori is the major risk factor related with gastric diseases. Flagellar motor is one of the most important virulence factors in H. pylori. However, the assembly mechanism of H. pylori flagellar motor is not fully understood yet. Previous report mainly described the overall structures of flagellum but had not focused on its specific components. Here, we focus on H. pylori flagellar C-ring protein FliY. We directly visualize the flagellar structures of H. pylori wild-type and FliY N-/C-terminal complementary strains by cryo-electron tomography and subtomogram averaging. Our results show that deletion of FliY or its C-terminal domain causes the loss of C ring, whereas deletion of FliY N-terminal does not affect C-ring assembly and flagellar structures. Our results provide direct evidence that C-ring protein FliY, especially its C-terminal domain, plays an indispensable role in H. pylori motor assembly and flagellar formation. This study will deepen our understanding about H. pylori pathogenesis.

Comprehensive structure and functional adaptations of the yeast nuclear pore complex

Nuclear pore complexes (NPCs) mediate the nucleocytoplasmic transport of macromolecules. Here we provide a structure of the isolated yeast NPC in which the inner ring is resolved by cryo-EM at sub-nanometer resolution to show how flexible connectors tie together different structural and functional layers. These connectors may be targets for phosphorylation and regulated disassembly in cells with an open mitosis. Moreover, some nucleoporin pairs and transport factors have similar interaction motifs, which suggests an evolutionary and mechanistic link between assembly and transport. We provide evidence for three major NPC variants that may foreshadow functional specializations at the nuclear periphery. Cryo-electron tomography extended these studies, providing a model of the in situ NPC with a radially expanded inner ring. Our comprehensive model reveals features of the nuclear basket and central transporter, suggests a role for the lumenal Pom152 ring in restricting dilation, and highlights structural plasticity that may be required for transport.

Structural basis for cytoplasmic dynein-1 regulation by Lis1

The lissencephaly 1 gene, LIS1, is mutated in patients with the neurodevelopmental disease lissencephaly. The Lis1 protein is conserved from fungi to mammals and is a key regulator of cytoplasmic dynein-1, the major minus-end-directed microtubule motor in many eukaryotes. Lis1 is the only dynein regulator known to bind directly to dynein's motor domain, and by doing so alters dynein's mechanochemistry. Lis1 is required for the formation of fully active dynein complexes, which also contain essential cofactors: dynactin and an activating adaptor. Here, we report the first high-resolution structure of the yeast dynein-Lis1 complex. Our 3.1 Å structure reveals, in molecular detail, the major contacts between dynein and Lis1 and between Lis1's ß-propellers. Structure-guided mutations in Lis1 and dynein show that these contacts are required for Lis1's ability to form fully active human dynein complexes and to regulate yeast dynein's mechanochemistry and in vivo function.

Evolution of Archaellum Rotation Involved Invention of a Stator Complex by Duplicating and Modifying a Core Component

Novelty in biology can arise from opportunistic repurposing of nascent characteristics of existing features. Understanding how this process happens at the molecular scale, however, suffers from a lack of case studies. The evolutionary emergence of rotary motors is a particularly clear example of evolution of a new function. The simplest of rotary motors is the archaellum, a molecular motor that spins a helical propeller for archaeal motility analogous to the bacterial flagellum. Curiously, emergence of archaellar rotation may have pivoted on the simple duplication and repurposing of a pre-existing component to produce a stator complex that anchors to the cell superstructure to enable productive rotation of the rotor component. This putative stator complex is composed of ArlF and ArlG, gene duplications of the filament component ArlB, providing an opportunity to study how gene duplication and neofunctionalization contributed to the radical innovation of rotary function. Toward understanding how this happened, we used electron cryomicroscopy to determine the structure of isolated ArlG filaments, the major component of the stator complex. Using a hybrid modeling approach incorporating structure prediction and validation, we show that ArlG filaments are open helices distinct to the closed helical filaments of ArlB. Curiously, further analysis reveals that ArlG retains a subset of the inter-protomer interactions of homologous ArlB, resulting in a superficially different assembly that nevertheless reflects the common ancestry of the two structures. This relatively simple mechanism to change quaternary structure was likely associated with the evolutionary neofunctionalization of the archaellar stator complex, and we speculate that the relative deformable elasticity of an open helix may facilitate elastic energy storage during the transmission of the discrete bursts of energy released by ATP hydrolysis to continuous archaellar rotation, allowing the inherent properties of a duplicated ArlB to be co-opted to fulfill a new role. Furthermore, agreement of diverse experimental evidence in our work supports recent claims to the power of new structure prediction techniques.

Bacteriophage Twort protein Gp168 is a β-clamp inhibitor by occupying the DNA sliding channel

Bacterial chromosome replication is mainly catalyzed by DNA polymerase III, whose beta subunits enable rapid processive DNA replication. Enabled by the clamp-loading complex, the two beta subunits form a ring-like clamp around DNA and keep the polymerase sliding along. Given the essential role of β-clamp, its inhibitors have been explored for antibacterial purposes. Similarly, β-clamp is an ideal target for bacteriophages to shut off host DNA synthesis during host takeover. The Gp168 protein of phage Twort is such an example, which binds to the β-clamp of Staphylococcus aureus and prevents it from loading onto DNA causing replication arrest. Here, we report a cryo-EM structure of the clamp–Gp168 complex at 3.2-Å resolution. In the structure of the complex, the Gp168 dimer occupies the DNA sliding channel of β-clamp and blocks its loading onto DNA, which represents a new inhibitory mechanism against β-clamp function. Interestingly, the key residues responsible for this interaction on the β-clamp are well conserved among bacteria. We therefore demonstrate that Gp168 is potentially a cross-species β-clamp inhibitor, as it forms complex with the Bacillus subtilis β-clamp. Our findings reveal an alternative mechanism for bacteriophages to inhibit β-clamp and provide a new strategy to combat bacterial drug resistance.

Application of transport-based metric for continuous interpolation between cryo-EM density maps

Cryogenic electron microscopy (cryo-EM) has become widely used for the past few years in structural biology, to collect single images of macromolecules "frozen in time". As this technique facilitates the identification of multiple conformational states adopted by the same molecule, a direct product of it is a set of 3D volumes, also called EM maps. To gain more insights on the possible mechanisms that govern transitions between different states, and hence the mode of action of a molecule, we recently introduced a bioinformatic tool that interpolates and generates morphing trajectories joining two given EM maps. This tool is based on recent advances made in optimal transport, that allow efficient evaluation of Wasserstein barycenters of 3D shapes. As the overall performance of the method depends on various key parameters, including the sensitivity of the regularization parameter, we performed various numerical experiments to demonstrate how MorphOT can be applied in different contexts and settings. Finally, we discuss current limitations and further potential connections between other optimal transport theories and the conformational heterogeneity problem inherent with cryo-EM data.

Bipartite binding and partial inhibition links DEPTOR and mTOR in a mutually antagonistic embrace

The mTORC1 kinase complex regulates cell growth, proliferation, and survival. Because mis-regulation of DEPTOR, an endogenous mTORC1 inhibitor, is associated with some cancers, we reconstituted mTORC1 with DEPTOR to understand its function. We find that DEPTOR is a unique partial mTORC1 inhibitor that may have evolved to preserve feedback inhibition of PI3K. Counterintuitively, mTORC1 activated by RHEB or oncogenic mutation is much more potently inhibited by DEPTOR. Although DEPTOR partially inhibits mTORC1, mTORC1 prevents this inhibition by phosphorylating DEPTOR, a mutual antagonism that requires no exogenous factors. Structural analyses of the mTORC1/DEPTOR complex showed DEPTOR's PDZ domain interacting with the mTOR FAT region, and the unstructured linker preceding the PDZ binding to the mTOR FRB domain. The linker and PDZ form the minimal inhibitory unit, but the N-terminal tandem DEP domains also significantly contribute to inhibition.

Validation, analysis and annotation of cryo-EM structures

The process of turning 2D micrographs into 3D atomic models of the imaged macromolecules has been under rapid development and scrutiny in the field of cryo-EM. Here, some important methods for validation at several stages in this process are described. Firstly, how Fourier shell correlation of two independent maps and phase randomization beyond a certain frequency address the assessment of map resolution is reviewed. Techniques for local resolution estimation and map sharpening are also touched upon. The topic of validating models which are either built de novo or based on a known atomic structure fitted into a cryo-EM map is then approached. Map-model comparison using Q-scores and Fourier shell correlation plots is used to assure the agreement of the model with the observed map density. The importance of annotating the model with B factors to account for the resolvability of individual atoms in the map is illustrated. Finally, the timely topic of detecting and validating water molecules and metal ions in maps that have surpassed ∼2 Å resolution is described.

Cryo-EM and antisense targeting of the 28-kDa frameshift stimulation element from the SARS-CoV-2 RNA genome

Drug discovery campaigns against COVID-19 are beginning to target the SARS-CoV-2 RNA genome. The highly conserved frameshift stimulation element (FSE), required for balanced expression of viral proteins, is a particularly attractive SARS-CoV-2 RNA target. Here we present a 6.9 Å resolution cryo-EM structure of the FSE (88 nucleotides, ~28 kDa), validated through an RNA nanostructure tagging method. The tertiary structure presents a topologically complex fold in which the 5' end is threaded through a ring formed inside a three-stem pseudoknot. Guided by this structure, we develop antisense oligonucleotides that impair FSE function in frameshifting assays and knock down SARS-CoV-2 virus replication in A549-ACE2 cells at 100 nM concentration.

Cryo-EM reveals conformational flexibility in apo DNA polymerase ζ

The translesion synthesis (TLS) DNA polymerases Rev1 and Polζ function together in DNA lesion bypass during DNA replication, acting as nucleotide inserter and extender polymerases, respectively. While the structural characterization of the Saccharomyces cerevisiae Polζ in its DNA-bound state has illuminated how this enzyme synthesizes DNA, a mechanistic understanding of TLS also requires probing conformational changes associated with DNA- and Rev1 binding. Here, we used single-particle cryo-electron microscopy to determine the structure of the apo Polζ holoenzyme. We show that compared with its DNA-bound state, apo Polζ displays enhanced flexibility that correlates with concerted motions associated with expansion of the Polζ DNA-binding channel upon DNA binding. We also identified a lysine residue that obstructs the DNA-binding channel in apo Polζ, suggesting a gating mechanism. The Polζ subunit Rev7 is a hub protein that directly binds Rev1 and is a component of several other protein complexes such as the shieldin DNA double-strand break repair complex. We analyzed the molecular interactions of budding yeast Rev7 in the context of Polζ and those of human Rev7 in the context of shieldin using a crystal structure of Rev7 bound to a fragment of the shieldin-3 protein. Overall, our study provides new insights into Polζ mechanism of action and the manner in which Rev7 recognizes partner proteins.

EvoSeg: Automated Electron Microscopy Segmentation through Random Forests and Evolutionary Optimization

Electron Microscopy Maps are key in the study of bio-molecular structures, ranging from borderline atomic level to the sub-cellular range. These maps describe the envelopes that cover possibly a very large number of proteins that form molecular machines within the cell. Within those envelopes, we are interested to find what regions correspond to specific proteins so that we can understand how they function, and design drugs that can enhance or suppress a process that they are involved in, along with other experimental purposes. A classic approach by which we can begin the exploration of map regions is to apply a segmentation algorithm. This yields a mask where each voxel in 3D space is assigned an identifier that maps it to a segment; an ideal segmentation would map each segment to one protein unit, which is rarely the case. In this work, we present a method that uses bio-inspired optimization, through an Evolutionary-Optimized Segmentation algorithm, to iteratively improve upon baseline segments obtained from a classical approach, called watershed segmentation. The cost function used by the evolutionary optimization is based on an ideal segmentation classifier trained as part of this development, which uses basic structural information available to scientists, such as the number of expected units, volume and topology. We show that a basic initial segmentation with the additional information allows our evolutionary method to find better segmentation results, compared to the baseline generated by the watershed.

Full-length de novo protein structure determination from cryo-EM maps using deep learning

MOTIVATION

Advances in microscopy instruments and image processing algorithms have led to an increasing number of cryo-EM maps. However, building accurate models for the EM maps at 3-5 Å resolution remains a challenging and time-consuming process. With the rapid growth of deposited EM maps, there is an increasing gap between the maps and reconstructed/modeled 3-dimensional (3D) structures. Therefore, automatic reconstruction of atomic-accuracy full-atomstructures fromEMmaps is pressingly needed.

RESULTS

We present a semi-automatic de novo structure determination method using a deep learningbased framework, named as DeepMM, which builds atomic-accuracy all-atom models from cryo-EM maps at near-atomic resolution. In our method, the main-chain and Cα positions as well as their amino acid and secondary structure types are predicted in the EM map using Densely Connected Convolutional Networks. DeepMM was extensively validated on 40 simulated maps at 5 Å resolution and 30 experimental maps at 2.6-4.8 Å resolution as well as an EMDB-wide data set of 2931 experimental maps at 2.6-4.9 Å resolution, and compared with state-of-the-art algorithms including RosettaES, MAINMAST, and Phenix. Overall, our DeepMM algorithm obtained a significant improvement over existing methods in terms of both accuracy and coverage in building full-length protein structures on all test sets, demonstrating the efficacy and general applicability of DeepMM.

AVAILABILITY

http://huanglab.phys.hust.edu.cn/DeepMM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predictable cholesterol binding sites in GPCRs lack consensus motifs

A rich diversity of transmembrane G protein-coupled receptors (GPCRs) are used by eukaryotes to sense physical and chemical signals. In humans alone, 800 GPCRs comprise the largest and most therapeutically targeted receptor class. Recent advances in GPCR structural biology have produced hundreds of GPCR structures solved by X-ray diffraction and increasingly, cryo-electron microscopy (cryo-EM). Many of these structures are stabilized by site-specific cholesterol binding, but it is unclear whether these interactions are a product of recurring cholesterol-binding motifs and if observed patterns of cholesterol binding differ by experimental technique. Here, we comprehensively analyze the location and composition of cholesterol binding sites in the current set of 473 human GPCR structural chains. Our findings establish that cholesterol binds similarly in cryo-EM and X-ray structures and show that 92% of cholesterol molecules on GPCR surfaces reside in predictable locations that lack discernable cholesterol-binding motifs.

Structural insights into integrin α5β1 opening by fibronectin ligand

Integrin α₅β₁ is a major fibronectin receptor critical for cell migration. Upon complex formation, fibronectin and α₅β₁ undergo conformational changes. While this is key for cell-tissue connections, its mechanism is unknown. Here, we report cryo-electron microscopy structures of native human α₅β₁ with fibronectin to 3.1-angstrom resolution, and in its resting state to 4.6-angstrom resolution. The α₅β₁-fibronectin complex revealed simultaneous interactions at the arginine-glycine-aspartate loop, the synergy site, and a newly identified binding site proximal to adjacent to metal ion-dependent adhesion site, inducing the translocation of helix α1 to secure integrin opening. Resting α₅β₁ adopts an incompletely bent conformation, challenging the model of integrin sharp bending inhibiting ligand binding. Our biochemical and structural analyses showed that affinity of α₅β₁ for fibronectin is increased with manganese ions (Mn²⁺) while adopting the half-bent conformation, indicating that ligand-binding affinity does not depend on conformation, and α₅β₁ opening is induced by ligand-binding.

Below 3 Å structure of apoferritin using a multipurpose TEM with a side entry cryoholder

Recently, the structural analysis of protein complexes by cryo-electron microscopy (cryo-EM) single particle analysis (SPA) has had great impact as a biophysical method. Many results of cryo-EM SPA are based on data acquired on state-of-the-art cryo-electron microscopes customized for SPA. These are currently only available in limited locations around the world, where securing machine time is highly competitive. One potential solution for this time-competitive situation is to reuse existing multi-purpose equipment, although this comes with performance limitations. Here, a multi-purpose TEM with a side entry cryo-holder was used to evaluate the potential of high-resolution SPA, resulting in a 3 Å resolution map of apoferritin with local resolution extending to 2.6 Å. This map clearly showed two positions of an aromatic side chain. Further, examination of optimal imaging conditions depending on two different multi-purpose electron microscope and camera combinations was carried out, demonstrating that higher magnifications are not always necessary or desirable. Since automation is effectively a requirement for large-scale data collection, and augmenting the multi-purpose equipment is possible, we expanded testing by acquiring data with SerialEM using a β-galactosidase test sample. This study demonstrates the possibilities of more widely available and established electron microscopes, and their applications for cryo-EM SPA.

SARM1 is a metabolic sensor activated by an increased NMN/NAD+ ratio to trigger axon degeneration

Axon degeneration is a central pathological feature of many neurodegenerative diseases. Sterile alpha and Toll/interleukin-1 receptor motif-containing 1 (SARM1) is a nicotinamide adenine dinucleotide (NAD+)-cleaving enzyme whose activation triggers axon destruction. Loss of the biosynthetic enzyme NMNAT2, which converts nicotinamide mononucleotide (NMN) to NAD+, activates SARM1 via an unknown mechanism. Using structural, biochemical, biophysical, and cellular assays, we demonstrate that SARM1 is activated by an increase in the ratio of NMN to NAD+ and show that both metabolites compete for binding to the auto-inhibitory N-terminal armadillo repeat (ARM) domain of SARM1. We report structures of the SARM1 ARM domain bound to NMN and of the homo-octameric SARM1 complex in the absence of ligands. We show that NMN influences the structure of SARM1 and demonstrate via mutagenesis that NMN binding is required for injury-induced SARM1 activation and axon destruction. Hence, SARM1 is a metabolic sensor responding to an increased NMN/NAD+ ratio by cleaving residual NAD+, thereby inducing feedforward metabolic catastrophe and axonal demise.

The myosin II coiled-coil domain atomic structure in its native environment

The atomic structure of the complete myosin tail within thick filaments isolated from Lethocerus indicus flight muscle is described and compared to crystal structures of recombinant, human cardiac myosin tail segments. Overall, the agreement is good with three exceptions: the proximal S2, in which the filament has heads attached but the crystal structure doesn't, and skip regions 2 and 4. At the head-tail junction, the tail α-helices are asymmetrically structured encompassing well-defined unfolding of 12 residues for one myosin tail, ∼4 residues of the other, and different degrees of α-helix unwinding for both tail α-helices, thereby providing an atomic resolution description of coiled-coil "uncoiling" at the head-tail junction. Asymmetry is observed in the nonhelical C termini; one C-terminal segment is intercalated between ribbons of myosin tails, the other apparently terminating at Skip 4 of another myosin tail. Between skip residues, crystal and filament structures agree well. Skips 1 and 3 also agree well and show the expected α-helix unwinding and coiled-coil untwisting in response to skip residue insertion. Skips 2 and 4 are different. Skip 2 is accommodated in an unusual manner through an increase in α-helix radius and corresponding reduction in rise/residue. Skip 4 remains helical in one chain, with the other chain unfolded, apparently influenced by the acidic myosin C terminus. The atomic model may shed some light on thick filament mechanosensing and is a step in understanding the complex roles that thick filaments of all species undergo during muscle contraction.

Fusion peptide priming reduces immune responses to HIV-1 envelope trimer base

Soluble "SOSIP"-stabilized envelope (Env) trimers are promising HIV-vaccine immunogens. However, they induce high-titer responses against the glycan-free trimer base, which is occluded on native virions. To delineate the effect on base responses of priming with immunogens targeting the fusion peptide (FP) site of vulnerability, here, we quantify the prevalence of trimer-base antibody responses in 49 non-human primates immunized with various SOSIP-stabilized Env trimers and FP-carrier conjugates. Trimer-base responses account for ∼90% of the overall trimer response in animals immunized with trimer only, ∼70% in animals immunized with a cocktail of SOSIP trimer and FP conjugate, and ∼30% in animals primed with FP conjugates before trimer immunization. Notably, neutralization breadth in FP-conjugate-primed animals correlates inversely with trimer-base responses. Our data provide methods to quantify the prevalence of trimer-base responses and reveal that FP-conjugate priming, either alone or as part of a cocktail, can reduce the trimer-base response and improve the neutralization outcome.

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

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

Structural interpretation of cryo-EM image reconstructions

The productivity of single-particle cryo-EM as a structure determination method has rapidly increased as many novel biological structures are being elucidated. The ultimate result of the cryo-EM experiment is an atomic model that should faithfully represent the computed image reconstruction. Although the principal approach of atomic model building and refinement from maps resembles that of the X-ray crystallographic methods, there are important differences due to the unique properties resulting from the 3D image reconstructions. In this review, we discuss the practiced work-flow from the cryo-EM image reconstruction to the atomic model. We give an overview of (i) resolution determination methods in cryo-EM including local and directional resolution variation, (ii) cryo-EM map contrast optimization including complementary map types that can help in identifying ambiguous density features, (iii) atomic model building and (iv) refinement in various resolution regimes including (v) their validation and (vi) discuss differences between X-ray and cryo-EM maps. Based on the methods originally developed for X-ray crystallography, the path from 3D image reconstruction to atomic coordinates has become an integral and important part of the cryo-EM structure determination work-flow that routinely delivers atomic models.

The C3/465 glycan hole cluster in BG505 HIV-1 envelope is the major neutralizing target involved in preventing mucosal SHIV infection

Stabilized HIV-1 envelope (Env) trimers elicit tier 2 autologous neutralizing antibody (nAb) responses in immunized animals. We previously demonstrated that BG505 SOSIP.664.T332N gp140 (BG505 SOSIP) immunization of rhesus macaques (RM) provided robust protection against autologous intra-vaginal simian-human immunodeficiency virus (SHIV) challenge that was predicted by high serum nAb titers. Here, we show that nAb in these protected RM targeted a glycan hole proximal to residue 465 in gp120 in all cases. nAb also targeted another glycan hole at residues 241/289 and an epitope in V1 at varying frequencies. Non-neutralizing antibodies directed at N611-shielded epitopes in gp41 were also present but were more prevalent in RM with low nAb titers. Longitudinal analysis demonstrated that nAb broadened in some RM during sequential immunization but remained focused in others, the latter being associated with increases in nAb titer. Thirty-eight monoclonal antibodies (mAbs) isolated from a protected RM with an exceptionally high serum neutralization titer bound to the trimer in ELISA, and four of the mAbs potently neutralized the BG505 Env pseudovirus (PV) and SHIV. The four neutralizing mAbs were clonally related and targeted the 465 glycan hole to varying degrees, mimicking the serum. The data demonstrate that the C3/465 glycan hole cluster was the dominant neutralization target in high titer protected RM, despite other co-circulating neutralizing and non-neutralizing specificities. The isolation of a neutralizing mAb family argues that clonotype expansion occurred during BG505 SOSIP immunization, leading to high titer, protective nAb and setting a desirable benchmark for HIV vaccines.

High-resolution mapping of the neutralizing and binding specificities of polyclonal sera post-HIV Env trimer vaccination

Mapping polyclonal serum responses is critical to rational vaccine design. However, most high-resolution mapping approaches involve isolating and characterizing individual antibodies, which incompletely defines the polyclonal response. Here we use two complementary approaches to directly map the specificities of the neutralizing and binding antibodies of polyclonal anti-HIV-1 sera from rabbits immunized with BG505 Env SOSIP trimers. We used mutational antigenic profiling to determine how all mutations in Env affected viral neutralization and electron microscopy polyclonal epitope mapping (EMPEM) to directly visualize serum Fabs bound to Env trimers. The dominant neutralizing specificities were generally only a subset of the more diverse binding specificities. Additional differences between binding and neutralization reflected antigenicity differences between virus and soluble Env trimer. Furthermore, we refined residue-level epitope specificity directly from sera, revealing subtle differences across sera. Together, mutational antigenic profiling and EMPEM yield a holistic view of the binding and neutralizing specificity of polyclonal sera.

Multiplexed Functional Assessment of Genetic Variants in CARD11

Genetic testing has increased the number of variants identified in disease genes, but the diagnostic utility is limited by lack of understanding variant function. CARD11 encodes an adaptor protein that expresses dominant-negative and gain-of-function variants associated with distinct immunodeficiencies. Here, we used a "cloning-free" saturation genome editing approach in a diploid cell line to simultaneously score 2,542 variants for decreased or increased function in the region of CARD11 associated with immunodeficiency. We also described an exon-skipping mechanism for CARD11 dominant-negative activity. The classification of reported clinical variants was sensitive (94.6%) and specific (88.9%), which rendered the data immediately useful for interpretation of seven coding and splicing variants implicated in immunodeficiency found in our clinic. This approach is generalizable for variant interpretation in many other clinically actionable genes, in any relevant cell type.

Structure of the shutdown state of myosin-2

Myosin-2 is essential for processes as diverse as cell division and muscle contraction. Dephosphorylation of its regulatory light chain (RLC) promotes an inactive, ‘shutdown’ state with the filament-forming tail folded onto the two heads¹, preventing filament formation and inactivating the motors². The mechanism by which this happens is obscure. Here we report a cryo-electron microscopy structure of shutdown smooth muscle myosin, with a resolution of 6 Å in the head region. A pseudo-atomic model, obtained by flexible fitting of crystal structures into the density and molecular dynamics simulations, describes interaction interfaces at the atomic level. The N-terminal extension of one RLC interacts with the tail and the other with the partner head, revealing how the RLCs stabilise the shutdown state in different ways and how their phosphorylation would allow myosin activation. Additional interactions between the three segments of the coiled coil, the motor domains and LCs stabilise the shutdown molecule. The structure of the lever in each head is competent to generate force upon activation. This shutdown structure is relevant to all myosin-2 isoforms and provides a framework for understanding their disease-causing mutations.

A 3.4-Å cryo-electron microscopy structure of the human coronavirus spike trimer computationally derived from vitrified NL63 virus particles

Human coronavirus NL63 (HCoV-NL63) is an enveloped pathogen of the family Coronaviridae that spreads worldwide and causes up to 10% of all annual respiratory diseases. HCoV-NL63 is typically associated with mild upper respiratory symptoms in children, elderly and immunocompromised individuals. It has also been shown to cause severe lower respiratory illness. NL63 shares ACE2 as a receptor for viral entry with SARS-CoV-1 and SARS-CoV-2. Here, we present the in situ structure of HCoV-NL63 spike (S) trimer at 3.4-Å resolution by single-particle cryo-EM imaging of vitrified virions without chemical fixative. It is structurally homologous to that obtained previously from the biochemically purified ectodomain of HCoV-NL63 S trimer, which displays a three-fold symmetric trimer in a single conformation. In addition to previously proposed and observed glycosylation sites, our map shows density at other sites, as well as different glycan structures. The domain arrangement within a protomer is strikingly different from that of the SARS-CoV-2 S and may explain their different requirements for activating binding to the receptor. This structure provides the basis for future studies of spike proteins with receptors, antibodies or drugs, in the native state of the coronavirus particles.

Molecular Architecture of the SARS-CoV-2 Virus

SARS-CoV-2 is an enveloped virus responsible for the COVID-19 pandemic. Despite recent advances in the structural elucidation of SARS-CoV-2 proteins, the detailed architecture of the intact virus remains to be unveiled. Here we report the molecular assembly of the authentic SARS-CoV-2 virus using cryoelectron tomography (cryo-ET) and subtomogram averaging (STA). Native structures of the S proteins in pre- and postfusion conformations were determined to average resolutions of 8.7-11 Å. Compositions of the N-linked glycans from the native spikes were analyzed by mass spectrometry, which revealed overall processing states of the native glycans highly similar to that of the recombinant glycoprotein glycans. The native conformation of the ribonucleoproteins (RNPs) and their higher-order assemblies were revealed. Overall, these characterizations revealed the architecture of the SARS-CoV-2 virus in exceptional detail and shed light on how the virus packs its ∼30-kb-long single-segmented RNA in the ∼80-nm-diameter lumen.

IM30 IDPs form a membrane-protective carpet upon super-complex disassembly

Members of the phage shock protein A (PspA) family, including the inner membrane-associated protein of 30 kDa (IM30), are suggested to stabilize stressed cellular membranes. Furthermore, IM30 is essential in thylakoid membrane-containing chloroplasts and cyanobacteria, where it is involved in membrane biogenesis and/or remodeling. While it is well known that PspA and IM30 bind to membranes, the mechanism of membrane stabilization is still enigmatic. Here we report that ring-shaped IM30 super-complexes disassemble on membranes, resulting in formation of a membrane-protecting protein carpet. Upon ring dissociation, the C-terminal domain of IM30 unfolds, and the protomers self-assemble on membranes. IM30 assemblies at membranes have been observed before in vivo and were associated with stress response in cyanobacteria and chloroplasts. These assemblies likely correspond to the here identified carpet structures. Our study defines the thus far enigmatic structural basis for the physiological function of IM30 and related proteins, including PspA, and highlights a hitherto unrecognized concept of membrane stabilization by intrinsically disordered proteins.

Full-length three-dimensional structure of the influenza A virus M1 protein and its organization into a matrix layer

Matrix proteins are encoded by many enveloped viruses, including influenza viruses, herpes viruses, and coronaviruses. Underneath the viral envelope of influenza virus, matrix protein 1 (M1) forms an oligomeric layer critical for particle stability and pH-dependent RNA genome release. However, high-resolution structures of full-length monomeric M1 and the matrix layer have not been available, impeding antiviral targeting and understanding of the pH-dependent transitions involved in cell entry. Here, purification and extensive mutagenesis revealed protein–protein interfaces required for the formation of multilayered helical M1 oligomers similar to those observed in virions exposed to the low pH of cell entry. However, single-layered helical oligomers with biochemical and ultrastructural similarity to those found in infectious virions before cell entry were observed upon mutation of a single amino acid. The highly ordered structure of the single-layered oligomers and their likeness to the matrix layer of intact virions prompted structural analysis by cryo-electron microscopy (cryo-EM). The resulting 3.4-Å–resolution structure revealed the molecular details of M1 folding and its organization within the single-shelled matrix. The solution of the full-length M1 structure, the identification of critical assembly interfaces, and the development of M1 assembly assays with purified proteins are crucial advances for antiviral targeting of influenza viruses.

Multi-subunit shells of matrix proteins line the interior of infectious influenza virus particles. In this study, biochemical purification of wild-type and mutant influenza M1 proteins allows the structural determination of an oligomer whose shape corresponds to that of infectious virions and suggests mechanisms for its formation and dismantling during infection.

A "Drug Sweeping" State of the TriABC Triclosan Efflux Pump from Pseudomonas aeruginosa

The structure of the TriABC inner membrane component of the triclosan/SDS-specific efflux pump from Pseudomonas aeruginosa was determined by cryoelectron microscopy to 4.5 Å resolution. The complete structure of the inner membrane transporter TriC of the resistance-nodulation-division (RND) superfamily was solved, including a partial structure of the fused periplasmic membrane fusion subunits, TriA and TriB. The substrate-free conformation of TriABC represents an intermediate step in efflux complex assembly before the engagement of the outer membrane channel. Structural analysis identified a tunnel network whose constriction impedes substrate efflux, indicating inhibition of TriABC in the unengaged state. Blind docking studies revealed binding to TriC at the same loci by substrates and bulkier non-substrates. Together with functional analyses, we propose that selective substrate translocation involves conformational gating at the tunnel narrowing that, together with conformational ordering of TriA and TriB, creates an engaged state capable of mediating substrate efflux.

The structure and symmetry of the radial spoke protein complex in Chlamydomonas flagella

The radial spoke is a key element in a transducer apparatus controlling the motility of eukaryotic cilia. The transduction biomechanics is a long-standing question in cilia biology. The radial spoke has three regions - a spoke head, a bifurcated neck and a stalk. Although the neck and the stalk are asymmetric, twofold symmetry of the head has remained controversial. In this work we used single particle cryo-electron microscopy (cryo-EM) analysis to generate a 3D structure of the whole radial spoke at unprecedented resolution. We show the head region at 15 Å (1.5 nm) resolution and confirm twofold symmetry. Using distance constraints generated by cross-linking mass spectrometry, we locate two components, RSP2 and RSP4, at the head and neck regions. Our biophysical analysis of isolated RSP4, RSP9, and RSP10 affirmed their oligomeric state. Our results enable us to redefine the boundaries of the regions and propose a model of organization of the radial spoke component proteins.

A 3.4-Å cryo-EM structure of the human coronavirus spike trimer computationally derived from vitrified NL63 virus particles

Human coronavirus NL63 (HCoV-NL63) is an enveloped pathogen of the family Coronaviridae that spreads worldwide and causes up to 10% of all annual respiratory diseases. HCoV-NL63 is typically associated with mild upper respiratory symptoms in children, elderly and immunocompromised individuals. It has also been shown to cause severe lower respiratory illness. NL63 shares ACE2 as a receptor for viral entry with SARS-CoV and SARS-CoV-2. Here we present the in situ structure of HCoV-NL63 spike (S) trimer at 3.4-Å resolution by single-particle cryo-EM imaging of vitrified virions without chemical fixative. It is structurally homologous to that obtained previously from the biochemically purified ectodomain of HCoV-NL63 S trimer, which displays a 3-fold symmetric trimer in a single conformation. In addition to previously proposed and observed glycosylation sites, our map shows density at other amino acid positions as well as differences in glycan structures. The domain arrangement within a protomer is strikingly different from that of the SARS-CoV-2 S and may explain their different requirements for activating binding to the receptor. This structure provides the basis for future studies of spike proteins with receptors, antibodies, or drugs, in the native state of the coronavirus particles.

Subdomain cryo-EM structure of nodaviral replication protein A crown complex provides mechanistic insights into RNA genome replication

For positive-strand RNA [(+)RNA] viruses, the major target for antiviral therapies is genomic RNA replication, which occurs at poorly understood membrane-bound viral RNA replication complexes. Recent cryoelectron microscopy (cryo-EM) of nodavirus RNA replication complexes revealed that the viral double-stranded RNA replication template is coiled inside a 30- to 90-nm invagination of the outer mitochondrial membrane, whose necked aperture to the cytoplasm is gated by a 12-fold symmetric, 35-nm diameter "crown" complex that contains multifunctional viral RNA replication protein A. Here we report optimizing cryo-EM tomography and image processing to improve crown resolution from 33 to 8.5 Å. This resolves the crown into 12 distinct vertical segments, each with 3 major subdomains: A membrane-connected basal lobe and an apical lobe that together comprise the ∼19-nm-diameter central turret, and a leg emerging from the basal lobe that connects to the membrane at ∼35-nm diameter. Despite widely varying replication vesicle diameters, the resulting two rings of membrane interaction sites constrain the vesicle neck to a highly uniform shape. Labeling protein A with a His-tag that binds 5-nm Ni-nanogold allowed cryo-EM tomography mapping of the C terminus of protein A to the apical lobe, which correlates well with the predicted structure of the C-proximal polymerase domain of protein A. These and other results indicate that the crown contains 12 copies of protein A arranged basally to apically in an N-to-C orientation. Moreover, the apical polymerase localization has significant mechanistic implications for template RNA recruitment and (-) and (+)RNA synthesis.

Mapping the immunogenic landscape of near-native HIV-1 envelope trimers in non-human primates

The induction of broad and potent immunity by vaccines is the key focus of research efforts aimed at protecting against HIV-1 infection. Soluble native-like HIV-1 envelope glycoproteins have shown promise as vaccine candidates as they can induce potent autologous neutralizing responses in rabbits and non-human primates. In this study, monoclonal antibodies were isolated and characterized from rhesus macaques immunized with the BG505 SOSIP.664 trimer to better understand vaccine-induced antibody responses. Our studies reveal a diverse landscape of antibodies recognizing immunodominant strain-specific epitopes and non-neutralizing neo-epitopes. Additionally, we isolated a subset of mAbs against an epitope cluster at the gp120-gp41 interface that recognize the highly conserved fusion peptide and the glycan at position 88 and have characteristics akin to several human-derived broadly neutralizing antibodies.

Modeling Epac1 interactions with the allosteric inhibitor AM-001 by co-solvent molecular dynamics

The exchange proteins activated by cAMP (EPAC) are implicated in a large variety of physiological processes and they are considered as promising targets for a wide range of therapeutic applications. Several recent reports provided evidence for the therapeutic effectiveness of the inhibiting EPAC1 activity cardiac diseases. In that context, we recently characterized a selective EPAC1 antagonist named AM-001. This compound was featured by a non-competitive mechanism of action but the localization of its allosteric site to EPAC1 structure has yet to be investigated. Therefore, we performed cosolvent molecular dynamics with the aim to identify a suitable allosteric binding site. Then, the docking and molecular dynamics were used to determine the binding of the AM-001 to the regions highlighted by cosolvent molecular dynamics for EPAC1. These analyses led us to the identification of a suitable allosteric AM-001 binding pocket at EPAC1. As a model validation, we also evaluated the binding poses of the available AM-001 analogues, with a different biological potency. Finally, the complex EPAC1 with AM-001 bound at the putative allosteric site was further refined by molecular dynamics. The principal component analysis led us to identify the protein motion that resulted in an inactive like conformation upon the allosteric inhibitor binding.

Cryo-electron Microscopy and Exploratory Antisense Targeting of the 28-kDa Frameshift Stimulation Element from the SARS-CoV-2 RNA Genome

Drug discovery campaigns against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) are beginning to target the viral RNA genome ^{1, 2} . The frameshift stimulation element (FSE) of the SARS-CoV-2 genome is required for balanced expression of essential viral proteins and is highly conserved, making it a potential candidate for antiviral targeting by small molecules and oligonucleotides ^3-6 . To aid global efforts focusing on SARS-CoV-2 frameshifting, we report exploratory results from frameshifting and cellular replication experiments with locked nucleic acid (LNA) antisense oligonucleotides (ASOs), which support the FSE as a therapeutic target but highlight difficulties in achieving strong inactivation. To understand current limitations, we applied cryogenic electron microscopy (cryo-EM) and the Ribosolve ⁷ pipeline to determine a three-dimensional structure of the SARS-CoV-2 FSE, validated through an RNA nanostructure tagging method. This is the smallest macromolecule (88 nt; 28 kDa) resolved by single-particle cryo-EM at subnanometer resolution to date. The tertiary structure model, defined to an estimated accuracy of 5.9 Å, presents a topologically complex fold in which the 5' end threads through a ring formed inside a three-stem pseudoknot. Our results suggest an updated model for SARS-CoV-2 frameshifting as well as binding sites that may be targeted by next generation ASOs and small molecules.

Electromechanical coupling in the hyperpolarization-activated K+ channel KAT1

Voltage-gated potassium (K_v) channels orchestrate electrical signaling and control cell volume by gating in response to either membrane depolarization or hyperpolarization. Yet, while all voltage-sensing domains transduce transmembrane electric field changes by a common mechanism involving the outward or inward translocation of gating charges^1–3, the general determinants of channel gating polarity remain poorly understood⁴. Here, we suggest a molecular mechanism for electromechanical coupling and gating polarity in non-domain-swapped K_v channels based on the cryo-EM structure of KAT1, the hyperpolarization-activated K_v channel from Arabidopsis thaliana. KAT1 displays a depolarized voltage sensor, which interacts with a closed pore domain directly via two interfaces and indirectly via an intercalated phospholipid. Functional evaluation of KAT1 structure-guided mutants at the sensor-pore interfaces suggests a mechanism in which direct interaction between the sensor and C-linker hairpin in the adjacent pore subunit is the primary determinant of gating polarity. We suggest that a ~5–7 Å inward motion of the S4 sensor helix can underlie a direct-coupling mechanism, driving a conformational reorientation of the C-linker and ultimately opening the activation gate formed by the S6 intracellular bundle. This direct-coupling mechanism contrasts with allosteric mechanisms proposed for hyperpolarization-activated cyclic nucleotide-gated (HCN) channels⁵, and may represent an unexpected link between depolarization and hyperpolarization-activated channels.

Practical Considerations for Atomistic Structure Modeling with Cryo-EM Maps

We describe common approaches to atomistic structure modeling with single particle analysis derived cryo-EM maps. Several strategies for atomistic model building and atomistic model fitting methods are discussed, including selection criteria and implementation procedures. In covering basic concepts and caveats, this short perspective aims to help facilitate active discussion between scientists at different levels with diverse backgrounds.

Structure of Microtubule-Trapped Human Kinesin-5 and Its Mechanism of Inhibition Revealed Using Cryoelectron Microscopy

Kinesin-5 motors are vital mitotic spindle components, and disruption of their function perturbs cell division. We investigated the molecular mechanism of the human kinesin-5 inhibitor GSK-1, which allosterically promotes tight microtubule binding. GSK-1 inhibits monomeric human kinesin-5 ATPase and microtubule gliding activities, and promotes the motor's microtubule stabilization activity. Using cryoelectron microscopy, we determined the 3D structure of the microtubule-bound motor-GSK-1 at 3.8 Å overall resolution. The structure reveals that GSK-1 stabilizes the microtubule binding surface of the motor in an ATP-like conformation, while destabilizing regions of the motor around the empty nucleotide binding pocket. Density corresponding to GSK-1 is located between helix-α4 and helix-α6 in the motor domain at its interface with the microtubule. Using a combination of difference mapping and protein-ligand docking, we characterized the kinesin-5-GSK-1 interaction and further validated this binding site using mutagenesis. This work opens up new avenues of investigation of kinesin inhibition and spindle perturbation.

2.7 Å cryo-EM structure of rotavirus core protein VP3, a unique capping machine with a helicase activity

In many viruses, including rotavirus (RV), the major pathogen of infantile gastroenteritis, capping of viral messenger RNAs is a pivotal step for efficient translation of the viral genome. In RV, VP3 caps the nascent transcripts synthesized from the genomic dsRNA segments by the RV polymerase VP1 within the particle core. Here, from cryo-electron microscopy, x-ray crystallography, and biochemical analyses, we show that VP3 forms a stable tetrameric assembly with each subunit having a modular domain organization, which uniquely integrates five distinct enzymatic steps required for capping the transcripts. In addition to the previously known guanylyl- and methyltransferase activities, we show that VP3 exhibits hitherto unsuspected RNA triphosphatase activity necessary for initiating transcript capping and RNA helicase activity likely required for separating the RNA duplex formed transiently during endogenous transcription. From our studies, we propose a new mechanism for how VP3 inside the virion core caps the nascent transcripts exiting from the polymerase.