Direct observation of cooperative protein structural dynamics of homodimeric hemoglobin from 100 ps to 10 ms with pump-probe X-ray solution scattering

Glycerol-slaved 1H-1H NMR cross-relaxation in quasi-native lysozyme

1H-1H nuclear cross-relaxation experiments have been carried out with lysozyme in variable glycerol viscosity to study intramolecular motion, self-diffusion, and isotropic rigid-body rotational tumbling at 298 K, pH 3.8. Dynamics of intramolecular 1H-1H cross-relaxation rates, the increase in internuclear spatial distances, and lateral and rotational diffusion coefficients all show fractional viscosity dependence with a power law exponent κ in the 0.17-0.83 range. The diffusion coefficient of glycerol Ds with the bulk viscosity itself is non-Stokesian, having a fractional viscosity dependence on the medium viscosity (Ds ∼ η-κ, κ ≈ 0.71). The concurrence and close similarity of the fractional viscosity dependence of glycerol diffusion on the one hand, and diffusion and intramolecular cross-relaxation rates of the protein on the other lead to infer that relaxation of glycerol slaves protein relaxations. Glycerol-transformed native lysozyme to a quasi-native state does not affect the conclusion that both global and internal fluctuations are slaved to glycerol relaxation.

Structural dynamics of protein-protein association involved in the light-induced transition of Avena sativa LOV2 protein

The Light-oxygen-voltage-sensing domain (LOV) superfamily, found in enzymes and signal transduction proteins, plays a crucial role in converting light signals into structural signals, mediating various biological mechanisms. While time-resolved spectroscopic studies have revealed the dynamics of the LOV-domain chromophore's electronic structures, understanding the structural changes in the protein moiety, particularly regarding light-induced dimerization, remains challenging. Here, we utilize time-resolved X-ray liquidography to capture the light-induced dimerization of Avena sativa LOV2. Our analysis unveils that dimerization occurs within milliseconds after the unfolding of the A'α and Jα helices in the microsecond time range. Notably, our findings suggest that protein-protein interactions (PPIs) among the β-scaffolds, mediated by helix unfolding, play a key role in dimerization. In this work, we offer structural insights into the dimerization of LOV2 proteins following structural changes in the A'α and Jα helices, as well as mechanistic insights into the protein-protein association process driven by PPIs.

Cooperative Protein Dynamics of Heterotetrameric Hemoglobin from Scapharca inaequivalvis

Hemoglobins achieve cooperative oxygen binding by diverse strategies based on different assemblies of globin subunits. Heterotetrameric hemoglobin from Scapharca inaequivalvis (HbII) consists of two AB-dimers, whose structure closely resembles that of homodimeric hemoglobin from the same organism (HbI). Herein, we investigated the structural dynamics of HbII following carbon monoxide (CO) dissociation using time-resolved resonance Raman (RR) spectroscopy. The observed spectra showed that the heme structure of the transient dissociated form of HbII was similar to that of HbI; however, the transition from the transient dissociated form to the equilibrium unligated form was faster for HbII than for HbI. Furthermore, the dependence of the time-resolved spectra on the yield of CO dissociation revealed that the transition became faster as the number of dissociated ligands increased from one to four. The positive correlation between the rate constants and number of dissociated ligands indicates that the structural transition of HbII following CO dissociation is cooperative.

BioCARS: Synchrotron facility for probing structural dynamics of biological macromolecules

A major goal in biomedical science is to move beyond static images of proteins and other biological macromolecules to the internal dynamics underlying their function. This level of study is necessary to understand how these molecules work and to engineer new functions and modulators of function. Stemming from a visionary commitment to this problem by Keith Moffat decades ago, a community of structural biologists has now enabled a set of x-ray scattering technologies for observing intramolecular dynamics in biological macromolecules at atomic resolution and over the broad range of timescales over which motions are functionally relevant. Many of these techniques are provided by BioCARS, a cutting-edge synchrotron radiation facility built under Moffat leadership and located at the Advanced Photon Source at Argonne National Laboratory. BioCARS enables experimental studies of molecular dynamics with time resolutions spanning from 100 ps to seconds and provides both time-resolved x-ray crystallography and small- and wide-angle x-ray scattering. Structural changes can be initiated by several methods-UV/Vis pumping with tunable picosecond and nanosecond laser pulses, substrate diffusion, and global perturbations, such as electric field and temperature jumps. Studies of dynamics typically involve subtle perturbations to molecular structures, requiring specialized computational techniques for data processing and interpretation. In this review, we present the challenges in experimental macromolecular dynamics and describe the current state of experimental capabilities at this facility. As Moffat imagined years ago, BioCARS is now positioned to catalyze the scientific community to make fundamental advances in understanding proteins and other complex biological macromolecules.

Single-Molecule X-ray Scattering Used to Visualize the Conformation Distribution of Biological Molecules via Single-Object Scattering Sampling

Biological macromolecules, the fundamental building blocks of life, exhibit dynamic structures in their natural environment. Traditional structure determination techniques often oversimplify these multifarious conformational spectra by capturing only ensemble- and time-averaged molecular structures. Addressing this gap, in this work, we extend the application of the single-object scattering sampling (SOSS) method to diverse biological molecules, including RNAs and proteins. Our approach, referred to as "Bio-SOSS", leverages ultrashort X-ray pulses to capture instantaneous structures. In Bio-SOSS, we employ two gold nanoparticles (AuNPs) as labels, which provide strong contrast in the X-ray scattering signal, to ensure precise distance determinations between labeled sites. We generated hypothetical Bio-SOSS images for RNAs, proteins, and an RNA-protein complex, each labeled with two AuNPs at specified positions. Subsequently, to validate the accuracy of Bio-SOSS, we extracted distances between these nanoparticle labels from the images and compared them with the actual values used to generate the Bio-SOSS images. Specifically, for a representative RNA (1KXK), the standard deviation in distance discrepancies between molecular dynamics snapshots and Bio-SOSS retrievals was found to be optimally around 0.2 Å, typically within 1 Å under practical experimental conditions at state-of-the-art X-ray free-electron laser facilities. Furthermore, we conducted an in-depth analysis of how various experimental factors, such as AuNP size, X-ray properties, and detector geometry, influence the accuracy of Bio-SOSS. This comprehensive investigation highlights the practicality and potential of Bio-SOSS in accurately capturing the diverse conformation spectrum of biological macromolecules, paving the way for deeper insights into their dynamic natures.

