Solvent Effects and Dynamic Averaging of 195Pt NMR Shielding in Cisplatin Derivatives

Dynamic and relativistic effects on Pt-Pt indirect spin-spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent-solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin-orbit Kohn-Sham (2c-ZKS) and four-component Dirac-Kohn-Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1-Pt2(NH3)4(Am)2-L2]n+ (Am = α-pyrrolidonate and pivalamidate; L = H2O, Cl-, and Br-) are studied. The computed Pt-Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N-C-O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X-) and aquahalo (L1 = H2O; L2 = X-) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Predicting 195 Pt NMR Chemical Shifts in Water-Soluble Inorganic/Organometallic Complexes with a Fast and Simple Protocol Combining Semiempirical Modeling and Machine Learning

Water-soluble Pt complexes are the key components in medicinal chemistry and catalysis. The well-known cisplatin family of anticancer drugs and industrial hydrosylilation catalysts are two leading examples. On the molecular level, the activity mechanisms of such complexes mostly involve changes in the Pt coordination sphere. Using ¹⁹⁵ Pt NMR spectroscopy for operando monitoring would be a valuable tool for uncovering the activity mechanisms; however, reliable approaches for the rapid correlation of Pt complex structure with ¹⁹⁵ Pt chemical shifts are very challenging and not available for everyday research practice. While NMR shielding is a response property, molecular 3D structure determines NMR spectra, as widely known, which allows us to build up 3D structure to ¹⁹⁵ Pt chemical shift correlations. Accordingly, we present a new workflow for the determination of lowest-energy configurational/conformational isomers based on the GFN2-xTB semiempirical method and prediction of corresponding chemical shifts with a Machine Learning (ML) model tuned for Pt complexes. The workflow was designed for the prediction of ¹⁹⁵ Pt chemical shifts of water-soluble Pt(II) and Pt(IV) anionic, neutral, and cationic complexes with halide, NO₂ ^- , (di)amino, and (di)carboxylate ligands with chemical shift values ranging from -6293 to 7090 ppm. The model offered an accuracy (normalized root-mean-square deviation/RMSD) of 1.08 %/145.02 ppm on the held-out test set.

Supramolecular Coronation of Platinum(II) Complexes by Macrocycles: Structure, Relativistic DFT Calculations, and Biological Effects

Platinum-based anticancer drugs are actively developed utilizing lipophilic ligands or drug carriers for the efficient penetration of biomembranes, reduction of side effects, and tumor targeting. We report the development of a supramolecular host-guest system built on cationic platinum(II) compounds bearing ligands anchored in the cavity of the macrocyclic host. The host-guest binding and hydrolysis process on the platinum core were investigated in detail by using NMR, MS, X-ray diffraction, and relativistic DFT calculations. The encapsulation process in cucurbit[7]uril unequivocally promotes the stability of hydrolyzed dicationic cis-[Pt^II(NH₃)₂(H₂O)(NH₂-R)]²⁺ compared to its trans isomer. Biological screening on the ovarian cancer lines A2780 and A2780/CP shows time-dependent toxicity. Notably, the reported complex and its β-cyclodextrin (β-CD) assembly achieve the same cellular uptake as cisplatin and cisplatin@β-CD, respectively, while maintaining a significantly lower toxicity profile.

Solvent effect on the 195Pt NMR properties in pyridonate-bridged PtIII dinuclear complex derivatives investigated by ab initio molecular dynamics and localized orbital analysis

