Structure is beauty, but not always truth.

Structural biology, as powerful as it is, can be misleading. We highlight four fundamental challenges: interpreting raw experimental data; accounting for motion; addressing the misleading nature of in vitro structures; and unraveling interactions between drugs and "anti-targets." Overcoming these challenges will amplify the impact of structural biology on drug discovery.

Bacteriophages inhibit and evade cGAS-like immune function in bacteria.

A fundamental strategy of eukaryotic antiviral immunity involves the cGAS enzyme, which synthesizes 2',3'-cGAMP and activates the effector STING. Diverse bacteria contain cGAS-like enzymes that produce cyclic oligonucleotides and induce anti-phage activity, known as CBASS. However, this activity has only been demonstrated through heterologous expression. Whether bacteria harboring CBASS antagonize and co-evolve with phages is unknown. Here, we identified an endogenous cGAS-like enzyme in Pseudomonas aeruginosa that generates 3',3'-cGAMP during phage infection, signals to a phospholipase effector, and limits phage replication. In response, phages express an anti-CBASS protein ("Acb2") that forms a hexamer with three 3',3'-cGAMP molecules and reduces phospholipase activity. Acb2 also binds to molecules produced by other bacterial cGAS-like enzymes (3',3'-cUU/UA/UG/AA) and mammalian cGAS (2',3'-cGAMP), suggesting broad inhibition of cGAS-based immunity. Upon Acb2 deletion, CBASS blocks lytic phage replication and lysogenic induction, but rare phages evade CBASS through major capsid gene mutations. Altogether, we demonstrate endogenous CBASS anti-phage function and strategies of CBASS inhibition and evasion.

Innate type 2 immunity controls hair follicle commensalism by Demodex mites.

Demodex mites are commensal parasites of hair follicles (HFs). Normally asymptomatic, inflammatory outgrowth of mites can accompany malnutrition, immune dysfunction, and aging, but mechanisms restricting Demodex outgrowth are not defined. Here, we show that control of mite HF colonization in mice required group 2 innate lymphoid cells (ILC2s), interleukin-13 (IL-13), and its receptor, IL-4Ra-IL-13Ra1. HF-associated ILC2s elaborated IL-13 that attenuated HFs and epithelial proliferation at anagen onset; in their absence, Demodex colonization led to increased epithelial proliferation and replacement of gene programs for repair by aberrant inflammation, leading to the loss of barrier function and HF exhaustion. Humans with rhinophymatous acne rosacea, an inflammatory condition associated with Demodex, had increased HF inflammation with decreased type 2 cytokines, consistent with the inverse relationship seen in mice. Our studies uncover a key role for skin ILC2s and IL-13, which comprise an immune checkpoint that sustains cutaneous integrity and restricts pathologic infestation by colonizing HF mites.

From structure to systems: high-resolution, quantitative genetic analysis of RNA polymerase II.

RNA polymerase II (RNAPII) lies at the core of dynamic control of gene expression. Using 53 RNAPII point mutants, we generated a point mutant epistatic miniarray profile (pE-MAP) comprising ∼60,000 quantitative genetic interactions in Saccharomyces cerevisiae. This analysis enabled functional assignment of RNAPII subdomains and uncovered connections between individual regions and other protein complexes. Using splicing microarrays and mutants that alter elongation rates in vitro, we found an inverse relationship between RNAPII speed and in vivo splicing efficiency. Furthermore, the pE-MAP classified fast and slow mutants that favor upstream and downstream start site selection, respectively. The striking coordination of polymerization rate with transcription initiation and splicing suggests that transcription rate is tuned to regulate multiple gene expression steps. The pE-MAP approach provides a powerful strategy to understand other multifunctional machines at amino acid resolution.

Recommendations for accelerating open preprint peer review to improve the culture of science.

