Measurement of atom resolvability in cryo-EM maps with Q-scores

Elucidating structural variability in p53 conformers using combinatorial refinement strategies and molecular dynamics

Low molecular weight proteins and protein assemblies can now be investigated using cryo-electron microscopy (EM) as a complement to traditional structural biology techniques. It is important, however, to not lose sight of the dynamic information inherent in macromolecules that give rise to their exquisite functionality. As computational methods continue to advance the field of biomedical imaging, so must strategies to resolve the minute details of disease-related entities. Here, we employed combinatorial modeling approaches to assess flexible properties among low molecular weight proteins (~100 kDa or less). Through a blend of rigid body refinement and simulated annealing, we determined new hidden conformations for wild type p53 monomer and dimer forms. Structures for both states converged to yield new conformers, each revealing good stereochemistry and dynamic information about the protein. Based on these insights, we identified fluid parts of p53 that complement the stable central core of the protein responsible for engaging DNA. Molecular dynamics simulations corroborated the modeling results and helped pinpoint the more flexible residues in wild type p53. Overall, the new computational methods may be used to shed light on other small protein features in a vast ensemble of structural data that cannot be easily delineated by other algorithms.

Elucidating the Architectural dynamics of MuB filaments in bacteriophage Mu DNA transposition

MuB is a non-specific DNA-binding protein and AAA+ ATPase that significantly influences the DNA transposition process of bacteriophage Mu, especially in target DNA selection for transposition. While studies have established the ATP-dependent formation of MuB filament as pivotal to this process, the high-resolution structure of a full-length MuB protomer and the underlying molecular mechanisms governing its oligomerization remain elusive. Here, we use cryo-EM to obtain a 3.4-Å resolution structure of the ATP(+)-DNA(+)-MuB helical filament, which encapsulates the DNA substrate within its axial channel. The structure categorizes MuB within the initiator clade of the AAA+ protein family and precisely locates the ATP and DNA binding sites. Further investigation into the oligomeric states of MuB show the existence of various forms of the filament. These findings lead to a mechanistic model where MuB forms opposite helical filaments along the DNA, exposing potential target sites on the bare DNA and then recruiting MuA, which stimulates MuB's ATPase activity and disrupts the previously formed helical structure. When this happens, MuB generates larger ring structures and dissociates from the DNA.

High-resolution yeast actin structures indicate the molecular mechanism of actin filament stiffening by cations

Actin filament assembly and the regulation of its mechanical properties are fundamental processes essential for eukaryotic cell function. Residue E167 in vertebrate actins forms an inter-subunit salt bridge with residue K61 of the adjacent subunit. Saccharomyces cerevisiae actin filaments are more flexible than vertebrate filaments and have an alanine at this position (A167). Substitution of this alanine for a glutamic acid (A167E) confers Saccharomyces cerevisiae actin filaments with salt-dependent stiffness similar to vertebrate actins. We developed an optimized cryogenic electron microscopy workflow refining sample preparation and vitrification to obtain near-atomic resolution structures of wild-type and A167E mutant Saccharomyces cerevisiae actin filaments. The difference between these structures allowed us to pinpoint the potential binding site of a filament-associated cation that controls the stiffness of the filaments in vertebrate and A167E Saccharomyces cerevisiae actins. Through an analysis of previously published high-resolution reconstructions of vertebrate actin filaments, along with a newly determined high-resolution vertebrate actin structure in the absence of potassium, we identified a unique peak near residue 167 consistent with the binding of a magnesium ion. Our findings show how magnesium can contribute to filament stiffening by directly bridging actin subunits and allosterically affecting the orientation of the DNase-I binding loop of actin, which plays a regulatory role in modulating actin filament stiffness and interactions with regulatory proteins.

Structural and mechanistic insights into Streptococcus pneumoniae NADPH oxidase

Nicotinamide adenine dinucleotide phosphate (NADPH) oxidases (NOXs) have a major role in the physiology of eukaryotic cells by mediating reactive oxygen species production. Evolutionarily distant proteins with the NOX catalytic core have been found in bacteria, including Streptococcus pneumoniae NOX (SpNOX), which is proposed as a model for studying NOXs because of its high activity and stability in detergent micelles. We present here cryo-electron microscopy structures of substrate-free and nicotinamide adenine dinucleotide (NADH)-bound SpNOX and of NADPH-bound wild-type and F397A SpNOX under turnover conditions. These high-resolution structures provide insights into the electron-transfer pathway and reveal a hydride-transfer mechanism regulated by the displacement of F397. We conducted structure-guided mutagenesis and biochemical analyses that explain the absence of substrate specificity toward NADPH and suggest the mechanism behind constitutive activity. Our study presents the structural basis underlying SpNOX enzymatic activity and sheds light on its potential in vivo function.

Molecular basis of FIGNL1 in dissociating RAD51 from DNA and chromatin

Maintaining genome integrity is an essential and challenging process. RAD51 recombinase, the central player of several crucial processes in repairing and protecting genome integrity, forms filaments on DNA. RAD51 filaments are tightly regulated. One of these regulators is FIGNL1, that prevents persistent RAD51 foci post-damage and genotoxic chromatin association in cells. The cryogenic electron microscopy structure of FIGNL1 in complex with RAD51 reveals that the FIGNL1 forms a non-planar hexamer and RAD51 N-terminus is enclosed in the FIGNL1 hexamer pore. Mutations in pore loop or catalytic residues of FIGNL1 render it defective in filament disassembly and are lethal in mouse embryonic stem cells. Our study reveals a unique mechanism for removing RAD51 from DNA and provides the molecular basis for FIGNL1 in maintaining genome stability.

Assembly of respiratory syncytial virus matrix protein lattice and its coordination with fusion glycoprotein trimers

