Asymmetric Environmental Effects on the Structure and Vibrations of

Aug 21, 2014 - Early efforts on the structural and electronic properties of cisplatin were contributed by Carloni et al. using density functional theo...
0 downloads 8 Views 2MB Size
Subscriber access provided by NATIONAL SUN YAT SEN UNIV

Article

Asymmetric Environmental Effects on the Structure and Vibrations of cis-[Pt(NH)Cl] in Condensed Phases 3

2

2

Chao Zhang, Emmanuel Baribefe Naziga, and Leonardo Guidoni J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/jp500865v • Publication Date (Web): 21 Aug 2014 Downloaded from http://pubs.acs.org on August 24, 2014

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Asymmetric Environmental Effects on the Structure and Vibrations of cis-[Pt(NH3)2Cl2] in Condensed Phases Chao Zhang,∗,†,¶ Emmanuel Baribefe Naziga,†,§ and Leonardo Guidoni∗,†,‡

1

Physics Department, Sapienza-Universita di Roma, P. le A. Moro 5, 00185, Rome, Italy, and Department of Physical and Chemical Sciences, University of L’Aquila, Via Vetoio, 67100, L’Aquila, Italy E-mail: [email protected]; [email protected]



To whom correspondence should be addressed Physics Department, Sapienza-Universita di Roma, P. le A. Moro 5, 00185, Rome, Italy ‡ Department of Physical and Chemical Sciences, University of L’Aquila, Via Vetoio, 67100, L’Aquila, Italy ¶ Current address: Institute of Physical Chemistry and Center for Computational Sciences, Johannes Gutenberg University Mainz, Staudinger Weg 7, D-55128, Mainz, Germamy § Current address: Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive, AB T1K 3M4, Lethbridge, Canada †

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract

2

3

We reported here the structural and vibrational properties of anti-cancer drug cis-

4

platin (cis-[Pt(NH3 )2 Cl2 ]) in gas phase, in solid phase and in aqueous solution using

5

DFT calculations, QM/MM molecular dynamics and effective normal modes analysis.

6

In contrast with the gas phase case, asymmetric hydrogen bonding environments are

7

found in both solid phase and aqueous solution. It is shown that the discrepancy of

8

the molecular geometry between previous gas phase calculations and the X-ray crystal

9

structure can be resolved by considering intermolecular hydrogen bonds in the calcu-

10

lations of solid phase. In addition, our simulations in solid phase and aqueous solution

11

reveal that asymmetric environmental effects lead to several spectral features observed

12

in experiments, such as the blue-shift in the N-H stretching region and the frequency

13

splitting of NH3 symmetric deformation modes. Furthermore, a similar decoupling and

14

localization of several vibrational modes of cisplatin is found in solid phase and aqueous

15

solution, in comparison to those of O−H stretching modes of water molecules in liquid

16

water (J. Phys. Chem. Lett., 2013, 4(19), pp 3245-3250).

17

KEYWORDS: Ab initio molecular dynamcis, vibrational spectroscopy, biomolecular sol-

18

vation, water dynamics.

19

Introduction

20

Cisplatin (cis-diamminedichloroplatinum(II) or cis-DDP) is a widely used anti-tumoral agent

21

in the therapy of various cancers since the discovery by Rosenberg and co-workers in 1960s. 1,2

22

Its structure consists of a square-planar central platinum(II) binding to two chlorine ligands

23

and two ammonia ligands in the cis-conformer. During the action, cisplatin molecule replace

24

the two chlorine ions with a pair of purine nitrogen atoms of DNA bases, inducing damages

25

and triggering cell death. 3,4

26

To understand better the mode of action, vibrational spectroscopies i.e. infra-red (IR)

27

and Raman, have been used to detect the intermediates, to measure the kinetics in drug 2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

28

delivery and to identify the binding modes. 5 One particular difficulty of applying vibrational

29

spectroscopy in biochemical systems is the interpretation of the experimental spectra. 6 To

30

this end, computational methods could offer great assistances. 7–14

31

Structural and vibrational properties of cisplatin have been addressed in several compu-

32

tational studies. Early efforts on the structural and electronic properties of cisplatin were

33

contributed by Carloni et al. using density functional theory (DFT) . 15 Later on, the bind-

34

ing of cisplatin to DNA was investigated by Quantum Mechanical/Molecular Mechanical

35

(QM/MM) Car-Parrinello molecular dynamics (MD). 16 Other computational works also in-

36

vestigated the hydrolysis process of cisplatin 17 and the structure of its hydration shell. 18

37

From the point of view of an isolated molecule, a comprehensive ab initio quantum chemical

38

analysis was given by Pavankumar et al. 19 A systematic evaluation of the performance of

39

different exchange-correlation functionals for cisplatin by Michalska and Wysokiński 20 con-

40

cluded that the MPW1PW/LanLDZ combination gave the best agreement with experiment

41

for structural and vibrational properties. Further discussions on the choice of basis sets and

42

effective core potentials were also reported by several groups. 21–23 Assessments on methods

43

for calculating Raman intensities was reported 24 and all-electron basis sets calculations in

44

gas phase were compared with newly measured Raman spectrum in solid phase. 25

45

So far, all reported vibrational studies of cisplatin are restricted to gas phase calculations,

46

although experimental spectra were collected in solid phase 25–29 and in aqueous solution. 30,31

47

Conceivably, the missing of environmental effects could be an important source for the ob-

48

served discrepancy between theoretical calculations and experimental measurements. 19,25

49

Moreover, it would be desirable to establish the structural and vibrational properties of

50

cisplatin in aqueous solution, where it acts as an anti-cancer drug. 3,4

51

In this work, we focused on structural and vibrational properties of cisplatin in condensed-

52

phase systems, i.e. solid phase and aqueous solution. The latter was treated by mixed

53

QM/MM Car-Parrinello MD, which includes explicit solvent and possible anharmonicity

54

effects at finite temperature. 32 This approach has been successfully used to calculate IR and

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

55

