Fate of cisplatin and its main hydrolysed forms in the presence of thiolates: a comprehensive computational and experimental study.

Interaction of platinum-based drugs with proteins containing sulphur amino acids is usually argued as one of the major reasons for the observed resistance to these drugs, mainly due to the deactivation of the native compounds by very efficient thiolation processes in the organism. In this work, we have investigated the detailed thermodynamics and kinetics of reaction between cisplatin cis-[PtCl2(NH3)2] and its major hydrolysed forms (monohydroxocisplatin cis-[PtCl(OH)(NH3)2] and monoaquacisplatin cis-[PtCl(H2O)(NH3)2]+) with various thiolates (methanethiolate, cysteine and glutathione) and methionine. We have used a demanding quantum chemistry approach at the MP2 and DFT levels of theory to determine the Gibbs free energies and the barrier of reactions of the most possible reaction paths. The substitution of the four ligands of the complexes studied here (Cl-, OH-, H2O and NH3) can either proceed by direct thiolations or bidentations. Our Raman spectroscopy measurements show that only two thiolations actually occur, although four are possible in principle. The reason could lie in the bidentation reactions eventually taking place after each thiolation, which is backed up by our computational results. The observed lability scale of the ligands under thiolate exposure was found to be in the following order H2O > Cl- ≈ NH3(trans) > NH3(cis) > OH-, the difference between ammine ligands being induced by a significant trans-labilization by thiolates. Finally, the S,N bidentation is shown to be preferred with respect to the S,O one.

Methionine oxidation by hydrogen peroxide in peptides and proteins: A theoretical and Raman spectroscopy study.

The oxidation of proteins results in their deterioration via the oxidation of reactive amino acids. Oxidation of the amino acid, methionine plays an important role during biological conditions of oxidative stress, and equally a role in protein stability. In this study the oxidation of the methionine residue using the tripeptide GlyMetGly with respect to hydrogen peroxide has been studied using both Raman spectroscopy and DFT calculations. Spectral modifications following the formation of methionine sulfoxide are shown with the appearance of the SO vibration whilst there is also the modification of the CS vibrations at approximately 700 cm-1. The changes in the intensity of the CS stretching band were used to calculate the kinetic rate constant as 7.9 ±â€¯0.6 × 10-3 dm3 mol-1 s-1. The energy barrier for the reaction. is determined both experimentally and using DFT calculations. The reaction of the dairy protein beta-lactoglobulin with hydrogen peroxide is equally studied using the same technique. The solvent accessible surface area of the methionine residues within the protein were also determined and a comparison of the reaction rate constant and the energy barriers of reaction for the oxidation of the tripeptide and for the protein respectively thus, provides information about the role of the protein environment in the oxidation process.

Nucleation of mercury sulfide by dealkylation.

Metal sulfide minerals are assumed to form naturally at ambient conditions via reaction of a metallic element with (poly)sulfide ions, usually produced by microbes in oxygen-depleted environments. Recently, the formation of mercury sulfide (ß-HgS) directly from linear Hg(II)-thiolate complexes (Hg(SR)2) in natural organic matter and in cysteine solutions was demonstrated under aerated conditions. Here, a detailed description of this non-sulfidic reaction is provided by computations at a high level of molecular-orbital theory. The HgS stoichiometry is obtained through the cleavage of the S-C bond in one thiolate, transfer of the resulting alkyl group (R') to another thiolate, and subsequent elimination of a sulfur atom from the second thiolate as a thioether (RSR'). Repetition of this mechanism leads to the formation of RS-(HgS)n-R chains which may self-assemble in parallel arrays to form cinnabar (α-HgS), or more commonly, quickly condense to four-coordinate metacinnabar (ß-HgS). The mechanistic pathway is thermodynamically favorable and its predicted kinetics agrees with experiment. The results provide robust theoretical support for the abiotic natural formation of nanoparticulate HgS under oxic conditions and in the absence of a catalyst, and suggest a new route for the (bio)synthesis of HgS nanoparticles with improved technological properties.

Chemical Forms of Mercury in Human Hair Reveal Sources of Exposure.