Respiratory syncytial virus (RSV) is an enveloped, filamentous, negative-strand RNA virus that causes significant respiratory illness worldwide. RSV vaccines are available, however there is still significant need for research to support the development of vaccines and therapeutics against RSV and related Mononegavirales viruses. Individual virions vary in size, with an average diameter of ~130 nm and ranging from ~500 nm to over 10 µm in length. Though the general arrangement of structural proteins in virions is known, we use cryo-electron tomography and sub-tomogram averaging to determine the molecular organization of RSV structural proteins. We show that the peripheral membrane-associated RSV matrix (M) protein is arranged in a packed helical-like lattice of M-dimers. We report that RSV F glycoprotein is frequently observed as pairs of trimers oriented in an anti-parallel conformation to support potential interactions between trimers. Our sub-tomogram averages indicate the positioning of F-trimer pairs is correlated with the underlying M lattice. These results provide insight into RSV virion organization and may aid in the development of RSV vaccines and anti-viral targets.

Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Factors affecting macromolecule orientations in thin films formed in cryo-EM

The formation of a vitrified thin film embedded with randomly oriented macromolecules is an essential prerequisite for cryogenic sample electron microscopy. Most commonly, this is achieved using the plunge-freeze method first described nearly 40 years ago. Although this is a robust method, the behaviour of different macromolecules shows great variation upon freezing and often needs to be optimized to obtain an isotropic, high-resolution reconstruction. For a macromolecule in such a film, the probability of encountering the air-water interface in the time between blotting and freezing and adopting preferred orientations is very high. 3D reconstruction using preferentially oriented particles often leads to anisotropic and uninterpretable maps. Currently, there are no general solutions to this prevalent issue, but several approaches largely focusing on sample preparation with the use of additives and novel grid modifications have been attempted. In this study, the effect of physical and chemical factors on the orientations of macromolecules was investigated through an analysis of selected well studied macromolecules, and important parameters that determine the behaviour of proteins on cryo-EM grids were revealed. These insights highlight the nature of the interactions that cause preferred orientations and can be utilized to systematically address orientation bias for any given macromolecule and to provide a framework to design small-molecule additives to enhance sample stability and behaviour.

Accurate model and ensemble refinement using cryo-electron microscopy maps and Bayesian inference

Converting cryo-electron microscopy (cryo-EM) data into high-quality structural models is a challenging problem of outstanding importance. Current refinement methods often generate unbalanced models in which physico-chemical quality is sacrificed for excellent fit to the data. Furthermore, these techniques struggle to represent the conformational heterogeneity averaged out in low-resolution regions of density maps. Here we introduce EMMIVox, a Bayesian inference approach to determine single-structure models as well as structural ensembles from cryo-EM maps. EMMIVox automatically balances experimental information with accurate physico-chemical models of the system and the surrounding environment, including waters, lipids, and ions. Explicit treatment of data correlation and noise as well as inference of accurate B-factors enable determination of structural models and ensembles with both excellent fit to the data and high stereochemical quality, thus outperforming state-of-the-art refinement techniques. EMMIVox represents a flexible approach to determine high-quality structural models that will contribute to advancing our understanding of the molecular mechanisms underlying biological functions.

Structure and function of the SIT1 proline transporter in complex with the COVID-19 receptor ACE2

Proline is widely known as the only proteogenic amino acid with a secondary amine. In addition to its crucial role in protein structure, the secondary amino acid modulates neurotransmission and regulates the kinetics of signaling proteins. To understand the structural basis of proline import, we solved the structure of the proline transporter SIT1 in complex with the COVID-19 viral receptor ACE2 by cryo-electron microscopy. The structure of pipecolate-bound SIT1 reveals the specific sequence requirements for proline transport in the SLC6 family and how this protein excludes amino acids with extended side chains. By comparing apo and substrate-bound SIT1 states, we also identify the structural changes that link substrate release and opening of the cytoplasmic gate and provide an explanation for how a missense mutation in the transporter causes iminoglycinuria.

Influence of the interaction between p53 and ZNF568 on mitochondrial oxidative phosphorylation

The tumor suppressor p53 plays important roles in suppressing the development and progression of cancer by responding to various stress signals. In addition, p53 can regulate the metabolic pathways of cancer cells by regulating energy metabolism and oxidative phosphorylation. Here, we present a mechanism for the interaction between p53 and ZNF568. Initially, we used X-ray crystallography to determine the irregular loop structure of the ZNF568 KRAB domain; this loop plays an important role in the interaction between p53 and ZNF568. In addition, Cryo-EM was used to examine how the p53 DBD and ZNF568 KRAB domains bind together. The function of ZNF568 on p53-mediated mitochondrial respiration was confirmed by measuring glucose consumption and lactate production. These findings show that ZNF568 can reduce p53-mediated mitochondrial respiratory activity by binding to p53 and inhibiting the transcription of SCO2. SIGNIFICANCE: ZNF568 can directly bind to the p53 DBD and transcriptionally regulate the SCO2 gene. SCO2 transcriptional regulation by interaction between ZNF568 and p53 may regulate the balance between mitochondrial respiration and glycolysis.

High-resolution single-particle imaging at 100-200 keV with the Gatan Alpine direct electron detector

