Investigation of the Interaction between the Large and Small Subunits of Potato ADP-Glucose Pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase), a key allosteric enzyme involved in higher plant starch biosynthesis, is composed of pairs of large (LS) and small subunits (SS). Current evidence indicates that the two subunit types play distinct roles in enzyme function. Recently the heterotetrameric structure of potato AGPase has been modeled. In the current study, we have applied the molecular mechanics generalized born surface area (MM-GBSA) method and identified critical amino acids of the potato AGPase LS and SS subunits that interact with each other during the native heterotetrameric structure formation. We have further shown the role of the LS amino acids in subunit-subunit interaction by yeast two-hybrid, bacterial complementation assay and native gel. Comparison of the computational results with the experiments has indicated that the backbone energy contribution (rather than the side chain energies) of the interface residues is more important in identifying critical residues. We have found that lateral interaction of the LS-SS is much stronger than the longitudinal one, and it is mainly mediated by hydrophobic interactions. This study will not only enhance our understanding of the interaction between the SS and the LS of AGPase, but will also enable us to engineer proteins to obtain better assembled variants of AGPase which can be used for the improvement of plant yield.     

Analysis of allosteric effector binding sites of potato ADP-glucose pyrophosphorylase through reverse genetics

ADP-glucose pyrophosphorylase (AGPase) is a key regulatory enzyme of bacterial glycogen and plant starch synthesis as it controls carbon flux via its allosteric regulatory behavior. Unlike the bacterial enzyme that is composed of a single subunit type, the plant AGPase is a heterotetrameric enzyme (alpha2beta2) with distinct roles for each subunit type. The large subunit (LS) is involved mainly in allosteric regulation through its interaction with the catalytic small subunit (SS). The LS modulates the catalytic activity of the SS by increasing the allosteric regulatory response of the hetero-oligomeric enzyme. To identify regions of the LS involved in binding of effector molecules, a reverse genetics approach was employed. A potato (Solanum tuberosum L.) AGPase LS down-regulatory mutant (E38A) was subjected to random mutagenesis using error-prone polymerase chain reaction and screened for the capacity to form an enzyme capable of restoring glycogen production in glgC(-) Escherichia coli. Dominant mutations were identified by their capacity to restore glycogen production when the LS containing only the second site mutations was co-expressed with the wild-type SS. Sequence analysis showed that most of the mutations were decidedly nonrandom and were clustered at conserved N- and C-terminal regions. Kinetic analysis of the dominant mutant enzymes indicated that the K(m) values for cofactor and substrates were comparable with the wild-type AGPase, whereas the affinities for activator and inhibitor were altered appreciably. These AGPase variants displayed increased resistance to P(i) inhibition and/or greater sensitivity toward 3-phosphoglyceric acid activation. Further studies of Lys-197, Pro-261, and Lys-420, residues conserved in AGPase sequences, by site-directed mutagenesis suggested that the effectors 3-phosphoglyceric acid and P(i) interact at two closely located binding sites.     

Insight into the 3D structure of ADP-glucose pyrophosphorylase from rice (Oryza sativa L.)

ADP-glucose pyrophosphorylase (E.C. 2.7.7.27; AGPase) is a key regulatory enzyme that catalyzes the rate-limiting step of starch biosynthesis in higher plants. AGPase consists of pair of small (SS) and large (LS) subunits thereby constituting a heterotetrameric structure. No crystal structure of the native heterotetrameric enzyme is available for any species, thus limiting the complete understanding of structure–function relationships of this enzyme. In this study, an attempt was made to deduce the heterotetrameric assembly of AGPase in rice. Homology modeling of the three-dimensional structure of the LS and SS was performed using the Swiss Model Server, and the models were evaluated and docked using GRAMM-X to obtain the stable heterodimer orientation (LS as receptor and SS as ligand) and then the heterotetrameric orientation. The initial heterotetrameric orientation was further refined using the RosettaDock Server. MD simulation of the representative heterodimer/tetramer was performed using NAMD, which indicated that the tail-to-tail interaction of LS and SS was more stable than the head-to-head orientation, and the heterotetramer energy was also minimized to ?767,011 kcal mol^?1. Subunit–subunit interaction studies were then carried out using the programs NACCESS and Dimplot. A total of 57 interface residues were listed in SS and 63 in LS. The residues plotted by Dimplot were similar to those listed by NACCESS. Multiple sequence alignment of the sequences of LS and SS from potato, maize and rice validated the interactions inferred in the study. RMSD of 1.093 Å was obtained on superimposition of the deduced heterotetramer on the template homo-tetramer (1YP2), showing the similarity between the two structures.     

