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

Structure of ATP synthase under strain during catalysis

ATP synthases are macromolecular machines consisting of an ATP-hydrolysis-driven F₁ motor and a proton-translocation-driven F_O motor. The F₁ and F_O motors oppose each other’s action on a shared rotor subcomplex and are held stationary relative to each other by a peripheral stalk. Structures of resting mitochondrial ATP synthases revealed a left-handed curvature of the peripheral stalk even though rotation of the rotor, driven by either ATP hydrolysis in F₁ or proton translocation through F_O, would apply a right-handed bending force to the stalk. We used cryoEM to image yeast mitochondrial ATP synthase under strain during ATP-hydrolysis-driven rotary catalysis, revealing a large deformation of the peripheral stalk. The structures show how the peripheral stalk opposes the bending force and suggests that during ATP synthesis proton translocation causes accumulation of strain in the stalk, which relaxes by driving the relative rotation of the rotor through six sub-steps within F₁, leading to catalysis.

CryoEM of mitochondrial ATP synthase frozen during rotary catalysis reveals dramatic conformational changes in the peripheral stalk subcomplex, which enable the enzyme’s efficient synthesis of ATP.

Novel sarbecovirus bispecific neutralizing antibodies with exceptional breadth and potency against currently circulating SARS-CoV-2 variants and sarbecoviruses

The emergence of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) variants of concern, including Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529) has aroused concerns over their increased infectivity and transmissibility, as well as decreased sensitivity to SARS-CoV-2-neutralizing antibodies (NAbs) and the current coronavirus disease 2019 (COVID-19) vaccines. Such exigencies call for the development of pan-sarbecovirus vaccines or inhibitors to combat the circulating SARS-CoV-2 NAb-escape variants and other sarbecoviruses. In this study, we isolated a broadly NAb against sarbecoviruses named GW01 from a donor who recovered from COVID-19. Cryo-EM structure and competition assay revealed that GW01 targets a highly conserved epitope in a wide spectrum of different sarbecoviruses. However, we found that GW01, the well-known sarbecovirus NAb S309, and the potent SARS-CoV-2 NAbs CC12.1 and REGN10989 only neutralize about 90% of the 56 tested currently circulating variants of SARS-CoV-2 including Omicron. Therefore, to improve efficacy, we engineered an IgG-like bispecific antibody GW01-REGN10989 (G9) consisting of single-chain antibody fragments (scFv) of GW01 and REGN10989. We found that G9 could neutralize 100% of NAb-escape mutants (23 out of 23), including Omicron variant, with a geometric mean (GM) 50% inhibitory concentration of 8.8 ng/mL. G9 showed prophylactic and therapeutic effects against SARS-CoV-2 infection of both the lung and brain in hACE2-transgenic mice. Site-directed mutagenesis analyses revealed that GW01 and REGN10989 bind to the receptor-binding domain in different epitopes and from different directions. Since G9 targets the epitopes for both GW01 and REGN10989, it was effective against variants with resistance to GW01 or REGN10989 alone and other NAb-escape variants. Therefore, this novel bispecific antibody, G9, is a strong candidate for the treatment and prevention of infection by SARS-CoV-2, NAb-escape variants, and other sarbecoviruses that may cause future emerging or re-emerging coronavirus diseases.

Building Biological Relevance Into Integrative Modelling of Macromolecular Assemblies

Recent advances in structural biophysics and integrative modelling methods now allow us to decipher the structures of large macromolecular assemblies. Understanding the dynamics and mechanisms involved in their biological function requires rigorous integration of all available data. We have developed a complete modelling pipeline that includes analyses to extract biologically significant information by consistently combining automated and interactive human-guided steps. We illustrate this idea with two examples. First, we describe the ryanodine receptor, an ion channel that controls ion flux across the cell membrane through transitions between open and closed states. The conformational changes associated with the transitions are small compared to the considerable system size of the receptor; it is challenging to consistently track these states with the available cryo-EM structures. The second example involves homologous recombination, in which long filaments of a recombinase protein and DNA catalyse the exchange of homologous DNA strands to reliably repair DNA double-strand breaks. The nucleoprotein filament reaction intermediates in this process are short-lived and heterogeneous, making their structures particularly elusive. The pipeline we describe, which incorporates experimental and theoretical knowledge combined with state-of-the-art interactive and immersive modelling tools, can help overcome these challenges. In both examples, we point to new insights into biological processes that arise from such interdisciplinary approaches.

Structural Characterization of the Full-Length Anti-CD20 Antibody Rituximab

Rituximab, a murine–human chimera, is the first monoclonal antibody (mAb) developed as a therapeutic agent to target CD20 protein. Its Fab domain and its interaction with CD20 have been extensively studied and high-resolution atomic models obtained by X-ray diffraction or cryo-electron microscopy are available. However, the structure of the full-length antibody is still missing as the inherent protein flexibility hampers the formation of well-diffracting crystals and the reconstruction of 3D microscope images. The global structure of rituximab from its dilute solution is here elucidated by small-angle X-ray scattering (SAXS). The limited data resolution achievable by this technique has been compensated by intensive computational modelling that led to develop a new and effective procedure to characterize the average mAb conformation as well as that of the single domains. SAXS data indicated that rituximab adopts an asymmetric average conformation in solution, with a radius of gyration and a maximum linear dimension of 52 Å and 197 Å, respectively. The asymmetry is mainly due to an uneven arrangement of the two Fab units with respect to the central stem (the Fc domain) and reflects in a different conformation of the individual units. As a result, the Fab elbow angle, which is a crucial determinant for antigen recognition and binding, was found to be larger (169°) in the more distant Fab unit than that in the less distant one (143°). The whole flexibility of the antibody has been found to strongly depend on the relative inter-domain orientations, with one of the Fab arms playing a major role. The average structure and the amount of flexibility has been studied in the presence of different buffers and additives, and monitored at increasing temperature, up to the complete unfolding of the antibody. Overall, the structural characterization of rituximab can help in designing next-generation anti-CD20 antibodies and finding more efficient routes for rituximab production at industrial level.

