Influence of pump laser fluence on ultrafast myoglobin structural dynamics.

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump-probe experiments for the investigation of electronic and nuclear changes during light-induced reactions. On timescales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer1,2. However, all ultrafast TR-SFX studies to date have employed such high pump laser energies that nominally several photons were absorbed per chromophore3-17. As multiphoton absorption may force the protein response into non-physiological pathways, it is of great concern18,19 whether this experimental approach20 allows valid conclusions to be drawn vis-à-vis biologically relevant single-photon-induced reactions18,19. Here we describe ultrafast pump-probe SFX experiments on the photodissociation of carboxymyoglobin, showing that different pump laser fluences yield markedly different results. In particular, the dynamics of structural changes and observed indicators of the mechanistically important coherent oscillations of the Fe-CO bond distance (predicted by recent quantum wavepacket dynamics21) are seen to depend strongly on pump laser energy, in line with quantum chemical analysis. Our results confirm both the feasibility and necessity of performing ultrafast TR-SFX pump-probe experiments in the linear photoexcitation regime. We consider this to be a starting point for reassessing both the design and the interpretation of ultrafast TR-SFX pump-probe experiments20 such that mechanistically relevant insight emerges.

Structural evidence for intermediates during O2 formation in photosystem II.

In natural photosynthesis, the light-driven splitting of water into electrons, protons and molecular oxygen forms the first step of the solar-to-chemical energy conversion process. The reaction takes place in photosystem II, where the Mn4CaO5 cluster first stores four oxidizing equivalents, the S0 to S4 intermediate states in the Kok cycle, sequentially generated by photochemical charge separations in the reaction center and then catalyzes the O-O bond formation chemistry1-3. Here, we report room temperature snapshots by serial femtosecond X-ray crystallography to provide structural insights into the final reaction step of Kok's photosynthetic water oxidation cycle, the S3→[S4]→S0 transition where O2 is formed and Kok's water oxidation clock is reset. Our data reveal a complex sequence of events, which occur over micro- to milliseconds, comprising changes at the Mn4CaO5 cluster, its ligands and water pathways as well as controlled proton release through the hydrogen-bonding network of the Cl1 channel. Importantly, the extra O atom Ox, which was introduced as a bridging ligand between Ca and Mn1 during the S2→S3 transition4-6, disappears or relocates in parallel with Yz reduction starting at approximately 700 µs after the third flash. The onset of O2 evolution, as indicated by the shortening of the Mn1-Mn4 distance, occurs at around 1,200 µs, signifying the presence of a reduced intermediate, possibly a bound peroxide.

Structures of the intermediates of Kok's photosynthetic water oxidation clock.

Inspired by the period-four oscillation in flash-induced oxygen evolution of photosystem II discovered by Joliot in 1969, Kok performed additional experiments and proposed a five-state kinetic model for photosynthetic oxygen evolution, known as Kok's S-state clock or cycle1,2. The model comprises four (meta)stable intermediates (S0, S1, S2 and S3) and one transient S4 state, which precedes dioxygen formation occurring in a concerted reaction from two water-derived oxygens bound at an oxo-bridged tetra manganese calcium (Mn4CaO5) cluster in the oxygen-evolving complex3-7. This reaction is coupled to the two-step reduction and protonation of the mobile plastoquinone QB at the acceptor side of PSII. Here, using serial femtosecond X-ray crystallography and simultaneous X-ray emission spectroscopy with multi-flash visible laser excitation at room temperature, we visualize all (meta)stable states of Kok's cycle as high-resolution structures (2.04-2.08 Å). In addition, we report structures of two transient states at 150 and 400 µs, revealing notable structural changes including the binding of one additional 'water', Ox, during the S2→S3 state transition. Our results suggest that one water ligand to calcium (W3) is directly involved in substrate delivery. The binding of the additional oxygen Ox in the S3 state between Ca and Mn1 supports O-O bond formation mechanisms involving O5 as one substrate, where Ox is either the other substrate oxygen or is perfectly positioned to refill the O5 position during O2 release. Thus, our results exclude peroxo-bond formation in the S3 state, and the nucleophilic attack of W3 onto W2 is unlikely.