Cerium Photocatalyst in Action: Structural Dynamics in the Presence of Substrate Visualized via Time-Resolved X-ray Liquidography

[Ce(III)Cl6]3-, with its earth-abundant metal element, is a promising photocatalyst facilitating carbon-halogen bond activation. Still, the structure of the reaction intermediate has yet to be explored. Here, we applied time-resolved X-ray liquidography (TRXL), which allows for direct observation of the structural details of reaction intermediates, to investigate the photocatalytic reaction of [Ce(III)Cl6]3-. Structural analysis of the TRXL data revealed that the excited state of [Ce(III)Cl6]3- has Ce-Cl bonds that are shorter than those of the ground state and that the Ce-Cl bond further contracts upon oxidation. In addition, this study represents the first application of TRXL to both photocatalyst-only and photocatalyst-and-substrate samples, providing insights into the substrate's influence on the photocatalyst's reaction dynamics. This study demonstrates the capability of TRXL in elucidating the reaction dynamics of photocatalysts under various conditions and highlights the importance of experimental determination of the structures of reaction intermediates to advance our understanding of photocatalytic mechanisms.

Emerging Time-Resolved X-Ray Diffraction Approaches for Protein Dynamics

Proteins guide the flows of information, energy, and matter that make life possible by accelerating transport and chemical reactions, by allosterically modulating these reactions, and by forming dynamic supramolecular assemblies. In these roles, conformational change underlies functional transitions. Time-resolved X-ray diffraction methods characterize these transitions either by directly triggering sequences of functionally important motions or, more broadly, by capturing the motions of which proteins are capable. To date, most successful have been experiments in which conformational change is triggered in light-dependent proteins. In this review, I emphasize emerging techniques that probe the dynamic basis of function in proteins lacking natively light-dependent transitions and speculate about extensions and further possibilities. In addition, I review how the weaker and more distributed signals in these data push the limits of the capabilities of analytical methods. Taken together, these new methods are beginning to establish a powerful paradigm for the study of the physics of protein function.

Projection to extract the perpendicular component (PEPC) method for extracting kinetics from time-resolved data

Time-resolved x-ray liquidography (TRXL) is a potent method for investigating the structural dynamics of chemical and biological reactions in the liquid phase. It has enabled the extraction of detailed structural aspects of various dynamic processes, the molecular structures of intermediates, and kinetics of reactions across a wide range of systems, from small molecules to proteins and nanoparticles. Proper data analysis is key to extracting the information of the kinetics and structural dynamics of the studied system encrypted in the TRXL data. In typical TRXL data, the signals from solute scattering, solvent scattering, and solute-solvent cross scattering are mixed in the q-space, and the solute kinetics and solvent hydrodynamics are mixed in the time domain, thus complicating the data analysis. Various methods developed so far generally require prior knowledge of the molecular structures of candidate species involved in the reaction. Because such information is often unavailable, a typical data analysis often involves tedious trial and error. To remedy this situation, we have developed a method named projection to extract the perpendicular component (PEPC), capable of removing the contribution of solvent kinetics from TRXL data. The resulting data then contain only the solute kinetics, and, thus, the solute kinetics can be easily determined. Once the solute kinetics is determined, the subsequent data analysis to extract the structural information can be performed with drastically improved convenience. The application of the PEPC method is demonstrated with TRXL data from the photochemistry of two molecular systems: [Au(CN)₂^-]₃ in water and CHI₃ in cyclohexane.

Extracting Kinetics and Thermodynamics of Molecules without Heavy Atoms via Time-Resolved Solvent Scattering Signals

Time-resolved X-ray liquidography (TRXL) has emerged as a powerful technique for studying the structural dynamics of small molecules and macromolecules in liquid solutions. However, TRXL has limited sensitivity for small molecules containing light atoms only, whose signal has lower contrast compared with the signal from solvent molecules. Here, we present an alternative approach to bypass this limitation by detecting the change in solvent temperature resulting from a photoinduced reaction. Specifically, we analyzed the heat dynamics of TRXL data obtained from p-hydroxyphenacyl diethyl phosphate (HPDP). This analysis enabled us to experimentally determine the number of intermediates and their respective enthalpy changes, which can be compared to theoretical enthalpies to identify the intermediates. This work demonstrates that TRXL can be used to uncover the kinetics and reaction pathways for small molecules without heavy atoms even if the scattering signal from the solute molecules is buried under the strong solvent scattering signal.

Sample-minimizing co-flow cell for time-resolved pump-probe X-ray solution scattering

A fundamental problem in biological sciences is understanding how macromolecular machines work and how the structural changes of a molecule are connected to its function. Time-resolved techniques are vital in this regard and essential for understanding the structural dynamics of biomolecules. Time-resolved small- and wide-angle X-ray solution scattering has the capability to provide a multitude of information about the kinetics and global structural changes of molecules under their physiological conditions. However, standard protocols for such time-resolved measurements often require significant amounts of sample, which frequently render time-resolved measurements impossible. A cytometry-type sheath co-flow cell, developed at the BioCARS 14-ID beamline at the Advanced Photon Source, USA, allows time-resolved pump-probe X-ray solution scattering measurements to be conducted with sample consumption reduced by more than ten times compared with standard sample cells and protocols. The comparative capabilities of the standard and co-flow experimental setups were demonstrated by studying time-resolved signals in photoactive yellow protein.