Allosteric regulation of the higher plant ADP-glucose pyrophosphorylase is a product of synergy between the two subunits

The higher plant ADP-glucose pyrophosphorylase (AGPase) is a heterotetramer consisting of two regulatory large subunits (LSs) and two catalytic small subunits (SSs). To further characterize the roles of these subunits in determining enzyme function, different combinations of wildtype LS (LWT) and variant forms (LUpReg1, LM345) were co-expressed with wildtype SS (SWT) and variant forms (STG-15 and Sdevo330) and their enzyme properties compared to those measured for the heterotetrameric wildtype enzyme and SS homotetrameric enzymes. Analysis of the allosteric regulatory properties of the various enzymes indicates that although the LS is required for optimal activation by 3-phosphoglyceric acid and resistance to Pi, the overall allosteric regulatory and kinetic properties are specified by both subunits. Our results show that the regulatory and kinetic properties of AGPase are not simply due to the LS modulating the properties of the SS but, instead, are a product of synergistic interaction between the two subunits.     

Catalytic implications of the higher plant ADP-glucose pyrophosphorylase large subunit

ADP-glucose pyrophosphorylase, a key regulatory enzyme of starch biosynthesis, is composed of a pair of catalytic small subunits (SSs) and a pair of catalytically disabled large subunits (LSs). The N-terminal region of the LS has been known to be essential for the allosteric regulatory properties of the heterotetrameric enzyme. To gain further insight on the role of this region and the LS itself in enzyme function, the six proline residues found in the N-terminal region of the potato tuber AGPase were subjected to scanning mutagenesis. The wildtype and various mutant heterotetramers were expressed using our newly developed host-vector system, purified, and their kinetic parameters assessed. While P(17)L, P(26)L, and P(55)L mutations only moderately affected the kinetic properties, P(52)L and P(66)L gave rise to significant and contrasting changes in allosteric properties: P(66)L enzyme displayed up-regulatory properties toward 3-PGA while the P(52)L enzyme had down-regulatory properties. Unlike the other mutants, however, various mutations at P(44) led to only moderate changes in regulatory properties, but had severely impaired catalytic rates, apparent substrate affinities, and responsiveness to metabolic effectors, indicating Pro-44 or the LS is essential for optimal catalysis and activation of the AGPase heterotetramer. The catalytic importance of the LS is further supported by photoaffinity labeling studies, which revealed that the LS binds ATP at the same efficiency as the SS. These results indicate that the LS, although considered having no catalytic activity, may mimic many of the catalytic events undertaken by the SS and, thereby, influences net catalysis of the heterotetrameric enzyme.     

Direct appraisal of the potato tuber ADP-glucose pyrophosphorylase large subunit in enzyme function by study of a novel mutant form

The higher plant ADP-glucose pyrophosphorylase is a heterotetramer consisting of two subunit types, which have evolved at different rates from a common ancestral gene. The potato tuber small subunit (SS) displays both catalytic and regulatory properties, whereas the exact role of the large subunit (LS), which contains substrate and effector binding sites, remains unresolved. We identified a mutation, S302N, which increased the solubility of the recombinant potato tuber LS and, in turn, enabling it to form a homotetrameric structure. The LS302N homotetramer possesses very little enzyme activity at a level 100-fold less than that seen for the unactivated SS homotetramer. Unlike the SS enzyme, however, the LS302N homotetramer enzyme is neither activated by the effector 3-phosphoglycerate nor inhibited by P(i). When combined with the catalytically silenced SS, S D143N, however, the LS302N-containing enzyme shows significantly enhanced catalytic activity and restored 3-PGA activation. This unmasking of catalytic and regulatory potential of the LS is conspicuously evident when the activities of the resurrected L(K41R.T51K.S302N) homotetramer are compared with its heterotetrameric form assembled with S D143N. Overall, these results indicate that the LS possesses catalytic and regulatory properties only when assembled with SS and that the net properties of the heterotetrameric enzyme is a product of subunit synergy.     

