Modeling of inhibitor–metalloenzyme interactions and selectivity using molecular mechanics grounded in quantum chemistry

Temperature effect on the structure and conformational fluctuations in two zinc knuckles from the mouse mammary tumor virus

Zinc fingers are small protein domains in which zinc plays a structural role, contributing to the stability of the zinc-peptide complex. Zinc fingers are structurally diverse and are present in proteins that perform a broad range of functions in various cellular processes, such as replication and repair, transcription and translation, metabolism and signaling, cell proliferation, and apoptosis. Zinc fingers typically function as interaction modules and bind to a wide variety of compounds, such as nucleic acids, proteins, and small molecules. In this study, we investigated the structural properties, in solution, of the proximal and distal zinc knuckles of the nucleocapsid (NC) protein from the mouse mammary tumor virus (MMTV) (MMTV NC). For this purpose, we performed a series of molecular dynamics simulations in aqueous solution at 300 K, 333 K, and 348 K. The temperature effect was evaluated in terms of root mean square deviation of the backbone atoms and root mean square fluctuation of the coordinating residue atoms. The stability of the zinc coordination sphere was analyzed based upon the time profile of the interatomic distances between the zinc ions and the chelator atoms. The results indicate that the hydrophobic character of the proximal zinc finger is dominant at 333 K. The low mobility of the coordinating residues suggests that the strong electrostatic effect exerted by the zinc ion on its coordinating residues is not influenced by the increase in temperature. The evolution of the structural parameters of the coordination sphere of the distal zinc finger at 300 K gives us a reasonable picture of the unfolding pathway, as proposed by Bombarda and coworkers (Bombarda et al., 2005), which can predict the binding order of the four conserved ligand-binding residues. Our results support the conclusion that the structural features can vary significantly between the two zinc knuckles of MMTV NC.

Aromatic Base Stacking in DNA: From ab initio Calculations to Molecular Dynamics Simulations

Abstract Aromatic stacking of nucleic acid bases is one of the key players in determining the structure and dynamics of nucleic acids. The arrangement of nucleic acid bases with extensive overlap of their aromatic rings gave rise to numerous often contradictory suggestions about the physical origins of stacking and the possible role of delocalized electrons in stacked aromatic π systems, leading to some confusion about the issue. The recent advance of computer hardware and software finally allowed the application of state of the art quantum-mechanical approaches with inclusion of electron correlation effects to study aromatic base stacking, now providing an ultimitate qualitative description of the phenomenon. Base stacking is determined by an interplay of the three most commonly encountered molecular interactions: dispersion attraction, electrostatic interaction, and short-range repulsion. Unusual (aromatic- stacking specific) energy contributions were in fact not evidenced and are not necessary to describe stacking. The currently used simple empirical potential form, relying on atom-centered constant point charges and Lennard-Jones van der Waals terms, is entirely able to reproduce the essential features of base stacking. Thus, we can conclude that base stacking is in principle one of the best described interactions in current molecular modeling and it allows to study base stacking in DNA using large-scale classical molecular dynamics simulations. Neglect of cooperativity of stacking appears to be the most serious approximation of the currently used force field form. This review summarizes recent developments in the field. It is written for an audience that is not necessarily expert in computational quantum chemistry and follows up on our previous contribution (Sponer et. al., J. Biomol. Struct. Dyn. 14, 117, (1997)). First, the applied methodology, its accuracy, and the physical nature of base stacking is briefly overviewed, including a comment on the accuracy of other molecular orbital methods and force fields. Then, base stacking is contrasted with hydrogen bonding, the other dominant force in nucleic acid structure. The sequence dependence and cooperativity of base stacking is commented on, and finally a brief introduction into recent progress in large-scale molecular dynamics simulations of nucleic acids is provided. Using four stranded DNA assemblies as an example, we demonstrate the efficacy of current molecular dynamics techniques that utilize refined and verified force fields in the study of stacking in nucleic acid molecules.

Modeling Structural Coordination and Ligand Binding in Zinc Proteins with a Polarizable Potential

