Structural evidence for ammonia tunneling across the (beta alpha)(8) barrel of the imidazole glycerol phosphate synthase bienzyme complex

Illuminating Protein Allostery by Chemically Accurate Contact Response Analysis (ChACRA)

Decoding allostery at the atomic level is essential for understanding the relationship between a protein's sequence, structure, and dynamics. Recently, we have shown that decomposing temperature responses of inter-residue contacts can reveal allosteric couplings and provide useful insight into the functional dynamics of proteins. The details of this Chemically Accurate Contact Response Analysis (ChACRA) are presented here along with its application to two well-known allosteric proteins. The first protein, IGPS, is a model of ensemble allostery that lacks clear structural differences between the active and inactive states. We show that the application of ChACRA reveals the experimentally identified allosteric coupling between effector and active sites of IGPS. The second protein, ATCase, is a classic example of allostery with distinct active and inactive structural states. Using ChACRA, we directly identify the most significant residue level interactions underlying the enzyme's cooperative behavior. Both test cases demonstrate the utility of ChACRA's unsupervised machine learning approach for dissecting allostery at the residue level.

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.

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.

Turning up the heat mimics allosteric signaling in imidazole-glycerol phosphate synthase

Allosteric drugs have the potential to revolutionize biomedicine due to their enhanced selectivity and protection against overdosage. However, we need to better understand allosteric mechanisms in order to fully harness their potential in drug discovery. In this study, molecular dynamics simulations and nuclear magnetic resonance spectroscopy are used to investigate how increases in temperature affect allostery in imidazole glycerol phosphate synthase. Results demonstrate that temperature increase triggers a cascade of local amino acid-to-amino acid dynamics that remarkably resembles the allosteric activation that takes place upon effector binding. The differences in the allosteric response elicited by temperature increase as opposed to effector binding are conditional to the alterations of collective motions induced by either mode of activation. This work provides an atomistic picture of temperature-dependent allostery, which could be harnessed to more precisely control enzyme function.

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.

A simple method to determine changes in the affinity between HisF and HisH in the Imidazole Glycerol Phosphate Synthase heterodimer

The bi-enzyme HisF-HisH heterodimer is part of the pathway that produces histidine and purines in bacteria and lower eukaryotes, but it is absent in mammals. This heterodimer has been largely studied probing the basis of the allosteric effects and the structural stability in proteins. It is also a potential target for antibacterial drugs. In this work, we developed a simple method to evaluate changes in the affinity between HisF and HisH in the heterodimer of the bacteria Thermotoga maritima. HisH contains a single tryptophan residue, which is exposed in the free protein, but buried in the heterodimer interface. Hence, the intrinsic fluorescence maximum of this residue changes to shorter wavelengths upon dimerization. Thus, we used the fluorescence intensity at this shorter wavelength to monitor heterodimer accumulation when HisH was combined with sub-stoichiometric HisF. Under conditions where the HisF-HisH heterodimer is in equilibrium with the free states of these enzymes, when [HisH] > [HisF], we deduced a linear function connecting [HisF-HisH] to [HisF], in which the slope depends on the heterodimer dissociation constant (Kd). Based on this equation, taking fluorescence intensities as proxies of the heterodimer and HisF concentrations, we experimentally determined the Kd at four different temperatures. These Kd values were compared to those evaluated using ITC. Both methods revealed an increase in the HisF and HisH binding affinity as the temperature increases. In spite of differences in their absolute values, the Kd determined using these methods presented an evident linear correlation. To demonstrate the effectiveness of the fluorescence method we determined the effect on the Kd caused by 12 single mutations in HisF. Coherently, this test singled out the only mutation in the binding interface. In brief, the method described here effectively probes qualitative effects on the Kd, can be carried out using common laboratory equipment and is scalable.

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.

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.

Effect of Non-Lethal Selection on Spontaneous Revertants of Frameshift Mutations: The Escherichia coli hisF Case

Microorganisms possess the potential to adapt to fluctuations in environmental parameters, and their evolution is driven by the continuous generation of mutations. The reversion of auxotrophic mutations has been widely studied; however, little is known about the reversion of frameshift mutations resulting in amino acid auxotrophy and on the structure and functioning of the protein encoded by the revertant mutated gene. The aims of this work were to analyze the appearance of reverse mutations over time and under different selective pressures and to investigate revertant enzymes' three-dimensional structures and their correlation with a different growth ability. Escherichia coli FB182 strain, carrying the hisF892 single nucleotide deletion resulting in histidine auxotrophy, was subjected to different selective pressures, and revertant mutants were isolated and characterized. The obtained results allowed us to identify different indels of different lengths located in different positions in the hisF gene, and relations with the incubation time and the selective pressure applied were observed. Moreover, the structure of the different mutant proteins was consistent with the respective revertant ability to grow in absence of histidine, highlighting a correlation between the mutations and the catalytic activity of the mutated HisF enzyme.

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.

