Atomic modeling of cryo-electron microscopy reconstructions--joint refinement of model and imaging parameters

Local heterogeneity analysis of crystallographic and cryo-EM maps using shell-approximation

In X-ray crystallography and cryo-EM, experimental maps can be heterogeneous, showing different level of details in different regions. In this work we interpret heterogeneity in terms of two parameters, assigned individually for each atom, combining the conventional atomic displacement parameter with the resolution of the atomic image in the map. We propose a local real-space procedure to estimate the values of these heterogeneity parameters, assuming that a fragment of the density map and atomic positions are given. The procedure is based on an analytic representation of the atomic image, as a function of the inhomogeneity parameters and atomic coordinates. In this article, we report the results of the tests both with maps simulated and those derived from experimental data. For simulated maps containing regions with different resolutions, the method determines the local map resolution around the atomic centers and the values of the displacement parameter with reasonable accuracy. For experimental maps, obtained as a Fourier synthesis of a given global resolution, estimated values of the local resolution are close to the global one, and the values of the estimated displacement parameters are close to the respective values of the closest atoms in the refined model. Shown successful applications of the proposed method to experimental crystallographic and cryo-EM maps can be seen as a practical proof of method.

Cryo-EM structure of adeno-associated virus 4 at 2.2 Å resolution

Adeno-associated virus (AAV) is the vector of choice for several approved gene-therapy treatments and is the basis for many ongoing clinical trials. Various strains of AAV exist (referred to as serotypes), each with their own transfection characteristics. Here, a high-resolution cryo-electron microscopy structure (2.2 Å) of AAV serotype 4 (AAV4) is presented. The receptor responsible for transduction of the AAV4 clade of AAV viruses (including AAV11, AAV12 and AAVrh32.33) is unknown. Other AAVs interact with the same cell receptor, adeno-associated virus receptor (AAVR), in one of two different ways. AAV5-like viruses interact exclusively with the polycystic kidney disease-like 1 (PKD1) domain of AAVR, while most other AAVs interact primarily with the PKD2 domain. A comparison of the present AAV4 structure with prior corresponding structures of AAV5, AAV2 and AAV1 in complex with AAVR provides a foundation for understanding why the AAV4-like clade is unable to interact with either PKD1 or PKD2 of AAVR. The conformation of the AAV4 capsid in variable regions I, III, IV and V on the viral surface appears to be sufficiently different from AAV2 to ablate binding with PKD2. Differences between AAV4 and AAV5 in variable region VII appear to be sufficient to exclude binding with PKD1.

Cross-Species Permissivity: Structure of a Goat Adeno-Associated Virus and Its Complex with the Human Receptor AAVR

Adeno-associated virus (AAV) is a small ssDNA satellite virus of high interest (in recombinant form) as a safe and effective gene therapy vector. AAV's human cell entry receptor (AAVR) contains polycystic kidney disease (PKD) domains bound by AAV. Seeking understanding of the spectrum of interactions, goat AAVGo.1 is investigated, because its host is the species most distant from human with reciprocal cross-species cell susceptibility. The structure of AAVGo.1, solved by cryo-EM to 2.9 Å resolution, is most similar to AAV5. Through ELISA (enzyme-linked immunosorbent assay) studies, it is shown that AAVGo.1 binds to human AAVR more strongly than do AAV2 or AAV5, and that it joins AAV5 in a class that binds exclusively to PKD domain 1 (PKD1), in contrast to other AAVs that interact primarily with PKD2. The AAVGo.1 cryo-EM structure of a complex with a PKD12 fragment of AAVR at 2.4 Å resolution shows PKD1 bound with minimal change in virus structure. There are only minor conformational adaptations in AAVR, but there is a near-rigid rotation of PKD1 with maximal displacement of the receptor domain by ~1 Å compared to PKD1 bound to AAV5. AAVGo.1 joins AAV5 as the second member of an emerging class of AAVs whose mode of receptor-binding is completely different from other AAVs, typified by AAV2. IMPORTANCE Adeno-associated virus (AAV) is a small ssDNA satellite parvovirus. As a recombinant vector with a protein shell encapsidating a transgene, recombinant AAV (rAAV) is a leading delivery vehicle for gene therapy, with two FDA-approved treatments and 150 clinical trials for 30 diseases. The human entry receptor AAVR has five PKD domains. To date, all serotypes, except AAV5, have interacted primarily with the second PKD domain, PKD2. Goat is the AAV host most distant from human with cross-species cell infectivity. AAVGo.1 is similar in structure to AAV5, the two forming a class with a distinct mode of receptor-binding. Within the two classes, binding interactions are mostly conserved, giving an indication of the latitude available in modulating delivery vectors.

