Resolution of structural heterogeneity in dynamic crystallography

Spin-Coupled Electron Densities of Iron-Sulfur Cluster Imaged by In Situ Serial Laue Diffraction

Iron-sulfur clusters are inorganic cofactors found in many proteins involved in fundamental biological processes including DNA processing. The prokaryotic DNA repair enzyme PhrB, a member of the protein family of cryptochromes and photolyases, carries a four-iron-four-sulfur cluster [4Fe4S] in addition to the catalytic cofactor flavin adenine dinucleotide (FAD) and a second pigment 6,7-dimethyl-8-ribityllumazine (DMRL). The light-induced redox reactions of this multi-cofactor protein complex were recently shown as two interdependent photoreductions of FAD and DMRL mediated by the [4Fe4S] cluster functioning as an electron cache to hold a fine balance of electrons. Here, we apply the more traditional temperature-scan cryo-trapping technique in protein crystallography and the newly developed technology of in situ serial Laue diffraction at room temperature. These diffraction methods in dynamic crystallography enable us to capture strong signals of electron density changes in the [4Fe4S] cluster that depict quantized electronic movements. The mixed valence layers of the [4Fe4S] cluster due to spin coupling and their dynamic responses to light illumination are observed directly in our difference maps between its redox states. These direct observations of the quantum effects in a protein bound iron-sulfur cluster have thus opened a window into the mechanistic understanding of metal clusters in biological systems.

Photoinduced isomerization sampling of retinal in bacteriorhodopsin

Photoisomerization of retinoids inside a confined protein pocket represents a critical chemical event in many important biological processes from animal vision, nonvisual light effects, to bacterial light sensing and harvesting. Light-driven proton pumping in bacteriorhodopsin entails exquisite electronic and conformational reconfigurations during its photocycle. However, it has been a major challenge to delineate transient molecular events preceding and following the photoisomerization of the retinal from noisy electron density maps when varying populations of intermediates coexist and evolve as a function of time. Here, I report several distinct early photoproducts deconvoluted from the recently observed mixtures in time-resolved serial crystallography. This deconvolution substantially improves the quality of the electron density maps, hence demonstrates that the all-trans retinal undergoes extensive isomerization sampling before it proceeds to the productive 13-cis configuration. Upon light absorption, the chromophore attempts to perform trans-to-cis isomerization at every double bond together with the stalled anti-to-syn rotations at multiple single bonds along its polyene chain. Such isomerization sampling pushes all seven transmembrane helices to bend outward, resulting in a transient expansion of the retinal binding pocket, and later, a contraction due to recoiling. These ultrafast responses observed at the atomic resolution support that the productive photoreaction in bacteriorhodopsin is initiated by light-induced charge separation in the prosthetic chromophore yet governed by stereoselectivity of its protein pocket. The method of a numerical resolution of concurrent events from mixed observations is also generally applicable.

Best practices for time-resolved serial synchrotron crystallography

