Compounds identified by virtual docking to a tetrameric EGFR extracellular domain can modulate Grb2 internalization.

BACKGROUND: Overexpression or mutation of the epidermal growth factor receptor (EGFR) potently enhances the growth of many solid tumors. Tumor cells frequently display resistance to mechanistically-distinct EGFR-directed therapeutic agents, making it valuable to develop therapeutics that work by additional mechanisms. Current EGFR-targeting therapeutics include antibodies targeting the extracellular domains, and small molecules inhibiting the intracellular kinase domain. Recent studies have identified a novel prone extracellular tetrameric EGFR configuration, which we identify as a potential target for drug discovery. METHODS: Our focus is on the prone EGFR tetramer, which contains a novel protein-protein interface involving extracellular domain III. This EGFR tetramer is computationally targeted for stabilization by small molecule ligand binding. This study performed virtual screening of a Life Chemicals, Inc. small molecule library of 345,232 drug-like compounds against a molecular dynamics simulation of protein-protein interfaces distinct to the novel tetramer. One hundred nine chemically diverse candidate molecules were selected and evaluated using a cell-based high-content imaging screen that directly assessed induced internalization of the EGFR effector protein Grb2. Positive hits were further evaluated for influence on phosphorylation of EGFR and its effector ERK1/2. RESULTS: Fourteen hit compounds affected internalization of Grb2, an adaptor responsive to EGFR activation. Most hits had limited effect on cell viability, and minimally influenced EGFR and ERK1/2 phosphorylation. Docked hit compound poses generally include Arg270 or neighboring residues, which are also involved in binding the effective therapeutic cetuximab, guiding further chemical optimization. CONCLUSIONS: These data suggest that the EGFR tetrameric configuration offers a novel cancer drug target.

A new model for allosteric regulation of phenylalanine hydroxylase: implications for disease and therapeutics.

The structural basis for allosteric regulation of phenylalanine hydroxylase (PAH), whose dysfunction causes phenylketonuria (PKU), is poorly understood. A new morpheein model for PAH allostery is proposed to consist of a dissociative equilibrium between two architecturally different tetramers whose interconversion requires a ∼90° rotation between the PAH catalytic and regulatory domains, the latter of which contains an ACT domain. This unprecedented model is supported by in vitro data on purified full length rat and human PAH. The conformational change is both predicted to and shown to render the tetramers chromatographically separable using ion exchange methods. One novel aspect of the activated tetramer model is an allosteric phenylalanine binding site at the intersubunit interface of ACT domains. Amino acid ligand-stabilized ACT domain dimerization follows the multimerization and ligand binding behavior of ACT domains present in other proteins in the PDB. Spectroscopic, chromatographic, and electrophoretic methods demonstrate a PAH equilibrium consisting of two architecturally distinct tetramers as well as dimers. We postulate that PKU-associated mutations may shift the PAH quaternary structure equilibrium in favor of the low activity assemblies. Pharmacological chaperones that stabilize the ACT:ACT interface can potentially provide PKU patients with a novel small molecule therapeutic.

Allostery and the dynamic oligomerization of porphobilinogen synthase.

The structural basis for allosteric regulation of porphobilinogen synthase (PBGS) is modulation of a quaternary structure equilibrium between octamer and hexamer (via dimers), which is represented schematically as 8mer ⇔ 2mer ⇔ 2mer∗⇔ 6mer∗. The "∗" represents a reorientation between two domains of each subunit that occurs in the dissociated state because it is sterically forbidden in the larger multimers. Allosteric effectors of PBGS are both intrinsic and extrinsic and are phylogenetically variable. In some species this equilibrium is modulated intrinsically by magnesium which binds at a site specific to the 8mer. In other species this equilibrium is modulated intrinsically by pH with the guanidinium group of an arginine being spatially equivalent to the allosteric magnesium ion. In humans, disease associated variants all shift the equilibrium toward the 6mer∗ relative to wild type. The 6mer∗ has a surface cavity that is not present in the 8mer and is proposed as a small molecule allosteric binding site. In silico and in vitro approaches have revealed species-specific allosteric PBGS inhibitors that stabilize the 6mer∗. Some of these inhibitors are drugs in clinical use leading to the hypothesis that extrinsic allosteric inhibition of human PBGS could be a mechanism for drug side effects.

