Aspartate 458 of human glutathione synthetase is important for cooperativity and active site structure.

Human glutathione synthetase (hGS) catalyzes the second ATP-dependent step in the biosynthesis of glutathione (GSH) and is negatively cooperative to the γ-glutamyl substrate. The hGS active site is composed of three highly conserved catalytic loops, notably the alanine rich A-loop. Experimental and computational investigations of the impact of mutation of Asp458 are reported, and thus the role of this A-loop residue on hGS structure, activity, negativity cooperativity and stability is defined. Several Asp458 hGS mutants (D458A, D458N and D458R) were constructed using site-directed mutagenesis and their activities determined (10%, 15% and 7% of wild-type hGS, respectively). The Michaelis-Menten constant (K(m)) was determined for all three substrates (glycine, GAB and ATP): glycine K(m) increased by 30-115-fold, GAB K(m) decreased by 8-17-fold, and the ATP K(m) was unchanged. All Asp458 mutants display a change in cooperativity from negative cooperativity to non-cooperative. All mutants show similar stability as compared to wild-type hGS, as determined by differential scanning calorimetry. The findings indicate that Asp458 is essential for hGS catalysis and that it impacts the allostery of hGS.

Investigation of uranyl nitrate ion pairs complexed with amide ligands using electrospray ionization ion trap mass spectrometry and density functional theory.

Ion populations formed from electrospray of uranyl nitrate solutions containing different amides vary depending on ligand nucleophilicity and steric crowding at the metal center. The most abundant species were ion pair complexes having the general formula [UO(2)(NO(3))(amide)(n=2,3)](+); however, singly charged complexes containing the amide conjugate base and reduced uranyl UO(2)(+) were also formed as were several doubly charged species. The formamide experiment produced the greatest diversity of species resulting from weaker amide binding, leading to dissociation and subsequent solvent coordination or metal reduction. Experiments using methyl formamide, dimethyl formamide, acetamide, and methyl acetamide produced ion pair and doubly charged complexes that were more abundant and less abundant complexes containing solvent or reduced uranyl. This pattern is reversed in the dimethylacetamide experiment, which displayed lower abundance doubly charged complexes, but augmented reduced uranyl complexes. DFT investigations of the tris-amide ion pair complexes showed that interligand repulsion distorts the amide ligands out of the uranyl equatorial plane and that complex stabilities do not increase with increasing amide nucleophilicity. Elimination of an amide ligand largely relieves the interligand repulsion, and the remaining amide ligands become closely aligned with the equatorial plane in the structures of the bis-amide ligands. The studies show that the phenomenological distribution of coordination complexes in a metal-ligand electrospray experiment is a function of both ligand nucleophilicity and interligand repulsion and that the latter factor begins exerting influence even in the case of relatively small ligands like the substituted methyl-formamide and methyl-acetamide ligands.

The role of the glycine triad in human glutathione synthetase.

Experimental kinetics and computational modeling of human glutathione synthetase (hGS) support the significant role of the G-loop glycine triad (G369, G370, G371) for activity of this ATP-grasp enzyme. Enzyme kinetic experiments indicate that G369V and G370V mutant hGS have little activity (<0.7 and 0.3%, respectively, versus wild-type hGS). However, G371V retains ∼13% of the activity of wild-type hGS. With respect to G-loop:A-loop interaction in hGS, mutations at Gly369 and Gly370 decrease ligand binding and prevent active site closure and protection. This research indicates that Gly369 and Gly370 have essential roles in hGS, while Gly371 has a lesser involvement. Implications for glycine-rich ensembles in other phosphate-binding enzymes are discussed.

Determination of the active site of Sphingobium chlorophenolicum 2,6-dichlorohydroquinone dioxygenase (PcpA).

