Increased flexibility as a strategy for cold adaptation: a comparative molecular dynamics study of cold- and warm-active uracil DNA glycosylase

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.     

Fluctuations between stabilizing and destabilizing electrostatic contributions of ion pairs in conformers of the c-Myc-Max leucine zipper

In solution proteins often exhibit backbone and side-chain flexibility. Yet electrostatic interactions in proteins are sensitive to motions. Hence, here we study the contribution of ion pairs toward protein stability in a range of conformers which sample the conformational space in solution. Specifically, we focus on the electrostatic contributions of ion pairs to the stability of each of the conformers in the NMR ensemble of the c-Myc-Max leucine zipper and to their average energy minimized structure. We compute the electrostatic contributions of inter- and intra-helical ion pairs and of an ion pair network. We find that the electrostatic contributions vary considerably among the 40 NMR conformers. Each ion pair, and the network, fluctuates between being stabilizing and being destabilizing. This fluctation reflects the variability in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other and with respect to other charged groups in the rest of the protein. Ion pair interactions in the c-Myc-Max leucine zipper in solution depend on the protein conformer which is analyzed. Hence, the overall stabilizing (or destabilizing) contribution of an ion pair is conformer population-dependent. This study indicates that free energy calculations performed using the continuum electrostatics methodology are sensitive to protein conformational details.     

Dynamic arrangement of ion pairs and individual contributions to the thermal stability of the cofactor-binding domain of glutamate dehydrogenase from Thermotoga maritima

The dynamics of a hyperthermophilic protein fragment in a water environment, as studied by performing molecular dynamics (MD) simulations at various temperatures, is compared to the dynamical behavior of a homologous mesophilic protein simulated under identical conditions. The effects on the stability of the spatial arrangement and mobility of the charged residues in solution were quantified by calculating free energy changes upon salt bridge formation in these proteins. Electrostatic free energy terms derived from a thermodynamic cycle were obtained by solving the linearized Poisson-Boltzmann equation for a series of protein conformations generated by MD simulations and placed subsequently in a continuum solvent medium. Our results show that the ion pairs are electrostatically stabilizing in most of the cases, but their individual contributions vary significantly. The greater contribution of the charged residues to the stability of the hyperthermophilic protein as compared with the mesophilic counterpart was evidenced only by the calculations that included conformations sampled at 343 and 373 K. The "dynamic" structure of the hyperthermophilic protein fragment simulated at elevated temperatures reveals an optimum placement of the ionizable residues within the protein structure as well as the role of their cooperative interactions in promoting thermal stability. The thermodynamic properties such as electrostatic free energy differences, configurational entropies, and specific heat capacities calculated in the dynamic context of the protein structure provided new insight into the mechanism of protein thermostabilization.     

Crystal structure of a family 4 uracil-DNA glycosylase from Thermus thermophilus HB8

Uracil-DNA glycosylase (UDG; EC 3.2.2.-) removes uracil from DNA to initiate DNA base excision repair. Since hydrolytic deamination of cytosine to uracil is one of the most frequent DNA-damaging events in all cells, UDG is an essential enzyme for maintaining the integrity of genomic information. For the first time, we report the crystal structure of a family 4 UDG from Thermus thermophilus HB8 (TthUDG) complexed with uracil, solved at 1.5 angstroms resolution. As opposed to UDG enzymes in its other families, TthUDG possesses a [4Fe-4S] cluster. This iron-sulfur cluster, which is distant from the active site, interacts with loop structures and has been suggested to be unessential to the activity but necessary for stabilizing the loop structures. In addition to the iron-sulfur cluster, salt-bridges and ion pairs on the molecular surface and the presence of proline on loops and turns is thought to contribute to the enzyme's thermostability. Despite very low levels of sequence identity with Escherichia coli and human UDGs (family 1) and E.coli G:T/U mismatch-specific DNA glycosylase (MUG) (family 2), the topology and order of secondary structures of TthUDG are similar to those of these distant relatives. Furthermore, the coordinates of the core structure formed by beta-strands are almost the same. Positive charge is distributed over the active-site groove, where TthUDG would bind DNA strands, as do UDG enzymes in other families. TthUDG recognizes uracil specifically in the same manner as does human UDG (family 1), rather than guanine in the complementary strand DNA, as does E.coli MUG (family 2). These results suggest that the mechanism by which family 4 UDGs remove uracils from DNA is similar to that of family 1 enzymes.     

Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation

