Biomedical subjects

Jean-Philip Piquemal

Publications and source records attributed to Jean-Philip Piquemal.

11 recordsLinked to original sources

PubMed2006

Generalization of the Gaussian electrostatic model: extension to arbitrary angular momentum, distributed multipoles, and speedup with reciprocal space methods.

The simulation of biological systems by means of current empirical force fields presents shortcomings due to their lack of accuracy, especially in the description of the nonbonded terms. We have previously introduced a force field based on density fitting termed the Gaussian electrostatic model-0 (GEM-0) J.-P. Piquemal et al. [J. Chem. Phys. 124, 104101 (2006)] that improves the description of the nonbonded interactions. GEM-0 relies on density fitting methodology to reproduce each contribution of the constrained space orbital variation (CSOV) energy decomposition scheme, by expanding the electronic density of the molecule in s-type Gaussian functions centered at specific sites. In the present contribution we extend the Coulomb and exchange components of the force field to auxiliary basis sets of arbitrary angular momentum. Since the basis functions with higher angular momentum have directionality, a reference molecular frame (local frame) formalism is employed for the rotation of the fitted expansion coefficients. In all cases the intermolecular interaction energies are calculated by means of Hermite Gaussian functions using the McMurchie-Davidson [J. Comput. Phys. 26, 218 (1978)] recursion to calculate all the required integrals. Furthermore, the use of Hermite Gaussian functions allows a point multipole decomposition determination at each expansion site. Additionally, the issue of computational speed is investigated by reciprocal space based formalisms which include the particle mesh Ewald (PME) and fast Fourier-Poisson (FFP) methods. Frozen-core (Coulomb and exchange-repulsion) intermolecular interaction results for ten stationary points on the water dimer potential-energy surface, as well as a one-dimensional surface scan for the canonical water dimer, formamide, stacked benzene, and benzene water dimers, are presented. All results show reasonable agreement with the corresponding CSOV calculated reference contributions, around 0.1 and 0.15 kcal/mol error for Coulomb and exchange, respectively. Timing results for single Coulomb energy-force calculations for (H(2)O)(n), n=64, 128, 256, 512, and 1024, in periodic boundary conditions with PME and FFP at two different rms force tolerances are also presented. For the small and intermediate auxiliaries, PME shows faster times than FFP at both accuracies and the advantage of PME widens at higher accuracy, while for the largest auxiliary, the opposite occurs.

Models, Chemical↗

PubMed2006

Towards accurate solvation dynamics of divalent cations in water using the polarizable amoeba force field: From energetics to structure.

Molecular dynamics simulations were performed using a modified amoeba force field to determine hydration and dynamical properties of the divalent cations Ca2+ and Mg2+. The extension of amoeba to divalent cations required the introduction of a cation specific parametrization. To accomplish this, the Thole polarization damping model parametrization was modified based on the ab initio polarization energy computed by a constrained space orbital variation energy decomposition scheme. Excellent agreement has been found with condensed phase experimental results using parameters derived from gas phase ab initio calculations. Additionally, we have observed that the coordination of the calcium cation is influenced by the size of the periodic water box, a recurrent issue in first principles molecular dynamics studies.

Amoeba↗

PubMed2006

Quantum mechanics/molecular mechanics electrostatic embedding with continuous and discrete functions.

A quantum mechanics/molecular mechanics (QM/MM) implementation that uses the Gaussian electrostatic model (GEM) as the MM force field is presented. GEM relies on the reproduction of electronic density by using auxiliary basis sets to calculate each component of the intermolecular interaction. This hybrid method has been used, along with a conventional QM/MM (point charges) method, to determine the polarization on the QM subsystem by the MM environment in QM/MM calculations on 10 individual H(2)O dimers and a Mg(2+)-H(2)O dimer. We observe that GEM gives the correct polarization response in cases when the MM fragment has a small charge, while the point charges produce significant over-polarization of the QM subsystem and in several cases present an opposite sign for the polarization contribution. In the case when a large charge is located in the MM subsystem, for example, the Mg(2+) ion, the opposite is observed at small distances. However, this is overcome by the use of a damped Hermite charge, which provides the correct polarization response.

Magnesium↗

PubMed2006

[Pb(H2O)]2+ and [Pb(OH)]+: four-component density functional theory calculations, correlated scalar relativistic constrained-space orbital variation energy decompositions, and topological analysis.

Within the scope of studying the molecular implications of the Pb(2+) cation in environmental and polluting processes, this paper reports Hartree-Fock and density functional theory (B3LYP) four-component relativistic calculations using an all-electron basis set applied to [Pb(H(2)O)](2+) and [Pb(OH)](+), two complexes expected to be found in the terrestrial atmosphere. It is shown that full-relativistic calculations validate the use of scalar relativistic approaches within the framework of density functional theory. [Pb(H(2)O)](2+) is found C(2v) at any level of calculations whereas [Pb(OH)](+) can be found bent or linear depending of the computational methodology used. When C(s) is found the barrier to inversion through the C(infinityv) structure is very low, and can be overcome at high enough temperature, making the molecule floppy. In order to get a better understanding of the bonding occurring between the Pb(2+) cation and the H(2)O and OH(-) ligands, natural bond orbital and atoms-in-molecule calculations have been performed. These approaches are supplemented by a topological analysis of the electron localization function. Finally, the description of these complexes is refined using constrained-space orbital variation complexation energy decompositions.

