Modeling Zn2+−Cysteinate Complexes in Proteins

How simple is too simple? Computational perspective on importance of second-shell environment for metal-ion selectivity

The metal-ion selectivity in biomolecules represents one of the most important phenomena in bioinorganic chemistry. The open question to what extent is the selectivity in the complex bioinorganic structures such as metallopeptides determined by the first-shell ligands of the metal ion is answered herein using six model peptides complexed with the set of divalent metal ions (Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), Zn(2+), Cd(2+), and Hg(2+)) and their various first-shell representations. By calculating the differences among the free energies of complexation of metal ions in these peptides and their model (truncated) systems it is quantitatively shown that the definition of the first shell is paramount to this discussion and revolves around the chemical nature of the binding site. Despite the vast conceivable diversity of peptidic structures, that suggest certain fluidity of this definition, major contributing factors are identified and assessed based on their importance for capturing metal-ion selectivity. These factors include soft/hard character of ligands and various non-covalent interactions in the vicinity of the binding site. The relative importance of these factors is considered and specific suggestions for effective construction of the models are made. The relationship of first-shell models and their corresponding parent peptides is discussed thoroughly, both with respect to their chemical similarity and potential disparity introduced by generally "non-alignable" conformational flexibility of the two systems. It is concluded that, in special cases, this disparity can be negligible and that heeding the chemical factors contributing to selectivity during construction of the model can successfully result in models that retain the affinity profile for various metal ions with high fidelity.

Synthesis, spectroscopic studies, antimicrobial activities and antitumor of a new monodentate V-shaped Schiff base and its transition metal complexes

Reaction of 4-aminoacetophenone and 4-bromobenzaldehyde in ethanol resulted in the formation of the monodentate V-shaped Schiff base (E)-1-(4-((4-bromo-benzylidene)amino)phenyl)ethanone (L). Interaction of L with different di- and trivalent metal ions revealed disubstituted derivatives. The ligand and its complexes were characterized by elemental analysis, mass, IR and NMR spectrometry. Biological activities of the ligand and complexes against the Escherchia coli and Staphylococcus aureus bacterias, and the two fungus Aspergillus flavus and Candida albicans were screened. The cytotoxicity of the compounds were checked as antitumor agents on liver carcinoma cell line (HepG2). They exhibited in vitro broad range of antitumor activities towards the cell line; the [ZnL2(H2O)2](NO3)2 complex was stronger antitumor towards HepG2 cell line as well as two breast cancer cell lines (MCF7 and T47D) relative to cis-platin.

Predicting the stability constants of metal-ion complexes from first principles

The most important experimental quantity describing the thermodynamics of metal-ion binding with various (in)organic ligands, or biomolecules, is the stability constant of the complex (β). In principle, it can be calculated as the free-energy change associated with the metal-ion complexation, i.e., its uptake from the solution under standard conditions. Because this process is associated with the interactions of charged species, large values of interaction and solvation energies are in general involved. Using the standard thermodynamic cycle (in vacuo complexation and solvation/desolvation of the reference state and of the resulting complexes), one usually subtracts values of several hundreds of kilocalories per mole to obtain final results on the order of units or tens of kilocalories per mole. In this work, we use density functional theory and Møller-Plesset second-order perturbation theory calculations together with the conductor-like screening model for realistic solvation to calculate the stability constants of selected complexes--[M(NH3)4](2+), [M(NH3)4(H2O)2](2+), [M(Imi)(H2O)5](2+), [M(H2O)3(His)](+), [M(H2O)4(Cys)], [M(H2O)3(Cys)], [M(CH3COO)(H2O)3](+), [M(CH3COO)(H2O)5](+), [M(SCH2COO)2](2-)--with eight divalent metal ions (Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), Zn(2+), Cd(2+), and Hg(2+)). Using the currently available computational protocols, we show that it is possible to achieve a relative accuracy of 2-4 kcal·mol(-1) (1-3 orders of magnitude in β). However, because most of the computed values are affected by metal- and ligand-dependent systematic shifts, the accuracy of the "absolute" (uncorrected) values is generally lower. For metal-dependent systematic shifts, we propose the specific values to be used for the given metal ion and current protocol. At the same time, we argue that ligand-dependent shifts (which cannot be easily removed) do not influence the metal-ion selectivity of the particular site, and therefore it can be computed to within 2 kcal·mol(-1) average accuracy. Finally, a critical discussion is presented that aims at potential caveats that one may encounter in theoretical predictions of the stability constants and highlights the perspective that theoretical calculations may become both competitive and complementary tools to experimental measurements.

