Co2+ Binding Cysteine and Selenocysteine: A DFT Study

Structure prediction of neutral physiological copper(II) compounds with l-cysteine and l-histidine

Bis(aminoacidato)copper(II) [Cu^II(aa)₂] coordination compounds are the physiological species of copper(II) amino acid compounds in blood plasma. Since there are no experimental data in the literature about the geometries that physiological Cu^II(aa)₂ could form with l-cysteine (Cys), that is, for bis(l-cysteinato)copper(II) [Cu(Cys)₂] and the ternary (l-histidinato)(l-cysteinato)copper(II) [Cu(His)(Cys)], this paper computationally examines the possible conformations that the two compounds could form with the Cys ligand having a protonated sulfur, as in the conventional zwitterion, which was determined to be prevailing in aqueous solution. These two amino acids can bind metals in a tridentate fashion and thus form many possible coordination patterns. Density functional calculations were performed for the conformational analyses in the gas phase and in implicitly modeled aqueous solution using a polarizable continuum model. Additionally, we examine which coordination mode, with thiol or thiolate group, is more stable. The Cys coordination via the amino N and carboxylato O atoms (a glycinato mode) is obtained as the most stable one in aqueous Cu(Cys)₂, and also in Cu(His)(Cys) when the His glycinato or histaminato mode combines with the intact thiol group. Whereas the conformers with N and thiol S as the copper(II) donor atoms are predicted to be the least stable, those with the Cu-N and Cu-S(thiolate) bonding (and protonated carboxylato group) are the most stable. The differences are explained by different covalent and ionic contributions of Cu-S(thiol) vs. Cu-S(thiolate). The study can contribute to the insight into formation and reactivity of the copper(II) cysteinato complexes in solution.

Theoretical Spectroscopic Studies on Chemical and Electronic Structures of Selenocysteine and Pyrrolysine

The chemical and electronic structures of the 21st and 22nd proteinogenic amino acid selenocysteine (Sec), pyrrolysine (Pyl), and their derivatives (deprotonated and protonated ions) were extensively characterized for the first time. Through the fragment based step-by-step research on their potential energy surface (PES), electronic energies of the most stable conformers of Sec, Pyl and the related ions were finally determined at the advanced CBS-QB3 and DSD-PBEP86-D3(BJ)/aug-cc-pVTZ levels, respectively, with the identification of many new low-energy conformers. The infrared spectra (IR) at 298 K of the most abundant conformers in different forms were scaled by comparison with the anharmonic frequency calculations and analyzed comparing with the experimental spectra of similar molecules. The characteristic soft X-ray spectra (including X-ray photoelectron spectra (XPS) and near-edge X-ray absorption fine-structure spectra (NEXAFS)) of the most stable conformers at 498 K were also simulated. In particular, the two possible protonated configurations of Pyl can be clearly distinguished by their different spectral features. Furthermore, a small binding energy intersection appeared around 293 eV at the C 1s edge between the canonical and protonated Pyl conformers, which is different from all the previous studies. This work thus filled the gap in our knowledge by providing detailed information on the chemical and electronic structures of Sec and Pyl and will be a useful guidance for future experimental research.

Use of cobalt(II) chelates of monothiol-containing ligands for the removal of nitric oxide

It is first reported herein that cobalt(II) complexes solution of monothiol-containing multidentate ligands are used to remove low concentration of nitric oxide (NO). These chelating ligands are water-soluble amines, alcohols or acids which containing at least one -SH group, include those of cysteine, mercaptosuccinic acid, mercaptoethanesulfonate, mercaptopropionic acid and the like. These -SH compounds when coordinated with cobalt ions, forming complexes are very effective for NO removal. The results indicate that the side group (methyl, carboxyl, carboxymethyl) on α-carbon atom of ligands contribute to the denitration of the chelate solution, whereas the substituents on sulfur atom of ligands deactivate the complexation system. In addition, we have found that several monothiol compounds with simple molecule structure and low cost exhibit good performance in denitration, and some of cobalt thiol complexes are more valuable in removing NO than ferrous thiol complexes.

