Structural Evolution of Primate Glutamate Dehydrogenase 2 as Revealed by In Silico Predictions and Experimentally Determined Structures

Glutamate dehydrogenase (GDH) interconverts glutamate to a-ketoglutarate and ammonia, interconnecting amino acid and carbohydrate metabolism. In humans, two functional GDH genes, GLUD1 and GLUD2, encode for hGDH1 and hGDH2, respectively. GLUD2 evolved from retrotransposition of the GLUD1 gene in the common ancestor of modern apes. These two isoenzymes are involved in the pathophysiology of human metabolic, neoplastic, and neurodegenerative disorders. The 3D structures of hGDH1 and hGDH2 have been experimentally determined; however, no information is available about the path of GDH2 structure changes during primate evolution. Here, we compare the structures predicted by the AlphaFold Colab method for the GDH2 enzyme of modern apes and their extinct primate ancestors. Also, we analyze the individual effect of amino acid substitutions emerging during primate evolution. Our most important finding is that the predicted structure of GDH2 in the common ancestor of apes was the steppingstone for the structural evolution of primate GDH2s. Two changes with a strong functional impact occurring at the first evolutionary step, Arg443Ser and Gly456Ala, had a destabilizing and stabilizing effect, respectively, making this step the most important one. Subsequently, GDH2 underwent additional modifications that fine-tuned its enzymatic properties to adapt to the functional needs of modern-day primate tissues.

The functional and structural characterization of Xanthomonas campestris pv. campestris core effector XopP revealed a new kinase activity

Exo70B1 is a protein subunit of the exocyst complex with a crucial role in a variety of cell mechanisms, including immune responses against pathogens. The calcium-dependent kinase 5 (CPK5) of Arabidopsis thaliana (hereafter Arabidopsis), phosphorylates AtExo70B1 upon functional disruption. We previously reported that, the Xanthomonas campestris pv. campestris effector XopP compromises AtExo70B1, while bypassing the host's hypersensitive response, in a way that is still unclear. Herein we designed an experimental approach, which includes biophysical, biochemical, and molecular assays and is based on structural and functional predictions, utilizing AplhaFold and DALI online servers, respectively, in order to characterize the in vivo XccXopP function. The interaction between AtExo70B1 and XccXopP was found very stable in high temperatures, while AtExo70B1 appeared to be phosphorylated at XccXopP-expressing transgenic Arabidopsis. XccXopP revealed similarities with known mammalian kinases and phosphorylated AtExo70B1 at Ser107, Ser111, Ser248, Thr309, and Thr364. Moreover, XccXopP protected AtExo70B1 from AtCPK5 phosphorylation. Together these findings show that XccXopP is an effector, which not only functions as a novel serine/threonine kinase upon its host target AtExo70B1 but also protects the latter from the innate AtCPK5 phosphorylation, in order to bypass the host's immune responses. Data are available via ProteomeXchange with the identifier PXD041405.

Structure of amino acid sequence-reversed wtRop protein: insights from atomistic molecular dynamics simulations

This study aims to the investigation of the advantages of designing new proteins presume upon a 'bias' sequence of amino acids, based on the reversed sequence of parent proteins, such as the retro ones. The structural simplicity of wtRop offers a very attractive model system to study these aspects. The current work is based on all-atom Molecular Dynamics (MD) simulations and corresponding experimental evidence on two different types of reversed wtRop protein, one with a fully reversed sequence of amino acids (rRop) and another with a partially reversed sequence (prRop), where only the five residues of the loop region (30ASP-34GLN) were not reversed. The exploration of the structure of the two retro proteins is performed highlighting the similarities and the differences with their parent protein, by employing various measures. Two models have been studied for both reversed proteins, a dimeric and a monomeric with the former one found to be more stable than the latter. Preferable equilibrium structures that the protein molecule can attain are explored, indicating the equilibration pathway. Simulation findings indicate a disruption of the α-helical structure and the appearance of additional secondary structures for both retro proteins. Reduced structural stability compared to their parent protein (wtRop) is also found. A corruption of the hydrophobic core is observed in the dimeric models. Furthermore, the simulations findings are consistent with the experimental characterization of prRop by circular dichroism spectroscopy (CD) which highlights an unstable, highly α-helical protein.Communicated by Ramaswamy H. Sarma.

Probing the conformational changes of in vivo overexpressed cell cycle regulator 6S ncRNA

