Symmetric versus asymmetric discretization of the integral equations in polarizable continuum solvation models

Synthesis, design, and optimization of a potent and selective series of pyridylpiperazines as promising antimalarial agents

An optimization of the pyridylpiperazine series against Plasmodium falciparum has been performed, exploring a structure-activity relationship carried out on the toluyl fragment of hit 1, a compound with low micromolar activity against Plasmodium falciparum discovered by high-throughput screening. After confirming the crucial role played by this aryl fragment in the antiplasmodial activity, the replacement of the ortho-methyl substituent of 1 by halogenated ones led to an improvement for four analogs, either in terms of potency, expected pharmacokinetics profile, or both. Further introduction of endocyclic nitrogens in this fragment identified two more optimized compounds, 20 and 23, which are expected to be much more metabolically stable than 1. Additional assessment of the cytotoxicity, Ligand Lipophilic Efficiency, potency against the chloroquine-resistant Dd2 strain and in silico ADMET predictions revealed a satisfactory profile for most compounds, ultimately identifying the four optimized compounds 7, 9, 20 and 23 as promising compounds for further lead optimization of this series against Plasmodium falciparum.

Iodine(I) pnictogenate complexes as Iodination reagents

Halogen(I) complexes are widely used as halogenation reagents and traditionally feature homoleptic stabilising Lewis bases, though the recent revitalisation of iodine(I) carboxylate chemistry has provided isolable examples of heteroleptic iodine(I) complexes. This work reports iodine(I) pnictogenate complexes stabilised by a Lewis base (L), Ph2P(O)O─I─L, synthesised via cation exchange from the silver(I) precursor, (Ph2P(O)OAg)n. The complexes were characterised in both solution (1H, 1H-15N HMBC, 31P) and the solid state, and supplemented computationally by DFT studies. Interestingly, these iodine(I) pnictogenates demonstrate a range of stabilities, and have been found to excel as iodination reagents in comparison to carbonyl hypoiodites, with comparable reactivity to the eponymous Barluenga's reagent in the iodination of antipyrine.

Hydroboration of vinylsilanes providing diversity-oriented hydrophobic building blocks for biofunctional molecules

Hydroboration of vinylsilanes with BH3 affords two silylethanol regioisomers. Herein, we investigated the regioisomeric ratio of hydroboration products from various vinylsilanes, focusing on the characteristic reaction profile. All investigated vinylsilanes afforded both regioisomers, and greater bulkiness increased the proportion of the Markovnikov products. The obtained silylethanols were used as hydrophobic building blocks for constructing nuclear progesterone receptor (PR) modulators. Notably, structural conversions from an α-isomer (silylethan-1-oxy derivative) to a β-isomer (2-silylethoxy derivative) caused complete activity-switching from a PR agonist to an antagonist. Our results indicate that silylethanols are useful for structural development, and vinylsilanes are a versatile source of hydrophobic building blocks for obtaining biofunctional molecules.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Pyridine Iodine(I) Cations: Kinetic Trapping as a Sulfonate Complexes

Seven pyridine iodine(I) sulfonate complexes were prepared and isolated at low temperatures and characterized by X-ray diffraction analysis. The inherently instable pyridine iodine(I) cations are stabilized by an oxygen of sulfonate anions via the I⋅⋅⋅O halogen bond. In these complexes, the iodine atom of the pyridine iodine(I) cation acts as an electron acceptor and the sulfonate oxygen as the electron donor. These complexes are stable enough in the crystalline state, yet decompose rapidly under ambient conditions, also being unstable in solution. The (pyridine)N-I bond lengths [2.140(3)-2.197(2) Å] and the I⋅⋅⋅O halogen bonds [2.345(6)-2.227(3) Å] are analogous to (imide)N-I⋅⋅⋅O-N-pyridine uncharged halogen-bonded complexes formed from N-haloimides and pyridine N-oxides, thus confirming the existence of elusive pyridine iodine(I) cation.

Controlling the repair mechanisms of oxetanes through functional group substitution

