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.

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.

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.

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 Å 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.

Electrically stimulated droplet injector for reduced sample consumption in serial crystallography

With advances in X-ray free-electron lasers (XFELs), serial femtosecond crystallography (SFX) has enabled the static and dynamic structure determination for challenging proteins such as membrane protein complexes. In SFX with XFELs, the crystals are typically destroyed after interacting with a single XFEL pulse. Therefore, thousands of new crystals must be sequentially introduced into the X-ray beam to collect full data sets. Because of the serial nature of any SFX experiment, up to 99% of the sample delivered to the X-ray beam during its "off-time" between X-ray pulses is wasted due to the intrinsic pulsed nature of all current XFELs. To solve this major problem of large and often limiting sample consumption, we report on improvements of a revolutionary sample-saving method that is compatible with all current XFELs. We previously reported 3D-printed injection devices coupled with gas dynamic virtual nozzles (GDVNs) capable of generating samples containing droplets segmented by an immiscible oil phase for jetting crystal-laden droplets into the path of an XFEL. Here, we have further improved the device design by including metal electrodes inducing electrowetting effects for improved control over droplet generation frequency to stimulate the droplet release to matching the XFEL repetition rate by employing an electrical feedback mechanism. We report the improvements in this electrically triggered segmented flow approach for sample conservation in comparison with a continuous GDVN injection using the microcrystals of lysozyme and 3-deoxy-D-manno-octulosonate 8-phosphate synthase and report the segmented flow approach for sample injection applied at the Macromolecular Femtosecond Crystallography instrument at the Linear Coherent Light Source for the first time.

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.

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

An indirect Fourier transform method is presented which describes a solution scattering profile from a reduced set of intensities. Equations are derived to fit the experimental profile using least squares and to calculate commonly used size and shape parameters directly from the reduced set of intensities, along with associated uncertainties. An analytical equation is derived enabling regularization of the real-space pair distribution function. Convenient software is provided to perform all described calculations.

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.fit_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.

Structural and biophysical properties of FopA, a major outer membrane protein of Francisella tularensis

Francisella tularensis is an extremely infectious pathogen and a category A bioterrorism agent. It causes the highly contagious zoonosis, Tularemia. Currently, FDA approved vaccines against tularemia are unavailable. F. tularensis outer membrane protein A (FopA) is a well-studied virulence determinant and protective antigen against tularemia. It is a major outer membrane protein (Omp) of F. tularensis. However, FopA-based therapeutic intervention is hindered due to lack of complete structural information for membrane localized mature FopA. In our study, we established recombinant expression, monodisperse purification, crystallization and X-ray diffraction (~6.5 Å) of membrane localized mature FopA. Further, we performed bioinformatics and biophysical experiments to unveil its structural organization in the outer membrane. FopA consists of 393 amino acids and has less than 40% sequence identity to known bacterial Omps. Using comprehensive sequence alignments and structure predictions together with existing partial structural information, we propose a two-domain organization for FopA. Circular dichroism spectroscopy and heat modifiability assay confirmed FopA has a β-barrel domain consistent with alphafold2’s prediction of an eight stranded β-barrel at the N-terminus. Small angle X-ray scattering (SAXS) and native-polyacrylamide gel electrophoresis revealed FopA purified in detergent micelles is predominantly dimeric. Molecular density derived from SAXS at 31 Å shows putative dimeric N-terminal β-barrels surrounded by detergent corona and connected to C-terminal domains via flexible linker. Disorder analysis predicts N- and C-terminal domains are interspersed by a long intrinsically disordered region and alphafold2 predicts this region to be largely unstructured. Taken together, we propose a dimeric, two-domain organization of FopA in the outer membrane: the N-terminal β-barrel is membrane embedded, provides dimerization interface and tethers to membrane extrinsic C-terminal domain via long flexible linker. Structure determination of membrane localized mature FopA is essential to understand its role in pathogenesis and develop anti-tularemia therapeutics. Our results pave the way towards it.

Synchronous RNA conformational changes trigger ordered phase transitions in crystals

