A direct interaction between NQO1 and a chemotherapeutic dimeric naphthoquinone

BACKGROUND

Multimeric naphthoquinones are redox-active compounds that exhibit antineoplastic, antiprotozoal, and antiviral activities. Due to their multimodal effect on perturbation of cellular oxidative state, these compounds hold great potential as therapeutic agents against highly proliferative neoplastic cells. In our previous work, we developed a series of novel dimeric naphthoquinones and showed that they were selectively cytotoxic to human acute myeloid leukemia (AML), breast and prostate cancer cell lines. We subsequently identified the oxidoreductase NAD(P)H dehydrogenase, quinone 1 (NQO1) as the major target of dimeric naphthoquinones and proposed a mechanism of action that entailed induction of a futile redox cycling.

RESULTS

Here, for the first time, we describe a direct physical interaction between the bromohydroxy dimeric naphthoquinone E6a and NQO1. Moreover, our studies reveal an extensive binding interface between E6a and the isoalloxazine ring of the flavin adenine dinucleotide (FAD) cofactor of NQO1 in addition to interactions with protein side chains in the active site. We also present biochemical evidence that dimeric naphthoquinones affect the redox state of the FAD cofactor of NQO1. Comparison of the mode of binding of E6a with those of other chemotherapeutics reveals unique characteristics of the interaction that can be leveraged in future drug optimization efforts.

CONCLUSION

The first structure of a dimeric naphthoquinone-NQO1 complex was reported, which can be used for design and synthesis of more potent next generation dimeric naphthoquinones to target NQO1 with higher affinity and specificity.

Thymine DNA glycosylase exhibits negligible affinity for nucleobases that it removes from DNA

Thymine DNA Glycosylase (TDG) performs essential functions in maintaining genetic integrity and epigenetic regulation. Initiating base excision repair, TDG removes thymine from mutagenic G·T mispairs caused by 5-methylcytosine (mC) deamination and other lesions including uracil (U) and 5-hydroxymethyluracil (hmU). In DNA demethylation, TDG excises 5-formylcytosine (fC) and 5-carboxylcytosine (caC), which are generated from mC by Tet (ten–eleven translocation) enzymes. Using improved crystallization conditions, we solved high-resolution (up to 1.45 Å) structures of TDG enzyme–product complexes generated from substrates including G·U, G·T, G·hmU, G·fC and G·caC. The structures reveal many new features, including key water-mediated enzyme–substrate interactions. Together with nuclear magnetic resonance experiments, the structures demonstrate that TDG releases the excised base from its tight product complex with abasic DNA, contrary to previous reports. Moreover, DNA-free TDG exhibits no significant binding to free nucleobases (U, T, hmU), indicating a K_d >> 10 mM. The structures reveal a solvent-filled channel to the active site, which might facilitate dissociation of the excised base and enable caC excision, which involves solvent-mediated acid catalysis. Dissociation of the excised base allows TDG to bind the beta rather than the alpha anomer of the abasic sugar, which might stabilize the enzyme–product complex.

A novel member of glycoside hydrolase family 30 subfamily 8 with altered substrate specificity

Endoxylanases classified into glycoside hydrolase family 30 subfamily 8 (GH30-8) are known to hydrolyze the hemicellulosic polysaccharide glucuronoxylan (GX) but not arabinoxylan or neutral xylooligosaccharides. This is owing to the specificity of these enzymes for the α-1,2-linked glucuronate (GA) appendage of GX. Limit hydrolysis of this substrate produces a series of aldouronates each containing a single GA substituted on the xylose penultimate to the reducing terminus. In this work, the structural and biochemical characterization of xylanase 30A from Clostridium papyrosolvens (CpXyn30A) is presented. This xylanase possesses a high degree of amino-acid identity to the canonical GH30-8 enzymes, but lacks the hallmark β8-α8 loop region which in part defines the function of this GH30 subfamily and its role in GA recognition. CpXyn30A is shown to have a similarly low activity on all xylan substrates, while hydrolysis of xylohexaose revealed a competing transglycosylation reaction. These findings are directly compared with the model GH30-8 enzyme from Bacillus subtilis, XynC. Despite its high sequence identity to the GH30-8 enzymes, CpXyn30A does not have any apparent specificity for the GA appendage. These findings confirm that the typically conserved β8-α8 loop region of these enzymes influences xylan substrate specificity but not necessarily β-1,4-xylanase function.

