Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

Tension-induced suppression of allosteric conformational changes coordinates kinesin-1 stepping

Kinesin-1 walks along microtubules by alternating ATP hydrolysis and movement of its two motor domains ("head"). The detached head preferentially binds to the forward tubulin-binding site after ATP binds to the microtubule-bound head, but the mechanism preventing premature microtubule binding while the partner head awaits ATP remains unknown. Here, we examined the role of the neck linker, the segment connecting two heads, in this mechanism. Structural analyses of the nucleotide-free head revealed a bulge just ahead of the neck linker's base, creating an asymmetric constraint on its mobility. While the neck linker can stretch freely backward, it must navigate around this bulge to extend forward. We hypothesized that increased neck linker tension suppresses premature binding of the tethered head, which was supported by molecular dynamics simulations and single-molecule fluorescence assays. These findings demonstrate a tension-dependent allosteric mechanism that coordinates the movement of two heads, where neck linker tension modulates the allosteric conformational changes rather than directly affecting the nucleotide state.

The Structure of the Drp1 Lattice on Membrane

Mitochondrial health relies on the membrane fission mediated by dynamin-related protein 1 (Drp1). Previous structural studies of Drp1 on remodeled membranes were hampered by heterogeneity, leaving a critical gap in the understanding of the mitochondrial fission mechanisms. Here we present a cryo-electron microscopy structure of full-length human Drp1 decorated on membrane tubules. Using the reconstruction of average subtracted tubular regions (RASTR) technique, we report that Drp1 forms a locally ordered lattice along the tubule without global helical symmetry. The filaments in the lattice are similar to dynamin rungs with conserved stalk interactions. Adjacent filaments are connected by GTPase domain interactions in a novel stacked conformation. We identified two states of the Drp1 lattice among the heterogenous dataset representing conformational changes around hinge 1. Additionally, we observed contact between Drp1 and membrane that can be assigned to the variable domain sequence. Together these structures revealed a putative mechanism by which Drp1 constricts mitochondria membranes in a stepwise, "ratchet" manner.

Dimerization of GAS2 mediates crosslinking of microtubules and F-actin

The spectraplakin family protein GAS2 was originally identified as a growth arrest-specific protein, and recent studies have revealed its involvement in multiple cellular processes. Its dual interaction with actin filaments and microtubules highlights its essential role in cytoskeletal organization, such as cell division, apoptosis, and possibly tumorigenesis. However, the structural basis of cytoskeletal dynamics regulation by GAS2 remains unclear. In this study, we present cryo-electron microscopy structures of the GAS2 type 3 calponin homology domain (CH3) in complex with F-actin at 2.8 Å resolution, thus solving the first type CH3 domain structure bound to F-actin and confirming its actin-binding activity. We also provide the first near-atomic resolution cryo-EM structure of the GAS2-GAR domain bound to microtubules and identify conserved microtubule-binding residues. Our biochemical experiments show that GAS2 promotes microtubule nucleation and polymerization, and that its C-terminal region is essential for dimerization, bundling of both F-actin and microtubules, and microtubule nucleation. As mutations leading to expression of C-terminally truncated GAS2 have been linked to hearing loss, these findings suggest that the disruption of GAS2-dependent cytoskeletal organisation could underlie auditory dysfunction.

Dissecting Large-Scale Structural Transitions in Membrane Transporters Using Advanced Simulation Technologies

Membrane transporters are integral membrane proteins that act as gatekeepers of the cell, controlling fundamental processes such as recruitment of nutrients and expulsion of waste material. At a basic level, transporters operate using the "alternating access model," in which transported substances are accessible from only one side of the membrane at a time. This model usually involves large-scale structural changes in the transporter, which often cannot be captured using unbiased, conventional molecular simulation techniques. In this article, we provide an overview of some of the major simulation techniques that have been applied to characterize the structural dynamics and energetics involved in the transition of membrane transporters between their functional states. After briefly introducing each technique, we discuss some of their advantages and limitations and provide some recent examples of their application to membrane transporters.

Unveiling Cas8 dynamics and regulation within a transposon-encoded Cascade-TniQ complex

The Vibrio cholerae Cascade-TniQ complex unveiled a new paradigm in biology, demonstrating that CRISPR-associated proteins can direct DNA transposition. Despite the tremendous potential of "knocking-in" genes at desired sites, the mechanisms underlying DNA binding and transposition remain elusive. In this system, a conformational change of the Cas8 protein is essential for DNA binding, yet how it occurs is unclear. Here, structural modeling and free energy simulations reconstruct the Cas8 helical bundle and reveal an open-closed conformational change that is key for the complex's function. We show that when Cascade-TniQ binds RNA, the Cas8 bundle changes conformation mediated by the interaction with the Cas7.1 protein. This interaction promotes the bundle's transition toward the open state, priming the complex for DNA binding. As the target DNA binds the guide RNA, the opening of the Cas8 bundle becomes more favorable, exposing positively charged residues and facilitating their interaction with DNA, which ultimately leads the DNA-binding process to completion. These outcomes provide a dynamic representation of a critical conformational change in one of the largest CRISPR systems and illustrate its role at critical steps of the Cascade-TniQ biophysical function, advancing our understanding of nucleic acid binding and transposition mechanisms.

Structural basis for the substrate recognition and transport mechanism of the human y+LAT1-4F2hc transporter complex

Heteromeric amino acid transporters (HATs), including y+LAT1-4F2hc complex, are responsible for transporting amino acids across membranes, and mutations in y+LAT1 cause lysinuric protein intolerance (LPI), a hereditary disorder characterized by defective cationic amino acid transport. The relationship between LPI and specific mutations in y+LAT1 has yet to be fully understood. In this study, we characterized the function of y+LAT1-4F2hc complex in mammalian cells and determined the cryo-EM structures of the human y+LAT1-4F2hc complex in two distinct conformations: the apo state in an inward-open conformation and the native substrate-bound state in an outward-open conformation. Structural analysis suggests that Asp243 in y+LAT1 plays a crucial role in coordination with sodium ion and substrate selectivity. Molecular dynamic (MD) simulations further revealed the different transport mechanism of cationic amino acids and neutral amino acids. These results provide important insights into the mechanisms of the substrate binding and working cycle of HATs.