Novel applications of perovskite oxide via catalytic peroxymonosulfate advanced oxidation in aqueous systems for trace L-cysteine detection

Perovskite oxides offer new opportunities in wastewater treatment via catalytic oxidation. Herein, we report a new application of perovskite oxides for biological detection via catalytic decolourisation and colorimetric determination. The presence of trace biomolecules in an aqueous system would interfere the decolourisation process of dyes, where the decolourisation rate is quantitatively correlated to the biomolecular concentration. In this work, trace L-cysteine (Cys) detection was demonstrated on the basis of a Ag-Ba_0.5Sr_0.5Co_0.75Fe_0.2O_3-δ (Ag-BSCF)/peroxymonosulfate/textile dye system. Thiol-containing cysteine can bind to Ag, Co and Fe atoms, therefore shielding the catalytic performance of the perovskite in degradation of dye solutions. Such a cost-effective biosensor system presents an excellent linear response in Cys concentration ranging from nM to μM.

Gas-phase interactions of organotin compounds with cysteine

The gas-phase interactions of cysteine with di-organotin and tri-organotin compounds have been studied by mass spectrometry experiments and quantum calculations. Positive-ion electrospray spectra show that the interaction of di- and tri-organotins with cysteine results in the formation of [(R)₂ Sn(Cys-H)]⁺ and [(R)₃ Sn(Cys)]⁺ ions, respectively. MS/MS spectra of [(R)₂ Sn(Cys-H)]⁺ complexes are characterized by numerous fragmentation processes, notably associated with elimination of NH₃ and (C,H₂ ,O₂ ). Several dissociation routes are characteristic of each given organic species. Upon collision, both the [(R)₃ Sn(Gly)]⁺ and [(R)₃ Sn(Cys)]⁺ complexes are associated with elimination of the intact amino acid, leading to the formation of [(R)₃ Sn]⁺ cation. But for the latter complex, two additional fragmentation processes are observed, associated with the elimination of NH₃ and C₃ H₄ O₂ S. Calculations indicate that the interaction between organotins and cysteine is predominantly electrostatic but also exhibits a considerable covalent character, which is slightly more pronounced in tri-organotin complexes. A preferred bidentate interaction of the type -η² -S-NH₂ , with sulfur and the amino group, is observed. As for the [(R)₃ Sn(Cys)]⁺ complexes, their stability is due to the combination of the hydrogen bond taking place between the amino group and the sulfur lone pair and the interaction between the carboxylic oxygen atom and the metal. Copyright © 2016 John Wiley & Sons, Ltd.

A RRKM study and a DFT assessment on gas-phase fragmentation of formamide-M(2+) (M = Ca, Sr)

A kinetic study of the unimolecular reactivity of formamide-M(2+) (M = Ca, Sr) systems was carried out by means of RRKM statistical theory using high-level DFT. The results predict M(2+), [M(NH2)](+) and [HCO](+) as the main products, together with an intermediate that could eventually evolve to produce [M(NH3)](2+) and CO, for high values of internal energy. In this framework, we also evaluated the influence of the external rotational energy on the reaction rate constants. In order to find a method to perform reliable electronic structure calculations for formamide-M(2+) (M = Ca, Sr) at a relatively low computational cost, an assessment of different methods was performed. In the first assessment twenty-one functionals, belonging to different DFT categories, and an MP2 wave function method using a small basis set were evaluated. CCSD(T)/cc-pWCVTZ single point calculations were used as reference. A second assessment has been performed on geometries and energies. We found BLYP/6-31G(d) and G96LYP/6-31+G(d,p) as the best performing methods, for formamide-Ca(2+) and formamide-Sr(2+), respectively. Furthermore, a detailed assessment was done on RRKM reactivity and G96LYP/6-31G(d) provided results in agreement with higher level calculations. The combination of geometrical, energetics and kinetics (RRKM) criteria to evaluate DFT functionals is rather unusual and provides an original assessment procedure. Overall, we suggest using G96LYP as the best performing functional with a small basis set for both systems.