The Catalytic Mechanism of Pdx2 Glutaminase Driven by a Cys-His-Glu Triad: A Computational Study

The catalytic mechanism of Pdx2 was studied with atomic detail employing the computational ONIOM hybrid QM/MM methodology. Pdx2 employs a Cys-His-Glu catalytic triad to deaminate glutamine to glutamate and ammonia - the source of the nitrogen of pyridoxal 5'-phosphate (PLP). This enzyme is, therefore, a rate-limiting step in the PLP biosynthetic pathway of Malaria and Tuberculosis pathogens that rely on this mechanism to obtain PLP. For this reason, Pdx2 is considered a novel and promising drug target to treat these diseases. The results obtained show that the catalytic mechanism of Pdx2 occurs in six steps that can be divided into four stages: (i) activation of Cys₈₇ , (ii) deamination of glutamine with the formation of the glutamyl-thioester intermediate, (iii) hydrolysis of the formed intermediate, and (iv) enzymatic turnover. The kinetic data available in the literature (19.1-19.5 kcal mol^-1 ) agree very well with the calculated free energy barrier of the hydrolytic step (18.2 kcal.mol^-11 ), which is the rate-limiting step of the catalytic process when substrate is readily available in the active site. This catalytic mechanism differs from other known amidases in three main points: i) it requires the activation of the nucleophile Cys₈₇ to a thiolate; ii) the hydrolysis occurs in a single step and therefore does not require the formation of a second tetrahedral reaction intermediate, as it is proposed, and iii) Glu₁₉₈ does not have a direct role in the catalytic process. Together, these results can be used for the synthesis of new transition state analogue inhibitors capable of inhibiting Pdx2 and impair diseases like Malaria and Tuberculosis.

Role of a high centrality residue in protein dynamics and thermal stability

Centralities determined from Residue Interaction Networks (RIN) in proteins have been used to predict aspects of their structure and dynamics. Here, we correlate the Eigenvector Centrality (Ec) with the rate constant for thermal denaturation (kden) of the HisF protein from Thermotoga maritima based on 12 single alanine substitution mutants. The molecular basis for this correlation was further explored by studying a mutant containing a replacement of a high Ec residue, Y182A, which displayed increased kden at 80 °C. The crystallographic structure of this mutant showed few changes, mostly in two flexible loops. The 1H-15N -HSQC showed only subtle changes of cross peak positions for residues located near the mutation site and scattered throughout the structure. However, the comparison of the RIN showed that Y182 is the vertex of a set of high centrality residues that spreads throughout the HisF structure, which is lacking in the mutant. Cross-correlation displacements of Cα calculated from a molecular dynamics simulation at different temperatures showed that the Y182A mutation reduced the correlated movements in the HisF structure above 70 °C. 1H-15N NMR chemical shift covariance using temperature as perturbation were consistent with these results. In conclusion the increase in temperature drives the structure of the mutant HisF-Y182A into a less connected state, richer in non-concerted motions, located predominantly in the C-terminal half of the protein where Y182 is placed. Conversely, wild-type HisF responds to increased temperature as a single unit. Hence the replacement of a high Ec residue alters the distribution of thermal energy through HisF structure.

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.

Residue-Level Allostery Propagates through the Effective Coarse-Grained Hessian

The long-ranged coupling between residues that gives rise to allostery in a protein is built up from short-ranged physical interactions. Computational tools used to predict this coupling and its functional relevance have relied on the application of graph theoretical metrics to residue-level correlations measured from all-atom molecular dynamics simulations. The short-ranged interactions that yield these long-ranged residue-level correlations are quantified by the effective coarse-grained Hessian. Here we compute an effective harmonic coarse-grained Hessian from simulations of a benchmark allosteric protein, IGPS, and demonstrate the improved locality of this graph Laplacian over two other connectivity matrices. Additionally, two centrality metrics are developed that indicate the direct and indirect importance of each residue at producing the covariance between the effector binding pocket and the active site. The residue importance indicated by these two metrics is corroborated by previous mutagenesis experiments and leads to unique functional insights; in contrast to previous computational analyses, our results suggest that fP76-hK181 is the most important contact for conveying direct allosteric paths across the HisF-HisH interface. The connectivity around fD98 is found to be important at affecting allostery through indirect means.

Development of a Force Field for the Simulation of Single-Chain Proteins and Protein-Protein Complexes

