Optimisation of the surface electrostatics as a strategy for cold adaptation of uracil-DNA N-glycosylase (UNG) from Atlantic cod (Gadus morhua)

Characterization of Heat-labile Uracil-DNA Glycosylase from Oncorhynchus mykiss and its Application for Carry-over Contamination Control in RT-qPCR

BACKGROUND

Heat-labile uracil-DNA glycosylase (HL-UDG) is commonly employed to eliminate carry-over contamination in DNA amplifications. However, the prevailing HL-UDG is markedly inactivated at 50°C, rendering it unsuitable for specific one-step RT-qPCR protocols utilizing reverse transcriptase at an optimal temperature of 42°C.

OBJECTIVE

This study aimed to explore novel HL-UDG with lower inactivation temperature and for recombinant expression.

METHODS

The gene encoding an HL-UDG was cloned from the cold-water fish rainbow trout (Oncorhynchus mykiss) and expressed in Escherichia coli with high yield. The thermostability of this enzyme and other enzymatic characteristics were thoroughly examined. The novel HL-UDG was then applied for controlling carry-over contamination in one-step RT-qPCR.

RESULTS

This recombinantly expressed truncated HL-UDG of rainbow trout (OmUDG) exhibited high amino acids similarity (84.1% identity) to recombinant Atlantic cod UDG (rcUDG) and was easily denatured at 40°C. The optimal pH of OmUDG was 8.0, and the optimal concentrations of both Na+ and K+ were 10 mM. Since its inactivation temperature was lower than that of rcUDG, the OmUDG could be used to eliminate carry-over contamination in one-step RT-qPCR with moderate reverse transcription temperature.

CONCLUSION

We successfully identified and recombinantly expressed a novel HL-UDG with an inactivation temperature of 40°C. It is suitable for eliminating carry-over contamination in one-step RT-qPCR.

A New Class of Uracil-DNA Glycosylase Inhibitors Active against Human and Vaccinia Virus Enzyme

Uracil-DNA glycosylases are enzymes that excise uracil bases appearing in DNA as a result of cytosine deamination or accidental dUMP incorporation from the dUTP pool. The activity of Family 1 uracil-DNA glycosylase (UNG) activity limits the efficiency of antimetabolite drugs and is essential for virulence in some bacterial and viral infections. Thus, UNG is regarded as a promising target for antitumor, antiviral, antibacterial, and antiprotozoal drugs. Most UNG inhibitors presently developed are based on the uracil base linked to various substituents, yet new pharmacophores are wanted to target a wide range of UNGs. We have conducted virtual screening of a 1,027,767-ligand library and biochemically screened the best hits for the inhibitory activity against human and vaccinia virus UNG enzymes. Although even the best inhibitors had IC₅₀ ≥ 100 μM, they were highly enriched in a common fragment, tetrahydro-2,4,6-trioxopyrimidinylidene (PyO3). In silico, PyO3 preferably docked into the enzyme's active site, and in kinetic experiments, the inhibition was better consistent with the competitive mechanism. The toxicity of two best inhibitors for human cells was independent of the presence of methotrexate, which is consistent with the hypothesis that dUMP in genomic DNA is less toxic for the cell than strand breaks arising from the massive removal of uracil. We conclude that PyO3 may be a novel pharmacophore with the potential for development into UNG-targeting agents.

A structurally conserved motif in γ-herpesvirus uracil-DNA glycosylases elicits duplex nucleotide-flipping

Efficient γ-herpesvirus lytic phase replication requires a virally encoded UNG-type uracil-DNA glycosylase as a structural element of the viral replisome. Uniquely, γ-herpesvirus UNGs carry a seven or eight residue insertion of variable sequence in the otherwise highly conserved minor-groove DNA binding loop. In Epstein-Barr Virus [HHV-4] UNG, this motif forms a disc-shaped loop structure of unclear significance. To ascertain the biological role of the loop insertion, we determined the crystal structure of Kaposi's sarcoma-associated herpesvirus [HHV-8] UNG (kUNG) in its product complex with a uracil-containing dsDNA, as well as two structures of kUNG in its apo state. We find the disc-like conformation is conserved, but only when the kUNG DNA-binding cleft is occupied. Surprisingly, kUNG uses this structure to flip the orphaned partner base of the substrate deoxyuridine out of the DNA duplex while retaining canonical UNG deoxyuridine-flipping and catalysis. The orphan base is stably posed in the DNA major groove which, due to DNA backbone manipulation by kUNG, is more open than in other UNG-dsDNA structures. Mutagenesis suggests a model in which the kUNG loop is pinned outside the DNA-binding cleft until DNA docking promotes rigid structuring of the loop and duplex nucleotide flipping, a novel observation for UNGs.

