Modelling vibrational coupling in DNA oligomers: a computational strategy combining QM and continuum solvation models

Effective fully polarizable QM/MM approaches to compute Raman and Raman Optical Activity spectra in aqueous solution

Raman and Raman Optical Activity (ROA) signals are amply affected by solvent effects, especially in the presence of strongly solute-solvent interactions such as Hydrogen Bonding (HB). In this work, we extend the fully atomistic polarizable Quantum Mechanics/Molecular Mechanics approach, based on the Fluctuating Charges and Fluctuating Dipoles force field to the calculation of Raman and ROA spectra. Such an approach is able to accurately describe specific HB interactions, by also accounting for anisotropic contributions due to the inclusion of fluctuating dipoles. To highlight the potentiality of the novel approach, Raman and ROA spectra of L-Serine and L-Cysteine dissolved in aqueous solution are computed and compared both with alternative theoretical approaches and experimental measurements.

UV-Resonance Raman Spectra of Systems in Complex Environments: A Multiscale Modeling Applied to Doxorubicin Intercalated into DNA

UV-Resonance Raman (RR) spectroscopy is a valuable tool to study the binding of drugs to biomolecular receptors. The extraction of information at the molecular level from experimental RR spectra is made much easier and more complete thanks to the use of computational approaches, specifically tuned to deal with the complexity of the supramolecular system. In this paper, we propose a protocol to simulate RR spectra of complex systems at different levels of sophistication, by exploiting a quantum mechanics/molecular mechanics (QM/MM) approach. The approach is challenged to investigate RR spectra of a widely used chemotherapy drug, doxorubicin (DOX) intercalated into a DNA double strand. The computed results show good agreement with experimental data, thus confirming the reliability of the computational protocol.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

A 4,4'-bis(2-benzoxazolyl)stilbene luminescent probe: assessment of aggregate formation through photophysics experiments and quantum-chemical calculations

A combination of experimental and quantum mechanical investigations is applied to the study of the optical features of 4,4'-bis(2-benzoxazolyl)stilbene (BBS) dissolved in solution or in a poly(l-lactic acid) (PLA) thermoplastic matrix at different concentrations. The experimental analyses allow the characterization of BBS solutions and dispersions in terms of absorption and emission features, along with the collection of some key parameters such as fluorescence quantum yield, anisotropy and lifetime, while the computational approach gives a detailed description of the photophysical behavior of BBS in the different environments. For the 10^-5 M BBS solution, the fluorescence spectra show the expected peaks at 425 and 455 nm of the non-interacting BBS molecules with a single fluorescence lifetime of 0.85 ns without revealing any aggregation phenomena, prevented by the short lifetime and fast diffusion rate of the monomer. Moreover, the calculated spectra are in excellent agreement with the experiments, thus showing the reliability of the computational approach. In time-resolved emission experiments (TRES) on more concentrated solutions (10^-4 M) and on BBS crystals, the presence of an excimer is revealed by the appearance of a broad peak around 540 nm, followed by the disappearance of the two main peaks at 460 nm on a time scale of about 10 ns. The computational analysis attributes this behavior to the formation of aggregates of different geometries. The BBS dispersions in PLA reveal the presence of different BBS architectures depending on the fluorophore content. Even at low concentrations, BBS is mainly dispersed as a monomer in the matrix, spheroid aggregates of about 800-900 nm in diameter are also present and the relevant fluorescence spectra arise from the combination of monomer and aggregate contributions. At higher concentrations, BBS starts forming crystals of a peculiar helicoidal shape, with a diameter of about 2 μm, variable length up to several hundreds of μm and emission spectra similar to those of isolated BBS crystals.

Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution

A polarizable QM/MM approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems is proposed and applied to aqueous solutions of (l)-methyl lactate and (S)-glycidol.

Effective Fully Polarizable QM/MM Approach To Model Vibrational Circular Dichroism Spectra of Systems in Aqueous Solution

We propose a methodology, based on the combination of classical Molecular Dynamics (MD) simulations with a fully polarizable Quantum Mechanical (QM)/Molecular Mechanics (MM)/Polarizable Continuum Model (PCM) Hamiltonian, to calculate Vibrational Circular Dichroism (VCD) spectra of chiral systems in aqueous solution. Polarization effects are included in the MM force field by exploiting an approach based on Fluctuating Charges (FQ). By performing the MD, the description of the solvating environment is enriched by taking into account the dynamical aspects of the solute-solvent interactions. On the other hand, the QM/FQ/PCM calculation of the VCD spectrum ensures an accurate description of the electronic density of the solute and a proper account for the specific interactions in solution. The application of our approach to (R)-methyloxirane and (l)-alanine in aqueous solution gives calculated spectra in remarkable agreement with their experimental counterparts and a substantial improvement with respect to the same spectra calculated with the PCM.

Intrinsically disordered proteins--relation to general model expressing the active role of the water environment

This work discusses the role of unstructured polypeptide chain fragments in shaping the protein's hydrophobic core. Based on the "fuzzy oil drop" model, which assumes an idealized distribution of hydrophobicity density described by the 3D Gaussian, we can determine which fragments make up the core and pinpoint residues whose location conflicts with theoretical predictions. We show that the structural influence of the water environment determines the positions of disordered fragments, leading to the formation of a hydrophobic core overlaid by a hydrophilic mantle. This phenomenon is further described by studying selected proteins which are known to be unstable and contain intrinsically disordered fragments. Their properties are established quantitatively, explaining the causative relation between the protein's structure and function and facilitating further comparative analyses of various structural models.

An investigation of the photophysical properties of minor groove bound and intercalated DAPI through quantum-mechanical and spectroscopic tools

The fluorescent probe 4',6-diamidino-2-phenylindole (DAPI) is a dye known to interact with polynucleotides in a non-univocal manner, both intercalation and minor groove binding modes being possible, and to specifically change its photophysical properties according to the different environments. To investigate this behavior, quantum-mechanical calculations using time-dependent density functional theory (TDDFT), coupled with polarizable continuum and/or atomistic models, were performed in combination with spectroscopic measurements of the probe in the different environments, ranging from a homogeneous solution to the minor groove or intercalation pockets of double stranded nucleic acids. According to our simulation, the electronic transition involves a displacement of the electron charge towards the external amidine groups and this feature makes the absorption energies very environment-sensitive while a much smaller sensitivity is seen in the fluorescence energies. Moreover, the calculations show that the DAPI molecule, when minor groove bound to the nucleic acid, presents both a reduced geometrical flexibility because of the rigid DNA pocket and a reduced polarization due to the very "apolar" microenvironment. All these effects can be used to better understand the observed enhancement of the fluorescence, which makes it an excellent marker for DNA.