Molecular dynamics flexible fitting: a practical guide to combine cryo-electron microscopy and X-ray crystallography

Quantitative Characterization of Domain Motions in Molecular Machines

The work of molecular machines such as the ribosome is accompanied by conformational changes, often characterized by relative motions of their domains. The method we have developed seeks to quantify these motions in a general way, facilitating comparisons of results obtained by different researchers. Typically there are multiple snapshots of a structure in the form of pdb coordinates resulting from flexible fitting of low-resolution density maps, from X-ray crystallography, or from molecular dynamics simulation trajectories. Our objective is to characterize the motion of each domain as a coordinate transformation using moments of inertia tensor, a method we developed earlier. What has been missing until now are ancillary tools that make this task practical, general, and biologically informative. We have provided a comprehensive solution to this task with a set of tools implemented on the VMD platform. These tools address the need for reproducible segmentation of domains, and provide a generalized description of their motions using principal axes of inertia. Although this methodology has been specifically developed for studying ribosome motion, it is applicable to any molecular machine.

Beta-Propiolactone Inactivation of Coxsackievirus A16 Induces Structural Alteration and Surface Modification of Viral Capsids

Beta-propiolactone (BPL) is an inactivating agent that is widely used in the vaccine industry. However, its effects on vaccine protein antigens and its mechanisms of action remain poorly understood. Here we present cryo-electron microscopy (cryo-EM) structures of BPL-treated coxsackievirus A16 (CVA16) mature virions and procapsids at resolutions of 3.9 Å and 6.5 Å, respectively. Notably, both particles were found to adopt an expanded conformation resembling the 135S-like uncoating intermediate, with characteristic features including an opened 2-fold channel, the externalization of the N terminus of VP1 capsid protein, and the absence of pocket factor. However, major neutralizing epitopes are very well preserved on these particles. Further biochemical analyses revealed that BPL treatment impairs the abilities of CVA16 particles to bind to the attachment receptor heparan sulfate and to a conformation-dependent monoclonal antibody in a BPL dose-dependent manner, indicating that BPL is able to modify surface-exposed amino acid residues. Taken together, our results demonstrate that BPL treatment may induce alteration of the overall structure and surface properties of a nonenveloped viral capsid, thus revealing a novel mode of action of BPL.IMPORTANCE Beta-propiolactone (BPL) is commonly used as an inactivating reagent to produce viral vaccines. It is recognized that BPL inactivates viral infectivity through modification of viral nucleic acids. However, its effect on viral proteins remains largely unknown. Here, we present high-resolution cryo-EM structures of BPL-treated coxsackievirus A16 (CVA16) mature virions and procapsids, which reveals an expanded overall conformation and characteristic features that are typical for the 135S-like uncoating intermediate. We further show that the BPL concentration affects the binding of inactivated CVA16 particles to their receptor/antibody. Thus, BPL treatment can alter the overall structure and surface properties of viral capsids, which may lead to antigenic and immunogenic variations. Our findings provide important information for future development of BPL-inactivated vaccines.

Cryo-EM structures of the triheteromeric NMDA receptor and its allosteric modulation

N-methyl-d-aspartate receptors (NMDARs) are heterotetrameric ion channels assembled as diheteromeric or triheteromeric complexes. Here, we report structures of the triheteromeric GluN1/GluN2A/GluN2B receptor in the absence or presence of the GluN2B-specific allosteric modulator Ro 25-6981 (Ro), determined by cryogenic electron microscopy (cryo-EM). In the absence of Ro, the GluN2A and GluN2B amino-terminal domains (ATDs) adopt "closed" and "open" clefts, respectively. Upon binding Ro, the GluN2B ATD clamshell transitions from an open to a closed conformation. Consistent with a predominance of the GluN2A subunit in ion channel gating, the GluN2A subunit interacts more extensively with GluN1 subunits throughout the receptor, in comparison with the GluN2B subunit. Differences in the conformation of the pseudo-2-fold-related GluN1 subunits further reflect receptor asymmetry. The triheteromeric NMDAR structures provide the first view of the most common NMDA receptor assembly and show how incorporation of two different GluN2 subunits modifies receptor symmetry and subunit interactions, allowing each subunit to uniquely influence receptor structure and function, thus increasing receptor complexity.

The First Extracellular Linker Is Important for Several Aspects of the Gating Mechanism of Human TRPA1 Channel

Transient receptor potential ankyrin 1 (TRPA1) is an excitatory ion channel involved in pain, inflammation and itching. This channel gates in response to many irritant and proalgesic agents, and can be modulated by calcium and depolarizing voltage. While the closed-state structure of TRPA1 has been recently resolved, also having its open state is essential for understanding how this channel works. Here we use molecular dynamics simulations combined with electrophysiological measurements and systematic mutagenesis to predict and explore the conformational changes coupled to the expansion of the presumptive channel's lower gate. We show that, upon opening, the upper part of the sensor module approaches the pore domain of an adjacent subunit and the conformational dynamics of the first extracellular flexible loop may govern the voltage-dependence of multimodal gating, thereby serving to stabilize the open state of the channel. These results are generally important in understanding the structure and function of TRPA1 and offer new insights into the gating mechanism of TRPA1 and related channels.

Structural insights into the functional cycle of the ATPase module of the 26S proteasome

In eukaryotic cells, the ubiquitin-proteasome system (UPS) is responsible for the regulated degradation of intracellular proteins. The 26S holocomplex comprises the core particle (CP), where proteolysis takes place, and one or two regulatory particles (RPs). The base of the RP is formed by a heterohexameric AAA⁺ ATPase module, which unfolds and translocates substrates into the CP. Applying single-particle cryo-electron microscopy (cryo-EM) and image classification to samples in the presence of different nucleotides and nucleotide analogs, we were able to observe four distinct conformational states (s1 to s4). The resolution of the four conformers allowed for the construction of atomic models of the AAA⁺ ATPase module as it progresses through the functional cycle. In a hitherto unobserved state (s4), the gate controlling access to the CP is open. The structures described in this study allow us to put forward a model for the 26S functional cycle driven by ATP hydrolysis.

Antibodyomics: bioinformatics technologies for understanding B-cell immunity to HIV-1

Numerous antibodies have been identified from HIV-1-infected donors that neutralize diverse strains of HIV-1. These antibodies may provide the basis for a B cell-mediated HIV-1 vaccine. However, it has been unclear how to elicit similar antibodies by vaccination. To address this issue, we have undertaken an informatics-based approach to understand the genetic and immunologic processes controlling the development of HIV-1-neutralizing antibodies. As DNA sequencing comprises the fastest growing database of biological information, we focused on incorporating next-generation sequencing of B-cell transcripts to determine the origin, maturation pathway, and prevalence of broadly neutralizing antibody lineages (Antibodyomics1, 2, 4, and 6). We also incorporated large-scale robotic analyses of serum neutralization to identify and quantify neutralizing antibodies in donor cohorts (Antibodyomics3). Statistical analyses furnish another layer of insight (Antibodyomics5), with physical characteristics of antibodies and their targets through molecular dynamics simulations (Antibodyomics7) and free energy perturbation analyses (Antibodyomics8) providing information-rich output. Functional interrogation of individual antibodies (Antibodyomics9) and synthetic antibody libraries (Antibodyomics10) also yields multi-dimensional data by which to understand and improve antibodies. Antibodyomics, described here, thus comprise resolution-enhancing tools, which collectively embody an information-driven discovery engine aimed toward the development of effective B cell-based vaccines.

