Modelling Studies of Crystalline PEEK

Abstract

The x-ray diffraction spectrum of the engineering thermoplastic poly (aryl-ether-ether-ketone), commonly known as PEEK, has been reported in the literature. The unit cell dimensions are consistent with two aryl units per unit cell, while the chemical repeat has three, suggesting that there is not a strict lateral register between chains. The polymer chain packing of crystalline PEEK has been investigated using lattice energy calculations and it is concluded that lateral chain ordering effects are rather unimportant in determining the crystal structure. The crystallographic equivalence of ether and ketone linkages arises due to essentially random chain stacking.     

First Principles Modeling of Dissociative Adsorption at Crystal Surfaces: Hydrogen on Pt(111)

Multiscale simulation has the potential of becoming the new modeling paradigm in chemical sciences. An important class of multiscale models involves the mapping of a finer scale model into an approximate surface that is used by a coarser scale model. As a specific example of this class we present the case of the adsorption dynamics of diatomic molecules on single crystal catalyst surfaces. The prototype system studied is the dissociative adsorption of H₂ on Pt(111). The finer scale model consists of density functional theory (DFT) periodic slab calculations that provide a small dataset for training an atomistic scale potential energy surface. The coarser scale model uses a semi-classical molecular dynamics (MD) algorithm to obtain the sticking coefficient as a function of the incident energy. Comparison to experimental data and published simulation work is presented. Finally, major challenges in multiscale modeling of chemical reactivity in coupled DFT/MD simulations are discussed, specifically the need for a systematic method of assessing the accuracy of the coarse graining process.     

Effect of temperature,Si/Al and metal loading on the positioning of Pt atoms in MOR-type zeolites

Metropolis Monte Carlo (MMC) and molecular dynamics (MD) simulations were performed to study the feasibility of synthesising Pt nanowires inside the pores of zeolites with the MOR framework, such as mordenite. Results show that the temperature, the silicon to aluminium ratio (Si/Al) of the zeolite and the Pt metal loading have an important combined effect on the positioning of Pt atoms inside the framework. MMC simulations indicate that low Si/Al and high metal loadings promote the positioning of Pt atoms inside the main pore channels, which are the places where metal nanowires could be formed. On the other hand, high Si/Al and low metal loadings favour the positioning of Pt atoms in the side pockets of the MOR structure. For Pt loadings of 5% (mol/mol), the guest metal atoms were positioned mainly in the main pore channels for all conditions studied. MD simulation results are in agreement with MMC simulations, showing that when Si/Al = ∞, the trend of single Pt atoms is to move towards the side pockets of the MOR framework, while for Si/Al = 5, the Pt atoms remain in the main pore channel of the structure. Temperature had the effect of increasing the mobility of Pt atoms through the MOR framework.     

Understanding the dual mechanism of bioactive peptides targeting the enzymes involved in Renin Angiotensin System (RAS): An in-silico approach

Abstract

Understanding the dual inhibition mechanism of food derivative peptides targeting the enzymes (Renin and Angiotensin Converting enzyme) in the Renin Angiotensin System. Two peptides RALP and WYT were reported to possess antihypertensive activity targeting both renin and ACE, and we have used molecular docking and molecular dynamics simulation, in order to understand the underlying mechanism. The selected peptides (RALP and WYT) from the series of peptides reported were docked to renin and ACE and two binding modes were selected based on the binding energy, interaction pattern and clusters of docking simulation. The enzyme-peptide complexes for renin and ACE (Renin/RALP_1,2; ACE/RALP_1,2; Renin/WYT_1,2 and ACE/WYT_1,2) were subjected to molecular dynamics simulation. Our results identified that the peptides inhibiting renin, tends to move out of the binding pockets (S1’ S2’) which is critical for potent binding and occupies the less important pockets (S4 and S3). This could possibly be the reason for its low potency. Whereas, the same peptides targeting ACE, tends to be intact in the pocket because of the metal ion coordination and there is an ample room to improve on its efficacy. Our results further pave way for the biochemist, medicinal chemist to design dual peptides targeting the RAS effectively.