Progressive assembly of multi-domain protein structures from cryo-EM density maps

Progress in cryo-electron microscopy has provided the potential for large-size protein structure determination. However, the success rate for solving multi-domain proteins remains low because of the difficulty in modelling inter-domain orientations. Here we developed domain enhanced modeling using cryo-electron microscopy (DEMO-EM), an automatic method to assemble multi-domain structures from cryo-electron microscopy maps through a progressive structural refinement procedure combining rigid-body domain fitting and flexible assembly simulations with deep-neural-network inter-domain distance profiles. The method was tested on a large-scale benchmark set of proteins containing up to 12 continuous and discontinuous domains with medium- to low-resolution density maps, where DEMO-EM produced models with correct inter-domain orientations (template modeling score (TM-score) >0.5) for 97% of cases and outperformed state-of-the-art methods. DEMO-EM was applied to the severe acute respiratory syndrome coronavirus 2 genome and generated models with average TM-score and root-mean-square deviation of 0.97 and 1.3 Å, respectively, with respect to the deposited structures. These results demonstrate an efficient pipeline that enables automated and reliable large-scale multi-domain protein structure modelling from cryo-electron microscopy maps.

Bacterial ribosome collision sensing by a MutS DNA repair ATPase paralogue

Ribosome stalling during translation is detrimental to cellular fitness, but how this is sensed and elicits recycling of ribosomal subunits and quality control of associated mRNA and incomplete nascent chains is poorly understood^1,2. Here we uncover Bacillus subtilis MutS2, a member of the conserved MutS family of ATPases that function in DNA mismatch repair³, as an unexpected ribosome-binding protein with an essential function in translational quality control. Cryo-electron microscopy analysis of affinity-purified native complexes shows that MutS2 functions in sensing collisions between stalled and translating ribosomes and suggests how ribosome collisions can serve as platforms to deploy downstream processes: MutS2 has an RNA endonuclease small MutS-related (SMR) domain, as well as an ATPase/clamp domain that is properly positioned to promote ribosomal subunit dissociation, which is a requirement both for ribosome recycling and for initiation of ribosome-associated protein quality control (RQC). Accordingly, MutS2 promotes nascent chain modification with alanine-tail degrons-an early step in RQC-in an ATPase domain-dependent manner. The relevance of these observations is underscored by evidence of strong co-occurrence of MutS2 and RQC genes across bacterial phyla. Overall, the findings demonstrate a deeply conserved role for ribosome collisions in mounting a complex response to the interruption of translation within open reading frames.

Energy landscape of the SARS-CoV-2 reveals extensive conformational heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As with any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the glycosylated spike protein, and identify the diversity of low-energy conformations in the vicinity of its open (so called 1RBD-up) state. The resulting atomic refinement reveal global and local molecular rearrangements that cannot be inferred from an average 1RBD-up cryo-EM model. Here we report varied degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these low-energy conformations. Five out of six binding sites we analyzed, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This apo-to-holo switching is reminiscent of a conformational preequilibrium. We found only one binding site, namely that of AB-C135 remains in apo state within all the sampled free energy-minimizing models, suggesting an induced fit mechanism for the docking of this antibody to the spike.

Insights into POT1 structural dynamics revealed by cryo-EM

Telomeres are protein-DNA complexes that protect the ends of linear eukaryotic chromosomes. Mammalian telomeric DNA consists of 5′-(TTAGGG)n-3′ double-stranded repeats, followed by up to several hundred bases of a 3′ single-stranded G-rich overhang. The G-rich overhang is bound by the shelterin component POT1 which interacts with TPP1, the component involved in telomerase recruitment. A previously published crystal structure of the POT1 N-terminal half bound to the high affinity telomeric ligand 5′-TTAGGGTTAG-3′ showed that the first six nucleotides, TTAGGG, are bound by the OB1 fold, while the adjacent OB2 binds the last four, TTAG. Here, we report two cryo-EM structures of full-length POT1 bound by the POT1-binding domain of TPP1. The structures differ in the relative orientation of the POT1 OB1 and OB2, suggesting that these two DNA-binding OB folds take up alternative conformations. Supporting DNA binding studies using telomeric ligands in which the OB1 and OB2 binding sites were spaced apart, show that POT1 binds with similar affinities to spaced or contiguous binding sites, suggesting plasticity in DNA binding and a role for the alternative conformations observed. A likely explanation is that the structural flexibility of POT1 enhances binding to the tandemly arranged telomeric repeats and hence increases telomere protection.

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

Cryo-EM, Protein Engineering, and Simulation Enable the Development of Peptide Therapeutics against Acute Myeloid Leukemia

Cryogenic electron microscopy (cryo-EM) has emerged as a viable structural tool for molecular therapeutics development against human diseases. However, it remains a challenge to determine structures of proteins that are flexible and smaller than 30 kDa. The 11 kDa KIX domain of CREB-binding protein (CBP), a potential therapeutic target for acute myeloid leukemia and other cancers, is a protein which has defied structure-based inhibitor design. Here, we develop an experimental approach to overcome the size limitation by engineering a protein double-shell to sandwich the KIX domain between apoferritin as the inner shell and maltose-binding protein as the outer shell. To assist homogeneous orientations of the target, disulfide bonds are introduced at the target-apoferritin interface, resulting in a cryo-EM structure at 2.6 Å resolution. We used molecular dynamics simulations to design peptides that block the interaction of the KIX domain of CBP with the intrinsically disordered pKID domain of CREB. The double-shell design allows for fluorescence polarization assays confirming the binding between the KIX domain in the double-shell and these interacting peptides. Further cryo-EM analysis reveals a helix-helix interaction between a single KIX helix and the best peptide, providing a possible strategy for developments of next-generation inhibitors.