The non-coding 6S RNA is a master regulator of the cell cycle in bacteria which binds to the RNA polymerase-σ70 holoenzyme during the stationary phase to inhibit transcription from the primary σ factor. Inhibition is reversed upon outgrowth from the stationary phase by synthesis of small product RNA transcripts (pRNAs). 6S and its complex with a pRNA were structurally characterized using Small Angle X-ray Scattering. The 3D models of 6S and 6S:pRNA complex presented here, demonstrate that the fairly linear and extended structure of 6S undergoes a major conformational change upon binding to pRNA. In particular, 6S:pRNA complex formation is associated with a compaction of the overall 6S size and an expansion of its central domain. Our structural models are consistent with the hypothesis that the resultant particle has a shape and size incompatible with binding to RNA polymerase-σ70. Overall, by use of an optimized in vivo methodological approach, especially useful for structural studies, our study considerably improves our understanding of the structural basis of 6S regulation by offering a mechanistic glimpse of the 6S transcriptional control.

Structure and Thermal Stability of wtRop and RM6 Proteins through All-Atom Molecular Dynamics Simulations and Experiments

In the current work we study, via molecular simulations and experiments, the folding and stability of proteins from the tertiary motif of 4-α-helical bundles, a recurrent motif consisting of four amphipathic α-helices packed in a parallel or antiparallel fashion. The focus is on the role of the loop region in the structure and the properties of the wild-type Rop (wtRop) and RM6 proteins, exploring the key factors which can affect them, through all-atom molecular dynamics (MD) simulations and supporting by experimental findings. A detailed investigation of structural and conformational properties of wtRop and its RM6 loopless mutation is presented, which display different physical characteristics even in their native states. Then, the thermal stability of both proteins is explored showing RM6 as more thermostable than wtRop through all studied measures. Deviations from native structures are detected mostly in tails and loop regions and most flexible residues are indicated. Decrease of hydrogen bonds with the increase of temperature is observed, as well as reduction of hydrophobic contacts in both proteins. Experimental data from circular dichroism spectroscopy (CD), are also presented, highlighting the effect of temperature on the structural integrity of wtRop and RM6. The central goal of this study is to explore on the atomic level how a protein mutation can cause major changes in its physical properties, like its structural stability.

α-Helices in the Type III Secretion Effectors: A Prevalent Feature with Versatile Roles

Type III Secretion Systems (T3SSs) are multicomponent nanomachines located at the cell envelope of Gram-negative bacteria. Their main function is to transport bacterial proteins either extracellularly or directly into the eukaryotic host cell cytoplasm. Type III Secretion effectors (T3SEs), latest to be secreted T3S substrates, are destined to act at the eukaryotic host cell cytoplasm and occasionally at the nucleus, hijacking cellular processes through mimicking eukaryotic proteins. A broad range of functions is attributed to T3SEs, ranging from the manipulation of the host cell's metabolism for the benefit of the bacterium to bypassing the host's defense mechanisms. To perform this broad range of manipulations, T3SEs have evolved numerous novel folds that are compatible with some basic requirements: they should be able to easily unfold, pass through the narrow T3SS channel, and refold to an active form when on the other side. In this review, the various folds of T3SEs are presented with the emphasis placed on the functional and structural importance of α-helices and helical domains.

Probing Protein Folding with Sequence-Reversed α-Helical Bundles

Recurrent protein folding motifs include various types of helical bundles formed by α-helices that supercoil around each other. While specific patterns of amino acid residues (heptad repeats) characterize the highly versatile folding motif of four-α-helical bundles, the significance of the polypeptide chain directionality is not sufficiently understood, although it determines sequence patterns, helical dipoles, and other parameters for the folding and oligomerization processes of bundles. To investigate directionality aspects in sequence-structure relationships, we reversed the amino acid sequences of two well-characterized, highly regular four-α-helical bundle proteins and studied the folding, oligomerization, and structural properties of the retro-proteins, using Circular Dichroism Spectroscopy (CD), Size Exclusion Chromatography combined with Multi-Angle Laser Light Scattering (SEC-MALS), and Small Angle X-ray Scattering (SAXS). The comparison of the parent proteins with their retro-counterparts reveals that while the α-helical character of the parents is affected to varying degrees by sequence reversal, the folding states, oligomerization propensities, structural stabilities, and shapes of the new molecules strongly depend on the characteristics of the heptad repeat patterns. The highest similarities between parent and retro-proteins are associated with the presence of uninterrupted heptad patterns in helical bundles sequences.

Crystal structure of glutamate dehydrogenase 2, a positively selected novel human enzyme involved in brain biology and cancer pathophysiology

INTRODUCTION

