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

Single-particle Cryo-EM and molecular dynamics simulations: A perfect match

Knowledge of the structure and dynamics of biomolecules is key to understanding the mechanisms underlying their biological functions. Single-particle cryo-electron microscopy (cryo-EM) is a powerful structural biology technique to characterize complex biomolecular systems. Here, we review recent advances of how Molecular Dynamics (MD) simulations are being used to increase and enhance the information extracted from cryo-EM experiments. We will particularly focus on the physics underlying these experiments, how MD facilitates structure refinement, in particular for heterogeneous and non-isotropic resolution, and how thermodynamic and kinetic information can be extracted from cryo-EM data.

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.

Four-dimensional imaging for cryo-electron microscopy experiments using molecular simulations and manifold learning

Elucidating protein conformational changes is essential because conformational changes are closely related to the functions of proteins. Cryo-electron microscopy (cryo-EM) experiment can be used to reconstruct protein conformational changes via a method that involves using the experimental data (two-dimensional protein images). In this study, a reconstruction method, referred to as the "four-dimensional imaging," was proposed. In our four-dimensional imaging technique, the protein conformational change was obtained using the two-dimensional protein images (the three-dimensional electron density maps used in previously proposed techniques were not used). The protein conformation for each two-dimensional protein image was obtained using our original protocol with molecular dynamics simulations. Using a manifold-learning technique and two-dimensional protein images, the protein conformations were arranged according to the conformational change of the protein. By arranging the protein conformations according to the arrangement of the protein images, four-dimensional imaging is constructed. A simulation for a cryo-EM experiment demonstrated the validity of our four-dimensional imaging technique.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces phosphatidylserine externalization, a key step in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements underlying groove opening and how lipids reorganize outside the closed groove remain unknown. Here we directly visualize how lipids associate at the closed groove of Ca2+-bound fungal nhTMEM16 in nanodiscs using cryo-EM. Functional experiments pinpoint lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryo-EM structure determination.

Structural basis of closed groove scrambling by a TMEM16 protein

Activation of Ca2+-dependent TMEM16 scramblases induces the externalization of phosphatidylserine, a key molecule in multiple signaling processes. Current models suggest that the TMEM16s scramble lipids by deforming the membrane near a hydrophilic groove, and that Ca2+ dependence arises from the different association of lipids with an open or closed groove. However, the molecular rearrangements involved in groove opening and of how lipids reorganize outside the closed groove remain unknown. Using cryogenic electron microscopy, we directly visualize how lipids associate at the closed groove of Ca2+-bound nhTMEM16 in nanodiscs. Functional experiments pinpoint the lipid-protein interaction sites critical for closed groove scrambling. Structural and functional analyses suggest groove opening entails the sequential appearance of two π-helical turns in the groove-lining TM6 helix and identify critical rearrangements. Finally, we show that the choice of scaffold protein and lipids affects the conformations of nhTMEM16 and their distribution, highlighting a key role of these factors in cryoEM structure determination.

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.

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.

Elucidating polymorphs of crystal structures by intensity-based hierarchical clustering analysis of multiple diffraction data sets