The accuracy of atomistic physics-based force fields for the simulation of biological macromolecules has typically been benchmarked experimentally using biophysical data from simple, often single-chain systems. In the case of proteins, the careful refinement of force field parameters associated with torsion-angle potentials and the use of improved water models have enabled a great deal of progress toward the highly accurate simulation of such monomeric systems in both folded and, more recently, disordered states. In living organisms, however, proteins constantly interact with other macromolecules, such as proteins and nucleic acids, and these interactions are often essential for proper biological function. Here, we show that state-of-the-art force fields tuned to provide an accurate description of both ordered and disordered proteins can be limited in their ability to accurately describe protein-protein complexes. This observation prompted us to perform an extensive reparameterization of one variant of the Amber protein force field. Our objective involved refitting not only the parameters associated with torsion-angle potentials but also the parameters used to model nonbonded interactions, the specification of which is expected to be central to the accurate description of multicomponent systems. The resulting force field, which we call DES-Amber, allows for more accurate simulations of protein-protein complexes, while still providing a state-of-the-art description of both ordered and disordered single-chain proteins. Despite the improvements, calculated protein-protein association free energies still appear to deviate substantially from experiment, a result suggesting that more fundamental changes to the force field, such as the explicit treatment of polarization effects, may simultaneously further improve the modeling of single-chain proteins and protein-protein complexes.

Light Regulation of Enzyme Allostery through Photo-responsive Unnatural Amino Acids

Imidazole glycerol phosphate synthase (ImGPS) is an allosteric bienzyme complex in which substrate binding to the synthase subunit HisF stimulates the glutaminase subunit HisH. To control this stimulation with light, we have incorporated the photo-responsive unnatural amino acids phenylalanine-4'-azobenzene (AzoF), o-nitropiperonyl-O-tyrosine (NPY), and methyl-o-nitropiperonyllysine (mNPK) at strategic positions of HisF. The light-mediated isomerization of AzoF at position 55 (fS55AzoF^E ↔ fS55AzoF^Z) resulted in a reversible 10-fold regulation of HisH activity. The light-mediated decaging of NPY at position 39 (fY39NPY → fY39) and of mNPK at position 99 (fK99mNPK → fK99) led to a 4- to 6-fold increase of HisH activity. Molecular dynamics simulations explained how the unnatural amino acids interfere with the allosteric machinery of ImGPS and revealed additional aspects of HisH stimulation in wild-type ImGPS. Our findings show that unnatural amino acids can be used as a powerful tool for the spatiotemporal control of a central metabolic enzyme complex by light.

Exploring Allosteric Pathways of a V-Type Enzyme with Dynamical Perturbation Networks

Elucidation of the allosteric pathways in proteins is a computational challenge that strongly benefits from combination of atomistic molecular dynamics (MD) simulations and coarse-grained analysis of the complex dynamical network of chemical interactions based on graph theory. Here, we introduce and assess the performances of the dynamical perturbation network analysis of allosteric pathways in a prototypical V-type allosteric enzyme. Dynamical atomic contacts obtained from MD simulations are used to weight the allosteric protein graph, which involves an extended network of contacts perturbed by the effector binding in the allosteric site. The outcome showed good agreement with previously reported theoretical and experimental extended studies and it provided recognition of new potential allosteric spots that can be exploited in future mutagenesis experiments. Overall, the dynamical perturbation network analysis proved to be a powerful computational tool, complementary to other network-based approaches that can assist the full exploitation of allosteric phenomena for advances in protein engineering and rational drug design.

Artificial Light Regulation of an Allosteric Bienzyme Complex by a Photosensitive Ligand

The artificial regulation of proteins by light is an emerging subdiscipline of synthetic biology. Here, we used this concept to photocontrol both catalysis and allostery within the heterodimeric enzyme complex imidazole glycerol phosphate synthase (ImGP-S). ImGP-S consists of the cyclase subunit HisF and the glutaminase subunit HisH, which is allosterically stimulated by substrate binding to HisF. We show that a light-sensitive diarylethene (1,2-dithienylethene, DTE)-based competitive inhibitor in its ring-open state binds with low micromolar affinity to the cyclase subunit and displaces its substrate from the active site. As a consequence, catalysis by HisF and allosteric stimulation of HisH are impaired. Following UV-light irradiation, the DTE ligand adopts its ring-closed state and loses affinity for HisF, restoring activity and allostery. Our approach allows for the switching of ImGP-S activity and allostery during catalysis and appears to be generally applicable for the light regulation of other multienzyme complexes.

First insights of peptidoglycan amidation in Gram-positive bacteria - the high-resolution crystal structure of Staphylococcus aureus glutamine amidotransferase GatD