Time-resolved studies of biomacromolecular crystals have been limited to systems involving only minute conformational changes within the same lattice. Ligand-induced changes greater than several angstroms, however, are likely to result in solid-solid phase transitions, which require a detailed understanding of the mechanistic interplay between conformational and lattice transitions. Here we report the synchronous behavior of the adenine riboswitch aptamer RNA in crystal during ligand-triggered isothermal phase transitions. Direct visualization using polarized video microscopy and atomic force microscopy shows that the RNA molecules undergo cooperative rearrangements that maintain lattice order, whose cell parameters change distinctly as a function of time. The bulk lattice order throughout the transition is further supported by time-resolved diffraction data from crystals using an X-ray free electron laser. The synchronous molecular rearrangements in crystal provide the physical basis for studying large conformational changes using time-resolved crystallography and micro/nanocrystals.

Time-resolved crystallography (TRX) is used for monitoring only small conformational changes of biomacromolecules within the same lattice. Here, the authors report the interplay between synchronous molecular rearrangements and lattice phase transitions in RNA crystals, providing the basis for the investigation of large conformational changes using TRX.

Segmented flow generator for serial crystallography at the European X-ray free electron laser

Serial femtosecond crystallography (SFX) with X-ray free electron lasers (XFELs) allows structure determination of membrane proteins and time-resolved crystallography. Common liquid sample delivery continuously jets the protein crystal suspension into the path of the XFEL, wasting a vast amount of sample due to the pulsed nature of all current XFEL sources. The European XFEL (EuXFEL) delivers femtosecond (fs) X-ray pulses in trains spaced 100 ms apart whereas pulses within trains are currently separated by 889 ns. Therefore, continuous sample delivery via fast jets wastes >99% of sample. Here, we introduce a microfluidic device delivering crystal laden droplets segmented with an immiscible oil reducing sample waste and demonstrate droplet injection at the EuXFEL compatible with high pressure liquid delivery of an SFX experiment. While achieving ~60% reduction in sample waste, we determine the structure of the enzyme 3-deoxy-D-manno-octulosonate-8-phosphate synthase from microcrystals delivered in droplets revealing distinct structural features not previously reported.

Membrane protein megahertz crystallography at the European XFEL

The world’s first superconducting megahertz repetition rate hard X-ray free-electron laser (XFEL), the European XFEL, began operation in 2017, featuring a unique pulse train structure with 886 ns between pulses. With its rapid pulse rate, the European XFEL may alleviate some of the increasing demand for XFEL beamtime, particularly for membrane protein serial femtosecond crystallography (SFX), leveraging orders-of-magnitude faster data collection. Here, we report the first membrane protein megahertz SFX experiment, where we determined a 2.9 Å-resolution SFX structure of the large membrane protein complex, Photosystem I, a > 1 MDa complex containing 36 protein subunits and 381 cofactors. We address challenges to megahertz SFX for membrane protein complexes, including growth of large quantities of crystals and the large molecular and unit cell size that influence data collection and analysis. The results imply that megahertz crystallography could have an important impact on structure determination of large protein complexes with XFELs.

The European X-ray free-electron laser (EuXFEL) in Hamburg is the first XFEL with a megahertz repetition rate. Here the authors present the 2.9 Å structure of the large membrane protein complex Photosystem I from T. elongatus that was determined at the EuXFEL.

3D printed droplet generation devices for serial femtosecond crystallography enabled by surface coating

