Correct Modeling of Cisplatin: a Paradigmatic Case

Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order

The unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.

In Vitro and In Silico Vibrational-Rotational Spectroscopic Characterization of the Next-Generation Refrigerant HFO-1123

Very short-lived substances have recently been proposed as replacements for hydrofluorocarbons (HFCs), in turn being used in place of ozone-depleting substances, in refrigerant applications. In this respect, hydro-fluoro-olefins (HFOs) are attracting particular interest because, due to their reduced global warming potential, they are supposed to be environmentally friendlier. Notwithstanding this feature, they represent a new class of compounds whose spectroscopic properties and reactivity need to be characterized to allow their atmospheric monitoring and to understand their environmental fate. In the present work, the structural, vibrational, and ro-vibrational properties of trifluorothene (HFO-1123, F2C = CHF) are studied by state-of-the-art quantum chemical calculations. The equilibrium molecular structure has an expected error within 2 mÅ and 0.2° for bond lengths and angles, respectively. This represents the first step toward the computation of highly accurate rotational constants for both the ground and first excited fundamental vibrational levels, which reproduce the available experimental data well within 0.1%. Centrifugal distortion parameters and vibrational-rotational coupling terms are computed as well and used to solve some conflicting experimental results. Simulation of the vibrational transition frequencies and intensities beyond the double harmonic approximation and up to three quanta of vibrational excitation provides insights into the couplings ruling the vibrational dynamics and guides the characterization of the gas-phase infrared spectrum experimentally recorded in the range of 200-5000 cm-1. The full characterization of the IR features is completed with the experimental determination of the absorption cross sections over the 400-5000 cm-1 region from which the radiative forcing and global warming potential of HFO-1123 are derived.

Accurate Quantum Chemical Spectroscopic Characterization of Glycolic Acid: A Route Toward its Astrophysical Detection

The first step to shed light on the abiotic synthesis of biochemical building blocks, and their further evolution toward biological systems, is the detection of the relevant species in astronomical environments, including earthlike planets. To this end, the species of interest need to be accurately characterized from structural, energetic, and spectroscopic viewpoints. This task is particularly challenging when dealing with flexible systems, whose spectroscopic signature is ruled by the interplay of small- and large-amplitude motions (SAMs and LAMs, respectively) and is further tuned by the conformational equilibrium. In such instances, quantum chemical (QC) calculations represent an invaluable tool for assisting the interpretation of laboratory measurements or even observations. In the present work, the role of QC results is illustrated with reference to glycolic acid (CH2OHCOOH), a molecule involved in photosynthesis and plant respiration and a precursor of oxalate in humans, which has been detected in the Murchison meteorite but not yet in the interstellar medium or in planetary atmospheres. In particular, the equilibrium structure of the lowest-energy conformer is derived by employing the so-called semiexperimental approach. Then, accurate yet cost-effective QC calculations relying on composite post-Hartree-Fock schemes and hybrid coupled-cluster/density functional theory approaches are used to predict the structural and ro-vibrational spectroscopic properties of the different conformers within the framework of the second-order vibrational perturbation theory. A purposely tailored discrete variable representation anharmonic approach is used to treat the LAMs related to internal rotations. The computed spectroscopic data, particularly those in the infrared region, complement the available experimental investigations, thus enhancing the possibility of an astronomical detection of this molecule.

The Spectroscopic Characterization of Halogenated Pollutants through the Interplay between Theory and Experiment: Application to R1122

