Brian N. Dominy

Analyzing the robustness of the MM/PBSA free energy calculation method: Application to DNA conformational transitions

The ability to predict and characterize free energy differences associated with conformational equilibria or the binding of biomolecules is vital to understanding the molecular basis of many important biological functions. As biological studies focus on larger molecular complexes and properties of the genome, proteome, and interactome, the development and characterization of efficient methods for calculating free energy becomes increasingly essential. The aim of this study is to examine the robustness of the end‐point free energy method termed the molecular mechanics Poisson‐Boltzmann solvent accessible surface area (MM/PBSA) method. Specifically, applications of MM/PBSA to the conformational equilibria of nucleic acid (NA) systems are explored. This is achieved by comparing A to B form DNA conformational free energy differences calculated using MM/PBSA with corresponding free energy differences determined with a more rigorous and time‐consuming umbrella sampling algorithm. In addition, the robustness of NA MM/PBSA calculations is also evaluated in terms of the sensitivity towards the choice of force field and the choice of solvent model used during conformational sampling. MM/PBSA calculations of the free energy difference between A‐form and B‐form DNA are shown to be in very close agreement with the PMF result determined using an umbrella sampling approach. Further, it is found that the MM/PBSA conformational free energy differences were also in agreement using either the CHARMM or AMBER force field. The influence of ionic strength on conformational stability was particularly insensitive to the choice of force field. Finally, it is also shown that the use of a generalized Born implicit solvent during conformational sampling results in free energy estimates that deviate slightly from those obtained using explicitly solvated MD simulations in these NA systems.

New Family of Deamination Repair Enzymes in Uracil-DNA Glycosylase Superfamily

DNA glycosylases play a major role in the repair of deaminated DNA damage. Previous investigations identified five families within the uracil-DNA glycosylase (UDG) superfamily. All enzymes within the superfamily studied thus far exhibit uracil-DNA glycosylase activity. Here we identify a new class of DNA glycosylases in the UDG superfamily that lacks UDG activity. Instead, these enzymes act as hypoxanthine-DNA glycosylases in vitro and in vivo. Molecular modeling and structure-guided mutational analysis allowed us to identify a unique catalytic center in this class of DNA glycosylases. Based on unprecedented biochemical properties and phylogenetic analysis, we propose this new class of DNA repair glycosylases that exists in bacteria, archaea, and eukaryotes as family 6 and designate it as the hypoxanthine-DNA glycosylase family. This study demonstrates the structural evolvability that underlies substrate specificity and catalytic flexibility in the evolution of enzymatic function.

Comparison of Solvation-Effect Methods for the Simulation of Peptide Interactions with a Hydrophobic Surface

In this study we investigated the interaction behavior between thirteen different small peptides and a hydrophobic surface using three progressively more complex methods of representing solvation effects: a united‐atom implicit solvation method [CHARMM 19 force field (C19) with Analytical Continuum Electrostatics (ACE)], an all‐atom implicit solvation method (C22 with GBMV), and an all‐atom explicit solvation method (C22 with TIP3P). The adsorption behavior of each peptide was characterized by the calculation of the potential of mean force as a function of peptide‐surface separation distance. The results from the C22/TIP3P model suggest that hydrophobic peptides exhibit relatively strong adsorption behavior, polar and positively‐charged peptides exhibit negligible to relatively weak favorable interactions with the surface, and negatively‐charged peptides strongly resist adsorption. Compared to the TIP3P model, the ACE and GBMV implicit solvent models predict much stronger attractions for the hydrophobic peptides as well as stronger repulsions for the negatively‐charged peptides on the CH3‐SAM surface. These comparisons provide a basis from which each of these implicit solvation methods may be reparameterized to provide closer agreement with explicitly represented solvation in simulations of peptide and protein adsorption to functionalized surfaces.

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.

Identification of Escherichia coli mismatch-specific uracil DNA glycosylase as a robust xanthine DNA glycosylase.

The gene for the mismatch-specific uracil DNA glycosylase (MUG) was identified in the Escherichia coli genome as a sequence homolog of the human thymine DNA glycosylase with activity against mismatched uracil base pairs. Examination of cell extracts led us to detect a previously unknown xanthine DNA glycosylase (XDG) activity in E. coli. DNA glycosylase assays with purified enzymes indicated the novel XDG activity is attributable to MUG. Here, we report a biochemical characterization of xanthine DNA glycosylase activity in MUG. The wild type MUG possesses more robust activity against xanthine than uracil and is active against all xanthine-containing DNA (C/X, T/X, G/X, A/X and single-stranded X). Analysis of potentials of mean force indicates that the double-stranded xanthine base pairs have a relatively narrow energetic difference in base flipping, whereas the tendency for uracil base flipping follows the order of C/U > G/U > T/U > A/U. Site-directed mutagenesis performed on conserved motifs revealed that Asn-140 and Ser-23 are important determinants for XDG activity in E. coli MUG. Molecular modeling and molecular dynamics simulations reveal distinct hydrogen-bonding patterns in the active site of E. coli MUG that account for the specificity differences between E. coli MUG and human thymine DNA glycosylase as well as that between the wild type MUG and the Asn-140 and Ser-23 mutants. This study underscores the role of the favorable binding interactions in modulating the specificity of DNA glycosylases.

Specificity and Catalytic Mechanism in Family 5 Uracil DNA Glycosylase

