CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase.

Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.

Protocol using similarity score and improved shrink-wrap algorithm for better convergence of phase-retrieval calculation in X-ray diffraction imaging.

In X-ray diffraction imaging (XDI), electron density maps of a targeted particle are reconstructed computationally from the diffraction pattern alone using phase-retrieval (PR) algorithms. However, the PR calculations sometimes fail to yield realistic electron density maps that approximate the structure of the particle. This occurs due to the absence of structure amplitudes at and near the zero-scattering angle and the presence of Poisson noise in weak diffraction patterns. Consequently, the PR calculation becomes a bottleneck for XDI structure analyses. Here, a protocol to efficiently yield realistic maps is proposed. The protocol is based on the empirical observation that realistic maps tend to yield low similarity scores, as suggested in our prior study [Sekiguchi et al. (2017), J. Synchrotron Rad. 24, 1024-1038]. Among independently and concurrently executed PR calculations, the protocol modifies all maps using the electron-density maps exhibiting low similarity scores. This approach, along with a new protocol for estimating particle shape, improved the probability of obtaining realistic maps for diffraction patterns from various aggregates of colloidal gold particles, as compared with PR calculations performed without the protocol. Consequently, the protocol has the potential to reduce computational costs in PR calculations and enable efficient XDI structure analysis of non-crystalline particles using synchrotron X-rays and X-ray free-electron laser pulses.

Similarity score for screening phase-retrieved maps in X-ray diffraction imaging - characterization in reciprocal space.

X-ray diffraction imaging (XDI) is utilized for visualizing the structures of non-crystalline particles in material sciences and biology. In the structural analysis, phase-retrieval (PR) algorithms are applied to the diffraction amplitude data alone to reconstruct the electron density map of a specimen particle projected along the direction of the incident X-rays. However, PR calculations may not lead to good convergence because of a lack of diffraction patterns in small-angle regions and Poisson noise in X-ray detection. Therefore, the PR calculation is still a bottleneck for the efficient application of XDI in the structural analyses of non-crystalline particles. For screening maps from hundreds of trial PR calculations, we have been using a score and measuring the similarity between a pair of retrieved maps. Empirically, probable maps approximating the particle structures gave a score smaller than a threshold value, but the reasons for the effectiveness of the score are still unclear. In this study, the score is characterized in terms of the phase differences between the structure factors of the retrieved maps, the usefulness of the score in screening the maps retrieved from experimental diffraction patterns is demonstrated, and the effective resolution of similarity-score-selected maps is discussed.

Coenzyme-binding pathway on glutamate dehydrogenase suggested from multiple-binding sites visualized by cryo-electron microscopy.

The structure of hexameric glutamate dehydrogenase (GDH) in the presence of the coenzyme nicotinamide adenine dinucleotide phosphate (NADP) was visualized using cryogenic transmission electron microscopy to investigate the ligand-binding pathways to the active site of the enzyme. Each subunit of GDH comprises one hexamer-forming core domain and one nucleotide-binding domain (NAD domain), which spontaneously opens and closes the active-site cleft situated between the two domains. In the presence of NADP, the potential map of GDH hexamer, assuming D3 symmetry, was determined at a resolution of 2.4 Å, but the NAD domain was blurred due to the conformational variety. After focused classification with respect to the NAD domain, the potential maps interpreted as NADP molecules appeared at five different sites in the active-site cleft. The subunits associated with NADP molecules were close to one of the four metastable conformations in the unliganded state. Three of the five binding sites suggested a pathway of NADP molecules to approach the active-site cleft for initiating the enzymatic reaction. The other two binding modes may rarely appear in the presence of glutamate, as demonstrated by the reaction kinetics. Based on the visualized structures and the results from the enzymatic kinetics, we discussed the binding modes of NADP to GDH in the absence and presence of glutamate.

Ultrastructure and fractal property of chromosomes in close-to-native yeast nuclei visualized using X-ray laser diffraction.