Raman spectra of (bio)molecules aqueous solutions. 7,8,11,12,33–39 In addition, we employed

56

the recently proposed effective normal modes analysis 11,40 to elucidate the finite temperature

57

picture of vibrational motions of cisplatin and to interpret corresponding vibrational spectra.

58

By comparing the condensed phase results with that of gas phase, we showed that both in

59

solid phase and aqueous solution share a similar asymmetric hydrogen bonding environment.

60

In solid state, such asymmetric environment originated from the asymmetric interactions

61

of cisplatin in the molecular crystal due to packing effects. On the contrary, in aqueous

62

solution, it is instead the hydrogen bonds (HBs) dynamics of the surrounding water molecules

63

which instantaneously create an asymmetric field. These effects lead to the blue-shift in the

64

N-H stretching region and the frequency splittings of NH3 symmetric deformation modes.

65

Subsequently, they provide an interpretation of the corresponding experimental spectra in

66

solid phase and aqueous solution.

67

This article is organized as follows: Section 2 provides the theoretical background and

68

computational protocols on the calculations of the IR and Raman spectra and vibrational

69

analyses at 0K and finite temperature; Section 3 presents our results and discussions on

70

asymmetric hydrogen bonding environments in solid phase and aqueous solution and their

71

effects on structural and vibrational features of cisplatin molecule as observed in experiments.

72

Conclusions and perspectives are given in Section 4.

73

Computational methods

74

Electronic structure calculations and QM/MM simulations

75

Cisplatin in gas phase. Calculations using localized basis sets were performed using Den-

76

sity Functional Theory (DFT) with Becke-Perdew (BP) exchange correlation functional. 41,42

77

We used Gaussian package 43 with the LanL2DZ basis set according to the previous work

78

from Wysokiński and coworkers. 20 In the Lan2DZ basis set, valence basis set and relativistic

79

effective core potential (ECP) were used for the platinum atom. 44 The IR absorption coeffi4

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

80

cient and Raman activity were subsequently obtained from these calculations. In plane wave

81

calculations we used the CPMD program. 45 A cubic box of size 25 a.u. and a plane wave

82

energy cut-off of 70 Ry were used. The BP exchange-correlation functional was employed in

83

combination with Troullier-Martins pseudopotentials. 46 The system was isolated using the

84

decoupling scheme by Martyna and Tuckerman. 47

85

Cisplatin in solid phase. Geometry optimizations of cisplatin in solid phase were ini-

86

tialized based on the alpha form of the crystal structure for non-hydrogen atoms. 48 Through-

87

out the context of this work, we discussed only the structural and vibrational properties of

88

alpha form which is energetically more stable than the beta form. Reader interested in

89

polymorphism in cisplatin crystals could refer to recent literatures on this issue. 49,50 All cal-

90

culations were performed using CPMD program 45 with BP exchange correlation functional

91

and plane wave cutoff of 70 Ry. K-point meshes for Brillouin zone sampling were constructed

92

using the Monkhorst-Pack 51 scheme with 2x2x2 and 4x4x4 grids. The inversion center of

93

two cisplatin molecules in the unit cell was kept during geometry optimizations.

94

Cisplatin in aqueous solution. QM/MM simulations of cisplatin in aqueous solu-

95

tion were performed using the CPMD/Gromos interface. 45,52 A single cisplatin molecule was

96

considered as the quantum system surrounded by 1606 TIP3P classical waters. 53 The sys-

97

tem was firstly equilibrated by classical MD simulations at constant pressure (1 atm.) and

98

temperature (300 K) using Nosé-Hoover thermostat 54 and Parrinello-Rahman barostat. 55

99

Parameterization of cisplatin in the classical MD simulations followed the recommended

100

AMBER procedure. 56 Most of the bond, angle and torsional parameters were adapted from

101

classical MD simulations of cisplatin bound to DNA. 57 For bond and angle parameters of

102

NH3 groups, standard AMBER values were adopted. 58

103

The last frame of classical MD simulations was taken for subsequent QM/MM simula-

104

tions. In QM/MM Car-Parrinello dynamics, we set the fictitious mass to 400 a.u and the

105

time step to 0.085 fs. Other settings for QM region were the same as for the gas phase

106

calculations. The QM/MM system was first equilibrated for 2 ps with Nosé-Hoover thermo-

5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 30

107

stat at 300 K. Then, the production run continued for 40 ps in NVE ensemble. The dipole

108

moment of cisplatin was sampled along the simulation every two time steps. Snapshots from

109

the trajectory were extracted every 10 time steps for analysis. For each snapshot, the polar-

110

izability tensor of the cisplatin within the field generated by the classical point charges, was

111

calculated using Gaussian 43 package, following the protocol described in ref. 38 For the pur-

112

pose of checking the convergence of the sampling, the QM/MM simulation was prolonged

113

for additional 40 ps. There are no significant differences in the shape of the VDOS (see

114

Fig. S1 in Supporting Information), although details of the intensity depend on the length

115

of samplings. This is in accord with the fact that 1 cm−1 spectral resolution in the Fourier

116

transform of the time correlation function needs only 30 ps MD simulations. 32

117

Calculation of IR spectra

118

IR absorption coefficient for ith normal vibration in gas phase at zero temperature is given

119

as 59 :

IiIR 120

121

122

123

=C



dM dQi

2

(1)

where C is a constant, M is the electric dipole moment of the system and Qi are the normal coordinates. IR absorption spectrum at finite temperature in aqueous solution can be calculated from the fluctuation-dissipation theorem: 60

I

IR

(ω) = C

Z



dte−iωt hM(t)M(0)i

(2)

−∞

124

The final line shape should be multiplied by ω 2 /kT as a quantum correction 61,62 to the

125

classical time-correlation function.

6

ACS Paragon Plus Environment

Page 7 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

126

Calculation of Raman spectra

127

For a gas phase molecule at zero temperature, the classical expression of Raman scattering

128