Humans are contaminated by mercury in different forms from different sources. In practice, contamination by methylmercury from fish consumption is assessed by measuring hair mercury concentration, whereas exposure to elemental and inorganic mercury from other sources is tested by analysis of blood or urine. Here, we show that diverse sources of hair mercury at concentrations as low as 0.5 ppm can be individually identified by specific coordination to C, N, and S ligands with high energy-resolution X-ray absorption spectroscopy. Methylmercury from seafood, ethylmercury used as a bactericide, inorganic mercury from dental amalgams, and exogenously derived atmospheric mercury bind in distinctive intermolecular configurations to hair proteins, as supported by molecular modeling. A mercury spike located by X-ray nanofluorescence on one hair strand could even be dated to removal of a single dental amalgam. Chemical forms of other known or putative toxic metals in human tissues could be identified by this approach with potential broader applications to forensic, energy, and materials science.

Protein reactivity with singlet oxygen: Influence of the solvent exposure of the reactive amino acid residues.

The singlet oxygen quenching rate constants were measured for three model proteins, bovine serum albumin, ß-lactoglobulin and lysozyme. The results were analyzed by comparing them with the corresponding singlet oxygen quenching rate constants for a series of tripeptides with the basic formula GlyAAGly where the central amino acid (AA) was the oxidizable amino acid, tryptophan, tyrosine, methionine and histidine. It was found that the reaction rate constant in proteins can be satisfactorily modelled by the sum of the individual contributions of the oxidizable AA residues corrected for the solvent accessible surface area (SASA) effects. The best results were obtained when the SASA of the AA residues were determined by averaging over molecular dynamics simulated trajectories of the proteins. The limits of this geometrical correction of the AA residue reactivity are also discussed.

Structure, Bonding, and Stability of Mercury Complexes with Thiolate and Thioether Ligands from High-Resolution XANES Spectroscopy and First-Principles Calculations.

We present results obtained from high energy-resolution L3-edge XANES spectroscopy and first-principles calculations for the structure, bonding, and stability of mercury(II) complexes with thiolate and thioether ligands in crystalline compounds, aqueous solution, and macromolecular natural organic matter (NOM). Core-to-valence XANES features that vary in intensity differentiate with unprecedented sensitivity the number and identity of Hg ligands and the geometry of the ligand environment. Post-Hartree-Fock XANES calculations, coupled with natural population analysis, performed on MP2-optimized Hg[(SR)2···(RSR)n] complexes show that the shape, position, and number of electronic transitions observed at high energy-resolution are directly correlated to the Hg and S (l,m)-projected empty densities of states and occupations of the hybridized Hg 6s and 5d valence orbitals. Linear two-coordination, the most common coordination geometry in mercury chemistry, yields a sharp 2p to 6s + 5d electronic transition. This transition varies in intensity for Hg bonded to thiol groups in macromolecular NOM. The intensity variation is explained by contributions from next-nearest, low-charge, thioether-type RSR ligands at 3.0-3.3 Å from Hg. Thus, Hg in NOM has two strong bonds to thiol S and k additional weak Hg···S contacts, or 2 + k coordination. The calculated stabilization energy is -5 kcal/mol per RSR ligand. Detection of distant ligands beyond the first coordination shell requires precise measurement of, and comparison to, spectra of reference compounds as well as accurate calculation of spectra for representative molecular models. The combined experimental and theoretical approaches described here for Hg can be applied to other closed-shell atoms, such as Ag(I) and Au(I). To facilitate further calculation of XANES spectra, experimental data, a new crystallographic structure of a key mercury thioether complex, Cartesian coordinates of the computed models, and examples of input files are provided as Supporting Information .

Empirical force field for cisplatin based on quantum dynamics data: case study of new parameterization scheme for coordination compounds.