Genome compaction and activity in the nucleus depend on spatiotemporal changes in the organization of chromatins in chromosomes. However, the direct imaging of the chromosome structures in the nuclei has been difficult and challenging. Herein, we directly visualized the structure of chromosomes in frozen-hydrated nuclei of budding yeast in the interphase using X-ray laser diffraction. The reconstructed projection electron density maps revealed inhomogeneous distributions of chromosomes, such as a 300 nm assembly and fibrous substructures in the elliptic-circular shaped nuclei of approximately 800 nm. In addition, from the diffraction patterns, we confirmed the absence of regular arrangements of chromosomes and chromatins with 400-20 nm spacing, and demonstrated that chromosomes were composed of self-similarly assembled substructural domains with an average radius of gyration of 58 nm and smooth surfaces. Based on these analyses, we constructed putative models to discuss the organization of 16 chromosomes, carrying DNA of 4.1 mm in 800 nm ellipsoid of the nucleus at the interphase. We anticipate the structural parameters on the fractal property of chromosomes and the experimental images to be a starting point for constructing more sophisticated 3D structural models of the nucleus.

Prediction of hydrophilic and hydrophobic hydration structure of protein by neural network optimized using experimental data.

The hydration structures of proteins, which are necessary for their folding, stability, and functions, were visualized using X-ray and neutron crystallography and transmission electron microscopy. However, complete visualization of hydration structures over the entire protein surface remains difficult. To compensate for this incompleteness, we developed a three-dimensional convolutional neural network to predict the probability distribution of hydration water molecules on the hydrophilic and hydrophobic surfaces, and in the cavities of proteins. The neural network was optimized using the distribution patterns of protein atoms around the hydration water molecules identified in the high-resolution X-ray crystal structures. We examined the feasibility of the neural network using water sites in the protein crystal structures that were not included in the datasets. The predicted distribution covered most of the experimentally identified hydration sites, with local maxima appearing in their vicinity. This computational approach will help to highlight the relevance of hydration structures to the biological functions and dynamics of proteins.

Red light-induced structure changes in phytochrome A from Pisum sativum.

Phytochrome A (phyA) is a photoreceptor protein of plants that regulates the red/far-red light photomorphogenic responses of plants essential for growth and development. PhyA, composed of approximately 1100 amino acid residues, folds into photosensory and output signaling modules. The photosensory module covalently binds phytochromobilin as a chromophore for photoreversible interconversion between inactive red light-absorbing (Pr) and active far-red light-absorbing (Pfr) forms to act as a light-driven phosphorylation enzyme. To understand the molecular mechanism in the initial process of photomorphogenic response, we studied the molecular structures of large phyA (LphyA) from Pisum sativum, which lacks the 52 residues in the N-terminal, by small-angle X-ray scattering combined with multivariate analyses applied to molecular models predicted from the scattering profiles. According to our analyses, Pr was in a dimer and had a four-leaf shape, and the subunit was approximated as a bent rod of 175 × 50 Å. The scattering profile of Pfr was calculated from that recorded for a mixture of Pr and Pfr under red-light irradiation by using their population determined from the absorption spectrum. The Pfr dimer exhibited a butterfly shape composed of subunits with a straight rod of 175 × 50 Å. The shape differences between Pr and Pfr indicated conformational changes in the Pr/Pfr interconversion which would be essential to the interaction with protein molecules involved in transcriptional control.

Common architectures in cyanobacteria Prochlorococcus cells visualized by X-ray diffraction imaging using X-ray free electron laser.

Visualization of intracellular structures and their spatial organization inside cells without any modification is essential to understand the mechanisms underlying the biological functions of cells. Here, we investigated the intracellular structure of cyanobacteria Prochlorococcus in the interphase by X-ray diffraction imaging using X-ray free-electron laser. A number of diffraction patterns from single cells smaller than 1 µm in size were collected with high signal-to-noise ratio with a resolution of up to 30 nm. From diffraction patterns, a set of electron density maps projected along the direction of the incident X-ray were retrieved with high reliability. The most characteristic structure found to be common among the cells was a C-shaped arrangement of 100-nm sized high-density spots, which surrounded a low-density area of 100 nm. Furthermore, a three-dimensional map reconstructed from the projection maps of individual cells was non-uniform, indicating the presence of common structures among cyanobacteria cells in the interphase. By referring to the fluorescent images for distributions of thylakoid membranes, nucleoids, and carboxysomes, we inferred and represented their spatial arrangements in the three-dimensional map. The arrangement allowed us to discuss the relevance of the intracellular organization to the biological functions of cyanobacteria.

