Local electric fields and molecular properties in heterogeneous environments through polarizable embedding.

In spectroscopies, the local field experienced by a molecule embedded in an environment will be different from the externally applied electromagnetic field, and this difference may significantly alter the response and transition properties of the molecule. The polarizable embedding (PE) model has previously been developed to model the local field contribution stemming from the direct molecule-environment coupling of the electromagnetic response properties of molecules in solution as well as in heterogeneous environments, such as proteins. Here we present an extension of this approach to address the additional effective external field effect, i.e., the manifestations of the environment polarization induced by the external field, which allows for the calculation of properties defined in terms of the external field. Within a response framework, we report calculations of the one- and two-photon absorption (1PA and 2PA, respectively) properties of PRODAN-methanol clusters as well as the fluorescent protein DsRed. Our results demonstrate the necessity of accounting for both the dynamical reaction field and effective external field contributions to the local field in order to reproduce full quantum chemical reference calculations. For the lowest π→π* transition in DsRed, inclusion of effective external field effects gives rise to a 1.9- and 3.5-fold reduction in the 1PA and 2PA cross-sections, respectively. The effective external field is, however, strongly influenced by the heterogeneity of the protein matrix, and the resulting effect can lead to either screening or enhancement depending on the nature of the transition under consideration.

[1]  L. Jensen,et al.  Microscopic and macroscopic polarization within a combined quantum mechanics and molecular mechanics model. , 2005, The Journal of chemical physics.

[2]  Jacopo Tomasi,et al.  Vibrational circular dichroism within the polarizable continuum model: a theoretical evidence of conformation effects and hydrogen bonding for (S)-(-)-3-butyn-2-ol in CCl4 solution , 2002 .

[3]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[4]  Jacopo Tomasi,et al.  The Cotton–Mouton effect of furan and its homologues in the gas phase, for the pure liquids and in solution , 2003 .

[5]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[6]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[7]  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.

[8]  S. Lukyanov,et al.  Fluorescent proteins from nonbioluminescent Anthozoa species , 1999, Nature Biotechnology.

[9]  D. M. Bishop,et al.  Effective polarizabilities and local field corrections for nonlinear optical experiments in condensed media , 1998 .

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

[11]  S. Boxer,et al.  Extreme electric fields power catalysis in the active site of ketosteroid isomerase , 2014, Science.

[12]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[13]  H. Ågren,et al.  Density-functional-theory study of the electric-field-induced second harmonic generation (EFISHG) of push–pull phenylpolyenes in solution , 2006 .

[14]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[15]  Xavier Saint Raymond,et al.  Elementary Introduction to the Theory of Pseudodifferential Operators , 2018 .

[16]  Olav Vahtras LoProp for Dalton , 2014 .

[17]  Rayleigh-Schrödinger perturbation theory for nonlinear Schrödinger equations with linear perturbation , 1993 .

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

[19]  Kenneth Ruud,et al.  Benchmarking two-photon absorption cross sections: performance of CC2 and CAM-B3LYP. , 2015, Physical chemistry chemical physics : PCCP.

[20]  Luca Frediani,et al.  Solvent effects on Raman optical activity spectra calculated using the polarizable continuum model. , 2006, The journal of physical chemistry. A.

[21]  D. Matyushov,et al.  Cavity field in liquid dielectrics , 2008 .

[22]  S. Pipolo,et al.  The cavity electromagnetic field within the polarizable continuum model of solvation. , 2014, The Journal of chemical physics.

[23]  H. Ågren,et al.  Relation between Nonlinear Optical Properties of Push-Pull Molecules and Metric of Charge Transfer Excitations. , 2015, Journal of chemical theory and computation.

[24]  Arieh Warshel,et al.  Energetics and Dynamics of Enzymatic Reactions , 2001 .

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

[26]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[27]  Jacopo Tomasi,et al.  On the Calculation of Infrared Intensities in Solution within the Polarizable Continuum Model , 2000 .

[28]  P. T. V. Duijnen,et al.  Molecular and Atomic Polarizabilities: Thole's Model Revisited , 1998 .

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

[30]  Jacopo Tomasi,et al.  Cavity field effects within a polarizable continuum model of solvation: Application to the calculation of electronic circular dichroism spectra of R‐(+)‐3‐methyl‐cyclopentanone , 2011 .

[31]  J. Tomasi,et al.  Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules , 2002 .

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

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

[34]  S. Lukyanov,et al.  Fluorescent proteins and their applications in imaging living cells and tissues. , 2010, Physiological reviews.

[35]  M. Head‐Gordon,et al.  Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes. , 2004, Journal of the American Chemical Society.

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

[37]  F. J. Luque,et al.  Generalized linear response approximation in discrete methods , 1997 .

[38]  F. Lipparini,et al.  Towards an accurate description of anharmonic infrared spectra in solution within the polarizable continuum model: reaction field, cavity field and nonequilibrium effects. , 2011, The Journal of chemical physics.

