The Allosteric Mechanism of Yeast Chorismate Mutase: A Dynamic Analysis

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.

Exploring movement and energy in human P-glycoprotein conformational rearrangement

Human P-glycoprotein (P-gp), a kind of ATP-Binding Cassette transporter, can export a diverse variety of anti-cancer drugs out of the tumor cell. Its overexpression is one of the main reasons for the multidrug resistance (MDR) of tumor cells. It has been confirmed that during the substrate transport process, P-gp experiences a large-scale structural rearrangement from the inward- to outward-facing states. However, the mechanism of how the nucleotide-binding domains (NBDs) control the transmembrane domains (TMDs) to open towards the periplasm in the outward-facing state has not yet been fully characterized. Herein, targeted molecular dynamics simulations were performed to explore the conformational rearrangement of human P-gp. The results show that the allosteric process proceeds in a coupled way, and first the transition is driven by the NBDs, and then transmitted to the cytoplasmic parts of TMDs, finally to the periplasmic parts. The trajectories show that besides the translational motions, the NBDs undergo a rotation movement, which mainly occurs in xy plane and ensures the formation of the correct ATP-binding pockets. The analyses on the interaction energies between the six structure segments (cICLs) from the TMDs and NBDs reveal that their subtle energy differences play an important role in causing the periplasmic parts of the transmembrane helices to separate from each other in the established directions and in appropriate amplitudes. This conclusion can explain the two experimental phenomena about human P-gp in some extent. These studies have provided a detailed exploration into human P-gp rearrangement process and given an energy insight into the TMD reorientation during P-gp transition.

Evolution of allosteric regulation in chorismate mutases from early plants

Plants, fungi, and bacteria synthesize the aromatic amino acids: l-phenylalanine, l-tyrosine, and l-tryptophan. Chorismate mutase catalyzes the branch point reaction of phenylalanine and tyrosine biosynthesis to generate prephenate. In Arabidopsis thaliana, there are two plastid-localized chorismate mutases that are allosterically regulated (AtCM1 and AtCM3) and one cytosolic isoform (AtCM2) that is unregulated. Previous analysis of plant chorismate mutases suggested that the enzymes from early plants (i.e. bryophytes/moss, lycophytes, and basal angiosperms) formed a clade distinct from the isoforms found in flowering plants; however, no biochemical information on these enzymes is available. To understand the evolution of allosteric regulation in plant chorismate mutases, we analyzed a basal lineage of plant enzymes homologous to AtCM1 based on sequence similarity. The chorismate mutases from the moss/bryophyte Physcomitrella patens (PpCM1 and PpCM2), the lycophyte Selaginella moellendorffii (SmCM), and the basal angiosperm Amborella trichopoda (AmtCM1 and AmtCM2) were characterized biochemically. Tryptophan was a positive effector for each of the five enzymes examined. Histidine was a weak positive effector for PpCM1 and AmtCM1. Neither tyrosine nor phenylalanine altered the activity of SmCM; however, tyrosine was a negative regulator of the other four enzymes. Phenylalanine down-regulates both moss enzymes and AmtCM2. The 2.0 Å X-ray crystal structure of PpCM1 in complex with the tryptophan identified the allosteric effector site and reveals structural differences between the R- (more active) and T-state (less active) forms of plant chorismate mutases. Molecular insight into the basal plant chorismate mutases guides our understanding of the evolution of allosteric regulation in these enzymes.

The "Hot-Solvent/Cold-Solute" Problem Revisited

The temperature steers the equilibrium and nonequilibrium conformational dynamics of macromolecules in solution. Therefore, corresponding molecular dynamics simulations require a strategy for temperature control which should guarantee that the experimental statistical ensemble is also sampled in silico. Several algorithms for temperature control have been proposed. All these thermostats interfere with the macromolecule's "natural" dynamics as given by the Newtonian mechanics. Furthermore, using a single thermostat for an inhomogeneous solute-solvent system can lead to stationary temperature gradients. To avoid this "hot solvent/cold solute" problem, two separate thermostats are frequently applied, one to the solute and one to the solvent. However, such a separate temperature control will perturb the dynamics of the macromolecule much more strongly than a global one and, therefore, can introduce large artifacts into its conformational dynamics. Based on the concept that an explicit solvent environment represents an ideal thermostat concerning the magnitude and time correlation of temperature fluctuations of the solute, we propose a temperature control strategy that, on the one hand, provides a homogeneous temperature distribution throughout the system together with the correct statistical ensemble for the solute molecule while, on the other hand, minimally perturbing its dynamics.