cross section of ith normal vibration over a solid angle dΩ is: 60 dσi (ω0 − ωi )4 =C dΩ ωi Bi



dα dQi

2

(3)

129

where (ω0 − ωi )4 is the density of states and Bi = 1 − e−¯hωi /kT is a temperature factor

130

which accounts for the intensity contribution of excited vibrational states for Raman Stokes

131

intensity, 63 and α is the polarizability tensor.

132

133

The corresponding Raman spectrum at finite temperature in aqueous solution can be obtained in a similar way to IR spectrum: 60,64 dσ = C(ω0 − ω)4 dΩ

134

135

Z

∞ −∞

dte−iωt [ǫ0 · α(t) · ǫ][ǫ0 · α(0) · ǫ]

(4)

where ǫ0 and ǫ are the unit vectors in which the electric vector of incident light and of scattered light are polarized along.

136

For spatially isotropic systems (liquids and gases) in the plan-polarization condition, i.e.

137

the direction of the incident beam, the polarization direction of this beam, and the direction

138

of observation are perpendicular to each other, the cross section can be rewritten as: dσ = C(ω0 − ω)4 dΩ

139

140

141

142

Z



dte−iωt hT rβ(0)β(t)i

(5)

−∞

where β is a traceless anisotropic part of the polarizability tensor α, β = α − α ¯ I (¯ α is the average of the polarizability tensor trace and I is the unit tensor). The final line shape should  hω ¯ hω ¯ coth 2kT as a quantum correction 38,65 to the classical time-correlation be multiplied by 2kT function.

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

143

Page 8 of 30

In practice, hT rβ(0)β(t)i can be calculated as a sum of scalar auto-correlation functions.

hT rβ(0)β(t)i = hβxx (0)βxx (t) + βyy (0)βyy (t) + βzz (0)βzz (t)i (6)

+2 (hβyx (0)βyz (t) + βzx (0)βzx (t) + βzy (0)βzy (t)i)

144

Vibrational analysis

145

Vibrational normals modes of a given system at zero temperature can be found by solving

146

the eigenvalue problem: 3N X

(7)

(hij − λk δij )νik = 0

i

1/2

147

where hij is the mass-weighted molecular Hessian, λk

148

the normal modes.

are the normal frequencies and νk

149

At finite temperature, vibrational analysis is not straightforward. Here we used the re-

150

cently proposed effective normal modes (ENM) analysis. 11,40 In this method, effective normal

151

modes νk can be extracted from its vibrational density of states (VDOSs) by minimizing the

152

following functional:

Ω(n) =

X k

153

154

1 2πkT

Z

dω|ω|2n P ν˙ k (ω) −



1 2πkT

Z

dω|ω|n P ν˙ k (ω)

2 !

(8)

with respect to νk . Here n is an integer constant and P ν˙ k is VDOS of νk . When n=2, the method equals to the standard normal modes analysis with an average Hessian matrix.

155

In practice, it is common to use internal coordinates such as bond lengths, bending

156

angles and dihedral angles, i.e. Decius coordinates, 59 instead of cartesian coordinates for

157

molecules. Furthermore, it is convenient to define another set of internal coordinates, i.e.

158

Pulay coordinates 66 as linear combinations of Decius coordinates. The advantage using

159

Pulay coordinates is that they reflect the the symmetry of the molecular functional groups

160

and resemble the actually localized vibrational modes for complicated cases.

161

For cisplatin molecule, we used C2v type Pulay coordinates for Pt square planar structure 8

ACS Paragon Plus Environment

Page 9 of 30

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 30

180

about 40◦ (N1) and 20◦ (N2) with respect to the gas phase structure (Fig. 2 inset) and

181

form asymmetric intermolecular HBs with neighbouring molecules (Fig. 2). We noted that

182

although the detailed orientation and the strength of HBs may depend on the exchange-

183

correlation functional used in the DFT calculations as well as by the difference of lattice

184

constants reported in the low- and the high-resolution X-ray crystal structures, the observed

185

asymmetric hydrogen bonding environment in solid phase is a general feature. This is also

186

supported by the high-resolution X-ray crystal structure of cisplatin where positions of hy-

187

drogen atoms were resolved from neutron powder diffraction. 49 Table 1: Comparison of the calculated and experimental bond lengths (Å) and bond angles(◦ ) of cisplatin.

r(P t − N 1) r(P t − N 2) r(P t − Cl1) r(P t − Cl2) 6 (N 1P tCl1) 6 (N 2P tCl2) 6 (N 1P tN 2) 6 (Cl1P tCl2)

gas 2.08 v.s. 2.31 v.s 82.3 v.s 99.7 95.8

Calc. solid 2.038 2.043 2.362 2.363 89.1 87.9 91.3 91.6

aqueous 2.04(4) v.s. 2.36(5) v.s 91(4) v.s. 85(3) 91(3)

Exp. 48 solid 2.05(4) 1.95(4) 2.328(9) 2.333(9) 92(1) 88.5(0.9) 87(2) 91.9(0.4)

Exp. 49 solid 2.049(3) 2.047(3) 2.3206(3) 2.3216(8) 89.18(9) 88.55(9) 90.62(12) 91.65(3)

188

For the vibrational properties, the full list of 27 normal frequencies of cisplatin in solid

189

phase and gas phase is reported in Table 2. We were aware that the cutoff for plane wave

190

basis sets and the choice of the exchange-correlation functional are important factors for

191

an accurate determination of vibrational frequencies. For instance, it is known that the

192

inclusion of Hartree-Fock exchange leads to a significant improvement on the calculated

193

vibrational spectra. 67 Regarding to this issue, we found that adopting a cutoff at 100 Ry

194

in gas phase calculations does not change the normal modes themselves compared to the

195

results using a cutoff at 70 Ry but slightly shifts the frequencies with no general tendency

196

(Table 2). As shown in Table 2, the calculated normal frequencies in solid phase are in good

197

agreement with the ones resolved from Raman spectroscopy under the same conditions. 25

198

This suggests that i) the level of theory used here is sufficient for describing the vibrational 10

