Joint quantum chemical and polarizable molecular mechanics investigation of formate complexes with penta- and hexahydrated Zn2+: Comparison between energetics of model bidentate,monodentate, and through-water Zn2+ binding modes and evaluation of nonadditivity effects

In order to gain an understanding of the energetics of polycoordinated Zn²⁺ binding to the formate anion (the end side chain of the Asp and Glu residues of proteins), we compare three competing binding modes in the presence of five and six water molecules: a, bidentate binding of Zn²⁺ to both formate oxygens; b, monodentate binding of Zn²⁺ to one formate oxygen; and c, through-water binding of Zn²⁺ to formate, in which the cation remains bound to its first-hydration shell waters and interacts with both formate oxygens through three water molecules. We also investigate a complex d, which is similar to c, in which formate is protonated into formic acid and one water molecule is deprotonated. The computations are carried out using the ab initio self-consistent field/MP2 with three basis sets of increasing size density functional theory, semiempirical AM1 and PM3, and the sum of interactions between fragments ab initio computed (SIBFA) molecular mechanics procedures. The summed energies of the isolated molecules making up the complexes disfavor tautomer d compared to a–c. On the other hand, the ab initio computations give the ordering of intermolecular interaction energies as d formic acid tautomer >b monodentate >a bidentate >c through-water. Whereas the first-order energy E₁ favors both inner-shell Zn²⁺ complexes with formate over the outer-shell complex, the polarization and the charge-transfer components of the second-order energy E₂ both favor the outer-shell complex over the inner-shell one, despite the increased separation between the cation and the highly polarizable formate ion. Energy balances including continuum solvation enthalpies produce an equilibration of complexes a–d. The preference favoring the monodentate complex over the bidentate one is consistent with other ab initio results for formate binding by a fully coordinated Zn²⁺ cation and with structural results from X-ray crystallography. The SIBFA results are consistent with the ab initio results, and the computed interaction energy values match the ab initio ones to within 3%. The effects of nonadditivity are analyzed in the ab initio, SIBFA, and semiempirical computations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1379–1390, 1999     

Complexes of thiomandelate and captopril mercaptocarboxylate inhibitors to metallo-beta-lactamase by polarizable molecular mechanics. Validation on model binding sites by quantum chemistry

Using the polarizable molecular mechanics method SIBFA, we have performed a search for the most stable binding modes of D- and L-thiomandelate to a 104-residue model of the metallo-beta-lactamase from B. fragilis, an enzyme involved in the acquired resistance of bacteria to antibiotics. Energy balances taking into account solvation effects computed with a continuum reaction field procedure indicated the D-isomer to be more stably bound than the L-one, conform to the experimental result. The most stably bound complex has the S(-) ligand bridging monodentately the two Zn(II) cations and one carboxylate O(-) H-bonded to the Asn193 side chain. We have validated the SIBFA energy results by performing additional SIBFA as well as quantum chemical (QC) calculations on small (88 atoms) model complexes extracted from the 104-residue complexes, which include the residues involved in inhibitor binding. Computations were done in parallel using uncorrelated (HF) as well as correlated (DFT, LMP2, MP2) computations, and the comparisons extended to corresponding captopril complexes (Antony et al., J Comput Chem 2002, 23, 1281). The magnitudes of the SIBFA intermolecular interaction energies were found to correctly reproduce their QC counterparts and their trends for a total of twenty complexes.     

Spectroscopy and structures of copper complexes with ethylenediamine and methyl-substituted derivatives

Copper complexes of ethylenediamine (en), N-methylethylenediamine (meen), N,N-dimethylethylenediamine (dmen), N,N,N'-trimethylethylenediamine (tren), and N,N,N',N'-tetramethylethylenediamine (tmen) are synthesized in laser-vaporization supersonic molecular beams and studied by pulsed-field ionization zero electron kinetic energy (ZEKE) and photoionization efficiency spectroscopies and second-order Moller-Plesset perturbation theory. Precise ionization energies and vibrational frequencies of Cu-en, -meen, and -dmen are measured from the ZEKE spectra, and ionization thresholds of Cu-tren and -tmen are estimated from the photoionization efficiency spectra. The measured vibrational modes span a frequency range of 35-1646 cm(-1) and include metal-ligand stretch and bend, hydrogen-bond stretch, and ligand-based torsion. A number of low-energy structures with Cu binding to one or two nitrogen atoms are predicted for each complex by the ab initio calculations. The combination of the spectroscopic measurements and ab initio calculations has identified a hydrogen-bond-stabilized monodentate structure for the Cu-en complex and bidentate cyclic structures for the methyl-substituted derivatives. The change of the Cu binding from the monodentate to the bidentate mode arises from the competition between copper coordination and hydrogen bonding.     