Structure determination of uracil-DNA N-glycosylase from Deinococcus radiodurans in complex with DNA

Uracil-DNA N-glycosylase (UNG) is a DNA-repair enzyme in the base-excision repair (BER) pathway which removes uracil from DNA. Here, the crystal structure of UNG from the extremophilic bacterium Deinococcus radiodurans (DrUNG) in complex with DNA is reported at a resolution of 1.35 Å. Prior to the crystallization experiments, the affinity between DrUNG and different DNA oligonucleotides was tested by electrophoretic mobility shift assays (EMSAs). As a result of this analysis, two 16 nt double-stranded DNAs were chosen for the co-crystallization experiments, one of which (16 nt AU) resulted in well diffracting crystals. The DNA in the co-crystal structure contained an abasic site (substrate product) flipped into the active site of the enzyme, with no uracil in the active-site pocket. Despite the high resolution, it was not possible to fit all of the terminal nucleotides of the DNA complex into electron density owing to disorder caused by a lack of stabilizing interactions. However, the DNA which was in contact with the enzyme, close to the active site, was well ordered and allowed detailed analysis of the enzyme-DNA interaction. The complex revealed that the interaction between DrUNG and DNA is similar to that in the previously determined crystal structure of human UNG (hUNG) in complex with DNA [Slupphaug et al. (1996). Nature (London), 384, 87-92]. Substitutions in a (here defined) variable part of the leucine loop result in a shorter loop (eight residues instead of nine) in DrUNG compared with hUNG; regardless of this, it seems to fulfil its role and generate a stabilizing force with the minor groove upon flipping out of the damaged base into the active site. The structure also provides a rationale for the previously observed high catalytic efficiency of DrUNG caused by high substrate affinity by demonstrating an increased number of long-range electrostatic interactions between the enzyme and the DNA. Interestingly, specific interactions between residues in the N-terminus of a symmetry-related molecule and the complementary DNA strand facing away from the active site were also observed which seem to stabilize the enzyme-DNA complex. However, the significance of this observation remains to be investigated. The results provide new insights into the current knowledge about DNA damage recognition and repair by uracil-DNA glycosylases.

Structural plasticity inMycobacterium tuberculosisuracil-DNA glycosylase (MtUng) and its functional implications

17 independent crystal structures of family I uracil-DNA glycosylase fromMycobacterium tuberculosis(MtUng) and its complexes with uracil and its derivatives, distributed among five distinct crystal forms, have been determined. Thermodynamic parameters of binding in the complexes have been measured using isothermal titration calorimetry. The two-domain protein exhibits open and closed conformations, suggesting that the closure of the domain on DNA binding involves conformational selection. Segmental mobility in the enzyme molecule is confined to a 32-residue stretch which plays a major role in DNA binding. Uracil and its derivatives can bind to the protein in two possible orientations. Only one of them is possible when there is a bulky substituent at the 5′ position. The crystal structures of the complexes provide a reasonable rationale for the observed thermodynamic parameters. In addition to providing fresh insights into the structure, plasticity and interactions of the protein molecule, the results of the present investigation provide a platform for structure-based inhibitor design.

Structural and biophysical analysis of interactions between cod and human uracil-DNA N-glycosylase (UNG) and UNG inhibitor (Ugi)

Uracil-DNA N-glycosylase from Atlantic cod (cUNG) shows cold-adapted features such as high catalytic efficiency, a low temperature optimum for activity and reduced thermal stability compared with its mesophilic homologue human UNG (hUNG). In order to understand the role of the enzyme-substrate interaction related to the cold-adapted properties, the structure of cUNG in complex with a bacteriophage encoded natural UNG inhibitor (Ugi) has been determined. The interaction has also been analyzed by isothermal titration calorimetry (ITC). The crystal structure of cUNG-Ugi was determined to a resolution of 1.9 Å with eight complexes in the asymmetric unit related through noncrystallographic symmetry. A comparison of the cUNG-Ugi complex with previously determined structures of UNG-Ugi shows that they are very similar, and confirmed the nucleotide-mimicking properties of Ugi. Biophysically, the interaction between cUNG and Ugi is very strong and shows a binding constant (Kb) which is one order of magnitude larger than that for hUNG-Ugi. The binding of both cUNG and hUNG to Ugi was shown to be favoured by both enthalpic and entropic forces; however, the binding of cUNG to Ugi is mainly dominated by enthalpy, while the entropic term is dominant for hUNG. The observed differences in the binding properties may be explained by an overall greater positive electrostatic surface potential in the protein-Ugi interface of cUNG and the slightly more hydrophobic surface of hUNG.