Human Inducible Hsp70: Structures, Dynamics, and Interdomain Communication from All-Atom Molecular Dynamics Simulations

The 70 kDa human heat shock protein is a major molecular chaperone involved in de novo folding of proteins in vivo and refolding of proteins under stress conditions. Hsp70 is related to several "misfolding diseases" and other major pathologies, such as cancer, and is a target for new therapies. Hsp70 is comprised of two main domains: an N-terminal nucleotide binding domain (NBD) and a C-terminal substrate protein binding domain (SBD). The chaperone function of Hsp70 is based on an allosteric mechanism. Binding of ATP in NBD decreases the affinity of the substrate for SBD, and hydrolysis of ATP is promoted by binding of polypeptide segments in the SBD. No complete structure of human Hsp70 is known. Here, we report two models of human Hsp70, constructed by homology with Saccharomyces cerevisiae cochaperone protein Hsp110 (open model) and with Escherichia coli 70 kDa DnaK (closed model) and relaxed for several tens to hundreds of nanoseconds by using all-atom molecular dynamics simulations in explicit solvent. We obtain two stable states, Hsp70 with SBD open and SBD closed, which agree with experimental and structural information for ATP-Hsp70 and ADP-Hsp70, respectively. The dynamics of the transition from the open to closed states is investigated with a coarse-grained model and normal-mode analysis. The results show that the conformational change between the two states can be represented by a relatively small number of collective modes which involved major conformational changes in the two domains. These modes provide a mechanistic representation of the communication between NBD and SBD and allow us to identify subdomains and residues that appear to have a critical role in the conformational change mechanism that guides the chaperoning cycle of Hsp70.

Conversion of aminodeoxychorismate synthase into anthranilate synthase with Janus mutations: mechanism of pyruvate elimination catalyzed by chorismate enzymes

The central importance of chorismate enzymes in bacteria, fungi, parasites, and plants combined with their absence in mammals makes them attractive targets for antimicrobials and herbicides. Two of these enzymes, anthranilate synthase (AS) and aminodeoxychorismate synthase (ADCS), are structurally and mechanistically similar. The first catalytic step, amination at C2, is common between them, but AS additionally catalyzes pyruvate elimination, aromatizing the aminated intermediate to anthranilate. Despite prior attempts, the conversion of a pyruvate elimination-deficient enzyme into an elimination-proficient one has not been reported. Janus, a bioinformatics method for predicting mutations required to functionally interconvert homologous enzymes, was employed to predict mutations to convert ADCS into AS. A genetic selection on a library of Janus-predicted mutations was performed. Complementation of an AS-deficient strain of Escherichia coli grown on minimal medium led to several ADCS mutants that allow growth in 6 days compared to 2 days for wild-type AS. The purified mutant enzymes catalyze the conversion of chorismate to anthranilate at rates that are ∼50% of the rate of wild-type ADCS-catalyzed conversion of chorismate to aminodeoxychorismate. The residues mutated do not contact the substrate. Molecular dynamics studies suggest that pyruvate elimination is controlled by the conformation of the C2-aminated intermediate. Enzymes that catalyze elimination favor the equatorial conformation, which presents the C2-H to a conserved active site lysine (Lys424) for deprotonation and maximizes stereoelectronic activation. Acid/base catalysis of pyruvate elimination was confirmed in AS and salicylate synthase by showing incorporation of a solvent-derived proton into the pyruvate methyl group and by solvent kinetic isotope effects on pyruvate elimination catalyzed by AS.

Structural evolution of differential amino acid effector regulation in plant chorismate mutases

Chorismate mutase converts chorismate into prephenate for aromatic amino acid biosynthesis. To understand the molecular basis of allosteric regulation in the plant chorismate mutases, we analyzed the three Arabidopsis thaliana chorismate mutase isoforms (AtCM1-3) and determined the x-ray crystal structures of AtCM1 in complex with phenylalanine and tyrosine. Functional analyses show a wider range of effector control in the Arabidopsis chorismate mutases than previously reported. AtCM1 is activated by tryptophan with phenylalanine and tyrosine acting as negative effectors; however, tryptophan, cysteine, and histidine activate AtCM3. AtCM2 is a nonallosteric form. The crystal structure of AtCM1 in complex with tyrosine and phenylalanine identifies differences in the effector sites of the allosterically regulated yeast enzyme and the other two Arabidopsis isoforms. Site-directed mutagenesis of residues in the effector site reveals key features leading to differential effector regulation in these enzymes. In AtCM1, mutations of Gly-213 abolish allosteric regulation, as observed in AtCM2. A second effector site position, Gly-149 in AtCM1 and Asp-132 in AtCM3, controls amino acid effector specificity in AtCM1 and AtCM3. Comparisons of chorismate mutases from multiple plants suggest that subtle differences in the effector site are conserved in different lineages and may lead to specialized regulation of this branch point enzyme.

