Influence of Outer-Shell Metal Ligands on the Structural and Electronic Properties of Horse Liver Alcohol Dehydrogenase Zinc Active Site

Formation of the iron-oxo hydroxylating species in the catalytic cycle of aromatic amino acid hydroxylases

The first part of the catalytic cycle of the pterin-dependent, dioxygen-using nonheme-iron aromatic amino acid hydroxylases, leading to the Fe(IV)=O hydroxylating intermediate, has been investigated by means of density functional theory. The starting structure in the present investigation is the water-free Fe-O(2) complex cluster model that represents the catalytically competent form of the enzymes. A model for this structure was obtained in a previous study of water-ligand dissociation from the hexacoordinate model complex of the X-ray crystal structure of the catalytic domain of phenylalanine hydroxylase in complex with the cofactor (6R)-L-erythro-5,6,7,8-tetrahydrobiopterin (BH(4)) (PAH-Fe(II)-BH(4)). The O-O bond rupture and two-electron oxidation of the cofactor are found to take place via a Fe-O-O-BH(4) bridge structure that is formed in consecutive radical reactions involving a superoxide ion, O(2)(-). The overall effective free-energy barrier to formation of the Fe(IV)=O species is calculated to be 13.9 kcal  mol(-1), less than 2 kcal  mol(-1) lower than that derived from experiment. The rate-limiting step is associated with a one-electron transfer from the cofactor to dioxygen, whereas the spin inversion needed to arrive at the quintet state in which the O-O bond cleavage is finalized, essentially proceeds without activation.

Insights into intramolecular Trp and His side-chain orientation and stereospecific π interactions surrounding metal centers: an investigation using protein metal-site mimicry in solution

Metal-binding scaffolds incorporating a Trp/His-paired epitope are instrumental in giving novel insights into the physicochemical basis of functional and mechanistic versatility conferred by the Trp-His interplay at a metal site. Herein, by coupling biometal site mimicry and (1)H and (13)C NMR spectroscopy experiments, modular constructs EDTA-(L-Trp, L-His) (EWH; EDTA=ethylenediamino tetraacetic acid) and DTPA-(L-Trp, L-His) (DWH; DTPA=diethylenetriamino pentaacetic acid) were employed to dissect the static and transient physicochemical properties of hydrophobic/hydrophilic aromatic interactive modes surrounding biometal centers. The binding feature and identities of the stoichiometric metal-bound complexes in solution were investigated by using (1)H and (13)C NMR spectroscopy, which facilitated a cross-validation of the carboxylate, amide oxygen, and tertiary amino groups as the primary ligands and indole as the secondary ligand, with the imidazole (Im) N3 nitrogen being weakly bound to metals such as Ca(2+) owing to a multivalency effect. Surrounding the metal centers, the stereospecific orientation of aromatic rings in the diastereoisomerism is interpreted with the Ca(2+)-EWH complex. With respect to perturbed Trp side-chain rotamer heterogeneity, drastically restricted Trp side-chain flexibility and thus a dynamically constrained rotamer interconversion due to π interactions is evident from the site-selective (13)C NMR spectroscopic signal broadening of the Trp indolyl C3 atom. Furthermore, effects of Trp side-chain fluctuation on indole/Im orientation were the subject of a 2D NMR spectroscopy study by using the Ca(2+)-bound state; a C-H2(indolyl)/C-H5(Im(+)) connectivity observed in the NOESY spectra captured direct evidence that the N-H1 of the Ca(2+)-Im(+) unit interacted with the pyrrole ring of the indole unit in Ca(2+)-bound EWH but not in DWH, which is assignable to a moderately static, anomalous, T-shaped, interplanar π(+)-π stacking alignment. Nevertheless, a comparative (13)C NMR spectroscopy study of the two homologous scaffolds revealed that the overall response of the indole unit arises predominantly from global attractions between the indole ring and the entire positively charged first coordination sphere. The study thus demonstrates the coordination-sphere/geometry dependence of the Trp/His side-chain interplay, and established that π interactions allow (13)C NMR spectroscopy to offer a new window for investigating Trp rotamer heterogeneity near metal-binding centers.

Mutation of outer-shell residues modulates metal ion co-ordination strength in a metalloenzyme