ACS Paragon Plus Environment

Page 11 of 30

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

199

properties of cisplatin molecule; ii) The anharmonicity, which usually leads to the red-shift

200

of the frequency, plays a less significant role in the case of cisplatin molecule, in particular

201

for high wavenumber modes. Moreover, throughout this work, choices of cutoff and the

202

exchange-correlation functional were the same in gas phase, solid phase and aqueous solution

203

calculations. The following results and conclusions will mainly concern spectral differences

204

between gas phase and condensed phases, therefore we expect a further benefit due to error

205

cancellation.

206

In solid phase, it is found that N−H stretching modes of cisplatin molecule (ν22 to ν27

207

in Table 2) are generally red-shifted with respect to gas phase. Since both calculations were

208

done with the harmonic approximation, forming of intermolecular HBs in solid phase is the

209

only reason for this red-shift. In addition, we found that stretching modes ν22 ↔ ν23 , ν24 ↔

210

ν25 ν26 ↔ ν27 become non-degenerate in frequency in contrast to that in gas phase. Certainly,

211

this comes from the fact that two ammonia groups form asymmetric number/strength of HBs

212

and the gas phase C2v symmetry of cisplatin molecules is broken in solid phase.

213

In terms of normal modes, it is found that Pt−N stretching modes (ν10 and ν11 ) are

214

decoupled and localized to individual Pt-N bonds, at variance with what was obtained in

215

gas phase, in which they distribute equally to each bond (Fig. 3). A similar pattern is found

216

also in Pt−Cl stretching modes (ν8 and ν9 ). Interestingly, the decoupling and localization of

217

Pt−N and Pt−Cl stretching modes of cisplatin in solid phase resembles that of O−H stretch-

218

ing modes of water molecules in water dimer and liquid water. 68 As revealed previously, this

219

phenomenon comes from the broken symmetry of the target molecules in asymmetric hydro-

220

gen bonding environment. 68 Moreover, it is worth to note that the degree of the decoupling

221

and localization indeed provides a way to gauge the details of HBs asymmetry.

222

In literature, the low frequency modes (ν1 to ν7 ) of cisplatin in solid phase were seldom

223

addressed, 19,20,24,25 probably because they are mostly Raman or IR inactive. It is found

224

that: 1) Modes δCl−P t−Cl and δN −P t−N are both blue-shifted (150 cm−1 → 161 cm−1 and 255

225

cm−1 → 290 cm−1 respectively). This is in agreement with the fact that angles of Cl−Pt−Cl

12

ACS Paragon Plus Environment

Page 12 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Gas Solid Aqueous a 25 Mode # Fundamental IR Raman Sym. BP/PW BPW91/LanLDZ BP/PW Exp. Pulay (n=2) 1 ρ′ a2 108(113) 111 186 163 2 δCl−P t−Cl Y a1 150(149) 135 161 162 155 3 ρ Y b1 151(165) 154 217 171 4 τN′ H3 Y b1 154(154) 157 263* 143* 5 τN H3 a2 136(148) 159 246* 140* 6 δN −P t−Cl Y b2 233(235) 228 224 210 184 7 δN −P t−N Y Y a1 255(249) 227 290 255 250 ′ 8 νP t−Cl Y Y b2 333(334) 318 304* 317 319* 9 νP t−Cl Y Y a1 345(345) 328 309* 323 348* Y Y b2 442(450) 469 512* 508 538* 10 νP′ t−N 11 νP t−N Y Y a1 479(480) 472 542* 524 550* ′ 12 θN H3 Y a2 690(701) 759 800* 724 795* 13 θN H3 Y Y b1 707(718) 777 830* 789 802* ′ 14 θN Y Y b 741(746) 811 838* 811 808* 2 H3 15 θN H3 Y Y a1 761(767) 843 865* 824 817* 16 δ N H3 Y Y b2 1195(1206) 1252 1272* 1295 1299* 17 δ N H3 Y Y a1 1202(1217) 1255 1318* 1316 1338* ′ 18 δ N H3 Y Y b2 1565(1576) 1631 1558* 1537 1595* ′ 19 δN Y Y a 1567(1580) 1639 1567* — 1597* 1 H3 ′ Y a2 1588(1602) 1663 1598* 1601 1602* 20 δ N H3 ′ 21 δ N H3 Y Y b1 1591(1608) 1668 1620* 1648 1612* 22 νN −H Y Y a1 3200(3214) 3228 3198* 3211 3214*# ′ 23 νN −H Y Y b2 3201(3216) 3228 3217* — 3237*# 24 νN′ −H Y Y b2 3365(3387) 3441 3274* 3287 3284*# 25 νN −H Y Y a1 3366(3388) 3442 3287* — 3290*# ′ Y a2 3419(3444) 3518 3306* 3309 3318*# 26 νN −H 27 νN −H Y Y b1 3420(3445) 3519 3350* — 3321*# ′ ′ Note: ν and ν (bond symmetric and asymmetric stretch), δ (bend or symmetric deformation), δ (degenerative deformation), θ and θ′ (in-phase and out-of-phase rock), τ and τ ′ (in-phase and out-of-phase twist), ρ and ρ′ (skeletal in-phase and out-of-phase deformation).*: Decoupled and localized modes (see Text for discussions). #: These modes are subject to a spurious red-shift because of the fictitious mass used in Car-Parrinello MD (see Text for discussions). For gas phase, the symmetries of the modes are also shown, where there are 5 belonged to a2 symmetry which are IR inactive. All the modes were sorted ascendingly in frequency with respect to the gas phase BPW91/LanLDZ calculation. Number in parenthesis calculated with using and increased plane wave cutoff of 100 Ry.

Table 2: The comparisons of normal frequencies of cisplatin molecule in gas phase, crystal and aqueous solution.

Page 13 of 30 The Journal of Physical Chemistry

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 14 of 30

Page 15 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

240

The Journal of Physical Chemistry