Allosteric inhibition of human porphobilinogen synthase.

Porphobilinogen synthase (PBGS) catalyzes the first common step in tetrapyrrole (e.g. heme, chlorophyll) biosynthesis. Human PBGS exists as an equilibrium of high activity octamers, low activity hexamers, and alternate dimer configurations that dictate the stoichiometry and architecture of further assembly. It is posited that small molecules can be found that inhibit human PBGS activity by stabilizing the hexamer. Such molecules, if present in the environment, could potentiate disease states associated with reduced PBGS activity, such as lead poisoning and ALAD porphyria, the latter of which is associated with human PBGS variants whose quaternary structure equilibrium is shifted toward the hexamer (Jaffe, E. K., and Stith, L. (2007) Am. J. Hum. Genet. 80, 329-337). Hexamer-stabilizing inhibitors of human PBGS were identified using in silico prescreening (docking) of approximately 111,000 structures to a hexamer-specific surface cavity of a human PBGS crystal structure. Seventy-seven compounds were evaluated in vitro; three provided 90-100% conversion of octamer to hexamer in a native PAGE mobility shift assay. Based on chemical purity, two (ML-3A9 and ML-3H2) were subjected to further evaluation of their effect on the quaternary structure equilibrium and enzymatic activity. Naturally occurring ALAD porphyria-associated human PBGS variants are shown to have an increased susceptibility to inhibition by both ML-3A9 and ML-3H2. ML-3H2 is a structural analog of amebicidal drugs, which have porphyria-like side effects. Data support the hypothesis that human PBGS hexamer stabilization may explain these side effects. The current work identifies allosteric ligands of human PBGS and, thus, identifies human PBGS as a medically relevant allosteric enzyme.

Shape shifting leads to small-molecule allosteric drug discovery.

Enzymes that regulate their activity by modulating an equilibrium of alternate, nonadditive, functionally distinct oligomeric assemblies (morpheeins) constitute a recently described mode of allostery. The oligomeric equilibrium for porphobilinogen synthase (PBGS) consists of high-activity octamers, low-activity hexamers, and two dimer conformations. A phylogenetically diverse allosteric site specific to hexamers is proposed as an inhibitor binding site. Inhibitor binding is predicted to draw the oligomeric equilibrium toward the low-activity hexamer. In silico docking enriched a selection from a small-molecule library for compounds predicted to bind to this allosteric site. In vitro testing of selected compounds identified one compound whose inhibition mechanism is species-specific conversion of PBGS octamers to hexamers. We propose that this strategy for inhibitor discovery can be applied to other proteins that use the morpheein model for allosteric regulation.

Probing the oligomeric assemblies of pea porphobilinogen synthase by analytical ultracentrifugation.

The enzyme porphobilinogen synthase (PBGS) can exist in different nonadditive homooligomeric assemblies, and under appropriate conditions, the distribution of these assemblies can respond to ligands such as metals or substrate. PBGS from most organisms was believed to be octameric until work on a rare allele of human PBGS revealed an alternate hexameric assembly, which is also available to the wild-type enzyme at elevated pH [Breinig, S., et al. (2003) Nat. Struct. Biol. 10, 757-763]. Herein, we establish that the distribution of pea PBGS quaternary structures also contains octamers and hexamers, using both sedimentation velocity and sedimentation equilibrium experiments. We report results in which the octamer dominates under purification conditions and discuss conditions that influence the octamer:hexamer ratio. As predicted by PBGS crystal structures from related organisms, in the absence of magnesium, the octameric assembly is significantly destabilized, and the oligomeric distribution is dominated largely by the hexameric assembly. Although the PBGS hexamer-to-octamer oligomeric rearrangement is well documented under some conditions, both assemblies are very stable (under AU conditions) in the time frame of our ultracentrifuge experiments.