Integrative Protein Assembly With LZerD and Deep Learning in CAPRI 47-55

We report the performance of the protein complex prediction approaches of our group and their results in CAPRI Rounds 47-55, excluding the joint CASP Rounds 50 and 54, as well as the special COVID-19 Round 51. Our approaches integrated classical pipelines developed in our group as well as more recently developed deep learning pipelines. In the cases of human group prediction, we surveyed the literature to find information to integrate into the modeling, such as assayed interface residues. In addition to any literature information, generated complex models were selected by a rank aggregation of statistical scoring functions, by generative model confidence, or by expert inspection. In these CAPRI rounds, our human group successfully modeled eight interfaces and achieved the top quality level among the submissions for all of them, including two where no other group did. We note that components of our modeling pipelines have become increasingly unified within deep learning approaches. Finally, we discuss several case studies that illustrate successful and unsuccessful modeling using our approaches.

In-cell Structure and Variability of Pyrenoid Rubisco

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) is a key enzyme in the global carbon cycle, catalyzing CO2 fixation during photosynthesis. To overcome Rubisco's inherent catalytic inefficiency, many photosynthetic organisms have evolved CO2-concentrating mechanisms. Central to these mechanisms is the pyrenoid, a protein-dense organelle within the chloroplast of eukaryotic algae, which increases the local concentration of CO2 around Rubisco and thereby enhances its catalytic efficiency. Although the structure of Rubisco has been extensively studied by in vitro methods such as X-ray crystallography and single particle cryo-EM, its native structure within the pyrenoid, its dynamics, and its association with binding partners remain elusive. Here, we investigate the structure of native pyrenoid Rubisco inside the green alga Chlamydomonas reinhardtii by applying cryo-electron tomography (cryo-ET) on cryo-focused ion beam (cryo-FIB) milled cells, followed by subtomogram averaging and 3D classification. Reconstruction at sub-nanometer resolution allowed accurate modeling and determination of a closed (activated) Rubisco conformation. Comparison to other reconstructed subsets revealed local variations at the complex active site and at the large subunit dimers interface, as well as association with binding proteins. The different structural subsets distribute stochastically within the pyrenoid. Taken together, these findings offer a comprehensive description of the structure, dynamics, and functional organization of Rubisco within the pyrenoid, providing valuable insights into its critical role in CO2 fixation.

The Structure of the Drp1 Lattice on Membrane

Mitochondrial health relies on the membrane fission mediated by dynamin-related protein 1 (Drp1). Previous structural studies of Drp1 on remodeled membranes were hampered by heterogeneity, leaving a critical gap in the understanding of the mitochondrial fission mechanisms. Here we present a cryo-electron microscopy structure of full-length human Drp1 decorated on membrane tubules. Using the reconstruction of average subtracted tubular regions (RASTR) technique, we report that Drp1 forms a locally ordered lattice along the tubule without global helical symmetry. The filaments in the lattice are similar to dynamin rungs with conserved stalk interactions. Adjacent filaments are connected by GTPase domain interactions in a novel stacked conformation. We identified two states of the Drp1 lattice among the heterogenous dataset representing conformational changes around hinge 1. Additionally, we observed contact between Drp1 and membrane that can be assigned to the variable domain sequence. Together these structures revealed a putative mechanism by which Drp1 constricts mitochondria membranes in a stepwise, "ratchet" manner.

A structural biology compatible file format for atomic force microscopy

Cryogenic electron microscopy (cryo-EM), X-ray crystallography, and nuclear magnetic resonance (NMR) contribute structural data that are interchangeable, cross-verifiable, and visualizable on common platforms, making them powerful tools for our understanding of protein structures. Unfortunately, atomic force microscopy (AFM) has so far failed to interface with these structural biology methods, despite the recent development of localization AFM (LAFM) that allows extracting high-resolution structural information from AFM data. Here, we build on LAFM and develop a pipeline that transforms AFM data into 3D-density files (.afm) that are readable by programs commonly used to visualize, analyze, and interpret structural data. We show that 3D-LAFM densities can serve as force fields to steer molecular dynamics flexible fitting (MDFF) to obtain structural models of previously unresolved states based on AFM observations in close-to-native environment. Besides, the .afm format enables direct 3D or 2D visualization and analysis of conventional AFM images. We anticipate that the file format will find wide usage and embed AFM in the repertoire of methods routinely used by the structural biology community, allowing AFM researchers to deposit data in repositories in a format that allows comparison and cross-verification with data from other techniques.

Structural insights into dynamics of the BMV TLS aminoacylation

Brome Mosaic Virus (BMV) utilizes a tRNA-like structure (TLS) within its 3' untranslated region to mimic host tRNA functions, aiding aminoacylation and viral replication. This study explores the structural dynamics of BMV TLS interacting with tyrosyl-tRNA synthetase (TyrRS) during aminoacylation. Using cryo-EM, we capture multiple states of the TLS-TyrRS complex, including unbound TLS, pre-1a, post-1a, and catalysis states, with resolutions of 4.6 Å, 3.5 Å, 3.7 Å, and 3.85 Å, respectively. These structural comparisons indicate dynamic changes in both TLS and TyrRS. Upon binding, TLS undergoes dynamic rearrangements, particularly with helices B3 and E pivoting, mediated by the unpaired A36 residue, ensuring effective recognition by TyrRS. The dynamic changes also include a more compact arrangement in the catalytic center of TyrRS and the insertion of 3' CCA end into the enzyme's active site, facilitating two-steps aminoacylation. Enzymatic assays further demonstrated the functional importance of TLS-TyrRS interactions, with mutations in key residues significantly impacting aminoacylation efficiency. Furthermore, Electrophoretic Mobility Shift Assay (EMSA) demonstrated that BMV TLS binds elongation factors EF1α and EF2, suggesting a multifaceted strategy to exploit host translational machinery. These findings not only enhance our knowledge of virus-host interactions but also offer potential targets for antiviral drug development.