Modeling difference x-ray scattering observations from an integral membrane protein within a detergent micelle

Time-resolved x-ray solution scattering (TR-XSS) is a sub-field of structural biology, which observes secondary structural changes in proteins as they evolve along their functional pathways. While the number of distinct conformational states and their rise and decay can be extracted directly from TR-XSS experimental data recorded from light-sensitive systems, structural modeling is more challenging. This step often builds from complementary structural information, including secondary structural changes extracted from crystallographic studies or molecular dynamics simulations. When working with integral membrane proteins, another challenge arises because x-ray scattering from the protein and the surrounding detergent micelle interfere and these effects should be considered during structural modeling. Here, we utilize molecular dynamics simulations to explicitly incorporate the x-ray scattering cross term between a membrane protein and its surrounding detergent micelle when modeling TR-XSS data from photoactivated samples of detergent solubilized bacteriorhodopsin. This analysis provides theoretical foundations in support of our earlier approach to structural modeling that did not explicitly incorporate this cross term and improves agreement between experimental data and theoretical predictions at lower x-ray scattering angles.

Light-induced protein structural dynamics in bacteriophytochrome revealed by time-resolved x-ray solution scattering

Bacteriophytochromes (BphPs) are photoreceptors that regulate a wide range of biological mechanisms via red light-absorbing (Pr)-to-far-red light-absorbing (Pfr) reversible photoconversion. The structural dynamics underlying Pfr-to-Pr photoconversion in a liquid solution phase are not well understood. We used time-resolved x-ray solution scattering (TRXSS) to capture light-induced structural transitions in the bathy BphP photosensory module of Pseudomonas aeruginosa. Kinetic analysis of the TRXSS data identifies three distinct structural species, which are attributed to lumi-F, meta-F, and Pr, connected by time constants of 95 μs and 21 ms. Structural analysis based on molecular dynamics simulations shows that the light activation of PaBphP accompanies quaternary structural rearrangements from an "II"-framed close form of the Pfr state to an "O"-framed open form of the Pr state in terms of the helical backbones. This study provides mechanistic insights into how modular signaling proteins such as BphPs transmit structural signals over long distances and regulate their downstream biological responses.

A setup for millisecond time-resolved X-ray solution scattering experiments at the CoSAXS beamline at the MAX IV Laboratory

The function of biomolecules is tightly linked to their structure, and changes therein. Time-resolved X-ray solution scattering has proven a powerful technique for interrogating structural changes and signal transduction in photoreceptor proteins. However, these only represent a small fraction of the biological macromolecules of interest. More recently, laser-induced temperature jumps have been introduced as a more general means of initiating structural changes in biomolecules. Here we present the development of a setup for millisecond time-resolved X-ray solution scattering experiments at the CoSAXS beamline, primarily using infrared laser light to trigger a temperature increase, and structural changes. We present results that highlight the characteristics of this setup along with data showing structural changes in lysozyme caused by a temperature jump. Further developments and applications of the setup are also discussed.

Reaction dynamics studied via femtosecond X-ray liquidography at X-ray free-electron lasers

X-ray free-electron lasers (XFELs) provide femtosecond X-ray pulses suitable for pump–probe time-resolved studies with a femtosecond time resolution. Since the advent of the first XFEL in 2009, recent years have witnessed a great number of applications with various pump–probe techniques at XFELs. Among these, time-resolved X-ray liquidography (TRXL) is a powerful method for visualizing structural dynamics in the liquid solution phase. Here, we classify various chemical and biological molecular systems studied via femtosecond TRXL (fs-TRXL) at XFELs, depending on the focus of the studied process, into (i) bond cleavage and formation, (ii) charge distribution and electron transfer, (iii) orientational dynamics, (iv) solvation dynamics, (v) coherent nuclear wavepacket dynamics, and (vi) protein structural dynamics, and provide a brief review on each category. We also lay out a plausible roadmap for future fs-TRXL studies for areas that have not been explored yet.

Femtosecond X-ray liquidography using X-ray free-electron lasers (XFELs) visualizes various aspects of reaction dynamics.

Reversible molecular motional switch based on circular photoactive protein oligomers exhibits unexpected photo-induced contraction

Molecular switches alterable between two stable states by environmental stimuli, such as light and temperature, offer the potential for controlling biological functions. Here, we report a circular photoswitchable protein complex made of multiple protein molecules that can rapidly and reversibly switch with significant conformational changes. The structural and photochromic properties of photoactive yellow protein (PYP) are harnessed to construct circular oligomer PYPs (coPYPs) of desired sizes. Considering the light-induced N-terminal protrusion of monomer PYP, we expected coPYPs would expand upon irradiation, but time-resolved X-ray scattering data reveal that the late intermediate has a pronounced light-induced contraction motion. This work not only provides an approach to engineering a novel protein-based molecular switch based on circular oligomers of well-known protein units but also demonstrates the importance of characterizing the structural dynamics of designed molecular switches.

Regulatory Switching by Concerted Motions on the Microsecond Time Scale of the Oxygen Sensor Protein FixL