Analytic modeling of inhomogeneous-resolution maps in cryo-electron microscopy and crystallography

Refinement of macromolecular atomic models versus experimental maps in crystallography and cryo-electron microscopy is a critical step in structure solution. For an appropriate comparison, model maps should mimic the imperfections in the experimental maps, mainly atomic disorder and limited resolution, which are often inhomogeneous over the molecular region. In the suggested method, these model maps are calculated as the sum of atomic contributions expressed through a specifically designed function describing a solitary spherical wave. Thanks to this function, atomic contributions are analytically expressed through their atomic displacement parameter and local resolution, a value now associated with each atom. Such a full analytic dependence of inhomogeneous-resolution map values on model parameters permits the refinement of all of these parameters together.

Toward Real Real-Space Refinement of Atomic Models

High-quality atomic models providing structural information are the results of their refinement versus diffraction data (reciprocal-space refinement), or versus experimental or experimentally based maps (real-space refinement). A proper real-space refinement can be achieved by comparing such a map with a map calculated from the atomic model. Similar to density distributions, the maps of a limited and even inhomogeneous resolution can also be calculated as sums of terms, known as atomic images, which are three-dimensional peaky functions surrounded by Fourier ripples. These atomic images and, consequently, the maps for the respective models, can be expressed analytically as functions of coordinates, atomic displacement parameters, and the local resolution. This work discusses the practical feasibility of such calculation for the real-space refinement of macromolecular atomic models.

Cryo-electron Microscopy of Adeno-associated Virus

Adeno-associated virus (AAV) has a single-stranded DNA genome encapsidated in a small icosahedrally symmetric protein shell with 60 subunits. AAV is the leading delivery vector in emerging gene therapy treatments for inherited disorders, so its structure and molecular interactions with human hosts are of intense interest. A wide array of electron microscopic approaches have been used to visualize the virus and its complexes, depending on the scientific question, technology available, and amenability of the sample. Approaches range from subvolume tomographic analyses of complexes with large and flexible host proteins to detailed analysis of atomic interactions within the virus and with small ligands at resolutions as high as 1.6 Å. Analyses have led to the reclassification of glycan receptors as attachment factors, to structures with a new-found receptor protein, to identification of the epitopes of antibodies, and a new understanding of possible neutralization mechanisms. AAV is now well-enough characterized that it has also become a model system for EM methods development. Heralding a new era, cryo-EM is now also being deployed as an analytic tool in the process development and production quality control of high value pharmaceutical biologics, namely AAV vectors.

