Biomedical subjects

Nohad Gresh

Publications and source records attributed to Nohad Gresh.

9 recordsLinked to original sources

PubMed2006

Towards a force field based on density fitting.

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.

Algorithms↗

PubMed2006

Development, validation, and applications of anisotropic polarizable molecular mechanics to study ligand and drug-receptor interactions.

A correct representation of intermolecular interaction energies is necessary for reliable drug-receptor docking studies. While ab initio quantum chemistry with extended basis sets is the most accurate tool for that purpose, its use is precluded for very large molecular complexes. This constitutes the incentive for the development of accurate molecular mechanics potentials, in which the first-order electrostatic, and the second-order polarization energy contributions, are of essential importance. In this paper, we review the most important steps in the development of anisotropic, polarizable molecular mechanics (APMM) procedures. Among these, we illustrate validation tests of the ab initio-grounded, polarizable molecular mechanics potential, SIBFA (Sum of Interactions Between Fragments Ab initio computed). These are done by comparisons with parallel quantum-chemical (QC) results on representative multiply hydrogen-bonded complexes and polycoordinated complexes of one, or of two, divalent metal cations. For both kinds of complexes, the need to reproduce the non-additivity of the QC interaction energies is emphasized. One difficulty arises upon handling flexible molecules, due to the need to account simultaneously and consistently for the onset of inter- and intra-molecular polarization and charge-transfer effects. A new approach in the context of SIBFA was recently developed towards this aim, and tested in two cases of conformation-dependent cation-ligand interaction energies. The first relates to the complexes formed between the mecapto-carboxamide anion, an essential building-block of several Zn-metalloenzyme inhibitors, and Zn(II). The second relates to the complexes of the tetra-anionic pyrophosphate anion, a key building-block of ATP and GTP, with one or two divalent Zn(II) cations used as a probe. In the domain of applications, two recent studies are then presented. The first is the docking of the captopril drug to the active site of the binuclear Zn(II)- beta-lactamase enzyme. The second is the complex of a non-hydrolyzable analog of ATP with the active site of a binuclear Mg(II)-dependent kinase. An extension to an open-shell cation, Cu(II), is finally presented. The encouraging results presented in this review show that APMM procedures could be used in large-scale studies of ligand and drug-receptor interactions.

Angiotensin-Converting Enzyme Inhibitors↗

PubMed2005

Design, synthesis, and biological evaluation of novel naphthoquinone derivatives with CDC25 phosphatase inhibitory activity.

CDC25 dual-specificity phosphatases are essential key regulators of eukaryotic cell cycle progression and the CDC25A and B isoforms are over-expressed in different tumors and related cancer cell lines. CDC25s are now considered to be interesting targets in the search for novel anticancer agents. We describe new compounds derived from vitamin K3 that inhibit CDC25B activity with IC50 values in the low micromolar range. These naphthoquinone derivatives also display antiproliferative activity on HeLa cells as expected for CDC25 inhibitors and inhibit cell growth in a clonogenic assay at submicromolar concentrations. They increase inhibitory tyrosine 15 phosphorylation of CDK and induce the cleavage of PARP, a hallmark of apoptosis.

Antineoplastic Agents↗

PubMed2005

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.

Mathematical Computing↗

PubMed2005

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.

Algorithms↗

PubMed2004

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.

Alanine↗

PubMed2004

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%.

Adenosine Triphosphate↗

PubMed2003

Inclusion of the ligand field contribution in a polarizable molecular mechanics: SIBFA-LF.

To account for the distortion of the coordination sphere that takes place in complexes containing open-shell metal cations such as Cu(II), we implemented, in sum of interactions between fragments ab initio computed (SIBFA) molecular mechanics, an additional contribution to take into account the ligand field splitting of the metal d orbitals. This term, based on the angular overlap model, has been parameterized for Cu(II) coordinated to oxygen and nitrogen ligands. The comparison of the results obtained from density functional theory computations on the one hand and SIBFA or SIBFA-LF on the other shows that SIBFA-LF gives geometric arrangements similar to those obtained from quantum mechanical computations. Moreover, the geometric improvement takes place without downgrading the energetic agreement obtained from SIBFA. The systems considered are Cu(II) interacting with six water molecules, four ammonia or four imidazoles, and four water plus two formate anions.

Journal Article↗

PubMed2002

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.

Algorithms↗