Domain Organization in Plant Blue-Light Receptor Phototropin2 of Arabidopsis thaliana Studied by Small-Angle X-ray Scattering.

Phototropin2 (phot2) is a blue-light (BL) receptor protein that regulates the BL-dependent activities of plants for efficient photosynthesis. Phot2 is composed of two light-oxygen-voltage sensing domains (LOV1 and LOV2) to absorb BL, and a kinase domain. Photo-activated LOV domains, especially LOV2, play a major role in photo-dependent increase in the phosphorylation activity of the kinase domain. The atomic details of the overall structure of phot2 and the intramolecular mechanism to convert BL energy to a phosphorylation signal remain unknown. We performed structural studies on the LOV fragments LOV1, LOV2, LOV2-linker, and LOV2-kinase, and full-length phot2, using small-angle X-ray scattering (SAXS). The aim of the study was to understand structural changes under BL irradiation and discuss the molecular mechanism that enhance the phosphorylation activity under BL. SAXS is a suitable technique for visualizing molecular structures of proteins in solution at low resolution and is advantageous for monitoring their structural changes in the presence of external physical and/or chemical stimuli. Structural parameters and molecular models of the recombinant specimens were obtained from SAXS profiles in the dark, under BL irradiation, and after dark reversion. LOV1, LOV2, and LOV2-linker fragments displayed minimal structural changes. However, BL-induced rearrangements of functional domains were noted for LOV2-kinase and full-length phot2. Based on the molecular model together with the absorption measurements and biochemical assays, we discuss the intramolecular interactions and domain motions necessary for BL-enhanced phosphorylation activity of phot2.

Assessment of Force Field Accuracy Using Cryogenic Electron Microscopy Data of Hyper-thermostable Glutamate Dehydrogenase.

Molecular dynamics (MD) simulations in biophysically relevant time scales of microseconds is a powerful tool for studying biomolecular processes, but results often display force field dependency. Therefore, assessment of force field accuracy using experimental data of biomolecules in solution is essential for simulation studies. Here, we propose the use of structural models obtained via cryo-electron microscopy (cryoEM), which provides biomolecular structures in vitreous ice mimicking the environment in solution. The accuracy of the AMBER (ff99SB-ILDN-NMR, ff14SB, ff15ipq, and ff15FB) and CHARMM (CHARMM22 and CHARMM36m) force fields was assessed by comparing their MD trajectories with the cryoEM data of thermostable hexameric glutamate dehydrogenase (GDH), which included a cryoEM map at a resolution of approximately 3 Å and structure models of subunits reflecting metastable conformations in domain motion occurring in GDH. In the assessment, we validated the force fields with respect to the reproducibility and stability of secondary structures and intersubunit interactions in the cryoEM data. Furthermore, we evaluated the force fields regarding the reproducibility of the energy landscape in the domain motion expected from the cryoEM data. As a result, among the six force fields, ff15FB and ff99SB-ILDN-NMR displayed good agreement with the experiment. The present study demonstrated the advantages of the high-resolution cryoEM map and suggested the optimal force field to reproduce experimentally observed protein structures.

Methods and application of coherent X-ray diffraction imaging of noncrystalline particles.