Gram-positive bacteria homeostasis and antibiotic resistance mechanisms are dependent on the intricate architecture of the cell wall, where amidated peptidoglycan plays an important role. The amidation reaction is carried out by the bi-enzymatic complex MurT-GatD, for which biochemical and structural information is very scarce. In this work, we report the first crystal structure of the glutamine amidotransferase member of this complex, GatD from Staphylococcus aureus, at 1.85 Å resolution. A glutamine molecule is found close to the active site funnel, hydrogen-bonded to the conserved R128. In vitro functional studies using 1H-NMR spectroscopy showed that S. aureus MurT-GatD complex has glutaminase activity even in the absence of lipid II, the MurT substrate. In addition, we produced R128A, C94A and H189A mutants, which were totally inactive for glutamine deamidation, revealing their essential role in substrate sequestration and catalytic reaction. GatD from S. aureus and other pathogenic bacteria share high identity to enzymes involved in cobalamin biosynthesis, which can be grouped in a new sub-family of glutamine amidotransferases. Given the ubiquitous presence of GatD, these results provide significant insights into the molecular basis of the so far undisclosed amidation mechanism, contributing to the development of alternative therapeutics to fight infections.

Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking

Understanding protein-protein interactions (PPIs) is quite important to elucidate crucial biological processes and even design compounds that interfere with PPIs with pharmaceutical significance. Protein-protein docking can afford the atomic structural details of protein-protein complexes, but the accurate prediction of the three-dimensional structures for protein-protein systems is still notoriously difficult due in part to the lack of an ideal scoring function for protein-protein docking. Compared with most scoring functions used in protein-protein docking, the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) methodologies are more theoretically rigorous, but their overall performance for the predictions of binding affinities and binding poses for protein-protein systems has not been systematically evaluated. In this study, we first evaluated the performance of MM/PBSA and MM/GBSA to predict the binding affinities for 46 protein-protein complexes. On the whole, different force fields, solvation models, and interior dielectric constants have obvious impacts on the prediction accuracy of MM/GBSA and MM/PBSA. The MM/GBSA calculations based on the ff02 force field, the GB model developed by Onufriev et al. and a low interior dielectric constant (εin = 1) yield the best correlation between the predicted binding affinities and the experimental data (rp = -0.647), which is better than MM/PBSA (rp = -0.523) and a number of empirical scoring functions used in protein-protein docking (rp = -0.141 to -0.529). Then, we examined the capability of MM/GBSA to identify the possible near-native binding structures from the decoys generated by ZDOCK for 43 protein-protein systems. The results illustrate that the MM/GBSA rescoring has better capability to distinguish the correct binding structures from the decoys than the ZDOCK scoring. Besides, the optimal interior dielectric constant of MM/GBSA for re-ranking docking poses may be determined by analyzing the characteristics of protein-protein binding interfaces. Considering the relatively high prediction accuracy and low computational cost, MM/GBSA may be a good choice for predicting the binding affinities and identifying correct binding structures for protein-protein systems.

Correlation of fitness landscapes from three orthologous TIM barrels originates from sequence and structure constraints

Sequence divergence of orthologous proteins enables adaptation to environmental stresses and promotes evolution of novel functions. Limits on evolution imposed by constraints on sequence and structure were explored using a model TIM barrel protein, indole-3-glycerol phosphate synthase (IGPS). Fitness effects of point mutations in three phylogenetically divergent IGPS proteins during adaptation to temperature stress were probed by auxotrophic complementation of yeast with prokaryotic, thermophilic IGPS. Analysis of beneficial mutations pointed to an unexpected, long-range allosteric pathway towards the active site of the protein. Significant correlations between the fitness landscapes of distant orthologues implicate both sequence and structure as primary forces in defining the TIM barrel fitness landscape and suggest that fitness landscapes can be translocated in sequence space. Exploration of fitness landscapes in the context of a protein fold provides a strategy for elucidating the sequence-structure-fitness relationships in other common motifs.

The TIM barrel fold is an evolutionarily conserved motif found in proteins with a variety of enzymatic functions. Here the authors explore the fitness landscape of the TIM barrel protein IGPS and uncover evolutionary constraints on both sequence and structure, accompanied by long range allosteric interactions.

Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex

It is important to identify hotspot residues that determine protein-protein interactions in interfaces of macromolecular complexes. We have applied a combination of ancestral sequence reconstruction and protein design to identify hotspots within imidazole glycerol phosphate synthase (ImGPS). ImGPS is a key metabolic enzyme complex, which links histidine and de novo purine biosynthesis and consists of the cyclase subunit HisF and the glutaminase subunit HisH. Initial fluorescence titration experiments showed that HisH from Zymomonas mobilis (zmHisH) binds with high affinity to the reconstructed HisF from the last universal common ancestor (LUCA-HisF) but not to HisF from Pyrobaculum arsenaticum (paHisF), which differ by 103 residues. Subsequent titration experiments with a reconstructed evolutionary intermediate linking LUCA-HisF and paHisF and inspection of the subunit interface of a contemporary ImGPS allowed us to narrow down the differences crucial for zmHisH binding to nine amino acids of HisF. Homology modeling and in silico mutagenesis studies suggested that at most two of these nine HisF residues are crucial for zmHisH binding. These computational results were verified by experimental site-directed mutagenesis, which finally enabled us to pinpoint a single amino acid residue in HisF that is decisive for high-affinity binding of zmHisH. Our work shows that the identification of protein interface hotspots can be very efficient when reconstructed proteins with different binding properties are included in the analysis. Proteins 2017; 85:312-321. © 2016 Wiley Periodicals, Inc.