Crystal structure of phosphotransacetylase from the methanogenic archaeon Methanosarcina thermophila.

Phosphotransacetylase (Pta) [EC 2.3.1.8] is ubiquitous in the carbon assimilation and energy-yielding pathways in anaerobic prokaryotes where it catalyzes the reversible transfer of the acetyl group from acetyl phosphate to CoA forming acetyl CoA and inorganic phosphate. The crystal structure of Pta from the methane-producing archaeon Methanosarcina thermophila, representing the first crystal structure of any Pta, was determined by multiwavelength anomalous diffraction at 2.7 A resolution. In solution and in the crystal, the enzyme forms a homodimer. Each monomer consists of two alpha/beta domains with a cleft along the domain boundary, which presumably contains the substrate binding sites. Comparison of the four monomers present in the asymmetric unit indicates substantial variations in the relative orientation of the two domains and the structure of the putative active site cleft. A search for structural homologs revealed the NADP(+)-dependent isocitrate and isopropylmalate dehydrogenases as the only homologs with a similar two-domain architecture.

wALADin benzimidazoles differentially modulate the function of porphobilinogen synthase orthologs.

The heme biosynthesis enzyme porphobilinogen synthase (PBGS) is a potential drug target in several human pathogens. wALADin1 benzimidazoles have emerged as species-selective PBGS inhibitors against Wolbachia endobacteria of filarial worms. In the present study, we have systematically tested wALADins against PBGS orthologs from bacteria, protozoa, metazoa, and plants to elucidate the inhibitory spectrum. However, the effect of wALADin1 on different PBGS orthologs was not limited to inhibition: several orthologs were stimulated by wALADin1; others remained unaffected. We demonstrate that wALADins allosterically modulate the PBGS homooligomeric equilibrium with inhibition mediated by favoring low-activity oligomers, while 5-aminolevulinic acid, Mg(2+), or K(+) stabilized high-activity oligomers. Pseudomonas aeruginosa PBGS could be inhibited or stimulated by wALADin1 depending on these factors and pH. We have defined the wALADin chemotypes responsible for either inhibition or stimulation, facilitating the design of tailored PBGS modulators for potential application as antimicrobial agents, herbicides, or drugs for porphyric disorders.

Environmental contaminants perturb fragile protein assemblies and inhibit normal protein function.

The molecular mechanisms whereby small molecules that contaminate our environment cause physiological effects are largely unknown, in terms of both targets and mechanisms. The essential human enzyme porphobilinogen synthase (HsPBGS, a.k.a. 5-aminolevulinate dehydratase, ALAD) functions in heme biosynthesis. HsPBGS catalytic activity is regulated allosterically via an equilibrium of inactive hexamers and active octamers, and we have shown that certain drugs and drug-like small molecules can inhibit HsPBGS in vitro by stabilizing the hexamer. Here we address whether components of the National Toxicology Program library of environmental contaminants can stabilize the HsPBGS hexamer and inhibit activity in vitro. Native polyacrylamide gel electrophoresis was used to screen the library (1,408 compounds) for components that alter the oligomeric distribution of HsPBGS. Freshly purchased samples of 37 preliminary hits were used to confirm the electrophoretic results and to determine the dose-dependence of the perturbation of oligomeric distribution. Seventeen compounds were identified which alter the oligomeric distribution toward the hexamer and also inhibit HsPBGS catalytic activity, including the most potent HsPBGS inhibitor yet characterized (Mutagen X, IC50 = 1.4 µM). PBGS dysfunction is associated with the inborn error of metabolism know as ALAD porphyria and with lead poisoning. The identified hexamer-stabilizing inhibitors could potentiate these diseases. Allosteric regulation of activity via an equilibrium of alternate oligomers has been proposed for many proteins. Based on the precedent set herein, perturbation of these oligomeric equilibria by small molecules (such as environmental contaminants) can be considered as a mechanism of toxicity.