The structure of apolipoprotein B100 from human low-density lipoprotein

Low-density lipoprotein (LDL) has a central role in lipid and cholesterol metabolism and is a key agent in the development and progression of atherosclerosis, the leading cause of mortality worldwide1,2. Apolipoprotein B100 (apoB100), one of the largest proteins in the genome, is the primary structural and functional component of LDL, yet its size and complex lipid associations have posed major challenges for structural studies3. Here we present the structure of apoB100 resolved to subnanometre resolution in most regions using an integrative approach of cryo-electron microscopy, AlphaFold24 and molecular-dynamics-based refinement5. The structure consists of a large globular N-terminal domain and an approximately 61-nm-long continuous amphipathic β-sheet that wraps around the LDL particle like a belt. Distributed quasi-symmetrically across the two sides of the β-belt are nine strategically located interstrand inserts that extend across the lipid surface to provide additional structural support through a network of long-range interactions. We further compare our structure to a comprehensive list of more than 200 intramolecular cross-links and find close agreement between the two. These results suggest a mechanism for how the various domains of apoB100 act in concert to maintain LDL shape and cohesion across a range of particle sizes. More generally, they advance our fundamental understanding of LDL synthesis, form and function, and will help to accelerate the design of potential therapeutics.

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.

3D variability analysis reveals a hidden conformational change controlling ammonia transport in human asparagine synthetase

Advances in X-ray crystallography and cryogenic electron microscopy (cryo-EM) offer the promise of elucidating functionally relevant conformational changes that are not easily studied by other biophysical methods. Here we show that 3D variability analysis (3DVA) of the cryo-EM map for wild-type (WT) human asparagine synthetase (ASNS) identifies a functional role for the Arg-142 side chain and test this hypothesis experimentally by characterizing the R142I variant in which Arg-142 is replaced by isoleucine. Support for Arg-142 playing a role in the intramolecular translocation of ammonia between the active site of the enzyme is provided by the glutamine-dependent synthetase activity of the R142 variant relative to WT ASNS, and MD simulations provide a possible molecular mechanism for these findings. Combining 3DVA with MD simulations is a generally applicable approach to generate testable hypotheses of how conformational changes in buried side chains might regulate function in enzymes.

A novel bispecific antibody targeting two overlapping epitopes in RBD improves neutralizing potency and breadth against SARS-CoV-2

SARS-CoV-2 has been evolving into a large number of variants, including the highly pathogenic Delta variant, and the currently prevalent Omicron subvariants with extensive evasion capability, which raises an urgent need to develop new broad-spectrum neutralizing antibodies. Herein, we engineer two IgG-(scFv)2 form bispecific antibodies with overlapping epitopes (bsAb1) or non-overlapping epitopes (bsAb2). Both bsAbs are significantly superior to the parental monoclonal antibodies in terms of their antigen-binding and virus-neutralizing activities against all tested circulating SARS-CoV-2 variants including currently dominant JN.1. The bsAb1 can efficiently neutralize all variants insensitive to parental monoclonal antibodies or the cocktail with IC50 lower than 20 ng/mL, even slightly better than bsAb2. Furthermore, the cryo-EM structures of bsAb1 in complex with the Omicron spike protein revealed that bsAb1 with overlapping epitopes effectively locked the S protein, which accounts for its conserved neutralization against Omicron variants. The bispecific antibody strategy engineered from overlapping epitopes provides a novel solution for dealing with viral immune evasion.

Advancing Cancer Immunotherapy through Engineering New PD-L1 Degraders: A Comprehensive Study from Small Molecules to PD-L1-Specific Peptide-Drug Conjugates

Despite the considerable achievements of antibodies targeting PD-1/PD-L1 in cancer immunotherapy, limitations in antitumor immune response and pharmacokinetics hinder their clinical adoption. Small molecules toward PD-L1 degradation signifies an innovative avenue to modulate PD-1/PD-L1 axis. Herein, we unveil a comprehensive engineering involving the development of new PD-L1 degraders based on the berberine (BBR) and palmatine (PMT) bioactive frameworks and explore their translational potential for cancer immunotherapy using a peptide-drug conjugate strategy. Chemical modifications at the O-9 position of PMT dramatically enhance the PD-L1 degradation capacity. Further conjugation of PMT degraders with an anti-PD-L1 peptide featuring disulfide linkers enables efficient GSH-specific prodrug activation, yielding synergistic immunotherapeutic benefits through both external PD-L1 blockade and internal PD-L1 degradation mechanisms. This work elucidates the compelling charm of the discovery and application of PD-L1 degraders, offering solutions to the challenges in advancing cancer immunotherapy in widespread clinics.

Non-averaged single-molecule tertiary structures reveal RNA self-folding through individual-particle cryo-electron tomography

Large-scale and continuous conformational changes in the RNA self-folding process present significant challenges for structural studies, often requiring trade-offs between resolution and observational scope. Here, we utilize individual-particle cryo-electron tomography (IPET) to examine the post-transcriptional self-folding process of designed RNA origami 6-helix bundle with a clasp helix (6HBC). By avoiding selection, classification, averaging, or chemical fixation and optimizing cryo-ET data acquisition parameters, we reconstruct 120 three-dimensional (3D) density maps from 120 individual particles at an electron dose of no more than 168 e-Å-2, achieving averaged resolutions ranging from 23 to 35 Å, as estimated by Fourier shell correlation (FSC) at 0.5. Each map allows us to identify distinct RNA helices and determine a unique tertiary structure. Statistical analysis of these 120 structures confirms two reported conformations and reveals a range of kinetically trapped, intermediate, and highly compacted states, demonstrating a maturation folding landscape likely driven by helix-helix compaction interactions.

