Systematic review of response criteria and endpoints in autoimmune hepatitis by the International Autoimmune Hepatitis Group

BACKGROUND & AIMS

Autoimmune hepatitis (AIH) has been well characterised and codified through the development of diagnostic criteria. These criteria have been adapted and simplified and are widely used in clinical practice. However, there is a need to update and precisely define the criteria for both treatment response and treatment.

METHODS

A systematic review was performed and a modified Delphi consensus process was used to identify and redefine the response criteria in autoimmune hepatitis.

RESULTS

The consensus process initiated by the International Autoimmune Hepatitis Group proposes that the term 'complete biochemical response' defined as 'normalization of serum transaminases and IgG below the upper limit of normal' be adopted to include a time point at 6 months after initiation of treatment. An insufficient response by 6 months was a failure to meet the above definition. Non-response was defined as '<50% decrease of serum transaminases within 4 weeks after initiation of treatment'. Remission is defined as liver histology with a Hepatitis Activity Index <4/18. Intolerance to treatment was agreed to stand for 'any adverse event possibly related to treatment leading to potential drug discontinuation'.

CONCLUSIONS

These definitions provide a simple and reproducible framework to define treatment response and non-response, irrespective of the therapeutic intervention. A consensus on endpoints is urgently required to set a global standard for the reporting of study results and to enable inter-study comparisons. Future prospective database studies are needed to validate these endpoints.

LAY SUMMARY

Consensus among international experts on response criteria and endpoints in autoimmune hepatitis is lacking. A consensus on endpoints is urgently required to set a global standard for the reporting of study results and to enable the comparison of results between clinical trials. Therefore, the International Autoimmune Hepatitis Group (IAIHG) herein presents a statement on 5 agreed response criteria and endpoints: complete biochemical response, insufficient response, non-response, remission, and intolerance to treatment, which can be used to guide future reporting.

Introducing an Artificial Deazaflavin Cofactor in Escherichia coli and Saccharomyces cerevisiae

Deazaflavin-dependent whole-cell conversions in well-studied and industrially relevant microorganisms such as Escherichia coli and Saccharomyces cerevisiae have high potential for the biocatalytic production of valuable compounds. The artificial deazaflavin FOP (FO-5′-phosphate) can functionally substitute the natural deazaflavin F₄₂₀ and can be synthesized in fewer steps, offering a solution to the limited availability of the latter due to its complex (bio)synthesis. Herein we set out to produce FOP in vivo as a scalable FOP production method and as a means for FOP-mediated whole-cell conversions. Heterologous expression of the riboflavin kinase from Schizosaccharomyces pombe enabled in vivo phosphorylation of FO, which was supplied by either organic synthesis ex vivo, or by a coexpressed FO synthase in vivo, producing FOP in E. coli as well as in S. cerevisiae. Through combined approaches of enzyme engineering as well as optimization of expression systems and growth media, we further improved the in vivo FOP production in both organisms. The improved FOP production yield in E. coli is comparable to the F₄₂₀ yield of native F₄₂₀-producing organisms such as Mycobacterium smegmatis, but the former can be achieved in a significantly shorter time frame. Our E. coli expression system has an estimated production rate of 0.078 μmol L^–1 h^–1 and results in an intracellular FOP concentration of about 40 μM, which is high enough to support catalysis. In fact, we demonstrate the successful FOP-mediated whole-cell conversion of ketoisophorone using E. coli cells. In S. cerevisiae, in vivo FOP production by SpRFK using supplied FO was improved through media optimization and enzyme engineering. Through structure-guided enzyme engineering, a SpRFK variant with 7-fold increased catalytic efficiency compared to the wild type was discovered. By using this variant in optimized media conditions, FOP production yield in S. cerevisiae was 20-fold increased compared to the very low initial yield of 0.24 ± 0.04 nmol per g dry biomass. The results show that bacterial and eukaryotic hosts can be engineered to produce the functional deazaflavin cofactor mimic FOP.

A Tailor-Made Deazaflavin-Mediated Recycling System for Artificial Nicotinamide Cofactor Biomimetics

Nicotinamide adenine dinucleotide (NAD) and its 2'-phosphorylated form NADP are crucial cofactors for a large array of biocatalytically important redox enzymes. Their high cost and relatively poor stability, however, make them less attractive electron mediators for industrial processes. Nicotinamide cofactor biomimetics (NCBs) are easily synthesized, are inexpensive, and are also generally more stable than their natural counterparts. A bottleneck for the application of these artificial hydride carriers is the lack of efficient cofactor recycling methods. Therefore, we engineered the thermostable F420:NADPH oxidoreductase from Thermobifida fusca (Tfu-FNO), by structure-inspired site-directed mutagenesis, to accommodate the unnatural N1 substituents of eight NCBs. The extraordinarily low redox potential of the natural cofactor F420H2 was then exploited to reduce these NCBs. Wild-type enzyme had detectable activity toward all selected NCBs, with K m values in the millimolar range and k cat values ranging from 0.09 to 1.4 min-1. Saturation mutagenesis at positions Gly-29 and Pro-89 resulted in mutants with up to 139 times higher catalytic efficiencies. Mutant G29W showed a k cat value of 4.2 s-1 toward 1-benzyl-3-acetylpyridine (BAP+), which is similar to the k cat value for the natural substrate NADP+. The best Tfu-FNO variants for a specific NCB were then used for the recycling of catalytic amounts of these nicotinamides in conversion experiments with the thermostable ene-reductase from Thermus scotoductus (TsOYE). We were able to fully convert 10 mM ketoisophorone with BAP+ within 16 h, using F420 or its artificial biomimetic FOP (FO-2'-phosphate) as an efficient electron mediator and glucose-6-phosphate as an electron donor. The generated toolbox of thermostable and NCB-dependent Tfu-FNO variants offers powerful cofactor regeneration biocatalysts for the reduction of several artificial nicotinamide biomimetics at both ambient and high temperatures. In fact, to our knowledge, this enzymatic method seems to be the best-performing NCB-recycling system for BNAH and BAPH thus far.