Background: Uracil DNA glycosylases are DNA repair enzymes involved in the removal of base damage. Results: Family 5 UDGb is a uracil, hypoxanthine, and xanthine DNA glycosylase. Conclusion: Family 5 UDGb adapts multiple catalytic amino acids for the excision of pyrimidine and purine deaminated DNA bases. Significance: Family 5 UDGb exemplifies functional diversity in enzyme superfamilies. UDGb belongs to family 5 of the uracil DNA glycosylase (UDG) superfamily. Here, we report that family 5 UDGb from Thermus thermophilus HB8 is not only a uracil DNA glycosyase acting on G/U, T/U, C/U, and A/U base pairs, but also a hypoxanthine DNA glycosylase acting on G/I, T/I, and A/I base pairs and a xanthine DNA glycosylase acting on all double-stranded and single-stranded xanthine-containing DNA. Analysis of potentials of mean force indicates that the tendency of hypoxanthine base flipping follows the order of G/I > T/I, A/I > C/I, matching the trend of hypoxanthine DNA glycosylase activity observed in vitro. Genetic analysis indicates that family 5 UDGb can also act as an enzyme to remove uracil incorporated into DNA through the existence of dUTP in the nucleotide pool. Mutational analysis coupled with molecular modeling and molecular dynamics analysis reveals that although hydrogen bonding to O2 of uracil underlies the UDG activity in a dissociative fashion, Tth UDGb relies on multiple catalytic residues to facilitate its excision of hypoxanthine and xanthine. This study underscores the structural and functional diversity in the UDG superfamily.

The evolution of cefotaximase activity in the TEM β-lactamase.

The development of a molecular-level understanding of drug resistance through β-lactamase is critical not only in designing newer-generation antibacterial agents but also in providing insight into the evolutionary mechanisms of enzymes in general. In the present study, we have evaluated the effect of four drug resistance mutations (A42G, E104K, G238S, and M182T) on the cefotaximase activity of the TEM-1 β-lactamase. Using computational methods, including docking and molecular mechanics calculations, we have been able to correctly identify the relative order of catalytic activities associated with these four single point mutants. Further analyses suggest that the changes in catalytic efficiency for mutant enzymes are correlated to structural changes within the binding site. Based on the energetic and structural analyses of the wild-type and mutant enzymes, structural rearrangement is suggested as a mechanism of evolution of drug resistance through TEM β-lactamase. The present study not only provides molecular-level insight into the effect of four drug resistance mutations on the structure and function of the TEM β-lactamase but also establishes a foundation for a future molecular-level analysis of complete evolutionary trajectory for this class of enzymes.

A structural determinant in the uracil DNA glycosylase superfamily for the removal of uracil from adenine/uracil base pairs

The uracil DNA glycosylase superfamily consists of several distinct families. Family 2 mismatch-specific uracil DNA glycosylase (MUG) from Escherichia coli is known to exhibit glycosylase activity on three mismatched base pairs, T/U, G/U and C/U. Family 1 uracil N-glycosylase (UNG) from E. coli is an extremely efficient enzyme that can remove uracil from any uracil-containing base pairs including the A/U base pair. Here, we report the identification of an important structural determinant that underlies the functional difference between MUG and UNG. Substitution of a Lys residue at position 68 with Asn in MUG not only accelerates the removal of uracil from mismatched base pairs but also enables the enzyme to gain catalytic activity on A/U base pairs. Binding and kinetic analysis demonstrate that the MUG-K68N substitution results in enhanced ground state binding and transition state interactions. Molecular modeling reveals that MUG-K68N, UNG-N123 and family 5 Thermus thermophiles UDGb-A111N can form bidentate hydrogen bonds with the N3 and O4 moieties of the uracil base. Genetic analysis indicates the gain of function for A/U base pairs allows the MUG-K68N mutant to remove uracil incorporated into the genome during DNA replication. The implications of this study in the origin of life are discussed.

Thermodynamic resolution: How do errors in modeled protein structures affect binding affinity predictions?

The present study addresses the effect of structural distortion, caused by protein modeling errors, on calculated binding affinities toward small molecules. The binding affinities to a total of 300 distorted structures based on five different protein–ligand complexes were evaluated to establish a broadly applicable relationship between errors in protein structure and errors in calculated binding affinities. Relatively accurate protein models (less than 2 Å RMSD within the binding site) demonstrate a 14.78 (±7.5)% deviation in binding affinity from that calculated by using the corresponding crystal structure. For structures of 2–3 Å, 3–4 Å, and >4 Å RMSD within the binding site, the error in calculated binding affinity increases to 20.8 (±5.98), 22.79 (±11.3), and 29.43 (±11.47)%, respectively. The results described here may be used in combination with other tools to evaluate the utility of modeled protein structures for drug development or other ligand‐binding studies. Proteins 2010.

Correlated Mutation in the Evolution of Catalysis in Uracil DNA Glycosylase Superfamily

Enzymes in Uracil DNA glycosylase (UDG) superfamily are essential for the removal of uracil. Family 4 UDGa is a robust uracil DNA glycosylase that only acts on double-stranded and single-stranded uracil-containing DNA. Based on mutational, kinetic and modeling analyses, a catalytic mechanism involving leaving group stabilization by H155 in motif 2 and water coordination by N89 in motif 3 is proposed. Mutual Information analysis identifies a complexed correlated mutation network including a strong correlation in the EG doublet in motif 1 of family 4 UDGa and in the QD doublet in motif 1 of family 1 UNG. Conversion of EG doublet in family 4 Thermus thermophilus UDGa to QD doublet increases the catalytic efficiency by over one hundred-fold and seventeen-fold over the E41Q and G42D single mutation, respectively, rectifying the strong correlation in the doublet. Molecular dynamics simulations suggest that the correlated mutations in the doublet in motif 1 position the catalytic H155 in motif 2 to stabilize the leaving uracilate anion. The integrated approach has important implications in studying enzyme evolution and protein structure and function.