Allosteric Communication Disrupted by a Small Molecule Binding to the Imidazole Glycerol Phosphate Synthase Protein-Protein Interface

Allosteric enzymes regulate a wide range of catalytic transformations, including biosynthetic mechanisms of important human pathogens, upon binding of substrate molecules to an orthosteric (or active) site and effector ligands at distant (allosteric) sites. We find that enzymatic activity can be impaired by small molecules that bind along the allosteric pathway connecting the orthosteric and allosteric sites, without competing with endogenous ligands. Noncompetitive allosteric inhibitors disrupted allostery in the imidazole glycerol phosphate synthase (IGPS) enzyme from Thermotoga maritima as evidenced by nuclear magnetic resonance, microsecond time-scale molecular dynamics simulations, isothermal titration calorimetry, and kinetic assays. The findings are particularly relevant for the development of allosteric antibiotics, herbicides, and antifungal compounds because IGPS is absent in mammals but provides an entry point to fundamental biosynthetic pathways in plants, fungi, and bacteria.

Dissecting Dynamic Allosteric Pathways Using Chemically Related Small-Molecule Activators

The allosteric mechanism of the heterodimeric enzyme imidazole glycerol phosphate synthase was studied in detail with solution nuclear magnetic resonance spectroscopy and molecular dynamics simulations. We studied IGPS in complex with a series of allosteric activators corresponding to a large range of catalytic rate enhancements (26- to 4,900-fold), in which ligand binding is entropically driven. Conformational flexibility on the millisecond timescale plays a crucial role in intersubunit communication. Carr-Purcell-Meiboom-Gill relaxation dispersion experiments probing Ile, Leu, and Val methyl groups reveal that the apo- and glutamine-mimicked complexes are static on the millisecond timescale. Domain-wide motions are stimulated in the presence of the allosteric activators. These studies, in conjunction with ligand titrations, demonstrate that the allosteric network is widely dispersed and varies with the identity of the effector. Furthermore, we find that stronger allosteric ligands create more conformational flexibility on the millisecond timescale throughout HisF. This domain-wide loosening leads to maximum catalytic activity.

Simulation Studies of Protein and Small Molecule Interactions and Reaction

Computational studies of protein and small molecule (protein-ligand/enzyme-substrate) interactions become more and more important in biological science and drug discovery. Computer modeling can provide molecular details of the processes such as conformational change, binding, and transportation of small molecules/proteins, which are not easily to be captured in experiments. In this chapter, we discussed simulation studies of both protein and small molecules from three aspects: conformation sampling, transportations of small molecules in enzymes, and enzymatic reactions involving small molecules. Both methodology developments and examples of simulation studies in this field were presented.

Solution NMR Spectroscopy for the Study of Enzyme Allostery

Allostery is a ubiquitous biological regulatory process in which distant binding sites within a protein or enzyme are functionally and thermodynamically coupled. Allosteric interactions play essential roles in many enzymological mechanisms, often facilitating formation of enzyme-substrate complexes and/or product release. Thus, elucidating the forces that drive allostery is critical to understanding the complex transformations of biomolecules. Currently, a number of models exist to describe allosteric behavior, taking into account energetics as well as conformational rearrangements and fluctuations. In the following Review, we discuss the use of solution NMR techniques designed to probe allosteric mechanisms in enzymes. NMR spectroscopy is unequaled in its ability to detect structural and dynamical changes in biomolecules, and the case studies presented herein demonstrate the range of insights to be gained from this valuable method. We also provide a detailed technical discussion of several specialized NMR experiments that are ideally suited for the study of enzymatic allostery.

The RNA-Binding Chaperone Hfq Is an Important Global Regulator of Gene Expression in Pasteurella multocida and Plays a Crucial Role in Production of a Number of Virulence Factors, Including Hyaluronic Acid Capsule

