Debashree Ghosh

Noncovalent Interactions in Extended Systems Described by the Effective Fragment Potential Method: Theory and Application to Nucleobase Oligomers

The implementation of the effective fragment potential (EFP) method within the Q-CHEM electronic structure package is presented. The EFP method is used to study noncovalent π-π and hydrogen-bonding interactions in DNA strands. Since EFP is a computationally inexpensive alternative to high-level ab initio calculations, it is possible to go beyond the dimers of nucleic acid bases and to investigate the asymptotic behavior of different components of the total interaction energy. The calculations demonstrated that the dispersion energy is a leading component in π-stacked oligomers of all sizes. Exchange-repulsion energy also plays an important role. The contribution of polarization is small in these systems, whereas the magnitude of electrostatics varies. Pairwise fragment interactions (i.e., the sum of dimer binding energies) were found to be a good approximation for the oligomer energy.

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk.

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymines VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 Å radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.

First-Principle Protocol for Calculating Ionization Energies and Redox Potentials of Solvated Molecules and Ions: Theory and Application to Aqueous Phenol and Phenolate

The effect of hydration on the lowest vertical ionization energy (VIE) of phenol and phenolate solvated in bulk water was characterized using the equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods (referred to as EOM/EFP) and determined experimentally by valence photoemission measurements using microjets and synchrotron radiation. The computed solvent-induced shifts in VIEs (ΔVIEs) are -0.66 and +5.72 eV for phenol and phenolate, respectively. Our best estimates of the absolute values of VIEs (7.9 and 7.7 eV for phenol and phenolate) agree reasonably well with the respective experimental values (7.8 ± 0.1 and 7.1 ± 0.1 eV). The EOM/EFP scheme was benchmarked against full EOM-IP-CCSD using microsolvated phenol and phenolate clusters. A protocol for calculating redox potentials with EOM/EFP was developed based on linear response approximation (LRA) of free energy determination. The oxidation potentials of phenol and phenolate calculated using LRA and EOM/EFP are 1.32 and 0.89 V, respectively; they agree well with experimental values.

A VUV Photoionization and Ab Initio Determination of the Ionization Energy of a Gas Phase Sugar (Deoxyribose).

The ionization energy of gas-phase deoxyribose was determined using tunable vacuum ultraviolet synchrotron radiation coupled to an effusive thermal source. Adiabatic and vertical ionization energies of the ground and first four excited states of α-pyranose, the structure that dominates in the gas phase, were calculated using high-level electronic structure methods. An appearance energy of 9.1(±0.05) eV was recorded, which agrees reasonably well with a theoretical value of 8.8 eV for the adiabatic ionization energy. A clear picture of the dissociative photoionization dynamics of deoxyribose emerges from the fragmentation pattern recorded using mass spectrometry and from ab initio molecular dynamics calculations. The experimental threshold 9.4 (±0.05) eV for neutral water elimination upon ionization is captured well in the calculations, and qualitative insights are provided by molecular orbital analysis and molecular dynamics snapshots along the reaction coordinate.

Effective fragment potential method in Q‐CHEM: A guide for users and developers

A detailed description of the implementation of the effective fragment potential (EFP) method in the Q‐CHEM electronic structure package is presented. The Q‐CHEM implementation interfaces EFP with standard quantum mechanical (QM) methods such as Hartree–Fock, density functional theory, perturbation theory, and coupled‐cluster methods, as well as with methods for electronically excited and open‐shell species, for example, configuration interaction, time‐dependent density functional theory, and equation‐of‐motion coupled‐cluster models. In addition to the QM/EFP functionality, a “fragment‐only” feature is also available (when the system is described by effective fragments only). To aid further developments of the EFP methodology, a detailed description of the C++ classes and EFP modules workflow is presented. The EFP input structure and EFP job options are described. To assist setting up and performing EFP calculations, a collection of Perl service scripts is provided. The precomputed EFP parameters for standard fragments such as common solvents are stored in Q‐CHEMs auxiliary library; they can be easily invoked, similar to specifying standard basis sets. The instructions for generating user‐defined EFP parameters are given. Fragments positions can be specified by their center of mass coordinates and Euler angles. The interface with the IQMOL and WEBMO software is also described.

Extension of the Effective Fragment Potential Method to Macromolecules