A unified view of "how allostery works"

The question of how allostery works was posed almost 50 years ago. Since then it has been the focus of much effort. This is for two reasons: first, the intellectual curiosity of basic science and the desire to understand fundamental phenomena, and second, its vast practical importance. Allostery is at play in all processes in the living cell, and increasingly in drug discovery. Many models have been successfully formulated, and are able to describe allostery even in the absence of a detailed structural mechanism. However, conceptual schemes designed to qualitatively explain allosteric mechanisms usually lack a quantitative mathematical model, and are unable to link its thermodynamic and structural foundations. This hampers insight into oncogenic mutations in cancer progression and biased agonists' actions. Here, we describe how allostery works from three different standpoints: thermodynamics, free energy landscape of population shift, and structure; all with exactly the same allosteric descriptors. This results in a unified view which not only clarifies the elusive allosteric mechanism but also provides structural grasp of agonist-mediated signaling pathways, and guides allosteric drug discovery. Of note, the unified view reasons that allosteric coupling (or communication) does not determine the allosteric efficacy; however, a communication channel is what makes potential binding sites allosteric.

Conformational transition pathway in the activation process of allosteric glucokinase

Glucokinase (GK) is a glycolytic enzyme that plays an important role in regulating blood glucose level, thus acting as a potentially attractive target for drug discovery in the treatment of diabetes of the young type 2 and persistent hyperinsulinemic hypoglycemia of infancy. To characterize the activation mechanism of GK from the super-open state (inactive state) to the closed state (active state), a series of conventional molecular dynamics (MD) and targeted MD (TMD) simulations were performed on this enzyme. Conventional MD simulation showed a specific conformational ensemble of GK when the enzyme is inactive. Seven TMD simulations depicted a reliably conformational transition pathway of GK from the inactive state to the active state, and the components important to the conformational change of GK were identified by analyzing the detailed structures of the TMD trajectories. In combination with the inactivation process, our findings showed that the whole conformational pathway for the activation-inactivation-activation of GK is a one-direction circulation, and the active state is less stable than the inactive state in the circulation. Additionally, glucose was demonstrated to gradually modulate its binding pose with the help of residues in the large domain and connecting region of GK during the activation process. Furthermore, the obtained energy barriers were used to explain the preexisting equilibrium and the slow binding kinetic process of the substrate by GK. The simulated results are in accordance with the recent findings from the mutagenesis experiments and kinetic analyses. Our observations reveal a complicated conformational process in the allosteric protein, resulting in new knowledge about the delicate mechanisms for allosteric biological macromolecules that will be useful in drug design for targeting allosteric proteins.

A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35 V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35 V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.

Despite recent progress in computational and experimental studies of dynamic regulation in protein kinases, a mechanistic understanding of long-range communication and mechanisms of mutation-induced signaling controlling kinase activity remains largely qualitative. In this study, we have performed a systematic modeling and analysis of allosteric activation in ABL and EGFR kinases at the increasing level of complexity - from catalytic domain to multi-domain regulatory complexes. The results of this study have revealed organizing structural and mechanistic principles of allosteric signaling in protein kinases. Although activation mechanisms in ABL and EGFR kinases have evolved through acquisition of structurally different regulatory complexes, we have found that long-range interdomain communication between common functional segments (αF-helix and αC-helix) may be important for allosteric activation. The results of study have revealed molecular signatures of activating cancer mutations and have shed the light on general mechanistic aspects of mutation-induced signaling in protein kinases. An advanced understanding and further characterization of molecular signatures of kinase mutations may aid in a better rationalization of mutational effects on clinical outcomes and facilitate molecular-based therapeutic strategies to combat kinase mutation-dependent tumorigenesis.

The energy landscape analysis of cancer mutations in protein kinases

The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling may involve a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. The presented study has demonstrated that activation cancer mutations may affect the thermodynamic equilibrium between kinase states by allosterically altering the distribution of locally frustrated sites and increasing the local frustration in the inactive form, while eliminating locally frustrated sites and restoring structural rigidity of the active form. The energy landsape analysis of protein kinases and the proposed role of locally frustrated sites in activation mechanisms may have useful implications for bioinformatics-based screening and detection of functional sites critical for allosteric regulation in complex biomolecular systems.

Coarse-grained representation of protein flexibility. Foundations, successes, and shortcomings

