1
|
Huang H, Zhao DX, Zhao J, Chen X, Liu C, Yang ZZ. Origin of Enantioselectivity in Engineered Cytochrome c-Catalyzed Carbon-Radical FePP Hydrolysis Revealed Using QM/MM (ABEEM Polarizable Force Field) and MD Simulations. J Phys Chem B 2024; 128:3807-3823. [PMID: 38605466 DOI: 10.1021/acs.jpcb.3c07158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2024]
Abstract
The origin of highly efficient asymmetric aminohydroxylation of styrene catalyzed by engineered cytochrome c is investigated by the developed Atom-Bond Electronegativity Equalization Method polarizable force field (ABEEM PFF), which is a combined outcome of electronic and steric effects. Model molecules were used to establish the charge parameters of the ABEEM PFF, for which the bond-stretching and angle-bending parameters were obtained by using a combination of modified Seminario and scan methods. The interactions between carbon-radical Fe-porphyrin (FePP) and waters are simulated by molecular dynamics, which shows a clear preference for the pre-R over the pre-S. This preference is attributed to the hydrogen-bond between the mutated 100S and 101P residues as well as van der Waals interactions, enforcing a specific conformation of the carbon-radical FePP complex within the binding pocket. Meanwhile, the hydrogen-bond between water and the nitrogen atom in the active intermediate dictates the stereochemical outcome. Quantum mechanics/molecular mechanics (QM/MM (ABEEM PFF)) and free-energy perturbation calculations elucidate that the 3RTS is characterized by sandwich-like structure among adjacent amino acid residues, which exhibits greater stability than crowed arrangement in 3STS and enables the R enantiomer to form more favorably. Thus, this study provides mechanistic insight into the catalytic reaction of hemoproteins.
Collapse
Affiliation(s)
- Hong Huang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| | - Jian Zhao
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| | - Xin Chen
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| |
Collapse
|
2
|
Wang Y, Liu L, Gao Y, Zhao J, Liu C, Gong L, Yang Z. Development of a QM/MM(ABEEM) method for the deprotonation of neutral and cation radicals in the G-tetrad and GGX(8-oxo-G) tetrad. Phys Chem Chem Phys 2023; 26:504-516. [PMID: 38084041 DOI: 10.1039/d3cp04357f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
The rapid deprotonation of G˙+ in the DNA strand impedes positive charge (hole) transfer, whereas the slow deprotonation rate of G˙+ in the G-tetrad makes it a more suitable carrier for hole conduction. The QM/MM(ABEEM) combined method, which involves the integration of QM and the ABEEM polarizable force field (ABEEM PFF), was developed to investigate the deprotonation of neutral and cation free radicals in the G-tetrad and GGX(8-oxo-G) tetrad (xanthine and 8-oxoguanine dual substituted G-tetrad). By incorporating valence-state electronegativity piecewise functions χ*(r) and implementing charge local conservation conditions, QM/MM(ABEEM) possesses the advantage of accurately simulating charge transfer and polarization effect during deprotonation. The activation energy calculated by the QM method of X˙ is the lowest among other bases in the GGX(8-oxo-G) tetrad, which is supported by the computation of the average electronegativity calculated by ABEEM PFF. By utilizing QM/MM(ABEEM) with a two-way free energy perturbation method, the deprotonation activation energy of X˙ in the GGX(8-oxo-G) tetrad is determined to be 33.0 ± 2.1 kJ mol-1, while that of G˙+ in the G-tetrad is 20.7 ± 0.6 kJ mol-1, consistent with the experimental measurement of 20 ± 1.0 kJ mol-1. These results manifest that X˙ in the GGX(8-oxo-G) tetrad exhibits a slower deprotonation rate than G˙+ in the G-tetrad, suggesting that the GGX(8-oxo-G) tetrad may serve as a more favorable hole transport carrier. Furthermore, the unequal average electronegativities of bases in the GGX(8-oxo-G) tetrad impede the deprotonation rate. This study provides a potential foundation for investigating the microscopic mechanism of DNA electronic devices.
Collapse
Affiliation(s)
- Yue Wang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Linlin Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Yue Gao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Jiayue Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Lidong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Zhongzhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| |
Collapse
|
3
|
Zhao J, Zhu ZW, Zhao DX, Yang ZZ. Atomic charges in molecules defined by molecular real space partition into atomic subspaces. Phys Chem Chem Phys 2023; 25:9020-9030. [PMID: 36928882 DOI: 10.1039/d2cp05428k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/10/2023]
Abstract
Atomic charge (AC), which is the charge distribution of a molecule, is an important property that is closely associated with structures, reactivities, and intra- and inter-molecular interactions among molecules. Several theoretical models or methods can be used to obtain the magnitudes of AC with different characteristics. These models can be classified into fuzzy-atoms models and models partitioning a molecule into individual atoms with sharp boundaries. The first category includes Mulliken, natural population analysis (NPA), Hirshfeld, Merz-Kollman-Singh (MK), CHELPG, the electronegativity equalization method (EEM), the atom-bond electronegativity equalization method (ABEEM), and atomic polar tensor (APT). The second category is derived from quantum chemical topology (QCT) and includes the quantum theory of atoms in molecules (QTAIM) and QCT analysis based on the potential acting on one electron in a molecule (PAEMQCT). Herein, after giving a bird's-eye view of the population methods of the first category, we specifically describe some features of the second category. We only present the basic framework of QCT for obtaining ACs from QTAIM and PAEMQCT and show their important characteristics. QCT establishes the basis of the following chemical concept: a molecule is spatially partitioned into individual atoms with sharp boundaries. The ACs from QTAIM are close to the atomic valence in chemistry, and ACs from PAEMQCT may be practically suitable for modeling intra- and inter-molecular interactions.
Collapse
Affiliation(s)
- Jian Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, Liaoning province, 116029, China. .,State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, Liaoning province, 116023, China
| | - Zun-Wei Zhu
- School of Materials Science and Engineering, Anyang Institute of Technology, Anyang, Henan province, 455000, China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, Liaoning province, 116029, China.
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, Liaoning province, 116029, China.
| |
Collapse
|
4
|
Liu C, Jiang H, Li Y, Xue B, Yao YY, Yang ZZ. Development of a QM/MM(ABEEM) method combined with a polarizable force field to investigate the excision reaction mechanism of damaged thymine. Phys Chem Chem Phys 2023; 25:3432-3448. [PMID: 36637033 DOI: 10.1039/d2cp05873a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
This paper focuses on the development of a quantum mechanics/molecular mechanics method using the ABEEM polarizable force field (QM/MM(ABEEM) method) to investigate the excision reaction mechanism of damaged thymine. This method does not simply combine the QM method with the polarizable force field. A valence electronegativity piecewise function with the distance between atoms as a variable is introduced to describe the atomic partial charges, and changes greatly during the reaction process. At the same time, the charge transfer effect is treated using the condition of local charge conservation. Compared with the traditional QM/MM method, the QM/MM(ABEEM) method can more accurately simulate the polarization effect and charge transfer effect in the reaction process. Focusing on the controversial problems of the excision of damaged bases, six reaction pathways were designed for monofunctional and difunctional deglycosylation of neutral bases and protonated bases. The results show that the QM/MM(ABEEM) method accurately simulates the polarization effect, charge transfer effect, activation energy and other properties of the reaction process. The process in which the active residue Asp activates the nucleophile H2O to attack the protonated base is the preferred path. The average activation energy and free activation energy of the protonated base are 7.00-14.00 kcal mol-1 lower than that of the neutral base. The study in this paper is helpful to understand the mechanism of repair enzymes in repairing bases.
Collapse
Affiliation(s)
- Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - He Jiang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Yue Li
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Bing Xue
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Yu-Ying Yao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.
| |
Collapse
|
5
|
Lu LN, Liu C, Yang ZZ, Zhao DX. Refined models of coordination between Al3+/Mg2+ and enzyme in molecular dynamics simulation in terms of ABEEM polarizable force field. J Mol Graph Model 2022; 114:108190. [DOI: 10.1016/j.jmgm.2022.108190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 04/04/2022] [Accepted: 04/05/2022] [Indexed: 10/18/2022]
|
6
|
Ren WS, Rahman A, Jiang KB, Deng H, Zhao YY, Zhang WJ, Liu K, Qian P, Guo H. Unraveling the Origins of Changing Product Specificity Properties of Arginine Methyltransferase PRMT7 by the E181D and E181D/Q329A Mutations through QM/MM MD and Free-Energy Simulations. J Chem Theory Comput 2022; 18:2631-2641. [PMID: 35316052 DOI: 10.1021/acs.jctc.1c01219] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Arginine methylations can regulate important biological processes and affect many cellular activities, and the enzymes that catalyze the methylations are protein arginine methyltransferases (PRMTs). The biological consequences of arginine methylations depend on the methylation states of arginine that are determined by the PRMT's product specificity. Nevertheless, it is still unclear how different PRMTs may generate different methylation states for the target proteins. PRMT7 is the only known member of type III PRMT that produces monomethyl arginine (MMA) product. Interestingly, its E181D and E181D/Q329A mutants can catalyze, respectively, the formation of asymmetrically dimethylated arginine (ADMA) and symmetrically dimethylated arginine (SDMA). The reasons as to why the mutants have the abilities to add the second methyl group and E181D (E181D/Q329A) has the unique product specificity in generating ADMA (SDMA) have not been understood. Here, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) and potential of mean force (PMF) free-energy simulations are performed for the E181D and E181D/Q329A mutants to understand the origin for their ability to generate, respectively, ADMA and SDMA. The simulations show that the free-energy barrier for adding the second methyl group to MMA in E181D (E181D/Q329A) to produce ADMA (SDMA) is considerably lower than the corresponding barriers in wild type and E181D/Q329A (wild type and E181D), consistent with experimental observations. Some important factors that contribute to the change of the activity and product specificity due to the E181D and E181D/Q329A mutations are identified based on the data from the simulations and analysis. It is shown that the transferable methyl group (from SAM) and Nη2 (the nitrogen atom that is methylated in the substrate MMA) can only form good near-attack conformations in the E181D reaction state for the methyl transfer (not in wild type and E181D/Q329A), while the transferable methyl group and Nη1 (the nitrogen atom that is not methylated in the substrate MMA) can only form good near-attack conformations in E181D/Q329A (not in wild type and E181D). The results suggest that the steric repulsions in the reaction state between the methyl group on MMA and active-site residues (e.g., Q329) and the release of such repulsions (e.g., from the Q329A mutation) may play an important role in generating specific near-attack conformations for the methyl transfer and controlling the product specificity for the mutants. The general principle identified in this work for PRMT7 is expected to be useful for understanding the activity and product specificity of other PRMTs as well.
Collapse
Affiliation(s)
- Wan-Sheng Ren
- Chemistry and Material Science Faculty, Shandong Agricultural University, Taian 271018, P. R. China
| | - Adua Rahman
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Kai-Bin Jiang
- Chemistry and Material Science Faculty, Shandong Agricultural University, Taian 271018, P. R. China
| | - Hao Deng
- Chemistry and Material Science Faculty, Shandong Agricultural University, Taian 271018, P. R. China
| | - Yuan-Yuan Zhao
- Chemistry and Material Science Faculty, Shandong Agricultural University, Taian 271018, P. R. China
| | - Wei-Jie Zhang
- Chemistry and Material Science Faculty, Shandong Agricultural University, Taian 271018, P. R. China
| | - Kedian Liu
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ping Qian
- Chemistry and Material Science Faculty, Shandong Agricultural University, Taian 271018, P. R. China
| | - Hong Guo
- Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
7
|
Zhang C, Zhao DX, Feng Y, Wang J, Yang ZZ. Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution. Phys Chem Chem Phys 2022; 24:4232-4250. [DOI: 10.1039/d1cp05676j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...
Collapse
|
8
|
P. Oliveira M, Hünenberger PH. Systematic optimization of a fragment-based force field against experimental pure-liquid properties considering large compound families: application to oxygen and nitrogen compounds. Phys Chem Chem Phys 2021; 23:17774-17793. [PMID: 34350931 PMCID: PMC8386690 DOI: 10.1039/d1cp02001c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 06/30/2021] [Indexed: 12/04/2022]
Abstract
The CombiFF approach is a workflow for the automated refinement of force-field parameters against experimental condensed-phase data, considering entire classes of organic molecules constructed using a fragment library via combinatorial isomer enumeration. One peculiarity of this approach is that it relies on an electronegativity-equalization scheme to account for induction effects within molecules, with values of the atomic hardness and electronegativity as electrostatic parameters, rather than the partial charges themselves. In a previous article [M. P. Oliveira, M. Andrey, S. R. Rieder, L. Kern, D. F. Hahn, S. Riniker, B. A. C. Horta and P. H. Hünenberger, J. Chem. Theory. Comput. 2020, 16, 7525], CombiFF was introduced and applied to calibrate a GROMOS-compatible united-atom force field for the saturated acyclic (halo-)alkane family. Here, this scheme is employed for the construction of a corresponding force field for saturated acyclic compounds encompassing eight common chemical functional groups involving oxygen and/or nitrogen atoms, namely: ether, aldehyde, ketone, ester, alcohol, carboxylic acid, amine, and amide. Monofunctional as well as homo-polyfunctional compounds are considered. A total of 1712 experimental liquid densities ρliq and vaporization enthalpies ΔHvap concerning 1175 molecules are used for the calibration (339 molecules) and validation (836 molecules) of the 102 non-bonded interaction parameters of the force field. Using initial parameter values based on the GROMOS 2016H66 parameter set, convergence is reached after five iterations. Given access to one processor per simulated system, this operation only requires a few days of wall-clock computing time. After optimization, the root-mean-square deviations from experiment are 29.9 (22.4) kg m-3 for ρliq and 4.1 (5.5) kJ mol-1 for ΔHvap for the calibration (validation) set. Thus, a very good level of agreement with experiment is achieved in terms of these two properties, although the errors are inhomogeneously distributed across the different chemical functional groups.
Collapse
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 5503
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 5503
| |
Collapse
|
9
|
Liu C, Lv C, Yao YY, Du X, Zhao DX, Yang ZZ. Water-Mediated Oxidation of Guanine by a Repair Enzyme: Simulation Using the ABEEM Polarizable Force Field. J Chem Theory Comput 2021; 17:3525-3538. [PMID: 34018392 DOI: 10.1021/acs.jctc.1c00107] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The recognition mechanism of oxidative damage in organisms has long been a research hotspot. Water is an important medium in the recognition process, but its specific role remains unknown. There is a need to develop a suitable force field that can adequately describe the electrostatic, hydrogen bond, and other interactions among the molecules in the complex system of the repair enzyme and oxidized base. The developing ABEEM polarizable force field (PFF) has been used to simulate the repaired enzyme hOGG1 and oxidized DNA (PDB ID: 1EBM) in a biological environment, and the corresponding results are better than those of the fixed-charge force fields OPLS/AA and AMBER OL15. 8-Oxo-G is recognized by Gln315 of hOGG1 mainly through hydrogen bonds mediated by continuous exchange of 2 water molecules. Phe319 and Cys253 are stacked on both sides of the π planes of bases to form sandwich structures. The charge polarization effect gives an important signal to drive the exchange of water molecules and maintains the recognition of oxidation bases by enzymes. The mediated main water molecule A and mediated auxiliary water molecule B together pull Gln315 to recognize 8-oxo-G by hydrogen bond interactions. Then, the charge polarization signal of solvent water molecule C with a large absolute charge causes the absolute charge of O atoms in water molecule A or B to increase by approximately 0.2 e, and water molecule A or B leaves Gln315 and 8-oxo-G. The other water molecule and water molecule C synergistically recognize 8-oxo-G with Gln315. Even though the water molecules between Gln315 and 8-oxo-G are removed, the MD simulation results show that water molecules appear between Gln315 and 8-oxo-G in a very short time (<2 ps). The dwell time of each water molecule is approximately 60 ps. The radial distribution function and dwell time support the correctness of the above mechanism. These polarization effects and hydrogen bonding interactions cannot be simulated by a fixed-charge force field.
Collapse
Affiliation(s)
- Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Change Lv
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Yu-Ying Yao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Xue Du
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| |
Collapse
|
10
|
Portelinha J, Duay SS, Yu SI, Heilemann K, Libardo MDJ, Juliano SA, Klassen JL, Angeles-Boza AM. Antimicrobial Peptides and Copper(II) Ions: Novel Therapeutic Opportunities. Chem Rev 2021; 121:2648-2712. [PMID: 33524257 DOI: 10.1021/acs.chemrev.0c00921] [Citation(s) in RCA: 48] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The emergence of new pathogens and multidrug resistant bacteria is an important public health issue that requires the development of novel classes of antibiotics. Antimicrobial peptides (AMPs) are a promising platform with great potential for the identification of new lead compounds that can combat the aforementioned pathogens due to their broad-spectrum antimicrobial activity and relatively low rate of resistance emergence. AMPs of multicellular organisms made their debut four decades ago thanks to ingenious researchers who asked simple questions about the resistance to bacterial infections of insects. Questions such as "Do fruit flies ever get sick?", combined with pioneering studies, have led to an understanding of AMPs as universal weapons of the immune system. This review focuses on a subclass of AMPs that feature a metal binding motif known as the amino terminal copper and nickel (ATCUN) motif. One of the metal-based strategies of hosts facing a pathogen, it includes wielding the inherent toxicity of copper and deliberately trafficking this metal ion into sites of infection. The sudden increase in the concentration of copper ions in the presence of ATCUN-containing AMPs (ATCUN-AMPs) likely results in a synergistic interaction. Herein, we examine common structural features in ATCUN-AMPs that exist across species, and we highlight unique features that deserve additional attention. We also present the current state of knowledge about the molecular mechanisms behind their antimicrobial activity and the methods available to study this promising class of AMPs.
Collapse
Affiliation(s)
- Jasmin Portelinha
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States
| | - Searle S Duay
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States.,Chemistry Department, Adamson University, 900 San Marcelino Street, Ermita, Manila 1000, Philippines
| | - Seung I Yu
- Department of Molecular and Cell Biology, University of Connecticut, 91 North Eagleville Road, Storrs, Connecticut 06269, United States
| | - Kara Heilemann
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States
| | - M Daben J Libardo
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States
| | - Samuel A Juliano
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States
| | - Jonathan L Klassen
- Department of Molecular and Cell Biology, University of Connecticut, 91 North Eagleville Road, Storrs, Connecticut 06269, United States
| | - Alfredo M Angeles-Boza
- Department of Chemistry, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States.,Institute of Material Science, University of Connecticut, 55 North Eagleville Road, Storrs, Connecticut 06269, United States
| |
Collapse
|
11
|
Oliveira MP, Andrey M, Rieder SR, Kern L, Hahn DF, Riniker S, Horta BAC, Hünenberger PH. Systematic Optimization of a Fragment-Based Force Field against Experimental Pure-Liquid Properties Considering Large Compound Families: Application to Saturated Haloalkanes. J Chem Theory Comput 2020; 16:7525-7555. [DOI: 10.1021/acs.jctc.0c00683] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Maurice Andrey
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Leyla Kern
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - David F. Hahn
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
12
|
Lu LN, Liu C, Yang ZZ. Systematic Parameterization and Simulation of Boronic Acid-β-Lactamase Aqueous Solution in Developing the ABEEMσπ Polarizable Force Field. J Phys Chem A 2020; 124:8614-8632. [PMID: 32910648 DOI: 10.1021/acs.jpca.0c06806] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Boronic acid, an inhibitor of β-lactamase, has begun to be applied to the treatment of biological infections and tumors. Scientists are working to develop new and more effective boronic acid. Molecular dynamics (MD) simulation provides a powerful auxiliary tool for drug design. However, the current force fields have no boron-related parameters. In this work, an atom-bond electronegativity equalization method at the σπ level (ABEEMσπ) polarizable force field (ABEEMσπ PFF) of boronic acid and β-lactamase has been developed to determine the potential functions and parameters. The interaction between boron and serine in β-lactamase is regarded as a bonded mode. The interaction between them is simulated by the Morse potential energy function, which is close to the experimental change of the stretching potential energy in a large range. The potential energy surfaces of the bond length, bond angle, and dihedral angle of boronic acid-β-lactamase have the same stability point and change trend as M06-2X/6-311G**. For 47 boronic acid-β-lactamase training molecules, the linear correlation coefficient (R) of the charge distribution between the ABEEMσπ PFF and HF/STO-3G is greater than 0.96. Attributed to the fact that the charge distribution of the ABEEMσπ PFF can fluctuate with the change of geometry and environment, the polarization effect and charge-transfer effect are well reflected. The binding ability of different boronic acids with the same β-lactamase is different. A total of 10 boronic acid-β-lactamase model molecules and 10 boronic acid-β-lactamase and water complexes are simulated. The order of binding energy of five large model molecules calculated by the ABEEMσπ PFF is consistent with that of the MP2 method. The binding energies of boronic acid-β-lactamase and water complexes are close to those of the MP2 method. The results of MD simulation of five aqueous boronic acid-β-lactamase complexes in the NVT ensemble verify the rationality of boron-related parameters of the ABEEMσπ PFF, which have a good application prospect. This study lays a solid theoretical foundation for further study of the inhibition of boronic acid on β-lactamase.
Collapse
Affiliation(s)
- Li-Nan Lu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| |
Collapse
|
13
|
Amin KS, Hu X, Salahub DR, Baldauf C, Lim C, Noskov S. Benchmarking polarizable and non-polarizable force fields for Ca2+–peptides against a comprehensive QM dataset. J Chem Phys 2020; 153:144102. [DOI: 10.1063/5.0020768] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Kazi S. Amin
- CMS – Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Xiaojuan Hu
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Dennis R. Salahub
- Department of Chemistry, CMS – Centre for Molecular Simulation, IQST – Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Carmay Lim
- Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan
- Department of Chemistry, National Tsing Hua University, Hsinchu 300, Taiwan
| | - Sergei Noskov
- CMS – Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
14
|
Shi H, Gong LD, Liu C, Lu LN, Yang ZZ. ABEEM/MM OH - Models for OH -(H 2O) n Clusters and Aqueous OH -: Structures, Charge Distributions, and Binding Energies. J Phys Chem A 2020; 124:5963-5978. [PMID: 32520555 DOI: 10.1021/acs.jpca.0c03941] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Based on the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), two fluctuating charge models of OH--water system were proposed. The difference between these two models is whether there is charge transfer between OH- and its first-shell water molecules. The structures, charge distributions, charge transfer, and binding energies of the OH-(H2O)n (n = 1-8, 10, 15, 23) clusters were studied by these two ABEEM/MM models, the OPLS/AA force field, the OPLS-SMOOTH/AA force field, and the QM methods. The results demonstrate that two ABEEM/MM models can search out all stable structures just as the QM methods, and the structures and charge distributions agree well with those from the QM calculations. The structures, the charge transfer, and the strength of hydrogen bonds in the first hydration shell are closely related to the coordination number of OH-. Molecular dynamics simulations on the aqueous OH- solution are performed at 298 and 278 K using ABEEM/MM-I model. The MD results show that the populations of three-, four-, and five-coordinated OH- are 29.6%, 67.1%, and 3.4% at 298 K, respectively, and those of two-, three-, four-, and five-coordinated OH- are 10.8%, 44.9%, 39.2%, and 4.9% at 278 K, respectively; the average hydrogen bond lengths and the hydrogen bond angle in the first shell increase with the temperature decreasing.
Collapse
Affiliation(s)
- Hua Shi
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.,School of Marine Science and Environment, Dalian Ocean University, Dalian 116023, People's Republic of China
| | - Li-Dong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Li-Nan Lu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| |
Collapse
|
15
|
Song LF, Sengupta A, Merz KM. Thermodynamics of Transition Metal Ion Binding to Proteins. J Am Chem Soc 2020; 142:6365-6374. [PMID: 32141296 DOI: 10.1021/jacs.0c01329] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Modeling the thermodynamics of a transition metal (TM) ion assembly be it in proteins or in coordination complexes affords us a better understanding of the assembly and function of metalloclusters in diverse application areas including metal organic framework design, TM-based catalyst design, the trafficking of TM ions in biological systems, and drug design in metalloprotein platforms. While the structural details of TM ions bound to metalloproteins are generally well understood via experimental and computational approaches, accurate studies describing the thermodynamics of TM ion binding are rare. Herein, we demonstrate that we can obtain accurate structural and absolute binding free energies of Co2+ and Ni2+ to the enzyme glyoxalase I using an optimized 12-6-4 (m12-6-4) potential. Critically, this model simultaneously reproduces the solvation free energy of the individual TM ions and reproduces the thermodynamics of TM ion-ligand coordination as well as the thermodynamics of TM ion binding to a protein active site unlike extant models. We find the incorporation of the thermodynamics associated with protonation state changes for the TM ion (un)binding to be crucial. The high accuracy of m12-6-4 potential in this study presents an accurate route to explore more complicated processes associated with TM cluster assembly and TM ion transport.
Collapse
|
16
|
Li H, Wang D, Zhao X, Lu LN, Liu C, Gong LD, Zhao DX, Yang ZZ. Reaction mechanism of NO with hydrolysates of NAMI-A: an MD simulation by combining the QM/MM(ABEEM) with the MD-FEP method. J Comput Chem 2019; 40:1141-1150. [PMID: 30375671 DOI: 10.1002/jcc.25734] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2018] [Revised: 09/27/2018] [Accepted: 09/28/2018] [Indexed: 11/06/2022]
Abstract
Nitrosylation reaction mechanisms of the hydrolysates of NAMI-A and hydrolysis reactions of ruthenium nitrosyl complexes were investigated in the triplet state and the singlet state. Activation free energies were calculated by combining the QM/MM(ABEEM) method with free energy perturbation theory, and the explicit solvent environment was simulated by an ABEEMσπ polarizable force field. Our results demonstrate that nitrosylation reactions of the hydrolysates of NAMI-A occur in both the triplet and the singlet states. The Ru-N-O angle of the triplet ruthenium nitrosyl complexes is in the range of 132.0°-138.2°. However, all the ruthenium nitrosyl complexes at the singlet state show an almost linear Ru-N-O angle. The nitrosylation reaction happens prior to the hydrolysis reaction for the first-step hydrolysates. The activation free energies of the nitrosylation reactions show that the H2 O-NO exchange reaction of [RuCl4 (Im)(H2 O)] in the singlet spin sate is the most likely one. Comparing with the activation free energies of the hydrolysis reactions of the ruthenium nitrosyl complexes, the results indicate that the rate of the DMSO-H2 O exchange reaction of [RuCl3 (NO)(Im)(DMSO)] is faster than that of [RuCl3 (H2 O)(Im)(DMSO)] in both the triplet spin state and the singlet spin state. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Hui Li
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China.,Department of Chemistry, Bohai University, Jinzhou 121013, China
| | - Di Wang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Xin Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Li-Nan Lu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Li-Dong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| |
Collapse
|
17
|
Jing Z, Liu C, Cheng SY, Qi R, Walker BD, Piquemal JP, Ren P. Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications. Annu Rev Biophys 2019; 48:371-394. [PMID: 30916997 DOI: 10.1146/annurev-biophys-070317-033349] [Citation(s) in RCA: 223] [Impact Index Per Article: 44.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Sara Y Cheng
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Brandon D Walker
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA; .,Sorbonne Université, CNRS, Laboratoire de Chimie Theórique, 75252 Paris CEDEX 05, France.,Institut Universitaire de France, 75005 Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| |
Collapse
|
18
|
Chu H, Cao L, Peng X, Li G. Polarizable force field development for lipids and their efficient applications in membrane proteins. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1312] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Huiying Chu
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| | - Liaoran Cao
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| | - Xiangda Peng
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| |
Collapse
|
19
|
Liu C, Li Y, Han BY, Gong LD, Lu LN, Yang ZZ, Zhao DX. Development of the ABEEMσπ Polarization Force Field for Base Pairs with Amino Acid Residue Complexes. J Chem Theory Comput 2017; 13:2098-2111. [DOI: 10.1021/acs.jctc.6b01206] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Cui Liu
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Yue Li
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Bing-Yu Han
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Li-Dong Gong
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Li-Nan Lu
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Zhong-Zhi Yang
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Dong-Xia Zhao
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| |
Collapse
|
20
|
Studies on hydrolysis mechanism of anticancer ruthenium drug ImH[trans-Ru(Im)2Cl4] via ABEEMσπ polarizable force field combined with QM and MD-FEP. Chem Res Chin Univ 2017. [DOI: 10.1007/s40242-017-6401-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
21
|
Li H, Sun TT, Zhang C, Liu L, Zhao D, Yang Z. QM/MM(ABEEM) Study on the Ligand Substitution Processes of Ruthenium(III) Complex NAMI-A. CHINESE J CHEM 2017. [DOI: 10.1002/cjoc.201600605] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Hui Li
- Chemistry and Chemical Engineering Faculty; Liaoning Normal University; Dalian Liaoning 116029 China
| | - Ting-Ting Sun
- Chemistry and Chemical Engineering Faculty; Liaoning Normal University; Dalian Liaoning 116029 China
| | - Chao Zhang
- Chemistry and Chemical Engineering Faculty; Liaoning Normal University; Dalian Liaoning 116029 China
| | - Linlin Liu
- Chemistry and Chemical Engineering Faculty; Liaoning Normal University; Dalian Liaoning 116029 China
| | - Dongxia Zhao
- Chemistry and Chemical Engineering Faculty; Liaoning Normal University; Dalian Liaoning 116029 China
| | - Zhongzhi Yang
- Chemistry and Chemical Engineering Faculty; Liaoning Normal University; Dalian Liaoning 116029 China
| |
Collapse
|
22
|
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
|
23
|
Fast calculation of molecular total energy with ABEEMσπ/MM method – For some series of organic molecules and peptides. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.03.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
24
|
Xu S, Zhao DX, Gong LD, Liu C, Yang ZZ. Search of the conformations of Val-dipeptide and Val-tripeptide by ab initio method and ABEEMσπ polarizable force field. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2014.10.073] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|