Parameterization of molecular complexes containing a metallic compound, such as cisplatin, is challenging due to the unconventional coordination nature of the bonds which involve platinum atoms. In this work, we develop a new methodology of parameterization for such compounds based on quantum dynamics (QD) calculations. We show that the coordination bonds and angles are more flexible than in normal covalent compounds. The influence of explicit solvent is also shown to be crucial to determine the flexibility of cisplatin in quantum dynamics simulations. Two empirical topologies of cisplatin were produced by fitting its atomic fluctuations against QD in vacuum and QD with explicit first solvation shell of water molecules respectively. A third topology built in a standard way from the static optimized structure was used for comparison. The later one leads to an excessively rigid molecule and exhibits much smaller fluctuations of the bonds and angles than QD reveals. It is shown that accounting for the high flexibility of cisplatin molecule is needed for adequate description of its first hydration shell. MD simulations with flexible QD-based topology also reveal a significant decrease of the barrier of passive diffusion of cisplatin accross the model lipid bilayer. These results confirm that flexibility of organometallic compounds is an important feature to be considered in classical molecular dynamics topologies. Proposed methodology based on QD simulations provides a systematic way of building such topologies.

Formation of Mercury Sulfide from Hg(II)-Thiolate Complexes in Natural Organic Matter.

Methylmercury is the environmental form of neurotoxic mercury that is biomagnified in the food chain. Methylation rates are reduced when the metal is sequestered in crystalline mercury sulfides or bound to thiol groups in macromolecular natural organic matter. Mercury sulfide minerals are known to nucleate in anoxic zones, by reaction of the thiol-bound mercury with biogenic sulfide, but not in oxic environments. We present experimental evidence that mercury sulfide forms from thiol-bound mercury alone in aqueous dark systems in contact with air. The maximum amount of nanoparticulate mercury sulfide relative to thiol-bound mercury obtained by reacting dissolved mercury and soil organic matter matches that detected in the organic horizon of a contaminated soil situated downstream from Oak Ridge, TN, in the United States. The nearly identical ratios of the two forms of mercury in field and experimental systems suggest a common reaction mechanism for nucleating the mineral. We identified a chemical reaction mechanism that is thermodynamically favorable in which thiol-bound mercury polymerizes to mercury-sulfur clusters. The clusters form by elimination of sulfur from the thiol complexes via breaking of mercury-sulfur bonds as in an alkylation reaction. Addition of sulfide is not required. This nucleation mechanism provides one explanation for how mercury may be immobilized, and eventually sequestered, in oxygenated surface environments.

A theoretical study of the unfolding pathway of reduced human serum albumin.

The unfolding of the reduced human serum albumin (HSA) was simulated by generating four molecular dynamics (MD) trajectories of 160 ns each at 350, 375, 400, and 425 K, respectively. A principal components analysis (PCA) was performed on the four trajectories. Based on this analysis, 17 representative protein conformers were identified and subsequently used to construct a sequence of partially unfolded structures. They were ordered according to their decreasing α-helix fractions. The structural evolution in this unfolding sequence was found to be continuous at global but also at local level supporting the hypothesis that the protein unfolding pathway is not significantly dependent on the simulation temperature. As a result, the α-helix fraction of the protein appears to be a good reaction coordinate for the unfolding process, as it was previously suggested by experiments. Based on this observation, two conformers in the unfolding sequence were predicted to be close to the equilibrium structure of reduced HSA thus providing a first theoretical model for this protein form.

Theoretical insights into the mechanism of redox switch in heat shock protein Hsp33.

Heat shock protein 33 (Hsp33) is activated in the presence of H2O2 by a very interesting redox switch based on a tetra-coordinated zinc-cysteine complex present in the fully reduced and inactive protein form. The oxidation of this zinc center by H2O2 induces formation of two S-S bridges and the zinc release followed by the protein unfolding. We report here a theoretical study of the step-by-step sequence of the overall process starting with the oxidation of the first cysteine residue and ending with the zinc release. Each reaction step is characterized by its Gibbs free energy barrier (∆G (‡)). It is predicted that the first reaction step consists in the oxidation of Cys263 by H2O2 which is by far the most reactive cysteine (∆G (‡) = 15.4 kcal mol(-1)). The next two reaction steps are the formation of the first S-S bridge between Cys263 and Cys266 (∆G (‡) = 13.6 kcal mol(-1)) and the oxidation of Cys231 by H2O2 (∆G (‡) = 20.4 kcal mol(-1)). It is then shown that the formation of the second S-S bridge (Cys231-Cys233) before the zinc release is most unlikely (∆G (‡) = 34.8 kcal mol(-1)). Instead, the release of zinc just after the oxidation of the third cysteine (Cys231) is shown to be thermodynamically (dissociation Gibbs free energy ∆G d = 6.0 kcal mol(-1)) and kinetically (reaction rate constant k d ≈ 10(6) s(-1)) favored. This result is in good agreement with the experimental data on the oxidation mechanism of Hsp33 zinc center available to date.