Modeling the IR Spectra of Aqueous Metal Carboxylate Complexes: Correlation between Bonding Geometry and Stretching Mode Wavenumber Shifts

A widely used principle is that shifts in the wavenumber of carboxylate stretching modes upon bonding with a metal center can be used to infer if the geometry of the bonding is monodentate or bidentate. We have tested this principle with ab initio modeling for aqueous metal carboxylate complexes and have shown that it does indeed hold. Modeling of the bonding of acetate and formate in aqueous solution to a range of cations was used to predict the infrared spectra of the metal‐carboxylate complexes, and the wavenumbers of the symmetric and antisymmetric vibrational modes are reported. Furthermore, we have shown that these shifts in wavenumber occur primarily due to how bonding with the metal changes the carboxylate C?O bond lengths and O‐C‐O angle.     

Carbamoylphosphine oxide complexes of trivalent lanthanide cations: role of counterions, ligand binding mode, and protonation investigated by quantum mechanical calculations

We present a quantum mechanical study of carbamoylphosphine oxide (CMPO) complexes of MX(3) (M(3+) = La(3+), Eu(3+), Yb(3+); X(-) = Cl(-), NO(3)(-)) with a systematic comparison of monodentate vs bidentate binding modes of CMPO. The per ligand interaction energies Delta E increase from La(3+) to Yb(3+) and are higher with Cl(-) than with NO(3)(-) as counterions, as a result of steric strain in the first coordination sphere with the bidentate anions. The energy difference between monodentate (via phosphoryl oxygen) and bidentate CMPO complexes is surprisingly small, compared to Delta E or to the binding energy of one solvent molecule. Protonation of uncomplexed CMPO takes place preferably at the phosphoryl oxygen O(P), while in the Eu(NO(3))(3)CMPOH(+) complex carbonyl (O(C)) protonation is preferred and O(P) is bonded to the metal. A comparison of uranyl and lanthanide nitrate complexes of CMPO shows that the interaction energies Delta E of the former are lower. Finally, the effect of grafting CMPO arms at the wide rim of a calix[4]arene platform is described. The results are important for our understanding of cation binding and extraction by potentially bidentate CMPO, diamide, and diphosphoryl types of ligands.     

Carboxylate binding to be(2+) in proteins and influence of the dielectric environment

To gain insight into the interaction of Be2+ ions with negatively charged protein residues, the free energy changes associated with the replacement of water molecules in the first hydration shell of with one and two acetate anions were computed for the gas phase reactions using ab initio methods at the MP2 and DFT-B3LYP computational levels. Both unidentate and bidentate modes of coordination of the carboxylate group with the Be2+ ion are considered. Continuum dielectric calculations were then performed to estimate the corresponding free energy changes in several environments of varying dielectric strength. Environments with dielectric constants of 2 and 4, which represent a protein interior, and 78, which corresponds to water, were used. It is found that the free energy changes for the substitution reactions decrease in magnitude with increasing dielectric strength, in agreement with similar results reported for Mg2+, Ca2+, and Zn2+ (Dudev et al. J. Phys. Chem. B 2000, 104, 3692). However, unlike Mg2+, Ca2+, and Zn2+, the free energy change for single-anion or concerted two-anion substitution reactions with remains negative and indicates the reactions are still favorable in the high dielectric aqueous environment. It is also found that the unidentate mode of binding is favored over the bidentate mode, and this is attributed, in part, to the introduction of hydrogen bonds between one carboxylate oxygen and a water molecule within the cluster when unidentate binding with Be2+ is involved.     

Effect of metal cations [Li+, Na+, K+, Be2+, Mg2+, and Ca2+] on the structure of 2‐(3′‐hydroxy‐2′‐pyridyl)benzoxazole: A theoretical investigation

Density functional theory calculations were performed at the B3LYP/6‐311++G(d,p) level to systematically explore the geometrical multiplicity and binding strength for the complexes formed by alkaline and alkaline earth metal cations, viz. Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺ (Mⁿ⁺, hereinafter), with 2‐(3′‐hydroxy‐2′‐pyridyl)benzoxazole. A total of 60 initial structures were designed and optimized, of which 51 optimized structures were found, which could be divided into two different types: monodentate complexes and bidentate complexes. In the cation‐heteroatom complex, bidentate binding is generally stronger than monodentate binding, and of which the bidentate binding with five‐membered ring structure has the strongest interaction. Energy decomposition revealed that the total binding energies mainly come from electrostatic interaction for alkaline metal ion complexes and orbital interaction energy for alkaline earth metal ion complex. In addition, the electron localization function analysis show that only the Be? O and Be? N bond are covalent character, and others are ionic character. © 2012 Wiley Periodicals, Inc.     