Intersystem crossing (ISC) plays a key role in the photolysis processes of oxetanes formed by benzophenone (BP)-like and thymine structures. In this work, we systematically explored the photophysical processes of oxetanes and ring-splitting products and investigated the effect of substituents on the repair mechanisms of oxetanes. The regioselectivity of oxetanes (head-to-head, HH and head-to-tail, HT) and the electron-donating and electron-withdrawing substituents, including CH₃, OCH₃ and NO₂, were considered. It was found that the substituents influence the ISC rates of these compounds more by changing their spin-orbit coupling (SOC) coefficients rather than energy gaps. The SOC coefficients of HH-oxetanes are more affected by these groups than HT-oxetanes and products, and they have greater ISC rates on the whole. Besides, the insertion of substituents can alter the radiative and nonradiative decay rates, thereby transforming the photoinduced cycloreversion mechanisms of oxetanes. The ring-splitting reactions of non-substituted oxetanes could occur via two pathways of singlet and triplet manifolds. Furthermore, oxetanes with NO₂ at the X site have the largest ISC rates but hardly undergo repair processes, while the introduction of electron-donating substituents can effectively promote the repair of oxetanes. The singlet ring-splitting reactions of HH-oxetanes are more inclined to occur after introducing CH₃ and OCH₃ at two sites. However, HT-oxeatnes with CH₃ are more likely to undergo triplet repair processes and OCH₃-substituted structures tend to originate cycloreversion in the singlet manifolds. Moreover, the introduction of CH₃ and OCH₃ at the Y site rather than the X site can more significantly accelerate the repair processes of HH-oxetanes. Contrarily, HT-oxetanes with electron-donating groups at the X site exhibit faster repair rates than those at the Y site. We hope this work can provide valuable insights into BP-like drugs and photosensitive DNA repair.

Fragment-Based Calculations of Enzymatic Thermochemistry Require Dielectric Boundary Conditions

Electronic structure calculations on enzymes require hundreds of atoms to obtain converged results, but fragment-based approximations offer a cost-effective solution. We present calculations on enzyme models containing 500-600 atoms using the many-body expansion, comparing to benchmarks in which the entire enzyme-substrate complex is described at the same level of density functional theory. When the amino acid fragments contain ionic side chains, the many-body expansion oscillates under vacuum boundary conditions but rapid convergence is restored using low-dielectric boundary conditions. This implies that full-system calculations in the gas phase are inappropriate benchmarks for assessing errors in fragment-based approximations. A three-body protocol retains sub-kilocalorie per mole fidelity with respect to a supersystem calculation, as does a two-body calculation combined with a full-system correction at a low-cost level of theory. These protocols pave the way for application of high-level quantum chemistry to large systems via rigorous, ab initio treatment of many-body polarization.

Carbonyl hypoiodites from pivalic and trimesic acid and their silver(I) intermediates

The first tris(O-I-N) carbonyl hypoiodites have been synthesised based on trimesic acid and pyridine or 4-methylpyridine, with their structures definitively confirmed by single crystal X-ray diffraction (SCXRD). The more soluble carbonyl hypoiodites based on pivalic acid have also been studied via NMR, SCXRD, and computational analyses, enabling the study of the direct silver(I) precursor and intermediates of the resulting carbonyl hypoiodites generated using a range of substituted pyridines.

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.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0-0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang-Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang-Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0-0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

Rethinking Uncaging: A New Antiaromatic Photocage Driven by a Gain of Resonance Energy

Photoactivatable compounds for example photoswitches or photolabile protecting groups (PPGs, photocages) for spatiotemporal light control, play a crucial role in different areas of research. For each application, parameters such as the absorption spectrum, solubility in the respective media and/or photochemical quantum yields for several competing processes need to be optimized. The design of new photochemical tools therefore remains an important task. In this study, we exploited the concept of excited-state-aromaticity, first described by N. Colin Baird in 1971, to investigate a new class of photocages, based on cyclic, ground-state-antiaromatic systems. Several thio- and nitrogen-functionalized compounds were synthesized, photochemically characterized and further optimized, supported by quantum chemical calculations. After choosing the optimal scaffold, which shows an excellent uncaging quantum yield of 28 %, we achieved a bathochromic shift of over 100 nm, resulting in a robust, well accessible, visible light absorbing, compact new photocage with a clean photoreaction and a high quantum product (ϵ⋅Φ) of 893 M-1  cm-1 at 405 nm.

Bioisosteric Modification of To042: Synthesis and Evaluation of Promising Use-Dependent Inhibitors of Voltage-Gated Sodium Channels