Evidence vs Expectations: How to validate your ligand in a protein structure

Determination of a protein-ligand complex structure is essential in many areas of structural biology. Details of the interactions between protein and a small molecule ligand often represent major findings from a crystal structure. Thorough validation of interpretation of such structural data is particularly important given high expectation of confirming prior experimental findings regarding targeted protein-ligand interaction. Modern methods of ligand validation are discussed and illustrated.

Structure of human apurinic/apyrimidinic endonuclease 1 with the essential Mg2+ cofactor

Apurinic/apyrimidinic endonuclease 1 (APE1) mediates the repair of abasic sites and other DNA lesions and is essential for base-excision repair and strand-break repair pathways. APE1 hydrolyzes the phosphodiester bond at abasic sites, producing 5'-deoxyribose phosphate and the 3'-OH primer needed for repair synthesis. It also has additional repair activities, including the removal of 3'-blocking groups. APE1 is a powerful enzyme that absolutely requires Mg2+, but the stoichiometry and catalytic function of the divalent cation remain unresolved for APE1 and for other enzymes in the DNase I superfamily. Previously reported structures of DNA-free APE1 contained either Sm3+ or Pb2+ in the active site. However, these are poor surrogates for Mg2+ because Sm3+ is not a cofactor and Pb2+ inhibits APE1, and their coordination geometry is expected to differ from that of Mg2+. A crystal structure of human APE1 was solved at 1.92 Å resolution with a single Mg2+ ion in the active site. The structure reveals ideal octahedral coordination of Mg2+ via two carboxylate groups and four water molecules. One residue that coordinates Mg2+ directly and two that bind inner-sphere water molecules are strictly conserved in the DNase I superfamily. This structure, together with a recent structure of the enzyme-product complex, inform on the stoichiometry and the role of Mg2+ in APE1-catalyzed reactions.

Structural and functional analysis of the pro-domain of human cathelicidin, LL-37

Cathelicidins form a family of small host defense peptides distinct from another class of cationic antimicrobial peptides, the defensins. They are expressed as large precursor molecules with a highly conserved pro-domain known as the cathelin-like domain (CLD). CLDs have high degrees of sequence homology to cathelin, a protein isolated from pig leukocytes and belonging to the cystatin family of cysteine protease inhibitors. In this report, we describe for the first time the X-ray crystal structure of the human CLD (hCLD) of the sole human cathelicidin, LL-37. The structure of the hCLD, determined at 1.93 Å resolution, shows the cystatin-like fold and is highly similar to the structure of the CLD of the pig cathelicidin, protegrin-3. We assayed the in vitro antibacterial activities of the hCLD, LL-37, and the precursor form, pro-cathelicidin (also known as hCAP18), and we found that the unprocessed protein inhibited the growth of Gram-negative bacteria with efficiencies comparable to that of the mature peptide, LL-37. In addition, the antibacterial activity of LL-37 was not inhibited by the hCLD intermolecularly, because exogenously added hCLD had no effect on the bactericidal activity of the mature peptide. The hCLD itself lacked antimicrobial function and did not inhibit the cysteine protease, cathepsin L. Our results contrast with previous reports of hCLD activity. A comparative structural analysis between the hCLD and the cysteine protease inhibitor stefin A showed why the hCLD is unable to function as an inhibitor of cysteine proteases. In this respect, the cystatin scaffold represents an ancestral structural platform from which proteins evolved divergently, with some losing inhibitory functions.

Visualizing ligand molecules in Twilight electron density

Three-dimensional models of protein structures determined by X-ray crystallography are based on the interpretation of experimentally derived electron-density maps. The real-space correlation coefficient (RSCC) provides an easily comprehensible, objective measure of the residue-based fit of atom coordinates to electron density. Among protein structure models, protein-ligand complexes are of special interest, given their contribution to understanding the molecular underpinnings of biological activity and to drug design. For consumers of such models, it is not trivial to determine the degree to which ligand-structure modelling is biased by subjective electron-density interpretation. A standalone script, Twilight, is presented for the analysis, visualization and annotation of a pre-filtered set of 2815 protein-ligand complexes deposited with the PDB as of 15 January 2012 with ligand RSCC values that are below a threshold of 0.6. It also provides simplified access to the visualization of any protein-ligand complex available from the PDB and annotated by the Uppsala Electron Density Server. The script runs on various platforms and is available for download at http://www.ruppweb.org/twilight/.