Cations, Divalent↗

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

Revisiting the geometry of nd10 (n+1)s0 [M(H2O)]p+ complexes using four-component relativistic DFT calculations and scalar relativistic correlated CSOV energy decompositions (M(p+) = Cu+, Zn2+, Ag+, Cd2+, Au+, Hg2+).

Hartree-Fock and DFT (B3LYP) nonrelativistic (scalar relativistic pseudopotentials for the metallic cation) and relativistic (molecular four-component approach coupled to an all-electron basis set) calculations are performed on a series of six nd10 (n+1)s0 [M(H2O)]p+ complexes to investigate their geometry, either planar C2v or nonplanar C(s). These complexes are, formally, entities originating from the complexation of a water molecule to a metallic cation: in the present study, no internal reorganization has been found, which ensures that the complexes can be regarded as a water molecule interacting with a metallic cation. For [Au(H2O)]+ and [Hg(H2O)]2+, it is observed that both electronic correlation and relativistic effects are required to recover the C(s) structures predicted by the four-component relativistic all-electron DFT calculations. However, including the zero-point energy corrections makes these shallow C(s) minima vanish and the systems become floppy. In all other systems, namely [Cu(H2O)]+, [Zn(H2O)]2+, [Ag(H2O)]+, and [Cd(H2O)]2+, all calculations predict a C2v geometry arising from especially flat potential energy surfaces related to the out-of-plane wagging vibration mode. In all cases, our computations point to the quasi-perfect transferability of the atomic pseudopotentials considered toward the molecular species investigated. A rationalization of the shape of the wagging potential energy surfaces (i.e., single well vs. double well) is proposed based on the Constrained Space Orbital Variation decompositions of the complexation energies. Any way of stabilizing the lowest unoccupied orbital of the metallic cation is expected to favor charge-transfer (from the highest occupied orbital(s) of the water ligand), covalence, and, consequently, C(s) structures. The CSOV complexation energy decompositions unambiguously reveal that such stabilizations are achieved by means of relativistic effects for [Au(H2O)]+, and, to a lesser extent, for [Hg(H2O)]2+. Such analyses allow to numerically quantify the rule of thumb known for Au+ which, once again, appears as a better archetype of a relativistic cation than Hg2+. This observation is reinforced due to the especially high contribution of the nonadditive correlation/relativity terms to the total complexation energy of [Au(H2O)]+.

Cadmium↗

PubMed2005

A CSOV study of the difference between HF and DFT intermolecular interaction energy values: the importance of the charge transfer contribution.

Intermolecular interaction energy decompositions using the Constrained Space Orbital Variation (CSOV) method are carried out at the Hartree-Fock level on the one hand and using DFT with usual GGA functionals on the other for a number of model complexes to analyze the role of electron correlation in the intermolecular stabilization energy. In addition to the overall stabilization, the results provide information on the variation, with respect to the computational level, of the different contributions to the interaction energy. The complexes studied are the water linear dimer, the N-methylformamide dimer, the nucleic acid base pairs, the benzene-methane and benzene-N2 van der Waals complexes, [Cu+ -(ImH)3]2, where "ImH" stands for the Imidazole ligand, and ImH-Zn++. The variation of the frozen core energy (the sum of the intermolecular electrostatic energy and the Pauli repulsion energy) calculated from the unperturbed orbitals of the interacting entities indicates that the intramolecular correlation contributions can be stabilizing as well as destabilizing, and that general trends can be derived from the results obtained using usual density functionals. The most important difference between the values obtained from HF and DFT computations concerns the charge transfer contribution, which, in most cases, undergoes the largest increase. The physical meaning of these results is discussed. The present work gives reference calculations that might be used to parametrize new correlated molecular mechanics potentials.

Algorithms↗

PubMed2005

Intermolecular electrostatic energies using density fitting.

A method is presented to calculate the electron-electron and nuclear-electron intermolecular Coulomb interaction energy between two molecules by separately fitting the unperturbed molecular electron density of each monomer. This method is based on the variational Coulomb fitting method which relies on the expansion of the ab initio molecular electron density in site-centered auxiliary basis sets. By expanding the electron density of each monomer in this way the integral expressions for the intermolecular electrostatic calculations are simplified, lowering the operation count as well as the memory usage. Furthermore, this method allows the calculation of intermolecular Coulomb interactions with any level of theory from which a one-electron density matrix can be obtained. Our implementation is initially tested by calculating molecular properties with the density fitting method using three different auxiliary basis sets and comparing them to results obtained from ab initio calculations. These properties include dipoles for a series of molecules, as well as the molecular electrostatic potential and electric field for water. Subsequently, the intermolecular electrostatic energy is tested by calculating ten stationary points on the water dimer potential-energy surface. Results are presented for electron densities obtained at four different levels of theory using two different basis sets, fitted with three auxiliary basis sets. Additionally, a one-dimensional electrostatic energy surface scan is performed for four different systems (H2O dimer, Mg2+-H2O, Cu+-H2O, and n-methyl-formamide dimer). Our results show a very good agreement with ab initio calculations for all properties as well as interaction energies.

Computer Simulation↗

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↗

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↗