Radiation damage and dose limits in serial synchrotron crystallography at cryo- and room temperatures.

Radiation damage limits the accuracy of macromolecular structures in X-ray crystallography. Cryogenic (cryo-) cooling reduces the global radiation damage rate and, therefore, became the method of choice over the past decades. The recent advent of serial crystallography, which spreads the absorbed energy over many crystals, thereby reducing damage, has rendered room temperature (RT) data collection more practical and also extendable to microcrystals, both enabling and requiring the study of specific and global radiation damage at RT. Here, we performed sequential serial raster-scanning crystallography using a microfocused synchrotron beam that allowed for the collection of two series of 40 and 90 full datasets at 2- and 1.9-Å resolution at a dose rate of 40.3 MGy/s on hen egg white lysozyme (HEWL) crystals at RT and cryotemperature, respectively. The diffraction intensity halved its initial value at average doses (D1/2) of 0.57 and 15.3 MGy at RT and 100 K, respectively. Specific radiation damage at RT was observed at disulfide bonds but not at acidic residues, increasing and then apparently reversing, a peculiar behavior that can be modeled by accounting for differential diffraction intensity decay due to the nonuniform illumination by the X-ray beam. Specific damage to disulfide bonds is evident early on at RT and proceeds at a fivefold higher rate than global damage. The decay modeling suggests it is advisable not to exceed a dose of 0.38 MGy per dataset in static and time-resolved synchrotron crystallography experiments at RT. This rough yardstick might change for proteins other than HEWL and at resolutions other than 2 Å.

Untangling the sequence of events during the S2 → S3 transition in photosystem II and implications for the water oxidation mechanism.

In oxygenic photosynthesis, light-driven oxidation of water to molecular oxygen is carried out by the oxygen-evolving complex (OEC) in photosystem II (PS II). Recently, we reported the room-temperature structures of PS II in the four (semi)stable S-states, S1, S2, S3, and S0, showing that a water molecule is inserted during the S2 → S3 transition, as a new bridging O(H)-ligand between Mn1 and Ca. To understand the sequence of events leading to the formation of this last stable intermediate state before O2 formation, we recorded diffraction and Mn X-ray emission spectroscopy (XES) data at several time points during the S2 → S3 transition. At the electron acceptor site, changes due to the two-electron redox chemistry at the quinones, QA and QB, are observed. At the donor site, tyrosine YZ and His190 H-bonded to it move by 50 µs after the second flash, and Glu189 moves away from Ca. This is followed by Mn1 and Mn4 moving apart, and the insertion of OX(H) at the open coordination site of Mn1. This water, possibly a ligand of Ca, could be supplied via a "water wheel"-like arrangement of five waters next to the OEC that is connected by a large channel to the bulk solvent. XES spectra show that Mn oxidation (τ of ∼350 µs) during the S2 → S3 transition mirrors the appearance of OX electron density. This indicates that the oxidation state change and the insertion of water as a bridging atom between Mn1 and Ca are highly correlated.

Hydroxyl radical mediated damage of proteins in low oxygen solution investigated using X-ray footprinting mass spectrometry.