An experimental and theoretical study of the amino acid side chain Raman bands in proteins.

The Raman spectra of a series of tripeptides with the basic formula GlyAAGly where the central amino acid (AA) was tryptophan, tyrosine, phenylalanine, glycine, methionine, histidine, lysine and leucine were measured in H2O. The theoretical Raman spectra obtained using density functional theory (DFT) calculations at the B3LYP/6-311+G(2df,2pd) level of theory allows a precise attribution of the vibrational bands. The experimental results show that there is a blue shift in the frequencies of several bands of the amino acid side chains in tripeptides compared to free amino acids, especially in the case of AAs containing aromatic rings. On the other hand, a very good agreement was found between the Raman bands of AA residues in tripeptides and those measured on three model proteins: bovine serum albumin, ß-lactoglobulin and lysozyme. The present analysis contributes to an unambiguous interpretation of the protein Raman spectra that is useful in monitoring the biological reactions involving AA side chains alteration.

A principal component analysis of the dynamics of subdomains and binding sites in human serum albumin.

The conformational dynamics of human serum albumin (HSA) was investigated by principal component analysis (PCA) applied to three molecular dynamics trajectories of 200 ns each. The overlap of the essential subspaces spanned by the first 10 principal components (PC) of different trajectories was about 0.3 showing that the PCA based on a trajectory length of 200 ns is not completely convergent for this protein. The contributions of the relative motion of subdomains and of the subdomains (internal) distortion to the first 10 PCs were found to be comparable. Based on the distribution of the first 3 PC, 10 protein conformers are identified showing relative root mean square deviations (RMSD) between 2.3 and 4.6 Å. The main PCs are found to be delocalized over the whole protein structure indicating that the motions of different protein subdomains are coupled. This coupling is considered as being related to the allosteric effects observed upon ligand binding to HSA. On the other hand, the first PC of one of the three trajectories describes a conformational transition of the protein domain I that is close to that experimentally observed upon myristate binding. This is a theoretical support for the older hypothesis stating that changes of the protein onformation favorable to binding can precede the ligand complexation. A detailed all atoms PCA performed on the primary Sites 1 and 2 confirms the multiconformational character of the HSA binding sites as well as the significant coupling of their motions.

Oxidation reactivity of zinc-cysteine clusters in metallothionein.

Evaluating the reactivity of the metal-thiolate clusters in metallothionein (MT) is a key step in understanding the biological functions of this protein. The effects of the metal clustering and protein environment on the thiolate reactivity with hydrogen peroxide (H(2)O(2)) were investigated by performing quantum theory calculations with chemical accuracy at two levels of complexity. At the first level, the reactivity with H(2)O(2) of a model system ([(Zn)(3)(MeS)(9)](3-), MeS is methanethiolate) of the ß domain cluster of MT was evaluated using density functional theory (DFT) with the mPW1PW91 functional. At the second level of complexity, the protein environment was included in the reactant system and the calculations were performed with the hybrid ONIOM method combining the DFT-mPW1PW91 and the semiempirical PM6 levels of theory. In these conditions, the energy barrier for the oxidation of the most reactive terminal thiolate was 21.5 kcal mol(-1). This is 3 kcal mol(-1) higher than that calculated for the terminal thiolate in the model system [(Zn)(3)(MeS)(9)](3-) and about 7 kcal mol(-1) higher than that obtained for the free thiolate. In spite of this rise of the energy barrier induced by the protein environment, the thiolate oxidation by H(2)O(2) is confirmed as a possible way for metal release from MT. On the other hand, the results suggest that the antioxidant role of MT in the living cell cannot be as important as that of glutathione (which bears a free thiol).