Mammalian glutamate dehydrogenase (hGDH1 in human cells) interconverts glutamate to α-ketoglutarate and ammonia while reducing NAD(P) to NAD(P)H. During primate evolution, humans and great apes have acquired hGDH2, an isoenzyme that underwent rapid evolutionary adaptation concomitantly with brain expansion, thereby acquiring unique catalytic and regulatory properties that permitted its function under conditions inhibitory to its ancestor hGDH1. Although the 3D-structures of GDHs, including hGDH1, have been determined, attempts to determine the hGDH2 structure were until recently unsuccessful. Comparison of the hGDH1/hGDH2 structures would enable a detailed understanding of their evolutionary differences. This work aimed at the determination of the hGDH2 crystal structure and the analysis of its functional implications. Recombinant hGDH2 was produced in the Spodoptera frugiperda ovarian cell line Sf21, using the Baculovirus expression system. Purification was achieved via a two-step chromatography procedure. hGDH2 was crystallized, X-ray diffraction data were collected using synchrotron radiation and the structure was determined by molecular replacement. The hGDH2 structure is reported at a resolution of 2.9 Å. The enzyme adopts a novel semi-closed conformation, which is an intermediate between known open and closed GDH1 conformations, differing from both. The structure enabled us to dissect previously reported biochemical findings and to structurally interpret the effects of evolutionary amino acid substitutions, including Arg470His, on ADP affinity. In conclusion, our data provide insights into the structural basis of hGDH2 properties, the functional evolution of hGDH isoenzymes, and open new prospects for drug design, especially for cancer therapeutics.

Probing the Structural Dynamics of the Catalytic Domain of Human Soluble Guanylate Cyclase

In the nitric oxide (NO) signaling pathway, human soluble guanylate cyclase (hsGC) synthesizes cyclic guanosine monophosphate (cGMP); responsible for the regulation of cGMP-specific protein kinases (PKGs) and phosphodiesterases (PDEs). The crystal structure of the inactive hsGC cyclase dimer is known, but there is still a lack of information regarding the substrate-specific internal motions that are essential for the catalytic mechanism of the hsGC. In the current study, the hsGC cyclase heterodimer complexed with guanosine triphosphate (GTP) and cGMP was subjected to molecular dynamics simulations, to investigate the conformational dynamics that have functional implications on the catalytic activity of hsGC. Results revealed that in the GTP-bound complex of the hsGC heterodimer, helix 1 of subunit α (α:h1) moves slightly inwards and comes close to helix 4 of subunit β (β:h4). This conformational change brings loop 2 of subunit β (β:L2) closer to helix 2 of subunit α (α:h2). Likewise, loop 2 of subunit α (α:L2) comes closer to helix 2 of subunit β (β:h2). These structural events stabilize and lock GTP within the closed pocket for cyclization. In the cGMP-bound complex, α:L2 detaches from β:h2 and establishes interactions with β:L2, which results in the loss of global structure compactness. Furthermore, with the release of pyrophosphate, the interaction between α:h1 and β:L2 weakens, abolishing the tight packing of the binding pocket. This study discusses the conformational changes induced by the binding of GTP and cGMP to the hsGC catalytic domain, valuable in designing new therapeutic strategies for the treatment of cardiovascular diseases.

Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands

Soluble guanylate cyclase (sGC) regulates numerous physiological processes. The β subunit Heme Nitric Oxide/Oxygen (HNOX) domain makes this protein sensitive to small gaseous ligands. The structural basis of the activation mechanism of sGC under the influence of ligands (NO, O₂, CO) is poorly understood. We examine the effect of different ligands on the human sGC HNOX domain. HNOX systems with gaseous ligands were generated and explored using Molecular Dynamics (MD). The distance between heme Fe2+ and histidine in the NO-ligated HNOX (NO-HNOX) system is larger compared to the O₂, CO systems. NO-HNOX rapidly adopts the conformation of the five-group metal coordination system. Loops α, β, γ and helix-f exhibit increased mobility and different hydrogen bond networks in NO-HNOX compared to the other systems. The removal of His from the Fe coordination sphere in NO-HNOX is assisted by interaction of the imidazole ring with the surrounding residues which in turn leads to the release of signaling helix-f and activation of the sGC enzyme. Insights into the conformational dynamics of a human sGC HNOX domain, especially for regions which are functionally critical for signal transduction, are valuable in the understanding of cardiovascular diseases.

Expression, purification and crystallization of a protein resulting from the inversion of the amino-acid sequence of a helical bundle

