1
|
Ghanavati R, Escobosa AC, Manz TA. An automated protocol to construct flexibility parameters for classical forcefields: applications to metal-organic frameworks. RSC Adv 2024; 14:22714-22762. [PMID: 39035129 PMCID: PMC11258866 DOI: 10.1039/d4ra01859a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2024] [Accepted: 06/18/2024] [Indexed: 07/23/2024] Open
Abstract
In this work, forcefield flexibility parameters were constructed and validated for more than 100 metal-organic frameworks (MOFs). We used atom typing to identify bond types, angle types, and dihedral types associated with bond stretches, angle bends, dihedral torsions, and other flexibility interactions. Our work used Manz's angle-bending and dihedral-torsion model potentials. For a crystal structure containing N atoms in its unit cell, the number of independent flexibility interactions is 3(N atoms - 1). Because the number of bonds, angles, and dihedrals is normally much larger than 3(N atoms - 1), these internal coordinates are redundant. To reduce (but not eliminate) this redundancy, our protocol prunes dihedral types in a way that preserves symmetry equivalency. Next, each dihedral type is classified as non-rotatable, hindered, rotatable, or linear. We introduce a smart selection method that identifies which particular torsion modes are important for each rotatable dihedral type. Then, we computed the force constants for all flexibility interactions together via LASSO regression (i.e., regularized linear least-squares fitting) of the training dataset. LASSO automatically identifies and removes unimportant forcefield interactions. For each MOF, the reference dataset was quantum-mechanically-computed in VASP via DFT with dispersion and included: (i) finite-displacement calculations along every independent atom translation mode, (ii) geometries randomly sampled via ab initio molecular dynamics (AIMD), (iii) the optimized ground-state geometry using experimental lattice parameters, and (iv) rigid torsion scans for each rotatable dihedral type. After training, the flexibility model was validated across geometries that were not part of the training dataset. For each MOF, we computed the goodness of fit (R-squared value) and the root-mean-squared error (RMSE) separately for the training and validation datasets. We compared flexibility models with and without bond-bond cross terms. Even without cross terms, the model yielded R-squared values of 0.910 (avg across all MOFs) ± 0.018 (st. dev.) for atom-in-material forces in the validation datasets. Our SAVESTEPS protocol should find widespread applications to parameterize flexible forcefields for material datasets. We performed molecular dynamics simulations using these flexibility parameters to compute heat capacities and thermal expansion coefficients for two MOFs.
Collapse
Affiliation(s)
- Reza Ghanavati
- Chemical & Materials Engineering, New Mexico State University Las Cruces NM 88001 USA
| | - Alma C Escobosa
- Chemical & Materials Engineering, New Mexico State University Las Cruces NM 88001 USA
| | - Thomas A Manz
- Chemical & Materials Engineering, New Mexico State University Las Cruces NM 88001 USA
| |
Collapse
|
2
|
Maity S, Gangan AS, Anshu A, V Valappil RR, Chakraborty B, Ramaniah LM, Srinivasan V. Pressure induced topochemical polymerization of solid acrylamide facilitated by anisotropic response of the hydrogen bond network. Phys Chem Chem Phys 2021; 23:9448-9456. [PMID: 33885052 DOI: 10.1039/d0cp04993j] [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 pressure induced polymerization of molecular solids is an appealing route to obtain pure, crystalline polymers without the need for radical initiators. Here, we report a detailed density functional theory (DFT) study of the structural and chemical changes that occur in defect free solid acrylamide, a hydrogen bonded crystal, when it is subjected to hydrostatic pressures. While our calculations are able to reproduce experimentally measured pressure dependent spectroscopic features in the 0-20 GPa range, our atomistic analysis predicts polymerization in acrylamide at a pressure of ∼23 GPa at 0 K albeit through large enthalpy barriers. Interestingly, we find that the two-dimensional hydrogen bond network in acrylamide templates topochemical polymerization by aligning the atoms through an anisotropic response at low pressures. This results not only in conventional C-C, but also unusual C-O polymeric linkages, as well as a new hydrogen bonded framework, with both N-HO and C-HO bonds. Using a simple model for thermal effects, we also show that at 300 K, higher pressures significantly accelerate the transformation into polymers by lowering the barrier. Thus, application of pressure offers an alternative route for topochemical polymerization when higher temperatures are undesirable.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal 462 066, India.
| | | | | | | | | | | | | |
Collapse
|
3
|
Rana RM, Rampogu S, Abid NB, Zeb A, Parate S, Lee G, Yoon S, Kim Y, Kim D, Lee KW. In Silico Study Identified Methotrexate Analog as Potential Inhibitor of Drug Resistant Human Dihydrofolate Reductase for Cancer Therapeutics. Molecules 2020; 25:molecules25153510. [PMID: 32752079 PMCID: PMC7435474 DOI: 10.3390/molecules25153510] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Revised: 07/28/2020] [Accepted: 07/30/2020] [Indexed: 11/16/2022] Open
Abstract
Drug resistance is a core issue in cancer chemotherapy. A known folate antagonist, methotrexate (MTX) inhibits human dihydrofolate reductase (hDHFR), the enzyme responsible for the catalysis of 7,8-dihydrofolate reduction to 5,6,7,8-tetrahydrofolate, in biosynthesis and cell proliferation. Structural change in the DHFR enzyme is a significant cause of resistance and the subsequent loss of MTX. In the current study, wild type hDHFR and double mutant (engineered variant) F31R/Q35E (PDB ID: 3EIG) were subject to computational study. Structure-based pharmacophore modeling was carried out for wild type (WT) and mutant (MT) (variant F31R/Q35E) hDHFR structures by generating ten models for each. Two pharmacophore models, WT-pharma and MT-pharma, were selected for further computations, and showed excellent ROC curve quality. Additionally, the selected pharmacophore models were validated by the Guner-Henry decoy test method, which yielded high goodness of fit for WT-hDHFR and MT-hDHFR. Using a SMILES string of MTX in ZINC15 with the selections of 'clean', in vitro and in vivo options, 32 MTX-analogs were obtained. Eight analogs were filtered out due to their drug-like properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment tests and Lipinski's Rule of five. WT-pharma and MT-pharma were further employed as a 3D query in virtual screening with drug-like MTX analogs. Subsequently, seven screening hits along with a reference compound (MTX) were subjected to molecular docking in the active site of WT- and MT-hDHFR. Through a clustering analysis and examination of protein-ligand interactions, one compound was found with a ChemPLP fitness score greater than that of MTX (reference compound). Finally, a simulation of molecular dynamics (MD) identified an MTX analog which exhibited strong affinity for WT- and MT-hDHFR, with stable RMSD, hydrogen bonds (H-bonds) in the binding site and the lowest MM/PBSA binding free energy. In conclusion, we report on an MTX analog which is capable of inhibiting hDHFR in wild type form, as well as in cases where the enzyme acquires resistance to drugs during chemotherapy treatment.
Collapse
Affiliation(s)
- Rabia Mukhtar Rana
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Shailima Rampogu
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Noman Bin Abid
- Division of Life Science and Applied Life Science (BK 21), College of Natural Sciences, Gyeongsang National University, Jinju 52828, Korea;
| | - Amir Zeb
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Shraddha Parate
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Gihwan Lee
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Sanghwa Yoon
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Yumi Kim
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Donghwan Kim
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
| | - Keun Woo Lee
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea; (R.M.R.); (S.R.); (A.Z.); (S.P.); (G.L.); (S.Y.); (Y.K.); (D.K.)
- Correspondence: ; Tel.: +82-55-772-1360
| |
Collapse
|
4
|
Abstract
ConspectusLiquids under confinement differ in behavior from their bulk counterparts and can acquire properties that are specific to the confined phase and linked to the nature and structure of the host matrix. While confined liquid water is not a new topic of research, the past few years have seen a series of intriguing novel features for water inside nanoscale pores. These unusual properties arise from the very specific nature of nanoporous materials, termed "soft porous crystals"; they combine large-scale flexibility with a heterogeneous internal surface. This creates a rich diversity of behavior for the adsorbed water, and the combination of different experimental characterization techniques along with computational chemistry at various scales is necessary to understand the phenomena observed and their microscopic origins. The range of systems of interest span the whole chemical range, from the inorganic (zeolites, imogolites) to the organic (microporous carbons, graphene, and its derivatives), and even encompass the hybrid organic-inorganic systems (such as metal-organic frameworks).The combination of large scale flexibility with the strong physisorption (or even chemisorption) of water can lead to unusual properties (belonging to the "metamaterials" category) and to novel phenomena. One striking example is the recent elucidation of the mechanism of negative hydration expansion in ZrW2O8, by which adsorption of ∼10 wt % water in the inorganic nonporous framework leads to large shrinkage of its volume. Another eye-catching case is the occurrence of multiple water adsorption-driven structural transitions in the MIL-53 family of materials: the specific interactions between water guest molecules and the host framework create behavior that has not been observed with any other adsorbate. Both are counterintuitive phenomena that have been elucidated by a combination of experimental in situ techniques and molecular simulation.Another important direction of research is the shift in the systems and phenomena studied, from physical adsorption toward studies of reactivity, hydrothermal stability, and the effect of confinement on aqueous phases more complex than pure water. There have been examples of water adsorption in highly flexible metal-organic frameworks being able to compete with the materials' coordination bonds, thereby limiting its hydrothermal stability, while tweaking the functional groups of the same framework can lead to increased stability while retaining the flexibility of the material. However, this additional complexity and tunability in the macroscopic behavior can occur from changes in the confined fluid rather than the material. Very recent studies have shown that aqueous solutions of high concentration (such as LiCl up to 20 mol L-1) confined in flexible nanoporous materials can have specific properties different from pure water and not entirely explained by osmotic effects. There, the strong ordering of the confined electrolyte competes with the structural flexibility of the framework to create an entirely new behavior for the {host, guest} system.
Collapse
Affiliation(s)
- François-Xavier Coudert
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France
| |
Collapse
|
5
|
Dürholt JP, Fraux G, Coudert FX, Schmid R. Ab Initio Derived Force Fields for Zeolitic Imidazolate Frameworks: MOF-FF for ZIFs. J Chem Theory Comput 2019; 15:2420-2432. [DOI: 10.1021/acs.jctc.8b01041] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Johannes P. Dürholt
- Computational Materials Chemistry group, Fakultät für Chemie und Biochemie, Ruhr-Universität Bochum, Bochum 44801, Germany
| | - Guillaume Fraux
- Chimie ParisTech, PSL University,
CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France
| | - François-Xavier Coudert
- Chimie ParisTech, PSL University,
CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France
| | - Rochus Schmid
- Computational Materials Chemistry group, Fakultät für Chemie und Biochemie, Ruhr-Universität Bochum, Bochum 44801, Germany
| |
Collapse
|
6
|
Rana RM, Rampogu S, Zeb A, Son M, Park C, Lee G, Yoon S, Baek A, Parameswaran S, Park SJ, Lee KW. In Silico Study Probes Potential Inhibitors of Human Dihydrofolate Reductase for Cancer Therapeutics. J Clin Med 2019; 8:jcm8020233. [PMID: 30754680 PMCID: PMC6406960 DOI: 10.3390/jcm8020233] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Revised: 02/05/2019] [Accepted: 02/08/2019] [Indexed: 01/10/2023] Open
Abstract
Dihydrofolate reductase (DHFR) is an essential cellular enzyme and thereby catalyzes thereduction of dihydrofolate to tetrahydrofolate (THF). In cancer medication, inhibition of humanDHFR (hDHFR) remains a promising strategy, as it depletes THF and slows DNA synthesis and cellproliferation. In the current study, ligand-based pharmacophore modeling identified and evaluatedthe critical chemical features of hDHFR inhibitors. A pharmacophore model (Hypo1) was generatedfrom known inhibitors of DHFR with a correlation coefficient (0.94), root mean square (RMS)deviation (0.99), and total cost value (125.28). Hypo1 was comprised of four chemical features,including two hydrogen bond donors (HDB), one hydrogen bond acceptor (HBA), and onehydrophobic (HYP). Hypo1 was validated using Fischer's randomization, test set, and decoy setvalidations, employed as a 3D query in a virtual screening at Maybridge, Chembridge, Asinex,National Cancer Institute (NCI), and Zinc databases. Hypo1-retrieved compounds were filtered byan absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment test andLipinski's rule of five, where the drug-like hit compounds were identified. The hit compounds weredocked in the active site of hDHFR and compounds with Goldfitness score was greater than 44.67(docking score for the reference compound), clustering analysis, and hydrogen bond interactionswere identified. Furthermore, molecular dynamics (MD) simulation identified three compounds asthe best inhibitors of hDHFR with the lowest root mean square deviation (1.2 Å to 1.8 Å), hydrogenbond interactions with hDHFR, and low binding free energy (-127 kJ/mol to -178 kJ/mol). Finally,the toxicity prediction by computer (TOPKAT) affirmed the safety of the novel inhibitors of hDHFRin human body. Overall, we recommend novel hit compounds of hDHFR for cancer and rheumatoidarthritis chemotherapeutics.
Collapse
Affiliation(s)
- Rabia Mukhtar Rana
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Shailima Rampogu
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Amir Zeb
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Minky Son
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Chanin Park
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Gihwan Lee
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Sanghwa Yoon
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Ayoung Baek
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Sarvanan Parameswaran
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| | - Seok Ju Park
- Department of Internal Medicine, College of Medicine, Busan Paik Hospital, Inje University,Busan 47392, Korea.
| | - Keun Woo Lee
- Division of Life Sciences, Division of Applied Life Science (BK21 Plus), Research Institute of NaturalScience (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Korea.
| |
Collapse
|
7
|
Millange F, Walton RI. MIL-53 and its Isoreticular Analogues: a Review of the Chemistry and Structure of a Prototypical Flexible Metal-Organic Framework. Isr J Chem 2018. [DOI: 10.1002/ijch.201800084] [Citation(s) in RCA: 49] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Franck Millange
- Département de Chimie; Université de Versailles-St-Quentin-en-Yvelines; 45 Avenue des États-Unis 78035 Versailles cedex France
| | - Richard I. Walton
- Department of Chemistry; University of Warwick; Gibbet Hill Road Coventry CV4 7AL United Kingdom
| |
Collapse
|
8
|
Akimbekov Z, Katsenis AD, Nagabhushana GP, Ayoub G, Arhangelskis M, Morris AJ, Friščić T, Navrotsky A. Experimental and Theoretical Evaluation of the Stability of True MOF Polymorphs Explains Their Mechanochemical Interconversions. J Am Chem Soc 2017; 139:7952-7957. [DOI: 10.1021/jacs.7b03144] [Citation(s) in RCA: 73] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Zamirbek Akimbekov
- Peter
A. Rock Thermochemistry Laboratory and NEAT ORU, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Athanassios D. Katsenis
- Department
of Chemistry, McGill University, 801 Sherbrooke Street West, Montreal, Canada H3A
0B8
| | - G. P. Nagabhushana
- Peter
A. Rock Thermochemistry Laboratory and NEAT ORU, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Ghada Ayoub
- Department
of Chemistry, McGill University, 801 Sherbrooke Street West, Montreal, Canada H3A
0B8
| | - Mihails Arhangelskis
- Department
of Chemistry, McGill University, 801 Sherbrooke Street West, Montreal, Canada H3A
0B8
| | - Andrew J. Morris
- Cavendish
Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, U.K
- Department
of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, U.K
| | - Tomislav Friščić
- Department
of Chemistry, McGill University, 801 Sherbrooke Street West, Montreal, Canada H3A
0B8
| | - Alexandra Navrotsky
- Peter
A. Rock Thermochemistry Laboratory and NEAT ORU, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
9
|
Banlusan K, Strachan A. First-principles study of elastic mechanical responses to applied deformation of metal-organic frameworks. J Chem Phys 2017. [DOI: 10.1063/1.4982356] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Kiettipong Banlusan
- School of Materials Engineering and Birck Nanotechnology Center, Purdue University, West Lafayette, Indiana 47907, USA
| | - Alejandro Strachan
- School of Materials Engineering and Birck Nanotechnology Center, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
10
|
Cockayne E. Thermodynamics of the Flexible Metal-Organic Framework Material MIL-53(Cr) From First Principles. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2017; 121:4312-4317. [PMID: 28690713 PMCID: PMC5497523 DOI: 10.1021/acs.jpcc.6b11692] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
We use first-principles density functional theory total energy and linear response phonon calculations to compute the Helmholtz and Gibbs free energy as a function of temperature, pressure, and cell volume in the flexible metal-organic framework material MIL-53(Cr) within the quasiharmonic approximation. GGA and metaGGA calculations were performed, each including empirical van der Waals (vdW) forces under the D2, D3, or D3(BJ) parameterizations. At all temperatures up to 500 K and pressures from -30 MPa to 30 MPa, two minima in the free energy versus volume are found, corresponding to the narrow pore (np) and large pore (lp) structures. Critical positive and negative pressures are identified, beyond which there is only one free energy minimum. While all results overestimated the stability of the np phase relative to the lp phase, the best overall agreement with experiment is found for the metaGGA PBEsol+RTPSS+U+J approach with D3 or D3(BJ) vdW forces. For these parameterizations, the calculated free energy barrier for the np-lp transition is only 3 to 6 kJ per mole of Cr4(OH)4(C8H4O4)4.
Collapse
Affiliation(s)
- Eric Cockayne
- Materials Measurement Science Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 USA
| |
Collapse
|
11
|
Parkes MV, Greathouse JA, Hart DB, Gallis DFS, Nenoff TM. Ab initio molecular dynamics determination of competitive O₂ vs. N₂ adsorption at open metal sites of M₂(dobdc). Phys Chem Chem Phys 2017; 18:11528-38. [PMID: 27063148 DOI: 10.1039/c6cp00768f] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The separation of oxygen from nitrogen using metal-organic frameworks (MOFs) is of great interest for potential pressure-swing adsorption processes for the generation of purified O2 on industrial scales. This study uses ab initio molecular dynamics (AIMD) simulations to examine for the first time the pure-gas and competitive gas adsorption of O2 and N2 in the M2(dobdc) (M = Cr, Mn, Fe) MOF series with coordinatively unsaturated metal centers. Effects of metal, temperature, and gas composition are explored. This unique application of AIMD allows us to study in detail the adsorption/desorption processes and to visualize the process of multiple guests competitively binding to coordinatively unsaturated metal sites of a MOF.
Collapse
Affiliation(s)
- Marie V Parkes
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185-0754, USA
| | - Jeffery A Greathouse
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185-0754, USA
| | - David B Hart
- Geochemistry Department, Sandia National Laboratories, Albuquerque, New Mexico 87185-0754, USA
| | - Dorina F Sava Gallis
- Nanoscale Sciences Department, Sandia National Laboratories, Albuquerque, New Mexico 87185-1415, USA
| | - Tina M Nenoff
- Physical Chemical and Nano Sciences Center, Sandia National Laboratories, Albuquerque, New Mexico 87185-1415, USA.
| |
Collapse
|
12
|
Fraux G, Coudert FX. Recent advances in the computational chemistry of soft porous crystals. Chem Commun (Camb) 2017. [DOI: 10.1039/c7cc03306k] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We highlight recent progress in the field of computational chemistry of nanoporous materials, focusing on methods and studies that address the extraordinary dynamic nature of these systems: the high flexibility of their frameworks, the large-scale structural changes upon external physical or chemical stimulation, and the presence of defects and disorder.
Collapse
Affiliation(s)
- Guillaume Fraux
- Chimie ParisTech
- PSL Research University
- CNRS
- Institut de Recherche de Chimie Paris
- 75005 Paris
| | | |
Collapse
|
13
|
Rogge SM, Wieme J, Vanduyfhuys L, Vandenbrande S, Maurin G, Verstraelen T, Waroquier M, Van Speybroeck V. Thermodynamic Insight in the High-Pressure Behavior of UiO-66: Effect of Linker Defects and Linker Expansion. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2016; 28:5721-5732. [PMID: 27594765 PMCID: PMC5006632 DOI: 10.1021/acs.chemmater.6b01956] [Citation(s) in RCA: 56] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/14/2016] [Revised: 07/25/2016] [Indexed: 05/19/2023]
Abstract
In this Article, we present a molecular-level understanding of the experimentally observed loss of crystallinity in UiO-66-type metal-organic frameworks, including the pristine UiO-66 to -68 as well as defect-containing UiO-66 materials, under the influence of external pressure. This goal is achieved by constructing pressure-versus-volume profiles at finite temperatures using a thermodynamic approach relying on ab initio derived force fields. On the atomic level, the phenomenon is reflected in a sudden drop in the number of symmetry operators for the crystallographic unit cell because of the disordered displacement of the organic linkers with respect to the inorganic bricks. For the defect-containing samples, a reduced mechanical stability is observed, however, critically depending on the distribution of these defects throughout the material, hence demonstrating the importance of judiciously characterizing defects in these materials.
Collapse
Affiliation(s)
- Sven M.
J. Rogge
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Jelle Wieme
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Louis Vanduyfhuys
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Steven Vandenbrande
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Guillaume Maurin
- Institut
Charles Gerhardt Montpellier, Université
Montpellier 2, Place
E. Bataillon, 34095 Montpellier cedex 05, France
| | - Toon Verstraelen
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Michel Waroquier
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Veronique Van Speybroeck
- Center
for Molecular Modeling (CMM), Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
- E-mail:
| |
Collapse
|
14
|
Coudert FX, Fuchs AH, Neimark AV. Adsorption deformation of microporous composites. Dalton Trans 2016; 45:4136-40. [PMID: 26600091 DOI: 10.1039/c5dt03978a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
We study here the behavior of flexible adsorbent materials, or soft porous crystals, when used in practical applications as nanostructured composites such as core-shell particles or mixed matrix membranes. Based on simple models and the well-established laws of elasticity, we demonstrate how the presence of a binder results in an attenuation of the adsorption-induced stress and deformation. In the case where the adsorbent undergoes adsorption-induced structural transitions, such as the gate opening phenomenon occurring in some metal-organic frameworks, we show that the presence of the binder will result in shifts of the adsorption-induced transition pressures.
Collapse
Affiliation(s)
- François-Xavier Coudert
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France.
| | - Alain H Fuchs
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France.
| | - Alexander V Neimark
- Department of Chemical and Biochemical Engineering, Rutgers University, Piscataway, New Jersey 08854, USA.
| |
Collapse
|
15
|
Haigis V, Coudert FX, Vuilleumier R, Boutin A, Fuchs AH. Hydrothermal Breakdown of Flexible Metal-Organic Frameworks: A Study by First-Principles Molecular Dynamics. J Phys Chem Lett 2015; 6:4365-4370. [PMID: 26722973 DOI: 10.1021/acs.jpclett.5b01926] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Flexible metal-organic frameworks, also known as soft porous crystals, have been proposed for a vast number of technological applications, because they respond by large changes in structure and properties to small external stimuli, such as adsorption of guest molecules and changes in temperature or pressure. While this behavior is highly desirable in applications such as sensing and actuation, their extreme flexibility can also be synonymous with decreased stability. In particular, their performance in industrial environments is limited by a lack of stability at elevated temperatures and in the presence of water. Here, we use first-principles molecular dynamics to study the hydrothermal breakdown of soft porous crystals. Focusing on the material MIL-53(Ga), we show that the weak point of the structure is the bond between the metal center and the organic linker and elucidate the mechanism by which water lowers the activation free energy for the breakdown. This allows us to propose strategies for the synthesis of MOFs with increased heat and water stability.
Collapse
Affiliation(s)
- Volker Haigis
- École Normale Supérieure, PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - François-Xavier Coudert
- PSL Research University, Chimie ParisTech - CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France
| | - Rodolphe Vuilleumier
- École Normale Supérieure, PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Anne Boutin
- École Normale Supérieure, PSL Research University, Département de Chimie, Sorbonne Universités - UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 rue Lhomond, 75005 Paris, France
| | - Alain H Fuchs
- PSL Research University, Chimie ParisTech - CNRS, Institut de Recherche de Chimie Paris, 75005 Paris, France
| |
Collapse
|
16
|
Cliffe MJ, Hill JA, Murray CA, Coudert FX, Goodwin AL. Defect-dependent colossal negative thermal expansion in UiO-66(Hf) metal-organic framework. Phys Chem Chem Phys 2015; 17:11586-92. [PMID: 25866163 DOI: 10.1039/c5cp01307k] [Citation(s) in RCA: 80] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Thermally-densified hafnium terephthalate UiO-66(Hf) is shown to exhibit the strongest isotropic negative thermal expansion (NTE) effect yet reported for a metal-organic framework (MOF). Incorporation of correlated vacancy defects within the framework affects both the extent of thermal densification and the magnitude of NTE observed in the densified product. We thus demonstrate that defect inclusion can be used to tune systematically the physical behaviour of a MOF.
Collapse
Affiliation(s)
- Matthew J Cliffe
- Department of Chemistry, University of Oxford, Inorganic Chemistry Laboratory, South Parks Road, Oxford OX1 3QR, UK.
| | | | | | | | | |
Collapse
|
17
|
Odoh SO, Cramer CJ, Truhlar DG, Gagliardi L. Quantum-Chemical Characterization of the Properties and Reactivities of Metal–Organic Frameworks. Chem Rev 2015; 115:6051-111. [DOI: 10.1021/cr500551h] [Citation(s) in RCA: 206] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Samuel O. Odoh
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|