wavenumber smaller than 2000 cm−1 (see Fig. S2 in Supporting Information).

241

Interestingly, we found that the first group of N-H stretching modes (ν22 and ν23 in

242

Table 2) have a blue-shift of 120 cm−1 in aqueous solution with respect to gas phase at 20K.

243

This leads to a reduction of about 70 cm−1 of the frequency gab between the first group

244

(ν22 and ν23 in Table 2) and the second group (ν24 and ν25 in Table 2) of N-H stretching

245

modes, with respect the gas phase value of 160 cm−1 . Such unexpected blue-shift has to

246

be associated with intermolecular HBs between chlorine and hydrogen atoms. It is found

247

that the Cl-H distance of cisplatin molecules in gas phase at 0K and 20K has two distinct

248

populations with ratio 1:2 (Fig. 4). Temperature effect alone help the rotation of NH3 group

249

and modulate the ratio between two populations. However, intramolecular HBs are strong

250

enough, therefore two populations remain distinct from each other. Only in the aqueous

251

solution, water molecules help to lower the rotational barrier of NH3 groups and to decrease

252

the strength of Cl-H HBs. It suggests that NH3 groups of cisplatin in aqueous solution

253

resemble the ammonia in gas phase with C3v symmetry. This consequently leads to the

254

observed blue-shift of the first group of N-H stretching modes (ν22 and ν23 ) and breaks the

255

gas phase C2v symmetry of cisplatin molecules in aqueous solution. It is worth to note

256

that these observations come from the VDOS of MD simulations and do not depend on the

257

procedure of the effective normal modes analysis.

258

Indeed, the blue-shift of the symmetric N-H stretching modes (ν22 and ν23 ) also appears

259

in solid phase. Experimental spectroscopy in solid phase reports that frequency gap between

260

the first group (ν22 and ν23 in Table 2) and the second group (ν24 and ν25 in Table 2) of

261

N-H stretching modes is reduced up to 76 cm−1 , 25 which is very close to the one calculated

262

in aqueous solution (70 cm−1 ). Considering that the calculated average bond distances and

263

angles of cisplatin in aqueous solution are very close to those reported in X-ray crystal struc-

264

ture (Table 1), it would not be surprising to see this spectral feature in both aqueous solution

265

and solid phase. Although a direct quantitative comparison cannot be done, we found the

266

calculated IR and Raman spectra of cisplatin in aqueous solution are in general agreement

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 16 of 30

Page 17 of 30

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 18 of 30

Page 19 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

294

neutron scattering study of cisplatin in solid phase assigned such splitting to be 28 cm−1 , 70

295

which can not be reproduced by gas phase calculations. On the contrary, our calculated IR

296

spectrum from time-correlation functional formalism fairly reproduces such splitting of NH3

297

symmetric deformation modes (Fig. 7).

298

It should be pointed out that in the long-time limit, the average Hessian and therefore

299

the normal modes of cisplatin in aqueous solution have C2v symmetry. However, in cases

300

where the environment is fluctuating between different conformations, the effective normal

301

modes analysis based on VDOS and the average Hessian analysis give different results. It has

302

to be considered that because of this fluctuating environment, the cisplatin in the condensed

303

phase does not have C2v symmetry anytime, although on average it has C2v symmetry.

304

This is indeed our case, since water is prevalently binding to one of the two NH3 groups

305

at a certain time. Indeed, when looking into the corresponding effective normal modes in

306

aqueous solution (Fig. 7 inset), it is found that the C2v symmetry is broken and modes are

307

decoupled and localized.

308

As demonstrated by the effective normal modes analysis in the case of liquid water, the

309

instantaneous asymmetric hydrogen bonding environment would lead to a frequency sepa-

310

ration of O-H stretching modes in a non-trivial manner. 68 Correspondingly, the observed

311

frequency splitting of NH3 symmetric deformation modes comes from a similar origin. Such

312

a phenomenon is not present in gas phase but is present in both our VDOS analysis and

313

in the experimental data. Therefore, our relatively short simulation catches the underlay-

314

ing physical picture and provides a rationale, although a quantitative determination of the

315

relationship between the degree of environmental asymmetry and the frequency separation

316

of cisplatin would require a large sample of QM/MM MD simulations, which is out of the

317

scope of the present work.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 20 of 30

Page 21 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

328

phase calculations and the X-ray crystal structure can be resolved by considering inter-

329

molecular HBs in the calculations of periodic crystal structure. In addition, our simulations

330

in condenses phase show that asymmetric environmental effects lead to the blue-shift in the

331

N-H stretching region and the frequency splitting of NH3 symmetric deformation modes.

332

Subsequently, it provides a physical explanation for these spectral features as observed in

333

experiments, which can not be revealed by gas phase calculations alone.

334

Furthermore, we noticed that the decoupling and localization of several vibrational modes

335

of cisplatin founded here in solid phase and aqueous solution, such as Pt-N and Pt-Cl stretch-

336

ing modes, resembles those of O−H stretching modes of water molecules in water dimer and

337

liquid water. 68,71 Because of the same gas phase C2v molecular symmetry for both cases,

338

we suggest that the decoupled and localized modes as well as the relationship between the

339

environmental asymmetry and the frequency separation could be a general characteristic for

340

this type of HBs forming molecules in condensed phases.

341

Acknowledgement

342

We gratefully acknowledge M. Martinez and D. Bovi for their assistance in the analysis of the

343

effective normal modes and W. Andreoni for providing us with the Platinum pseudopoten-

344

tial. C.Z. thanks for the interesting discussions on the instantaneous asymmetry with T. D.

345

Kühne and R. Z. Khaliullin. Computational resources were supplied by CASPUR, CINECA,

346

and the Caliban-HPC centre at the University of L’Aquila. C.Z. acknowledges the Gauss

347

Center for Supercomputing (GCS) for providing computing time through the John von Neu-

348

mann Institute for Computing (NIC) on the GCS share of the supercomputer JUQUEEN at

349

Jülich Supercomputing Center (JSC). L.G. acknowledges funding provided by the European