Uranyl complexes with diamide ligands: a quantum mechanics study of chelating properties in the gas phase

We report a quantum mechanical study on the complexes of UO(2)(2+) with diamide ligands L of malonamide and succinamide type, respectively, forming 6- and 7-chelate rings in their bidentate coordination to uranium. The main aims are to (i) assess how strong the chelate effect is (i.e., the preference for bi- versus monodentate binding modes of L), (ii) compare these ligands as a function of the chelate ring size, and (iii) assess the role of neutralizing counterions. For this purpose, we consider UO(2)L(2+), UO(2)L(2)(2+), UO(2)L(3)(2+), and UO(2)X(2)L type complexes with X(-) = Cl(-) versus NO(3)(-). Hartree-Fock and DFT calculations lead to similar trends and reveal the importance of saturation and steric repulsions ("strain") in the first coordination sphere. In the unsaturated UO(2)L(2+), UO(2)L(2)(2+), and UO(2)Cl(2)L complexes, the 7-ring chelate is preferred over the 6-ring chelate, and bidentate coordination is preferred over the monodentate one. However, in the saturated UO(2)(NO(3))(2)L complexes, the 6- and 7-chelating ligands have similar binding energies, and for a given ligand, the mono- and bidentate binding modes are quasi-isoenergetic. These conclusions are confirmed by the calculations of free energies of complexation in the gas phase. In condensed phases, the monodentate form of UO(2)X(2)L complexes should be further stabilized by coordination of additional ligands, as well as by interactions (e.g., hydrogen bonding) of the "free" carbonyl oxygen, leading to an enthalpic preference for this form, compared to the bidentate one. We also considered an isodesmic reaction exchanging one bidentate ligand L with two monoamide analogues, which reveals that the latter are clearly preferred (by 23-14 kcal/mol at the HF level and 24-12 kcal/mol at the DFT level). Thus, in the gas phase, the studied bidentate ligands are enthalpically disfavored, compared to bis-monodentate analogues. The contrast with trends observed in solution hints at the importance of "long range" forces (e.g., second shell interactions) and entropy effects on the chelate effect in condensed phases.     

Dispersion interactions of carbohydrates with condensate aromatic moieties: theoretical study on the CH-π interaction additive properties

In this article we present the first systematic study of the additive properties (i.e. degree of additivity) of the carbohydrate-aromatic moiety CH-π dispersion interaction. The additive properties were studied on the β-D-glucopyranose, β-D-mannopyranose and α-L-fucopyranose complexes with the naphthalene molecule by comparing the monodentate (single CH-π) and bidentate (two CH-π) complexes. All model complexes were optimized using the DFT-D approach, at the BP/def2-TZVPP level of theory. The interaction energies were refined using single point calculations at highly correlated ab initio methods at the CCSD(T)/CBS level, calculated as E + (E(CCSD(T))-E(MP2))(Small Basis). Bidentate complexes show very strong interactions in the range from -10.79 up to -7.15 and -8.20 up to -6.14 kcal mol(-1) for the DFT-D and CCSD(T)/CBS level, respectively. These values were compared with the sum of interaction energies of the appropriate monodentate carbohydrate-naphthalene complexes. The comparison reveals that the bidentate complex interaction energy is higher (interaction is weaker) than the sum of monodentate complex interaction energies. Bidentate complex interaction energy corresponds to 2/3 of the sum of the appropriate monodentate complex interaction energies (averaging over all modeled carbohydrate complexes). The observed interaction energies were also compared with the sum of interaction energies of the corresponding previously published carbohydrate-benzene complexes. Also in this case the interaction energy of the bidentate complex was higher (i.e. weaker interaction) than the sum of interaction energies of the corresponding benzene complexes. However, the obtained difference is lower than before, while the bidentate complex interaction energy corresponds to 4/5 of the sum of interaction energy of the benzene complexes, averaged over all structures. The mentioned comparison might aid protein engineering efforts where amino acid residues phenylalanine or tyrosine are to be replaced by a tryptophan and can help to predict the changes in the interactions. The observed results also show that DFT-D correctly describes the CH-π interaction energy and their additive properties in comparison to CCSD(T)/CBS calculated interaction energies. Thus, the DFT-D approach might be used for calculation of larger complexes of biological interest, where dispersion interaction plays an important role.     