As the second most abundant cation in human body, zinc is vital for the structures and functions of many proteins. Zinc-containing matrix metalloproteinases (MMPs) have been widely investigated as potential drug targets in a range of diseases ranging from cardiovascular disorders to cancers. However, it remains a challenge in theoretical studies to treat zinc in proteins with classical mechanics. In this study, we examined Zn(2+) coordination with organic compounds and protein side chains using a polarizable atomic multipole based electrostatic model. We find that polarization effect plays a determining role in Zn(2+) coordination geometry in both matrix metalloproteinase (MMP) complexes and in zinc-finger proteins. In addition, the relative binding free energies of selected inhibitors binding with MMP13 have been estimated and compared with experimental results. While not directly interacting with the small molecule inhibitors, the permanent and polarizing field of Zn(2+) exerts a strong influence on the relative affinities of the ligands. The simulation results also reveal the polarization effect on binding is ligand dependent and thus difficult to be incorporated into fixed-charge models implicitly.

The reaction mechanism of type I phosphomannose isomerases: new information from inhibition and polarizable molecular mechanics studies

Type I phosphomannose isomerases (PMIs) are zinc-dependent metalloenzymes involved in the reversible isomerization of D-mannose 6-phosphate (M6P) and D-fructose 6-phosphate (F6P). 5-Phospho-D-arabinonohydroxamic acid (5PAH), an inhibitor endowed with nanomolar affinity for yeast (Type I) and Pseudomonas aeruginosa (Type II) PMIs (Roux et al., Biochemistry 2004; 43:2926-2934), strongly inhibits human (Type I) PMI (for which we report an improved expression and purification procedure), as well as Escherichia coli (Type I) PMI. Its K(i) value of 41 nM for human PMI is the lowest value ever reported for an inhibitor of PMI. 5-Phospho-D-arabinonhydrazide, a neutral analogue of the reaction intermediate 1,2-cis-enediol, is about 15 times less efficient at inhibiting both enzymes, in accord with the anionic nature of the postulated high-energy reaction intermediate. Using the polarizable molecular mechanics, sum of interactions between fragments ab initio computed (SIBFA) procedure, computed structures of the complexes between Candida albicans (Type I) PMI and the cyclic substrate β-D-mannopyranose 6-phosphate (β-M6P) and between the enzyme and the high-energy intermediate analogue inhibitor 5PAH are reported. Their analysis allows us to identify clearly the nature of each individual active site amino acid and to formulate a hypothesis for the overall mechanism of the reaction catalyzed by Type I PMIs, that is, the ring-opening and isomerization steps, respectively. Following enzyme-catalyzed ring-opening of β-M6P by zinc-coordinated water and Gln111 ligands, Lys136 is identified as the probable catalytic base involved in proton transfer between the two carbon atoms C1 and C2 of the substrate D-mannose 6-phosphate.

Synthesis and evaluation of non-hydrolyzable D-mannose 6-phosphate surrogates reveal 6-deoxy-6-dicarboxymethyl-D-mannose as a new strong inhibitor of phosphomannose isomerases

Non-hydrolyzable d-mannose 6-phosphate analogues in which the phosphate group was replaced by a phosphonomethyl, a dicarboxymethyl, or a carboxymethyl group were synthesized and kinetically evaluated as substrate analogues acting as potential inhibitors of type I phosphomannose isomerases (PMIs) from Saccharomyces cerevisiae and Escherichia coli. While 6-deoxy-6-phosphonomethyl-d-mannose and 6-deoxy-6-carboxymethyl-D-mannose did not inhibit the enzymes significantly, 6-deoxy-6-dicarboxymethyl-D-mannose appeared as a new strong competitive inhibitor of both S. cerevisiae and E. coli PMIs with K(m)/K(i) ratios of 28 and 8, respectively. We thus report the first malonate-based inhibitor of an aldose-ketose isomerase to date. Phosphonomethyl mimics of the 1,2-cis-enediolate high-energy intermediate postulated for the isomerization reaction catalyzed by PMIs were also synthesized but behave as poor inhibitors of PMIs. A polarizable molecular mechanics (SIBFA) study was performed on the complexes of d-mannose 6-phosphate and two of its analogues with PMI from Candida albicans, an enzyme involved in yeast infection homologous to S. cerevisiae and E. coli PMIs. It shows that effective binding to the catalytic site occurs with retention of the Zn(II)-bound water molecule. Thus the binding of the hydroxyl group on C1 of the ligand to Zn(II) should be water-mediated. The kinetic study reported here also suggests the dianionic character of the phosphate surrogate as a likely essential parameter for strong binding of the inhibitor to the enzyme active site.

