High-resolution structure of phosphoketolase from Bifidobacterium longum determined by cryo-EM single-particle analysis.

In bifidobacteria, phosphoketolase (PKT) plays a key role in the central hexose fermentation pathway called "bifid shunt." The three-dimensional structure of PKT from Bifidobacterium longum with co-enzyme thiamine diphosphate (ThDpp) was determined at 2.1 Å resolution by cryo-EM single-particle analysis using 196,147 particles to build up the structural model of a PKT octamer related by D4 symmetry. Although the cryo-EM structure of PKT was almost identical to the X-ray crystal structure previously determined at 2.2 Å resolution, several interesting structural features were observed in the cryo-EM structure. Because this structure was solved at relatively high resolution, it was observed that several amino acid residues adopt multiple conformations. Among them, Q546-D547-H548-N549 (the QN-loop) demonstrate the largest structural change, which seems to be related to the enzymatic function of PKT. The QN-loop is at the entrance to the substrate binding pocket. The minor conformer of the QN-loop is similar to the conformation of the QN-loop in the crystal structure. The major conformer is located further from ThDpp than the minor conformer. Interestingly, the major conformer in the cryo-EM structure of PKT resembles the corresponding loop structure of substrate-bound Escherichia coli transketolase. That is, the minor and major conformers may correspond to "closed" and "open" states for substrate access, respectively. Moreover, because of the high-resolution analysis, many water molecules were observed in the cryo-EM structure of PKT. Structural features of the water molecules in the cryo-EM structure are discussed and compared with water molecules observed in the crystal structure.

Protocol for Retrieving Three-Dimensional Biological Shapes for a Few XFEL Single-Particle Diffraction Patterns.

X-ray free-electron laser (XFEL) scattering promises to probe single biomolecular complexes without crystallization, enabling the study of biomolecular structures under near-physiological conditions at room temperature. However, such structural determination of biomolecules is extremely challenging thus far. In addition to the large numbers of diffraction patterns required, the orientation of each diffraction pattern needs to be accurately estimated and the missing phase information needs to be recovered for three-dimensional (3D) structure reconstruction. Given the current limitations to the amount and resolution of the data available from single-particle XFEL scattering experiments, we propose an alternative approach to find plausible 3D biological shapes from a limited number of diffraction patterns to serve as a starting point for further analyses. In our proposed strategy, small sets of input (e.g., five) XFEL diffraction patterns were matched against a library of diffraction patterns simulated from 1628 electron microscopy (EM) models to find potential matching 3D models that are consistent with the input diffraction patterns. This approach was tested for three example cases: EMD-3457 (Thermoplasma acidophilum 20S proteasome), EMD-5141 (Escherichia coli 70S ribosome complex), and EMD-5152 (budding yeast Nup84 complex). We observed that choosing the best strategy to define matching regions on the diffraction patterns is critical for identifying correctly matching diffraction patterns. While increasing the number of input diffraction patterns improved the matches in some cases, we found that the resulting matches are more dependent on the uniqueness or complexity of the shape as captured in the individual input diffraction patterns and the availability of a similar 3D biological shape in the search library. The protocol could be useful for finding candidate models for a limited amount of low-resolution data, even when insufficient for reconstruction, performing a quick exploration of new data upon collection, and the analysis of the conformational heterogeneity of the particle of interest as captured within the diffraction patterns.

Computational Protocol for Assessing the Optimal Pixel Size to Improve the Accuracy of Single-particle Cryo-electron Microscopy Maps.

Cryo-electron microscopy (cryo-EM) single-particle analysis has come a long way in achieving atomic-level resolution when imaging biomolecules. To obtain the best possible three-dimensional (3D) structure in cryo-EM, many parameters have to be carefully considered. Here we address the often-overlooked parameter of the pixel size, which describes the magnification of the image produced by the experiment. While efforts are made to refine and validate this parameter in the analysis of cryo-EM experimental data, there is no systematic protocol in place. Since the pixel size parameter can have an impact on the resolution and accuracy of a cryo-EM map, and the atomic resolution 3D structure models derived from it, we propose a computational protocol to estimate the appropriate pixel size parameter. In our protocol, we fit and refine atomic structures against cryo-EM maps at multiple pixel sizes. The resulting fitted and refined structures are evaluated using the GOAP (generalized orientation-dependent, all-atom statistical potential) score, which we found to perform better than other commonly used functions, such as Molprobity and the correlation coefficient from refinement. Finally, we describe the efficacy of this protocol in retrieving appropriate pixel sizes for several examples; simulated data based on yeast elongation factor 2 and experimental data from Gro-EL chaperone, beta-galactosidase, and the TRPV1 ion channel.

