Molecular dynamics simulations of the adsorption and diffusion behavior of pure and mixed alkanes in silicalite

The adsorption and diffusion of mixed hydrocarbon components in silicalite have been studied using molecular dynamic simulation methods. We have investigated the effect of molecular loadings and temperature on the diffusional behavior of both pure and mixed alkane components. For binary mixtures with components of similar sizes, molecular diffusional behavior in the channels was noticed to be reversed as loading is increased. This behavior was noticeably absent for components of different sizes in the mixture. Methane molecules in the methane/propane mixture have the highest diffusion coefficients across the entire loading range. Binary mixtures containing ethane molecules prove more difficult to separate compared to other binary components. In the ternary mixture, however, ethane molecules diffuse much faster at 400 K in the channel with a tendency to separate out quickly from other components. © 2005 Elsevier Inc. All rights reserved.

Design and characterization of zeolite membranes for fluid separation

Experimental and theoretical methods have been used to study zeolite structures, properties and applications as membranes for separation purposes. Thin layers of silicalite-1 and Na-LTA zeolites have been synthesised onto carbon-graphite supports using a hydrothermal synthesis procedure. The separation behaviour of the composite membranes was characterized by gas permeation studies of pure, binary and ternary mixtures of methane, ethane and propane. The influence of temperature and feed gas mixture composition on the separation and selectivity performance of the membranes was also investigated. It was found that the silicalite-1 composite membranes synthesised onto the 4 hour oxidized carbon-graphite supports showed the most promising separation behaviour of all the composite membranes investigated. Molecular simulation methods were used to gain an understanding of how hydrocarbon molecules behave both within the pores and on the surfaces of silicalite-1, mordenite and LTA zeolites. Molecular dynamic simulations were used to investigate the influence of temperature and molecular loadings on the diffusional behaviour of hydrocarbons in zeolites. Both hydroxylated (surface termination with hydroxyl groups) and non-hydroxylated silicalite-1 and Na-mordenite surfaces were generated. For both zeolites the most stable surfaces correspond to the {010} surface. For the silicalite-1 {010} surface the adsorption of hydrocarbons and molecular water onto the hydroxylated surface showed a favourable exothermic adsorption process compared to adsorption on the non-hydroxylated surface. With the Na-mordenite {010} surface the adsorption of hydrocarbons onto both the hydroxylated and non-hydroxylated surfaces had a combination of favourable and non-favourable adsorption energies, while the adsorption of molecular water onto both types of surface was found to be a favourable adsorption process.

Interdependent lateral interactions, hydrophobicity and acid strength and their influence on the catalytic activity of nanoporous sulfonic acid silicas

A series of propylsulfonic (MCM-SOH) and octyl co-functionalised propylsulfonic (MCM-Oc-SOH) catalysts have been prepared by post modification of MCM-41 with mercaptopropyltrimethoxysilane (MPTS) to achieve SOH surface coverages spanning the range 0.12-1 monolayer. Within the MCM-Oc-SOH series, samples with submonolayer MPTS coverages were further grafted with octyltrimethoxysilane to cap bare hydroxyl sites and tune the hydrophobicity of the support. For the MCM-SO H series NH calorimetry revealed acid strength increases as a function of sulfonic acid loading, with -ΔH(NH ) increasing from 87 to 118 kJ mol. In contrast, MCM-Oc-SOH exhibits a dramatic enhancement of acid strength for submonolayer SOH coverages, with -ΔH(NH ) found to increase to 103 kJ mol. In line with these acid strength measurements the per-site activity of the MCM-SOH series in the esterification of butanol with acetic acid was found to increase with SOH content. Incorporation of octyl groups further promotes esterification activity of all the samples within the MCM-Oc-SOH series, such that the turn over frequency of the sample with the lowest loading of SOH more than doubles. Molecular dynamic simulations indicate that the interaction of isolated sulfonic acid groups with the pore walls is the primary cause of the decrease in acid strength and activity of submonolayer samples within the MCM-SOH series. Incorporation of octyl groups results in a combination of increased hydrophobicity and lateral interactions between adjacent sulfonic acid head groups, resulting in a striking enhancement of acid strength and esterification activity. © 2010 The Royal Society of Chemistry.

Development of potent and selective tissue transglutaminase inhibitors:their effect on TG2 function and application in pathological conditions

