Computational Approach for Studying Optical Properties of DNA Systems in Solution.

In this paper we present a study of the methodological aspects regarding calculations of optical properties for DNA systems in solution. Our computational approach will be built upon a fully polarizable QM/MM/Continuum model within a damped linear response theory framework. In this approach the environment is given a highly advanced description in terms of the electrostatic potential through the polarizable embedding model. Furthermore, bulk solvent effects are included in an efficient manner through a conductor-like screening model. With the aim of reducing the computational cost we develop a set of averaged partial charges and distributed isotropic dipole-dipole polarizabilities for DNA suitable for describing the classical region in ground-state and excited-state calculations. Calculations of the UV-spectrum of the 2-aminopurine optical probe embedded in a DNA double helical structure are presented. We show that inclusion of polarizabilities in the embedding potential stemming from the DNA double helix is of crucial importance, while the water cluster surrounding the DNA system is well represented using a continuum approach.

[1]  A. Klamt.  Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena , 1995 .

[2]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[3]  Thanh N. Truong,et al.  Optimized atomic radii for quantum dielectric continuum solvation models , 1995 .

[4]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[5]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[6]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[7]  Jacopo Tomasi,et al.  Excited states and solvatochromic shifts within a nonequilibrium solvation approach: A new formulation of the integral equation formalism method at the self-consistent field, configuration interaction, and multiconfiguration self-consistent field level , 1998 .

[8]  D. M. Bishop,et al.  Near-resonant absorption in the time-dependent self-consistent field and multiconfigurational self-consistent field approximations , 2001 .

[9]  Benedetta Mennucci,et al.  Toward a Unified Modeling of Environment and Bridge-Mediated Contributions to Electronic Energy Transfer: A Fully Polarizable QM/MM/PCM Approach. , 2012, Journal of chemical theory and computation.

[10]  Jógvan Magnus Haugaard Olsen,et al.  Averaged Solvent Embedding Potential Parameters for Multiscale Modeling of Molecular Properties. , 2016, Journal of chemical theory and computation.

[11]  J. Kelly,et al.  [Ru(TAP)2(dppz)]2+: a DNA intercalating complex, which luminesces strongly in water and undergoes photo-induced proton-coupled electron transfer with guanosine-5'-monophosphate. , 2004, Dalton transactions.

[12]  D C Ward,et al.  Fluorescence studies of nucleotides and polynucleotides. I. Formycin, 2-aminopurine riboside, 2,6-diaminopurine riboside, and their derivatives. , 1969, The Journal of biological chemistry.

[13]  F. Javier Luque,et al.  Towards a molecular dynamics consensus view of B-DNA flexibility , 2008, Nucleic acids research.

[14]  Qiang Cui,et al.  Combining implicit solvation models with hybrid quantum mechanical/molecular mechanical methods: A critical test with glycine , 2002 .

[15]  Mark S. Gordon,et al.  Self‐consistent molecular orbital methods. XXIII. A polarization‐type basis set for second‐row elements , 1982 .

[16]  Jim A. Thomas,et al.  Kinetically inert transition metal complexes that reversibly bind to DNA. , 2003, Chemical Society reviews.

[17]  Jacob Kongsted,et al.  Solvation Effects on Electronic Transitions: Exploring the Performance of Advanced Solvent Potentials in Polarizable Embedding Calculations. , 2011, Journal of chemical theory and computation.

[18]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[19]  Jacob Kongsted,et al.  Polarizable embedding with a multiconfiguration short-range density functional theory linear response method. , 2015, The Journal of chemical physics.

[20]  Jacob Kongsted,et al.  Damped Response Theory in Combination with Polarizable Environments: The Polarizable Embedding Complex Polarization Propagator Method. , 2014, Journal of chemical theory and computation.

[21]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[22]  U. Rothlisberger,et al.  Insights into intrastrand cross-link lesions of DNA from QM/MM molecular dynamics simulations. , 2012, Journal of the American Chemical Society.

[23]  Patrick Norman,et al.  A perspective on nonresonant and resonant electronic response theory for time-dependent molecular properties. , 2011, Physical chemistry chemical physics : PCCP.

[24]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[25]  Timothy Clark,et al.  Efficient diffuse function‐augmented basis sets for anion calculations. III. The 3‐21+G basis set for first‐row elements, Li–F , 1983 .

[26]  Christian Griesinger,et al.  2-Aminopurine incorporation perturbs the dynamics and structure of DNA. , 2010, Angewandte Chemie.

[27]  Pär Söderhjelm,et al.  How accurate can a force field become? A polarizable multipole model combined with fragment-wise quantum-mechanical calculations. , 2009, The journal of physical chemistry. A.

[28]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[29]  M. Surin,et al.  DNA Electronic Circular Dichroism on the Inter-Base Pair Scale: An Experimental-Theoretical Case Study of the AT Homo-Oligonucleotide. , 2015, The journal of physical chemistry letters.

[30]  Vincenzo Barone,et al.  Linear Response Theory and Electronic Transition Energies for a Fully Polarizable QM/Classical Hamiltonian. , 2012, Journal of chemical theory and computation.

[31]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