Techniques, tools and best practices for ligand electron-density analysis and results from their application to deposited crystal structures

As a result of substantial instrumental automation and the continuing improvement of software, crystallographic studies of biomolecules are conducted by non-experts in increasing numbers. While improved validation almost ensures that major mistakes in the protein part of structure models are exceedingly rare, in ligand-protein complex structures, which in general are most interesting to the scientist, ambiguous ligand electron density is often difficult to interpret and the modelled ligands are generally more difficult to properly validate. Here, (i) the primary technical reasons and potential human factors leading to problems in ligand structure models are presented; (ii) the most common categories of building errors or overinterpretation are classified; (iii) a few instructive and specific examples are discussed in detail, including an electron-density-based analysis of ligand structures that do not contain any ligands; (iv) means of avoiding such mistakes are suggested and the implications for database validity are discussed and (v) a user-friendly software tool that allows non-expert users to conveniently inspect ligand density is provided.

Structure of extracellular signal-regulated kinase 2 in complex with ATP and ADP

Extracellular signal-regulated kinases 1 and 2 (ERK1 and ERK2) are members of the mitogen-activated protein (MAP) kinase family. Constitutive activation of the ERK proteins contributes to the development and progression of numerous human tumors. Thus, ERK1 and ERK2 are promising targets for the design and the development of anticancer drugs. The detailed structural analysis of ERK complexed with ATP can provide valuable information for the design of new ligands that can bind in the ATP-binding pocket and inhibit ERK activity. In this study, the structures of apo-form ERK2 and of its complexes with the substrate ATP and the product ADP were determined. Comparison with the structural homolog cyclin-dependent kinase 2 reveals differences in the way that the ATP binding to the protein is mediated by magnesium. Only minor conformational changes are identified that occur upon substrate binding, and these are limited to the active-site residues.

On the variability of experimental data in macromolecular crystallography

Experimental errors as determined by data-processing algorithms in macromolecular crystallography are compared with the direct error estimates obtained by a multiple crystal data-collection protocol. It is found that several-fold error inflation is necessary to account for crystal-to-crystal variation. It is shown that similar error inflation is observed for data collected from multiple sections of the same crystal, indicating non-uniform crystal growth as one of the likely sources of additional data variation. Other potential sources of error inflation include differential X-ray absorption for different reflections and variation of unit-cell parameters. The underestimation of the experimental errors is more severe in lower resolution shells and for reflections characterized by a higher signal-to-noise ratio. These observations partially account for the gap between the expected and the observed R values in macromolecular crystallography.

Rapid Crystallization of L-Alanine on Engineered Surfaces using Metal-Assisted and Microwave-Accelerated Evaporative Crystallization

This study demonstrates the application of metal-assisted and microwave-accelerated evaporative crystallization (MA-MAEC) technique to rapid crystallization of L-alanine on surface engineered silver nanostructures. In this regard, silver island films (SIFs) were modified with hexamethylenediamine (HMA), 1-undecanethiol (UDET), and 11-mercaptoundecanoic acid (MUDA), which introduced -NH(2), -CH(3) and -COOH functional groups to SIFs, respectively. L-Alanine was crystallized on these engineered surfaces and blank SIFs at room temperature and using MA-MAEC technique. Significant improvements in crystal size, shape, and quality were observed on HMA-, MUDA- and UDET-modified SIFs at room temperature (crystallization time = 144, 40 and 147 min, respectively) as compared to those crystals grown on blank SIFs. Using the MA-MAEC technique, the crystallization time of L-alanine on engineered surfaces were reduced to 17 sec for microwave power level 10 (i.e., duty cycle 100%) and 7 min for microwave power level 1 (duty cycle 10%). Raman spectroscopy and powder x-ray diffraction (XRD) measurements showed that L-Alanine crystals grown on engineered surfaces using MA-MAEC technique had identical characteristic peaks of L-alanine crystals grown using traditional evaporative crystallization.