Representation of Zn(II) complexes in polarizable molecular mechanics. Further refinements of the electrostatic and short-range contributions. Comparisons with parallel ab initio computations

We present refinements of the SIBFA molecular mechanics procedure to represent the intermolecular interaction energies of Zn(II). The two first-order contributions, electrostatic (E(MTP)), and short-range repulsion (E(rep)), are refined following the recent developments due to Piquemal et al. (Piquemal et al. J Phys Chem A 2003, 107, 9800; and Piquemal et al., submitted). Thus, E(MTP) is augmented with a penetration component, E(pen), which accounts for the effects of reduction in electronic density of a given molecular fragment sensed by another interacting fragment upon mutual overlap. E(pen) is fit in a limited number of selected Zn(II)-mono-ligated complexes so that the sum of E(MTP) and E(pen) reproduces the Coulomb contribution E(c) from an ab initio Hartree-Fock energy decomposition procedure. Denoting by S, the overlap matrix between localized orbitals on the interacting monomers, and by R, the distance between their centroids, E(rep) is expressed by a S(2)/R term now augmented with an S(2)/R(2) one. It is calibrated in selected monoligated Zn(II) complexes to fit the corresponding exchange repulsion E(exch) from ab initio energy decomposition, and no longer as previously the difference between (E(c) + E(exch)) and E(MTP). Along with the reformulation of the first-order contributions, a limited recalibration of the second-order contributions was carried out. As in our original formulation (Gresh, J Comput Chem 1995, 16, 856), the Zn(II) parameters for each energy contribution were calibrated to reproduce the radial behavior of its ab initio HF counterpart in monoligated complexes with N, O, and S ligands. The SIBFA procedure was subsequently validated by comparisons with parallel ab initio computations on several Zn(II) polyligated complexes, including binuclear Zn(II) complexes as in models for the Gal4 and beta-lactamase metalloproteins. The largest relative error with respect to the RVS computations is 3%, and the ordering in relative energies of competing structures reproduced even though the absolute numerical values of the ab initio interaction energies can be as large as 1220 kcal/mol. A term-to-term identification of the SIBFA contributions to their ab initio counterparts remained possible even for the largest sized complexes.     

General molecular mechanics method for transition metal carboxylates and its application to the multiple coordination modes in mono- and dinuclear Mn(II) complexes

A general molecular mechanics method is presented for modeling the symmetric bidentate, asymmetric bidentate, and bridging modes of metal-carboxylates with a single parameter set by using a double-minimum M-O-C angle-bending potential. The method is implemented within the Molecular Operating Environment (MOE) with parameters based on the Merck molecular force field although, with suitable modifications, other MM packages and force fields could easily be used. Parameters for high-spin d (5) manganese(II) bound to carboxylate and water plus amine, pyridyl, imidazolyl, and pyrazolyl donors are developed based on 26 mononuclear and 29 dinuclear crystallographically characterized complexes. The average rmsd for Mn-L distances is 0.08 A, which is comparable to the experimental uncertainty required to cover multiple binding modes, and the average rmsd in heavy atom positions is around 0.5 A. In all cases, whatever binding mode is reported is also computed to be a stable local minimum. In addition, the structure-based parametrization implicitly captures the energetics and gives the same relative energies of symmetric and asymmetric coordination modes as density functional theory calculations in model and "real" complexes. Molecular dynamics simulations show that carboxylate rotation is favored over "flipping" while a stochastic search algorithm is described for randomly searching conformational space. The model reproduces Mn-Mn distances in dinuclear systems especially accurately, and this feature is employed to illustrate how MM calculations on models for the dimanganese active site of methionine aminopeptidase can help determine some of the details which may be missing from the experimental structure.     

Interaction of neutral and zwitterionic glycine with Zn2+ in gas phase: ab initio and SIBFA molecular mechanics calculations