CHARMM-GUI MDFF/xMDFF Utilizer for Molecular Dynamics Flexible Fitting Simulations in Various Environments

X-ray crystallography and cryo-electron microscopy are two popular methods for the structure determination of biological molecules. Atomic structures are derived through the fitting and refinement of an initial model into electron density maps constructed by both experiments. Two computational approaches, MDFF and xMDFF, have been developed to facilitate this process by integrating the experimental data with molecular dynamics simulation. However, the setup of an MDFF/xMDFF simulation requires knowledge of both experimental and computational methods, which is not straightforward for nonexpert users. In addition, sometimes it is desirable to include realistic environments, such as explicit solvent and lipid bilayers during the simulation, which poses another challenge even for expert users. To alleviate these difficulties, we have developed MDFF/xMDFF Utilizer in CHARMM-GUI that helps users to set up an MDFF/xMDFF simulation. The capability of MDFF/xMDFF Utilizer is greatly enhanced by integration with other CHARMM-GUI modules, including protein structure manipulation, a diverse set of lipid types, and all-atom CHARMM and coarse-grained PACE force fields. With this integration, various simulation environments are available for MDFF Utilizer (vacuum, implicit/explicit solvent, and bilayers) and xMDFF Utilizer (vacuum and solution). In this work, three examples are shown to demonstrate the usage of MDFF/xMDFF Utilizer.

Insights into Adenovirus Uncoating from Interactions with Integrins and Mediators of Host Immunity

Human adenoviruses are large (150 MDa) nonenveloped double-stranded DNA (dsDNA) viruses that cause acute respiratory, gastrointestinal and ocular infections. Despite these disease associations, adenovirus has aided basic and clinical research efforts through studies of its association with cells and as a target of host antiviral responses. This review highlights the knowledge of adenovirus disassembly and nuclear transport gleaned from structural, biophysical and functional analyses of adenovirus interactions with soluble and membrane-associated host molecules.

Asymmetric cryo-EM structure of the canonical Allolevivirus Qβ reveals a single maturation protein and the genomic ssRNA in situ

Single-stranded (ss) RNA viruses infect all domains of life. To date, for most ssRNA virions, only the structures of the capsids and their associated protein components have been resolved to high resolution. Qβ, an ssRNA phage specific for the conjugative F-pilus, has a T = 3 icosahedral lattice of coat proteins assembled around its 4,217 nucleotides of genomic RNA (gRNA). In the mature virion, the maturation protein, A2, binds to the gRNA and is required for adsorption to the F-pilus. Here, we report the cryo-electron microscopy (cryo-EM) structures of Qβ with and without symmetry applied. The icosahedral structure, at 3.7-Å resolution, resolves loops not previously seen in the published X-ray structure, whereas the asymmetric structure, at 7-Å resolution, reveals A2 and the gRNA. A2 contains a bundle of α-helices and replaces one dimer of coat proteins at a twofold axis. The helix bundle binds gRNA, causing denser packing of RNA in its proximity, which asymmetrically expands the surrounding coat protein shell to potentially facilitate RNA release during infection. We observe a fixed pattern of gRNA organization among all viral particles, with the major and minor grooves of RNA helices clearly visible. A single layer of RNA directly contacts every copy of the coat protein, with one-third of the interactions occurring at operator-like RNA hairpins. These RNA-coat interactions stabilize the tertiary structure of gRNA within the virion, which could further provide a roadmap for capsid assembly.

Architecture of the Human Mitochondrial Iron-Sulfur Cluster Assembly Machinery

Fe-S clusters, essential cofactors needed for the activity of many different enzymes, are assembled by conserved protein machineries inside bacteria and mitochondria. As the architecture of the human machinery remains undefined, we co-expressed in Escherichia coli the following four proteins involved in the initial step of Fe-S cluster synthesis: FXN42-210 (iron donor); [NFS1]·[ISD11] (sulfur donor); and ISCU (scaffold upon which new clusters are assembled). We purified a stable, active complex consisting of all four proteins with 1:1:1:1 stoichiometry. Using negative staining transmission EM and single particle analysis, we obtained a three-dimensional model of the complex with ∼14 Å resolution. Molecular dynamics flexible fitting of protein structures docked into the EM map of the model revealed a [FXN42-210]24·[NFS1]24·[ISD11]24·[ISCU]24 complex, consistent with the measured 1:1:1:1 stoichiometry of its four components. The complex structure fulfills distance constraints obtained from chemical cross-linking of the complex at multiple recurring interfaces, involving hydrogen bonds, salt bridges, or hydrophobic interactions between conserved residues. The complex consists of a central roughly cubic [FXN42-210]24·[ISCU]24 sub-complex with one symmetric ISCU trimer bound on top of one symmetric FXN42-210 trimer at each of its eight vertices. Binding of 12 [NFS1]2·[ISD11]2 sub-complexes to the surface results in a globular macromolecule with a diameter of ∼15 nm and creates 24 Fe-S cluster assembly centers. The organization of each center recapitulates a previously proposed conserved mechanism for sulfur donation from NFS1 to ISCU and reveals, for the first time, a path for iron donation from FXN42-210 to ISCU.

Chemotaxis cluster 1 proteins form cytoplasmic arrays in Vibrio cholerae and are stabilized by a double signaling domain receptor DosM

Nearly all motile bacterial cells use a highly sensitive and adaptable sensory system to detect changes in nutrient concentrations in the environment and guide their movements toward attractants and away from repellents. The best-studied bacterial chemoreceptor arrays are membrane-bound. Many motile bacteria contain one or more additional, sometimes purely cytoplasmic, chemoreceptor systems. Vibrio cholerae contains three chemotaxis clusters (I, II, and III). Here, using electron cryotomography, we explore V. cholerae's cytoplasmic chemoreceptor array and establish that it is formed by proteins from cluster I. We further identify a chemoreceptor with an unusual domain architecture, DosM, which is essential for formation of the cytoplasmic arrays. DosM contains two signaling domains and spans the two-layered cytoplasmic arrays. Finally, we present evidence suggesting that this type of receptor is important for the structural stability of the cytoplasmic array.

Structural and mechanistic insights into human splicing factor SF3b complex derived using an integrated approach guided by the cryo-EM density maps