In the last decade, halogenated ethenes have seen an increasing interest for different applications; in particular, in refrigeration, air-conditioning and heat pumping. At the same time, their adverse effects as atmospheric pollutants require environmental monitoring, especially by remote sensing spectroscopic techniques. For this purpose, an accurate characterization of the spectroscopic fingerprint—in particular, those of relevance for rotational–vibrational spectroscopy—of the target molecules is strongly needed. This work provides an integrated computational–theoretical investigation on R1122 (2-Chloro-1,1-difluoro-ethylene, ClHC=CF2), a compound widely employed as a key intermediate in different chemical processes. State-of-the-art quantum chemical calculations relying on CCSD(T)-based composite schemes and hybrid CCSD(T)/DFT approaches are used to obtain an accurate prediction of the structural, rotational and vibrational spectroscopic properties. In addition, the equilibrium geometry is obtained by exploiting the semi-experimental method. The theoretical predictions are used to guide the analysis of the experimentally recorded gas-phase infrared spectrum, which is assigned in the 400–6500 cm−1 region. Furthermore, absorption cross sections are accurately determined over the same spectral range. Finally, by using the obtained spectroscopic data, a first estimate of the global warming potential of R1122 vibrational spectra is obtained.

Systematic assessment of DFT methods for geometry optimization of mononuclear platinum-containing complexes

Platinum is used extensively as a catalyst for a wide variety of chemical reactions, though its scarcity and price present limitations to expansions of its use. To understand the origin of platinum's versatility-with the goals of both improving the efficiency of existing catalysts and mimicking its reactivity with more abundant metals-the mechanisms of platinum-catalyzed chemical reactions must be understood via structural and spectroscopic characterization of these catalysts under operando conditions. Such data, typically consisting of complex mixtures of species, often prove challenging to interpret, inviting the aid of chemical theory. DFT calculations in particular have proven successful at predicting structural and spectroscopic parameters of transition metal species, though a thorough investigation of how these methods perform for platinum-based complexes has yet to be undertaken. Herein, we evaluated the performance of geometry optimization for five commonly used functionals (BP86, PBE, B3LYP, PBE0, and TPSSh) in combination with various ligand basis sets, relativistic approximations, and solvation and dispersion models. We applied these DFT methods to a training set of 14 platinum-containing complexes with varying sizes, oxidation states, and number and type of ligands and determined that the best-performing method was the PBE0 functional together with the def2-TZVP basis set for the ligand atoms, the ZORA relativistic approximation, and solvation and dispersion corrections. The ability of this DFT methodology to accurately predict metrical parameters was confirmed using two case studies, most notably by comparing the DFT optimized geometry of a previously uncharacterized complex to newly collected EXAFS data, which showed excellent agreement.

Accurate Biomolecular Structures by the Nano-LEGO Approach: Pick the Bricks and Build Your Geometry

The determination of accurate equilibrium molecular structures plays a fundamental role for understanding many physical–chemical properties of molecules, ranging from the precise evaluation of the electronic structure to the analysis of the role played by dynamical and environmental effects in tuning their overall behavior. For small semi-rigid systems in the gas phase, state-of-the-art quantum chemical computations rival the most sophisticated experimental (from, for example, high-resolution spectroscopy) results. For larger molecules, more effective computational approaches must be devised. To this end, we have further enlarged the compilation of available semi-experimental (SE) equilibrium structures, now covering the most important fragments containing H, B, C, N, O, F, P, S, and Cl atoms collected in the new SE100 database. Next, comparison with geometries optimized by methods rooted in the density functional theory showed that the already remarkable results delivered by PW6B95 and, especially, rev-DSDPBEP86 functionals can be further improved by a linear regression (LR) approach. Use of template fragments (taken from the SE100 library) together with LR estimates for the missing interfragment parameters paves the route toward accurate structures of large molecules, as witnessed by the very small deviations between computed and experimental rotational constants. The whole approach has been implemented in a user-friendly tool, termed nano-LEGO, and applied to a number of demanding case studies.

Advances in Our Understanding of the Molecular Mechanisms of Action of Cisplatin in Cancer Therapy