Broad ultra-potent neutralization of SARS-CoV-2 variants by monoclonal antibodies specific to the tip of RBD

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) continue to wreak havoc across the globe. Higher transmissibility and immunologic resistance of VOCs bring unprecedented challenges to epidemic extinguishment. Here we describe a monoclonal antibody, 2G1, that neutralizes all current VOCs and has surprising tolerance to mutations adjacent to or within its interaction epitope. Cryo-electron microscopy structure showed that 2G1 bound to the tip of receptor binding domain (RBD) of spike protein with small contact interface but strong hydrophobic effect, which resulted in nanomolar to sub-nanomolar affinities to spike proteins. The epitope of 2G1 on RBD partially overlaps with angiotensin converting enzyme 2 (ACE2) interface, which enables 2G1 to block interaction between RBD and ACE2. The narrow binding epitope but high affinity bestow outstanding therapeutic efficacy upon 2G1 that neutralized VOCs with sub-nanomolar half maximal inhibitory concentration in vitro. In SARS-CoV-2, Beta or Delta variant-challenged transgenic mice and rhesus macaque models, 2G1 protected animals from clinical illness and eliminated viral burden, without serious impact to animal safety. Mutagenesis experiments suggest that 2G1 is potentially capable of dealing with emerging SARS-CoV-2 variants in the future. This report characterized the therapeutic antibodies specific to the tip of spike against SARS-CoV-2 variants and highlights the potential clinical applications as well as for developing vaccine and cocktail therapy.

Interface refinement of low- to medium-resolution Cryo-EM complexes using HADDOCK2.4

A wide range of cellular processes requires the formation of multimeric protein complexes. The rise of cryo-electron microscopy (cryo-EM) has enabled the structural characterization of these protein assemblies. The density maps produced can, however, still suffer from limited resolution, impeding the process of resolving structures at atomic resolution. In order to solve this issue, monomers can be fitted into low- to medium-resolution maps. Unfortunately, the models produced frequently contain atomic clashes at the protein-protein interfaces (PPIs), as intermolecular interactions are typically not considered during monomer fitting. Here, we present a refinement approach based on HADDOCK2.4 to remove intermolecular clashes and optimize PPIs. A dataset of 14 cryo-EM complexes was used to test eight protocols. The best-performing protocol, consisting of a semi-flexible simulated annealing refinement with centroid restraints on the monomers, was able to decrease intermolecular atomic clashes by 98% without significantly deteriorating the quality of the cryo-EM density fit.

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

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

Human Learning for Molecular Simulations: The Collective Variables Dashboard in VMD

The Collective Variables Dashboard is a software tool for real-time, seamless exploration of molecular structures and trajectories in a customizable space of collective variables. The Dashboard arises from the integration of the Collective Variables Module (also known as Colvars) with the visualization software VMD, augmented with a fully discoverable graphical interface offering interactive workflows for the design and analysis of collective variables. Typical use cases include a priori design of collective variables for enhanced sampling and free energy simulations as well as analysis of any type of simulation or collection of structures in a collective variable space. A combination of those cases commonly occurs when preliminary simulations, biased or unbiased, reveal that an optimized set of collective variables is necessary to improve sampling in further simulations. Then the Dashboard provides an efficient way to intuitively explore the space of likely collective variables, validate them on existing data, and use the resulting collective variable definitions directly in further biased simulations using the Collective Variables Module. Visualization of biasing energies and forces is proposed to help analyze or plan biased simulations. We illustrate the use of the Dashboard on two applications: discovering coordinates to describe ligand unbinding from a protein binding site and designing volume-based variables to bias the hydration of a transmembrane pore.

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

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

Glycosylation and Serological Reactivity of an Expression-enhanced SARS-CoV-2 Viral Spike Mimetic

Extensive glycosylation of viral glycoproteins is a key feature of the antigenic surface of viruses and yet glycan processing can also be influenced by the manner of their recombinant production. The low yields of the soluble form of the trimeric spike (S) glycoprotein from SARS-CoV-2 has prompted advances in protein engineering that have greatly enhanced the stability and yields of the glycoprotein. The latest expression-enhanced version of the spike incorporates six proline substitutions to stabilize the prefusion conformation (termed SARS-CoV-2 S HexaPro). Although the substitutions greatly enhanced expression whilst not compromising protein structure, the influence of these substitutions on glycan processing has not been explored. Here, we show that the site-specific N-linked glycosylation of the expression-enhanced HexaPro resembles that of an earlier version containing two proline substitutions (2P), and that both capture features of native viral glycosylation. However, there are site-specific differences in glycosylation of HexaPro when compared to 2P. Despite these discrepancies, analysis of the serological reactivity of clinical samples from infected individuals confirmed that both HexaPro and 2P protein are equally able to detect IgG, IgA, and IgM responses in all sera analysed. Moreover, we extend this observation to include an analysis of glycan engineered S protein, whereby all N-linked glycans were converted to oligomannose-type and conclude that serological activity is not impacted by large scale changes in glycosylation. These observations suggest that variations in glycan processing will not impact the serological assessments currently being performed across the globe.