The morpheein model of allostery: evaluating proteins as potential morpheeins.

An equilibrium mixture of alternate quaternary structure assemblies can form a basis for allostery. The morpheein model of allostery is a concerted dissociative model that describes an equilibrium of alternate quaternary structure assemblies whose architectures are dictated by alternate conformations in the dissociated state. Kinetic and biophysical anomalies that suggest that the morpheein model of allostery applies for a given protein of interest are briefly described. Two methods are presented for evaluating proteins as potential morpheeins. One is a subunit interchange method that uses chromatography, dialysis, and mass spectroscopy to monitor changes in multimer composition. The other is a two-dimensional native gel electrophoresis method to monitor ligand-induced changes in an equilibrium of alternate multimeric assemblies.

Diverse clinical compounds alter the quaternary structure and inhibit the activity of an essential enzyme.

An in vitro evaluation of the Johns Hopkins Clinical Compound Library demonstrates that certain drugs can alter the quaternary structure of an essential human protein. Human porphobilinogen synthase (HsPBGS) is an essential enzyme involved in heme biosynthesis; it exists as an equilibrium of high-activity octamers, low-activity hexamers, and alternate dimer configurations that dictate the stoichiometry and architecture of further assembly. Decreased HsPBGS activity is implicated in toxicities associated with lead poisoning and 5-aminolevulinate dehydratase (ALAD) porphyria, the latter of which involves hexamer-favoring HsPBGS variants. A medium-throughput native PAGE mobility-shift screen coupled with evaluation of hits as HsPBGS inhibitors revealed 12 drugs that stabilize the HsPBGS hexamer and inhibit HsPBGS activity in vitro. A detailed characterization of these effects is presented. Drug inhibition of HsPBGS in vivo by inducing hexamer formation would constitute an unprecedented mechanism for side effects. We suggest that small-molecule perturbation of quaternary structure equilibria be considered as a general mechanism for drug action and side effects.

Expanding the Concepts in Protein Structure-Function Relationships and Enzyme Kinetics: Teaching using Morpheeins.

A morpheein is a homo-oligomeric protein that can exist as an ensemble of physiologically significant and functionally different alternate quaternary assemblies. Morpheeins exist in nature and utilize conformational equilibria between different tertiary structures to form distinct oligomers as a means of regulating their function. Notably, alternate morpheeins are not misfolded forms of a protein; they are differently assembled native states that contain alternate subunit conformations. Transitions between alternate morpheein assemblies involve oligomer dissociation, conformational change in the dissociated state, and reassembly to a different oligomer. These transitions occur in response to the protein's environment, e.g., effector molecules, and represent a new model of allosteric regulation. The unique features of morpheeins are being revealed through detailed characterization of the prototype enzyme, porphobilinogen synthase, which exists in a dynamic equilibrium of a high activity octamer, a low activity hexamer, and two dimer conformations. Morpheeins are likely far more common than previously appreciated. There are, however, both intellectual and experimental barriers to recognizing proteins as morpheeins. These barriers derive from the way we were taught and continue to teach about protein folding, protein purification, protein structure-function relationships, and enzyme kinetics. This article explores some of these limitations and encourages incorporation of morpheeins into both introductory and advanced biochemistry classes.

Kinetics and thermodynamics of the interchange of the morpheein forms of human porphobilinogen synthase.

