Search

The excited state proton transfer (ESPT) reaction from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid (HPTS or pyranine) to an acetate molecule has been investigated in explicit aqueous solution via excited state ab initio molecular dynamics simulations based on hybrid quantum/molecular mechanics (QM/MM) potentials. In all the trajectories, the direct proton transfer has been observed in the excited state within 1 ps. We find that the initial structural configuration extracted from the ground state distribution strongly affects the ESPT kinetics. Indeed, the relative orientation of the proton donor-acceptor pair and the presence of a water molecule hydrogen bonded to the phenolic acid group of the pyranine are the key factors to facilitate the ESPT. Furthermore, we analyze the vibrational fingerprints of the ESPT reaction, reproducing the blue shift of the acetate CO stretching (COac), from 1666 to 1763 cm-1 testifying the transformation of acetate to acetic acid. Finally, our findings suggest that the acetate CC stretching (CCac) is also sensitive to the progress of the ESPT reaction. The CCac stretching is indeed ruled by the two vibrational modes (928 and 1426 cm-1), that in the excited state are alternately activated when the proton is shared or bound to the donor/acceptor, respectively.

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(ii) complex, [Ru(dcbpy)2(NCS)2]4- [dcbpy = (4,4′-dicarboxy-2,2′-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C2 symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies. This journal is

In the search for new therapeutic strategies to contrast SARS-CoV-2, we here studied the interaction of polydatin (PD) and resveratrol (RESV)—two natural stilbene polyphenols with manifold, well known biological activities—with Spike, the viral protein essential for virus entry into host cells, and ACE2, the angiotensin-converting enzyme present on the surface of multiple cell types (including respiratory epithelial cells) which is the main host receptor for Spike binding. Molecular Docking simulations evidenced that both compounds can bind Spike, ACE2 and the ACE2:Spike complex with good affinity, although the interaction of PD appears stronger than that of RESV on all the investigated targets. Preliminary biochemical assays revealed a significant inhibitory activity of the ACE2:Spike recognition with a dose-response effect only in the case of PD.

We present electronic structure methods to unveil the non-radiative pathways of photoinduced charge transfer (CT) reactions that play a main role in photophysics and light harvesting technologies. A prototypical π-stacked molecular complex consisting of an electron donor (1-chloronaphthalene, 1ClN) and an electron acceptor (tetracyanoethylene, TCNE) was investigated in dichloromethane solution for this purpose. The characterization of TCNE:π:1ClN in both its equilibrium ground and photoinduced low-lying CT electronic states was performed by using a reliable and accurate theoretical–computational methodology exploiting ab initio molecular dynamics simulations. The structural and vibrational time evolution of key vibrational modes is found to be in excellent agreement with femtosecond stimulated Raman spectroscopy experiments [R. A. Mathies et al., J. Phys. Chem. A, 2018, 122, 14, 3594], unveiling a correlation between vibrational fingerprints and electronic properties. The evaluation of nonadiabatic coupling matrix elements along generalized normal modes has made possible the interpretation on the molecular scale of the activation of nonradiative relaxation pathways towards the ground electronic state. In particular, two low frequency vibrational modes such as the out of plane bending and dimer breathing and the TCNE central C Created by potrace 1.16, written by Peter Selinger 2001-2019 C stretching play a prominent role in relaxation phenomena from the electronic CT state to the ground state one., We present electronic structure methods to unveil the non-radiative pathways of photoinduced charge transfer (CT) reactions that play a main role in photophysics and light harvesting technologies.

In this work, we simulate the excited state proton transfer (ESPT) reaction involving the pyranine photoacid and an acetate molecule as proton acceptor, connected by a bridge water molecule. We employ ab initio molecular dynamics combined with an hybrid quantum/molecular mechanics (QM/MM) framework. Furthermore, a time-resolved vibrational analysis based on the wavelet-transform allows one to identify two low frequency vibrational modes that are fingerprints of the ESPT event: a ring wagging and ring breathing. Their composition suggests their key role in optimizing the structure of the proton donor-acceptor couple and promoting the ESPT event. We find that the choice of the QM/MM partition dramatically affects the photoinduced reactivity of the system. The QM subspace was gradually extended including the water molecules directly interacting with the pyranine-water-acetate system. Indeed, the ESPT reaction takes place when the hydrogen bond network around the reactive system is taken into account at full QM level.

We present a novel time-resolved vibrational analysis for studying photoinduced nuclear relaxation. Generalized modes velocities are defined from ab initio molecular dynamics and wavelet transformed, providing the time localization of vibrational signals in the electronic excited state. The photoexcited pyranine in aqueous solution is presented as a case study. The transient and sequential activation of the simulated vibrational signals is in good agreement with vibrational dynamics obtained from femtosecond stimulated Raman spectroscopy data.

Advances in time resolved vibrational spectroscopy techniques provided a new stimulus for understanding the transient molecular dynamics triggered by the electronic excitation. The detailed interpretation of such time dependent spectroscopic signals is a challenging task from both experimental and theoretical points of view. We simulated and analyzed the transient photorelaxation of the pyranine photoacid in aqueous solution, with special focus on structural parameters and low frequency skeleton modes that are possibly preparatory for the photoreaction occurring at later time, as suggested by experimental spectroscopic studies. To this aim, we adopted an accurate computational protocol that combines excited state ab-initio molecular dynamics within a hybrid quantum mechanical/molecular mechanics framework and a time resolved vibrational analysis based on the wavelet transform. According to our results, the main nuclear relaxation on the excited potential energy surface is completed in about 500 fs, in agreement with experimental data. The rearrangement of C-C bonds occurs according to a complex vibrational dynamics, showing oscillatory patterns that are out of phase and modulated by modes below 200 cm-1. We also analyzed in both the ground and the excited state the evolution of some structural parameters involved in excited state proton transfer reaction, namely those involving the pyranine and the water molecule hydrogen bonded to the phenolic O-H group. Both the hydrogen bond distance and the intermolecular orientation are optimized in the excited state, resulting in a tighter proton donor-acceptor couple. Indeed, we found evidence that collective low frequency skeleton modes, such as the out of plane wagging at 108 cm-1and the deformation at 280 cm-1are photoactivated since the ultrafast part of the relaxation and modulate the pyranine-water molecule rearrangement, favouring the step preparatory for the photoreactivity.

Synthetic receptors for biomacromolecules lack the supramolecular self-assembly behavior typical of biological systems. Here we propose a new method for the preparation of protein imprinted polymers based on the specific interaction of a peptide multi-functional block with a protein target. This peptide block contains a protein-binding peptide domain, a polymerizable moiety at the C-terminus and an environment-sensitive fluorescent molecule at the N-terminus. The method relies on a preliminary step consisting of peptide/protein supramolecular assembly, followed by copolymerization with the most common acrylate monomers (acrylamide, acrylic acid and bis-acrylamide) to produce a protein imprinted hydrogel polymer. Such a peptide block can function as an active assistant recognition element to improve affinity, and guarantees its effective polymerization at the protein/cavity interface, allowing for proper placement of a dye. As a proof of concept, we chose Bovine Serum Albumin (BSA) as the protein target and built the peptide block around a BSA binding dodecapeptide, with an allyl group as the polymerizable moiety and a dansyl molecule as the responsive dye. Compared to conventional approaches these hydrogels showed higher affinity (more than 45%) and imprinted sensitivity (about twenty fold) to the target, with a great BSA selectivity with respect to ovalbumin (alpha = 1.25) and lysozyme (alpha = 6.02). Upon protein binding, computational and experimental observations showed a blue shift of the emission peak (down to 440 nm) and an increase of fluorescence emission (twofold) and average lifetime (Delta(tau) = 4.3 ns). Such an approach generates recognition cavities with controlled chemical information and represents an a priori method for self-responsive materials. Provided a specific peptide and minimal optimization conditions are used, such a method could be easily implemented for any protein target.

Among polyphenols, trans-resveratrol (tRES) and trans-polydatin (tPD) exert multiple biological effects, particularly antioxidant and antiproliferative. In this work, we have investigated the interaction of tPD with three cancer-related DNA sequences able to form G-quadruplex (G4) structures, as well as with a model duplex, and compared its behaviour with tRES. Interestingly, fluorescence analysis evidenced the ability of tPD to bind all the studied DNA systems, similarly to tRES, with tRES displaying a higher ability to discriminate G4 over duplex with respect to tPD. However, neither tRES nor tPD produced significant conformational changes of the analyzed DNA upon binding, as determined by CD-titration analysis. Computational analysis and biological data confirmed the biophysical results: indeed, molecular docking evidenced the stronger interaction of tRES with the promoter of c-myc oncogene, and immunoblotting assays revealed a reduction of c-myc expression, more effective for tRES than tPD. Furthermore, in vitro assays on melanoma cells proved that tPD was able to significantly reduce telomerase activity, and inhibit cell proliferation, with tRES producing higher effects than tPD.

We simulated for the first time the oxygen-oxygen bond formation in a synthetic calcium-tetra manganese complex recently developed by Zhang and co-workers. In spite of promising structural similarities to the native oxygen evolving complex (OEC) in Photosystem II, several uncertainties on the mimic stability in water and on its potential catalytic activity still persist. Here, we characterized at density functional theory level the electronic and structural features of the Sn states of the complex, along with the oxygen-oxygen bond formation reaction, proposing a reasonable model for the hydrate complex. As a main finding, both the synthetic compound and the natural OEC show very close energetic barriers for the oxo-oxyl coupling process, suggesting that key electronic features of the natural OEC reactivity are well reproduced. This result strongly encourages the use of this synthetic complex in combination with other molecular assemblies for the design of successful artificial leaves.

We simulated the intrinsic reaction path of the Green Fluorescent Protein (GFP) proton shuttle in both the ground state (S0) and first singlet excited state (S1), accounting for the main energetic and steric effects of the protein in a convenient model including the chromophore, the crystallographic water, and the residues directly involved in the proton transfer event. We adopted density functional theory (DFT) and time-dependent density functional theory (TD-DFT) levels to define the potential energy surfaces of the two electronic states, and we compared results obtained by the Damped Velocity Verlet and the Hessian-based Predictor-Corrector integrators of the intrinsic reaction coordinate, which gave a comparable and consistent picture of the mechanism. We show that, at S1, the GFP proton transfer becomes favored, with respect to S0, as suggested by the experimental evidence. As an important finding, this change is strictly related to the rearrangement of the hydrogen bond network composing the reaction path, which, in S1, relaxes to a tighter and planar configuration, as a consequence of the photoinduced relaxation in the GFP chromophore structure, thus prompting more effectively for the proton shuttle. Therefore, we give an unprecedented direct proof of the key role played by the photoinduced structural relaxation of the GFP on the chromophore photoacidity, validating, in particular, the hypothesis of Fang and co-workers [Nature 2009, 462, 200].

The growing production of biodiesel as a promising alternative and renewable fuel led as the main problem the dramatic increase of its by-product: glycerol. Different strategies for glycerol derivatization have been reported so far, some more efficient or sustainable than others. Herein, we report a very promising and eco-friendly transformation of glycerol in nontoxic solvents and chemicals (i.e., solketal, ketals), proposing three new families of Fe(III) compounds capable of catalysing glycerol acetalization with unpublished turn over frequencies (TOFs), and adhering most of the principles of green chemistry. The comparison between the activity of complexes of formula [FeCl3(1-R)] (1-R = substituted pyridinimine), [FeCl(2-R,R′)] (2-R,R′ = substituted O,O′-deprotonated salens) and their corresponding simple salts reveals that the former are extremely convenient because they are able to promote solketal formation with excellent TOFs, up to 105 h–1. Satisfactory performances were shown with respect to the entire range of substrates, with results being competitive to those reported in the literature so far. Moreover, the experimental activity was supported by an accurate and complete ab initio study, which disclosed the fundamental role of iron(III) as Lewis acid in promoting the catalytic activity. The unprecedented high activity and the low loading of the catalyst, combined with the great availability and the good eco-toxicological profile of iron, foster future applications of this catalytic process for the sustainable transformation of an abundant by-product in a variety of chemicals.

The detailed knowledge of excited state proton transfer mechanisms in complex environments is of paramount importance in chemistry. However, the definition of an effective reaction coordinate and the understanding of the driving force of the reaction can be difficult from both the experimental and the theoretical points of view. Here we analyzed by theoretical approaches the mechanism and the driving forces of the excited state proton transfer reaction occurring between the 7-hydroxy-4-(trifluoromethyl)coumarin photoacid and the 1-methylimidazole base molecules in toluene solution. In particular, we compared the intrinsic and the dynamical reaction pathways, obtained by integrating the reaction coordinate, and by performing ab initio simulations of molecular dynamics, respectively. Time-dependent density functional theory and polarizable solvation continuum models were adopted to define the excited state potential energy surface. Results were analyzed by means of the D(CT) electronic density based index. Our findings suggest that the reaction coordinate is mainly composed of several intra- and intermolecular modes of the reactants. An analysis of both the intrinsic coordinate and the dynamical results shows that the charge transfer induced by electronic excitation of the coumarin molecule is the main proton transfer driving force. With regards to the methodological validation, the combination of ab initio molecular dynamics with time-dependent density functional theory appears to be feasible and reliable to study excited state proton transfer reactions in the condensed phase.

A new time-resolved vibrational analysis unveils the mechanism of an excited state proton shuttle in green fluorescent protein., We simulated an excited state proton transfer in green fluorescent protein by excited state ab initio dynamics, and examined the reaction mechanism in both the time and the frequency domain through a multi resolution wavelet analysis. This original approach allowed us, for the first time, to directly compare the trends of photoactivated vibrations to femtosecond stimulated Raman spectroscopy results, and to give an unequivocal interpretation of the role played by low frequency modes in promoting the reaction. We could attribute the main driving force of the reaction to an important photoinduced softening of the ring–ring orientational motion of the chromophore, thus permitting the tightening of the hydrogen bond network and the opening of the reaction pathway. We also found that both the chromophore (in terms of its inter-ring dihedral angle and phenolic C–O and imidazolinone C–N bond distances) and its pocket (in terms of the inter-molecular oxygen’s dihedral angle of the chromophore pocket) relaxations are modulated by low frequency (about 120 cm–1) modes involving the oxygen atoms of the network. This is in agreement with the femtosecond Raman spectroscopy findings in the time-frequency domain. Moreover, the rate in proximity to the Franck Condon region involves a picosecond time scale, with a significant influence from fluctuations of nearby hydrogen bonded residues such as His148. This approach opens a new scenario with ab initio simulations as routinely used tools to understand photoreactivity and the results of advanced time resolved spectroscopy techniques.

The mechanism of base to base intermolecular proton shuttling occurring in the excited state proton transfer reaction between 7-hydroxy-4-(trifluoromethyl)coumarin (CouOH) and concentrated 1-methylimidazole base (1-MeId) in toluene solution is disclosed here by means of a computational approach based on Density Functional Theory (DFT) and Time Dependent DFT (TD-DFT). These methods allow us to characterize both the ground and excited state potential energy surfaces along the proton shuttling coordinate, and to assess the nature of the emitting species in the presence of an excess of 1-MeId. As a result, the tautomerism of CouOH is found to be photo-activated and, from a mechanistic point of view, the calculations clearly show that the overall driving force of the entire shuttling is the coumarin photoacidity, which is responsible for both the first proton transfer event and the strengthening of the following chain mechanism of base to base proton hopping.

In-chain donor/acceptor block copolymers comprised of alternating electron rich/poor moieties are emerging as promising semiconducting chromophores for use in organic photovoltaic devices. The mobilities of charge carriers in these materials are experimentally probed using gated organic field-effect transistors to quantify electron and hole mobilities, but a mechanistic understanding of the relevant charge diffusion pathways is lacking. To elucidate the mechanisms of electron and hole transport following excitation to optically accessible low-lying valence states, we utilize mean-field quantum electronic dynamics in the TDDFT formalism to explicitly track the evolution of these photo-accessible states. From the orbital pathway traversed in the dynamics, p- and n-type conductivities can be distinguished. The electronic dynamics of the studied polymers show the time-resolved transitions between the initial photoexcited state, a tightly-bound excitonic state that is dark to the ground state, and a partially charge separated state indicated by long-lived, out-of-phase charge oscillations along the polymer backbone. The frequency of these charge oscillations yields an insight into the characteristic mobilities of charge carriers in these materials. When the barycenters of the electron and hole densities are followed during the dynamics, a pseudo-classical picture for the translation of charge carrier densities along the polymer backbone emerges that clarifies a crucial aspect in the design of efficient organic photovoltaic materials.

The excited state proton transfer (ESPT) reaction taking place between 7-hydroxy-4-(trifluorometyl)coumarin and 1-methylimidazole, recently experimentally characterized, has been here considered as a case study to illustrate the possibility of using theoretical approaches rooted in density functional theory (DFT) and time-dependent DFT (TD-DFT) for the description of complex reactions occurring at the excited state. In particular, beside identifying all stable species occurring at the ground and excited state during the ESPT reaction, a quantitative characterization of their photophysical properties, such as absorption and emission, is obtained by properly including solvent effects. More interestingly, a computational protocol enabling one to locate possible reaction pathways for the ESPT is here proposed. This protocol is based on the use of density based indices purposely developed to characterize the properties of vertical and relaxed excited states which allow one to discriminate the most favorable reaction paths on potential energy surfaces that are in the case of ESPT intrinsically very flat and difficult to characterize based on sole energy criteria, thus opening a new scenario for the description of photoinduced proton transfer reactions.

A new descriptor tool enabling to qualitatively identify excited-state potential energy regions with high decay probability is here disclosed and applied to analyze a photoinduced ring-opening reaction already well characterized from the experimental and theoretical point of view. The analysis based on such descriptor allows one to highlight a high probability of excited-state deactivation within the Franck–Condon region, in agreement with experiments, and to qualitatively indentify the main mechanisms providing efficient pathways of photoreactivity.

The photophysical signature of the tautomeric species of the asymmetric (N,N-dimethylanilino)-1,3-diketone molecule are investigated using approaches rooted in density functional theory (DFT) and time-dependent DFT (TD-DFT). In particular, since this molecule, in the excited state, can undergo proton transfer reactions coupled to intramolecular charge transfer events, the different radiative and nonradiative channels are investigated by making use of different density-based indexes. The use of these tools, together with the analysis of both singlet and triplet potential energy surfaces, provide new insights into excited-state reactivity allowing one to rationalize the experimental findings including different behavior of the molecule as a function of solvent polarity.

In spite of the detailed information provided by advanced time-resolved spectroscopy, the understanding of the excited-state proton transfer (ESPT) reactivity remains difficult to obtain at molecular level. In this work we studied three photoacids showing different strength: the 8-hydroxy-1,3,6-pyrenetrisulfonate weak photoacid, the N-methyl-6-hydroxyquinolinium strong photoacid and the phenol-carboxyether dipicolinium cyanine (QCy9) superphotoacid, focusing on the intermolecular ESPT toward a solvent molecule or a base molecule in aqueous solution. To this aim, the ground and the first singlet excited-state potential energy surfaces of the three systems were characterized by means of the time-dependent density functional theory and a hybrid implicit/explicit model of the solvent. Main structural and photophysical features of the photoacids were assessed and satisfactorily compared with the experimental data. Energy profiles along the PT coordinate were analyzed in both the electronic states. We reproduced many important features of the photoacidity experimentally observed. The results suggest that the relative strength is mainly due to the different extent of charge transfer caused by the electronic transition in proximity of the acid group. Remarkably, we found that even in the case of the strongest photoacid (QCy9), showing a ESPT rate as rapid as to escape the solvent dynamics control, the PT is modulated and supported by the first solvation shell of the proton-accepting molecule. However, a complete understanding of this fascinating field needs the full account for the electronic and the molecular dynamics in play at different timescales.

In the present work, the general liquid optimized boundary (GLOB) model, which is a discrete/continuum approach for molecular simulations of liquids and solutions, has been extended to include a pressure coupling algorithm based on an extended phase-space scheme. To this end, a definition of the instantaneous pressure for a microscopic system has been derived from the minimum energy principle for the Helmholtz free energy. Applications to a pure liquid, such as water, and an aqueous solution of myoglobin are presented and the results are compared to those obtained using standard periodic boundary conditions (PBC).

In this work, the GLOB model, an effective and reliable computational approach well suited for ab initio and QM/MM molecular dynamics simulations of complex molecular systems in solution, has been applied to study two representative open-shell systems, the cobalt(II) ion and the glycine radical in aqueous solution, with special reference to their structural and magnetic properties. The main structural features of the solvent cage around the cobalt ion and the hydrogen bonding patterns around the neutral and zwitterionic forms of the glycine radical have been investigated in some detail. The general good agreement with experiments supports the use of the present model to investigate more challenging and biological/technological relevant open-shell systems.

Semi-empirical correction to density functional theory for dispersion (DFT-D) has been implemented for energies, analytical gradients, and Hessians in order to explore potential energy surfaces by means of a complete set of first-principle methods. The impact of non-bonding interactions on structures, binding energies and zero-point energy contributions as well as on ab initio molecular dynamics trajectories have been investigated for the well known case of benzene dimer. While the static results are in remarkable agreement with the most sophisticated post-Hartree–Fock approaches, the low cost of DFT-D allows to unravel dynamical aspects too, which are mandatory for situations ruled by weak interactions.

A recently developed extended Lagrangian model employing localized basis functions and nonperiodic boundary conditions (GLOB/ADMP) was applied to the radicals issuing from the homolytic breaking of the C(alpha)-H(alpha) bond of glycine in aqueous solution at different pH values. Although the modifications of the structure and the magnetic properties of these species induced by the solvent are qualitatively reproduced by a static discrete-continuum model, magnetic parameters are further tuned by short-time dynamical effects (solute vibrations and solvent librations). The results delivered by GLOB/ADMP simulations for both hyperfine tensors and g-tensors are in remarkable agreement with their experimental counterparts, allowing a reliable disentanglement of the overall observables into well-defined contributions. The dominant role of out-of-plane vibrations in determining hyperfine splittings is confirmed and quantified, together with the remarkable sensitivity of the gyromagnetic tensor to bond lengths and valence angles defining the NC(alpha)C' moiety. Together with their specific interest for the title radical, our results suggest some interesting trends for other biologically significant radicals and point out the need of extending magneto-structural relationships to dynamical aspects.

The recent advances in the development and application of density-based indexes for the description of the nature and the quantification of the extent of charge transfer associated with a given electronic transition are here reviewed. Starting from the basic definition of the indexes, a brief overview of their potential as indicators of potentially problematic cases in the description of charge transfer excitations using Time Dependent Density Functional Theory (TD-DFT) will be first given together with their possible application for comparing TD-DFT results to post Hartree–Fock (post-HF) calculations. After this methodological part, several examples of the application of density-based indexes to describe, from a quantitative and qualitative point of view, the charge transfer character (for instance in push–pull systems) or to map excited state reaction pathways (for instance in the case of Excited State Proton Transfer reactions) will be given to exemplify the insights that these indexes may bring to the description and design of new compounds of potential technological relevance.

The characterization of fundamental vibrational levels from the analysis of ab-initio dynamics is discussed. Results from test calculations are compared with those obtained by a second-order perturbative procedure to include anharmonicity exploiting the same ab-initio potential. The influence of the single molecular modes on vibrationally averaged properties calculated on the fly during the dynamics is also analyzed. Formaldehyde, both in the gas phase and in acetonitrile solution, and the vinyl radical in the gas phase have been chosen as test molecules, while the hyperfine coupling constants of the open shell system is computed as an example of molecular properties. Results allow for a quantitative comparison with experimental data.

This paper describes the extension of our computational strategy for pK predictions of small molecules to large solutes. The basic computational tool results from the coupling of quantum mechanical methods rooted in the density functional theory with the most recent version of the Polarizable Continuum Model. However, a third level is introduced, which includes solute regions far from the reactive center, which are described at a simplified level. This partition, together with the recent implementation of fast cavity generation, powerful iterative solvers, and fast multipole technology, allows us to tackle solutes of the dimension of a small protein. The problems and perspectives of this methodology are analyzed with special reference to the behavior of different Polarizable Continuum Model versions on the challenging playground represented by the pK’s of the different histidine residues occurring in the human prion protein.

In this paper, we report on a general approach for the detection of a specific tumoural biomarker directly in serum. Such detection is made possible using a protein-binding peptide selected through an improved phage display technique and then conjugated to engineered microparticles (MPs). Protein biomarkers represent an unlimited source of information for non-invasive diagnostic and prognostic tests; MP-based assays are becoming largely used in manipulation of soluble biomarkers, but their direct use in serum is hampered by the complex biomolecular environment. Our technique overcomes the current limitations as it produces a selective MP—engineered with an antifouling layer—that ‘captures’ the relevant protein staying impervious to the background. Our system succeeds in fishing-out the human tumour necrosis factor alpha directly in serum with a high selectivity degree. Our method could have great impact in soluble protein manipulation and detection for a wide variety of diagnostic applications.

The excited state intramolecular proton transfer (ESIPT) reaction taking place within 2-(2-hydroxyphenyl) benzoxazole (HBT) and two recently experimentally characterised napthalimide derivatives – known as N-1 and N-4 – has been investigated in order to identify and test a possible protocol for the description and complete mechanistic and electronic characterisation of the reaction at the excited state. This protocol is based on density functional theory (DFT), time-dependent density functional theory (TD-DFT) and a recently proposed electron density based index (DCT). This method is able to identify all stable species involved in the reaction; discriminate between possible reaction pathways over potential energy surfaces (PES), which are intrinsically very flat and difficult to characterise, and quantitatively measure the excited state charge transfer character throughout the reaction. The photophysical properties of the molecules (i.e. absorption and emission wavelength) are also quantitatively determined via the implicit inclusion of solvent effects in the case of toluene and, more polar, tetrahydrofuran (THF). The accuracy obtained with this protocol then opens up the possibility of the ab initio design of molecules exhibiting ESIPT for tailored applications such as highly selective molecular sensors.

We individuate a photoinduced electron transfer (PeT) as a quenching mechanism affecting rhodamine B photophysics in solvent. The PeT involves an electron transfer from the carboxylate group to the xanthene ring of rhodamine B. This is finely modulated by the subtle balance of coulombic and non-classical interactions between the carboxyphenyl and xanthene rings, also mediated by the solvent. We propose the use of an electronic density based index, the so called DCT index, as a new tool to assess and quantify the nature of the excited states involved in non-radiative decays near the region of their intersection. In the present case, this analysis allows us to gain insight on the interconversion process from the bright state to the dark state responsible for the quenching of rhodamine B fluorescence. Our findings encourage the use of density based indices to study the processes affecting excited state reactions that are characterized by a drastic change in the excitation nature, in order to rationalize the photophysical behavior of complex molecular systems.

The detailed interpretation of time-resolved spectroscopic signals in terms of the molecular rearrangement during a photoreaction or a photophysical event is one of the most important challenges of both experimental and theoretical chemistry. Here we simulate a time-resolved fluorescence spectrum of a dye in aqueous solution, the N-methyl-6-oxyquinolinium betaine, and analyze it in terms of far IR and THz frequency contributions, providing a direct connection to specific molecular motions. To obtain this result, we build up an innovative and general approach based on excited state ab-initio molecular dynamics and a wavelet-based time-dependent frequency analysis of nonstationary signals. We obtain a nice agreement with key parameters of the solvent dynamics, such as the total Stokes shift and the Stokes shift relaxation times. As an important finding, we observe a strong change of specific solute-solvent interactions upon the electronic excitation, with the migration of about 1.5 water molecules from the first solvation shell toward the bulk. In spite of this event, the Stokes shift dynamics is ruled by collective solvent motions in the THz and far IR, which guide and modulate the strong rearrangement of the dye microsolvation. By the relaxation of THz and IR contributions to the emission signal, we can follow and understand in detail the molecularity of the process. The protocol presented here is, in principle, transferable to other time-resolved spectroscopic techniques.

Photophysical signatures, namely absorption and emission energies, lifetime and quantum yields, have been computed through TD-DFT approaches and compared with experimental counterparts for the Rhodamine 110 dye in aqueous solution. Thanks to a new protocol of analysis, based on the use of very promising electronic based indices, it has been possible to investigate the interplay between Rhodamine 110 dye's structure, degree of charge transfer upon excitation, and fluorescence signatures. This combined analysis is very promising to support the understanding of charge transfer based mechanisms affecting dyes photophysics.

Solvent effects on the conformational preferences of homopolypeptides constituted by α-aminoisobutyric acid (polyAib) have been investigated coupling the polarizable continuum model (PCM) either to a quantum mechanical (PBE0/6-31G(d)) or to a molecular−mechanical (Amber) representation of the solute. The results allow a deeper insight into the influence of the solvent on the 310/α-helix equilibrium in polyAib and show that PCM/Amber is a useful tool to study the conformational preferences of large peptides in condensed phase. As a matter of fact, all the PCM/Amber results are close to their PCM/PBE0 counterparts, except for some underestimation of the absolute solvation energies in polar solvents.

Possible radical reaction products issuing from H-atom addition to cytosine have been characterized and analyzed by means of a comprehensive quantum mechanical approach including density functional computations (B3LYP), together with simulation of the solvent by the polarizable continuum model (PCM), and averaging of spectroscopic properties over the most important vibrational motions. The hyperfine couplings of the semirigid 5,6-dihydrocytos-6yl radical computed at the optimized geometry are in good agreement with their experimental counterparts. On the other hand, vibrational averaging is mandatory for obtaining an effectively planar structure for the 5,6-dihydrocytos-5yl radical with the consequent equivalence of beta-hydrogens. Finally, only proper consideration of environmental effects restores the agreement between computed and experimental couplings for the base anion protonated at N3.

An approximate method, recently proposed to include in continuum solvation models the effects of electronic charge lying outside the solute cavity, has been adapted and implemented in the framework of the polarizable continuum model (PCM). This formulation exploits all the features already developed for the other PCM versions; it provides molecular free energies, gradients and second derivatives with respect to nuclear coordinates. The performances of this method have been tested in comparison with other PCM versions, in particular, we examined the reliability of this technique to reproduce actual volume charge distribution effects, compared to traditional procedures based on Gauss’ Law.

We have implemented the polarizable continuum model for solute–solvent interactions into a plane-wave (PW) Car–Parrinello code. A comparison between PW and Gaussian-type orbitals (GTO) basis sets implementations, concerning solvation energies and dipole moments of a number of small systems, is reported. All the computations have been performed using the BPW91 exchange-correlation functional. GTO calculation has been performed using the 6-311G(d,p) basis set for geometry optimizations, and the aug-cc-pVTZ basis for energies and dipole moments. PW calculations were carried out on the Gaussian optimized geometries, up to an energy cutoff of 30 Ry. Excellent agreement between PW and GTO calculations is obtained, with average deviations of solute–solvent interaction energies lower than 0.5 kcal/mol.

The conformational behavior of β-alanine zwitterion in aqueous solution has been investigated by post-Hartree−Fock and density functional methods using the polarizable continuum model (PCM) for the description of the solvent. Our results show that the PCM is able to reproduce the delicate balance between intramolecular hydrogen bridges and preferential solvation of isolated charged groups in determining the overall conformational preferences without any need of explicit solvent molecules. From a methodological point of view, the Hartree−Fock method is not sufficiently reliable for quantitative predictions of conformer equilibria, whereas approaches including electron correlation (DFT, MP2, MP4) perform remarkably well when using medium size basis sets including diffuse functions. Vibrational averaging arising from the torsional motion around the central CC bond leads to vicinal NMR coupling constants (J) in remarkable agreement with experiments when using an extended Karplus-type equation for the torsiona...

We present a new implementation of the polarizable continuum model, which allows the computation of solvation energies and their gradients for large systems. This is obtained by a balanced two-term iterative procedure free from any N 3 computational bottleneck and from any N 2 storage requirement, N being the number of finite elements used to describe the cavity occupied by the solute. Comparison with the results obtained by the standard matrix-inversion algorithm for small and medium-sized systems shows that the new algorithm couples a fully satisfactory numerical accuracy with remarkable computational efficiency.

A recently developed quantum mechanical approach devoted to the study of unstable species in solution was applied to the radicals resulting from the homolytic breaking of the Cα−Hα bond of glycine in aqueous solution at different pH values.The computational protocol includes density functional calculations, simulation of the solvent by a mixed discrete-continuum approach, and averaging of spectroscopic properties over the most important vibrational motions. In vacuo computations provide reliable results for the zwitterionic form when using hybrid Hartree−Fock/density-functional methods with purposely tailored basis sets. Under the same conditions, disappointing results are obtained for the magnetic properties of neutral and, especially, anionic forms. Although the modifications of the structure and the magnetic properties of these species induced by the solvent are well reproduced by either a continuum model or a supermolecule approach, quantitative results can be obtained only by a mixed discrete-continu...

Although lifetimes and quantum yields of widely used fluorophores are often largely characterized, a systematic approach providing a rationale of their photophysical behavior on a quantitative basis is still a challenging goal. Here we combine methods rooted in the time-dependent density functional theory and fluorescence lifetime imaging microscopy to accurately determine and analyze fluorescence signatures (lifetime, quantum yield, and band peaks) of several commonly used rhodamine and pyronin dyes. We show that the radiative lifetime of rhodamines can be correlated to the charge transfer from the phenyl toward the xanthene moiety occurring upon the S(0) ← S(1) de-excitation, and to the xanthene/phenyl relative orientation assumed in the S(1) minimum structure, which in turn is variable upon the amino and the phenyl substituents. These findings encourage the synergy of experiment and theory as unique tool to design finely tuned fluorescent probes, such those conceived for modern optical sensors.

Most spectroscopic parameters are influenced by nuclear dynamics and by the chemical environment. However, proper inclusion of these effects still represents a challenge in computational spectroscopy studies. In many cases, a route coupling satisfactory accuracy with reasonable computational costs consists in the integration of DFT-based methods to compute spectroscopic parameters, with ab initio molecular dynamics simulations to sample from the classical phase space of the system. Here, we discuss the application of this approach in two case studies of remarkable practical interest, namely the simulation of the absorption spectrum of 9-methyladenine, an adenine nucleoside model; and the prediction of electron spin resonance parameters for nitroxyl radicals, the prototypical spin probes. In both cases, the accuracy of the results increases significantly when the subtle interplay of intra-molecular dynamics and solvent effects is introduced.

Herein, we propose a conceptually innovative approach to investigating reaction mechanisms. This study demonstrates the importance of considering explicitly the effects of large amplitude motions, aside from the intrinsic reaction coordinate, when tuning the free energy landscape of reaction pathways. We couple the path collective variables method with DFT-based enhanced sampling simulations to characterize the associative mechanism of the hydrolysis of the methyl phosphate dianion in solution. Importantly, energetics and mechanistic differences are observed when passing from the potential to the free energy surface.

In this work, we consider a new combination of vibrational analysis and normal-like mode decomposition of Debye-Waller factors of solvated ions entirely based on molecular dynamics data. Such a novel time-dependent analysis procedure provides a direct link between x-ray absorption fine structure parameters and normal mode contributions for an ion-solvent system. The potentialities of such a methodology rely on two fundamental aspects which distinguish it from already available tools. First, a general vibrational analysis that does not require any Gaussian or harmonic model for describing atomic fluctuations in liquids. Second, a very accurate sampling of the short range motions around the structural probe via the recently developed atom centered density matrix propagation/general liquid optimized boundary method. This novel molecular dynamics methodology is based on an integrated ab initio/classical potential using localized basis functions and nonperiodic boundary conditions. As a case study we have chosen the Zn(II) ion in aqueous solution. The consistency of our results and the observed good agreement with experiments show how the key support to advanced structural techniques from molecular dynamics can be further expanded and investigated.

Radical species of the nucleic acid bases (NABs) play a primary role in radiation damage and electron transfer processes of DNA. In the present work, a novel QM/MM model for simulating molecular systems in the liquid phase has been applied to study the microsolvation of the uracil anion radical in comparison with the parent molecule. Moreover, we have analyzed the vertical detachment energy as a function of structural intramolecular and intermolecular fluctuations in order to find molecular modes that may possibly promote the electron transfer process.

Application of a new integrated computational approach for two widely used nitroxide spin probes allows to show unequivocally that proper account of stereoelectronic, environmental, and dynamical effects leads to magnetic properties in quantitative agreement with experimental results without the need of any empirical parameter. Together with their specific interest, our results point out, in our opinion, the importance of developing and validating computational approaches able to switch on and off different effects, including environmental and dynamical ones, in order to evaluate their specific role in determining the overall experimental outcome.

ESR experimental data combined with quantum mechanical calculations and first-principle molecular dynamics simulations have allowed the unequivocal identification and characterization of the canonical form of the uracil anion radical as the prevalent tautomer in aqueous solution, contrary to both gas-phase results and recently reported results in solution based on thermodynamic analyses. The present study confirms the effectiveness of spectroscopic/theoretical investigations of short lived molecular species in different environments.

Interpretation of structural properties and dynamic behavior of molecules in solution is of fundamental importance to understand their stability, chemical reactivity, and catalytic action. While information can be gained, in principle, by a variety of spectroscopic techniques, the interpretation of the rich indirect information that can be inferred from the analysis of experimental spectra is seldom straightforward because of the subtle interplay of several different effects, whose specific role is not easy to separate and evaluate. In such a complex scenario, theoretical studies can be very helpful at two different levels: (i) supporting and complementing experimental results to determine the structure of the target molecule starting from its spectral properties; (ii) dissecting and evaluating the role of different effects in determining the observed spectroscopic properties. This is the reason why computational spectroscopy is rapidly evolving from a highly specialized research field into a versatile an...

A novel molecular dynamics methodology, which is based on an integrated ab initio /classical potential using localized basis functions and non-periodic boundary conditions, has been applied to study the microsolvation of the Zn(II) ion in aqueous solution. The metal ion, along with its first solvent shell, was treated at DFT B3LYP level, while the remaining bulk solvent was modeled with an effective TIP3P water model. The solvent structural arrangement around the ion has been analyzed and the characteristic Zn–O distance favourably compared to recent ND, XRD and EXAFS experiments, as well as to sophisticated cluster calculations.