Discovery of potent, selective, and orally active carboxylic acid based inhibitors of matrix metalloproteinase-13

The matrix metalloproteinase enzyme MMP-13 plays a key role in the degradation of type II collagen in cartilage and bone in osteoarthritis (OA). An effective MMP-13 inhibitor would therefore be a novel disease modifying therapy for the treatment of arthritis. Our efforts have resulted in the discovery of a series of carboxylic acid inhibitors of MMP-13 that do not significantly inhibit the related MMP-1 (collagenase-1) or tumor necrosis factor-alpha (TNF-alpha) converting enzyme (TACE). It has previously been suggested (but not proven) that inhibition of the latter two enzymes could lead to side effects. A promising carboxylic acid lead 9 was identified and a convergent synthesis developed. This paper describes the optimization of 9 and the identification of a compound 24f for further development. Compound 24f is a subnanomolar inhibitor of MMP-13 (IC(50) value 0.5 nM and K(i) of 0.19 nM) having no activity against MMP-1 or TACE (IC(50) of >10000 nM). Furthermore, in a rat model of MMP-13-induced cartilage degradation, 24f significantly reduced proteoglycan release following oral dosing at 30 mg/kg (75% inhibition, p < 0.05) and at 10 mg/kg (40% inhibition, p < 0.05).

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

Competition between protein ligands and cytoplasmic inorganic anions for the metal cation: a DFT/CDM study

Many of the essential metalloproteins are located in the cell, whose cytoplasmic fluid contains several small inorganic anions, such as Cl-, NO2-, NO3-, H2PO4-, and SO4(2-), that play an indispensable role in determining the cell's volume, regulating the cell's pH, signal transduction, muscle contraction, as well as cell growth and metabolism. However, the physical principles governing the competition between these abundant, intracellular anions and protein or nucleic acid residues in binding to cytoplasmic metal cations such as Na+, K+, Mg2+, and Ca2+ are not well understood; hence, we have delineated the physicochemical basis for this competition using density functional theory in conjunction with the continuum dielectric method. The results show that the metal cation can bind to its target protein against a high background concentration of inorganic anions because (i) desolvating a negatively charged Asp/Glu carboxylate in a protein cavity costs much less than desolvating an inorganic anion in aqueous solution and (ii) the metal-binding site acts as a polydentate ligand that uses all its ligating entities to bind the metal cation either directly or indirectly. The results also show that the absolute hydration free energy of the "alien" anion as well as the net charge and relative solvent exposure of the metal-binding protein cavity are the key factors governing the competition between protein and inorganic ligands for a given cytoplasmic metal cation. Increasing the net negative charge of the protein cavity, while decreasing the number of available amide groups for metal binding, protects the metal-bound ligands from being dislodged by cellular anions, thus revealing a "protective" role for carboxylate groups in a protein cavity, in addition to their role in high affinity metal-binding.

Binding of 5-phospho-D-arabinonohydroxamate and 5-phospho-D-arabinonate inhibitors to zinc phosphomannose isomerase from Candida albicans studied by polarizable molecular mechanics and quantum mechanics

Type I phosphomannose isomerase (PMI) is a Zn-dependent metalloenzyme involved in the isomerization of D-fructose 6-phosphate to D-mannose 6-phosphate. One of our laboratories has recently designed and synthesized 5-phospho-D-arabinonohydroxamate (5PAH), an inhibitor endowed with a nanomolar affinity for PMI (Roux et al., Biochemistry 2004, 43, 2926). By contrast, the 5-phospho-D-arabinonate (5PAA), in which the hydroxamate moiety is replaced by a carboxylate one, is devoid of inhibitory potency. Subsequent biochemical studies showed that in its PMI complex, 5PAH binds Zn(II) through its hydroxamate moiety rather than through its phosphate. These results have stimulated the present theoretical investigation in which we resort to the SIBFA polarizable molecular mechanics procedure to unravel the structural and energetical aspects of 5PAH and 5PAA binding to a 164-residue model of PMI. Consistent with the experimental results, our theoretical studies indicate that the complexation of PMI by 5PAH is much more favorable than by 5PAA, and that in the 5PAH complex, Zn(II) ligation by hydroxamate is much more favorable than by phosphate. Validations by parallel quantum-chemical computations on model of the recognition site extracted from the PMI-inhibitor complexes, and totaling up to 140 atoms, showed the values of the SIBFA intermolecular interaction energies in such models to be able to reproduce the quantum-chemistry ones with relative errors < 3%. On the basis of the PMI-5PAH SIBFA energy-minimized structure, we report the first hypothesis of a detailed view of the active site of the zinc PMI complexed to the high-energy intermediate analogue inhibitor, which allows us to identify active site residues likely involved in the proton transfer between the two adjacent carbons of the substrates.