2,6-Dichlorohydroquinone 1,2-dioxygenase (PcpA) from Sphingobium chlorophenolicum ATCC 39723 is a member of a class of Fe(II)-containing hydroquinone dioxygenases that is involved in the mineralization of the pollutant pentachlorophenol. This enzyme has not been extensively characterized, despite its interesting ring-cleaving activity and use of Fe(II), which are reminiscent of the well-known extradiol catechol dioxygenases. On the basis of limited sequence homology to the extradiol catechol dioxygenases, the residues ligating the Fe(II) center were originally proposed to be H159, H227, and E276 (Xu et al. in Biochemistry 38:7659-7669, 1999). However, PcpA has higher sequence homology to a newly reported, crystallographically characterized zinc metalloenzyme that has a similar predicted fold. We generated a homology model of the structure of PcpA based upon the structure of this zinc metalloenzyme. The homology model predicts that the tertiary structure of PcpA differs significantly from that of the extradiol dioxygenases, and that the residues ligating the Fe(II) are H11, H227, and E276. This structural model was tested by mutating each of H11, H159, H227, and E276 to alanine. An additional residue that is predicted to lie near the active site and is conserved among PcpA, its closest homologues, and the extradiol dioxygenases, Y266, was mutated to phenylalanine. Of these mutants, only H159A retained significant activity, thus confirming the active-site location predicted by the homology-based structural model. The model provides an important basis for understanding the origin of the unique function of PcpA.

Bonding and structure of copper nitrenes.

Copper nitrenes are of interest as intermediates in the catalytic aziridination of olefins and the amination of C-H bonds. However, despite advances in the isolation and study of late-transition-metal multiply bonded complexes, a bona fide structurally characterized example of a terminal copper nitrene has, to our knowledge, not been reported. In anticipation of such a report, terminal copper nitrenes are studied from a computational perspective. The nitrene complexes studied here are of the form (beta-diketiminate)Cu(NPh). Density functional theory (DFT), complete active space self-consistent-field (CASSCF) electronic structure techniques, and hybrid quantum mechanical/molecular mechanical (QM/MM) methods are employed to study such species. While DFT methods indicate that a triplet (S = 1) is the ground state, CASSCF calculations indicate that a singlet (S = 0) is the ground state, with only a small energy gap between the singlet and triplet. Moreover, the ground-state (open-shell) singlet copper nitrene is found to be highly multiconfigurational (i.e., biradical) and to possess a bent geometry about the nitrene nitrogen, contrasting with the linear nitrene geometry of the triplet copper nitrenes. CASSCF calculations also reveal the existence of a closed-shell singlet state with some degree of multiple bonding character for the copper-nitrene bond.

Thermodynamic and structural features of aqueous Ce(III).

With a single f-electron, Ce(III) is the simplest test case for benchmarking the thermodynamic and structural properties of hydrated Ln(III) against varying density functionals and reaction field models, in addition to determining the importance of multiconfigurational character in their wave functions. Here, the electronic structure of Ce(H2O)x(H 2O)y(3+) (x = 8, 9; y = 0, 12-14) has been examined using DFT and CASSCF calculations. The latter confirmed that the wave function of octa- and nona-aqua Ce(III) is well-described by a single configuration. Benchmarking was performed for density functionals, reaction field cavity types, and solvation reactions against the experimental free energy of hydration, DeltaG(hyd)(Ce(3+)). The UA0, UAKS, Pauling, and UFF polarized continuum model cavities displayed different performance, depending on whether one or two hydration shells were examined, and as a function of the size of the metal basis set. These results were essentially independent of the density functional employed. Using these benchmarks, the free energy for water exchange between CN = 8 and CN = 9, for which no experimental data are available, was estimated to be approximately -4 kcal/mol.

Electronic structure studies on deprotonation of dithiophosphinic acids in water clusters.

We report herein a computational study of proton transfer reactions between dithiophosphinic acids (HAs) and water clusters using B3LYP and MP2 methods. The ground-state and transition-state structures of HA-(H(2)O)(n) (n = 1, 2, 3) cluster complexes have been calculated. The influence of water molecules on energy barrier heights of proton transfer reactions has been examined in the gas phase and solution for bis[o-(trifluoromethyl)phenyl]- and bis(2,4,4-trimethylpentyl)dithiophosphinic acids (HA1 and HA2, respectively). Gas-phase calculations indicate that electron-withdrawing substituents and trifluoromethyl groups in the ortho position favor deprotonation of HA1 when three water molecules are included in the cluster. This suggests that at least three water molecules are necessary to solvate the abstracted proton in the presence of the anion. In the case of HA2, the electron-donating groups favor the reverse proton transfer reaction, namely, protonation of dithiophosphinate anion. Bulk solvent effects have been modeled for aqueous and organic media with the CPCM model. The calculated results show that polar solvents can lower the activation energy for less energetically stable transition states that have more localized charges.

A molecular modeling study on the enantioselectivity of aryl alkyl ketone reductions by a NADPH-dependent carbonyl reductase.