[39]  J. Tainer,et al.  Crystallographic structures of Discosoma red fluorescent protein with immature and mature chromophores: linking peptide bond trans-cis isomerization and acylimine formation in chromophore maturation. , 2005, Biochemistry.

[40]  P. Åstrand,et al.  Local electric field factors by a combined charge-transfer and point–dipole interaction model , 2015 .

[41]  Jógvan Magnus Haugaard Olsen,et al.  Molecular-Level Insight into the Spectral Tuning Mechanism of the DsRed Chromophore. , 2012, The journal of physical chemistry letters.

[42]  János G. Ángyán,et al.  Common theoretical framework for quantum chemical solvent effect theories , 1992 .

[43]  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 .

[44]  Jing Kong,et al.  Lennard–Jones parameters for the combined QM/MM method using the B3LYP/6‐31G*/AMBER potential , 2005, J. Comput. Chem..

[45]  S. Boxer,et al.  Measuring electric fields and noncovalent interactions using the vibrational stark effect. , 2015, Accounts of chemical research.

[46]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[47]  M. Drobizhev,et al.  Color hues in red fluorescent proteins are due to internal quadratic Stark effect. , 2009, The journal of physical chemistry. B.

[48]  Jacopo Tomasi,et al.  An attempt to bridge the gap between computation and experiment for nonlinear optical properties: Macroscopic susceptibilities in solution , 2000 .

[49]  Arieh Warshel,et al.  Apparent NAC effect in chorismate mutase reflects electrostatic transition state stabilization. , 2003, Journal of the American Chemical Society.

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

[51]  A. Warshel,et al.  Electrostatic basis for enzyme catalysis. , 2006, Chemical reviews.

[52]  Ove Christiansen,et al.  Multiphoton transition moments and absorption cross sections in coupled cluster response theory employing variational transition moment functionals , 1998 .

[53]  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 .

[54]  H. A. Lorentz (1) Anfangsgründe der Maxwellschen Theorie, verknüpft mit der Elektronentheorie (2) The Theory of Electrons, and its Applications to the Phenomena of Light and Radiant Heat , 2009, Nature.

[55]  S. Tretiak,et al.  Enhanced Two‐Photon Absorption of Organic Chromophores: Theoretical and Experimental Assessments , 2008 .

[56]  Jacob Kongsted,et al.  Molecular Properties through Polarizable Embedding , 2011 .

[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]  Kenneth Ruud,et al.  Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone. , 2015, The Journal of chemical physics.

[60]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[61]  Jacopo Tomasi,et al.  Theoretical Approach to the Calculation of Vibrational Raman Spectra in Solution within the Polarizable Continuum Model , 2001 .

[62]  Jacopo Tomasi,et al.  On the Calculation of Local Field Factors for Microscopic Static Hyperpolarizabilities of Molecules in Solution with the Aid of Quantum-Mechanical Methods , 1998 .

[63]  S. Pipolo,et al.  The cavity electromagnetic field within the polarizable continuum model of solvation: An application to the real-time time dependent density functional theory , 2014 .

[64]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[65]  B. Honig,et al.  Electrostatic interactions in membranes and proteins. , 1986, Annual review of biophysics and biophysical chemistry.

[66]  S. Benkovic,et al.  Probing the Electrostatics of Active Site Microenvironments along the Catalytic Cycle for Escherichia coli Dihydrofolate Reductase , 2014, Journal of the American Chemical Society.

[67]  M. Drobizhev,et al.  Two-photon absorption properties of fluorescent proteins , 2011, Nature Methods.

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

[69]  H. Nakatsuji,et al.  Color tuning in photofunctional proteins. , 2011, Chemphyschem : a European journal of chemical physics and physical chemistry.

[70]  J. Tomasi,et al.  Calculation of nonlinear optical susceptibilities of pure liquids within the Polarizable Continuum Model: the effect of the macroscopic nonlinear polarization at the output frequency , 2003 .

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

[72]  H. Nakamura,et al.  Roles of electrostatic interaction in proteins , 1996, Quarterly Reviews of Biophysics.

[73]  Roberto Cammi,et al.  Quantum mechanical polarizable continuum model approach to the Kerr effect of pure liquids. , 2005, The journal of physical chemistry. B.

[74]  J. R. Carl,et al.  Atom dipole interaction model for molecular polarizability. Application to polyatomic molecules and determination of atom polarizabilities , 1972 .

[75]  S. Hammes-Schiffer,et al.  Catalytic efficiency of enzymes: a theoretical analysis. , 2013, Biochemistry.

[76]  M. Drobizhev,et al.  Absolute two-photon absorption spectra and two-photon brightness of orange and red fluorescent proteins. , 2009, The journal of physical chemistry. B.

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

[78]  L. Onsager Electric Moments of Molecules in Liquids , 1936 .