Psychrophilic enzymes: from folding to function and biotechnology

Psychrophiles thriving permanently at near-zero temperatures synthesize cold-active enzymes to sustain their cell cycle. Genome sequences, proteomic, and transcriptomic studies suggest various adaptive features to maintain adequate translation and proper protein folding under cold conditions. Most psychrophilic enzymes optimize a high activity at low temperature at the expense of substrate affinity, therefore reducing the free energy barrier of the transition state. Furthermore, a weak temperature dependence of activity ensures moderate reduction of the catalytic activity in the cold. In these naturally evolved enzymes, the optimization to low temperature activity is reached via destabilization of the structures bearing the active site or by destabilization of the whole molecule. This involves a reduction in the number and strength of all types of weak interactions or the disappearance of stability factors, resulting in improved dynamics of active site residues in the cold. These enzymes are already used in many biotechnological applications requiring high activity at mild temperatures or fast heat-inactivation rate. Several open questions in the field are also highlighted.

Uracil-DNA glycosylase: Structural, thermodynamic and kinetic aspects of lesion search and recognition

Uracil appears in DNA as a result of cytosine deamination and by incorporation from the dUTP pool. As potentially mutagenic and deleterious for cell regulation, uracil must be removed from DNA. The major pathway of its repair is initiated by uracil-DNA glycosylases (UNG), ubiquitously found enzymes that hydrolyze the N-glycosidic bond of deoxyuridine in DNA. This review describes the current understanding of the mechanism of uracil search and recognition by UNG. The structure of UNG proteins from several species has been solved, revealing a specific uracil-binding pocket located in a DNA-binding groove. DNA in the complex with UNG is highly distorted to allow the extrahelical recognition of uracil. Thermodynamic studies suggest that UNG binds with appreciable affinity to any DNA, mainly due to the interactions with the charged backbone. The increase in the affinity for damaged DNA is insufficient to account for the exquisite specificity of UNG for uracil. This specificity is likely to result from multistep lesion recognition process, in which normal bases are rejected at one or several pre-excision stages of enzyme-substrate complex isomerization, and only uracil can proceed to enter the active site in a catalytically competent conformation. Search for the lesion by UNG involves random sliding along DNA alternating with dissociation-association events and partial eversion of undamaged bases for initial sampling.

Structure of uracil-DNA N-glycosylase (UNG) from Vibrio cholerae: mapping temperature adaptation through structural and mutational analysis

The crystal structure of Vibrio cholerae uracil-DNA N-glycosylase (vcUNG) has been determined to 1.5 A resolution. Based on this structure, a homology model of Aliivibrio salmonicida uracil-DNA N-glycosylase (asUNG) was built. A previous study demonstrated that asUNG possesses typical cold-adapted features compared with vcUNG, such as a higher catalytic efficiency owing to increased substrate affinity. Specific amino-acid substitutions in asUNG were suggested to be responsible for the increased substrate affinity and the elevated catalytic efficiency by increasing the positive surface charge in the DNA-binding region. The temperature adaptation of these enzymes has been investigated using structural and mutational analyses, in which mutations of vcUNG demonstrated an increased substrate affinity that more resembled that of asUNG. Visualization of surface potentials revealed a more positive potential for asUNG compared with vcUNG; a modelled double mutant of vcUNG had a potential around the substrate-binding region that was more like that of asUNG, thus rationalizing the results obtained from the kinetic studies.

Convergent evolution of novel protein function in shrew and lizard venom

How do proteins evolve novel functions? To address this question, we are studying the evolution of a mammalian toxin, the serine protease BLTX [1], from the salivary glands of the North American shrew Blarina brevicauda. Here, we examine the molecular changes responsible for promoting BLTX toxicity. First, we show that regulatory loops surrounding the BLTX active site have evolved adaptively via acquisition of small insertions and subsequent accelerated sequence evolution. Second, these mutations introduce a novel chemical environment into the catalytic cleft of BLTX. Third, molecular-dynamic simulations show that the observed changes create a novel chemical and physical topology consistent with increased enzyme catalysis. Finally, we show that a toxic serine protease from the Mexican beaded lizard (GTX) [2] has evolved convergently through almost identical functional changes. Together, these results suggest that the evolution of toxicity might be predictable-arising via adaptive structural modification of analogous labile regulatory loops of an ancestral serine protease-and thus might aid in the identification of other toxic proteins.