Signal transduction proteins perceive external stimuli in their sensor module and regulate the biological activities of the effector module, allowing cellular adaptation in response to environmental changes. FixL is a dimeric heme protein kinase that senses the oxygen level in plant root nodules to regulate the transcription of nitrogen fixation genes via the phosphorylation of its cognate transcriptional activator. Dissociation of oxygen from the heme induces conformational changes in the protein, converting it from the inactive form for phosphorylation to the active form. However, how FixL undergoes conformational change to regulate kinase activity upon oxygen dissociation remains poorly understood. Here we report time-resolved ultraviolet resonance Raman spectra showing conformational changes for FixL from Sinorhizobium meliloti. We observed spectral changes with a time constant of about 3 μs, which were oxygen-specific. Furthermore, we found that the conformational changes in the sensor and kinase domains are coupled, enabling allosteric control of kinase activity. Our results demonstrate that concerted structural changes on the microsecond time scale serve as the regulatory switch in FixL.

Ultrafast coherent motion and helix rearrangement of homodimeric hemoglobin visualized with femtosecond X-ray solution scattering

Ultrafast motion of molecules, particularly the coherent motion, has been intensively investigated as a key factor guiding the reaction pathways. Recently, X-ray free-electron lasers (XFELs) have been utilized to elucidate the ultrafast motion of molecules. However, the studies on proteins using XFELs have been typically limited to the crystalline phase, and proteins in solution have rarely been investigated. Here we applied femtosecond time-resolved X-ray solution scattering (fs-TRXSS) and a structure refinement method to visualize the ultrafast motion of a protein. We succeeded in revealing detailed ultrafast structural changes of homodimeric hemoglobin involving the coherent motion. In addition to the motion of the protein itself, the time-dependent change of electron density of the hydration shell was tracked. Besides, the analysis on the fs-TRXSS data of myoglobin allows for observing the effect of the oligomeric state on the ultrafast coherent motion.

Femtosecond time-resolved X-ray solution scattering (fs-TRXSS) measurements provide information on the structural dynamics of proteins in solution. Here, the authors present a structure refinement method for the analysis of fs-TRXSS data and use it to characterise the ultrafast structural changes of homodimeric haemoglobin.

Femtosecond X-ray spectroscopy of haem proteins

We discuss our recently reported femtosecond (fs) X-ray emission spectroscopy results on the ligand dissociation and recombination in nitrosylmyoglobin (MbNO) in the context of previous studies on ferrous haem proteins. We also present a preliminary account of femtosecond X-ray absorption studies on MbNO, pointing to the presence of more than one species formed upon photolysis.

Effect of the abolition of intersubunit salt bridges on allosteric protein structural dynamics

A salt bridge, one of the representative structural factors established by non-covalent interactions, plays a crucial role in stabilizing the structure and regulating the protein function, but its role in dynamic processes has been elusive. Here, to scrutinize the structural and functional roles of the salt bridge in the process of performing the protein function, we investigated the effects of salt bridges on the allosteric structural transition of homodimeric hemoglobin (HbI) by applying time-resolved X-ray solution scattering (TRXSS) to the K30D mutant, in which the interfacial salt bridges of the wild type (WT) are abolished. The TRXSS data of K30D are consistent with the kinetic model that requires one monomer intermediate in addition to three structurally distinct dimer intermediates (I1, I2, and I3) observed in WT and other mutants. The kinetic and structural analyses show that K30D has an accelerated biphasic transition from I2 to I3 by more than nine times compared to WT and lacks significant structural changes in the transition from R-like I2 to T-like I3 observed in WT, unveiling that the loss of the salt bridges interrupts the R-T allosteric transition of HbI. Besides, the correlation between the bimolecular CO recombination rates in K30D, WT, and other mutants reveals that the bimolecular CO recombination is abnormally decelerated in K30D, indicating that the salt bridges also affect the cooperative ligand binding in HbI. These comparisons of the structural dynamics and kinetics of K30D and WT show that the interfacial salt bridges not only assist the physical connection of two subunits but also play a critical role in the global structural signal transduction of one subunit to the other subunit via a series of well-organized structural transitions.

Energy Transfer across Nonpolar and Polar Contacts in Proteins: Role of Contact Fluctuations

Molecular dynamics simulations of the villin headpiece subdomain HP36 have been carried out to examine relations between rates of vibrational energy transfer across non-covalently bonded contacts and equilibrium structural fluctuations, with focus on van der Waals contacts. Rates of energy transfer across van der Waals contacts vary inversely with the variance of the contact length, with the same constant of proportionality for all nonpolar contacts of HP36. A similar relation is observed for hydrogen bonds, but the proportionality depends on contact pairs, with hydrogen bonds stabilizing the α-helices all exhibiting the same constant of proportionality, one that is distinct from those computed for other polar contacts. Rates of energy transfer across van der Waals contacts are found to be up to 2 orders of magnitude smaller than rates of energy transfer across polar contacts.

Real-time tracking of protein unfolding with time-resolved x-ray solution scattering

The correct folding of proteins is of paramount importance for their function, and protein misfolding is believed to be the primary cause of a wide range of diseases. Protein folding has been investigated with time-averaged methods and time-resolved spectroscopy, but observing the structural dynamics of the unfolding process in real-time is challenging. Here, we demonstrate an approach to directly reveal the structural changes in the unfolding reaction. We use nano- to millisecond time-resolved x-ray solution scattering to probe the unfolding of apomyoglobin. The unfolding reaction was triggered using a temperature jump, which was induced by a nanosecond laser pulse. We demonstrate a new strategy to interpret time-resolved x-ray solution scattering data, which evaluates ensembles of structures obtained from molecular dynamics simulations. We find that apomyoglobin passes three states when unfolding, which we characterize as native, molten globule, and unfolded. The molten globule dominates the population under the conditions investigated herein, whereas native and unfolded structures primarily contribute before the laser jump and 30 μs after it, respectively. The molten globule retains much of the native structure but shows a dynamic pattern of inter-residue contacts. Our study demonstrates a new strategy to directly observe structural changes over the cause of the unfolding reaction, providing time- and spatially resolved atomic details of the folding mechanism of globular proteins.

Integrating solvation shell structure in experimentally driven molecular dynamics using x-ray solution scattering data