Developments in direct electron detector technology have played a pivotal role in enabling high-resolution structural studies by cryo-EM at 200 and 300 keV. Yet, theory and recent experiments indicate advantages to imaging at 100 keV, energies for which the current detectors have not been optimized. In this study, we evaluated the Gatan Alpine detector, designed for operation at 100 and 200 keV. Compared to the Gatan K3, Alpine demonstrated a significant DQE improvement at these energies, specifically a ∼ 4-fold improvement at Nyquist at 100 keV. In single-particle cryo-EM experiments, Alpine datasets yielded better than 2 Å resolution reconstructions of apoferritin at 120 and 200 keV on a ThermoFisher Scientific (TFS) Glacios microscope fitted with a non-standard SP-Twin lens. We also achieved a ∼ 3.2 Å resolution reconstruction of a 115 kDa asymmetric protein complex, proving Alpine's effectiveness with complex biological samples. In-depth analysis revealed that Alpine reconstructions are comparable to K3 reconstructions at 200 keV, and remarkably, reconstruction from Alpine at 120 keV on a TFS Glacios surpassed all but the 300 keV data from a TFS Titan Krios with GIF/K3. Additionally, we show Alpine's capability for high-resolution data acquisition and screening on lower-end systems by obtaining ∼ 3 Å resolution reconstructions of apoferritin and aldolase at 100 keV and detailed 2D averages of a 55 kDa sample using a side-entry cryo holder. Overall, we show that Gatan Alpine performs well with the standard 200 keV imaging systems and may potentially capture the benefits of lower accelerating voltages, bringing smaller sized particles within the scope of cryo-EM.

Hub stability in the calcium calmodulin-dependent protein kinase II

The calcium calmodulin protein kinase II (CaMKII) is a multi-subunit ring assembly with a central hub formed by the association domains. There is evidence for hub polymorphism between and within CaMKII isoforms, but the link between polymorphism and subunit exchange has not been resolved. Here, we present near-atomic resolution cryogenic electron microscopy (cryo-EM) structures revealing that hubs from the α and β isoforms, either standalone or within an β holoenzyme, coexist as 12 and 14 subunit assemblies. Single-molecule fluorescence microscopy of Venus-tagged holoenzymes detects intermediate assemblies and progressive dimer loss due to intrinsic holoenzyme lability, and holoenzyme disassembly into dimers upon mutagenesis of a conserved inter-domain contact. Molecular dynamics (MD) simulations show the flexibility of 4-subunit precursors, extracted in-silico from the β hub polymorphs, encompassing the curvature of both polymorphs. The MD explains how an open hub structure also obtained from the β holoenzyme sample could be created by dimer loss and analysis of its cryo-EM dataset reveals how the gap could open further. An assembly model, considering dimer concentration dependence and strain differences between polymorphs, proposes a mechanism for intrinsic hub lability to fine-tune the stoichiometry of αβ heterooligomers for their dynamic localization within synapses in neurons.