Three analogues of To042, a tocainide-related lead compound recently reported for the treatment of myotonia, were synthesized and evaluated in vitro as skeletal muscle sodium channel blockers possibly endowed with enhanced use-dependent behavior. Patch-clamp experiments on hNav1.4 expressed in HEK293 cells showed that N-[(naphthalen-1-yl)methyl]-4-[(2,6-dimethyl)phenoxy]butan-2-amine, the aryloxyalkyl bioisostere of To042, exerted a higher use-dependent block than To042 thus being able to preferentially block the channels in over-excited membranes while preserving healthy tissue function. It also showed the lowest active transport across BBB according to the results of P-glycoprotein (P-gp) interacting activity evaluation and the highest cytoprotective effect on HeLa cells. Quantum mechanical calculations and dockings gave insights on the most probable conformation of the aryloxyalkyl bioisostere of To042 in solution and the target residues involved in the binding, respectively. Both approaches indicated the conformations that might be adopted in both the unbound and bound state of the ligand. Overall, N-[(naphthalen-1-yl)methyl]-4-[(2,6-dimethyl)phenoxy]butan-2-amine exhibits an interesting toxico-pharmacological profile and deserves further investigation.

Electronic Spectra of C60 Films Using Screened Range Separated Hybrid Functionals

We study computationally the electronic spectra of C₆₀ thin films using the recently developed density functional theory (DFT) framework combining a screened range separated hybrid (SRSH) functional with a polarizable continuum model (PCM). The SRSH-PCM approach achieves excellent correspondence between the frontier orbital's energy levels and the ionization potential and electron affinity of the molecular system at the condensed phase and consequently leads to high quality electronic excitation energies when used in time-dependent DFT calculations. Our calculated excited states reproduce the experimentally main reported spectral peaks at the 3.6-4.6 eV energy range and when addressing excitonic effects also reproduce the red-shifted spectral feature. Notably, we analyze the low-lying peak at 2.7 eV and associate it to an excitonic state.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces "natural charge-transfer orbitals" that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I^-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with the experiment, and only one-third of the water molecules in the I^-(H₂O)₉₆ simulation cell need to be described with LR-TDDFT to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.

Densely Functionalized 2-Methylideneazetidines: Evaluation as Antibacterials

Twenty-two novel, variously substituted nitroazetidines were designed as both sulfonamide and urethane vinylogs possibly endowed with antimicrobial activity. The compounds under study were obtained following a general procedure recently developed, starting from 4-nitropentadienoates deriving from a common β-nitrothiophenic precursor. While being devoid of any activity against fungi and Gram-negative bacteria, most of the title compounds performed as potent antibacterial agents on Gram-positive bacteria (E. faecalis and three strains of S. aureus), with the most potent congener being the 1-(4-chlorobenzyl)-3-nitro-4-(p-tolyl)azetidine 22, which displayed potency close to that of norfloxacin, the reference antibiotic (minimum inhibitory concentration values 4 and 1-2 μg/mL, respectively). Since 22 combines a relatively efficient activity against Gram-positive bacteria and a cytotoxicity on eucharyotic cells only at 4-times higher concentrations (inhibiting concentration on 50% of the cultured eukaryotic cells: 36 ± 10 μM, MIC: 8.6 μM), it may be considered as a promising hit compound for the development of a new series of antibacterials selectively active on Gram-positive pathogens. The relatively concise synthetic route described herein, based on widely available starting materials, could feed further structure-activity relationship studies, thus allowing for the fine investigation and optimization of the toxico-pharmacological profile.

Xanthenylacetic Acid Derivatives Effectively Target Lysophosphatidic Acid Receptor 6 to Inhibit Hepatocellular Carcinoma Cell Growth

Despite the increasing incidence of hepatocellular carcinoma (HCC) worldwide, current pharmacological treatments are still unsatisfactory. We have previously shown that lysophosphatidic acid receptor 6 (LPAR6) supports HCC growth and that 9-xanthenylacetic acid (XAA) acts as an LPAR6 antagonist inhibiting HCC growth without toxicity. Here, we synthesized four novel XAA derivatives, (±)-2-(9H-xanthen-9-yl)propanoic acid (compound 4 - MC9), (±)-2-(9H-xanthen-9-yl)butanoic acid (compound 5 - MC6), (±)-2-(9H-xanthen-9-yl)hexanoic acid (compound 7 - MC11), and (±)-2-(9H-xanthen-9-yl)octanoic acid (compound 8 - MC12, sodium salt) by introducing alkyl groups of increasing length at the acetic α-carbon atom. Two of these compounds were characterized by X-ray powder diffraction and quantum mechanical calculations, while molecular docking simulations suggested their enantioselectivity for LPAR6. Biological data showed anti-HCC activity for all XAA derivatives, with the maximum effect observed for MC11. Our findings support the view that increasing the length of the alkyl group improves the inhibitory action of XAA and that enantioselectivity can be exploited for designing novel and more effective XAA-based LPAR6 antagonists.