Pre-mRNA splicing in eukaryotes is performed by the spliceosome, a highly complex macromolecular machine. SF3b is a multi-protein complex which recognizes the branch point adenosine of pre-mRNA as part of a larger U2 snRNP or U11/U12 di-snRNP in the dynamic spliceosome machinery. Although a cryo-EM map is available for human SF3b complex, the structure and relative spatial arrangement of all components in the complex are not yet known. We have recognized folds of domains in various proteins in the assembly and generated comparative models. Using an integrative approach involving structural and other experimental data, guided by the available cryo-EM density map, we deciphered a pseudo-atomic model of the closed form of SF3b which is found to be a "fuzzy complex" with highly flexible components and multiplicity of folds. Further, the model provides structural information for 5 proteins (SF3b10, SF3b155, SF3b145, SF3b130 and SF3b14b) and localization information for 4 proteins (SF3b10, SF3b145, SF3b130 and SF3b14b) in the assembly for the first time. Integration of this model with the available U11/U12 di-snRNP cryo-EM map enabled elucidation of an open form. This now provides new insights on the mechanistic features involved in the transition between closed and open forms pivoted by a hinge region in the SF3b155 protein that also harbors cancer causing mutations. Moreover, the open form guided model of the 5' end of U12 snRNA, which includes the branch point duplex, shows that the architecture of SF3b acts as a scaffold for U12 snRNA: pre-mRNA branch point duplex formation with potential implications for branch point adenosine recognition fidelity.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

Molecular dynamics-based refinement and validation for sub-5 Å cryo-electron microscopy maps

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.

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

To understand the roles that proteins and other large molecules play inside cells, it is important to determine their structures. One of the techniques that researchers can use to do this is called cryo-electron microscopy (cryo-EM), which rapidly freezes molecules to fix them in position before imaging them in fine detail.

The cryo-EM images are like maps that show the approximate position of atoms. These images must then be processed in order to build a three-dimensional model of the protein that shows how its atoms are arranged relative to each other. One computational approach called Molecular Dynamics Flexible Fitting (MDFF) works by flexibly fitting possible atomic structures into cryo-EM maps. Although this approach works well with relatively undetailed (or ‘low resolution’) cryo-EM images, it struggles to handle the high-resolution cryo-EM maps now being generated.

Singharoy, Teo, McGreevy et al. have now developed two MDFF methods – called cascade MDFF and resolution exchange MDFF – that help to resolve atomic models of biological molecules from cryo-EM images. Each method can refine poorly guessed models into ones that are consistent with the high-resolution experimental images. The refinement is achieved by interpreting a range of images that starts with a ‘fuzzy’ image. The contrast of the image is then progressively improved until an image is produced that has a resolution that is good enough to almost distinguish individual atoms.

The method works because each cryo-EM image shows not just one, but a collection of atomic structures that the molecule can take on, with the fuzzier parts of the image representing the more flexible parts of the molecule. By taking into account this flexibility, the large-scale features of the protein structure can be determined first from the fuzzier images, and increasing the contrast of the images allows smaller-scale refinements to be made to the structure.

The MDFF tools have been designed to be easy to use and are available to researchers at low cost through cloud computing platforms. They can now be used to unravel the structure of many different proteins and protein complexes including those involved in photosynthesis, respiration and protein synthesis.

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

Improved cryoEM-Guided Iterative Molecular Dynamics--Rosetta Protein Structure Refinement Protocol for High Precision Protein Structure Prediction

Many excellent methods exist that incorporate cryo-electron microscopy (cryoEM) data to constrain computational protein structure prediction and refinement. Previously, it was shown that iteration of two such orthogonal sampling and scoring methods – Rosetta and molecular dynamics (MD) simulations – facilitated exploration of conformational space in principle. Here, we go beyond a proof-of-concept study and address significant remaining limitations of the iterative MD–Rosetta protein structure refinement protocol. Specifically, all parts of the iterative refinement protocol are now guided by medium-resolution cryoEM density maps, and previous knowledge about the native structure of the protein is no longer necessary. Models are identified solely based on score or simulation time. All four benchmark proteins showed substantial improvement through three rounds of the iterative refinement protocol. The best-scoring final models of two proteins had sub-Ångstrom RMSD to the native structure over residues in secondary structure elements. Molecular dynamics was most efficient in refining secondary structure elements and was thus highly complementary to the Rosetta refinement which is most powerful in refining side chains and loop regions.

Structure of the human 26S proteasome at a resolution of 3.9 Å

Protein degradation in eukaryotic cells is performed by the Ubiquitin-Proteasome System (UPS). The 26S proteasome holocomplex consists of a core particle (CP) that proteolytically degrades polyubiquitylated proteins, and a regulatory particle (RP) containing the AAA-ATPase module. This module controls access to the proteolytic chamber inside the CP and is surrounded by non-ATPase subunits (Rpns) that recognize substrates and deubiquitylate them before unfolding and degradation. The architecture of the 26S holocomplex is highly conserved between yeast and humans. The structure of the human 26S holocomplex described here reveals previously unidentified features of the AAA-ATPase heterohexamer. One subunit, Rpt6, has ADP bound, whereas the other five have ATP in their binding pockets. Rpt6 is structurally distinct from the other five Rpt subunits, most notably in its pore loop region. For Rpns, the map reveals two main, previously undetected, features: the C terminus of Rpn3 protrudes into the mouth of the ATPase ring; and Rpn1 and Rpn2, the largest proteasome subunits, are linked by an extended connection. The structural features of the 26S proteasome observed in this study are likely to be important for coordinating the proteasomal subunits during substrate processing.

The bacteriophage ϕ29 tail possesses a pore-forming loop for cell membrane penetration

Most bacteriophages are tailed bacteriophages with an isometric or a prolate head attached to a long contractile, long non-contractile, or short non-contractile tail. The tail is a complex machine that plays a central role in host cell recognition and attachment, cell wall and membrane penetration, and viral genome ejection. The mechanisms involved in the penetration of the inner host cell membrane by bacteriophage tails are not well understood. Here we describe structural and functional studies of the bacteriophage ϕ29 tail knob protein gene product 9 (gp9). The 2.0 Å crystal structure of gp9 shows that six gp9 molecules form a hexameric tube structure with six flexible hydrophobic loops blocking one end of the tube before DNA ejection. Sequence and structural analyses suggest that the loops in the tube could be membrane active. Further biochemical assays and electron microscopy structural analyses show that the six hydrophobic loops in the tube exit upon DNA ejection and form a channel that spans the lipid bilayer of the membrane and allows the release of the bacteriophage genomic DNA, suggesting that cell membrane penetration involves a pore-forming mechanism similar to that of certain non-enveloped eukaryotic viruses. A search of other phage tail proteins identified similar hydrophobic loops, which indicates that a common mechanism might be used for membrane penetration by prokaryotic viruses. These findings suggest that although prokaryotic and eukaryotic viruses use apparently very different mechanisms for infection, they have evolved similar mechanisms for breaching the cell membrane.