Spectroscopic studies and biological activity of some transition metal complexes of unusual Schiff base

Unusual Schiff base ligand, 4-ethanimidoyl-6-[(1E)-N-(2-hydroxy-4-methylphenyl)ethanimidoyl]benzene-1,3-diol, L, was synthesized via catalytic process involving the interaction of some metal ions with a macrocyclic Schiff base (MSB). The transition metal derivatives [ML(H(2)O)(4)](NO(3))(3), M=Cr(III) and Fe(III), [NiL(H(2)O)(4)](NO(3))(2), [ML(H(2)O)(2)](NO(3))(2), M=Zn(II) and Cd(II), [Cl(2)Pd(μ-Cl)(2)PdL], [PtL(Cl)(2)] and [PtL(Cl)(4)] were also synthesized from the corresponding metal species with L. The Schiff bases and complexes were characterized by elemental analysis, mass spectrometry, IR and (1)H NMR spectroscopy. The crystal structure of L was determined by X-ray analysis. The spectroscopic studies revealed a variety of structure arrangements for the complexes. The biological activities of L and metal complexes against the Escherchia coli as Gram-negative bacteria and Staphylococcus aureus as Gram-positive bacteria, and the two fungus Aspergillus flavus and Candida albicans were screened. The cytotoxicity of [PtL(Cl)(2)] complex, a cis-platin analogous, was checked as an antitumor agent on two breast cancer cell lines (MCF7 and T47D) and human liver carcinoma cell line (HepG2).

A short survey on protein blocks

Protein structures are classically described in terms of secondary structures. Even if the regular secondary structures have relevant physical meaning, their recognition from atomic coordinates has some important limitations such as uncertainties in the assignment of boundaries of helical and β-strand regions. Further, on an average about 50% of all residues are assigned to an irregular state, i.e., the coil. Thus different research teams have focused on abstracting conformation of protein backbone in the localized short stretches. Using different geometric measures, local stretches in protein structures are clustered in a chosen number of states. A prototype representative of the local structures in each cluster is generally defined. These libraries of local structures prototypes are named as "structural alphabets". We have developed a structural alphabet, named Protein Blocks, not only to approximate the protein structure, but also to predict them from sequence. Since its development, we and other teams have explored numerous new research fields using this structural alphabet. We review here some of the most interesting applications.

Proton transfer function of carbonic anhydrase: Insights from QM/MM simulations

Recent QM/MM analyses of proton transfer function of human carbonic anhydrase II (CAII) are briefly reviewed. The topics include a preliminary analysis of nuclear quadrupole coupling constant calculations for the zinc ion and more detailed analyses of microscopic pK(a) of the zinc-bound water and free energy profile for the proton transfer. From a methodological perspective, our results emphasize that performing sufficient sampling is essential to the calculation of all these quantities, which reflects the well solvated nature of CAII active site. From a mechanistic perspective, our analyses highlight the importance of electrostatics in shaping the energetics and kinetics of proton transfer in CAII for its function. We argue that once the pK(a) for the zinc-bound water is modulated to be in the proper range (approximately 7.0), proton transfer through a relatively well solvated cavity towards/from the protein surface (His64) does not require any major acceleration. Therefore, although structural details like the length of the water wire between the donor and acceptor groups still may make a non-negligible contribution, our computational results and the framework of analysis suggest that the significance of such "fine-tuning" is likely secondary to the modulation of pK(a) of the zinc-bound water. We encourage further experimental analysis with mutation of (charged) residues not in the immediate neighborhood of the zinc ion to quantitatively test this electrostatics based framework; in particular, Phi analysis based on these mutations may shed further light into the relative importance of the classical Grotthus mechanism and the "proton hole" pathway that we have proposed recently for CAII.