Computational protocol for predicting the binding affinities of zinc containing metalloprotein-ligand complexes

Zinc is one of the most important metal ions found in proteins performing specific functions associated with life processes. Coordination geometry of the zinc ion in the active site of the metalloprotein-ligand complexes poses a challenge in determining ligand binding affinities accurately in structure-based drug design. We report here an all atom force field based computational protocol for estimating rapidly the binding affinities of zinc containing metalloprotein-ligand complexes, considering electrostatics, van der Waals, hydrophobicity, and loss in conformational entropy of protein side chains upon ligand binding along with a nonbonded approach to model the interactions of the zinc ion with all the other atoms of the complex. We examined the sensitivity of the binding affinity predictions to the choice of Lennard-Jones parameters, partial atomic charges, and dielectric treatments adopted for system preparation and scoring. The highest correlation obtained was R2 = 0.77 (r = 0.88) for the predicted binding affinity against the experiment on a heterogenous dataset of 90 zinc containing metalloprotein-ligand complexes consisting of five unique protein targets. Model validation and parameter analysis studies underscore the robustness and predictive ability of the scoring function. The high correlation obtained suggests the potential applicability of the methodology in designing novel ligands for zinc-metalloproteins. The scoring function has been web enabled for free access at www.scfbio-iitd.res.in/software/drugdesign/bapplz.jsp as BAPPL-Z server (Binding Affinity Prediction of Protein-Ligand complexes containing Zinc metal ions).

An ab initio investigation of zinc chloro complexes

A series of geometry, frequency, and energy calculations of chloroaquazinc(II) complexes were carried out at up to the MP2/6-31+G* level. A thorough examination of all species up to and including hexacoordinate species, and with up to six chlorides, was carried out. The structures of the complexes are compared with experimental data where available. The solution chemistry of zinc(II) in the presence of chloride is discussed, and Raman spectra of zinc perchlorate with increasing amount of chloride are presented.

Molecular docking study and development of an empirical binding free energy model for phosphodiesterase 4 inhibitors

In the present work, several computational methodologies were combined to develop a model for the prediction of PDE4B inhibitors' activity. The adequacy of applying the ligand docking approach, keeping the enzyme rigid, to the study of a series of PDE4 inhibitors was confirmed by a previous molecular dynamics analysis of the complete enzyme. An exhaustive docking procedure was performed to identify the most probable binding modes of the ligands to the enzyme, including the active site metal ions and the surrounding structural water molecules. The enzyme-inhibitor interaction enthalpies, refined by using the semiempirical molecular orbital approach, were combined with calculated solvation free energies and entropy considerations in an empirical free energy model that enabled the calculation of binding free energies that correlated very well with experimentally derived binding free energies. Our results indicate that both the inclusion of the structural water molecules close to the ions in the binding site and the use of a free energy model with a quadratic dependency on the ligand free energy of solvation are important aspects to be considered for molecular docking investigations involving the PDE4 enzyme family.

Complexes of thiomandelate and captopril mercaptocarboxylate inhibitors to metallo-β-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.

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.

CHARMM fluctuating charge force field for proteins: II protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model