Structural basis for the evolution and antibody evasion of SARS-CoV-2 BA.2.86 and JN.1 subvariants

The Omicron subvariants of SARS-CoV-2, especially for BA.2.86 and JN.1, have rapidly spread across multiple countries, posing a significant threat in the ongoing COVID-19 pandemic. Distinguished by 34 additional mutations on the Spike (S) protein compared to its BA.2 predecessor, the implications of BA.2.86 and its evolved descendant, JN.1 with additional L455S mutation in receptor-binding domains (RBDs), are of paramount concern. In this work, we systematically examine the neutralization susceptibilities of SARS-CoV-2 Omicron subvariants and reveal the enhanced antibody evasion of BA.2.86 and JN.1. We also determine the cryo-EM structures of the trimeric S proteins from BA.2.86 and JN.1 in complex with the host receptor ACE2, respectively. The mutations within the RBDs of BA.2.86 and JN.1 induce a remodeling of the interaction network between the RBD and ACE2. The L455S mutation of JN.1 further induces a notable shift of the RBD-ACE2 interface, suggesting the notably reduced binding affinity of JN.1 than BA.2.86. An analysis of the broadly neutralizing antibodies possessing core neutralizing epitopes reveals the antibody evasion mechanism underlying the evolution of Omicron BA.2.86 subvariant. In general, we construct a landscape of evolution in virus-receptor of the circulating Omicron subvariants.

Outer membrane protein assembly mediated by BAM-SurA complexes

The outer membrane is a formidable barrier that protects Gram-negative bacteria against environmental threats. Its integrity requires the correct folding and insertion of outer membrane proteins (OMPs) by the membrane-embedded β-barrel assembly machinery (BAM). Unfolded OMPs are delivered to BAM by the periplasmic chaperone SurA, but how SurA and BAM work together to ensure successful OMP delivery and folding remains unclear. Here, guided by AlphaFold2 models, we use disulphide bond engineering in an attempt to trap SurA in the act of OMP delivery to BAM, and solve cryoEM structures of a series of complexes. The results suggest that SurA binds BAM at its soluble POTRA-1 domain, which may trigger conformational changes in both BAM and SurA that enable transfer of the unfolded OMP to the BAM lateral gate for insertion into the outer membrane. Mutations that disrupt the interaction between BAM and SurA result in outer membrane assembly defects, supporting the key role of SurA in outer membrane biogenesis.

Tumor-agnostic cancer therapy using antibodies targeting oncofetal chondroitin sulfate

Molecular similarities between embryonic and malignant cells can be exploited to target tumors through specific signatures absent in healthy adult tissues. One such embryonic signature tumors express is oncofetal chondroitin sulfate (ofCS), which supports disease progression and dissemination in cancer. Here, we report the identification and characterization of phage display-derived antibody fragments recognizing two distinct ofCS epitopes. These antibody fragments show binding affinity to ofCS in the low nanomolar range across a broad selection of solid tumor types in vitro and in vivo with minimal binding to normal, inflamed, or benign tumor tissues. Anti-ofCS antibody drug conjugates and bispecific immune cell engagers based on these targeting moieties disrupt tumor progression in animal models of human and murine cancers. Thus, anti-ofCS antibody fragments hold promise for the development of broadly effective therapeutic and diagnostic applications targeting human malignancies.

Potassium dependent structural changes in the selectivity filter of HERG potassium channels

The fine tuning of biological electrical signaling is mediated by variations in the rates of opening and closing of gates that control ion flux through different ion channels. Human ether-a-go-go related gene (HERG) potassium channels have uniquely rapid inactivation kinetics which are critical to the role they play in regulating cardiac electrical activity. Here, we exploit the K+ sensitivity of HERG inactivation to determine structures of both a conductive and non-conductive selectivity filter structure of HERG. The conductive state has a canonical cylindrical shaped selectivity filter. The non-conductive state is characterized by flipping of the selectivity filter valine backbone carbonyls to point away from the central axis. The side chain of S620 on the pore helix plays a central role in this process, by coordinating distinct sets of interactions in the conductive, non-conductive, and transition states. Our model represents a distinct mechanism by which ion channels fine tune their activity and could explain the uniquely rapid inactivation kinetics of HERG.

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria that encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo, the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We used molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo. Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed the overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceum model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

Bayesian reweighting of biomolecular structural ensembles using heterogeneous cryo-EM maps with the cryoENsemble method

Cryogenic electron microscopy (cryo-EM) has emerged as a powerful method for the determination of structures of complex biological molecules. The accurate characterisation of the dynamics of such systems, however, remains a challenge. To address this problem, we introduce cryoENsemble, a method that applies Bayesian reweighting to conformational ensembles derived from molecular dynamics simulations to improve their agreement with cryo-EM data, thus enabling the extraction of dynamics information. We illustrate the use of cryoENsemble to determine the dynamics of the ribosome-bound state of the co-translational chaperone trigger factor (TF). We also show that cryoENsemble can assist with the interpretation of low-resolution, noisy or unaccounted regions of cryo-EM maps. Notably, we are able to link an unaccounted part of the cryo-EM map to the presence of another protein (methionine aminopeptidase, or MetAP), rather than to the dynamics of TF, and model its TF-bound state. Based on these results, we anticipate that cryoENsemble will find use for challenging heterogeneous cryo-EM maps for biomolecular systems encompassing dynamic components.

Structure of adenylyl cyclase 5 in complex with Gβγ offers insights into ADCY5-related dyskinesia