Communicated by Ramaswamy H. Sarma     

A Tight-binding Molecular Dynamics Study of the Dissociation of Boron Clusters in c-Si

Boron clusters with silicon self-interstitials have been implicated in the phenomenon of transient enhanced diffusion (TED) following ion implantation of boron and subsequent annealing steps. This paper explores possible dissolution mechanisms for boron-interstitial clusters during the simulation of a typical annealing process. Using tight-binding molecular dynamics (TBMD) and employing a Goodwin-Skinner-Pettifor sp-based TB model, we have been able to observe the complete dissolution of a B 4 I 4 cluster into the surrounding crystalline silicon matrix. Many unsuccessful attempts to observe dissolution are also presented, highlighting the effect of cluster stability, temperature and the role of vacancies in cluster dissolution. Though we can make no unambiguous statements on the definitive dissolution mechanism of boron-defect clusters based on one successful dissolution event, we can hint at key events that appear to be important, such as the diffusion of self-interstitials (presciently predicted by Pelaz et al. ), the "stranding" of boron atoms in their wake, and the importance of mobile boron-self-interstitial (B-I) pairs. The intrinsic diffusivity of boron in a c-Si lattice and its retardation of the diffusivity of Si self-interstitials is also discussed.     

Molecular modeling of layered double hydroxide intercalated with benzoate,modeling and experiment

The structure of Zn4Al2 Layered Double Hydroxide intercalated with benzencarboxylate (C6H5COO-) was solved using molecular modeling combined with experiment (X-ray powder diffraction, IR spectroscopy, TG measurements). Molecular modeling revealed the arrangement of guest molecules, layer stacking, water content and water location in the interlayer space of the host structure. Molecular modeling using empirical force field was carried out in Cerius(2) modeling environment. Results of modeling were confronted with experiment that means comparing the calculated and measured diffraction pattern and comparing the calculated water content with the thermogravimetric value. Good agreement has been achieved between calculated and measured basal spacing: d(calc) = 15.3 A and d(exp) = 15.5 A. The number of water molecules per formula unit (6H2O per Zn4Al2(OH)12) obtained by modeling (i.e., corresponding to the energy minimum) agrees with the water content estimated by thermogravimetry. The long axis of guest molecules are almost perpendicular to the LDH layers, anchored to the host layers via COO- groups. Mutual orientation of benzoate ring planes in the interlayer space keeps the parquet arrangement. Water molecules are roughly arranged in planes adjacent to host layers together with COO- groups.     

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins—ion channels and transporters—which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.     

Computer Simulation of Surface Segregation

Abstract

Computer simulation techniques can now reliably model the surface structure and energies of inorganic solids. We present recent work which has been directed at modelling the segregation of impurities. The techniques are based on energy minimisation where the forces between the atoms are described by the Born model of solids. We initially show that magnesium and calcium segregation to the prism surfaces of alumina reproduce the available experimental data within the uncertainties of the experimental technique. However, reliable experimental data is not always available, and hence we show using the examples of calcium segregation to zirconia surfaces can give insights to the surface structure and provide predictions for experiment to test. Secondly, the simulations can be verified by comparison with morphologies. This is illustrated by comparison of the experimentally determined morphology of calcite with magnesium and lithium impurities with those calculated. Next extensions to the energy minimisation methods are described using dynamical techniques based on lattice and molecular dynamics. Finally, we describe the situation where there is phase separation and the structure and stability of the interface is governed by the epitaxial relations between the underlying oxide and the impurity oxide.     

Generation of initial molecular dynamics configurations in arbitrary geometries and in parallel

