Publisher Correction: Structural basis of ligand recognition at the human MT1 melatonin receptor.

Change history: In this Letter, the rotation signs around 90°, 135° and 15° were missing and in the HTML, Extended Data Tables 2 and 3 were the wrong tables; these errors have been corrected online.

Structural basis of ligand recognition at the human MT1 melatonin receptor.

Melatonin (N-acetyl-5-methoxytryptamine) is a neurohormone that maintains circadian rhythms1 by synchronization to environmental cues and is involved in diverse physiological processes2 such as the regulation of blood pressure and core body temperature, oncogenesis, and immune function3. Melatonin is formed in the pineal gland in a light-regulated manner4 by enzymatic conversion from 5-hydroxytryptamine (5-HT or serotonin), and modulates sleep and wakefulness5 by activating two high-affinity G-protein-coupled receptors, type 1A (MT1) and type 1B (MT2)3,6. Shift work, travel, and ubiquitous artificial lighting can disrupt natural circadian rhythms; as a result, sleep disorders affect a substantial population in modern society and pose a considerable economic burden7. Over-the-counter melatonin is widely used to alleviate jet lag and as a safer alternative to benzodiazepines and other sleeping aids8,9, and is one of the most popular supplements in the United States10. Here, we present high-resolution room-temperature X-ray free electron laser (XFEL) structures of MT1 in complex with four agonists: the insomnia drug ramelteon11, two melatonin analogues, and the mixed melatonin-serotonin antidepressant agomelatine12,13. The structure of MT2 is described in an accompanying paper14. Although the MT1 and 5-HT receptors have similar endogenous ligands, and agomelatine acts on both receptors, the receptors differ markedly in the structure and composition of their ligand pockets; in MT1, access to the ligand pocket is tightly sealed from solvent by extracellular loop 2, leaving only a narrow channel between transmembrane helices IV and V that connects it to the lipid bilayer. The binding site is extremely compact, and ligands interact with MT1 mainly by strong aromatic stacking with Phe179 and auxiliary hydrogen bonds with Asn162 and Gln181. Our structures provide an unexpected example of atypical ligand entry for a non-lipid receptor, lay the molecular foundation of ligand recognition by melatonin receptors, and will facilitate the design of future tool compounds and therapeutic agents, while their comparison to 5-HT receptors yields insights into the evolution and polypharmacology of G-protein-coupled receptors.

Fitting high-resolution electron density maps from atomic models to solution scattering data.

Solution scattering techniques, such as small- and wide-angle X-ray scattering (SWAXS), provide valuable insights into the structure and dynamics of biological macromolecules in solution. In this study, we present an approach to accurately predict solution X-ray scattering profiles at wide angles from atomic models by generating high-resolution electron density maps. Our method accounts for the excluded volume of bulk solvent by calculating unique adjusted atomic volumes directly from the atomic coordinates. This approach eliminates the need for one of the free fitting parameters commonly used in existing algorithms, resulting in improved accuracy of the calculated SWAXS profile. An implicit model of the hydration shell is generated that uses the form factor of water. Two parameters, namely the bulk solvent density and the mean hydration shell contrast, are adjusted to best fit the data. Results using eight publicly available SWAXS profiles show high-quality fits to the data. In each case, the optimized parameter values show small adjustments demonstrating that the default values are close to the true solution. Disabling parameter optimization produces significantly more accurate predicted scattering profiles compared to the leading software. The algorithm is computationally efficient, comparable to the leading software and up to 10 times faster for large molecules. The algorithm is encoded in a command line script called denss.pdb2mrc.py and is available open source as part of the DENSS v1.7.0 software package. In addition to improving the ability to compare atomic models to experimental SWAXS data, these developments pave the way for increasing the accuracy of modeling algorithms using SWAXS data and decreasing the risk of overfitting.

Structure of the full-length glucagon class B G-protein-coupled receptor.