The role of surface wetting properties and their impact on the performance of 3D printed microfluidic droplet generation devices for serial femtosecond crystallography (SFX) are reported. SFX is a novel crystallography method enabling structure determination of proteins at room temperature with atomic resolution using X-ray free-electron lasers (XFELs). In SFX, protein crystals in their mother liquor are delivered and intersected with a pulsed X-ray beam using a liquid jet injector. Owing to the pulsed nature of the X-ray beam, liquid jets tend to waste the vast majority of injected crystals, which this work aims to overcome with the delivery of aqueous protein crystal suspension droplets segmented by an oil phase. For this purpose, 3D printed droplet generators that can be easily customized for a variety of XFEL measurements have been developed. The surface properties, in particular the wetting properties of the resist materials compatible with the employed two-photon printing technology, have so far not been characterized extensively, but are crucial for stable droplet generation. This work investigates experimentally the effectiveness and the long-term stability of three different surface treatments on photoresist films and glass as models for our 3D printed droplet generator and the fused silica capillaries employed in the other fluidic components of an SFX experiment. Finally, the droplet generation performance of an assembly consisting of the 3D printed device and fused silica capillaries is examined. Stable and reproducible droplet generation was achieved with a fluorinated surface coating which also allowed for robust downstream droplet delivery. Experimental XFEL diffraction data of crystals formed from the large membrane protein complex photosystem I demonstrate the full compatibility of the new injection method with very fragile membrane protein crystals and show that successful droplet generation of crystal-laden aqueous droplets intersected by an oil phase correlates with increased crystal hit rates.

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.

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.

X-ray Emission Spectroscopy at X-ray Free Electron Lasers: Limits to Observation of the Classical Spectroscopic Response for Electronic Structure Analysis

X-ray free electron lasers (XFELs) provide ultrashort intense X-ray pulses suitable to probe electron dynamics but can also induce a multitude of nonlinear excitation processes. These affect spectroscopic measurements and interpretation, particularly for upcoming brighter XFELs. Here we identify and discuss the limits to observing classical spectroscopy, where only one photon is absorbed per atom for a Mn²⁺ in a light element (O, C, H) environment. X-ray emission spectroscopy (XES) with different incident photon energies, pulse intensities, and pulse durations is presented. A rate equation model based on sequential ionization and relaxation events is used to calculate populations of multiply ionized states during a single pulse and to explain the observed X-ray induced spectral lines shifts. This model provides easy estimation of spectral shifts, which is essential for experimental designs at XFELs and illustrates that shorter X-ray pulses will not overcome sequential ionization but can reduce electron cascade effects.

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.

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.

Femtosecond structural dynamics drives the trans/cis isomerization in photoactive yellow protein

A variety of organisms have evolved mechanisms to detect and respond to light, in which the response is mediated by protein structural changes after photon absorption. The initial step is often the photoisomerization of a conjugated chromophore. Isomerization occurs on ultrafast time scales and is substantially influenced by the chromophore environment. Here we identify structural changes associated with the earliest steps in the trans-to-cis isomerization of the chromophore in photoactive yellow protein. Femtosecond hard x-ray pulses emitted by the Linac Coherent Light Source were used to conduct time-resolved serial femtosecond crystallography on photoactive yellow protein microcrystals over a time range from 100 femtoseconds to 3 picoseconds to determine the structural dynamics of the photoisomerization reaction.

The structure of the PanD/PanZ protein complex reveals negative feedback regulation of pantothenate biosynthesis by coenzyme A

Coenzyme A (CoA) is an ubiquitous and essential cofactor, synthesized from the precursor pantothenate. Vitamin biosynthetic pathways are normally tightly regulated, including the pathway from pantothenate to CoA. However, no regulation of pantothenate biosynthesis has been identified. We have recently described an additional component in the pantothenate biosynthetic pathway, PanZ, which promotes the activation of the zymogen, PanD, to form aspartate α-decarboxylase (ADC) in a CoA-dependent manner. Here we report the structure of PanZ in complex with PanD, which reveals the structural basis for the CoA dependence of this interaction and activation. In addition, we show that PanZ acts as a CoA-dependent inhibitor of ADC catalysis. This inhibitory effect can effectively regulate the biosynthetic pathway to pantothenate, and thereby also regulate CoA biosynthesis. This represents a previously unobserved mode of metabolic regulation whereby a cofactor-utilizing protein negatively regulates the biosynthesis of the same cofactor.

•

Structure of the PanD-PanZ.AcCoA complex is reported at a resolution of 1.6 Å

•

Binding of AcCoA to PanZ is required to form the PanZ/PanD interface

•

PanZ.AcCoA activates PanD via selection of a reactive conformation of PanD

•

PanZ.AcCoA inhibits the activated enzyme, regulating pantothenate biosynthesis