An ab initio molecular dynamics investigation of the solvent effect (water) on the structural parameters, ¹⁹⁵Pt NMR spin-spin coupling constants (SSCCs) and chemical shifts of a series of pyridonate-bridged Pt^III dinuclear complexes is performed using Kohn-Sham (KS) Car-Parrinello molecular dynamics (CPMD) and relativistic hybrid KS NMR calculations. The indirect solvent effect (via structural changes) has a dramatic effect on the ¹J_PtPt SSCCs. The complexes exhibit a strong trans influence in solution, where the Pt-Pt bond lengthens with increasing axial ligand σ-donor strength. In the diaqua complex, where the solvent effect is more pronounced, the SSCCs averaged for CPMD configurations with explicit plus implicit solvation agree much better with the experimental data, while the calculations for static geometry and CPMD unsolvated configurations show large deviations with respect to experiment. The combination of CPMD with hybrid KS NMR calculations provides a much more realistic computational model that reproduces the large magnitudes of ¹J_PtPt and ¹⁹⁵Pt chemical shifts. An analysis of ¹J_PtPt in terms of localized and canonical orbitals shows that the SSCCs are driven by changes in the s-character of the natural atomic orbitals of Pt atoms, which affect the 'Fermi contact' mechanism.

1H and 195Pt NMR prediction for inclusion compounds formed by cisplatin and oxidized carbon nanostructures

Prediction of NMR chemical shifts can assist experimentalists in the characterization of drug delivery systems based on carbon nanocomposites. Chemical shifts are strongly correlated to the nucleus position and its chemical neighborhood. Therefore, to predict structures and NMR properties of complex chemical models, choosing a more consistent theoretical level capable of providing more realistic results and moderate computational demand is a major challenge. In this work, we predicted the NMR spectra of inclusion compounds formed by cisplatin (cDDP) and an oxidized carbon nanotube (CNTox) and nanocone (CNCox) considered by specialists as potential drug delivery systems. The ¹⁹⁵Pt NMR chemical shifts calculated at the DFT level with the new relativistic NMR-DKH basis set were −2314 ppm and −2192 ppm for cDDP@CNTox and cDDP@CNCox complexes, respectively, which are both high-field shifted relative to the free cDDP (−2110 ppm). ¹H NMR chemical shifts are also sensitive to the inclusion process. The H (NH₃) signals are found on average at +4.3 (cDDP), −5.1 (cDDP@CNTox) and +6.6 ppm (cDDP@CNCox). Interestingly, despite the similar inclusion modes in CNTox and CNCox cavities, the ¹H NMR shifts were in opposite directions. A possible reason might be the higher stability of cDDP@CNTox (ΔE_F = −19.9 kcal mol⁻¹) than that of cDDP@CNCox (ΔE_F = −5.7 kcal mol⁻¹), which suggests a short guest–host contact in the former and consequently, a more efficient shielding of hydrogen atoms due to the electron-rich carbon structure. These results may be helpful as comparison data in the NMR spectra assignment in solution and the inclusion compounds' structural elucidation.

¹⁹⁵Pt NMR chemical shifts of the cisplatin molecule in the inclusion complex formed by oxidized carbon nanotubes.

First-Principles Calculation of 1H NMR Chemical Shifts of Complex Metal Polyhydrides: The Essential Inclusion of Relativity and Dynamics

1H NMR spectroscopy has become an important technique for the characterization of transition-metal hydride complexes, whose metal-bound hydrides are often difficult to locate by X-ray diffraction. In this regard, the accurate prediction of 1H NMR chemical shifts provides a useful, but challenging, strategy to help in the interpretation of the experimental spectra. In this work, we establish a density-functional-theory protocol that includes relativistic, solvent, and dynamic effects at a high level of theory, allowing us to report an accurate and reliable interpretation of 1H NMR hydride chemical shifts of iridium polyhydride complexes. In particular, we have studied in detail the hydride chemical shifts of the [Ir6(IMe)8(CO)2H14]2+ complex in order to validate previous assignments. The computed 1H NMR chemical shifts are strongly dependent on the relativistic treatment, the choice of the DFT exchange-correlation functional, and the conformational dynamics. By combining a fully relativistic four-component electronic-structure treatment with ab initio molecular dynamics, we were able to reliably model both the terminal and bridging hydride chemical shifts and to show that two NMR hydride signals were inversely assigned in the experiment.

Al27 NMR chemical shift of Al(OH)4 - calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H₂O)₆ ³⁺, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)₄ ^- requires an understanding of how these errors cancel in order to define the limits of accurately predicting Al27 chemical shielding in Al(OH)₄ ^-. In this work, we estimate the absolute shielding of the Al27 nucleus in Al(OH)₄ ^- at the coupled cluster level (515.1 ± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H₂O)₆ ³⁺ and Al(OH)₄ ^- ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the Al27 chemical shift of Al(OH)₄ ^- from first principles.

