Situs: A package for docking crystal structures into low-resolution maps from electron microscopy

Systematic comparison of molecular conformations of H+,K+-ATPase reveals an important contribution of the A-M2 linker for the luminal gating

Gastric H⁺,K⁺-ATPase, an ATP-driven proton pump responsible for gastric acidification, is a molecular target for anti-ulcer drugs. Here we show its cryo-electron microscopy (EM) structure in an E2P analog state, bound to magnesium fluoride (MgF), and its K⁺-competitive antagonist SCH28080, determined at 7 Å resolution by electron crystallography of two-dimensional crystals. Systematic comparison with other E2P-related cryo-EM structures revealed that the molecular conformation in the (SCH)E2·MgF state is remarkably distinguishable. Although the azimuthal position of the A domain of the (SCH)E2·MgF state is similar to that in the E2·AlF (aluminum fluoride) state, in which the transmembrane luminal gate is closed, the arrangement of transmembrane helices in the (SCH)E2·MgF state shows a luminal-open conformation imposed on by bound SCH28080 at its luminal cavity, based on observations of the structure in the SCH28080-bound E2·BeF (beryllium fluoride) state. The molecular conformation of the (SCH)E2·MgF state thus represents a mixed overall structure in which its cytoplasmic and luminal half appear to be independently modulated by a phosphate analog and an antagonist bound to the respective parts of the enzyme. Comparison of the molecular conformations revealed that the linker region connecting the A domain and the transmembrane helix 2 (A-M2 linker) mediates the regulation of luminal gating. The mechanistic rationale underlying luminal gating observed in H⁺,K⁺-ATPase is consistent with that observed in sarcoplasmic reticulum Ca²⁺-ATPase and other P-type ATPases and is most likely conserved for the P-type ATPase family in general.

High-resolution modeling of protein structures based on flexible fitting of low-resolution structural data