Potent-selective peptidomimetic inhibitors of tissue transglutaminase (TG2) were developed through a combination of protein-ligand docking and molecular dynamic techniques. Derivatives of these inhibitors were made with the aim of specific TG2 targeting to the intra- and extracellular space. A cell-permeable fluorescently labeled derivative enabled detection of in situ cellular TG2 activity in human umbilical cord endothelial cells and TG2-transduced NIH3T3 cells, which could be enhanced by treatment of cells with ionomycin. Reaction of TG2 with this fluorescent inhibitor in NIH3T3 cells resulted in loss of binding of TG2 to cell surface syndecan-4 and inhibition of translocation of the enzyme into the extracellular matrix, with a parallel reduction in fibronectin deposition. In human umbilical cord endothelial cells, this same fluorescent inhibitor also demonstrated a reduction in fibronectin deposition, cell motility, and cord formation in Matrigel. Use of the same inhibitor in a mouse model of hypertensive nephrosclerosis showed over a 40% reduction in collagen deposition.

S100A4 downregulates filopodia formation through increased dynamic instability

Cell migration requires the initial formation of cell protrusions, lamellipodia and/or filopodia, the attachment of the leading lamella to extracellular cues and the formation and efficient recycling of focal contacts at the leading edge. The small calcium binding EF-hand protein S100A4 has been shown to promote cell motility but the direct molecular mechanisms responsible remain to be elucidated. In this work, we provide new evidences indicating that elevated levels of S100A4 affect the stability of filopodia and prevent the maturation of focal complexes. Increasing the levels of S100A4 in a rat mammary benign tumor derived cell line results in acquired cellular migration on the wound healing scratch assay. At the cellular levels, we found that high levels of S100A4 induce the formation of many nascent filopodia, but that only a very small and limited number of those can stably adhere and mature, as opposed to control cells, which generate fewer protrusions but are able to maintain these into more mature projections. This observation was paralleled by the fact that S100A4 overexpressing cells were unable to establish stable focal adhesions. Using different truncated forms of the S100A4 proteins that are unable to bind to myosin IIA, our data suggests that this newly identified functions of S100A4 is myosin-dependent, providing new understanding on the regulatory functions of S100A4 on cellular migration.

Liquid-phase adsorption of n-paraffins on type 5A molecular sieves

As a basis for the commercial separation of normal paraffins a detailed study has been made of factors affecting the adsorption of binary liquid mixtures of high molecular weight normal paraffins (C12, C16, and C20) from isooctane on type 5A molecular sieves. The literature relating to molecular sieve properties and applications, and to liquid-phase adsorption of high molecular weight normal paraffin compounds by zeolites, was reviewed. Equilibrium isotherms were determined experimentally for the normal paraffins under investigation at temperatures of 303oK, 323oK and 343oK and showed a non-linear, favourable- type of isotherm. A higher equilibrium amount was adsorbed with lower molecular weight normal paraffins. An increase in adsorption temperature resulted in a decrease in the adsorption value. Kinetics of adsorption were investigated for the three normal paraffins at different temperatures. The effective diffusivity and the rate of adsorption of each normal paraffin increased with an increase in temperature in the range 303 to 343oK. The value of activation energy was between 2 and 4 kcal/mole. The dynamic properties of the three systems were investigated over a range of operating conditions (i.e. temperature, flow rate, feed concentration, and molecular sieve size in the range 0.032 x 10-3 to 2 x 10-3m) with a packed column. The heights of adsorption zones calculated by two independent equations (one based on a constant width, constant velocity and adsorption zone and the second on a solute material balance within the adsorption zone) agreed within 3% which confirmed the validity of using the mass transfer zone concept to provide a simple design procedure for the systems under study. The dynamic capacity of type 5A sieves for n-eicosane was lower than for n-hexadecane and n-dodecane corresponding to a lower equilibrium loading capacity and lower overall mass transfer coefficient. The values of individual external, internal, theoretical and experimental overall mass transfer coefficient were determined. The internal resistance was in all cases rate-controlling. A mathematical model for the prediction of dynamic breakthrough curves was developed analytically and solved from the equilibrium isotherm and the mass transfer rate equation. The experimental breakthrough curves were tested against both the proposed model and a graphical method developed by Treybal. The model produced the best fit with mean relative percent deviations of 26, 22, and 13% for the n-dodecane, n-hexadecane, and n-eicosane systems respectively.

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Complexity and chirality indices for molecular informatics::differential geometry approach

Novel molecular complexity measures are designed based on the quantum molecular kinematics. The Hamiltonian matrix constructed in a quasi-topological approximation describes the temporal evolution of the modelled electronic system and determined the time derivatives for the dynamic quantities. This allows to define the average quantum kinematic characteristics closely related to the curvatures of the electron paths, particularly, the torsion reflecting the chirality of the dynamic system. A special attention has been given to the computational scheme for this chirality measure. The calculations on realistic molecular systems demonstrate reasonable behaviour of the proposed molecular complexity indices.