Automated multiconformer model building for X-ray crystallography and cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift toward modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior Rfree and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g., Coot) and fit can be further improved by refinement using standard pipelines (e.g., Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Structural Requirements for Reverse Transcription by a Diversity-generating Retroelement

Diversity-generating retroelements (DGRs) create massive protein sequence variation in ecologically diverse microbes. Variation occurs during reverse transcription of a protein-encoding RNA template coupled to misincorporation at adenosines. In the prototypical Bordetella bacteriophage DGR, the template must be surrounded by upstream and downstream RNA segments for cDNA synthesis by the reverse transcriptase bRT and associated protein Avd. The function of the surrounding RNA was unknown. Cryo-EM revealed that this RNA enveloped bRT and lay over barrel-shaped Avd, forming an intimate ribonucleoprotein (RNP). An abundance of essential interactions between RNA structural elements and bRT-Avd precisely positioned an RNA homoduplex for initiation of cDNA synthesis by cis -priming. Our results explain how the surrounding RNA primes cDNA synthesis, promotes processivity, terminates polymerization, and strictly limits mutagenesis to select proteins through mechanisms that are likely conserved in DGRs from distant taxa.

Principles of peptide selection by the transporter associated with antigen processing

Our ability to fight pathogens relies on major histocompatibility complex class I (MHC-I) molecules presenting diverse antigens on the surface of diseased cells. The transporter associated with antigen processing (TAP) transports nearly the entire repertoire of antigenic peptides into the endoplasmic reticulum for MHC-I loading. How TAP transports peptides specific for MHC-I is unclear. In this study, we used cryo-EM to determine a series of structures of human TAP, both in the absence and presence of peptides with various sequences and lengths. The structures revealed that peptides of eight or nine residues in length bind in a similarly extended conformation, despite having little sequence overlap. We also identified two peptide-anchoring pockets on either side of the transmembrane cavity, each engaging one end of a peptide with primarily main chain atoms. Occupation of both pockets results in a global conformational change in TAP, bringing the two halves of the transporter closer together to prime it for isomerization and ATP hydrolysis. Shorter peptides are able to bind to each pocket separately but are not long enough to bridge the cavity to bind to both simultaneously. Mutations that disrupt hydrogen bonds with the N and C termini of peptides almost abolish MHC-I surface expression. Our findings reveal that TAP functions as a molecular caliper that selects peptides according to length rather than sequence, providing antigen diversity for MHC-I presentation.

An interaction network in the polymerase active site is a prerequisite for Watson-Crick base pairing in Pol γ

The replication accuracy of DNA polymerase gamma (Pol γ) is essential for mitochondrial genome integrity. Mutation of human Pol γ arginine-853 has been linked to neurological diseases. Although not a catalytic residue, Pol γ arginine-853 mutants are void of polymerase activity. To identify the structural basis for the disease, we determined a crystal structure of the Pol γ mutant ternary complex with correct incoming nucleotide 2'-deoxycytidine 5'-triphosphate (dCTP). Opposite to the wild type that undergoes open-to-closed conformational changes when bound to a correct nucleotide that is essential for forming a catalytically competent active site, the mutant complex failed to undergo the conformational change, and the dCTP did not base pair with its Watson-Crick complementary templating residue. Our studies revealed that arginine-853 coordinates an interaction network that aligns the 3'-end of primer and dCTP with the catalytic residues. Disruption of the network precludes the formation of Watson-Crick base pairing and closing of the active site, resulting in an inactive polymerase.

CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase

Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.

Platform-directed allostery and quaternary structure dynamics of SAMHD1 catalysis

SAMHD1 regulates cellular nucleotide homeostasis, controlling dNTP levels by catalysing their hydrolysis into 2'-deoxynucleosides and triphosphate. In differentiated CD4+ macrophage and resting T-cells SAMHD1 activity results in the inhibition of HIV-1 infection through a dNTP blockade. In cancer, SAMHD1 desensitizes cells to nucleoside-analogue chemotherapies. Here we employ time-resolved cryogenic-EM imaging and single-particle analysis to visualise assembly, allostery and catalysis by this multi-subunit enzyme. Our observations reveal how dynamic conformational changes in the SAMHD1 quaternary structure drive the catalytic cycle. We capture five states at high-resolution in a live catalytic reaction, revealing how allosteric activators support assembly of a stable SAMHD1 tetrameric core and how catalysis is driven by the opening and closing of active sites through pairwise coupling of active sites and order-disorder transitions in regulatory domains. This direct visualisation of enzyme catalysis dynamics within an allostery-stabilised platform sets a precedent for mechanistic studies into the regulation of multi-subunit enzymes.

Uncovering Protein Ensembles: Automated Multiconformer Model Building for X-ray Crystallography and Cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift towards modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior R f r e e and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g. Coot) and fit can be further improved by refinement using standard pipelines (e.g. Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Substrate binding plasticity revealed by Cryo-EM structures of SLC26A2

SLC26A2 is a vital solute carrier responsible for transporting essential nutritional ions, including sulfate, within the human body. Pathogenic mutations within SLC26A2 give rise to a spectrum of human diseases, ranging from lethal to mild symptoms. The molecular details regarding the versatile substrate-transporter interactions and the impact of pathogenic mutations on SLC26A2 transporter function remain unclear. Here, using cryo-electron microscopy, we determine three high-resolution structures of SLC26A2 in complexes with different substrates. These structures unveil valuable insights, including the distinct features of the homodimer assembly, the dynamic nature of substrate binding, and the potential ramifications of pathogenic mutations. This structural-functional information regarding SLC26A2 will advance our understanding of cellular sulfate transport mechanisms and provide foundations for future therapeutic development against various human diseases.

Human polymerase theta helicase positions DNA microhomologies for double-strand break repair

DNA double-strand breaks occur in all human cells on a daily basis and must be repaired with high fidelity to minimize genomic instability1. Deficiencies in high-fidelity DNA repair by homologous recombination lead to dependence on DNA polymerase theta, which identifies DNA microhomologies in 3' single-stranded DNA overhangs and anneals them to initiate error-prone double-strand break repair. The resulting genomic instability is associated with numerous cancers, thereby making this polymerase an attractive therapeutic target2,3. However, despite the biomedical importance of polymerase theta, the molecular details of how it initiates DNA break repair remain unclear4,5. Here we present cryo-electron microscopy structures of the polymerase theta helicase domain bound to microhomology-containing DNA, revealing DNA-induced rearrangements of the helicase that enable DNA repair. Our structures show that DNA-bound helicase dimers facilitate a microhomology search that positions 3' single-stranded DNA ends in proximity to align complementary base pairs and anneal DNA microhomology. We define the molecular determinants that enable the polymerase theta helicase domain to identify and pair DNA microhomologies to initiate mutagenic DNA repair, providing mechanistic insights into therapeutic targeting of these interactions.

Machine learning approaches to cryoEM density modification differentially affect biomacromolecule and ligand density quality

The application of machine learning to cryogenic electron microscopy (cryoEM) data analysis has added a valuable set of tools to the cryoEM data processing pipeline. As these tools become more accessible and widely available, the implications of their use should be assessed. We noticed that machine learning map modification tools can have differential effects on cryoEM densities. In this perspective, we evaluate these effects to show that machine learning tools generally improve densities for biomacromolecules while generating unpredictable results for ligands. This unpredictable behavior manifests both in quantitative metrics of map quality and in qualitative investigations of modified maps. The results presented here highlight the power and potential of machine learning tools in cryoEM, while also illustrating some of the risks of their unexamined use.

Dynamic inter-domain transformations mediate the allosteric regulation of human 5, 10-methylenetetrahydrofolate reductase

5,10-methylenetetrahydrofolate reductase (MTHFR) commits folate-derived one-carbon units to generate the methyl-donor S-adenosyl-L-methionine (SAM). Eukaryotic MTHFR appends to the well-conserved catalytic domain (CD) a unique regulatory domain (RD) that confers feedback inhibition by SAM. Here we determine the cryo-electron microscopy structures of human MTHFR bound to SAM and its demethylated product S-adenosyl-L-homocysteine (SAH). In the active state, with the RD bound to a single SAH, the CD is flexible and exposes its active site for catalysis. However, in the inhibited state the RD pocket is remodelled, exposing a second SAM-binding site that was previously occluded. Dual-SAM bound MTHFR demonstrates a substantially rearranged inter-domain linker that reorients the CD, inserts a loop into the active site, positions Tyr404 to bind the cofactor FAD, and blocks substrate access. Our data therefore explain the long-distance regulatory mechanism of MTHFR inhibition, underpinned by the transition between dual-SAM and single-SAH binding in response to cellular methylation status.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

Removal of Pseudomonas type IV pili by a small RNA virus

The retractile type IV pilus (T4P) is important for virulence of the opportunistic human pathogen Pseudomonas aeruginosa. The single-stranded RNA (ssRNA) phage PP7 binds to T4P and is brought to the cell surface through pilus retraction. Using fluorescence microscopy, we discovered that PP7 detaches T4P, which impairs cell motility and restricts the pathogen's virulence. Using cryo-electron microscopy, mutagenesis, optical trapping, and Langevin dynamics simulation, we resolved the structure of PP7, T4P, and the PP7/T4P complex and showed that T4P detachment is driven by the affinity between the phage maturation protein and its bound pilin, plus the pilus retraction force and speed, and pilus bending. Pilus detachment may be widespread among other ssRNA phages and their retractile pilus systems and offers new prospects for antibacterial prophylaxis and therapeutics.

Purification, identification and Cryo-EM structure of prostatic acid phosphatase in human semen

Prostatic acid phosphatase (PAP) is a glycoprotein that plays a crucial role in the hydrolysis of phosphate ester present in prostatic exudates. It is a well-established indicator for prostate cancer due to its elevated serum levels in disease progression. Despite its abundance in semen, PAP's influence on male fertility has not been extensively studied. In our study, we report a significantly optimized method for purifying human endogenous PAP, achieving remarkably high efficiency and active protein recovery rate. This achievement allowed us to better analyze and understand the PAP protein. We determined the cryo-electron microscopic (Cryo-EM) structure of prostatic acid phosphatase in its physiological state for the first time. Our structural and gel filtration analysis confirmed the formation of a tight homodimer structure of human PAP. This functional homodimer displayed an elongated conformation in the cryo-EM structure compared to the previously reported crystal structure. Additionally, there was a notable 5-degree rotation in the angle between the α domain and α/β domain of each monomer. Through structural analysis, we revealed three potential glycosylation sites: Asn94, Asn220, and Asn333. These sites contained varying numbers and forms of glycosyl units, suggesting sugar moieties influence PAP function. Furthermore, we found that the active sites of PAP, His44 and Asp290, are located between the two protein domains. Overall, our study not only provide an optimized approach for PAP purification, but also offer crucial insights into its structural characteristics. These findings lay the groundwork for further investigations into the physiological function and potential therapeutic applications of this important protein.

Automated model building and protein identification in cryo-EM maps

Interpreting electron cryo-microscopy (cryo-EM) maps with atomic models requires high levels of expertise and labour-intensive manual intervention in three-dimensional computer graphics programs1,2. Here we present ModelAngelo, a machine-learning approach for automated atomic model building in cryo-EM maps. By combining information from the cryo-EM map with information from protein sequence and structure in a single graph neural network, ModelAngelo builds atomic models for proteins that are of similar quality to those generated by human experts. For nucleotides, ModelAngelo builds backbones with similar accuracy to those built by humans. By using its predicted amino acid probabilities for each residue in hidden Markov model sequence searches, ModelAngelo outperforms human experts in the identification of proteins with unknown sequences. ModelAngelo will therefore remove bottlenecks and increase objectivity in cryo-EM structure determination.

Structural basis of prostaglandin efflux by MRP4

Multidrug resistance protein 4 (MRP4) is a broadly expressed ATP-binding cassette transporter that is unique among the MRP subfamily for transporting prostanoids, a group of signaling molecules derived from unsaturated fatty acids. To better understand the basis of the substrate selectivity of MRP4, we used cryogenic-electron microscopy to determine six structures of nanodisc-reconstituted MRP4 at various stages throughout its transport cycle. Substrate-bound structures of MRP4 in complex with PGE1, PGE2 and the sulfonated-sterol DHEA-S reveal a common binding site that accommodates a diverse set of organic anions and suggest an allosteric mechanism for substrate-induced enhancement of MRP4 ATPase activity. Our structure of a catalytically compromised MRP4 mutant bound to ATP-Mg2+ is outward-occluded, a conformation previously unobserved in the MRP subfamily and consistent with an alternating-access transport mechanism. Our study provides insights into the endogenous function of this versatile efflux transporter and establishes a basis for MRP4-targeted drug design.

High-resolution cryo-EM of the human CDK-activating kinase for structure-based drug design

Rational design of next-generation therapeutics can be facilitated by high-resolution structures of drug targets bound to small-molecule inhibitors. However, application of structure-based methods to macromolecules refractory to crystallization has been hampered by the often-limiting resolution and throughput of cryogenic electron microscopy (cryo-EM). Here, we use high-resolution cryo-EM to determine structures of the CDK-activating kinase, a master regulator of cell growth and division, in its free and nucleotide-bound states and in complex with 15 inhibitors at up to 1.8 Å resolution. Our structures provide detailed insight into inhibitor interactions and networks of water molecules in the active site of cyclin-dependent kinase 7 and provide insights into the mechanisms contributing to inhibitor selectivity, thereby providing the basis for rational design of next-generation therapeutics. These results establish a methodological framework for the use of high-resolution cryo-EM in structure-based drug design.

Efficient Sampling of Cavity Hydration in Proteins with Nonequilibrium Grand Canonical Monte Carlo and Polarizable Force Fields

Prediction of the hydration levels of protein cavities and active sites is important to both mechanistic analysis and ligand design. Due to the unique microscopic environment of these buried water molecules, a polarizable model is expected to be crucial for an accurate treatment of protein internal hydration in simulations. Here we adapt a nonequilibrium candidate Monte Carlo approach for conducting grand canonical Monte Carlo simulations with the Drude polarizable force field. The GPU implementation enables the efficient sampling of internal cavity hydration levels in biomolecular systems. We also develop an enhanced sampling approach referred to as B-walking, which satisfies detailed balance and readily combines with grand canonical integration to efficiently calculate quantitative binding free energies of water to protein cavities. Applications of these developments are illustrated in a solvent box and the polar ligand binding site in trypsin. Our simulation results show that including electronic polarization leads to a modest but clear improvement in the description of water position and occupancy compared to the crystal structure. The B-walking approach enhances the range of water sampling in different chemical potential windows and thus improves the accuracy of water binding free energy calculations.

Structural basis for plasmid restriction by SMC JET nuclease

DNA loop-extruding SMC complexes play crucial roles in chromosome folding and DNA immunity. Prokaryotic SMC Wadjet (JET) complexes limit the spread of plasmids through DNA cleavage, yet the mechanisms for plasmid recognition are unresolved. We show that artificial DNA circularization renders linear DNA susceptible to JET nuclease cleavage. Unlike free DNA, JET cleaves immobilized plasmid DNA at a specific site, the plasmid-anchoring point, showing that the anchor hinders DNA extrusion but not DNA cleavage. Structures of plasmid-bound JetABC reveal two presumably stalled SMC motor units that are drastically rearranged from the resting state, together entrapping a U-shaped DNA segment, which is further converted to kinked V-shaped cleavage substrate by JetD nuclease binding. Our findings uncover mechanical bending of residual unextruded DNA as molecular signature for plasmid recognition and non-self DNA elimination. We moreover elucidate key elements of SMC loop extrusion, including the motor direction and the structure of a DNA-holding state.

Tertiary folds of the SL5 RNA from the 5' proximal region of SARS-CoV-2 and related coronaviruses

Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron microscopy (cryo-EM) and computational modeling with biochemically determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2), resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T's "arms." Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4 to 6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across these human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9 to 8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4 to 9.0 Å resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities and notable differences, with implications for potential protein-binding modes and therapeutic targets.