Structure of a group II intron in complex with its reverse transcriptase

Bacterial group II introns are large catalytic RNAs related to nuclear spliceosomal introns and eukaryotic retrotransposons. They self-splice, yielding mature RNA, and integrate into DNA as retroelements. A fully active group II intron forms a ribonucleoprotein complex comprising the intron ribozyme and an intron-encoded protein that performs multiple activities including reverse transcription, in which intron RNA is copied into the DNA target. Here we report cryo-EM structures of an endogenously spliced Lactococcus lactis group IIA intron in its ribonucleoprotein complex form at 3.8-Å resolution and in its protein-depleted form at 4.5-Å resolution, revealing functional coordination of the intron RNA with the protein. Remarkably, the protein structure reveals a close relationship between the reverse transcriptase catalytic domain and telomerase, whereas the active splicing center resembles the spliceosomal Prp8 protein. These extraordinary similarities hint at intricate ancestral relationships and provide new insights into splicing and retromobility.

Structure of the intact ATM/Tel1 kinase

The ataxia-telangiectasia mutated (ATM) protein is an apical kinase that orchestrates the multifaceted DNA-damage response. Normally, ATM kinase is in an inactive, homodimer form and is transformed into monomers upon activation. Besides a conserved kinase domain at the C terminus, ATM contains three other structural modules, referred to as FAT, FATC and N-terminal helical solenoid. Here we report the first cryo-EM structure of ATM kinase, which is an intact homodimeric ATM/Tel1 from Schizosaccharomyces pombe. We show that two monomers directly contact head-to-head through the FAT and kinase domains. The tandem N-terminal helical solenoid tightly packs against the FAT and kinase domains. The structure suggests that ATM/Tel1 dimer interface and the consecutive HEAT repeats inhibit the binding of kinase substrates and regulators by steric hindrance. Our study provides a structural framework for understanding the mechanisms of ATM/Tel1 regulation as well as the development of new therapeutic agents.

QwikMD - Integrative Molecular Dynamics Toolkit for Novices and Experts

The proper functioning of biomolecules in living cells requires them to assume particular structures and to undergo conformational changes. Both biomolecular structure and motion can be studied using a wide variety of techniques, but none offers the level of detail as do molecular dynamics (MD) simulations. Integrating two widely used modeling programs, namely NAMD and VMD, we have created a robust, user-friendly software, QwikMD, which enables novices and experts alike to address biomedically relevant questions, where often only molecular dynamics simulations can provide answers. Performing both simple and advanced MD simulations interactively, QwikMD automates as many steps as necessary for preparing, carrying out, and analyzing simulations while checking for common errors and enabling reproducibility. QwikMD meets also the needs of experts in the field, increasing the efficiency and quality of their work by carrying out tedious or repetitive tasks while enabling easy control of every step. Whether carrying out simulations within the live view mode on a small laptop or performing complex and large simulations on supercomputers or Cloud computers, QwikMD uses the same steps and user interface. QwikMD is freely available by download on group and personal computers. It is also available on the cloud at Amazon Web Services.

Architecture of the Yeast Mitochondrial Iron-Sulfur Cluster Assembly Machinery: THE SUB-COMPLEX FORMED BY THE IRON DONOR, Yfh1 PROTEIN, AND THE SCAFFOLD, Isu1 PROTEIN

The biosynthesis of Fe-S clusters is a vital process involving the delivery of elemental iron and sulfur to scaffold proteins via molecular interactions that are still poorly defined. We reconstituted a stable, functional complex consisting of the iron donor, Yfh1 (yeast frataxin homologue 1), and the Fe-S cluster scaffold, Isu1, with 1:1 stoichiometry, [Yfh1]24·[Isu1]24 Using negative staining transmission EM and single particle analysis, we obtained a three-dimensional reconstruction of this complex at a resolution of ∼17 Å. In addition, via chemical cross-linking, limited proteolysis, and mass spectrometry, we identified protein-protein interaction surfaces within the complex. The data together reveal that [Yfh1]24·[Isu1]24 is a roughly cubic macromolecule consisting of one symmetric Isu1 trimer binding on top of one symmetric Yfh1 trimer at each of its eight vertices. Furthermore, molecular modeling suggests that two subunits of the cysteine desulfurase, Nfs1, may bind symmetrically on top of two adjacent Isu1 trimers in a manner that creates two putative [2Fe-2S] cluster assembly centers. In each center, conserved amino acids known to be involved in sulfur and iron donation by Nfs1 and Yfh1, respectively, are in close proximity to the Fe-S cluster-coordinating residues of Isu1. We suggest that this architecture is suitable to ensure concerted and protected transfer of potentially toxic iron and sulfur atoms to Isu1 during Fe-S cluster assembly.

Secondary Interaction Interfaces with PCNA Control Conformational Switching of DNA Polymerase PolB from Polymerization to Editing

Replicative DNA polymerases (Pols) frequently possess two distinct DNA processing activities: DNA synthesis (polymerization) and proofreading (3'-5' exonuclease activity). The polymerase and exonuclease reactions are performed alternately and are spatially separated in different protein domains. Thus, the growing DNA primer terminus has to undergo dynamic conformational switching between two distinct functional sites on the polymerase. Furthermore, the transition from polymerization (pol) mode to exonuclease (exo) mode must occur in the context of a DNA Pol holoenzyme, wherein the polymerase is physically associated with processivity factor proliferating cell nuclear antigen (PCNA) and primer-template DNA. The mechanism of this conformational switching and the role that PCNA plays in it have remained obscure, largely due to the dynamic nature of ternary Pol/PCNA/DNA assemblies. Here, we present computational models of ternary assemblies for archaeal polymerase PolB. We have combined all available structural information for the binary complexes with electron microscopy data and have refined atomistic models for ternary PolB/PCNA/DNA assemblies in pol and exo modes using molecular dynamics simulations. In addition to the canonical PIP-box/interdomain connector loop (IDCL) interface of PolB with PCNA, contact analysis of the simulation trajectories revealed new secondary binding interfaces, distinct between the pol and exo states. Using targeted molecular dynamics, we explored the conformational transition from pol to exo mode. We identified a hinge region between the thumb and palm domains of PolB that is critical for conformational switching. With the thumb domain anchored onto the PCNA surface, the neighboring palm domain executed rotational motion around the hinge, bringing the core of PolB down toward PCNA to form a new interface with the clamp. A helix from PolB containing a patch of arginine residues was involved in the binding, locking the complex in the exo mode conformation. Together, these results provide a structural view of how the transition between the pol and exo states of PolB is coordinated through PCNA to achieve efficient proofreading.

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.