Flexibility is the key magnitude to understand the variety of functions of proteins. Unfortunately, its experimental study is quite difficult, and in fact, most experimental procedures are designed to reduce flexibility and allow a better definition of the structure. Theoretical approaches have become then the alternative but face serious timescale problems, since many biologically relevant deformation movements happen in a timescale that is far beyond the possibility of current atomistic models. In this complex scenario, coarse-grained simulation methods have emerged as a powerful and inexpensive alternative. Along this chapter, we will review these coarse-grained methods, and explain their physical foundations and their range of applicability.

Conformational sampling, catalysis, and evolution of the bacterial phosphotriesterase

To efficiently catalyze a chemical reaction, enzymes are required to maintain fast rates for formation of the Michaelis complex, the chemical reaction and product release. These distinct demands could be satisfied via fluctuation between different conformational substates (CSs) with unique configurations and catalytic properties. However, there is debate as to how these rapid conformational changes, or dynamics, exactly affect catalysis. As a model system, we have studied bacterial phosphotriesterase (PTE), which catalyzes the hydrolysis of the pesticide paraoxon at rates limited by a physical barrier-either substrate diffusion or conformational change. The mechanism of paraoxon hydrolysis is understood in detail and is based on a single, dominant, enzyme conformation. However, the other aspects of substrate turnover (substrate binding and product release), although possibly rate-limiting, have received relatively little attention. This work identifies "open" and "closed" CSs in PTE and dominant structural transition in the enzyme that links them. The closed state is optimally preorganized for paraoxon hydrolysis, but seems to block access to/from the active site. In contrast, the open CS enables access to the active site but is poorly organized for hydrolysis. Analysis of the structural and kinetic effects of mutations distant from the active site suggests that remote mutations affect the turnover rate by altering the conformational landscape.

Evolutionary basis for the coupled-domain motions in Thermus thermophilus leucyl-tRNA synthetase

Aminoacyl-tRNA synthetases are multidomain proteins that catalyze the covalent attachment of amino acids to their cognate transfer RNA. Various domains of an aminoacyl-tRNA synthetase perform their specific functions in a highly coordinated manner to maintain high accuracy in protein synthesis in cells. The coordination of their function, therefore, requires communication between domains. In this study we explored the relevance of enzyme motion in domain-domain communications. Specifically, we attempted to probe whether the communication between distantly located domains of a multidomain protein is accomplished through a coordinated movement of structural elements. We investigated the collective motion in Thermus thermophilus leucyl-tRNA synthetase by studying the low frequency normal modes. We identified the mode that best described the experimentally observed conformational changes of T. thermophilus leucyl-tRNA synthetase upon substrate binding and analyzed the correlated and anticorrelated motions between different domains. Furthermore, we used statistical coupling analysis to explore if the amino acid pairs and/or clusters whose motions are thermally coupled have also coevolved. Our study demonstrates that a small number of residues belong to the category whose coupled thermal motions correspond to evolutionary coupling as well. These residue clusters constitute a distinguished set of interacting networks that are sparsely distributed in the domain interface. Residues of these networking clusters are within van der Waals contact, and we suggest that they are critical in the propagation of long range mechanochemical motions in T. thermophilus leucyl-tRNA synthetase.

Allosteric communication in myosin V: from small conformational changes to large directed movements

The rigor to post-rigor transition in myosin, a consequence of ATP binding, plays an essential role in the Lymn–Taylor functional cycle because it results in the dissociation of the actomyosin complex after the powerstroke. On the basis of the X-ray structures of myosin V, we have developed a new normal mode superposition model for the transition path between the two states. Rigid-body motions of the various subdomains and specific residues at the subdomain interfaces are key elements in the transition. The allosteric communication between the nucleotide binding site and the U50/L50 cleft is shown to result from local changes due to ATP binding, which induce large amplitude motions that are encoded in the structure of the protein. The triggering event is the change in the interaction of switch I and the P-loop, which is stabilized by ATP binding. The motion of switch I, which is a relatively rigid element of the U50 subdomain, leads directly to a partial opening of the U50/L50 cleft; the latter is expected to weaken the binding of myosin to actin. The calculated transition path demonstrates the nature of the subdomain coupling and offers an explanation for the mutual exclusion of ATP and actin binding. The mechanism of the uncoupling of the converter from the motor head, an essential part of the transition, is elucidated. The origin of the partial untwisting of the central β-sheet in the rigor to post-rigor transition is described.