Microscopic imaging techniques have been developed to visualize events occurring in biological cells. Coherent X-ray diffraction imaging is one of the techniques applicable to structural analyses of cells and organelles, which have never been crystallized. In the experiment, a single noncrystalline particle is illuminated by an X-ray beam with almost complete spatial coherence. The structure of the particle projected along the direction of the beam is, in principle, retrieved from a finely recorded diffraction pattern alone by using iterative phase-retrieval algorithms. Here, we describe fundamental theory and experimental methods of coherent X-ray diffraction imaging and the recent application in structural studies of noncrystalline specimens by using X-rays available at Super Photon Ring of 8-Gev and SPring-8 Angstrom Compact Free Electron Laser in Japan.

Energy landscape of domain motion in glutamate dehydrogenase deduced from cryo-electron microscopy.

Analysis of the conformational changes of protein is important to elucidate the mechanisms of protein motions correlating with their function. Here, we studied the spontaneous domain motion of unliganded glutamate dehydrogenase from Thermococcus profundus using cryo-electron microscopy and proposed a novel method to construct free-energy landscape of protein conformations. Each subunit of the homo-hexameric enzyme comprises nucleotide-binding domain (NAD domain) and hexamer-forming core domain. A large active-site cleft is situated between the two domains and varies from open to close according to the motion of a NAD domain. A three-dimensional map reconstructed from all cryo-electron microscopy images displayed disordered volumes of NAD domains, suggesting that NAD domains in the collected images adopted various conformations in domain motion. Focused classifications on NAD domain of subunits provided several maps of possible conformations in domain motion. To deduce what kinds of conformations appeared in EM images, we developed a novel analysis method that describe the EM maps as a linear combination of representative conformations appearing in a 200-ns molecular dynamics simulation as reference. The analysis enabled us to estimate the appearance frequencies of the representative conformations, which illustrated a free-energy landscape in domain motion. In the open/close domain motion, two free-energy basins hindered the direct transformation from open to closed state. Structure models constructed for representative EM maps in classifications demonstrated the correlation between the energy landscape and conformations in domain motion. Based on the results, the domain motion in glutamate dehydrogenase and the analysis method to visualize conformational changes and free-energy landscape were discussed. DATABASE: The EM maps of the four conformations were deposited to Electron Microscopy Data Bank (EMDB) as accession codes EMD-9845 (open), EMD-9846 (half-open1), EMD-9847 (half-open2), and EMD-9848 (closed), respectively. In addition, the structural models built for the four conformations were deposited to the Protein Data Bank (PDB) as accession codes 6JN9 (open), 6JNA (half-open1), 6JNC (half-open2), and 6JND (closed), respectively.

Dipole-dipole interactions between tryptophan side chains and hydration water molecules dominate the observed dynamic stokes shift of lysozyme.

The fluorescence intensity of tryptophan residues in hen egg-white lysozyme was measured up to 500 ps after the excitation by irradiation pulses at 290 nm. From the time-dependent variation of fluorescence intensity in a wavelength range of 320-370 nm, the energy relaxation in the dynamic Stokes shift was reconstructed as the temporal variation in wavenumber of the estimated fluorescence maximum. The relaxation was approximated by two exponential curves with decay constants of 1.2 and 26.7 ps. To interpret the relaxation, a molecular dynamics simulation of 75 ns was conducted for lysozyme immersed in a water box. From the simulation, the energy relaxation in the electrostatic interactions of each tryptophan residue was evaluated by using a scheme derived from the linear response theory. Dipole-dipole interactions between each of the Trp62 and Trp123 residues and hydration water molecules displayed an energy relaxation similar to that experimentally observed regarding time constants and magnitudes. The side chains of these residues were partly or fully exposed to the solvent. In addition, by inspecting the variation in dipole moments of the hydration water molecules around lysozyme, it was suggested that the observed relaxation could be attributed to the orientational relaxation of hydration water molecules participating in the hydrogen-bond network formed around each of the two tryptophan residues.

Molecular shape under far-red light and red light-induced association of Arabidopsis phytochrome B.