In macromolecular structure determination using X-ray diffraction from multiple crystals, the presence of different structures (structural polymorphs) necessitates the classification of the diffraction data for appropriate structural analysis. Hierarchical clustering analysis (HCA) is a promising technique that has so far been used to extract isomorphous data, mainly for single-structure determination. Although in principle the use of HCA can be extended to detect polymorphs, the absence of a reference to define the threshold used to group the isomorphous data sets (the `isomorphic threshold') poses a challenge. Here, unit-cell-based and intensity-based HCAs have been applied to data sets for apo trypsin and inhibitor-bound trypsin that were mixed post data acquisition to investigate the efficacy of HCA in classifying polymorphous data sets. Single-step intensity-based HCA successfully classified polymorphs with a certain `isomorphic threshold'. In data sets for several samples containing an unknown degree of structural heterogeneity, polymorphs could be identified by intensity-based HCA using the suggested `isomorphic threshold'. Polymorphs were also detected in single crystals using data collected using the continuous helical scheme. These findings are expected to facilitate the determination of multiple structural snapshots by exploiting automated data collection and analysis.

Conformational heterogeneity and probability distributions from single-particle cryo-electron microscopy

Single-particle cryo-electron microscopy (cryo-EM) is a technique that takes projection images of biomolecules frozen at cryogenic temperatures. A major advantage of this technique is its ability to image single biomolecules in heterogeneous conformations. While this poses a challenge for data analysis, recent algorithmic advances have enabled the recovery of heterogeneous conformations from the noisy imaging data. Here, we review methods for the reconstruction and heterogeneity analysis of cryo-EM images, ranging from linear-transformation-based methods to nonlinear deep generative models. We overview the dimensionality-reduction techniques used in heterogeneous 3D reconstruction methods and specify what information each method can infer from the data. Then, we review the methods that use cryo-EM images to estimate probability distributions over conformations in reduced subspaces or predefined by atomistic simulations. We conclude with the ongoing challenges for the cryo-EM community.

Structural and functional studies of Arabidopsis thaliana glutamate dehydrogenase isoform 2 demonstrate enzyme dynamics and identify its calcium binding site

Glutamate dehydrogenase (GDH) is an enzyme at the crossroad of plant nitrogen and carbon metabolism. GDH catalyzes the conversion of 2-oxoglutarate into glutamate (2OG → Glu), utilizing ammonia as cosubstrate and NADH as coenzyme. The GDH reaction is reversible, meaning that the NAD+-dependent reaction (Glu → 2OG) releases ammonia. In Arabidopsis thaliana, three GDH isoforms exist, AtGDH1, AtGDH2, and AtGDH3. The subject of this work is AtGDH2. Previous reports have suggested that enzymes homologous to AtGDH2 contain a calcium-binding EF-hand motif located in the coenzyme binding domain. Here, we show that while AtGDH2 indeed does bind calcium, the binding occurs elsewhere and the region predicted to be the EF-hand motif has a completely different structure. As the true calcium binding site is > 20 Å away from the active site, it seems to play a structural, rather than catalytic role. We also performed comparative kinetic characterization of AtGDH1 and AtGDH2 using spectroscopic methods and isothermal titration calorimetry, to note that the isoenzymes generally exhibit similar behavior, with calcium having only a minor effect. However, the spatial and temporal changes in the gene expression profiles of the three AtGDH genes point to AtGDH2 as the most prevalent isoform.

NAD(P)H binding configurations revealed by time-resolved fluorescence and two-photon absorption

NADH and NADPH play key roles in the regulation of metabolism. Their endogenous fluorescence is sensitive to enzyme binding, allowing changes in cellular metabolic state to be determined using fluorescence lifetime imaging microscopy (FLIM). However, to fully uncover the underlying biochemistry, the relationships between their fluorescence and binding dynamics require greater understanding. Here we accomplish this through time- and polarization-resolved fluorescence and polarized two-photon absorption measurements. Two lifetimes result from binding of both NADH to lactate dehydrogenase and NADPH to isocitrate dehydrogenase. The composite fluorescence anisotropy indicates the shorter (1.3-1.6 ns) decay component to be accompanied by local motion of the nicotinamide ring, pointing to attachment solely via the adenine moiety. For the longer lifetime (3.2-4.4 ns), the nicotinamide conformational freedom is found to be fully restricted. As full and partial nicotinamide binding are recognized steps in dehydrogenase catalysis, our results unify photophysical, structural, and functional aspects of NADH and NADPH binding and clarify the biochemical processes that underlie their contrasting intracellular lifetimes.

Equilibrium Dynamics of a Biomolecular Complex Analyzed at Single-amino Acid Resolution by Cryo-electron Microscopy

The biological function of macromolecular complexes depends not only on large-scale transitions between conformations, but also on small-scale conformational fluctuations at equilibrium. Information on the equilibrium dynamics of biomolecular complexes could, in principle, be obtained from local resolution (LR) data in cryo-electron microscopy (cryo-EM) maps. However, this possibility had not been validated by comparing, for a same biomolecular complex, LR data with quantitative information on equilibrium dynamics obtained by an established solution technique. In this study we determined the cryo-EM structure of the minute virus of mice (MVM) capsid as a model biomolecular complex. The LR values obtained correlated with crystallographic B factors and with hydrogen/deuterium exchange (HDX) rates obtained by mass spectrometry (HDX-MS), a gold standard for determining equilibrium dynamics in solution. This result validated a LR-based cryo-EM approach to investigate, with high spatial resolution, the equilibrium dynamics of biomolecular complexes. As an application of this approach, we determined the cryo-EM structure of two mutant MVM capsids and compared their equilibrium dynamics with that of the wild-type MVM capsid. The results supported a previously suggested linkage between mechanical stiffening and impaired equilibrium dynamics of a virus particle. Cryo-EM is emerging as a powerful approach for simultaneously acquiring information on the atomic structure and local equilibrium dynamics of biomolecular complexes.

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.

Rational identification of a high catalytic efficiency leucine dehydrogenase and process development for efficient synthesis of l-phenylglycine

Enzymatic asymmetric synthesis of chiral amino acids has great industrial potential. However, the low catalytic efficiency of high-concentration substrates limits their industrial application. Herein, using a combination of substrate catalytic efficiency prediction based on "open to closed" conformational change and substrate specificity prediction, a novel leucine dehydrogenase (TsLeuDH), with high substrate catalytic efficiency toward benzoylformic acid (BFA) for producing l-phenylglycine (l-Phg), was directly identified from 4695 putative leucine dehydrogenases in a public database. The specific activity of TsLeuDH was determined to be as high as 4253.8 U mg^-1 . Through reaction process optimization, a high-concentration substrate (0.7 m) was efficiently and completely converted within 90 min in a single batch, without any external coenzyme addition. Moreover, a continuous flow-feeding approach was designed using gradient control of the feed rate to reduce substrate accumulation. Finally, the highest overall substrate concentration of up to 1.2 m BFA could be aminated to l-Phg with conversion of >99% in 3 h, demonstrating that this new combination of enzyme process development is promising for large-scale application of l-Phg.

Weight average approaches for predicting dynamical properties of biomolecules

Recent advances in atomistic molecular dynamics (MD) simulations of biomolecules allow us to explore their conformational spaces widely, observing large-scale conformational fluctuations or transitions between distinct structures. To reproduce or refine experimental data using MD simulations, structure ensembles, which are characterized by multiple structures and their statistical weights on the rugged free-energy landscapes, are often used. Here, we summarize weight average approaches for various experimental measurements. Weight average approaches are now applied to hybrid quantum mechanics/molecular mechanics MD simulations to predict fast vibrational motions in a protein with a high accuracy for better understanding of molecular functions from atomic structures.

Multiscale simulations of large complexes in conjunction with cryo-EM analysis

The cellular environment is highly crowded with most proteins and RNA/DNA forming homomeric and heteromeric complexes. Essential questions regarding how these complexes switch between functional, rest, and abnormal states with regulators or modifications remain challenging and complicated. Here, we review the recent progress integrating cryoelectron microscopy and multiscale molecular modeling to understand the dynamics and function-related mechanism in protein-RNA/DNA complexes, protein-protein complexes/assemblies, and membrane protein complexes. One future direction of multiscale simulations will be to interpret the large complex multibody regulation in assembly-induced function enhancement in conjunction with advanced atomic resolution structural-biology techniques and specialized computing architectures.

A Bayesian approach to extracting free-energy profiles from cryo-electron microscopy experiments

Cryo-electron microscopy (cryo-EM) extracts single-particle density projections of individual biomolecules. Although cryo-EM is widely used for 3D reconstruction, due to its single-particle nature it has the potential to provide information about a biomolecule's conformational variability and underlying free-energy landscape. However, treating cryo-EM as a single-molecule technique is challenging because of the low signal-to-noise ratio (SNR) in individual particles. In this work, we propose the cryo-BIFE method (cryo-EM Bayesian Inference of Free-Energy profiles), which uses a path collective variable to extract free-energy profiles and their uncertainties from cryo-EM images. We test the framework on several synthetic systems where the imaging parameters and conditions were controlled. We found that for realistic cryo-EM environments and relevant biomolecular systems, it is possible to recover the underlying free energy, with the pose accuracy and SNR as crucial determinants. We then use the method to study the conformational transitions of a calcium-activated channel with real cryo-EM particles. Interestingly, we recover not only the most probable conformation (used to generate a high-resolution reconstruction of the calcium-bound state) but also a metastable state that corresponds to the calcium-unbound conformation. As expected for turnover transitions within the same sample, the activation barriers are on the order of [Formula: see text]. We expect our tool for extracting free-energy profiles from cryo-EM images to enable more complete characterization of the thermodynamic ensemble of biomolecules.

Coming of Age: Cryo-Electron Tomography as a Versatile Tool to Generate High-Resolution Structures at Cellular/Biological Interfaces

Over the last few years, cryo electron microscopy has become the most important method in structural biology. While 80% of deposited maps are from single particle analysis, electron tomography has grown to become the second most important method. In particular sub-tomogram averaging has matured as a method, delivering structures between 2 and 5 Å from complexes in cells as well as in vitro complexes. While this resolution range is not standard, novel developments point toward a promising future. Here, we provide a guide for the workflow from sample to structure to gain insight into this emerging field.

Regulation of the activity of maize glutamate dehydrogenase by ammonium and potassium

Glutamate dehydrogenase (GDH) is an important enzyme in ammonium metabolism, the activity of which is regulated by multiple factors. In this study, we investigate the effects of ammonium and potassium on the activity of maize GDH. Our results show that both ammonium and potassium play multiple roles in regulating the activity of maize GDH, with the specific roles depending on the concentration of potassium. Together with the structural information of GDH, we propose models for the substrate inhibition of ammonium, and the elimination of substrate inhibition by potassium. These models are supported by the analysis of statistic thermodynamics. We also analyze the binding sites of ammonium and potassium on maize GDH, and the conformational changes of maize GDH. The findings provide insight into the regulation of maize GDH activity by ammonium and potassium and reveal the importance of the dose and ratio of nitrogen and potassium in crop cultivation.

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.

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.

Structural Studies of Glutamate Dehydrogenase (Isoform 1) From Arabidopsis thaliana, an Important Enzyme at the Branch-Point Between Carbon and Nitrogen Metabolism

Glutamate dehydrogenase (GDH) releases ammonia in a reversible NAD(P)⁺-dependent oxidative deamination of glutamate that yields 2-oxoglutarate (2OG). In current perception, GDH contributes to Glu homeostasis and plays a significant role at the junction of carbon and nitrogen assimilation pathways. GDHs are members of a superfamily of ELFV (Glu/Leu/Phe/Val) amino acid dehydrogenases and are subdivided into three subclasses, based on coenzyme specificity: NAD⁺-specific, NAD⁺/NADP⁺ dual-specific, and NADP⁺-specific. We determined in this work that the mitochondrial AtGDH1 isozyme from A. thaliana is NAD⁺-specific. Altogether, A. thaliana expresses three GDH isozymes (AtGDH1-3) targeted to mitochondria, of which AtGDH2 has an extra EF-hand motif and is stimulated by calcium. Our enzymatic assays of AtGDH1 established that its sensitivity to calcium is negligible. In vivo the AtGDH1-3 enzymes form homo- and heterohexamers of varied composition. We solved the crystal structure of recombinant AtGDH1 in the apo-form and in complex with NAD⁺ at 2.59 and 2.03 Å resolution, respectively. We demonstrate also that both in the apo form and in 1:1 complex with NAD⁺, it forms D ₃-symmetric homohexamers. A subunit of AtGDH1 consists of domain I, which is involved in hexamer formation and substrate binding, and of domain II which binds coenzyme. Most of the subunits in our crystal structures, including those in NAD⁺ complex, are in open conformation, with domain II forming a large (albeit variable) angle with domain I. One of the subunits of the AtGDH1-NAD⁺ hexamer contains a serendipitous 2OG molecule in the active site, causing a dramatic (∼25°) closure of the domains. We provide convincing evidence that the N-terminal peptide preceding domain I is a mitochondrial targeting signal, with a predicted cleavage site for mitochondrial processing peptidase (MPP) at Leu17-Leu18 that is followed by an unexpected potassium coordination site (Ser27, Ile30). We also identified several MPD [(+/-)-2-methyl-2,4-pentanediol] binding sites with conserved sequence. Although AtGDH1 is insensitive to MPD in our assays, the observation of druggable sites opens a potential for non-competitive herbicide design.