Flavin-tag: A Facile Method for Site-Specific Labeling of Proteins with a Flavin Fluorophore

Site-specific protein labeling methods are highly valuable tools for research and applications. We present a new protein labeling method that allows covalent attachment of a chromo- and fluorogenic flavin (FMN) to any targeted protein using a short flavinylation peptide-tag. We show that this peptide can be as short as 7 residues and can be located at the N-terminus, C-terminus, or in internal regions of the target protein. Analogous to kinase-catalyzed phosphorylation, the flavin is covalently attached via a stable phosphothreonyl linkage. The site-specific covalent tethering of FMN is accomplished by using a bacterial flavin transferase. The covalent coupling of FMN was shown to work in Escherichia coli and Saccharomyces cerevisiae cells and could be performed in vitro, rendering the "Flavin-tag" method a powerful tool for the selective decoration of proteins with a biocompatible redox-active fluorescent chromophore.

Thermostable D-amino acid decarboxylases derived from Thermotoga maritima diaminopimelate decarboxylase

Diaminopimelate decarboxylases (DAPDCs) are highly selective enzymes that catalyze the common final step in different lysine biosynthetic pathways, i.e. the conversion of meso-diaminopimelate (DAP) to L-lysine. We examined the modification of the substrate specificity of the thermostable decarboxylase from Thermotoga maritima with the aim to introduce activity with 2-aminopimelic acid (2-APA) since its decarboxylation leads to 6-aminocaproic acid (6-ACA), a building block for the synthesis of nylon-6. Structure-based mutagenesis of the distal carboxylate binding site resulted in a set of enzyme variants with new activities toward different D-amino acids. One of the mutants (E315T) had lost most of its activity toward DAP and primarily acted as a 2-APA decarboxylase. We next used computational modeling to explain the observed shift in catalytic activities of the mutants. The results suggest that predictive computational protocols can support the redesign of the catalytic properties of this class of decarboxylating PLP-dependent enzymes.

Molecular characterization of the waxy locus in sorghum

A comparison of approximately 4.5 kb of nucleotide sequence from the waxy locus (the granule-bound starch synthase I [GBSS I] locus) from a waxy line, BTxARG1, and a non-waxy line, QL39, revealed an extremely high level of sequence conservation. Among a total of 24 nucleotide differences and 9 indels, only 2 nucleotide changes resulted in altered amino acid residues. Protein folding prediction software suggested that one of the amino acid changes (Glu to His) may result in an altered protein structure, which may explain the apparently inactive GBSS I present in BTxARG1. This SNP was not found in the second waxy line, RTx2907, which does not produce GBSS I, and no other SNPs or indels were found in the approximately 4 kb of sequence obtained from RTx2907. Using one indel, the waxy locus was mapped to sorghum chromosome SBI-10, which is syntenous to maize chromosome 9; the waxy locus has been mapped to this maize chromosome. The distribution of indels in a diverse set of sorghum germplasm suggested that there are two broad types of non-waxy GBSS I alleles, each type comprising several alleles, and that the two waxy alleles in BTxARG1 and RTx2907 have evolved from one of the non-waxy allele types. The Glu/His polymorphism was found only in BTxARG1 and derived lines and has potential as a perfect marker for the BTxARG1 source of the waxy allele at the GBSS I locus. The indels correctly predicted the non-waxy phenotype in approximately 65% of diverse sorghum germplasm. The indels co-segregated perfectly with phenotype in two sorghum populations derived from crosses between a waxy and a non-waxy sorghum line, correctly identifying heterozygous lines. Thus, these indel markers or sequence-based SNP markers can be used to follow waxy alleles in sorghum breeding programs in selected pedigrees.

Resistance gene analogues in sugarcane and sorghum and their association with quantitative trait loci for rust resistance