Peer review is an important part of the scientific process, but traditional peer review at journals is coming under increased scrutiny for its inefficiency and lack of transparency. As preprints become more widely used and accepted, they raise the possibility of rethinking the peer-review process. Preprints are enabling new forms of peer review that have the potential to be more thorough, inclusive, and collegial than traditional journal peer review, and to thus fundamentally shift the culture of peer review toward constructive collaboration. In this Consensus View, we make a call to action to stakeholders in the community to accelerate the growing momentum of preprint sharing and provide recommendations to empower researchers to provide open and constructive peer review for preprints.

Systematic functional prioritization of protein posttranslational modifications.

Protein function is often regulated by posttranslational modifications (PTMs), and recent advances in mass spectrometry have resulted in an exponential increase in PTM identification. However, the functional significance of the vast majority of these modifications remains unknown. To address this problem, we compiled nearly 200,000 phosphorylation, acetylation, and ubiquitination sites from 11 eukaryotic species, including 2,500 newly identified ubiquitylation sites for Saccharomyces cerevisiae. We developed methods to prioritize the functional relevance of these PTMs by predicting those that likely participate in cross-regulatory events, regulate domain activity, or mediate protein-protein interactions. PTM conservation within domain families identifies regulatory "hot spots" that overlap with functionally important regions, a concept that we experimentally validated on the HSP70 domain family. Finally, our analysis of the evolution of PTM regulation highlights potential routes for neutral drift in regulatory interactions and suggests that only a fraction of modification sites are likely to have a significant biological role.

Synthetic group A streptogramin antibiotics that overcome Vat resistance.

Natural products serve as chemical blueprints for most antibiotics in clinical use. The evolutionary process by which these molecules arise is inherently accompanied by the co-evolution of resistance mechanisms that shorten the clinical lifetime of any given class of antibiotics1. Virginiamycin acetyltransferase (Vat) enzymes are resistance proteins that provide protection against streptogramins2, potent antibiotics against Gram-positive bacteria that inhibit the bacterial ribosome3. Owing to the challenge of selectively modifying the chemically complex, 23-membered macrocyclic scaffold of group A streptogramins, analogues that overcome the resistance conferred by Vat enzymes have not been previously developed2. Here we report the design, synthesis, and antibacterial evaluation of group A streptogramin antibiotics with extensive structural variability. Using cryo-electron microscopy and forcefield-based refinement, we characterize the binding of eight analogues to the bacterial ribosome at high resolution, revealing binding interactions that extend into the peptidyl tRNA-binding site and towards synergistic binders that occupy the nascent peptide exit tunnel. One of these analogues has excellent activity against several streptogramin-resistant strains of Staphylococcus aureus, exhibits decreased rates of acetylation in vitro, and is effective at lowering bacterial load in a mouse model of infection. Our results demonstrate that the combination of rational design and modular chemical synthesis can revitalize classes of antibiotics that are limited by naturally arising resistance mechanisms.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2.

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Molecular evidence of widespread benzimidazole drug resistance in Ancylostoma caninum from domestic dogs throughout the USA and discovery of a novel ß-tubulin benzimidazole resistance mutation.