Cold adaptation of zinc metalloproteases in the thermolysin family from deep sea and arctic sea ice bacteria revealed by catalytic and structural properties and molecular dynamics: new insights into relationship between conformational flexibility and hydrogen bonding

Increased conformational flexibility is the prevailing explanation for the high catalytic efficiency of cold-adapted enzymes at low temperatures. However, less is known about the structural determinants of flexibility. We reported two novel cold-adapted zinc metalloproteases in the thermolysin family, vibriolysin MCP-02 from a deep sea bacterium and vibriolysin E495 from an Arctic sea ice bacterium, and compared them with their mesophilic homolog, pseudolysin from a terrestrial bacterium. Their catalytic efficiencies, k(cat)/K(m) (10-40 degrees C), followed the order pseudolysin < MCP-02 < E495 with a ratio of approximately 1:2:4. MCP-02 and E495 have the same optimal temperature (T(opt), 57 degrees C, 5 degrees C lower than pseudolysin) and apparent melting temperature (T(m) = 64 degrees C, approximately 10 degrees C lower than pseudolysin). Structural analysis showed that the slightly lower stabilities resulted from a decrease in the number of salt bridges. Fluorescence quenching experiments and molecular dynamics simulations showed that the flexibilities of the proteins were pseudolysin < MCP-02 < E495, suggesting that optimization of flexibility is a strategy for cold adaptation. Molecular dynamics results showed that the ordinal increase in flexibility from pseudolysin to MCP-02 and E495, especially the increase from MCP-02 to E495, mainly resulted from the decrease of hydrogen-bond stability in the dynamic structure, which was due to the increase in asparagine, serine, and threonine residues. Finally, a model for the cold adaptation of MCP-02 and E495 was proposed. This is the first report of the optimization of hydrogen-bonding dynamics as a strategy for cold adaptation and provides new insights into the structural basis underlying conformational flexibility.

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k(cat) and reduced K(m) relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K(m)) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation.

Comparative unfolding studies of psychrophilic and mesophilic uracil DNA glycosylase: MD simulations show reduced thermal stability of the cold-adapted enzyme

Uracil DNA glycosylase (UDG) is a DNA repair enzyme involved in the base excision repair (BER) pathway, removing misincorporated uracil from the DNA strand. The native and mutant forms of Atlantic cod and human UDG have previously been characterized in terms of kinetic and thermodynamic properties as well as the determination of several crystal structures. This data shows that the cold-adapted enzyme is more catalytically efficient but at the same time less resistant to heat compared to its warm-active counterpart. In this study, the structure-function relationship is further explored by means of comparative molecular dynamics (MD) simulations at three different temperatures (375, 400 and 425K) to gain a deeper insight into the structural features responsible for the reduced thermostability of the cold-active enzyme. The simulations show that there are distinct structural differences in the unfolding pathway between the two homologues, particularly evident in the N- and C-terminals. Distortion of the mesophilic enzyme is initiated simultaneously in the N- and C-terminal, while the C-terminal part plays a key role for the stability of the psychrophilic enzyme. The simulations also show that at certain temperatures the cold-adapted enzyme unfolds faster than the warm-active homologues in accordance with the lower thermal stability found experimentally.

Cold-adapted enzymes

By far the largest proportion of the Earth's biosphere is comprised of organisms that thrive in cold environments (psychrophiles). Their ability to proliferate in the cold is predicated on a capacity to synthesize cold-adapted enzymes. These enzymes have evolved a range of structural features that confer a high level of flexibility compared to thermostable homologs. High flexibility, particularly around the active site, is translated into low-activation enthalpy, low-substrate affinity, and high specific activity at low temperatures. High flexibility is also accompanied by a trade-off in stability, resulting in heat lability and, in the few cases studied, cold lability. This review addresses the structure, function, and stability of cold-adapted enzymes, highlighting the challenges for immediate and future consideration. Because of the unique properties of cold-adapted enzymes, they are not only an important focus in extremophile biology, but also represent a valuable model for fundamental research into protein folding and catalysis.

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included alpha-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein.

Overexpression, purification, crystallization and preliminary X-ray analysis of uracil N-glycosylase from Mycobacterium tuberculosis in complex with a proteinaceous inhibitor

