Preparation of Hydrophobic Metal-Organic Frameworks via Plasma Enhanced Chemical Vapor Deposition of Perfluoroalkanes for the Removal of Ammonia

Plasma enhanced chemical vapor deposition (PECVD) of perfluoroalkanes has long been studied for tuning the wetting properties of surfaces. For high surface area microporous materials, such as metal-organic frameworks (MOFs), unique challenges present themselves for PECVD treatments. Herein the protocol for development of a MOF that was previously unstable to humid conditions is presented. The protocol describes the synthesis of Cu-BTC (also known as HKUST-1), the treatment of Cu-BTC with PECVD of perfluoroalkanes, the aging of materials under humid conditions, and the subsequent ammonia microbreakthrough experiments on milligram quantities of microporous materials. Cu-BTC has an extremely high surface area (~1,800 m²/g) when compared to most materials or surfaces that have been previously treated by PECVD methods. Parameters such as chamber pressure and treatment time are extremely important to ensure the perfluoroalkane plasma penetrates to and reacts with the inner MOF surfaces. Furthermore, the protocol for ammonia microbreakthrough experiments set forth here can be utilized for a variety of test gases and microporous materials.     

Molecular simulations of physical and chemical adsorption under gas and liquid environments using force field- and quantum mechanics-based methods

Here we review our simulations of adsorption on metal–organic frameworks (MOFs) and platinum (Pt) catalysts, focusing on the modelling methods required to understand these two very different systems. MOFs are porous, crystalline materials with large surface areas, which are promising for a variety of adsorption applications. We review our simulations of gas uptake in PCN-53 (porous coordination network) as well as gas storage in MOFs functionalised with metal alkoxide sites. While fluid–solid interactions in both systems can be modelled quite well using algebraic force fields, the alkoxide sites in the functionalised MOFs require specialised versions, in order to describe the stronger adsorption energies. We discuss grand canonical Monte Carlo (GCMC) simulations of both systems. Pt is a common catalyst, and simulations have proven quite useful for providing molecular level details to understand its functionality. This involves understanding adsorption phenomena, which often requires quantum mechanical calculations. We describe our periodic boundary condition density functional theory (DFT) simulations of Pt-catalysed NO oxidation, focusing on adsorbate geometries and coverage effects. Finally, we describe one of the current ‘grand challenges’ in molecular simulations of adsorption, modelling catalytic activity in aqueous phase, which requires a combination of algebraic force fields, DFT and GCMC.     

Motional timescale predictions by molecular dynamics simulations: Case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin‐lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.     

Water activity coefficients in aqueous amino acid solutions by molecular dynamics simulation: 1. Force field development

New force fields for molecular dynamics (MD) simulation of aqueous zwitterionic amino acid simulations were developed. These were especially designed to calculate activity coefficient of water in amino acid solutions with high accuracy. For example, aqueous solutions of the following amino acids were considered: glycine, alanine, α-aminobutyric acid, α-aminovalerianic acid, valine and leucine. The force fields were obtained by quantum chemical calculations using B3LYP/6-31G and MP2/6-311(d,p) model theories in combination with the Merz–Kollmann–Singh scheme. To further increase the accuracy of the force field, a polarised continuum was considered in all quantum chemical calculations. Water activity coefficients obtained from MD using different all-purpose literature force fields, namely, OPLS, AMBER ff03 and GROMOS 53A6 as well as experimental data are compared with the results utilising the new force field. The new force field is shown to give better results compared with experimental data than existing force fields.     

Simulating the vapour–liquid equilibria of 1,4-dioxane

1,4-dioxane, a cyclic ether, is an emerging contaminant which is difficult to remove from water with conventional water treatment methods and resistant to biodegradation. Once a reliable force field is developed for 1,4-dioxane, molecular simulation techniques can be useful to study alternative adsorbents for its removal. For this purpose, we carried out Monte Carlo simulations in a constant volume Gibbs Ensemble to generate a force field which is capable of predicting the vapour–liquid coexistence curve and critical data of 1,4-dioxane. Results are given in comparison with experimental data and results from simulations with other force fields. Liquid densities and critical temperature are predicted in excellent agreement with experimental data using the new force field. At high temperatures, predicted vapour densities are in good agreement with experimental data, however, at lower temperatures the predicted vapour densities deviate about an order of magnitude from the experimental values. The critical density is slightly underestimated with our new force field. However, overall, the results of simulations with the new parameters give much better agreement with experimental data compared to the results obtained using other force fields.     

Tetracycline and its analogues as inhibitors of amyloid fibrils: searching for a geometrical pharmacophore by theoretical investigation of their conformational behavior in aqueous solution