A comparative study of semiempirical, ab initio, and DFT methods in evaluating metal-ligand bond strength, proton affinity, and interactions between first and second shell ligands in Zn-biomimetic complexes

Although theoretical methods are now available which give very accurate results, often comparable to the experimental ones, modeling chemical or biological interesting systems often requires less demanding and less accurate theoretical methods, mainly due to computer limitations. Therefore, it is crucial to know the precision of such less reliable methods for relevant models and data. This has been done in this work for small zinc-active site models including O- (H(2)O and OH(-)) and N-donor (NH(3) and imidazole) ligands. Calculations using a number of quantum mechanical methods were carried out to determine their precision for geometries, coordination number relative stability, metal-ligand bond strengths, proton affinities, and interaction energies between first and second shell ligands. We have found that obtaining chemical accuracy can be as straightforward as HF geometry optimization with a double-zeta plus polarization basis followed by a B3LYP energy calculation with a triple-zeta quality basis set including diffuse and polarization functions. The use of levels as low as PM3 geometry optimization followed by a B3LYP single-point energy calculation with a double-zeta quality basis including polarization functions already yields useful trends in bond length, proton affinities or bond dissociation energies, provided that appropriate caution is taken with the optimized structures. The reliability of these levels of calculation has been successfully demonstrated for real biomimetic cases.

Metal active site elasticity linked to activation of homocysteine in methionine synthases

Enzymes possessing catalytic zinc centers perform a variety of fundamental processes in nature, including methyl transfer to thiols. Cobalamin-independent (MetE) and cobalamin-dependent (MetH) methionine synthases are two such enzyme families. Although they perform the same net reaction, transfer of a methyl group from methyltetrahydrofolate to homocysteine (Hcy) to form methionine, they display markedly different catalytic strategies, modular organization, and active site zinc centers. Here we report crystal structures of zinc-replete MetE and MetH, both in the presence and absence of Hcy. Structural investigation of the catalytic zinc sites of these two methyltransferases reveals an unexpected inversion of zinc geometry upon binding of Hcy and displacement of an endogenous ligand in both enzymes. In both cases a significant movement of the zinc relative to the protein scaffold accompanies inversion. These structures provide new information on the activation of thiols by zinc-containing enzymes and have led us to propose a paradigm for the mechanism of action of the catalytic zinc sites in these and related methyltransferases. Specifically, zinc is mobile in the active sites of MetE and MetH, and its dynamic nature helps facilitate the active site conformational changes necessary for thiol activation and methyl transfer.

Catalytic mechanism of guanine deaminase: an ONIOM and molecular dynamics study

The catalytic mechanism of Bacillus subtilis guanine deaminase (bGD), a Zn metalloenzyme, has been investigated by a combination of quantum mechanical calculations using the multilayered ONIOM method and molecular dynamics simulations. In contrast to a previously proposed catalytic mechanism, which requires the bound guanine to assume a rare tautomeric state, the ONIOM calculations showed that the active-site residues of the enzyme do not affect the tautomeric state of guanine, and consequently the bound guanine is a tautomer that is the most abundant in aqueous solution. Two residues, Glutamate 55 and Aspartate 114, were found to play important roles in proton shuttling in the reaction. The proposed reaction path is initiated by proton transfer from a Zn-bound water to protonate Asp114. This process may be quite complex and rather dynamic in nature, as revealed by the molecular dynamics (MD) simulations, whereby another water may bridge the Zn-bound water and Asp114, which then is eliminated by positioning of guanine in the active site. The binding of guanine stabilizes protonated Asp114 by hydrogen bond formation. Asp114 can then transfer its proton to the N3 of the bound guanine, facilitating the nucleophilic attack on C2 of the guanine by the Zn-bound hydroxide to form a tetrahedral intermediate. This occurs with a rather low barrier. Glu55 then transfers a proton from the Zn-hydroxide to the amino group of the reaction intermediate and, at this point, the C2-N2 bond has lengthened by 0.2 A compared to guanine, making C2-N2 bond cleavage more facile. The C2-N2 bond breaks forming ammonia, with an energy barrier of approximately 8.8 kcal/mol. Ammonia leaves the active site, and xanthine is freed by the cleavage of the Zn-O2 bond, with a barrier approximately 8.4 kcal/mol. Along this reaction path, the highest barrier comes from C2-N2 bond cleavage, while the barrier from the cleavage of the Zn-O2 bond is slightly smaller. The Zn-O2 bond can be broken without the assistance of water during the release of xanthine.

