Combined quantum‐mechanical molecular mechanics calculations with NWChem and AMBER: Excited state properties of green fluorescent protein chromophore analogue in aqueous solution

Combined quantum mechanical molecular mechanics (QM/MM) calculations have become a popular methodology for efficient and accurate description of large molecular systems. In this work we introduce our development of a QM/MM framework based on two well‐known codes—NWChem and AMBER. As an initial application area we are focused on excited state properties of small molecules in an aqueous phase using an analogue of the green fluorescent protein (GFP) chromophore as a particular test case. Our approach incorporates high level coupled cluster theory for the analysis of excited states providing a reliable theoretical analysis of effects of an aqueous solvation environment on the photochemical properties of the GFP chromophore. Using a systematic approach, which involves comparison of gas phase and aqueous phase results for different protonation states and conformations, we resolve existing uncertainties regarding the theoretical interpretation of experimental data. We observe that the impact of aqueous environment on charged states generally results in blue shifts of the absorption spectra, but the magnitude of the effect is sensitive to both protonation state and conformation and can be rationalized based on charge movement into the area of higher/lower external electrostatic potentials. At neutral pH levels the experimentally observed absorption signal is most likely coming from the phenol protonated form. Our results also show that the high level electron correlated method is essential for a proper description of excited states of GFP. © 2017 Wiley Periodicals, Inc.

[1]  Peter J Tonge,et al.  Ground state isomerization of a model green fluorescent protein chromophore , 2003, FEBS letters.

[2]  M. Zimmer,et al.  The Role of the Tight-Turn, Broken Hydrogen Bonding, Glu222 and Arg96 in the Post-translational Green Fluorescent Protein Chromophore Formation. , 2008, Chemical physics.

[3]  Anna I Krylov,et al.  On the Nature of an Extended Stokes Shift in the mPlum Fluorescent Protein. , 2015, The journal of physical chemistry. B.

[4]  Weitao Yang,et al.  Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface , 2000 .

[5]  R. A. Jockusch,et al.  Deactivation pathways of an isolated green fluorescent protein model chromophore studied by electronic action spectroscopy. , 2009, Journal of the American Chemical Society.

[6]  Alexander V. Nemukhin,et al.  Modeling of the structure and electronic spectra of green fluorescent protein chromophore , 2008 .

[7]  Piotr Piecuch,et al.  Left‐eigenstate completely renormalized equation‐of‐motion coupled‐cluster methods: Review of key concepts, extension to excited states of open‐shell systems, and comparison with electron‐attached and ionized approaches , 2009 .

[8]  Ville R. I. Kaila,et al.  Electrostatic spectral tuning mechanism of the green fluorescent protein. , 2013, Physical chemistry chemical physics : PCCP.

[9]  Evgeny Epifanovsky,et al.  Quantum Chemical Benchmark Studies of the Electronic Properties of the Green Fluorescent Protein Chromophore. 1. Electronically Excited and Ionized States of the Anionic Chromophore in the Gas Phase. , 2009, Journal of chemical theory and computation.

[10]  Karol Kowalski,et al.  Extension of the method of moments of coupled-cluster equations to excited states: The triples and quadruples corrections to the equation-of-motion coupled-cluster singles and doubles energies , 2002 .

[11]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[12]  Debashree Ghosh,et al.  Extension of the Effective Fragment Potential Method to Macromolecules. , 2016, The journal of physical chemistry. B.

[13]  Anna I Krylov,et al.  Quantum chemistry behind bioimaging: insights from ab initio studies of fluorescent proteins and their chromophores. , 2012, Accounts of chemical research.

[14]  I. Polyakov,et al.  Potential Energy Landscape of the Electronic States of the GFP Chromophore in Different Protonation Forms: Electronic Transition Energies and Conical Intersections. , 2010, Journal of chemical theory and computation.

[15]  Volkhard Helms,et al.  Low-lying electronic excitations of the green fluorescent protein chromophore , 2000 .

[16]  Maodu Chen,et al.  Photoabsorption of green and red fluorescent protein chromophore anions in vacuo. , 2007, Biophysical chemistry.

[17]  John A Tainer,et al.  Understanding GFP posttranslational chemistry: structures of designed variants that achieve backbone fragmentation, hydrolysis, and decarboxylation. , 2006, Journal of the American Chemical Society.

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

[19]  Marat Valiev,et al.  Probing microhydration effect on the electronic structure of the GFP chromophore anion: Photoelectron spectroscopy and theoretical investigations. , 2015, The Journal of chemical physics.

[20]  Evgeny Epifanovsky,et al.  Quantum Chemical Benchmark Studies of the Electronic Properties of the Green Fluorescent Protein Chromophore: 2. Cis-Trans Isomerization in Water. , 2009, Journal of chemical theory and computation.

[21]  B. Reid,et al.  Chromophore formation in green fluorescent protein. , 1997, Biochemistry.

[22]  S J Remington,et al.  Structural basis for dual excitation and photoisomerization of the Aequorea victoria green fluorescent protein. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[23]  T. Martínez,et al.  Conical intersection dynamics in solution: the chromophore of Green Fluorescent Protein. , 2004, Faraday discussions.

