Significance of the Protein Interface Configuration for Allostery in Imidazole Glycerol Phosphate Synthase

Catalytic Effects of Active Site Conformational Change in the Allosteric Activation of Imidazole Glycerol Phosphate Synthase

Imidazole glycerol phosphate synthase (IGPS) is a class-I glutamine amidotransferase (GAT) that hydrolyzes glutamine. Ammonia is produced and transferred to a second active site, where it reacts with N1-(5'-phosphoribosyl)-formimino-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) to form precursors to purine and histidine biosynthesis. Binding of PrFAR over 25 Å away from the active site increases glutaminase efficiency by ∼4500-fold, primarily altering the glutamine turnover number. IGPS has been the focus of many studies on allosteric communication; however, atomic details for how the glutamine hydrolysis rate increases in the presence of PrFAR are lacking. We present a density functional theory study on 237-atom active site cluster models of IGPS based on crystallized structures representing the inactive and allosterically active conformations and investigate the multistep reaction leading to thioester formation and ammonia production. The proposed mechanism is supported by similar, well-studied enzyme mechanisms, and the corresponding energy profile is consistent with steady-state kinetic studies of PrFAR + IGPS. Additional active site models are constructed to examine the relationship between active site structural change and transition-state stabilization via energy decomposition schemes. The results reveal that the inactive IGPS conformation does not provide an adequately formed oxyanion hole structure and that repositioning of the oxyanion strand relative to the substrate is vital for a catalysis-competent oxyanion hole, with or without the hVal51 dihedral flip. These findings are valuable for future endeavors in modeling the IGPS allosteric mechanism by providing insight into the atomistic changes required for rate enhancement that can inform suitable reaction coordinates for subsequent investigations.

Connected Component Analysis of Dynamical Perturbation Contact Networks

Describing protein dynamical networks through amino acid contacts is a powerful way to analyze complex biomolecular systems. However, due to the size of the systems, identifying the relevant features of protein-weighted graphs can be a difficult task. To address this issue, we present the connected component analysis (CCA) approach that allows for fast, robust, and unbiased analysis of dynamical perturbation contact networks (DPCNs). We first illustrate the CCA method as applied to a prototypical allosteric enzyme, the imidazoleglycerol phosphate synthase (IGPS) enzyme from Thermotoga maritima bacteria. This approach was shown to outperform the clustering methods applied to DPCNs, which could not capture the propagation of the allosteric signal within the protein graph. On the other hand, CCA reduced the DPCN size, providing connected components that nicely describe the allosteric propagation of the signal from the effector to the active sites of the protein. By applying the CCA to the IGPS enzyme in different conditions, i.e., at high temperature and from another organism (yeast IGPS), and to a different enzyme, i.e., a protein kinase, we demonstrated how CCA of DPCNs is an effective and transferable tool that facilitates the analysis of protein-weighted networks.

MDiGest: A Python package for describing allostery from molecular dynamics simulations

Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.

Rational design, production and in vitro analysis of photoxenoproteins

In synthetic biology, the artificial control of proteins by light is of growing interest since it enables the spatio-temporal regulation of downstream molecular processes. This precise photocontrol can be established by the site-directed incorporation of photo-sensitive non-canonical amino acids (ncAAs) into proteins, which generates so-called photoxenoproteins. Photoxenoproteins can be engineered using ncAAs that facilitate the irreversible activation or reversible regulation of their activity upon irradiation. In this chapter, we provide a general outline of the engineering process based on the current methodological state-of-the-art to obtain artificial photocontrol in proteins using the ncAAs o-nitrobenzyl-O-tyrosine as example for photocaged ncAAs (irreversible), and phenylalanine-4'-azobenzene as example for photoswitchable ncAAs (reversible). We thereby focus on the initial design as well as the production and characterization of photoxenoproteins in vitro. Finally, we outline the analysis of photocontrol under steady-state and non-steady-state conditions using the allosteric enzyme complexes imidazole glycerol phosphate synthase and tryptophan synthase as examples.

Time Evolution of the Millisecond Allosteric Activation of Imidazole Glycerol Phosphate Synthase

Deciphering the molecular mechanisms of enzymatic allosteric regulation requires the structural characterization of functional states and also their time evolution toward the formation of the allosterically activated ternary complex. The transient nature and usually slow millisecond time scale interconversion between these functional states hamper their experimental and computational characterization. Here, we combine extensive molecular dynamics simulations, enhanced sampling techniques, and dynamical networks to describe the allosteric activation of imidazole glycerol phosphate synthase (IGPS) from the substrate-free form to the active ternary complex. IGPS is a heterodimeric bienzyme complex whose HisH subunit is responsible for hydrolyzing glutamine and delivering ammonia for the cyclase activity in HisF. Despite significant advances in understanding the underlying allosteric mechanism, essential molecular details of the long-range millisecond allosteric activation of IGPS remain hidden. Without using a priori information of the active state, our simulations uncover how IGPS, with the allosteric effector bound in HisF, spontaneously captures glutamine in a catalytically inactive HisH conformation, subsequently attains a closed HisF:HisH interface, and finally forms the oxyanion hole in HisH for efficient glutamine hydrolysis. We show that the combined effector and substrate binding dramatically decreases the conformational barrier associated with oxyanion hole formation, in line with the experimentally observed 4500-fold activity increase in glutamine hydrolysis. The allosteric activation is controlled by correlated time-evolving dynamic networks connecting the effector and substrate binding sites. This computational strategy tailored to describe millisecond events can be used to rationalize the effect of mutations on the allosteric regulation and guide IGPS engineering efforts.

Metal-responsive regulation of enzyme catalysis using genetically encoded chemical switches

Dynamic control over protein function is a central challenge in synthetic biology. To address this challenge, we describe the development of an integrated computational and experimental workflow to incorporate a metal-responsive chemical switch into proteins. Pairs of bipyridinylalanine (BpyAla) residues are genetically encoded into two structurally distinct enzymes, a serine protease and firefly luciferase, so that metal coordination biases the conformations of these enzymes, leading to reversible control of activity. Computational analysis and molecular dynamics simulations are used to rationally guide BpyAla placement, significantly reducing experimental workload, and cell-free protein synthesis coupled with high-throughput experimentation enable rapid prototyping of variants. Ultimately, this strategy yields enzymes with a robust 20-fold dynamic range in response to divalent metal salts over 24 on/off switches, demonstrating the potential of this approach. We envision that this strategy of genetically encoding chemical switches into enzymes will complement other protein engineering and synthetic biology efforts, enabling new opportunities for applications where precise regulation of protein function is critical.

A guide to designing photocontrol in proteins: methods, strategies and applications

Light is essential for various biochemical processes in all domains of life. In its presence certain proteins inside a cell are excited, which either stimulates or inhibits subsequent cellular processes. The artificial photocontrol of specifically proteins is of growing interest for the investigation of scientific questions on the organismal, cellular and molecular level as well as for the development of medicinal drugs or biocatalytic tools. For the targeted design of photocontrol in proteins, three major methods have been developed over the last decades, which employ either chemical engineering of small-molecule photosensitive effectors (photopharmacology), incorporation of photoactive non-canonical amino acids by genetic code expansion (photoxenoprotein engineering), or fusion with photoreactive biological modules (hybrid protein optogenetics). This review compares the different methods as well as their strategies and current applications for the light-regulation of proteins and provides background information useful for the implementation of each technique.

Distinct allosteric pathways in imidazole glycerol phosphate synthase from yeast and bacteria

Understanding the relationship between protein structures and their function is still an open question that becomes very challenging when allostery plays an important functional role. Allosteric proteins, in fact, exploit different ranges of motions (from sidechain local fluctuations to long-range collective motions) to effectively couple distant binding sites, and of particular interest is whether allosteric proteins of the same families with similar functions and structures also necessarily share the same allosteric mechanisms. Here, we compared the early dynamics initiating the allosteric communication of a prototypical allosteric enzyme from two different organisms, i.e., the imidazole glycerol phosphate synthase (IGPS) enzymes from the thermophilic bacteria and the yeast, working at high and room temperatures, respectively. By combining molecular dynamics simulations and network models derived from graph theory, we found rather distinct early allosteric dynamics in the IGPS from the two organisms, involving significatively different allosteric pathways in terms of both local and collective motions. Given the successful prediction of key allosteric residues in the bacterial IGPS, whose mutation disrupts its allosteric communication, the outcome of this study paves the way for future experimental studies on the yeast IGPS that could foster therapeutic applications by exploiting the control of IGPS enzyme allostery.

Molecular photoswitches in aqueous environments

Molecular photoswitches enable dynamic control of processes with high spatiotemporal precision, using light as external stimulus, and hence are ideal tools for different research areas spanning from chemical biology to smart materials. Photoswitches are typically organic molecules that feature extended aromatic systems to make them responsive to (visible) light. However, this renders them inherently lipophilic, while water-solubility is of crucial importance to apply photoswitchable organic molecules in biological systems, like in the rapidly emerging field of photopharmacology. Several strategies for solubilizing organic molecules in water are known, but there are not yet clear rules for applying them to photoswitchable molecules. Importantly, rendering photoswitches water-soluble has a serious impact on both their photophysical and biological properties, which must be taken into consideration when designing new systems. Altogether, these aspects pose considerable challenges for successfully applying molecular photoswitches in aqueous systems, and in particular in biologically relevant media. In this review, we focus on fully water-soluble photoswitches, such as those used in biological environments, in both in vitro and in vivo studies. We discuss the design principles and prospects for water-soluble photoswitches to inspire and enable their future applications.

Computational design and experimental characterization of a photo-controlled mRNA-cap guanine-N7 methyltransferase

The spatial and temporal control of gene expression at the post-transcriptional level is essential in eukaryotic cells and developing multicellular organisms. In recent years optochemical and optogenetic tools have enabled the manipulation and investigation of many steps in the involved processes. However, examples for light-mediated control of eukaryotic mRNA processing and the responsible enzymes are still rare. In particular, methylation of the 5′ cap of mRNA is required for ribosome assembly, and the responsible guanine-N7 methyltransferase (MTase) from E. cuniculi (Ecm1) proved suitable for activating translation. Here, we report on a photoswitchable MTase obtained by bridging the substrate-binding cleft of Ecm1 with a tetra-ortho-methoxy-azobenzene. This azobenzene derivative is characterized by efficient trans-to-cis isomerization using red light at 615 nm. Starting from a cysteine-free Ecm1 variant (ΔCys), we used a computational approach to identify suitable conjugation sites for the azobenzene moiety. We created and characterized the four best-ranked variants, each featuring two appropriately positioned cysteines close to the substrate-binding cleft. Conjugating and crosslinking the azobenzene between C149/C155 in a designed Ecm1 variant (VAR3-Az) enabled light-dependent modulation of the MTase activity and showed a 50% higher activity for the cis form than the trans-form of the azobenzene conjugated to VAR3-Az.

Guided by computational design, we engineered a light-dependent 5' cap guanine-N7 methyltransferase by bridging the substrate-binding cleft with an azobenzene derivative.

Molecular basis for the allosteric activation mechanism of the heterodimeric imidazole glycerol phosphate synthase complex

Imidazole glycerol phosphate synthase (HisFH) is a heterodimeric bienzyme complex operating at a central branch point of metabolism. HisFH is responsible for the HisH-catalyzed hydrolysis of glutamine to glutamate and ammonia, which is then used for a cyclase reaction by HisF. The HisFH complex is allosterically regulated but the underlying mechanism is not well understood. Here, we elucidate the molecular basis of the long range, allosteric activation of HisFH. We establish that the catalytically active HisFH conformation is only formed when the substrates of both HisH and HisF are bound. We show that in this conformation an oxyanion hole in the HisH active site is established, which rationalizes the observed 4500-fold allosteric activation compared to the inactive conformation. In solution, the inactive and active conformations are in a dynamic equilibrium and the HisFH turnover rates correlate with the population of the active conformation, which is in accordance with the ensemble model of allostery.

The allosteric regulation of the bienzyme complex imidazole glycerol phosphate synthase (HisFH) remains to be elucidated. Here, the authors provide structural insights into the dynamic allosteric mechanism by which ligand binding to the cyclase and glutaminase active sites of HisFH regulate enzyme activation.

Computational investigations of allostery in aromatic amino acid biosynthetic enzymes

Allostery, in which binding of ligands to remote sites causes a functional change in the active sites, is a fascinating phenomenon observed in enzymes. Allostery can occur either with or without significant conformational changes in the enzymes, and the molecular basis of its mechanism can be difficult to decipher using only experimental techniques. Computational tools for analyzing enzyme sequences, structures, and dynamics can provide insights into the allosteric mechanism at the atomic level. Combining computational and experimental methods offers a powerful strategy for the study of enzyme allostery. The aromatic amino acid biosynthesis pathway is essential in microorganisms and plants. Multiple enzymes involved in this pathway are sensitive to feedback regulation by pathway end products and are known to use allostery to control their activities. To date, four enzymes in the aromatic amino acid biosynthesis pathway have been computationally investigated for their allosteric mechanisms, including 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase, anthranilate synthase, chorismate mutase, and tryptophan synthase. Here we review the computational studies and findings on the allosteric mechanisms of these four enzymes. Results from these studies demonstrate the capability of computational tools and encourage future computational investigations of allostery in other enzymes of this pathway.