Investigation of Nascent Base Pair and Polymerase Behavior in the Presence of Mismatches in DNA Polymerase I Using Molecular Dynamics

Incorporation of N7-Platinated Guanines into Thermus Aquaticus (Taq) DNA Polymerase: Atomistic Insights from Molecular Dynamics Simulations

In this work, we elucidated some key aspects of the mechanism of action of the cisplatin anticancer drug, cis-[Pt(NH₃)₂Cl₂], involving direct interactions with free nucleotides. A comprehensive in silico molecular modeling analysis was conducted to compare the interactions of Thermus aquaticus (Taq) DNA polymerase with three distinct N7-platinated deoxyguanosine triphosphates: [Pt(dien)(N7-dGTP)] (1), cis-[Pt(NH₃)₂Cl(N7-dGTP)] (2), and cis-[Pt(NH₃)₂(H₂O)(N7-dGTP)] (3) {dien = diethylenetriamine; dGTP = 5'-(2'-deoxy)-guanosine-triphosphate}, using canonical dGTP as a reference, in the presence of DNA. The goal was to elucidate the binding site interactions between Taq DNA polymerase and the tested nucleotide derivatives, providing valuable atomistic insights. Unbiased molecular dynamics simulations (200 ns for each complex) with explicit water molecules were performed on the four ternary complexes, yielding significant findings that contribute to a better understanding of experimental results. The molecular modeling highlighted the crucial role of a specific α-helix (O-helix) within the fingers subdomain, which facilitates the proper geometry for functional contacts between the incoming nucleotide and the DNA template needed for incorporation into the polymerase. The analysis revealed that complex 1 exhibits a much lower affinity for Taq DNA polymerase than complexes 2-3. The affinities of cisplatin metabolites 2-3 for Taq DNA polymerase were found to be quite similar to those of natural dGTP, resulting in a lower incorporation rate for complex 1 compared to complexes 2-3. These findings could have significant implications for the cisplatin mechanism of action, as the high intracellular availability of free nucleobases might promote the competitive incorporation of platinated nucleotides over direct cisplatin attachment to DNA. The study's insights into the incorporation of platinated nucleotides into the Taq DNA polymerase active site suggest that the role of platinated nucleotides in the cisplatin mechanism of action may have been previously underestimated.

Single-molecule Taq DNA polymerase dynamics

Taq DNA polymerase functions at elevated temperatures with fast conformational dynamics-regimes previously inaccessible to mechanistic, single-molecule studies. Here, single-walled carbon nanotube transistors recorded the motions of Taq molecules processing matched or mismatched template-deoxynucleotide triphosphate pairs from 22° to 85°C. By using four enzyme orientations, the whole-enzyme closures of nucleotide incorporations were distinguished from more rapid, 20-μs closures of Taq's fingers domain testing complementarity and orientation. On average, one transient closure was observed for every nucleotide binding event; even complementary substrate pairs averaged five transient closures between each catalytic incorporation at 72°C. The rate and duration of the transient closures and the catalytic events had almost no temperature dependence, leaving all of Taq's temperature sensitivity to its rate-determining open state.

Molecular Dynamics Simulations Reveal Differentiated Context-Dependent Conformational Dynamics of Two Proteins of the Same Family

The Arabidopsis pyrabactin resistant 1 (PYR1)-like family of proteins (PYLs) are receptors of abscisic acid (ABA), an essential small signaling molecule in plants. Here, we report a comparative molecular dynamics (MD) study on two PYL members, PYR1 and PYL10, which, despite their highly similar sequences and structures, have been suggested to belong to two different subclasses of PYLs, one being dimeric and relying on binding to ABA to inhibit downstream type 2C protein phosphatases (PP2Cs) and the other being monomeric and able to constitutively inhibit downstream PP2Cs without ABA. MD simulations have been carried out on these proteins in various monomeric or complexation states. Analyses of the simulations unambiguously confirm that ABA has large effects on the conformational dynamics of PYR1 but not PYL10, whereas a downstream PP2C has much larger effects on PYL10 than on PYR1. The differentiated effects are consistent with the functional differences between the two proteins. Potential of mean forces (PMFs) calculated by umbrella sampling showed that binding to ABA strengthens the PYR1-PP2C complex, increasing the PMF change for dissociation from 7.5 to 12.0 kcal mol^-1. On the other hand, the same PMF change for an apo-PYL10-PP2C complex was computed to be 9.5 kcal mol^-1, suggesting stronger binding in apo-PYL10-PP2C than in apo-PYR1-PP2C. Several specific sequence features that may contribute to the functional differentiation between PYR1 and PYL10 are suggested based on the intersubunit residue-residue contacts occurred in the simulations.