Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach

Study of absorption and emission spectra of substituted terthiophene compounds by methods of theoretical chemistry

Terthiophene derivatives attract interest due to their prospective applications in optoelectronic or sensor devices. Due to their nontoxicity they can be considered as suitable candidates in biological applications. Supramolecular organization of the matter is one of the most interesting topics in contemporary materials science. Amphiphilic chromophores based on substituted terthiophenes are capable of self-assembly into supramolecular architectures. In this work, we aim at simulation of the spectral properties of terthiophene with oligo(ethylene oxide) substituents by the methods of quantum chemistry and molecular dynamics (MD). The potential energy surface (PES) of this molecule was determined by the methods of density functional theory (DFT) for the ground state and time-dependent density-functional theory (TD-DFT) for the excited state. MD simulations in water than revealed the most frequented molecular conformations in both these states. Absorption and fluorescence spectra were determined for all these conformations, including the surrounding water molecules, using TD-DFT and averaged over the conformation space to obtain the final absorption and fluorescence spectrum. The calculated spectra were compared with their experimental counterparts and the differences were discussed in context of the supramolecular structure revealed by confocal microscopy. In spite of its simplicity, this approach provides a satisfactory approximation of absorption and fluorescent spectra of these molecules obtained by computational methods.

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.

Storage Stability of Atheroglitatide, an Echogenic Liposomal Formulation of Pioglitazone Targeted to Advanced Atheroma with a Fibrin-Binding Peptide

We have conducted a stability study of a complex liposomal pharmaceutical product, Atheroglitatide (AGT), stored at three temperatures, 4, 24, and 37 °C, for up to six months. The six parameters measured were functions of liposomal integrity (size and number), drug payload (loading efficiency), targeting peptide integrity (conjugation efficiency and specific avidity), and echogenicity (ultrasound-dependent controlled drug release), which were considered most relevant to the product's intended use. At 4 °C, liposome diameter trended upward, indicative of aggregation, while liposome number per mg lipid and echogenicity trended downward. At 24 °C, peptide conjugation efficiency (CE) and targeting efficiency (TE, specific avidity) trended downward. At 37 °C, CE and drug (pioglitazone) loading efficiency trended downward. At 4 °C, the intended storage temperature, echogenicity, and liposome size reached their practical tolerance limits at 6 months, fixing the product expiration at that point. Arrhenius analysis of targeting peptide CE and drug loading efficiency decay at the higher temperatures indicated complete stability of these characteristics at 4 °C. The results of this study underscore the storage stability challenges presented by complex nanopharmaceutical formulations.

An Expedited Route to Optical and Electronic Properties at Finite Temperature via Unsupervised Learning

Electronic properties and absorption spectra are the grounds to investigate molecular electronic states and their interactions with the environment. Modeling and computations are required for the molecular understanding and design strategies of photo-active materials and sensors. However, the interpretation of such properties demands expensive computations and dealing with the interplay of electronic excited states with the conformational freedom of the chromophores in complex matrices (i.e., solvents, biomolecules, crystals) at finite temperature. Computational protocols combining time dependent density functional theory and ab initio molecular dynamics (MD) have become very powerful in this field, although they require still a large number of computations for a detailed reproduction of electronic properties, such as band shapes. Besides the ongoing research in more traditional computational chemistry fields, data analysis and machine learning methods have been increasingly employed as complementary approaches for efficient data exploration, prediction and model development, starting from the data resulting from MD simulations and electronic structure calculations. In this work, dataset reduction capabilities by unsupervised clustering techniques applied to MD trajectories are proposed and tested for the ab initio modeling of electronic absorption spectra of two challenging case studies: a non-covalent charge-transfer dimer and a ruthenium complex in solution at room temperature. The K-medoids clustering technique is applied and is proven to be able to reduce by ∼100 times the total cost of excited state calculations on an MD sampling with no loss in the accuracy and it also provides an easier understanding of the representative structures (medoids) to be analyzed on the molecular scale.

Directing the Morphology, Packing, and Properties of Chiral Metal–Organic Frameworks by Cation Exchange**