Ancylostoma caninum is an important zoonotic gastrointestinal nematode of dogs worldwide and a close relative of human hookworms. We recently reported that racing greyhound dogs in the USA are infected with A. caninum that are commonly resistant to multiple anthelmintics. Benzimidazole resistance in A. caninum in greyhounds was associated with a high frequency of the canonical F167Y(TTC>TAC) isotype-1 ß-tubulin mutation. In this work, we show that benzimidazole resistance is remarkably widespread in A. caninum from domestic dogs across the USA. First, we identified and showed the functional significance of a novel benzimidazole isotype-1 ß-tubulin resistance mutation, Q134H(CAA>CAT). Several benzimidazole resistant A. caninum isolates from greyhounds with a low frequency of the F167Y(TTC>TAC) mutation had a high frequency of a Q134H(CAA>CAT) mutation not previously reported from any eukaryotic pathogen in the field. Structural modeling predicted that the Q134 residue is directly involved in benzimidazole drug binding and that the 134H substitution would significantly reduce binding affinity. Introduction of the Q134H substitution into the C. elegans ß-tubulin gene ben-1, by CRISPR-Cas9 editing, conferred similar levels of resistance as a ben-1 null allele. Deep amplicon sequencing on A. caninum eggs from 685 hookworm positive pet dog fecal samples revealed that both mutations were widespread across the USA, with prevalences of 49.7% (overall mean frequency 54.0%) and 31.1% (overall mean frequency 16.4%) for F167Y(TTC>TAC) and Q134H(CAA>CAT), respectively. Canonical codon 198 and 200 benzimidazole resistance mutations were absent. The F167Y(TTC>TAC) mutation had a significantly higher prevalence and frequency in Western USA than in other regions, which we hypothesize is due to differences in refugia. This work has important implications for companion animal parasite control and the potential emergence of drug resistance in human hookworms.

A single inactivating amino acid change in the SARS-CoV-2 NSP3 Mac1 domain attenuates viral replication in vivo.

Despite unprecedented efforts, our therapeutic arsenal against SARS-CoV-2 remains limited. The conserved macrodomain 1 (Mac1) in NSP3 is an enzyme exhibiting ADP-ribosylhydrolase activity and a possible drug target. To determine the role of Mac1 catalytic activity in viral replication, we generated recombinant viruses and replicons encoding a catalytically inactive NSP3 Mac1 domain by mutating a critical asparagine in the active site. While substitution to alanine (N40A) reduced catalytic activity by ~10-fold, mutations to aspartic acid (N40D) reduced activity by ~100-fold relative to wild-type. Importantly, the N40A mutation rendered Mac1 unstable in vitro and lowered expression levels in bacterial and mammalian cells. When incorporated into SARS-CoV-2 molecular clones, the N40D mutant only modestly affected viral fitness in immortalized cell lines, but reduced viral replication in human airway organoids by 10-fold. In mice, the N40D mutant replicated at >1000-fold lower levels compared to the wild-type virus while inducing a robust interferon response; all animals infected with the mutant virus survived infection. Our data validate the critical role of SARS-CoV-2 NSP3 Mac1 catalytic activity in viral replication and as a promising therapeutic target to develop antivirals.

Accurate positioning of functional residues with robotics-inspired computational protein design.

SignificanceComputational protein design promises to advance applications in medicine and biotechnology by creating proteins with many new and useful functions. However, new functions require the design of specific and often irregular atom-level geometries, which remains a major challenge. Here, we develop computational methods that design and predict local protein geometries with greater accuracy than existing methods. Then, as a proof of concept, we leverage these methods to design new protein conformations in the enzyme ketosteroid isomerase that change the protein's preference for a key functional residue. Our computational methods are openly accessible and can be applied to the design of other intricate geometries customized for new user-defined protein functions.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

A counter-enzyme complex regulates glutamate metabolism in Bacillus subtilis.

Multi-enzyme assemblies composed of metabolic enzymes catalyzing sequential reactions are being increasingly studied. Here, we report the discovery of a 1.6 megadalton multi-enzyme complex from Bacillus subtilis composed of two enzymes catalyzing opposite ('counter-enzymes') rather than sequential reactions: glutamate synthase (GltAB) and glutamate dehydrogenase (GudB), which make and break glutamate, respectively. In vivo and in vitro studies show that the primary role of complex formation is to inhibit the activity of GudB. Using cryo-electron microscopy, we elucidated the structure of the complex and the molecular basis of inhibition of GudB by GltAB. The complex exhibits unusual oscillatory progress curves and is necessary for both planktonic growth, in glutamate-limiting conditions, and for biofilm growth, in glutamate-rich media. The regulation of a key metabolic enzyme by complexing with its counter enzyme may thus enable cell growth under fluctuating glutamate concentrations.