Tetracycline (TC) and its derivatives have recently been proposed as a new class of antagonists in prion diseases as they prevent the aggregation of prion protein peptides and their acquisition of protease resistance in vitro and in vivo. Looking for relationships between conformational flexibility and biological activity, we searched for a geometrical pharmacophore by investigating, in aqueous solution, the conformational behavior of 15 TCs in both the zwitterionic and the anionic forms. For TC similar conformational flexibility was found for the two forms and two main conformational families were detected, an extended and a folded conformation characterized by different intramolecular hydrogen-bond networks. On comparing the Molecular Mechanics results with the ab initio ones and the experimental evidence, it can be seen that the conformational behavior of TC is reasonably well predicted by the MM2 force field, whereas the conformational energies provided by the Amber force field are unreliable. The conformational analysis of the other TC derivatives was then performed by the MM2 force field. As a result, their conformational behavior was similar to that observed for TC itself. Despite the hydronaphthacene moiety's conformational flexibility, no geometrical pharmacophore was found among the TCs, i.e. properties other than geometrical ones should play a crucial role in determining their anti-fibrillogenic ability.     

Metal–Organic Framework Composites for Theragnostics and Drug Delivery Applications

Among a plethora of nano-sized therapeutics, metal-organic frameworks (MOFs) have been some of the most investigated novel materials for, predominantly, cancer drug delivery applications. Due to their large drug uptake capacities and slow-release mechanisms, MOFs are desirable drug delivery vehicles that protect and transport sensitive drug molecules to target sites. The inclusion of other guest materials into MOFs to make MOF-composite materials has added further functionality, from externally triggered drug release to improved pharmacokinetics and diagnostic aids. MOF-composites are synthetically versatile and can include examples such as magnetic nanoparticles in MOFs for MRI image contrast and polymer coatings that improve the blood-circulation time. From synthesis to applications, this review will consider the main developments in MOF-composite chemistry for biomedical applications and demonstrate the potential of these novel agents in nanomedicine. It is concluded that, although vast synthetic progress has been made in the field, it requires now to develop more biomedical expertise with a focus on rational model selection, a major comparative toxicity study, and advanced targeting techniques.     

Magnetic Field–Suppressed Lithium Dendrite Growth for Stable Lithium‐Metal Batteries

Lithium metal is the most attractive anode material due to its extremely high specific capacity, minimum potential, and low density. However, uncontrollable growth of lithium dendrite results in severe safety and cycling stability concerns, which hinders the application in next generation secondary batteries. In this paper, a new and facile method imposing a magnetic field to lithium metal anodes is proposed. That is, the lithium ions suffering Lorentz force due to the electromagnetic fields are put into spiral motion causing magnetohydrodynamics (MHD) effect. This MHD effect can effectively promote mass transfer and uniform distribution of lithium ions to suppress the dendrite growth as well as obtain uniform and compact lithium deposition. The results show that the lithium metal electrodes within the magnetic field exhibit excellent cycling and rate performance in a symmetrical battery. Additionally, full batteries using limited lithium metal as anodes and commercial LiFePO₄ as cathodes show improved performance within the magnetic field. In summary, a new and facile strategy of suppressing lithium dendrites using the MHD effect by imposing a magnetic field is proposed, which may be generalized to other advanced alkali metal batteries.     

Development of an AMBER-compatible transferable force field for poly(ethylene glycol) ethers (glymes)

An all-atom force field consistent with the general AMBER force field (GAFF) format for poly(ethylene glycol) dimethyl ether (diglyme or G2) was developed by fitting to experimental liquid densities and dielectric constants. Not surprisingly, the new force field gives excellent agreement with experimental liquid phase densities and dielectric constants over a wide temperature range. Other dynamic and thermodynamic properties of liquid G2 such as its self-diffusion coefficient, shear viscosity, and vaporization enthalpy were also calculated and compared to experimental data. For all of the properties studied, the performance of the proposed new force field is better than that of the standard GAFF force field. The force field parameters were transferred to model two other poly(ethylene glycol) ethers: monoglyme (G1) and tetraglyme (G4). The predictive ability of the modified force field for G1 and G4 was significantly better than that of the original GAFF force field. The proposed force field provides an alternative option for the simulation of mixtures containing glymes using GAFF-compatible force fields, particularly for electrochemical applications. The accuracy of a previously published force field based on the OPLS-AA format and the accuracies of two modified versions of that force field were also examined for G1, G2, and G4. It was found that the original OPLS-AA force field is superior to the modified versions of it, and that it has a similar accuracy to the proposed new GAFF-compatible force field.

Open image in new window

Graphical abstract Transferability of an AMBER-compatible force field parameterized for G2 to other glymes

     

Why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution when comparing with the GROMOS 43A1 force field?