The human glucagon receptor, GCGR, belongs to the class B G-protein-coupled receptor family and plays a key role in glucose homeostasis and the pathophysiology of type 2 diabetes. Here we report the 3.0 Å crystal structure of full-length GCGR containing both the extracellular domain and transmembrane domain in an inactive conformation. The two domains are connected by a 12-residue segment termed the stalk, which adopts a ß-strand conformation, instead of forming an α-helix as observed in the previously solved structure of the GCGR transmembrane domain. The first extracellular loop exhibits a ß-hairpin conformation and interacts with the stalk to form a compact ß-sheet structure. Hydrogen-deuterium exchange, disulfide crosslinking and molecular dynamics studies suggest that the stalk and the first extracellular loop have critical roles in modulating peptide ligand binding and receptor activation. These insights into the full-length GCGR structure deepen our understanding of the signalling mechanisms of class B G-protein-coupled receptors.

Snapshot of an oxygen intermediate in the catalytic reaction of cytochrome c oxidase.

Cytochrome c oxidase (CcO) reduces dioxygen to water and harnesses the chemical energy to drive proton translocation across the inner mitochondrial membrane by an unresolved mechanism. By using time-resolved serial femtosecond crystallography, we identified a key oxygen intermediate of bovine CcO. It is assigned to the PR-intermediate, which is characterized by specific redox states of the metal centers and a distinct protein conformation. The heme a3 iron atom is in a ferryl (Fe4+ = O2-) configuration, and heme a and CuB are oxidized while CuA is reduced. A Helix-X segment is poised in an open conformational state; the heme a farnesyl sidechain is H-bonded to S382, and loop-I-II adopts a distinct structure. These data offer insights into the mechanism by which the oxygen chemistry is coupled to unidirectional proton translocation.

Ab initio electron density determination directly from solution scattering data.

Using a novel iterative structure factor retrieval algorithm, here I show that electron density can be directly calculated from solution scattering data without modeling. The algorithm was validated with experimental data from 12 different biological macromolecules. This approach avoids many of the assumptions limiting the resolution and accuracy of modeling algorithms by explicitly calculating electron density. This algorithm can be applied to a wide variety of molecular systems.

Crystal structure of CO-bound cytochrome c oxidase determined by serial femtosecond X-ray crystallography at room temperature.

Cytochrome c oxidase (CcO), the terminal enzyme in the electron transfer chain, translocates protons across the inner mitochondrial membrane by harnessing the free energy generated by the reduction of oxygen to water. Several redox-coupled proton translocation mechanisms have been proposed, but they lack confirmation, in part from the absence of reliable structural information due to radiation damage artifacts caused by the intense synchrotron radiation. Here we report the room temperature, neutral pH (6.8), damage-free structure of bovine CcO (bCcO) in the carbon monoxide (CO)-bound state at a resolution of 2.3 Å, obtained by serial femtosecond X-ray crystallography (SFX) with an X-ray free electron laser. As a comparison, an equivalent structure was obtained at a resolution of 1.95 Å, from data collected at a synchrotron light source. In the SFX structure, the CO is coordinated to the heme a3 iron atom, with a bent Fe-C-O angle of ∼142°. In contrast, in the synchrotron structure, the Fe-CO bond is cleaved; CO relocates to a new site near CuB, which, in turn, moves closer to the heme a3 iron by ∼0.38 Å. Structural comparison reveals that ligand binding to the heme a3 iron in the SFX structure is associated with an allosteric structural transition, involving partial unwinding of the helix-X between heme a and a3, thereby establishing a communication linkage between the two heme groups, setting the stage for proton translocation during the ensuing redox chemistry.

Structural consequences of transforming growth factor beta-1 activation from near-therapeutic X-ray doses.

Dissociation of transforming growth factor beta-1 (TGFß-1) from the inhibitory protein latency-associated peptide (LAP) can occur from low doses of X-ray irradiation of the LAP-TGFß-1 complex, resulting in the activation of TGFß-1, and can have health-related consequences. Using the tools and knowledge developed in the study of radiation damage in the crystallographic setting, small-angle X-ray scattering (SAXS) and complementary techniques suggest an activation process that is initiated but not driven by the initial X-ray exposure. LAP is revealed to be extended when not bound to TGFß-1 and has a different structural conformation compared to the bound state. These studies pave the way for the structural understanding of systems impacted at therapeutic X-ray doses and show the potential impact of radiation damage studies beyond their original intent.

