Cryo_fit: Democratization of flexible fitting for cryo-EM

Beta-arrestin 1 mediated Src activation via Src SH3 domain revealed by cryo-electron microscopy

Beta-arrestins (βarrs) are key regulators and transducers of G-protein coupled receptor signaling; however, little is known of how βarrs communicate with their downstream effectors. Here, we use cryo-electron microscopy to elucidate how βarr1 recruits and activates non-receptor tyrosine kinase Src. βarr1 binds Src SH3 domain via two distinct sites: a polyproline site in the N-domain and a non-proline site in the central crest region. At both sites βarr1 interacts with the aromatic surface of SH3 which is critical for Src autoinhibition, suggesting that βarr1 activates Src by SH3 domain displacement. Binding of SH3 to the central crest region induces structural rearrangements in the β-strand V, finger, and middle loops of βarr1 and interferes with βarr1 coupling to the receptor core potentially impacting receptor desensitization and downstream signaling.

One-Sentence Summary

Beta-arrestin 1 uses two distinct sites to bind the SH3 domain of Src and drive relief of Src autoinhibition.

Dissemination of pathogenic bacteria is reinforced by a MARTX toxin effector duet

Multiple bacterial genera take advantage of the multifunctional autoprocessing repeats-in-toxin (MARTX) toxin to invade host cells. Secretion of the MARTX toxin by Vibrio vulnificus, a deadly opportunistic pathogen that causes primary septicemia, the precursor of sepsis, is a major driver of infection; however, the molecular mechanism via which the toxin contributes to septicemia remains unclear. Here, we report the crystal and cryo-electron microscopy (EM) structures of a toxin effector duet comprising the domain of unknown function in the first position (DUF1)/Rho inactivation domain (RID) complexed with human targets. These structures reveal how the duet is used by bacteria as a potent weapon. The data show that DUF1 acts as a RID-dependent transforming NADase domain (RDTND) that disrupts NAD+ homeostasis by hijacking calmodulin. The cryo-EM structure of the RDTND-RID duet complexed with calmodulin and Rac1, together with immunological analyses in vitro and in mice, provide mechanistic insight into how V. vulnificus uses the duet to suppress ROS generation by depleting NAD(P)+ and modifying Rac1 in a mutually-reinforcing manner that ultimately paralyzes first line immune responses, promotes dissemination of invaders, and induces sepsis. These data may allow development of tools or strategies to combat MARTX toxin-related human diseases.

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.

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

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

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

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

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge-driven or artificial intelligence-driven and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary-structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models that contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly while refining them against the density maps. We introduce such an adaptive decision-making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two prototypical protein systems, adenylate kinase and carbon monoxide dehydrogenase. For these test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2-3 Å map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that are inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from the chimpanzee adenovirus. The overhead for decision making remains low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install R-MDFF on different computing environments, from local Linux-based workstations to high-performance computing environments.

Gentle and fast all-atom model refinement to cryo-EM densities via a maximum likelihood approach

Better detectors and automated data collection have generated a flood of high-resolution cryo-EM maps, which in turn has renewed interest in improving methods for determining structure models corresponding to these maps. However, automatically fitting atoms to densities becomes difficult as their resolution increases and the refinement potential has a vast number of local minima. In practice, the problem becomes even more complex when one also wants to achieve a balance between a good fit of atom positions to the map, while also establishing good stereochemistry or allowing protein secondary structure to change during fitting. Here, we present a solution to this challenge using a maximum likelihood approach by formulating the problem as identifying the structure most likely to have produced the observed density map. This allows us to derive new types of smooth refinement potential-based on relative entropy-in combination with a novel adaptive force scaling algorithm to allow balancing of force-field and density-based potentials. In a low-noise scenario, as expected from modern cryo-EM data, the relative-entropy based refinement potential outperforms alternatives, and the adaptive force scaling appears to aid all existing refinement potentials. The method is available as a component in the GROMACS molecular simulation toolkit.

Release of frustration drives corneal amyloid disaggregation by brain chaperone

TGFBI-related corneal dystrophy (CD) is characterized by the accumulation of insoluble protein deposits in the corneal tissues, eventually leading to progressive corneal opacity. Here we show that ATP-independent amyloid-β chaperone L-PGDS can effectively disaggregate corneal amyloids in surgically excised human cornea of TGFBI-CD patients and release trapped amyloid hallmark proteins. Since the mechanism of amyloid disassembly by ATP-independent chaperones is unknown, we reconstructed atomic models of the amyloids self-assembled from TGFBIp-derived peptides and their complex with L-PGDS using cryo-EM and NMR. We show that L-PGDS specifically recognizes structurally frustrated regions in the amyloids and releases those frustrations. The released free energy increases the chaperone's binding affinity to amyloids, resulting in local restructuring and breakage of amyloids to protofibrils. Our mechanistic model provides insights into the alternative source of energy utilized by ATP-independent disaggregases and highlights the possibility of using these chaperones as treatment strategies for different types of amyloid-related diseases.

