Aqueous Stability of Metal-Organic Frameworks Using ReaxFF-Based Metadynamics Simulations.

J Phys Chem B

Department of Materials Science and Engineering, The Pennsylvania State University, University Park, Pennsylvania 16802, United States.

Published: July 2023

Aqueous stability is a critical property for the application of metal-organic framework (MOF) materials in humid conditions. The sampling of the free energy surface for a water reaction is challenging due to a lack of a reactive force field. Here, we developed a ReaxFF force field for simulating the reaction of zeolitic imidazole frameworks (ZIFs) with water. We carried out metadynamics simulations based on ReaxFF to study the reaction of water with a few different types of MOFs. We also conducted an experimental water immersion test and characterized the XRD, TG, and gas adsorption properties of the MOFs before and after the immersion test. By considering the energy barrier for a hydrolysis reaction, the simulation results are in good agreement with the experiments. MOFs with open structures and large pores are found to be unstable in metadynamics simulations, where the water molecule can attack or bond with the metallic node relatively easily. In contrast, it is more difficult for water to attack the Zn atom in the ZnN tetrahedral structure of ZIFs. We also found that ZIFs with the -NO functional groups have higher water stability. Discrepancies between the metadynamics simulation and gas adsorption experiments have been accounted for from the phase/crystallinity change of the structure reflected in the X-ray diffraction and thermogravimetry analysis of the MOF samples.

Download full-text PDF

Source
http://dx.doi.org/10.1021/acs.jpcb.3c00563DOI Listing

Publication Analysis

Top Keywords

metadynamics simulations
12
aqueous stability
8
force field
8
immersion test
8
gas adsorption
8
water
7
stability metal-organic
4
metal-organic frameworks
4
frameworks reaxff-based
4
metadynamics
4

Similar Publications

Theoretical study of the formation of HO by lytic polysaccharide monooxygenases: the reaction mechanism depends on the type of reductant.

Chem Sci

January 2025

State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry College of Chemistry and Chemical Engineering, Xiamen University Xiamen 361005 China

Lytic polysaccharide monooxygenases (LPMOs) are a unique group of monocopper enzymes that exhibit remarkable ability to catalyze the oxidative cleavage of recalcitrant carbohydrate substrates, such as cellulose and chitin, by utilizing O or HO as the oxygen source. One of the key challenges in understanding the catalytic mechanism of LPMOs lies in deciphering how they activate dioxygen using diverse reductants. To shed light on this intricate process, we conducted in-depth investigations using quantum mechanical/molecular mechanical (QM/MM) metadynamics simulations, molecular dynamics (MD) simulations, and density functional theory (DFT) calculations.

View Article and Find Full Text PDF

Ligand-based cheminformatics and free energy-inspired molecular simulations for prioritizing and optimizing G-protein coupled receptor kinase-6 (GRK6) inhibitors in multiple myeloma treatment.

Comput Biol Chem

January 2025

Drug Discovery and Development Laboratory (DDD Lab), Department of Pharmaceutical Technology, Jadavpur University, Kolkata 700032, India. Electronic address:

Multiple myeloma (MM) is the second most frequently diagnosed hematological malignancy, presenting limited treatment options with no curative potential and significant drug resistance. Recent studies involving genetic knockdown established the crucial role of GRK6 in upholding the viability of MM cells, emphasizing the need to identify potential inhibitors. Computational exploration of GRK6 inhibitors has not been attempted previously.

View Article and Find Full Text PDF

Computational Explorations of Th First Hydrolysis Reaction Constants: Insights from Molecular Dynamics and Density Functional Theory Calculations.

J Phys Chem A

January 2025

Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of Ministry of Education, Tsinghua University, Beijing 100084, China.

The fundamental hydrolysis behavior of tetravalent actinide cations (An) with a high charge is crucial for understanding their solution chemistry, particularly in nuclear fuel reprocessing and environmental behavior. Using Th as a reference of the An series, this work employed both the periodic model and the cluster model to calculate the first hydrolysis reaction constant (p) of the Th aqua ion and conducted a detailed evaluation of these approaches. In the periodic model, molecular dynamics (AIMD) simulations of Th in the explicit solvation environment are conducted, using metadynamics and constrained molecular dynamics to calculate p values.

View Article and Find Full Text PDF

JAK3 Inhibitors: Covalent and Noncovalent Interactions of a Cyanamide Group Investigated by Multiscale Free-Energy Simulations.

J Chem Inf Model

January 2025

Medicinal Chemistry and Drug Design Technologies Department, Chiesi Farmaceutici S.p.A., Largo F. Belloli 11/A, 43122 Parma, Italy.

Janus kinase type 3 (JAK3), an emerging target for treating autoimmune diseases, possesses a front pocket cysteine that is targeted by covalent modifiers, best represented by the marketed drug ritlecitinib (). Recently, 2,3-dihydro-1-inden-1-ylcyanamides have been developed as novel JAK3 inhibitors. Among them, the -(6-(7-pyrrolo[2,3-]pyrimidin-4-yl)-2,3-dihydro-1-inden-1-yl)cyanamide inhibitor () and its methylated analogue (), while being potent inhibitors, displayed different mechanisms of action (covalent vs noncovalent) and binding modes (Casimiro-Garcia et al.

View Article and Find Full Text PDF

Cannabinoid receptor 1 (CB1R) has been extensively studied as a potential therapeutic target for various conditions, including pain management, obesity, emesis, and metabolic syndrome. Unlike orthosteric agonists such as Δ-tetrahydrocannabinol (THC), cannabidiol (CBD) has been identified as a negative allosteric modulator (NAM) of CB1R, among its other pharmacological targets. Previous computational and structural studies have proposed various binding sites for CB1R NAMs.

View Article and Find Full Text PDF

Want AI Summaries of new PubMed Abstracts delivered to your In-box?

Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!