In the past few decades, prediction of macromolecular structures beyond the native conformation has been aided by the development of molecular dynamics (MD) protocols aimed at exploration of the energetic landscape of proteins. Yet, the computed structures do not always agree with experimental observables, calling for further development of the MD strategies to bring the computations and experiments closer together. Here, we report a scalable, efficient MD simulation approach that incorporates an x-ray solution scattering signal as a driving force for the conformational search of stable structural configurations outside of the native basin. We further demonstrate the importance of inclusion of the hydration layer effect for a precise description of the processes involving large changes in the solvent exposed area, such as unfolding. Utilization of the graphics processing unit allows for an efficient all-atom calculation of scattering patterns on-the-fly, even for large biomolecules, resulting in a speed-up of the calculation of the associated driving force. The utility of the methodology is demonstrated on two model protein systems, the structural transition of lysine-, arginine-, ornithine-binding protein and the folding of deca-alanine. We discuss how the present approach will aid in the interpretation of dynamical scattering experiments on protein folding and association.

A tool for visualizing protein motions in time-resolved crystallography

Time-resolved serial femtosecond crystallography (TR-SFX) at an x-ray free electron laser enables protein structural changes to be imaged on time-scales from femtoseconds to seconds. It can, however, be difficult to grasp the nature and timescale of global protein motions when structural changes are not isolated near a single active site. New tools are, therefore, needed to represent the global nature of electron density changes and their correlation with modeled protein structural changes. Here, we use TR-SFX data from bacteriorhodopsin to develop and validate a method for quantifying time-dependent electron density changes and correlating them throughout the protein. We define a spherical volume of difference electron density about selected atoms, average separately the positive and negative electron difference densities within each volume, and walk this spherical volume through all atoms within the protein. By correlating the resulting difference electron density amplitudes with time, our approach facilitates an initial assessment of the number and timescale of structural intermediates and highlights quake-like motions on the sub-picosecond timescale. This tool also allows structural models to be compared with experimental data using theoretical difference electron density changes calculated from refined resting and photo-activated structures.

Effect of Occluded Ligand Migration on the Kinetics and Structural Dynamics of Homodimeric Hemoglobin

Small molecules such as molecular oxygen, nitric oxide, and carbon monoxide play important roles in life, and many proteins require the transport of small molecules to and from the bulk solvent for their function. Ligand migration within a protein molecule is expected to be closely related to the overall structural changes of the protein, but the detailed and quantitative connection remains elusive. For example, despite numerous studies, how occluded ligand migration affects the kinetics and structural dynamics of the R-T transition remains unclear. To shed light on this issue, we chose homodimeric hemoglobin (HbI) with the I114F mutation (I114F), which is known to interfere with ligand migration between the primary and secondary docking sites, and studied its kinetics and structural dynamics using time-resolved X-ray solution scattering. The kinetic analysis shows that I114F has three structurally distinct intermediates (I₁, I₂, and I₃) as in the wild type (WT), but its geminate CO recombination occurs directly from I₁ without the path via I₂ observed in WT. Moreover, the structural transitions, which involve ligand migration (the transitions from I₁ to I₂ and from I₃ to the initial state), are decelerated compared to WT. The structural analysis revealed that I114F involves generally smaller structural changes in all three intermediates compared to WT.

Model Reconstruction from Small-Angle X-Ray Scattering Data Using Deep Learning Methods

Small-angle X-ray scattering (SAXS) method is widely used in investigating protein structures in solution, but high-quality 3D model reconstructions are challenging. We present a new algorithm based on a deep learning method for model reconstruction from SAXS data. An auto-encoder for protein 3D models was trained to compress 3D shape information into vectors of a 200-dimensional latent space, and the vectors are optimized using genetic algorithms to build 3D models that are consistent with the scattering data. The program has been tested with experimental SAXS data, demonstrating the capacity and robustness of accurate model reconstruction. Furthermore, the model size information can be optimized using this algorithm, enhancing the automation in model reconstruction directly from SAXS data. The program was implemented using Python with the TensorFlow framework, with source code and webserver available from http://liulab.csrc.ac.cn/decodeSAXS.

•

A convolutional neural network auto-encoder framework for 3D models is developed

•

The auto-encoder compresses protein shape information to 200 parameters

•

Accurate 3D models (both shape and radius) can be reconstructed from 1D SAXS data

Revisiting the nitrite reductase activity of hemoglobin with differential pulse voltammetry

Nitric oxide (NO) is an omnipresent signalling molecule in all vertebrates. NO modulates blood flow and neural activity. Nitrite anion is one of the most important sources of NO. Nitrite is reduced to NO by various physiological mechanisms including reduction by hemoglobin in vascular system. In this study, nitrite reductase activity (NRA) of hemoglobin is reported using cyclic voltammetry (CV) and differential pulse voltammetry (DPV) in a wide potential window from +0.3 V to -1.3 V (vs. Ag/AgCl). To the best of our knowledge, a detailed look into NRA of hemoglobin is proposed here for the first time. Our results indicated two different regimes for reduction of nitrite by hemoglobin in its Fe(II) and Fe(I) states. Both reactions showed a reversible behaviour in the time scale of the experiments. The first reduction displayed a normal redox behaviour, while the latter one had the characteristics of a catalytic electro-reduction/oxidation. The reduction in Fe(II) state was selected as a tool for comparing the NRA of hemoglobin (Hb) and hemoglobin-S (Hb-S) under native-like conditions in a didodecyldimethyl ammonium bromide (DDAB) liquid crystal film. These investigations lay the prospects and guidelines for understanding the direct electrochemistry of hemoglobin utilizing a simplified mediator-free platform.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Broad-Band Ultraviolet CD Spectroscopy of Ultrafast Peptide Backbone Conformational Dynamics