StarMap: a user-friendly workflow for Rosetta-driven molecular structure refinement

Cryogenic electron microscopy (cryo-EM) data represent density maps of macromolecular systems at atomic or near-atomic resolution. However, building and refining 3D atomic models by using data from cryo-EM maps is not straightforward and requires significant hands-on experience and manual intervention. We recently developed StarMap, an easy-to-use interface between the popular structural display program ChimeraX and Rosetta, a powerful molecular modeling engine. StarMap offers a general approach for refining structural models of biological macromolecules into cryo-EM density maps by combining Monte Carlo sampling with local density-guided optimization, Rosetta-based all-atom refinement and real-space B-factor calculations in a straightforward workflow. StarMap includes options for structural symmetry, local refinements and independent model validation. The overall quality of the refinement and the structure resolution is then assessed via analytical outputs, such as magnification calibration (pixel size calibration) and Fourier shell correlations. Z-scores reported by StarMap provide an easily interpretable indicator of the goodness of fit for each residue and can be plotted to evaluate structural models and improve local residue refinements, as well as to identify flexible regions and potentially functional sites in large macromolecular complexes. The protocol requires general computer skills, without the need for coding expertise, because most parts of the workflow can be operated by clicking tabs within the ChimeraX graphical user interface. Time requirements for the model refinement depend on the size and quality of the input data; however, this step can typically be completed within 1 d. The analytical parts of the workflow are completed within minutes.

Open Structural Data in Precision Medicine

Three-dimensional protein structural data at the molecular level are pivotal for successful precision medicine. Such data are crucial not only for discovering drugs that act to block the active site of the target mutant protein but also for clarifying to the patient and the clinician how the mutations harbored by the patient work. The relative paucity of structural data reflects their cost, challenges in their interpretation, and lack of clinical guidelines for their utilization. Rapid technological advancements in experimental high-resolution structural determination increasingly generate structures. Computationally, modeling algorithms, including molecular dynamics simulations, are becoming more powerful, as are compute-intensive hardware, particularly graphics processing units, overlapping with the inception of the exascale era. Accessible, freely available, and detailed structural and dynamical data can be merged with big data to powerfully transform personalized pharmacology. Here we review protein and emerging genome high-resolution data, along with means, applications, and examples underscoring their usefulness in precision medicine. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Unusual nucleosome formation and transcriptome influence by the histone H3mm18 variant

Histone H3mm18 is a non-allelic H3 variant expressed in skeletal muscle and brain in mice. However, its function has remained enigmatic. We found that H3mm18 is incorporated into chromatin in cells with low efficiency, as compared to H3.3. We determined the structures of the nucleosome core particle (NCP) containing H3mm18 by cryo-electron microscopy, which revealed that the entry/exit DNA regions are drastically disordered in the H3mm18 NCP. Consistently, the H3mm18 NCP is substantially unstable in vitro. The forced expression of H3mm18 in mouse myoblast C2C12 cells markedly suppressed muscle differentiation. A transcriptome analysis revealed that the forced expression of H3mm18 affected the expression of multiple genes, and suppressed a group of genes involved in muscle development. These results suggest a novel gene expression regulation system in which the chromatin landscape is altered by the formation of unusual nucleosomes with a histone variant, H3mm18, and provide important insight into understanding transcription regulation by chromatin.

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.

EMDA: A Python package for Electron Microscopy Data Analysis

An open-source Python library EMDA for cryo-EM map and model manipulation is presented with a specific focus on validation. The use of several functionalities in the library is presented through several examples. The utility of local correlation as a metric for identifying map-model differences and unmodeled regions in maps, and how it is used as a metric of map-model validation is demonstrated. The mapping of local correlation to individual atoms, and its use to draw insights on local signal variations are discussed. EMDA’s likelihood-based map overlay is demonstrated by carrying out a superposition of two domains in two related structures. The overlay is carried out first to bring both maps into the same coordinate frame and then to estimate the relative movement of domains. Finally, the map magnification refinement in EMDA is presented with an example to highlight the importance of adjusting the map magnification in structural comparison studies.

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.

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.

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

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

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

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

Structure of the far-red light utilizing photosystem I of Acaryochloris marina