The interaction of Zn²⁺ with glycine (Gly) in the gas phase is studied by a combination of ab initio and molecular mechanics techniques. The structures and energetics of the various isomers of the Gly–Zn²⁺ complex are first established via high‐level ab initio calculations. Two low‐energy isomers are characterized: one in which the metal ion interacts with the carboxylate end of zwitterionic glycine, and another in which it chelates the amino nitrogen and the carbonyl oygen of neutral glycine. These calculations lead to the first accurate value of the gas‐phase affinity of glycine for Zn²⁺. Ab initio calculations were also used to evaluate the performance of various implementations of the SIBFA force field. To assess the extent of transferability of the distributed multipoles and polarizabilities used in the SIBFA computations, two approaches are followed. In the first, approach (a), these quantities are extracted from the ab initio Hartree–Fock wave functions of glycine or its zwitterion in its entirety, and for each individual Zn²⁺‐binding conformation. In the second, approach (b), they are assembled from the appropriate constitutive fragments, namely methylamine and formic acid for neutral glycine, and protonated methylamine and formate for the zwitterion; they undergo the appropriate vector or matrix rotation to be assembled in the conformation studied. The values of the Zn²⁺–glycine interaction energies are compared to those resulting from ab initio SCF and MP2 computations using both the all‐electron 6‐311+G(2d,2p) basis set and an effective core potential together with the valence CEP 4‐31G(2d) basis set. Approach (a) values closely reproduce the ab initio ones, both in terms of the total interaction energies and of the individual components. Approach (b) can provide a similar match to ab initio interaction energies as does approach (a), provided that the two constitutive Gly building blocks are considered as separate entities having mutual interactions that are computed simultaneously with those occurring with Zn²⁺. Thus, the supermolecule is treated as a three‐body rather than a two‐body system. These results indicate that the current implementation of the SIBFA force field should be adequate to undertake accurate studies on zinc metallopeptides. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 963–973, 2000     

Coordination number of zinc ions in the phosphotriesterase active site by molecular dynamics and quantum mechanics

We have run several molecular dynamics (MD) simulations on zinc-containing phosphotriesterase (PTE) with two bound substrates, sarin and paraoxon, and with the substrate analog diethyl 4-methylbenzylphosphonate. A standard nonbonded model was employed to treat the zinc ions with the commonly used charge of +2. In all the trajectories, we observed a tightly bound water (TBW) molecule in the active site that was coordinated to the less buried zinc ion. The phosphoryl oxygen of the substrate/inhibitor was found to be coordinated to the same zinc ion so that, considering all ligands, the less buried zinc was hexa-coordinated. The hexa-coordination of this zinc ion was not seen in the deposited X-ray pdb files for PTE. Several additional MD simulations were then performed using different charges (+1, +1.5) on the zinc ions, along with ab initio and density functional theory (DFT) calculations, to evaluate the following possibilities: the crystal diffraction data were not correctly interpreted; the hexa-coordinated zinc ion in PTE is only present in solution and not in the crystal; and the hexa-coordinated zinc ion in PTE is an artifact of the force field used. A charge of +1.5 leads to a coordination number (CN) of 5 on both zinc ions, which is consistent with the results from ab initio and DFT calculations and with the latest high resolution X-ray crystal structure. The commonly used charge of +2 produces a CN of 6 on the less buried zinc. The CN on the more buried zinc ion is 5 when the substrate/inhibitor is present in the simulation, and increases to 6 when the substrate/inhibitor is removed prior to the simulation. The results of both of the MD and quantum mechanical calculations lead to the conclusion that the zinc ions in the PTE active site are both penta-coordinated, and that the MD simulations performed with the charge of +2 overestimate the CN of the zinc ions in the PTE active site. The overall protein structures in the simulations remain unaffected by the change in zinc charge from +2 to +1.5. The results also suggest that the charge +1.5 is the most appropriate for the molecular dynamics simulations on zinc-containing PTE when a nonbonded model is used and no global thermodynamic conclusion is sought. We also show that the standard nonbonded model is not able to properly treat the CN and energy at the same time. A preliminary, promising charge-transfer model is discussed with the use of the zinc charge of +1.5.     

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.     

Correlation of infrared spectra of zinc(II) carboxylates with their structures