Ten things I `hate' about refinement

Macromolecular refinement is an optimization process that aims to produce the most likely macromolecular structural model in the light of experimental data. As such, macromolecular refinement is one of the most complex optimization problems in wide use. Macromolecular refinement programs have to deal with the complex relationship between the parameters of the atomic model and the experimental data, as well as a large number of types of prior knowledge about chemical structure. This paper draws attention to areas of unfinished business in the field of macromolecular refinement. In it, we describe ten refinement topics that we think deserve attention and discuss directions leading to macromolecular refinement software that would make the best use of modern computer resources to meet the needs of structural biologists of the twenty-first century.

Adeno-Associated Virus (AAV-DJ)-Cryo-EM Structure at 1.56 Å Resolution

Adeno-associated virus is the leading viral vector for gene therapy. AAV-DJ is a recombinant variant developed for tropism to the liver. The AAV-DJ structure has been determined to 1.56 Å resolution through cryo-electron microscopy (cryo-EM). Only apoferritin is reported in preprints at 1.6 Å or higher resolution, and AAV-DJ nearly matches the highest resolutions ever attained through X-ray diffraction of virus crystals. However, cryo-EM has the advantage that most of the hydrogens are clear, improving the accuracy of atomic refinement, and removing ambiguity in hydrogen bond identification. Outside of secondary structures where hydrogen bonding was predictable a priori, the networks of hydrogen bonds coming from direct observation of hydrogens and acceptor atoms are quite different from those inferred even at 2.8 Å resolution. The implications for understanding viral assembly mean that cryo-EM will likely become the favored approach for high resolution structural virology.

Structure of the gene therapy vector, adeno-associated virus with its cell receptor, AAVR

Adeno-associated virus (AAV) vectors are preeminent in emerging clinical gene therapies. Generalizing beyond the most tractable genetic diseases will require modulation of cell specificity and immune neutralization. Interactions of AAV with its cellular receptor, AAVR, are key to understanding cell-entry and trafficking with the rigor needed to engineer tissue-specific vectors. Cryo-electron tomography shows ordered binding of part of the flexible receptor to the viral surface, with distal domains in multiple conformations. Regions of the virus and receptor in close physical proximity can be identified by cross-linking/mass spectrometry. Cryo-electron microscopy with a two-domain receptor fragment reveals the interactions at 2.4 Å resolution. AAVR binds between AAV's spikes on a plateau that is conserved, except in one clade whose structure is AAVR-incompatible. AAVR's footprint overlaps the epitopes of several neutralizing antibodies, prompting a re-evaluation of neutralization mechanisms. The structure provides a roadmap for experimental probing and manipulation of viral-receptor interactions.

Advances in image processing for single-particle analysis by electron cryomicroscopy and challenges ahead

Electron cryomicroscopy (cryoEM) is essential for the study and functional understanding of non-crystalline macromolecules such as proteins. These molecules cannot be imaged using X-ray crystallography or other popular methods. CryoEM has been successfully used to visualize macromolecular complexes such as ribosomes, viruses, and ion channels. Determination of structural models of these at various conformational states leads to insight on how these molecules function. Recent advances in imaging technology have given cryoEM a scientific rebirth. As a result of these technological advances image processing and analysis have yielded molecular structures at atomic resolution. Nevertheless there continue to be challenges in image processing, and in this article we will touch on the most essential in order to derive an accurate three-dimensional model from noisy projection images. Traditional approaches, such as k-means clustering for class averaging, will be provided as background. We will then highlight new approaches for each image processing subproblem, including a 3D reconstruction method for asymmetric molecules using just two projection images and deep learning algorithms for automated particle picking.

Anti-diabetic drug binding site in a mammalian KATP channel revealed by Cryo-EM

Sulfonylureas are anti-diabetic medications that act by inhibiting pancreatic K_ATP channels composed of SUR1 and Kir6.2. The mechanism by which these drugs interact with and inhibit the channel has been extensively investigated, yet it remains unclear where the drug binding pocket resides. Here, we present a cryo-EM structure of a hamster SUR1/rat Kir6.2 channel bound to a high-affinity sulfonylurea drug glibenclamide and ATP at 3.63 Å resolution, which reveals unprecedented details of the ATP and glibenclamide binding sites. Importantly, the structure shows for the first time that glibenclamide is lodged in the transmembrane bundle of the SUR1-ABC core connected to the first nucleotide binding domain near the inner leaflet of the lipid bilayer. Mutation of residues predicted to interact with glibenclamide in our model led to reduced sensitivity to glibenclamide. Our structure provides novel mechanistic insights of how sulfonylureas and ATP interact with the K_ATP channel complex to inhibit channel activity.

Structure of the calcium-dependent type 2 secretion pseudopilus

Many Gram-negative bacteria use type 2 secretion systems (T2SS) to secrete proteins involved in virulence and adaptation. Transport of folded proteins via T2SS nanomachines requires the assembly of inner membrane-anchored fibers called pseudopili. Although efficient pseudopilus assembly is essential for protein secretion, structure-based functional analyses are required to unravel the mechanistic link between these processes. Here, we report an atomic model for a T2SS pseudopilus from Klebsiella oxytoca, obtained by fitting the NMR structure of its calcium-bound subunit PulG into the ~ 5 Å resolution cryo-electron microscopy (cryoEM) reconstruction of assembled fibers. This structure reveals the comprehensive network of inter-subunit contacts and unexpected features, including a disordered central region of the PulG helical stem, and highly flexible C-terminal residues on the fiber surface. NMR, mutagenesis and functional analyses highlight the key role of calcium in PulG folding and stability. Fiber disassembly in the absence of calcium provides a basis for pseudopilus length control, essential for protein secretion, and supports the Archimedes' screw model for T2S mechanism.

The 2.8 Å Electron Microscopy Structure of Adeno-Associated Virus-DJ Bound by a Heparinoid Pentasaccharide

Atomic structures of adeno-associated virus (AAV)-DJ, alone and in complex with fondaparinux, have been determined by cryoelectron microscopy at 3 Å resolution. The gene therapy vector, AAV-DJ, is a hybrid of natural serotypes that was previously derived by directed evolution, selecting for hepatocyte entry and resistance to neutralization by human serum. The structure of AAV-DJ differs from that of parental serotypes in two regions where neutralizing antibodies bind, so immune escape appears to have been the primary driver of AAV-DJ's directed evolution. Fondaparinux is an analog of cell surface heparan sulfate to which several AAVs bind during entry. Fondaparinux interacts with viral arginines at a known heparin binding site, without the large conformational changes whose presence was controversial in low-resolution imaging of AAV2-heparin complexes. The glycan density suggests multi-modal binding that could accommodate sequence variation and multivalent binding along a glycan polymer, consistent with a role in attachment, prior to more specific interactions with a receptor protein mediating entry.

Cryo-EM structure of the ATP-sensitive potassium channel illuminates mechanisms of assembly and gating

K_ATP channels are metabolic sensors that couple cell energetics to membrane excitability. In pancreatic β-cells, channels formed by SUR1 and Kir6.2 regulate insulin secretion and are the targets of antidiabetic sulfonylureas. Here, we used cryo-EM to elucidate structural basis of channel assembly and gating. The structure, determined in the presence of ATP and the sulfonylurea glibenclamide, at ~6 Å resolution reveals a closed Kir6.2 tetrameric core with four peripheral SUR1s each anchored to a Kir6.2 by its N-terminal transmembrane domain (TMD0). Intricate interactions between TMD0, the loop following TMD0, and Kir6.2 near the proposed PIP₂ binding site, and where ATP density is observed, suggest SUR1 may contribute to ATP and PIP₂ binding to enhance Kir6.2 sensitivity to both. The SUR1-ABC core is found in an unusual inward-facing conformation whereby the two nucleotide binding domains are misaligned along a two-fold symmetry axis, revealing a possible mechanism by which glibenclamide inhibits channel activity.

DOI:http://dx.doi.org/10.7554/eLife.24149.001

The hormone insulin reduces blood sugar levels by encouraging fat, muscle and other body cells to take up sugar. When blood sugar levels rise following a meal, cells within the pancreas known as beta cells should release insulin. In people with diabetes, the beta cells fail to release insulin, meaning that the high blood sugar levels are not corrected.

When blood sugar levels are high, beta cells generate more energy in the form of ATP molecules. The increased level of ATP causes channels called ATP-sensitive potassium (K_ATP) channels in the membrane of the cell to close. This triggers a cascade of events that leads to the release of insulin.

Some treatments for diabetes alter how the K_ATP channels work. For example, a widely prescribed medication called glibenclamide (also known as glyburide in the United States) stimulates the release of insulin by preventing the flow of potassium through K_ATP channels. It remains unknown exactly how ATP and glibenclamide interact with the channel’s molecular structure to stop the flow of potassium ions.

K_ATP channels are made up of two proteins called SUR1 and Kir6.2. To investigate the structure of the K_ATP channel, Martin et al. purified channels made of the hamster form of the SUR1 protein and the mouse form of Kir6.2, which each closely resemble their human counterparts. The channels were purified in the presence of ATP and glibenclamide and were then rapidly frozen to preserve their structure, which allowed them to be visualized individually using electron microscopy. By analyzing the images taken from many channels, Martin et al. constructed a highly detailed, three-dimensional map of the K_ATP channel. The structure revealed by this map shows how SUR1 and Kir6.2 work together and provides insight into how ATP and glibenclamide interact with the channel to block the flow of potassium, and hence stimulate the release of insulin.

An important next step will be to improve the structure to more clearly identify where ATP and glibenclamide bind to the K_ATP channel. It will also be important to study the structures of channels that are bound to other regulatory molecules. This will help researchers to fully understand how K_ATP channels located throughout the body operate under healthy and diseased conditions. This knowledge will aid in the design of more effective drugs to treat several devastating diseases caused by defective K_ATP channels.

DOI:http://dx.doi.org/10.7554/eLife.24149.002

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

The influence of frame alignment with dose compensation on the quality of single particle reconstructions

As direct electron detection devices in cryo-electron microscopy become ubiquitous, the field is now ripe for new developments in image analysis techniques that take advantage of their increased SNR coupled with their high-throughput frame collection abilities. In approaching atomic resolution of native-like biomolecules, the accurate extraction of structural locations and orientations of side-chains from frames depends not only on the electron dose that a sample receives but also on the ability to accurately estimate the CTF. Here we use a new 2.8Å resolution structure of a recombinant gene therapy virus, AAV-DJ with Arixtra, imaged on an FEI Titan Krios with a DE-20 direct electron detector to probe new metrics including relative side-chain density and ResLog analysis for optimizing the compensation of electron beam damage and to characterize the factors that are limiting the resolution of the reconstruction. The influence of dose compensation on the accuracy of CTF estimation and particle classifiability are also presented. We show that rigorous dose compensation allows for better particle classifiability and greater recovery of structural information from negatively charged, electron-sensitive side-chains, resulting in a more accurate macromolecular model.

Biogenesis and structure of a type VI secretion membrane core complex

Bacteria share their ecological niches with other microbes. The bacterial type VI secretion system is one of the key players in microbial competition, as well as being an important virulence determinant during bacterial infections. It assembles a nano-crossbow-like structure in the cytoplasm of the attacker cell that propels an arrow made of a haemolysin co-regulated protein (Hcp) tube and a valine-glycine repeat protein G (VgrG) spike and punctures the prey's cell wall. The nano-crossbow is stably anchored to the cell envelope of the attacker by a membrane core complex. Here we show that this complex is assembled by the sequential addition of three type VI subunits (Tss)-TssJ, TssM and TssL-and present a structure of the fully assembled complex at 11.6 Å resolution, determined by negative-stain electron microscopy. With overall C5 symmetry, this 1.7-megadalton complex comprises a large base in the cytoplasm. It extends in the periplasm via ten arches to form a double-ring structure containing the carboxy-terminal domain of TssM (TssMct) and TssJ that is anchored in the outer membrane. The crystal structure of the TssMct-TssJ complex coupled to whole-cell accessibility studies suggest that large conformational changes induce transient pore formation in the outer membrane, allowing passage of the attacking Hcp tube/VgrG spike.

GPU-accelerated analysis and visualization of large structures solved by molecular dynamics flexible fitting

Hybrid structure fitting methods combine data from cryo-electron microscopy and X-ray crystallography with molecular dynamics simulations for the determination of all-atom structures of large biomolecular complexes. Evaluating the quality-of-fit obtained from hybrid fitting is computationally demanding, particularly in the context of a multiplicity of structural conformations that must be evaluated. Existing tools for quality-of-fit analysis and visualization have previously targeted small structures and are too slow to be used interactively for large biomolecular complexes of particular interest today such as viruses or for long molecular dynamics trajectories as they arise in protein folding. We present new data-parallel and GPU-accelerated algorithms for rapid interactive computation of quality-of-fit metrics linking all-atom structures and molecular dynamics trajectories to experimentally-determined density maps obtained from cryo-electron microscopy or X-ray crystallography. We evaluate the performance and accuracy of the new quality-of-fit analysis algorithms vis-à-vis existing tools, examine algorithm performance on GPU-accelerated desktop workstations and supercomputers, and describe new visualization techniques for results of hybrid structure fitting methods.

Finding the right fit: chiseling structures out of cryo-electron microscopy maps

Cryo-electron microscopy is a central tool for studying the architecture of macromolecular complexes at subnanometer resolution. Interpretation of an electron microscopy map requires its computational integration with data about the structure's components from all available sources, notably atomic models. Selecting a protocol for EM density-guided integrative structural modeling depends on the resolution and quality of the EM map as well as the available complimentary datasets. Here, we review rigid, flexible, and de novo integrative fitting into EM maps and provide guidelines and considerations for the design of modeling experiments. Finally, we discuss efforts towards establishing unified criteria for map and model assessment and validation.

Conformational States of macromolecular assemblies explored by integrative structure calculation

A detailed description of macromolecular assemblies in multiple conformational states can be very valuable for understanding cellular processes. At present, structural determination of most assemblies in different biologically relevant conformations cannot be achieved by a single technique and thus requires an integrative approach that combines information from multiple sources. Different techniques require different computational methods to allow efficient and accurate data processing and analysis. Here, we summarize the latest advances and future challenges in computational methods that help the interpretation of data from two techniques—mass spectrometry and three-dimensional cryo-electron microscopy (with focus on alignment and classification of heterogeneous subtomograms from cryo-electron tomography). We evaluate how new developments in these two broad fields will lead to further integration with atomic structures to broaden our picture of the dynamic behavior of assemblies in their native environment.

The joys and perils of flexible fitting

While performing their functions, biological macromolecules often form large, dynamically changing macromolecular assemblies. Only a relatively small number of such assemblies have been accessible to the atomic-resolution techniques X-ray crystallography and NMR. Electron microscopy in conjunction with image reconstruction has become the preferred alternative for revealing the structures of such macromolecular complexes. However, for most assemblies the achievable resolution is too low to allow accurate atomic modeling directly from the data. Yet, useful models often can be obtained by fitting atomic models of individual components into a low-resolution reconstruction of the entire assembly. Several algorithms for achieving optimal fits in this context were developed recently, many allowing considerable degrees of flexibility to account for binding-induced conformational changes of the assembly components. This chapter describes the advantages and potential pitfalls of these methods and puts them into perspective with alternative approaches such as iterative modular fitting of rigid-body domains.

ResLog plots as an empirical metric of the quality of cryo-EM reconstructions

Compared to the field of X-ray crystallography, the field of single particle three-dimensional electron microscopy has few reliable metrics for assessing the quality of 3D reconstructions. New metrics are needed that can determine whether a given 3D reconstruction accurately reflects the structure of the particles from which it was derived or instead depicts a plausible though incorrect structure due to coarse misalignment of particles. Here an empirical procedure is presented for differentiating between a reconstruction with well-aligned particles and a reconstruction with grossly misclassified particles. For a given dataset, 3D reconstructions are computed from subsets of particles with decreasing numbers of particles contributing to the reconstruction. A plot of inverse resolution vs. the logarithm of the number of particles (a "ResLog" plot) provides metrics for the reliability of the reconstruction and the overall quality of the dataset and processing. Specifically, the y-intercept of a regression line provides a measure of the relative accuracy of the particle alignment and classification, and the slope is an indicator of the overall data quality including the imaging conditions and processing steps. ResLog plots can also be used to optimize conditions for data collection and reconstruction parameters. Although resolution estimates can vary by method of calculation, ResLog-derived parameters are consistent whether calculated by Fourier shell correlation or Fourier neighbor correlation, or a new coordinate-based metric that serves as a yardstick for structures where atomic coordinates are available. ResLog plots could become part of a standard set of parameters to be included in 3D reconstruction reports.

Electron microscopy analysis of a disaccharide analog complex reveals receptor interactions of adeno-associated virus

Mechanistic studies of macromolecular complexes often feature X-ray structures of complexes with bound ligands. The attachment of adeno-associated virus (AAV) to cell surface glycosaminoglycans (GAGs) is an example that has not proven amenable to crystallography, because the binding of GAG analogs disrupts lattice contacts. The interactions of AAV with GAGs are of interest in mediating the cell specificity of AAV-based gene therapy vectors. Previous electron microscopy led to differing conclusions on the exact binding site and the existence of large ligand-induced conformational changes in the virus. Conformational changes are expected during cell entry, but it has remained unclear whether the electron microscopy provided evidence of their induction by GAG-binding. Taking advantage of automated data collection, careful processing and new methods of structure refinement, the structure of AAV-DJ complexed with sucrose octasulfate is determined by electron microscopy difference map analysis to 4.8Å resolution. At this higher resolution, individual sulfate groups are discernible, providing a stereochemical validation of map interpretation, and highlighting interactions with two surface arginines that have been implicated in genetic studies. Conformational changes induced by the SOS are modest and limited to the loop most directly interacting with the ligand. While the resolution attainable will depend on sample order and other factors, there are an increasing number of macromolecular complexes that can be studied by cryo-electron microscopy at resolutions beyond 5Å, for which the approaches used here could be used to characterize the binding of inhibitors and other small molecule effectors when crystallography is not tractable.