Experimental and Theoretical Characterization of Ultrafast Water-Soluble Photochromic Photoacids

UV-visible transient absorption spectroscopy and quantum mechanical simulations are combined to elucidate the photochemical mechanism of two metastable merocyanine/spiropyran photoacids, 2-[(E)-2-(2-hydroxyphenyl)ethenyl]-3,3-dimethyl-1-(3-sulfopropyl)-3H-indol-1-ium (phenylhydroxy-MCH) and 2-[(E)-2-(1H-indazol-7-yl)ethenyl]-3-(3-sulfopropyl)-1,3-benzothiazol-3-ium (indazole-MCH). Transient absorption spectra demonstrate that trans-acid isomerization to the cis form results in deprotonation on a picosecond time scale. Ring closure to form spiropyran follows promptly from the appropriate conformation or follows at longer time delays (≫3.5 ns) following a barrier crossing for single-bond isomerization to the appropriate conformation. Consistent with the results of Berton et al. [ Chem. Sci. 2020, 11, 8457-8468] , we find that cis-phenylhydroxy-MCH is a stronger acid than trans-phenylhydroxy-MCH. The decrease in pK_a upon isomerization is further investigated to benchmark quantum chemical methods for their accuracy. Calculations were performed with nine levels of theory including continuum solvent models and explicit water. The calculations are not sufficient to describe the ΔpK_a following isomerization of these photoacids, and more work is necessary to properly evaluate the physical basis for the acidity of the cis photoacids.

How Well Does a Solvated Octa-acid Capsule Shield the Embedded Chromophore? A Computational Analysis Based on an Anisotropic Dielectric Continuum Model

The optical properties of chromophores embedded in a water-solvated dimer of octa-acid that forms a molecular-shaped capsule are investigated. In particular, we address the anisotropic dielectric environment that appears to blue-shift excitation energies compared to the free aqueous chromophores. Recently we reported that using an effective scalar dielectric constant ε ≈ 3 appears to reproduce the measured spectra of the embedded coumarins, suggesting that the capsule provides a significant, albeit not perfect, screening of the aqueous dielectric environment. Here, we report absorption energies using a theoretical treatment that includes continuum solvation affected by an anisotropic dielectric function reflecting the high-dielectric environment outside of the capsule and the low-dielectric region within. We report time-dependent density functional theory calculations using a range-separated functional with the Poisson boundary conditions that model the anisotropic dielectric environment. Our calculations find that the anisotropic environment due to the water-solvated hydrophobic capsule is equivalent to a homogeneous effective dielectric constant of ≈3. The calculated values also appear to reproduce measured absorption of the embedded coumarin, where we study the effect of the hydrophobic capsule on the excited state.

A theoretical study on electronic spectra of a novel series of metal substituted boron subphthalocyanine chloride

Boron subphthalocyanine chloride has been extensively studied by experimentalists and computational chemists due to its unique optical and electronic properties. It has been practical to modify the optical and physical properties of subphthalocyanine through axial, peripheral, and center substitutions or ring expansion. However, there have been few investigations on the substitution of central boron atom. In the present work, a new metal-substituted (center substitution of boron atom) series of boron subphthalocyanine chloride (metal = Fe, Co, Ni, Cu, and Zn) are theoretically designed utilizing modern density functional theory. The optimized results of this series in gas phase and with polarizable continuum model show that they may be chemically stable, and the predicted order of the stability of MSubPC is Fe>Cu>Ni>Co>Zn. Also, this new series of MSubPC molecules all becomes more non-planar and has much smaller dipole moments, which imply that they may be feasible for blend with organic acceptors. The HOMO-LUMO energy gaps of MSubPC (M=Co, Ni, Cu) are smaller than that of subPC. Furthermore, the wavelength of simulated absorption peaks of ZnSubPC and NiSubPC is red-shifted with respect to prototype subPC molecule in the visible region, and FeSubPC has noticeably stronger absorption strength than subPC because its excitation involves more orbital transitions and d electrons. The work here shows a new way to design photoelectric materials based on subphthalocyanine with center metal substitution.