A fluctuating charge (FQ) force field is applied to molecular dynamics simulations for six small proteins in explicit polarizable solvent represented by the TIP4P-FQ potential. The proteins include 1FSV, 1ENH, 1PGB, 1VII, 1H8K, and 1CRN, representing both helical and beta-sheet secondary structural elements. Constant pressure and temperature (NPT) molecular dynamics simulations are performed on time scales of several nanoseconds, the longest simulations yet reported using explicitly polarizable all-atom empirical potentials (for both solvent and protein) in the condensed phase. In terms of structure, the FQ force field allows deviations from native structure up to 2.5 A (with a range of 1.0 to 2.5 A). This is commensurate to the performance of the CHARMM22 nonpolarizable model and other currently existing polarizable models. Importantly, secondary structural elements maintain native structure in general to within 1 A (both helix and beta-strands), again in good agreement with the nonpolarizable case. In qualitative agreement with QM/MM ab initio dynamics on crambin (Liu et al. Proteins 2001, 44, 484), there is a sequence dependence of average condensed phase atomic charge for all proteins, a dependence one would anticipate considering the differing chemical environments around individual atoms; this is a subtle quantum mechanical feature captured in the FQ model but absent in current state-of-the-art nonpolarizable models. Furthermore, there is a mutual polarization of solvent and protein in the condensed phase. Solvent dipole moment distributions within the first and second solvation shells around the protein display a shift towards higher dipole moments (increases on the order of 0.2-0.3 Debye) relative to the bulk; protein polarization is manifested via the enhanced condensed phase charges of typical polar atoms such as backbone carbonyl oxygens, amide nitrogens, and amide hydrogens. Finally, to enlarge the sample set of proteins, gas-phase minimizations and 1 ps constant temperature simulations are performed on various-sized proteins to compare to earlier work by Kaminsky et al. (J Comp Chem 2002, 23, 1515). The present work establishes the feasibility of applying a fully polarizable force field for protein simulations and demonstrates the approach employed in extending the CHARMM force field to include these effects.

Intramolecular interaction energies in model alanine and glycine tetrapeptides. Evaluation of anisotropy, polarization, and correlation effects. A parallel ab initio HF/MP2, DFT, and polarizable molecular mechanics study

An extension of the SIBFA polarizable molecular mechanics procedure to flexible oligopeptides is reported. The procedure is evaluated by computing the relative conformational energies, deltaE(conf), of the alanine tetrapeptide in 10 representative conformations, which were originally derived by Beachy et al. (J Am Chem Soc 1997, 119, 5908) to benchmark molecular mechanics procedures with respect to ab initio computations. In the present study, a particular emphasis is on the separable nature of the components of the energy and the particular impact of the polarization energy component on deltaE(conf). We perform comparisons with respect to single-point HF, DFT, LMP2, and MP2 computations done at the SIBFA-derived energy minima. Such comparisons are made first for the 10 conformers derived from phi/psi torsional angle energy-minimization (the rigid rotor approach), and, in a second step, after allowing additional relaxation of the C(alpha) centered valence angles. In both series of energy-minimization, the SIBFA deltaE(conf) compared best with the LMP2 results using the 6-311G** basis set, the rms being 1.3 kcal/mol. In the absence of the polarization component, the rms is 3.5 kcal/mol. In both series of minimizations, the magnitudes of deltaE(conf), computed as differences with respect to the most stable conformer taken as energy zero, decrease along the series: HF > DFT > LMP2 > SIBFA > MP2, indicative of increasing stabilization of the most highly folded conformers.

Conformation-dependent intermolecular interaction energies of the triphosphate anion with divalent metal cations. Application to the ATP-binding site of a binuclear bacterial enzyme. A parallel quantum chemical and polarizable molecular mechanics investigation

We have explored the conformation-dependent interaction energy of the triphosphate moiety, a key constituent of ATP and GTP, with a closed-shell divalent cation, Zn2+, used as a probe. This was done using the SIBFA polarizable molecular mechanics procedure. We have resorted to a previously developed approach in which triphosphate is built out from its elementary constitutive fragments, and the intramolecular, interfragment, interaction energies are computed simultaneously with their intermolecular interactions with the divalent cation. This approach has enabled reproduction of the values of the intermolecular interaction energies from ab initio quantum-chemistry with relative errors <3%. It was extended to the complex of a nonhydrolyzable analog of ATP with the active site of a bacterial enzyme having two Mg2+ cations as cofactors. We obtained following energy-minimization a very close overlap of the ATP analog over its position from X-ray crystallography. For models of the ATP analog-enzyme complex encompassing up to 169 atoms, the values of the SIBFA interaction energies were found to match their DFT counterparts with relative errors of <2%.

Inhibition of Type I and Type II Phosphomannose Isomerases by the Reaction Intermediate Analogue 5-Phospho-d-Arabinonohydroxamic Acid Supports a Catalytic Role for the Metal Cofactor