A Combined ONIOM Quantum Chemical−Molecular Dynamics Study of Zinc−Uracil Bond Breaking in Yeast Cytosine Deaminase

A QM/MM method that combines ONIOM quantum chemistry and molecular dynamics is developed and applied to a step in the deamination of cytosine to uracil in yeast cytosine deaminase (yCD). A two-layer ONIOM calculation is used for the reaction complex, with an inner part treated at a high level for the chemical reaction (bond breaking) and a middle part treated at a lower level for relevant protein residues that are frozen in the quantum optimization. An outer layer (protein and solvent) is treated using MD. Configurations for the entire system are generated using MD and optimized with ONIOM. The method permits the use of high-level quantum calculations along with sufficient configurational sampling to approximate the potential of mean force for certain bond-breaking reactions. A previously proposed reaction mechanism for deamination (Sklenak, S.; Yao, L. S.; Cukier, R. I.; Yan, H. G. J. Am. Chem. Soc. 2004, 126, 14879) requires breaking the bond between a catalytic zinc and the O4 of uracil in order to permit product release. Using an ONIOM approach, direct bond cleavage was found to be energetically unfavorable. In the work presented here, the combined ONIOM MD method is used to show that the barrier for bond cleavage is small, approximately 3 kcal/mol, and, consequently, should not be the rate-limiting step in the reaction.

An ab initio investigation of zinc chloro complexes

A series of geometry, frequency, and energy calculations of chloroaquazinc(II) complexes were carried out at up to the MP2/6-31+G* level. A thorough examination of all species up to and including hexacoordinate species, and with up to six chlorides, was carried out. The structures of the complexes are compared with experimental data where available. The solution chemistry of zinc(II) in the presence of chloride is discussed, and Raman spectra of zinc perchlorate with increasing amount of chloride are presented.

Cys(x)His(y)-Zn2+ interactions: possibilities and limitations of a simple pairwise force field

In zinc proteins, the Zn2+ cation frequently binds with a tetrahedral coordination to cysteine and histidine side chains. We examine the possibilities and limitations of a classical, pairwise force field for molecular dynamics of such systems. Hartree Fock and density functional calculations are used to obtain geometries, charge distributions, and association energies of side chain analogues bound to Zn2+. Both ionized and neutral cysteines are considered. Two parameterizations are obtained, then tested and compared through molecular dynamics simulations of two small, homologous proteins in explicit solvent: Protein Kinase C and the Cysteine Rich Domain (CRD) of Raf, which have two Cys3His-Zn2+ groups each. The lack of explicit polarizability and charge transfer in the force field leads to poor accuracy for the association energies, and to parameters--including the zinc charge, that depend on the number of bound cysteines and their protonation state. Nevertheless, the structures sampled with the best parameterization are in good overall agreement with experiment, and have zinc coordination geometries compatible with related structures in the Cambridge Structural Database and the Protein Data Bank. Non-optimized parameters lead to poorer structures. This suggests that while a simple force field is not appropriate for processes involving exchange between water and amino acids in the zinc coordination sphere (e.g. protein unfolding), it can be useful for equilibrium simulations of stable Cys3His zinc fingers.

A Molecular Dynamics Exploration of the Catalytic Mechanism of Yeast Cytosine Deaminase

Yeast cytosine deaminase (yCD), a zinc metalloenzyme of significant biomedical interest, is investigated by a series of molecular dynamics simulations in its free form and complexed with its reactant (cytosine), product (uracil), several reaction intermediates, and an intermediate analogue. Quantum chemical calculations, used to construct a model for the catalytic Zn ion with its ligands (two cysteines, a histidine, and one water) show, by comparison with crystal structure data, that the cysteines are deprotonated and the histidine is monoprotonated. The simulations suggest that Glu64 plays a critical role in the catalysis by yCD. The rotation of the Glu64 side-chain carboxyl group that can be protonated or deprotonated permits it to act as a proton shuttle between the Zn-bound water and cytosine and subsequent reaction intermediates. Free energy methods are used to obtain the barriers for these rotations, and they are sufficiently small to permit rotation on a nanosecond time scale. In the course of the reaction, cytosine reorients to a geometry to favor nucleophilic attack by a Zn-bound hydroxide. A stable position for a reaction product, ammonia, was located in the active site, and the free energy of exchange with a water molecule was evaluated. The simulations also reveal small motions of the C-terminus and the loop that contains Phe114 that may be important for reactant binding and product release.