The Gram-negative bacterium Pasteurella multocida is the causative agent of a number of economically important animal diseases, including avian fowl cholera. Numerous P. multocida virulence factors have been identified, including capsule, lipopolysaccharide (LPS), and filamentous hemagglutinin, but little is known about how the expression of these virulence factors is regulated. Hfq is an RNA-binding protein that facilitates riboregulation via interaction with small noncoding RNA (sRNA) molecules and their mRNA targets. Here, we show that a P. multocida hfq mutant produces significantly less hyaluronic acid capsule during all growth phases and displays reduced in vivo fitness. Transcriptional and proteomic analyses of the hfq mutant during mid-exponential-phase growth revealed altered transcript levels for 128 genes and altered protein levels for 78 proteins. Further proteomic analyses of the hfq mutant during the early exponential growth phase identified 106 proteins that were produced at altered levels. Both the transcript and protein levels for genes/proteins involved in capsule biosynthesis were reduced in the hfq mutant, as were the levels of the filamentous hemagglutinin protein PfhB2 and its secretion partner LspB2. In contrast, there were increased expression levels of three LPS biosynthesis genes, encoding proteins involved in phosphocholine and phosphoethanolamine addition to LPS, suggesting that these genes are negatively regulated by Hfq-dependent mechanisms. Taken together, these data provide the first evidence that Hfq plays a crucial role in regulating the global expression of P. multocida genes, including the regulation of key P. multocida virulence factors, capsule, LPS, and filamentous hemagglutinin.

Protein Allostery and Conformational Dynamics

The functions of many proteins are regulated through allostery, whereby effector binding at a distal site changes the functional activity (e.g., substrate binding affinity or catalytic efficiency) at the active site. Most allosteric studies have focused on thermodynamic properties, in particular, substrate binding affinity. Changes in substrate binding affinity by allosteric effectors have generally been thought to be mediated by conformational transitions of the proteins or, alternatively, by changes in the broadness of the free energy basin of the protein conformational state without shifting the basin minimum position. When effector binding changes the free energy landscape of a protein in conformational space, the change affects not only thermodynamic properties but also dynamic properties, including the amplitudes of motions on different time scales and rates of conformational transitions. Here we assess the roles of conformational dynamics in allosteric regulation. Two cases are highlighted where NMR spectroscopy and molecular dynamics simulation have been used as complementary approaches to identify residues possibly involved in allosteric communication. Perspectives on contentious issues, for example, the relationship between picosecond-nanosecond local and microsecond-millisecond conformational exchange dynamics, are presented.

Case study on the evolution of hetero-oligomer interfaces based on the differences in paralogous proteins

We addressed the evolutionary trace of hetero-oligomer interfaces by comparing the structures of paralogous proteins; one of them is a monomer or homo-oligomer and the other is a hetero-oligomer. We found different trends in amino acid conservation pattern and hydrophobicity between homo-oligomer and hetero-oligomer. The degree of amino acid conservation in the interface of homo-oligomer has no obvious difference from that in the surface, whereas the degree of conservation is much higher in the interface of hetero-oligomer. The interface of homo-oligomer has a few very conserved residue positions, whereas the residue conservation in the interface of hetero-oligomer tends to be higher. In addition, the interface of hetero-oligomer has a tendency of being more hydrophobic compared with the one in homo-oligomer. We conjecture that these differences are related to the inherent symmetry in homo-oligomers that cannot exist in hetero-oligomers. Paucity of the structural data precludes statistical tests of these tendencies, yet the trend can be applied to the prediction of the interface of hetero-oligomer. We obtained putative interfaces of the subunits in CPSF (cleavage and polyadenylation specificity factor), one of the human pre-mRNA 3′-processing complexes. The locations of predicted interface residues were consistent with the known experimental data.

DynaFace: Discrimination between Obligatory and Non-obligatory Protein-Protein Interactions Based on the Complex's Dynamics

Protein-protein interfaces have been evolutionarily-designed to enable transduction between the interacting proteins. Thus, we hypothesize that analysis of the dynamics of the complex can reveal details about the nature of the interaction, and in particular whether it is obligatory, i.e., persists throughout the entire lifetime of the proteins, or not. Indeed, normal mode analysis, using the Gaussian network model, shows that for the most part obligatory and non-obligatory complexes differ in their decomposition into dynamic domains, i.e., the mobile elements of the protein complex. The dynamic domains of obligatory complexes often mix segments from the interacting chains, and the hinges between them do not overlap with the interface between the chains. In contrast, in non-obligatory complexes the interface often hinges between dynamic domains, held together through few anchor residues on one side of the interface that interact with their counterpart grooves in the other end. In automatic analysis, 117 of 139 obligatory (84.2%) and 203 of 246 non-obligatory (82.5%) complexes are correctly classified by our method: DynaFace. We further use DynaFace to predict obligatory and non-obligatory interactions among a set of 300 putative protein complexes. DynaFace is available at: http://safir.prc.boun.edu.tr/dynaface.

Biosynthesis of Histidine