The metal ion co-ordination sites of many metalloproteins have been characterized by a variety of spectroscopic techniques and small-molecule model systems, revealing many important insights into the structural determinants of metal ion co-ordination. However, our understanding of this fundamentally and practically important phenomenon remains frustratingly simplistic; in many proteins it is essentially impossible to predict metal ion specificity and the effects of remote ‘outer-shell’ residues on metal ion co-ordination strength are also poorly defined. This is exemplified by our inability to explain why metalloenzymes with identical metal ion co-ordination spheres, such as the closely related orthologues of bacterial PTE (phosphotriesterase) from Agrobacterium radiobacter and Pseudomonas diminuta, display different metal ion specificity and co-ordination strength. In the present study, we present a series of PTE variants that all possess identical metal ion co-ordination spheres, yet display large differences in their metal ion co-ordination strength. Using measurement of the rates of metal ion dissociation from the active site alongside analysis of structural data obtained through X-ray crystallography, we show that ‘outer-shell’ residues provide essential support for the metal ion ligands, in effect buttressing them in their optimal orientation. Remote mutations appear to modulate metal ion interactions by increasing or decreasing the stabilizing effects of these networks. The present study therefore provides a description of how the greater protein fold can be modified to ‘tune’ the strength of metal ion co-ordination and metal ion specificity, as well as reinforcing the concept of proteins as ensembles of conformational states with unique structures and biochemical properties.

Energies, Geometries, and Charge Distributions of Zn Molecules, Clusters, and Biocenters from Coupled Cluster, Density Functional, and Neglect of Diatomic Differential Overlap Models

We present benchmark databases of Zn-ligand bond distances, bond angles, dipole moments, and bond dissociation energies for Zn-containing small molecules and Zn coordination compounds with H, CH3, C2H5, NH3, O, OH, H2O, F, Cl, S, and SCH3 ligands. The test set also includes clusters with Zn-Zn bonds. In addition, we calculated dipole moments and binding energies for Zn centers in coordination environments taken from zinc metalloenzyme X-ray structures, representing both structural and catalytic zinc centers. The benchmark values are based on relativistic-core coupled cluster calculations. These benchmark calculations are used to test the predictions of four density functionals, namely B3LYP and the more recently developed M05-2X, M06, and M06-2X levels of theory, and six semiempirical methods, including neglect of diatomic differential overlap (NDDO) calculations incorporating the new PM3 parameter set for Zn called ZnB, developed by Brothers and co-workers, and the recent PM6 parametrization of Stewart. We found that the best DFT method to reproduce dipole moments and dissociation energies of our Zn compound database is M05-2X, which is consistent with a previous study employing a much smaller and less diverse database and a much larger set of density functionals. Here we show that M05-2X geometries and single-point coupled cluster calculations with M05-2X geometries can also be used as benchmarks for larger compounds, where coupled cluster optimization is impractical, and in particular we use this strategy to extend the geometry, binding energy, and dipole moment databases to additional molecules, and we extend the tests involving crystal-site coordination compounds to two additional proteins. We find that the most predictive NDDO methods for our training set are PM3 and MNDO/d. Notably, we also find large errors in B3LYP for the coordination compounds based on experimental X-ray geometries.

Protonation state and substrate binding to B2 metallo-beta-lactamase CphA from Aeromonas hydrofila

The zinc enzymes metallo beta-lactamases counteract the beneficial action of beta-lactam antibiotics against bacterial infections, by hydrolyzing their beta-lactam rings. To understand structure/function relationships on a representative member of this class, the B2 M beta L CphA, we have investigated the H-bond pattern at the Zn enzymatic active site and substrate binding mode by molecular simulation methods. Extensive QM calculations at the DFT-BLYP level on eleven models of the protein active site, along with MD simulations of the protein in aqueous solution, allow us to propose two plausible protonation states for the unbound enzyme, which are probably in equilibrium. Docking procedures along with MD simulations and QM calculations suggest that in the complex between the enzyme and its substrate (biapenem), the latter is stable in only one of the two protonation states, in addition it exhibits two different binding modes, of which only one agrees with previous proposals. In both cases, the substrate is polarized as in aqueous solution. We conclude that addressing mechanistic issues on this class of enzymes requires a careful procedure to assign protonation states and substrate docking modes.

Computational Study of the Phosphoryl Transfer Catalyzed by a Cyclin-Dependent Kinase