Factors Governing the Protonation State of Zn-Bound Histidine in Proteins: A DFT/CDM Study

We have performed systematic theoretical studies to elucidate the factors governing the His protonation/deprotonation state in Zn-binding sites, especially those containing the ubiquitous Zn-His-Asp/Glu triad. Specifically, we have addressed the following three questions: (1) How does the transfer of the Zn-bound His imidazole proton to the second-shell Asp/Glu carboxylate oxygen depend on the composition of the other first-shell ligands and the solvent accessibility of the metal-binding site? (2) Can any second-shell ligand with a proton acceptor group such as the backbone carbonyl oxygen also act as a proton acceptor? (3) What is the effect of the Asp/Glu in the Zn-His-Asp/Glu triad on the Zn-bound water protonation state? To address these questions, we used a combination of quantum mechanical and continuum dielectric methods to compute the free energies for deprotonating a Zn-bound imidazole/water in various Zn complexes. The calculations show that whether the Zn-bound His is protonated or deprotonated depends on (1) the solvent accessibility of the metal-binding site, and (2) the Lewis acid ability of Zn, which is indirectly determined by both the first- and the second-shell Zn ligands. The calculations also show that the effect of the Zn-His-Asp/Glu interaction on the nucleophilicity of the Zn-bound water depends on the solvent accessibility of the catalytic Zn site. Furthermore, they show that the Asp/Glu side chain in the Zn-His-Asp/Glu triad can increase the negative charge of its partner, His, and create an anionic hole that may stabilize a cation in buried cavities, provided that the Zn complex is cationic/neutral. The findings of this work are in accord with available experimental data.

A combined experimental and theoretical study of divalent metal ion selectivity and function in proteins: application to E. coli ribonuclease H1

Structural and thermodynamic aspects of alkaline earth metal dication (Mg(2+), Ca(2+), Sr(2+), Ba(2+)) binding to E. coli ribonuclease H1 (RNase H1) have been investigated using both experimental and theoretical methods. The various metal-binding modes of the enzyme were explored using classical molecular dynamics simulations, and relative binding free energies were subsequently evaluated by free energy simulations. The trends in the free energies of model systems based on the simulation structures were subsequently verified using a combination of density functional theory and continuum dielectric methods. The calculations provide a physical basis for the experimental results and suggest plausible role(s) for the metal cation and the catalytically important acidic residues in protein function. Magnesium ion indirectly activates water attack of the phosphorus atom by freeing one of the active site carboxylate residues, D70, to act as a general base through its four first-shell water molecules, which prevent D70 from binding directly to Mg(2+). Calcium ion, on the other hand, inhibits enzyme activity by preventing D70 from acting as a general base through bidentate interactions with both carboxylate oxygen atoms of D70. These additional interactions to D70, in addition to the D10 and E48 monodentate interactions found for Mg(2+), enable Ca(2+) to bind tighter than the other divalent ions. However, a bare Mg(2+) ion with two or less water molecules in the first shell could bind directly to the three active-site carboxylates, in particular D70, thus inhibiting enzymatic activity. The present analyses and results could be generalized to other members of the RNase H family that possess the same structural fold and show similar metal-binding site and Mg(2+)-dependent activity.

Modeling zinc in biomolecules with the self consistent charge-density functional tight binding (SCC-DFTB) method: applications to structural and energetic analysis

