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.

Non-Markov state model of peptide dynamics

A hidden Markov state model has been applied to classical molecular dynamics simulated small peptide in explicit water. The methodology allows increasing the time resolution of the model and describe the dynamics with the precision of 0.3 ps (comparing to 6 ps for the standard methodology). It also permits the investigation of the mechanisms of transitions between the conformational states of the peptide. The detailed description of one of such transitions for the studied molecule is presented. © 2012 Elsevier B.V. All rights reserved.

Molecular mechanisms of salt effects on carbon nanotube dispersions in an organic solvent (N-methyl-2-pyrrolidone)

We consider the effects of salt (sodium iodide) on pristine carbon nanotube (CNT) dispersions in an organic solvent, N-methyl-2-pyrrolidone (NMP). We investigate the molecular-scale mechanisms of ion interactions with the nanotube surface and we show how the microscopic ion-surface interactions affect the stability of CNT dispersions in NMP. In our study we use a combination of fully atomistic Molecular Dynamics simulations of sodium and iodide ions at the CNT-NMP interface with direct experiments on the CNT dispersions. In the experiments we analyze the effects of salt on the stability of the dispersions by photoluminescence (PL) and optical absorption spectroscopy of the samples as well as by visual inspection. By fully atomistic Molecular Dynamics simulations we investigate the molecular-scale mechanisms of sodium and iodide ion interactions with the nanotube surface. Our simulations reveal that both ions are depleted from the CNT surface in the CNT-NMP dispersions mainly due to the two reasons: (1) there is a high energy penalty for the ion partial desolvation at the CNT surface; (2) NMP molecules form a dense solvation layer at the CNT surface that prevents ions to come close to the CNT surface. As a result, an increase of the salt concentration increases the "osmotic" stress in the CNT-NMP system and, thus, decreases the stability of the CNT dispersions in NMP. Direct experiments confirm the simulation results: addition of NaI salt into the NMP dispersions of pristine CNTs leads to precipitation of CNTs (bundle formation) even at very small salt concentration (∼10 -3 mol L -1). In line with the simulation predictions, the effect increases with the increase of the salt concentration. Overall, our results show that dissolved salt ions have strong effects on the stability of CNT dispersions. Therefore, it is possible to stimulate the bundle formation in the CNT-NMP dispersions and regulate the overall concentration of nanotubes in the dispersions by changing the NaI concentration in the solvent. © 2012 The Royal Society of Chemistry.

Multiscale modelling:approaches and challenges