Four-component relativistic 31P NMR calculations for trans-platinum(ii) complexes: importance of the solvent and dynamics in spectral simulations

We report a combined experimental-theoretical study on the 31P NMR chemical shift for a number of trans-platinum(ii) complexes. Validity and reliability of the 31P NMR chemical shift calculations are examined by comparing with the experimental data. A successful computational protocol for the accurate prediction of the 31P NMR chemical shifts was established for trans-[PtCl2(dma)PPh3] (dma = dimethylamine) complexes. The reliability of the computed values is shown to be critically dependent on the level of relativistic effects (two-component vs. four component), choice of density functionals, dynamical averaging, and solvation effects. Snapshots obtained from ab initio molecular dynamics simulations were used to identify those solvent molecules which show the largest interactions with the platinum complex, through inspection by using the non-covalent interaction program. We observe satisfactory accuracy from the full four-component matrix Dirac-Kohn-Sham method (mDKS) based on the Dirac-Coulomb Hamiltonian, in conjunction with the KT2 density functional, and dynamical averaging with explicit solvent molecules.

195Pt NMR and Molecular Dynamics Simulation Study of the Solvation of [PtCl6]2- in Water-Methanol and Water-Dimethoxyethane Binary Mixtures

The experimental ¹⁹⁵Pt NMR chemical shift, δ(¹⁹⁵Pt), of the [PtCl₆]^2- anion dissolved in binary mixtures of water and a fully miscible organic solvent is extremely sensitive to the composition of the mixture at room temperature. Significantly nonlinear δ(¹⁹⁵Pt) trends as a function of solvent composition are observed in mixtures of water-methanol, or ethylene glycol, 2-methoxyethanol, and 1,2-dimethoxyethane (DME). The extent of the deviation from linearity of the δ(¹⁹⁵Pt) trend depends strongly on the nature of the organic component in these solutions, which broadly suggests preferential solvation of the [PtCl₆]^2- anion by the organic molecule. This simplistic interpretation is based on an accepted view pertaining to monovalent cations in similar binary solvent mixtures. To elucidate these phenomena in detail, classical molecular dynamics computer simulations were performed for [PtCl₆]^2- in water-methanol and water-DME mixtures using the anionic charge scaling approach to account for the effect of electronic dielectric screening. Our simulations suggest that the simplistic model of preferential solvation of [PtCl₆]^2- by the organic component as inferred from nonlinear δ(¹⁹⁵Pt) trends is not entirely accurate, particularly for water-DME mixtures. The δ(¹⁹⁵Pt) trend in these mixtures levels off for high DME mole fractions, which results from apparent preferential location of [PtCl₆]^2- anions at the borders of water-rich regions or clusters within these inherently micro-heterogeneous mixtures. By contrast in water-methanol mixtures, apparently less pronounced mixed solvent micro-heterogeneity is found, suggesting the experimental δ(¹⁹⁵Pt) trend is consistent with a more moderate preferential solvation of [PtCl₆]^2- anions. This finding underlines the important role of solvent-solvent interactions and micro-heterogeneity in determining the solvation environment of [PtCl₆]^2- anions in binary solvent mixtures, probed by highly sensitive ¹⁹⁵Pt NMR. The notion that preferential solvation of [PtCl₆]^2- results primarily from competing ion-solvent interactions as generally assumed for monatomic ions, may not be appropriate in general.

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and ¹⁷O NMR quadrupolar relaxation rates are investigated in the series of complexes UO₂²⁺, UO₂(CO₃)₃^4-, and (UO₂)₃(CO₃)₆^6-. Car-Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO₂²⁺ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl ¹⁷O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO₂(CO₃)₃^4- and (UO₂)₃(CO₃)₆^6- complexes show a solvent dependence of the terminal carbonate ¹⁷O and ¹³C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin-orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique ¹³C chemical shifts is correctly reproduced for (UO₂)₃(CO₃)₆^6-. NMR relaxation rate data support large ¹⁷O peak widths, but remain below those noted in the experimental literature. Comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