Parameters for the zinc ion have been developed in the self-consistent charge density functional tight-binding (SCC-DFTB) framework. The approach was tested against B3LYP calculations for a range of systems, including small molecules that contain the typical coordination environment of zinc in biological systems (cysteine, histidine, glutamic/aspartic acids, and water) and active site models for a number of enzymes such as alcohol dehydrogenase, carbonic anhydrase, and aminopeptidase. The SCC-DFTB approach reproduces structural and energetic properties rather reliably (e.g., total and relative ligand binding energies and deprotonation energies of ligands and barriers for zinc-assisted proton transfers), as compared with B3LYP/6-311+G** or MP2/6-311+G** calculations.

First-second shell interactions in metal binding sites in proteins: a PDB survey and DFT/CDM calculations

The role of the second shell in the process of metal binding and selectivity in metalloproteins has been elucidated by combining Protein Data Bank (PDB) surveys of Mg, Mn, Ca, and Zn binding sites with density functional theory/continuum dielectric methods (DFT/CDM). Peptide backbone groups were found to be the most common second-shell ligand in Mg, Mn, Ca, and Zn binding sites, followed (in decreasing order) by Asp/Glu, Lys/Arg, Asn/Gln, and Ser/Thr side chains. Aromatic oxygen- or nitrogen-containing side chains (Tyr, His, and Trp) and sulfur-containing side chains (Cys and Met) are seldom found in the second coordination layer. The backbone and Asn/Gln side chain are ubiquitous in the metal second coordination layer as their carbonyl oxygen and amide hydrogen can act as a hydrogen-bond acceptor and donor, respectively, and can therefore partner practically every first-shell ligand. The second most common outer-shell ligand, Asp/Glu, predominantly hydrogen bonds to a metal-bound water or Zn-bound histidine and polarizes the H-O or H-N bond. In certain cases, a second-shell Asp/Glu could affect the protonation state of the metal ligand. It could also energetically stabilize a positively charged metal complex more than a neutral ligand such as the backbone and Asn/Gln side chain. As for the first shell, the second shell is predicted to contribute to the metal selectivity of the binding site by discriminating between metal cations of different ionic radii and coordination geometries. The first-shell-second-shell interaction energies decay rapidly with increasing solvent exposure of the metal binding site. They are less favorable but are of the same order of magnitude as compared to the respective metal-first-shell interaction energies. Altogether, the results indicate that the structure and properties of the second shell are dictated by those of the first layer. The outer shell is apparently designed to stabilize/protect the inner-shell and complement/enhance its properties.

Orbital symmetry as a tool for understanding the bonding in Krossing's cation

The geometric and electronic structure of Krossing's cation, Ag(eta(2)-P(4))(2)(+), which shows an unexpected planar coordination environment at the metal center and D(2)(h) symmetry both in solution and in the solid state, have been investigated using density functional theory and orbital-symmetry-based energy decomposition. The analysis reveals that the contribution from electrostatic interactions to the bond energy is greater than that of orbital interactions. Partitioning of the latter term into the irreducible representations shows that, in addition to the 5s orbital, 5p orbitals of silver act as acceptor orbitals for electron donation from sigma(P-P) orbitals (a(1)(g), b(1)(u)) and n(P) orbitals (b(3)(u)). Back-donation from the 4d(10) closed shell of Ag into sigma orbitals of the pnictogen cages (b(2)(g)) is also important. However, this contribution is shown not to determine the D(2)(h) structure, contradicting conclusions from the pioneering study of the title cation (J. Am. Chem.Soc. 2001, 123, 4603). The contributions from the irreducible representations to the stabilizing orbital interactions in the D(2)(h) structure and in its D(2)(d)-symmetric conformer are analogous, indicating that the planar coordination environment at the metal center in Ag(eta(2)-P(4))(2)(+) is induced by intermolecular rather than by intramolecular interactions. Because ethylene coordination to a metal ion is an elementary reaction step in industrial processes, the bonding in Ag(C(2)H(4))(2)(+) has been analyzed as well and compared to that in Krossing's cation. Surprisingly, similar contributions to the bond energies and an involvement of metal 4d and 5p orbitals have been found, whereas a recent atoms in molecules analysis suggested that the metal-ligand interactions in silver(I) olefin complexes fundamentally differ from those in tetrahedro P(4) complexes. The only qualitative difference between the bonding patterns in Ag(eta(2)-P(4))(2)(+) and Ag(C(2)H(4))(2)(+) is the negligible energy contribution from the b(3)(u) irreducible representation in the ethylene complex because a respective symmetry-adapted linear combination of ligand orbitals is not available.