Complex formation of Sn(II) with L-cysteine: an IR, DTA/TGA and DFT investigation

The novel complex of Sn(II) with L-cysteine (L-H2Cys) has been synthesized and characterized by elemental analysis, TGA and IR spectroscopy. Vibrational assignment and DFT/PBE0/def2-TZVP ab initio simulation give evidence of cysteine molecule being coordinated to Sn(II) as three-dentate chelating N,O,S-donor ligand. The four Perdew density functionals TPSS, PBE0, PBE, TPSSh have been tested to provide consistency of simulated and experimental IR spectra, the best result is provided by unweighted Hartree-Fock density functionals (PBE, TPSS). On the contrary, the Hartree-Fock weighted functionals (PBE0, TPPSh) provide the most accurate geometry optimization. Unharmonic frequencies are obtained via ab initio vibrational self-consistent field (PT2-VSCF) calculations at DFT/TPSS/Def2-TZVP level, the vibrational assignment of IR spectra has been carried out.

Rational design of [Co(acacen)L2]+ inhibitors of protein function

Cobalt(III) Schiff base complexes, such as [Co(acacen)L(2)](+), inhibit the function of Zn(II)-dependent proteins through dissociative exchange of the axial ligands with key histidine residues of the target protein. Consequently the efficacy of these compounds depends strongly on the lability of the axial ligands. A series of [Co(acacen)L(2)](+) complexes with various axial ligands was investigated using DFT to determine the kinetics and thermodynamics of ligand exchange and hydrolysis. Results showed excellent agreement with experimental data, indicating that axial ligand lability is determined by several factors: pK(a) of the axial ligand, the kinetic barrier to ligand dissociation, and the relative thermodynamic stability of the complexes before and after exchange. Hammett plots were constructed to determine if the kinetics and thermodynamics of exchange can be modulated by the addition of an electron-withdrawing group (EWG) to either the axial ligand itself or to the equatorial acacen ligand. Results predict that addition of an EWG to the axial ligand will shift the kinetics and thermodynamics so as to promote axial ligand exchange, while addition of an EWG to acacen will decrease axial ligand lability. These investigations will aid in the design of the next generation of [Co(acacen)L(2)](2+), allowing researchers to develop new, more effective inhibitors.

Structure of dipeptides having N-terminal selenocysteine residues: a DFT study in gas and aqueous phase

Over the last few decades, dipeptides as well as their analogues have served as important model systems for the computational studies concerning the structure of protein and energetics of protein folding. Here, we present a density functional structural study on a set of seven dipeptides having N-terminal selenocysteine residues (the component in the C-terminus is varied with seven different combinations viz. Ala, Phe, Glu, Thr, Asn, Arg and Sec) in gas and simulated aqueous phase using a polarizable continuum model (PCM). The molecular geometries of the dipeptides are fully optimized at B3LYP/6-311++G(d,p) level and subsequent frequency calculations confirm them as true minima. The effects of solvation and identity of the varying C-terminal residue on the energetics, structural features of the peptide planes, values of the ψ and ф dihedrals, geometry around the α-carbon atoms and theoretically predicted vibrational spectra of the dipeptides are investigated. Two types of intramolecular H-bonds, namely N…H-N and O…H-C, are found to play important roles in influencing the planarity of the peptide planes and geometry around the α-carbon atoms of the dipeptides. The identity of the varying C-terminal residue influences the values of ф, planarity of the peptide planes and geometry around the C₇ α-carbon atoms while the solvation effects are evident on the values of bond lengths and bond angles of the amide planes.

