Search

The ongoing integration of quantum chemistry, statistical mechanics, and artificial intelligence is paving the route toward more effective and accurate strategies for the investigation of the spectroscopic properties of medium-to-large size chromo- phores in condensed phases. In this context we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available computational tools. In this paper we report our latest developments with specific reference to unsupervised atomistic simulations employing non periodic boundary conditions (NPBC) followed by clustering of the trajectories employing optimized feature spaces. Next accurate variational computations are performed for a representative point of each cluster, whereas intracluster fluctuations are taken into account by a cheap yet reliable perturbative approach. A number of methodological improvements have been introduced including, e.g., more realistic reaction field effects at the outer boundary of the simulation sphere, automatic definition of the feature space by continuous perception of solute−solvent interactions, full account of polarization and charge transfer in the first solvation shell, and inclusion of vibronic contributions. After its validation, this new approach has been applied to the challenging case of solvatochromic effects on the UV−vis spectra of a prototypical nitroxide radical (TEMPO) in different solvents. The reliability, effectiveness, and robustness of the new platform is demonstrated by the remarkable agreement with experiment of the results obtained through an unsupervised approach characterized by a strongly reduced computational cost as compared to that of conventional quantum mechanics and molecular mechanics models without any accuracy reduction.

The local structure of the hydration of pyrazole has been analysed through static and dynamical microsolvation models described by quantum mechanical methods. Then, a reliable classical force field of pyrazole has been obtained on the basis of the quantum mechanical results and the dynamical properties of aqueous pyrazole solutions have been studied by molecular dynamics simulations. Finally, the structure of pyrazole-water solutions at different concentrations has been investigated by energy dispersive X-ray diffraction and experimental results have been compared to calculations. This comparison provides both a tool for interpretation of experiments and a way to validate the computational protocol.

The ultraviolet-visible absorption spectrum of camptothecin (CPT) has been been recorded in aqueous solution at pH 5.3, where the equilibrium among the different CPT forms is shifted toward the lactonic one. Time-dependent density functional theory (TD-DFT) computations lead to a remarkable reproduction of the experimental spectrum only upon addition of explicit water molecules in interaction with specific moieties of the camptothecin molecule. Molecular dynamics (MD) simulations enforcing boundary periodic conditions for CPT embedded with 865 water molecules, with a force field derived from DFT computations, show that the experimental spectrum is due to the contributions of CPT molecules with different solvation patterns. A similar solvent effect is observed for several CPT derivatives, including the clinically relevant SN-38 and topotecan drugs. The quantitative agreement between TD-DFT/MD computations and experimental data allow us to identify specific spectroscopic signatures diagnostic of the drug environment and to develop procedures that can be used to monitor the drug-DNA/protein interaction.

Multiscale methods combining quantum mechanics and molecular mechanics (QM/MM) have become the most suitable and effective strategies for the investigation of the spectroscopic properties of medium-to-large size chromophores in condensed phases. In this context, we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available computational tools. In this paper, we report our latest developments with specific reference to a general protocol based on atomistic simulations, carried out under nonperiodic boundary conditions (NPBC). In particular, we add to our in house MD engine a new efficient treatment of mean field electrostatic contributions to energy and forces, together with the capability of performing the simulations either in the canonical (NVT) or in the isothermal-isobaric (NPT) ensemble. Next, we provide convincing evidence that the NBPC approach enhanced by specific tweaks for rigid body propagation, allows for the simulation of solute-solvent systems with a minimum number of degrees of freedom and large integration time step. After its validation, this new approach is applied to the challenging case of solvatochromic effects on the electron paramagnetic resonance (EPR) spectrum of a prototypical nitroxide radical. To this end, we propose and validate also an automated protocol to extract and weight simulation snapshots, making use of a continuous description of the strength of solute-solvent hydrogen bridges. While further developments are being worked on, the effectiveness of our approach, even in its present form, is demonstrated by the accuracy of the results obtained through an unsupervised approach characterized by a strongly reduced computational cost as compared to that of conventional QM/MM models, without any appreciable deterioration of accuracy.