The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg²⁺) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg²⁺ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg²⁺ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg²⁺-excess concentrations and comparison with fluorescence titration measurements shows that on average 20–30% of the Mg²⁺ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.     

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. Figure Electrostatic isosurfaces of cod uracil DNA glycosylase in complex with double stranded DNA     

A pragmatic approach to structure based calculation of coupled proton and electron transfer in proteins

The coupled motion of electrons and protons occurs in many proteins. Using appropriate tools for calculation, the three-dimensional protein structure can show how each protein modulates the observed electron and proton transfer reactions. Some of the assumptions and limitations involved in calculations that rely on continuum electrostatics to calculate the energy of charges in proteins are outlined. Approaches that mix molecular mechanics and continuum electrostatics are described. Three examples of the analysis of reactions in photosynthetic reaction centers are given: comparison of the electrochemistry of hemes in different sites; analysis of the role of the protein in stabilizing the early charge separated state in photosynthesis; and calculation of the proton uptake and protein motion coupled to the electron transfer from the primary (Q(A)) to secondary (Q(B)) quinone. Different mechanisms for stabilizing intra-protein charged cofactors are highlighted in each reaction.     

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).     

Heat-shock proteins associated with base excision repair enzymes in HeLa cells

Two enzymes of base excision repair (BER), uracil DNA glycosylase (UDG) and DNA polymerase beta (beta pol), from HeLa cells co-eluted from Superose 12 FPLC columns. The UDG was completely displaced from 150-180-kDa fractions to 30- 70-kDa fractions by brief treatment with 0.5 N NaCl, pH 3.0, as expected when protein-protein associations are disrupted, but beta pol was not displaced by this treatment. UDG was not essential to the presence of beta pol in the 150-180-kDa enzyme complex. beta pol and UDG apparently reside in separate but co-eluting structures. Immunoaffinity chromatography showed that the association of UDG and beta pol was accounted for by attachment in common to DNA and that the association was abolished by eliminating DNA. Evidence for base excision repairosomes containing UDG and beta pol in protein-protein assemblies was not found. However, UDG and human AP endonuclease (HAP1) were associated with HSP70 and HSP27, which are present in 150-180-kDa and 30-70-kDa proteins of cell sonicates. The association of HSPs with BER enzymes was confirmed by hydroxyl radical protein-protein footprinting and immunoaffinity tests. The association of HSPs and BER enzymes is a novel finding. HSP binding may account for the presence of BER enzymes in the two large size class fractions and HSPs may have functional roles in BER.     

Insights from xanthine and uracil DNA glycosylase activities of bacterial and human SMUG1: switching SMUG1 to UDG

Single-strand-selective monofunctional uracil DNA glycosylase (SMUG1) belongs to Family 3 of the uracil DNA glycosylase (UDG) superfamily. Here, we report that a bacterial SMUG1 ortholog in Geobacter metallireducens (Gme) and the human SMUG1 enzyme are not only UDGs but also xanthine DNA glycosylases (XDGs). In addition, mutational analysis and molecular dynamics (MD) simulations of Gme SMUG1 identify important structural determinants in conserved motifs 1 and 2 for XDG and UDG activities. Mutations at M57 (M57L) and H210 (H210G, H210M, and H210N), both of which are involved in interactions with the C2 carbonyl oxygen in uracil or xanthine, cause substantial reductions in XDG and UDG activities. Increased selectivity is achieved in the A214R mutant of Gme SMUG1, which corresponds to a position involved in base flipping. This mutation results in an activity profile resembling a human SMUG1-like enzyme as exemplified by the retention of UDG activity on mismatched base pairs and weak XDG activity. MD simulations indicate that M57L increases the flexibility of the motif 2 loop region and specifically A214, which may account for the reduced catalytic activity. G60Y completely abolishes XDG and UDG activity, which is consistent with a modeled structure in which G60Y blocks the entry of either xanthine or uracil to the base binding pocket. Most interestingly, a proline substitution at the G63 position switches the Gme SMUG1 enzyme to an exclusive UDG as demonstrated by the uniform excision of uracil in both double-stranded and single-stranded DNA and the complete loss of XDG activity. MD simulations indicate that a combination of a reduced free volume and altered flexibility in the active-site loops may underlie the dramatic effects of the G63P mutation on the activity profile of SMUG1. This study offers insights on the important role that modulation of conformational flexibility may play in defining specificity and catalytic efficiency.     