Cryogenic electron microscopy and tomography reveal imperfect icosahedral symmetry in alphaviruses

Alphaviruses are spherical, enveloped RNA viruses with two-layered icosahedral architecture. The structures of many alphaviruses have been studied using cryogenic electron microscopy (cryo-EM) reconstructions, which impose icosahedral symmetry on the viral particles. Using cryogenic electron tomography (cryo-ET), we revealed a polarized symmetry defect in the icosahedral lattice of Chikungunya virus (CHIKV) in situ, similar to the late budding particles, suggesting the inherent imperfect symmetry originates from the final pinch-off of assembled virions. We further demonstrated this imperfect symmetry is also present in in vitro purified CHIKV and Mayaro virus, another arthritogenic alphavirus. We employed a subparticle-based single-particle analysis protocol to circumvent the icosahedral imperfection and boosted the resolution of the structure of the CHIKV to ∼3 Å resolution, which revealed detailed molecular interactions between glycoprotein E1-E2 heterodimers in the transmembrane region and multiple lipid-like pocket factors located in a highly conserved hydrophobic pocket. This complementary use of in situ cryo-ET and single-particle cryo-EM approaches provides a more precise structural description of near-icosahedral viruses and valuable insights to guide the development of structure-based antiviral therapies against alphaviruses.