The biosynthesis of histidine in Escherichia coli and Salmonella typhimurium has been an important model system for the study of relationships between the flow of intermediates through a biosynthetic pathway and the control of the genes encoding the enzymes that catalyze the steps in a pathway. This article provides a comprehensive review of the histidine biosynthetic pathway and enzymes, including regulation of the flow of intermediates through the pathway and mechanisms that regulate the amounts of the histidine biosynthetic enzymes. In addition, this article reviews the structure and regulation of the histidine (his) biosynthetic operon, including transcript processing, Rho-factor-dependent "classical" polarity, and the current model of his operon attenuation control. Emphasis is placed on areas of recent progress. Notably, most of the enzymes that catalyze histidine biosynthesis have recently been crystallized, and their structures have been determined. Many of the histidine biosynthetic intermediates are unstable, and the histidine biosynthetic enzymes catalyze some chemically unusual reactions. Therefore, these studies have led to considerable mechanistic insight into the pathway itself and have provided deep biochemical understanding of several fundamental processes, such as feedback control, allosteric interactions, and metabolite channeling. Considerable recent progress has also been made on aspects of his operon regulation, including the mechanism of pp(p)Gpp stimulation of his operon transcription, the molecular basis for transcriptional pausing by RNA polymerase, and pathway evolution. The progress in these areas will continue as sophisticated new genomic, metabolomic, proteomic, and structural approaches converge in studies of the histidine biosynthetic pathway and mechanisms of control of his biosynthetic genes in other bacterial species.

Crystal structures capture three states in the catalytic cycle of a pyridoxal phosphate (PLP) synthase

PLP synthase (PLPS) is a remarkable single-enzyme biosynthetic pathway that produces pyridoxal 5'-phosphate (PLP) from glutamine, ribose 5-phosphate, and glyceraldehyde 3-phosphate. The intact enzyme includes 12 synthase and 12 glutaminase subunits. PLP synthesis occurs in the synthase active site by a complicated mechanism involving at least two covalent intermediates at a catalytic lysine. The first intermediate forms with ribose 5-phosphate. The glutaminase subunit is a glutamine amidotransferase that hydrolyzes glutamine and channels ammonia to the synthase active site. Ammonia attack on the first covalent intermediate forms the second intermediate. Glyceraldehyde 3-phosphate reacts with the second intermediate to form PLP. To investigate the mechanism of the synthase subunit, crystal structures were obtained for three intermediate states of the Geobacillus stearothermophilus intact PLPS or its synthase subunit. The structures capture the synthase active site at three distinct steps in its complicated catalytic cycle, provide insights into the elusive mechanism, and illustrate the coordinated motions within the synthase subunit that separate the catalytic states. In the intact PLPS with a Michaelis-like intermediate in the glutaminase active site, the first covalent intermediate of the synthase is fully sequestered within the enzyme by the ordering of a generally disordered 20-residue C-terminal tail. Following addition of ammonia, the synthase active site opens and admits the Lys-149 side chain, which participates in formation of the second intermediate and PLP. Roles are identified for conserved Asp-24 in the formation of the first intermediate and for conserved Arg-147 in the conversion of the first to the second intermediate.

Weighted Implementation of Suboptimal Paths (WISP): An Optimized Algorithm and Tool for Dynamical Network Analysis

Allostery can occur by way of subtle cooperation among protein residues (e.g., amino acids) even in the absence of large conformational shifts. Dynamical network analysis has been used to model this cooperation, helping to computationally explain how binding to an allosteric site can impact the behavior of a primary site many ångstroms away. Traditionally, computational efforts have focused on the most optimal path of correlated motions leading from the allosteric to the primary active site. We present a program called Weighted Implementation of Suboptimal Paths (WISP) capable of rapidly identifying additional suboptimal pathways that may also play important roles in the transmission of allosteric signals. Aside from providing signal redundancy, suboptimal paths traverse residues that, if disrupted through pharmacological or mutational means, could modulate the allosteric regulation of important drug targets. To demonstrate the utility of our program, we present a case study describing the allostery of HisH-HisF, an amidotransferase from T. maritima thermotiga. WISP and its VMD-based graphical user interface (GUI) can be downloaded from http://nbcr.ucsd.edu/wisp.

A general method for artificial metalloenzyme formation through strain-promoted azide-alkyne cycloaddition

Strain-promoted azide-alkyne cycloaddition (SPAAC) can be used to generate artificial metalloenzymes (ArMs) from scaffold proteins containing a p-azido-L-phenylalanine (Az) residue and catalytically active bicyclononyne-substituted metal complexes. The high efficiency of this reaction allows rapid ArM formation when using Az residues within the scaffold protein in the presence of cysteine residues or various reactive components of cellular lysate. In general, cofactor-based ArM formation allows the use of any desired metal complex to build unique inorganic protein materials. SPAAC covalent linkage further decouples the native function of the scaffold from the installation process because it is not affected by native amino acid residues; as long as an Az residue can be incorporated, an ArM can be generated. We have demonstrated the scope of this method with respect to both the scaffold and cofactor components and established that the dirhodium ArMs generated can catalyze the decomposition of diazo compounds and both Si-H and olefin insertion reactions involving these carbene precursors.