N-Acetylfarnesylcysteine is a novel class of peroxisome proliferator-activated receptor γ ligand with partial and full agonist activity in vitro and in vivo

The thiazolidedione (TZD) class of drugs is clinically approved for the treatment of type 2 diabetes. The therapeutic actions of TZDs are mediated via activation of peroxisome proliferator-activated receptor γ (PPARγ). Despite their widespread use, concern exists regarding the safety of currently used TZDs. This has prompted the development of selective PPARγ modulators (SPPARMs), compounds that promote glucose homeostasis but with reduced side effects due to partial PPARγ agonism. However, this also results in partial agonism with respect to PPARγ target genes promoting glucose homeostasis. Using a gene expression-based screening approach we identified N-acetylfarnesylcysteine (AFC) as both a full and partial agonist depending on the PPARγ target gene (differential SPPARM). AFC activated PPARγ as effectively as rosiglitazone with regard to Adrp, Angptl4, and AdipoQ, but was a partial agonist of aP2, a PPARγ target gene associated with increased adiposity. Induction of adipogenesis by AFC was also attenuated compared with rosiglitazone. Reporter, ligand binding assays, and dynamic modeling demonstrate that AFC binds and activates PPARγ in a unique manner compared with other PPARγ ligands. Importantly, treatment of mice with AFC improved glucose tolerance similar to rosiglitazone, but AFC did not promote weight gain to the same extent. Finally, AFC had effects on adipose tissue remodeling similar to those of rosiglitazone and had enhanced antiinflammatory effects. In conclusion, we describe a new approach for the identification of differential SPPARMs and have identified AFC as a novel class of PPARγ ligand with both full and partial agonist activity in vitro and in vivo.

Holo-Ni(II)HpNikR is an asymmetric tetramer containing two different nickel-binding sites

The metalloregulatory protein NikR from Helicobacter pylori (HpNikR) is a master regulator of gene expression which both activates and represses specific genes in response to nickel availability. Here, we report the first crystal structure (at 2.37 Å resolution) of Ni(II)HpNikR prepared directly from the holo protein. The protein contains four nickel ions located in two distinct coordination environments. Two nickel ions are bound to sites in a four-coordinate square-planar geometry as predicted on the basis of the structures of NikR from Escherichia coli and Pyrococcus horikoshii . The remaining two nickel ions are bound to sites with unexpected 5- or 6-coordination geometries which were previously thought to be involved in nickel incorporation into the protein. The nickel with 5-/6-coordination geometry utilizes three histidines from two separate monomeric HpNikR units along with two or three water molecules as ligands. The spatial location of the nickel in the 5-/6-coordinate site is within approximately 5 Å of the expected site if a 4-coordinate square-planar geometry occurred. Two of the histidines that participate as ligands in the 5-/6-coordinate site would also participate as ligands if the 4-coordinate site was occupied, making it impossible for both sites to be occupied simultaneously. DFT calculations show that the 5-/6-coordinate geometries are energetically favorable when the local protein environment is included in the calculations. The presence of two distinct coordination environments in HpNikR is suggested to be related to the specificity and binding affinity of this transcription factor for DNA.

Ligand bound structures of a glycosyl hydrolase family 30 glucuronoxylan xylanohydrolase

Xylanases of glycosyl hydrolase family 30 (GH30) have been shown to cleave β-1,4 linkages of 4-O-methylglucuronoxylan (MeGX(n)) as directed by the position along the xylan chain of an α-1,2-linked 4-O-methylglucuronate (MeGA) moiety. Complete hydrolysis of MeGX(n) by these enzymes results in singly substituted aldouronates having a 4-O-methylglucuronate moiety linked to a xylose penultimate from the reducing terminal xylose and some number of xylose residues toward the nonreducing terminus. This novel mode of action distinguishes GH30 xylanases from the more common xylanase families that cleave MeGX(n) in accessible regions. To help understand this unique biochemical function, we have determined the structure of XynC in its native and ligand-bound forms. XynC structure models derived from diffraction data of XynC crystal soaks with the simple sugar glucuronate (GA) and the tetrameric sugar 4-O-methyl-aldotetrauronate resulted in models containing GA and 4-O-methyl-aldotriuronate, respectively. Each is observed in two locations within XynC surface openings. Ligand coordination occurs within the XynC catalytic substrate binding cleft and on the structurally fused side β-domain, demonstrating a substrate targeting role for this putative carbohydrate binding module. Structural data reveal that GA acts as a primary functional appendage for recognition and hydrolysis of the MeGX(n) polymer by the protein. This work compares the structure of XynC with a previously reported homologous enzyme, XynA, from Erwinia chrysanthemi and analyzes the ligand binding sites. Our results identify the molecular interactions that define the unique function of XynC and homologous GH30 enzymes.