The nine different membrane-anchored adenylyl cyclase isoforms (AC1-9) in mammals are stimulated by the heterotrimeric G protein, Gαs, but their response to Gβγ regulation is isoform specific. In the present study, we report cryo-electron microscope structures of ligand-free AC5 in complex with Gβγ and a dimeric form of AC5 that could be involved in its regulation. Gβγ binds to a coiled-coil domain that links the AC transmembrane region to its catalytic core as well as to a region (C1b) that is known to be a hub for isoform-specific regulation. We confirmed the Gβγ interaction with both purified proteins and cell-based assays. Gain-of-function mutations in AC5 associated with human familial dyskinesia are located at the interface of AC5 with Gβγ and show reduced conditional activation by Gβγ, emphasizing the importance of the observed interaction for motor function in humans. We propose a molecular mechanism wherein Gβγ either prevents dimerization of AC5 or allosterically modulates the coiled-coil domain, and hence the catalytic core. As our mechanistic understanding of how individual AC isoforms are uniquely regulated is limited, studies such as this may provide new avenues for isoform-specific drug development.

Molecular basis for Gβγ-mediated activation of phosphoinositide 3-kinase γ

The conversion of phosphatidylinositol 4,5-bisphosphate to phosphatidylinositol 3,4,5-triphosphate by phosphoinositide 3-kinase γ (PI3Kγ) is critical for neutrophil chemotaxis and cancer metastasis. PI3Kγ is activated by Gβγ heterodimers released from G protein-coupled receptors responding to extracellular signals. Here we determined cryo-electron microscopy structures of Sus scrofa PI3Kγ-human Gβγ complexes in the presence of substrates/analogs, revealing two Gβγ binding sites: one on the p110γ helical domain and another on the p101 C-terminal domain. Comparison with PI3Kγ alone reveals conformational changes in the kinase domain upon Gβγ binding that are similar to Ras·GTP-induced changes. Assays of variants perturbing the Gβγ binding sites and interdomain contacts altered by Gβγ binding suggest that Gβγ recruits the enzyme to membranes and allosterically regulates activity via both sites. Studies of zebrafish neutrophil migration align with these findings, paving the way for in-depth investigation of Gβγ-mediated activation mechanisms in this enzyme family and drug development for PI3Kγ.

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.

GOLEM: Automated and Robust Cryo-EM-Guided Ligand Docking with Explicit Water Molecules

A detailed understanding of ligand-protein interaction is essential for developing rational drug-design strategies. In recent years, technological advances in cryo-electron microscopy (cryo-EM) brought a new era to the structural determination of biological macromolecules and assemblies at high resolution, marking cryo-EM as a promising tool for studying ligand-protein interactions. However, even in high-resolution cryo-EM results, the densities for the bound small-molecule ligands are often of lower quality due to their relatively dynamic and flexible nature, frustrating their accurate coordinate assignment. To address the challenge of ligand modeling in cryo-EM maps, here we report the development of GOLEM (Genetic Optimization of Ligands in Experimental Maps), an automated and robust ligand docking method that predicts a ligand's pose and conformation in cryo-EM maps. GOLEM employs a Lamarckian genetic algorithm to perform a hybrid global/local search for exploring the ligand's conformational, orientational, and positional space. As an important feature, GOLEM explicitly considers water molecules and places them at optimal positions and orientations. GOLEM takes into account both molecular energetics and the correlation with the cryo-EM maps in its scoring function to optimally place the ligand. We have validated GOLEM against multiple cryo-EM structures with a wide range of map resolutions and ligand types, returning ligand poses in excellent agreement with the densities. As a VMD plugin, GOLEM is free of charge and accessible to the community. With these features, GOLEM will provide a valuable tool for ligand modeling in cryo-EM efforts toward drug discovery.

MDFF_NM: Improved Molecular Dynamics Flexible Fitting into Cryo-EM Density Maps with a Multireplica Normal Mode-Based Search

Molecular Dynamics Flexible Fitting (MDFF) is a widely used tool to refine high-resolution structures into cryo-EM density maps. Despite many successful applications, MDFF is still limited by its high computational cost, overfitting, accuracy, and performance issues due to entrapment within wrong local minima. Modern ensemble-based MDFF tools have generated promising results in the past decade. In line with these studies, we present MDFF_NM, a stochastic hybrid flexible fitting algorithm combining Normal Mode Analysis (NMA) and simulation-based flexible fitting. Initial tests reveal that, besides accelerating the fitting process, MDFF_NM increases the diversity of fitting routes leading to the target, uncovering ensembles of conformations in closer agreement with experimental data. The potential integration of MDFF_NM with other existing methods and integrative modeling pipelines is also discussed.

SANS reveals lipid-dependent oligomerization of an intramembrane aspartyl protease from H. volcanii

Reactions that occur within the lipid membrane involve, at minimum, ternary complexes among the enzyme, substrate, and lipid. For many systems, the impact of the lipid in regulating activity or oligomerization state is poorly understood. Here, we used small-angle neutron scattering (SANS) to structurally characterize an intramembrane aspartyl protease (IAP), a class of membrane-bound enzymes that use membrane-embedded aspartate residues to hydrolyze transmembrane segments of biologically relevant substrates. We focused on an IAP ortholog from the halophilic archaeon Haloferax volcanii (HvoIAP). HvoIAP purified in n-dodecyl-β-D-maltoside (DDM) fractionates on size-exclusion chromatography (SEC) as two fractions. We show that, in DDM, the smaller SEC fraction is consistent with a compact HvoIAP monomer. Molecular dynamics flexible fitting conducted on an AlphaFold2-generated monomer produces a model in which loops are compact alongside the membrane-embedded helices. In contrast, SANS data collected on the second SEC fraction indicate an oligomer consistent with an elongated assembly of discrete HvoIAP monomers. Analysis of in-line SEC-SANS data of the HvoIAP oligomer, the first such experiment to be conducted on a membrane protein at Oak Ridge National Lab (ORNL), shows a diversity of elongated and spherical species, including one consistent with the tetrameric assembly reported for the Methanoculleus marisnigri JR1 IAP crystal structure not observed previously in solution. Reconstitution of monomeric HvoIAP into bicelles increases enzyme activity and results in the assembly of HvoIAP into a species with similar dimensions as the ensemble of oligomers isolated from DDM. Our study reveals lipid-mediated HvoIAP self-assembly and demonstrates the utility of in-line SEC-SANS in elucidating oligomerization states of small membrane proteins.

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.