[24]  Noel S. Hush,et al.  Solvent effects on the electronic spectra of transition metal complexes. , 2000, Chemical reviews.

[25]  S J Remington,et al.  Structural basis of spectral shifts in the yellow-emission variants of green fluorescent protein. , 1998, Structure.

[26]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[27]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[28]  Nadia Rega,et al.  On the optical absorption of the anionic GFP chromophore in vacuum, solution, and protein. , 2013, Physical chemistry chemical physics : PCCP.

[29]  Alexander V Nemukhin,et al.  Role of zwitterions in kindling fluorescent protein photochemistry. , 2015, The journal of physical chemistry. B.

[30]  Benedetta Mennucci,et al.  Modeling environment effects on spectroscopies through QM/classical models. , 2013, Physical chemistry chemical physics : PCCP.

[31]  Steen Brøndsted Nielsen,et al.  Chromophores of the green fluorescent protein studied in the gas phase , 2002 .

[32]  Debashree Ghosh,et al.  Effect of Solvation on Electron Detachment and Excitation Energies of a Green Fluorescent Protein Chromophore Variant. , 2016, The journal of physical chemistry. B.

[33]  José M. Lluch,et al.  Photoinduced proton transfer from the green fluorescent protein chromophore to a water molecule: analysis of the transfer coordinate , 2004 .

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

[35]  B. A. Lindquist,et al.  Photodynamics in complex environments: ab initio multiple spawning quantum mechanical/molecular mechanical dynamics. , 2009, The journal of physical chemistry. B.

[36]  Alexander V Nemukhin,et al.  Electronic Excitations of the Chromophore from the Fluorescent Protein asFP595 in Solutions. , 2006, Journal of chemical theory and computation.

[37]  R. Parr Density-functional theory of atoms and molecules , 1989 .

[38]  Notker Rösch,et al.  Quantum chemical modeling of structure and absorption spectra of the chromophore in green fluorescent proteins , 1998 .

[39]  Anna I Krylov,et al.  Insight into the common mechanism of the chromophore formation in the red fluorescent proteins: the elusive blue intermediate revealed. , 2012, Journal of the American Chemical Society.

[40]  J A McCammon,et al.  Shedding light on the dark and weakly fluorescent states of green fluorescent proteins. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[41]  William W. Ward,et al.  SPECTRAL PERTURBATIONS OF THE AEQUOREA GREEN‐FLUORESCENT PROTEIN , 1982 .

[42]  Karol Kowalski,et al.  New coupled-cluster methods with singles, doubles, and noniterative triples for high accuracy calculations of excited electronic states. , 2004, The Journal of chemical physics.

[43]  Piotr Piecuch,et al.  Benchmarking the completely renormalised equation-of-motion coupled-cluster approaches for vertical excitation energies , 2015 .

[44]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[45]  John A Tainer,et al.  Structural evidence for an enolate intermediate in GFP fluorophore biosynthesis. , 2006, Journal of the American Chemical Society.

[46]  Marco Garavelli,et al.  Solvent effects on the vibrational activity and photodynamics of the green fluorescent protein chromophore: a quantum-chemical study. , 2005, Journal of the American Chemical Society.

[47]  Roger Y. Tsien,et al.  Improved green fluorescence , 1995, Nature.

[48]  M. Zimmer,et al.  Green fluorescent protein (GFP): applications, structure, and related photophysical behavior. , 2002, Chemical reviews.

[49]  Benedetta Mennucci,et al.  Chromophore-protein coupling beyond nonpolarizable models: understanding absorption in green fluorescent protein. , 2015, Journal of chemical theory and computation.

[50]  John F. Stanton,et al.  The equation of motion coupled‐cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties , 1993 .

[51]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[52]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[53]  Peter J Tonge,et al.  Light-driven decarboxylation of wild-type green fluorescent protein. , 2003, Journal of the American Chemical Society.

[54]  R Y Tsien,et al.  Wavelength mutations and posttranslational autoxidation of green fluorescent protein. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[55]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

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

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

[58]  S. Inouye,et al.  Chemical nature of the light emitter of the Aequorea green fluorescent protein. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[59]  W. M. Westler,et al.  Chemical structure of the hexapeptide chromophore of the Aequorea green-fluorescent protein. , 1993, Biochemistry.

[60]  Sean C. Smith,et al.  Bond selection in the photoisomerization reaction of anionic green fluorescent protein and kindling fluorescent protein chromophore models. , 2008, Journal of the American Chemical Society.

[61]  Notker Rösch,et al.  PROTONATION EFFECTS ON THE CHROMOPHORE OF GREEN FLUORESCENT PROTEIN. QUANTUM CHEMICAL STUDY OF THE ABSORPTION SPECTRUM , 1997 .

[62]  Renato Zenobi,et al.  Absorption of the green fluorescent protein chromophore anion in the gas phase studied by a combination of FTICR mass spectrometry with laser-induced photodissociation spectroscopy , 2011 .