With recent developments in X-ray sources, instrumentation and data-analysis tools, time-resolved crystallographic experiments, which were originally the preserve of a few expert groups, are becoming simpler and can be carried out at more radiation sources, and are thus increasingly accessible to a growing user base. However, these experiments are just that: discrete experiments, not just `data collections'. As such, careful planning and consideration of potential pitfalls is required to enable a successful experiment. Here, some of the key factors that should be considered during the planning and execution of a time-resolved structural study are outlined, with a particular focus on synchrotron-based experiments.

Moving beyond static snapshots: Protein dynamics and the Protein Data Bank

Proteins are the molecular machines of living systems. Their dynamics are an intrinsic part of their evolutionary selection in carrying out their biological functions. Although the dynamics are more difficult to observe than a static, average structure, we are beginning to observe these dynamics and form sound mechanistic connections between structure, dynamics, and function. This progress is highlighted in case studies from myoglobin and adenylate kinase to the ribosome and molecular motors where these molecules are being probed with a multitude of techniques across many timescales. New approaches to time-resolved crystallography are allowing simple “movies” to be taken of proteins in action, and new methods of mapping the variations in cryo-electron microscopy are emerging to reveal a more complete description of life’s machines. The results of these new methods are aided in their dissemination by continual improvements in curation and distribution by the Protein Data Bank and their partners around the world.

Dimer Asymmetry and Light Activation Mechanism in Brucella Blue-Light Sensor Histidine Kinase

Bacteria employ two-component systems (TCSs) to sense and respond to changes in their surroundings. At the core of the TCS signaling pathway is the multidomain sensor histidine kinase, where the enzymatic activity of its output domain is allosterically controlled by the input signal perceived by the sensor domain.

The ability to sense and respond to environmental cues is essential for adaptation and survival in living organisms. In bacteria, this process is accomplished by multidomain sensor histidine kinases that undergo autophosphorylation in response to specific stimuli, thereby triggering downstream signaling cascades. However, the molecular mechanism of allosteric activation is not fully understood in these important sensor proteins. Here, we report the full-length crystal structure of a blue light photoreceptor LOV histidine kinase (LOV-HK) involved in light-dependent virulence modulation in the pathogenic bacterium Brucella abortus. Joint analyses of dark and light structures determined in different signaling states have shown that LOV-HK transitions from a symmetric dark structure to a highly asymmetric light state. The initial local and subtle structural signal originated in the chromophore-binding LOV domain alters the dimer asymmetry via a coiled-coil rotary switch and helical bending in the helical spine. These amplified structural changes result in enhanced conformational flexibility and large-scale rearrangements that facilitate the phosphoryl transfer reaction in the HK domain.

Ultrafast Structural Changes Decomposed from Serial Crystallographic Data

Direct visualization of electronic and molecular events during biochemical reactions is essential to mechanistic insights. This Letter presents an in-depth analysis of the serial crystallographic data sets collected by Barends and Schlichting et al. ( Science 2015 , 350 , 445 ) that probe the ligand photodissociation in carbonmonoxy myoglobin. This analysis reveals electron density changes caused by the formation of high-spin 3d atomic orbitals of the heme iron upon photolysis and their dynamic behaviors within the first few picoseconds. The heme iron is found popping out of and recoiling back into the heme plane in succession. These findings provide long-awaited visual validations for previous works using ultrafast spectroscopy and molecular dynamics simulations. Electron density variations are also found largely in the solvent during the first period of a low-frequency oscillation. This work demonstrates the importance of the analytical methods in detecting and isolating weak, transient signals of electronic changes arising from chemical reactions.

Structural basis of molecular logic OR in a dual-sensor histidine kinase

Signal detection and integration by sensory proteins constitute the critical molecular events as living organisms respond to changes in a complex environment. Many sensory proteins adopt a modular architecture that integrates the perception of distinct chemical or physical signals and the generation of a biological response in the same protein molecule. Currently, how signal perception and integration are achieved in such a modular, often dimeric, framework remains elusive. Here, we report a dynamic crystallography study on the tandem sensor domains of a dual-sensor histidine kinase PPHK (phosphorylation-responsive photosensitive histidine kinase) that operates a molecular logic OR, by which the output kinase activity is modulated by a phosphorylation signal and a light signal. A joint analysis of ∼170 crystallographic datasets probing different signaling states shows remarkable dimer asymmetry as PPHK responds to the input signals and transitions from one state to the other. Supported by mutational data and structural analysis, these direct observations reveal the working mechanics of the molecular logic OR in PPHK, where the light-induced bending of a long signaling helix at the dimer interface is counteracted by the ligand-induced structural changes from a different sensor domain. We propose that the logic OR of PPHK, together with an upstream photoreceptor, implements a "long-pass" red light response distinct from those accomplished by classical phytochromes.

Time-resolved crystallography reveals allosteric communication aligned with molecular breathing

A comprehensive understanding of protein function demands correlating structure and dynamic changes. Using time-resolved serial synchrotron crystallography, we visualized half-of-the-sites reactivity and correlated molecular-breathing motions in the enzyme fluoroacetate dehalogenase. Eighteen time points from 30 milliseconds to 30 seconds cover four turnover cycles of the irreversible reaction. They reveal sequential substrate binding, covalent-intermediate formation, setup of a hydrolytic water molecule, and product release. Small structural changes of the protein mold and variations in the number and placement of water molecules accompany the various chemical steps of catalysis. Triggered by enzyme-ligand interactions, these repetitive changes in the protein framework’s dynamics and entropy constitute crucial components of the catalytic machinery.

Data-driven challenges and opportunities in crystallography

Structural biology is in the midst of a revolution fueled by faster and more powerful instruments capable of delivering orders of magnitude more data than their predecessors. This increased pace in data gathering introduces new experimental and computational challenges, frustrating real-time processing and interpretation of data and requiring long-term solutions for data archival and retrieval. This combination of challenges and opportunities is driving the exploration of new areas of structural biology, including studies of macromolecular dynamics and the investigation of molecular ensembles in search of a better understanding of conformational landscapes. The next generation of instruments promises to yield even greater data rates, requiring a concerted effort by institutions, centers and individuals to extract meaning from every bit and make data accessible to the community at large, facilitating data mining efforts by individuals or groups as analysis tools improve.

Structural and dynamical description of the enzymatic reaction of a phosphohexomutase

Enzymes are known to adopt various conformations at different points along their catalytic cycles. Here, we present a comprehensive analysis of 15 isomorphous, high resolution crystal structures of the enzyme phosphoglucomutase from the bacterium Xanthomonas citri. The protein was captured in distinct states critical to function, including enzyme-substrate, enzyme-product, and enzyme-intermediate complexes. Key residues in ligand recognition and regions undergoing conformational change are identified and correlated with the various steps of the catalytic reaction. In addition, we use principal component analysis to examine various subsets of these structures with two goals: (1) identifying sites of conformational heterogeneity through a comparison of room temperature and cryogenic structures of the apo-enzyme and (2) a priori clustering of the enzyme-ligand complexes into functionally related groups, showing sensitivity of this method to structural features difficult to detect by traditional methods. This study captures, in a single system, the structural basis of diverse substrate recognition, the subtle impact of covalent modification, and the role of ligand-induced conformational change in this representative enzyme of the α-D-phosphohexomutase superfamily.

Time-resolved x-ray crystallography capture of a slow reaction tetrahydrofolate intermediate

Time-resolved crystallography is a powerful technique to elucidate molecular mechanisms at both spatial (angstroms) and temporal (picoseconds to seconds) resolutions. We recently discovered an unusually slow reaction at room temperature that occurs on the order of days: the in crystalline reverse oxidative decay of the chemically labile (6S)-5,6,7,8-tetrahydrofolate in complex with its producing enzyme Escherichia coli dihydrofolate reductase. Here, we report the critical analysis of a representative dataset at an intermediate reaction time point. A quinonoid-like intermediate state lying between tetrahydrofolate and dihydrofolate features a near coplanar geometry of the bicyclic pterin moiety, and a tetrahedral sp ³ C6 geometry is proposed based on the apparent mFo-DFc omit electron densities of the ligand. The presence of this intermediate is strongly supported by Bayesian difference refinement. Isomorphous Fo-Fo difference map and multi-state refinement analyses suggest the presence of end-state ligand populations as well, although the putative intermediate state is likely the most populated. A similar quinonoid intermediate previously proposed to transiently exist during the oxidation of tetrahydrofolate was confirmed by polarography and UV-vis spectroscopy to be relatively stable in the oxidation of its close analog tetrahydropterin. We postulate that the constraints on the ligand imposed by the interactions with the protein environment might be the origin of the slow reaction observed by time-resolved crystallography.

Microfocus diffraction from different regions of a protein crystal: structural variations and unit-cell polymorphism

Real macromolecular crystals can be non-ideal in a myriad of ways. This often creates challenges for structure determination, while also offering opportunities for greater insight into the crystalline state and the dynamic behavior of macromolecules. To evaluate whether different parts of a single crystal of a dynamic protein, EutL, might be informative about crystal and protein polymorphism, a microfocus X-ray synchrotron beam was used to collect a series of 18 separate data sets from non-overlapping regions of the same crystal specimen. A principal component analysis (PCA) approach was employed to compare the structure factors and unit cells across the data sets, and it was found that the 18 data sets separated into two distinct groups, with large R values (in the 40% range) and significant unit-cell variations between the members of the two groups. This categorization mapped the different data-set types to distinct regions of the crystal specimen. Atomic models of EutL were then refined against two different data sets obtained by separately merging data from the two distinct groups. A comparison of the two resulting models revealed minor but discernable differences in certain segments of the protein structure, and regions of higher deviation were found to correlate with regions where larger dynamic motions were predicted to occur by normal-mode molecular-dynamics simulations. The findings emphasize that large spatially dependent variations may be present across individual macromolecular crystals. This information can be uncovered by simultaneous analysis of multiple partial data sets and can be exploited to reveal new insights about protein dynamics, while also improving the accuracy of the structure-factor data ultimately obtained in X-ray diffraction experiments.

Real-space analysis of radiation-induced specific changes with independent component analysis

A method of analysis is presented that allows for the separation of specific radiation-induced changes into distinct components in real space. The method relies on independent component analysis (ICA) and can be effectively applied to electron density maps and other types of maps, provided that they can be represented as sets of numbers on a grid. Here, for glucose isomerase crystals, ICA was used in a proof-of-concept analysis to separate temperature-dependent and temperature-independent components of specific radiation-induced changes for data sets acquired from multiple crystals across multiple temperatures. ICA identified two components, with the temperature-independent component being responsible for the majority of specific radiation-induced changes at temperatures below 130 K. The patterns of specific temperature-independent radiation-induced changes suggest a contribution from the tunnelling of electron holes as a possible explanation. In the second case, where a group of 22 data sets was collected on a single thaumatin crystal, ICA was used in another type of analysis to separate specific radiation-induced effects happening on different exposure-level scales. Here, ICA identified two components of specific radiation-induced changes that likely result from radiation-induced chemical reactions progressing with different rates at different locations in the structure. In addition, ICA unexpectedly identified the radiation-damage state corresponding to reduced disulfide bridges rather than the zero-dose extrapolated state as the highest contrast structure. The application of ICA to the analysis of specific radiation-induced changes in real space and the data pre-processing for ICA that relies on singular value decomposition, which was used previously in data space to validate a two-component physical model of X-ray radiation-induced changes, are discussed in detail. This work lays a foundation for a better understanding of protein-specific radiation chemistries and provides a framework for analysing effects of specific radiation damage in crystallographic and cryo-EM experiments.

Photoactivation mechanism of a carotenoid-based photoreceptor

Photoprotection is essential for efficient photosynthesis. Cyanobacteria have evolved a unique photoprotective mechanism mediated by a water-soluble carotenoid-based photoreceptor known as orange carotenoid protein (OCP). OCP undergoes large conformational changes in response to intense blue light, and the photoactivated OCP facilitates dissipation of excess energy via direct interaction with allophycocyanins at the phycobilisome core. However, the structural events leading up to the OCP photoactivation remain elusive at the molecular level. Here we present direct observations of light-induced structural changes in OCP captured by dynamic crystallography. Difference electron densities between the dark and illuminated states reveal widespread and concerted atomic motions that lead to altered protein-pigment interactions, displacement of secondary structures, and domain separation. Based on these crystallographic observations together with site-directed mutagenesis, we propose a molecular mechanism for OCP light perception, in which the photochemical property of a conjugated carbonyl group is exploited. We hypothesize that the OCP photoactivation starts with keto-enol tautomerization of the essential 4-keto group in the carotenoid, which disrupts the strong hydrogen bonds between the bent chromophore and the protein moiety. Subsequent structural changes trapped in the crystal lattice offer a high-resolution glimpse of the initial molecular events as OCP begins to transition from the orange-absorbing state to the active red-absorbing state.

Graphene-based microfluidics for serial crystallography

Microfluidic strategies to enable the growth and subsequent serial crystallographic analysis of micro-crystals have the potential to facilitate both structural characterization and dynamic structural studies of protein targets that have been resistant to single-crystal strategies. However, adapting microfluidic crystallization platforms for micro-crystallography requires a dramatic decrease in the overall device thickness. We report a robust strategy for the straightforward incorporation of single-layer graphene into ultra-thin microfluidic devices. This architecture allows for a total material thickness of only ∼1 μm, facilitating on-chip X-ray diffraction analysis while creating a sample environment that is stable against significant water loss over several weeks. We demonstrate excellent signal-to-noise in our X-ray diffraction measurements using a 1.5 μs polychromatic X-ray exposure, and validate our approach via on-chip structure determination using hen egg white lysozyme (HEWL) as a model system. Although this work is focused on the use of graphene for protein crystallography, we anticipate that this technology should find utility in a wide range of both X-ray and other lab on a chip applications.

Angular-split/temporal-delay approach to ultrafast protein dynamics at XFELs

X-ray crystallography promises direct insights into electron-density changes that lead to and arise from structural changes such as electron and proton transfer and the formation, rupture and isomerization of chemical bonds. The ultrashort pulses of hard X-rays produced by free-electron lasers present an exciting opportunity for capturing ultrafast structural events in biological macromolecules within femtoseconds after photoexcitation. However, shot-to-shot fluctuations, which are inherent to the very process of self-amplified spontaneous emission (SASE) that generates the ultrashort X-ray pulses, are a major source of noise that may conceal signals from structural changes. Here, a new approach is proposed to angularly split a single SASE pulse and to produce a temporal delay of picoseconds between the split pulses. These split pulses will allow the probing of two distinct states before and after photoexcitation triggered by a laser pulse between the split X-ray pulses. The split pulses originate from a single SASE pulse and share many common properties; thus, noise arising from shot-to-shot fluctuations is self-canceling. The unambiguous interpretation of ultrafast structural changes would require diffraction data at atomic resolution, as these changes may or may not involve any atomic displacement. This approach, in combination with the strategy of serial crystallography, offers a solution to study ultrafast dynamics of light-initiated biochemical reactions or biological processes at atomic resolution.

Conformational heterogeneity of the Pfr chromophore in plant and cyanobacterial phytochromes

Phytochromes are biological photoreceptors that can be reversibly photoconverted between a dark and photoactivated state. The underlying reaction sequences are initiated by the photoisomerization of the tetrapyrrole cofactor, which in plant and cyanobacterial phytochromes are a phytochromobilin (PΦB) and a phycocyanobilin (PCB), respectively. The transition between the two states represents an on/off-switch of the output module activating or deactivating downstream physiological processes. In addition, the photoactivated state, i.e., Pfr in canonical phytochromes, can be thermally reverted to the dark state (Pr). The present study aimed to improve our understanding of the specific reactivity of various PΦB- and PCB-binding phytochromes in the Pfr state by analysing the cofactor structure by vibrational spectroscopic techniques. Resonance Raman (RR) spectroscopy revealed two Pfr conformers (Pfr-I and Pfr-II) forming a temperature-dependent conformational equilibrium. The two sub-states—found in all phytochromes studied, albeit with different relative contributions—differ in structural details of the C-D and A-B methine bridges. In the Pfr-I sub-state the torsion between the rings C and D is larger by ca. 10° compared to Pfr-II. This structural difference is presumably related to different hydrogen bonding interactions of ring D as revealed by time-resolved IR spectroscopic studies of the cyanobacterial phytochrome Cph1. The transitions between the two sub-states are evidently too fast (i.e., nanosecond time scale) to be resolved by NMR spectroscopy which could not detect a structural heterogeneity of the chromophore in Pfr. The implications of the present findings for the dark reversion of the Pfr state are discussed.

Towards time-resolved serial crystallography in a microfluidic device

Serial methods for crystallography have the potential to enable dynamic structural studies of protein targets that have been resistant to single-crystal strategies. The use of serial data-collection strategies can circumvent challenges associated with radiation damage and repeated reaction initiation. This work utilizes a microfluidic crystallization platform for the serial time-resolved Laue diffraction analysis of macroscopic crystals of photoactive yellow protein (PYP). Reaction initiation was achieved via pulsed laser illumination, and the resultant electron-density difference maps clearly depict the expected pR(1)/pR(E46Q) and pR(2)/pR(CW) states at 10 µs and the pB1 intermediate at 1 ms. The strategies presented here have tremendous potential for extension to chemical triggering methods for reaction initiation and for extension to dynamic, multivariable analyses.

Molecular dynamics simulation of triclinic lysozyme in a crystal lattice

Molecular dynamics simulations of crystals can enlighten interpretation of experimental X-ray crystallography data and elucidate structural dynamics and heterogeneity in biomolecular crystals. Furthermore, because of the direct comparison against experimental data, they can inform assessment of molecular dynamics methods and force fields. We present microsecond scale results for triclinic hen egg-white lysozyme in a supercell consisting of 12 independent unit cells using four contemporary force fields (Amber ff99SB, ff14ipq, ff14SB, and CHARMM 36) in crystalline and solvated states (for ff14SB only). We find the crystal simulations consistent across multiple runs of the same force field and robust to various solvent equilibration schemes. However, convergence is slow compared with solvent simulations. All the tested force fields reproduce experimental structural and dynamic properties well, but Amber ff14SB maintains structure and reproduces fluctuations closest to the experimental model: its average backbone structure differs from the deposited structure by 0.37Å; by contrast, the average backbone structure in solution differs from the deposited by 0.65Å. All the simulations are affected by a small progressive deterioration of the crystal lattice, presumably due to imperfect modeling of hydrogen bonding and other crystal contact interactions; this artifact is smallest in ff14SB, with average lattice positions deviating by 0.20Å from ideal. Side-chain disorder is surprisingly low with fewer than 30% of the nonglycine or alanine residues exhibiting significantly populated alternate rotamers. Our results provide helpful insight into the methodology of biomolecular crystal simulations and indicate directions for future work to obtain more accurate energy models for molecular dynamics.

Chemical selectivity in structure determination by the time dependent analysis of in situ XRPD data: a clear view of Xe thermal behavior inside a MFI zeolite

PSD/PCA analysis of MED data allowed to enhance the chemical selectivity in X-ray powder diffraction and to obtain Xe substructure into MFI zeolite.

Dynamic Crystallography Reveals Early Signalling Events in Ultraviolet Photoreceptor UVR8

Arabidopsis thaliana UVR8 (AtUVR8) is a long-sought-after photoreceptor that undergoes dimer dissociation in response to UV-B light. Crystallographic and mutational studies have identified two crucial tryptophan residues for UV-B responses in AtUVR8. However, the mechanism of UV-B perception and structural events leading up to dimer dissociation remain elusive at the molecular level. We applied dynamic crystallography to capture light-induced structural events in photoactive AtUVR8 crystals. Here we report two intermediate structures at 1.67Å resolution. At the epicenter of UV-B signaling, concerted motions associated with Trp285/Trp233 lead to ejection of a water molecule, which weakens an intricate network of hydrogen bonds and salt bridges at the dimer interface. Partial opening of the β-propeller structure due to thermal relaxation of conformational strains originating in the epicenter further disrupts the dimer interface and leads to dimer dissociation. These dynamic crystallographic observations provide structural insights into the photo-perception and signaling mechanism of UVR8.

In situ serial Laue diffraction on a microfluidic crystallization device

Renewed interest in room-temperature diffraction has been prompted by the desire to observe structural dynamics of proteins as they function. Serial crystallography, an experimental strategy that aggregates small pieces of data from a large uniform pool of crystals, has been demonstrated at synchrotrons and X-ray free-electron lasers. This work utilizes a microfluidic crystallization platform for serial Laue diffraction from macroscopic crystals and proposes that a collection of small slices of Laue data from many individual crystals is a realistic solution to the difficulties in dynamic studies of irreversible biochemical reactions.

Singular value decomposition as a tool for background corrections in time-resolved XFEL scattering data

The development of new X-ray light sources, XFELs, with unprecedented time and brilliance characteristics has led to the availability of very large datasets with high time resolution and superior signal strength. The chaotic nature of the emission processes in such sources as well as entirely novel detector demands has also led to significant challenges in terms of data analysis. This paper describes a heuristic approach to datasets where spurious background contributions of a magnitude similar to (or larger) than the signal of interest prevents conventional analysis approaches. The method relies on singular-value decomposition of no-signal subsets of acquired datasets in combination with model inputs and appears generally applicable to time-resolved X-ray diffuse scattering experiments.

Reaction trajectory revealed by a joint analysis of protein data bank

Structural motions along a reaction pathway hold the secret about how a biological macromolecule functions. If each static structure were considered as a snapshot of the protein molecule in action, a large collection of structures would constitute a multidimensional conformational space of an enormous size. Here I present a joint analysis of hundreds of known structures of human hemoglobin in the Protein Data Bank. By applying singular value decomposition to distance matrices of these structures, I demonstrate that this large collection of structural snapshots, derived under a wide range of experimental conditions, arrange orderly along a reaction pathway. The structural motions along this extensive trajectory, including several helical transformations, arrive at a reverse engineered mechanism of the cooperative machinery (Ren, companion article), and shed light on pathological properties of the abnormal homotetrameric hemoglobins from α-thalassemia. This method of meta-analysis provides a general approach to structural dynamics based on static protein structures in this post genomics era.