Myosins are molecular motor proteins that interact with actin filaments to perform a wide range of cellular functions. They use the universal energy storage molecule adenosine triphosphate (ATP). The functional cycle involves myosin binding to actin, a “powerstroke” leading to directed movement, and myosin release in preparation for the next step. A fundamental question concerns the mechanism by which the local structural changes due to ATP binding, hydrolysis, and products release can generate the large myosin changes of conformation required for this cycle. Here, we focus on the rigor to post-rigor transition of myosin V, which results in the release of myosin from actin. Starting from the X-ray structures of the two states, we have used the optimal superposition of normal modes to determine the transition path. The path shows the allosteric mechanism by which ATP binding leads to the opening of the U50/L50 cleft, the essential step in the unbinding of myosin from actin. More generally, the new normal-mode superposition model can be useful for describing large-amplitude conformational transitions encoded in protein structures by evolution.

Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C(alpha)-root mean square deviations between two end structures of 7.13 A and 3.55 A, respectively). Target conformations are reached by root mean-square deviations of 2.27 A and 1.90 A for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-A root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.

Insights into protein flexibility: The relationship between normal modes and conformational change upon protein-protein docking

Understanding protein interactions has broad implications for the mechanism of recognition, protein design, and assigning putative functions to uncharacterized proteins. Studying protein flexibility is a key component in the challenge of describing protein interactions. In this work, we characterize the observed conformational change for a set of 20 proteins that undergo large conformational change upon association (>2 A Calpha RMSD) and ask what features of the motion are successfully reproduced by the normal modes of the system. We demonstrate that normal modes can be used to identify mobile regions and, in some proteins, to reproduce the direction of conformational change. In 35% of the proteins studied, a single low-frequency normal mode was found that describes well the direction of the observed conformational change. Finally, we find that for a set of 134 proteins from a docking benchmark that the characteristic frequencies of normal modes can be used to predict reliably the extent of observed conformational change. We discuss the implications of the results for the mechanics of protein recognition.

Direct observation in solution of a preexisting structural equilibrium for a mutant of the allosteric aspartate transcarbamoylase

Many signaling and metabolic pathways rely on the ability of some of the proteins involved to undergo a substrate-induced transition between at least two structural states. Among the various models put forward to account for binding and activity curves of those allosteric proteins, the Monod, Wyman, and Changeux model for allostery theory has certainly been the most influential, although a central postulate, the preexisting equilibrium between the low-activity, low-affinity quaternary structure and the high-activity, high-affinity quaternary structure states in the absence of substrates, has long awaited direct experimental substantiation. Upon substrate binding, allosteric Escherichia coli aspartate transcarbamoylase adopts alternate quaternary structures, stabilized by a set of interdomain and intersubunit interactions, which are readily differentiated by their solution x-ray scattering curves. Disruption of a salt link, which is observed only in the low-activity, low-affinity quaternary structure, between Lys-143 of the regulatory chain and Asp-236 of the catalytic chain yields a mutant enzyme that is in a reversible equilibrium between at least two states in the absence of ligand, a major tenet of the Monod, Wyman, and Changeux model. By using this mutant as a magnifying glass of the structural effect of ligand binding, a comparative analysis of the binding of carbamoyl phosphate (CP) and analogs points out the crucial role of the amine group of CP in facilitating the transition toward the high-activity, high-affinity quaternary state. Thus, the cooperative binding of aspartate in aspartate transcarbamoylase appears to result from the combination of the preexisting quaternary structure equilibrium with local changes induced by CP binding.

Conformational transition pathway in the allosteric process of human glucokinase

Glucokinase (GK) is an important enzyme for regulating blood glucose levels and a potentially attractive target for diabetes of the young type 2 and persistent hyperinsulinemic hypoglycemia of infancy. To characterize the conformational transition of GK from the closed state to the superopen state, a series of conventional molecular dynamics (MD) and target MD (TMD) simulations were performed on both the wild-type enzyme and its mutants. Two 10-ns conventional MD simulations showed that, although the allosteric site of GK is approximately 20 A away from the active site, the activator is able to enhance the activity of the enzyme through conformational restriction. Fourteen TMD simulations on GK and five of its mutants revealed a reliably conformational transition pathway. The overall conformational transition includes three stages, and three likely stable intermediate states were identified by free energy scanning for the snapshots throughout the pathway. The conformational transition feature revealed by our TMD simulations rationalized several important mutagenesis and kinetic data. Remarkably, the TMD simulations predicted that Y61S, I159A, A201R, V203E, and V452S mutations, which have not been investigated so far, may facilitate the opening process of GK. These predictions also have been verified by mutagenesis and kinetic analyses in this study. These observations are beneficial to understanding the mechanism of GK regulation and designing the compounds for treating metabolic diseases.