Cryo-Cooling Effect on DHFR Crystal Studied by Replica-Exchange Molecular Dynamics Simulations.

Cryo-cooling is routinely performed before x-ray diffraction image collection to reduce the damage to crystals due to ionizing radiation. It has been suggested that although backbone structures are usually very similar between room temperature and cryo-temperature, cryo-cooling may hamper biologically relevant dynamics. In this study, the crystal of Escherichia coli dihydrofolate reductase is studied with replica-exchange molecular dynamics simulation, and the results are compared with the crystal structure determined at cryo-temperature and room temperature with the time-averaged ensemble method. Although temperature dependence of unit cell compaction and root mean-square fluctuation of Cα is found in accord with experiment, it is found that the protein structure at low temperature can be more heterogeneous than the ensemble of structures reported by using the time-averaged ensemble method, encouraging further development of the time-averaged ensemble method and indicating that data should be examined carefully to avoid overinterpretation of one average structure.

Computational investigation of the conformational dynamics in Tom20-mitochondrial presequence tethered complexes.

The translocase of the outer membrane (TOM) mediates the membrane permeation of mitochondrial matrix proteins. Tom20 is a subunit of the TOM complex and binds to the N-terminal region (ie, presequence) in mitochondrial matrix precursor proteins. Previous experimental studies indicated that the presequence recognition by Tom20 was achieved in a dynamic-equilibrium among multiple bound states of the α-helical presequence. Accordingly, the co-crystallization of Tom20 and a presequence peptide required a disulfide-bond cross-linking. A 3-residue spacer sequence (XAG) was inserted between the presequence and the anchoring Cys residue at the C-terminus to not disturb the movement of the presequence peptide in the binding site of Tom20. Two crystalline forms were obtained according to Ala or Tyr at the X position of the spacer sequence, which may reflect the dynamic-equilibrium of the presequence. Here, we have performed replica-exchange molecular dynamics (REMD) simulations to study the effect of disulfide-bond linker and single amino acid difference in the spacer region of the linker on the conformational dynamics of Tom20-presequence complex. Free energy and network analyses of the REMD simulations were compared against previous simulations of non-tethered system. We concluded that the disulfide-bond tethering did not strongly affect the conformational ensemble of the presequence peptide in the complex. Further investigation showed that the choice of Ala or Tyr at the X position did not affect the most distributions of the conformational ensemble of the presequence. The present study provides a rational basis for the disulfide-bond tethering to study the dynamics of weakly binding complexes.

Searching for 3D structural models from a library of biological shapes using a few 2D experimental images.

BACKGROUND: Advancements in biophysical experimental techniques have pushed the limits in terms of the types of phenomena that can be characterized, the amount of data that can be produced and the resolution at which we can visualize them. Single particle techniques such as Electron Microscopy (EM) and X-ray free electron laser (XFEL) scattering require a large number of 2D images collected to resolve three-dimensional (3D) structures. In this study, we propose a quick strategy to retrieve potential 3D shapes, as low-resolution models, from a few 2D experimental images by searching a library of 2D projection images generated from existing 3D structures. RESULTS: We developed the protocol to assemble a non-redundant set of 3D shapes for generating the 2D image library, and to retrieve potential match 3D shapes for query images, using EM data as a test. In our strategy, we disregard differences in volume size, giving previously unknown structures and conformations a greater number of 3D biological shapes as possible matches. We tested the strategy using images from three EM models as query images for searches against a library of 22750 2D projection images generated from 250 random EM models. We found that our ability to identify 3D shapes that match the query images depends on how complex the outline of the 2D shapes are and whether they are represented in the search image library. CONCLUSIONS: Through our computational method, we are able to quickly retrieve a 3D shape from a few 2D projection images. Our approach has the potential for exploring other types of 2D single particle structural data such as from XFEL scattering experiments, for providing a tool to interpret low-resolution data that may be insufficient for 3D reconstruction, and for estimating the mixing of states or conformations that could exist in such experimental data.

