首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The coordination chemistry of the zinc ion in the active site of alcohol dehydrogenase has been studied by the ab initio Hartree–Fock method. Geometry optimizations were performed using analytical gradients and basis sets of double-zeta quality. Correlation effects were included at the MP 2 level. The active site was modeled by Zn(HS)2XL(H2O)0–2, where X denotes ammonia or imidazole and L denotes water, methanol, ethanol, or the corresponding aldehydes or anions. It is shown that with uncharged L-ligands the four-coordinate complexes are about 20, 17, and 40kJ/mol more stable than are the corresponding three-, five-, and six-coordinate complexes, respectively. If the L-ligand is negatively charged, only the four-coordinate complexes are stable. These results suggest that the active-site zinc ion in alcohol dehydrogenase prefers a coordination number of four during the catalytic reaction, especially when the nonprotein ligand is negatively charged. Ligand exchange at the zinc ion is likely to proceed by an associative mechanism with intermittent formation of a five-coordinate complex. The results lend no support to mechanistic proposals attributing an important catalytic role to a negatively charged five-coordinate hydroxide or alkoxide ligand. © 1994 John Wiley & Sons, Inc.  相似文献   

2.
The reactions of nickel(II), copper(II), and zinc(II) acetate salts with a potentially tetradentate biphenyl-bridged bis(pyrrole-2-yl-methyleneamine) ligand yielded three complexes with different coordination geometries. X-ray crystal structural analysis reveals that in the nickel(II) complex each nickel is five-coordinate, distorted trigonal bipyramid. In the copper(II) complex, each copper is four-coordinate, between square planar and tetrahedral. In the zinc(II) complex, each zinc is four-coordinate with a distorted tetrahedral geometry and the molar ratio of the zinc and ligand is 1 : 2.  相似文献   

3.
Neutral transition metal complexes of different alpha,omega-dimethyltripyrrins TrpyMX with M = Cu(II) and Zn(II) have been prepared with a variety of anionic halogeno and pseudohalogeno ligands X, and have been studied with respect to coordination modes and structural distortion. Only four- and five-coordinate species have been observed throughout the series. All four-coordinate species display unstrained, but distorted tetrahedral or strained and distorted square-planar coordination environments for zinc(II) and copper(II) species, respectively, thus following the expectations from simple ligand field arguments. Five-coordinate species do not form easily and were observed either in donor solvents or in the solid as 1D coordination polymers with distorted trigonal-bipyramidal coordination and different topologies.  相似文献   

4.
High-energy X-ray scattering (HEXS) has been used to understand the coordination environment of the uranyl ion in a perchlorate solution. Assuming the two coordinating oxo ligands bound to U(VI) are represented in a peak in the pair distribution function (PDF) at 1.766(1) A, integration of the peak intensity is used to quantify the charge located on the oxygens. The dioxo ligands are essentially neutral, as predicted by numerous published calculations, with a charge of -16.4(8) electrons. The peak in the PDF at 2.420(1) A is consistent with equatorial ligating waters. The intensity of this peak is inconsistent with an integral coordination number and is used to propose a solution equilibrium of five and four waters coordinating to the uranyl(VI) ion favoring the five-coordinate species. This equilibrium is then used to experimentally determine that five-coordinate uranyl is 1.19+/-0.42 kcal/mol more stable than its four-coordinate counterpart under the conditions of the experiment. Further peaks in the Fourier transform of the scattering data at 4.50, 7, and 8.7 A are attributed to uranium-solvent correlations.  相似文献   

