QM/MM calculation of solvent effects on absorption spectra of guanine

Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited‐state calculations were extracted from ground‐state molecular dynamics (MD) simulations using the self‐consistent‐charge density functional tight binding (SCC‐DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited‐state calculations used time‐dependent density functional theory (TDDFT) and the DFT‐based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto‐N7H and keto‐N9H guanine, with particular focus on solvent effects in the low‐energy spectrum of the keto‐N9H tautomer. When compared with the vertical excitation energies of gas‐phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC‐DFTB‐based rather than B3LYP‐based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas‐phase and solution spectra, typically ca. 0.1–0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent‐induced structural changes and from electrostatic solute–solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Calculating absorption shifts for retinal proteins: computational challenges

Rhodopsins can modulate the optical properties of their chromophores over a wide range of wavelengths. The mechanism for this spectral tuning is based on the response of the retinal chromophore to external stress and the interaction with the charged, polar, and polarizable amino acids of the protein environment and is connected to its large change in dipole moment upon excitation, its large electronic polarizability, and its structural flexibility. In this work, we investigate the accuracy of computational approaches for modeling changes in absorption energies with respect to changes in geometry and applied external electric fields. We illustrate the high sensitivity of absorption energies on the ground-state structure of retinal, which varies significantly with the computational method used for geometry optimization. The response to external fields, in particular to point charges which model the protein environment in combined quantum mechanical/molecular mechanical (QM/MM) applications, is a crucial feature, which is not properly represented by previously used methods, such as time-dependent density functional theory (TDDFT), complete active space self-consistent field (CASSCF), and Hartree-Fock (HF) or semiempirical configuration interaction singles (CIS). This is discussed in detail for bacteriorhodopsin (bR), a protein which blue-shifts retinal gas-phase excitation energy by about 0.5 eV. As a result of this study, we propose a procedure which combines structure optimization or molecular dynamics simulation using DFT methods with a semiempirical or ab initio multireference configuration interaction treatment of the excitation energies. Using a conventional QM/MM point charge representation of the protein environment, we obtain an absorption energy for bR of 2.34 eV. This result is already close to the experimental value of 2.18 eV, even without considering the effects of protein polarization, differential dispersion, and conformational sampling.     

Investigation of the S<Subscript>0</Subscript>S<Subscript>1</Subscript> excitation in bacteriorhodopsin with the ONIOM(MO:MM) hybrid method

 We have investigated the S₀ and S₁ electronic states in bacteriorhodopsin using a variety of QM/MM levels. The decomposition of the calculated excitation energies into electronic and electrostatic components shows that the interaction of the chromophore with the protein electric field increases the excitation energy, while polarization effects are negligible. Therefore, the experimentally observed reduction in excitation energy from solution phase to protein environment (the Opsin shift) does not come from the electrostatic interaction with the protein environment, but from either the interaction ofthe chromophore with the solvent or counter ion, or structural effects. Our high-level ONIOM(TD– B3LYP:Amber) calculation predicts the excitation energy within 8 kcal/mol from experiment, the discrepancy probably being caused by the neglect of polarization of the protein environment. In addition, we have shown that the level of optimization is extremely critical for the calculation of accurate excitation energies in bacteriorhodopsin. Received: 13 October 2001 / Accepted: 6 September 2002 / Published online: 3 February 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 222nd National Meeting of the American Chemical Society, 2001 Correspondence to: K. Morokuma e-mail: morokuma@emory.edu     

Coupled-cluster studies of the lowest excited states of the 11-cis-retinal chromophore

The first few excited states of the 11-cis-retinal (PSB11) chromophore have been studied at the coupled-cluster approximative singles and doubles (CC2) level using triple-zeta quality basis sets augmented with double sets of polarisation functions. The two lowest vertical excitation energies of 2.14 and 3.21 eV are in good agreement with recently reported experimental values of 2.03 and 3.18 eV obtained in molecular beam measurements. Calculations at the time-dependent density functional theory (TDDFT) level using the B3LYP hybrid functional yield vertical excitation energies of 2.34 and 3.10 eV for the two lowest states. Zero-point vibrational energy (ZPVE) corrections of -0.09 and -0.17 eV were deduced from the harmonic vibrational frequencies for the ground and excited states calculated at the density functional theory (DFT) and TDDFT level, respectively, using the B3LYP hybrid functional.     

Electronic spectra of nitroethylene

A systematic study of the electronic excited states of nitroethylene (C₂H₃NO₂) was carried out using the approximate coupled‐cluster singles‐and‐doubles approach with the resolution of the identity (RI‐CC2), the time dependent density functional theory with the CAMB3LYP functional (TDDFT/CAMB3LYP) and the DFT multireference configuration interaction (DFT/MRCI) method. Vertical transition energies and optical oscillator strengths were computed for a maximum of 20 singlet transitions. Semiclassical simulations of the ultraviolet (UV) spectra were performed at the RI‐CC2 and DFT/MRCI levels. The main features in the UV spectrum were assigned to a weak n‐π* transition, and two higher energy π_CC+O‐π* bands. These characteristics are common to molecules containing NO₂ groups. Simulated spectra are in good agreement with the experimental spectrum. The energy of the bands in the DFT/MRCI simulation agrees quite well with the experiment, although it overestimates the band intensities. RI‐CC2 produced intensities comparable to the experiment, but the bands were blue shifted. A strong π_CC+O‐π* band, not previously measured, was found in the 8–9 eV range. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.     

Insight into the common mechanism of the chromophore formation in the red fluorescent proteins: the elusive blue intermediate revealed

Understanding the chromophore maturation process in fluorescent proteins is important for the design of proteins with improved properties. Here, we present the results of electronic structure calculations identifying the nature of a blue intermediate, a key species in the process of the red chromophore formation in DsRed, TagRFP, fluorescent timers, and PAmCherry. The chromophore of the blue intermediate has a structure in which the π-system of the imidazole ring is extended by the acylimine bond, which can be represented by the model N-[(5-hydroxy-1H-imidazole-2yl)methylidene]acetamide (HIMA) compound. Ab initio and QM/MM calculations of the isolated model and protein-bound (mTagBFP) chromophores identify the anionic form of HIMA as the only structure that has absorption that is consistent with the experiment and is stable in the protein binding pocket. The anion and zwitterion are the only protonation forms of HIMA whose absorption (421 and 414 nm, or 2.95 and 3.00 eV) matches the experimental spectrum of the blue form in DsRed (the absorption maximum is 408 nm or 3.04 eV) and mTagBFP (400 nm or 3.10 eV). The QM/MM optimization of the protein-bound anionic form results in a structure that is close to the X-ray one, whereas the zwitterionic chromophore is unstable in the protein binding pocket and undergoes prompt proton transfer. The computed excitation energy of the protein-bound anionic form of the mTagBFP-like chromophore (3.04 eV) agrees with the experimental absorption spectrum of the protein. The DsRed-like chromophore formation in red fluorescent proteins is revisited on the basis of ab initio results and verified by directed mutagenesis revealing a key role of the amino acid residue 70, which is the second after the chromophore tripeptide, in the formation process.     

Benchmark studies on the building blocks of DNA. 1. Superiority of coupled cluster methods in describing the excited states of nucleobases in the Franck-Condon region

Equation of motion excitation energy coupled-cluster (EOMEE-CC) methods including perturbative triple excitations have been used to set benchmark results for the excitation energy and oscillator strength of the building units of DNA, i.e., cytosine, guanine, adenine and thymine. In all cases the lowest twelve transitions have been considered including valence and Rydberg ones. Triple-ζ basis sets with diffuse functions have been used and the results are compared to CC2, CASPT2, TDDFT, and DFT/MRCI results from the literature. The results clearly show that it is only the EOMEE-CCSD(T) that is capable of providing accuracy of about 0.1 eV. EOMEE-CCSD systematically overshoots the energy of all types of transitions by 0.1-0.3 eV, whereas CC2 is surprisingly accurate for ππ* transitions but fails (often badly) for nπ* and Rydberg transitions. DFT and CASPT2 seem to give reliable results for the lowest transition, but the error increases fast with the excitation level. The differences in the excitation energies often change the energy ordering of the states, which should even influence the conclusions of excited state dynamics obtained with these approximate methods. The results call for further benchmark calculations on larger building blocks of DNA (nucleosides, basis pairs) at the CCSD(T) level.     

Modeling of phytochrome absorption spectra

Phytochromes constitute one of the six well‐characterized families of photosensory proteins in Nature. From the viewpoint of computational modeling, however, phytochromes have been the subject of much fewer studies than most other families of photosensory proteins, which is likely a consequence of relevant high‐resolution structural data becoming available only in recent years. In this work, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are used to calculate UV‐vis absorption spectra of Deinococcus radiodurans bacteriophytochrome. We investigate how the choice of QM/MM methodology affects the resulting spectra and demonstrate that QM/MM methods can reproduce the experimental absorption maxima of both the Q and Soret bands with an accuracy of about 0.15 eV. Furthermore, we assess how the protein environment influences the intrinsic absorption of the bilin chromophore, with particular focus on the Q band underlying the primary photochemistry of phytochromes. © 2013 Wiley Periodicals, Inc.