In the method of X-ray footprinting mass spectrometry (XFMS), proteins at micromolar concentration in solution are irradiated with a broadband X-ray source, and the resulting hydroxyl radical modifications are characterized using liquid chromatography mass spectrometry to determine sites of solvent accessibility. These data are used to infer structural changes in proteins upon interaction with other proteins, folding, or ligand binding. XFMS is typically performed under aerobic conditions; dissolved molecular oxygen in solution is necessary in many, if not all, the hydroxyl radical modifications that are generally reported. In this study we investigated the result of X-ray induced modifications to three different proteins under aerobic versus low oxygen conditions, and correlated the extent of damage with dose calculations. We observed a concentration-dependent protecting effect at higher protein concentration for a given X-ray dose. For the typical doses used in XFMS experiments there was minimal X-ray induced aggregation and fragmentation, but for higher doses we observed formation of covalent higher molecular weight oligomers, as well as fragmentation, which was affected by the amount of dissolved oxygen in solution. The higher molecular weight products in the form of dimers, trimers, and tetramers were present in all sample preparations, and, upon X-ray irradiation, these oligomers became non-reducible as seen in SDS-PAGE. The results provide an important contribution to the large body of X-ray radiation damage literature in structural biology research, and will specifically help inform the future planning of XFMS, and well as X-ray crystallography and small-angle X-ray scattering experiments.

Macromolecular X-ray structure determination using weak, single-wavelength anomalous data.

We describe a likelihood-based method for determining the substructure of anomalously scattering atoms in macromolecular crystals that allows successful structure determination by single-wavelength anomalous diffraction (SAD) X-ray analysis with weak anomalous signal. With the use of partial models and electron density maps in searches for anomalously scattering atoms, testing of alternative values of parameters and parallelized automated model-building, this method has the potential to extend the applicability of the SAD method in challenging cases.

Protein structural ensembles are revealed by redefining X-ray electron density noise.

To increase the power of X-ray crystallography to determine not only the structures but also the motions of biomolecules, we developed methods to address two classic crystallographic problems: putting electron density maps on the absolute scale of e(-)/Å(3) and calculating the noise at every point in the map. We find that noise varies with position and is often six to eight times lower than thresholds currently used in model building. Analyzing the rescaled electron density maps from 485 representative proteins revealed unmodeled conformations above the estimated noise for 45% of side chains and a previously hidden, low-occupancy inhibitor of HIV capsid protein. Comparing the electron density maps in the free and nucleotide-bound structures of three human protein kinases suggested that substrate binding perturbs distinct intrinsic allosteric networks that link the active site to surfaces that recognize regulatory proteins. These results illustrate general approaches to identify and analyze alternative conformations, low-occupancy small molecules, solvent distributions, communication pathways, and protein motions.

Mycobacterium tuberculosis FtsX extracellular domain activates the peptidoglycan hydrolase, RipC.

Bacterial growth and cell division are coordinated with hydrolysis of the peptidoglycan (PG) layer of the cell wall, but the mechanisms of regulation of extracellular PG hydrolases are not well understood. Here we report the biochemical, structural, and genetic analysis of the Mycobacterium tuberculosis homolog of the transmembrane PG-hydrolase regulator, FtsX. The purified FtsX extracellular domain binds the PG peptidase Rv2190c/RipC N-terminal segment, causing a conformational change that activates the enzyme. Deletion of ftsEX and ripC caused similar phenotypes in Mycobacterium smegmatis, as expected for genes in a single pathway. The crystal structure of the FtsX extracellular domain reveals an unprecedented fold containing two lobes connected by a flexible hinge. Mutations in the hydrophobic cleft between the lobes reduce RipC binding in vitro and inhibit FtsX function in M. smegmatis. These studies suggest how FtsX recognizes RipC and support a model in which a conformational change in FtsX links the cell division apparatus with PG hydrolysis.

Goniometer-based femtosecond crystallography with X-ray free electron lasers.

The emerging method of femtosecond crystallography (FX) may extend the diffraction resolution accessible from small radiation-sensitive crystals and provides a means to determine catalytically accurate structures of acutely radiation-sensitive metalloenzymes. Automated goniometer-based instrumentation developed for use at the Linac Coherent Light Source enabled efficient and flexible FX experiments to be performed on a variety of sample types. In the case of rod-shaped Cpl hydrogenase crystals, only five crystals and about 30 min of beam time were used to obtain the 125 still diffraction patterns used to produce a 1.6-Å resolution electron density map. For smaller crystals, high-density grids were used to increase sample throughput; 930 myoglobin crystals mounted at random orientation inside 32 grids were exposed, demonstrating the utility of this approach. Screening results from cryocooled crystals of ß2-adrenoreceptor and an RNA polymerase II complex indicate the potential to extend the diffraction resolution obtainable from very radiation-sensitive samples beyond that possible with undulator-based synchrotron sources.