Acaryochloris marina is one of the cyanobacterial species that can use far-red light to drive photochemical reactions for oxygenic photosynthesis. Here, we report the structure of A. marina photosystem I (PSI) reaction center, determined by cryo-electron microscopy at 2.58 Å resolution. The structure reveals an arrangement of electron carriers and light-harvesting pigments distinct from other type I reaction centers. The paired chlorophyll, or special pair (also referred to as P740 in this case), is a dimer of chlorophyll d and its epimer chlorophyll d'. The primary electron acceptor is pheophytin a, a metal-less chlorin. We show the architecture of this PSI reaction center is composed of 11 subunits and we identify key components that help explain how the low energy yield from far-red light is efficiently utilized for driving oxygenic photosynthesis.

Mitochondrial sorting and assembly machinery operates by β-barrel switching

The mitochondrial outer membrane contains so-called β-barrel proteins, which allow communication between the cytosol and the mitochondrial interior^1-3. Insertion of β-barrel proteins into the outer membrane is mediated by the multisubunit mitochondrial sorting and assembly machinery (SAM, also known as TOB)^4-6. Here we use cryo-electron microscopy to determine the structures of two different forms of the yeast SAM complex at a resolution of 2.8-3.2 Å. The dimeric complex contains two copies of the β-barrel channel protein Sam50-Sam50a and Sam50b-with partially open lateral gates. The peripheral membrane proteins Sam35 and Sam37 cap the Sam50 channels from the cytosolic side, and are crucial for the structural and functional integrity of the dimeric complex. In the second complex, Sam50b is replaced by the β-barrel protein Mdm10. In cooperation with Sam50a, Sam37 recruits and traps Mdm10 by penetrating the interior of its laterally closed β-barrel from the cytosolic side. The substrate-loaded SAM complex contains one each of Sam50, Sam35 and Sam37, but neither Mdm10 nor a second Sam50, suggesting that Mdm10 and Sam50b function as placeholders for a β-barrel substrate released from Sam50a. Our proposed mechanism for dynamic switching of β-barrel subunits and substrate explains how entire precursor proteins can fold in association with the mitochondrial machinery for β-barrel assembly.

Cryo-EM structure of native human uromodulin, a zona pellucida module polymer

Assembly of extracellular filaments and matrices mediating fundamental biological processes such as morphogenesis, hearing, fertilization, and antibacterial defense is driven by a ubiquitous polymerization module known as zona pellucida (ZP) "domain". Despite the conservation of this element from hydra to humans, no detailed information is available on the filamentous conformation of any ZP module protein. Here, we report a cryo-electron microscopy study of uromodulin (UMOD)/Tamm-Horsfall protein, the most abundant protein in human urine and an archetypal ZP module-containing molecule, in its mature homopolymeric state. UMOD forms a one-start helix with an unprecedented 180-degree twist between subunits enfolded by interdomain linkers that have completely reorganized as a result of propeptide dissociation. Lateral interaction between filaments in the urine generates sheets exposing a checkerboard of binding sites to capture uropathogenic bacteria, and UMOD-based models of heteromeric vertebrate egg coat filaments identify a common sperm-binding region at the interface between subunits.

Structural basis for the inhibition of cGAS by nucleosomes

The cyclic guanosine monophosphate-adenosine monophosphate synthase (cGAS) senses invasion of pathogenic DNA and stimulates inflammatory signaling, autophagy, and apoptosis. Organization of host DNA into nucleosomes was proposed to limit cGAS autoinduction, but the underlying mechanism was unknown. Here, we report the structural basis for this inhibition. In the cryo-electron microscopy structure of the human cGAS-nucleosome core particle (NCP) complex, two cGAS monomers bridge two NCPs by binding the acidic patch of the histone H2A-H2B dimer and nucleosomal DNA. In this configuration, all three known cGAS DNA binding sites, required for cGAS activation, are repurposed or become inaccessible, and cGAS dimerization, another prerequisite for activation, is inhibited. Mutating key residues linking cGAS and the acidic patch alleviates nucleosomal inhibition. This study establishes a structural framework for why cGAS is silenced on chromatinized self-DNA.

Advances in methods for atomic resolution macromolecular structure determination

Recent technical advances have dramatically increased the power and scope of structural biology. New developments in high-resolution cryo-electron microscopy, serial X-ray crystallography, and electron diffraction have been especially transformative. Here we highlight some of the latest advances and current challenges at the frontiers of atomic resolution methods for elucidating the structures and dynamical properties of macromolecules and their complexes.

Practical Considerations for Atomistic Structure Modeling with Cryo-EM Maps

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

Advances in Structure Modeling Methods for Cryo-Electron Microscopy Maps

Cryo-electron microscopy (cryo-EM) has now become a widely used technique for structure determination of macromolecular complexes. For modeling molecular structures from density maps of different resolutions, many algorithms have been developed. These algorithms can be categorized into rigid fitting, flexible fitting, and de novo modeling methods. It is also observed that machine learning (ML) techniques have been increasingly applied following the rapid progress of the ML field. Here, we review these different categories of macromolecule structure modeling methods and discuss their advances over time.