Regioselectivity enhancement in synthesis of [70]fullerene derivatives by introduction of a branched structure

Regio-purity in [70]fullerene derivatives is of great importance to improve the power conversion efficiencies of organic photovoltaics. We found that the introduction of a branched structure to [70]PCBM enhanced the yield of α-isomers. The effect of the steric group in the reaction mechanism was theoretically investigated and a difference in the activation energies within specific pathways was revealed.

SE2 reaction in noncarbon system: Metal-halide catalysis for dehydrogenation of ammonia borane

An electrophilic substitution (SE) reaction of BN isosteres has been investigated for the dehydrogenation of ammonia borane (AB) by metal chlorides (MCl2) using various ab initio calculations. In contrast to the typical SE reaction occurring at the carbon atom, the nitrogen atom in AB serves as the reaction center for the SE reaction with the boron moiety as the leaving group when the MCl2 approaches the AB. The SE2 backside reaction is favored as a trigger step for the dehydrogenation of AB by the MCl2 The SE2 reaction is found for 3d-transition-metal chlorides (e.g., FeCl2, CoCl2, NiCl2, CuCl2, and ZnCl2), while PdCl2 leads to the dehydrogenation of AB by a direct B-H σ-bond activation, similar to most organometallic catalysts. Interestingly, the polymerization of AB promoted by MCl2 can be explained with the similar SE2 mechanism, and the dehydrogenation of the BN derivative 3-methyl-1,2-BN-cyclopentane (CBN) bearing a carbon backbone ring also follows the SE2 reaction. In particular, the experimental observation that the use of metal-chloride catalysis decreases the by-products obtained during the hydrogenation of AB can be explained by our mechanism involving the SE2 reaction. This work is helpful for the development of novel metal-halide catalysts for practical hydrogen storage materials, including the BN moiety.

Tuned range separated hybrid functionals for solvated low bandgap oligomers

The description of charge transfer excitations has long been a challenge to time dependent density functional theory. The recently developed concept of "optimally tuned range separated hybrid (OT-RSH) functionals" has proven to describe charge transfer excitations accurately in many cases. However, describing solvated or embedded systems is yet a challenge. This challenge is not only computational but also conceptual, because the tuning requires identifying a specific orbital, typically the highest occupied one of the molecule under study. For solvated molecules, this orbital may be delocalized over the solvent. We here demonstrate that one way of overcoming this problem is to use a locally projected self-consistent field diagonalization on an absolutely localized molecular orbital expansion. We employ this approach to determine ionization energies and the optical gap of solvated oligothiophenes, i.e., paradigm low gap systems that are of relevance in organic electronics. Dioxane solvent molecules are explicitly represented in our calculations, and the ambiguities of straightforward parameter tuning in solution are elucidated. We show that a consistent estimate of the optimal range separated parameter (ω) at the limit of bulk solvation can be obtained by gradually extending the solvated system. In particular, ω is influenced by the solvent beyond the first coordination sphere. For determining ionization energies, a considerable number of solvent molecules on the first solvation shell must be taken into account. We demonstrate that accurately calculating optical gaps of solvated systems using OT-RSH can be done in three steps: (i) including the chemical environment when determining the range-separation parameter, (ii) taking into account the screening due to the solvent, and (iii) using realistic molecular geometries.

Bioinspired design of redox-active ligands for multielectron catalysis: effects of positioning pyrazine reservoirs on cobalt for electro- and photocatalytic generation of hydrogen from water