Automated structural analysis of Sporobolomyces salmonicolor carbonyl reductase (SSCR) indicates that the two largest potential receptor sites are in the vicinity of the nicotinamide reductant. The largest receptor site is a scalene triangle with sides of approximately 8 A by 9 A by 13 A, which is narrow in width; one corner is surrounded by hydrophilic residues that can favorably bond with the ketone oxygen. Docking aryl alkyl ketones shows a distinct preference for binding to the largest receptor site, and for conformations that place the carbonyl oxygen of the substrate in the hydrophilic corner of the largest receptor site. Favorable docking conformations for aryl alkyl ketones fall into two low-energy ensembles. These conformational ensembles are distinguished by the positions of the substituents, presenting either the Si- or Re-face of the ketone to the nicotinamide reductant. For the ketones investigated here, there is a correspondence between the major enantiomer of the alcohol obtained from the reduction of the ketone and the conformer found to have the most stable interaction energy with the receptor site in all cases. The receptor site modeling, docking simulations, molecular dynamics, and enzyme-substrate geometry optimizations lead to a model for understanding the enantioselectivity of this NADPH-dependent carbonyl reductase.

Catalytic loop motion in human glutathione synthetase: A molecular modeling approach.

Conformational changes of three flexible loops (G, A, and S) in human glutathione synthetase (hGS) arise to accommodate the substrates inside the active site. The crystal structure of hGS, a member of the ATP-grasp superfamily, has been reported only for the product-enzyme complex. To study the function of the hGS loops, molecular dynamics simulations are performed on three different conformational models: unbound enzyme, reactant-enzyme, and product-enzyme complex of hGS. The conformational changes among the three models are analyzed and the roles of the loops during the catalytic process are described. The modeled structures of hGS show that the central portions of the G- and A-loop have a double role in the reactant complex conformation: they bind the substrates and simultaneously interact with each other through an extensive network of hydrogen bonds. The present study proposes that these favorable loop-ligand and loop-loop interactions are required for opening and closing of the active site of hGS. Additionally, this research identifies important amino acid residues and explains their function within the catalytic loops of hGS.

Oxygen atom transfer energetics: assessment of the effect of method and solvent.

Several density functional methods, the semiempirical methods AM1 and PM3, Hartree-Fock, and Gaussian3 theories were applied to compute the oxygen atom transfer enthalpies for 14 X/XO couples (inorganic and organic systems, charged and neutral species, light and heavy main group element containing molecules). The calculated reaction enthalpies were compared to available experimental data. The G3 method alone was found to perform within the experimental error, while the popular B3LYP and BLYP functionals provided inadequate results. Solvent effects were estimated for 19 neutral and anionic X/XO couples by using the conductor-like polarizable continuum model and several cavity models coupled with the B3LYP/6-31++G(2d,2p) level of theory. Surprisingly, the magnitude of the aqueous solvent correction was found to vary significantly for different solute cavity models, occasionally giving larger errors than the gas-phase calculation.

The butterfly dimer [(tBu3SiO)Cr]2(mu-OSitBu3)2 and its oxidative cleavage to (tBu3SiO)2Cr(=N-N=CPh2)2 and (tBu3SiO)2Cr=N(2,6-Ph2-C6H3).