A cyclin-dependent kinase, Cdk2, catalyzes the transfer of the gamma-phosphate from ATP to a threonine or serine residue of its polypeptide substrates. Here, we investigate aspects of the reaction mechanism of Cdk2 by gas-phase density functional calculations, classical molecular dynamics, and Car-Parrinello QM/MM simulations. We focus on the role of the conserved Asp127 and on the nature of the phosphoryl transfer reaction mechanism catalyzed by Cdk2. Our findings suggest that Asp127 is active in its deprotonated form by assisting the formation of the near-attack orientation of the substrate serine or threonine. Therefore, the residue does not act as a general base during the catalysis. The mechanism for the phosphoryl transfer is a single SN2-like concerted step, which shows a phosphorane-like transition state geometry. Although the resulting reaction mechanism is in agreement with a previous density functional study of the same catalytic reaction mechanism (Cavalli et al., Chem. Comm. 2003, 1308-1309), the reaction barrier is considerably lower when QM/MM calculations are performed, as in this study ( approximately 42 kcal mol(-1) QM vs. approximately 24 kcal mol(-1) QM/MM); this indicates that important roles for the catalysis are played by the protein environment and solvent waters. Because of the high amino acid sequence conservation among the whole family of cyclin-dependent kinases (CDKs), these results could be general for the CDK family.

The biochemistry of drug metabolism--an introduction: Part 2. Redox reactions and their enzymes

This review continues a general presentation of the metabolism of drugs and other xenobiotics started in a recent issue of Chemistry & Biodiversity. This Part 2 presents the numerous oxidoreductases involved, their nomenclature, relevant biochemical properties, catalytic mechanisms, and the very diverse reactions they catalyze. Many medicinally, environmentally, and toxicologically relevant examples are presented and discussed. Cytochromes P450 occupy a majority of the pages of Part 2, but a large number of relevant oxidoreductases are also considered, e.g., flavin-containing monooxygenases, amine oxidases, molybdenum hydroxylases, peroxidases, and the innumerable dehydrogenases/reductases.

Dissecting the General Physicochemical Properties of Noncovalent Interactions Involving Tyrosine Side Chain as a Second-Shell Ligand in Biomolecular Metal-Binding Site Mimetics: An Experimental Study Combining Fluorescence,13C NMR Spectroscopy and ESI Mass Spectrometry

Detailed physicochemical features inherent in the dynamic cation-pi interactions of aromatic amino acid side chains in the secondary coordination spheres around metal ions were extracted and mapped by intrinsic tyrosine fluorescence titration experiments with two homologous, artificially engineered metal-binding scaffolds which mimic metal-binding sites in metalloproteins. A newly formulated method for the treatment of fluorescence titration data allows straightforward assessment of both the magnitudes and properties of metal-chelation-assisted cation-aromatic interactions (K2) underlying a proposed two-step metallosupramolecular association process. The unprecedented linear platform-motif correlations between the two contrasting scaffolds in their changes in tyrosine fluorescence on binding of 3d metal cations help to elucidate the properties of general cation-arene recognition corresponding to the metal-responsive characteristics of the second-shell Tyr residue surrounding the metal-binding sites in the supramolecular context, and thereby define a new noncovalent design principle for metal-ion recognition in aqueous solution. As supported by NMR spectroscopic and ESI-MS analyses and molecular mechanics force field calculations, the systematic study exemplifies the concept of using steady-state tyrosine fluorescence as a powerful tool for comprehensive descriptions of cation-pi interactions in the extended environment of a metal-binding site. We established that the physicochemical properties pertaining to indirect metal-arene interactions are highly dependent on the electronic properties of the metal ions. This work suggests that second-shell cation-pi interactions may play more diverse roles, including modulation of structure, reactivity, and function of metal-binding sites, than the previously well-established direct cation-pi interactions involving hard cations (e.g., alkali metal ions). Moreover, such a study will continue to complement theoretical predications and/or the early experimental investigations in organic solvents.

Quantum chemical studies of protein structure

Quantum chemical methods now permit the prediction of many spectroscopic observables in proteins and related model systems, in addition to electrostatic properties, which are found to be in excellent accord with those determined from experiment. I discuss the developments over the past decade in these areas, including predictions of nuclear magnetic resonance chemical shifts, chemical shielding tensors, scalar couplings and hyperfine (contact) shifts, the isomer shifts and quadrupole splittings in Mössbauer spectroscopy, molecular energies and conformations, as well as a range of electrostatic properties, such as charge densities, the curvatures, Laplacians and Hessians of the charge density, electrostatic potentials, electric field gradients and electrostatic field effects. The availability of structure/spectroscopic correlations from quantum chemistry provides a basis for using numerous spectroscopic observables in determining aspects of protein structure, in determining electrostatic properties which are not readily accessible from experiment, as well as giving additional confidence in the use of these techniques to investigate questions about chemical bonding and chemical reactions.

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.