Mechanism of bleomycin suicide: a Car-Parrinello molecular dynamics investigation

Using first-principles molecular dynamics simulations (Car-Parrinello method) we investigated the possible reaction pathways for decay of the active bleomycin-Fe(III)-OOH complex, so-called bleomycin suicide. The theoretical model of activated bleomycin contains the whole metal bonding domain of the bleomycin ligand. Simulations performed both in a vacuum and in water show that a facile decaying process involves a homolytic O-O bond cleavage with an almost simultaneous hydrogen atom abstraction. The formation of an intra- or intermolecular hydrogen bond appears to be crucial for the decay of the activated bleomycin. We did not observe any evidence of heterolytic cleavage of the O-O bond of the Fe(III)-OOH species.     

Car-Parrinello molecular dynamics simulations and EPR property calculations on aqueous ubisemiquinone radical anion

Car-Parrinello molecular dynamics (CPMD) simulations have been performed on ubisemiquinone radical anion in aqueous solution. The different types of hydrogen bonding formed between the semiquinone and the solvent were studied in terms of frequency and directionality, in comparison with the parent benzosemiquinone radical anion. The EPR parameters (g-tensors and hyperfine coupling constants) were obtained from quantum chemical property calculations performed on snapshots along the MD trajectory, and the contributions of different solvation effects to the EPR parameters have been evaluated. The influence of the anion’s conformational behaviour was examined, including the orientation-dependent effects of hyperconjugation on side-chain hyperfine coupling.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.     

Formation and decay of tetrazane derivatives--a Car-Parrinello molecular dynamics study

The complications during flight 510 of the Ariane Project were ascribed to problems in the upper stage engine that employs the bipropellant monomethylhydrazine (MMH) and nitrogen tetroxide (NTO). This has led to the question what conditions or reactions possibly cause an uncontrolled behaviour in the combustion process of MMH/NTO. We use first-principles molecular dynamics to investigate the reactions of the hypergolic mixture in different chemical situations. It was possible to observe the ultrafast redox reaction between the reactants on the timescale of an unconstrained simulation. We show that electrostatic attraction is crucial for the understanding of this reaction. Besides a cold reaction preceding the ignition, a reaction path leading to the highly reactive compound dimethyltetrazane could be identified.     

On-the-fly localization of electronic orbitals in Car-Parrinello molecular dynamics

The ab initio molecular-dynamics formalism of Car and Parrinello is extended to preserve the locality of the orbitals. The supplementary term in the Lagrangian does not affect the nuclear dynamics, but ensures "on the fly" localization of the electronic orbitals within a periodic supercell in the Gamma-point approximation. The relationship between the resulting equations of motion and the formation of a gauge-invariant Lagrangian combined with a gauge-fixing procedure is briefly discussed. The equations of motion can be used to generate a very stable and easy to implement numerical integration algorithm. It is demonstrated that this algorithm can be used to compute the trajectory of the maximally localized orbitals, known as Wannier orbitals, in ab initio molecular dynamics with only a modest increase in the overall computer time. In the present paper, the new method is implemented within the generalized gradient approximation to Kohn-Sham density-functional theory employing plane wave basis sets and atomic pseudopotentials. In the course of the presentation, we briefly discuss how the present approach can be combined with localized basis sets to design fast linear scaling ab initio molecular-dynamics methods.     

Solvation structure of hydroxyl radical by Car-Parrinello molecular dynamics

Car-Parrinello molecular dynamics simulations of a hydroxyl radical in liquid water have been performed. Structural and dynamical properties of the solvated structure have been studied in details. The partial atom-atom radial distribution functions for the hydrated hydroxyl do not show drastic differences with the radial distribution functions for liquid water. The OH is found to be a more active hydrogen bond donor and acceptor than the water molecule, but the accepted hydrogen bonds are much weaker than for the hydroxide OH- ion. The first solvation shell of the OH is less structured than the water's one and contains a considerable fraction of water molecules that are not hydrogen bonded to the hydroxyl. Part of them are found to come closer to the solvated radical than the hydrogen bonded molecules do. The lifetime of the hydrogen bonds accepted by the hydroxyl is found to be shorter than the hydrogen bond lifetime in water. A hydrogen transfer between a water molecule and the OH radical has been observed, though it is a much rarer event than a proton transfer between water and an OH- ion. The velocity autocorrelation power spectrum of the hydroxyl hydrogen shows the properties both of the OH radical in clusters and of the OH- ion in liquid.     

Anharmonic infrared and Raman spectra in Car-Parrinello molecular dynamics simulations

The infrared and Raman spectra of naphthalene crystal with inclusion of anharmonic effects have been calculated by adopting the generalized variational density functional perturbation theory in the framework of Car-Parrinello molecular dynamics simulations. The computational approach has been generalized for cells of arbitrary shape. The intermolecular interactions have been analyzed with and without the van der Waals corrections, showing the importance of such interactions in the naphthalene crystal to reproduce the structural, dynamical, and spectroscopic properties.     