Structural and Functional Characterization of Aerobactin Synthetase IucA from a Hypervirulent Pathotype of Klebsiella pneumoniae.

Iron is a vital mineral nutrient required by virtually all life forms to prosper; pathogenic bacteria are no exception. Despite the abundance of iron within the human host, highly regulated iron physiology can result in exceedingly low levels of iron bioavailable to prospective invading bacteria. To combat this scarcity of iron, many pathogenic bacteria have acquired specific and efficient iron acquisition systems, which allow them to thrive in iron-deficient host environments. One of the more prominent bacterial iron acquisition systems involves the synthesis, secretion, and reuptake of small-molecule iron chelators known as siderophores. Aerobactin, a citrate-hydroxamate siderophore originally isolated nearly 50 years ago, is produced by a number of pathogenic Gram-negative bacteria. Aerobactin has recently been demonstrated to play a pivotal role in mediating the enhanced virulence of a particularly invasive pathotype of Klebsiella pneumoniae (hvKP). Toward further understanding of this key virulence factor, we report the structural and functional characterization of aerobactin synthetase IucA from a strain of hvKP. The X-ray crystal structures of unliganded and ATP-bound forms of IucA were solved, forming the foundation of our structural analysis. Small angle X-ray scattering (SAXS) data suggest that, unlike its closest structurally characterized homologues, IucA adopts a tetrameric assembly in solution. Finally, we employed activity assays to investigate the substrate specificity and determine the apparent steady-state kinetic parameters of IucA.

A hybrid NMR/SAXS-based approach for discriminating oligomeric protein interfaces using Rosetta.

Oligomeric proteins are important targets for structure determination in solution. While in most cases the fold of individual subunits can be determined experimentally, or predicted by homology-based methods, protein-protein interfaces are challenging to determine de novo using conventional NMR structure determination protocols. Here we focus on a member of the bet-V1 superfamily, Aha1 from Colwellia psychrerythraea. This family displays a broad range of crystallographic interfaces none of which can be reconciled with the NMR and SAXS data collected for Aha1. Unlike conventional methods relying on a dense network of experimental restraints, the sparse data are used to limit conformational search during optimization of a physically realistic energy function. This work highlights a new approach for studying minor conformational changes due to structural plasticity within a single dimeric interface in solution.

The accurate assessment of small-angle X-ray scattering data.

Small-angle X-ray scattering (SAXS) has grown in popularity in recent times with the advent of bright synchrotron X-ray sources, powerful computational resources and algorithms enabling the calculation of increasingly complex models. However, the lack of standardized data-quality metrics presents difficulties for the growing user community in accurately assessing the quality of experimental SAXS data. Here, a series of metrics to quantitatively describe SAXS data in an objective manner using statistical evaluations are defined. These metrics are applied to identify the effects of radiation damage, concentration dependence and interparticle interactions on SAXS data from a set of 27 previously described targets for which high-resolution structures have been determined via X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. The studies show that these metrics are sufficient to characterize SAXS data quality on a small sample set with statistical rigor and sensitivity similar to or better than manual analysis. The development of data-quality analysis strategies such as these initial efforts is needed to enable the accurate and unbiased assessment of SAXS data quality.

Structural conservation of an ancient tRNA sensor in eukaryotic glutaminyl-tRNA synthetase.