The phosphomannose isomerases (PMI) comprise three families of proteins: type I, type II, and type III PMIs. Members of all three families catalyze the reversible isomerization of D-mannose 6-phosphate (M6P) and D-fructose 6-phosphate (F6P) but share little or no sequence identity. Because (1) PMIs are essential for the survival of several microorganisms, including yeasts and bacteria, and (2) the PMI enzymes from several pathogens do not share significant sequence identity to the human protein, PMIs have been considered as potential therapeutic targets. Elucidation of the catalytic and regulatory mechanisms of the different types of PMIs is strongly needed for rational species-specific drug design. To date, inhibition and crystallographic studies of all PMIs are still largely unexplored. As part of our research program on aldose-ketose isomerases, we report in this paper the evaluation of two new inhibitors of type I and type II PMIs from baker's yeast and Pseudomonas aeruginosa, respectively. We found that 5-phospho-D-arabinonohydroxamic acid (5PAH), which is the most potent inhibitor of phosphoglucose isomerase (PGI), is by far the best inhibitor ever reported of both type I and type II PMI-catalyzed isomerization of M6P to F6P. 5PAH, which has an inhibition constant at least 3 orders of magnitude smaller than that of previously reported PMI inhibitors, may be the first high-energy intermediate analogue inhibitor of the enzymes. We also tested the related molecule 5-phospho-D-arabinonate (5PAA), which is a strong competitive inhibitor of PGI, and found that it does not inhibit either PMI. All together, our results are consistent with a catalytic role for the metal cofactor in PMI activity.

Electrostatic properties in the catalytic site of papain: A possible regulatory mechanism for the reactivity of the ion pair

We present an analysis of the electrostatic properties in the catalytic site of papain (EC 3.4.22.2), an archetype enzyme of the C1 cysteine proteinase family, and we investigate their possible role in the formation, stabilization and regulation of the Cys25((-))...His159((+)) catalytic ion pair. The electrostatic properties were computed using a reassociation method based in multicentered multipolar expansions obtained from ab initio quantum calculations of overlapping protein fragments. Solvent effects were introduced by coupling the use of multicentered multipolar expansions to two continuum boundary element methods to solve the Poisson and the linearized Poisson-Boltzmann equations. The electrostatic profile found in the proton transfer region of papain showed that this enzyme has a well-defined electrostatic environment to favor the formation and stabilization of the catalytic ion pair. The papain catalytic site electrostatic profile can be considered as an electrostatic fingerprint of the papain family with the following characteristics: (i) the presence of a net electric field highly aligned in the (Cys25)-SG-->(His159)-ND1 direction; (ii) the electrostatic profile has a saddle-point character; (iii) it is basically a local environmental effect. Furthermore, our analysis describes a possible regulatory mechanism (the E(SG-->ND1) attenuation effect) controlling the ion pair reactivity and permits to infer the Asp57 acidic residue as the most probable candidate to act as the electrostatic modulator.

Coordination and thermodynamics of stable Zn(II) complexes in the gas phase

Ab initio molecular orbital methods in combination with DFT calculations were used to study the structural and thermodynamic properties of 17 complexes containing zinc cation and four first-shell ligands as models of active site of metalloenzymes (e.g. angiotensin converting enzyme, thermolysin). The geometry of the complexes was relaxed by complete optimization by ab initio molecular orbital methods at Hertree-Fock level with 3-21G* basis set. Following single point calculation with tight SCF criteria at the B3LYP level with 6-311+G(2d,p) basis set was used to calculate accurate interaction enthalpies. The structure and thermodynamics of optimized complexes are discussed from the point of view of their biological importance.

Binding of D- and L-captopril inhibitors to metallo-beta-lactamase studied by polarizable molecular mechanics and quantum mechanics