Chemoenzymatic Syntheses of Fluorine-18-Labeled Disaccharides from [18F] FDG Yield Potent Sensors of Living Bacteria In Vivo.

Chemoenzymatic techniques have been applied extensively to pharmaceutical development, most effectively when routine synthetic methods fail. The regioselective and stereoselective construction of structurally complex glycans is an elegant application of this approach that is seldom applied to positron emission tomography (PET) tracers. We sought a method to dimerize 2-deoxy-[18F]-fluoro-d-glucose ([18F]FDG), the most common tracer used in clinical imaging, to form [18F]-labeled disaccharides for detecting microorganisms in vivo based on their bacteria-specific glycan incorporation. When [18F]FDG was reacted with ß-d-glucose-1-phosphate in the presence of maltose phosphorylase, the α-1,4- and α-1,3-linked products 2-deoxy-[18F]-fluoro-maltose ([18F]FDM) and 2-deoxy-2-[18F]-fluoro-sakebiose ([18F]FSK) were obtained. This method was further extended with the use of trehalose (α,α-1,1), laminaribiose (ß-1,3), and cellobiose (ß-1,4) phosphorylases to synthesize 2-deoxy-2-[18F]fluoro-trehalose ([18F]FDT), 2-deoxy-2-[18F]fluoro-laminaribiose ([18F]FDL), and 2-deoxy-2-[18F]fluoro-cellobiose ([18F]FDC). We subsequently tested [18F]FDM and [18F]FSK in vitro, showing accumulation by several clinically relevant pathogens including Staphylococcus aureus and Acinetobacter baumannii, and demonstrated their specific uptake in vivo. Both [18F]FDM and [18F]FSK were stable in human serum with high accumulation in preclinical infection models. The synthetic ease and high sensitivity of [18F]FDM and [18F]FSK to S. aureus including methicillin-resistant (MRSA) strains strongly justify clinical translation of these tracers to infected patients. Furthermore, this work suggests that chemoenzymatic radiosyntheses of complex [18F]FDG-derived oligomers will afford a wide array of PET radiotracers for infectious and oncologic applications.

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray-derived conformational ensembles.

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme's reaction cycle. To examine positioning and motions and test catalytic proposals, we collected "room temperature" X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI's reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI's general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure-function relationships, providing the basis for a new era of "ensemble-function" interrogation of enzymes.

Integration of software tools for integrative modeling of biomolecular systems.

Integrative modeling computes a model based on varied types of input information, be it from experiments or prior models. Often, a type of input information will be best handled by a specific modeling software package. In such a case, we desire to integrate our integrative modeling software package, Integrative Modeling Platform (IMP), with software specialized to the computational demands of the modeling problem at hand. After several attempts, however, we have concluded that even in collaboration with the software's developers, integration is either impractical or impossible. The reasons for the intractability of integration include software incompatibilities, differing modeling logic, the costs of collaboration, and academic incentives. In the integrative modeling software ecosystem, several large modeling packages exist with often redundant tools. We reason, therefore, that the other development groups have similarly concluded that the benefit of integration does not justify the cost. As a result, modelers are often restricted to the set of tools within a single software package. The inability to integrate tools from distinct software negatively impacts the quality of the models and the efficiency of the modeling. As the complexity of modeling problems grows, we seek to galvanize developers and modelers to consider the long-term benefit that software interoperability yields. In this article, we formulate a demonstrative set of software standards for implementing a model search using tools from independent software packages and discuss our efforts to integrate IMP and the crystallography suite Phenix within the Bayesian modeling framework.

Assessment of the nucleotide modifications in the high-resolution cryo-electron microscopy structure of the Escherichia coli 50S subunit.

