Production of Dihydrogen Using Ammonia Borane as Reagent and Pyrazole as Catalyst.

Theoretical chemistry (DLPNO-CCSD(T)/def2-TZVP//M06-2x/aug-cc-pVDZ) was used to design a system based on ammonia boranes catalyzed by pyrazoles with the aim of producing dihydrogen, nowadays of high interest as clean fuel. The reactivity of ammonia borane and cyclotriborazane were investigated, including catalytic activation through 1H-pyrazole, 4-methoxy-1H-pyrazole, and 4-nitro-1H-pyrazole. The results point toward a catalytic cycle by which, at the same time, ammonia borane can initially store and then, through catalysis, produce dihydrogen and amino borane. Subsequently, amino borane can trimerize to form cyclotriborazane that, in presence of the same catalyst, can also produce dihydrogen. This study proposes therefore a consistent progress in using environmentally sustainable (metal free) catalysts to efficiently extract dihydrogen from small B-N bonded molecules.

Binding of Nitriles and Isonitriles to V(III) and Mo(III) Complexes: Ligand vs Metal Controlled Mechanism.

The synthesis and structures of nitrile complexes of V(N[tBu]Ar)3, 2 (Ar = 3,5-Me2C6H3), are described. Thermochemical and kinetic data for their formation were determined by variable temperature Fourier transform infrared (FTIR), calorimetry, and stopped-flow techniques. The extent of back-bonding from metal to coordinated nitrile indicates that electron donation from the metal to the nitrile plays a less prominent role for 2 than for the related complex Mo(N[tBu]Ar)3, 1. Kinetic studies reveal similar rate constants for nitrile binding to 2, but the activation parameters depend critically on the nature of R in RCN. Activation enthalpies range from 2.9 to 7.2 kcal·mol-1, and activation entropies from -9 to -28 cal·mol-1·K-1 in an opposing manner. Density functional theory (DFT) calculations provide a plausible explanation supporting the formation of a π-stacking interaction between a pendant arene of the metal anilide of 2 and the arene substituent on the incoming nitrile in favorable cases. Data for ligand binding to 1 do not exhibit this range of activation parameters and are clustered in a small area centered at ΔH‡ = 5.0 kcal·mol-1 and ΔS‡ = -26 cal·mol-1·K-1. Computational studies are in agreement with the experimental data and indicate a stronger dependence on electronic factors associated with the change in spin state upon ligand binding to 1.

Revealing the Molecular Interactions between Human ACE2 and the Receptor Binding Domain of the SARS-CoV-2 Wild-Type, Alpha and Delta Variants.

After a sudden and first spread of the pandemic caused by the novel SARS-CoV-2 (Severe Acute Respiratory Syndrome-Coronavirus 2) wild-type strain, mutants have emerged which have been associated with increased infectivity, inducing surges in the contagions. The first of the so-called variants of concerns, was firstly isolated in the United Kingdom and later renamed Alpha variant. Afterwards, in the middle of 2021, a new variant appeared called Delta. The latter is characterized by the presence of point mutations in the Spike protein of SARS-CoV-2, especially in the Receptor Binding Domain (RBD). When in its active conformation, the RBD can interact with the human receptor Angiotensin-Converting Enzyme 2 (ACE2) to allow the entry of the virions into cells. In this contribution, by using extended all-atom molecular dynamic simulations, complemented with machine learning post-processing, we analyze the changes in the molecular interaction network induced by these different strains in comparison with the wild-type. On one hand, although relevant variations are evidenced, only limited changes in the global stability indicators and in the flexibility profiles have been observed. On the other hand, key differences were obtained by tracking hydrophilic and hydrophobic molecular interactions, concerning both positioning at the ACE2/RBD interface and formation/disruption dynamic behavior.

Interaction of a 1,3-Dicarbonyl Toxin with Ru(II)-Biimidazole Complexes for Luminescence Sensing: A Spectroscopic and Photochemical Experimental Study Rationalized by Time-Dependent Density Functional Theory Calculations.

A family of ruthenium(II) complexes containing one 2,2'-biimidazole (bim) ligand and two polypyridyl (NN) ligands has been prepared and their photophysical and photochemical features have been tested in the presence of tenuazonic acid (TeA), a widespread food and feed mycotoxin of current concern. While not tested in in vivo studies, TeA and other secondary metabolites of Alternaria fungi are suspected to exert adverse effects on the human health, so sensors and rapid analytical procedures are required. It is well-known that 1,3-dicarbonyl compounds such as TeA are relatively easy to deprotonate (the pKa of TeA is 3.5), yielding an enolate anion stabilized by resonance. The chelating and hydrogen-donor features of bim allow simultaneous binding to the metal core and to the target ß-diketonate delocalized anion. Such a binding induces changes in the blue absorption (40 nm bathochromic shift), red luminescence intensity (>75% quenching), and triplet lifetime (0.2 µs decrease) of the Ru(NN)2(bim)2+ luminophore. Moreover, we have computationally rationalized, by time-dependent density functional theory, the structure of the different adducts of Ru-bim complexes with TeA and the electronic nature of the spectral absorption bands and their change upon the addition of TeA.

Atomistic-Level Description of the Covalent Inhibition of SARS-CoV-2 Papain-like Protease.

Inhibition of the papain-like protease (PLpro) of SARS-CoV-2 has been demonstrated to be a successful target to prevent the spreading of the coronavirus in the infected body. In this regard, covalent inhibitors, such as the recently proposed VIR251 ligand, can irreversibly inactivate PLpro by forming a covalent bond with a specific residue of the catalytic site (Cys111), through a Michael addition reaction. An inhibition mechanism can therefore be proposed, including four steps: (i) ligand entry into the protease pocket; (ii) Cys111 deprotonation of the thiol group by a Brønsted-Lowry base; (iii) Cys111-S- addition to the ligand; and (iv) proton transfer from the protonated base to the covalently bound ligand. Evaluating the energetics and PLpro conformational changes at each of these steps could aid the design of more efficient and selective covalent inhibitors. For this aim, we have studied by means of MD simulations and QM/MM calculations the whole mechanism. Regarding the first step, we show that the inhibitor entry in the PLpro pocket is thermodynamically favorable only when considering the neutral Cys111, that is, prior to the Cys111 deprotonation. For the second step, MD simulations revealed that His272 would deprotonate Cys111 after overcoming an energy barrier of ca. 32 kcal/mol (at the QM/MM level), but implying a decrease of the inhibitor stability inside the protease pocket. This information points to a reversible Cys111 deprotonation, whose equilibrium is largely shifted toward the neutral Cys111 form. Although thermodynamically disfavored, if Cys111 is deprotonated in close proximity to the vinylic carbon of the ligand, then covalent binding takes place in an irreversible way (third step) to form the enolate intermediate. Finally, due to Cys111-S- negative charge redistribution over the bound ligand, proton transfer from the initially protonated His272 is favored, finally leading to an irreversibly modified Cys111 and a restored His272. These results elucidate the selectivity of Cys111 to enable formation of a covalent bond, even if a weak proton acceptor is available, as His272.

π-Bridge Substitution in DASAs: The Subtle Equilibrium between Photochemical Improvements and Thermal Control*.

Donor-acceptor Stenhouse adducts (DASAs) are playing an outstanding role as innovative and versatile photoswitches. Until now, all the efforts have been spent on modifying the donor and acceptor moieties to modulate the absorption energy and improve the cyclization and reversion kinetics. However, there is a strong dependence on specific structural modifications and a lack of predictive behavior, mostly owing to the complex photoswitching mechanism. Here, by means of a combined experimental and theoretical study, the effect of chemical modification of the π-bridge linking the donor and acceptor moieties is systematically explored, revealing the significant impact on the absorption, photocyclization, and relative stability of the open form. In particular, a position along the π-bridge is found to be the most suited to redshift the absorption while preserving the cyclization. However, thermal back-reaction to the initial isomer is blocked. These effects are explained in terms of an increased acceptor capability offered by the π-bridge substituent that can be modulated. This strategy opens the path toward derivatives with infra-red absorption and a potential anchoring point for further functionalization.

Microscopic interactions between ivermectin and key human and viral proteins involved in SARS-CoV-2 infection.

The identification of chemical compounds able to bind specific sites of the human/viral proteins involved in the SARS-CoV-2 infection cycle is a prerequisite to design effective antiviral drugs. Here we conduct a molecular dynamics study with the aim to assess the interactions of ivermectin, an antiparasitic drug with broad-spectrum antiviral activity, with the human Angiotensin-Converting Enzyme 2 (ACE2), the viral 3CLpro and PLpro proteases, and the viral SARS Unique Domain (SUD). The drug/target interactions have been characterized in silico by describing the nature of the non-covalent interactions found and by measuring the extent of their time duration along the MD simulation. Results reveal that the ACE2 protein and the ACE2/RBD aggregates form the most persistent interactions with ivermectin, while the binding with the remaining viral proteins is more limited and unspecific.

E/Z Molecular Photoswitches Activated by Two-Photon Absorption: Comparison between Different Families.

Nonlinear optical techniques as two-photon absorption (TPA) have raised relevant interest within the last years due to the capability to excite chromophores with photons of wavelength equal to only half of the corresponding one-photon absorption energy. At the same time, its probability being proportional to the square of the light source intensity, it allows a better spatial control of the light-induced phenomenon. Although a consistent number of experimental studies focus on increasing the TPA cross section, very few of them are devoted to the study of photochemical phenomena induced by TPA. Here, we show a design strategy to find suitable E/Z photoswitches that can be activated by TPA. A theoretical approach is followed to predict the TPA cross sections related to different excited states of various photoswitches' families, finally concluding that protonated Schiff-bases (retinal)-like photoswitches outperform compared to the others. The donor-acceptor substitution effect is therefore rationalized for the successful TPA activatable photoswitch, in order to maximize its properties, finally also forecasting a possible application in optogenetics. Some experimental measurements are also carried out to support our conclusions.

Molecular Basis of SARS-CoV-2 Infection and Rational Design of Potential Antiviral Agents: Modeling and Simulation Approaches.

The emergence in late 2019 of the coronavirus SARS-CoV-2 has resulted in the breakthrough of the COVID-19 pandemic that is presently affecting a growing number of countries. The development of the pandemic has also prompted an unprecedented effort of the scientific community to understand the molecular bases of the virus infection and to propose rational drug design strategies able to alleviate the serious COVID-19 morbidity. In this context, a strong synergy between the structural biophysics and molecular modeling and simulation communities has emerged, resolving at the atomistic level the crucial protein apparatus of the virus and revealing the dynamic aspects of key viral processes. In this Review, we focus on how in silico studies have contributed to the understanding of the SARS-CoV-2 infection mechanism and the proposal of novel and original agents to inhibit the viral key functioning. This Review deals with the SARS-CoV-2 spike protein, including the mode of action that this structural protein uses to entry human cells, as well as with nonstructural viral proteins, focusing the attention on the most studied proteases and also proposing alternative mechanisms involving some of its domains, such as the SARS unique domain. We demonstrate that molecular modeling and simulation represent an effective approach to gather information on key biological processes and thus guide rational molecular design strategies.

The role of CO2 detachment in fungal bioluminescence: thermally vs. excited state induced pathways.

Different fungi lineages are known to emit light on Earth, mainly in tropical climates. Although the preparation of bioluminescent cell-free extracts allowed one to characterize the enzymatic requirements, the molecular mechanism underlying luminescence is still largely unknown and is based on the experimental putative assumption that a high-energy intermediate should be formed by reaction with O2 and formation of an endoperoxide. Here, we aim at determining, through state-of-the-art multiconfigurational quantum chemistry, the full mechanistic landscape leading from the endoperoxide to the emitting species, envisaging different possible pathways and proposing their viability. Especially, thermal CO2 detachment followed by excited-state peroxide opening (thermal-chemiluminescence) can compete with a parallel pathway, i.e., first excited-state endoperoxide opening, followed by CO2 detachment on the same excited-state (excited state-chemiluminescence). Clear differences in the energy supplies, as well as the possibility to directly populate the emitting species from the intersection seam between ground and excited states, land credence to a kinetically efficient thermal-chemiluminescent pathway, establishing for the first time a detailed description of fungal bioluminescence.

Trans-to-cis photoisomerization of cyclocurcumin in different environments rationalized by computational photochemistry.

Cyclocurcumin is a turmeric component that has attracted much less attention compared to the well-known curcumin. In spite of the less deep characterization of its properties, cyclocurcumin has shown promising anticancer effects when used in combination with curcumin. Especially, due to its peculiar molecular structure, cyclocurcumin can be regarded as an almost ideal photoswitch, whose capabilities can also be exploited for relevant biological applications. Here, by means of state-of-the-art computational methods for electronic excited-state calculations (TD-DFT, MS-CASPT2, and XMS-CASPT2), we analyze in detail the absorption and photoisomerization pathways leading from the more stable trans isomer to the cis one. The different molecular surroundings, taken into account by means of the electrostatic solvent effect and compared with available experimental data, have been found to be critical in describing the fate of irradiated cyclocurcumin: when in non-polar environments, an excited state barrier prevents photoisomerization and favours fluorescence, whereas in polar solvents, an almost barrierless path results in a striking decrease of fluorescence, opening the way toward a crossing region with the ground state and thus funneling the photoproduction of the cis isomer.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases.

An extensive theoretical characterization of the singlet excited state manifold of the five canonical DNA/RNA nucleobases (thymine, cytosine, uracil, adenine and guanine) in gas-phase is carried out with time-dependent density functional theory (TD-DFT) and restricted active space second-order perturbation theory (RASPT2) approaches. Both ground state and excited state absorptions are analyzed and compared between these different theoretical approaches, assessing the performance of the hybrid B3LYP and CAM-B3LYP (long-range corrected) functionals with respect to the RASPT2 reference. By comparing the TD-DFT estimates with our reference for high-lying excited states, we are able to narrow down specific energetic windows where TD-DFT may be safely employed to qualitatively reproduce the excited state absorption (ESA) signals registered in non-linear and time-resolved spectroscopy for monitoring photoinduced phenomena. Our results show a qualitative agreement between the RASPT2 reference and the B3LYP computed ESAs of pyrimidines in the near-IR/Visible spectral probing window while for purines the agreement is limited to the near-IR ESAs, with generally larger discrepancies obtained with the CAM-B3LYP functional. This outcome paves the way for appropriate application of cost-effective TD-DFT approaches to simulate linear and non-linear spectroscopies of realistic multichromophoric DNA/RNA systems with biological and nanotechnological relevance.

Competing ultrafast photoinduced electron transfer and intersystem crossing of [Re(CO) 3 (Dmp)(His124)(Trp122)] + in Pseudomonas aeruginosa azurin: a nonadiabatic dynamics study.

We present a computational study of sub-picosecond nonadiabatic dynamics in a rhenium complex coupled electronically to a tryptophan (Trp) side chain of Pseudomonas aeruginosa azurin, a prototypical protein used in the study of electron transfer in proteins. To gain a comprehensive understanding of the photoinduced processes in this system, we have carried out vertical excitation calculations at the TDDFT level of theory as well as nonadiabatic dynamics simulations using the surface hopping including arbitrary couplings (SHARC) method coupled to potential energy surfaces represented with a linear vibronic coupling model. The results show that the initial photoexcitation populates both singlet metal-to-ligand charge transfer (MLCT) and singlet charge-separated (CS) states, where in the latter an electron was transferred from the Trp amino acid to the complex. Subsequently, a complex mechanism of simultaneous intersystem crossing and electron transfer leads to the sub-picosecond population of triplet MLCT and triplet CS states. These results confirm the assignment of the sub-ps time constants of previous experimental studies and constitute the first computational evidence for the ultrafast formation of the charge-separated states in Re-sensitized azurin.

Charge-Transfer versus Charge-Separated Triplet Excited States of [ReI (dmp)(CO)3 (His124)(Trp122)]+ in Water and in Modified Pseudomonas aeruginosa Azurin Protein.

A computational investigation of the triplet excited states of a rhenium complex electronically coupled with a tryptophan side chain and bound to an azurin protein is presented. In particular, by using high-level molecular modeling, evidence is provided for how the electronic properties of the excited-state manifolds strongly depend on coupling with the environment. Indeed, only upon explicitly taking into account the protein environment can two stable triplet states of metal-to-ligand charge transfer or charge-separated nature be recovered. In addition, it is also demonstrated how the rhenium complex plus tryptophan system in an aqueous environment experiences too much flexibility, which prevents the two chromophores from being electronically coupled. This occurrence disables the formation of a charge-separated state. The successful strategy requires a multiscale approach of combining molecular dynamics and quantum chemistry. In this context, the strategy used to parameterize the force fields for the electronic triplet states of the metal complex is also presented.

Correction: Different hydrogen bonding environments of the retinal protonated Schiff base control the photoisomerization in channelrhodopsin-2.

Correction for 'Different hydrogen bonding environments of the retinal protonated Schiff base control the photoisomerization in channelrhodopsin-2' by Yanan Guo et al., Phys. Chem. Chem. Phys., 2018, 20, 27501-27509.

The effect of solvent relaxation in the ultrafast time-resolved spectroscopy of solvated benzophenone.

Benzophenone (BP) despite its relatively simple molecular structure is a paradigmatic sensitizer, featuring both photocatalytic and photobiological effects due to its rather complex photophysical properties. In this contribution we report an original theoretical approach to model realistic, ultra-fast spectroscopy data, which requires describing intra- and intermolecular energy and structural relaxation. In particular we explicitly simulate time-resolved pump-probe spectra using a combination of state-of-the art hybrid quantum mechanics/molecular mechanics dynamics to treat relaxation and vibrational effects. The comparison with experimental transient absorption data demonstrates the efficiency and accuracy of our approach. Furthermore the explicit inclusion of the solvent, water for simulation and methanol for experiment, allows us, despite the inherent different behavior of the two, to underline the role played by the H-bonding relaxation in the first hundreds of femtoseconds after optical excitation. Finally we predict for the first time the two-dimensional electronic spectrum (2DES) of BP taking into account the vibrational effects and hence modelling partially symmetric and asymmetric ultrafast broadening.

Different hydrogen bonding environments of the retinal protonated Schiff base control the photoisomerization in channelrhodopsin-2.

The first event of the channelrhodopsin-2 (ChR2) photocycle, i.e. trans-to-cis photoisomerization, is studied by means of quantum mechanics/molecular mechanics, taking into account the flexible retinal environment in the ground state. By treating the chromophore at the ab initio multiconfigurational level of theory, we can rationalize the experimental findings based on pump-probe spectroscopy, explaining the different and more complex scenario found for ChR2 in comparison to other rhodopsins. In particular, we find that depending on the hydrogen bonding pattern, different excited states are involved, hence making it possible to suggest one pattern as the most productive. Moreover, after photoisomerization the structure of the first photocycle intermediate, P5001, is characterized by simulating the infrared spectrum and compared to available experimental data. This was obtained by extensive molecular dynamics, where the chromophore is described by a semi-empirical method based on density functional theory. The results clearly identify which counterion is responsible for accepting the proton from the retinal Schiff base: the side chain of the glutamic acid E123.

Deciphering the photosensitization mechanisms of hypericin towards biological membranes.

Resorting to state-of-the art molecular modeling and simulation techniques we provide full characterization of the photophysical properties of the naturally occurring hypericin chromophore, currently used in photodynamic therapy. In particular, we reveal the different photophysical pathways leading to intersystem-crossing and hence, triplet manifold population that is necessary for the subsequent production of singlet oxygen. In particular we identify an extended region of quasi-degeneracy between the first singlet excited state and three triplet state surfaces. This energetic factor allows the occurrence of intersystem-crossing even in the presence of a relatively small spin-orbit coupling. Furthermore, thanks to extended all-atom molecular dynamics simulations we provide insight into the interaction of hypericin with lipid bilayers. We demonstrate the formation of stable interactions with the membrane and, in particular, the penetration of hypericin into its hydrophobic core. This organization allows a spatial overlap between hypericin and the lipid oxidizable double bond pointing towards the production of singlet oxygen in close spatial proximity to its reactant, hence favoring photosensitization.

Ab initio molecular dynamics of thiophene: the interplay of internal conversion and intersystem crossing.

The fast and slow components of the relaxation of photoexcited thiophene have been investigated by means of SHARC (surface hopping including arbitrary couplings) molecular dynamics based on multiconfiguration electronic structure calculations. Triplet states are included to ascertain their role in the relaxation process. After thiophene is excited to the S1 state, ultrafast dynamics (τfast = 96 fs) initiates a ring opening due to cleavage of a carbon sulfur bond and simultaneous ring puckering. This time constant is in agreement with previous experimental and theoretical studies. The subsequent dynamics of the open-ring structures is characterized by the interplay of internal conversion and intersystem crossing. For the open-ring structures, the S0, S1, T1 and T2 states are nearly degenerate and the spin-orbit couplings are large. The underlying potential energy surface is flat and long-lived open-ring structures in the singlet as well as in the triplet states are formed. Both the participation of triplet states and the shape of the energy surface explain the experimentally observed slow ring closure in the ground state.

Optomechanical Control of Quantum Yield in Trans-Cis Ultrafast Photoisomerization of a Retinal Chromophore Model.

The quantum yield of a photochemical reaction is one of the most fundamental quantities in photochemistry, as it measures the efficiency of the transduction of light energy into chemical energy. Nature has evolved photoreceptors in which the reactivity of a chromophore is enhanced by its molecular environment to achieve high quantum yields. The retinal chromophore sterically constrained inside rhodopsin proteins represents an outstanding example of such a control. In a more general framework, mechanical forces acting on a molecular system can strongly modify its reactivity. Herein, we show that the exertion of tensile forces on a simplified retinal chromophore model provokes a substantial and regular increase in the trans-to-cis photoisomerization quantum yield in a counterintuitive way, as these extension forces facilitate the formation of the more compressed cis photoisomer. A rationale for the mechanochemical effect on this photoisomerization mechanism is also proposed.