Fifty-four different sugarcane resistance gene analogue (RGA) sequences were isolated, characterized, and used to identify molecular markers linked to major disease-resistance loci in sugarcane. Ten RGAs were identified from a sugarcane stem expressed sequence tag (EST) library; the remaining 44 were isolated from sugarcane stem, leaf, and root tissue using primers designed to conserved RGA motifs. The map location of 31 of the RGAs was determined in sugarcane and compared with the location of quantitative trait loci (QTL) for brown rust resistance. After 2 years of phenotyping, 3 RGAs were shown to generate markers that were significantly associated with resistance to this disease. To assist in the understanding of the complex genetic structure of sugarcane, 17 of the 31 RGAs were also mapped in sorghum. Comparative mapping between sugarcane and sorghum revealed syntenic localization of several RGA clusters. The 3 brown rust associated RGAs were shown to map to the same linkage group (LG) in sorghum with 2 mapping to one region and the third to a region previously shown to contain a major rust-resistance QTL in sorghum. These results illustrate the value of using RGAs for the identification of markers linked to disease resistance loci and the value of simultaneous mapping in sugarcane and sorghum.

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class ( Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue ( Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.

Identifications of two different mechanisms for sorghum midge resistance through QTL mapping

Sorghum midge is the one of the most damaging insect pests of grain sorghum production worldwide. At least three different mechanisms are involved in midge resistance. The genetic bases of these mechanisms, however, are poorly understood. In this study, for the first time, quantitative trait loci associated with two of the mechanisms of midge resistance, antixenosis and antibiosis, were identified in an RI (recombinant inbred) population from the cross of sorghum lines ICSV745 x 90562. Two genetic regions located on separate linkage groups were found to be associated with antixenosis and explained 12% and 15%, respectively, of the total variation in egg numbers/spikelet laid in a cage experiment. One region was significantly associated with antibiosis and explained 34.5% of the variation of the difference of egg and pupal counts in the RI population. The identification of genes for different mechanisms of midge resistance will be particularly useful for exploring new sources of midge resistance and for gene pyramiding of different mechanisms for increased security in sorghum breeding through marker-assisted selection.

Direct sequencing of RAPD fragments using 3'-extended oligonucleotide primers and dye terminator cycle-sequencing

Random amplified polymorphic DNA (RAPD) markers are used widely to develop high resolution genetic maps and for genome fingerprinting. Typically, single oligomers of approximately 10 nucleotides are used to PCR amplify characteristic RAPD marker fragments. We describe an efficient method for the direct end-sequencing of gel-purified RAPD fragments using one primer from a set of four 3'-terminal extended (A, T, C or G) oligonucleotides, identical to the RAPD primer but for the single nucleotide extension. Strand-specific DNA sequence could be independently read from each of the RAPD fragments without recourse to strand separation or fragment cloning. Informative RAPD fragments could be readily converted into mapped STS or SCAR loci using this technology. The 3'-extended primers may also be used to amplify independent genomic RAPD markers.

Nucleation in protein crystallization

In this paper protein crystallization is regarded as a process starting with phase separation in a two-component system. The nucleation time of a lysozyme solution is measured by recording the NMR spectra of crystallizing solutions as a function of time. It is found that there is an appreciable induction time before the first nuclei appear in the solution and that this induction time depends strongly on the degree of supersaturation due to the protein concentration at a given ionic strength or due to the temperature. From the experimental data it is evident that (at least for lysozyme) crystallization under the prevailing experimental conditions is a transient process with an induction time and not a steady-state process.

A 200 bp region of the pea ENOD12 promoter is sufficient for nodule-specific and nod factor induced expression

ENOD12 is one of the first nodulin genes expressed upon inoculation with Rhizobium and also purified Nod factors are able to induce ENOD12 expression. The ENOD12 gene family in pea (Pisum sativum) has two members. A cDNA clone representing PsENOD12A [26] and a PsENOD12B genomic clone [7] have been previously described. The isolation and characterization of a PsENOD12A genomic clone is presented in this paper. By using a Vicia hirsuta-Agrobacterium rhizogenes transformation system it is shown that both genes have a similar expression pattern in transgenic V. hirsuta root nodules. Promoter analyses of both PsENOD12 promoters showed that the 200 bp immediately upstream of the transcription start are sufficient to direct nodule-specific and Nod factor-induced gene expression.

Microhydrodynamics simulation of protein crystallization. I. Static calculations

A computer simulation method is proposed to study the effects of hydrodynamic interactions on protein crystallization. It is a combination of Stokesian dynamics and continuum hydrodynamics and is referred to as "microhydrodynamics." The method is checked against analytical expressions for Stokes drag and diffusion coefficients for unit spheres. For a number of protein molecules the diffusion coefficients have been calculated and compared with experimental values. It is shown that the method works well for stationary calculations. Using dynamical calculations interacting protein molecules will be simulated to study the events in the early stages of protein crystallization.

An optimal strategy for X-ray data collection on macromolecular crystals with position-sensitive detectors

X-ray data collection on macromolecular crystals is preferably done with minimum exposure time and high completeness. A Fortran procedure – DCS – has been written in the environment of the MADNES program to predict the completeness of data before the start of actual data collection. In addition, the program can check and refine the data-collection strategy and suggest optimal settings and rotation ranges for one or more crystals in different orientations that give highest completeness in minimum exposure time, thus extending the life of the crystal. The method has been tested on previously collected data as well as on new structures. Since the program has been in full use in this laboratory, the completeness of collected data has improved significantly.

Mesoscopic theories for protein crystal growth

A computer-simulation method is proposed for studying the hydrodynamic interactions of rigid protein molecules. It is a combination of Stokes dynamics and continuum hydrodynamics. The Stokes equations of motion for the protein molecules, the creeping-flow equation for the solvent together with the no-slip boundary conditions give a complete representation of the system. The resulting three-dimensional boundary-value problem can be rewritten in a two-dimensional form (without any loss of information) considering the surfaces of the particles only. Then, by solving the equations on discrete surface elements, the so-called mobility matrix is determined in which all hydrodynamic interactions are included. Finally, after calculation of the conservative forces and the stochastic force, the new velocities of the protein molecules can be determined. The simulation method can be applied to arbitrary particle shapes. It can also handle arbitrary flow fields, and the effects of applying a flow field to the system can be studied. From analysis of the trajectories, information can be gained on the kinetics and thermodynamics in the early stages of the crystallization process.

Crystal structure of a porcine pancreatic phospholipase A2 mutant. A large conformational change caused by the F63V point mutation

The highly homologous bovine and porcine pancreatic phospholipase A2 (85% amino acid residue identity) show a large conformational difference in the loop from residue 59 to 71. In bovine phospholipase A2 residues 59 to 66 adopt an alpha-helix conformation, while residues 67 to 71 are in a surface loop. Residues 59 to 66 in the porcine enzyme have a random coil conformation, and residues 67 to 71 form a short 3(10)-helix. It has been suggested that most probably this conformational difference is caused by the substitution Val63 (bovine) to Phe63 (porcine) in the otherwise invariant loop 59 to 70. To test this hypothesis, a mutant porcine phospholipase A2 was constructed in which residue Phe63 was replaced by a Val. The activity of this F63V mutant towards aggregated substrates was about half the activity of wild-type porcine phospholipase A2, but significantly different from that of the bovine enzyme. The affinity for zwitterionic interfaces was found to be intermediate between porcine and bovine phospholipase. The mutation did not have any effect on the stability of the enzyme towards denaturation by guanidine.HCl. The F63V mutant was crystallized in space group P2(1)2(1)2(1) with cell dimensions a = 79.88 A, b = 65.23 A, c = 52.62 A, with two molecules per asymmetric unit. Its three-dimensional structure was solved by molecular replacement methods, and refined to a crystallographic R-factor of 17.6% for all data between 10 and 2.2 A resolution. In one molecule the 58 to 71 loop is in very weak density, suggesting a high degree of disorder or flexibility. The conformation of the same loop in the other molecule could be determined unambiguously. It shows a conformation which resembles more that of bovine phospholipase A2 than that of porcine phospholipase. It is concluded that indeed the single F63V substitution causes a dramatic conformational change.

Crystal structure of the reduced form of p-hydroxybenzoate hydroxylase refined at 2.3 A resolution

The crystal structure of the reduced form of the enzyme p-hydroxybenzoate hydroxylase from Pseudomonas fluorescens, complexed with its substrate p-hydroxybenzoate, has been obtained by protein X-ray crystallography. Crystals of the reduced form were prepared by soaking crystals of the oxidized enzyme-substrate complex in deaerated mother liquor containing 300-400 mM NADPH. A rapid bleaching of the crystals indicated the reduction of the enzyme-bound FAD by NADPH. This was confirmed by single crystal spectroscopy. X-ray data to 2.3 A were collected on oscillation films using a rotating anode generator as an X-ray source. After data processing and reduction, restrained least squares refinement using the 1.9 A structure of the oxidized enzyme-substrate complex as a starting model, yielded a crystallographic R-factor of 14.8% for 11,394 reflections. The final model of the reduced complex contains 3,098 protein atoms, the FAD molecule, the substrate p-hydroxybenzoate and 322 solvent molecules. The structures of the oxidized and reduced forms of the enzyme-substrate complex were found to be very similar. The root-mean-square discrepancy for all atoms between both structures was 0.38 A. The flavin ring is almost completely planar in the final model, although it was allowed to bend or twist during refinement. The observed angle between the benzene and the pyrimidine ring is 2 degrees. This value should be compared with observed values of 10 degrees for the oxidized enzyme-substrate complex and 19 degrees for the enzyme-product complex. The position of the substrate is virtually unaltered with respect to its position in the oxidized enzyme. No trace of a bound NADP+ or NADPH molecule was found.

Site-directed mutagenesis and X-ray crystallography of two phospholipase A2 mutants: Y52F and Y73F

Tyr52 and Tyr73 are conserved amino acid residues throughout all vertebrate phospholipases A2. They are part of an extended hydrogen bonding system that links the N-terminal alpha-NH3(+)-group to the catalytic residues His48 and Asp99. These tyrosines were replaced by phenylalanines in a porcine pancreatic phospholipase A2 mutant, in which residues 62-66 had been deleted (delta 62-66PLA2). The mutations did not affect the catalytic properties of the enzyme, nor the folding kinetics. The stability against denaturation by guanidine hydrochloride was decreased, however. To analyse how the enzyme compensates for the loss of the tyrosine hydroxyl group, the X-ray structures of the delta Y52F and delta Y73F mutants were determined. After crystallographic refinement the final crystallographic R-factors were 18.1% for the delta Y52F mutant (data between 7 and 2.3 A resolution) and 19.1% for the delta Y73F mutant (data between 7 and 2.4 A resolution). No conformational changes occurred in the mutants compared with the delta 62-66PLA2, but an empty cavity formed at the site of the hydroxyl group of the former tyrosine. In both mutants the Asp99 side chain loses one of its hydrogen bonds and this might explain the observed destabilization.

Structure determination of quinoprotein methylamine dehydrogenase from Thiobacillus versutus

The crystal structure of quinoprotein methylamine dehydrogenase from Thiobacillus versutus (EC 1.4.99.3, Mr = 123,500) has been solved to 2.25 A resolution. The crystals of space group P3(1)21 (a = b = 129.8, c = 104.3 A) contain half a tetrameric enzyme molecule in the asymmetric unit, with a solvent content of ca 70%. The procedure used to solve this structure involved multiple isomorphous-replacement phasing, complemented by phase extension using solvent flattening, and phase combination with partial-model phases. The use of solvent flattening was essential to generate good quality electron density maps into which initial models were built. These partial models were refined using molecular-dynamics procedures. Refined model phases were then combined with solvent-flattening phases to generate improved electron density distributions. In the absence of an amino-acid sequence for this enzyme, the current 2.25 A resolution electron density map was interpreted to provide a model for the complete molecule. The crystallographic R factor for this model, which lacks any water molecules, is 28.6% for data between 6.0 and 2.25 A resolution.

Structure of an engineered porcine phospholipase A2 with enhanced activity at 2.1 A resolution. Comparison with the wild-type porcine and Crotalus atrox phospholipase A2

The crystal structure of an engineered phospholipase A2 with enhanced activity has been refined to an R-factor of 18.6% at 2.1 A resolution using a combination of molecular dynamics refinement by the GROMOS package and least-squares refinement by TNT. This mutant phospholipase was obtained previously by deleting residues 62 to 66 in porcine pancreatic phospholipase A2, and changing Asp59 to Ser, Ser60 to Gly and Asn67 to Tyr. The refined structure allowed a detailed comparison with wild-type porcine and Crotalus atrox phospholipase A2. The conformation of the deletion region appears to be intermediate between that in those two enzymes. The residues in the active center are virtually the same. An internal hydrophobic area occupied by Phe63 in the wild-type porcine phospholipase A2 is kept as conserved as possible by local rearrangement of neighboring atoms. In the mutant structure, this hydrophobic pocket is now occupied by the disulfide bond between residues 61 and 91. A detailed description of the second binding site for a calcium ion in this enzyme is given.

X-ray structure of phospholipase A2 complexed with a substrate-derived inhibitor

Phospholipases A2 play a part in a number of physiologically important cellular processes such as inflammation, blood platelet aggregation and acute hypersensitivity. These processes are all initiated by the release of arachidonic acid from cell membranes which is catalysed by intracellular phospholipases A2 and followed by conversion of arachidonic acid to prostaglandins, leukotrienes or thromboxanes. An imbalance in the production of these compounds can lead to chronic inflammatory diseases such as rheumatoid arthritis and asthma. Inhibitors of phospholipase A2 might therefore act to reduce the effects of inflammation, so structural information about the binding of phospholipase A2 to its substrates could be helpful in the design of therapeutic drugs. The three-dimensional structure is not known for any intracellular phospholipase A2, but these enzymes share significant sequence homology with secreted phospholipases, for which some of the structures have been determined. Here we report the structure of a complex between an extracellular phospholipase A2 and a competitively inhibiting substrate analogue, which reveals considerable detail about the interaction and suggests a mechanism for catalysis by this enzyme.

Analysis of the active site of the flavoprotein p-hydroxybenzoate hydroxylase and some ideas with respect to its reaction mechanism

The flavoprotein p-hydroxybenzoate hydroxylase has been studied extensively by biochemical techniques by others and in our laboratory by X-ray crystallography. As a result of the latter investigations, well-refined crystal structures are known of the enzyme complexed (i) with its substrate p-hydroxybenzoate and (ii) with its reaction product 3,4-dihydroxybenzoate and (iii) the enzyme with reduced FAD. Knowledge of these structures and the availability of the three-dimensional structure of a model compound for the reactive flavin 4a-hydroperoxide intermediate has allowed a detailed analysis of the reaction with oxygen. In the model of this reaction intermediate, fitted to the active site of p-hydroxybenzoate hydroxylase, all possible positions of the distal oxygen were surveyed by rotating this oxygen about the single bond between the C4a and the proximal oxygen. It was found that the distal oxygen is free to sweep an arc of about 180 degrees in the active site. The flavin 4a-peroxide anion, which is formed after reaction of molecular oxygen with reduced FAD, might accept a proton from an active-site water molecule or from the hydroxyl group of the substrate. The position of the oxygen to be transferred with respect to the substrate appears to be almost ideal for nucleophilic attack of the substrate onto this oxygen. The oxygen is situated above the 3-position of the substrate where the substitution takes place, at an angle of about 60 degrees with the aromatic plane, allowing strong interactions with the pi electrons of the substrate. Polarization of the peroxide oxygen-oxygen bond by the enzyme may enhance the reactivity of flavin 4a-peroxide.

The coenzyme analogue adenosine 5-diphosphoribose displaces FAD in the active site of p-hydroxybenzoate hydroxylase. An x-ray crystallographic investigation

p-Hydroxybenzoate hydroxylase (PHBH) is an NADPH-dependent enzyme. To locate the NADPH binding site, the enzyme was crystallized under anaerobic conditions in the presence of the substrate p-hydroxybenzoate, the coenzyme analogue adenosine 5-diphosphoribose (ADPR), and sodium dithionite. This yielded colorless crystals that were suitable for X-ray analysis. Diffraction data were collected up to 2.7-A resolution. A difference Fourier between data from these colorless crystals and data from yellow crystals of the enzyme-substrate complex showed that in the colorless crystals the flavin ring was absent. The adenosine 5'-diphosphate moiety, which is the common part between FAD and ADPR, was still present. After restrained least-squares refinement of the enzyme-substrate complex with the riboflavin omitted from the model, additional electron density appeared near the pyrophosphate, which indicated the presence of an ADPR molecule in the FAD binding site of PHBH. The complete ADPR molecule was fitted to the electron density, and subsequent least-squares refinement resulted in a final R factor of 16.8%. Replacement of bound FAD by ADPR was confirmed by equilibrium dialysis, where it was shown that ADPR can effectively remove FAD from the enzyme under mild conditions in 0.1 M potassium phosphate buffer, pH 8.0. The empty pocket left by the flavin ring is filled by solvent, leaving the architecture of the active site and the binding of the substrate largely unaffected.

Crystal structure of the p-hydroxybenzoate hydroxylase-substrate complex refined at 1.9 A resolution. Analysis of the enzyme-substrate and enzyme-product complexes

Using synchrotron radiation, the X-ray diffraction intensities of crystals of p-hydroxy-benzoate hydroxylase, complexed with the substrate p-hydroxybenzoate, were measured to a resolution of 1.9 A. Restrained least-squares refinement alternated with rebuilding in electron density maps yielded an atom model of the enzyme-substrate complex with a crystallographic R-factor of 15.6% for 31,148 reflections between 6.0 and 1.9 A. A total of 330 solvent molecules was located. In the final model, only three residues have deviating phi-psi angle combinations. One of them, the active site residue Arg44, has a well-defined electron density and may be strained to adopt this conformation for efficient catalysis. The mode of binding of FAD is distinctly different for the different components of the coenzyme. The adenine ring is engaged in three water-mediated hydrogen bonds with the protein, while making only one direct hydrogen bond with the enzyme. The pyrophosphate moiety makes five water-mediated versus three direct hydrogen bonds. The ribityl and ribose moieties make only direct hydrogen bonds, in all cases, except one, with side-chain atoms. The isoalloxazine ring also makes only direct hydrogen bonds, but virtually only with main-chain atoms. The conformation of FAD in p-hydroxybenzoate hydroxylase is strikingly similar to that in glutathione reductase, while the riboflavin-binding parts of these two enzymes have no structural similarity at all. The refined 1.9 A structure of the p-hydroxybenzoate hydroxylase-substrate complex was the basis of further refinement of the 2.3 A structure of the enzyme-product complex. The result was a final R-factor of 16.7% for 14,339 reflections between 6.0 and 2.3 A and an improved geometry. Comparison between the complexes indicated only small differences in the active site region, where the product molecule is rotated by 14 degrees compared with the substrate in the enzyme-substrate complex. During the refinements of the enzyme-substrate and enzyme-product complexes, the flavin ring was allowed to bend or twist by imposing planarity restraints on the benzene and pyrimidine ring, but not on the flavin ring as a whole. The observed angle between the benzene ring and the pyrimidine ring was 10 degrees for the enzyme-substrate complex and 19 degrees for the enzyme-product complex. Because of the high temperature factors of the flavin ring in the enzyme-product complex, the latter value should be treated with caution. Six out of eight peptide residues near the flavin ring are oriented with their nitrogen atom pointing towards the ring.(ABSTRACT TRUNCATED AT 400 WORDS)

Structure of quinoprotein methylamine dehydrogenase at 2.25 A resolution

The three-dimensional structure of quinoprotein methylamine dehydrogenase from Thiobacillus versutus has been determined at 2.25 A resolution by a combination of multiple isomorphous replacement, phase extension by solvent flattening and partial structure phasing using molecular dynamics refinement. In the resulting map, the polypeptide chain for both subunits could be followed and an X-ray sequence was established. The tetrameric enzyme, made up of two heavy (H) and two light (L) subunits, is a flat parallellepiped with overall dimensions of approximately 76 x 61 x 45 A. The H subunit, comprising 370 residues, is made up of two distinct segments: the first 31 residues form an extension which embraces one of the L subunits; the remaining residues are found in a disc-shaped domain. This domain is formed by a circular arrangement of seven topologically identical four-stranded antiparallel beta-sheets, with approximately 7-fold symmetry. In spite of distinct differences, this arrangement is reminiscent of the structure found in influenza virus neuraminidase. The L subunit consists of 121 residues, out of which 53 form a beta-sheet scaffold of a central three-stranded antiparallel sheet flanked by two shorter two-stranded antiparallel sheets. The remaining residues are found in segments of irregular structure. This subunit is stabilized by six disulphide bridges, plus two covalent bridges involving the quinone co-factor and residues 57 and 107 of this subunit. The active site is located in a channel at the interface region between the H and L subunits, and the electron density in this part of the molecule suggests that the co-factor of this enzyme is not pyrrolo quinoline quinone (PQQ) itself, but might be instead a precursor of PQQ.

Enhanced activity and altered specificity of phospholipase A2 by deletion of a surface loop

Protein engineering and x-ray crystallography have been used to study the role of a surface loop that is present in pancreatic phospholipases but is absent in snake venom phospholipases. Removal of residues 62 to 66 from porcine pancreatic phospholipase A2 does not change the binding constant for micelles significantly, but it improves catalytic activity up to 16 times on micellar (zwitterionic) lecithin substrates. In contrast, the decrease in activity on negatively charged substrates is greater than fourfold. A crystallographic study of the mutant enzyme shows that the region of the deletion has a well-defined structure that differs from the structure of the wild-type enzyme. No structural changes in the active site of the enzyme were detected.

X-ray structure of lipoamide dehydrogenase from Azotobacter vinelandii determined by a combination of molecular and isomorphous replacement techniques

The crystal structure of lipoamide dehydrogenase from Azotobacter vinelandii has been determined by a combination of molecular replacement and isomorphous replacement techniques yielding eventually a good-quality 2.8 A electron density map. Initially, the structure determination was attempted by molecular replacement procedures alone using a model of human glutathione reductase, which has 26% sequence identity with this bacterial dehydrogenase. The rotation function yielded the correct orientation of the model structure both when the glutathione reductase dimer and monomer were used as starting model. The translation function could not be solved, however. Consequently, data for two heavy-atom derivatives were collected using the Hamburg synchotron facilities. The derivatives had several sites in common, which was presumably a major reason why the electron density map obtained by isomorphous information alone was of poor quality. Application of solvent flattening procedures cleaned up the map considerably, however, showing clearly the outline of the lipoamide dehydrogenase dimer, which has a molecular weight of 100,000. Application of the "phased translation function", which combines the phase information of both isomorphous and molecular replacement, led to an unambiguous determination of the position of the model structure in the lipoamide dehydrogenase unit cell. The non-crystallographic 2-fold axis of the dimer was optimized by several cycles of constrained-restrained least-squares refinement and subsequently used for phase improvement by 2-fold density averaging. After ten cycles at 3.5 A, the resolution was gradually extended to 2.8 A in another 140 cycles. The 2.8 A electron density distribution obtained in this manner was of much improved quality and allowed building of an atomic model of A. vinelandii lipoamide dehydrogenase. It appears that in the orthorhombic crystals used each dimer is involved in contacts with eight surrounding dimers, leaving unexplained why the crystals are rather fragile. Contacts between subunits within one dimer, which are quite extensive, can be divided into two regions separated by a cavity. In one of the contact regions, the level of sequence identity with glutathione reductase is very low but it is quite high in the other. The folding of the polypeptide chain in each subunit is quite similar to that of glutathione reductase, as is the extended conformation of the co-enzyme FAD.(ABSTRACT TRUNCATED AT 400 WORDS)

The influence of purification and protein heterogeneity on the crystallization of p-hydroxybenzoate hydroxylase

The structure of the enzyme p-hydroxybenzoate hydroxylase was determined to a resolution of 0.25 nm [Wierenga et al. (1979) J. Mol. Biol. 131, 53-73] with crystals belonging to space group C222(1). Subsequently it was impossible to repeat the growth of this crystal form and only poor quality tetragonal crystals could be obtained. We have thoroughly investigated this problem and found that Cibacron-blue-purified enzyme appears to be heterogeneous with respect to aggregation state and Cys-116 oxidation. Most importantly, it could be firmly established that C222(1) crystals can only be grown from purely dimeric p-hydroxybenzoate hydroxylase possessing an intact SH group. Ion-exchange chromatography on DEAE-Sepharose can successfully remove those forms of the enzyme which impede successful crystallization. Sulfite and dithiothreitol improve crystallization by dissociating the enzyme oligomers into dimers; sulfite especially gives excellent results.

Crystal structure of bovine pancreatic phospholipase A2 covalently inhibited by p-bromo-phenacyl-bromide

Bovine pancreatic phospholipase A2 covalently inhibited by p-bromo-phenacyl-bromide was crystallized from 50% (v/v) 2-methyl-2,4-pentanediol. The space group was P3(1)21 with cell dimensions a = b = 46.73 A and c = 102.5 A (1 A = 0.1 nm). Diffraction data were collected by oscillation photography from one single crystal of dimensions 0.2 mm x 0.2 mm x 0.2 mm. The crystal structure was determined to a resolution of 2.5 A by crystallographic refinement of a starting model, which consisted of native bovine pancreatic phospholipase A2 positioned and oriented in the P3(1)21 cell as in the bovine pro-phospholipase A2. The crystallographic R-factor decreased from 0.378 to 0.197 after 70 refinement cycles. For the greater part the three-dimensional structure was very similar to that of native phospholipase. The inhibitor group shows up clearly. However, as in solution, there is no calcium ion bound any more in the active site, and this causes a significant conformational change in the loop from residue 59 to 73. This loop is remote from the calcium binding site. Interestingly, this is the same loop that also shows different conformations in other phospholipase A2 molecules. The inhibitor molecule has hydrophobic interactions with Phe5 and Cys45. Rational design of specific and potent inhibitors of phospholipase A2 catalysis is discussed on the basis of the present three-dimensional structure.

Molecular modeling reveals the possible importance of a carbonyl oxygen binding pocket for the catalytic mechanism of p-hydroxybenzoate hydroxylase

p-Hydroxybenzoate hydroxylase catalyzes the hydroxylation of an aromatic substrate and uses flavin as a cofactor. The reaction probably occurs via a flavin 4a-hydroperoxide intermediate. In this study the crystal structure of 4a,5-epoxyethano-3-methyl-4a,5-dihydrolumiflavin, an analogue of the flavin 4a-hydroperoxide intermediate, was fitted to the active site in the crystal structure of the p-hydroxybenzoate hydroxylase-3,4-dihydroxybenzoate complex. This model of an important catalytic intermediate fitted very well in the active site of p-hydroxybenzoate hydroxylase. The most striking result was that whereas with the normal flavin, the 0-4 of the flavin ring makes only poor hydrogen bonds with the protein, with the flavin 4a-hydroperoxide analogue, the same 0-4 makes strong hydrogen bonds with the NH groups of Gly-46 and Val-47. These two NH groups form a carbonyl oxygen binding pocket which has a geometry almost identical to the oxyanion hole found in several proteases. The possible consequences of this model for the reaction mechanism of p-hydroxybenzoate hydroxylase are discussed.

Crystal structure of p-hydroxybenzoate hydroxylase complexed with its reaction product 3,4-dihydroxybenzoate

Crystals of the flavin-containing enzyme p-hydroxybenzoate hydroxylase (PHBHase) complexed with its reaction product were investigated in order to obtain insight into the catalytic cycle of this enzyme involving two substrates and two cofactors. PHBHase was crystallized initially with its substrate, p-hydroxybenzoate and the substrate was then converted into the product 3,4-dihydroxybenzoate by allowing the catalytic reaction to proceed in the crystals. In addition, crystals were soaked in mother liquor containing a high concentration of this product. Data up to 2.3 A (1 A = 0.1 nm) were collected by the oscillation method and the structure of the enzyme product complex was refined by alternate restrained least-squares procedures and model building by computer graphics techniques. A total of 273 solvent molecules could be located, four of them being presumably sulfate ions. The R-factor for 14,339 reflections between 6.0 A and 2.3 A is 19.3%. The 3-hydroxyl group of the product introduced by the enzyme is clearly visible in the electron density, showing unambiguously which carbon atom of the substrate is hydroxylated. A clear picture of the hydroxylation site is obtained. The plane of the product is rotated 21 degrees with respect to the plane of the substrate in the current model of enzyme-substrate complex. The 4-hydroxyl group of the product is hydrogen bonded to the hydroxyl group of Tyr201, its carboxyl group is interacting with the side-chains of Tyr222, Arg214 and Ser212, while the newly introduced 3-hydroxyl group makes a hydrogen bond with the backbone carbonyl oxygen of Pro293.

Purification, crystallization and preliminary X-ray investigation of quinoprotein methylamine dehydrogenase from Thiobacillus versutus

The enzyme methylamine dehydrogenase or primary-amine:(acceptor) oxidoreductase (deaminating) (EC 1.4.99.3) was purified from the bacterium Thiobacillus versutus to homogeneity, as judged by polyacrylamide gel electrophoresis. The native enzyme has a Mr of 123 500 and contains four subunits arranged in a alpha 2 beta 2 configuration, the light and heavy subunits having a Mr of 12900 and 47500 respectively. The isoelectric point is 3.9. The purified enzyme was crystallized from 37--42% saturated ammonium sulphate in 0.1 M sodium acetate buffer, pH 5.0. The space group is P3(1)21 or P3(2)21, with one alpha 2 beta 2 molecule in the asymmetric unit. The cell dimensions are: a = b = 13.01 nm; c = 10.40 nm. The X-ray diffraction pattern extends to at least 0.25-nm resolution.

Heat-labile enterotoxin of Escherichia coli. Characterization of different crystal forms

Heat-labile enterotoxin (LT) was obtained in large quantities (several-gram amounts) and great purity from Escherichia coli C600 carrying the LT-coding multicopy plasmid EWD299. By growing this strain on a medium that allows high cell densities in the early stationary phase, we increased the net LT production per milliliter by a factor of 200, compared to natural porcine enterotoxigenic E. coli. Adsorption and redesorption on Controlled-Pore Glass usually resulted in a 50-100-fold purification of LT in one simple step, which was followed by established purification procedures. LT shows a natural tendency to form large crystals, which, however, are disordered. After numerous trials, conditions were found which virtually eliminated the disorder. Much better crystals were obtained by employing CdCl2 and KF as coprecipitating agents. CdCl2 yielded crystals which did not contain A subunits as judged by electrophoresis of dissolved crystals. Adding KF results in beautifully shaped crystals which diffracted beyond 2 A and are suitable for a high resolution structure determination.