Uracil N-glycosylase is an enzyme which initiates the pathway of uracil-excision repair of DNA. The enzyme from Mycobacterium tuberculosis was co-expressed with a proteinaceous inhibitor from Bacillus subtilis phage and was crystallized in monoclinic space group C2, with unit-cell parameters a = 201.14, b = 64.27, c = 203.68 A, beta = 109.7 degrees. X-ray data from the crystal have been collected for structure analysis.

Conformation and stability of elastase from Atlantic cod, Gadus morhua

Metal binding and conformational stability characteristics of psychrophilic elastase (ACE) from Atlantic cod (Gadus morhua) has been investigated. Chelation to Ca(2+) was found to be important for maintaining the biologically active conformation and for the thermal stability of the enzyme. However, presence of metal ions such as Zn(2+), Fe(3+) and Cu(2+) was found to inhibit its hydrolytic activity and so did the chelating agent EDTA. Both pH and guanidinium chloride induced denaturation of the enzyme was followed by monitoring the changes in the tryptophan fluorescence. ACE exhibited a simple two-state unfolding pattern in both acidic and basic conditions with the midpoint of transition at pH values 4.08 and 10.29, respectively. Guanidinium chloride and heat induced denaturation of the enzyme was investigated at two pH values, 5.50 and 8.00, wherein the enzyme possesses similar tertiary structure but differ in its hydrolytic activity. Guanidinium chloride induced denaturation indicated that the enzyme unfolds with a C(m) of 1.53 M at pH 8.0 and a DeltaG(H2O) of 6.91 kJ mol(-1) (28.65 J mol(-1) residue(-1)) which is the lowest reported for psychrophilic enzymes investigated till-date. However, at pH 5.50, DeltaG(H2O) value is slightly lowered by 0.65 kJ mol(-1) consistent with the observed increase in the apparent quenching constant obtained with acrylamide. On the other hand, increase in T(m) by 38.45 degrees C was observed for the enzyme at acid pH (5.50) in comparison to the heat induced unfolding at pH 8.0. The increase in the apparent T(m) has been attributed to the possible weak intermolecular association of the enzyme molecules at moderately high temperatures that is favoured by the increase in the accessible surface area / dynamics under acidic conditions. The stability characteristics of ACE have been compared with the available data for mesophilic porcine pancreatic elastase and possible mechanism for the low temperature adaptation of ACE has been proposed.

Cold-active β-Galactosidase from Arthrobacter sp. C2-2 Forms Compact 660kDa Hexamers: Crystal Structure at 1.9Å Resolution

The X-ray structure of cold-active beta-galactosidase (isoenzyme C-2-2-1) from an Antarctic bacterium Arthrobacter sp. C2-2 was solved at 1.9A resolution. The enzyme forms 660 kDa hexamers with active sites opened to the central cavity of the hexamer and connected by eight channels with exterior solvent. To our best knowledge, this is the first cold-active beta-galactosidase with known structure and also the first known beta-galactosidase structure in the form of compact hexamers. The hexamer organization regulates access of substrates and ligands to six active sites and this unique packing, present also in solution, raises questions about its purpose and function. This enzyme belongs to glycosyl hydrolase family 2, similarly to Escherichia coli beta-galactosidase, forming tetramers necessary for its enzymatic function. However, we discovered significant differences between these two enzymes affecting the ability of tetramer/hexamer formation and complementation of the active site. This structure reveals new insights into the cold-adaptation mechanisms of enzymatic pathways of extremophiles.

Increased Flexibility as a Strategy for Cold Adaptation

Uracil DNA glycosylase (UDG) is a DNA repair enzyme in the base excision repair pathway and removes uracil from the DNA strand. Atlantic cod UDG (cUDG), which is a cold-adapted enzyme, has been found to be up to 10 times more catalytically active in the temperature range 15-37 degrees C as compared with the warm-active human counterpart. The increased catalytic activity of cold-adapted enzymes as compared with their mesophilic homologues are partly believed to be caused by an increase in the structural flexibility. However, no direct experimental evidence supports the proposal of increased flexibility of cold-adapted enzymes. We have used molecular dynamics simulations to gain insight into the structural flexibility of UDG. The results from these simulations show that an important loop involved in DNA recognition (the Leu(272) loop) is the most flexible part of the cUDG structure and that the human counterpart has much lower flexibility in the Leu(272) loop. The flexibility in this loop correlates well with the experimental k(cat)/K(m) values. Thus, the data presented here add strong support to the idea that flexibility plays a central role in adaptation to cold environments.