Cys(x)His(y)-Zn2+ interactions: thiol vs. thiolate coordination

In zinc proteins, the Zn2+ cation frequently binds with a tetrahedral coordination to cysteine and histidine side chains, for example, in many DNA-binding proteins, where it plays primarily a structural role. We examine the possibility of thiolate protonation in Cys(x)His(y)-Zn2+ groups, both in proteins and in solution, through a combination of theoretical calculations and database analysis. Seventy-five percent of the thiolate-coordinated zincs in the Cambridge Structural Database are tetrahedral, while di-alkanethiol coordination always involves five or more ligands. Ab initio quantum calculations are performed on (ethanethiol/thiolate)(3)imidazole-Zn2+ complexes in vacuum, yielding geometries and gas phase basicities. Protonating one (respectively two) thiolates increases the Zn-S(thiol) distance by 0.4 A (respectively 0.3 A), providing a structural marker for protonation. The stabilities of the complexes in solution are compared by combining the gas phase basicities with continuum dielectric solvation calculations. In a continuum solvent with permittivity epsilon = 4, 20, or 80, one of three thiolates is predicted to be protonated at neutral pH. By extension, Cys4-Zn2+ groups are expected to be protonated in the same conditions. In contrast, most Cys3His and Cys4 geometries in the Protein Data Bank (PDB) appear consistent with all-thiolate Zn2+ coordination. This apparent discrepancy is resolved by two recent surveys of zinc protein structures, which suggest that these all-thiolate sites are stabilized by charged and polar groups nearby in the protein, thus overcoming their intrinsic instability. However, the experimental resolution is not sufficient in all the PDB structures to rule out a thiol/thiolate mixture, and protonated thiolates may occur in some proteins not solved at high resolution or not represented in the PDB, as suggested by recent mass spectrometry experiments; this possibility should be allowed for in X-ray structure refinement.

On the competition of the purine bases, functionalities of peptide side chains, and protecting agents for the coordination sites of dicationic cisplatin derivatives

The Pt-L bond energies of simple triammineplatinum(II) complexes, [Pt(NH(3))(3)L](2+), with oxygen-, nitrogen-, and sulfur-containing donor ligands L have been predicted and rationalized using density functional theory. The ligands L have been chosen as models for functionalities of peptide side chains, for sulfur-containing protecting agents, and for adenine and guanine sites of the DNA as the ultimate target of platinum anticancer drugs. Calculation of the Pt-L bond energy in [Pt(NH(3))(3)L](2+) reveals that the soft metal center of triammineplatinum(II) prefers N ligands over S ligands. This remarkable result has been discussed in light of several interpretations of the hard and soft acids and bases principle. The concept of orbital-symmetry-based energy decomposition has been employed for the determination of the contributions from sigma and pi orbital interactions, electrostatics, and intramolecular hydrogen bonding to the Pt-L bond energy. The calculations show that considerable differences in the bond energies of the triammineplatinum(II) complexes with N-heterocycles such as 1-methylimidazole, 9-methyladenine, and 9-methylguanine arise from electrostatics rather than from orbital interactions. Surprisingly, the net stabilization by hydrogen bonding between the (Pt)N-H group and the oxygen of 9-methylguanine is as weak as the intramolecular hydrogen bond in the aqua complex [Pt(NH(3))(3)(H(2)O)](2+), challenging the common hypothesis that DNA-active anticancer drugs require carrier ligands with N-H functionalities because of their hydrogen-bonding ability. The influence of a polarizable environment on the stability of the complexes has been investigated systematically with the dependence of the dielectric constant epsilon. With increasing epsilon, the complexes with S-containing ligands are more strongly stabilized than the complexes of the N-containing heterocycles. At epsilon = 78.4, the dielectric constant of water, 9-methylguanine remains the only purine derivative investigated which is competitive to neutral sulfur ligands. These findings are particularly important for a rationalization of the results from recent experimental studies on the competition of biological donor ligands L for coordination with the metal center of [Pt(dien)L](2+) (dien = 1,5-diamino 3-azapentane).