350

Research Council project n. 240624.

351

Supporting Information Available

352

Supporting information includes: A) The effects on C3v Pulay coordinates in effective normal 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 30

353

modes analysis of cisplatin. B) The convergence issue of VDOS of cisplatin in aqueous

354

solution; C) The effects of fictitious mass in Car-Parrinello MD on the vibrational spectra of

355

cisplatin; C) The residence time of water molecules around NH3 groups of cisplatin.

356

material is available free of charge via the Internet at http://pubs.acs.org/.

357

References

358

359

360

361

362

363

364

365

366

367

368

369

This

(1) Rosenberg, B.; Vancamp, L.; Krigas, T. Inhibition of Cell Division in Escherichia Coli by Electrolysis Products from a Platinum Electrode. Nature 1965, 205, 698–699. (2) Rosenberg, B.; Camp, L. V.; Trosko, E. J.; Mansour, V. H. Platinum Compounds: a New Class of Potent Antitumour Agents. Nature 1969, 222, 385–386. (3) Boulikas, T.; Vougiouka, M. Cisplatin and Platinum Drugs at the Molecular Level. Oncol. Rep. 2003, 10, 1663–1682. (4) Siddik, Z. H. Cisplatin: Mode of Cytotoxic Action and Molecular Basis of Resistance. Oncogene 2003, 22, 7265–79. (5) Wartewig, S.; Neubert, R. H. H. Pharmaceutical Applications of Mid-IR and Raman spectroscopy. Adv. Drug Deliv. Rev. 2005, 57, 1144–1170. (6) Movasaghi, Z.; Rehman, S.; Rehman, I. U. Raman Spectroscopy of Biological Tissues. Appl. Spectrosc. Rev. 2007, 42, 493–541.

370

(7) Schmitz, M.; Tavan, P. Vibrational Spectra from Atomic Fluctuations in Dynamics Sim-

371

ulations. II. Solvent-induced Frequency Fluctuations at Femtosecond Time Resolution.

372

J. Chem. Phys. 2004, 121, 12247–12258.

373

(8) Schmitz, M.; Tavan, P. Vibrational Spectra from Atomic Fluctuations in Dynamics

374

Simulations. I. Theory, Limitations, and a Sample Application. J. Chem. Phys. 2004,

375

121, 12233–12246. 22

ACS Paragon Plus Environment

Page 23 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

376

(9) Gaigeot, M.-P.; Martinez, M.; Vuilleumier, R. Infrared Spectroscopy in the Gas and

377

Liquid Phase from First Principle Molecular Dynamics Simulations: Application to

378

Small Peptides. Mol. Phys. 2007, 105, 2857–2878.

379

380

(10) Meier, R. J. Calculating the Vibrational Spectra of Molecules: an Introduction for Experimentalists with Contemporary Examples. Vib. Spectro. 2007, 43, 26–37.

381

(11) Bovi, D.; Mezzetti, A.; Vuilleumier, R.; Gaigeot, M.-P.; Chazallon, B.; Spezia, R.;

382

Guidoni, L. Environmental Effects on Vibrational Properties of Carotenoids: Experi-

383

ments and Calculations on Peridinin. Phys. Chem. Chem. Phys. 2011, 13, 20954–20964.

384

(12) Montagna, M.; Sterpone, F.; Guidoni, L. Structural and Spectroscopic Properties of

385

Water around Small Hydrophobic Solutes. J. Phys. Chem. B 2012, 116, 11695–11700.

386

(13) Mathias, G.; Baer, M. Generalized Normal Coordinates for the Vibrational Analysis of

387

Molecular Dynamics Simulations. J. Chem. Theory and Comput. 2011, 7, 2028–2039.

388

(14) VandeVondele, J.; Tröster, P.; Tavan, P.; Mathias, G. Vibrational Spectra of Phosphate

389

Ions in Aqueous Solution Probed by First-Principles Molecular Dynamics. J. Phys.

390

Chem. A 2012, 116, 2466–2474.

391

(15) Carloni, P.; Andreoni, W.; Hutter, J.; Curioni, A.; Giannozzi, P.; Parrinello, M. Struc-

392

ture and Bonding in Cisplatin and Other Pt(II) Complexes. Chem. Phys. Lett. 1995,

393

234, 50–56.

394

(16) Spiegel, K.; Rothlisberger, U.; Carloni, P. Cisplatin Binding to DNA Oligomers from

395

Hybrid Car-Parrinello/Molecular Dynamics Simulations. J. Phys. Chem. B 2004, 108,

396

2699–2707.

397

398

(17) Lau, J. K.-C.; Ensing, B. Hydrolysis of Cisplatin—a First-Principles Metadynamics Study. Phys. Chem. Chem. Phys. 2010, 12, 10348–10355.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 30

399

(18) Melchior, A.; Martínez, J.; Pappalardo, R.; Marcos, E. Hydration of Cisplatin Studied

400

by an Effective Ab Initio Pair Potential Including Solute-Solvent Polarization. J. Chem.

401

Theory and Comput. 2013, 9, 4562–4573.

402

(19) Pavankumar, P. N. V.; Seetharamulu, P.; Yao, S.; D.Saxe, J.; Reddy, D. G.;

403

Hausheer, F. H. Comprehensive Ab Initio Quantum Mechanical and Molecular Or-

404

bital(MO) Analysis of Cisplatin: Structure, Bonding, Charge Density, and Vibrational

405

Frequencies. J. Comp. Chem. 1999, 20(3), 365–382.

406

(20) Wysokiński, R.; Michalska, D. The Performance of Different Density Functional Meth-

407

ods in the calculation of Molecular Sturctures and Vibrational Spectra of Platinum (II)

408

Antitumor Drugs: Cisplatin and Carboplatin. J. Comp. Chem. 2001, 22, 901–912.

409

(21) Fiuza, S.; Amado, A.; Marques, M.; de Carvalho, L. B. Use of Effective Core Potential

410