5.
Electrospray ionization of appropriate precursors is used to deliver [Fe (III)-heme] (+) and [Fe (II)-hemeH] (+) ions as naked species in the gas phase where their ion chemistry has been examined by Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. In the naked, four-coordinate [Fe (II)-hemeH] (+) and [Fe (III)-heme] (+) ions, the intrinsic reactivity of iron(II)- and iron(III)-hemes is revealed free from any influence due to axial ligand, counterion, or solvent effects. Ligand (L) addition and ligand transfer equilibria with a series of selected neutrals are attained when [Fe (II)-hemeH] (+), corresponding to protonated Fe (II)-heme, is allowed to react in the FT-ICR cell. A Heme Cation Basicity (HCB) ladder for the various ligands toward [Fe (II)-hemeH] (+), corresponding to -Delta G degrees for the process [Fe (II)-hemeH] (+) + L --> [Fe (II)-hemeH(L)] (+) and named HCB (II), can thus be established. The so-obtained HCB (II) values are compared with the corresponding HCB (III) values for [Fe (III)-heme] (+). In spite of pronounced differences displayed by various ligands, NO shows a quite similar HCB of about 67 kJ mol (-1) at 300 K toward both ions, estimated to correspond to a binding energy of 124 kJ mol (-1). Density Functional Theory (DFT) computations confirm the experimental results, yielding very similar values of NO binding energies to [Fe (II)-hemeH] (+) and [Fe (III)-heme] (+), equal to 140 and 144 kJ mol (-1), respectively. The kinetic study of the NO association reaction supports the equilibrium HCB data and reveals that the two species share very close rate constant values both for the forward and for the reverse reaction. These gas phase results diverge markedly from the kinetics and thermodynamic behavior of NO binding to iron(II)- and iron(III)-heme proteins and model complexes in solution. The requisite of either a very labile or a vacant coordination site on iron for a facile addition of NO to occur, suggested to explain the bias for typically five-coordinate iron(II) species in solution, is fully supported by the present work.  相似文献   

6.
For complexes of transition metals (manganese, iron, cobalt, nickel) with monodentate ligands, equilibrium metal-ligand distances and ligand bond energies in the first and second coordination spheres have been calculated by the CNDO method. Some effects of ligand bond energies in different coordination spheres are analyzed. These effects significantly differ between the first and second coordination spheres. In the first sphere, the ligand bond energy is mainly determined by the nature of the central ion and the type of donor atom of the ligand, but weakly depends on the structure of the ligand. Conversely, in the second coordination sphere, the ligand bond energy weakly depends on the nature of the central ion and the type of donor atom, but considerably depends on the structure of the ligands in the first coordination sphere. In the second coordination sphere, ligand binding is determined by ligand interactions with both the central ion and the ligands of the first sphere. In the general case, when strong specific interactions between ligands are absent, the energetics of the second sphere is determined by the size of the inner-spheric ligands, which may be considered to be a specific steric effect. V. I. Vernadskii Institute of Geochemistry and Analytical Chemistry. Translated fromZhurnal Strukturnoi Khimii, Vol. 36, No. 2, pp. 370–374, March–April, 1995. Translated from L. Smolina  相似文献   

7.
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.  相似文献   

8.
侯廷军  章威  徐筱杰 《化学学报》2002,60(2):221-227
采用基于线性响应近似的自由能计算方法计算了一类hydroxamate抑制剂和MMP-2的绝对结合自由能。计算中,催化锌离子和MMP-2以及配体之间采用了非键模型。分子动力学模拟结果显示,采用非键模型时,催化Zn离子采用五配位的形式,但配位键的形式和初始结构比较有很大的差别。通过拟合,分别得到了单参数、双参数以及三参数的自由能预测模型,其中,含有常数校正项的三参数模型具有最佳的预测能力,预测自由能和实际自由能之间平均绝对误差仅为2.38kJ/mol。  相似文献   