A computational pre-processing tool for generating initial configurations of molecules for molecular dynamics simulations in geometries described by a mesh of unstructured arbitrary polyhedra is described. The mesh is divided into separate zones and each can be filled with a single crystal lattice of atoms. Each zone is filled by creating an expanding cube of crystal unit cells, initiated from an anchor point for the lattice. Each unit cell places the appropriate atoms for the user-specified crystal structure and orientation. The cube expands until the entire zone is filled with the lattice; zones with concave and disconnected volumes may be filled. When the mesh is spatially decomposed into portions for distributed parallel processing, each portion may be filled independently, meaning that the entire molecular system never needs to fit onto a single processor, allowing very large systems to be created. The computational time required to fill a zone with molecules scales linearly with the number of cells in the zone for a fixed number of molecules, and better than linearly with the number of molecules for a fixed number of mesh cells. Our tool, molConfig, has been implemented in the open source C++ code OpenFOAM.     

Molecular Dynamics Study of the Diffusion of Hydrogen in Solid Palladium

Abstract

The diffusion of hydrogen atoms in solid palladium is studied via a molecular dynamics simulation. In this calculation, the palladium atoms are fixed on the sites of the fcc lattice and the hydrogen atoms are initially put on the O (octahedral) sites. Through the present molecular dynamics simulation, the diffusion constant and activation energy are calculated for several different concentrations of hydrogen atoms. We find that the hydrogen atoms show jump motions between the O sites which lead to super diffusion in the solid palladium. We have also obtained the temperature dependence and concentration dependence of the vacancy factor and jump correlation functions.     

Atom ordering in cuboctahedral Ni-Al nanoalloys

Energy calculations have been carried out on high-symmetry cuboctahedral Ni-Al nanoalloy clusters, of varying composition, with the interatomic interactions modelled by the Gupta many-body potential. Relaxations of cuboctahedral fragments cut from the bulk lattice of Ni₃Al, with 13-561 atoms, were undertaken, as were relaxations of high symmetry clusters with 55 and 147 atoms. The lowest energy isomers were found to be dominated by three factors: the tendency toward mixing due to the favourable energy of mixing, Δ_mixE; the size difference between nickel and aluminium; and the higher cohesive and surface energy of nickel compared to aluminium. The latter two factors favour Al-segregation to the surface. The most stable Ni:Al composition approaches 3:1 for larger clusters.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.     

A QM/MM study of the complexes formed by aluminum and iron with serum transferrin at neutral and acidic pH

Serum transferrin (sTf) transports iron in serum and internalizes in cells via receptor mediated endocytosis. Additionally, sTf has been identified as the predominant aluminum carrier in serum. Some questions remain unclear about the exact mechanism for the metal release or whether the aluminum and iron show the same binding mode during the entire process. In the present work, simulation techniques at quantum and atomic levels have been employed in order to gain access into a molecular level understanding of the metal-bound sTf complex, and to describe the binding of Al(III) and Fe(III) ions to sTf. First, hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations were carried out in order to analyze the dynamics of the aluminum-loaded complex, taking into account the different pH conditions in blood and into the cell. Moreover, the complexes formed by transferrin with Al(III) and Fe(III) were optimized with high level density functional theory (DFT)/MM methods. All these results indicate that the interaction mode of Al(III) and Fe(III) with sTf change upon different pH conditions, and that the coordination of Al(III) and Fe(III) is not equivalent during the metal intake, transport and release processes. Our results emphasize the importance of the pH on the metal binding and release mechanism and suggest that Al(III) can follow the iron pathway to get access into cells, although once there, it may show a different binding mode, leading to a different mechanism for its release.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

The thermal properties of binary structure sI clathrate hydrate from molecular dynamics simulation

In this work, we present temperature dependence of lattice parameter and normalised lattice parameter in the atmospheric pressure and 120 bar and also pressure dependence of unit cell volume and normalised unit cell volume at 150 and 250?K for variety guests with different size, polarity and guest–host hydrogen bonding capability such as trimethylene oxide (TMO), ethylene oxide (EO), formaldehyde (FA), cyclobutane (CB), cyclopropane (CP) and ethane (Et) in the large cages with CH₄ in small cages of sI clathrate hydrates by molecular dynamics simulations. The obtained values of lattice parameters for the guest species are compatible with the experimental values. These clathrate hydrates are simulated with TIP4P/ice four-site water potential. Herein, isobaric thermal expansivity and isothermal compressibility are calculated at a temperature range of 50–250?K and a wide pressure range. These structural properties have been compared for guests which they are isoelectronic and have similar masses but with different size and polarity. We use molecular dynamics simulations to relate microscopic guest properties, like guest–host hydrogen bonding to macroscopic sI clathrate hydrate properties. The temperature dependence of thermodynamic properties such as constant-volume and constant-pressure heat capacity is presented in the atmospheric pressure for these guest species.     

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Density functional study on size-dependent structures, stabilities, electronic and magnetic properties of Au n M (M = Al and Si, n?=?1–9) clusters: comparison with pure gold clusters

