A simple polarizable continuum solvation model for electrolyte solutions

Isomerization of bioactive acylhydrazones triggered by light or thiols

The acylhydrazone unit is well represented in screening databases used to find ligands for biological targets, and numerous bioactive acylhydrazones have been reported. However, potential E/Z isomerization of the C=N bond in these compounds is rarely examined when bioactivity is assayed. Here we analysed two ortho-hydroxylated acylhydrazones discovered in a virtual drug screen for modulators of N-methyl-D-aspartate receptors and other bioactive hydroxylated acylhydrazones with structurally defined targets reported in the Protein Data Bank. We found that ionized forms of these compounds, which are populated under laboratory conditions, photoisomerize readily and the isomeric forms have markedly different bioactivity. Furthermore, we show that glutathione, a tripeptide involved with cellular redox balance, catalyses dynamic E⇄Z isomerization of acylhydrazones. The ratio of E to Z isomers in cells is determined by the relative stabilities of the isomers regardless of which isomer was applied. We conclude that E/Z isomerization may be a common feature of the bioactivity observed with acylhydrazones and should be routinely analysed.

State-specific solvation for restricted active space spin-flip (RAS-SF) wave functions based on the polarizable continuum formalism

The restricted active space spin-flip (RAS-SF) formalism is a particular form of single-reference configuration interaction that can describe some forms of strong correlation at a relatively low cost and which has recently been formulated for the description of charge-transfer excited states. Here, we introduce both equilibrium and nonequilibrium versions of a state-specific solvation correction for vertical transition energies computed using RAS-SF wave functions, based on the framework of a polarizable continuum model (PCM). Ground-state polarization is described using the solvent's static dielectric constant and in the nonequilibrium solvation approach that polarization is modified upon vertical excitation using the solvent's optical dielectric constant. Benchmark calculations are reported for well-studied models of photo-induced charge transfer, including naphthalene dimer, C₂H₄⋯C₂F₄, pentacene dimer, and perylene diimide (PDI) dimer, several of which are important in organic photovoltaic applications. For the PDI dimer, we demonstrate that the charge-transfer character of the excited states is enhanced in the presence of a low-dielectric medium (static dielectric constant ɛ₀ = 3) as compared to a gas-phase calculation (ɛ₀ = 1). This stabilizes mechanistic traps for singlet fission and helps to explain experimental singlet fission rates. We also examine the effects of nonequilibrium solvation on charge-separated states in an intramolecular singlet fission chromophore, where we demonstrate that the energetic ordering of the states changes as a function of solvent polarity. The RAS-SF + PCM methodology that is reported here provides a framework to study charge-separated states in solution and in photovoltaic materials.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

The role of substituted pyridine Schiff bases as ancillary ligands in the optical properties of a new series of fac-rhenium(i) tricarbonyl complexes: a theoretical view