Single-particle XFEL 3D reconstruction of ribosome-size particles based on Fourier slice matching: requirements to reach subnanometer resolution.

Three-dimensional (3D) structures of biomolecules provide insight into their functions. Using X-ray free-electron laser (XFEL) scattering experiments, it was possible to observe biomolecules that are difficult to crystallize, under conditions that are similar to their natural environment. However, resolving 3D structure from XFEL data is not without its challenges. For example, strong beam intensity is required to obtain sufficient diffraction signal and the beam incidence angles to the molecule need to be estimated for diffraction patterns with significant noise. Therefore, it is important to quantitatively assess how the experimental conditions such as the amount of data and their quality affect the expected resolution of the resulting 3D models. In this study, as an example, the restoration of 3D structure of ribosome from two-dimensional diffraction patterns created by simulation is shown. Tests are performed using the diffraction patterns simulated for different beam intensities and using different numbers of these patterns. Guidelines for selecting parameters for slice-matching 3D reconstruction procedures are established. Also, the minimum requirements for XFEL experimental conditions to obtain diffraction patterns for reconstructing molecular structures to a high-resolution of a few nanometers are discussed.

Gaussian mixture model for coarse-grained modeling from XFEL.

We explore the advantage of Gaussian mixture model (GMM) for interpretation of single particle diffraction patterns from X-ray free electron laser (XFEL) experiments. GMM approximates a biomolecular shape by the superposition of Gaussian distributions. As the Fourier transformation of GMM can be quickly performed, we can efficiently simulate XFEL diffraction patterns from approximated structure models. We report that the resolution that GMM can accurately reproduce is proportional to the cubic root of the number of Gaussians used in the modeling. This behavior can be attributed to the correspondence between the number of adjustable parameters in GMM and the amount of sampling points in diffraction space. Furthermore, GMMs can successfully be used to perform angular assignment and to detect conformational variation. These results demonstrate that GMMs serve as useful coarse-grained models for hybrid approach in XFEL single particle experiments.

Hybrid Methods for Macromolecular Modeling by Molecular Mechanics Simulations with Experimental Data.

Hybrid approaches for the modeling of macromolecular complexes that combine computational molecular mechanics simulations with experimental data are discussed. Experimental data for biological molecular structures are often low-resolution, and thus, do not contain enough information to determine the atomic positions of molecules. This is especially true when the dynamics of large macromolecules are the focus of the study. However, computational modeling can complement missing information. Significant increase in computational power, as well as the development of new modeling algorithms allow us to model structures of biological macromolecules reliably, using experimental data as references. We review the basics of molecular mechanics approaches, such as atomic model force field, and coarse-grained models, molecular dynamics simulation and normal mode analysis and describe how they could be used for flexible fitting hybrid modeling with experimental data, especially from cryo-EM and SAXS.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics.

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.

Flexible fitting to cryo-EM density map using ensemble molecular dynamics simulations.

Flexible fitting is a computational algorithm to derive a new conformational model that conforms to low-resolution experimental data by transforming a known structure. A common application is against data from cryo-electron microscopy to obtain conformational models in new functional states. The conventional flexible fitting algorithms cannot derive correct structures in some cases due to the complexity of conformational transitions. In this study, we show the importance of conformational ensemble in the refinement process by performing multiple fittings trials using a variety of different force constants. Application to simulated maps of Ca2+ ATPase and diphtheria toxin as well as experimental data of release factor 2 revealed that for these systems, multiple conformations with similar agreement with the density map exist and a large number of fitting trials are necessary to generate good models. Clustering analysis can be an effective approach to avoid over-fitting models. In addition, we show that an automatic adjustment of the biasing force constants during the fitting process, implemented as replica-exchange scheme, can improve the success rate. © 2017 Wiley Periodicals, Inc.

Three-dimensional reconstruction for coherent diffraction patterns obtained by XFEL.