Modular assembly of the principal microtubule nucleator γ-TuRC

The gamma-tubulin ring complex (γ-TuRC) is the principal microtubule nucleation template in vertebrates. Recent cryo-EM reconstructions visualized the intricate quaternary structure of the γ-TuRC, containing more than thirty subunits, raising fundamental questions about γ-TuRC assembly and the role of actin as an integral part of the complex. Here, we reveal the structural mechanism underlying modular γ-TuRC assembly and identify a functional role of actin in microtubule nucleation. During γ-TuRC assembly, a GCP6-stabilized core comprising GCP2-3-4-5-4-6 is expanded by stepwise recruitment, selective stabilization and conformational locking of four pre-formed GCP2-GCP3 units. Formation of the lumenal bridge specifies incorporation of the terminal GCP2-GCP3 unit and thereby leads to closure of the γ-TuRC ring in a left-handed spiral configuration. Actin incorporation into the complex is not relevant for γ-TuRC assembly and structural integrity, but determines γ-TuRC geometry and is required for efficient microtubule nucleation and mitotic chromosome alignment in vivo.

Comprehensive structure and functional adaptations of the yeast nuclear pore complex

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

Using cryo-EM to uncover mechanisms of bacterial transcriptional regulation

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

Computational Modeling of the Virucidal Inhibition Mechanism for Broad-Spectrum Antiviral Nanoparticles and HPV16 Capsid Segments

Solid core nanoparticles (NPs) coated with sulfonated ligands that mimic heparan sulfate proteoglycans (HSPGs) can exhibit virucidal activity against many viruses that utilize HSPG interactions with host cells for the initial stages of infection. How the interactions of these NPs with large capsid segments of HSPG-interacting viruses lead to their virucidal activity has been unclear. Here, we describe the interactions between sulfonated NPs and segments of the human papilloma virus type 16 (HPV16) capsids using atomistic molecular dynamics simulations. The simulations demonstrate that the NPs primarily bind at the interfaces of two HPV16 capsid proteins. After equilibration, the distances and angles between capsid proteins in the capsid segments are larger for the systems in which the NPs bind at the interfaces of capsid proteins. Over time, NP binding can lead to breaking of contacts between two neighboring proteins. The revealed mechanism of NPs targeting the interfaces between pairs of capsid proteins can be utilized for designing new generations of virucidal materials and contribute to the development of new broad-spectrum non-toxic virucidal materials.

Plasticity within the barrel domain of BamA mediates a hybrid-barrel mechanism by BAM

In Gram-negative bacteria, the biogenesis of β-barrel outer membrane proteins is mediated by the β-barrel assembly machinery (BAM). The mechanism employed by BAM is complex and so far- incompletely understood. Here, we report the structures of BAM in nanodiscs, prepared using polar lipids and native membranes, where we observe an outward-open state. Mutations in the barrel domain of BamA reveal that plasticity in BAM is essential, particularly along the lateral seam of the barrel domain, which is further supported by molecular dynamics simulations that show conformational dynamics in BAM are modulated by the accessory proteins. We also report the structure of BAM in complex with EspP, which reveals an early folding intermediate where EspP threads from the underside of BAM and incorporates into the barrel domain of BamA, supporting a hybrid-barrel budding mechanism in which the substrate is folded into the membrane sequentially rather than as a single unit.

ChAdOx1 interacts with CAR and PF4 with implications for thrombosis with thrombocytopenia syndrome

Vaccines derived from chimpanzee adenovirus Y25 (ChAdOx1), human adenovirus type 26 (HAdV-D26), and human adenovirus type 5 (HAdV-C5) are critical in combatting the severe acute respiratory coronavirus 2 (SARS-CoV-2) pandemic. As part of the largest vaccination campaign in history, ultrarare side effects not seen in phase 3 trials, including thrombosis with thrombocytopenia syndrome (TTS), a rare condition resembling heparin-induced thrombocytopenia (HIT), have been observed. This study demonstrates that all three adenoviruses deployed as vaccination vectors versus SARS-CoV-2 bind to platelet factor 4 (PF4), a protein implicated in the pathogenesis of HIT. We have determined the structure of the ChAdOx1 viral vector and used it in state-of-the-art computational simulations to demonstrate an electrostatic interaction mechanism with PF4, which was confirmed experimentally by surface plasmon resonance. These data confirm that PF4 is capable of forming stable complexes with clinically relevant adenoviruses, an important step in unraveling the mechanisms underlying TTS.

Integrative structural modeling of macromolecular complexes using Assembline

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

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles

In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H₂O)₁₂₈) and ≈6.8 s ((H₂O)₂₅₆) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.

Processive cleavage of substrate at individual proteolytic active sites of the Lon protease complex

The Lon protease is the prototype of a family of proteolytic machines with adenosine triphosphatase modules built into a substrate degradation chamber. Lon is known to degrade protein substrates in a processive fashion, cutting a protein chain processively into small peptides before commencing cleavages of another protein chain. Here, we present structural and biochemical evidence demonstrating that processive substrate degradation occurs at each of the six proteolytic active sites of Lon, which forms a deep groove that partially encloses the substrate polypeptide chain by accommodating only the unprimed residues and permits processive cleavage in the C-to-N direction. We identify a universally conserved acidic residue at the exit side of the binding groove indispensable for the proteolytic activity. This noncatalytic residue likely promotes processive proteolysis by carboxyl-carboxylate interactions with cleaved intermediates. Together, these results uncover a previously unrecognized mechanism for processive substrate degradation by the Lon protease.

Mapping the catalytic conformations of an assembly-line polyketide synthase module

[Figure: see text].