The role of the dielectric barrier in narrow biological channels: a novel composite approach to modeling single-channel currents

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.     

Identification of a prototypical single-stranded uracil DNA glycosylase from Listeria innocua

A recent phylogenetic study on UDG superfamily estimated a new clade of family 3 enzymes (SMUG1-like), which shares a lower homology with canonic SMUG1 enzymes. The enzymatic properties of the newly found putative DNA glycosylase are unknown. To test the potential UDG activity and evaluate phylogenetic classification, we isolated one SMUG1-like glycosylase representative from Listeria innocua (Lin). A biochemical screening of DNA glycosylase activity in vitro indicates that Lin SMUG1-like glycosylase is a single-strand selective uracil DNA glycosylase. The UDG activity on DNA bubble structures provides clue to its physiological significance in vivo. Mutagenesis and molecular modeling analyses reveal that Lin SMUG1-like glycosylase has similar functional motifs with SMUG1 enzymes; however, it contains a distinct catalytic doublet S67-S68 in motif 1 that is not found in any families in the UDG superfamily. Experimental investigation shows that the S67M-S68N double mutant is catalytically more active than either S67M or S68N single mutant. Coupled with mutual information analysis, the results indicate a high degree of correlation in the evolution of SMUG1-like enzymes. This study underscores the functional and catalytic diversity in the evolution of enzymes in UDG superfamily.     

Proteolytic degradation of the nuclear isoform of uracil-DNA glycosylase occurs during the S phase of the cell cycle

Uracil-DNA glycosylases are enzymes that remove uracil from DNA and initiate base-excision repair. These enzymes play a key role in maintaining genomic integrity by reducing the mutagenic events caused by G:C to A:T transition mutations. The recent finding that a family of RNA editing enzymes (APOBECs) can deaminate cytosine in DNA has raised the interest in these base-excision repair enzymes. This research focuses on the regulation of the nuclear isoform of uracil-DNA glycosylase, a 36000 Da protein that contains a unique 44 amino acid N-terminus. In synchronized HeLa cells, UDG1A protein levels decrease to barely detectable levels during the S phase of the cell cycle. Immunoblot analysis of immunoprecipitated or affinity-isolated UDG1A reveals ubiquitin-conjugated UDG1A when proteolysis is inhibited using N-acetyl-leu-leu-norleu-al or MG132, inhibitors of proteosomal dependent protein degradation. Transient transfection experiments, with histidine-tagged ubiquitin, were used to confirm that endogenous UDG1A is ubiquitinated in vivo. Addition of the nuclear export inhibitor, leptomycin B, prevents ubiquitination and degradation of UDG1A. This indicates that translocation from the nucleus may be a step in UDG1A turnover. Finally, UDG1A protein degradation is prevented when cells are incubated with the cyclin-dependent kinase inhibitor, roscovitine. These results suggest that protein phosphorylation and/or nuclear export participate in the post-translational regulation of UDG1A protein levels.     

On the role of electrostatic interactions in the design of protein-protein interfaces

Here, the methods of continuum electrostatics are used to investigate the contribution of electrostatic interactions to the binding of four protein-protein complexes; barnase-barstar, human growth hormone and its receptor, subtype N9 influenza virus neuraminidase and the NC41 antibody, the Ras binding domain (RBD) of kinase cRaf and a Ras homologue Rap1A. In two of the four complexes electrostatics are found to strongly oppose binding (hormone-receptor and neuraminidase-antibody complexes), in one case the net effect is close to zero (barnase-barstar) and in one case electrostatics provides a significant driving force favoring binding (RBD-Rap1A). In order to help understand the wide range of electrostatic contributions that were calculated, the electrostatic free energy was partitioned into contributions of individual charged and polar residues, salt bridges and networks involving salt bridges and hydrogen bonds. Although there is no one structural feature that accounts for the differences between the four interfaces, the extent to which the desolvation of buried charges is compensated by the formation of hydrogen bonds and ion pairs appears to be an important factor. Structural features that are correlated with contribution of an individual residue to stability are also discussed. These include partial burial of a charged group in the free monomer, the formation of networks involving charged and polar amino acids, and the formation of partially exposed ion-pairs. The total electrostatic contribution to binding is found to be inversely correlated with buried total and non-polar surface area. This suggests that different interfaces can be designed to exploit electrostatic and hydrophobic forces in very different ways.     

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.     