Earlier studies have found that the occurrence of inverse sequence identity in proteins is not indicative of three-dimensional similarity, but rather leads to different folds or unfolded proteins. Short helices, however, frequently keep their conformations when their sequences are inverted. To explore the impact of sequence inversion on long helices, revRM6, with the inverse amino-acid sequence relative to RM6, a highly stable variant of the ColE1 Rop protein, was engineered. RM6 is a highly regular four-α-helical bundle that serves as a model system for protein-folding studies. Here, the crystallization and preliminary crystallographic characterization of revRM6 are reported. The protein was overexpressed in Escherichia coli, purified to homogeneity and crystallized. The crystals belonged to space group P41212, with unit-cell parameters a = b = 44.98, c = 159.74 Å, and diffracted to a resolution of 3.45 Å.

Molecular Dynamics Simulations of BcZBP, A Deacetylase from Bacillus cereus: Active Site Loops Determine Substrate Accessibility and Specificity

BcZBP is an LmbE-like, homohexameric, zinc-dependent deacetylase from the opportunistic pathogen Bacillus cereus with three, thus far uncharacterized, homologues in B. anthracis. Although its specific substrate is still unknown, the enzyme has been shown to preferentially deacetylate N-acetylglucosamine and diacetylchitobiose via an active site based on a zinc-binding motif of the type HXDDXnH. In the crystal structure, the active site is located at a deep and partially blocked cleft formed at the interface between monomers related by the molecular 3-fold axis, although the major, in structural terms, building block of the enzyme is not the trimer, but the intertwined dimer. Here, we report results from a 50 ns molecular dynamics simulation of BcZBP in explicit solvent with full electrostatics and show that (i) the view of the intertwined dimer as the major structural and functional building block of this class of hexameric enzymes is possibly an oversimplification of the rather complex dynamics observed in the simulation, (ii) the most mobile (with respect to their atomic fluctuations) parts of the structure coincide with three surface loops surrounding the active site, and (iii) these mobile loops define the active site's accessibility, and may be implicated in the determination of the enzyme's specificity.

HrpG and HrpV proteins from the Type III secretion system of Erwinia amylovora form a stable heterodimer

Bacterial type III secretion systems (T3SSs) are specialized multicomponent nanomachines that mediate the transport of proteins either to extracellular locations or directly into eukaryotic host cell cytoplasm. Erwinia amylovora, the main agent of rosaceous plants fireblight disease, employs an Hrp/Hrc1 T3SS to accomplish its pathogenesis. The regulatory network that controls the activation of this T3SS is largely unknown in E. amylovora. However, in Pseudomonas syringae pathovars, the HrpG/HrpV complex has been shown to directly regulate the activity of transcription factor HrpS and consequently the upregulation of the Hrp/Hrc1 T3SS related genes. In this work, we report the successful recombinant production and purification of a stable E. amylovora HrpG/HrpV complex, using pPROpET, a bicistronic expression vector. Furthermore, we present the first solution structure of this complex based on small-angle X-ray scattering data.

The benefit of the European User Community from transnational access to national radiation facilities

Transnational access (TNA) to national radiation sources is presently provided via programmes of the European Commission by BIOSTRUCT-X and CALIPSO with a major benefit for scientists from European countries. Entirely based on scientific merit, TNA allows all European scientists to realise synchrotron radiation experiments for addressing the Societal Challenges promoted in HORIZON2020. In addition, by TNA all European users directly take part in the development of the research infrastructure of facilities. The mutual interconnection of users and facilities is a strong prerequisite for future development of the research infrastructure of photon science. Taking into account the present programme structure of HORIZON2020, the European Synchrotron User Organization (ESUO) sees considerable dangers for the continuation of this successful collaboration in the future.

Structure determination through homology modelling and torsion-angle simulated annealing: application to a polysaccharide deacetylase from Bacillus cereus

The structure of BC0361, a polysaccharide deacetylase from Bacillus cereus, has been determined using an unconventional molecular-replacement procedure. Tens of putative models of the C-terminal domain of the protein were constructed using a multitude of homology-modelling algorithms, and these were tested for the presence of signal in molecular-replacement calculations. Of these, only the model calculated by the SAM-T08 server gave a consistent and convincing solution, but the resulting model was too inaccurate to allow phase determination to proceed to completion. The application of slow-cooling torsion-angle simulated annealing (started from a very high temperature) drastically improved this initial model to the point of allowing phasing through cycles of model building and refinement to be initiated. The structure of the protein is presented with emphasis on the presence of a C(α)-modified proline at its active site, which was modelled as an α-hydroxy-L-proline.

Phylogenetic analysis of a gene cluster encoding an additional, rhizobial-like type III secretion system that is narrowly distributed among Pseudomonas syringae strains