Redeployment of automated MrBUMP search-model identification for map fitting in cryo-EM

In crystallography, the phase problem can often be addressed by the careful preparation of molecular-replacement search models. This has led to the development of pipelines such as MrBUMP that can automatically identify homologous proteins from an input sequence and edit them to focus on the areas that are most conserved. Many of these approaches can be applied directly to cryo-EM to help discover, prepare and correctly place models (here called cryo-EM search models) into electrostatic potential maps. This can significantly reduce the amount of manual model building that is required for structure determination. Here, MrBUMP is repurposed to fit automatically obtained PDB-derived chains and domains into cryo-EM maps. MrBUMP was successfully able to identify and place cryo-EM search models across a range of resolutions. Methods such as map segmentation are also explored as potential routes to improved performance. Map segmentation was also found to improve the effectiveness of the pipeline for higher resolution (<8 Å) data sets.

Lessons Learned from Responsive Molecular Dynamics Studies of the COVID-19 Virus

Over the past 18 months, the need to perform atomic detail molecular dynamics simulations of the SARS-CoV-2 virion, its spike protein, and other structures related to the viral infection cycle has led biomedical researchers worldwide to urgently seek out all available biomolecular structure information, appropriate molecular modeling and simulation software, and the necessary computing resources to conduct their work. We describe our experiences from several COVID-19 research collaborations and the challenges they presented in terms of our molecular modeling software development and support efforts, our laboratory's local computing environment, and our scientists' use of non-traditional HPC hardware platforms such as public clouds for large scale parallel molecular dynamics simulations.

Multiple nanocages of a cyanophage small heat shock protein with icosahedral and octahedral symmetries

The structures of a cyanophage small heat shock protein (sHSP) were determined as octahedrons of 24-mers and 48-mers and as icosahedrons of 60-mers. An N-terminal deletion construct of an 18 kDa sHSP of Synechococcus sp. phage S-ShM2 crystallized as a 24-mer and its structure was determined at a resolution of 7 Å. The negative stain electron microscopy (EM) images showed that the full-length protein is a mixture of a major population of larger and a minor population of smaller cage-like particles. Their structures have been determined by electron cryomicroscopy 3D image reconstruction at a resolution of 8 Å. The larger particles are 60-mers with icosahedral symmetry and the smaller ones are 48-mers with octahedral symmetry. These structures are the first of the viral/phage origin and the 60-mer is the largest and the first icosahedral assembly to be reported for sHSPs.

Complete three-dimensional structures of the Lon protease translocating a protein substrate

Lon is an evolutionarily conserved proteolytic machine carrying out a wide spectrum of biological activities by degrading misfolded damaged proteins and specific cellular substrates. Lon contains a large N-terminal domain and forms a hexameric core of fused adenosine triphosphatase and protease domains. Here, we report two complete structures of Lon engaging a substrate, determined by cryo–electron microscopy to 2.4-angstrom resolution. These structures show a multilayered architecture featuring a tensegrity triangle complex, uniquely constructed by six long N-terminal helices. The interlocked helix triangle is assembled on the top of the hexameric core to spread a web of six globular substrate-binding domains. It serves as a multipurpose platform that controls the access of substrates to the AAA+ ring, provides a ruler-based mechanism for substrate selection, and acts as a pulley device to facilitate unfolding of the translocated substrate. This work provides a complete framework for understanding the structural mechanisms of Lon.

SMARCAD1 is an ATP-dependent histone octamer exchange factor with de novo nucleosome assembly activity

The adenosine 5′-triphosphate (ATP)–dependent chromatin remodeler SMARCAD1 acts on nucleosomes during DNA replication, repair, and transcription, but despite its implication in disease, information on its function and biochemical activities is scarce. Chromatin remodelers use the energy of ATP hydrolysis to slide nucleosomes, evict histones, or exchange histone variants. Here, we show that SMARCAD1 transfers the entire histone octamer from one DNA segment to another in an ATP-dependent manner but is also capable of de novo nucleosome assembly from histone octamer because of its ability to simultaneously bind all histones. We present a low-resolution cryo–electron microscopy structure of SMARCAD1 in complex with a nucleosome and show that the adenosine triphosphatase domains engage their substrate unlike any other chromatin remodeler. Our biochemical and structural data provide mechanistic insights into SMARCAD1-induced nucleosome disassembly and reassembly.

CryoFold: determining protein structures and data-guided ensembles from cryo-EM density maps

Cryo-electron microscopy (EM) requires molecular modeling to refine structural details from data. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins that cannot be resolved by cryo-EM. Here, we introduce CryoFold - a pipeline of molecular dynamics simulations that determines ensembles of protein structures directly from sequence by integrating density data of varying sparsity at 3-5 Å resolution with coarse-grained topological knowledge of the protein folds. We present six examples showing its broad applicability for folding proteins between 72 to 2000 residues, including large membrane and multi-domain systems, and results from two EMDB competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations, unseen by traditional methods, are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python GUI.

Molecular mechanism of sugar transport in plants unveiled by structures of glucose/H+ symporter STP10

Sugars are essential sources of energy and carbon and also function as key signalling molecules in plants. Sugar transport proteins (STP) are proton-coupled symporters responsible for uptake of glucose from the apoplast into plant cells. They are integral to organ development in symplastically isolated tissues such as seed, pollen and fruit. Additionally, STPs play a vital role in plant responses to stressors such as dehydration and prevalent fungal infections like rust and mildew. Here we present a structure of Arabidopsis thaliana STP10 in the inward-open conformation at 2.6 Å resolution and a structure of the outward-occluded conformation at improved 1.8 Å resolution, both with glucose and protons bound. The two structures describe key states in the STP transport cycle. Together with molecular dynamics simulations that establish protonation states and biochemical analysis, they pinpoint structural elements, conserved in all STPs, that clarify the basis of proton-to-glucose coupling. These results advance our understanding of monosaccharide uptake, which is essential for plant organ development, and set the stage for bioengineering strategies in crops.