Microfluidic sorting of protein nanocrystals by size for X-ray free-electron laser diffraction

The advent and application of the X-ray free-electron laser (XFEL) has uncovered the structures of proteins that could not previously be solved using traditional crystallography. While this new technology is powerful, optimization of the process is still needed to improve data quality and analysis efficiency. One area is sample heterogeneity, where variations in crystal size (among other factors) lead to the requirement of large data sets (and thus 10-100 mg of protein) for determining accurate structure factors. To decrease sample dispersity, we developed a high-throughput microfluidic sorter operating on the principle of dielectrophoresis, whereby polydisperse particles can be transported into various fluid streams for size fractionation. Using this microsorter, we isolated several milliliters of photosystem I nanocrystal fractions ranging from 200 to 600 nm in size as characterized by dynamic light scattering, nanoparticle tracking, and electron microscopy. Sorted nanocrystals were delivered in a liquid jet via the gas dynamic virtual nozzle into the path of the XFEL at the Linac Coherent Light Source. We obtained diffraction to ∼4 Å resolution, indicating that the small crystals were not damaged by the sorting process. We also observed the shape transforms of photosystem I nanocrystals, demonstrating that our device can optimize data collection for the shape transform-based phasing method. Using simulations, we show that narrow crystal size distributions can significantly improve merged data quality in serial crystallography. From this proof-of-concept work, we expect that the automated size-sorting of protein crystals will become an important step for sample production by reducing the amount of protein needed for a high quality final structure and the development of novel phasing methods that exploit inter-Bragg reflection intensities or use variations in beam intensity for radiation damage-induced phasing. This method will also permit an analysis of the dependence of crystal quality on crystal size.

Serial femtosecond X-ray diffraction of enveloped virus microcrystals

Serial femtosecond crystallography (SFX) using X-ray free-electron lasers has produced high-resolution, room temperature, time-resolved protein structures. We report preliminary SFX of Sindbis virus, an enveloped icosahedral RNA virus with ∼700 Å diameter. Microcrystals delivered in viscous agarose medium diffracted to ∼40 Å resolution. Small-angle diffuse X-ray scattering overlaid Bragg peaks and analysis suggests this results from molecular transforms of individual particles. Viral proteins undergo structural changes during entry and infection, which could, in principle, be studied with SFX. This is an important step toward determining room temperature structures from virus microcrystals that may enable time-resolved studies of enveloped viruses.

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

A set of quantitative techniques is suggested for assessing SAXS data quality. These are applied in the form of a script, SAXStats, to a test set of 27 proteins, showing that these techniques are more sensitive than manual assessment of data quality.

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.

Comparing chemistry to outcome: the development of a chemical distance metric, coupled with clustering and hierarchal visualization applied to macromolecular crystallography

Many bioscience fields employ high-throughput methods to screen multiple biochemical conditions. The analysis of these becomes tedious without a degree of automation. Crystallization, a rate limiting step in biological X-ray crystallography, is one of these fields. Screening of multiple potential crystallization conditions (cocktails) is the most effective method of probing a proteins phase diagram and guiding crystallization but the interpretation of results can be time-consuming. To aid this empirical approach a cocktail distance coefficient was developed to quantitatively compare macromolecule crystallization conditions and outcome. These coefficients were evaluated against an existing similarity metric developed for crystallization, the C6 metric, using both virtual crystallization screens and by comparison of two related 1,536-cocktail high-throughput crystallization screens. Hierarchical clustering was employed to visualize one of these screens and the crystallization results from an exopolyphosphatase-related protein from Bacteroides fragilis, (BfR192) overlaid on this clustering. This demonstrated a strong correlation between certain chemically related clusters and crystal lead conditions. While this analysis was not used to guide the initial crystallization optimization, it led to the re-evaluation of unexplained peaks in the electron density map of the protein and to the insertion and correct placement of sodium, potassium and phosphate atoms in the structure. With these in place, the resulting structure of the putative active site demonstrated features consistent with active sites of other phosphatases which are involved in binding the phosphoryl moieties of nucleotide triphosphates. The new distance coefficient, CD_coeff, appears to be robust in this application, and coupled with hierarchical clustering and the overlay of crystallization outcome, reveals information of biological relevance. While tested with a single example the potential applications related to crystallography appear promising and the distance coefficient, clustering, and hierarchal visualization of results undoubtedly have applications in wider fields.