Predicting the X-ray lifetime of protein crystals.

Radiation damage is a major cause of failure in macromolecular crystallography experiments. Although it is always best to evenly illuminate the entire volume of a homogeneously diffracting crystal, limitations of the available equipment and imperfections in the sample often require a more sophisticated targeting strategy, involving microbeams smaller than the crystal, and translations of the crystal during data collection. This leads to a highly inhomogeneous distribution of absorbed X-rays (i.e., dose). Under these common experimental conditions, the relationship between dose and time is nonlinear, making it difficult to design an experimental strategy that optimizes the radiation damage lifetime of the crystal, or to assign appropriate dose values to an experiment. We present, and experimentally validate, a predictive metric diffraction-weighted dose for modeling the rate of decay of total diffracted intensity from protein crystals in macromolecular crystallography, and hence we can now assign appropriate "dose" values to modern experimental setups. Further, by taking the ratio of total elastic scattering to diffraction-weighted dose, we show that it is possible to directly compare potential data-collection strategies to optimize the diffraction for a given level of damage under specific experimental conditions. As an example of the applicability of this method, we demonstrate that by offsetting the rotation axis from the beam axis by 1.25 times the full-width half maximum of the beam, it is possible to significantly extend the dose lifetime of the crystal, leading to a higher number of diffracted photons, better statistics, and lower overall radiation damage.

Hyaloperonospora arabidopsidis ATR1 effector is a repeat protein with distributed recognition surfaces.

The in planta association of the Hyaloperonospora arabidopsidis effector ATR1 with the cognate Arabidopsis thaliana RPP1 immune receptor activates a disease-resistance signaling pathway that inhibits pathogen growth. To define the molecular events specifying effector recognition by RPP1, we determined the crystal structure of ATR1 and assayed in planta the effects of surface polymorphisms that are critical to activating plant immunity. ATR1 adopts an elongated, all-helical, two-domain, seahorse-like structure with an overall architecture unlike any previously described fold. Structural comparisons highlight a tandemly duplicated, five-helix motif in the C-terminal domain that creates a structural framework for rapid diversification. Identification and mapping of critical recognition sites suggest that ATR1 detection by the RPP1 resistance protein is mediated by several distinct protein surfaces that allow the effectors to escape recognition through diverse surface polymorphisms. ATR1 gain-of-recognition mutants demonstrate that multiple amino acid substitutions are necessary for recognition and that surface polymorphisms exert additive effects. These results suggest that ATR1 is a modular repeat protein belonging to an ancient family of oomycete effectors that rapidly evolves to escape host detection and adopt diverse virulence functions.

Accessing protein conformational ensembles using room-temperature X-ray crystallography.

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.

Deep residual networks for crystallography trained on synthetic data.

The use of artificial intelligence to process diffraction images is challenged by the need to assemble large and precisely designed training data sets. To address this, a codebase called Resonet was developed for synthesizing diffraction data and training residual neural networks on these data. Here, two per-pattern capabilities of Resonet are demonstrated: (i) interpretation of crystal resolution and (ii) identification of overlapping lattices. Resonet was tested across a compilation of diffraction images from synchrotron experiments and X-ray free-electron laser experiments. Crucially, these models readily execute on graphics processing units and can thus significantly outperform conventional algorithms. While Resonet is currently utilized to provide real-time feedback for macromolecular crystallography users at the Stanford Synchrotron Radiation Lightsource, its simple Python-based interface makes it easy to embed in other processing frameworks. This work highlights the utility of physics-based simulation for training deep neural networks and lays the groundwork for the development of additional models to enhance diffraction collection and analysis.

Crystallographic data processing for free-electron laser sources.