A morpheein is a homo-oligomeric protein that can adopt different nonadditive quaternary assemblies (morpheein forms) with different functionalities. The human porphobilinogen synthase (PBGS) morpheein forms are a high activity octamer, a low activity hexamer, and two structurally distinct dimer conformations. Conversion between hexamer and octamer involves dissociation to dimers, conformational change at the dimer level, followed by association to the alternate assembly. The current work promotes an alternative and novel view of the physiologically relevant dimeric structures, which are derived from the crystal structures, but are distinct from the asymmetric units of their crystal forms. Using a well characterized heteromeric system (WT+F12L; Tang, L. et al. (2005) J. Biol. Chem. 280, 15786-15793), extensive study of the human PBGS morpheein reequilibration process now reveals that the intervening dimers do not dissociate to monomers. The morpheein equilibria of wild type (WT) human PBGS are found to respond to changes in pH, PBGS concentration, and substrate turnover. Notably, the WT enzyme is predominantly an octamer at neutral pH, but increasing pH results in substantial conversion to lower order oligomers. Most significantly, the free energy of activation for the conversion of WT+F12L human PBGS heterohexamers to hetero-octamers is determined to be the same as that for the catalytic conversion of substrate to product by the octamer, remarkably suggesting a common rate-limiting step for both processes, which is postulated to be the opening/closing of the active site lid.

Steady-state kinetic analysis of phosphotransacetylase from Methanosarcina thermophila.

Phosphotransacetylase (EC 2.3.1.8) catalyzes the reversible transfer of the acetyl group from acetyl phosphate to coenzyme A (CoA), forming acetyl-CoA and inorganic phosphate. A steady-state kinetic analysis of the phosphotransacetylase from Methanosarcina thermophila indicated that there is a ternary complex kinetic mechanism rather than a ping-pong kinetic mechanism. Additionally, inhibition patterns of products and a nonreactive substrate analog suggested that the substrates bind to the enzyme in a random order. Dynamic light scattering revealed that the enzyme is dimeric in solution.

Structural and functional studies suggest a catalytic mechanism for the phosphotransacetylase from Methanosarcina thermophila.

Phosphotransacetylase (EC 2.3.1.8) catalyzes reversible transfer of the acetyl group from acetyl phosphate to coenzyme A (CoA), forming acetyl-CoA and inorganic phosphate. Two crystal structures of phosphotransacetylase from the methanogenic archaeon Methanosarcina thermophila in complex with the substrate CoA revealed one CoA (CoA1) bound in the proposed active site cleft and an additional CoA (CoA2) bound at the periphery of the cleft. The results of isothermal titration calorimetry experiments are described, and they support the hypothesis that there are distinct high-affinity (equilibrium dissociation constant [KD], 20 microM) and low-affinity (KD, 2 mM) CoA binding sites. The crystal structures indicated that binding of CoA1 is mediated by a series of hydrogen bonds and extensive van der Waals interactions with the enzyme and that there are fewer of these interactions between CoA2 and the enzyme. Different conformations of the protein observed in the crystal structures suggest that domain movements which alter the geometry of the active site cleft may contribute to catalysis. Kinetic and calorimetric analyses of site-specific replacement variants indicated that there are catalytic roles for Ser309 and Arg310, which are proximal to the reactive sulfhydryl of CoA1. The reaction is hypothesized to proceed through base-catalyzed abstraction of the thiol proton of CoA by the adjacent and invariant residue Asp316, followed by nucleophilic attack of the thiolate anion of CoA on the carbonyl carbon of acetyl phosphate. We propose that Arg310 binds acetyl phosphate and orients it for optimal nucleophilic attack. The hypothesized mechanism proceeds through a negatively charged transition state stabilized by hydrogen bond donation from Ser309.

Structural and kinetic analyses of arginine residues in the active site of the acetate kinase from Methanosarcina thermophila.