Phytochrome B (phyB) is a plant photoreceptor protein that regulates various photomorphogenic responses to optimize plant growth and development. PhyB exists in two photoconvertible forms: a red light-absorbing (Pr) and a far-red light-absorbing (Pfr) form. Therefore, to understand the mechanism of phototransformation, the structural characterization of full-length phyB in these two forms is necessary. Here, we report the molecular structure of Arabidopsis thaliana phyB in Pr form and the molecular properties of the Pfr form determined by small-angle X-ray scattering coupled with size-exclusion chromatography. In solution, the Pr form associated as a dimer with a radius of gyration of 50 Å. The molecular shape was a crossed shape, in which the orientation of the photosensory modules differed from that in the crystal structure of dimeric photosensory module. PhyB exhibited structural reversibility in the Pfr-to-Pr phototransformation and thermal reversion from Pfr to Pr in the dark. In addition, Pfr only exhibited nonspecific association, which distinguished molecular properties of Pfr form from those of the inactive Pr form.

Diffraction apparatus and procedure in tomography X-ray diffraction imaging for biological cells at cryogenic temperature using synchrotron X-ray radiation.

X-ray diffraction imaging is a technique for visualizing the structure of biological cells. In X-ray diffraction imaging experiments using synchrotron radiation, cryogenic conditions are necessary in order to reduce radiation damage in the biological cells. Frozen-hydrated biological specimens kept at cryogenic temperatures are also free from drying and bubbling, which occurs in wet specimens under vacuum conditions. In a previous study, the diffraction apparatus KOTOBUKI-1 [Nakasako et al. (2013), Rev. Sci. Instrum. 84, 093705] was constructed for X-ray diffraction imaging at cryogenic temperatures by utilizing a cryogenic pot, which is a cooling device developed in low-temperature physics. In this study a new cryogenic pot, suitable for tomography experiments, has been developed. The pot can rotate a biological cell over an angular range of ±170° against the direction of the incident X-ray beam. Herein, the details and the performance of the pot and miscellaneous devices are reported, along with established experimental procedures including specimen preparation. The apparatus has been used in tomography experiments for visualizing the three-dimensional structure of a Cyanidioschyzon merolae cell with an approximate size of 5 µm at a resolution of 136 nm. Based on the experimental results, the necessary improvements for future experiments and the resolution limit achievable under experimental conditions within a maximum tolerable dose are discussed.

Classification of ab initio models of proteins restored from small-angle X-ray scattering.

In structure analyses of proteins in solution by using small-angle X-ray scattering (SAXS), the molecular models are restored by using ab initio molecular modeling algorithms. There can be variation among restored models owing to the loss of phase information in the scattering profiles, averaging with regard to the orientation of proteins against the direction of the incident X-ray beam, and also conformational fluctuations. In many cases, a representative molecular model is obtained by averaging models restored in a number of ab initio calculations, which possibly provide nonrealistic models inconsistent with the biological and structural information about the target protein. Here, a protocol for classifying predicted models by multivariate analysis to select probable and realistic models is proposed. In the protocol, each structure model is represented as a point in a hyper-dimensional space describing the shape of the model. Principal component analysis followed by the clustering method is applied to visualize the distribution of the points in the hyper-dimensional space. Then, the classification provides an opportunity to exclude nonrealistic models. The feasibility of the protocol was examined through the application to the SAXS profiles of four proteins.

Growth of Cuprous Oxide Particles in Liquid-Phase Synthesis Investigated by X-ray Laser Diffraction.

Cuprous oxide (Cu2O) particles obtained by surfactant-assisted liquid-phase synthesis have cuboid shapes but the internal structures are difficult to be visualized by electron microscopy. Herein, we investigated the internal structures of numerous individual Cu2O particles with submicrometer dimensions by X-ray diffraction imaging (XDI) using X-ray free-electron laser (XFEL) pulses. The reconstructed two-dimensional electron density maps, which displayed inhomogeneous internal structures, were divided into five classes characterized by the positions and shapes of high and low electron density areas. Further analysis of the maps in each class by a manifold learning algorithm revealed that the internal structures of Cu2O particles varied in correlation with total electron density while retaining the characteristics within each class. On the basis of the analyses, we proposed a growth mechanism to yield the inhomogeneity in the internal structures of Cu2O particles in surfactant-mediated liquid-phase synthesis.