9.
Cobalamin-dependent methionine synthase (MetH) is an important metalloenzyme responsible for the biosynthesis of methionine. It catalyzes methyl transfer from N(5)-methyl-tetrahydrofolate to homocysteine (Hcy) by using a zinc ion to activate the Hcy substrate. Density functional theory (B3LYP) calculations on the active-site model in gas phase and in a polarized continuum model were performed to study the Zn coordination changes from the substrate-unbound state to the substrate-bound state. The protein effect on the Zn(2+) coordination exchange was further investigated by ONIOM (B3LYP:AMBER)-ME and EE calculations. The Zn(2+)-coordination exchange is found to be highly unfavorable in the gas phase with a high barrier and endothermicity. In the water solution, the reaction becomes exothermic and the reaction barrier is drastically decreased to about 10.0 kcal/mol. A considerable protein effect on the coordination exchange was also found; the reaction is even more exothermic and occurs without barrier. The enzyme was suggested to constrain the zinc coordination sphere in the reactant state (Hcy-unbound state) more than that in the product state (Hcy-bound state), which promotes ligation of the Hcy substrate. Molecular dynamics simulations using molecular mechanics (MM) and PM3/MM potentials suggest a correlation between the flexibility of the Zn(2+)-binding site and regulation of the enzyme function. Directed in silico mutations of selected residues in the active site were also performed. Our studies support a dissociative mechanism starting with the Zn-O(Asn234) bond breaking followed by the Zn-S((Hcy)) bond formation; the proposed associative mechanism for the Zn(2+)-coordination exchange is not supported.  相似文献   

10.
The structure of [Cu(aq)]2+ has been investigated by using full multiple-scattering theoretical (MXAN) analysis of the copper K-edge X-ray absorption (XAS) spectrum and density functional theory (DFT) to test both ideal Td and square-planar four-coordinate, five-coordinate square-pyramidal, and six-coordinate octahedral [Cu(aq)]2+ models. The best fit was an elongated five-coordinate square pyramid with four Cu-O(eq) bonds (2 x 1.98 +/- 0.03 A and 2 x 1.95 +/- 0.03 A) and a long Cu-O(ax) bond (2.35 +/- 0.05 A). The four equatorial ligands were D2d-distorted from the mean equatorial plane by +/-(17 +/- 4) degrees, so that the overall symmetry of [Cu(H2O)5]2+ is C2v. The four-coordinate MXAN fit was nearly as good, but the water ligands (4 x 1.96 +/- 0.02 A) migrated +/-(13 +/- 4) degrees from the mean equatorial plane, making the [Cu(H2O)4]2+ model again D2d-distorted. Spectroscopically calibrated DFT calculations were carried out on the C2v elongate square-pyramidal and D2d-distorted four-coordinate MXAN copper models, providing comparative electronic structures of the experimentally observed geometries. These calculations showed 0.85e spin on Cu(II) and 0.03e electron spin on each of the four equatorial water oxygens. All covalent bonding was restricted to the equatorial plane. In the square-pyramidal model, the electrostatic Cu-O(ax) bond was worth only 96.8 kJ mol(-1), compared to 304.6 kJ mol(-1) for each Cu-O(eq) bond. Both MXAN and DFT showed the potential well of the axial bond to be broad and flat, allowing large low-energy excursions. The irregular geometry and D2d-distorted equatorial ligand set sustained by unconstrained [Cu(H2O)5]2+ warrants caution in drawing conclusions regarding structural preferences from small molecule crystal structures and raises questions about the site-structural basis of the rack-induced bonding hypothesis of blue copper proteins. Further, previously neglected protein folding thermodynamic consequences of the rack-bonding hypothesis indicate an experimental disconfirmation.  相似文献   

11.
High-level ab initio quantum chemical calculations, at the CP-dG2thaw level of theory, are reported for coordination of Na+ to a wide assortment of small organic and inorganic ligands. The ligands range in size from H to C6H6, and include 22 of the ligands for which precise relative sodium ion binding free energies have been determined by recent Fourier transform ion cyclotron resonance and guided ion beam studies. Agreement with the relative experimental values is excellent (+/-1.1 kJ mol(-1)), and agreement with the absolute scale (obtained when these relative values are pegged to the CH3NH2 "anchor" value measured in a high-pressure mass spectrometric study) is only marginally poorer, with CP-dG2thaw values exceeding the absolute experimental DeltaG(298) values by an average of 2.1 kJ mol(-1). The excellent agreement between experiment and the CP-dG2thaw technique also suggests that the additional 97 ligands surveyed here (which, in many cases, are not readily susceptible to laboratory investigation) can also be reliably fitted to the existing experimental scale. However, while CP-dG2thaw and the experimental ladder are in close accord, a small set of higher level ab initio calculations on sodium ion/ligand complexes (including several values obtained here using the W1 protocol) suggests that the CP-dG2thaw values are themselves too low by approximately 2.5 kJ mol(-1), thereby implying that the accepted laboratory values are typically 4.6 kJ mol(-1) too low. The present work also highlights the importance of Na+/ligand binding energy determinations (whether by experimental or theoretical approaches) on a case-by-case basis: trends in increasing binding energy along homologous series of compounds are not reliably predictable, nor are binding site preferences or chelating tendencies in polyfunctional compounds.  相似文献   

