Brian N. Dominy

One of the anthrax receptors, tumor endothelial marker 8 (TEM8), is reported to be a potential anticancer target due to its over‐expression during tumor angiogenesis. To extend our BioLayer Interferometry study in PA‐TEM8 binding, we present a computational approach to reveal the role of an integral metal ion on receptor structure and binding thermodynamics. We estimated the interaction energy between PA and TEM8 using computer simulation. Consistent with experimental study, computational results indicate the metal ion in TEM8 contributes significantly to the binding affinity, and PA‐TEM8 binding is more favorable in the presence of Mg²⁺ than Ca²⁺. Further, computational analysis suggests that the differences in PA‐TEM8 binding affinity are comparable to the closely related integrin proteins. The conformation change, which linked to changes in activity of integrins, was not found in TEM8. In the present of Mg²⁺, TEM8 remains in a conformation analogous to an integrin open (high‐affinity) conformation.

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.

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.

Examining Electrostatic Influences on Base-Flipping: A Comparison of TIP3P and GB Solvent Models - Volume 13 Issue 1 - Allyn R. Brice, Brian N. Dominy

Modeling the Influence of Salt on the Hydrophobic Effect and Protein Fold Stability - Volume 13 Issue 1 - Mihir S. Date, Brian N. Dominy

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. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011