The far-UV spectral window widely used for the conformational analysis of biomolecules is not easily covered with broad-band lasers. This has made it difficult to use circular dichroism (CD) spectroscopy to directly follow fast structure changes. By combining transient CD spectroscopy in the deep-UV with thioamide substitution, we demonstrate a method to overcome this difficulty. We investigated a dipeptide whose two carbonyl oxygen atoms were replaced by sulfur, red-shifting the strong lowest-lying ππ* transitions into the more accessible 250-370 nm spectral window. Coupling of the two thioamide units cannot be resolved by achiral 2D-UV spectroscopy, but it gives rise to a pronounced bisignate CD spectrum. The transient CD spectra reveal weakening of this coupling in the electronically excited state, where conformational constraints are released. Our results show that direct local probing of fast backbone conformational change via CD spectroscopy is possible in combination with site-selective thio substitution in peptides and proteins.

SVD-aided non-orthogonal decomposition (SANOD) method to exploit prior knowledge of spectral components in the analysis of time-resolved data

Analysis of time-resolved data typically involves discriminating noise against the signal and extracting time-independent components and their time-dependent contributions. Singular value decomposition (SVD) serves this purpose well, but the extracted time-independent components are not necessarily the physically meaningful spectra directly representing the actual dynamic or kinetic processes but rather a mathematically orthogonal set necessary for constituting the physically meaningful spectra. Converting the orthogonal components into physically meaningful spectra requires subsequent posterior analyses such as linear combination fitting (LCF) and global fitting (GF), which takes advantage of prior knowledge about the data but requires that all components are known or satisfactory components are guessed. Since in general not all components are known, they have to be guessed and tested via trial and error. In this work, we introduce a method, which is termed SVD-aided Non-Orthogonal Decomposition (SANOD), to circumvent trial and error. The key concept of SANOD is to combine the orthogonal components from SVD with the known prior knowledge to fill in the gap of the unknown signal components and to use them for LCF. We demonstrate the usefulness of SANOD via applications to a variety of cases.

Do Osmolytes Impact the Structure and Dynamics of Myoglobin?

Osmolytes are small organic compounds that can affect the stability of proteins in living cells. The mechanism of osmolytes' protective effects on protein structure and dynamics has not been fully explained, but in general, two possibilities have been suggested and examined: a direct interaction of osmolytes with proteins (water replacement hypothesis), and an indirect interaction (vitrification hypothesis). Here, to investigate these two possible mechanisms, we studied myoglobin-osmolyte systems using FTIR, UV-vis, CD, and femtosecond IR pump-probe spectroscopy. Interestingly, noticeable changes are observed in both the lifetime of the CO stretch of CO-bound myoglobin and the spectra of UV-vis, CD, and FTIR upon addition of the osmolytes. In addition, the temperature-dependent CD studies reveal that the protein's thermal stability depends on molecular structure, hydrogen-bonding ability, and size of osmolytes. We anticipate that the present experimental results provide important clues about the complicated and intricate mechanism of osmolyte effects on protein structure and dynamics in a crowded cellular environment.

Protein Structural Dynamics of Wild-Type and Mutant Homodimeric Hemoglobin Studied by Time-Resolved X-Ray Solution Scattering

The quaternary transition between the relaxed (R) and tense (T) states of heme-binding proteins is a textbook example for the allosteric structural transition. Homodimeric hemoglobin (HbI) from Scapharca inaequivalvis is a useful model system for investigating the allosteric behavior because of the relatively simple quaternary structure. To understand the cooperative transition of HbI, wild-type and mutants of HbI have been studied by using time-resolved X-ray solution scattering (TRXSS), which is sensitive to the conformational changes. Herein, we review the structural dynamics of HbI investigated by TRXSS and compare the results of TRXSS with those of other techniques.

Insulin hexamer dissociation dynamics revealed by photoinduced T-jumps and time-resolved X-ray solution scattering

The structural dynamics of insulin hexamer dissociation were studied by the photoinduced temperature jump technique and monitored by time-resolved X-ray scattering. The process of hexamer dissociation was found to involve several transient intermediates, including an expanded hexamer and an unstable tetramer. Our findings provide insights into the mechanisms of protien-protein association.

Probing Cytochrome c Folding Transitions upon Phototriggered Environmental Perturbations Using Time-Resolved X-ray Scattering

Direct tracking of protein structural dynamics during folding-unfolding processes is important for understanding the roles of hierarchic structural factors in the formation of functional proteins. Using cytochrome c (cyt c) as a platform, we investigated its structural dynamics during folding processes triggered by local environmental changes (i.e., pH or heme iron center oxidation/spin/ligation states) with time-resolved X-ray solution scattering measurements. Starting from partially unfolded cyt c, a sudden pH drop initiated by light excitation of a photoacid caused a structural contraction in microseconds, followed by active site restructuring and unfolding in milliseconds. In contrast, the reduction of iron in the heme via photoinduced electron transfer did not affect conformational stability at short timescales (<1 ms), despite active site coordination geometry changes. These results demonstrate how different environmental perturbations can change the nature of interaction between the active site and protein conformation, even within the same metalloprotein, which will subsequently affect the folding structural dynamics.

Transient isomers in the photodissociation of bromoiodomethane

The photochemistry of halomethanes is fascinating for the complex cascade reactions toward either the parent or newly synthesized molecules. Here, we address the structural rearrangement of photodissociated CH₂IBr in methanol and cyclohexane, probed by time-resolved X-ray scattering in liquid solution. Upon selective laser cleavage of the C-I bond, we follow the reaction cascade of the two geminate geometrical isomers, CH₂I-Br and CH₂Br-I. Both meta-stable isomers decay on different time scales, mediated by solvent interaction, toward the original parent molecule. We observe the internal rearrangement of CH₂Br-I to CH₂I-Br in cyclohexane by extending the time window up to 3 μs. We track the photoproduct kinetics of CH₂Br-I in methanol solution where only one isomer is observed. The effect of the polarity of solvent on the geminate recombination pathways is discussed.