[32]  Hui Li,et al.  The polarizable continuum model (PCM) interfaced with the fragment molecular orbital method (FMO) , 2006, J. Comput. Chem..

[33]  V. Barone,et al.  Quantum Calculation of Molecular Energies and Energy Gradients in Solution by a Conductor Solvent Model , 1998 .

[34]  Thanh N. Truong,et al.  A new method for incorporating solvent effect into the classical, ab initio molecular orbital and density functional theory frameworks for arbitrary shape cavity , 1995 .

[35]  Kazuo Kitaura,et al.  Extending the power of quantum chemistry to large systems with the fragment molecular orbital method. , 2007, The journal of physical chemistry. A.

[36]  D. M. Bishop,et al.  Nonlinear response theory with relaxation: the first-order hyperpolarizability. , 2005, The Journal of chemical physics.

[37]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[38]  F. J. Luque,et al.  Frontiers in molecular dynamics simulations of DNA. , 2012, Accounts of chemical research.

[39]  F. Jensen.  Unifying General and Segmented Contracted Basis Sets. Segmented Polarization Consistent Basis Sets. , 2014, Journal of chemical theory and computation.

[40]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[41]  Poul Jørgensen,et al.  On the Efficiency of Algorithms for Solving Hartree-Fock and Kohn-Sham Response Equations. , 2011, Journal of chemical theory and computation.

[42]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[43]  Jógvan Magnus Haugaard Olsen,et al.  PERI-CC2: A Polarizable Embedded RI-CC2 Method. , 2012, Journal of chemical theory and computation.

[44]  Kurt V. Mikkelsen,et al.  On the importance of excited state dynamic response electron correlation in polarizable embedding methods , 2012, J. Comput. Chem..

[45]  Jan H. Jensen,et al.  FragIt: A Tool to Prepare Input Files for Fragment Based Quantum Chemical Calculations , 2012, PloS one.

[46]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[47]  Hui Li,et al.  Note: FixSol solvation model and FIXPVA2 tessellation scheme. , 2012, The Journal of chemical physics.

[48]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[49]  X. Assfeld,et al.  A QM/MM study of the absorption spectrum of harmane in water solution and interacting with DNA: the crucial role of dynamic effects. , 2013, The journal of physical chemistry. B.

[50]  Jacob Kongsted,et al.  The polarizable embedding coupled cluster method. , 2011, Journal of Chemical Physics.

[51]  Hui Li,et al.  Quantum mechanical/molecular mechanical/continuum style solvation model: linear response theory, variational treatment, and nuclear gradients. , 2009, The Journal of chemical physics.

[52]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[53]  Hongbo Zhu,et al.  QuanPol: A full spectrum and seamless QM/MM program , 2013, J. Comput. Chem..

[54]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .

[55]  Jan H. Jensen,et al.  Effective fragment molecular orbital method: a merger of the effective fragment potential and fragment molecular orbital methods. , 2010, The journal of physical chemistry. A.

[56]  M. Gordon,et al.  Polarization energy gradients in combined quantum mechanics, effective fragment potential, and polarizable continuum model calculations. , 2007, The Journal of chemical physics.

[57]  Luca Frediani,et al.  The Dalton quantum chemistry program system , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[58]  F. Jensen,et al.  Electrostatic Potential of Insulin: Exploring the Limitations of Density Functional Theory and Force Field Methods. , 2013, Journal of chemical theory and computation.

[59]  K. Ruud,et al.  Multiconfigurational self-consistent field linear response for the polarizable continuum model: Theory and application to ground and excited-state polarizabilities of para-nitroaniline in solution , 2003 .

[60]  Roland Lindh,et al.  Local properties of quantum chemical systems: the LoProp approach. , 2004, The Journal of chemical physics.

[61]  Jógvan Magnus Haugaard Olsen,et al.  Accuracy of Protein Embedding Potentials: An Analysis in Terms of Electrostatic Potentials. , 2015, Journal of chemical theory and computation.

[62]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[63]  Jacopo Tomasi,et al.  An Integrated Effective Fragment—Polarizable Continuum Approach to Solvation: Theory and Application to Glycine , 2002 .

[64]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[65]  D. Case,et al.  A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA , 2009, Nucleic acids research.

[66]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[67]  Jacob Kongsted,et al.  The multi-configuration self-consistent field method within a polarizable embedded framework. , 2013, The Journal of chemical physics.

[68]  Jacob Kongsted,et al.  Excited States in Solution through Polarizable Embedding , 2010 .

[69]  B. Nordén,et al.  Electronic Transition Moments of 2-Aminopurine , 1997 .

[70]  J. Widom,et al.  Role of DNA sequence in nucleosome stability and dynamics , 2001, Quarterly Reviews of Biophysics.

[71]  Trygve Helgaker,et al.  Recent advances in wave function-based methods of molecular-property calculations. , 2012, Chemical reviews.

[72]  Anita C Jones,et al.  2-aminopurine as a fluorescent probe of DNA conformation and the DNA–enzyme interface , 2015, Quarterly Reviews of Biophysics.

[73]  Luca Frediani,et al.  Excitation energies in solution: the fully polarizable QM/MM/PCM method. , 2011, The journal of physical chemistry. B.