Modeling Electronic Circular Dichroism within the Polarizable Embedding Approach

Importance of Polarizable Embedding for Computing Optical Rotation: The Case of Camphor in Ethanol

Solvation effects on optical rotation are notoriously challenging to model for computational chemistry, as the specific rotatory power of a molecule can vary wildly going from apolar to polar or even protic solvents. To address such a problem, we present a polarizable embedding implementation of an electric and magnetic response property based on density functional theory and the AMOEBA polarizable force field, and apply such an implementation to the study of the optical rotation of camphor in ethanol. By comparing a continuum model, and electrostatic and polarizable embedding QM/MM models, we observe that accounting for the environment's polarization gives rise to not only a different quantitative prediction, in very good agreement with experiments for the QM/AMOEBA model, but also to a very different qualitative picture, with the values of the optical rotation computed along a classical molecular dynamics trajectory with electrostatic embedding being statistically uncorrelated to the ones obtained with the polarizable description.

Close-Up Look at Electronic Spectroscopic Signatures of Common Pharmaceuticals in Solution

Simulating electronic properties and spectral signals requires robust computational approaches that need tuning with the system's peculiarities. In this paper, we test implicit and fully atomistic solvation models for the calculation of UV-vis and electronic circular dichroism (ECD) spectra of two pharmaceutically relevant molecules, namely, (2S)-captopril and (S)-naproxen, dissolved in aqueous solution. Room temperature molecular dynamics simulations reveal that these two drugs establish strong contacts with the surrounding solvent molecules via hydrogen bonds. Such specific interactions, which play a major role in the spectral response and are neglected in implicit approaches, are further characterized and quantified with natural bond orbital methods. Our calculations show that simulated spectra, and especially ECD, are in good agreement with experiments solely when conformational and configurational dynamics, mutual polarization, and solute-solvent repulsion effects are considered.

Polarizable Embedding Complex Polarization Propagator in Four- and Two-Component Frameworks

Explicit embedding methods combined with the complex polarization propagator (CPP) enable the modeling of spectroscopy for increasingly complex systems with a high density of states. We present the first derivation and implementation of the CPP in four- and exact-two-component (X2C) polarizable embedding (PE) frameworks. We denote the developed methods PE-4c-CPP and PE-X2C-CPP, respectively. We illustrate the methods by estimating the solvent effect on ultraviolet-visible (UV-vis) and X-ray atomic absorption (XAS) spectra of [Rh(H₂O)₆]³⁺ and [Ir(H₂O)₆]³⁺ immersed in aqueous solution. We moreover estimate solvent effects on UV-vis spectra of a platinum complex that can be photochemically activated (in water) to kill cancer cells. Our results clearly show that the inclusion of the environment is required: UV-vis and (to a lesser degree) XAS spectra can become qualitatively different from vacuum calculations. Comparison of PE-4c-CPP and PE-X2C-CPP methods shows that X2C essentially reproduces the solvent effect obtained with the 4c methods.

Natural and magnetic circular dichroism spectra of nucleosides: effect of the dynamics and environment

Chiroptical spectroscopic methods are excellent tools to study structure and interactions of biomolecules. However, their sensitivity to different structural aspects varies. To understand the dependence of absorption, electronic and magnetic circular dichroism (ECD, MCD) intensities on the structure, dynamics and environment, we measured and simulated spectra of nucleosides and other nucleic acid model components. The conformation space was explored by molecular dynamics (MD), the electronic spectra were generated using time dependent density functional theory (TDDFT). The sum over state (SOS) method was employed for MCD. The results show that accounting for the dynamics is crucial for reproduction of the experiment. While unpolarized absorption spectroscopy is relatively indifferent, ECD reflects the conformation and geometry dispersion more. MCD spectra provide variable response dependent on the wavelength and structural change. In general, MCD samples the structure more locally than ECD. Simple computational tests suggest that the optical spectroscopies coupled with the computational tools provide useful information about nucleic acid components, including base pairing and stacking.

Characterization of β-turns by electronic circular dichroism spectroscopy: a coupled molecular dynamics and time-dependent density functional theory computational study

Electronic circular dichroism is one of the most used spectroscopic techniques for peptide and protein structural characterization. However, while valuable experimental spectra exist for α-helix, β-sheet and random coil secondary structures, previous studies showed important discrepancies for β-turns, limiting their use as a reference for structural studies. In this paper, we simulated circular dichroism spectra for the best-characterized β-turns in peptides, namely types I, II, I' and II'. In particular, by combining classical molecular dynamics simulations and state-of-the-art quantum time-dependent density functional theory (with the polarizable embedding multiscale model) computations, two common electronic circular dichroism patterns were found for couples of β-turn types (namely, type I/type II' and type II/type I'), at first for a minimal di-peptide model (Ace-Ala-Ala-NHMe), but also for all sequences tested with non-aromatic residues in the central positions. On the other hand, as expected, aromatic substitution causes important perturbations to the previously found ECD patterns. Finally, by applying suitable approximations, these patterns were subsequently rationalized based on the exciton chirality rule. All these results provide useful predictions and pave the way for a possible experimental characterization of β-turns based on circular dichroism spectroscopy.

The shape of the electronic circular dichroism spectrum of (2,6-dimethylphenyl)(phenyl)methanol: interplay between conformational equilibria and vibronic effects

Analysis of the interplay between conformational equilibria, solvent effects and vibronic contributions in the ECD spectra.