Cyclophilin A stabilizes the HIV-1 capsid through a novel non-canonical binding site

The host cell factor cyclophilin A (CypA) interacts directly with the HIV-1 capsid and regulates viral infectivity. Although the crystal structure of CypA in complex with the N-terminal domain of the HIV-1 capsid protein (CA) has been known for nearly two decades, how CypA interacts with the viral capsid and modulates HIV-1 infectivity remains unclear. We determined the cryoEM structure of CypA in complex with the assembled HIV-1 capsid at 8-Å resolution. The structure exhibits a distinct CypA-binding pattern in which CypA selectively bridges the two CA hexamers along the direction of highest curvature. EM-guided all-atom molecular dynamics simulations and solid-state NMR further reveal that the CypA-binding pattern is achieved by single-CypA molecules simultaneously interacting with two CA subunits, in different hexamers, through a previously uncharacterized non-canonical interface. These results provide new insights into how CypA stabilizes the HIV-1 capsid and is recruited to facilitate HIV-1 infection.

The host cell factor cyclophilin A (CypA) interacts with the HIV-1 capsid and regulates infectivity. Here the authors combine cryo-EM, solid-state NMR and all-atom MD simulations, identifying an interaction interface between CypA and the HIV capsid that stabilizes the viral capsid and regulates infectivity.

Structural Basis for Recognition of Human Enterovirus 71 by a Bivalent Broadly Neutralizing Monoclonal Antibody

Enterovirus 71 (EV71) is the main pathogen responsible for hand, foot and mouth disease with severe neurological complications and even death in young children. We have recently identified a highly potent anti-EV71 neutralizing monoclonal antibody, termed D5. Here we investigated the structural basis for recognition of EV71 by the antibody D5. Four three-dimensional structures of EV71 particles in complex with IgG or Fab of D5 were reconstructed by cryo-electron microscopy (cryo-EM) single particle analysis all at subnanometer resolutions. The most critical EV71 mature virion-Fab structure was resolved to a resolution of 4.8 Å, which is rare in cryo-EM studies of virus-antibody complex so far. The structures reveal a bivalent binding pattern of D5 antibody across the icosahedral 2-fold axis on mature virion, suggesting that D5 binding may rigidify virions to prevent their conformational changes required for subsequent RNA release. Moreover, we also identified that the complementary determining region 3 (CDR3) of D5 heavy chain directly interacts with the extremely conserved VP1 GH-loop of EV71, which was validated by biochemical and virological assays. We further showed that D5 is indeed able to neutralize a variety of EV71 genotypes and strains. Moreover, D5 could potently confer protection in a mouse model of EV71 infection. Since the conserved VP1 GH-loop is involved in EV71 binding with its uncoating receptor, the scavenger receptor class B, member 2 (SCARB2), the broadly neutralizing ability of D5 might attribute to its inhibition of EV71 from binding SCARB2. Altogether, our results elucidate the structural basis for the binding and neutralization of EV71 by the broadly neutralizing antibody D5, thereby enhancing our understanding of antibody-based protection against EV71 infection.

Structural insights into the mechanism of activation of the TRPV1 channel by a membrane-bound tarantula toxin

Venom toxins are invaluable tools for exploring the structure and mechanisms of ion channels. Here, we solve the structure of double-knot toxin (DkTx), a tarantula toxin that activates the heat-activated TRPV1 channel. We also provide improved structures of TRPV1 with and without the toxin bound, and investigate the interactions of DkTx with the channel and membranes. We find that DkTx binds to the outer edge of the external pore of TRPV1 in a counterclockwise configuration, using a limited protein-protein interface and inserting hydrophobic residues into the bilayer. We also show that DkTx partitions naturally into membranes, with the two lobes exhibiting opposing energetics for membrane partitioning and channel activation. Finally, we find that the toxin disrupts a cluster of hydrophobic residues behind the selectivity filter that are critical for channel activation. Collectively, our findings reveal a novel mode of toxin-channel recognition that has important implications for the mechanism of thermosensation.

3DIANA: 3D Domain Interaction Analysis: A Toolbox for Quaternary Structure Modeling

Electron microscopy (EM) is experiencing a revolution with the advent of a new generation of Direct Electron Detectors, enabling a broad range of large and flexible structures to be resolved well below 1 nm resolution. Although EM techniques are evolving to the point of directly obtaining structural data at near-atomic resolution, for many molecules the attainable resolution might not be enough to propose high-resolution structural models. However, accessing information on atomic coordinates is a necessary step toward a deeper understanding of the molecular mechanisms that allow proteins to perform specific tasks. For that reason, methods for the integration of EM three-dimensional maps with x-ray and NMR structural data are being developed, a modeling task that is normally referred to as fitting, resulting in the so called hybrid models. In this work, we present a novel application—3DIANA—specially targeted to those cases in which the EM map resolution is medium or low and additional experimental structural information is scarce or even lacking. In this way, 3DIANA statistically evaluates proposed/potential contacts between protein domains, presents a complete catalog of both structurally resolved and predicted interacting regions involving these domains and, finally, suggests structural templates to model the interaction between them. The evaluation of the proposed interactions is computed with DIMERO, a new method that scores physical binding sites based on the topology of protein interaction networks, which has recently shown the capability to increase by 200% the number of domain-domain interactions predicted in interactomes as compared to previous approaches. The new application displays the information at a sequence and structural level and is accessible through a web browser or as a Chimera plugin at http://3diana.cnb.csic.es.

Nascent peptide assists the ribosome in recognizing chemically distinct small molecules

Regulation of gene expression in response to the changing environment is critical for cell survival. For instance, binding of macrolide antibiotics to the ribosome promote the translation arrest at the leader ORFs ermCL and ermBL necessary for inducing antibiotic resistance genes ermC and ermB. Cladinose-containing macrolides, like erythromycin (ERY), but not ketolides e.g., telithromycin (TEL), arrest translation of ermCL, while either ERY or TEL stall ermBL translation. How the ribosome distinguishes between chemically similar small molecules is unknown.

We show that single amino acid changes in the leader peptide switch the specificity of recognition of distinct molecules, triggering gene activation in response to only ERY, only TEL, to both antibiotics, or preventing stalling altogether. Thus, the ribosomal response to chemical signals can be modulated by minute changes in the nascent peptide, suggesting that protein sequences could have been optimized for rendering translation sensitive to environmental cues.

Single-particle cryo-electron microscopy of macromolecular complexes

Recent technological breakthroughs in image acquisition have enabled single-particle cryo-electron microscopy (cryo-EM) to achieve near-atomic resolution structural information for biological complexes. The improvements in image quality coupled with powerful computational methods for sorting distinct particle populations now also allow the determination of compositional and conformational ensembles, thereby providing key insights into macromolecular function. However, the inherent instability and dynamic nature of biological assemblies remain a tremendous challenge that often requires tailored approaches for successful implementation of the methodology. Here, we briefly describe the fundamentals of single-particle cryo-EM with an emphasis on covering the breadth of techniques and approaches, including low- and high-resolution methods, aiming to illustrate specific steps that are crucial for obtaining structural information by this method.