Rapid purification of the potato ADP-glucose pyrophosphorylase by polyhistidine-mediated chromatography

In an attempt to obtain facile methods to purify the heterotetrameric ADP-glucose pyrophosphorylase (AGPase), polyhistidine tags were attached to either the large (LS) or small (SS) subunits of this oligomeric enzyme. The addition of polyhistidine tag to the N-terminus of the LS or SS and co-expression with its unmodified counterpart subunit resulted in substantial induction of enzyme activity. In contrast, attachment of a polyhistidine-containing peptide through the use of a commercially available pET vector or addition of polyhistidine tags to the C-terminal ends of either subunit resulted in poor expression and/or production of enzyme activity. Preliminary experiment showed that these polyhistidine N-terminal-tagged enzymes interacted with Ni-NTA-agarose, indicating that immobilized metal affinity chromatography (IMAC) would be useful for efficient purification of the heterotetrameric AGPases. When ion-exchange chromatography step was employed prior to the IMAC, the polyhistidine-tagged AGPases were purified to near homogeneity. Comparison of kinetic parameters between AGPases with and without the polyhistidine tags revealed that attachment of the polyhistidine did not alter the allosteric and catalytic properties of the enzymes. These results indicate that polyhistidine tags will be useful for the rapid purification of preparative amounts of AGPases for biochemical and physical studies.     

Determinants of human glucokinase activation and implications for small molecule allosteric control

Glucokinase (GK) is an enzyme that catalyzes the ATP-dependent phosphorylation of glucose to form glucose-6-phosphate, and it is a tightly regulated checkpoint in glucose homeostasis. GK is known to undergo substantial conformational changes upon glucose binding. The monomeric enzyme possesses a highly exotic kinetic activity profile with an unusual sigmoidal dependence on glucose concentration. In this interdisciplinary study, which draws on small angle X-ray scattering (SAXS) integrated with 250?ns of atomistic molecular dynamics (MD) simulations and experimental glucose binding thermodynamics, we reveal that the critical regulation of this glucose sensor is due to a solvent controlled “switch”. We demonstrate that the “solvent switch” is driven by specific protein structural dynamics, which leads to an enzyme structure that has a much more favorable solvation energy than most of the protein ensemble. These findings uncover the physical workings of an agile and flexible protein scaffold, which derives its long-range allosteric control through specific regions with favorable solvation energy. The physiological framework presented herein provides insights that have direct implications for the design of small molecule GK activators as anti-diabetes therapeutics as well as for understanding how proteins can be designed to have built-in regulatory functions via solvation energy dynamics.     

Molecular dynamics simulations and rigid body (TLS) analysis of aspartate carbamoyltransferase: evidence for an uncoupled R state.

In the R form of ATCase complexed with the bisubstrate analogue, N-(phosphonacetyl)-L-aspartate, large temperature factors are reported for the allosteric domains of the regulatory chains. We studied the conformational flexibility of the holoenzyme with molecular dynamics simulations and rigid body (TLS) analysis. The results of the molecular dynamics simulations suggest that, although local atomic fluctuations account for the temperature factors of the catalytic and zinc domains, they do not account for the large temperature factors of the allosteric regions. However, the temperature factors of the allosteric domains can be satisfactorily analyzed using a rigid body model. The simulations and rigid body analysis support the idea that the allosteric regions are mechanically uncoupled from the rest of the enzyme in the PALA structure. Implications of this uncoupling for allosteric regulation are discussed.     

Structural modulation of factor VIIa by full-length tissue factor (TF1-263): implication of novel interactions between EGF2 domain and TF