A new view on crystal harvesting

X-ray crystallography typically requires the mounting of crystals, which can make the sample difficult to manipulate when it is small and the microscope objective is close to the crystallization plate. By simply moving the objective to the bottom of a clear crystallization plate (inverting the normal view), crystals were able to be manipulated and harvested from wells having a 0.9 mm diameter and 5.0 mm depth. The mounting system enabled the structural solution of the 187 amino acid N-terminal domain of Saccharomyces cerevisiae glutaminyl-tRNA synthetase from crystals that appeared during high-throughput screening but proved recalcitrant to scale-up and optimization. While not a general mounting solution, the simple expedient of removing the objective lens from the area where manipulation and harvesting occur greatly facilitates the manual, or even automated, process.

The structure of yeast glutaminyl-tRNA synthetase and modeling of its interaction with tRNA

Eukaryotic glutaminyl-tRNA synthetase (GlnRS) contains an appended N-terminal domain (NTD) whose precise function is unknown. Although GlnRS structures from two prokaryotic species are known, no eukaryotic GlnRS structure has been reported. Here we present the first crystallographic structure of yeast GlnRS, finding that the structure of the C-terminal domain is highly similar to Escherichia coli GlnRS but that 214 residues, including the NTD, are crystallographically disordered. We present a model of the full-length enzyme in solution, using the structures of the C-terminal domain, and the isolated NTD, with small-angle X-ray scattering data of the full-length molecule. We proceed to model the enzyme bound to tRNA, using the crystallographic structures of GatCAB and GlnRS-tRNA complex from bacteria. We contrast the tRNA-bound model with the tRNA-free solution state and perform molecular dynamics on the full-length GlnRS-tRNA complex, which suggests that tRNA binding involves the motion of a conserved hinge in the NTD.

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.

Small angle X-ray scattering as a complementary tool for high-throughput structural studies

Structural crystallography and nuclear magnetic resonance (NMR) spectroscopy are the predominant techniques for understanding the biological world on a molecular level. Crystallography is constrained by the ability to form a crystal that diffracts well and NMR is constrained to smaller proteins. Although powerful techniques, they leave many soluble, purified structurally uncharacterized protein samples. Small angle X-ray scattering (SAXS) is a solution technique that provides data on the size and multiple conformations of a sample, and can be used to reconstruct a low-resolution molecular envelope of a macromolecule. In this study, SAXS has been used in a high-throughput manner on a subset of 28 proteins, where structural information is available from crystallographic and/or NMR techniques. These crystallographic and NMR structures were used to validate the accuracy of molecular envelopes reconstructed from SAXS data on a statistical level, to compare and highlight complementary structural information that SAXS provides, and to leverage biological information derived by crystallographers and spectroscopists from their structures. All the ab initio molecular envelopes calculated from the SAXS data agree well with the available structural information. SAXS is a powerful albeit low-resolution technique that can provide additional structural information in a high-throughput and complementary manner to improve the functional interpretation of high-resolution structures.

Col2-GFP reporter mouse--a new tool to study skeletal development

Transgenic mice were generated that harbor a Col2-GFP reporter that marks chondrocytes and their immediate precursors during skeletal development. Cells engaged in chondrogenesis were identified by conventional fluorescence microscopy and confocal optical sectioning within their native environments in live embryos and in thick tissue slices. The use of these mice offers a novel approach for studying the role of chondrocytes in skeletal development.

Epithelial cell dynamics in rabbit cecum and proximal colon P1