Mechanism of inactivation of Escherichia coli aspartate aminotransferase by (S)-4-amino-4,5-dihydro-2-furancarboxylic acid

As a potential drug to treat neurological diseases, the mechanism-based inhibitor (S)-4-amino-4,5-dihydro-2-furancarboxylic acid (S-ADFA) has been found to inhibit the γ-aminobutyric acid aminotransferase (GABA-AT) reaction. To circumvent the difficulties in structural studies of a S-ADFA-enzyme complex using GABA-AT, l-aspartate aminotransferase (l-AspAT) from Escherichia coli was used as a model PLP-dependent enzyme. Crystal structures of the E. coli aspartate aminotransferase with S-ADFA bound to the active site were obtained via cocrystallization at pH 7.5 and 8. The complex structures suggest that S-ADFA inhibits the transamination reaction by forming adducts with the catalytic lysine 246 via a covalent bond while producing 1 equiv of pyridoxamine 5'-phosphate (PMP). Based on the structures, formation of the K246-S-ADFA adducts requires a specific initial binding configuration of S-ADFA in the l-AspAT active site, as well as deprotonation of the ε-amino group of lysine 246 after the formation of the quinonoid and/or ketimine intermediate in the overall inactivation reaction.

A conserved protein interaction interface on the type 5 G protein beta subunit controls proteolytic stability and activity of R7 family regulator of G protein signaling proteins

Regulators of G protein signaling (RGS) proteins of the R7 subfamily limit signaling by neurotransmitters in the brain and by light in the retina. They form obligate complexes with the Gβ5 protein that are subject to proteolysis to control their abundance and alter signaling. The mechanisms that regulate this proteolysis, however, remain unclear. We used genetic screens to find mutations in Gβ5 that selectively destabilize one of the R7 RGS proteins in Caenorhabditis elegans. These mutations cluster at the binding interface between Gβ5 and the N terminus of R7 RGS proteins. Equivalent mutations within mammalian Gβ5 allowed the interface to still bind the N-terminal DEP domain of R7 RGS proteins, and mutant Gβ5-R7 RGS complexes initially formed in cells but were then rapidly degraded by proteolysis. Molecular dynamics simulations suggest the mutations weaken the Gβ5-DEP interface, thus promoting dynamic opening of the complex to expose determinants of proteolysis known to exist on the DEP domain. We propose that conformational rearrangements at the Gβ5-DEP interface are key to controlling the stability of R7 RGS protein complexes.

Consolidation of glycosyl hydrolase family 30: a dual domain 4/7 hydrolase family consisting of two structurally distinct groups

In this work glycosyl hydrolase (GH) family 30 (GH30) is analyzed and shown to consist of its currently classified member sequences as well as several homologous sequence groups currently assigned within family GH5. A large scale amino acid sequence alignment and a phylogenetic tree were generated and GH30 groups and subgroups were designated. A partial rearrangement in the GH30 defining side-associated β-domain contributes to the differentiation of two major groups that contain up to eight subgroups. For this CAZy family of Clan A enzymes the dual domain fold is conserved, suggesting that it may be a requirement for evolved function. This work redefines GH family 30 and serves as a guide for future efforts regarding enzymes classified within this family.

Percentile-based spread: a more accurate way to compare crystallographic models

The comparison of biomacromolecular crystal structures is traditionally based on the root-mean-square distance between corresponding atoms. This measure is sensitive to the presence of outliers, which inflate it disproportionately to their fraction. An alternative measure, the percentile-based spread (p.b.s.), is proposed and is shown to represent the average variation in atomic positions more adequately. It is discussed in the context of isomorphous crystal structures, conformational changes and model ensembles generated by repetitive automated rebuilding.