Prediction of 195 Pt NMR of photoactivable diazido- and azine-Pt(IV) anticancer agents by DFT computational protocols

¹⁹⁵ Pt NMR chemical shifts for a series of large-sized photoactivable anticancer diazido-Pt(IV), homopiperizine-Pt(IV) and multifunctional azine-Pt(IV) complexes hardly to be probed experimentally and by sophisticated four-component and two-component relativistic calculations are predicted with high accuracy by density functional theory computational protocols. The calculated ¹⁹⁵ Pt NMR chemical shifts constitute a crucial descriptor for making highly predictive one-parameter quantitative structure activity relationships models that help in designing photoactivable Pt(IV)-based antitumor agents with high cytotoxicity and selectivity. Copyright © 2016 John Wiley & Sons, Ltd.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.

3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One

Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.

Alternative role of cisplatin in DNA damage – theoretical studies on the influence of excess electrons on the cisplatin–DNA complex

Interaction of excess electrons with cisplatin–DNA generates highly reactive Pt-containing species towards C–H abstraction, which strongly contributes to DNA damage.

Structure, solvent, and relativistic effects on the NMR chemical shifts in square-planar transition-metal complexes: assessment of DFT approaches

The role of various factors (geometry, solvent, relativistic treatment, functional) influencing the relativistic NMR chemical shift calculations for square-planar transition-metal complexes.

Accurate prediction of 195Pt NMR chemical shifts for a series of Pt(ii) and Pt(iv) antitumor agents by a non-relativistic DFT computational protocol

Exhaustive benchmark DFT calculations reveal that the non-relativistic GIAO-PBE0/SARC-ZORA(Pt)∪6-31+G(d)(E) computational protocol predicts accurate ¹⁹⁵Pt NMR chemical shifts for a wide range of square planar Pt(ii) and octahedral Pt(iv) anticancer agents.

Quadrupolar NMR Spin Relaxation Calculated Using Ab Initio Molecular Dynamics: Group 1 and Group 17 Ions in Aqueous Solution

Electric field gradient (EFG) fluctuations for the monoatomic ions (7)Li(+), (23)Na(+), (35)Cl(-), (81)Br(-), and (127)I(-) in aqueous solution are studied using Car-Parrinello ab initio molecular dynamics (aiMD) simulations based on density functional theory. EFG calculations are typically performed with 1024 ion-solvent configurations from the aiMD simulation, using the Zeroth Order Regular Approximation (ZORA) relativistic Hamiltonian. Autocorrelation functions for the spherical EFG tensor elements are computed, transformed into the corresponding spectral densities (under the extreme narrowing condition), and subsequently converted into NMR quadrupolar relaxation rates for the ions. The relaxation rates are compared with experimental data. The order of magnitude is correctly predicted by the simulations. The computational protocol is tested in detail for (81)Br(-).

Probing isotope shifts in 103Rh and 195Pt NMR spectra with density functional theory

Zero-point vibrationally averaged (rg(0)) structures were computed at the PBE0/SDD/6-31G* level for the [Pt(35)Cln(37)Cl5-n(H2(18)O)](-) (n = 0-5), cis-Pt(35)Cln(37)Cl4-n(H2(18)O)(H2(16)O) (n = 0-4), fac-[Pt(35)Cln(37)Cl3-n(H2(18)O)(H2(16)O)2](+) (n = 0-3), [Pt(35)Cln(37)Cl5-n((16/18)OH)](2-) (n = 0-5), cis-[Pt(35)Cln(37)Cl4-n((16/18)OH)2](2-) (n = 0-4), fac-[Pt(35)Cln(37)Cl3-n((16/18)OH)3](2-) (n = 0-3), cis-[Pt(35)Cln(37)Cl2-n((16/18)OH)4](2-) (n = 0-2), [Pt(35)Cln(37)Cl1-n((16/18)OH)5](2-) (n = 0-1), [Rh(35)Cln(37)Cl5-n(H2O)](2-) (n = 0-5), cis-[Rh(35)Cln(37)Cl4-n(H2O)2](-) (n = 0-4), and fac-Rh(35)Cln(37)Cl3-n(H2O)3 (n = 0-3) isotopologues and isotopomers. Magnetic shielding constants, computed at the ZORA-SO/PW91/QZ4P/TZ2P level, were used to evaluate the corresponding (35/37)Cl isotope shifts on the (195)Pt and (103)Rh NMR spectra, which are known experimentally. While the observed effects are reproduced reasonably well computationally in terms of qualitative trends and the overall order of magnitude (ca. 1 ppm), quantitative agreement with experiment is not yet achieved. Only small changes in M-Cl and M-O bonds upon isotopic substitution, on the order of femtometers, are necessary to produce the observed isotope shifts.