Multiscale systems that are characterized by a great range of spatial–temporal scales arise widely in many scientific domains. These range from the study of protein conformational dynamics to multiphase processes in, for example, granular media or haemodynamics, and from nuclear reactor physics to astrophysics. Despite the diversity in subject areas and terminology, there are many common challenges in multiscale modelling, including validation and design of tools for programming and executing multiscale simulations. This Theme Issue seeks to establish common frameworks for theoretical modelling, computing and validation, and to help practical applications to benefit from the modelling results. This Theme Issue has been inspired by discussions held during two recent workshops in 2013: ‘Multiscale modelling and simulation’ at the Lorentz Center, Leiden (http://www.lorentzcenter.nl/lc/web/2013/569/info.php3?wsid=569&venue=Snellius), and ‘Multiscale systems: linking quantum chemistry, molecular dynamics and microfluidic hydrodynamics’ at the Royal Society Kavli Centre. The objective of both meetings was to identify common approaches for dealing with multiscale problems across different applications in fluid and soft matter systems. This was achieved by bringing together experts from several diverse communities.

Molecular alignment relaxation in polymer optical fibers for sensing applications

A systematic study of annealing behavior of drawn PMMA fibers was performed. Annealing dynamics were investigated under different environmental conditions by fiber longitudinal shrinkage monitoring. The shrinkage process was found to follow a stretched exponential decay function revealing the heterogeneous nature of the underlying molecular dynamics. The complex dependence of the fiber shrinkage on initial degree of molecular alignment in the fiber, annealing time and temperature was investigated and interpreted. Moreover, humidity was shown to have a profound effect on the annealing process, which was not recognized previously. Annealing was also shown to have considerable effect on the fiber mechanical properties associated with the relaxation of molecular alignment in the fiber. The consequences of fiber annealing for the climatic stability of certain polymer optical fiber-based sensors are discussed, emphasizing the importance of fiber controlled pre-annealing with respect to the foreseeable operating conditions.

Theoretical study of chloroperoxidase catalyzed chlorination of beta-cyclopentanedione and role of water in the chlorination mechanism

Chloroperoxidase (CPO) is a potential biocatalyst for use in asymmetric synthesis. The mechanisms of CPO catalysis are therefore of interest. The halogenation reaction, one of several chemical reactions that CPO catalyzes, is not fully understood and is the subject of this dissertation. The mechanism by which CPO catalyzes halogenation is disputed. It has been postulated that halogenation of substrates occurs at the active site. Alternatively, it has been proposed that hypochlorous acid, produced at the active site via oxidation of chloride, is released prior to reaction, so that halogenation occurs in solution. The free-solution mechanism is supported by the observation that halogenation of most substrates often occurs non-stereospecifically. On the other hand, the enzyme-bound mechanism is supported by the observation that some large substrates undergo halogenation stereospecifically. The major purpose of this research is to compare chlorination of the substrate β-cyclopentanedione in the two environments. One study was of the reaction with limited hydration because such a level of hydration is typical of the active site. For this work, a purely quantum mechanical approach was used. To model the aqueous environment, the limited hydration environment approach is not appropriate. Instead, reaction precursor conformations were obtained from a solvated molecular dynamics simulation, and reaction of potentially reactive molecular encounters was modeled with a hybrid quantum mechanical/molecular mechanical approach. Extensive work developing parameters for small molecules was pre-requisite for the molecular dynamics simulation. It is observed that a limited and optimized (active-site-like) hydration environment leads to a lower energetic barrier than the fully solvated model representative of the aqueous environment at room temperature, suggesting that the stable water network near the active site is likely to facilitate the chlorination mechanism. The influence of the solvent environment on the reaction barrier is critical. It is observed that stabilization of the catalytic water by other solvent molecules lowers the barrier for keto-enol tautomerization. Placement of water molecules is more important than the number of water molecules in such studies. The fully-solvated model demonstrates that reaction proceeds when the instantaneous dynamical water environment is close to optimal for stabilizing the transition state.

Stable algorithm for event detection in event-driven particle dynamics : logical states

Acknowledgments The authors gratefully acknowledge the support of the German Research Foundation (DFG) through the Cluster of Excellence ‘Engineering of Advanced Materials’ at the University of Erlangen-Nuremberg and through Grant Po 472/25.

An allosteric role for receptor activity-modifying proteins in defining GPCR pharmacology

G protein-coupled receptors are allosteric proteins that control transmission of external signals to regulate cellular response. Although agonist binding promotes canonical G protein signalling transmitted through conformational changes, G protein-coupled receptors also interact with other proteins. These include other G protein-coupled receptors, other receptors and channels, regulatory proteins and receptor-modifying proteins, notably receptor activity-modifying proteins (RAMPs). RAMPs have at least 11 G protein-coupled receptor partners, including many class B G protein-coupled receptors. Prototypic is the calcitonin receptor, with altered ligand specificity when co-expressed with RAMPs. To gain molecular insight into the consequences of this protein–protein interaction, we combined molecular modelling with mutagenesis of the calcitonin receptor extracellular domain, assessed in ligand binding and functional assays. Although some calcitonin receptor residues are universally important for peptide interactions (calcitonin, amylin and calcitonin gene-related peptide) in calcitonin receptor alone or with receptor activity-modifying protein, others have RAMP-dependent effects, whereby mutations decreased amylin/calcitonin gene-related peptide potency substantially only when RAMP was present. Remarkably, the key residues were completely conserved between calcitonin receptor and AMY receptors, and between subtypes of AMY receptor that have different ligand preferences. Mutations at the interface between calcitonin receptor and RAMP affected ligand pharmacology in a RAMP-dependent manner, suggesting that RAMP may allosterically influence the calcitonin receptor conformation. Supporting this, molecular dynamics simulations suggested that the calcitonin receptor extracellular N-terminal domain is more flexible in the presence of receptor activity-modifying protein 1. Thus, RAMPs may act in an allosteric manner to generate a spectrum of unique calcitonin receptor conformational states, explaining the pharmacological preferences of calcitonin receptor-RAMP complexes. This provides novel insight into our understanding of G protein-coupled receptor-protein interaction that is likely broadly applicable for this receptor class.

DEVELOPMENT OF A CONDENSED PHASE VELOCITY MAP IMAGING APPARATUS: DISSOCIATION PHOTODYNAMICS OF NITROGEN DIOXIDE AT 355 NM

The photochemistry of the polar regions of Earth, as well as the interstellar medium, is driven by the effect of ultraviolet radiation on ice surfaces and on the materials trapped within them. While the area of ice photochemistry is vast and much research has been completed, it has only recently been possible to study the dynamics of these processes on a microscopic level. One of the leading techniques for studying photoreaction dynamics is Velocity Map Imaging (VMI). This technique has been used extensively to study several types of reaction dynamics processes. Although the majority of these studies have utilized molecular beams as the main medium for reactants, new studies showed the versatility of the technique when applied to molecular dynamics of molecules adsorbed on metal surfaces. Herein the development of a velocity map imaging apparatus capable of studying the photochemistry of condensed phase materials is described. The apparatus is used to study of the photo-reactivity of NO2 condensed within argon matrices to illustrate its capabilities. A doped ice surface is formed by condensing Ar and NO2 gas onto a sapphire rod which is cooled using a helium compressor to 20 K. The matrix is irradiated using an Nd:YAG laser at 355 nm, and the resulting NO fragment is state-selectively ionized using an excimer-pumped dye laser. In all, we are able to detect transient photochemically generated species and can collect information on their quantum state and kinetic energy distribution. It is found that the REMPI spectra changes as different sections of the dissociating cloud are probed. The rotational and translational energy populations are found to be bimodal with a low temperature component roughly at the temperature of the matrix, and a second component with much higher temperature, the rotational temperature showing a possible population inversion, and the translational temperature of 100-200 K. The low temperature translational component is found to dominate at long delay times between dissociation and ionization, while at short time delays the high temperature component plays a larger role. The velocity map imaging technique allows for the detection of both the axial and radial components of the translational energy. The distribution of excess energy over the rotational, electronic and translational states of the NO photofragments provides evidence for collisional quenching of the fragments in the Ar-matrix prior to their desorption.

A study of electrochemical transport and diffuse charge dynamics using a Langevin equation

This thesis aims to develop new numerical and computational tools to study electrochemical transport and diffuse charge dynamics at small scales. Previous efforts at modeling electrokinetic phenomena at scales where the noncontinuum effects become significant have included continuum models based on the Poisson-Nernst-Planck equations and atomic simulations using molecular dynamics algorithms. Neither of them is easy to use or conducive to electrokinetic transport modeling in strong confinement or over long time scales. This work introduces a new approach based on a Langevin equation for diffuse charge dynamics in nanofluidic devices, which incorporates features from both continuum and atomistic methods. The model is then extended to include steric effects resulting from finite ion size, and applied to the phenomenon of double layer charging in a symmetric binary electrolyte between parallel-plate blocking electrodes, between which a voltage is applied. Finally, the results of this approach are compared to those of the continuum model based on the Poisson-Nernst-Planck equations.

A influência da incorporação de ácido α-eleosteárico, purificado e isolado do óleo de tungue, na dinâmica molecular de lipossomos

Este trabalho descreve o isolamento e purificação do ácido α-eleosteárico (α-ESA) a partir do óleo de tungue e sua caracterização por espectroscopia de infravermelho com transformada de Fourier (FTIR), cromatografia gasosa acoplada com espectrometria de massas (GC-MS) e espectroscopia de ressonância magnética nuclear (RMN) de 1H e 13C. O α-ESA apresenta atividades biológicas (antitumorais, anti-inflamatórias e antioxidantes), tornando-se importante compreender sua interação com membranas lipídicas. Assim, este trabalho também descreve resultados referentes ao efeito da incorporação de α-ESA na dinâmica molecular de lipossomos compostos por fosfatidilcolina. O sistema lipossomal puro e contendo α-ESA foi caracterizado através do uso de espectroscopia de UV-visível, FTIR, RMN e calorimetria de varredura diferencial (DSC). Como resultados da purificação do α-ESA, obtivemos uma pureza de 95,9% utilizando acetona como solvente de recristalização em detrimento dos 92,2% em solução etanólica. Na incorporação em lipossomos, observou-se uma maior interação do α-ESA com a parte polar, de interface e os primeiros metilenos da região apolar da fosfatidilcolina. Além disso, α-ESA apresentou um efeito de redução da fluidez de lipossomos. Os resultados contribuem para a geração de conhecimento para o desenvolvimento de novos sistemas farmacológicos.

Host-guest system of 4-nerolidylcatechol in 2-hydroxypropyl-beta-cyclodextrin: preparation, characterization and molecular modeling

The interaction of 4-nerolidylcatechol (4-NRC), a potent antioxidant agent, and 2-hydroxypropyl-beta-cyclodextrin (HP-beta-CD) was investigated by the solubility method using Fourier transform infrared (FTIR) methods in addition to UV-Vis, (1)H-nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. The inclusion complexes were prepared using grinding, kneading and freeze-drying methods. According to phase solubility studies in water a B(S)-type diagram was found, displaying a stoichiometry complexation of 2:1 (drug:host) and stability constant of 6494 +/- A 837 M(-1). Stoichiometry was established by the UV spectrophotometer using Job's plot method and, also confirmed by molecular modeling. Data from (1)H-NMR, and FTIR, experiments also provided formation evidence of an inclusion complex between 4-NRC and HP-beta-CD. 4-NRC complexation indeed led to higher drug solubility and stability which could probably be useful to improve its biological properties and make it available to oral administration and topical formulations.

A THEORETICAL STUDY OF INTERACTION OF NANOPARTICLES WITH BIOMOLECULE

Many types of materials at nanoscale are currently being used in everyday life. The production and use of such products based on engineered nanomaterials have raised concerns of the possible risks and hazards associated with these nanomaterials. In order to evaluate and gain a better understanding of their effects on living organisms, we have performed first-principles quantum mechanical calculations and molecular dynamics simulations. Specifically, we will investigate the interaction of nanomaterials including semiconducting quantum dots and metallic nanoparticles with various biological molecules, such as dopamine, DNA nucleobases and lipid membranes. Firstly, interactions of semiconducting CdSe/CdS quantum dots (QDs) with the dopamine and the DNA nucleobase molecules are investigated using similar quantum mechanical approach to the one used for the metallic nanoparticles. A variety of interaction sites are explored. Our results show that small-sized Cd4Se4 and Cd4S4 QDs interact strongly with the DNA nucleobase if a DNA nucleobase has the amide or hydroxyl chemical group. These results indicate that these QDs are suitable for detecting subcellular structures, as also reported by experiments. The next two chapters describe a preparation required for the simulation of nanoparticles interacting with membranes leading to accurate structure models for the membranes. We develop a method for the molecular crystalline structure prediction of 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC), 1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE) and cyclic di-amino acid peptide using first-principles methods. Since an accurate determination of the structure of an organic crystal is usually an extremely difficult task due to availability of the large number of its conformers, we propose a new computational scheme by applying knowledge of symmetry, structural chemistry and chemical bonding to reduce the sampling size of the conformation space. The interaction of metal nanoparticles with cell membranes is finally carried out by molecular dynamics simulations, and the results are reported in the last chapter. A new force field is developed which accurately describes the interaction forces between the clusters representing small-sized metal nanoparticles and the lipid bilayer molecules. The permeation of nanoparticles into the cell membrane is analyzed together with the RMSD values of the membrane modeled by a lipid bilayer. The simulation results suggest that the AgNPs could cause the same amount of deformation as the AuNPs for the dysfunction of the membrane.