Structural basis for the inhibition mechanism of the DNA polymerase holoenzyme from mpox virus

There are three key components at the core of the mpox virus (MPXV) DNA polymerase holoenzyme: DNA polymerase F8, processivity factors A22, and the Uracil-DNA glycosylase E4. The holoenzyme is recognized as a vital antiviral target because MPXV replicates in the cytoplasm of host cells. Nucleotide analogs such as cidofovir and cytarabine (Ara-C) have shown potential in curbing MPXV replication and they also display promise against other poxviruses. However, the mechanism behind their inhibitory effects remains unclear. Here, we present the cryo-EM structure of the DNA polymerase holoenzyme F8/A22/E4 bound with its competitive inhibitor Ara-C-derived cytarabine triphosphate (Ara-CTP) at an overall resolution of 3.0 Å and reveal its inhibition mechanism. Ara-CTP functions as a direct chain terminator in proximity to the deoxycytidine triphosphate (dCTP)-binding site. The extra hydrogen bond formed with Asn665 makes it more potent in binding than dCTP. Asn665 is conserved among eukaryotic B-family polymerases.

Structures of human γδ T cell receptor-CD3 complex

Gamma delta (γδ) T cells, a unique T cell subgroup, are crucial in various immune responses and immunopathology1-3. The γδ T cell receptor (TCR), which is generated by γδ T cells, recognizes a diverse range of antigens independently of the major histocompatibility complex2. The γδ TCR associates with CD3 subunits, initiating T cell activation and holding great potential in immunotherapy4. Here we report the structures of two prototypical human Vγ9Vδ2 and Vγ5Vδ1 TCR-CD3 complexes5,6, revealing two distinct assembly mechanisms that depend on Vγ usage. The Vγ9Vδ2 TCR-CD3 complex is monomeric, with considerable conformational flexibility in the TCRγ-TCRδ extracellular domain and connecting peptides. The length of the connecting peptides regulates the ligand association and T cell activation. A cholesterol-like molecule wedges into the transmembrane region, exerting an inhibitory role in TCR signalling. The Vγ5Vδ1 TCR-CD3 complex displays a dimeric architecture, whereby two protomers nestle back to back through the Vγ5 domains of the TCR extracellular domains. Our biochemical and biophysical assays further corroborate the dimeric structure. Importantly, the dimeric form of the Vγ5Vδ1 TCR is essential for T cell activation. These findings reveal organizing principles of the γδ TCR-CD3 complex, providing insights into the unique properties of γδ TCR and facilitating immunotherapeutic interventions.

Single-particle Cryo-EM and molecular dynamics simulations: A perfect match

Knowledge of the structure and dynamics of biomolecules is key to understanding the mechanisms underlying their biological functions. Single-particle cryo-electron microscopy (cryo-EM) is a powerful structural biology technique to characterize complex biomolecular systems. Here, we review recent advances of how Molecular Dynamics (MD) simulations are being used to increase and enhance the information extracted from cryo-EM experiments. We will particularly focus on the physics underlying these experiments, how MD facilitates structure refinement, in particular for heterogeneous and non-isotropic resolution, and how thermodynamic and kinetic information can be extracted from cryo-EM data.

Vimentin filaments integrate low-complexity domains in a complex helical structure

Intermediate filaments (IFs) are integral components of the cytoskeleton. They provide cells with tissue-specific mechanical properties and are involved in numerous cellular processes. Due to their intricate architecture, a 3D structure of IFs has remained elusive. Here we use cryo-focused ion-beam milling, cryo-electron microscopy and tomography to obtain a 3D structure of vimentin IFs (VIFs). VIFs assemble into a modular, intertwined and flexible helical structure of 40 α-helices in cross-section, organized into five protofibrils. Surprisingly, the intrinsically disordered head domains form a fiber in the lumen of VIFs, while the intrinsically disordered tails form lateral connections between the protofibrils. Our findings demonstrate how protein domains of low sequence complexity can complement well-folded protein domains to construct a biopolymer with striking mechanical strength and stretchability.

Chromatin image-driven modelling

The challenge of modelling the spatial conformation of chromatin remains an open problem. While multiple data-driven approaches have been proposed, each has limitations. This work introduces two image-driven modelling methods based on the Molecular Dynamics Flexible Fitting (MDFF) approach: the force method and the correlational method. Both methods have already been used successfully in protein modelling. We propose a novel way to employ them for building chromatin models directly from 3D images. This approach is termed image-driven modelling. Additionally, we introduce the initial structure generator, a tool designed to generate optimal starting structures for the proposed algorithms. The methods are versatile and can be applied to various data types, with minor modifications to accommodate new generation imaging techniques.

A potent Henipavirus cross-neutralizing antibody reveals a dynamic fusion-triggering pattern of the G-tetramer

The Hendra and Nipah viruses (HNVs) are highly pathogenic pathogens without approved interventions for human use. In addition, the interaction pattern between the attachment (G) and fusion (F) glycoproteins required for virus entry remains unclear. Here, we isolate a panel of Macaca-derived G-specific antibodies that cross-neutralize HNVs via multiple mechanisms. The most potent antibody, 1E5, confers adequate protection against the Nipah virus challenge in female hamsters. Crystallography demonstrates that 1E5 has a highly similar binding pattern to the receptor. In cryo-electron microscopy studies, the tendency of 1E5 to bind to the upper or lower heads results in two distinct quaternary structures of G. Furthermore, we identify the extended outer loop β1S2-β1S3 of G and two pockets on the apical region of fusion (F) glycoprotein as the essential sites for G-F interactions. This work highlights promising drug candidates against HNVs and contributes deeper insights into the viruses.