Over the last few years, luminescent Re(i) tricarbonyl complexes have been increasingly proposed as fluorophores suitable for fluorescence microscopy to visualize biological structures and cells. In this sense, incorporating an asymmetrical pyridine Schiff base (PSB) as the ancillary ligand strongly modifies the staining and luminescent properties of Re(i) tricarbonyl complexes. In this work, we analyzed two series of Re(i) tricarbonyl complexes with their respective PSB ligands: (1) fac-[Re(CO)3(2,2'-bpy)(PSB)]1+ and (2) fac-[Re(CO)3(4,4'-bis(ethoxycarbonyl)-2,2'-bpy)(PSB)]1+, where the PSB exhibits substitutions at positions 4 or 6 in the phenolic ring with methyl or halogen substituents. Thus, we performed computational relativistic DFT and TDDFT studies to determine their optical properties. The ten complexes analyzed showed absorption in the visible light range. Furthermore, our analyses, including zero-field splitting (ZFS), allowed us to determine that the low-lying excited state locates below the 3LLCT states. Interestingly, seven of the ten analyzed complexes, whose corresponding PSB harbors an intramolecular hydrogen bond (IHB), exhibited luminescent emission that could be suitable for biological purposes: large Stokes shift, emission in the range 600-700 nm and τ in the order of 10-2 to 10-3 s. Conversely, the three complexes lacking the IHB due to two halogen substituents in the corresponding PSB showed a predicted emission with the lowest triplet excited state energy entering the NIR region. The main differences in the complexes' photophysical behavior have been explained by the energy gap law and time-resolved luminescence. These results emphasize the importance of choosing suitable substituents at the 4 and 6 positions in the phenolic ring of the PSB, which determine the presence of the IHB since they modulate the luminescence properties of the Re(i) core. Therefore, this study could predict Re(i) tricarbonyl complexes' properties, considering the desired emission features for biological and other applications.

Do water's electrons care about electrolytes?

Ions have a profound effect on the geometrical structure of liquid water and an aqueous environment is known to change the electronic structure of ions. Here we combine photoelectron spectroscopy measurements from liquid microjets with molecular dynamical and quantum chemical calculations to address the reverse question, to what extent do ions affect the electronic structure of liquid water? We study aqueous solutions of sodium iodide (NaI) over a wide concentration range, from nearly pure water to 8 M solutions, recording spectra in the 5 to 60 eV binding energy range to include all water valence and the solute Na+ 2p, I- 4d, and I- 5p orbital ionization peaks. We observe that the electron binding energies of the solute ions change only slightly as a function of electrolyte concentration, less than 150 ± 60 meV over an ∼8 M range. Furthermore, the photoelectron spectrum of liquid water is surprisingly mildly affected as we transform the sample from a dilute aqueous salt solution to a viscous, crystalline-like phase. The most noticeable spectral changes are a negative binding energy shift of the water 1b2 ionizing transition (up to -370 ± 60 meV) and a narrowing of the flat-top shape water 3a1 ionization feature (up to 450 ± 90 meV). A novel computationally efficient technique is introduced to calculate liquid-state photoemission spectra using small clusters from molecular dynamics (MD) simulations embedded in dielectric continuum. This theoretical treatment captured the characteristic positions and structures of the aqueous photoemission peaks, reproducing the experimentally observed narrowing of the water 3a1 feature and weak sensitivity of the water binding energies to electrolyte concentration. The calculations allowed us to attribute the small binding energy shifts to ion-induced disruptions of intermolecular electronic interactions. Furthermore, they demonstrate the importance of considering concentration-dependent screening lengths for a correct description of the electronic structure of solvated systems. Accounting for electronic screening, the calculations highlight the minimal effect of electrolyte concentration on the 1b1 binding energy reference, in accord with the experiments. This leads us to a key finding that the isolated, lowest-binding-energy, 1b1, photoemission feature of liquid water is a robust energetic reference for aqueous liquid microjet photoemission studies.

A new set of atomic radii for accurate estimation of solvation free energy by Poisson-Boltzmann solvent model

The Poisson-Boltzmann implicit solvent (PB) is widely used to estimate the solvation free energies of biomolecules in molecular simulations. An optimized set of atomic radii (PB radii) is an important parameter for PB calculations, which determines the distribution of dielectric constants around the solute. We here present new PB radii for the AMBER protein force field to accurately reproduce the solvation free energies obtained from explicit solvent simulations. The presented PB radii were optimized using results from explicit solvent simulations of the large systems. In addition, we discriminated PB radii for N- and C-terminal residues from those for nonterminal residues. The performances using our PB radii showed high accuracy for the estimation of solvation free energies at the level of the molecular fragment. The obtained PB radii are effective for the detailed analysis of the solvation effects of biomolecules.

Progress in developing Poisson-Boltzmann equation solvers

This review outlines the recent progress made in developing more accurate and efficient solutions to model electrostatics in systems comprised of bio-macromolecules and nano-objects, the last one referring to objects that do not have biological function themselves but nowadays are frequently used in biophysical and medical approaches in conjunction with bio-macromolecules. The problem of modeling macromolecular electrostatics is reviewed from two different angles: as a mathematical task provided the specific definition of the system to be modeled and as a physical problem aiming to better capture the phenomena occurring in the real experiments. In addition, specific attention is paid to methods to extend the capabilities of the existing solvers to model large systems toward applications of calculations of the electrostatic potential and energies in molecular motors, mitochondria complex, photosynthetic machinery and systems involving large nano-objects.