Mechanism of the cis-[Pt(1R,2R-DACH)(H2O)2]2+ intrastrand binding to the double-stranded (pGpG)·(CpC) dinucleotide in aqueous solution: a computational DFT study

A mechanism of the intrastrand 1,2-cross-link formation between the double-stranded pGpG·CpC dinucleotide (ds(pGpG)) and fully aquated oxaliplatin cis-[Pt(DACH)(H2O)2](2+) (DACH = cyclohexane-1R,2R-diamine) is presented. All structures of the reaction pathways including the transition states (TSs) were fully optimized in water solvent using DFT methodology with dispersion corrections. Both 5' → 3' and 3' → 5' binding directions were considered. In the first step there is a slight kinetic preference for 5'-guanine (5'G) monoadduct formation with an activation Gibbs free energy of 18.7 kcal/mol since the N7 center of the 5'G base is fully exposed to the solvent. On the other hand, the N7 atom of 3'-guanine (3'G) is sterically shielded by 5'G. The lowest energy path for formation of the 3'G monoadduct with an activation barrier of 19.3 kcal/mol is connected with a disruption of the 'DNA-like' structure of ds(pGpG). Monoadduct formation is the rate-determining process. The second step, chelate formation, is kinetically preferred in the 3' → 5' direction. The whole process of the platination is exergonic by up to -18.8 kcal/mol. Structural changes of ds(pGpG), charge transfer effects, and the influence of platination on the G·C base pair interaction strengths are also discussed in detail.

Computational study and molecular orbital analysis of NMR shielding, spin-spin coupling, and electric field gradients of azido platinum complexes

(195)Pt, (14)N, and (15)N NMR data for five azido (N(3)(-)) complexes are studied using relativistic density functional theory (DFT). Good agreement with experiment is obtained for Pt and N chemical shifts as well as Pt-N J-coupling constants. Calculated (14)N electric field gradients (EFGs) reflect experimentally observed trends for the line broadening of azido (14)N NMR signals. A localized molecular orbital analysis of the nitrogen EFGs and chemical shifts is performed to explain some interesting trends seen experimentally and in the first-principles calculations: (i) (14)N NMR signals for the Pt-coordinating (N(α)) nuclei in the azido ligands are much broader than for the central (N(β)) or terminal (N(γ)) atoms. The N(β) signals are particularly narrow; (ii) compared to N(γ), the N(α) nuclei are particularly strongly shielded; (iii) N(β) nuclei have much larger chemical shifts than N(α) and N(γ) ; and (iv) The Pt-N(α) J-coupling constants are small in magnitude when considering the formal sp hybridization of N(α). It is found that for N(α) a significant shielding reduction due to formation of the dative N(α)-Pt bond is counterbalanced by an increased shielding from spin-orbit (SO) coupling originating at Pt. Upon coordination, the strongly delocalized π system of free azide localizes somewhat on N(β) and N(γ). This effect, along with rehybridization at N(α) upon bond formation with Pt, is shown to cause a deshielding of N(γ) relative to N(α) and a strong increase of the EFG at N(α). The large 2p character of the azide σ bonds is responsible for the particularly high N(β) chemical shifts. The nitrogen s-character of the Pt-N(α) bond is low, which is the reason for the small J-coupling. Similar bonding situations are likely to be found in azide complexes with other transition metals.

Probing platinum azido complexes by 14N and 15N NMR spectroscopy

Metal azido complexes are of general interest due to their high energetic properties, and platinum azido complexes in particular because of their potential as photoactivatable anticancer prodrugs. However, azido ligands are difficult to probe by NMR spectroscopy due to the quadrupolar nature of (14)N and the lack of scalar (1)H coupling to enhance the sensitivity of the less abundant (15)N by using polarisation transfer. In this work, we report (14)N and (15)N NMR spectroscopic studies of cis,trans,cis-[Pt(N(3))(2)(OH)(2)(NH(3))] (1) and trans,trans,trans-[Pt(N(3))(2)(OH)(2)(X)(Y)], where X=Y=NH(3) (2); X=NH(3), Y=py (3) (py=pyridine); X=Y=py (4); and selected Pt(II) precursors. These studies provide the first (15)N NMR data for azido groups in coordination complexes. We discuss one- and three-bond J((15)N,(195)Pt) couplings for azido and am(m)ine ligands. The (14)N(α) (coordinated azido nitrogen) signal in the Pt(IV) azido complexes is extremely broad (W(1/2)≈2124 Hz for 4) in comparison to other metal azido complexes, attributable to a highly asymmetrical electric field gradient at the (14)N(α) atom. Through the use of anti-ringing pulse sequences, the (14)N NMR spectra, which show resolution of the broad (14)N(α) peak, were obtained rapidly (e.g., 1.5 h for 10 mM 4). The linewidths of the (14)N(α) signals correlate with the viscosity of the solvent. For (15) N-enriched samples, it is possible to detect azido (15)N resonances directly, which will allow photoreactions to be followed by 1D (15)N NMR spectroscopy. The T(1) relaxation times for 3 and 4 were in the range 5.7-120 s for (15)N, and 0.9-11.3 ms for (14)N. Analysis of the (1)J((15)N,(195)Pt) coupling constants suggests that an azido ligand has a moderately strong trans influence in octahedral Pt(IV) complexes, within the series 2-pic<py<NH(3)<Cl(-)<N(3)(-)<NO(2)(-)<SCN(-) (2-pic=2-methylpyridine). In addition, an axial Cl(-) appears to weaken an equatorial Pt(IV)-NH(3) bond to a greater extent than an axial OH(-) ligand.

Multinuclear solid-state NMR of square-planar platinum complexes — Cisplatin and related systems

Multinuclear solid-state nuclear magnetic resonance (SSNMR) experiments have been performed on cisplatin and four related square-planar compounds. The wideband uniform rate smooth truncation – Carr–Purcell–Meiboom–Gill (WURST–CPMG) pulse sequence was utilized in NMR experiments to acquire 195Pt, 14N, and 35Cl ultra-wideline NMR spectra of high quality. Standard Hahn-echo and magic-angle spinning 195Pt NMR experiments are also performed to refine extracted chemical shielding (CS) tensor parameters. Platinum magnetic shielding (MS) tensor orientations are calculated using both plane-wave density functional theory (DFT) and standard DFT methods. The tensor orientations are shown to be highly constrained by molecular symmetry elements, but also influenced to some degree by intermolecular interactions. 14N WURST–CPMG experiments were performed on three compounds and electric field gradient (EFG) parameters (the quadrupolar coupling constant, CQ, and the asymmetry parameter, ηQ) are reported. First principles calculations of the 14N EFG tensor parameters and orientations and affirm their dependence on the local hydrogen bonding environment. 35Cl WURST–CPMG experiments on cisplatin and transplatin are reported, using two different static magnetic fields to extract EFG and CS tensor parameters, and 35Cl EFG tensor magnitudes and orientations are predicted using first principles calculations. Transverse (T2) relaxation data for all nuclei are used to investigate heteronuclear dipolar relaxation mechanisms, as well as the nature of the local hydrogen bonding environments.