To circumvent the difficulty of directly solving high-resolution biomolecular structures, low-resolution structural data from Cryo-electron microscopy (EM) and small angle solution X-ray scattering (SAXS) are increasingly used to explore multiple conformational states of biomolecular assemblies. One promising venue to obtain high-resolution structural models from low-resolution data is via data-constrained flexible fitting. To this end, we have developed a new method based on a coarse-grained Cα-only protein representation, and a modified form of the elastic network model (ENM) that allows large-scale conformational changes while maintaining the integrity of local structures including pseudo-bonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified ENM energy with an EM/SAXS-fitting score and a collision energy that penalizes steric collisions. Unlike some previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled with accuracy. This method is also highly efficient in computing time. We have demonstrated our method using adenylate kinase as a test case which undergoes a large open-to-close conformational change. The EM-fitting method is available at a web server (http://enm.lobos.nih.gov), and the SAXS-fitting method is available as a pre-compiled executable upon request.

Structural and biochemical characterization of the Cop9 signalosome CSN5/CSN6 heterodimer

The Cop9 signalosome complex (CSN) regulates the functional cycle of the major E3 ubiquitin ligase family, the cullin RING E3 ubiquitin ligases (CRLs). Activated CRLs are covalently modified by the ubiquitin-like protein Nedd8 (neural precursor cell expressed developmentally down-regulated protein 8). CSN serves an essential role in myriad cellular processes by reversing this modification through the isopeptidase activity of its CSN5 subunit. CSN5 alone is inactive due to an auto-inhibited conformation of its catalytic domain. Here we report the molecular basis of CSN5 catalytic domain activation and unravel a molecular hierarchy in CSN deneddylation activity. The association of CSN5 and CSN6 MPN (for Mpr1/Pad1 N-terminal) domains activates its isopeptidase activity. The CSN5/CSN6 module, however, is inefficient in CRL deneddylation, indicating a requirement of further elements in this reaction such as other CSN subunits. A hybrid molecular model of CSN5/CSN6 provides a structural framework to explain these functional observations. Docking this model into a published CSN electron density map and using distance constraints obtained from cross-linking coupled to mass-spectrometry, we find that the C-termini of the CSN subunits could form a helical bundle in the centre of the structure. They likely play a key scaffolding role in the spatial organization of CSN and precise positioning of the dimeric MPN catalytic core.

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.

Automated target segmentation and real space fast alignment methods for high-throughput classification and averaging of crowded cryo-electron subtomograms

MOTIVATION

Cryo-electron tomography allows the imaging of macromolecular complexes in near living conditions. To enhance the nominal resolution of a structure it is necessary to align and average individual subtomograms each containing identical complexes. However, if the sample of complexes is heterogeneous, it is necessary to first classify subtomograms into groups of identical complexes. This task becomes challenging when tomograms contain mixtures of unknown complexes extracted from a crowded environment. Two main challenges must be overcomed: First, classification of subtomograms must be performed without knowledge of template structures. However, most alignment methods are too slow to perform reference-free classification of a large number of (e.g. tens of thousands) of subtomograms. Second, subtomograms extracted from crowded cellular environments, contain often fragments of other structures besides the target complex. However, alignment methods generally assume that each subtomogram only contains one complex. Automatic methods are needed to identify the target complexes in a subtomogram even when its shape is unknown.

RESULTS

In this article, we propose an automatic and systematic method for the isolation and masking of target complexes in subtomograms extracted from crowded environments. Moreover, we also propose a fast alignment method using fast rotational matching in real space. Our experiments show that, compared with our previously proposed fast alignment method in reciprocal space, our new method significantly improves the alignment accuracy for highly distorted and especially crowded subtomograms. Such improvements are important for achieving successful and unbiased high-throughput reference-free structural classification of complexes inside whole-cell tomograms.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structural characterization of Staphylococcus aureus biotin protein ligase and interaction partners: an antibiotic target

The essential metabolic enzyme biotin protein ligase (BPL) is a potential target for the development of new antibiotics required to combat drug-resistant pathogens. Staphylococcus aureus BPL (SaBPL) is a bifunctional protein, possessing both biotin ligase and transcription repressor activities. This positions BPL as a key regulator of several important metabolic pathways. Here, we report the structural analysis of both holo- and apo-forms of SaBPL using X-ray crystallography. We also present small-angle X-ray scattering data of SaBPL in complex with its biotin-carboxyl carrier protein substrate as well as the SaBPL:DNA complex that underlies repression. This has revealed the molecular basis of ligand (biotinyl-5'-AMP) binding and conformational changes associated with catalysis and repressor function. These data provide new information to better understand the bifunctional activities of SaBPL and to inform future strategies for antibiotic discovery.

Fibrinogen species as resolved by HPLC-SAXS data processing within the UltraScan Solution Modeler (US-SOMO) enhanced SAS module

Fibrinogen is a large heterogeneous aggregation/degradation-prone protein playing a central role in blood coagulation and associated pathologies, whose structure is not completely resolved. When a high-molecular-weight fraction was analyzed by size-exclusion high-performance liquid chromatography/small-angle X-ray scattering (HPLC-SAXS), several composite peaks were apparent and because of the stickiness of fibrinogen the analysis was complicated by severe capillary fouling. Novel SAS analysis tools developed as a part of the UltraScan Solution Modeler (US-SOMO; http://somo.uthscsa.edu/), an open-source suite of utilities with advanced graphical user interfaces whose initial goal was the hydrodynamic modeling of biomacromolecules, were implemented and applied to this problem. They include the correction of baseline drift due to the accumulation of material on the SAXS capillary walls, and the Gaussian decomposition of non-baseline-resolved HPLC-SAXS elution peaks. It was thus possible to resolve at least two species co-eluting under the fibrinogen main monomer peak, probably resulting from in-column degradation, and two others under an oligomers peak. The overall and cross-sectional radii of gyration, molecular mass and mass/length ratio of all species were determined using the manual or semi-automated procedures available within the US-SOMO SAS module. Differences between monomeric species and linear and sideways oligomers were thus identified and rationalized. This new US-SOMO version additionally contains several computational and graphical tools, implementing functionalities such as the mapping of residues contributing to particular regions of P(r), and an advanced module for the comparison of primary I(q) versus q data with model curves computed from atomic level structures or bead models. It should be of great help in multi-resolution studies involving hydrodynamics, solution scattering and crystallographic/NMR data.

Structural characterization of the bacteriophage T7 tail machinery

Most bacterial viruses need a specialized machinery, called "tail," to inject their genomes inside the bacterial cytoplasm without disrupting the cellular integrity. Bacteriophage T7 is a well characterized member of the Podoviridae family infecting Escherichia coli, and it has a short noncontractile tail that assembles sequentially on the viral head after DNA packaging. The T7 tail is a complex of around 2.7 MDa composed of at least four proteins as follows: the connector (gene product 8, gp8), the tail tubular proteins gp11 and gp12, and the fibers (gp17). Using cryo-electron microscopy and single particle image reconstruction techniques, we have determined the precise topology of the tail proteins by comparing the structure of the T7 tail extracted from viruses and a complex formed by recombinant gp8, gp11, and gp12 proteins. Furthermore, the order of assembly of the structural components within the complex was deduced from interaction assays with cloned and purified tail proteins. The existence of common folds among similar tail proteins allowed us to obtain pseudo-atomic threaded models of gp8 (connector) and gp11 (gatekeeper) proteins, which were docked into the corresponding cryo-EM volumes of the tail complex. This pseudo-atomic model of the connector-gatekeeper interaction revealed the existence of a common molecular architecture among viruses belonging to the three tailed bacteriophage families, strongly suggesting that a common molecular mechanism has been favored during evolution to coordinate the transition between DNA packaging and tail assembly.

Targeted conformational search with map-restrained self-guided Langevin dynamics: application to flexible fitting into electron microscopic density maps

We present a map-restrained self-guided Langevin dynamics (MapSGLD) simulation method for efficient targeted conformational search. The targeted conformational search represents simulations under restraints defined by experimental observations and/or by user specified structural requirements. Through map-restraints, this method provides an efficient way to maintain substructures and to set structure targets during conformational searching. With an enhanced conformational searching ability of self-guided Langevin dynamics, this approach is suitable for simulating large-scale conformational changes, such as the formation of macromolecular assemblies and transitions between different conformational states. Using several examples, we illustrate the application of this method in flexible fitting of atomic structures into density maps derived from cryo-electron microscopy.

NONUNIFORM FOURIER TRANSFORMS FOR RIGID-BODY AND MULTI-DIMENSIONAL ROTATIONAL CORRELATIONS

The task of evaluating correlations is central to computational structural biology. The rigid-body correlation problem seeks the rigid-body transformation (R, t), R ∈ SO(3), t ∈ ℝ³ that maximizes the correlation between a pair of input scalar-valued functions representing molecular structures. Exhaustive solutions to the rigid-body correlation problem take advantage of the fast Fourier transform to achieve a speedup either with respect to the sought translation or rotation. We present PFcorr, a new exhaustive solution, based on the non-equispaced SO(3) Fourier transform, to the rigid-body correlation problem; unlike previous solutions, ours achieves a combination of translational and rotational speedups without requiring equispaced grids. PFcorr can be straightforwardly applied to a variety of problems in protein structure prediction and refinement that involve correlations under rigid-body motions of the protein. Additionally, we show how it applies, along with an appropriate flexibility model, to analogs of the above problems in which the flexibility of the protein is relevant.

Computational methods for constructing protein structure models from 3D electron microscopy maps

Protein structure determination by cryo-electron microscopy (EM) has made significant progress in the past decades. Resolutions of EM maps have been improving as evidenced by recently reported structures that are solved at high resolutions close to 3Å. Computational methods play a key role in interpreting EM data. Among many computational procedures applied to an EM map to obtain protein structure information, in this article we focus on reviewing computational methods that model protein three-dimensional (3D) structures from a 3D EM density map that is constructed from two-dimensional (2D) maps. The computational methods we discuss range from de novo methods, which identify structural elements in an EM map, to structure fitting methods, where known high resolution structures are fit into a low-resolution EM map. A list of available computational tools is also provided.

A molecule-centered method for accelerating the calculation of hydrodynamic interactions in Brownian dynamics simulations containing many flexible biomolecules

Inclusion of hydrodynamic interactions (HIs) is essential in simulations of biological macromolecules that treat the solvent implicitly if the macromolecules are to exhibit correct translational and rotational diffusion. The present work describes the development and testing of a simple approach aimed at allowing more rapid computation of HIs in coarse-grained Brownian dynamics simulations of systems that contain large numbers of flexible macromolecules. The method combines a complete treatment of intramolecular HIs with an approximate treatment of the intermolecular HIs which assumes that the molecules are effectively spherical; all of the HIs are calculated at the Rotne-Prager-Yamakawa level of theory. When combined with Fixman's Chebyshev polynomial method for calculating correlated random displacements, the proposed method provides an approach that is simple to program but sufficiently fast that it makes it computationally viable to include HIs in large-scale simulations. Test calculations performed on very coarse-grained models of the pyruvate dehydrogenase (PDH) E2 complex and on oligomers of ParM (ranging in size from 1 to 20 monomers) indicate that the method reproduces the translational diffusion behavior seen in more complete HI simulations surprisingly well; the method performs less well at capturing rotational diffusion but its discrepancies diminish with increasing size of the simulated assembly. Simulations of residue-level models of two tetrameric protein models demonstrate that the method also works well when more structurally detailed models are used in the simulations. Finally, test simulations of systems containing up to 1024 coarse-grained PDH molecules indicate that the proposed method rapidly becomes more efficient than the conventional BD approach in which correlated random displacements are obtained via a Cholesky decomposition of the complete diffusion tensor.

Cross-validation in cryo-EM-based structural modeling

Single-particle cryo-EM is a powerful approach to determine the structure of large macromolecules and assemblies thereof in many cases at subnanometer resolution. It has become popular to refine or flexibly fit atomic models into density maps derived from cryo-EM experiments. These density maps are typically significantly lower in resolution than electron density maps obtained from X-ray diffraction experiments, such that the number of parameters that need to be determined is much larger than the number of experimental observables. Overfitting and misinterpretation of the density, thus, become a serious problem. For diffraction data, a cross-validation approach was introduced almost 20 y ago; however, no such approach has been described yet for structure refinement against cryo-EM density maps, although the overfitting problem is, because of the lower resolution, significantly larger. We present a cross-validation approach for real-space refinement against cryo-EM density maps in analogy to cross-validation typically used in crystallography. Our approach is able to detect overfitting and allows for optimizing the choice of restraints used in the refinement. The approach is shown on three protein structures with simulated data and experimental data of the rotavirus double-layer particle. Because cross-validation requires splitting the dataset into at least two independent sets, we further present an approach to quantify correlations between the structure factor sets. This analysis is also helpful for other cross-validation applications, such as refinements against diffraction data or 3D reconstructions of cryo-EM density maps.

Large terminase conformational change induced by connector binding in bacteriophage T7

During bacteriophage morphogenesis DNA is translocated into a preformed prohead by the complex formed by the portal protein, or connector, plus the terminase, which are located at an especial prohead vertex. The terminase is a powerful motor that converts ATP hydrolysis into mechanical movement of the DNA. Here, we have determined the structure of the T7 large terminase by electron microscopy. The five terminase subunits assemble in a toroid that encloses a channel wide enough to accommodate dsDNA. The structure of the complete connector-terminase complex is also reported, revealing the coupling between the terminase and the connector forming a continuous channel. The structure of the terminase assembled into the complex showed a different conformation when compared with the isolated terminase pentamer. To understand in molecular terms the terminase morphological change, we generated the terminase atomic model based on the crystallographic structure of its phage T4 counterpart. The docking of the threaded model in both terminase conformations showed that the transition between the two states can be achieved by rigid body subunit rotation in the pentameric assembly. The existence of two terminase conformations and its possible relation to the sequential DNA translocation may shed light into the molecular bases of the packaging mechanism of bacteriophage T7.

Details of ssDNA annealing revealed by an HSV-1 ICP8-ssDNA binary complex

Infected cell protein 8 (ICP8) from herpes simplex virus 1 was first identified as a single-strand (ss) DNA-binding protein. It is essential for, and abundant during, viral replication. Studies in vitro have shown that ICP8 stimulates model replication reactions, catalyzes annealing of complementary ssDNAs and, in combination with UL12 exonuclease, will catalyze ssDNA annealing homologous recombination. DNA annealing and strand transfer occurs within large oligomeric filaments of ssDNA-bound ICP8. We present the first 3D reconstruction of a novel ICP8-ssDNA complex, which seems to be the basic unit of the DNA annealing machine. The reconstructed volume consists of two nonameric rings containing ssDNA stacked on top of each other, corresponding to a molecular weight of 2.3 MDa. Fitting of the ICP8 crystal structure suggests a mechanism for the annealing reaction catalyzed by ICP8, which is most likely a general mechanism for protein-driven DNA annealing.

Enhanced sampling and overfitting analyses in structural refinement of nucleic acids into electron microscopy maps

Flexible fitting computational algorithms are often useful to interpret low-resolution maps of many macromolecular complexes generated by electron microscopy (EM) imaging. One such atomistic simulation technique is molecular dynamics flexible fitting (MDFF), which has been widely applied to generate structural models of large ribonucleoprotein assemblies such as the ribosome. We have previously shown that MDFF simulations of globular proteins are sensitive to the resolution of the target EM map and the strength of restraints used to preserve the secondary structure elements during fitting (Vashisth, H.; et al. Structure 2012, 20, 1453-1462). In this work, we aim to systematically examine the quality of structural models of various nucleic acids obtained via MDFF by varying the map resolution and the strength of structural restraints. We also demonstrate how an enhanced conformational sampling technique for proteins, temperature-accelerated molecular dynamics (TAMD), can be combined with MDFF for the structural refinement of nucleic acids in EM maps. Finally, we also demonstrate application of TAMD-assisted MDFF (TAMDFF) on a RNA/protein complex and suggest that TAMDFF is a viable strategy for enhanced conformational fitting in target maps of ribonucleoprotein complexes.

Computer Simulations of the Bacterial Cytoplasm

Ever since the pioneering work of Minton, it has been recognized that the highly crowded interior of biological cells has the potential to cause dramatic changes to both the kinetics and thermodynamics of protein folding and association events relative to behavior that might be observed in dilute solution conditions. One very productive way to explore the effects of crowding on protein behavior has been to use macromolecular crowding agents that exclude volume without otherwise strongly interacting with the protein under study. An alternative, complementary approach to understanding the potential differences between behavior in vivo and in vitro is to develop simulation models that explicitly attempt to model intracellular environments at the molecular scale, and that thereby can be used to directly monitor biophysical behavior in conditions that accurately mimic those encountered in vivo. It is with studies of this type that the present review will be concerned. We review in detail four published studies that have attempted to simulate the structure and dynamics of the bacterial cytoplasm and that have each explored different biophysical aspects of the cellular interior. While each of these studies has yielded important new insights, there are important questions that remain to be resolved in terms of determining the relative contributions made by energetic and hydrodynamic interactions to the diffusive behavior of macromolecules and to the thermodynamics of protein folding and associations in vivo. Some possible new directions for future generation simulation models of the cytoplasm are outlined.

Developing advanced X-ray scattering methods combined with crystallography and computation

The extensive use of small angle X-ray scattering (SAXS) over the last few years is rapidly providing new insights into protein interactions, complex formation and conformational states in solution. This SAXS methodology allows for detailed biophysical quantification of samples of interest. Initial analyses provide a judgment of sample quality, revealing the potential presence of aggregation, the overall extent of folding or disorder, the radius of gyration, maximum particle dimensions and oligomerization state. Structural characterizations include ab initio approaches from SAXS data alone, and when combined with previously determined crystal/NMR, atomistic modeling can further enhance structural solutions and assess validity. This combination can provide definitions of architectures, spatial organizations of protein domains within a complex, including those not determined by crystallography or NMR, as well as defining key conformational states of a protein interaction. SAXS is not generally constrained by macromolecule size, and the rapid collection of data in a 96-well plate format provides methods to screen sample conditions. This includes screening for co-factors, substrates, differing protein or nucleotide partners or small molecule inhibitors, to more fully characterize the variations within assembly states and key conformational changes. Such analyses may be useful for screening constructs and conditions to determine those most likely to promote crystal growth of a complex under study. Moreover, these high throughput structural determinations can be leveraged to define how polymorphisms affect assembly formations and activities. This is in addition to potentially providing architectural characterizations of complexes and interactions for systems biology-based research, and distinctions in assemblies and interactions in comparative genomics. Thus, SAXS combined with crystallography/NMR and computation provides a unique set of tools that should be considered as being part of one's repertoire of biophysical analyses, when conducting characterizations of protein and other macromolecular interactions.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Structural and biochemical studies of RIG-I antiviral signaling

Retinoic acid-inducible gene I (RIG-I) is an important pattern recognition receptor that detects viral RNA and triggers the production of type-I interferons through the downstream adaptor MAVS (also called IPS-1, CARDIF, or VISA). A series of structural studies have elaborated some of the mechanisms of dsRNA recognition and activation of RIG-I. Recent studies have proposed that K63-linked ubiquitination of, or unanchored K63-linked polyubiquitin binding to RIG-I positively regulates MAVS-mediated antiviral signaling. Conversely phosphorylation of RIG-I appears to play an inhibitory role in controlling RIG-I antiviral signal transduction. Here we performed a combined structural and biochemical study to further define the regulatory features of RIG-I signaling. ATP and dsRNA binding triggered dimerization of RIG-I with conformational rearrangements of the tandem CARD domains. Full length RIG-I appeared to form a complex with dsRNA in a 2:2 molar ratio. Compared with the previously reported crystal structures of RIG-I in inactive state, our electron microscopic structure of full length RIG-I in complex with blunt-ended dsRNA, for the first time, revealed an exposed active conformation of the CARD domains. Moreover, we found that purified recombinant RIG-I proteins could bind to the CARD domain of MAVS independently of dsRNA, while S8E and T170E phosphorylation-mimicking mutants of RIG-I were defective in binding E3 ligase TRIM25, unanchored K63-linked polyubiquitin, and MAVS regardless of dsRNA. These findings suggested that phosphorylation of RIG inhibited downstream signaling by impairing RIG-I binding with polyubiquitin and its interaction with MAVS.

Using enhanced sampling and structural restraints to refine atomic structures into low-resolution electron microscopy maps

For a variety of problems in structural biology, low-resolution maps generated by electron microscopy imaging are often interpreted with the help of various flexible-fitting computational algorithms. In this work, we systematically analyze the quality of final models of various proteins obtained via molecular dynamics flexible fitting (MDFF) by varying the map-resolution, strength of structural restraints, and the steering forces. We find that MDFF can be extended to understand conformational changes in lower-resolution maps if larger structural restraints and lower steering forces are used to prevent overfitting. We further show that the capabilities of MDFF can be extended by combining it with an enhanced conformational sampling method, temperature-accelerated molecular dynamics (TAMD). Specifically, either TAMD can be used to generate better starting configurations for MDFF fitting or TAMD-assisted MDFF (TAMDFF) can be performed to accelerate conformational search in atomistic simulations.

Modeling, docking, and fitting of atomic structures to 3D maps from cryo-electron microscopy

Electron microscopy (EM) and image analysis offer an effective approach for determining the three-dimensional structure of macromolecular complexes. The versatility of these methods means that molecular species not normally amenable to other structural methods, e.g., X-ray crystallography and NMR spectroscopy, can be analyzed. However, the resolution of EM structures is often too low to provide an atomic model directly by chain tracing. Instead, a combination of modeling and fitting can be an effective way to analyze the EM structure at an atomic level, thus allowing localization of subunits or evaluation of conformational changes. Here we describe the steps involved in this process: building a homology model, fitting this model to an EM map, and using computational methods for docking of additional domains to the model. As an example, we illustrate the methods using an integral membrane protein, CopA, which functions to pump copper across the membrane in an ATP-dependent manner. In this example, we build a homology model based on the published atomic coordinates for a related calcium pump from sarcoplasmic reticulum (SERCA). After fitting this homology model to a 17 Å resolution EM map, computational software is used to dock a metal-binding domain (MBD) that is unique to the copper pump. Although this software identifies a number of plausible interfaces for docking, the constraints of the EM map steer us to select a unique solution. Thus, the synergy of these two methods allows us to describe both the location of the unknown MBD relative to the other cytoplasmic domains and the atomic details of the domain interface.

Tubular crystals and helical arrays: structural determination of HIV-1 capsid assemblies using iterative helical real-space reconstruction

Helical structures are important in many different life forms and are well-suited for structural studies by cryo-EM. A unique feature of helical objects is that a single projection image contains all the views needed to perform a three-dimensional (3D) crystallographic reconstruction. Here, we use HIV-1 capsid assemblies to illustrate the detailed approaches to obtain 3D density maps from helical objects. Mature HIV-1 particles contain a conical- or tubular-shaped capsid that encloses the viral RNA genome and performs essential functions in the virus life cycle. The capsid is composed of capsid protein (CA) oligomers which are helically arranged on the surface. The N-terminal domain (NTD) of CA is connected to its C-terminal domain (CTD) through a flexible hinge. Structural analysis of two- and three-dimensional crystals provided molecular models of the capsid protein (CA) and its oligomer forms. We determined the 3D density map of helically assembled HIV-1 CA hexamers at 16 Å resolution using an iterative helical real-space reconstruction method. Docking of atomic models of CA-NTD and CA-CTD dimer into the electron density map indicated that the CTD dimer interface is retained in the assembled CA. Furthermore, molecular docking revealed an additional, novel CTD trimer interface.

Applications of small-angle X-ray scattering to biomacromolecular solutions

Small-angle scattering of X-rays (SAXS) is an established method for low-resolution structural characterization of biological macromolecules in solution. Being complementary to the high resolution methods (X-ray crystallography and NMR), SAXS is often used in combination with them. The technique provides overall three-dimensional structures using ab initio reconstructions and hybrid modeling, and allows one to quantitatively characterize equilibrium mixtures as well as flexible systems. Recent progress in SAXS instrumentation, most notably, high brilliance synchrotron sources, has paved the way for high throughput automated SAXS studies allowing screening of external conditions (pH, temperature, ligand binding etc.). The modern approaches for SAXS data analysis are presented in this review including rapid characterization of macromolecular solutions in terms of low-resolution shapes, validation of high-resolution models in close-to-native conditions, quaternary structure analysis of complexes and quantitative description of the oligomeric composition in mixtures. Practical aspects of SAXS as a standalone tool and its combinations with other structural, biophysical or bioinformatics methods are reviewed. The capabilities of the technique are illustrated by a selection of recent applications for the studies of biological molecules. Future perspectives on SAXS and its potential impact to structural molecular biology are discussed.

ESCRT-III CHMP2A and CHMP3 form variable helical polymers in vitro and act synergistically during HIV-1 budding

The endosomal sorting complex required for transport-III (ESCRT-III) proteins are essential for budding of some enveloped viruses, for the formation of intraluminal vesicles at the endosome and for the abscission step of cytokinesis. ESCRT-III proteins form polymers that constrict membrane tubes, leading to fission. We have used electron cryomicroscopy to determine the molecular organization of pleiomorphic ESCRT-III CHMP2A-CHMP3 polymers. The three-dimensional reconstruction at 22 Å resolution reveals a helical organization of filaments of CHMP molecules organized in a head-to-tail fashion. Protease susceptibility experiments indicate that polymerization is achieved via conformational changes that increase the protomer stability. Combinatorial siRNA knockdown experiments indicate that CHMP3 contributes synergistically to HIV-1 budding, and the CHMP3 contribution is ~ 10-fold more pronounced in concert with CHMP2A than with CHMP2B. This is consistent with surface plasmon resonance affinity measurements that suggest sequential CHMP4B-CHMP3-CHMP2A recruitment while showing that both CHMP2A and CHMP2B interact with CHMP4B, in agreement with their redundant functions in HIV-1 budding. Our data thus indicate that the CHMP2A-CHMP3 polymer observed in vitro contributes to HIV-1 budding by assembling on CHMP4B polymers.

Cryo-EM structure of gastric H+,K+-ATPase with a single occupied cation-binding site

Gastric H(+),K(+)-ATPase is responsible for gastric acid secretion. ATP-driven H(+) uptake into the stomach is efficiently accomplished by the exchange of an equal amount of K(+), resulting in a luminal pH close to 1. Because of the limited free energy available for ATP hydrolysis, the stoichiometry of transported cations is thought to vary from 2H(+)/2K(+) to 1H(+)/1K(+) per hydrolysis of one ATP molecule as the luminal pH decreases, although direct evidence for this hypothesis has remained elusive. Here, we show, using the phosphate analog aluminum fluoride (AlF) and a K(+) congener (Rb(+)), the 8-Å resolution structure of H(+),K(+)-ATPase in the transition state of dephosphorylation, (Rb(+))E2~AlF, which is distinct from the preceding Rb(+)-free E2P state. A strong density located in the transmembrane cation-binding site of (Rb(+))E2~AlF highly likely represents a single bound Rb(+) ion, which is clearly different from the Rb(+)-free E2AlF or K(+)-bound (K(+))E2~AlF structures. Measurement of radioactive (86)Rb(+) binding suggests that the binding stoichiometry varies depending on the pH, and approximately half of the amount of Rb(+) is bound under acidic crystallization conditions compared with at a neutral pH. These data represent structural and biochemical evidence for the 1H(+)/1K(+)/1ATP transport mode of H(+),K(+)-ATPase, which is a prerequisite for generation of the 10(6)-fold proton gradient in terms of thermodynamics. Together with the released E2P-stabilizing interaction between the β subunit's N terminus and the P domain observed in the (Rb(+))E2~AlF structure, we propose a refined vectorial transport model of H(+),K(+)-ATPase, which must prevail against the highly acidic state of the gastric lumen.

Specificity of an anti-capsid antibody associated with Hepatitis B Virus-related acute liver failure

Previously, the livers of patients suffering from acute liver failure (ALF), a potentially fatal syndrome arising from infection by Hepatitis B Virus (HBV), were found to contain massive amounts of an antibody specific for the core antigen (HBcAg) capsid. We have used cryo-electron microscopy and molecular modeling to define its epitope. HBV capsids are icosahedral shells with 25Å-long dimeric spikes, each a 4-helix bundle, protruding from the contiguous "floor". Of the anti-HBcAg antibodies previously characterized, most bind around the spike tip while one binds to the floor. The ALF-associated antibody binds tangentially to a novel site on the side of the spike. This epitope is conformational. The Fab binds with high affinity to its principal determinants but has lower affinities for quasi-equivalent variants. The highest occupancy site is on one side of a spike, with no detectable binding to the corresponding site on the other side. Binding of one Fab per dimer was also observed by analytical ultracentrifugation. The Fab did not bind to the e-antigen dimer, a non-assembling variant of capsid protein. These findings support the propositions that antibodies with particular specificities may correlate with different clinical expressions of HBV infection and that antibodies directed to particular HBcAg epitopes may be involved in ALF pathogenesis.

Interaction of the HIV-1 intasome with transportin 3 protein (TNPO3 or TRN-SR2)

Transportin 3 (TNPO3 or TRN-SR2) has been shown to be an important cellular factor for early steps of lentiviral replication. However, separate studies have implicated distinct mechanisms for TNPO3 either through its interaction with HIV-1 integrase or capsid. Here we have carried out a detailed biophysical characterization of TNPO3 and investigated its interactions with viral proteins. Biophysical analyses including circular dichroism, analytical ultracentrifugation, small-angle x-ray scattering, and homology modeling provide insight into TNPO3 architecture and indicate that it is highly structured and exists in a monomer-dimer equilibrium in solution. In vitro biochemical binding assays argued against meaningful direct interaction between TNPO3 and the capsid cores. Instead, TNPO3 effectively bound to the functional intasome but not to naked viral DNA, suggesting that TNPO3 can directly engage the HIV-1 IN tetramer prebound to the cognate DNA. Mass spectrometry-based protein footprinting and site-directed mutagenesis studies have enabled us to map several interacting amino acids in the HIV-1 IN C-terminal domain and the cargo binding domain of TNPO3. Our findings provide important information for future genetic analysis to better understand the role of TNPO3 and its interacting partners for HIV-1 replication.

The hexameric structure of a conjugative VirB4 protein ATPase provides new insights for a functional and phylogenetic relationship with DNA translocases

VirB4 proteins are ATPases essential for pilus biogenesis and protein transport in type IV secretion systems. These proteins contain a motor domain that shares structural similarities with the motor domains of DNA translocases, such as the VirD4/TrwB conjugative coupling proteins and the chromosome segregation pump FtsK. Here, we report the three-dimensional structure of full-length TrwK, the VirB4 homologue in the conjugative plasmid R388, determined by single-particle electron microscopy. The structure consists of a hexameric double ring with a barrel-shaped structure. The C-terminal half of VirB4 proteins shares a striking structural similarity with the DNA translocase TrwB. Docking the atomic coordinates of the crystal structures of TrwB and FtsK into the EM map revealed a better fit for FtsK. Interestingly, we have found that like TrwB, TrwK is able to bind DNA with a higher affinity for G4 quadruplex structures than for single-stranded DNA. Furthermore, TrwK exerts a dominant negative effect on the ATPase activity of TrwB, which reflects an interaction between the two proteins. Our studies provide new insights into the structure-function relationship and the evolution of these DNA and protein translocases.

Structural determination of the phosphorylation domain of the ryanodine receptor

The ryanodine receptor (RyR) is a large, homotetrameric sarcoplasmic reticulum membrane protein that is essential for Ca(2+) cycling in both skeletal and cardiac muscle. Genetic mutations in RyR1 are associated with severe conditions including malignant hyperthermia (MH) and central core disease. One phosphorylation site (Ser 2843) has been identified in a segment of RyR1 flanked by two RyR motifs, which are found exclusively in all RyR isoforms as closely associated tandem (or paired) motifs, and are named after the protein itself. These motifs also contain six known MH mutations. In this study, we designed, expressed and purified the tandem RyR motifs, and show that this domain contains a putative binding site for the Ca(2+)/calmodulin-dependent protein kinase β isoform. We present a 2.2 Å resolution crystal structure of the RyR domain revealing a two-fold, symmetric, extended four-helix bundle stabilized by a β sheet. Using mathematical modelling, we fit our crystal structure within a tetrameric electron microscopy (EM) structure of native RyR1, and propose that this domain is localized in the RyR clamp region, which is absent in its cousin protein inositol 1,4,5-trisphosphate receptor.

Autocovariance Structures for Radial Averages in Small Angle X-Ray Scattering Experiments

Small-angle X-ray scattering (SAXS) is a technique for obtaining low-resolution structural information about biological macromolecules, by exposing a dilute solution to a high-intensity X-ray beam and capturing the resulting scattering pattern on a two-dimensional detector. The two-dimensional pattern is reduced to a one-dimensional curve through radial averaging; that is, by averaging across annuli on the detector plane. Subsequent analysis of structure relies on these one-dimensional data. This paper reviews the technique of SAXS and investigates autocorrelation structure in the detector plane and in the radial averages. Across a range of experimental conditions and molecular types, spatial autocorrelation in the detector plane is present and is well-described by a stationary kernel convolution model. The corresponding autocorrelation structure for the radial averages is non-stationary. Implications of the autocorrelation structure for inference about macromolecular structure are discussed.

Cryo-electron microscopy modeling by the molecular dynamics flexible fitting method

The increasing power and popularity of cryo-electron microscopy (cryo-EM) in structural biology brought about the development of so-called hybrid methods, which permit the interpretation of cryo-EM density maps beyond their nominal resolution in terms of atomic models. The Cryo-EM Modeling Challenge 2010 is the first community effort to bring together developers of hybrid methods as well as cryo-EM experimentalists. Participating in the challenge, the molecular dynamics flexible fitting (MDFF) method was applied to a number of cryo-EM density maps. The results are described here with special emphasis on the use of symmetry-based restraints to improve the quality of atomic models derived from density maps of symmetric complexes; on a comparison of the stereochemical quality of atomic models resulting from different hybrid methods; and on application of MDFF to electron crystallography data.

Real-space refinement with DireX: from global fitting to side-chain improvements

Single-particle cryo-electron microscopy (cryo-EM) has become an important tool to determine the structure of large biomolecules and assemblies thereof. However, the achievable resolution varies considerably over a wide range of about 3.5-20 Å. The interpretation of these intermediate- to low-resolution density maps in terms of atomic models is a big challenge and an area of active research. Here, we present our real-space structure refinement program DireX, which was developed primarily for cryo-EM-derived density maps. The basic principle and its main features are described. DireX employs Deformable Elastic Network (DEN) restraints to reduce overfitting by decreasing the effective number of degrees of freedom used in the refinement. Missing or reduced density due to flexible parts of the protein can lead to artifacts in the structure refinement, which is addressed through the concept of restrained grouped occupancy refinement. Furthermore, we describe the performance of DireX in the 2010 Cryo-EM Modeling Challenge, where we chose six density maps of four different proteins provided by the Modeling Challenge exemplifying typical refinement results at a large resolution range from 3 to 23 Å.

Macromolecular structure modeling from 3D EM using VolRover 2.0

We review tools for structure identification and model-based refinement from three-dimensional electron microscopy implemented in our in-house software package, VOLROVER 2.0. For viral density maps with icosahedral symmetry, we segment the capsid, polymeric, and monomeric subunits using techniques based on automatic symmetry detection and multidomain fast marching. For large biomolecules without symmetry information, we again use our multidomain fast-marching method with manual or fit-based multiseeding to segment meaningful substructures. In either case, we subject the resulting segmented subunit to secondary structure detection when the EM resolution is sufficiently high, and rigid-body structure fitting when the corresponding X-ray structure is available. Secondary structure elements are identified by three techniques: our earlier volume-based and boundary-based skeletonization methods as well as a new method, currently in development, based on solving the grassfire flow equation. For rigid-body fitting, we adapt our earlier fast Fourier-based correlation scheme F2Dock. Our reported segmentation, secondary structure elements identification, and rigid-body fitting techniques, implemented in VOLROVER 2.0 are applied to the PSB 2011 cryo-EM modeling challenge data, and our results are briefly compared to similar results submitted from other research groups. The comparisons show that our techniques are equally capable of segmenting relatively accurate subunits from a viral or protein assembly, and that high segmentation quality leads in turn to higher-quality results of secondary structure elements identification and correlation-based rigid-body fitting. © 2012 Wiley Periodicals, Inc. Biopolymers 97: 709-731, 2012.

Comparison of Segger and other methods for segmentation and rigid-body docking of molecular components in cryo-EM density maps

Segmentation and docking are useful methods for the discovery of molecular components in electron cryo-microscopy (cryo-EM) density maps of macromolecular complexes. In this article, we describe the segmentation and docking methods implemented in Segger. For 11 targets posted in the 2010 cryo-EM challenge, we segmented the regions corresponding to individual molecular components using Segger. We then used the segmented regions to guide rigid-body docking of individual components. Docking results were evaluated by comparing the docked components with published structures, and by calculation of several scores, such as atom inclusion, density occupancy, and geometry clash. The accuracy of the component segmentation using Segger and other methods was assessed by comparing segmented regions with docked components.

Structure of the RNA claw of the DNA packaging motor of bacteriophage Φ29

Bacteriophage DNA packaging motors translocate their genomic DNA into viral heads, compacting it to near-crystalline density. The Bacillus subtilis phage ϕ29 has a unique ring of RNA (pRNA) that is an essential component of its motor, serving as a scaffold for the packaging ATPase. Previously, deletion of a three-base bulge (18-CCA-20) in the pRNA A-helix was shown to abolish packaging activity. Here, we solved the structure of this crucial bulge by nuclear magnetic resonance (NMR) using a 27mer RNA fragment containing the bulge (27b). The bulge actually involves five nucleotides (17-UCCA-20 and A100), as U17 and A100 are not base paired as predicted. Mutational analysis showed these newly identified bulge residues are important for DNA packaging. The bulge introduces a 33–35° bend in the helical axis, and inter-helical motion around this bend appears to be restricted. A model of the functional 120b pRNA was generated using a 27b NMR structure and the crystal structure of the 66b prohead-binding domain. Fitting this model into a cryo-EM map generated a pentameric pRNA structure; five helices projecting from the pRNA ring resemble an RNA claw. Biochemical analysis suggested that this shape is important for coordinated motor action required for DNA translocation.

High precision alignment of cryo-electron subtomograms through gradient-based parallel optimization

BACKGROUND

Cryo-electron tomography emerges as an important component for structural system biology. It not only allows the structural characterization of macromolecular complexes, but also the detection of their cellular localizations in near living conditions. However, the method is hampered by low resolution, missing data and low signal-to-noise ratio (SNR). To overcome some of these difficulties and enhance the nominal resolution one can align and average a large set of subtomograms. Existing methods for obtaining the optimal alignments are mostly based on an exhaustive scanning of all but discrete relative rigid transformations (i.e. rotations and translations) of one subtomogram with respect to the other.

RESULTS

In this paper, we propose gradient-guided alignment methods based on two popular subtomogram similarity measures, a real space as well as a Fourier-space constrained score. We also propose a stochastic parallel refinement method that increases significantly the efficiency for the simultaneous refinement of a set of alignment candidates. We estimate that our stochastic parallel refinement is on average about 20 to 40 fold faster in comparison to the standard independent refinement approach. Results on simulated data of model complexes and experimental structures of protein complexes show that even for highly distorted subtomograms and with only a small number of very sparsely distributed initial alignment seeds, our combined methods can accurately recover true transformations with a substantially higher precision than the scanning based alignment methods.

CONCLUSIONS

Our methods increase significantly the efficiency and accuracy for subtomogram alignments, which is a key factor for the systematic classification of macromolecular complexes in cryo-electron tomograms of whole cells.

Structural determinants of the supramolecular organization of G protein-coupled receptors in bilayers

The G protein-coupled receptor (GPCR) rhodopsin self-assembles into supramolecular structures in native bilayers, but the structural determinants of receptor oligomerization are not known. We carried out multiple self-assembly coarse-grained molecular dynamics (CGMD) simulations of model membranes containing up to 64 molecules of the visual receptor rhodopsin over time scales reaching 100 μs. The simulations show strong preferential interaction modes between receptors. Two primary modes of receptor-receptor interactions are consistent with umbrella sampling/potential of mean force (PMF) calculations as a function of the distance between a pair of receptors. The preferential interfaces, involving helices (H) 1/8, 4/5 and 5, present no energy barrier to forming a very stable receptor dimer. Most notably, the PMFs show that the preferred rhodopsin dimer exists in a tail-to-tail conformation, with the interface comprising transmembrane H1/H2 and amphipathic H8 at the extracellular and cytoplasmic surfaces, respectively. This dimer orientation is in line with earlier electron microscopy, X-ray, and cross-linking experiments of rhodopsin and other GPCRs. Less stable interfaces, involving H4 and H6, have a free energy barrier for desolvation (delipidation) of the interfaces and appear to be designed to stabilize "lubricated" (i.e., lipid-coated) dimers. The overall CGMD strategy used here is general and can be applied to study the homo- and heterodimerization of GPCRs and other transmembrane proteins. Systematic extension of the work will deepen our understanding of the forces involved in the membrane organization of integral membrane proteins.

Fitting multimeric protein complexes into electron microscopy maps using 3D Zernike descriptors

A novel computational method for fitting high-resolution structures of multiple proteins into a cryoelectron microscopy map is presented. The method named EMLZerD generates a pool of candidate multiple protein docking conformations of component proteins, which are later compared with a provided electron microscopy (EM) density map to select the ones that fit well into the EM map. The comparison of docking conformations and the EM map is performed using the 3D Zernike descriptor (3DZD), a mathematical series expansion of three-dimensional functions. The 3DZD provides a unified representation of the surface shape of multimeric protein complex models and EM maps, which allows a convenient, fast quantitative comparison of the three-dimensional structural data. Out of 19 multimeric complexes tested, near native complex structures with a root-mean-square deviation of less than 2.5 Å were obtained for 14 cases while medium range resolution structures with correct topology were computed for the additional 5 cases.

Outcome of the first electron microscopy validation task force meeting

This Meeting Review describes the proceedings and conclusions from the inaugural meeting of the Electron Microscopy Validation Task Force organized by the Unified Data Resource for 3DEM (http://www.emdatabank.org) and held at Rutgers University in New Brunswick, NJ on September 28 and 29, 2010. At the workshop, a group of scientists involved in collecting electron microscopy data, using the data to determine three-dimensional electron microscopy (3DEM) density maps, and building molecular models into the maps explored how to assess maps, models, and other data that are deposited into the Electron Microscopy Data Bank and Protein Data Bank public data archives. The specific recommendations resulting from the workshop aim to increase the impact of 3DEM in biology and medicine.

Structural and functional insight into how the Plasmodium falciparum VAR2CSA protein mediates binding to chondroitin sulfate A in placental malaria

Malaria is a major global health problem. Pregnant women are susceptible to infection regardless of previously acquired immunity. Placental malaria is caused by parasites capable of sequestering in the placenta. This is mediated by VAR2CSA, a parasite antigen that interacts with chondroitin sulfate A (CSA). One vaccine strategy is to block this interaction with VAR2CSA-specific antibodies. It is a priority to define a small VAR2CSA fragment that can be used in an adhesion blocking vaccine. In this, the obvious approach is to define regions of VAR2CSA involved in receptor binding. It has been shown that full-length recombinant VAR2CSA binds specifically to CSA with nanomolar affinity, and that the CSA-binding site lies in the N-terminal part of the protein. In this study we define the minimal binding region by truncating VAR2CSA and analyzing CSA binding using biosensor technology. We show that the core CSA-binding site lies within the DBL2X domain and parts of the flanking interdomain regions. This is in contrast to the idea that single domains do not possess the structural requirements for specific CSA binding. Small-angle x-ray scattering measurements enabled modeling of VAR2CSA and showed that the CSA-binding DBL2X domain is situated in the center of the structure. Mutating classic sulfate-binding sites in VAR2CSA, along with testing dependence of ionic interactions, suggest that the CSA binding is not solely dependent on the sulfated CSA structure. Based on these novel PfEMP1 structure-function studies, we have constructed a small VAR2CSA antigen that has the capacity to induce highly adhesion-blocking antibodies.