In all organisms, aminoacyl tRNA synthetases covalently attach amino acids to their cognate tRNAs. Many eukaryotic tRNA synthetases have acquired appended domains, whose origin, structure and function are poorly understood. The N-terminal appended domain (NTD) of glutaminyl-tRNA synthetase (GlnRS) is intriguing since GlnRS is primarily a eukaryotic enzyme, whereas in other kingdoms Gln-tRNA(Gln) is primarily synthesized by first forming Glu-tRNA(Gln), followed by conversion to Gln-tRNA(Gln) by a tRNA-dependent amidotransferase. We report a functional and structural analysis of the NTD of Saccharomyces cerevisiae GlnRS, Gln4. Yeast mutants lacking the NTD exhibit growth defects, and Gln4 lacking the NTD has reduced complementarity for tRNA(Gln) and glutamine. The 187-amino acid Gln4 NTD, crystallized and solved at 2.3 Å resolution, consists of two subdomains, each exhibiting an extraordinary structural resemblance to adjacent tRNA specificity-determining domains in the GatB subunit of the GatCAB amidotransferase, which forms Gln-tRNA(Gln). These subdomains are connected by an apparent hinge comprised of conserved residues. Mutation of these amino acids produces Gln4 variants with reduced affinity for tRNA(Gln), consistent with a hinge-closing mechanism proposed for GatB recognition of tRNA. Our results suggest a possible origin and function of the NTD that would link the phylogenetically diverse mechanisms of Gln-tRNA(Gln) synthesis.

Reconstruction of 3D density from solution scattering.

Ab initio modeling methods have proven to be powerful means of interpreting solution scattering data. In the absence of atomic models, or complementary to them, ab initio modeling approaches can be used for generating low-resolution particle envelopes using only solution scattering profiles. Recently, a new ab initio reconstruction algorithm has been introduced to the scientific community, called DENSS. DENSS is unique among ab initio modeling algorithms in that it solves the inverse scattering problem, i.e., the 1D scattering intensities are directly used to determine the 3D particle density. The reconstruction of particle density has several advantages over conventional uniform density modeling approaches, including the ability to reconstruct a much wider range of particle types and the ability to visualize low-resolution density fluctuations inside the particle envelope. In this chapter we will discuss the theory behind this new approach, how to use DENSS, and how to interpret the results. Several examples with experimental and simulated data will be provided.

Fitting high-resolution electron density maps from atomic models to solution scattering data.