Direct Observation of Insulin Association Dynamics with Time-Resolved X-ray Scattering

Biological functions frequently require protein-protein interactions that involve secondary and tertiary structural perturbation. Here we study protein-protein dissociation and reassociation dynamics in insulin, a model system for protein oligomerization. Insulin dimer dissociation into monomers was induced by a nanosecond temperature-jump (T-jump) of ∼8 °C in aqueous solution, and the resulting protein and solvent dynamics were tracked by time-resolved X-ray solution scattering (TRXSS) on time scales of 10 ns to 100 ms. The protein scattering signals revealed the formation of five distinguishable transient species during the association process that deviate from simple two-state kinetics. Our results show that the combination of T-jump pump coupled to TRXSS probe allows for direct tracking of structural dynamics in nonphotoactive proteins.

SVD-aided pseudo principal-component analysis: A new method to speed up and improve determination of the optimum kinetic model from time-resolved data

Determination of the optimum kinetic model is an essential prerequisite for characterizing dynamics and mechanism of a reaction. Here, we propose a simple method, termed as singular value decomposition-aided pseudo principal-component analysis (SAPPA), to facilitate determination of the optimum kinetic model from time-resolved data by bypassing any need to examine candidate kinetic models. We demonstrate the wide applicability of SAPPA by examining three different sets of experimental time-resolved data and show that SAPPA can efficiently determine the optimum kinetic model. In addition, the results of SAPPA for both time-resolved X-ray solution scattering (TRXSS) and transient absorption (TA) data of the same protein reveal that global structural changes of protein, which is probed by TRXSS, may occur more slowly than local structural changes around the chromophore, which is probed by TA spectroscopy.

Ultrafast X-Ray Crystallography and Liquidography

Time-resolved X-ray diffraction provides direct information on three-dimensional structures of reacting molecules and thus can be used to elucidate structural dynamics of chemical and biological reactions. In this review, we discuss time-resolved X-ray diffraction on small molecules and proteins with particular emphasis on its application to crystalline (crystallography) and liquid-solution (liquidography) samples. Time-resolved X-ray diffraction has been used to study picosecond and slower dynamics at synchrotrons and can now access even femtosecond dynamics with the recent arrival of X-ray free-electron lasers.

Interfacial water effect on cooperativity and signal communication in Scapharca dimeric hemoglobin

Cooperativity is important in controlling the biological functions of allosteric proteins.

Identifying and Visualizing Macromolecular Flexibility in Structural Biology

Structural biology comprises a variety of tools to obtain atomic resolution data for the investigation of macromolecules. Conventional structural methodologies including crystallography, NMR and electron microscopy often do not provide sufficient details concerning flexibility and dynamics, even though these aspects are critical for the physiological functions of the systems under investigation. However, the increasing complexity of the molecules studied by structural biology (including large macromolecular assemblies, integral membrane proteins, intrinsically disordered systems, and folding intermediates) continuously demands in-depth analyses of the roles of flexibility and conformational specificity involved in interactions with ligands and inhibitors. The intrinsic difficulties in capturing often subtle but critical molecular motions in biological systems have restrained the investigation of flexible molecules into a small niche of structural biology. Introduction of massive technological developments over the recent years, which include time-resolved studies, solution X-ray scattering, and new detectors for cryo-electron microscopy, have pushed the limits of structural investigation of flexible systems far beyond traditional approaches of NMR analysis. By integrating these modern methods with powerful biophysical and computational approaches such as generation of ensembles of molecular models and selective particle picking in electron microscopy, more feasible investigations of dynamic systems are now possible. Using some prominent examples from recent literature, we review how current structural biology methods can contribute useful data to accurately visualize flexibility in macromolecular structures and understand its important roles in regulation of biological processes.

Ultrafast anisotropic protein quake propagation after CO photodissociation in myoglobin

"Protein quake" denotes the dissipation of excess energy across a protein, in response to a local perturbation such as the breaking of a chemical bond or the absorption of a photon. Femtosecond time-resolved small- and wide-angle X-ray scattering (TR-SWAXS) is capable of tracking such ultrafast protein dynamics. However, because the structural interpretation of the experiments is complicated, a molecular picture of protein quakes has remained elusive. In addition, new questions arose from recent TR-SWAXS data that were interpreted as underdamped oscillations of an entire protein, thus challenging the long-standing concept of overdamped global protein dynamics. Based on molecular-dynamics simulations, we present a detailed molecular movie of the protein quake after carbon monoxide (CO) photodissociation in myoglobin. The simulations suggest that the protein quake is characterized by a single pressure peak that propagates anisotropically within 500 fs across the protein and further into the solvent. By computing TR-SWAXS patterns from the simulations, we could interpret features in the reciprocal-space SWAXS signals as specific real-space dynamics, such as CO displacement and pressure wave propagation. Remarkably, we found that the small-angle data primarily detect modulations of the solvent density but not oscillations of the bare protein, thereby reconciling recent TR-SWAXS experiments with the notion of overdamped global protein dynamics.

Femtosecond X-ray solution scattering reveals that bond formation mechanism of a gold trimer complex is independent of excitation wavelength