Electron Tomography: A Three-Dimensional Analytic Tool for Hard and Soft Materials Research

Three-dimensional (3D) structural analysis is essential to understand the relationship between the structure and function of an object. Many analytical techniques, such as X-ray diffraction, neutron spectroscopy, and electron microscopy imaging, are used to provide structural information. Transmission electron microscopy (TEM), one of the most popular analytic tools, has been widely used for structural analysis in both physical and biological sciences for many decades, in which 3D objects are projected into two-dimensional (2D) images. In many cases, 2D-projection images are insufficient to understand the relationship between the 3D structure and the function of nanoscale objects. Electron tomography (ET) is a technique that retrieves 3D structural information from a tilt series of 2D projections, and is gradually becoming a mature technology with sub-nanometer resolution. Distinct methods to overcome sample-based limitations have been separately developed in both physical and biological science, although they share some basic concepts of ET. This review discusses the common basis for 3D characterization, and specifies difficulties and solutions regarding both hard and soft materials research. It is hoped that novel solutions based on current state-of-the-art techniques for advanced applications in hybrid matter systems can be motivated.

The Structure of Immature Virus-Like Rous Sarcoma Virus Gag Particles Reveals a Structural Role for the p10 Domain in Assembly

The polyprotein Gag is the primary structural component of retroviruses. Gag consists of independently folded domains connected by flexible linkers. Interactions between the conserved capsid (CA) domains of Gag mediate formation of hexameric protein lattices that drive assembly of immature virus particles. Proteolytic cleavage of Gag by the viral protease (PR) is required for maturation of retroviruses from an immature form into an infectious form. Within the assembled Gag lattices of HIV-1 and Mason-Pfizer monkey virus (M-PMV), the C-terminal domain of CA adopts similar quaternary arrangements, while the N-terminal domain of CA is packed in very different manners. Here, we have used cryo-electron tomography and subtomogram averaging to study in vitro-assembled, immature virus-like Rous sarcoma virus (RSV) Gag particles and have determined the structure of CA and the surrounding regions to a resolution of ∼8 Å. We found that the C-terminal domain of RSV CA is arranged similarly to HIV-1 and M-PMV, whereas the N-terminal domain of CA adopts a novel arrangement in which the upstream p10 domain folds back into the CA lattice. In this position the cleavage site between CA and p10 appears to be inaccessible to PR. Below CA, an extended density is consistent with the presence of a six-helix bundle formed by the spacer-peptide region. We have also assessed the affect of lattice assembly on proteolytic processing by exogenous PR. The cleavage between p10 and CA is indeed inhibited in the assembled lattice, a finding consistent with structural regulation of proteolytic maturation.

IMPORTANCE

Retroviruses first assemble into immature virus particles, requiring interactions between Gag proteins that form a protein layer under the viral membrane. Subsequently, Gag is cleaved by the viral protease enzyme into separate domains, leading to rearrangement of the virus into its infectious form. It is important to understand how Gag is arranged within immature retroviruses, in order to understand how virus assembly occurs, and how maturation takes place. We used the techniques cryo-electron tomography and subtomogram averaging to obtain a detailed structural picture of the CA domains in immature assembled Rous sarcoma virus Gag particles. We found that part of Gag next to CA, called p10, folds back and interacts with CA when Gag assembles. This arrangement is different from that seen in HIV-1 and Mason-Pfizer monkey virus, illustrating further structural diversity of retroviral structures. The structure provides new information on how the virus assembles and undergoes maturation.

Integrative Modeling of Macromolecular Assemblies from Low to Near-Atomic Resolution

While conventional high-resolution techniques in structural biology are challenged by the size and flexibility of many biological assemblies, recent advances in low-resolution techniques such as cryo-electron microscopy (cryo-EM) and small angle X-ray scattering (SAXS) have opened up new avenues to define the structures of such assemblies. By systematically combining various sources of structural, biochemical and biophysical information, integrative modeling approaches aim to provide a unified structural description of such assemblies, starting from high-resolution structures of the individual components and integrating all available information from low-resolution experimental methods. In this review, we describe integrative modeling approaches, which use complementary data from either cryo-EM or SAXS. Specifically, we focus on the popular molecular dynamics flexible fitting (MDFF) method, which has been widely used for flexible fitting into cryo-EM maps. Second, we describe hybrid molecular dynamics, Rosetta Monte-Carlo and minimum ensemble search (MES) methods that can be used to incorporate SAXS into pseudoatomic structural models. We present concise descriptions of the two methods and their most popular alternatives, along with select illustrative applications to protein/nucleic acid assemblies involved in DNA replication and repair.

Sinorhizobium meliloti Phage ΦM9 Defines a New Group of T4 Superfamily Phages with Unusual Genomic Features but a Common T=16 Capsid

Relatively little is known about the phages that infect agriculturally important nitrogen-fixing rhizobial bacteria. Here we report the genome and cryo-electron microscopy structure of the Sinorhizobium meliloti-infecting T4 superfamily phage ΦM9. This phage and its close relative Rhizobium phage vB_RleM_P10VF define a new group of T4 superfamily phages. These phages are distinctly different from the recently characterized cyanophage-like S. meliloti phages of the ΦM12 group. Structurally, ΦM9 has a T=16 capsid formed from repeating units of an extended gp23-like subunit that assemble through interactions between one subunit and the adjacent E-loop insertion domain. Though genetically very distant from the cyanophages, the ΦM9 capsid closely resembles that of the T4 superfamily cyanophage Syn9. ΦM9 also has the same T=16 capsid architecture as the very distant phage SPO1 and the herpesviruses. Despite their overall lack of similarity at the genomic and structural levels, ΦM9 and S. meliloti phage ΦM12 have a small number of open reading frames in common that appear to encode structural proteins involved in interaction with the host and which may have been acquired by horizontal transfer. These proteins are predicted to encode tail baseplate proteins, tail fibers, tail fiber assembly proteins, and glycanases that cleave host exopolysaccharide.

IMPORTANCE Despite recent advances in the phylogenetic and structural characterization of bacteriophages, only a small number of phages of plant-symbiotic nitrogen-fixing soil bacteria have been studied at the molecular level. The effects of phage predation upon beneficial bacteria that promote plant growth remain poorly characterized. First steps in understanding these soil bacterium-phage dynamics are genetic, molecular, and structural characterizations of these groups of phages. The T4 superfamily phages are among the most complex phages; they have large genomes packaged within an icosahedral head and a long, contractile tail through which the DNA is delivered to host cells. This phylogenetic and structural study of S. meliloti-infecting T4 superfamily phage ΦM9 provides new insight into the diversity of this family. The comparison of structure-related genes in both ΦM9 and S. meliloti-infecting T4 superfamily phage ΦM12, which comes from a completely different lineage of these phages, allows the identification of host infection-related factors.