Solution scattering techniques, such as small and wide-angle X-ray scattering (SWAXS), provide valuable insights into the structure and dynamics of biological macromolecules in solution. In this study, we present an approach to accurately predict solution X-ray scattering profiles at wide angles from atomic models by generating high-resolution electron density maps. Our method accounts for the excluded volume of bulk solvent by calculating unique adjusted atomic volumes directly from the atomic coordinates. This approach eliminates the need for a free fitting parameter commonly used in existing algorithms, resulting in improved accuracy of the calculated SWAXS profile. An implicit model of the hydration shell is generated which uses the form factor of water. Two parameters, namely the bulk solvent density and the mean hydration shell contrast, are adjusted to best fit the data. Results using eight publicly available SWAXS profiles show high quality fits to the data. In each case, the optimized parameter values show small adjustments demonstrating that the default values are close to the true solution. Disabling parameter optimization results in a significant improvement of the calculated scattering profiles compared to the leading software. The algorithm is computationally efficient, showing more than tenfold reduction in execution time compared to the leading software. The algorithm is encoded in a command line script called denss.pdb2mrc.py and is available open source as part of the DENSS v1.7.0 software package (https://github.com/tdgrant1/denss). In addition to improving the ability to compare atomic models to experimental SWAXS data, these developments pave the way for increasing the accuracy of modeling algorithms utilizing SWAXS data while decreasing the risk of overfitting.

RNA structures and dynamics with Å resolution revealed by x-ray free-electron lasers.

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free-electron laser sources to reveal the formation and ready identification of angstrom-scale features in structured and unstructured RNAs. Previously unrecognized structural signatures of RNA secondary and tertiary structures are identified through wide-angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base-paired intermediate to assume a triple-helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. This method may help to rapidly characterize and identify structural elements in nucleic acids in both equilibrium and time-resolved experiments.

RNA structures and dynamics with Å resolution revealed by x-ray free electron lasers.

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free electron laser sources to reveal the formation and ready identification of Å scale features in structured and unstructured RNAs. New structural signatures of RNA secondary and tertiary structures are identified through wide angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base paired intermediate to assume a triple helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. In addition to understanding how RNA triplexes form and thereby function as dynamic signaling elements, this new method can vastly increase the rate of structure determination for these biologically essential, but mostly uncharacterized macromolecules.

Room-temperature structural studies of SARS-CoV-2 protein NendoU with an X-ray free-electron laser.

NendoU from SARS-CoV-2 is responsible for the virus's ability to evade the innate immune system by cleaving the polyuridine leader sequence of antisense viral RNA. Here we report the room-temperature structure of NendoU, solved by serial femtosecond crystallography at an X-ray free-electron laser to 2.6 Å resolution. The room-temperature structure provides insight into the flexibility, dynamics, and other intrinsic properties of NendoU, with indications that the enzyme functions as an allosteric switch. Functional studies examining cleavage specificity in solution and in crystals support the uridine-purine cleavage preference, and we demonstrate that enzyme activity is fully maintained in crystal form. Optimizing the purification of NendoU and identifying suitable crystallization conditions set the benchmark for future time-resolved serial femtosecond crystallography studies. This could advance the design of antivirals with higher efficacy in treating coronaviral infections, since drugs that block allosteric conformational changes are less prone to drug resistance.

Modular droplet injector for sample conservation providing new structural insight for the conformational heterogeneity in the disease-associated NQO1 enzyme.

Droplet injection strategies are a promising tool to reduce the large amount of sample consumed in serial femtosecond crystallography (SFX) measurements at X-ray free electron lasers (XFELs) with continuous injection approaches. Here, we demonstrate a new modular microfluidic droplet injector (MDI) design that was successfully applied to deliver microcrystals of the human NAD(P)H:quinone oxidoreductase 1 (NQO1) and phycocyanin. We investigated droplet generation conditions through electrical stimulation for both protein samples and implemented hardware and software components for optimized crystal injection at the Macromolecular Femtosecond Crystallography (MFX) instrument at the Stanford Linac Coherent Light Source (LCLS). Under optimized droplet injection conditions, we demonstrate that up to 4-fold sample consumption savings can be achieved with the droplet injector. In addition, we collected a full data set with droplet injection for NQO1 protein crystals with a resolution up to 2.7 Å, leading to the first room-temperature structure of NQO1 at an XFEL. NQO1 is a flavoenzyme associated with cancer, Alzheimer's and Parkinson's disease, making it an attractive target for drug discovery. Our results reveal for the first time that residues Tyr128 and Phe232, which play key roles in the function of the protein, show an unexpected conformational heterogeneity at room temperature within the crystals. These results suggest that different substates exist in the conformational ensemble of NQO1 with functional and mechanistic implications for the enzyme's negative cooperativity through a conformational selection mechanism. Our study thus demonstrates that microfluidic droplet injection constitutes a robust sample-conserving injection method for SFX studies on protein crystals that are difficult to obtain in amounts necessary for continuous injection, including the large sample quantities required for time-resolved mix-and-inject studies.

Describing small-angle scattering profiles by a limited set of intensities.

Small-angle scattering (SAS) probes the size and shape of particles at low resolution through the analysis of the scattering of X-rays or neutrons passing through a solution of particles. One approach to extracting structural information from SAS data is the indirect Fourier transform (IFT). The IFT approach parameterizes the real-space pair distribution function [P(r)] of a particle using a set of basis functions, which simultaneously determines the scattering profile [I(q)] using corresponding reciprocal-space basis functions. This article presents an extension of an IFT algorithm proposed by Moore [J. Appl. Cryst. (1980), 13, 168-175] which used a trigonometric series to describe the basis functions, where the real-space and reciprocal-space basis functions are Fourier mates. An equation is presented relating the Moore coefficients to the intensities of the SAS profile at specific positions, as well as a series of new equations that describe the size and shape parameters of a particle from this distinct set of intensity values. An analytical real-space regularizer is derived to smooth the P(r) curve and ameliorate systematic deviations caused by series termination. Regularization is commonly used in IFT methods though not described in Moore's original approach, which is particularly susceptible to such effects. The algorithm is provided as a script, denss.f it_data.py, as part of the DENSS software package for SAS, which includes both command line and interactive graphical interfaces. Results of the program using experimental data show that it is as accurate as, and often more accurate than, existing tools.