A molecular dynamics investigation of the titration of a trivalent aqueous ion

We carried out a series of molecular dynamics simulations of the hydrolysis of a model trivalent metal ion in aqueous solution. We use a dissociative model for water and examine the spontaneous speciation of M³⁺ into M(OH) _n ^(3-n)+ (n =1,4) both in neutral solution and as a function of added protons and hydroxide ions. The species distributions in neutral solution correspond reasonably well with those expected for real trivalent metal ions at neutral pH. However, the change in the species distributions as a function of either added protons or hydroxide ions is much less than expected with very large concentrations of protons or hydroxide ions required to shift the species equilibria in either direction. The influence of added protons and hydroxide ions on the species distributions appears to be proportional to the average charge of the hydrolysis couples, being highest for the 3+/2+ couple and lowest for the 1+/0 and 0/1- couples. Proton exchange rates vary with proton/hydroxide ion concentration giving a minimum at intermediate values ([H+]≈ 0.166) with increasing rates at both lower and higher pH.     

Revisiting the structure of (LiCH3)n aggregates using Car-Parrinello molecular dynamics

The theoretical study of (LiMe)(n) aggregates using Car-Parrinello molecular dynamics was undertaken. With respect to a quantum chemical static treatment, this approach furnishes supplementary information about the structural parameters. Equilibrium structures are indeed stable to ca. 300 K, provided the methyl groups in the aggregates are considered to rotate essentially freely. The Li-C distance depends on the coordination number of Li and not so much on the degree of aggregation. Finally, above 650 K, the cubic LiCH(3) tetramer (which is energetically favored) undergoes an entropy-driven rearrangement to a planar structure.     

Distorted five-fold coordination of Cu2+(aq) from a Car-Parrinello molecular dynamics simulation

The solvation shell structure and dynamics of a single Cu2+ ion in a periodic box with 32 water molecules under ambient conditions has been investigated using Car-Parrinello molecular dynamics simulations in a time-window of 18 ps. Five-fold coordination with four equidistant equatorial water molecules at 2.00 A and one axial water molecule at 2.45 A from the Cu2+ ion is found. A "hole" without water molecules is found on the opposite side of the axial water. The ion-water bonding character for the equatorial water molecules is different from that of the axial water molecules, as shown by a localized orbital analysis of the electronic structure. Moreover, the calculated OD stretching vibrational band for the equatorial water molecules lies ca. 175 cm-1 below the axial-water band, in good agreement with experimental data. The equatorial-water band lies below, and the axial-water band above, the pure liquid D2O band, also in agreement with experimental data.     

Hydration of Y3+ ion: a Car-Parrinello molecular dynamics study

The solvation shell structure of Y3+ and the dynamics of the hydrated ion in an aqueous solution of 0.8 M YCl3 are studied in two conditions with and without an excess proton by using first principles molecular dynamics method. We find that the first solvation shell around Y3+ contains eight water molecules forming a square antiprism as expected from x-ray absorption near edge structure in both the conditions we examined. A detailed analysis relying upon localized orbitals reveals that the complexation of water molecules with yttrium cation leads to a substantial amount of charge redistribution particularly on the oxygen atoms, giving rise to the chemical shifts of approximately -20 ppm in 17O nuclear magnetic resonance relative to the computed nuclear shieldings of the bulk water.     

Car-Parrinello molecular dynamics study of the rearrangement of the valeramide radical cation

Car-Parrinello molecular dynamics (CPMD) studies of neutral (1) and ionized (1 (+.)) valeramide are performed with the aim of providing a rationalization for the unusual temperature effect on the dissociation pattern of 1(+.) observed in mass spectrometric experiments. According to CPMD simulations of neutral valeramide 1 performed at approximately 500 K, the conformation with the fully relaxed carbon backbone predominates (96 %). Conformational changes involving folding of the carbon backbone into conformers that would allow intramolecular H transfers are predicted not to take place spontaneously at this temperature because of the barrier heights associated with these transitions (3.5 and 6.9 kcal mol(-1)), which cannot be overcome by thermal motion alone. For 1(+.), CPMD simulations performed at approximately 300 K reveal a substantial stability of a conformation in which the carbon backbone is fully relaxed; no reaction is observed even after 7 ps. However, when conformers with already folded carbon-backbones are used as initial geometries in the CPMD simulations, the gamma-hydrogen migration (McLafferty rearrangement resulting in C(3)H(6)) is already completed within 2 ps. For this important process, the free activation energy associated with both a required conformational change and the subsequent H transfer equals 4.5 kcal mol(-1), while for the formally related delta-H shift (which eventually gives rise to the elimination of C(2)H(4)/C(2)H(5.)) it amounts to 7.0 kcal mol(-1). Since the barriers associated with conformational changes are energetically more demanding than those of the corresponding hydrogen transfers, 1(+.) is essentially trapped by conformational barriers and long-lived at approximately 300 K. At elevated temperatures (500 K), the preferred reaction (within 7.3 ps) in the CPMD simulation corresponds to the McLafferty rearrangement. The estimated free activation energy associated with this process amounts to 2.5 kcal mol(-1), while the free activation energy for the delta-H transfer equals 4.4 kcal mol(-1). This relatively small free activation energy for the McLafferty rearrangement might cause dissociation of a substantial fraction of 1(+.) prior to the time-delayed mass selection, which would reduce the C3/C2 ratio in the experiments conducted with metastable ions that have a lifetime in the order of some micros at a source temperature of 500 K.     