Mechanism of CRISPR-RNA guided recognition of DNA targets in Escherichia coli

In bacteria and archaea, short fragments of foreign DNA are integrated into Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) loci, providing a molecular memory of previous encounters with foreign genetic elements. In Escherichia coli, short CRISPR-derived RNAs are incorporated into a multi-subunit surveillance complex called Cascade (CRISPR-associated complex for antiviral defense). Recent structures of Cascade capture snapshots of this seahorse-shaped RNA-guided surveillance complex before and after binding to a DNA target. Here we determine a 3.2 Å x-ray crystal structure of Cascade in a new crystal form that provides insight into the mechanism of double-stranded DNA binding. Molecular dynamic simulations performed using available structures reveal functional roles for residues in the tail, backbone and belly subunits of Cascade that are critical for binding double-stranded DNA. Structural comparisons are used to make functional predictions and these predictions are tested in vivo and in vitro. Collectively, the results in this study reveal underlying mechanisms involved in target-induced conformational changes and highlight residues important in DNA binding and protospacer adjacent motif recognition.

Re-analysis of cryoEM data on HCV IRES bound to 40S subunit of human ribosome integrated with recent structural information suggests new contact regions between ribosomal proteins and HCV RNA

In this study, we combine available high resolution structural information on eukaryotic ribosomes with low resolution cryo-EM data on the Hepatitis C Viral RNA (IRES) human ribosome complex. Aided further by the prediction of RNA-protein interactions and restrained docking studies, we gain insights on their interaction at the residue level. We identified the components involved at the major and minor contact regions, and propose that there are energetically favorable local interactions between 40S ribosomal proteins and IRES domains. Domain II of the IRES interacts with ribosomal proteins S5 and S25 while the pseudoknot and the downstream domain IV region bind to ribosomal proteins S26, S28 and S5. We also provide support using UV cross-linking studies to validate our proposition of interaction between the S5 and IRES domains II and IV. We found that domain IIIe makes contact with the ribosomal protein S3a (S1e). Our model also suggests that the ribosomal protein S27 interacts with domain IIIc while S7 has a weak contact with a single base RNA bulge between junction IIIabc and IIId. The interacting residues are highly conserved among mammalian homologs while IRES RNA bases involved in contact do not show strict conservation. IRES RNA binding sites for S25 and S3a show the best conservation among related viral IRESs. The new contacts identified between ribosomal proteins and RNA are consistent with previous independent studies on RNA-binding properties of ribosomal proteins reported in literature, though information at the residue level is not available in previous studies.

Single-Particle Cryo-EM at Crystallographic Resolution

Until only a few years ago, single-particle electron cryo-microscopy (cryo-EM) was usually not the first choice for many structural biologists due to its limited resolution in the range of nanometer to subnanometer. Now, this method rivals X-ray crystallography in terms of resolution and can be used to determine atomic structures of macromolecules that are either refractory to crystallization or difficult to crystallize in specific functional states. In this review, I discuss the recent breakthroughs in both hardware and software that transformed cryo-microscopy, enabling understanding of complex biomolecules and their functions at atomic level.

Molecular dynamics simulations of large macromolecular complexes

Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

The importance of dynamics in integrative modeling of supramolecular assemblies

Revealing the atomistic architecture of supramolecular complexes is a fundamental step toward a deeper understanding of cellular functioning. To date, this formidable task is facilitated by an emerging array of integrative modeling approaches that combine experimental data from different sources. One major challenge these methods have to face is the treatment of the dynamic rearrangements of the individual subunits upon assembly. While this flexibility can be sampled at different levels, integrating native dynamic determinants with available experimental inputs can provide an effective way to reveal the molecular recognition mechanisms at the basis of supramolecular assembly.

Atomic-accuracy models from 4.5-Å cryo-electron microscopy data with density-guided iterative local refinement

We describe a general approach for refining protein structure models on the basis of cryo-electron microscopy maps with near-atomic resolution. The method integrates Monte Carlo sampling with local density-guided optimization, Rosetta all-atom refinement and real-space B-factor fitting. In tests on experimental maps of three different systems with 4.5-Å resolution or better, the method consistently produced models with atomic-level accuracy largely independently of starting-model quality, and it outperformed the molecular dynamics-based MDFF method. Cross-validated model quality statistics correlated with model accuracy over the three test systems.

Decrypting the structural, dynamic, and energetic basis of a monomeric kinesin interacting with a tubulin dimer in three ATPase states by all-atom molecular dynamics simulation

We have employed molecular dynamics (MD) simulation to investigate, with atomic details, the structural dynamics and energetics of three major ATPase states (ADP, APO, and ATP state) of a human kinesin-1 monomer in complex with a tubulin dimer. Starting from a recently solved crystal structure of ATP-like kinesin-tubulin complex by the Knossow lab, we have used flexible fitting of cryo-electron-microscopy maps to construct new structural models of the kinesin-tubulin complex in APO and ATP state, and then conducted extensive MD simulations (total 400 ns for each state), followed by flexibility analysis, principal component analysis, hydrogen bond analysis, and binding free energy analysis. Our modeling and simulation have revealed key nucleotide-dependent changes in the structure and flexibility of the nucleotide-binding pocket (featuring a highly flexible and open switch I in APO state) and the tubulin-binding site, and allosterically coupled motions driving the APO to ATP transition. In addition, our binding free energy analysis has identified a set of key residues involved in kinesin-tubulin binding. On the basis of our simulation, we have attempted to address several outstanding issues in kinesin study, including the possible roles of β-sheet twist and neck linker docking in regulating nucleotide release and binding, the structural mechanism of ADP release, and possible extension and shortening of α4 helix during the ATPase cycle. This study has provided a comprehensive structural and dynamic picture of kinesin's major ATPase states, and offered promising targets for future mutational and functional studies to investigate the molecular mechanism of kinesin motors.

Structural and biophysical methods to analyze clock function and mechanism

Structural approaches have provided insight into the mechanisms of circadian clock oscillators. This review focuses upon the myriad structural methods that have been applied to the molecular architecture of cyanobacterial circadian proteins, their interactions with each other, and the mechanism of the KaiABC posttranslational oscillator. X-ray crystallography and solution NMR were deployed to gain an understanding of the three-dimensional structures of the three proteins KaiA, KaiB, and KaiC that make up the inner timer in cyanobacteria. A hybrid structural biology approach including crystallography, electron microscopy, and solution scattering has shed light on the shapes of binary and ternary Kai protein complexes. Structural studies of the cyanobacterial oscillator demonstrate both the strengths and the limitations of the divide-and-conquer strategy. Thus, investigations of complexes involving domains and/or peptides have afforded valuable information into Kai protein interactions. However, high-resolution structural data are still needed at the level of complexes between the 360-kDa KaiC hexamer that forms the heart of the clock and its KaiA and KaiB partners.