The bacterial Zn2+ metallo-beta-lactamase from B. fragilis is a zinc-enzyme with two potential metal ion binding sites. It cleaves the lactam ring of antibiotics, thus contributing to the acquired resistance of bacteria against antibiotics. The present study bears on the binuclear form of the enzyme. We compare several possible binding modes of captopril, a mercaptocarboxamide inhibitor of several zinc-metalloenzymes. Two diastereoisomers of captopril were considered, with either a D- or an L-proline residue. We have used the polarizable molecular mechanics procedure SIBFA (Sum of Interactions Between Fragments ab initio computed). Two beta-lactamase models were considered, encompassing 104 and 188 residues, respectively. The energy balances included the inter and intramolecular interaction energies as well as the contribution from solvation computed using a continuum reaction field procedure. The thiolate ion of the inhibitor is binding to both metal ions, expelling the bridging solvent molecule from the uncomplexed enzyme. Different competing binding modes of captopril were considered, either where the inhibitor binds in a monodentate mode to the zinc cations only with its thiolate ion, or in bidentate modes involving additional zinc binding by its carboxylate or ketone carbonyl groups. The additional coordination by the inhibitor's carboxylate or carbonyl group always occurs at the zinc ion, which is bound by a histidine, a cysteine, and an aspartate side chain. For both diastereomers, the energy balances favor monodentate binding of captopril via S-. The preference over bidentate binding is small. The interaction energies were recomputed in model sites restricted to captopril, the Zn2+ cations, and their coordinating end side chains from beta-lactamase (98 atoms). The interaction energies and their ranking among competing arrangements were consistent with those computed by ab initio HF and DFT procedures.

Polarizable force fields

Standard force fields used in biomolecular computing describe electrostatic interactions in terms of fixed, usually atom-centered, charges. Real physical systems, however, polarize substantially when placed in a high-dielectric medium such as water--or even when a strongly charged system approaches a neutral body in the gas phase. Such polarization strongly affects the geometry and energetics of molecular recognition. First introduced more than 20 years ago, polarizable force fields seek to account for appropriate variations in charge distribution with dielectric environment. Over the past five years, an accelerated pace of development of such force fields has taken place on systems ranging from liquid water to metalloenzymes. Noteworthy progress has been made in better understanding the capabilities and limitations of polarizable models for water and in the formulation and utilization of complete specifically parameterized polarizable force fields for peptides and proteins.

Structural differences of matrix metalloproteinases. Homology modeling and energy minimization of enzyme-substrate complexes

Matrix metalloproteinases are extracellular enzymes taking part in the remodeling of extracellular matrix. The structures of the catalytic domain of MMP1, MMP3, MMP7 and MMP8 are known, but structures of enzymes belonging to this family still remain to be determined. A general approach to the homology modeling of matrix metalloproteinases, exemplified by the modeling of MMP2, MMP9, MMP12 and MMP14 is described. The models were refined using an energy minimization procedure developed for matrix metalloproteinases. This procedure includes incorporation of parameters for zinc and calcium ions in the AMBER 4.1 force field, applying a non-bonded approach and a full ion charge representation. Energy minimization of the apoenzymes yielded structures with distorted active sites, while reliable three-dimensional structures of the enzymes containing a substrate in active site were obtained. The structural differences between the eight enzyme-substrate complexes were studied with particular emphasis on the active site, and possible sites for obtaining selectivity among the MMP's are discussed. Differences in the P1' pocket are well-documented and have been extensively exploited in inhibitor design. The present work indicates that selectivity could be further improved by considering the P2 pocket as well.

Theoretical Investigations of the Influence of Pressure on the Selectivity of the Michael Addition of Diphenylmethaneamine to Stereogenic Crotonates

SIBFA (sum of interactions between fragments ab initio computed) molecular mechanics systematics has been modified to take into account the effect of pressure on intra- and intermolecular energies. The van der Waals radii are related to the pressure, using Bridgman results on the variation of crystal volume, on one hand, and the relation between the volume of an atom and its van der Waals radius on the other. This procedure produces a decrease of the volume of the systems considered. The modified systematics is used for the study of the conformation at 1 atm and 15 kbar of two stereogenic crotonates and of the complexes formed by these two molecules with the diphenylmethaneamine and the three solvent molecules present in the experiment. The results obtained show that in the case of NMECC 1a the diastereoselectivity induced by high pressure and by the presence of methanol proceeds from an important stabilization of the pro-R reactive complex in which the crotonate has a stacked-transoid conformation. This stabilization is mainly due to intermolecular interactions. In the case of the second crotonate considered, NMCC 1b, our results indicate that pressure induces a stabilization of the pro-R and pro-S complexes having the axial conformation for which the reaction exhibits little diastereoselectivity in qualitative agreement with experimental data. This study tends to show that it is possible to account theoretically for the influence of pressure on molecular conformation and/or complex structure, using a molecular mechanics method that is able to take into account the variation of volumes of the different entities present in the system studied.