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

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

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

Stabilization of Ribosomal RNA of the Small Subunit by Spermidine in Staphylococcus aureus

Cryo-electron microscopy is now used as a method of choice in structural biology for studying protein synthesis, a process mediated by the ribosome machinery. In order to achieve high-resolution structures using this approach, one needs to obtain homogeneous and stable samples, which requires optimization of ribosome purification in a species-dependent manner. This is especially critical for the bacterial small ribosomal subunit that tends to be unstable in the absence of ligands. Here, we report a protocol for purification of stable 30 S from the Gram-positive bacterium Staphylococcus aureus and its cryo-EM structures: in presence of spermidine at a resolution ranging between 3.4 and 3.6 Å and in its absence at 5.3 Å. Using biochemical characterization and cryo-EM, we demonstrate the importance of spermidine for stabilization of the 30 S via preserving favorable conformation of the helix 44.

Mechanism of Rad26-assisted rescue of stalled RNA polymerase II in transcription-coupled repair

Transcription-coupled repair is essential for the removal of DNA lesions from the transcribed genome. The pathway is initiated by CSB protein binding to stalled RNA polymerase II. Mutations impairing CSB function cause severe genetic disease. Yet, the ATP-dependent mechanism by which CSB powers RNA polymerase to bypass certain lesions while triggering excision of others is incompletely understood. Here we build structural models of RNA polymerase II bound to the yeast CSB ortholog Rad26 in nucleotide-free and bound states. This enables simulations and graph-theoretical analyses to define partitioning of this complex into dynamic communities and delineate how its structural elements function together to remodel DNA. We identify an allosteric pathway coupling motions of the Rad26 ATPase modules to changes in RNA polymerase and DNA to unveil a structural mechanism for CSB-assisted progression past less bulky lesions. Our models allow functional interpretation of the effects of Cockayne syndrome disease mutations.

The Viral Class II Membrane Fusion Machinery: Divergent Evolution from an Ancestral Heterodimer

A key step during the entry of enveloped viruses into cells is the merger of viral and cell lipid bilayers. This process is driven by a dedicated membrane fusion protein (MFP) present at the virion surface, which undergoes a membrane–fusogenic conformational change triggered by interactions with the target cell. Viral MFPs have been extensively studied structurally, and are divided into three classes depending on their three-dimensional fold. Because MFPs of the same class are found in otherwise unrelated viruses, their intra-class structural homology indicates horizontal gene exchange. We focus this review on the class II fusion machinery, which is composed of two glycoproteins that associate as heterodimers. They fold together in the ER of infected cells such that the MFP adopts a conformation primed to react to specific clues only upon contact with a target cell, avoiding premature fusion in the producer cell. We show that, despite having diverged in their 3D fold during evolution much more than the actual MFP, the class II accompanying proteins (AP) also derive from a distant common ancestor, displaying an invariant core formed by a β-ribbon and a C-terminal immunoglobulin-like domain playing different functional roles—heterotypic interactions with the MFP, and homotypic AP/AP contacts to form spikes, respectively. Our analysis shows that class II APs are easily identifiable with modern structural prediction algorithms, providing useful information in devising immunogens for vaccine design.

Structure of native HIV-1 cores and their interactions with IP6 and CypA

The viral capsid plays essential roles in HIV replication and is a major platform engaging host factors. To overcome challenges in study native capsid structure, we used the perfringolysin O to perforate the membrane of HIV-1 particles, thus allowing host proteins and small molecules to access the native capsid while improving cryo–electron microscopy image quality. Using cryo–electron tomography and subtomogram averaging, we determined the structures of native capsomers in the presence and absence of inositol hexakisphosphate (IP6) and cyclophilin A and constructed an all-atom model of a complete HIV-1 capsid. Our structures reveal two IP6 binding sites and modes of cyclophilin A interactions. Free energy calculations substantiate the two binding sites at R18 and K25 and further show a prohibitive energy barrier for IP6 to pass through the pentamer. Our results demonstrate that perfringolysin O perforation is a valuable tool for structural analyses of enveloped virus capsids and interactions with host cell factors.

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

[Figure: see text].

Double Proton Transfer during a Novel Tertiary α-Ketol Rearrangement in Ketol-Acid Reductoisomerase: A Water-Mediated, Metal-Catalyzed, Base-Induced Mechanism

(KARI) catalyzes the conversion of (S)-2-acetolactate or (S)-2-aceto-2-hydroxybutyrate to 2,3-dihydroxy-3-alkylbutyrate, the second step in the biosynthesis of branched chain amino acids (BCAAs). Because the BCAA biosynthetic pathway is present in bacteria, plants, and fungi, but absent in animals, it is an excellent target for the development of new-generation antibiotics and herbicides. Nevertheless, the mechanism of the KARI-catalyzed reaction has not yet been fully solved. In this study, we used iterative molecular dynamics (MD) flexible fitting-Rosetta techniques to optimize the three-dimensional solution structure of archaea KARI from Sulfolobus solfataricus (Sso-KARI) determined from cryo-electron microscopy. On the basis of the structure of the Sso-KARI/2Mg²⁺/NADH/(S)-2-acetolactate complex, we deciphered the catalytic mechanism of the KARI-mediated reaction through hybrid quantum mechanics/molecular mechanics MD simulations in conjunction with umbrella sampling. With an activation energy of only 6.06 kcal/mol, a water-mediated, metal-catalyzed, base-induced (WMMCBI) mechanism was preferred for deprotonation of the tertiary OH group of (S)-2-acetolactate in Sso-KARI. The WMMCBI mechanism for double proton transfer occurred within a proton wire route with two steps involving the formation of hydroxide: (i) Glu233 served as a general base to deprotonate the Mg²⁺-bound water, forming a hydroxide-coordinated Mg²⁺ ion; (ii) this hydroxide ion acted as a strong base that rapidly deprotonated the ternary OH group of the substrate. In contrast, the direct deprotonation of the substrate by Glu233 was kinetically unfavorable. This mechanism suggests a novel approach for designing catalysts for deprotonation and provides clues for the development of new-generation antibiotics and herbicides.

A Tool for Segmentation of Secondary Structures in 3D Cryo-EM Density Map Components Using Deep Convolutional Neural Networks