Tissue factor (TF)-mediated factor VII (FVII) activation and a subsequent proteolytic TF-FVIIa binary complex formation is the key step initiating the coagulation cascade, with implications in various homeostatic and pathologic scenarios. TF binding allosterically modifies zymogen-like free FVIIa to its highly catalytically active form. As a result of unresolved crystal structure of the full-length TF_1-263-FVIIa binary complex and free FVIIa, allosteric alterations in FVIIa following its binding to full-length TF and the consequences of these on function are not entirely clear. The present study aims to map and identify structural alterations in FVIIa and TF resulting from full-length TF binding to FVIIa and the key events responsible for enhanced FVIIa activity in coagulation. We constructed the full-length TF_1-263-FVIIa membrane bound complex using computational modeling and subjected it to molecular dynamics (MD) simulations. MD simulations showed that TF alters the structure of each domain of FVIIa and these combined alterations contribute to enhanced TF-FVIIa activity. Detailed, domain-wise investigation revealed several new non-covalent interactions between TF and FVIIa that were not found in the truncated soluble TF-FVIIa crystal structure. The structural modulation of each FVIIa domain imparted by TF indicated that both inter and intra-domain communication is crucial for allosteric modulation of FVIIa. Our results suggest that these newly formed interactions can provide additional stability to the protease domain and regulate its activity profile by governing catalytic triad (CT) orientation and localization. The unexplored newly formed interactions between EGF2 and TF provides a possible explanation for TF-induced allosteric activation of FVIIa.     

Ligand‐induced formation of transient dimers of mammalian 12/15‐lipoxygenase: A key to allosteric behavior of this class of enzymes?

Mammalian lipoxygenases (LOXs) have been implicated in cellular defense response and are important for physiological homeostasis. Since their discovery, LOXs have been believed to function as monomeric enzymes that exhibit allosteric properties. In aqueous solutions, the rabbit 12/15-LOX is mainly present as hydrated monomer but changes in the local physiochemical environment suggested a monomer-dimer equilibrium. Because the allosteric character of the enzyme can hardly be explained using a single ligand binding-site model, we proposed that the binding of allosteric effectors may shift the monomer-dimer equilibrium toward dimer formation. To test this hypothesis, we explored the impact of an allosteric effector [13(S)-hydroxyoctadeca-9(Z),11(E)-dienoic acid] on the structural properties of rabbit 12/15-LOX by small-angle X-ray scattering. Our data indicate that the enzyme undergoes ligand-induced dimerization in aqueous solution, and molecular dynamics simulations suggested that LOX dimers may be stable in the presence of substrate fatty acids. These data provide direct structural evidence for the existence of LOX dimers, where two noncovalently linked enzyme molecules might work in unison and, therefore, such mode of association might be related to the allosteric character of 12/15-LOX. Introduction of negatively charged residues (W181E + H585E and L183E + L192E) at the intermonomer interface disturbs the hydrophobic dimer interaction of the wild-type LOX, and this structural alteration may lead to functional distortion of mutant enzymes.     

Crystal, solution and in silico structural studies of dihydrodipicolinate synthase from the common grapevine

Dihydrodipicolinate synthase (DHDPS) catalyzes the rate limiting step in lysine biosynthesis in bacteria and plants. The structure of DHDPS has been determined from several bacterial species and shown in most cases to form a homotetramer or dimer of dimers. However, only one plant DHDPS structure has been determined to date from the wild tobacco species, Nicotiana sylvestris (Blickling et al. (1997) J. Mol. Biol. 274, 608-621). Whilst N. sylvestris DHDPS also forms a homotetramer, the plant enzyme adopts a 'back-to-back' dimer of dimers compared to the 'head-to-head' architecture observed for bacterial DHDPS tetramers. This raises the question of whether the alternative quaternary architecture observed for N. sylvestris DHDPS is common to all plant DHDPS enzymes. Here, we describe the structure of DHDPS from the grapevine plant, Vitis vinifera, and show using analytical ultracentrifugation, small-angle X-ray scattering and X-ray crystallography that V. vinifera DHDPS forms a 'back-to-back' homotetramer, consistent with N. sylvestris DHDPS. This study is the first to demonstrate using both crystal and solution state measurements that DHDPS from the grapevine plant adopts an alternative tetrameric architecture to the bacterial form, which is important for optimizing protein dynamics as suggested by molecular dynamics simulations reported in this study.     