Calculations for the Conformational and Vibrational Study of Platinum (II) Anticancer

411

Drugs. cis-Diamminedichloroplatinum (II) as a Case Study. J. Phys. Chem. A 2008,

412

112, 3253–3259.

413

(22) de Berredo, R.; Jorge, F. E. All-Electron Double Zeta Basis Sets for Platinum: Esti-

414

mating Scalar Relativistic Effects on Platinum(II) Anticancer Drugs. J. Mol. Struct.:

415

THEOCHEM 2010, 961, 107–112.

416

(23) Paschoal, D.; Marcial, B. L.; Lopes, J. F.; Almeida, W. B. D.; Santos, H. F. D. The

417

Role of the Basis Set and the Level of Quantum Mechanical Theory in the Prediction

418

of the Structure and Reactivity of Cisplatin. J. Comput. Chem. 2012, 33, 2292–2302.

419

420

421

(24) Michalska, D.; Wysokiński, R. The Prediction of Raman Spectra of Platinum (II) Anticancer Drugs by Density Functional Theory. Chem. Phys. Lett. 2005, 403, 211–217. (25) Amado, A. M.; Fiuza, S. M.; Marques, M. P. M.; de Carvalho, L. A. E. B.

422

Conformational and Vibrational Study of Platinum (II) Anticancer Drugs:

423

diamminedichloroplatinum (II) as a Case Study. J. Chem. Phys. 2007, 127, 185104. 24

ACS Paragon Plus Environment

Cis-

Page 25 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

424

(26) Adams, D. M.; Chatt, J.; Gerratt, J.; Westland, A. D. Far-infrared Spectra of Some

425

Square-Planar Complexes [PtX, L,] (X = C1 or Br). The Influence of L upon the

426

Platinum-Halogen Stretching Frequency and Its Relation to the Trans-Effect. J. Chem.

427

Soc. 1964, FEB, 734–739.

428

(27) Nakamoto, K.; McCarthy, P. J.; Fujita, J.; Condrate, R. A.; Behnke, G. T. Infrared

429

Studies of Ligand-Ligand Interaction in Dihalogenodiammineplatinum(II) Complexes.

430

Inorg. Chem. 1965, 4, 36–43.

431

432

(28) Pouchert, C. The Aldrich Library of FT-IR spectra; Aldrich Chem. Co.: Milwaukee, Wis., 1997.

433

(29) Tripisciano, C.; Costa, S.; Kalenczuk, R.; Borowiak-Palen, E. Cisplatin Filled Multi-

434

walled Carbon Nanotubes–a Novel Molecular Hybrid of Anticancer Drug Container.

435

Eur. Phys. J. B 2010, 75, 141–146.

436

437

438

439

(30) Yan, X.; Gemeinhart, R. A. Cisplatin Delivery from Poly(Acrylic Acid-Co-Methyl Methacrylate) Microparticles. J. Control. Release 2005, 106, 198–208. (31) Casolaro, M.; Cini, R.; Bello, B. D.; Ferrali, M.; Maellaro, E. Cisplatin/Hydrogel Complex in Cancer Therapy. Biomacromolecules 2009, 10, 944–949.

440

(32) Schmitz, M.; Tavan, P. In Modern Methods for Theoretical Physical Chemistry of

441

Biopolymers; Starikov, E. B., Lewis, J. P., Tanaka, S., Eds.; Elsevier: Amsterdam,

442

2006.

443

(33) Rousseau, R.; Kleinschmidt, V.; Schmitt, U. W.; Marx, D. Assigning Protonation Pat-

444

terns in Water Networks in Bacteriorhodopsin Based on Computed IR Spectra. Angew.

445

Chem. Int. Ed. 2004, 43, 4804–4807.

446

(34) Klähn, M.; Schlitter, J.; Gerwert, K. Theoretical IR Spectroscopy based on QM/MM

25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

447

Calculations Provides Changes in Charge Distribution, Bond Lengths, and Bond Angles

448

of the GTP Ligand Induced by the Ras-Protein. Biophys. J. 2005, 88, 3829–3844.

449

(35) Yang, S.; Cho, M. IR spectra of N-methylacetamide in Water Predicted by Com-

450

bined Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations.

451

J. Chem. Phys. 2005, 123, 134503.

452

(36) Kinnaman, C. S.; Cremeens, M. E.; Romesberg, F. E.; Corcelli, S. A. Infrared Line

453

Shape of an Alpha-Carbon Deuterium-Labeled Amino Acid. J. Am. Chem. Soc. 2006,

454

128, 13334–13335.

455

(37) Mroginski, M. A.; Mark, F.; Thiel, W.; Hildebrandt, P. Quantum Mechanics/Molecular

456

Mechanics Calculation of the Raman Spectra of the Phycocyanobilin Chromophore in

457

Alpha-C-Phycocyanin. Biophys. J. 2007, 93, 1885–1894.

458

459

(38) Miani, A.; Raugei, S.; Carloni, P.; Helfand, M. S. Structure and Raman Spectrum of Clavulanic Acid in Aqueous Solution. J. Phys. Chem. B. 2007, 111, 2621–2630.

460

(39) Tanzi, L.; Ramondo, F.; Guidoni, L. Vibrational Spectra of Water Solutions of Azoles

461

from QM/MM Calculations: Effects of Solvation. J. Phys. Chem. A 2012, 116, 10160–

462

10171.

463

(40) Martinez, M.; Gaigeot, M. P.; Borgis, D.; Vuilleumier, R. Extracting Effective Normal

464

Modes from Equilibrium Dynamics at Finite Temperature. J. Chem. Phys. 2006, 125,

465

144106.

466

467

468

469

470

(41) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098–3100. (42) Perdew, J. P. Density-Functional Approximation for the Correlation Energy of the Inhomogeneous Electron Gas. Phys. Rev. B 1986, 33, 8822–8824. (43) Frisch, M. J. et al. Gaussian 03 ; 2004. 26

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

471

(44) Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations.

472