Although cryo-electron microscopy (cryo-EM) has been successfully used to derive atomic structures for many proteins, it is still challenging to derive atomic structures when the resolution of cryo-EM density maps is in the medium resolution range, such as 5-10 Å. Detection of protein secondary structures, such as helices and β-sheets, from cryo-EM density maps provides constraints for deriving atomic structures from such maps. As more deep learning methodologies are being developed for solving various molecular problems, effective tools are needed for users to access them. We have developed an effective software bundle, DeepSSETracer, for the detection of protein secondary structure from cryo-EM component maps in medium resolution. The bundle contains the network architecture and a U-Net model trained with a curriculum and gradient of episodic memory (GEM). The bundle integrates the deep neural network with the visualization capacity provided in ChimeraX. Using a Linux server that is remotely accessed by Windows users, it takes about 6 s on one CPU and one GPU for the trained deep neural network to detect secondary structures in a cryo-EM component map containing 446 amino acids. A test using 28 chain components of cryo-EM maps shows overall residue-level F1 scores of 0.72 and 0.65 to detect helices and β-sheets, respectively. Although deep learning applications are built on software frameworks, such as PyTorch and Tensorflow, our pioneer work here shows that integration of deep learning applications with ChimeraX is a promising and effective approach. Our experiments show that the F1 score measured at the residue level is an effective evaluation of secondary structure detection for individual classes. The test using 28 cryo-EM component maps shows that DeepSSETracer detects β-sheets more accurately than Emap2sec+, with a weighted average residue-level F1 score of 0.65 and 0.42, respectively. It also shows that Emap2sec+ detects helices more accurately than DeepSSETracer with a weighted average residue-level F1 score of 0.77 and 0.72 respectively.

Integration of Mass Spectrometry Data for Structural Biology

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

Poor Person's pH Simulation of Membrane Proteins

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

Prediction of Protein Complex Structure Using Surface-Induced Dissociation and Cryo-Electron Microscopy

A variety of techniques involving the use of mass spectrometry (MS) have been developed to obtain structural information on proteins and protein complexes. One example of these techniques, surface-induced dissociation (SID), has been used to study the oligomeric state and connectivity of protein complexes. Recently, we demonstrated that appearance energies (AE) could be extracted from SID experiments and that they correlate with structural features of specific protein-protein interfaces. While SID AE provides some structural information, the AE data alone are not sufficient to determine the structures of the complexes. For this reason, we sought to supplement the data with computational modeling, through protein-protein docking. In a previous study, we demonstrated that the scoring of structures generated from protein-protein docking could be improved with the inclusion of SID data; however, this work relied on knowledge of the correct tertiary structure and only built full complexes for a few cases. Here, we performed docking using input structures that require less prior knowledge, using homology models, unbound crystal structures, and bound+perturbed crystal structures. Using flexible ensemble docking (to build primarily subcomplexes from an ensemble of backbone structures), the RMSD₁₀₀ of all (15/15) predicted structures using the combined Rosetta, cryo-electron microscopy (cryo-EM), and SID score was less than 4 Å, compared to only 7/15 without SID and cryo-EM. Symmetric docking (which used symmetry to build full complexes) resulted in predicted structures with RMSD₁₀₀ less than 4 Å for 14/15 cases with experimental data, compared to only 5/15 without SID and cryo-EM. Finally, we also developed a confidence metric for which all (26/26) proteins flagged as high confidence were accurately predicted.

Energy Landscape of the SARS-CoV-2 Reveals Extensive Conformational Heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As for any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the spike protein conformations, and identify molecular rearrangements along the most-likely conformational path in the vicinity of the open (so called 1RBD-up) state. The resulting global and local atomic refinements reveal larger movements than those expected by comparing the reported 1RBD-up and 1RBD-down cryo-EM models. Here we report greater degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these conformations along the minimum free-energy pathway. A local analysis of seven key binding pockets reveals that six out them, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, and protein-glycan interaction sites, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This is reminiscent of a conformational pre-equilibrium. We found only one binding pocket, namely antibody AB-C135 to remain closed along the entire minimum free energy path, suggesting an induced fit mechanism for this enzyme.

Molecular determinants of Ebola nucleocapsid stability from molecular dynamics simulations

Ebola virus (EBOV) is a human pathogen with the ability to cause hemorrhagic fever and bleeding diathesis in hosts. The life cycle of EBOV depends on its nucleocapsid. The Ebola nucleocapsid consists of a helical assembly of nucleoproteins (NPs) encapsidating single-stranded viral RNA (ssRNA). Knowledge of the molecular determinants of Ebola nucleocapsid stability is essential for the development of therapeutics against EBOV. However, large degrees of freedom associated with the Ebola nucleocapsid helical assembly pose a computational challenge, thereby limiting the previous simulation studies to the level of monomers. In the present work, we have performed all atom molecular dynamics (MD) simulations of the helical assembly of EBOV nucleoproteins in the absence and presence of ssRNA. We found that ssRNA is essential for maintaining structural integrity of the nucleocapsid. Other molecular determinants observed to stabilize the nucleocapsid include NP-RNA and NP-NP interactions and ion distributions. Additionally, the structural and dynamical behavior of the nucleocapsid monomer depends on its position in the helical assembly. NP monomers present on the longitudinal edges of the helical tube are more exposed, flexible, and have weaker NP-NP interactions than those residing in the center. This work provides key structural features stabilizing the nucleocapsid that may serve as therapeutic targets.

Conformational distortion in a fibril-forming oligomer arrests alpha-Synuclein fibrillation and minimizes its toxic effects

The fibrillation pathway of alpha-Synuclein, the causative protein of Parkinson's disease, encompasses transient, heterogeneous oligomeric forms whose structural understanding and link to toxicity are not yet understood. We report that the addition of the physiologically-available small molecule heme at a sub-stoichiometric ratio to either monomeric or aggregated α-Syn, targets a His50 residue critical for fibril-formation and stabilizes the structurally-heterogeneous populations of aggregates into a minimally-toxic oligomeric state. Cryo-EM 3D reconstruction revealed a 'mace'-shaped structure of this monodisperse population of oligomers, which is comparable to a solid-state NMR Greek key-like motif (where the core residues are arranged in parallel in-register sheets with a Greek key topology at the C terminus) that forms the fundamental unit/kernel of protofilaments. Further structural analyses suggest that heme binding induces a distortion in the Greek key-like architecture of the mace oligomers, which impairs their further appending into protofilaments and fibrils. Additionally, our study reports a novel mechanism of prevention as well as reclamation of amyloid fibril formation by blocking an inter-protofilament His50 residue using a small molecule.