Energetics of the cleft closing transition and the role of electrostatic interactions in conformational rearrangements of the glutamate receptor ligand binding domain

The ionotropic glutamate receptors are localized in the pre- and postsynaptic membrane of neurons in the brain. Activation by the principal excitatory neurotransmitter glutamate allows the ligand binding domain to change conformation, communicating opening of the channel for ion conduction. The free energy of the GluR2 S1S2 ligand binding domain (S1S2) closure transition was computed using a combination of thermodynamic integration and umbrella sampling modeling methods. A path that involves lowering the charge on E705 was chosen to clarify the role of this binding site residue. A continuum electrostatics approach in S1S2 is used to show E705, located in the ligand binding cleft, stabilizes the closed conformation of S1S2 via direct interactions with other protein residues, not through the ligand. In the closed conformation, in the absence of a ligand, S1S2 is somewhat more closed than what has been reported in X-ray structures. A semiopen conformation has been identified which is characterized by disruption of a single cross-cleft interaction and differs only slightly in energy from the fully closed S1S2. The fully open S1S2 conformation exhibits a wide energy well and shares structural similarity with the apo S1S2 crystal structure. Hybrid continuum electrostatics/MD calculations along the chosen closure transition pathway reveal solvation energies, and electrostatic interaction energies between two lobes of the protein increase the relative energetic difference between the open and closed conformational states. By analyzing the role of several cross-cleft contacts as well as other binding site residues, we demonstrate how S1S2 interactions facilitate formation of the closed conformation of the GluR2 ligand binding domain.     

Mutational effects on O(6)-methylguanine-DNA methyltransferase from hyperthermophile: contribution of ion-pair network to protein thermostability

Ion pairs have been considered to be general stabilizing factors in hyperthermophilic proteins, but the present experimental data cannot fully explain how ion pairs and ion-pair networks contribute to the stability. In this paper, we show experimental evidence that not all of the internal ion pairs contribute to the thermal and thermodynamic stability, using O(6)-methylguanine-DNA methyltransferase from Thermococcus kodakaraensis KOD1 (Tk-MGMT) as a model protein. Of three mutants in which an inter-helical ion pair was disrupted, only one mutant (E93A) was shown to be destabilized. Delta G of E93A was lower by approximately 4 kJ mol(-1) than that of the wild type, and E93A unfolded one order of magnitude faster than did the wild type and other variants. Glu 93 has unique properties in forming an ion-pair network that bridges the N- and C-terminal domains and connects three helices in the protein interior.     

Linear free energy correlations for enzymatic base flipping: how do damaged base pairs facilitate specific recognition?

To efficiently maintain their genomic integrity, DNA repair glycosylases must exhibit high catalytic specificity for their cognate damaged bases using an extrahelical recognition mechanism. One possible contribution to specificity is the weak base pairing and inherent instability of damaged sites which may lead to increased extrahelicity of the damaged base and enhanced recognition of these sites. This model predicts that the binding affinity of the enzyme should increase as the thermodynamic stability of the lesion base pair decreases, because less work is required to extrude the base into its active site. We have tested this hypothesis with uracil DNA glycosylase (UDG) by constructing a series of DNA duplexes containing a single uracil (U) opposite a variety of bases (X) that formed from zero to three hydrogen bonds with U. Linear free energy (LFE) relationships were observed that correlated UDG binding affinity with the entropy and enthalpy of duplex melting, and the dynamic accessibility of the damaged site to chemical oxidation. These LFEs indicate that the increased conformational freedom of the damaged site brought about by enthalpic destabilization of the base pair promotes the formation of extrahelical states that enhance specific recognition by as much as 3000-fold. However, given the small stability differences between normal base pairs and U.A or U.G base pairs, relative base pair stability contributes little to the >10(6)-fold discrimination of UDG for uracil sites in cellular DNA. In contrast, the intrinsic instability of other more egregious DNA lesions may contribute significantly to the specificity of other DNA repair enzymes that bind to extrahelical bases.