The Green Fluorescent Protein (GFP) is a widely studied chemical system both for its large amount of applications and the complexity of the excited state proton transfer responsible of the change in the protonation state of the chromophore. A detailed investigation on the structure of the chromophore environment and the influence of chromophore form (either neutral or anionic) on it is of crucial importance to understand how these factors could potentially influence the protein function. In this study, we perform a detailed computational investigation based on the analysis of ab-initio molecular dynamics simulations, to disentangle the main structural quantities determining the fine balance in the chromophore environment. We found that specific hydrogen bonds interactions directly involving the chromophore (or not), are correlated to quantities, such as the volume of the cavity in which the chromophore is embedded and that it is importantly affected by the chromophore protonation state. The cross-correlation analysis performed on some of these hydrogen bonds and the cavity volume, demonstrates a direct correlation among them and we also identified the ones specifically involved in this correlation. We also found that specific interactions among residues far in the space are correlated, demonstrating the complexity of the chromophore environment and that many structural quantities have to be taken into account to properly describe and understand the main factors tuning the active site of the protein. From an overall evaluation of the results obtained in this work, it is shown that the residues which a priori are perceived to be spectators play instead an important role in both influencing the chromophore environment (cavity volume) and its dynamics (cross-correlations among spatially distant residues).

Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical–computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute–solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural–vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.

Light induces non-equilibrium time evolving molecular phenomena. The computational modeling of photo-induced processes in large systems, embedded in complex environments (i.e., solutions, proteins, materials), demands for a quantum and statistical mechanic treatment to achieve the required accuracy in the description of both the excited-state energy potentials and the choice of the initial conditions for dynamical simulations. On the other hand, the theoretical investigation on the atomistic scale of times and sizes of the ultrafast photo-induced reactivity and non-equilibrium relaxation dynamics right upon excitation requests tailored computational protocols. These methods often exploit hierarchic computation schemes, where a large part of the degrees of freedom are required to be treated explicitly to achieve the right accuracy. Additionally, part of the explicit system needs to be treated at ab initio level, where density functional theory, using hybrid functionals, represents a good compromise between accuracy and computational cost, when proton transfers, non-covalent interactions, and hydrogen bond dynamics play important roles. Thus, the modeling strategies presented in this review stress the importance of hierarchical quantum/molecular mechanics with effective non-periodic boundary conditions and efficient phase-sampling schemes to achieve chemical accuracy in ultrafast time-resolved spectroscopy and photo-induced phenomena. These approaches can allow explicit and accurate treatment of molecule/environment interactions, including also the electrostatic and dispersion forces of the bulk. At the same time, the specificities of the different case studies of photo-induced phenomena in solutions and biological environments are highlighted and discussed, with special attention to the computational and modeling challenges.

We develop a theoretical–computational protocol to simulate Joule heating on nanocomposite materials, allowing to investigate morphology, voltage and concentration effects.