VESPER: global and local cryo-EM map alignment using local density vectors

An increasing number of density maps of biological macromolecules have been determined by cryo-electron microscopy (cryo-EM) and stored in the public database, EMDB. To interpret the structural information contained in EM density maps, alignment of maps is an essential step for structure modeling, comparison of maps, and for database search. Here, we developed VESPER, which captures the similarity of underlying molecular structures embedded in density maps by taking local gradient directions into consideration. Compared to existing methods, VESPER achieved substantially more accurate global and local alignment of maps as well as database retrieval.

Molecular Dynamics to Predict Cryo-EM: Capturing Transitions and Short-Lived Conformational States of Biomolecules

Single-particle cryogenic electron microscopy (cryo-EM) has revolutionized the field of the structural biology, providing an access to the atomic resolution structures of large biomolecular complexes in their near-native environment. Today's cryo-EM maps can frequently reach the atomic-level resolution, while often containing a range of resolutions, with conformationally variable regions obtained at 6 Å or worse. Low resolution density maps obtained for protein flexible domains, as well as the ensemble of coexisting conformational states arising from cryo-EM, poses new challenges and opportunities for Molecular Dynamics (MD) simulations. With the ability to describe the biomolecular dynamics at the atomic level, MD can extend the capabilities of cryo-EM, capturing the conformational variability and predicting biologically relevant short-lived conformational states. Here, we report about the state-of-the-art MD procedures that are currently used to refine, reconstruct and interpret cryo-EM maps. We show the capability of MD to predict short-lived conformational states, finding remarkable confirmation by cryo-EM structures subsequently solved. This has been the case of the CRISPR-Cas9 genome editing machinery, whose catalytically active structure has been predicted through both long-time scale MD and enhanced sampling techniques 2 years earlier than cryo-EM. In summary, this contribution remarks the ability of MD to complement cryo-EM, describing conformational landscapes and relating structural transitions to function, ultimately discerning relevant short-lived conformational states and providing mechanistic knowledge of biological function.

Structural insights into photosystem II assembly

Biogenesis of photosystem II (PSII), nature's water-splitting catalyst, is assisted by auxiliary proteins that form transient complexes with PSII components to facilitate stepwise assembly events. Using cryo-electron microscopy, we solved the structure of such a PSII assembly intermediate from Thermosynechococcus elongatus at 2.94 Å resolution. It contains three assembly factors (Psb27, Psb28 and Psb34) and provides detailed insights into their molecular function. Binding of Psb28 induces large conformational changes at the PSII acceptor side, which distort the binding pocket of the mobile quinone (Q_B) and replace the bicarbonate ligand of non-haem iron with glutamate, a structural motif found in reaction centres of non-oxygenic photosynthetic bacteria. These results reveal mechanisms that protect PSII from damage during biogenesis until water splitting is activated. Our structure further demonstrates how the PSII active site is prepared for the incorporation of the Mn₄CaO₅ cluster, which performs the unique water-splitting reaction.

The SARS-CoV-2 Spike variant D614G favors an open conformational state

The COVID-19 (coronavirus disease 2019) pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the "D-form" to the "G-form") that carried an amino acid substitution D614G in its "Spike" protein. The G-form is more infectious in vitro and is associated with increased viral loads in the upper airways. To gain insight into the molecular-level underpinnings of these characteristics, we used microsecond all-atom simulations. We show that changes in the protein energetics favor a higher population of infection-capable states in the G-form through release of asymmetry present in the D-form inter-protomer interactions. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive owing to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are critical for vaccine design.

Structural insights into the interplay of protein biogenesis factors with the 70S ribosome

Bacterial co-translational N-terminal methionine excision, an early event of nascent polypeptide chain processing, is mediated by two enzymes: peptide deformylase (PDF) and methionine aminopeptidase (MetAP). Trigger factor (TF), the only ribosome-associated bacterial chaperone, offers co-translational chaperoning assistance. Here, we present two high-resolution cryoelectron microscopy structures of tRNA-bound E. coli ribosome complexes showing simultaneous binding of PDF and TF, in the absence (3.4 Å) and presence of MetAP (4.1 Å). These structures establish molecular details of the interactions of the factors with the ribosome, and thereby reveal the structural basis of nascent chain processing. Our results suggest that simultaneous binding of all three factors is not a functionally favorable mechanism of nascent chain processing. Strikingly, an unusual structural distortion of the 70S ribosome, potentially driven by binding of multiple copies of MetAP, is observed when MetAP is incubated with a pre-formed PDF-TF-bound ribosome complex.

Molecular Dynamics Flexible Fitting: All You Want to Know About Resolution Exchange