DeepQs: Local quality assessment of cryo-EM density map by deep learning map-model fit score

Cryogenic electron microscopy maps are valuable for determining macromolecule structures. A proper quality assessment method is essential for cryo-EM map selection or revision. This article presents DeepQs, a novel approach to estimate local quality for 3D cryo-EM density maps, using a deep-learning algorithm based on map-model fit score. DeepQs is a parameter-free method for users and incorporates structural information between map and its related atomic model into well-trained models by deep learning. More specifically, the DeepQs approach leverages the interplay between map and atomic model through predefined map-model fit score, Q-score. DeepQs can get close results to the ground truth map-model fit scores with only cryo-EM map as input. In experiments, DeepQs demonstrates the lowest root mean square error with standard method Fourier shell correlation metric and high correlation with map-model fit score, Q-score, when compared with other local quality estimation methods in high-resolution dataset (<=5 Å). DeepQs can also be applied to evaluate the quality of the post-processed maps. In both cases, DeepQs runs faster by using GPU acceleration. Our program is available at http://www.csbio.sjtu.edu.cn/bioinf/DeepQs for academic use.

Community recommendations on cryoEM data archiving and validation

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for the deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and the resulting consensus recommendations. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

Rapid simulation of glycoprotein structures by grafting and steric exclusion of glycan conformer libraries

Most membrane proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, shielding potentially large fractions of protein surface. High glycan conformational freedom hinders complete structural elucidation of glycoproteins. Computer simulations may be used to model glycosylated proteins but require hundreds of thousands of computing hours on supercomputers, thus limiting routine use. Here, we describe GlycoSHIELD, a reductionist method that can be implemented on personal computers to graft realistic ensembles of glycan conformers onto static protein structures in minutes. Using molecular dynamics simulation, small-angle X-ray scattering, cryoelectron microscopy, and mass spectrometry, we show that this open-access toolkit provides enhanced models of glycoprotein structures. Focusing on N-cadherin, human coronavirus spike proteins, and gamma-aminobutyric acid receptors, we show that GlycoSHIELD can shed light on the impact of glycans on the conformation and activity of complex glycoproteins.

High-resolution single-particle imaging at 100-200 keV with the Gatan Alpine direct electron detector

