We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER–HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single–point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo-bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization.

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

Proton transfer reactions are ubiquitous in chemistry, especially in aqueous solutions. We investigate photoinduced proton transfer between the photoacid 8-hydroxypyrene-1,3,6trisulfonate (HPTS) and water using fast fluorescence spectroscopy and ab initio molecular dynamics simulations. Photoexcitation causes rapid proton release from the HPTS hydroxyl. Previous experiments on HPTS/water described the progress from photoexcitation to proton diffusion using kinetic equations with two time constants. The shortest time constant has been interpreted as protonated and photoexcited HPTS evolving into an "associated" state, where the proton is "shared" between the HPTS hydroxyl and an originally hydrogen bonded water. The longer time constant has been interpreted as indicating evolution to a "solvent separated" state where the shared proton undergoes long distance diffusion. In this work, we refine the previous experimental results using very pure HPTS. We then use excited state ab initio molecular dynamics to elucidate the detailed molecular mechanism of aqueous excited state proton transfer in HPTS. We find that the initial excitation results in rapid rearrangement of water, forming a strong hydrogen bonded network (a "water wire") around HPTS. HPTS then deprotonates in ≤3 ps, resulting in a proton that migrates back and forth along the wire before localizing on a single water molecule. We find a near linear relationship between the emission wavelength and proton-HPTS distance over the simulated time scale, suggesting that the emission wavelength can be used as a ruler for the proton distance. Our simulations reveal that the "associated" state corresponds to a water wire with a mobile proton and that the diffusion of the proton away from this water wire (to a generalized "solventseparated" state) corresponds to the longest experimental time constant.

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) was investigated using fast fluorescence measurements. The closely related molecule, 8-methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as protonated (P), associated I (AI), associated II (AII), and deprotonated (D). We find that the simple kinetic scheme of P → AI → AII → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into AI, which mainly goes on to deprotonation (D), but AI can also feed into AII. AII can return to AI or decay to the ground state, but does not deprotonate within experimental error. Quantum chemistry and excited state QM/MM Born–Oppenheimer molecular dynamics simulations indicate that AI and AII are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial, and equatorial, respectively.

The search for prostate cancer biomarkers has received increased attention and several DNA repair related enzymes have been linked to this dysfunction. Here we report a targeted search for single nucleotide polymorphisms (SNPs) and functional impact characterization of human ALKBH family dioxygenases related to prostate cancer. Our results uncovered a SNP of ALKBH7, rs7540, which is associated with prostate cancer disease in a statistically significantly manner in two separate cohorts, and maintained in African American men. Comparisons of molecular dynamics (MD) simulations on the wild-type and variant protein structures indicate that the resulting alteration in the enzyme induces a significant structural change that reduces ALKBH7’s ability to bind its cosubstrate. Experimental spectroscopy studies with purified proteins validate our MD predictions and corroborate the conclusion that this cancer-associated mutation affects productive cosubstrate binding in ALKBH7.

The aggregation of amyloid fibrils can lead to various diseases including Alzheimer’s, Parkinson’s disease, and transmissible spongiform encephalopathy. Amyloid fibrils can develop from a variety of proteins in the body as they misfold into a primarily β-sheet structure and aggregate. Human lysozyme has been shown to have far reaching effects in the human healtha homologous enzyme, hen egg-white lysozyme, has been shown to denature to a primarily β-sheet structure at low pH and high alcohol content solution. We have studied these systems in atomic-level detail with a combination of constant pH and microsecond long molecular dynamics simulation in explicit solvent, which cumulatively total over 10 μs of simulation time. These studies have allowed us to determine two potential unfolding pathways depending on the protonation state of a key glutamic acid residue as well as the effect of solution dynamics and pH on the unfolding process.

DNA synthesis, carried out by DNA polymerases, requires balancing speed and accuracy for faithful replication of the genome. High fidelity DNA polymerases contain a 3′–5′ exonuclease domain that can remove misincorporated nucleotides on the 3′ end of the primer strand, a process called proofreading. The E. coli replicative polymerase, DNA polymerase III, has spatially separated (~55 Å apart) polymerase and exonuclease subunits. Here, we report on the dynamics of E. coli DNA polymerase III proofreading in the presence of its processivity factor, the β2-sliding clamp, at varying base pair termini using single-molecule FRET. We find that the binding kinetics do not depend on the base identity at the termini, indicating a tolerance for DNA mismatches. Further, our single-molecule data and MD simulations show two previously unobserved features: (1) DNA Polymerase III is a highly dynamic protein that adopts multiple conformational states while bound to DNA with matched or mismatched ends, and (2) an exonuclease-deficient DNA polymerase III has reduced conformational flexibility. Overall, our single-molecule experiments provide high time-resolution insight into a mechanism that ensures high fidelity DNA replication to maintain genome integrity.

Translesion DNA synthesis is an essential process that helps resume DNA replication at forks stalled near bulky adducts on the DNA. Benzo[a]pyrene (B[a]P) is a polycyclic aromatic hydrocarbon (PAH) that can be metabolically activated to benzo[a]pyrene diol epoxide (BPDE), which then can react with DNA to form carcinogenic DNA adducts. Here, we have used single-molecule florescence resonance energy transfer (smFRET) experiments, classical molecular dynamics simulations, and nucleotide incorporation assays to investigate the mechanism by which the model Y-family polymerase, Dpo4, bypasses a (+)-cis-B[a]P-N 2-dG adduct in DNA. Our data show that when (+)-cis-B[a]P-N 2-dG is the templating base, the B[a]P moiety is in a non-solvent exposed conformation stacked within the DNA helix, where it effectively blocks nucleotide incorporation across the adduct by Dpo4. However, when the media contains a small amount of dimethyl sulfoxide (DMSO), the adduct is able to move to a solvent-exposed conformation, which enables error-prone DNA replication past the adduct. When the primer terminates across from the adduct position, the addition of DMSO leads to the formation of an insertion complex capable of accurate nucleotide incorporation.