Allosteric Inhibitors Have Distinct Effects,but Also Common Modes of Action,in the HCV Polymerase

The RNA-dependent RNA polymerase from the Hepatitis C Virus (gene product NS5B) is a validated drug target because of its critical role in genome replication. There are at least four distinct allosteric sites on the polymerase to which several small molecule inhibitors bind. In addition, numerous crystal structures have been solved with different allosteric inhibitors bound to the polymerase. However, the molecular mechanisms by which these small molecules inhibit the enzyme have not been fully elucidated. There is evidence that allosteric inhibitors alter the intrinsic motions and distribution of conformations sampled by the enzyme. In this study we use molecular dynamics simulations to understand the structural and dynamic changes that result when inhibitors are bound at three different allosteric binding sites on the enzyme. We observe that ligand binding at each site alters the structure and dynamics of NS5B in a distinct manner. Nonetheless, our studies also highlight commonalities in the mechanisms of action of the different inhibitors. Each inhibitor alters the conformational states sampled by the enzyme, either by rigidifying the enzyme and preventing transitions between functional conformational states or by destabilizing the enzyme and preventing functionally relevant conformations from being adequately sampled. By illuminating the molecular mechanisms of allosteric inhibition, these studies delineate the intrinsic functional properties of the enzyme and pave the way for designing novel and more effective polymerase inhibitors. This information may also be important to understand how allosteric regulation occurs in related viral polymerases and other enzymes.     

Energetics for displacing a single chain from the surface of microcrystalline cellulose into the active site of Acidothermus cellulolyticus Cel5A

A series of molecular mechanics calculations were used to analyze the energetics for moving a single polysaccharide chain from the surface of microcrystalline cellulose into the binding cleft of the Cel5A cellulase from Acidothermus cellulolyticus. A build-up procedure was used to model the placement of a 12-residue oligosaccharide chain along the surface of the enzyme, using as a guide the four residues of the tetrasaccharide substrate co-crystallized with the protein in the crystallographic structure determination. The position of this 12-residue oligosaccharide was used to orient the enzyme properly above two different surfaces of cellulose 1beta, the (1,0,0) and the (1,1,0) faces of the crystal. Constrained molecular dynamics simulations were then used to pull a target chain directly below the enzyme up out of the crystal surface and into the binding groove. The energetics for this process were favorable for both faces, with the step face being more favorable than the planar face, implying that this surface could be hydrolyzed more readily.     

Dynamic energy correlation analysis of E. coli aspartokinase III and alteration of allosteric regulation by manipulating energy transduction pathways

Conformational change associated with allosteric regulation in a protein is ultimately driven by energy transformation. However, little is known about the latter process. In this work, we combined steered molecular dynamics simulations and sequence conservation analysis to investigate the conformational changes and energy transformation in the allosteric enzyme aspartokinase III (AK III) from Escherichia coli. Correlation analysis of energy change at residue level indicated significant transformation between electrostatic energy and dihedral angle energy during the allosteric regulation. Key amino acid residues located in the corresponding energy transduction pathways were identified by dynamic energy correlation analysis. To verify their functions, residues with a high energy correlation in the pathways were altered and their effects on allosteric regulation of AKIII were determined. This study sheds new insights into energy transformation during allosteric regulation of AK III and proposes a strategy to identify key residues that are involved in intramolecular energy transduction and thus in driving the allosteric process.     

Dissecting mutational allosteric effects in alkaline phosphatases associated with different Hypophosphatasia phenotypes: An integrative computational investigation

