Thermal stabilities of DNA duplexes containing Gua (), α- () or β-anomer of formamidopyrimidine-N7-9-hydroxy-aflatoxin B₁ () differ markedly (T: a < g < b ), but the underlying molecular origin of this experimentally observed phenomenon is yet to be identified and determined. Here, by employing explicit-solvent molecular dynamics simulations coupled with free-energy calculations using a combined linear-interaction-energy/linear-response-approximation approach, we explain the quantitative differences in T m in terms of three structural features (bulkiness, order, and compactness) and three energetical contributions (non-polar, electrostatic, and preorganized-electrostatic), and thus advance the current understanding of the relationships between structures, free energies, and thermal stabilities of DNA double helices.
View Article and Find Full Text PDFHuman X-family DNA polymerases β (Polβ) and λ (Polλ) catalyze the nucleotidyl-transfer reaction in the base excision repair pathway of the cellular DNA damage response. Using empirical valence bond and free-energy perturbation simulations, we explore the feasibility of various mechanisms for the deprotonation of the 3'-OH group of the primer DNA strand, and the subsequent formation and cleavage of P-O bonds in four Polβ, two truncated Polλ (tPolλ), and two tPolλ Loop1 mutant (tPolλΔL1) systems differing in the initial X-ray crystal structure and nascent base pair. The average calculated activation free energies of 14, 18, and 22 kcal mol for Polβ, tPolλ, and tPolλΔL1, respectively, reproduce the trend in the observed catalytic rate constants.
View Article and Find Full Text PDFTunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures.
View Article and Find Full Text PDFWe carried out free-energy calculations and transient kinetic experiments for the insertion of the right (dC) and wrong (dA) nucleotides by wild-type (WT) and six mutant variants of human DNA polymerase β (Pol β). Since the mutated residues in the point mutants, I174S, I260Q, M282L, H285D, E288K, and K289M, were not located in the Pol β catalytic site, we assumed that the WT and its point mutants share the same dianionic phosphorane transition-state structure of the triphosphate moiety of deoxyribonucleotide 5'-triphosphate (dNTP) substrate. On the basis of this assumption, we have formulated a thermodynamic cycle for calculating relative dNTP insertion efficiencies, Ω = (k(pol)/K(D))(mut)/(k(pol)/K(D))(WT) using free-energy perturbation (FEP) and linear interaction energy (LIE) methods.
View Article and Find Full Text PDFWe show how a restricted reaction surface can be used to facilitate the calculation of biologically important contributions of active site geometries and dynamics to DNA polymerase fidelity. Our analysis, using human DNA polymerase beta (pol β), is performed within the framework of an electrostatic linear free energy response (EFER) model. The structure, dynamics, and energetics of pol β-DNA-dNTP interactions are computed between two points on the multidimensional reaction free energy surface.
View Article and Find Full Text PDF