We show that metal–organic frameworks, based on tetrahedral pyridyl ligands, can be used as a morphological and structural template to form a series of isostructural crystals having different metal ions and properties. An iterative crystal‐to‐crystal conversion has been demonstrated by consecutive cation exchanges. The primary manganese‐based crystals are characterized by an uncommon space group (P622). The packing includes chiral channels that can mediate the cation exchange, as indicated by energy‐dispersive X‐ray spectroscopy on microtome‐sectioned crystals. The observed cation exchange is in excellent agreement with the Irving–Williams series (Mn<Fe<Co<Ni< Cu>Zn) associated with the relative stability of the resulting coordination nodes. Furthermore, we demonstrate how the metal cation controls the optical and magnetic properties. The crystals maintain their morphology, allowing a quantitative comparison of their properties at both the ensemble and single‐crystal level.

Development, Validation, and Pilot Application of a Generalized Fluctuating Charge Model for Computational Spectroscopy in Solution

A general approach enforcing nonperiodic boundary conditions for the computation of spectroscopic properties in solution has been improved including an effective description of charge-transfer contributions and coordination number adjustment for explicit solvent molecules. Both contributions are obtained from a continuous description of intermolecular hydrogen bonds, which has been employed also for an effective clustering of molecular dynamics trajectories. Fine tuning of the model has been performed for several water clusters, and then its efficiency and reliability have been demonstrated by computing the absorption spectra of different creatinine tautomers in aqueous solution.

How Hydrogen Bonding Amplifies Isomeric Differences in Pyridones toward Strong Changes in Acidity and Tautomerism

Steric hindrance of hydration and hydrogen bond enhancement by localized charges have been identified as key factors for the massive chemical differences between the hydroxypyridine/pyridone isomers in aqueous solution. While all isomers occur mainly in the hydroxypyridine form in the gas phase, they differ by more than 3 orders of magnitude both in their acidity and tautomeric equilibrium constants upon hydration. By monitoring the electronic and solvation structures as a function of the protonation state and the O^– substitution position on the pyridine ring, the amplification of the isomeric differences in aqueous solution has been investigated. Near-edge X-ray absorption fine structure (NEXAFS) measurements at the N K-edge served as the probe of the chemical state. The combination of molecular dynamics simulations, complete active space self-consistent field (CASSCF), and time-dependent density functional theory (TD-DFT) spectral calculations contributes to unraveling the principles of tautomerism and acidity in multiple biochemical systems based on tautomerism.

Combining multilevel Hartree–Fock and multilevel coupled cluster approaches with molecular mechanics: a study of electronic excitations in solutions

We present the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable.

Structural and Photophysical Properties of Various Polypyridyl Ligands: A Combined Experimental and Computational Study

Covalent triazine frameworks (CTFs) with polypyridyl ligands are very promising supports to anchor photocatalytic complexes. Herein, we investigate the photophysical properties of a series of ligands which vary by the extent of the aromatic system, the nitrogen content and their topologies to aid in selecting interesting building blocks for CTFs. Interestingly, some linkers have a rotational degree of freedom, allowing both a trans and cis structure, where only the latter allows anchoring. Therefore, the influence of the dihedral angle on the UV-Vis spectrum is studied. The photophysical properties are investigated by a combined computational and experimental study. Theoretically, both static and molecular dynamics simulations are performed to deduce ground- and excited state properties based on density functional theory (DFT) and time-dependent DFT. The position of the main absorption peak shifts towards higher wavelengths for an increased size of the π-system and a higher π-electron deficiency. We found that the position of the main absorption peak among the different ligands studied in this work can amount to 271 nm; which has a significant impact on the photophysical properties of the ligands. This broad range of shifts allows modulation of the electronic structure by varying the ligands and may help in a rational design of efficient photocatalysts.

Charge Anisotropy of Nitrogen: Where Chemical Intuition Fails