Acetate kinase catalyzes transfer of the gamma-phosphate of ATP to acetate. The only crystal structure reported for acetate kinase is the homodimeric enzyme from Methanosarcina thermophila containing ADP and sulfate in the active site (Buss, K. A., Cooper, D. C., Ingram-Smith, C., Ferry, J. G., Sanders, D. A., and Hasson, M. S. (2001) J. Bacteriol. 193, 680-686). Here we report two new crystal structure of the M. thermophila enzyme in the presence of substrate and transition state analogs. The enzyme co-crystallized with the ATP analog adenosine 5'-[gamma-thio]triphosphate contained AMP adjacent to thiopyrophosphate in the active site cleft of monomer B. The enzyme co-crystallized with ADP, acetate, Al(3+), and F(-) contained a linear array of ADP-AlF(3)-acetate in the active site cleft of monomer B. Together, the structures clarify the substrate binding sites and support a direct in-line transfer mechanism in which AlF(3) mimics the meta-phosphate transition state. Monomers A of both structures contained ADP and sulfate, and the active site clefts were closed less than in monomers B, suggesting that domain movement contributes to catalysis. The finding that His(180) was in close proximity to AlF(3) is consistent with a role for stabilization of the meta-phosphate that is in agreement with a previous report indicating that this residue is essential for catalysis. Residue Arg(241) was also found adjacent to AlF(3), consistent with a role for stabilization of the transition state. Kinetic analyses of Arg(241) and Arg(91) replacement variants indicated that these residues are essential for catalysis and also indicated a role in binding acetate.

Characterization of the acetate binding pocket in the Methanosarcina thermophila acetate kinase.

Acetate kinase catalyzes the reversible magnesium-dependent synthesis of acetyl phosphate by transfer of the ATP gamma-phosphoryl group to acetate. Inspection of the crystal structure of the Methanosarcina thermophila enzyme containing only ADP revealed a solvent-accessible hydrophobic pocket formed by residues Val(93), Leu(122), Phe(179), and Pro(232) in the active site cleft, which identified a potential acetate binding site. The hypothesis that this was a binding site was further supported by alignment of all acetate kinase sequences available from databases, which showed strict conservation of all four residues, and the recent crystal structure of the M. thermophila enzyme with acetate bound in this pocket. Replacement of each residue in the pocket produced variants with K(m) values for acetate that were 7- to 26-fold greater than that of the wild type, and perturbations of this binding pocket also altered the specificity for longer-chain carboxylic acids and acetyl phosphate. The kinetic analyses of variants combined with structural modeling indicated that the pocket has roles in binding the methyl group of acetate, influencing substrate specificity, and orienting the carboxyl group. The kinetic analyses also indicated that binding of acetyl phosphate is more dependent on interactions of the phosphate group with an unidentified residue than on interactions between the methyl group and the hydrophobic pocket. The analyses also indicated that Phe(179) is essential for catalysis, possibly for domain closure. Alignments of acetate kinase, propionate kinase, and butyrate kinase sequences obtained from databases suggested that these enzymes have similar catalytic mechanisms and carboxylic acid substrate binding sites.

Expression, purification, crystallization and preliminary X-ray analysis of phosphotransacetylase from Methanosarcina thermophila.

Phosphotransacetylase (Pta) from the anaerobic archaeon Methanosarcina thermophila has been heterologously expressed in a soluble form which facilitated crystallization using the hanging-drop vapor-diffusion method with ammonium sulfate as a precipitant. This is the first report of the crystallization of any Pta. While the M. thermophila Pta has high sequence identity to Ptas from other organisms, it has no homology to any previously crystallized proteins. The protein crystallized in space group I4(1), with unit-cell parameters a = b = 114.8, c = 127.8 A, alpha = beta = gamma = 90 degrees. The crystals diffracted to 2.5 A resolution using Cu Kalpha radiation. The enzyme had previously been reported to exist as a monomer; however, the self-rotation function showed the presence of a non-crystallographic symmetry axis at psi = 90, phi = 90, kappa = 180 degrees, suggesting oligomerization. Dynamic light-scattering analysis supported a dimeric state for Pta in solution.