The Perturbed Matrix Method (PMM) approach to be used in combination with Molecular Dynamics (MD) trajectories (MD-PMM) has been recoded from scratch, improved in several aspects, and implemented in the Gaussian suite of programs for allowing a user-friendly and yet flexible tool to estimate quantum chemistry observables in complex systems in condensed phases. Particular attention has been devoted to a description of rigid and flexible quantum centers together with powerful essential dynamics and clustering approaches. The default implementation is fully black-box and does not require any external action concerning both MD and PMM sections. At the same time, fine-tuning of different parameters and use of external data are allowed in all the steps of the procedure. Two specific systems (Tyrosine and Uridine) have been reinvestigated with the new version of the code in order to validate the implementation, check the performances, and illustrate some new features.

Computational modeling involving Quantum Mechanics (QM) and Molecular Mechanics (MM) calculations are widely utilized to unveil the atomic-molecular properties that underpin their inherent characteristic features. The choice over the either of the QM and MM methods or a multiscale composite approach is driven by the target property of interest, and of course, the molecular size. Often, tailor-made schemes need to be devised to match the specific study purpose. Herein, we provide a perspective of these approaches addressing their effectiveness in terms of the delicate balance between the accuracy and computational feasibility. We focus on representative examples to highlight how different approaches can be fruitfully exploited for modeling the conformational landscape, and possibly, the spectroscopic behavior of biochemical molecules, especially amino acids building blocks.

The low-lying valence electronic excitations of the topotecan anticancer drug, in two stable lactone forms, have been addressed in infinite dilute aqueous solution by combining time-dependent density functional theory calculations with nanoseconds time-scale classical molecular dynamics simulations at 298 K. The effects of the surrounding and fluctuating classical environment on the investigated topotecan forms are included in a perturbed electronic Hamiltonian, which is computed, and then diagonalized, at each frame stored during the molecular dynamics sampling in explicit solution. Current results clearly indicate that, at moderately acidic and physiological conditions, the valence UV-vis absorption spectra of topotecan drug are strongly affected by the surrounding dielectric media and by its perturbing trajectory as arising from finite-temperature fluctuations and supramolecular interactions. Furthermore, the extension of the proposed computational study to hydrated topotecan complexes in liquid water shows that all of the experimentally detected UV-vis spectroscopic features in solution are accurately reproduced only when direct solute-solvent intermolecular interactions are also explicitly taken into account in our simulating scenario. Finally, the present investigation opens up a chance regarding the computational prediction of the UV-vis absorption spectra of topotecan interacting, in silico, with the topoisomerase-DNA binary complex in physiological conditions (i.e., water dilute solution, room temperature).

Molecular dynamics simulations of the Hg2+ ion in aqueous solution have been carried out using an effective two-body potential derived from quantum mechanical calculations. A stable heptacoordinated structure of the Hg2+ first hydration shell has been observed and confirmed by extended X-ray absorption fine structure (EXAFS) experimental data. The structural properties of the Hg2+ hydration shells have been investigated using radial and angular distribution functions, while the dynamical behavior has been discussed in terms of reorientational correlation functions, mean residence times of water molecules in the first and second hydration shells, and self-diffusion coefficients. The effect of water-water interactions on the Hg2+ hydration properties has been evaluated using the SPC/E and TIP5P water models.

A quite unexpected sevenfold coordination of the hydrated Hg(II) complex in aqueous solution is revealed by an extensive study combining X-ray absorption spectroscopy (XAS) and quantum mechanics/molecular dynamics (QM/MD) calculations. As a matter of fact, the generally accepted octahedral solvation of Hg(II) ion cannot be reconciled with XAS results. Next, refined QM computations point out the remarkable stability of a heptacoordinated structure with C 2 symmetry, and long-time MD simulations by new interaction potentials including many-body effects reveal that the hydrated complex has a quite flexible structure, corresponding for most of the time to heptacoordinated species. This picture is fully consistent with X-ray absorption near-edge structure experimental data which unambiguously show the preference for a sevenfold instead of a sixfold coordination.