12.
The periplasmic nitrate reductase NAP belongs to the DMSO reductase family that regroups molybdoenzymes housing a bis-molybdopterin cofactor as the active site. Several forms of the Mo(V) state, an intermediate redox state in the catalytic cycle of the enzyme, have been evidenced by EPR spectroscopy under various conditions, but their structure and catalytic relevance are not fully understood. On the basis of structural data available from the literature, we built several models that reproduce the first coordination sphere of the molybdenum cofactor and used DFT methods to make magneto-structural correlations on EPR-detected species. "High-g" states, which are the most abundant Mo(V) species, are characterized by a low-anisotropy g tensor and a high g(min) value. We assign this signature to a six-sulfur coordination sphere in a pseudotrigonal prismatic geometry with a partial disulfide bond. The "very high-g" species is well described with a sulfido ion as the sixth ligand. The "low-g" signal can be successfully associated to a Mo(V) sulfite-oxidase-type active site with only one pterin moiety coordinated to the molybdenum ion with an oxo or sulfido axial ligand. For all these species we investigate their catalytic activity using a thermodynamic point of view on the molybdenum coordination sphere. Beyond the periplasmic nitrate reductase case, this work provides useful magneto-structural correlations to characterize EPR-detected species in mononuclear molybdoenzymes.  相似文献   

13.
This paper describes the application of PROLEADS to the flexible docking of ligands into crystallographically derived enzyme structures that are assumed to be rigid. PROLEADS uses a Tabu search methodology to perform the flexible search and an empirically derived estimate of the binding affinity to drive the docking process. The paper tests the extent to which the assumption of a rigid enzyme compromises the accuracy of the results. All-pairs docking experiments are performed for three enzymes (thrombin, thermolysin and influenza virus neuraminidase) based on six or more ligand-enzyme crystal structures for each enzyme. In 76% of the cases, PROLEADS can successfully identify the correct ligand conformation as the lowest energy configuration when the enzyme structure is derived from that ligand's crystal structure, but the methodology only docks 49% of the cases successfully when the ligand is docked against enzyme crystal structures derived from other ligands. Small movements in the enzyme structure lead to an under-prediction in the energy of the correct binding mode by up to 14 kJ/mol and in some cases this under-prediction can lead to the native mode not being recognised as the lowest energy solution. The type of movements responsible for mis-docking are: the movement of sidechains as a result of changes in C position; the movement of sidechains without changes in C position; the movement of flexible portions of main chains to facilitate the formation of hydrogen bonds; and the movement of metal atoms bound to the enzyme active site. The work illustrates that the assumption of a rigid active site can lead to errors in identification of the correct binding mode and the assessment of binding affinity, even for enzymes which show relatively small shift in atomic positions from one ligand to the next. A good docking code, such as PROLEADS, can usually dock successfully if there is induced fit in relatively rigid enzymes but there remains the need to develop improved strategies for dealing with enzyme flexibility. The work implies that treatments of enzyme flexibility which focus only on sidechain rotations will not deal with the critical shifts responsible for mis-docking of ligands in thrombin, thermolysin and neuraminidase. The paper demonstrates the utility of all pairs docking experiments as a method of assessing the effectiveness of docking methodologies in dealing with enzyme flexibility.  相似文献   