Collision induced dissociation of doubly-charged ions: Coulomb explosion vs. neutral loss in [Ca(urea)]2+ gas phase unimolecular reactivity via chemical dynamics simulations

In this paper we report different theoretical approaches to study the gas-phase unimolecular dissociation of the doubly-charged cation [Ca(urea)](2+), in order to rationalize recent experimental findings. Quantum mechanical plus molecular mechanical (QM/MM) direct chemical dynamics simulations were used to investigate collision induced dissociation (CID) and rotational-vibrational energy transfer for Ar + [Ca(urea)](2+) collisions. For the picosecond time-domain of the simulations, both neutral loss and Coulomb explosion reactions were found and the differences in their mechanisms elucidated. The loss of neutral urea subsequent to collision with Ar occurs via a shattering mechanism, while the formation of two singly-charged cations follows statistical (or almost statistical) dynamics. Vibrational-rotational energy transfer efficiencies obtained for trajectories that do not dissociate during the trajectory integration were used in conjunction with RRKM rate constants to approximate dissociation pathways assuming complete intramolecular vibrational energy redistribution (IVR) and statistical dynamics. This statistical limit predicts, as expected, that at long time the most stable species on the potential energy surface (PES) dominate. These results, coupled with experimental CID from which both neutral loss and Coulomb explosion products were obtained, show that the gas phase dissociation of this ion occurs by multiple mechanisms leading to different products and that reactivity on the complicated PES is dynamically complex.

On the mechanism and rate of gold incorporation into thiol-dependent flavoreductases

NADPH-dependent flavoreductases are important drug targets. During their enzymatic cycle thiolates and selenolates that have high affinity for transition metals are generated. Auranofin (AF), a gold-containing compound, is classified by the World Health Organization as an antirheumatic agent and it is indicated as the scaffold for the development of new anticancer and antiparasitic drugs. AF inhibits selenocysteine-containing flavoreductases (thioredoxin reductase and thioredoxin glutathione reductase) more effectively than non Se-containing ones (glutathione reductase); this preference has been ascribed to the high affinity of selenium for gold. We solved the 3D structure of the Se-containing Thioredoxin Glutathione Reductase from the human parasite Schistosoma mansoni complexed with Au and our results challenge this view: we believe that the relative velocity of the reaction rather than the relative affinity, depends on the presence of Sec residues, which appear to dictate AF selectivity.

Conformational preferences and pK(a) value of selenocysteine residue

The conformational preferences of the L-selenocysteine (Sec) dipeptides with selenol and selenolate groups (Ac-Sec-NHMe and Ac-Sec(-) -NHMe, respectively) and the apparent (i.e., macroscopic) pK(a) value of the Sec residue have been studied using the dispersion-corrected density functionals M06-2X and B2PLYP-D with the implicit solvation method in the gas phase and in water. In the gas phase, the backbone-to-backbone and/or side chain-to-backbone hydrogen bonds are found to contribute in stabilizing the most preferred conformations for the Sec and Sec(-) residues, as seen for the Cys and Cys(-) residues. However, the polyproline II-like conformations prevail over the conformations with the backbone-to-backbone hydrogen bonds in water because of the weakened hydrogen bonds by the favorable direct interactions between the backbone CO and HN groups and water molecules. The Sec and Sec(-) residues are found to adopt more various conformations than the Cys and Cys(-) residues in water, although the most preferred conformations of the neutral and/or anionic forms of the two residues are similar each other in the gas phase and in water. Using the statistically weighted free energies of the Sec and Sec(-) dipeptides in the gas phase and their solvation free energies, the pK(a) value of the Sec residue is estimated to be 5.47 at 25°C, which is in good agreement with the experimental value of 5.43 ± 0.02. It is found that the lower pK(a) value of the selenol side chain for the Sec residue by ∼3 units than the thiol side chain for the Cys residue is ascribed to the higher gas-phase acidity of the Sec residue.