The three-dimensional (3D) structural analysis of single particles using an X-ray free-electron laser (XFEL) is a new structural biology technique that enables observations of molecules that are difficult to crystallize, such as flexible biomolecular complexes and living tissue in the state close to physiological conditions. In order to restore the 3D structure from the diffraction patterns obtained by the XFEL, computational algorithms are necessary as the orientation of the incident beam with respect to the sample needs to be estimated. A program package for XFEL single-particle analysis based on the Xmipp software package, that is commonly used for image processing in 3D cryo-electron microscopy, has been developed. The reconstruction program has been tested using diffraction patterns of an aerosol nanoparticle obtained by tomographic coherent X-ray diffraction microscopy.

Thermodynamic properties of water molecules in the presence of cosolute depend on DNA structure: a study using grid inhomogeneous solvation theory.

In conditions that mimic those of the living cell, where various biomolecules and other components are present, DNA strands can adopt many structures in addition to the canonical B-form duplex. Previous studies in the presence of cosolutes that induce molecular crowding showed that thermal stabilities of DNA structures are associated with the properties of the water molecules around the DNAs. To understand how cosolutes, such as ethylene glycol, affect the thermal stability of DNA structures, we investigated the thermodynamic properties of water molecules around a hairpin duplex and a G-quadruplex using grid inhomogeneous solvation theory (GIST) with or without cosolutes. Our analysis indicated that (i) cosolutes increased the free energy of water molecules around DNA by disrupting water-water interactions, (ii) ethylene glycol more effectively disrupted water-water interactions around Watson-Crick base pairs than those around G-quartets or non-paired bases, (iii) due to the negative electrostatic potential there was a thicker hydration shell around G-quartets than around Watson-Crick-paired bases. Our findings suggest that the thermal stability of the hydration shell around DNAs is one factor that affects the thermal stabilities of DNA structures under the crowding conditions.

Hybrid approach for structural modeling of biological systems from X-ray free electron laser diffraction patterns.

We present a new hybrid approach for structural modeling using X-ray free electron laser (XFEL) diffraction patterns from non-crystalline biological samples. Reconstruction of a 3D structure requires a large number of diffraction patterns; however, in the current XFEL experiments with biological systems, the analysis often relies on a small number of 2D diffraction patterns. In this study, we explore the strategies to identify plausible 3D structural models by combining the 2D analysis of such diffraction patterns with computational modeling (normal mode analysis or molecular dynamics simulations). As the first step toward such hybrid modeling, we established a protocol to assess the agreement between the model structure and the target XFEL diffraction pattern and showed that XFEL data can be used to study the conformational transitions of biological molecules. We tested the proposed algorithms using data of three biomolecular complexes of different sizes (elongation factor 2, CCM virus, and ribosome) and examined the experimental conditions that are required to perform such studies, in particular the XFEL beam intensity requirements. The results indicate that the current beam intensity is close to a strength that enables us to study conformational transitions of macromolecules, such as ribosomes. The proposed algorithm can be combined with molecular mechanics approaches, such as molecular dynamics simulations and normal mode analysis, to generate a large number of candidate structures to perform hybrid structural modeling.

Packing interface energetics in different crystal forms of the λ Cro dimer.

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site-directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson-Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ∼5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures.

Structure determination using high-order spatial correlations in single-particle X-ray scattering.

Single-particle imaging using X-ray free-electron lasers (XFELs) is a promising technique for observing nanoscale biological samples under near-physiological conditions. However, as the sample's orientation in each diffraction pattern is unknown, advanced algorithms are required to reconstruct the 3D diffraction intensity volume and subsequently the sample's density model. While most approaches perform 3D reconstruction via determining the orientation of each diffraction pattern, a correlation-based approach utilizes the averaged spatial correlations of diffraction intensities over all patterns, making it well suited for processing experimental data with a poor signal-to-noise ratio of individual patterns. Here, a method is proposed to determine the 3D structure of a sample by analyzing the double, triple and quadruple spatial correlations in diffraction patterns. This ab initio method can reconstruct the basic shape of an irregular unsymmetric 3D sample without requiring any prior knowledge of the sample. The impact of background and noise on correlations is investigated and corrected to ensure the success of reconstruction under simulated experimental conditions. Additionally, the feasibility of using the correlation-based approach to process incomplete partial diffraction patterns is demonstrated. The proposed method is a variable addition to existing algorithms for 3D reconstruction and will further promote the development and adoption of XFEL single-particle imaging techniques.

