Revealing intrinsic changes of DNA induced by spore photoproduct lesion through computer simulation.

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.

Reaction Path-Force Matching in Collective Variables: Determining Ab Initio QM/MM Free Energy Profiles by Fitting Mean Force.

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.

Trafficking of ciliary G protein-coupled receptors.

In the last decade highly conserved cellular appendages called cilia have enjoyed a renewed interest from basic, biomedical scientists, and clinicians alike. This interest has grown upon the elucidation that cilia throughout the body serve as important sensory and signaling centers in both development and adult homeostasis. Furthermore, the identification of several rare genetic disorders associated with cilia dysfunction has broadened the field. However, even though their potential role in human health and disease is now recognized many basic questions about their functions remain. This chapter seeks to explore the trafficking of cilia-specific G protein-coupled receptors (GPCRs) and discusses several model systems in which this has been explored. We open the chapter by briefly discussing cilia and GPCRs then begin discussing some aspects of rhodopsin trafficking, arguably the most well studied of cilia GPCRs. We continue with sections on neuronal cilia and olfactory cilia receptor trafficking. Finally, we conclude with the emerging area of dynamic ciliary GPCR trafficking and speculate about future directions and some of the questions that remain for ciliary GPCRs.