The correlation of the infrared spectra of zinc(II) carboxylates with their structures was investigated in the paper. The complexes with different modes of the carboxylate binding, from chelating, through bridging (syn-syn, syn-anti, monatomic), ionic to monodentate were used for the study, namely [Zn(C6H5CHCHCOO)2(H2O)2] (I) with chelating carboxylate group (C6H5CHCHCOO=cinnamate), [Zn2(C6H5COO)4(pap)2] (II) with syn-syn bridging carboxylate (C6H5COO=benzoate; pap=papaverine), [Zn(C6H5CHCHCOO)2(mpcm)]n (III) with syn-anti carboxylate bridge (mpcm=methyl-3-pyridylcarbamate), [Zn(C5H4NCOO)2(H2O)4] (IV) with ionic carboxylate group (C5H4NCOO=nicotinate), [Zn(C6H5COO)2(pcb)2]n (V) with monodentate carboxylate coordination (pcb=3-pyridylcarbinol) and [Zn3(C6H5COO)6(nia)2] (VI) with syn-syn and monatomic carboxylate bridges (nia=nicotinamide). First, the mode of the carboxylate binding was assigned from the infrared spectra using the magnitude of the separation between the carboxylate stretches, Deltaexp=nuas(COO-)-nus(COO-). Then the values Deltaexp were compared with those calculated from structural data of the carboxylate anion (Deltacalc). The conclusions about the carboxylate binding which resulted from the Delta values, were confronted with the crystal structure of the complexes. The limitations and recommendations were formulated to assign the mode of the carboxylate binding from the infrared spectra. The dependence of the Deltaexp values on the magnitudes of Zn-O-C angles in bidentate carboxylate coordination was observed.     

Synthesis, thermal and infrared spectroscopic studies of hydrazinated mixed metal fumarates

A series of hydrazinated complexes of mixed nickel manganese zinc ferrous fumarates were synthesized from aqueous metal chlorides solution and sodium fumarate–hydrazine hydrate mixture. These complexes were characterized by chemical analysis, infrared spectroscopy, thermogravimetry-differential scanning calorimetry and isothermal mass loss studies. The hydrazine ligand in these complexes shows bidentate bridging nature while, carboxylate ion displays monodentate behaviour. Thermal decomposition studies indicate 3-step decomposition, involving 2-step dehydrazination and 1-step decarboxylation. The thermal decomposition residues were characterized by X-ray diffractometry, transmission electron microscopy and infrared spectroscopy.     

Mapping the interaction energy surfaces of cyclic alkanes: evaluating the transferability of an ab initio based potential model

Detailed interaction energy maps are computed for symmetric cyclopropane and tetrahedrane dimer systems using ab initio methods. Interaction energies of cubane and cyclohexane dimers are also reported. The global minimum energy structures of cyclopropane and tetrahedrane systems are both D(3d) structures with energies of -1.850 and -2.171 kcal mol(-1). The ability of NIPE potential model, based on ab initio nonbonding data of neopentane (N), isobutane (I), propane (P), ethane (E) and all their combinations to predict the pair interaction energies of these strained cyclic hydrocarbons is also investigated. The difference between the energies predicted by NIPE and those obtained from the ab initio calculations increases with ring strain In general, NIPE values are in close agreement with the ab initio results for alkane ring structures having low ring strain.     

Polarizable interaction potential for water from coupled cluster calculations. I. Analysis of dimer potential energy surface

A six-dimensional interaction potential for the water dimer has been fitted to ab initio interaction energies computed at 2510 dimer configurations. These energies were obtained by combining the supermolecular second-order energies extrapolated to the complete basis set limit from up to quadruple-zeta quality basis sets with the contribution from the coupled-cluster method including single, double, and noniterative triple excitations computed in a triple-zeta quality basis set. All basis sets were augmented by diffuse functions and supplemented by midbond functions. The energies have been fitted using an analytic form with the induction component represented by a polarizable term, making the potential directly transferable to clusters and the bulk phase. Geometries and energies of stationary points on the potential surface agree well with the results of high-level ab initio geometry optimizations.     

Insight into the complex and dynamic process of activation of matrix metalloproteinases.

Matrix metalloproteinases (MMPs) are important hydrolytic enzymes with profound physiological and pathological functions in living organisms. MMPs are produced in their inactive zymogenic forms, which are subsequently proteolytically activated in an elaborate set of events. The propeptide in the zymogen blocks the active site, with a cysteine side-chain thiolate from this propeptide achieving coordination with the catalytically important zinc ion in the active site. Molecular dynamics simulations, ab initio calculations, and wet chemistry experiments presented herein argue for the critical importance of a protonation event at the coordinated thiolate as a prerequisite for the departure of the propeptide from the active site. Furthermore, a catalytically important glutamate is shown to coordinate transiently to the active-site zinc ion to "mask" the positive potential of the zinc ion and lower the energy barrier for dissociation of the protonated cysteine side chain from the zinc ion. In addition, a subtle conformational change by the propeptide is needed in the course of zymogen activation. These elaborate processes take place in concert in the activation process of MMPs, and the insight into these processes presented herein sheds light on a highly regulated physiological process with profound consequences for eukaryotic organisms.