The optimal combination of force field and water model is an essential problem that is able to increase molecular dynamics simulation quality for different types of proteins and peptides. In this work, an attempt has been made to explore the problem by studying H1 peptide using four different models based on different force fields, water models and electrostatic schemes. The driving force for H1 peptide conformation transition and the reason why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution while the GROMOS 43A1 force field can do were investigated by temperature replica exchange molecular dynamics simulation (T-REMD). The simulation using the GROMOS 43A1 force field preferred to adopt a β-hairpin structure, which was in good agreement with the several other simulations and the experimental evidences. However, the simulation using the OPLS-AA force field has a significant difference from the simulations with the GROMOS 43A1 force field simulation. The results show that the driving force in H1 peptide conformation transition is solvent exposure of its hydrophobic residues. However, the subtle balances between residue-residue interactions and residue-solvent interaction are disrupted by using the OPLS-AA force field, which induced the reduction in the number of residue-residue contact. Similar solvent exposure of the hydrophobic residues is observed for all the conformations sampled using the OPLS-AA force field. For H1 peptide which exhibits large solvent exposure of the hydrophobic residues, the GROMOS 43A1 force field with the SPC water model can provide more accurate results.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Molecular simulations in metal–organic frameworks for diverse potential applications

As a new class of intriguing nanoporous materials, metal–organic frameworks (MOFs) have been considered for diverse potential applications. The number of MOFs synthesised to date is extremely large; thus, experimental testing alone is economically expensive and practically formidable. With rapidly growing computational resources, molecular simulation has become an indispensable tool to characterise, screen and design MOFs. Our research group has conducted comprehensive simulation studies in MOFs for carbon capture, hydrocarbon separation, alcohol adsorption and biofuel purification, water treatment, bio- and chiral separation and drug loading; furthermore, mechanical moduli, structural transition and thermal conductivity have also been examined. These systematic simulation studies are summarised in this review to demonstrate that simulation at a molecular level can secure the quantitative interpretation of experimental observation, provide the microscopic insight from bottom-up and facilitate the development of new MOFs.     

基于金属有机骨架的生物复合物抗逆性研究进展

金属有机骨架 (Metal-organic frameworks,MOFs) 由金属离子/簇和有机配体通过自组装形成,在催化、传感、能源和生物医药等领域有十分广泛的应用。近年来,金属有机骨架和生物活性物质结合引起了许多研究者的关注。金属有机骨架具有超大孔容、高比表面积、多样性结构组成等优势,使其可作为固定载体保护生物活性物质免受外界不良微环境的影响,包括高温、高压、有机溶剂等,从而提高生物活性物质的抗逆性。文中从金属有机骨架作为保护层提高不同种类生物活性物质抗逆性的角度进行了综述,并贯穿介绍了基于金属有机骨架的生物复合物的合成策略,旨在为相关领域的研究人员提供参考,促进基于金属有机骨架的生物复合物的实际应用。     

Concept and prototype of protein-ligand docking simulator with force feedback technology

A novel concept for a protein-ligand docking simulator using Virtual Reality (VR) technologies, in particular the tactile sense technology, was designed and a prototype was developed. Most conventional docking simulators are based on numerical differential calculations of the total energy between a protein and a ligand. However, the basic concept of our method differs from that of conventional simulators. Our design utilizes the force between a ligand and a protein instead of the total energy. The most characteristic function of the system is its ability to enable the user to 'touch' and sense the electrostatic potential field of a protein molecule. The user can scan the surface of a protein using a globular probe, which is given an electrostatic charge, and is controlled by a force feedback device. The electrostatic force between the protein and the probe is calculated in real time and immediately fed back into the force feedback device. The user can easily search interactively for positions where the probe is strongly attracted to the force field. Such positions can be regarded as candidate sites where functional groups of ligands corresponding to the probe can bind to the target protein. Certain limitations remain; for example, only twenty protein atoms can be used to generate the electrostatic field. Furthermore, the system can only use globular probes, preventing drug molecules or small chemical groups from being simulated. These limitations are the result of our insufficient computer resources. However, our prototype system has the potential to become a novel application method as well as being applicable to conventional VR technologies, especially to force feedback technologies.     

A new pairwise folding potential based on improved decoy generation and side-chain packing

A new force field for pairwise residue interactions as a function of C(alpha) to C(alpha) distances is presented. The force field was developed through the solution of a linear programming formulation with large sets of constraints. The constraints are based on the construction of >80,000 low-energy decoys for a set of proteins and requiring the decoy energies for each protein system to be higher than the native conformation of that particular protein. The generation of a robust force field was facilitated by the use of a novel decoy generation process, which involved the rational selection of proteins to add to the training set and included a significant energy minimization of the decoys. The force field was tested on a large set of decoys for various proteins not included in the training set and shown to perform well compared with a leading force field in identifying the native conformation for these proteins.     

Temperature and urea induced denaturation of the TRP‐cage mini protein TC5b: A simulation study consistent with experimental observations