The large intestine of mammals has long been viewed as an osmoregulatory organ, and evidence suggests that fluid and solute transport mechanisms within the intestine are heterogeneous, varying depending on the particular segment involved. Variations in function are often matched by morphological correlates, but despite the widespread use of rabbit large intestine as an experimental model, there is a lack of knowledge about the cellular makeup and dynamics in the colonic mucosal epithelium. The presence of mitotic figures and immunohistochemical localization of proliferating cell nuclear antigen (PCNA) were used to identify the proliferative zone(s). Cellular migration patterns were determined through the use of the thymidine analog 5-bromo-2-deoxyuridine (BrdU) over a 24-, 48-, and 72-hr period. Apoptotic nuclei were identified utilizing terminal deoxynucleotidyl transferase d-UTP nick-end labeling (TUNEL). Both cecum and the initial portion of the proximal colon (P1) exhibited a proliferative zone at or near the crypt base, and migration proceeded upwards toward the surface epithelium lining the intestinal lumen, where apoptosis occurred Turnover time of crypt columnar cells was determined to be about 3 days; that of mucous cells was estimated to be about 5 weeks. Rabbit cecum and proximal colon P1 are similar in their cellular morphology and epithelial cell kinetics. In both, the major proliferative zone is located at or near the crypt base, from which crypt columnar cells migrate toward the lumenal surface epithelium over a period of 3 days. Goblet cell turnover rate is much slower than that of columnar cells.

Col2-GFP reporter marks chondrocyte lineage and chondrogenesis during mouse skeletal development

Mice were generated in which a Col2-GFP transgene serves as a reporter for the chondrocyte lineage and for chondrogenesis in live embryos and newborn pups. Cells actively engaged in chondrogenesis were identified by confocal optical sectioning within their native environments in embryos and in thick tissue slices. Chondrocytes exhibiting GFP fluorescence were purified from rib cages by high-speed cell sorting of crude cell suspensions. Intensity of fluorescence correlated with biosynthesis of procollagen II in these cells. The use of these mice and their cells provides a novel approach for studying chondrocyte differentiation and chondrogenesis during skeletal development.

Proliferation of goblet cells and vacuolated cells in the rabbit distal colon

Previous studies of colonic epithelial cell kinetics in mice and rats revealed a pattern similar to small intestine, where basally located stem cells proliferate, differentiating as they migrate towards the surface epithelium. Vacuolated and goblet cells are assumed to co-migrate at the same rate. The present study indicates that rabbit distal colon has more complicated epithelial cell kinetics. The zone of proliferation was detected immunohistochemically using proliferating cell nuclear antigen (PCNA) and confirmed with the use of colchicine to arrest dividing cells in metaphase. Migrating cells were tracked from the zero-hour position (PCNA labeling, mitosis) to positions 24, 48, 72 hrs by monitoring cell migration with the thymidine analog 5-Bromo-2-Deoxyuridine (BrdU). PCNA revealed a major proliferative zone in the upper third of the crypt column and the presence of mitotic figures after colchicine corroborated these results. Differentiated vacuolated cell proliferation was detected at three crypt sites: base, middle, and top of the crypt, while columnar cells arose from a population of dividing cells at the top of the crypt. Turnover of columnar and vacuolated cells occurred within 72 hrs. Goblet cells exhibited maximal proliferation at the crypt base and migrated at a much slower rate than the other cell types. In rabbit distal colon, populations of proliferating cells exist at multiple levels of the crypt column. Vacuolated and goblet cells differ in their labeling indices and migration rates, suggesting that the two cell types arise and migrate independently.

Epidermal growth factor modifies cell cycle control in A431 human squamous carcinoma cells damaged by ionizing radiation

Epidermal growth factor (EGF) has been shown to radiosensitize A431 and other human squamous carcinoma cells with high numbers of surface EGF receptors. In this study of the mechanistic basis of EGF-induced radiosensitization, both EGF and ionizing radiation caused G1 phase arrests in cycling A431 cells, but only radiation caused a G2-M arrest. However, EGF enhanced the magnitude of this G2-M arrest, suggesting an interaction of signaling pathways involved in cellular responses to EGF and radiation damage. EGF and radiation also uniquely perturbed cyclin A and B1 mRNA levels during the time of maximum radiation-induced G2-M arrest. The effects of EGF on G2-M events probably originated in cells in G1. It is possible that aberrant EGF signal transduction in human squamous carcinoma cells may be exploited as a novel strategy for radiotherapy.