Dynamic complexity of chaotic transitions in high-dimensional classical dynamics:Leu-Enkephalin folding

Leu-Enkephalin in explicit water is simulated using classical molecular dynamics. A ß-turn transition is investigated by calculating the topological complexity (in the "computational mechanics" framework [J. P. Crutchfield and K. Young, Phys. Rev. Lett., 63, 105 (1989)]) of the dynamics of both the peptide and the neighbouring water molecules. The complexity of the atomic trajectories of the (relatively short) simulations used in this study reflect the degree of phase space mixing in the system. It is demonstrated that the dynamic complexity of the hydrogen atoms of the peptide and almost all of the hydrogens of the neighbouring waters exhibit a minimum precisely at the moment of the ß-turn transition. This indicates the appearance of simplified periodic patterns in the atomic motion, which could correspond to high-dimensional tori in the phase space. It is hypothesized that this behaviour is the manifestation of the effect described in the approach to molecular transitions by Komatsuzaki and Berry [T. Komatsuzaki and R.S. Berry, Adv. Chem. Phys., 123, 79 (2002)], where a "quasi-regular" dynamics at the transition is suggested. Therefore, for the first time, the less chaotic character of the folding transition in a realistic molecular system is demonstrated. © Springer-Verlag Berlin Heidelberg 2006.

Modular construction of dynamic nucleodendrimers

Isoguanosine-containing dendritic small molecules self-assemble into decameric nucleodendrimers as observed by 1D NMR spectroscopy, 2D DOSY, and mass spectrometry. In particular, apolar building blocks readily form pentameric structures in acetonitrile while the presence of alkali metals promotes the formation of stable decameric assemblies with a preference for cesium ions. Remarkably, co-incubation of guanosine and isoguanosine-containing nucleodendrons results in the formation of decameric structures in absence of added salts. Further analysis of the mixture indicated that guanosine derivatives facilitate the formation, but are not involved in decameric structures; a process reminiscent of molecular crowding. This molecular system provides a powerful canvas for the rapid and modular assembly of polyfunctional dendritic macromolecules. © 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Complexity of classical dynamics of molecular systems. I. Methodology

Methods for the calculation of complexity have been investigated as a possible alternative for the analysis of the dynamics of molecular systems. “Computational mechanics” is the approach chosen to describe emergent behavior in molecular systems that evolve in time. A novel algorithm has been developed for symbolization of a continuous physical trajectory of a dynamic system. A method for calculating statistical complexity has been implemented and tested on representative systems. It is shown that the computational mechanics approach is suitable for analyzing the dynamic complexity of molecular systems and offers new insight into the process.

Molecular dynamics simulation of brittle fracture in silicon

The fracture process involves converting potential energy from a strained body into surface energy, thermal energy, and the energy needed to create lattice defects. In dynamic fracture, energy is also initially converted into kinetic energy. This paper uses molecular dynamics (MD) to simulate brittle frcture in silicon and determine how energy is converted from potential energy (strain energy) into other forms.

Molecular Dynamics implementation of BN2D or 'Mercedes Benz' water model

Two-dimensional 'Mercedes Benz' (MB) or BN2D water model (Naim, 1971) is implemented in Molecular Dynamics. It is known that the MB model can capture abnormal properties of real water (high heat capacity, minima of pressure and isothermal compressibility, negative thermal expansion coefficient) (Silverstein et al., 1998). In this work formulas for calculating the thermodynamic, structural and dynamic properties in microcanonical (NVE) and isothermal-isobaric (NPT) ensembles for the model from Molecular Dynamics simulation are derived and verified against known Monte Carlo results. The convergence of the thermodynamic properties and the system's numerical stability are investigated. The results qualitatively reproduce the peculiarities of real water making the model a visually convenient tool that also requires less computational resources, thus allowing simulations of large (hydrodynamic scale) molecular systems. We provide the open source code written in C/C++ for the BN2D water model implementation using Molecular Dynamics.

Light scattering study of human serum albumin in pre-denaturation:relation to dynamic transition in water at 42°C

Protein functional motions are ultimately connected to water dynamics. The goal of this study is to link the conformational dynamics of albumin to a dynamic transition taking place at ∼ 42°C in water. We report the results of dynamic light scattering measurements of albumin aqueous solution in the temperature interval 20-65°C. The processing of the experimental data produced the temperature dependence of the macromolecular hydrodynamic radius. We demonstrate that the growth of the macromolecular size in this temperature range can be divided into two stages that are connected to the dynamical properties of water. © 2012 Elsevier B.V. All rights reserved.

Behavioral and molecular genetics of reading-related AM and FM detection thresholds

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h2 = 0.20) and FM (h2 = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.