[63]  R. Tsien,et al.  green fluorescent protein , 2020, Catalysis from A to Z.

[64]  J U Andersen,et al.  Absorption spectrum of the green fluorescent protein chromophore anion in vacuo. , 2001, Physical review letters.

[65]  Sriram Krishnamoorthy,et al.  Active-space completely-renormalized equation-of-motion coupled-cluster formalism: Excited-state studies of green fluorescent protein, free-base porphyrin, and oligoporphyrin dimer. , 2010, The Journal of chemical physics.

[66]  Jun Zeng,et al.  Electronic excitations of green fluorescent proteins: Protonation states of chromophore model compound in solutions , 2005, J. Comput. Chem..

[67]  Claudia Filippi,et al.  Absorption Spectrum of the Green Fluorescent Protein Chromophore: A Difficult Case for ab Initio Methods? , 2009, Journal of chemical theory and computation.

[68]  John A Tainer,et al.  Understanding GFP chromophore biosynthesis: controlling backbone cyclization and modifying post-translational chemistry. , 2005, Biochemistry.

[69]  R. Wachter,et al.  Chromogenic cross-link formation in green fluorescent protein. , 2007, Accounts of chemical research.

[70]  Ross C. Walker,et al.  An overview of the Amber biomolecular simulation package , 2013 .

[71]  S. Meech,et al.  Excited state reactions in fluorescent proteins. , 2009, Chemical Society reviews.

[72]  S J Remington,et al.  Crystal structure and photodynamic behavior of the blue emission variant Y66H/Y145F of green fluorescent protein. , 1997, Biochemistry.

[73]  Rodney J. Bartlett,et al.  The equation-of-motion coupled-cluster method: Excitation energies of Be and CO , 1989 .

[74]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[75]  K. Solntsev,et al.  Solvatochromism of the green fluorescence protein chromophore and its derivatives. , 2006, Journal of the American Chemical Society.

[76]  Konstantin A Lukyanov,et al.  Photoinduced Chemistry in Fluorescent Proteins: Curse or Blessing? , 2017, Chemical reviews.

[77]  Teodoro Laino,et al.  Relationship between structure and optical properties in green fluorescent proteins: a quantum mechanical study of the chromophore environment , 2004 .

[78]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[79]  H. Nakatsuji,et al.  Electronic excitations of the green fluorescent protein chromophore in its protonation states: SAC/SAC‐CI study , 2003, J. Comput. Chem..

[80]  Massimo Olivucci,et al.  Origin, nature, and fate of the fluorescent state of the green fluorescent protein chromophore at the CASPT2//CASSCF resolution. , 2004, Journal of the American Chemical Society.

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

[82]  Anna I Krylov,et al.  First-principles characterization of the energy landscape and optical spectra of green fluorescent protein along the A→I→B proton transfer route. , 2013, Journal of the American Chemical Society.

[83]  Haruki Niwa,et al.  Fluorescent properties of model chromophores of tyrosine-66 substituted mutants of Aequorea green fluorescent protein (GEP) , 1998 .

[84]  M. J. Cormier,et al.  Primary structure of the Aequorea victoria green-fluorescent protein. , 1992, Gene.

[85]  Jordan Miles,et al.  Resonantly Enhanced Multiphoton Ionization Spectrum of the Neutral Green Fluorescent Protein Chromophore. , 2014, The journal of physical chemistry letters.

[86]  P. Tonge,et al.  Probing the ground state structure of the green fluorescent protein chromophore using Raman spectroscopy. , 2000, Biochemistry.

[87]  A. Krylov,et al.  Effect of protein environment on electronically excited and ionized states of the green fluorescent protein chromophore. , 2011, The journal of physical chemistry. B.

[88]  Hetal N. Patel,et al.  Reaction progress of chromophore biogenesis in green fluorescent protein. , 2006, Journal of the American Chemical Society.

[89]  Miquel Moreno,et al.  Unveiling how an archetypal fluorescent protein operates: theoretical perspective on the ultrafast excited state dynamics of GFP variant S65T/H148D. , 2015, The journal of physical chemistry. B.

[90]  O. Christiansen,et al.  Photoabsorption studies of neutral green fluorescent protein model chromophores in vacuo. , 2009, Physical chemistry chemical physics : PCCP.

[91]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

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

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

[94]  Riccardo Basosi,et al.  Properties of the emitting state of the green fluorescent protein resolved at the CASPT2//CASSCF/CHARMM level. , 2005, Journal of the American Chemical Society.

[95]  Jeffrey R. Reimers,et al.  Solvent effects on molecular spectra. III. Absorption to and emission from the lowest singlet (n,π*) state of dilute pyrimidine in water , 1993 .

[96]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules , 1971 .

[97]  Notker Rösch,et al.  Absorption spectra of the GFP chromophore in solution: comparison of theoretical and experimental results , 2001 .

[98]  Andreas W. Götz,et al.  An extensible interface for QM/MM molecular dynamics simulations with AMBER , 2014, J. Comput. Chem..