KBDOCK 2013: a spatial classification of 3D protein domain family interactions

Comparing, classifying and modelling protein structural interactions can enrich our understanding of many biomolecular processes. This contribution describes Kbdock (http://kbdock.loria.fr/), a database system that combines the Pfam domain classification with coordinate data from the PDB to analyse and model 3D domain–domain interactions (DDIs). Kbdock can be queried using Pfam domain identifiers, protein sequences or 3D protein structures. For a given query domain or pair of domains, Kbdock retrieves and displays a non-redundant list of homologous DDIs or domain–peptide interactions in a common coordinate frame. Kbdock may also be used to search for and visualize interactions involving different, but structurally similar, Pfam families. Thus, structural DDI templates may be proposed even when there is little or no sequence similarity to the query domains.

Importance of hydrophobic cavities in allosteric regulation of formylglycinamide synthetase: insight from xenon trapping and statistical coupling analysis

Formylglycinamide ribonucleotide amidotransferase (FGAR-AT) is a 140 kDa bi-functional enzyme involved in a coupled reaction, where the glutaminase active site produces ammonia that is subsequently utilized to convert FGAR to its corresponding amidine in an ATP assisted fashion. The structure of FGAR-AT has been previously determined in an inactive state and the mechanism of activation remains largely unknown. In the current study, hydrophobic cavities were used as markers to identify regions involved in domain movements that facilitate catalytic coupling and subsequent activation of the enzyme. Three internal hydrophobic cavities were located by xenon trapping experiments on FGAR-AT crystals and further, these cavities were perturbed via site-directed mutagenesis. Biophysical characterization of the mutants demonstrated that two of these three voids are crucial for stability and function of the protein, although being ∼20 Å from the active centers. Interestingly, correlation analysis corroborated the experimental findings, and revealed that amino acids lining the functionally important cavities form correlated sets (co-evolving residues) that connect these regions to the amidotransferase active center. It was further proposed that the first cavity is transient and allows for breathing motion to occur and thereby serves as an allosteric hotspot. In contrast, the third cavity which lacks correlated residues was found to be highly plastic and accommodated steric congestion by local adjustment of the structure without affecting either stability or activity.

Solution NMR and computational methods for understanding protein allostery

Allosterism is an essential biological regulatory mechanism. In enzymes, allosteric regulation results in an activation or inhibition of catalytic turnover. The mechanisms by which this is accomplished are unclear and vary significantly depending on the enzyme. It is commonly the case that a metabolite binds to the enzyme at a site distant from the catalytic site, yet its binding is coupled to and sensed by the active site. This coupling can manifest in changes in structure, dynamics, or both at the active site. These interactions between the allosteric and active site, which are often quite distant from one another, involve numerous atoms as well as complex conformational rearrangements of the protein secondary and tertiary structure. Interrogation of this complex biological phenomenon necessitates multiple experimental approaches. In this article, we outline a combined solution NMR spectroscopic and computational approach using molecular dynamics and network models to uncover mechanistic aspects of allostery in the enzyme imidazole glycerol phosphate synthase.

Exploring residue component contributions to dynamical network models of allostery

Allosteric regulation in biological systems is of considerable interest given the vast number of proteins that exhibit such behavior. Network models obtained from molecular dynamics simulations have been shown to be powerful tools for the analysis of allostery. In this work, different coarse-grain residue representations (nodes) are used together with a dynamical network model to investigate models of allosteric regulation. This model assumes that allosteric signals are dependent on positional correlations of protein substituents, as determined through molecular dynamics simulations, and uses correlated motion to generate a signaling weight between two given nodes. We examine four types of network models using different node representations in Cartesian coordinates: the (i) residue alpha-carbons, (ii) sidechain center of mass, (iii) backbone center of mass, and the entire (iv) residue center of mass. All correlations are filtered by a dynamic contact map that defines the allowable interactions between nodes based on physical proximity. We apply the four models to imidazole glycerol phosphate synthase (IGPS), which provides a well-studied experimental framework in which allosteric communication is known to persist across disparate protein domains (e.g. a protein dimer interface). IGPS is modeled as a network of nodes and weighted edges. Optimal allosteric pathways are traced using the Floyd Warshall algorithm for weighted networks, and community analysis (a form of hierarchical clustering) is performed using the Girvan-Newman algorithm. Our results show that dynamical information encoded in the residue center of mass must be included in order to detect residues that are experimentally known to play a role in allosteric communication for IGPS. More broadly, this new method may be useful for predicting pathways of allosteric communication for any biomolecular system in atomic detail.