Metastable intermediates as stepping stones on the maturation pathways of viral capsids

As they mature, many capsids undergo massive conformational changes that transform their stability, reactivity, and capacity for DNA. In some cases, maturation proceeds via one or more intermediate states. These structures represent local minima in a rich energy landscape that combines contributions from subunit folding, association of subunits into capsomers, and intercapsomer interactions. We have used scanning calorimetry and cryo-electron microscopy to explore the range of capsid conformations accessible to bacteriophage HK97. To separate conformational effects from those associated with covalent cross-linking (a stabilization mechanism of HK97), a cross-link-incompetent mutant was used. The mature capsid Head I undergoes an endothermic phase transition at 60°C in which it shrinks by 7%, primarily through changes in its hexamer conformation. The transition is reversible, with a half-life of ~3 min; however, >50% of reverted capsids are severely distorted or ruptured. This observation implies that such damage is a potential hazard of large-scale structural changes such as those involved in maturation. Assuming that the risk is lower for smaller changes, this suggests a rationalization for the existence of metastable intermediates: that they serve as stepping stones that preserve capsid integrity as it switches between the radically different conformations of its precursor and mature states.

Large-scale conformational changes are widespread in virus maturation and infection processes. These changes are accompanied by the release of conformational free energy as the virion (or fusogenic glycoprotein) switches from a precursor state to its mature state. Each state corresponds to a local minimum in an energy landscape. The conformational changes in capsid maturation are so radical that the question arises of how maturing capsids avoid being torn apart. Offering proof of principle, severe damage is inflicted when a bacteriophage HK97 capsid reverts from the (nonphysiological) state that it enters when heated past 60°C. We suggest that capsid proteins have been selected in part by the criterion of being able to avoid sustaining collateral damage as they mature. One way of achieving this—as with the HK97 capsid—involves breaking the overall transition down into several smaller steps in which the risk of damage is reduced.

Structure of the immature HIV-1 capsid in intact virus particles at 8.8 Å resolution

Human immunodeficiency virus type 1 (HIV-1) assembly proceeds in two stages. First, the 55 kilodalton viral Gag polyprotein assembles into a hexameric protein lattice at the plasma membrane of the infected cell, inducing budding and release of an immature particle. Second, Gag is cleaved by the viral protease, leading to internal rearrangement of the virus into the mature, infectious form. Immature and mature HIV-1 particles are heterogeneous in size and morphology, preventing high-resolution analysis of their protein arrangement in situ by conventional structural biology methods. Here we apply cryo-electron tomography and sub-tomogram averaging methods to resolve the structure of the capsid lattice within intact immature HIV-1 particles at subnanometre resolution, allowing unambiguous positioning of all α-helices. The resulting model reveals tertiary and quaternary structural interactions that mediate HIV-1 assembly. Strikingly, these interactions differ from those predicted by the current model based on in vitro-assembled arrays of Gag-derived proteins from Mason-Pfizer monkey virus. To validate this difference, we solve the structure of the capsid lattice within intact immature Mason-Pfizer monkey virus particles. Comparison with the immature HIV-1 structure reveals that retroviral capsid proteins, while having conserved tertiary structures, adopt different quaternary arrangements during virus assembly. The approach demonstrated here should be applicable to determine structures of other proteins at subnanometre resolution within heterogeneous environments.

Flexibility in the insulin receptor ectodomain enables docking of insulin in crystallographic conformation observed in a hormone-bound microreceptor

Insulin binding to the insulin receptor (IR) is the first key step in initiating downstream signaling cascades for glucose homeostasis in higher organisms. The molecular details of insulin recognition by IR are not yet completely understood, but a picture of hormone/receptor interactions at one of the epitopes (Site 1) is beginning to emerge from recent structural evidence. However, insulin-bound structures of truncated IR suggest that crystallographic conformation of insulin cannot be accommodated in the full IR ectodomain due to steric overlap of insulin with the first two type III fibronectin domains (F1 and F2), which are contributed to the insulin binding-pocket by the second subunit in the IR homodimer. A conformational change in the F1-F2 pair has thus been suggested. In this work, we present an all-atom structural model of complex of insulin and the IR ectodomain, where no structural overlap of insulin with the receptor domains (F1 and F2) is observed. This structural model was arrived at by flexibly fitting parts of our earlier insulin/IR all-atom model into the simulated density maps of crystallized constructs combined with conformational sampling from apo-IR solution conformations. Importantly, our experimentally-consistent model helps rationalize yet unresolved Site.

Molecular architecture and mechanism of the anaphase-promoting complex

The ubiquitination of cell cycle regulatory proteins by the anaphase-promoting complex/cyclosome (APC/C) controls sister chromatid segregation, cytokinesis and the establishment of the G1 phase of the cell cycle. The APC/C is an unusually large multimeric cullin-RING ligase. Its activity is strictly dependent on regulatory coactivator subunits that promote APC/C-substrate interactions and stimulate its catalytic reaction. Because the structures of many APC/C subunits and their organization within the assembly are unknown, the molecular basis for these processes is poorly understood. Here, from a cryo-electron microscopy reconstruction of a human APC/C-coactivator-substrate complex at 7.4 Å resolution, we have determined the complete secondary structural architecture of the complex. With this information we identified protein folds for structurally uncharacterized subunits, and the definitive location of all 20 APC/C subunits within the 1.2 MDa assembly. Comparison with apo APC/C shows that the coactivator promotes a profound allosteric transition involving displacement of the cullin-RING catalytic subunits relative to the degron-recognition module of coactivator and APC10. This transition is accompanied by increased flexibility of the cullin-RING subunits and enhanced affinity for UBCH10-ubiquitin, changes which may contribute to coactivator-mediated stimulation of APC/C E3 ligase activity.

Zooming in on disordered systems: neutron reflection studies of proteins associated with fluid membranes

Neutron reflectometry (NR) is an emerging experimental technique for the structural characterization of proteins interacting with fluid bilayer membranes under conditions that mimic closely the cellular environment. Thus, cellular processes can be emulated in artificial systems and their molecular basis studied by adding cellular components one at a time in a well-controlled environment while the resulting structures, or structural changes in response to external cues, are monitored with neutron reflection. In recent years, sample environments, data collection strategies and data analysis were continuously refined. The combination of these improvements increases the information which can be obtained from NR to an extent that enables structural characterization of protein-membrane complexes at a length scale that exceeds the resolution of the measurement by far. Ultimately, the combination of NR with molecular dynamics (MD) simulations can be used to cross-validate the results of the two techniques and provide atomic-scale structural models. This review discusses these developments in detail and demonstrates how they provide new windows into relevant biomedical problems. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.