Increased redox-sensitive green fluorescent protein reduction potential in the endoplasmic reticulum following glutathione-mediated dimerization.

As the endoplasmic reticulum (ER) is the compartment where disulfide bridges in secreted and cell surface proteins are formed, the disturbance of its redox state has profound consequences, yet regulation of ER redox potential remains poorly understood. To monitor the ER redox state in live cells, several fluorescence-based sensors have been developed. However, these sensors have yielded results that are inconsistent with each other and with earlier non-fluorescence-based studies. One particular green fluorescent protein (GFP)-based redox sensor, roGFP1-iL, could detect oxidizing changes in the ER despite having a reduction potential significantly lower than that previously reported for the ER. We have confirmed these observations and determined the mechanisms by which roGFP1-iL detects oxidizing changes. First, glutathione mediates the formation of disulfide-bonded roGFP1-iL dimers with an intermediate excitation fluorescence spectrum resembling a mixture of oxidized and reduced monomers. Second, glutathione facilitates dimerization of roGFP1-iL, which shifted the equilibrium from oxidized monomers to dimers, thereby increasing the molecule's reduction potential compared with that of a dithiol redox buffer. We conclude that the glutathione redox couple in the ER significantly increased the reduction potential of roGFP1-iL in vivo by facilitating its dimerization while preserving its ratiometric nature, which makes it suitable for monitoring oxidizing and reducing changes in the ER with a high degree of reliability in real time. The ability of roGFP1-iL to detect both oxidizing and reducing changes in ER and its dynamic response in glutathione redox buffer between approximately -190 and -130 mV in vitro suggests a range of ER redox potentials consistent with those determined by earlier approaches that did not involve fluorescent sensors.

Advancing cryo-electron microscopy data analysis through accelerated simulation-based flexible fitting approaches.

Flexible fitting based on molecular dynamics simulation is a technique for structure modeling from cryo-EM data. It has been utilized for nearly two decades, and while cryo-EM resolution has improved significantly, it remains a powerful approach that can provide structural and dynamical insights that are not directly accessible from experimental data alone. Molecular dynamics simulations provide a means to extract atomistic details of conformational changes that are encoded in cryo-EM data and can also assist in improving the quality of structural models. Additionally, molecular dynamics simulations enable the characterization of conformational heterogeneity in cryo-EM data. We will summarize the advancements made in these techniques and highlight recent developments in this field.

Molecular size dependence on achievable resolution from XFEL single-particle 3D reconstruction.

Single-particle analysis using x-ray free-electron lasers (XFELs) is a novel method for obtaining structural information of samples in a state close to nature. In particular, it is suitable for observing the inner structure of large biomolecules by taking advantage of the high transmittance of x-rays. However, systematic studies on the resolution achievable for large molecules are lacking. In this study, the molecular size dependence of the resolution of a three-dimensional (3D) structure resulting from XFEL single-particle reconstruction is evaluated using synthetic data. Evidently, 3D structures of larger molecules can be restored with higher detail (defined relative to the molecular sizes) than smaller ones; however, reconstruction with high absolute resolution (defined in nm-1) is challenging. Our results provide useful information for the experimental design of 3D structure reconstruction using coherent x-ray diffraction patterns of single-particles.

MDSPACE: Extracting Continuous Conformational Landscapes from Cryo-EM Single Particle Datasets Using 3D-to-2D Flexible Fitting based on Molecular Dynamics Simulation.

This article presents an original approach for extracting atomic-resolution landscapes of continuous conformational variability of biomolecular complexes from cryo electron microscopy (cryo-EM) single particle images. This approach is based on a new 3D-to-2D flexible fitting method, which uses molecular dynamics (MD) simulation and is embedded in an iterative conformational-landscape refinement scheme. This new approach is referred to as MDSPACE, which stands for Molecular Dynamics simulation for Single Particle Analysis of Continuous Conformational hEterogeneity. The article describes the MDSPACE approach and shows its performance using synthetic and experimental datasets.