The [Au(CN)2 (-)]3 trimer in water experiences a strong van der Waals interaction between the d(10) gold atoms due to large relativistic effect and can serve as an excellent model system to study the bond formation process in real time. The trimer in the ground state (S0) exists as a bent structure without the covalent bond between the gold atoms, and upon the laser excitation, one electron in the antibonding orbital goes to the bonding orbital, thereby inducing the formation of a covalent bond between gold atoms. This process has been studied by various time-resolved techniques, and most of the interpretation on the structure and dynamics converge except that the structure of the first intermediate (S1) has been debated due to different interpretations between femtosecond optical spectroscopy and femtosecond X-ray solution scattering. Recently, the excitation wavelength of 267 nm employed in our previous scattering experiment was suggested as the culprit for misinterpretation. Here, we revisited this issue by performing femtosecond X-ray solution scattering with 310 nm excitation and compared the results with our previous study employing 267 nm excitation. The data show that a linear S1 structure is formed within 500 fs regardless of excitation wavelength and the structural dynamics observed at both excitation wavelengths are identical to each other within experimental errors.

Application of advanced X-ray methods in life sciences

BACKGROUND

Synchrotron radiation (SR) sources provide diverse X-ray methods for the investigation of structure-function relationships in biological macromolecules.

SCOPE OF REVIEW

Recent developments in SR sources and in the X-ray tools they offer for life sciences are reviewed. Specifically, advances in macromolecular crystallography, small angle X-ray solution scattering, X-ray absorption and fluorescence spectroscopy, and imaging are discussed with examples.

MAJOR CONCLUSIONS

SR sources offer a range of X-ray techniques that can be used in a complementary fashion in studies of biological systems at a wide range of resolutions from atomic to cellular scale. Emerging applications of X-ray techniques include the characterization of disordered proteins, noncrystalline and nonequilibrium systems, elemental imaging of tissues, cells and organs, and detection of time-resolved changes in molecular structures.

GENERAL SIGNIFICANCE

X-ray techniques are in the center of hybrid approaches that are used to gain insight into complex problems relating to biomolecular mechanisms, disease and possible therapeutic solutions. This article is part of a Special Issue entitled "Science for Life". Guest Editors: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Cooperative protein structural dynamics of homodimeric hemoglobin linked to water cluster at subunit interface revealed by time-resolved X-ray solution scattering

Homodimeric hemoglobin (HbI) consisting of two subunits is a good model system for investigating the allosteric structural transition as it exhibits cooperativity in ligand binding. In this work, as an effort to extend our previous study on wild-type and F97Y mutant HbI, we investigate structural dynamics of a mutant HbI in solution to examine the role of well-organized interfacial water cluster, which has been known to mediate intersubunit communication in HbI. In the T72V mutant of HbI, the interfacial water cluster in the T state is perturbed due to the lack of Thr72, resulting in two less interfacial water molecules than in wild-type HbI. By performing picosecond time-resolved X-ray solution scattering experiment and kinetic analysis on the T72V mutant, we identify three structurally distinct intermediates (I1, I2, and I3) and show that the kinetics of the T72V mutant are well described by the same kinetic model used for wild-type and F97Y HbI, which involves biphasic kinetics, geminate recombination, and bimolecular CO recombination. The optimized kinetic model shows that the R-T transition and bimolecular CO recombination are faster in the T72V mutant than in the wild type. From structural analysis using species-associated difference scattering curves for the intermediates, we find that the T-like deoxy I3 intermediate in solution has a different structure from deoxy HbI in crystal. In addition, we extract detailed structural parameters of the intermediates such as E-F distance, intersubunit rotation angle, and heme-heme distance. By comparing the structures of protein intermediates in wild-type HbI and the T72V mutant, we reveal how the perturbation in the interfacial water cluster affects the kinetics and structures of reaction intermediates of HbI.

Tracking reaction dynamics in solution by pump–probe X-ray absorption spectroscopy and X-ray liquidography (solution scattering)

TRXL and TRXAS are powerful techniques for real-time probing of structural and electronic dynamics of photoinduced reactions in solution phase.

Ligand and interfacial dynamics in a homodimeric hemoglobin

The structural dynamics of dimeric hemoglobin (HbI) from Scapharca inaequivalvis in different ligand-binding states is studied from atomistic simulations on the μs time scale. The intermediates are between the fully ligand-bound (R) and ligand-free (T) states. Tertiary structural changes, such as rotation of the side chain of Phe97, breaking of the Lys96-heme salt bridge, and the Fe-Fe separation, are characterized and the water dynamics along the R-T transition is analyzed. All these properties for the intermediates are bracketed by those determined experimentally for the fully ligand-bound and ligand-free proteins, respectively. The dynamics of the two monomers is asymmetric on the 100 ns timescale. Several spontaneous rotations of the Phe97 side chain are observed which suggest a typical time scale of 50-100 ns for this process. Ligand migration pathways include regions between the B/G and C/G helices and, if observed, take place in the 100 ns time scale.

Identifying the major intermediate species by combining time-resolved X-ray solution scattering and X-ray absorption spectroscopy

Identifying the intermediate species along a reaction pathway is a first step towards a complete understanding of the reaction mechanism, but often this task is not trivial. There has been a strong on-going debate: which of the three intermediates, the CHI2 radical, the CHI2-I isomer, and the CHI2(+) ion, is the dominant intermediate species formed in the photolysis of iodoform (CHI3)? Herein, by combining time-resolved X-ray liquidography (TRXL) and time-resolved X-ray absorption spectroscopy (TR-XAS), we present strong evidence that the CHI2 radical is dominantly formed from the photolysis of CHI3 in methanol at 267 nm within the available time resolution of the techniques (∼20 ps for TRXL and ∼100 ps for TR-XAS). The TRXL measurement, conducted using the time-slicing scheme, detected no CHI2-I isomer within our signal-to-noise ratio, indicating that, if formed, the CHI2-I isomer must be a minor intermediate. The TR-XAS transient spectra measured at the iodine L1 and L3 edges support the same conclusion. The present work demonstrates that the application of these two complementary time-resolved X-ray methods to the same system can provide a detailed understanding of the reaction mechanism.