For more than half a century computer simulations were developed and employed to study ensemble properties of a wide variety of atomic and molecular systems with tremendous success. Nowadays, a selection of force-fields is available that describe the interactions in such systems. A key feature of force-fields is an adequate description of the electrostatic potential (ESP). Several force-fields model the ESP via point charges positioned at the atom centers. A major shortcoming of this approach, its inability to model anisotropies in the ESP, can be mitigated using additional charge sites. It has been shown that nitrogen is the most problematic element abundant in many polymers as well as large molecules of biological origin. To tackle this issue, small organic molecules containing a single nitrogen atom were studied. In performing rigorous scans of the surroundings of these nitrogen atoms, positions where a single extra charge can enhance the ESP description the most were identified. Significant improvements are found for ammonia, amines, and amides. Interestingly, the optimal location for the extra charge does not correlate with the chemically intuitive position of the nitrogen lone pair. In fact, the placement of an extra charge in the lone-pair location does not lead to significant improvements in most cases.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes

The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.

Quantum Confinement Effects on Solvatochromic Shifts of Molecular Solutes

We demonstrate the pivotal role of quantum mechanics density confinement effects on solvatochromic shifts. In particular, by resorting to a quantum mechanics/molecular mechanics (QM/MM) approach capable of accounting for confinement effects we successfully reproduce vacuo-to-water solvatochromic shifts for dark n → π* and bright π → π* transitions of acrolein and dark n → π* transitions of pyridine and pyrimidine without the need of including explicit water molecules in the QM portion. Remarkably, our approach is also able to dissect the effects of the single forces acting on the solute-solvent couple and allows for a rationalization of the experimental findings in terms of physicochemical quantities.

The importance of nuclear quantum effects in spectral line broadening of optical spectra and electrostatic properties in aromatic chromophores

In this work, we examine the importance of nuclear quantum effects on capturing the line broadening and vibronic structure of optical spectra. We determine the absorption spectra of three aromatic molecules indole, pyridine, and benzene using time dependent density functional theory with several molecular dynamics sampling protocols: force-field based empirical potentials, ab initio simulations, and finally path-integrals for the inclusion of nuclear quantum effects. We show that the absorption spectrum for all these chromophores are similarly broadened in the presence of nuclear quantum effects regardless of the presence of hydrogen bond donor or acceptor groups. We also show that simulations incorporating nuclear quantum effects are able to reproduce the heterogeneous broadening of the absorption spectra even with empirical force fields. The spectral broadening associated with nuclear quantum effects can be accounted for by the broadened distribution of chromophore size as revealed by a particle in the box model. We also highlight the role that nuclear quantum effects have on the underlying electronic structure of aromatic molecules as probed by various electrostatic properties.

Solvent effects for vertical absorption and emission processes in solution using a self-consistent state specific method based on constrained equilibrium thermodynamics

A self-consistent state specific (SS) method in the framework of TDDFT is presented to account for solvent effects on absorption and emission processes for molecules in solution. In these processes, the initial state is an equilibrium state, while the polarization of the solvent is in nonequilibrium with the electron density of the solute in the final state. Nonequilibrium solvation free energy is calculated based on a novel nonequilibrium solvation model with constrained equilibrium manipulation. The bulk solvent effects are considered using the polarizable continuum method (PCM), where the solvent-solute interaction is described with a reaction field. Molecular orbitals and orbital energies in the presence of the reaction field corresponding to the excited state are employed and the response of the solvent is not included in the TDDFT calculations. A self-consistent procedure is designed to obtain the excited state reaction field. The equations based on this new nonequilibrium solvation model in the framework of the self-consistent SS-PCM/TDDFT method for calculation of vertical absorption and emission energies are presented and implemented in the Q-Chem package. Vertical absorption and emission energies for several small molecules in solution using the newly developed code are calculated and compared with available experimental data and the results of other theoretical studies. Solvent shifts of absorption and emission energies are reasonably reproduced with this approach. The new model is a promising approach to study nonequilibrium absorption and emission processes in solution.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?

Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach

The electronic absorption spectra of pyridine and nicotine in aqueous solution have been computed using a multistep approach. The computational protocol consists in studying the solute solvation with accurate molecular dynamics simulations, characterizing the hydrogen bond interactions, and calculating electronic transitions for a series of configurations extracted from the molecular dynamics trajectories with a polarizable QM/MM scheme based on the fluctuating charge model. Molecular dynamics simulations and electronic transition calculations have been performed on both pyridine and nicotine. Furthermore, the contributions of solute vibrational effect on electronic absorption spectra have been taken into account in the so called vertical gradient approximation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.