Developments in direct electron detector technology have played a pivotal role in enabling high-resolution structural studies by cryo-EM at 200 and 300 keV. Yet, theory and recent experiments indicate advantages to imaging at 100 keV, energies for which the current detectors have not been optimized. In this study, we evaluated the Gatan Alpine detector, designed for operation at 100 and 200 keV. Compared to the Gatan K3, Alpine demonstrated a significant DQE improvement at these voltages, specifically a ~4-fold improvement at Nyquist at 100 keV. In single-particle cryo-EM experiments, Alpine datasets yielded better than 2 Å resolution reconstructions of apoferritin at 120 and 200 keV on a ThermoFisher Scientific (TFS) Glacios microscope. We also achieved a ~3.2 Å resolution reconstruction for a 115 kDa asymmetric protein complex, proving its effectiveness with complex biological samples. In-depth analysis revealed that Alpine reconstructions are comparable to K3 reconstructions at 200 keV, and remarkably, reconstruction from Alpine at 120 keV on a TFS Glacios surpassed all but the 300 keV data from a TFS Titan Krios with GIF/K3. Additionally, we show Alpine's capability for high-resolution data acquisition and screening on lower-end systems by obtaining ~3 Å resolution reconstructions of apoferritin and aldolase at 100 keV and detailed 2D averages of a 55 kDa sample using a side-entry cryo holder. Overall, we show that Gatan Alpine performs well with the standard 200 keV imaging systems and may potentially capture the benefits of lower accelerating voltages, possibly bringing smaller sized particles within the scope of cryo-EM.

Mobile barrier mechanisms for Na+-coupled symport in an MFS sugar transporter

While many 3D structures of cation-coupled transporters have been determined, the mechanistic details governing the obligatory coupling and functional regulations still remain elusive. The bacterial melibiose transporter (MelB) is a prototype of major facilitator superfamily transporters. With a conformation-selective nanobody, we determined a low-sugar affinity inward-facing Na+-bound cryoEM structure. The available outward-facing sugar-bound structures showed that the N- and C-terminal residues of the inner barrier contribute to the sugar selectivity. The inward-open conformation shows that the sugar selectivity pocket is also broken when the inner barrier is broken. Isothermal titration calorimetry measurements revealed that this inward-facing conformation trapped by this nanobody exhibited a greatly decreased sugar-binding affinity, suggesting the mechanisms for substrate intracellular release and accumulation. While the inner/outer barrier shift directly regulates the sugar-binding affinity, it has little or no effect on the cation binding, which is supported by molecular dynamics simulations. Furthermore, the hydron/deuterium exchange mass spectrometry analyses allowed us to identify dynamic regions; some regions are involved in the functionally important inner barrier-specific salt-bridge network, which indicates their critical roles in the barrier switching mechanisms for transport. These complementary results provided structural and dynamic insights into the mobile barrier mechanism for cation-coupled symport.

CryoEM structures of the human CLC-2 voltage-gated chloride channel reveal a ball-and-chain gating mechanism

CLC-2 is a voltage-gated chloride channel that contributes to electrical excitability and ion homeostasis in many different tissues. Among the nine mammalian CLC homologs, CLC-2 is uniquely activated by hyperpolarization, rather than depolarization, of the plasma membrane. The molecular basis for the divergence in polarity of voltage gating among closely related homologs has been a long-standing mystery, in part because few CLC channel structures are available. Here, we report cryoEM structures of human CLC-2 at 2.46 - 2.76 Å, in the presence and absence of the selective inhibitor AK-42. AK-42 binds within the extracellular entryway of the Cl--permeation pathway, occupying a pocket previously proposed through computational docking studies. In the apo structure, we observed two distinct conformations involving rotation of one of the cytoplasmic C-terminal domains (CTDs). In the absence of CTD rotation, an intracellular N-terminal 15-residue hairpin peptide nestles against the TM domain to physically occlude the Cl--permeation pathway. This peptide is highly conserved among species variants of CLC-2 but is not present in other CLC homologs. Previous studies suggested that the N-terminal domain of CLC-2 influences channel properties via a "ball-and-chain" gating mechanism, but conflicting data cast doubt on such a mechanism, and thus the structure of the N-terminal domain and its interaction with the channel has been uncertain. Through electrophysiological studies of an N-terminal deletion mutant lacking the 15-residue hairpin peptide, we support a model in which the N-terminal hairpin of CLC-2 stabilizes a closed state of the channel by blocking the cytoplasmic Cl--permeation pathway.

Community recommendations on cryoEM data archiving and validation: Outcomes of a wwPDB/EMDB workshop on cryoEM data management, deposition and validation

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and consensus recommendations resulting from the workshop. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

A structural vista of phosducin-like PhLP2A-chaperonin TRiC cooperation during the ATP-driven folding cycle

Proper cellular proteostasis, essential for viability, requires a network of chaperones and cochaperones. ATP-dependent chaperonin TRiC/CCT partners with cochaperones prefoldin (PFD) and phosducin-like proteins (PhLPs) to facilitate folding of essential eukaryotic proteins. Using cryoEM and biochemical analyses, we determine the ATP-driven cycle of TRiC-PFD-PhLP2A interaction. PhLP2A binds to open apo-TRiC through polyvalent domain-specific contacts with its chamber's equatorial and apical regions. PhLP2A N-terminal H3-domain binding to subunits CCT3/4 apical domains displace PFD from TRiC. ATP-induced TRiC closure rearranges the contacts of PhLP2A domains within the closed chamber. In the presence of substrate, actin and PhLP2A segregate into opposing chambers, each binding to positively charged inner surface residues from CCT1/3/6/8. Notably, actin induces a conformational change in PhLP2A, causing its N-terminal helices to extend across the inter-ring interface to directly contact a hydrophobic groove in actin. Our findings reveal an ATP-driven PhLP2A structural rearrangement cycle within the TRiC chamber to facilitate folding.

2.7 Å cryo-EM structure of human telomerase H/ACA ribonucleoprotein