Cisplatin and other platinum-based chemotherapeutic drugs have been used extensively for the treatment of human cancers such as bladder, blood, breast, cervical, esophageal, head and neck, lung, ovarian, testicular cancers, and sarcoma. Cisplatin is commonly administered intravenously as a first-line chemotherapy for patients suffering from various malignancies. Upon absorption into the cancer cell, cisplatin interacts with cellular macromolecules and exerts its cytotoxic effects through a series of biochemical mechanisms by binding to Deoxyribonucleic acid (DNA) and forming intra-strand DNA adducts leading to the inhibition of DNA synthesis and cell growth. Its primary molecular mechanism of action has been associated with the induction of both intrinsic and extrinsic pathways of apoptosis resulting from the production of reactive oxygen species through lipid peroxidation, activation of various signal transduction pathways, induction of p53 signaling and cell cycle arrest, upregulation of pro-apoptotic genes/proteins, and down-regulation of proto-oncogenes and anti-apoptotic genes/proteins. Despite great clinical outcomes, many studies have reported substantial side effects associated with cisplatin monotherapy, while others have shown substantial drug resistance in some cancer patients. Hence, new formulations and several combinational therapies with other drugs have been tested for the purpose of improving the clinical utility of cisplatin. Therefore, this review provides a comprehensive understanding of its molecular mechanisms of action in cancer therapy and discusses the therapeutic approaches to overcome cisplatin resistance and side effects.

Cyanide Bridged Platinum-Iron Complexes as Cisplatin Prodrug Systems: Design and Computational Study

The potential role of cyanide-bridged platinum-iron complexes as an anti-cancer Pt(IV) prodrug is studied. We present design principles of a dual-function prodrug that can upon reduction dissociate and release concurrently six cisplatin units and a ferricyanide anion per prodrug unit. The prodrug molecule is a unique complex of hepta metal centers consisting of a ferricyanide core with six Pt(IV) centers each bonded to the Fe(III) core through a cyano ligand. The functionality of the prodrug is addressed through density functional theory (DFT) calculations.

Accuracy Meets Interpretability for Computational Spectroscopy by Means of Hybrid and Double-Hybrid Functionals

Accuracy and interpretability are often seen as the devil and holy grail in computational spectroscopy and their reconciliation remains a primary research goal. In the last few decades, density functional theory has revolutionized the situation, paving the way to reliable yet effective models for medium size molecules, which could also be profitably used by non-specialists. In this contribution we will compare the results of some widely used hybrid and double hybrid functionals with the aim of defining the most suitable recipe for all the spectroscopic parameters of interest in rotational and vibrational spectroscopy, going beyond the rigid rotor/harmonic oscillator model. We will show that last-generation hybrid and double hybrid functionals in conjunction with partially augmented double- and triple-zeta basis sets can offer, in the framework of second order vibrational perturbation theory, a general, robust, and user-friendly tool with unprecedented accuracy for medium-size semi-rigid molecules.

Metal-Mediated Morphology Regulation of Self-Assembled Double-Hydrophilic Block Copolymers

We report a nanoscale morphology-regulation strategy of self-assembled double-hydrophilic block copolymers with square planar Pt^II compounds. The selective coordination of Pt^II on the chelating blocks of poly(acrylic acid)-b-poly(ethylene glycol) (PAA-b-PEG) induced the self-association of metal-chelated unimers by the known cohesive force of Pt^II. The block-length variation of PAA with constant PEG led to the shape transition from normal core/shell and crew-cut spheres to anisotropic pearl-string structures. On the other hand, Pt^II adsorption on PEG blocks by extensive hydrogen bonding can further modify the molecular geometry of metal-chelated unimers by decreasing the volume of hydrophilic segments, eventually leading to the shape transition to vesicular structures. This result was well correlated to the structural constraint of PEG conformation estimated by the quantitative ¹H NMR analysis. The vesicles also exhibited the enclosing nature for the fluorescent guest molecules, which demonstrated the promising potential for the encapsulating delivery vehicle.

Influence of DFT Functionals and Solvation Models on the Prediction of Far-Infrared Spectra of Pt-Based Anticancer Drugs: Why Do Different Complexes Require Different Levels of Theory?