Treatment of CrCl2(THF)2 with NaOSitBu3 afforded the butterfly dimer [(tBu3SiO)Cr]2(mu-OSitBu3)2 (1(2)), whose d(CrCr) of 2.658(31) A and magnetism were indicative of strong antiferromagnetic coupling. A Boltzmann distribution of low-energy 1A1, 3B1, 5A1, 7B1, and 9A1 states obtained from calculations on [(HO)2Cr]2(muOH)2 (1'(2)) were used to provide a reasonable fit of the mu(eff) vs T data. Cleavage of 1(2) with various L (L = 4-picoline, p-tolunitrile, tBuCN, tBuNC, Ph2CO, and PMe3) generated (tBu3SiO)2CrL2 (1-L2). The dimer was oxidatively severed by Ph2CN2 to give (tBu3SiO)2Cr(N2CPh2)2 (2) and by RN3 at 23 degrees C to afford (silox)2Cr=NR (3-R) for bulky R (adamantyl (Ad), 2,6-iPr2-C6H3, 2,4,6-Me3-C6H2 = Mes, 2,6-Ph2-C6H3) and (tBu3SiO)2Cr(=NR)2 (4-R) for smaller substituents (R = 1-Naph, 2-Anth). X-ray structural studies were conducted on 1(2), square planar 1-(OCPh2)2, pseudo-Td 2 and pseudo-trigonal 3-(2,6-Ph2-C6H3), whose S = 1 ground state was discussed on the basis of calculations of (H3SiO)2Cr=NPh (3' '-Ph).

The correlation-consistent composite approach: application to the G3/99 test set.

The correlation-consistent composite approach (ccCA), an ab initio composite technique for computing atomic and molecular energies, recently has been shown to successfully reproduce experimental data for a number of systems. The ccCA is applied to the G3/99 test set, which includes 223 enthalpies of formation, 88 adiabatic ionization potentials, 58 adiabatic electron affinities, and 8 adiabatic proton affinities. Improvements on the original ccCA formalism include replacing the small basis set quadratic configuration interaction computation with a coupled cluster computation, employing a correction for scalar relativistic effects, utilizing the tight-d forms of the second-row correlation-consistent basis sets, and revisiting the basis set chosen for geometry optimization. With two types of complete basis set extrapolation of MP2 energies, ccCA results in an almost zero mean deviation for the G3/99 set (with a best value of -0.10 kcal mol(-1)), and a 0.96 kcal mol(-1) mean absolute deviation, which is equivalent to the accuracy of the G3X model chemistry. There are no optimized or empirical parameters included in the computation of ccCA energies. Except for a few systems to be discussed, ccCA performs as well as or better than Gn methods for most systems containing first-row atoms, while for systems containing second-row atoms, ccCA is an improvement over Gn model chemistries.

Stability studies of transition-metal linkage isomers using quantum mechanical methods. Groups 11 and 12 transition metals.

Several hypotheses to elucidate the linkage isomer preference of the thiocyanate (SCN(-)) ion have been offered. For complexes with small coordination numbers (i.e., 1 and 2) and groups 11 (Cu-triad) and 12 (Zn-triad) metals, different levels of theory and a variety of basis sets have been employed to study linkage isomerism. Similar results are obtained for all density functionals tested, pure and hybrid. Overall, good agreement, vis-à-vis experimentally identified linkage isomers, is achieved for ab initio techniques, whereas semiempirical quantum mechanical methods show a bias toward S-ligated isomers. Despite the seeming ease for the a priori prediction of the most stable thiocyanate isomers using acid/base principles, this research highlights the sensitivity of quantitative calculations of transition-metal linkage isomerism to the choice of basis set and electron correlation, particularly with post-Hartree-Fock treatments.

A T-shaped three-coordinate nickel(I) carbonyl complex and the geometric preferences of three-coordinate d9 complexes.

A three-coordinate diketiminate-nickel(I) complex with a carbonyl ligand has been characterized using EPR and IR spectroscopies and X-ray crystallography. The T geometry (bending from the sterically favored C(2)(v)() structure) contrasts with that of isosteric d(9) copper(II) complexes. DFT calculations on a truncated model reproduce experimental geometries, implying that the geometric differences are electronic in nature. Analysis of the charge distribution in the complexes shows that the geometry of the three-coordinate d(9) complexes is affected by differential charge donation of the ligands to the metal center.

Function of conserved residues of human glutathione synthetase: implications for the ATP-grasp enzymes.

Glutathione synthetase is an enzyme that belongs to the glutathione synthetase ATP-binding domain-like superfamily. It catalyzes the second step in the biosynthesis of glutathione from gamma-glutamylcysteine and glycine in an ATP-dependent manner. Glutathione synthetase has been purified and sequenced from a variety of biological sources; still, its exact mechanism is not fully understood. A variety of structural alignment methods were applied and four highly conserved residues of human glutathione synthetase (Glu-144, Asn-146, Lys-305, and Lys-364) were identified in the binding site. The function of these was studied by experimental and computational site-directed mutagenesis. The three-dimensional coordinates for several human glutathione synthetase mutant enzymes were obtained using molecular mechanics and molecular dynamics simulation techniques, starting from the reported crystal structure of human glutathione synthetase. Consistent with circular dichroism spectroscopy, our results showed no major changes to overall enzyme structure upon residue mutation. However, semiempirical calculations revealed that ligand binding is affected by these mutations. The key interactions between conserved residues and ligands were detected and found to be essential for enzymatic activity. Particularly, the negatively charged Glu-144 residue plays a major role in catalysis.