Potentials for K to Au Including the Outermost Core Orbitals. J. Chem. Phys. 1985,

473

82, 299–310.

474

475

476

477

(45) Hutter, J.; Alavi, A.; Deutch, T.; Bernasconi, M.; Goedecker, S.; Marx, D.; Tuckerman, M.; Parrinello, M. CPMD; 1995. (46) Troullier, M.; Martins, J. L. Efficient Pseudopotentials for Plane-Wave Calculations. Phys. Rev. B 1991, 43, 1993–2006.

478

(47) Martyna, G.; Tuckerman, M. A Reciprocal Space Based Method for Treating Long

479

Range Interactions in Ab Initio and Force-Field-Based Calculations in Clusters. J.

480

Chem. Phys. 1999,

481

482

(48) Milburn, G. H. W.; Truter, M. R. The Crystal Structures of Cis- and TransDichlorodiammineplatinum(II). J. Chem. Soc. A 1966, 1609–1616.

483

(49) Ting, V. P.; Schmidtmann, M.; Wilson, C. C.; Weller, M. T. Cisplatin: Polymorphism

484

and Structural Insights into an Important Chemotherapeutic Drug. Angew. Chem. Int.

485

Ed. 2010, 49, 9408–9411.

486

487

488

489

(50) Marques, M. P. M.; Valero, R.; Parker, S. F.; Tomkinson, J.; de Carvalho, L. A. E. B. Polymorphism in Cisplatin Anticancer Drug. J. Phys. Chem. B 2013, 117, 6421–9. (51) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1972, 13, 5188–5192.

490

(52) Laio, A.; VandeVondele, J.; Röthlisberger, U. A Hamiltonian Electrostatic Coupling

491

Scheme for Hybrid Car–Parrinello Molecular Dynamics Simulations. J. Chem. Phy.

492

2002, 116, 6941–6947.

493

(53) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L.

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

494

Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys.

495

1983, 79, 926–935.

496

497

498

499

(54) Nosé, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511–519. (55) Parrinello, M.; Rahman, A. Crystal Structure and Pair Potentials: A MolecularDynamics Study. Phys. Rev. Lett. 1980, 45, 1196–1199.

500

(56) Weiner, S. J.; Kollman, P. A.; Case, D. A.; Singh, U. C.; Ghio, C.; Alagona, G.;

501

Profeta, S. J.; Weiner, P. K. A New Force Field for Molecular Mechanical Simulation

502

of Nucliec Acids and Proteins. J. Am. Chem. Soc. 1984, 106, 765–784.

503

(57) Scheeff, E. D.; Briggs, J. M.; Howell, S. B. Molecular Modeling of the Intrastrand

504

Guanine-Guanine DNA Adducts Produced by Cisplatin and Oxaliplatin. Mol. Phar-

505

macol. 1999, 56, 633–643.

506

507

508

509

(58) Ponder, J. W.; Case, D. A. Force Fields for Protein Simulations. Adv. Protein Chem. 2003, 66, 27–85. (59) Wilson, E. B. J.; Decius, J. C. Molecular vibrations: The Theory of Infrared and Raman Vibrational Spectra; McGraw-Hill: New York, 1955.

510

(60) McQuarrie, D. A. Statistical Mechanics; Harper & Row: New York, 1976.

511

(61) Gaigeot, M. P.; Sprik, M. Ab initio Molecular Dynamics Computation of the Infrared

512

Spectrum of Aqueous Uracil. J. Phys. Chem. B 2003, 107, 10344–10358.

513

(62) Ramirez, R.; Lopez-Ciudad, T.; Kumar, P.; Marx, D. Quantum Corrections to Classical

514

Time-Correlation Functions: Hydrogen Bonding and Anharmonic Floppy Modes. J.

515

Chem. Phys. 2004, 121, 3973–3983.

516

517

(63) Shuker, R.; Gammon, R. W. Raman-Scattering Selection-Rule Breaking and the Density of States in Amorphous Materials. Phys. Rev. Lett. 1970, 25, 222–225. 28

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

518

519

520

521

The Journal of Physical Chemistry

(64) Gordon, R. G. Relations between Raman Spectroscopy and Nuclear Spin Relaxation. J. Chem. Phys. 1965, 42, 3658–3665. (65) Skinner, J. L. Semiclassical Approximation to Golden Rule Rate Constans. J. Chem. Phys. 1997, 107, 8717–8718.

522

(66) Pulay, P.; Fogarasi, G.; Pang, F.; Boggs, J. E. Systematic Ab Initio Calculation of

523

Molecular Geometries, Force Constans, and Dipole Moment Derivatives. J. Am. Chem.

524

Soc. 1979, 101, 2550–2560.

525

(67) Neugebauer, J.; Hess, B. A. Fundamental Vibrational Frequencies of Small Polyatomic

526

Molecules from Density-functional Calculations and Vibrational Perturbation Theory.

527

J. Chem. Phys. 2003, 118, 7215–7225.

528

(68) Zhang, C.; Khaliullin, R.; Bovi, D.; Guidoni, L.; Kühne, T. Vibrational Signature of

529

Water Molecules in Asymmetric Hydrogen Bonding Environments. J. Phys. Chem. Lett.

530

2013, 4, 3245–3250.

531

532

(69) Ong, S. W.; Tok, E. S.; Kang, H. C. Vibrational frequencies in Car-Parrinello molecular dynamics. Phys. Chem. Chem. Phys. 2010, 12, 14960–14966.

533

(70) Carvalho, L. A. E. B. D.; Marques, M. P. M.; Martin, C.; Parker, S. F.; Tomkinson, J.

534

Inelastic Neutron Scattering Study of PtII Complexes Displaying Anticancer Properties.

535

ChemPhysChem 2011, 12, 1334–1341.

536

537

(71) Kühne, T. D.; Khaliullin, R. Z. Electronic Signature of the Instantaneous Asymmetry in the First Coordination Shell of Liquid Water. Nat. Commun. 2013, 4, 1450.

29

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 30 of 30