Conformational changes in the Niemann-Pick type C1 protein NCR1 drive sterol translocation

The membrane protein Niemann-Pick type C1 (NPC1, named NCR1 in yeast) is central to sterol homeostasis in eukaryotes. Saccharomyces cerevisiae NCR1 is localized to the vacuolar membrane, where it is suggested to carry sterols across the protective glycocalyx and deposit them into the vacuolar membrane. However, documentation of a vacuolar glycocalyx in fungi is lacking, and the mechanism for sterol translocation has remained unclear. Here, we provide evidence supporting the presence of a glycocalyx in isolated S. cerevisiae vacuoles and report four cryo-EM structures of NCR1 in two distinct conformations, named tense and relaxed. These two conformations illustrate the movement of sterols through a tunnel formed by the luminal domains, thus bypassing the barrier presented by the glycocalyx. Based on these structures and on comparison with other members of the Resistance-Nodulation-Division (RND) superfamily, we propose a transport model that links changes in the luminal domains with a cycle of protonation and deprotonation within the transmembrane region of the protein. Our model suggests that NPC proteins work by a generalized RND mechanism where the proton motive force drives conformational changes in the transmembrane domains that are allosterically coupled to luminal/extracellular domains to promote sterol transport.

Cryo-EM structure of the human Asc-1 transporter complex

The Alanine-Serine-Cysteine transporter 1 (Asc-1 or SLC7A10) forms a crucial heterodimeric transporter complex with 4F2hc (SLC3A2) through a covalent disulfide bridge. This complex enables the sodium-independent transport of small neutral amino acids, including L-Alanine (L-Ala), Glycine (Gly), and D-Serine (D-Ser), within the central nervous system (CNS). D-Ser and Gly are two key endogenous glutamate co-agonists that activate N-methyl-d-aspartate (NMDA) receptors by binding to the allosteric site. Mice deficient in Asc-1 display severe symptoms such as tremors, ataxia, and seizures, leading to early postnatal death. Despite its physiological importance, the functional mechanism of the Asc-1-4F2hc complex has remained elusive. Here, we present cryo-electron microscopy (cryo-EM) structures of the human Asc-1-4F2hc complex in its apo state, D-Ser bound state, and L-Ala bound state, resolved at 3.6 Å, 3.5 Å, and 3.4 Å, respectively. Through detailed structural analysis and transport assays, we uncover a comprehensive alternating access mechanism that underlies conformational changes in the complex. In summary, our findings reveal the architecture of the Asc-1 and 4F2hc complex and provide valuable insights into substrate recognition and the functional cycle of this essential transporter complex.

Potent human neutralizing antibodies against Nipah virus derived from two ancestral antibody heavy chains

Nipah virus (NiV) is a World Health Organization priority pathogen and there are currently no approved drugs for clinical immunotherapy. Through the use of a naïve human phage-displayed Fab library, two neutralizing antibodies (NiV41 and NiV42) targeting the NiV receptor binding protein (RBP) were identified. Following affinity maturation, antibodies derived from NiV41 display cross-reactivity against both NiV and Hendra virus (HeV), whereas the antibody based on NiV42 is only specific to NiV. Results of immunogenetic analysis reveal a correlation between the maturation of antibodies and their antiviral activity. In vivo testing of NiV41 and its mature form (41-6) show protective efficacy against a lethal NiV challenge in hamsters. Furthermore, a 2.88 Å Cryo-EM structure of the tetrameric RBP and antibody complex demonstrates that 41-6 blocks the receptor binding interface. These findings can be beneficial for the development of antiviral drugs and the design of vaccines with broad spectrum against henipaviruses.

Differences between bacteria and eukaryotes in clamp loader mechanism, a conserved process underlying DNA replication

Clamp loaders are pentameric ATPases that place circular sliding clamps onto DNA, where they function in DNA replication and genome integrity. The central activity of a clamp loader is the opening of the ring-shaped sliding clamp and the subsequent binding to primer-template (p/t)-junctions. The general architecture of clamp loaders is conserved across all life, suggesting that their mechanism is retained. Recent structural studies of the eukaryotic clamp loader replication factor C (RFC) revealed that it functions using a crab-claw mechanism, where clamp opening is coupled to a massive conformational change in the loader. Here we investigate the clamp loading mechanism of the Escherichia coli clamp loader at high resolution using cryo-electron microscopy. We find that the E. coli clamp loader opens the clamp using a crab-claw motion at a single pivot point, whereas the eukaryotic RFC loader uses motions distributed across the complex. Furthermore, we find clamp opening occurs in multiple steps, starting with a partly open state with a spiral conformation, and proceeding to a wide open clamp in a surprising planar geometry. Finally, our structures in the presence of p/t-junctions illustrate how the clamp closes around p/t-junctions and how the clamp loader initiates release from the loaded clamp. Our results reveal mechanistic distinctions in a macromolecular machine that is conserved across all domains of life.

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semi-permeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo, the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We use molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo. Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceum model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

Structural insights into the C-terminus of the histone-lysine N-methyltransferase NSD3 by small-angle X-ray scattering