Molecular basis for ATPase-powered substrate translocation by the Lon AAA+ protease

The Lon AAA+ (adenosine triphosphatases associated with diverse cellular activities) protease (LonA) converts ATP-fuelled conformational changes into sufficient mechanical force to drive translocation of a substrate into a hexameric proteolytic chamber. To understand the structural basis for the substrate translocation process, we determined the cryo-electron microscopy (cryo-EM) structure of Meiothermus taiwanensis LonA (MtaLonA) in a substrate-engaged state at 3.6 Å resolution. Our data indicate that substrate interactions are mediated by the dual pore loops of the ATPase domains, organized in spiral staircase arrangement from four consecutive protomers in different ATP-binding and hydrolysis states. However, a closed AAA+ ring is maintained by two disengaged ADP-bound protomers transiting between the lowest and highest position. This structure reveals a processive rotary translocation mechanism mediated by LonA-specific nucleotide-dependent allosteric coordination among the ATPase domains, which is induced by substrate binding.

Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump

The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun's lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.

Integration of Mass Spectrometry Data for Structural Biology

Mass spectrometry (MS) is increasingly being used to probe the structure and dynamics of proteins and the complexes they form with other macromolecules. There are now several specialized MS methods, each with unique sample preparation, data acquisition, and data processing protocols. Collectively, these methods are referred to as structural MS and include cross-linking, hydrogen-deuterium exchange, hydroxyl radical footprinting, native, ion mobility, and top-down MS. Each of these provides a unique type of structural information, ranging from composition and stoichiometry through to residue level proximity and solvent accessibility. Structural MS has proved particularly beneficial in studying protein classes for which analysis by classic structural biology techniques proves challenging such as glycosylated or intrinsically disordered proteins. To capture the structural details for a particular system, especially larger multiprotein complexes, more than one structural MS method with other structural and biophysical techniques is often required. Key to integrating these diverse data are computational strategies and software solutions to facilitate this process. We provide a background to the structural MS methods and briefly summarize other structural methods and how these are combined with MS. We then describe current state of the art approaches for the integration of structural MS data for structural biology. We quantify how often these methods are used together and provide examples where such combinations have been fruitful. To illustrate the power of integrative approaches, we discuss progress in solving the structures of the proteasome and the nuclear pore complex. We also discuss how information from structural MS, particularly pertaining to protein dynamics, is not currently utilized in integrative workflows and how such information can provide a more accurate picture of the systems studied. We conclude by discussing new developments in the MS and computational fields that will further enable in-cell structural studies.

Validation, analysis and annotation of cryo-EM structures

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

Reconstruction of Three-Dimensional Conformations of Bacterial ClpB from High-Speed Atomic-Force-Microscopy Images

ClpB belongs to the cellular disaggretase machinery involved in rescuing misfolded or aggregated proteins during heat or other cellular shocks. The function of this protein relies on the interconversion between different conformations in its native condition. A recent high-speed-atomic-force-microscopy (HS-AFM) experiment on ClpB from Thermus thermophilus shows four predominant conformational classes, namely, open, closed, spiral, and half-spiral. Analyses of AFM images provide only partial structural information regarding the molecular surface, and thus computational modeling of three-dimensional (3D) structures of these conformations should help interpret dynamical events related to ClpB functions. In this study, we reconstruct 3D models of ClpB from HS-AFM images in different conformational classes. We have applied our recently developed computational method based on a low-resolution representation of 3D structure using a Gaussian mixture model, combined with a Monte-Carlo sampling algorithm to optimize the agreement with target AFM images. After conformational sampling, we obtained models that reflect conformational variety embedded within the AFM images. From these reconstructed 3D models, we described, in terms of relative domain arrangement, the different types of ClpB oligomeric conformations observed by HS-AFM experiments. In particular, we highlighted the slippage of the monomeric components around the seam. This study demonstrates that such details of information, necessary for annotating the different conformational states involved in the ClpB function, can be obtained by combining HS-AFM images, even with limited resolution, and computational modeling.

Structures of tmRNA and SmpB as they transit through the ribosome

In bacteria, trans-translation is the main rescue system, freeing ribosomes stalled on defective messenger RNAs. This mechanism is driven by small protein B (SmpB) and transfer-messenger RNA (tmRNA), a hybrid RNA known to have both a tRNA-like and an mRNA-like domain. Here we present four cryo-EM structures of the ribosome during trans-translation at resolutions from 3.0 to 3.4 Å. These include the high-resolution structure of the whole pre-accommodated state, as well as structures of the accommodated state, the translocated state, and a translocation intermediate. Together, they shed light on the movements of the tmRNA-SmpB complex in the ribosome, from its delivery by the elongation factor EF-Tu to its passage through the ribosomal A and P sites after the opening of the B1 bridges. Additionally, we describe the interactions between the tmRNA-SmpB complex and the ribosome. These explain why the process does not interfere with canonical translation.

Large-Scale Molecular Dynamics Simulations of Cellular Compartments

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.

Poor Person's pH Simulation of Membrane Proteins

pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and V_o domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations.

Deactivation blocks proton pathways in the mitochondrial complex I