A theoretical study of hydration effects on the prototropic tautomerism of selenouracils

The prototropic tautomerism of 2-, 4-selenouracil and 2,4-diselenouracil has been studied using density functional theory (DFT) methods, at the B3LYP/6-311 + G(3df,2p)//B3LYP/6-31G(d,p) level. The relative stability order of selenouracil tautomers does not resemble that of uracil tautomers, but it is similar to that of thiouracils, even though the energy gaps between the different tautomers of selenouracils are smaller than for thiouracils. The tautomerism activation barriers are high enough as to conclude that only the oxo-selenone or the diselenone structures should be found in the gas phase. The specific interaction with one water molecule reduces these barriers by a half, but still the oxo-selenone form is always the most stable tautomer. The addition of a second water molecule has a relatively small effect, as well as bulk effects, evaluated by means of a continuum-polarized model. For isolated 2- and 4-selenouracils, the more favorable tautomerization process corresponds to a hydrogen transfer towards the selenium atom, the activation barriers for transfer towards the oxygen atom being much higher. This situation changes when specific and bulk effects are included, and the latter process becomes the more favorable one. For 2,4-diselenouracil the more favorable tautomerization, in the gas phase, corresponds to the H shift from N1 to the Se atom at C2, while solvation effects favor the transfer from N3 to the Se atom at C4.

Structures and fragmentations of cobalt(II)-cysteine complexes in the gas phase

The electronebulization of a cobalt(II)/cysteine(Cys) mixture in water/methanol (50/50) produced mainly cobalt-cationized species. Three main groups of the Co-cationized species can be distinguished in the ESI-MS spectrum: (1) the cobalt complexes including the cysteine amino acid only (they can be singly charged, for example, [Co(Cys)n- H]+ with n = 1-3 or doubly charged such as [Co + (Cys)2]2+); (2) the cobalt complexes with methanol: [Co(CH3OH)n- H]+ with n = 1-3, [Co(CH3OH)4]2+; and (3) the complexes with the two different types of ligands: [Co(Cys)(CH3OH) - H]+. Only the singly charged complexes were observed. Collision-induced dissociation (CID) products of the [Co(Cys)2]2+, [Co(Cys)2 - H]+ and [Co(Cys) - H]+ complexes were studied as a function of the collision energy, and mechanisms for the dissociation reactions are proposed. These were supported by the results of deuterium labelling experiments and by density functional theory calculations. Since [Co(Cys) - H]+ was one of the main product ions obtained upon the CID of [Co(Cys)2]2+ and of [Co(Cys)2 - H]+ under low-energy conditions, the fragmentation pathways of [Co(Cys) - H]+ and the resulting product ion structures were studied in detail. The resulting product ion structures confirmed the high affinity of cobalt(II) for the sulfur atom of cysteine.

A Coupled Car-Parrinello Molecular Dynamics and EXAFS Data Analysis Investigation of Aqueous Co2+

We have studied the microscopic solvation structure of Co(2+) in liquid water by means of density functional theory (DFT)-based Car-Parrinello molecular dynamics (CPMD) simulations and extended X-ray absorption fine structure (EXAFS) data analysis. The effect of the number of explicit water molecules in the simulation box on the first and second hydration shell structures has been considered. Classical molecular dynamics simulations, using an effective two-body potential for Co(2+)-water interactions, were also performed to show box size effects in a larger range. We have found that the number of explicit solvent molecules has a marginal role on the first solvation shell structural parameters, whereas larger boxes may be necessary to provide a better description of the second solvation shell. Car-Parrinello simulations were determined to provide a reliable description of structural and dynamical properties of Co(2+) in liquid water. In particular, they seem to describe both the first and second hydration shells correctly. The EXAFS signal was reconstructed from Car-Parrinello simulations. Good agreement between the theoretical and experimental signals was observed, thus strengthening the microscopic picture of the Co(2+) solvation properties obtained using first-principle simulations.