NSD3 is a member of six H3K36-specific histone lysine methyltransferases in metazoans. Its overexpression or mutation is implicated in developmental defects and oncogenesis. Aside from the well-characterized catalytic SET domain, NSD3 has multiple clinically relevant potential chromatin-binding motifs, such as the proline-tryptophan-tryptophan-proline (PWWP), the plant homeodomain (PHD), and the adjacent Cys-His-rich domain located at the C-terminus. The crystal structure of the individual domains is available, and this structural knowledge has allowed the designing of potential inhibitors, but the intrinsic flexibility of larger constructs has hindered the characterization of mutual domain conformations. Here, we report the first structural characterization of the NSD3 C-terminal region comprising the PWWP2, SET, and PHD4 domains, which has been achieved at a low resolution in solution by small-angle X-ray scattering (SAXS) data on two multiple-domain NSD3 constructs complemented with size-exclusion chromatography and advanced computational modeling. Structural models predicted by machine learning have been validated in direct space, by comparison with the SAXS-derived molecular envelope, and in reciprocal space, by reproducing the experimental SAXS profile. Selected models have been refined by SAXS-restrained molecular dynamics. This study shows how SAXS data can be used with advanced computational modeling techniques to achieve a detailed structural characterization and sheds light on how NSD3 domains are interconnected in the C-terminus.

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.

Structural basis for regulated assembly of the mitochondrial fission GTPase Drp1

Mitochondrial fission is a critical cellular event to maintain organelle function. This multistep process is initiated by the enhanced recruitment and oligomerization of dynamin-related protein 1 (Drp1) at the surface of mitochondria. As such, Drp1 is essential for inducing mitochondrial division in mammalian cells, and homologous proteins are found in all eukaryotes. As a member of the dynamin superfamily of proteins (DSPs), controlled Drp1 self-assembly into large helical polymers stimulates its GTPase activity to promote membrane constriction. Still, little is known about the mechanisms that regulate correct spatial and temporal assembly of the fission machinery. Here we present a cryo-EM structure of a full-length Drp1 dimer in an auto-inhibited state. This dimer reveals two key conformational rearrangements that must be unlocked through intramolecular rearrangements to achieve the assembly-competent state observed in previous structures. This structural insight provides understanding into the mechanism for regulated self-assembly of the mitochondrial fission machinery.

Understanding the cell: Future views of structural biology

Determining the structure and mechanisms of all individual functional modules of cells at high molecular detail has often been seen as equal to understanding how cells work. Recent technical advances have led to a flush of high-resolution structures of various macromolecular machines, but despite this wealth of detailed information, our understanding of cellular function remains incomplete. Here, we discuss present-day limitations of structural biology and highlight novel technologies that may enable us to analyze molecular functions directly inside cells. We predict that the progression toward structural cell biology will involve a shift toward conceptualizing a 4D virtual reality of cells using digital twins. These will capture cellular segments in a highly enriched molecular detail, include dynamic changes, and facilitate simulations of molecular processes, leading to novel and experimentally testable predictions. Transferring biological questions into algorithms that learn from the existing wealth of data and explore novel solutions may ultimately unveil how cells work.

DEMO-EM2: assembling protein complex structures from cryo-EM maps through intertwined chain and domain fitting

The breakthrough in cryo-electron microscopy (cryo-EM) technology has led to an increasing number of density maps of biological macromolecules. However, constructing accurate protein complex atomic structures from cryo-EM maps remains a challenge. In this study, we extend our previously developed DEMO-EM to present DEMO-EM2, an automated method for constructing protein complex models from cryo-EM maps through an iterative assembly procedure intertwining chain- and domain-level matching and fitting for predicted chain models. The method was carefully evaluated on 27 cryo-electron tomography (cryo-ET) maps and 16 single-particle EM maps, where DEMO-EM2 models achieved an average TM-score of 0.92, outperforming those of state-of-the-art methods. The results demonstrate an efficient method that enables the rapid and reliable solution of challenging cryo-EM structure modeling problems.

Observation of Bacteriophage Ultrastructure by Cryo-Electron Microscopy

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

Calcium-gated potassium channel blockade via membrane-facing fenestrations

Quaternary ammonium blockers were previously shown to bind in the pore to block both open and closed conformations of large-conductance calcium-activated potassium (BK and MthK) channels. Because blocker entry was assumed through the intracellular entryway (bundle crossing), closed-pore access suggested that the gate was not at the bundle crossing. Structures of closed MthK, a Methanobacterium thermoautotrophicum homolog of BK channels, revealed a tightly constricted intracellular gate, leading us to investigate the membrane-facing fenestrations as alternative pathways for blocker access directly from the membrane. Atomistic free energy simulations showed that intracellular blockers indeed access the pore through the fenestrations, and a mutant channel with narrower fenestrations displayed no closed-state TPeA block at concentrations that blocked the wild-type channel. Apo BK channels display similar fenestrations, suggesting that blockers may use them as access paths into closed channels. Thus, membrane fenestrations represent a non-canonical pathway for selective targeting of specific channel conformations, opening novel ways to selectively drug BK channels.

Differences in clamp loader mechanism between bacteria and eukaryotes

Clamp loaders are pentameric ATPases that place circular sliding clamps onto DNA, where they function in DNA replication and genome integrity. The central activity of a clamp loader is the opening of the ring-shaped sliding clamp, and the subsequent binding to primer-template (p/t)-junctions. The general architecture of clamp loaders is conserved across all life, suggesting that their mechanism is retained. Recent structural studies of the eukaryotic clamp loader Replication Factor C (RFC) revealed that it functions using a crab-claw mechanism, where clamp opening is coupled to a massive conformational change in the loader. Here we investigate the clamp loading mechanism of the E. coli clamp loader at high resolution using cryo-electron microscopy (cryo-EM). We find that the E. coli clamp loader opens the clamp using a crab-claw motion at a single pivot point, whereas the eukaryotic RFC loader uses motions distributed across the complex. Furthermore, we find clamp opening occurs in multiple steps, starting with a partly open state with a spiral conformation, and proceeding to a wide open clamp in a surprising planar geometry. Finally, our structures in the presence of p/t-junctions illustrate how clamp closes around p/t-junctions and how the clamp loader initiates release from the loaded clamp. Our results reveal mechanistic distinctions in a macromolecular machine that is conserved across all domains of life.

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

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