Background

The central role of Type III secretion systems (T3SS) in bacteria-plant interactions is well established, yet unexpected findings are being uncovered through bacterial genome sequencing. Some Pseudomonas syringae strains possess an uncharacterized cluster of genes encoding putative components of a second T3SS (T3SS-2) in addition to the well characterized Hrc1 T3SS which is associated with disease lesions in host plants and with the triggering of hypersensitive response in non-host plants. The aim of this study is to perform an in silico analysis of T3SS-2, and to compare it with other known T3SSs.

Results

Based on phylogenetic analysis and gene organization comparisons, the T3SS-2 cluster of the P. syringae pv. phaseolicola strain is grouped with a second T3SS found in the pNGR234b plasmid of Rhizobium sp. These additional T3SS gene clusters define a subgroup within the Rhizobium T3SS family. Although, T3SS-2 is not distributed as widely as the Hrc1 T3SS in P. syringae strains, it was found to be constitutively expressed in P. syringae pv phaseolicola through RT-PCR experiments.

Conclusions

The relatedness of the P. syringae T3SS-2 to a second T3SS from the pNGR234b plasmid of Rhizobium sp., member of subgroup II of the rhizobial T3SS family, indicates common ancestry and/or possible horizontal transfer events between these species. Functional analysis and genome sequencing of more rhizobia and P. syringae pathovars may shed light into why these bacteria maintain a second T3SS gene cluster in their genome.

Purification, crystallization and preliminary X-ray diffraction analysis of the C-terminal fragment of the MvfR protein from Pseudomonas aeruginosa

The LysR-type transcriptional regulator MvfR plays a critical role in Pseudomonas aeruginosa pathogenicity via the transcriptional regulation of multiple quorum-sensing-regulated virulence factors. The protein also controls pathogenic type VI secretion loci. MvfRC87, a 242-residue C-terminal segment of MvfR, was produced in Escherichia coli, purified and crystallized. X-ray diffraction data were collected using synchrotron radiation and crystallographic parameters were determined.

The rotation-coupled sliding of EcoRV

It has been proposed that certain type II restriction enzymes (REs), such as EcoRV, track the helical pitch of DNA as they diffuse along DNA, a so-called rotation-coupled sliding. As of yet, there is no direct experimental observation of this phenomenon, but mounting indirect evidence gained from single-molecule imaging of RE-DNA complexes support the hypothesis. We address this issue by conjugating fluorescent labels of varying size (organic dyes, proteins and quantum dots) to EcoRV, and by fusing it to the engineered Rop protein scRM6. Single-molecule imaging of these modified EcoRVs sliding along DNA provides us with their linear diffusion constant (D(1)), revealing a significant size dependency. To account for the dependence of D(1) on the size of the EcoRV label, we have developed four theoretical models describing different types of motion along DNA and find that our experimental results are best described by rotation-coupled sliding of the protein. The similarity of EcoRV to other type II REs and DNA binding proteins suggests that this type of motion could be widely preserved in other biological contexts.

In silico analysis reveals multiple putative type VI secretion systems and effector proteins in Pseudomonas syringae pathovars

Type VI secretion systems (T6SS) of Gram-negative bacteria form injectisomes that have the potential to translocate effector proteins into eukaryotic host cells. In silico analysis of the genomes in six Pseudomonas syringae pathovars revealed that P. syringae pv. tomato DC3000, pv. tabaci ATCC 11528, pv. tomato T1 and pv. oryzae 1-6 each carry two putative T6SS gene clusters (HSI-I and HSI-II; HSI: Hcp secretion island), whereas pv. phaseolicola 1448A and pv. syringae B728 each carry one. The pv. tomato DC3000 HSI-I and pv. tomato T1 HSI-II possess a highly similar organization and nucleotide sequence, whereas the pv. tomato DC3000, pv. oryzae 1-6 and pv. tabaci 11528 HSI-II are more divergent. Putative effector orthologues vary in number among the strains examined. The Clp-ATPases and IcmF orthologues form distinct phylogenetic groups: the proteins from pv. tomato DC3000, pv. tomato T1, pv. oryzae and pv. tabaci 11528 from HSI-II group together with most orthologues from other fluorescent pseudomonads, whereas those from pv. phaseolicola, pv. syringae, pv. tabaci, pv. tomato T1 and pv. oryzae from HSI-I group closer to the Ralstonia solanacearum and Xanthomonas orthologues. Our analysis suggests multiple independent acquisitions and possible gene attrition/loss of putative T6SS genes by members of P. syringae.