1
|
Moreno Martinez D, Guillaumont D, Guilbaud P. Force Field Parameterization of Actinyl Molecular Cations Using the 12-6-4 Model. J Chem Inf Model 2022; 62:2432-2445. [PMID: 35537184 DOI: 10.1021/acs.jcim.2c00153] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, a set of 12-6-4 force fields (FFs) parameters were developed for the actinyl molecular cations, AnO2n+ (n = 1, 2), from uranium to plutonium for classical molecular dynamics (MD) for four water models: TIP3P, SPC/E, OPC3, and TIP4Pew. Such a non-bonded potential model taking into account the induced dipole between the metallic center and the surrounding molecules has shown better performances for various cations than the classic 12-6 non-bonded potentials. The parametrization method proposed elsewhere for metallic cations has been extended to these molecular cations. In contrast to the actinyl 12-6 FFs from the literature, the new models reproduce correctly both solvation and thermodynamic properties, thanks to the inclusion of the induced dipole term (C4). The transferability of such force fields was assessed by performing MD simulations of carbonato actinyl species, which are highly implicated in actinide migration or actinide extraction from seawater. A highly satisfying agreement was found when comparing the EXAFS signals computed from our MD simulation to the experimental ones. The set of FFs developed here opens new possibilities for the study of actinide chemistry.
Collapse
|
2
|
Macchiagodena M, Pagliai M, Andreini C, Rosato A, Procacci P. Upgrading and Validation of the AMBER Force Field for Histidine and Cysteine Zinc(II)-Binding Residues in Sites with Four Protein Ligands. J Chem Inf Model 2019; 59:3803-3816. [PMID: 31385702 DOI: 10.1021/acs.jcim.9b00407] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
We developed and validated a novel force field in the context of the AMBER parameterization for the simulation of zinc(II)-binding proteins. The proposed force field assumes nonbonded spherical interactions between the central zinc(II) and the coordinating residues. A crucial innovative aspect of our approach is to account for the polarization effects of the cation by redefining the atomic charges of the coordinating residues and an adjustment of Lennard-Jones parameters of Zn-interacting atoms to reproduce mean distance distributions. The optimal transferable parametrization was obtained by performing accurate quantum mechanical calculations on a training set of high-quality protein structures, encompassing the most common folds of zinc(II) sites. The addressed sites contain a zinc(II) ion tetra-coordinated by histidine and cysteine residues and represent about 70% of all physiologically relevant zinc(II) sites in the Protein Data Bank. Molecular dynamics simulations with explicit solvent, carried out on several zinc(II)-binding proteins not included in the training set, show that our model for zinc(II) sites preserves the tetra-coordination of the metal site with remarkable stability, yielding zinc(II)-X mean distances similar to experimental data. Finally, the model was tested by evaluating the zinc(II)-binding affinities, using the alchemical free energy perturbation approach. The calculated dissociation constants correlate satisfactorily with the experimental counterpart demonstrating the validity and transferability of the proposed parameterization for zinc(II)-binding proteins.
Collapse
Affiliation(s)
- Marina Macchiagodena
- Dipartimento di Chimica "Ugo Schiff" , Università degli Studi di Firenze , Via della Lastruccia 3 , 50019 Sesto Fiorentino , Italy
| | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff" , Università degli Studi di Firenze , Via della Lastruccia 3 , 50019 Sesto Fiorentino , Italy
| | - Claudia Andreini
- Dipartimento di Chimica "Ugo Schiff" , Università degli Studi di Firenze , Via della Lastruccia 3 , 50019 Sesto Fiorentino , Italy.,Magnetic Resonance Center (CERM)-Università degli Studi di Firenze , Via L. Sacconi 6 , 50019 Sesto Fiorentino , Italy
| | - Antonio Rosato
- Dipartimento di Chimica "Ugo Schiff" , Università degli Studi di Firenze , Via della Lastruccia 3 , 50019 Sesto Fiorentino , Italy.,Magnetic Resonance Center (CERM)-Università degli Studi di Firenze , Via L. Sacconi 6 , 50019 Sesto Fiorentino , Italy
| | - Piero Procacci
- Dipartimento di Chimica "Ugo Schiff" , Università degli Studi di Firenze , Via della Lastruccia 3 , 50019 Sesto Fiorentino , Italy
| |
Collapse
|
3
|
Nedjoua D, Krallafa AM. Temperature effect on the structure and conformational fluctuations in two zinc knuckles from the mouse mammary tumor virus. Comput Biol Chem 2018; 74:86-93. [PMID: 29567490 DOI: 10.1016/j.compbiolchem.2018.03.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2017] [Revised: 10/28/2017] [Accepted: 03/07/2018] [Indexed: 11/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Drici Nedjoua
- LCPM, Department of Chemistry, University of Oran 1 Ahmed Benbella, PO Box 1524, El m'naouer, Oran, 31000, Algeria.
| | - Abdelghani Mohamed Krallafa
- LCPM, Department of Chemistry, University of Oran 1 Ahmed Benbella, PO Box 1524, El m'naouer, Oran, 31000, Algeria.
| |
Collapse
|
4
|
Abstract
Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
5
|
Falcón‐León MP, Tapia‐Benavides AR, Tlahuext H, Galán‐Vidal C, Suarez‐Castillo OR, Tlahuextl M. The Effect of Zn
II
Coordination on the Addition of 2‐(Aminomethyl)benzimidazole to Acrylonitrile. Eur J Inorg Chem 2014. [DOI: 10.1002/ejic.201402346] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Martha P. Falcón‐León
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Hidalgo, Carr. Pachuca‐Tulancingo km 4.5, Hidalgo, México CP 42184, http://www.uaeh.edu.mx
| | - Antonio R. Tapia‐Benavides
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Hidalgo, Carr. Pachuca‐Tulancingo km 4.5, Hidalgo, México CP 42184, http://www.uaeh.edu.mx
| | - Hugo Tlahuext
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Morelos, Av. Universidad 1001, Morelos, Mexico CP 62209
| | - Carlos Galán‐Vidal
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Hidalgo, Carr. Pachuca‐Tulancingo km 4.5, Hidalgo, México CP 42184, http://www.uaeh.edu.mx
| | - Oscar R. Suarez‐Castillo
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Hidalgo, Carr. Pachuca‐Tulancingo km 4.5, Hidalgo, México CP 42184, http://www.uaeh.edu.mx
| | - Margarita Tlahuextl
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Hidalgo, Carr. Pachuca‐Tulancingo km 4.5, Hidalgo, México CP 42184, http://www.uaeh.edu.mx
| |
Collapse
|
6
|
Lu X, Gaus M, Elstner M, Cui Q. Parametrization of DFTB3/3OB for magnesium and zinc for chemical and biological applications. J Phys Chem B 2014; 119:1062-82. [PMID: 25178644 PMCID: PMC4306495 DOI: 10.1021/jp506557r] [Citation(s) in RCA: 108] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
We report the parametrization of
the approximate density functional
theory, DFTB3, for magnesium and zinc for chemical and biological
applications. The parametrization strategy follows that established
in previous work that parametrized several key main group elements
(O, N, C, H, P, and S). This 3OB set of parameters can thus be used
to study many chemical and biochemical systems. The parameters are
benchmarked using both gas-phase and condensed-phase systems. The
gas-phase results are compared to DFT (mostly B3LYP), ab initio (MP2 and G3B3), and PM6, as well as to a previous DFTB parametrization
(MIO). The results indicate that DFTB3/3OB is particularly successful
at predicting structures, including rather complex dinuclear metalloenzyme
active sites, while being semiquantitative (with a typical mean absolute
deviation (MAD) of ∼3–5 kcal/mol) for energetics. Single-point
calculations with high-level quantum mechanics (QM) methods generally
lead to very satisfying (a typical MAD of ∼1 kcal/mol) energetic
properties. DFTB3/MM simulations for solution and two enzyme systems
also lead to encouraging structural and energetic properties in comparison
to available experimental data. The remaining limitations of DFTB3,
such as the treatment of interaction between metal ions and highly
charged/polarizable ligands, are also discussed.
Collapse
Affiliation(s)
- Xiya Lu
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | | | | | | |
Collapse
|
7
|
Devereux M, Gresh N, Piquemal JP, Meuwly M. A supervised fitting approach to force field parametrization with application to the SIBFA polarizable force field. J Comput Chem 2014; 35:1577-91. [DOI: 10.1002/jcc.23661] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Revised: 04/14/2014] [Accepted: 05/25/2014] [Indexed: 01/07/2023]
Affiliation(s)
- Mike Devereux
- Department of Chemistry; University of Basel; Klingelbergstr 80 CH 4056 Switzerland
| | - Nohad Gresh
- Chemistry & Biology, Nucleo(s)tides & Immunology for Therapy (CBNIT), UMR 8601, CNRS, UFR Biomedicale; Paris France
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, UPMC, Sorbonne Université, Campus de Jussieu; 4 place Jussieu Paris France
| | - Markus Meuwly
- Department of Chemistry; University of Basel; Klingelbergstr 80 CH 4056 Switzerland
| |
Collapse
|
8
|
Chaudret R, Gresh N, Narth C, Lagardère L, Darden TA, Cisneros GA, Piquemal JP. S/G-1: an ab initio force-field blending frozen Hermite Gaussian densities and distributed multipoles. Proof of concept and first applications to metal cations. J Phys Chem A 2014; 118:7598-612. [PMID: 24878003 DOI: 10.1021/jp5051657] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We demonstrate as a proof of principle the capabilities of a novel hybrid MM'/MM polarizable force field to integrate short-range quantum effects in molecular mechanics (MM) through the use of Gaussian electrostatics. This lead to a further gain in accuracy in the representation of the first coordination shell of metal ions. It uses advanced electrostatics and couples two point dipole polarizable force fields, namely, the Gaussian electrostatic model (GEM), a model based on density fitting, which uses fitted electronic densities to evaluate nonbonded interactions, and SIBFA (sum of interactions between fragments ab initio computed), which resorts to distributed multipoles. To understand the benefits of the use of Gaussian electrostatics, we evaluate first the accuracy of GEM, which is a pure density-based Gaussian electrostatics model on a test Ca(II)-H2O complex. GEM is shown to further improve the agreement of MM polarization with ab initio reference results. Indeed, GEM introduces nonclassical effects by modeling the short-range quantum behavior of electric fields and therefore enables a straightforward (and selective) inclusion of the sole overlap-dependent exchange-polarization repulsive contribution by means of a Gaussian damping function acting on the GEM fields. The S/G-1 scheme is then introduced. Upon limiting the use of Gaussian electrostatics to metal centers only, it is shown to be able to capture the dominant quantum effects at play on the metal coordination sphere. S/G-1 is able to accurately reproduce ab initio total interaction energies within closed-shell metal complexes regarding each individual contribution including the separate contributions of induction, polarization, and charge-transfer. Applications of the method are provided for various systems including the HIV-1 NCp7-Zn(II) metalloprotein. S/G-1 is then extended to heavy metal complexes. Tested on Hg(II) water complexes, S/G-1 is shown to accurately model polarization up to quadrupolar response level. This opens up the possibility of embodying explicit scalar relativistic effects in molecular mechanics thanks to the direct transferability of ab initio pseudopotentials. Therefore, incorporating GEM-like electron density for a metal cation enable the introduction of nonambiguous short-range quantum effects within any point-dipole based polarizable force field without the need of an extensive parametrization.
Collapse
Affiliation(s)
- Robin Chaudret
- Sorbonne Universités , UPMC, Univ Paris 06, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005 Paris, France
| | | | | | | | | | | | | |
Collapse
|
9
|
Gkionis K, Kruse H, Platts JA, Mládek A, Koča J, Šponer J. Ion Binding to Quadruplex DNA Stems. Comparison of MM and QM Descriptions Reveals Sizable Polarization Effects Not Included in Contemporary Simulations. J Chem Theory Comput 2014; 10:1326-40. [DOI: 10.1021/ct4009969] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Konstantinos Gkionis
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Holger Kruse
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - James A. Platts
- School
of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| | - Jaroslav Koča
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- CEITEC
- Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| |
Collapse
|
10
|
Li P, Merz KM. Taking into Account the Ion-induced Dipole Interaction in the Nonbonded Model of Ions. J Chem Theory Comput 2014; 10:289-297. [PMID: 24659926 PMCID: PMC3960013 DOI: 10.1021/ct400751u] [Citation(s) in RCA: 255] [Impact Index Per Article: 25.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Metal ions exist in almost half of the proteins in the protein databank and they serve as structural, electron-transfer and catalytic elements in the metabolic processes of organisms. Molecular Dynamics (MD) simulation is a powerful tool that provides information about biomolecular systems at the atomic level. Coupled with the growth in computing power, algorithms like the Particle Mesh Ewald (PME) method have become the accepted standard when dealing with long-range interactions in MD simulations. The nonbonded model of metal ions consists of an electrostatic plus 12-6 Lennard Jones (LJ) potential and is used largely because of its speed relative to more accurate models. In previous work we found that ideal parameters do not exist that reproduce several experimental properties for M(II) ions simultaneously using the nonbonded model coupled with the PME method due to the underestimation of metal ion-ligand interactions. Via a consideration of the nature of the nonbonded model, we proposed that the observed error largely arises from overlooking charge-induced dipole interactions. The electrostatic plus 12-6 LJ potential model works reasonably well for neutral systems but does struggle with more highly charged systems. In the present work we designed and parameterized a new nonbonded model for metal ions by adding a 1/r4 term to the 12-6 model. We call it the 12-6-4 LJ-type nonbonded model due to its mathematical construction. Parameters were determined for 16 +2 metal ions for the TIP3P, SPC/E and TIP4PEW water models. The final parameters reproduce the experimental hydration free energies (HFE), ion-oxygen distances (IOD) in the first solvation shell and coordination numbers (CN) accurately for the metal ions investigated. Preliminary tests on MgCl2 at different concentrations in aqueous solution and Mg2+--nucleic acid systems show reasonable results suggesting that the present parameters can work in mixed systems. The 12-6-4 LJ-type nonbonded model is readily adopted into standard force fields like AMBER, CHARMM and OPLS-AA with only a modest computational overhead. The new nonbonded model doesn't consider charge-transfer effects explicitly and, hence, may not suitable for the simulation of systems where charge-transfer effects play a decisive role.
Collapse
Affiliation(s)
- Pengfei Li
- 2328 New Physics Building, PO Box 118435, University of Florida, Gainesville, Florida 32611-8435, And Department of Chemistry, Michigan State University, East Lansing, Michigan 48824
| | - Kenneth M. Merz
- 2328 New Physics Building, PO Box 118435, University of Florida, Gainesville, Florida 32611-8435, And Department of Chemistry, Michigan State University, East Lansing, Michigan 48824
| |
Collapse
|
11
|
Hofmann FD, Devereux M, Pfaltz A, Meuwly M. Toward force fields for atomistic simulations of iridium-containing complexes. J Comput Chem 2014; 35:18-29. [PMID: 24155105 DOI: 10.1002/jcc.23460] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2013] [Revised: 08/09/2013] [Accepted: 09/08/2013] [Indexed: 12/25/2022]
Abstract
The structural and energetic characterization of metal complexes is important in catalysis and photochemical applications. Unraveling their modes-of-action can be greatly assisted by computation, which typically is restricted to computationally demanding methods including electronic structure calculations with density functional theory. Here, we present an empirical force field based on valence bond theory applicable to a range of octahedral Ir(III) complexes with different coordinating ligands, including iridium complexes with a chiral P,N ligand. Using an approach applicable to metal-containing complexes in general, it is shown that with one common parametrization 85% of the 116 diastereomers--all within 21 kcal/mol of the lowest energy conformation of each series--can be correctly ranked. For neutral complexes, all diastereomers are ranked correctly. This helps to identify the most relevant diastereomers which, if necessary, can be further investigated by more demanding computational methods. Furthermore, if one specific complex is considered, the root mean square deviation between reference data from electronic structure calculations and the force field is ≈1 kcal/mol. This, together with the possibility to carry out explicit simulations in solution paves the way for an atomistic understanding of iridium-containing complexes in catalysis.
Collapse
Affiliation(s)
- Franziska D Hofmann
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056, Basel, Switzerland
| | | | | | | |
Collapse
|
12
|
Zhu T, Xiao X, Ji C, Zhang JZH. A New Quantum Calibrated Force Field for Zinc-Protein Complex. J Chem Theory Comput 2013; 9:1788-98. [PMID: 26587635 DOI: 10.1021/ct301091z] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A quantum calibrated polarizable-charge transfer force field (QPCT) has been proposed to accurately describe the interaction dynamics of zinc-protein complexes. The parameters of the QPCT force field were calibrated by quantum chemistry calculation and capture the polarization and charge transfer effect. QPCTs are validated by molecular dynamic simulation of the hydration shell of the zinc ion, five proteins containing the most common zinc-binding sites (ZnCys2His2, ZnCys3His1, ZnCys4, Zn2Cys6), as well as protein-ligand binding energy in zinc protein MMP3. The calculated results show excellent agreement with the experimental measurement and with results from QM/MM simulation, demonstrating that QPCT is accurate enough to maintain the correct structural integrity of the zinc binding pocket and provide accurate interaction dynamics of the zinc-residue complex. The current approach can also be extended to the study of interaction dynamics of other metal-containing proteins by recalibrating the corresponding parameters to the specific complexes.
Collapse
Affiliation(s)
- Tong Zhu
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China
| | - Xudong Xiao
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China
| | - Changge Ji
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China
| | - John Z H Zhang
- Center for Laser and Computational Biophysics, State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Theoretical and Computational Science, Institutes for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
13
|
Suárez D, Rayón VM, Díaz N, Valdés H. Ab Initio Benchmark Calculations on Ca(II) Complexes and Assessment of Density Functional Theory Methodologies. J Phys Chem A 2011; 115:11331-43. [DOI: 10.1021/jp205101z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- Dimas Suárez
- Departamento de Química Física y Analítica, Facultad de Química Universidad de Oviedo, 33007 Oviedo, Spain
| | - Víctor M. Rayón
- Departamento de Química Física y Química Inorgánica, Facultad de Ciencias Universidad de Valladolid, 47005 Valladolid, Spain
| | - Natalia Díaz
- Departamento de Química Física y Analítica, Facultad de Química Universidad de Oviedo, 33007 Oviedo, Spain
| | - Haydée Valdés
- Departamento de Química Física y Analítica, Facultad de Química Universidad de Oviedo, 33007 Oviedo, Spain
| |
Collapse
|
14
|
Wu R, Lu Z, Cao Z, Zhang Y. A Transferable Non-bonded Pairwise Force Field to Model Zinc Interactions in Metalloproteins. J Chem Theory Comput 2011; 7:433-443. [PMID: 21552372 PMCID: PMC3087386 DOI: 10.1021/ct100525r] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Herein we introduce a novel practical strategy to overcome the well-known challenge of modeling the divalent zinc cation in metalloproteins. The main idea is to design short-long effective functions (SLEF) to describe charge interactions between the zinc ion and all other atoms. This SLEF approach has the following desired features: (1). It is pairwise, additive and compatible with widely used atomic pair-wise force fields for modeling biomolecules; (2). It only changes interactions between the zinc ion and other atoms, and does not affect force field parameters that model other interactions in the system; (3). It is a non-bonded model that is inherently capable to describe different zinc ligands and coordination modes. By optimizing two SLEF parameters as well as zinc vdW parameters through force matching based on Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, we have successfully developed the first SLEF force field (SLEF1) to describe zinc interactions. Extensive molecular dynamics simulations of seven zinc enzyme systems with different coordination ligands and distinct chelation modes (4-,5-,6-fold), including the binuclear zinc active site, yielded zinc coordination numbers and binding distances in good agreement with the corresponding crystal structures as well as ab initio QM/MM MD results. This not only demonstrates the transferability and adequacy of the new SLEF1 force field in describing a variety of zinc proteins, but also indicates that this novel SLEF approach is a promising direction to explore for improving force field description of metal ion interactions.
Collapse
Affiliation(s)
- Ruibo Wu
- Department of Chemistry, New York University, New York, NY 10003 USA
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Zhenyu Lu
- Department of Chemistry, New York University, New York, NY 10003 USA
| | - Zexing Cao
- Department of Chemistry and State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003 USA
| |
Collapse
|
15
|
Keyes T, Napoleon RL. Extending Classical Molecular Theory with Polarization. J Phys Chem B 2010; 115:522-31. [DOI: 10.1021/jp105595q] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Tom Keyes
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Raeanne L. Napoleon
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
16
|
Devereux M, van Severen MC, Parisel O, Piquemal JP, Gresh N. Role of Cation Polarization in holo- and hemi-Directed [Pb(H2O)n]2+ Complexes and Development of a Pb2+ Polarizable Force Field. J Chem Theory Comput 2010; 7:138-47. [DOI: 10.1021/ct1004005] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Mike Devereux
- Université Paris Descartes, Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques, UMR 8601 CNRS, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex06, France; UPMC, Université Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France; and CNRS, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005 Paris, France
| | - Marie-Céline van Severen
- Université Paris Descartes, Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques, UMR 8601 CNRS, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex06, France; UPMC, Université Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France; and CNRS, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005 Paris, France
| | - Olivier Parisel
- Université Paris Descartes, Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques, UMR 8601 CNRS, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex06, France; UPMC, Université Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France; and CNRS, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Université Paris Descartes, Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques, UMR 8601 CNRS, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex06, France; UPMC, Université Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France; and CNRS, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005 Paris, France
| | - Nohad Gresh
- Université Paris Descartes, Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques, UMR 8601 CNRS, UFR Biomédicale, 45 rue des Saints-Pères, 75270 Paris Cedex06, France; UPMC, Université Paris 06, UMR 7616, Laboratoire de Chimie Théorique, Case Courrier 137, 4 Place Jussieu, F-75005 Paris, France; and CNRS, UMR 7616, Laboratoire de Chimie Théorique, case courrier 137, 4 place Jussieu, F-75005 Paris, France
| |
Collapse
|
17
|
Wu JC, Piquemal JP, Chaudret R, Reinhardt P, Ren P. Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field. J Chem Theory Comput 2010; 6:2059-2070. [PMID: 21116445 DOI: 10.1021/ct100091j] [Citation(s) in RCA: 125] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn(2+) are derived from gas-phase ab initio calculation of Zn(2+)-water dimer. The Thole-based dipole polarization is adjusted based on the Constrained Space Orbital Variations (CSOV) calculation while the Symmetry Adapted Perturbation Theory (SAPT) approach is also discussed. The vdW parameters of Zn(2+) have been obtained by comparing the AMOEBA Zn(2+)-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn(2+) solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive Reduced Variational Space (RVS) and Electron Localization Function (ELF) computations in order to unravel the nature of the bonding in Zn(2+)(H(2)O)(n) (n=1,6) complexes and to analyze the charge transfer contribution to the complexes. Results show that the importance of charge transfer decreases as the size of Zn-water cluster grows due to anticooperativity and to changes in the nature of the metal-ligand bonds. Induction could be dominated by polarization when the system approaches condensed-phase and the covelant effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn(2+) in bulk water, one should therefore avoid over-fitting to the ab initio charge transfer energy of Zn(2+)-water dimer. Indeed, in order to avoid overestimation of condensed-phase many-body effects, which is crucial to the transferability of polarizable molecular dynamics, charge transfer should not be included within the classical polarization contribution and should preferably be either incorporated in to the pairwise van der Waals contribution or treated explicitly.
Collapse
Affiliation(s)
- Johnny C Wu
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712-1062, USA
| | | | | | | | | |
Collapse
|
18
|
Gresh N, Audiffren N, Piquemal JP, de Ruyck J, Ledecq M, Wouters J. Analysis of the Interactions Taking Place in the Recognition Site of a Bimetallic Mg(II)−Zn(II) Enzyme, Isopentenyl Diphosphate Isomerase. A Parallel Quantum-Chemical and Polarizable Molecular Mechanics Study. J Phys Chem B 2010; 114:4884-95. [DOI: 10.1021/jp907629k] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Nohad Gresh
- Laboratoire de Pharmacochimie Moléculaire et Cellulaire, U648 INSERM, UFR Biomédicale, Université Paris Descartes, 45 rue des Saints-Pères, 75006 Paris, France, Centre Informatique National de l’Enseignement Supérieur, 950, rue de Saint Priest, 34097 Montpellier, France, Laboratoire de Chimie Théorique, Centre National de la Recherche Scientifique, UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France, and Laboratoire de Chimie Biologique Structurale, FUNDP, 61 Rue de
| | - Nicole Audiffren
- Laboratoire de Pharmacochimie Moléculaire et Cellulaire, U648 INSERM, UFR Biomédicale, Université Paris Descartes, 45 rue des Saints-Pères, 75006 Paris, France, Centre Informatique National de l’Enseignement Supérieur, 950, rue de Saint Priest, 34097 Montpellier, France, Laboratoire de Chimie Théorique, Centre National de la Recherche Scientifique, UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France, and Laboratoire de Chimie Biologique Structurale, FUNDP, 61 Rue de
| | - Jean-Philip Piquemal
- Laboratoire de Pharmacochimie Moléculaire et Cellulaire, U648 INSERM, UFR Biomédicale, Université Paris Descartes, 45 rue des Saints-Pères, 75006 Paris, France, Centre Informatique National de l’Enseignement Supérieur, 950, rue de Saint Priest, 34097 Montpellier, France, Laboratoire de Chimie Théorique, Centre National de la Recherche Scientifique, UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France, and Laboratoire de Chimie Biologique Structurale, FUNDP, 61 Rue de
| | - Jerome de Ruyck
- Laboratoire de Pharmacochimie Moléculaire et Cellulaire, U648 INSERM, UFR Biomédicale, Université Paris Descartes, 45 rue des Saints-Pères, 75006 Paris, France, Centre Informatique National de l’Enseignement Supérieur, 950, rue de Saint Priest, 34097 Montpellier, France, Laboratoire de Chimie Théorique, Centre National de la Recherche Scientifique, UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France, and Laboratoire de Chimie Biologique Structurale, FUNDP, 61 Rue de
| | - Marie Ledecq
- Laboratoire de Pharmacochimie Moléculaire et Cellulaire, U648 INSERM, UFR Biomédicale, Université Paris Descartes, 45 rue des Saints-Pères, 75006 Paris, France, Centre Informatique National de l’Enseignement Supérieur, 950, rue de Saint Priest, 34097 Montpellier, France, Laboratoire de Chimie Théorique, Centre National de la Recherche Scientifique, UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France, and Laboratoire de Chimie Biologique Structurale, FUNDP, 61 Rue de
| | - Johan Wouters
- Laboratoire de Pharmacochimie Moléculaire et Cellulaire, U648 INSERM, UFR Biomédicale, Université Paris Descartes, 45 rue des Saints-Pères, 75006 Paris, France, Centre Informatique National de l’Enseignement Supérieur, 950, rue de Saint Priest, 34097 Montpellier, France, Laboratoire de Chimie Théorique, Centre National de la Recherche Scientifique, UMR 7616, Université Pierre et Marie Curie, 4 place Jussieu, 75252 Paris Cedex 05, France, and Laboratoire de Chimie Biologique Structurale, FUNDP, 61 Rue de
| |
Collapse
|
19
|
De Courcy B, Gresh N, Piquemal JP. Importance of lone pair interactions/redistribution in hard and soft ligands within the active site of alcohol dehydrogenase Zn-metalloenzyme: Insights from electron localization function. Interdiscip Sci 2009; 1:55-60. [DOI: 10.1007/s12539-008-0027-0] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2008] [Revised: 11/18/2008] [Accepted: 12/03/2008] [Indexed: 11/29/2022]
|
20
|
Picot D, Ohanessian G, Frison G. Thermodynamic Stability Versus Kinetic Lability of ZnS4Core. Chem Asian J 2008; 5:1445-54. [DOI: 10.1002/asia.200900624] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|