The effects of temperature and urea denaturation (6M urea) on the dominant structures of the 20‐residue Trp‐cage mini‐protein TC5b are investigated by molecular dynamics simulations of the protein at different temperatures in aqueous and in 6M urea solution using explicit solvent degrees of freedom and the GROMOS force‐field parameter set 45A3. In aqueous solution at 278 K, TC5b is stable throughout the 20 ns of MD simulation and the trajectory structures largely agree with the NMR‐NOE atom–atom distance data available. Raising the temperature to 360 K and to 400 K, the protein denatures within 22 ns and 3 ns, showing that the denaturation temperature is well below 360 K using the GROMOS force field. This is 40–90 K lower than the denaturation temperatures observed in simulations using other much used protein force fields. As the experimental denaturation temperature is about 315 K, the GROMOS force field appears not to overstabilize TC5b, as other force fields and the use of continuum solvation models seem to do. This feature may directly stem from the GROMOS force‐field parameter calibration protocol, which primarily involves reproduction of condensed phase thermodynamic quantities such as energies, densities, and solvation free energies of small compounds representative for protein fragments. By adding 6M urea to the solution, the onset of denaturation is observed in the simulation, but is too slow to observe a particular side‐chain side‐chain contact (Trp6‐Ile4) that was experimentally observed to be characteristic for the denatured state. Interestingly, using temperature denaturation, the process is accelerated and the experimental data are reproduced.     

Increasing the precision of comparative models with YASARA NOVA--a self-parameterizing force field

One of the conclusions drawn at the CASP4 meeting in Asilomar was that applying various force fields during refinement of template-based models tends to move predictions in the wrong direction, away from the experimentally determined coordinates. We have derived an all-atom force field aimed at protein and nucleotide optimization in vacuo (NOVA), which has been specifically designed to avoid this problem. NOVA resembles common molecular dynamics force fields but has been automatically parameterized with two major goals: (i) not to make high resolution X-ray structures worse and (ii) to improve homology models built by WHAT IF. Force-field parameters were not required to be physically correct; instead, they were optimized with random Monte Carlo moves in force-field parameter space, each one evaluated by simulated annealing runs of a 50-protein optimization set. Errors inherent to the approximate force-field equation could thus be canceled by errors in force-field parameters. Compared with the optimization set, the force field did equally well on an independent validation set and is shown to move in silico models closer to reality. It can be applied to modeling applications as well as X-ray and NMR structure refinement. A new method to assign force-field parameters based on molecular trees is also presented. A NOVA server is freely accessible at http://www.yasara.com/servers     

MD simulation of subtilisin BPN' in a crystal environment.

In this paper we present a molecular dynamics (MD) simulation of subtilisin BPN' in a crystalline environment containing four protein molecules and solvent. Conformational and dynamic properties of the molecules are compared with each other and with respect to the X-ray structure to test the validity of the force field. The agreement between simulated and experimental structure using the GROMOS force field is better than that obtained in the literature using other force fields for protein crystals. The overall shape of the molecule is well preserved, as is the conformation of alpha-helices and beta-strands. Structural differences are mainly found in loop regions. Solvent networks found in the X-ray structure were reproduced by the simulation, which was unbiased with respect to the crystalline hydration structure. These networks seem to play an important role in the stability of the protein; evidence of this is found in the structure of the active site. The weak ion binding site in the X-ray structure of subtilisin BPN' is occupied by a monovalent ion. When a calcium ion is placed in the initial structure, three peptide ligands are replaced by 5 water ligands, whereas a potassium ion retains (in part) its original ligands. Existing force fields yield a reliable method to probe local structure and short-time dynamics of proteins, providing an accuracy of about 0.1 nm.     

Fabrication of Ordered Macro‐Microporous Single‐Crystalline MOF and Its Derivative Carbon Material for Supercapacitor

Because of their good performance in diffusion‐limited processes, ordered macro‐microporous single‐crystalline metal‐organic frameworks (MOFs) have potential for use in various fields. However, there are still very few reports of the synthesis of such MOFs. A general synthesis methodology for ordered macro‐microporous single‐crystalline MOFs is highly desired. Here, a novel strategy is reported for synthesizing single‐crystalline ordered macro‐microporous MOFs by monodentate‐ligand‐induced in situ crystallization within a 3D ordered hard template in a double‐solvent system. A space‐confined growth model is proposed to clarify the shaping effect of the template; the role of the monodentate ligand is also analyzed. Moreover, a carbon material derived from the macro‐microporous MOF inherits the ordered interconnected macroporous structure. The improved diffusion and lower resistance, as well as the structural robustness, endow the derivative carbon material with superior rate performance and excellent cycling stability when prepared as electrodes for a supercapacitor. It is anticipated that the method will provide new paths to the synthesis of such macro‐microporous materials for applications in energy‐related fields and beyond.