About the structural role of disulfide bridges in serum albumins: evidence from protein simulated unfolding.

The role of the 17 disulfide (S-S) bridges in preserving the native conformation of human serum albumin (HSA) is investigated by performing classical molecular dynamics (MD) simulations on protein structures with intact and, respectively, reduced S-S bridges. The thermal unfolding simulations predict a clear destabilization of the protein secondary structure upon reduction of the S-S bridges as well as a significant distortion of the tertiary structure that is revealed by the changes in the protein native contacts fraction. The effect of the S-S bridges reduction on the protein compactness was tested by calculating Gibbs free energy profiles with respect to the protein gyration radius. The theoretical results obtained using the OPLS-AA and the AMBER ff03 force fields are in agreement with the available experimental data. Beyond the validation of the simulation method, the results here reported provide new insights into the mechanism of the protein reductive/oxidative unfolding/folding processes. It is predicted that in the native conformation of the protein, the thiol (-SH) groups belonging to the same reduced S-S bridge are located in potential wells that maintain them in contact. The -SH pairs can be dispatched by specific conformational transitions of the peptide chain located in the neighborhood of the cysteine residues.

Oxidation of zinc-thiolate complexes of biological interest by hydrogen peroxide: a theoretical study.

Zinc-thiolate complexes play a major structural and functional role in the living cell. Their stability is directly related to the thiolate reactivity toward reactive oxygen species naturally present in the cell. Oxidation of some zinc-thiolate complexes has a functional role, as is the case of zinc finger redox switches. Herein, we report a theoretical investigation on the oxidation of thiolate by hydrogen peroxide in zinc finger cores of CCCC, CCHC, and CCHH kinds containing either cysteine or histidine residues. In the case of the CCCC core, the calculated energy barrier for the oxidation to sulfenate of the complexed thiolate was found to be 16.0 kcal mol(-1), which is 2 kcal mol(-1) higher than that for the free thiolate. The energy barrier increases to 19.3 and 22.2 kcal mol(-1) for the monoprotonated and diprotonated CCCC cores, respectively. Substitution of cysteine by histidine also induces an increase in the magnitude of the reaction energy barrier: It becomes 20.0 and 20.9 kcal mol(-1) for the CCCH and CCHH cores, respectively. It is concluded that the energy barrier for the oxidation of zinc fingers is strictly dependent on the type of ligands coordinated to zinc and on the protonation state of the complex. These changes in the thiolate reactivity can be explained by the lowering of the nucleophilicity of complexed sulfur and by the internal reorganization of the complex (changes in the metal-ligand distances) upon oxidation. The next reaction steps subsequent to sulfenate formation are also considered. The oxidized thiolate (sulfenate) is predicted to dissociate very fast: For all complexes, the calculated dissociation energy barrier is lower than 3 kcal mol(-1). It is also shown that the dissociated sulfenic acid can interact with a free thiolate to form a sulfur-sulfur (SS) bridge in a reaction that is predicted to be quasi-diffusion limited. The interesting biological consequences of the modulation of thiolate reactivity by the chemical composition of the zinc finger cores are discussed.

Free energy calculations on disulfide bridges reduction in proteins by combining ab initio and molecular mechanics methods.

Free energy profiles were calculated for the reduction of the four disulfide bridges in lysozyme by tris(2-carboxyethyl)phosphine (TCEP). The computational method combines high-precision density functional theory (DFT) calculations performed on the core of the reactant system with classical mechanical free energy evaluations based on the sampling of the configuration space of reaction environment. The predicted reaction energy barriers are in satisfactory agreement with experimental data, proving that the present method provides a reliable description of the mechanism of reaction. The role of the protein environment in this mechanism is further emphasized by analyzing the different contributions to the free energy profiles. It is shown that the protein environment affects the reaction by three factors: polarizability, steric hindrance of the reactant site, and S-S bridge distortion due to structural constraints. The corresponding effects are quantitatively evaluated, and the results are discussed in connection with the current two-step reaction model for the reduction of S-S bridges in proteins.