The density functional theory (DFT) method has been employed to systematically investigate the geometrical structures, relative stabilities, and electronic and magnetic properties of Au(n)M (M = Al and Si, n = 1-9) clusters for clarifying the effect of Al(Si) modulation on the gold nanostructures. Of all the clusters studied, the most stable configurations adopt a three-dimensional structure for Au(n)Al at n = 4-8 and Au(n)Si at n = 3-9, while for pure gold systems, no three-dimensional lowest energy structures are obtained. Through a careful analysis of the fragmentation energy, second-order difference of energy, HOMO-LUMO energy gap, and magnetic moment as a function of cluster size, an odd-even alternative phenomenon has been observed. The results show that the clusters with even-number valence electrons have a higher relative stability, but lower magnetic moments. Furthermore, Al(Si) doping is found to enhance the stabilities of gold frameworks. In addition, the charge analysis has been given to understand the different effects of individual doped atom on electronic properties and compared further.     

Multishelled Gold Nanowires

Abstract

The current miniaturization of electronic devices raises many questions about the properties of various materials at nanometre-scales. Recent molecular dynamics computer simulations have shown that small finite nanowires of gold exist as multishelled structures of lasting stability. These classical simulations are based on a well-tested embedded atom potential. Molecular dynamics simulation studies of metallic nanowires should help in developing methods for their fabrication, such as electron-beam litography and scanning tunneling microscopy.     

Molecular dynamics simulations of void defects in the energetic material HMX

A molecular dynamics (MD) simulation was carried out to characterize the dynamic evolution of void defects in crystalline octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine (HMX). Different models were constructed with the same concentration of vacancies (10 %) to discuss the size effects of void. Energetic ground state properties were determined by annealing simulations. The void formation energy per molecule removed was found to be 55–63 kcal/mol^?1, and the average binding energy per molecule was between 32 and 34 kcal/mol^?1 according to the change in void size. Voids with larger size had lower formation energy. Local binding energies for molecules directly on the void surface decreased greatly compared to those in defect-free lattice, and then gradually increased until the distance away from the void surface was around 10 Å. Analysis of 1 ns MD simulations revealed that the larger the void size, the easier is void collapse. Mean square displacements (MSDs) showed that HMX molecules that had collapsed into void present liquid structure characteristics. Four unique low-energy conformers were found for HMX molecules in void: two whose conformational geometries corresponded closely to those found in HMX polymorphs and two, additional, lower energy conformers that were not seen in the crystalline phases. The ratio of different conformers changed with the simulated temperature, in that the ratio of α conformer increased with the increase in temperature.     

A Molecular Dynamics Computer Simulation Performance Comparison of Java Versus C

Abstract

The relative performance of molecular dynamics (MD) computer simulations of fluids written in ANSI C is compared to that achieved by a comparable program written in Java. The performance of the Java program is shown to be dependent upon its runtime environment. The Java Runtime Environment (JRE) from the Java Development Kit (JDK) 1.2 provides a Just-In-Time (JIT) compiler option on Solaris and Windows 95 platforms which decreases the execution time by approximately 4–10× compared to the standard Java interpreter. The compiled Java implementation of the MD computer simulation runs between 30–100% slower. depending on the platform, compared to the equivalent C implementation. The stability of the two simulations, as measured by conservation of energy is shown to be identical to within ～ 1% over 10⁵ time steps.