The electron transport chain of mitochondria is initiated by the respiratory complex I that converts chemical energy into a proton motive force to power synthesis of adenosine triphosphate. On a chemical level, complex I catalyzes elementary electron and proton transfer processes that couple across large molecular distances of >300 Å. However, under low oxygen concentrations, the respiratory chain operates in reverse mode and produces harmful reactive oxygen species. To avoid cell damage, the mitochondrial complex I transitions into a deactive state that inhibits turnover by molecular principles that remain elusive. By combining large-scale molecular simulations with cryo-electron microscopy data, we show here that complex I deactivation blocks the communication between proton pumping and redox modules by conformational and hydration changes.

Cellular respiration is powered by membrane-bound redox enzymes that convert chemical energy into an electrochemical proton gradient and drive the energy metabolism. By combining large-scale classical and quantum mechanical simulations with cryo-electron microscopy data, we resolve here molecular details of conformational changes linked to proton pumping in the mammalian complex I. Our data suggest that complex I deactivation blocks water-mediated proton transfer between a membrane-bound quinone site and proton-pumping modules, decoupling the energy-transduction machinery. We identify a putative gating region at the interface between membrane domain subunits ND1 and ND3/ND4L/ND6 that modulates the proton transfer by conformational changes in transmembrane helices and bulky residues. The region is perturbed by mutations linked to human mitochondrial disorders and is suggested to also undergo conformational changes during catalysis of simpler complex I variants that lack the “active”-to-“deactive” transition. Our findings suggest that conformational changes in transmembrane helices modulate the proton transfer dynamics by wetting/dewetting transitions and provide important functional insight into the mammalian respiratory complex I.

Zinc binding alters the conformational dynamics and drives the transport cycle of the cation diffusion facilitator YiiP

YiiP is a secondary transporter that couples Zn2+ transport to the proton motive force. Structural studies of YiiP from prokaryotes and Znt8 from humans have revealed three different Zn2+ sites and a conserved homodimeric architecture. These structures define the inward-facing and outward-facing states that characterize the archetypal alternating access mechanism of transport. To study the effects of Zn2+ binding on the conformational transition, we use cryo-EM together with molecular dynamics simulation to compare structures of YiiP from Shewanella oneidensis in the presence and absence of Zn2+. To enable single-particle cryo-EM, we used a phage-display library to develop a Fab antibody fragment with high affinity for YiiP, thus producing a YiiP/Fab complex. To perform MD simulations, we developed a nonbonded dummy model for Zn2+ and validated its performance with known Zn2+-binding proteins. Using these tools, we find that, in the presence of Zn2+, YiiP adopts an inward-facing conformation consistent with that previously seen in tubular crystals. After removal of Zn2+ with high-affinity chelators, YiiP exhibits enhanced flexibility and adopts a novel conformation that appears to be intermediate between inward-facing and outward-facing states. This conformation involves closure of a hydrophobic gate that has been postulated to control access to the primary transport site. Comparison of several independent cryo-EM maps suggests that the transition from the inward-facing state is controlled by occupancy of a secondary Zn2+ site at the cytoplasmic membrane interface. This work enhances our understanding of individual Zn2+ binding sites and their role in the conformational dynamics that govern the transport cycle.

The role of membrane destabilisation and protein dynamics in BAM catalysed OMP folding

The folding of β-barrel outer membrane proteins (OMPs) in Gram-negative bacteria is catalysed by the β-barrel assembly machinery (BAM). How lateral opening in the β-barrel of the major subunit BamA assists in OMP folding, and the contribution of membrane disruption to BAM catalysis remain unresolved. Here, we use an anti-BamA monoclonal antibody fragment (Fab1) and two disulphide-crosslinked BAM variants (lid-locked (LL), and POTRA-5-locked (P5L)) to dissect these roles. Despite being lethal in vivo, we show that all complexes catalyse folding in vitro, albeit less efficiently than wild-type BAM. CryoEM reveals that while Fab1 and BAM-P5L trap an open-barrel state, BAM-LL contains a mixture of closed and contorted, partially-open structures. Finally, all three complexes globally destabilise the lipid bilayer, while BamA does not, revealing that the BAM lipoproteins are required for this function. Together the results provide insights into the role of BAM structure and lipid dynamics in OMP folding.

Integrative Structural Biology in the Era of Accurate Structure Prediction

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

Structures of rhodopsin in complex with G-protein-coupled receptor kinase 1

G-protein-coupled receptor (GPCR) kinases (GRKs) selectively phosphorylate activated GPCRs, thereby priming them for desensitization¹. Although it is unclear how GRKs recognize these receptors^2-4, a conserved region at the GRK N terminus is essential for this process^5-8. Here we report a series of cryo-electron microscopy single-particle reconstructions of light-activated rhodopsin (Rho*) bound to rhodopsin kinase (GRK1), wherein the N terminus of GRK1 forms a helix that docks into the open cytoplasmic cleft of Rho*. The helix also packs against the GRK1 kinase domain and stabilizes it in an active configuration. The complex is further stabilized by electrostatic interactions between basic residues that are conserved in most GPCRs and acidic residues that are conserved in GRKs. We did not observe any density for the regulator of G-protein signalling homology domain of GRK1 or the C terminus of rhodopsin. Crosslinking with mass spectrometry analysis confirmed these results and revealed dynamic behaviour in receptor-bound GRK1 that would allow the phosphorylation of multiple sites in the receptor tail. We have identified GRK1 residues whose mutation augments kinase activity and crosslinking with Rho*, as well as residues that are involved in activation by acidic phospholipids. From these data, we present a general model for how a small family of protein kinases can recognize and be activated by hundreds of different GPCRs.

Structure and transport mechanism of P5B-ATPases