The Jahn-Teller effect of the Cr(2+) ion in aqueous solution: Ab initio QM/MM molecular dynamics simulations

The hydration structure of Cr(2+) has been studied using molecular dynamics (MD) simulations including three-body corrections and combined ab initio quantum mechanical/molecular mechanical (QM/MM) MD simulations at the Hartree-Fock level. The structural properties are determined in terms of radial distribution functions, coordination numbers, and several angle distributions. The mean residence time was evaluated for describing ligand exchange processes in the second hydration shell. The Jahn-Teller distorted octahedral [Cr(H(2)O)(6)](2+) complex was pronounced in the QM/MM MD simulation. The first-shell distances of Cr(2+) are in the range of 1.9-2.8 A, which are slightly larger than those observed in the cases of Cu(2+) and Ti(3+). No first-shell water exchange occurred during the simulation time of 35 ps. Several water-exchange processes were observed in the second hydration shell with a mean residence time of 7.3 ps.     

Hyperfine coupling tensors of the benzosemiquinone radical anion from Car-Parrinello molecular dynamics.

Based on Car-Parrinello ab initio molecular dynamics simulations of the benzosemiquinone radical anion in both aqueous solution and the gas phase, density functional calculations provide the currently most refined EPR hyperfine coupling (HFC) tensors of semiquinone nuclei and solvent protons. For snapshots taken at regular intervals from the molecular dynamics trajectories, cluster models with different criteria for inclusion of water molecules and an additional continuum solvent model are used to analyse the HFCs. These models provide a detailed picture of the effects of dynamics and of different intermolecular interactions on the spin-density distribution and HFC tensors. Comparison with static calculations allows an assessment of the importance of dynamical effects, and of error compensation in static DFT calculations. Solvent proton HFCs depend characteristically on the position relative to the semiquinone radical anion. A point-dipolar model works well for in-plane hydrogen-bonded protons but deviates from the quantum chemical values for out-of-plane hydrogen bonding.     

Car-Parrinello molecular dynamics study of anharmonic systems: a Mannich base in solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schr?dinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.     

Molecular structure and dynamics of off-center Cu(2+) ions and strongly coupled Cu(2+)-Cu(2+) pairs in BaF(2) crystals: electron paramagnetic resonance and electron spin relaxation studies

X-band and Q-band electron paramagnetic resonance (EPR) spectra of Cu(2+) in BaF(2) crystal were recorded in the temperature range of 4.2-200 K. Spin-Hamiltonian parameters of single Cu(2+) complexes and of Cu(2+)-Cu(2+) pairs were derived and discussed. A special attention was paid to the dimeric species. Their molecular ground state configuration was found as having antiferromagnetic intradimer coupling with the singlet-triplet splitting J=-35 cm(-1). The zero-field splitting being D=0.0365 cm(-1) at 4.2 K increases with temperature as an effect of thermal population of excited dimer configurations. Electron spin echo (ESE) method was used for measurements of electron spin lattice and phase relaxation. The spin-lattice relaxation data show that except for coupling to the host lattice phonons the Cu(2+) ions are involved in local mode motions with energy of 82 cm(-1). Phase relaxation (ESE dephasing) of single Cu(2+) ions is due to spin diffusion at low temperatures. This relaxation is hampered for temperatures higher than 30 K due to the triplet state population of neighboring Cu(2+)-Cu(2+) dimers, which disturb dipolar coupling between Cu(2+) ions. For higher temperatures the relaxation is dominated by Raman T(1) processes. Fourier transform ESE spectrum displays dipolar Cu-F splitting which allowed determination of the off-center shift of Cu(2+) as delta(s)=0.132 nm. The dynamical effects observed in EPR spectra and in electron spin relaxation both for single Cu(2+) ions and Cu(2+)-Cu(2+) pairs are discussed as due to jumps between six off-center positions in the crystal unit cell and jumps between various dimer configurations.