14.
An extensive set of organozinc iodides, useful for Negishi-type cross-coupling reactions, are investigated as respective cations after formal loss of iodide in the gas phase. Firstly, two new alkylzinc compounds derived from Tyrosine (Tyr) and Tryptophan (Trp) are closely examined. Secondly, the influence of specific protecting groups on the subtle balance between intra- and intermolecular coordination of zinc in these reagents is probed through trifluoroacetyl (TFA)-derivatized alkylzinc compounds. Finally, the influence of the strongly coordinating bidentate ligand N,N,N',N'-tetramethylethylenediamine (TMEDA) on the structure of alkylzinc cations is further explored in order to better understand the stability of the respective complexes towards water. A combination of electrospray (ESI)-MS/MS, accurate ion mass measurements, infrared multiple-photon dissociation (IRMPD) spectroscopy and computational modelling allowed the full characterisation of all dimethylformamide (DMF)-solvated and TMEDA-coordinated alkylzinc cations in the gas phase. The calculations indicate that the zinc cation in gas-phase alkylzinc-DMF or TMEDA-complex ions preferentially adopts a tetrahedral coordination sphere with four ligands. Additionally, conformers with only three binding partners bound to zinc but with effectively combined hydrogen-bond interactions are also found. Collision induced dissociation (CID) patterns demonstrate that the zinc-DMF interaction in tetrahedral four-coordinate mono-DMF-zinc complex ions as well as the interaction between TMEDA and zinc in the corresponding complex ions is even stronger than typical covalent bonds. In most cases, all major features of the IRMPD spectra are consistent with only a single major isomer, allowing secured identification and assignment.  相似文献   

15.
Papain-like cysteine proteases are ubiquitous proteolytic enzymes. The protonated His199/deprotonated Cys29 ion pair (cathepsin B numbering) in the active site is essential for their proper functioning. The presence of this ion pair stands in contrast to the corresponding intrinsic residue p K a values, indicating a strong influence of the enzyme environment. In the present work we show by molecular dynamics simulations on quantum mechanical/molecular mechanical (QM/MM) potentials that the ion pair is stabilized by a complex hydrogen bond network which comprises several amino acids situated in the active site of the enzyme and 2-4 water molecules. QM/MM reaction path computations for the proton transfer from His199 to the thiolate of the Cys29 moiety indicate that the ion pair is about 32-36 kJ mol (-1) more stable than the neutral form if the whole hydrogen bonding network is active. Without any hydrogen bonding network the ion pair is predicted to be significantly less stable than the neutral form. QM/MM charge deletion analysis and QM model calculations are used to quantify the stabilizing effect of the active-site residues and the L1 helix in favor of the zwitterionic form. The active-site water molecules contribute about 30 kJ mol (-1) to the overall stabilization. Disruption of the hydrogen bonding network upon substrate binding is expected to enhance the nucleophilic reactivity of the thiolate.  相似文献   

16.
A dataset of 82 protein–ligand complexes of known 3D structure and binding constant Ki was analysed to elucidate the important factors that determine the strength of protein–ligand interactions. The following parameters were investigated: the number and geometry of hydrogen bonds and ionic interactions between the protein and the ligand, the size of the lipophilic contact surface, the flexibility of the ligand, the electrostatic potential in the binding site, water molecules in the binding site, cavities along the protein–ligand interface and specific interactions between aromatic rings. Based on these parameters, a new empirical scoring function is presented that estimates the free energy of binding for a protein–ligand complex of known 3D structure. The function distinguishes between buried and solvent accessible hydrogen bonds. It tolerates deviations in the hydrogen bond geometry of up to 0.25 Å in the length and up to 30 °Cs in the hydrogen bond angle without penalizing the score. The new energy function reproduces the binding constants (ranging from 3.7 × 10-2 M to 1 × 10-14 M, corresponding to binding energies between -8 and -80 kJ/mol) of the dataset with a standard deviation of 7.3 kJ/mol corresponding to 1.3 orders of magnitude in binding affinity. The function can be evaluated very fast and is therefore also suitable for the application in a 3D database search or de novo ligand design program such as LUDI. The physical significance of the individual contributions is discussed.  相似文献   