In human cells, P5B-ATPases execute the active export of physiologically important polyamines such as spermine from lysosomes to the cytosol, a function linked to a palette of disorders. Yet, the overall shape of P5B-ATPases and the mechanisms of polyamine recognition, uptake and transport remain elusive. Here we describe a series of cryo-electron microscopy structures of a yeast homolog of human ATP13A2-5, Ypk9, determined at resolutions reaching 3.4 Å, and depicting three separate transport cycle intermediates, including spermine-bound conformations. Surprisingly, in the absence of cargo, Ypk9 rests in a phosphorylated conformation auto-inhibited by the N-terminus. Spermine uptake is accomplished through an electronegative cleft lined by transmembrane segments 2, 4 and 6. Despite the dramatically different nature of the transported cargo, these findings pinpoint shared principles of transport and regulation among the evolutionary related P4-, P5A- and P5B-ATPases. The data also provide a framework for analysis of associated maladies, such as Parkinson’s disease.

In human cells, P5B‐ATPases execute export of spermine from lysosomes to the cytosol, but the mechanisms of spermine recognition, uptake and transport remain elusive. Here the authors present cryo‐EM structures of a yeast homolog of human ATP13A2‐5, Ypk9, which depict three separate transport cycle intermediates, including spermine‐bound conformations

Structural basis for VIPP1 oligomerization and maintenance of thylakoid membrane integrity

Vesicle-inducing protein in plastids 1 (VIPP1) is essential for the biogenesis and maintenance of thylakoid membranes, which transform light into life. However, it is unknown how VIPP1 performs its vital membrane-remodeling functions. Here, we use cryo-electron microscopy to determine structures of cyanobacterial VIPP1 rings, revealing how VIPP1 monomers flex and interweave to form basket-like assemblies of different symmetries. Three VIPP1 monomers together coordinate a non-canonical nucleotide binding pocket on one end of the ring. Inside the ring's lumen, amphipathic helices from each monomer align to form large hydrophobic columns, enabling VIPP1 to bind and curve membranes. In vivo mutations in these hydrophobic surfaces cause extreme thylakoid swelling under high light, indicating an essential role of VIPP1 lipid binding in resisting stress-induced damage. Using cryo-correlative light and electron microscopy (cryo-CLEM), we observe oligomeric VIPP1 coats encapsulating membrane tubules within the Chlamydomonas chloroplast. Our work provides a structural foundation for understanding how VIPP1 directs thylakoid biogenesis and maintenance.

Systemic cancer therapy with engineered adenovirus that evades innate immunity

Oncolytic virus therapy is a cancer treatment modality that has the potential to improve outcomes for patients with currently incurable malignancies. Although intravascular delivery of therapeutic viruses provides access to disseminated tumors, this delivery route exposes the virus to opsonizing and inactivating factors in the blood, which limit the effective therapeutic virus dose and contribute to activation of systemic toxicities. When human species C adenovirus HAdv-C5 is delivered intravenously, natural immunoglobulin M (IgM) antibodies and coagulation factor X rapidly opsonize HAdv-C5, leading to virus sequestration in tissue macrophages and promoting infection of liver cells, triggering hepatotoxicity. Here, we showed that natural IgM antibody binds to the hypervariable region 1 (HVR1) of the main HAdv-C5 capsid protein hexon. Using compound targeted mutagenesis of hexon HVR1 loop and other functional sites that mediate virus-host interactions, we engineered and obtained a high-resolution cryo-electron microscopy structure of an adenovirus vector, Ad5-3M, which resisted inactivation by blood factors, avoided sequestration in liver macrophages, and failed to trigger hepatotoxicity after intravenous delivery. Systemic delivery of Ad5-3M to mice with localized or disseminated lung cancer led to viral replication in tumor cells, suppression of tumor growth, and prolonged survival. Thus, compound targeted mutagenesis of functional sites in the virus capsid represents a generalizable approach to tailor virus interactions with the humoral and cellular arms of the immune system, enabling generation of "designer" viruses with improved therapeutic properties.

Cryo-Electron Microscopy and Mass Analysis of Oligolysine-Coated DNA Nanostructures

Cationic coatings can enhance the stability of synthetic DNA objects in low ionic strength environments such as physiological fluids. Here, we used single-particle cryo-electron microscopy (cryo-EM), pseudoatomic model fitting, and single-molecule mass photometry to study oligolysine and polyethylene glycol (PEG)-oligolysine-coated multilayer DNA origami objects. The coatings preserve coarse structural features well on a resolution of multiple nanometers but can also induce deformations such as twisting and bending. Higher-density coatings also led to internal structural deformations in the DNA origami test objects, in which a designed honeycomb-type helical lattice was deformed into a more square-lattice-like pattern. Under physiological ionic strength, where the uncoated objects disassembled, the coated objects remained intact but they shrunk in the helical direction and expanded in the direction perpendicular to the helical axis. Helical details like major/minor grooves and crossover locations were not discernible in cryo-EM maps that we determined of DNA origami coated with oligolysine and PEG-oligolysine, whereas these features were visible in cryo-EM maps determined from the uncoated reference objects. Blunt-ended double-helical interfaces remained accessible underneath the coating and may be used for the formation of multimeric DNA origami assemblies that rely on stacking interactions between blunt-ended helices. The ionic strength requirements for forming multimers from coated DNA origami differed from those needed for uncoated objects. Using single-molecule mass photometry, we found that the mass of coated DNA origami objects prior to and after incubation in low ionic strength physiological conditions remained unchanged. This finding indicated that the coating effectively prevented strand dissociation but also that the coating itself remained stable in place. Our results validate oligolysine coatings as a powerful stabilization method for DNA origami but also reveal several potential points of failure that experimenters should watch to avoid working with false premises.