In recent years, owing to the advances in instrumentation, cryo-EM has emerged as the go-to tool for obtaining high-resolution structures of biomolecular systems. However, building three-dimensional atomic structures of biomolecules from these high-resolution maps remains a concern for the traditional map-guided structure-determination schemes. Recently, we developed a computational tool, Resolution Exchange Molecular Dynamics Flexible Fitting (ReMDFF) to address this problem by re-refining a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution (Wang et al., J Struct Biol 204(2):319-328, 2018). In this chapter, we present a step-by-step outline for preparing, executing, and analyzing ReMDFF refinements of simple proteins and multimeric complexes. The structure determination of carbon monoxide dehydrogenase and Mg²⁺-channel CorA are employed as case studies. All scripts are provided via GitHub (Vant, Resolution exchange molecular dynamics flexible fitting (ReMDFF) all you want to know about flexible fitting, 2019, https://github.com/jvant/ReMDFF_Singharoy_Group.git ).

C-terminal troponin-I residues trap tropomyosin in the muscle thin filament blocked-state

Tropomyosin and troponin regulate muscle contraction by participating in a macromolecular scale steric-mechanism to control myosin-crossbridge - actin interactions and consequently contraction. At low-Ca²⁺, the C-terminal 30% of troponin subunit-I (TnI) is proposed to trap tropomyosin in a position on thin filaments that sterically interferes with myosin-binding, thus causing muscle relaxation. In contrast, at high-Ca²⁺, inhibition is released after the C-terminal domains dissociate from F-actin-tropomyosin as its component switch-peptide domain binds to the N-lobe of troponin-C (TnC). Recent, paradigm-shifting, cryo-EM reconstructions by the Namba group have revealed density attributed to TnI along cardiac muscle thin filaments at both low- and high-Ca²⁺ concentration. Modeling the reconstructions showed expected high-Ca²⁺ hydrophobic interactions of the TnI switch-peptide and TnC. However, under low-Ca²⁺ conditions, sparse interactions of TnI and tropomyosin, and in particular juxtaposition of non-polar switch-peptide residues and charged tropomyosin amino acids in the published model seem difficult to reconcile with an expected steric-blocking conformation. This anomaly is likely due to inaccurate fitting of tropomyosin into the cryo-EM volume. In the current study, the low-Ca²⁺ cryo-EM volume was fitted with a more accurate tropomyosin model and representation of cardiac TnI. Our results show that at low-Ca²⁺ a cluster of hydrophobic residues at the TnI switch-peptide and adjacent H₄ helix (Ala149, Ala151, Met 154, Leu159, Gly160, Ala161, Ala163, Leu167, Leu169, Ala171, Leu173) draw-in tropomyosin surface residues (Ile143, Ile146, Ala151, Ile154), presumably attracting the entire tropomyosin cable to its myosin-blocking position on actin. The modeling confirms that neighboring TnI "inhibitory domain" residues (Arg145, Arg148) bind to thin filaments at actin residue Asp25, as previously suggested. ClusPro docking of TnI residues 137-184 to actin-tropomyosin, including the TnI inhibitory-domain, switch-peptide and Helix H₄, verified the modeled configuration. Our residue-to-residue contact-mapping of the TnI-tropomyosin association lends itself to experimental validation and functional localization of disease-bearing mutations.

Refinement of a cryo-EM structure of hERG: Bridging structure and function

The human-ether-a-go-go-related gene (hERG) encodes the voltage-gated potassium channel (KCNH2 or Kv11.1, commonly known as hERG). This channel plays a pivotal role in the stability of phase 3 repolarization of the cardiac action potential. Although a high-resolution cryo-EM structure is available for its depolarized (open) state, the structure surprisingly did not feature many functionally important interactions established by previous biochemical and electrophysiology experiments. Using molecular dynamics flexible fitting (MDFF), we refined the structure and recovered the missing functionally relevant salt bridges in hERG in its depolarized state. We also performed electrophysiology experiments to confirm the functional relevance of a novel salt bridge predicted by our refinement protocol. Our work shows how refinement of a high-resolution cryo-EM structure helps to bridge the existing gap between the structure and function in the voltage-sensing domain (VSD) of hERG.

Negative-Stain Electron Microscopy Reveals Dramatic Structural Rearrangements in Ni-Fe-S-Dependent Carbon Monoxide Dehydrogenase/Acetyl-CoA Synthase

The Ni-Fe-S-containing A-cluster of acetyl-coenzyme A (CoA) synthase (ACS) assembles acetyl-CoA from carbon monoxide (CO), a methyl group (CH3+), and CoA. To accomplish this feat, ACS must bind CoA and interact with two other proteins that contribute the CO and CH3+, respectively: CO dehydrogenase (CODH) and corrinoid Fe-S protein (CFeSP). Previous structural data show that, in the model acetogen Moorella thermoacetica, domain 1 of ACS binds to CODH such that a 70-Å-long internal channel is created that allows CO to travel from CODH to the A-cluster. The A-cluster is largely buried and is inaccessible to CFeSP for methylation. Here we use electron microscopy to capture multiple snapshots of ACS that reveal previously uncharacterized domain motion, forming extended and hyperextended structural states. In these structural states, the A-cluster is accessible for methylation by CFeSP.

Hybrid methods for combined experimental and computational determination of protein structure

Knowledge of protein structure is paramount to the understanding of biological function, developing new therapeutics, and making detailed mechanistic hypotheses. Therefore, methods to accurately elucidate three-dimensional structures of proteins are in high demand. While there are a few experimental techniques that can routinely provide high-resolution structures, such as x-ray crystallography, nuclear magnetic resonance (NMR), and cryo-EM, which have been developed to determine the structures of proteins, these techniques each have shortcomings and thus cannot be used in all cases. However, additionally, a large number of experimental techniques that provide some structural information, but not enough to assign atomic positions with high certainty have been developed. These methods offer sparse experimental data, which can also be noisy and inaccurate in some instances. In cases where it is not possible to determine the structure of a protein experimentally, computational structure prediction methods can be used as an alternative. Although computational methods can be performed without any experimental data in a large number of studies, inclusion of sparse experimental data into these prediction methods has yielded significant improvement. In this Perspective, we cover many of the successes of integrative modeling, computational modeling with experimental data, specifically for protein folding, protein-protein docking, and molecular dynamics simulations. We describe methods that incorporate sparse data from cryo-EM, NMR, mass spectrometry, electron paramagnetic resonance, small-angle x-ray scattering, Förster resonance energy transfer, and genetic sequence covariation. Finally, we highlight some of the major challenges in the field as well as possible future directions.

Permeability of the HIV-1 capsid to metabolites modulates viral DNA synthesis

Reverse transcription, an essential event in the HIV-1 life cycle, requires deoxynucleotide triphosphates (dNTPs) to fuel DNA synthesis, thus requiring penetration of dNTPs into the viral capsid. The central cavity of the capsid protein (CA) hexamer reveals itself as a plausible channel that allows the passage of dNTPs into assembled capsids. Nevertheless, the molecular mechanism of nucleotide import into the capsid remains unknown. Employing all-atom molecular dynamics (MD) simulations, we established that cooperative binding between nucleotides inside a CA hexamer cavity results in energetically favorable conditions for passive translocation of dNTPs into the HIV-1 capsid. Furthermore, binding of the host cell metabolite inositol hexakisphosphate (IP₆) enhances dNTP import, while binding of synthesized molecules like benzenehexacarboxylic acid (BHC) inhibits it. The enhancing effect on reverse transcription by IP₆ and the consequences of interactions between CA and nucleotides were corroborated using atomic force microscopy, transmission electron microscopy, and virological assays. Collectively, our results provide an atomistic description of the permeability of the HIV-1 capsid to small molecules and reveal a novel mechanism for the involvement of metabolites in HIV-1 capsid stabilization, nucleotide import, and reverse transcription.

This study shows that the HIV-1 capsid protein, in addition to its structural role, regulates reverse transcription, an essential metabolic process of the virus, by mediating the import of nucleotides. In addition, host cell metabolites such as inositol phosphates are recruited by the capsid to regulate viral DNA synthesis.

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. Here, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

Molecular Architecture of the SARS-CoV-2 Virus

SARS-CoV-2 is an enveloped virus responsible for the COVID-19 pandemic. Despite recent advances in the structural elucidation of SARS-CoV-2 proteins, the detailed architecture of the intact virus remains to be unveiled. Here we report the molecular assembly of the authentic SARS-CoV-2 virus using cryoelectron tomography (cryo-ET) and subtomogram averaging (STA). Native structures of the S proteins in pre- and postfusion conformations were determined to average resolutions of 8.7-11 Å. Compositions of the N-linked glycans from the native spikes were analyzed by mass spectrometry, which revealed overall processing states of the native glycans highly similar to that of the recombinant glycoprotein glycans. The native conformation of the ribonucleoproteins (RNPs) and their higher-order assemblies were revealed. Overall, these characterizations revealed the architecture of the SARS-CoV-2 virus in exceptional detail and shed light on how the virus packs its ∼30-kb-long single-segmented RNA in the ∼80-nm-diameter lumen.

Polymerization and editing modes of a high-fidelity DNA polymerase are linked by a well-defined path

Proofreading by replicative DNA polymerases is a fundamental mechanism ensuring DNA replication fidelity. In proofreading, mis-incorporated nucleotides are excised through the 3'-5' exonuclease activity of the DNA polymerase holoenzyme. The exonuclease site is distal from the polymerization site, imposing stringent structural and kinetic requirements for efficient primer strand transfer. Yet, the molecular mechanism of this transfer is not known. Here we employ molecular simulations using recent cryo-EM structures and biochemical analyses to delineate an optimal free energy path connecting the polymerization and exonuclease states of E. coli replicative DNA polymerase Pol III. We identify structures for all intermediates, in which the transitioning primer strand is stabilized by conserved Pol III residues along the fingers, thumb and exonuclease domains. We demonstrate switching kinetics on a tens of milliseconds timescale and unveil a complete pol-to-exo switching mechanism, validated by targeted mutational experiments.

Assembly intermediates of orthoreovirus captured in the cell

Traditionally, molecular assembly pathways for viruses are inferred from high resolution structures of purified stable intermediates, low resolution images of cell sections and genetic approaches. Here, we directly visualise an unsuspected 'single shelled' intermediate for a mammalian orthoreovirus in cryo-preserved infected cells, by cryo-electron tomography of cellular lamellae. Particle classification and averaging yields structures to 5.6 Å resolution, sufficient to identify secondary structural elements and produce an atomic model of the intermediate, comprising 120 copies each of protein λ1 and σ2. This λ1 shell is 'collapsed' compared to the mature virions, with molecules pushed inwards at the icosahedral fivefolds by ~100 Å, reminiscent of the first assembly intermediate of certain prokaryotic dsRNA viruses. This supports the supposition that these viruses share a common ancestor, and suggests mechanisms for the assembly of viruses of the Reoviridae. Such methodology holds promise for dissecting the replication cycle of many viruses.

Structures of human pannexin 1 reveal ion pathways and mechanism of gating

Pannexin 1 (PANX1) is an ATP-permeable channel with critical roles in a variety of physiological functions such as blood pressure regulation1, apoptotic cell clearance2 and human oocyte development3. Here we present several structures of human PANX1 in a heptameric assembly at resolutions of up to 2.8 angström, including an apo state, a caspase-7-cleaved state and a carbenoxolone-bound state. We reveal a gating mechanism that involves two ion-conducting pathways. Under normal cellular conditions, the intracellular entry of the wide main pore is physically plugged by the C-terminal tail. Small anions are conducted through narrow tunnels in the intracellular domain. These tunnels connect to the main pore and are gated by a long linker between the N-terminal helix and the first transmembrane helix. During apoptosis, the C-terminal tail is cleaved by caspase, allowing the release of ATP through the main pore. We identified a carbenoxolone-binding site embraced by W74 in the extracellular entrance and a role for carbenoxolone as a channel blocker. We identified a gap-junction-like structure using a glycosylation-deficient mutant, N255A. Our studies provide a solid foundation for understanding the molecular mechanisms underlying the channel gating and inhibition of PANX1 and related large-pore channels.

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

Cryo-EM structure of the complete and ligand-saturated insulin receptor ectodomain

Glucose homeostasis and growth essentially depend on the hormone insulin engaging its receptor. Despite biochemical and structural advances, a fundamental contradiction has persisted in the current understanding of insulin ligand-receptor interactions. While biochemistry predicts two distinct insulin binding sites, 1 and 2, recent structural analyses have resolved only site 1. Using a combined approach of cryo-EM and atomistic molecular dynamics simulation, we present the structure of the entire dimeric insulin receptor ectodomain saturated with four insulin molecules. Complementing the previously described insulin-site 1 interaction, we present the first view of insulin bound to the discrete insulin receptor site 2. Insulin binding stabilizes the receptor ectodomain in a T-shaped conformation wherein the membrane-proximal domains converge and contact each other. These findings expand the current models of insulin binding to its receptor and of its regulation. In summary, we provide the structural basis for a comprehensive description of ligand-receptor interactions that ultimately will inform new approaches to structure-based drug design.

The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State

The COVID-19 pandemic underwent a rapid transition with the emergence of a SARS-CoV-2 variant that carried the amino acid substitution D614G in the Spike protein that became globally prevalent. The G-form is both more infectious in vitro and associated with increased viral loads in infected people. To gain insight into the mechanism underlying these distinctive characteristics, we employed multiple replicas of microsecond all-atom simulations to probe the molecular-level impact of this substitution on Spike closed and open states. The open state enables Spike interactions with its human cellular receptor, ACE2. Here we show that changes in the inter-protomer energetics due to the D614G substitution favor a higher population of infection-capable (open) states. The inter-protomer interactions between S1 and S2 subunits in the open state of the D-form are asymmetric. This asymmetry is resolved in the G-form due to the release of tensile hydrogen bonds resulting in an increased population of open conformations. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies.

Cryo-EM and Molecular Docking Shows Myosin Loop 4 Contacts Actin and Tropomyosin on Thin Filaments

The motor protein myosin drives muscle and nonmuscle motility by binding to and moving along actin of thin filaments. Myosin binding to actin also modulates interactions of the regulatory protein, tropomyosin, on thin filaments, and conversely tropomyosin affects myosin binding to actin. Insight into this reciprocity will facilitate a molecular level elucidation of tropomyosin regulation of myosin interaction with actin in muscle contraction, and in turn, promote better understanding of nonmuscle cell motility. Indeed, experimental approaches such as fiber diffraction, cryoelectron microscopy, and three-dimensional reconstruction have long been used to define regulatory interaction of tropomyosin and myosin on actin at a structural level. However, their limited resolution has not proven sufficient to determine tropomyosin and myosin contacts at an atomic-level and thus to fully substantiate possible functional contributions. To overcome this deficiency, we have followed a hybrid approach by performing new cryogenic electron microscopy reconstruction of myosin-S1-decorated F-actin-tropomyosin together with atomic scale protein-protein docking of tropomyosin to the EM models. Here, cryo-EM data were derived from filaments reconstituted with α1-actin, cardiac αα-tropomyosin, and masseter muscle β-myosin complexes; masseter myosin, which shares sequence identity with β-cardiac myosin-heavy chain, was used because of its stability in vitro. The data were used to build an atomic model of the tropomyosin cable that fits onto the actin filament between the tip of the myosin head and a cleft on the innermost edge of actin subunits. The docking and atomic scale fitting showed multiple discrete interactions of myosin loop 4 and acidic residues on successive 39-42 residue-long tropomyosin pseudorepeats. The contacts between S1 and tropomyosin on actin appear to compete with and displace ones normally found between actin and tropomyosin on myosin-free thin filaments in relaxed muscle, thus restructuring the filament during myosin-induced activation.

Structural mechanism of helicase loading onto replication origin DNA by ORC-Cdc6

The loading of the core Mcm2-7 helicase onto origin DNA is essential for the formation of replication forks and genomic stability. Here, we report two cryo-electron microscopy (cryo-EM) structures that capture helicase loader–helicase complexes just prior to DNA insertion. These pre-loading structures, combined with a computational simulation of the dynamic transition from the pre-loading state to the loaded state, provide crucial insights into the mechanism required for topologically linking the helicase to DNA. The helicase loading system is highly conserved from yeast to human, which means that the molecular principles described here for the yeast system are likely applicable to the human system.

DNA replication origins serve as sites of replicative helicase loading. In all eukaryotes, the six-subunit origin recognition complex (Orc1-6; ORC) recognizes the replication origin. During late M-phase of the cell-cycle, Cdc6 binds to ORC and the ORC–Cdc6 complex loads in a multistep reaction and, with the help of Cdt1, the core Mcm2-7 helicase onto DNA. A key intermediate is the ORC–Cdc6–Cdt1–Mcm2-7 (OCCM) complex in which DNA has been already inserted into the central channel of Mcm2-7. Until now, it has been unclear how the origin DNA is guided by ORC–Cdc6 and inserted into the Mcm2-7 hexamer. Here, we truncated the C-terminal winged-helix-domain (WHD) of Mcm6 to slow down the loading reaction, thereby capturing two loading intermediates prior to DNA insertion in budding yeast. In “semi-attached OCCM,” the Mcm3 and Mcm7 WHDs latch onto ORC–Cdc6 while the main body of the Mcm2-7 hexamer is not connected. In “pre-insertion OCCM,” the main body of Mcm2-7 docks onto ORC–Cdc6, and the origin DNA is bent and positioned adjacent to the open DNA entry gate, poised for insertion, at the Mcm2–Mcm5 interface. We used molecular simulations to reveal the dynamic transition from preloading conformers to the loaded conformers in which the loading of Mcm2-7 on DNA is complete and the DNA entry gate is fully closed. Our work provides multiple molecular insights into a key event of eukaryotic DNA replication.

Cryo-EM structure of islet amyloid polypeptide fibrils reveals similarities with amyloid-β fibrils

Amyloid deposits consisting of fibrillar islet amyloid polypeptide (IAPP) in pancreatic islets are associated with beta-cell loss and have been implicated in type 2 diabetes (T2D). Here, we applied cryo-EM to reconstruct densities of three dominant IAPP fibril polymorphs, formed in vitro from synthetic human IAPP. An atomic model of the main polymorph, built from a density map of 4.2-Å resolution, reveals two S-shaped, intertwined protofilaments. The segment 21-NNFGAIL-27, essential for IAPP amyloidogenicity, forms the protofilament interface together with Tyr37 and the amidated C terminus. The S-fold resembles polymorphs of Alzheimer's disease (AD)-associated amyloid-β (Aβ) fibrils, which might account for the epidemiological link between T2D and AD and reports on IAPP-Aβ cross-seeding in vivo. The results structurally link the early-onset T2D IAPP genetic polymorphism (encoding Ser20Gly) with the AD Arctic mutation (Glu22Gly) of Aβ and support the design of inhibitors and imaging probes for IAPP fibrils.

Cryo-EM structure of islet amyloid polypeptide fibrils reveals similarities with amyloid-β fibrils

Amyloid deposits consisting of fibrillar islet amyloid polypeptide (IAPP) in pancreatic islets are associated with beta-cell loss and have been implicated in type 2 diabetes (T2D). Here, we applied cryo-EM to reconstruct densities of three dominant IAPP fibril polymorphs, formed in vitro from synthetic human IAPP. An atomic model of the main polymorph, built from a density map of 4.2-Å resolution, reveals two S-shaped, intertwined protofilaments. The segment 21-NNFGAIL-27, essential for IAPP amyloidogenicity, forms the protofilament interface together with Tyr37 and the amidated C terminus. The S-fold resembles polymorphs of Alzheimer's disease (AD)-associated amyloid-β (Aβ) fibrils, which might account for the epidemiological link between T2D and AD and reports on IAPP-Aβ cross-seeding in vivo. The results structurally link the early-onset T2D IAPP genetic polymorphism (encoding Ser20Gly) with the AD Arctic mutation (Glu22Gly) of Aβ and support the design of inhibitors and imaging probes for IAPP fibrils.

Using NMR Chemical Shifts and Cryo-EM Density Restraints in Iterative Rosetta-MD Protein Structure Refinement

Cryo-EM has become one of the prime methods for protein structure elucidation, frequently yielding density maps with near-atomic or medium resolution. If protein structures cannot be deduced unambiguously from the density maps, computational structure refinement tools are needed to generate protein structural models. We have previously developed an iterative Rosetta-MDFF protocol that used cryo-EM densities to refine protein structures. Here we show that, in addition to cryo-EM densities, incorporation of other experimental restraints into the Rosetta-MDFF protocol further improved refined structures. We used NMR chemical shift (CS) data integrated with cryo-EM densities in our hybrid protocol in both the Rosetta step and the molecular dynamics (MD) simulations step. In 15 out of 18 cases for all MD rounds, the refinement results obtained when density maps and NMR chemical shift data were used in combination outperformed those of density map-only refinement. Notably, the improvement in refinement was highest when medium and low-resolution density maps were used. With our hybrid method, the RMSDs of final models obtained were always better than the RMSDs obtained by our previous protocol with just density refinement for both medium (6.9 Å) and low (9 Å) resolution maps. For all the six test proteins with medium resolution density maps (6.9 Å), the final refined structure RMSDs were lower for the hybrid method than for the cryo-EM only refinement. The final refined RMSDs were less than 1.5 Å when our hybrid protocol was used with 4 Å density maps. For four out of the six proteins the final RMSDs were even less than 1 Å. This study demonstrates that by using a combination of cryo-EM and NMR restraints, it is possible to refine structures to atomic resolution, outperforming single restraint refinement. This hybrid protocol will be a valuable tool when only low-resolution cryo-EM density data and NMR chemical shift data are available to refine structures.

Cylindrical Similarity Measurement for Helices in Medium-Resolution Cryo-Electron Microscopy Density Maps

Cryo-electron microscopy (cryo-EM) density maps at medium resolution (5-10 Å) reveal secondary structural features such as α-helices and β-sheets, but they lack the side chain details that would enable a direct structure determination. Among the more than 800 entries in the Electron Microscopy Data Bank (EMDB) of medium-resolution density maps that are associated with atomic models, a wide variety of similarities can be observed between maps and models. To validate such atomic models and to classify structural features, a local similarity criterion, the F1 score, is proposed and evaluated in this study. The F1 score is theoretically normalized to a range from zero to one, providing a local measure of cylindrical agreement between the density and atomic model of a helix. A systematic scan of 30,994 helices (among 3,247 protein chains modeled into medium-resolution density maps) reveals an actual range of observed F1 scores from 0.171 to 0.848, suggesting that the cylindrical fit of the current data is well stratified by the proposed measure. The best (highest) F1 scores tend to be associated with regions that exhibit high and spatially homogeneous local resolution (between 5 Å and 7.5 Å) in the helical density. The proposed F1 scores can be used as a discriminative classifier for validation studies and as a ranking criterion for cryo-EM density features in databases.

MemBlob database and server for identifying transmembrane regions using cryo-EM maps

SUMMARY

The identification of transmembrane helices in transmembrane proteins is crucial, not only to understand their mechanism of action but also to develop new therapies. While experimental data on the boundaries of membrane-embedded regions are sparse, this information is present in cryo-electron microscopy (cryo-EM) density maps and it has not been utilized yet for determining membrane regions. We developed a computational pipeline, where the inputs of a cryo-EM map, the corresponding atomistic structure, and the potential bilayer orientation determined by TMDET algorithm of a given protein result in an output defining the residues assigned to the bulk water phase, lipid interface and the lipid hydrophobic core. Based on this method, we built a database involving published cryo-EM protein structures and a server to be able to compute this data for newly obtained structures.

AVAILABILITY AND IMPLEMENTATION

http://memblob.hegelab.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Resveratrol as a nontoxic excipient stabilizes insulin in a bioactive hexameric form

Insulin aggregation is the leading cause of considerable reduction in the amount of active drug molecules in liquid formulations manufactured for diabetes management. Phenolic compounds, such as phenol and m-cresol, are routinely used to stabilize insulin in a hexameric form during its commercial preparation. However, long term usage of commercial insulin results in various adverse secondary responses, for which toxicity of the phenolic excipients is primarily responsible. In this study we aimed to find out a nontoxic insulin stabilizer. To that end, we have selected resveratrol, a natural polyphenol, as a prospective nontoxic insulin stabilizer because of its structural similarity with commercially used phenolic compounds. Atomic force microscopy visualization of resveratrol-treated human insulin revealed that resveratrol has a unique ability to arrest hINS in a soluble oligomeric form having discrete spherical morphology. Most importantly, resveratrol-treated insulin is nontoxic for HepG2 cells and it effectively maintains low blood glucose in a mouse model. Cryo-electron microscopy revealed 3D morphology of resveratrol-stabilized insulin that strikingly resembles crystal structures of insulin hexamer formulated with m-cresol. Significantly, we found that, in a condition inductive to amyloid fibrillation at physiological pH, resveratrol is capable of stabilizing insulin more efficiently than m-cresol. Thus, this study describes resveratrol as an effective nontoxic natural molecule that can be used for stabilizing insulin in a bioactive oligomeric form during its commercial formulation.

Mechanism of ribosome shutdown by RsfS in Staphylococcus aureus revealed by integrative structural biology approach

For the sake of energy preservation, bacteria, upon transition to stationary phase, tone down their protein synthesis. This process is favored by the reversible binding of small stress-induced proteins to the ribosome to prevent unnecessary translation. One example is the conserved bacterial ribosome silencing factor (RsfS) that binds to uL14 protein onto the large ribosomal subunit and prevents its association with the small subunit. Here we describe the binding mode of Staphylococcus aureus RsfS to the large ribosomal subunit and present a 3.2 Å resolution cryo-EM reconstruction of the 50S-RsfS complex together with the crystal structure of uL14-RsfS complex solved at 2.3 Å resolution. The understanding of the detailed landscape of RsfS-uL14 interactions within the ribosome shed light on the mechanism of ribosome shutdown in the human pathogen S. aureus and might deliver a novel target for pharmacological drug development and treatment of bacterial infections.

Measurement of atom resolvability in cryo-EM maps with Q-scores

Cryogenic electron microscopy (cryo-EM) maps are now at the point where resolvability of individual atoms can be achieved. However, resolvability is not necessarily uniform throughout the map. We introduce a quantitative parameter to characterize the resolvability of individual atoms in cryo-EM maps, the map Q-score. Q-scores can be calculated for atoms in proteins, nucleic acids, water, ligands and other solvent atoms, using models fitted to or derived from cryo-EM maps. Q-scores can also be averaged to represent larger features such as entire residues and nucleotides. Averaged over entire models, Q-scores correlate very well with the estimated resolution of cryo-EM maps for both protein and RNA. Assuming the models they are calculated from are well fitted to the map, Q-scores can be used as a measure of resolvability in cryo-EM maps at various scales, from entire macromolecules down to individual atoms. Q-score analysis of multiple cryo-EM maps of the same proteins derived from different laboratories confirms the reproducibility of structural features from side chains down to water and ion atoms.

Structural basis of transcription inhibition by the DNA mimic protein Ocr of bacteriophage T7

Bacteriophage T7 infects Escherichia coli and evades the host restriction/modification system. The Ocr protein of T7 was shown to exist as a dimer mimicking DNA and to bind to host restriction enzymes, thus preventing the degradation of the viral genome by the host. Here we report that Ocr can also inhibit host transcription by directly binding to bacterial RNA polymerase (RNAP) and competing with the recruitment of RNAP by sigma factors. Using cryo electron microscopy, we determined the structures of Ocr bound to RNAP. The structures show that an Ocr dimer binds to RNAP in the cleft, where key regions of sigma bind and where DNA resides during transcription synthesis, thus providing a structural basis for the transcription inhibition. Our results reveal the versatility of Ocr in interfering with host systems and suggest possible strategies that could be exploited in adopting DNA mimicry as a basis for forming novel antibiotics.

Using MemBlob to Analyze Transmembrane Regions Based on Cryo-EM Maps

Transmembrane proteins include membrane channels, pores, and receptors and, as such, comprise an important part of the proteome, yet our knowledge about them is much less complete than about soluble, globular proteins. An important aspect of transmembrane protein structure is their exact position within the lipid bilayer, a feature hard to investigate experimentally at the atomic level. Here we describe MemBlob, a novel approach utilizing difference electron density maps obtained by cryo-EM studies of transmembrane proteins. The idea behind is that the nonprotein part of such maps carries information on the exact localization of the membrane mimetics used in the experiment and can be used to extract the positional information of the protein within the membrane. MemBlob uses a structural model of the protein and an experimental electron density map to provide an estimation of the surface residues interacting with the membrane.

Structure and dynamics of the E. coli chemotaxis core signaling complex by cryo-electron tomography and molecular simulations

To enable the processing of chemical gradients, chemotactic bacteria possess large arrays of transmembrane chemoreceptors, the histidine kinase CheA, and the adaptor protein CheW, organized as coupled core-signaling units (CSU). Despite decades of study, important questions surrounding the molecular mechanisms of sensory signal transduction remain unresolved, owing especially to the lack of a high-resolution CSU structure. Here, we use cryo-electron tomography and sub-tomogram averaging to determine a structure of the Escherichia coli CSU at sub-nanometer resolution. Based on our experimental data, we use molecular simulations to construct an atomistic model of the CSU, enabling a detailed characterization of CheA conformational dynamics in its native structural context. We identify multiple, distinct conformations of the critical P4 domain as well as asymmetries in the localization of the P3 bundle, offering several novel insights into the CheA signaling mechanism.

Structure and mechanism of bactericidal mammalian perforin-2, an ancient agent of innate immunity

Perforin-2 (MPEG1) is thought to enable the killing of invading microbes engulfed by macrophages and other phagocytes, forming pores in their membranes. Loss of perforin-2 renders individual phagocytes and whole organisms significantly more susceptible to bacterial pathogens. Here, we reveal the mechanism of perforin-2 activation and activity using atomic structures of pre-pore and pore assemblies, high-speed atomic force microscopy, and functional assays. Perforin-2 forms a pre-pore assembly in which its pore-forming domain points in the opposite direction to its membrane-targeting domain. Acidification then triggers pore formation, via a 180° conformational change. This novel and unexpected mechanism prevents premature bactericidal attack and may have played a key role in the evolution of all perforin family proteins.

The structure of the endogenous ESX-3 secretion system

The ESX (or Type VII) secretion systems are protein export systems in mycobacteria and many Gram-positive bacteria that mediate a broad range of functions including virulence, conjugation, and metabolic regulation. These systems translocate folded dimers of WXG100-superfamily protein substrates across the cytoplasmic membrane. We report the cryo-electron microscopy structure of an ESX-3 system, purified using an epitope tag inserted with recombineering into the chromosome of the model organism Mycobacterium smegmatis. The structure reveals a stacked architecture that extends above and below the inner membrane of the bacterium. The ESX-3 protomer complex is assembled from a single copy of the EccB₃, EccC₃, and EccE₃ and two copies of the EccD₃ protein. In the structure, the protomers form a stable dimer that is consistent with assembly into a larger oligomer. The ESX-3 structure provides a framework for further study of these important bacterial transporters.