Time-Resolved Vibrational Analysis of Excited State Ab Initio Molecular Dynamics to Understand Photorelaxation: The Case of the Pyranine Photoacid in Aqueous Solution.
Ontology highlight
ABSTRACT:
SUBMITTER: Chiariello MG
PROVIDER: S-EPMC8011922 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description: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.
| S-EPMC8279639 | biostudies-literature
Project description:We present an on-the-fly ab initio semiclassical study of vibrational energy levels of glycine, calculated by Fourier transform of the wavepacket correlation function. It is based on a multiple coherent states approach integrated with monodromy matrix regularization for chaotic dynamics. All four lowest-energy glycine conformers are investigated by means of single-trajectory semiclassical spectra obtained upon classical evolution of on-the-fly trajectories with harmonic zero-point energy. For the most stable conformer I, direct dynamics trajectories are also run for each vibrational mode with energy equal to the first harmonic excitation. An analysis of trajectories evolved up to 50 000 atomic time units demonstrates that, in this time span, conformers II and III can be considered as isolated species, while conformers I and IV show a pretty facile interconversion. Therefore, previous perturbative studies based on the assumption of isolated conformers are often reliable but might be not completely appropriate in the case of conformer IV and conformer I for which interconversion occurs promptly.
| S-EPMC5472367 | biostudies-literature
Project description: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.
| S-EPMC5890789 | biostudies-literature
Project description:The interplay between light absorption and the molecular environment has a central role in the observed photophysics of a wide range of photoinduced chemical and biological phenomena. The understanding of the interplay between vibrational and electronic transitions is the focus of this work, since it can provide a rationale to tune the optical properties of charge transfer (CT) materials used for technological applications. A clear description of these processes poses a nontrivial challenge from both the theoretical and experimental points of view, where the main issue is how to accurately describe and probe drastic changes in the electronic structure and the ultrafast molecular relaxation and dynamics. In this work we focused on the intermolecular CT reaction that occurs upon photon absorption in a π-stacked model system in dichloromethane solution, in which the 1-chloronaphthalene (1ClN) acts as the electron donor and tetracyanoethylene (TCNE) is the electron acceptor. Density functional theory calculations have been carried out to characterize both the ground-state properties and more importantly the low-lying CT electronic transition, and excellent agreement with recently available experimental results [Mathies, R. A.; et al. J. Phys. Chem. A 2018, 122 (14), 3594] was obtained. The minima of the ground state and first singlet excited state have been accurately characterized in terms of spatial arrangements and vibrational Raman frequencies, and the CT natures of the first two low-lying electronic transitions in the absorption spectra have been addressed and clarified too. Finally, by modeling the possible coordination sites of the TCNE electron acceptor with respect to monovalent ions (Na+, K+) in an implicit solution of acetonitrile, we find that TCNE can accommodate a counterion in two different arrangements, parallel and orthogonal to the C═C axis, leading to the formation of a contact ion pair. The nature of the counterion and its relative position entail structural modifications of the TCNE radical anion, mainly the central C═C and C≡N bonds, compared to the isolated case. An important red shift of the C═C stretching frequency was observed when the counterion is orthogonal to the double bond, to a greater extent for Na+. On the contrary, in the second case, where the counterion ion lies along the internuclear C═C axis, we find that K+ polarizes the electron density of the double bond more, resulting in a greater red shift than with Na+.
| S-EPMC9574931 | biostudies-literature
Project description:The majority of macromolecular crystal structures are determined using the method of molecular replacement, in which known related structures are rotated and translated to provide an initial atomic model for the new structure. A theoretical understanding of the signal-to-noise ratio in likelihood-based molecular replacement searches has been developed to account for the influence of model quality and completeness, as well as the resolution of the diffraction data. Here we show that, contrary to current belief, molecular replacement need not be restricted to the use of models comprising a substantial fraction of the unknown structure. Instead, likelihood-based methods allow a continuum of applications depending predictably on the quality of the model and the resolution of the data. Unexpectedly, our understanding of the signal-to-noise ratio in molecular replacement leads to the finding that, with data to sufficiently high resolution, fragments as small as single atoms of elements usually found in proteins can yield ab initio solutions of macromolecular structures, including some that elude traditional direct methods.
| S-EPMC5389281 | biostudies-literature
Project description:Crystal structure determination requires solving the phase problem. This can be accomplished using ab initio direct methods for small molecules and macromolecules at resolutions higher than 1.2 Å, whereas macromolecular structure determination at lower resolution requires either molecular replacement using a homologous structure or experimental phases using a derivative such as covalent labeling (for example selenomethionine or mercury derivatization) or heavy-atom soaking (for example iodide ions). Here, a case is presented in which crystals were obtained from a 30.8 kDa protein sample and yielded a 1.6 Å resolution data set with a unit cell that could accommodate approximately 8 kDa of protein. Thus, it was unclear what had been crystallized. Molecular replacement with pieces of homologous proteins and attempts at iodide ion soaking failed to yield a solution. The crystals could not be reproduced. Sequence-independent molecular replacement using the structures available in the Protein Data Bank also failed to yield a solution. Ultimately, ab initio structure solution proved successful using the program ARCIMBOLDO, which identified two α-helical elements and yielded interpretable maps. The structure was the C-terminal dimerization domain of the intended target from Mycobacterium smegmatis. This structure is presented as a user-friendly test case in which an unknown protein fragment could be determined using ARCIMBOLDO.
| S-EPMC6130419 | biostudies-literature
Project description:The fluorescent base guanine analog, 8-vinyl-deoxyguanosine (8vdG), is studied in solution using a combination of optical spectroscopies, notably femtosecond fluorescence upconversion and quantum chemical calculations, based on time-dependent density functional theory (TD-DFT) and including solvent effect by using a mixed discrete-continuum model. In all investigated solvents, the fluorescence is very long lived (3-4 ns), emanating from a stable excited state minimum with pronounced intramolecular charge-transfer character. The main non-radiative decay channel features a sizeable energy barrier and it is affected by the polarity and the H-bonding properties of the solvent. Calculations provide a picture of dynamical solvation effects fully consistent with the experimental results and show that the photophysical properties of 8vdG are modulated by the orientation of the vinyl group with respect to the purine ring, which in turn depends on the solvent. These findings may have importance for the understanding of the fluorescence properties of 8vdG when incorporated in a DNA helix.
| S-EPMC7071107 | biostudies-literature
Project description:We built a full-dimensional analytical potential energy surface of the ground electronic state of Li₂H from ca. 20,000 ab initio multi-reference configuration interaction calculations, including core⁻valence correlation effects. The surface is flexible enough to accurately describe the three dissociation channels: Li (2s ²S) + LiH (¹Σ⁺), Li₂ (¹Σg⁺) + H (1s ²S) and 2Li (2s ²S) + H (1s ²S). Using a local fit of this surface, we calculated pure (J = 0) vibrational states of Li₂H up to the barrier to linearity (ca. 3400 cm-1 above the global minimum) using a vibrational self-consistent field/virtual state configuration interaction method. We found 18 vibrational states below this barrier, with a maximum of 6 quanta in the bending mode, which indicates that Li₂H could be spectroscopically observable. Moreover, we show that some of these vibrational states are highly correlated already ca. 1000 cm-1 below the height of the barrier. We hope these calculations can help the assignment of experimental spectra. In addition, the first low-lying excited states of each B₁, B₂ and A₂ symmetry of Li₂H were characterized.
| S-EPMC6337435 | biostudies-literature
Project description:Ab initio molecular dynamics simulations (AIMD) are systematically performed to study the Vibrational Energy Redistribution (VER) in solid nitromethane (NM) by combining normal mode decomposition and short-time Fourier transform technique. After the selective excitations of all fourteen intramolecular vibrational modes above 400 cm-1, four three-dimensional (3D) excitation and detected vibrational spectra are obtained. The evolution of the kinetic energy proportion of all vibrations are also given and discussed quantitatively. These results show that, as the daughter modes, NO2 symmetric stretches, CH3 stretches and bends are usually excited quickly and relatively conspicuously compared with the other vibrations. Interestingly, we found that, although the stretching vibration of the CN bond which is a bridge between the methyl and nitro group can not respond immediately to the selective excitations, it always accumulates the vibrational energy slowly and steadily. Then, the underlying mechanisms are discussed based on the response of vibrational modes in both the time and frequency domain. As a result, we found that anharmonic transfers following symmetry rules which involve the couplings assisted by the overtones and rotations, as well as the transfers among the adjacent modes, play important roles in the VER of solid NM.
| S-EPMC8695419 | biostudies-literature
Project description:Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).
| S-EPMC3797183 | biostudies-literature