Computational modeling was applied to far-infrared (FIR) spectra of Pt-based anticancer drugs to study the hydrolysis of these important molecules. Here, we present a study that investigates the influence of different factors-basis sets on non-Pt atoms, relativistic effective core potentials (RECPs) on the Pt atom, density functional theory (DFT) functionals, and solvation models-on the prediction of FIR spectra of two Pt-based anticancer drugs, cisplatin and carboplatin. Geometry optimizations and frequency calculations were performed with a range of functionals (PBE, PBE0, M06-L, and M06-2X), Dunning's correlation-consisted basis sets (VDZ, VTZ, aVDZ, and aVTZ), RECPs (VDZ-pp, VTZ-pp, aVDZ-pp, and aVTZ-pp), and solvation models (IEFPCM, CPCM, and SMD). The best combination of the basis set/DFT functional/solvation model was identified for each anticancer drug by comparing with experimentally available FIR spectra. Different combinations were established for cisplatin and carboplatin, which was rationalized by means of the partial atomic charge scheme, ChelpG, that was utilized to study the charge transfer between the Pt ion and ligands in both cisplatin and carboplatin.

Diving for Accurate Structures in the Ocean of Molecular Systems with the Help of Spectroscopy and Quantum Chemistry

The prediction and interpretation of structural properties are the starting points for a deep understanding of thermochemistry, kinetics, and spectroscopic signatures of molecular systems. To give an example, detailed knowledge of the conformational behavior of the main building blocks of biomolecules in the gas phase (i.e., without the perturbing effect of the environment) is a mandatory prerequisite toward the understanding of the role played by different interactions in determining the biological activity in terms of structure-activity relationships. The first step to take is an unambiguous definition of molecular structure. We address the so-called Born-Oppenheimer equilibrium structure, which is defined in a rigorous manner and isotopically independent, and the target accuracy. For the latter, we aim at so-called "spectroscopic" accuracy, which implies uncertainties of a few milliangstroms for bond lengths and smaller than a tenth of degree for angles. If on one side the continuous enhancements of the experimental techniques give access to new and unprecedented spectroscopic determinations, on the other side they require increasing efforts for an unbiased interpretation and analysis. Among the pieces of information, accurate molecular structures play a particularly important role. Indeed, there is a strong relationship between the experimental outcome and the electronic structure of the system. Spectroscopic techniques, in particular those exploited in the gas phase, are therefore accurate and reliable sources for structural information. However, it is seldom straightforward to derive molecular structures directly from the experimental information. Indeed, even in the favorable case of investigations in the gas phase, vibrational effects are always present, and disentangling their contributions requires collection of information for all vibrational modes, a nearly impossible task. To overcome these limitations, joint theory-spectroscopy strategies can be identified, which are referred to as "top-down" and "bottom-up". The first approach, denoted as the semiexperimental approach, relies on extracting from experimental outcomes the equilibrium structure by using quantum-chemical computations to recover vibrational effects. The bottom-up approach consists in verifying the computed equilibrium geometry by means of a comparison between calculated and experimental spectroscopic parameters that probe structural characteristics. In this contribution, we try to review the most important challenges in accurate molecular structure determinations, with particular emphasis on the "solution" provided by a joint theoretical-experimental approach and on the current state of the art. Starting from the illustration of different strategies, we proceed by addressing the increasing complexity in the derivation of equilibrium geometries: we start from the construction of a database of accurate structures, we then face the problem of extending the dimension of the systems amenable to accurate structural determinations, and finally we move to the challenge of understanding the nature of intermolecular interactions.

Correct Modeling of Cisplatin: a Paradigmatic Case

Quantum chemistry is a useful tool in modern approaches to drug and material design, but only when the adopted model reflects a correct physical picture. Paradigmatic is the case of cis‐diaminodichloroplatinum(II), cis‐[Pt(NH₃)₂Cl₂], for which the correct simulation of the structural and vibrational properties measured experimentally still remains an open question. By using this molecule as a proof of concept, it is shown that state‐of‐the‐art quantum chemical calculations and a simple model, capturing the basic physical flavors, a cis‐[Pt(NH₃)₂Cl₂] dimer, can provide the accuracy required for interpretative purposes. The present outcomes have fundamental implications for benchmark studies aiming at assessing the accuracy of a given computational protocol.