Cysteine oxidation by the superoxide radical: a theoretical study.

The cysteine residue oxidation by the superoxide radical in the gas phase and in aqueous solution is studied using the integrated molecular orbital+molecular orbital (IMOMO) method combining the quadratic configuration interaction [QCISD(T)] and density functional (DFT) methods. The molecular environment effects are systematically investigated by considering two alternative directions of attack of the superoxide radical on the thiol and two different cysteine residue conformations. It is found that hydrogen bonding and the electrostatic interactions between the superoxide radical and cysteine side chain significantly affect the reaction energy barrier, as compared to that derived for the simple thiol model methanethiol. Among the two possible reaction channels, the one involving the sulfinyl radical formation is predicted to be the dominant channel in aqueous solution. In a highly hydrophobic environment the thiyl radical formation channel becomes the main cysteine oxidation channel.

Protein S-S bridge reduction: a Raman and computational study of lysozyme interaction with TCEP.

The role of protein structure in the reactivity of the four disulfide (S-S) bridges of lysozyme was studied using Raman spectroscopy and molecular modelling. The experimental kinetics of S-S bridge reduction by tris-2-carboxyethyl phosphine (TCEP) was obtained by monitoring the protein S-S Raman bands. The kinetics are heterogeneous and were fitted using two apparent reaction rate constants. Kinetic measurements performed at different pH values indicate only moderate charge effects. The two intrinsic reaction rate constants derived for the neutral TCEP species were 0.45 and 0.052 mol(-1) s(-1), respectively. The molecular dynamics simulation of the reactants encounter shows that the accessibility of the lysozyme S-S bridges by TCEP decreases in the following order: cys30-cys115 > cys6-cys127 > cys64-cys80 > cys76-cys94. This simulation also illustrates the reaction mechanism which consists of a local unfolding followed by the reduction of the exposed S-S bridge. The Gibbs free energy for local unfolding was evaluated by comparing the actual reaction rate constant with that of a model system containing a fully exposed S-S bridge (oxidized glutathione). These values corresponding to the fast- and slow-reaction rate-constants were 8.5 and 13.8 kJ mol(-1), respectively. On the other hand, Raman measurements, as well as the molecular dynamics simulations, strongly suggest that the protein global unfolding following S-S bridge cleavage has only limited effects in stabilizing the reaction products.

Reductive unfolding of serum albumins uncovered by Raman spectroscopy.

The reductive unfolding of bovine serum albumin (BSA) and human serum albumin (HSA) induced by dithiothreitol (DTT) is investigated using Raman spectroscopy. The resolution of the S-S Raman band into both protein and oxidized DTT contributions provides a reliable basis for directly monitoring the S-S bridge exchange reaction. The related changes in the protein secondary structure are identified by analyzing the protein amide I Raman band. For the reduction of one S-S bridge of BSA, a mean Gibbs free energy of -7 kJ mol(-1) is derived by studying the reaction equilibrium. The corresponding value for the HSA S-S bridge reduction is -2 kJ mol(-1). The reaction kinetics observed via the S-S or amide I Raman bands are identical giving a reaction rate constant of (1.02 +/- 0.11) M(-1) s(-1) for BSA. The contribution of the conformational Gibbs free energy to the overall Gibbs free energy of reaction is further estimated by combining experimental data with ab initio calculations.

Mechanism of thiol oxidation by the superoxide radical.

In spite of the large quantity of experimental work that deals with the oxidation of thiols by superoxide, the mechanism of this reaction is still controversial. The ab initio molecular orbital calculations reported here predict that the main reaction pathway includes the formation of a three-electron-bonded adduct followed by the elimination of the hydroxide anion, giving the sulfinyl radical as the reaction product. The alternative reaction pathway consisting of hydrogen atom transfer from the thiol to the protonated superoxide radical involves a reaction energy barrier that is significantly higher. The difference between the two reaction energy barriers is clearly beyond the expected computational uncertainty. The systematic scanning of the potential energy surface reveals no other competitive reaction pathways. The present results provide a useful basis for the interpretation of the complex experimental data related to thiol oxidation by superoxide radical in a biological environment.