Post-transcriptional ribosomal RNA (rRNA) modifications are present in all organisms, but their exact functional roles and positions are yet to be fully characterized. Modified nucleotides have been implicated in the stabilization of RNA structure and regulation of ribosome biogenesis and protein synthesis. In some instances, rRNA modifications can confer antibiotic resistance. High-resolution ribosome structures are thus necessary for precise determination of modified nucleotides' positions, a task that has previously been accomplished by X-ray crystallography. Here, we present a cryo-electron microscopy (cryo-EM) structure of the Escherichia coli 50S subunit at an average resolution of 2.2 Å as an additional approach for mapping modification sites. Our structure confirms known modifications present in 23S rRNA and additionally allows for localization of Mg2+ ions and their coordinated water molecules. Using our cryo-EM structure as a testbed, we developed a program for assessment of cryo-EM map quality. This program can be easily used on any RNA-containing cryo-EM structure, and an associated Coot plugin allows for visualization of validated modifications, making it highly accessible.

Effects of α-tubulin acetylation on microtubule structure and stability.

Acetylation of K40 in α-tubulin is the sole posttranslational modification to mark the luminal surface of microtubules. It is still controversial whether its relationship with microtubule stabilization is correlative or causative. We have obtained high-resolution cryo-electron microscopy (cryo-EM) reconstructions of pure samples of αTAT1-acetylated and SIRT2-deacetylated microtubules to visualize the structural consequences of this modification and reveal its potential for influencing the larger assembly properties of microtubules. We modeled the conformational ensembles of the unmodified and acetylated states by using the experimental cryo-EM density as a structural restraint in molecular dynamics simulations. We found that acetylation alters the conformational landscape of the flexible loop that contains αK40. Modification of αK40 reduces the disorder of the loop and restricts the states that it samples. We propose that the change in conformational sampling that we describe, at a location very close to the lateral contacts site, is likely to affect microtubule stability and function.

Mix-and-inject XFEL crystallography reveals gated conformational dynamics during enzyme catalysis.

How changes in enzyme structure and dynamics facilitate passage along the reaction coordinate is a fundamental unanswered question. Here, we use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL), ambient-temperature X-ray crystallography, computer simulations, and enzyme kinetics to characterize how covalent catalysis modulates isocyanide hydratase (ICH) conformational dynamics throughout its catalytic cycle. We visualize this previously hypothetical reaction mechanism, directly observing formation of a thioimidate covalent intermediate in ICH microcrystals during catalysis. ICH exhibits a concerted helical displacement upon active-site cysteine modification that is gated by changes in hydrogen bond strength between the cysteine thiolate and the backbone amide of the highly strained Ile152 residue. These catalysis-activated motions permit water entry into the ICH active site for intermediate hydrolysis. Mutations at a Gly residue (Gly150) that modulate helical mobility reduce ICH catalytic turnover and alter its pre-steady-state kinetic behavior, establishing that helical mobility is important for ICH catalytic efficiency. These results demonstrate that MISC can capture otherwise elusive aspects of enzyme mechanism and dynamics in microcrystalline samples, resolving long-standing questions about the connection between nonequilibrium protein motions and enzyme catalysis.

State of the structure address on MET receptor activation by HGF.

The MET receptor tyrosine kinase (RTK) and its cognate ligand hepatocyte growth factor (HGF) comprise a signaling axis essential for development, wound healing and tissue homeostasis. Aberrant HGF/MET signaling is a driver of many cancers and contributes to drug resistance to several approved therapeutics targeting other RTKs, making MET itself an important drug target. In RTKs, homeostatic receptor signaling is dependent on autoinhibition in the absence of ligand binding and orchestrated set of conformational changes induced by ligand-mediated receptor dimerization that result in activation of the intracellular kinase domains. A fundamental understanding of these mechanisms in the MET receptor remains incomplete, despite decades of research. This is due in part to the complex structure of the HGF ligand, which remains unknown in its full-length form, and a lack of high-resolution structures of the complete MET extracellular portion in an apo or ligand-bound state. A current view of HGF-dependent MET activation has evolved from biochemical and structural studies of HGF and MET fragments and here we review what these findings have thus far revealed.