Telomerase is a ribonucleoprotein (RNP) enzyme that extends telomeric repeats at eukaryotic chromosome ends to counterbalance telomere loss caused by incomplete genome replication. Human telomerase is comprised of two distinct functional lobes tethered by telomerase RNA (hTR): a catalytic core, responsible for DNA extension; and a Hinge and ACA (H/ACA) box RNP, responsible for telomerase biogenesis. H/ACA RNPs also have a general role in pseudouridylation of spliceosomal and ribosomal RNAs, which is critical for the biogenesis of the spliceosome and ribosome. Much of our structural understanding of eukaryotic H/ACA RNPs comes from structures of the human telomerase H/ACA RNP. Here we report a 2.7 Å cryo-electron microscopy structure of the telomerase H/ACA RNP. The significant improvement in resolution over previous 3.3 Å to 8.2 Å structures allows us to uncover new molecular interactions within the H/ACA RNP. Many disease mutations are mapped to these interaction sites. The structure also reveals unprecedented insights into a region critical for pseudouridylation in canonical H/ACA RNPs. Together, our work advances understanding of telomerase-related disease mutations and the mechanism of pseudouridylation by eukaryotic H/ACA RNPs.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

The conformational landscape of human transthyretin revealed by cryo-EM

Transthyretin (TTR) is a natively tetrameric thyroxine transporter found in blood and cerebrospinal fluid whose misfolding and aggregation causes transthyretin amyloidosis. A rational drug design campaign identified the small molecule tafamidis (Vyndaqel/Vyndamax) as an effective stabilizer of the native TTR fold, and this aggregation inhibitor is regulatory agency-approved for the treatment of TTR amyloidosis. Despite 50 years of structural studies on TTR and this triumph of structure-based drug design, there remains a notable dearth of structural information available to understand ligand binding allostery and amyloidogenic TTR unfolding intermediates. We used single-particle cryo-electron microscopy (cryo-EM) to investigate the conformational landscape of this 55 kiloDalton tetramer in the absence and presence of one or two ligands, revealing inherent asymmetries in the tetrameric architecture and previously unobserved conformational states. These findings provide critical mechanistic insights into negatively cooperative ligand binding and the structural pathways responsible for TTR amyloidogenesis. This study underscores the capacity of cryo-EM to provide new insights into protein structures that have been historically considered too small to visualize and to identify pharmacological targets suppressed by the confines of the crystal lattice, opening uncharted territory in structure-based drug design.

ChemEM: Flexible Docking of Small Molecules in Cryo-EM Structures

Cryo-electron microscopy (cryo-EM), through resolution advancements, has become pivotal in structure-based drug discovery. However, most cryo-EM structures are solved at 3-4 Å resolution, posing challenges for small-molecule docking and structure-based virtual screening due to issues in the precise positioning of ligands and the surrounding side chains. We present ChemEM, a software package that employs cryo-EM data for the accurate docking of one or multiple ligands in a protein-binding site. Validated against a highly curated benchmark of high- and medium-resolution cryo-EM structures and the corresponding high-resolution controls, ChemEM displayed impressive performance, accurately placing ligands in all but one case, often surpassing cryo-EM PDB-deposited solutions. Even without including the cryo-EM density, the ChemEM scoring function outperformed the well-established AutoDock Vina score. Using ChemEM, we illustrate that valuable information can be extracted from maps at medium resolution and underline the utility of cryo-EM structures for drug discovery.

Exploiting the full potential of cryo-EM maps

The Coulomb potential maps generated by electron microscopy (EM) experiments contain not only information about the position but also about the charge state of the atom. This feature of EM maps allows the identification of specific ions and the protonation state of amino acid side chains in the sample. Here, we summarize qualitative observations of charges in EM maps, discuss the difficulties in interpreting the charge in Coulomb potential maps with respect to distinguishing it from radiation damage, and outline considerations to implement the correct charge in fitting algorithms.

Genome-wide mapping and cryo-EM structural analyses of the overlapping tri-nucleosome composed of hexasome-hexasome-octasome moieties

The nucleosome is a fundamental unit of chromatin in which about 150 base pairs of DNA are wrapped around a histone octamer. The overlapping di-nucleosome has been proposed as a product of chromatin remodeling around the transcription start site, and previously found as a chromatin unit, in which about 250 base pairs of DNA continuously bind to the histone core composed of a hexamer and an octamer. In the present study, our genome-wide analysis of human cells suggests another higher nucleosome stacking structure, the overlapping tri-nucleosome, which wraps about 300-350 base-pairs of DNA in the region downstream of certain transcription start sites of actively transcribed genes. We determine the cryo-electron microscopy (cryo-EM) structure of the overlapping tri-nucleosome, in which three subnucleosome moieties, hexasome, hexasome, and octasome, are associated by short connecting DNA segments. Small angle X-ray scattering and coarse-grained molecular dynamics simulation analyses reveal that the cryo-EM structure of the overlapping tri-nucleosome may reflect its structure in solution. Our findings suggest that nucleosome stacking structures composed of hexasome and octasome moieties may be formed by nucleosome remodeling factors around transcription start sites for gene regulation.

EMDB-the Electron Microscopy Data Bank

The Electron Microscopy Data Bank (EMDB) is the global public archive of three-dimensional electron microscopy (3DEM) maps of biological specimens derived from transmission electron microscopy experiments. As of 2021, EMDB is managed by the Worldwide Protein Data Bank consortium (wwPDB; wwpdb.org) as a wwPDB Core Archive, and the EMDB team is a core member of the consortium. Today, EMDB houses over 30 000 entries with maps containing macromolecules, complexes, viruses, organelles and cells. Herein, we provide an overview of the rapidly growing EMDB archive, including its current holdings, recent updates, and future plans.

Improving resolution and resolvability of single-particle cryoEM structures using Gaussian mixture models

Cryogenic electron microscopy is widely used in structural biology, but its resolution is often limited by the dynamics of the macromolecule. Here we developed a refinement protocol based on Gaussian mixture models that integrates particle orientation and conformation estimation and improves the alignment for flexible domains of protein structures. We demonstrated this protocol on multiple datasets, resulting in improved resolution and resolvability, locally and globally, by visual and quantitative measures.