A processing pipeline for diffraction data acquired using the `serial crystallography' methodology with a free-electron laser source is described with reference to the crystallographic analysis suite CrystFEL and the pre-processing program Cheetah. A detailed analysis of the nature and impact of indexing ambiguities is presented. Simulations of the Monte Carlo integration scheme, which accounts for the partially recorded nature of the diffraction intensities, are presented and show that the integration of partial reflections could be made to converge more quickly if the bandwidth of the X-rays were to be increased by a small amount or if a slight convergence angle were introduced into the incident beam.

Molecular functions of the histone acetyltransferase chaperone complex Rtt109-Vps75.

Histone acetylation and nucleosome remodeling regulate DNA damage repair, replication and transcription. Rtt109, a recently discovered histone acetyltransferase (HAT) from Saccharomyces cerevisiae, functions with the histone chaperone Asf1 to acetylate lysine K56 on histone H3 (H3K56), a modification associated with newly synthesized histones. In vitro analysis of Rtt109 revealed that Vps75, a Nap1 family histone chaperone, could also stimulate Rtt109-dependent acetylation of H3K56. However, the molecular function of the Rtt109-Vps75 complex remains elusive. Here we have probed the molecular functions of Vps75 and the Rtt109-Vps75 complex through biochemical, structural and genetic means. We find that Vps75 stimulates the kcat of histone acetylation by approximately 100-fold relative to Rtt109 alone and enhances acetylation of K9 in the H3 histone tail. Consistent with the in vitro evidence, cells lacking Vps75 showed a substantial reduction (60%) in H3K9 acetylation during S phase. X-ray structural, biochemical and genetic analyses of Vps75 indicate a unique, structurally dynamic Nap1-like fold that suggests a potential mechanism of Vps75-dependent activation of Rttl09. Together, these data provide evidence for a multifunctional HAT-chaperone complex that acetylates histone H3 and deposits H3-H4 onto DNA, linking histone modification and nucleosome assembly.

Basis for a ubiquitin-like protein thioester switch toggling E1-E2 affinity.

Ubiquitin-like proteins (UBLs) are conjugated by dynamic E1-E2-E3 enzyme cascades. E1 enzymes activate UBLs by catalysing UBL carboxy-terminal adenylation, forming a covalent E1 throught UBL thioester intermediate, and generating a thioester-linked E2 throught UBL product, which must be released for subsequent reactions. Here we report the structural analysis of a trapped UBL activation complex for the human NEDD8 pathway, containing NEDD8's heterodimeric E1 (APPBP1-UBA3), two NEDD8s (one thioester-linked to E1, one noncovalently associated for adenylation), a catalytically inactive E2 (Ubc12), and MgATP. The results suggest that a thioester switch toggles E1-E2 affinities. Two E2 binding sites depend on NEDD8 being thioester-linked to E1. One is unmasked by a striking E1 conformational change. The other comes directly from the thioester-bound NEDD8. After NEDD8 transfer to E2, reversion to an alternate E1 conformation would facilitate release of the E2 throught NEDD8 thioester product. Thus, transferring the UBL's thioester linkage between successive conjugation enzymes can induce conformational changes and alter interaction networks to drive consecutive steps in UBL cascades.

Large language models generate functional protein sequences across diverse families.

Deep-learning language models have shown promise in various biotechnological applications, including protein design and engineering. Here we describe ProGen, a language model that can generate protein sequences with a predictable function across large protein families, akin to generating grammatically and semantically correct natural language sentences on diverse topics. The model was trained on 280 million protein sequences from >19,000 families and is augmented with control tags specifying protein properties. ProGen can be further fine-tuned to curated sequences and tags to improve controllable generation performance of proteins from families with sufficient homologous samples. Artificial proteins fine-tuned to five distinct lysozyme families showed similar catalytic efficiencies as natural lysozymes, with sequence identity to natural proteins as low as 31.4%. ProGen is readily adapted to diverse protein families, as we demonstrate with chorismate mutase and malate dehydrogenase.