Shot-by-shot characterization of focused X-ray free electron laser pulses.

X-ray free electron lasers (XFEL) provide intense and almost coherent X-ray pulses. They are used for various experiments investigating physical and chemical properties in materials and biological science because of their complete coherence, high intensity, and very short pulse width. In XFEL experiments, specimens are irradiated by XFEL pulses focused by mirror optics. The focused pulse is too intense to measure its coherence by placing an X-ray detector on the focal spot. Previously, a method was proposed for evaluating the coherence of focused pulses from the visibility of the diffraction intensity of colloidal particles by the speckle visibility spectroscopy (SVS). However, the visibility cannot be determined exactly because the diffraction intensity is integrated into each finite size detector pixel. Here, we propose a method to evaluate the coherence of each XFEL pulse by using SVS in combination with a theory for exact sampling of the diffraction pattern and a technique of multiplying the diffraction data by a Gaussian masks, which reduces the influence of data missing in small-angle regions due to the presence of a direct beamstop. We also introduce a method for characterizing the shot-by-shot size of each XFEL pulse by analysing the X-ray irradiated area.

Blue light-excited LOV1 and LOV2 domains cooperatively regulate the kinase activity of full-length phototropin2 from Arabidopsis.

Phototropin2 (phot2) is a blue-light (BL) receptor that regulates BL-dependent activities for efficient photosynthesis in plants. phot2 comprises two BL-receiving light-oxygen-voltage-sensing domains (LOV1 and LOV2) and a kinase domain. BL-excited LOV2 is thought to be primarily responsible for the BL-dependent activation of the kinase. However, the molecular mechanisms by which small BL-induced conformational changes in the LOV2 domain are transmitted to the kinase remain unclear. Here, we used full-length wild-type and mutant phot2 proteins from Arabidopsis to study their molecular properties in the dark and under BL irradiation. Phosphorylation assays and absorption measurements indicated that the LOV1 domain assists the thermal relaxation of BL-excited LOV2 and vice versa. Using small-angle X-ray scattering and electron microscopy, we observed that phot2 forms a dimer and has a rod shape with a maximum length of 188 Å and a radius of gyration of 44 Å. Under BL, phot2 displayed large conformational changes that bent the rod shape. By superimposing the crystal structures of the LOV1 dimer, LOV2, and a homology model of the kinase to the observed changes, we inferred that the BL-dependent change consisted of positional shifts of both LOV2 and the kinase relative to LOV1. Furthermore, phot2 mutants lacking the photocycle in LOV1 or LOV2 still exhibited conformational changes under BL, suggesting that LOV1 and LOV2 cooperatively contribute to the conformational changes that activate the kinase. These results suggest that BL-activated LOV1 contributes to the kinase activity of phot2. We discuss the possible intramolecular interactions and signaling mechanisms in phot2.

Influences of lone-pair electrons on directionality of hydrogen bonds formed by hydrophilic amino acid side chains in molecular dynamics simulation.

The influence of lone-pair electrons on the directionality of hydrogen bonds that are formed by oxygen and nitrogen atoms in the side chains of nine hydrophilic was investigated using molecular dynamics simulations. The simulations were conducted using two types of force fields; one incorporated lone-pair electrons placed at off-atom sites and the other did not. The density distributions of the hydration water molecules around the oxygen and nitrogen atoms were calculated from the simulation trajectories, and were compared with the empirical hydration distribution functions, which were constructed from a large number of hydration water molecules found in the crystal structures of proteins. Only simulations using the force field explicitly incorporating lone-pair electrons reproduced the directionality of hydrogen bonds that is observed in the empirical distribution functions for the deprotonated oxygen and nitrogen atoms in the sp 2-hybridization. The amino acids that include such atoms are functionally important glutamate, aspartate, and histidine. Therefore, a set of force field that incorporates lone-pair electrons as off-atom charge sites would be effective for considering hydrogen bond formation by these amino acids in molecular dynamics simulation studies.