17.
Ferredoxin:thioredoxin reductase (FTR) catalyzes the reduction of the disulfide in thioredoxin in two one-electron steps using an active site comprising a [4Fe-4S] in close proximity to a redox active disulfide. M?ssbauer spectroscopy has been used to investigate the ligation and electronic properties of the [4Fe-4S] cluster in as-prepared FTR which has the active-site disulfide intact and in the N-ethylmaleimide (NEM)-modified form which provides a stable analogue of the one-electron-reduced heterodisulfide intermediate and has one of the cysteines of the active-site disulfide alkylated with NEM. The results reveal novel site-specific cluster chemistry involving weak interaction of the active-site disulfide with a unique Fe site of the [4Fe-4S]2+ cluster in the resting enzyme and cleavage of the active-site disulfide with concomitant coordination of one of the cysteines to yield a [4Fe-4S]3+ cluster with a five-coordinate Fe site ligated by two cysteine residues in the NEM-modified enzyme. The results provide molecular-level insight into the catalytic mechanism of FTR and other Fe-S-cluster-containing disulfide reductases, and suggest a possible mechanism for the reductive cleavage of S-adenosylmethionine by the radical SAM family of Fe-S enzymes.  相似文献   

18.
Binding of dioxygen to a non-heme enzyme has been modeled using the ONIOM combined quantum mechanical/molecular mechanical (QM/MM) method. For the present system, isopenicillin N synthase (IPNS), binding of dioxygen is stabilized by 8-10 kcal/mol for a QM:MM (B3LYP:Amber) protein model compared to a quantum mechanical model of the active site only. In the protein system, the free energy change of O2 binding is close to zero. Two major factors consistently stabilize O2 binding. The first effect, evaluated at the QM level, originates from a change in coordination geometry of the iron center. The active-site model artificially favors the deoxy state (O2 not bound) because it allows too-large rearrangements of the five-coordinate iron site. This error is corrected when the protein is included. The corresponding effect on binding energies is 3-6 kcal/mol, depending on the coordination mode of O2 (side-on or end-on). The second major factor that stabilizes O2 binding is van der Waals interactions between dioxygen and the surrounding enzyme. These interactions, 3-4 kcal/mol at the MM level, are neglected in models that include only the active site. Polarization of the active site by surrounding amino acids does not have a significant effect on the binding energy in the present system.  相似文献   

19.
The method of mass spectrometry with bombardment by fast atoms has been applied to the study of a number of pentacyanide complexes of iron(II) and iron(III). We have investigated the effect of the valency state of the metal and of the nature of the sixth ligand in the coordination sphere in pentacyanides of iron on the process of expelling ligands from the coordination sphere and on oxidation-reduction reactions. It is shown that the activation of small molecules by substituting a sixth ligand in the coordination sphere of pentacyanide complexes of iron(III) should be more effective than in the case of complexes of iron(II).Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 25, No. 6, pp. 733–737, November–December, 1989.  相似文献   

20.
The solvation number and conformation of N,N-dimethylacrylamide (DMAA) in the coordination sphere of the cobalt(II) ion in solution were studied, and compared with those of N,N-dimethylpropionamide (DMPA) by means of FT-Raman and FT-IR spectroscopy. Both solvents are present as either the planar cis or nonplanar staggered conformer in equilibrium, and the former is more stable in the bulk. As these solvents solvate the metal ion through the carbonyl O atom of the acryl (DMAA) or propionyl (DMPA) group, the solvation structure around the metal ion is highly congested to reduce the solvation number and/or to lead to a conformational geometry change of solvent. It turns out that the solvation number of the cobalt(II) ion is 4 for both DMAA and DMPA at 298 K, and that DMPA changes its conformation upon solvation, whereas DMAA hardly changes. The enthalpy of conformational change DeltaH degrees for DMPA is 5 kJ mol(-1) in the bulk, and is -9 kJ mol(-1) in the coordination sphere of the cobalt(II) ion. On the other hand, the DeltaH degrees value for DMAA is 9 kJ mol(-1) in the bulk.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号