MOLECULAR MODELING OF EPON 862-DETDA / CARBON COMPOSITES

The thermoset epoxy resin EPON 862, coupled with the DETDA hardening agent, are utilized as the polymer matrix component in many graphite (carbon fiber) composites. Because it is difficult to experimentally characterize the interfacial region, computational molecular modeling is a necessary tool for understanding the influence of the interfacial molecular structure on bulk-level material properties. The purpose of this research is to investigate the many possible variables that may influence the interfacial structure and the effect they will have on the mechanical behavior of the bulk level composite. Molecular models are established for EPON 862-DETDA polymer in the presence of a graphite surface. Material characteristics such as polymer mass-density, residual stresses, and molecular potential energy are investigated near the polymer/fiber interface. Because the exact degree of crosslinking in these thermoset systems is not known, many different crosslink densities (degrees of curing) are investigated. It is determined that a region exists near the carbon fiber surface in which the polymer mass density is different than that of the bulk mass density. These surface effects extend ~10 Å into the polymer from the center of the outermost graphite layer. Early simulations predict polymer residual stress levels to be higher near the graphite surface. It is also seen that the molecular potential energy in polymer atoms decreases with increasing crosslink density. New models are then established in order to investigate the interface between EPON 862-DETDA polymer and graphene nanoplatelets (GNPs) of various atomic thicknesses. Mechanical properties are extracted from the models using Molecular Dynamics techniques. These properties are then implemented into micromechanics software that utilizes the generalized method of cells to create representations of macro-scale composites. Micromechanics models are created representing GNP doped epoxy with varying number of graphene layers and interfacial polymer crosslink densities. The initial micromechanics results for the GNP doped epoxy are then taken to represent the matrix component and are re-run through the micromechanics software with the addition of a carbon fiber to simulate a GNP doped epoxy/carbon fiber composite. Micromechanics results agree well with experimental data, and indicate GNPs of 1 to 2 atomic layers to be highly favorable. The effect of oxygen bonded to the surface of the GNPs is lastly investigated. Molecular Models are created for systems with varying graphene atomic thickness, along with different amounts of oxygen species attached to them. Models are created for graphene containing hydroxyl groups only, epoxide groups only, and a combination of epoxide and hydroxyl groups. Results show models of oxidized graphene to decrease in both tensile and shear modulus. Attaching only epoxide groups gives the best results for mechanical properties, though pristine graphene is still favored.