Browsing by Author "Hege, Mellisa"

Now showing 1 - 3 of 3
  • Thumbnail Image
    Item
    Molecular dynamics simulations of spore photoproduct containing DNA systems
    (2023-05) Hege, Mellisa; Pu, Jingzhi; Blacklock, Brenda; Georgiadis, Millie; Deng, Yongming
    Bacterial endospores have been a topic of research interest over the last several decades given their high resistance to ultraviolet (UV) damage. Unlike vegetative bacterial cells, which form cyclobutane pyrimidine dimers (CPD) and pyrimidine 6-4 pyrimidone photoproducts (6-4PPs) as the major product upon UV irradiation, endospore bacteria form a spore photoproduct (5-(R-thyminyl)-5,6-dihydrothymine or SP) as the major product. Vegetative bacteria cells are subject to regular cell activities and processes such as division and deoxyribonucleic acid (DNA) replication, which are prone to damage from UV exposure. However, in endospores, which have a largely anhydrous inner environment, the DNA remains dormant when bound to spore-specific small acid-soluble proteins (SASP) and dipicolinic acid, making spores highly resistant to radiation, heat, desiccation, and chemical harm. During spore germination, SP lesions in DNA are repaired by a distinctive repair enzyme, spore photoproduct lyase (SPL). In this thesis, molecular dynamics (MD) simulations were carried out to (i) examine how the formation of the SP lesion in DNA affects the global and local structural properties of duplex DNA and (ii) study how this lesion is recognized and repaired in endospore. The first part of this work was focused on designing and developing a structurally and dynamically stable model for dinucleotide SP molecule (TpT), which was subsequently used as an SP patch incorporated into duplex DNA. Computationally, this requires modifications of the bond and nonbonded force field parameters. The stability of the patch and developed parameters was tested via solution-phase MD simulations for the SP lesion incorporated within the B-DNA dodecamer duplex (PDB 463B). The second part involved applying the new SP patch to simulate the crystallographic structure of the DNA oligomer containing SP lesions. Solution-phase MD simulations were performed for the SP-containing DNA oligomers (modeled based on PDB 4M94) and compared to the simulations of the native structure (PDB 4M95). Our analysis of the MD trajectories revealed a range of SP-induced structural and dynamical changes, including the weakened hydrogen bonds at the SP sites, increased DNA bending, and distinct conformational stability and distribution. In the third part of this thesis project, we carried out MD simulations of SP-containing DNA bound with SASPs to examine how the DNA interacts differently with SASP in the presence and absence of the SP lesion. The simulation results suggested decreased electrostatic and hydrogen bonding interactions between SASP and the damaged DNA at the SP site compared to the undamaged DNA-protein complex. In addition, decreased helicity percentage was observed in the SASPs that directly interact with the SP lesion. The last part of this this thesis work focused on the SP-dimer flipping mechanism, as the lesion is likely flipped out to its extrahelical state to be recognized and repaired by SPL. Using steered molecular dynamic (SMD) simulations and a pseudo-dihedral angle reaction coordinate, we obtained possible SP flipping pathways both in the presence and absence of SASP. Collectively, these simulation results lend new perspectives toward understanding the unique behavior of the SP lesion within the DNA duplex and the nucleoprotein complex. They also provide new insights into how the SP lesion is efficiently recognized and repaired during spore germination.
  • Thumbnail Image
    Item
    Reaction Path-Force Matching in Collective Variables: Determining Ab Initio QM/MM Free Energy Profiles by Fitting Mean Force
    (American Chemical Society, 2021) Kim, Bryant; Snyder, Ryan; Nagaraju, Mulpuri; Zhou, Yan; Ojeda-May, Pedro; Keeton, Seth; Hege, Mellisa; Shao, Yihan; Pu, Jingzhi; Chemistry and Chemical Biology, School of Science
    First-principles determination of free energy profiles for condensed-phase chemical reactions is hampered by the daunting costs associated with configurational sampling on ab initio quantum mechanical/molecular mechanical (AI/MM) potential energy surfaces. Here, we report a new method that enables efficient AI/MM free energy simulations through mean force fitting. In this method, a free energy path in collective variables (CVs) is first determined on an efficient reactive aiding potential. Based on the configurations sampled along the free energy path, correcting forces to reproduce the AI/MM forces on the CVs are determined through force matching. The AI/MM free energy profile is then predicted from simulations on the aiding potential in conjunction with the correcting forces. Such cycles of correction-prediction are repeated until convergence is established. As the instantaneous forces on the CVs sampled in equilibrium ensembles along the free energy path are fitted, this procedure faithfully restores the target free energy profile by reproducing the free energy mean forces. Due to its close connection with the reaction path-force matching (RP-FM) framework recently introduced by us, we designate the new method as RP-FM in collective variables (RP-FM-CV). We demonstrate the effectiveness of this method on a type-II solution-phase SN2 reaction, NH3 + CH3Cl (the Menshutkin reaction), simulated with an explicit water solvent. To obtain the AI/MM free energy profiles, we employed the semiempirical AM1/MM Hamiltonian as the base level for determining the string minimum free energy pathway, along which the free energy mean forces are fitted to various target AI/MM levels using the Hartree-Fock (HF) theory, density functional theory (DFT), and the second-order Møller-Plesset perturbation (MP2) theory as the AI method. The forces on the bond-breaking and bond-forming CVs at both the base and target levels are obtained by force transformation from Cartesian to redundant internal coordinates under the Wilson B-matrix formalism, where the linearized FM is facilitated by the use of spline functions. For the Menshutkin reaction tested, our FM treatment greatly reduces the deviations on the CV forces, originally in the range of 12-33 to ∼2 kcal/mol/Å. Comparisons with the experimental and benchmark AI/MM results, tests of the new method under a variety of simulation protocols, and analyses of the solute-solvent radial distribution functions suggest that RP-FM-CV can be used as an efficient, accurate, and robust method for simulating solution-phase chemical reactions.
  • Thumbnail Image
    Item
    Revealing intrinsic changes of DNA induced by spore photoproduct lesion through computer simulation
    (Elsevier, 2023-05) Hege, Mellisa; Li, Lei; Pu, Jingzhi; Chemistry and Chemical Biology, School of Science
    In bacterial endospores, a cross-linked thymine dimer, 5-thyminyl-5,6-dihydrothymine, commonly referred to as the spore photoproduct (SP), is found as the dominant DNA photo lesion under UV radiation. During spore germination, SP is faithfully repaired by the spore photoproduct lyase (SPL) for normal DNA replication to resume. Despite this general mechanism, the exact way in which SP modifies the duplex DNA structure so that the damaged site can be recognized by SPL to initiate the repair process is still unclear. A previous X-ray crystallographic study, which used a reverse transcriptase as a DNA host template, captured a protein-bound duplex oligonucleotide containing two SP lesions; the study showed shortened hydrogen bonds between the AT base pairs involved in the lesions and widened minor grooves near the damaged sites. However, it remains to be determined whether the results accurately reflect the conformation of SP-containing DNA (SP-DNA) in its fully hydrated pre-repair form. To uncover the intrinsic changes in DNA conformation caused by SP lesions, we performed molecular dynamics (MD) simulations of SP-DNA duplexes in aqueous solution, using the nucleic acid portion of the previously determined crystal structure as a template. After MD relaxation, our simulated SP-DNAs showed weakened hydrogen bonds at the damaged sites compared to those in the undamaged DNA. Our analyses of the MD trajectories revealed a range of local and global structural distortions of DNA induced by SP. Specifically, the SP region displays a greater tendency to adopt an A-like-DNA conformation, and curvature analysis revealed an increase in the global bending compared to the canonical B-DNA. Although these SP-induced DNA conformational changes are relatively minor, they may provide a sufficient structural basis for SP to be recognized by SPL during the lesion repair process.