The effective fragment potential (EFP) approach, which can be described as a nonempirical polarizable force field, affords an accurate first-principles treatment of noncovalent interactions in extended systems. EFP can also describe the effect of the environment on the electronic properties (e.g., electronic excitation energies and ionization and electron-attachment energies) of a subsystem via the QM/EFP (quantum mechanics/EFP) polarizable embedding scheme. The original formulation of the method assumes that the system can be separated, without breaking covalent bonds, into closed-shell fragments, such as solvent and solute molecules. Here, we present an extension of the EFP method to macromolecules (mEFP). Several schemes for breaking a large molecule into small fragments described by EFP are presented and benchmarked. We focus on the electronic properties of molecules embedded into a protein environment and consider ionization, electron-attachment, and excitation energies (single-point calculations only). The model systems include chromophores of green and red fluorescent proteins surrounded by several nearby amino acid residues and phenolate bound to the T4 lysozyme. All mEFP schemes show robust performance and accurately reproduce the reference full QM calculations. For further applications of mEFP, we recommend either the scheme in which the peptide is cut along the Cα-C bond, giving rise to one fragment per amino acid, or the scheme with two cuts per amino acid, along the Cα-C and Cα-N bonds. While using these fragmentation schemes, the errors in solvatochromic shifts in electronic energy differences (excitation, ionization, electron detachment, or electron-attachment) do not exceed 0.1 eV. The largest error of QM/mEFP against QM/EFP (no fragmentation of the EFP part) is 0.06 eV (in most cases, the errors are 0.01-0.02 eV). The errors in the QM/molecular mechanics calculations with standard point charges can be as large as 0.3 eV.

Toward Understanding the Redox Properties of Model Chromophores from the Green Fluorescent Protein Family: An Interplay between Conjugation, Resonance Stabilization, and Solvent Effects

The redox properties of model chromophores from the green fluorescent protein family are characterized computationally using density functional theory with a long-range corrected functional, the equation-of-motion coupled-cluster method, and implicit solvation models. The analysis of electron-donating abilities of the chromophores reveals an intricate interplay between the size of the chromophore, conjugation, resonance stabilization, presence of heteroatoms, and solvent effects. Our best estimates of the gas-phase vertical/adiabatic detachment energies of the deprotonated (i.e., anionic) model red, green, and blue chromophores are 3.27/3.15, 2.79/2.67, and 2.75/2.35 eV, respectively. Vertical/adiabatic ionization energies of the respective protonated (i.e., neutral) species are 7.64/7.35, 7.38/7.15, and 7.70/7.32 eV, respectively. The standard reduction potentials (E(red)(0)) of the anionic (Chr•/Chr–) and neutral (Chr+•/Chr) model chromophores in acetonitrile are 0.34/1.40 V (red), 0.22/1.24 V (green), and −0.12/1.02 V (blue), suggesting, counterintuitively, that the red chromophore is more difficult to oxidize than the green and blue ones (in both neutral and deprotonated forms). The respective redox potentials in water follow a similar trend but are more positive than the acetonitrile values.

Effect of Solvation on Electron Detachment and Excitation Energies of a Green Fluorescent Protein Chromophore Variant

Hybrid quantum mechanics/molecular mechanics (QM/MM) is applied to the fluorinated green fluorescent protein (GFP) chromophore (DFHBDI) in its deprotonated form to understand the solvatochromic shifts in its vertical detachment energy (VDE) and vertical excitation energy (VEE). This variant of the GFP chromophore becomes fluorescent in an RNA environment and has a wide range of applications in biomedical and biochemical fields. From microsolvation studies, we benchmark (with respect to full QM) the accuracy of our QM/MM calculations with effective fragment potential (EFP) as the MM method of choice. We show that while the solvatochromic shift in the VEE is minimal (0.1 eV blue shift) and its polarization component is only 0.03 eV, the effect of the solvent on the VDE is quite large (3.85 eV). We also show by accurate calculations on the solvatochromic shift of the VDE that polarization accounts for ∼0.23 eV and therefore cannot be neglected. The effect of the counterions on the VDE of the deprotonated chromophore in solvation is studied in detail, and a charge-smearing scheme is suggested for charged chromophores.

Effects of the benzoxazole group on green fluorescent protein chromophore crystal structure and solid state photophysics

Four benzoxazole-substituted GFP chromophores that differ by the length of their alkyl chain (from C1 to C12) were synthesized. In solution, the four compounds showed identical spectroscopic behavior, emitting blue light with moderate quantum yield. In the solid state, the butyl, pentyl and dodecyl derivatives strongly emitted orange light, while the methyl derivative was only weakly emissive. Based on the X-ray data and DFT calculations, emission in the solid state was explained by the formation of excimers. A very unusual “hot-dog”-type excimer was found for the dodecyl derivative, in which two overlapping chromophores are separated by an alkyl chain.

Feasibility of Ionization-Mediated Pathway for Ultraviolet-Induced Melanin Damage

Melanin is the pigment found in human skin that is responsible for both photoprotection and photodamage. Recently there have been reports that greater photodamage of DNA occurs when cells containing melanin are irradiated with ultraviolet (UV) radiation, thus suggesting that the photoproducts of melanin cause DNA damage. Photoionization processes have also been implicated in the photodegradation of melanin. However, not much is known about the oxidation potential of melanin and its monomers. In this work we calculate the ionization energies of monomers, dimers, and few oligomers of eumelanin to estimate the threshold energy required for the ionization of eumelanin. We find that this threshold is within the UV-B region for eumelanin. We also look at the charge and spin distributions of the various ionized states of the monomers that are formed to understand which of the ionization channels might favor monomerization from a covalent dimer.