Hypophosphatasia (HPP) is a rare inherited disorder characterized by defective bone mineralization and is highly variable in its clinical phenotype. The disease occurs due to various loss-of-function mutations in ALPL, the gene encoding tissue-nonspecific alkaline phosphatase (TNSALP). In this work, a data-driven and biophysics-based approach is proposed for the large-scale analysis of ALPL mutations-from nonpathogenic to severe HPPs. By using a pipeline of synergistic approaches including sequence-structure analysis, network modeling, elastic network models and atomistic simulations, we characterized allosteric signatures and effects of the ALPL mutations on protein dynamics and function. Statistical analysis of molecular features computed for the ALPL mutations showed a significant difference between the control, mild and severe HPP phenotypes. Molecular dynamics simulations coupled with protein structure network analysis were employed to analyze the effect of single-residue variation on conformational dynamics of TNSALP dimers, and the developed machine learning model suggested that the topological network parameters could serve as a robust indicator of severe mutations. The results indicated that the severity of disease-associated mutations is often linked with mutation-induced modulation of allosteric communications in the protein. This study suggested that ALPL mutations associated with mild and more severe HPPs can exert markedly distinct effects on the protein stability and long-range network communications. By linking the disease phenotypes with dynamic and allosteric molecular signatures, the proposed integrative computational approach enabled to characterize and quantify the allosteric effects of ALPL mutations and role of allostery in the pathogenesis of HPPs.     

Structural, kinetic and computational investigation of Vitis vinifera DHDPS reveals new insight into the mechanism of lysine-mediated allosteric inhibition

Lysine is one of the most limiting amino acids in plants and its biosynthesis is carefully regulated through inhibition of the first committed step in the pathway catalyzed by dihydrodipicolinate synthase (DHDPS). This is mediated via a feedback mechanism involving the binding of lysine to the allosteric cleft of DHDPS. However, the precise allosteric mechanism is yet to be defined. We present a thorough enzyme kinetic and thermodynamic analysis of lysine inhibition of DHDPS from the common grapevine, Vitis vinifera (Vv). Our studies demonstrate that lysine binding is both tight (relative to bacterial DHDPS orthologs) and cooperative. The crystal structure of the enzyme bound to lysine (2.4 Å) identifies the allosteric binding site and clearly shows a conformational change of several residues within the allosteric and active sites. Molecular dynamics simulations comparing the lysine-bound (PDB ID 4HNN) and lysine free (PDB ID 3TUU) structures show that Tyr132, a key catalytic site residue, undergoes significant rotational motion upon lysine binding. This suggests proton relay through the catalytic triad is attenuated in the presence of lysine. Our study reveals for the first time the structural mechanism for allosteric inhibition of DHDPS from the common grapevine.     

Dynamical probing of allosteric control in nuclear receptors

The dynamical behavior of the nuclear receptor LXR/RXR heterodimer was investigated with molecular dynamics simulations. The simulations reveal correlated motion between residues across the dimer interface that depends significantly on occupation of the ligand binding sites of the monomers. These results are broadly consistent with the observed experimental behavior of the dimers, where structural perturbation is thought to be a key element in signal transduction. Our results provide dynamical support for this model of allosteric control.     

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α₂ and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α₄-α₅ KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).     

Clarifying allosteric control of flap conformations in the 1TW7 crystal structure of HIV-1 protease

The 1TW7 crystal structure of HIV-1 protease shows the flaps placed wider and more open than what is seen in other examples of the semi-open, apo form. It has been proposed that this might be experimental evidence of allosteric control, because crystal packing creates contacts to the "elbow region" of the protease, which may cause deformation of the flaps. Recent dynamics simulations have shown that the conformation seen in 1TW7 relaxes into the typical semi-open conformation in the absence of the crystal contacts, definitively showing that the crystal contacts cause the deformation (Layten et al., J Am Chem Soc 2006;128:13360-13361). However, this does not prove or disprove allosteric modulation at the elbow. In this study, we have conducted additional simulations, supplemented with experimental testing, to further probe the possibility of 1TW7 providing an example of allosteric control of the flap region. We show that the contacts are unstable and do not restrict the conformational sampling of the flaps. The deformation seen in the 1TW7 crystal structure is simply opportunistic crystal packing and not allosteric control.