Mononuclear metalloenzymes in nature can function in cooperation with precisely positioned redox-active organic cofactors in order to carry out multielectron catalysis. Inspired by the finely tuned redox management of these bioinorganic systems, we present the design, synthesis, and experimental and theoretical characterization of a homologous series of cobalt complexes bearing redox-active pyrazines. These donor moieties are locked into key positions within a pentadentate ligand scaffold in order to evaluate the effects of positioning redox non-innocent ligands on hydrogen evolution catalysis. Both metal- and ligand-centered redox features are observed in organic as well as aqueous solutions over a range of pH values, and comparison with analogs bearing redox-inactive zinc(ii) allows for assignments of ligand-based redox events. Varying the geometric placement of redox non-innocent pyrazine donors on isostructural pentadentate ligand platforms results in marked effects on observed cobalt-catalyzed proton reduction activity. Electrocatalytic hydrogen evolution from weak acids in acetonitrile solution, under diffusion-limited conditions, reveals that the pyrazine donor of axial isomer 1-Co behaves as an unproductive electron sink, resulting in high overpotentials for proton reduction, whereas the equatorial pyrazine isomer complex 2-Co is significantly more active for hydrogen generation at lower voltages. Addition of a second equatorial pyrazine in complex 3-Co further minimizes overpotentials required for catalysis. The equatorial derivative 2-Co is also superior to its axial 1-Co congener for electrocatalytic and visible-light photocatalytic hydrogen generation in biologically relevant, neutral pH aqueous media. Density functional theory calculations (B3LYP-D2) indicate that the first reduction of catalyst isomers 1-Co, 2-Co, and 3-Co is largely metal-centered while the second reduction occurs at pyrazine. Taken together, the data establish that proper positioning of non-innocent pyrazine ligands on a single cobalt center is indeed critical for promoting efficient hydrogen catalysis in aqueous media, akin to optimally positioned redox-active cofactors in metalloenzymes. In a broader sense, these findings highlight the significance of electronic structure considerations in the design of effective electron-hole reservoirs for multielectron transformations.

Unraveling the Mechanism of Photoinduced Charge Transfer in Carotenoid-Porphyrin-C60 Molecular Triad

Photoinduced charge transfer (CT) plays a central role in biologically significant systems and in applications that harvest solar energy. We investigate the relationship of CT kinetics and conformation in a molecular triad. The triad, consisting of carotenoid, porphyrin, and fullerene is structurally flexible and able to acquire significantly varied conformations under ambient conditions. With an integrated approach of quantum calculations and molecular dynamics simulations, we compute the rate of CT at two distinctive conformations. The linearly extended conformation, in which the donor (carotenoid) and the acceptor (fullerene) are separated by nearly 50 Å, enables charge separation through a sequential CT process. A representative bent conformation that is entropically dominant, however, attenuates the CT, although the donor and the acceptor are spatially closer. Our computed rate of CT at the linear conformation is in good agreement with measured values. Our work provides unique fundamental understanding of the photoinduced CT process in the molecular triad.

Molecular structure, spectroscopy, and photoinduced kinetics in trinuclear cyanide bridged complex in solution: a first-principles perspective

We investigate the molecular structure of the solvated complex, [(NC)6Fe-Pt(NH3)4-Fe(CN)6](4-), and related dinuclear and mononuclear model complexes using first-principles calculations. Mixed nuclear complexes in both solution and crystal phases were widely studied as models for charge transfer (CT) reactions using advanced spectroscopical and electrochemical tools. In contrast to earlier interpretations, we find that the most stable gas phase and solvated geometries are substantially different from the crystal phase geometry, mainly due to variance in the underlying oxidation numbers of the metal centers. Specifically, in the crystal phase a Pt(IV) metal center resulting from Fe ← Pt backward electron transfers is stabilized by an octahedral ligand field, whereas in the solution phase a Pt(II) metal complex that prefers a square planar ligand field forms a CT salt by bridging to the iron complexes through long-range electrostatic interactions. The different geometry is shown to be consistent with spectroscopical data and measured CT rates of the solvated complex. Interestingly, we find that the experimentally indicated photoinduced process in the solvated complex is of backward CT (Fe ← Pt).

Quantum chemical study of the enzymatic repair of T(6-4)C/C(6-4)T UV-photolesions by DNA photolyases

Several strategies have evolved to repair one of the abundant UV radiation-induced damages caused to DNA, namely the mutagenic pyrimidine (6-4) pyrimidone photolesions. DNA (6-4)-photolyases are enzymes repairing these lesions by a photoinitiated electron transfer. An important aspect of a possible repair mechanism is its generality and transferability to different (6-4) lesions. Therefore, previously suggested mechanisms for the repair of the T(6-4)T lesion are here transferred to the T(6-4)C and C(6-4)T lesions and investigated theoretically using quantum chemical methods. Despite the different functional groups of the pyrimidine bases involved, a general valid molecular mechanism was identified, in which the initial step is an electron transfer coupled to a proton transfer from the protonated HIS365 to the N3(') nitrogen of the 3(') pyrimidine, followed by an intramolecular OH/NH2 transfer in one concerted step, which does not require an oxetane/azetidine or isolated water/ammonia intermediate.