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.

Interdigitating organic bilayers direct the short interlayer spacing in hybrid organic–inorganic layered vanadium oxide nanostructures

Layered metal oxides provide a single-step route to sheathed superlattices of atomic layers of a variety of inorganic materials, where the interlayer spacing and overall layered structure forms the most critical feature in the nanomaterials’ growth and application in electronics, health, and energy storage. We use a combination of computer simulations and experiments to describe the atomic-scale structure, dynamics and energetics of alkanethiol-intercalated layered vanadium oxide-based nanostructures. Molecular dynamics (MD) simulations identify the unusual substrate-constrained packing of the alkanethiol surfactant chains along each V2O5 (010) face that combines with extensive interdigitation between chains on opposing faces to maximize three-dimensional packing in the interlayer regions. The findings are supported by high resolution electron microscopy analyses of synthesized alkanethiol-intercalated vanadium oxide nanostructures, and the preference for this new interdigitated model is clarified using a large set of MD simulations. This dependency stresses the importance of organic–inorganic interactions in layered material systems, the control of which is central to technological applications of flexible hybrid nanomaterials.

A first principles analysis of the effect of hydrogen concentration in hydrogenated amorphous silicon on the formation of strained Si-Si bonds and the optical and mobility gaps

In this paper, we use a model of hydrogenated amorphous silicon generated from molecular dynamics with density functional theory calculations to examine how the atomic geometry and the optical and mobility gaps are influenced by mild hydrogen oversaturation. The optical and mobility gaps show a volcano curve as the hydrogen content varies from undersaturation to mild oversaturation, with largest gaps obtained at the saturation hydrogen concentration. At the same time, mid-gap states associated with dangling bonds and strained Si-Si bonds disappear at saturation but reappear at mild oversaturation, which is consistent with the evolution of optical gap. The distribution of Si-Si bond distances provides the key to the change in electronic properties. In the undersaturation regime, the new electronic states in the gap arise from the presence of dangling bonds and strained Si-Si bonds, which are longer than the equilibrium Si-Si distance. Increasing hydrogen concentration up to saturation reduces the strained bonds and removes dangling bonds. In the case of mild oversaturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structure shows that the extra hydrogen atoms form a bridge between neighbouring silicon atoms, thus increasing the Si-Si distance and increasing disorder in the sample.

The (Un)Folding of Multidomain Proteins Through the Lens of Single-molecule Force-spectroscopy and Computer Simulation

Proteins are specialized molecules that catalyze most of the reactions that can sustain life, and they become functional by folding into a specific 3D structure. Despite their importance, the question, "how do proteins fold?" - first pondered in in the 1930's - is still listed as one of the top unanswered scientific questions as of 2005, according to the journal Science. Answering this question would provide a foundation for understanding protein function and would enable improved drug targeting, efficient biofuel production, and stronger biomaterials. Much of what we currently know about protein folding comes from studies on small, single-domain proteins, which may be quite different from the folding of large, multidomain proteins that predominate the proteomes of all organisms.

In this thesis I will discuss my work to fill this gap in understanding by studying the unfolding and refolding of large, multidomain proteins using the powerful combination of single-molecule force-spectroscopy experiments and molecular dynamic simulations.

The three model proteins studied - Luciferase, Protein S, and Streptavidin - lend insight into the inter-domain dependence for unfolding and the subdomain stabilization of binding ligands, and ultimately provide new insight into atomistic details of the intermediate states along the folding pathway.

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.

Interactions of Cytochrome P450 3A4 with Phospholipid Bilayer Nanodiscs

Thesis (Ph.D.)--University of Washington, 2016-08

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

Coarse Graining to Invesitigate Peptide-Lipid Interactions

Experimental characterization of molecular details is challenging, and although single molecule experiments have gained prominence, oligomer characterization remains largely unexplored. The ability to monitor the time evolution of individual molecules while they self assemble is essential in providing mechanistic insights about biological events. Molecular dynamics (MD) simulations can fill the gap in knowledge between single molecule experiments and ensemble studies like NMR, and are increasingly used to gain a better understanding of microscopic properties. Coarse-grained (CG) models aid in both exploring longer length and time scale molecular phenomena, and narrowing down the key interactions responsible for significant system characteristics. Over the past decade, CG techniques have made a significant impact in understanding physicochemical processes. However, the realm of peptide-lipid interfacial interactions, primarily binding, partitioning and folding of amphipathic peptides, remains largely unexplored compared to peptide folding in solution. The main drawback of existing CG models is the inability to capture environmentally sensitive changes in dipolar interactions, which are indigenous to protein folding, and lipid dynamics. We have used the Drude oscillator approach to incorporate structural polarization and dipolar interactions in CG beads to develop a minimalistic peptide model, WEPPROM (Water Explicit Polarizable PROtein Model), and a lipid model WEPMEM (Water Explicit Polarizable MEmbrane Model). The addition of backbone dipolar interactions in a CG model for peptides enabled us to achieve alpha-beta secondary structure content de novo, without any added bias. As a prelude to studying amphipathic peptide-lipid membrane interactions, the balance between hydrophobicity and backbone dipolar interactions in driving ordered peptide aggregation in water and at a hydrophobic-hydrophilic interface, was explored. We found that backbone dipole interactions play a crucial role in driving ordered peptide aggregation, both in water and at hydrophobic-hydrophilic interfaces; while hydrophobicity is more relevant for aggregation in water. A zwitterionic (POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and an anionic lipid (POPS: 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) are used as model lipids for WEPMEM. The addition of head group dipolar interactions in lipids significantly improved structural, dynamic and dielectric properties of the model bilayer. Using WEPMEM and WEPPROM, we studied membrane-induced peptide folding of a cationic antimicrobial peptide with anticancer activity, SVS-1. We found that membrane-induced peptide folding is driven by both (a) cooperativity in peptide self interaction and (b) cooperativity in membrane-peptide interactions. The dipolar interactions between the peptide and the lipid head-groups contribute to stabilizing folded conformations. The role of monovalent ion size and peptide concentration in driving lipid domain formation in anionic/zwitterionic lipid mixtures was also investigated. Our study suggest monovalent ion size to be a crucial determinant of interaction with lipid head groups, and hence domain formation in lipid mixtures. This study reinforces the role of dipole interactions in protein folding, lipid membrane properties, membrane induced peptide folding and lipid domain formation. Therefore, the models developed in this thesis can be used to explore a multitude of biomolecular processes, both at longer time-scales and larger system sizes.

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.

Nonlinear breathing modes at a defect

Recent molecular dynamics (MD) simulations of Cubero et al (1999) of a DNA duplex containing the 'rogue' base difluorotoluene (F) in place of a thymine (T) base show that breathing events can occur on the nanosecond timescale, whereas breathing events in a normal DNA duplex take place on the microsecond timescale. The main aim of this paper is to analyse a nonlinear Klein-Gordon lattice model of the DNA duplex including both nonlinear interactions between opposing bases and a defect in the interaction at one lattice site; each of which can cause localisation of energy. Solutions for a breather mode either side of the defect are derived using multiple-scales asymptotics and are pieced together across the defect to form a solution which includes the effects of the nonlinearity and the defect. We consider defects in the inter-chain interactions and in the along chain interactions. In most cases we find in-phase breather modes and/or out-of-phase breather modes, with one case displaying a shifted mode.

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.

Computational and experimental study of the behavior of cyano-based ionic liquids in aqueous solution

The solvation of cyano- (CN-) based ionic liquids (ILs) and their capacity to establish hydrogen bonds (H-bonds) with water was studied by means of experimental and computational approaches. Experimentally, water activity data were measured for aqueous solutions of ILs based on 1-butyl-3-methylimidazolium ([BMIM](+)) cation combined with one of the following anions: thiocyanate ([SCN](-)), dicyanamide ([DCA](-)), or tricyanomethanide ([TCM](-)), and of 1-ethyl-3-methylimidazolium tetracyanoborate ([EMIM][TCB]). From the latter data, water activity coefficients were estimated showing that [BMIM][SCN] and [BMIM][DCA], unlike [BMIM][TCM] and [EMIM][TCB], are able to establish favorable interactions with water. Computationally, the conductor like screening model for real solvents (COSMO-RS) was used to estimate the water activity coefficients which compare well with the experimental ones. From the COSMO-RS results, it is suggested that the polarity of each ion composing the ILs has a strong effect on the solvation phenomena. Furthermore, classical molecular dynamics (MD) simulations were performed for obtaining an atomic level picture of the local molecular neighborhood of the different species. From the experimental and computational data it is showed that increasing the number of CN groups in the ILs' anions does not enhance their ability to establish H-bonds with water but decreases their polarities, being [BMIM][DCA] and [BMIM][SCN] the ones presenting higher propensity to interact.

Urban dynamics evaluation of Envigado city

This article presents the methodological conceptualization and the main results of a System Dynamics model, which main objective is to support the housing policies in the city of Envigado -- The used methodology developed a scenario-based model to emulate the approximate evolution of the housing demand and supply for the city, using a scenario of the Gross Domestic Product (GDP) and a housing authorization strategy as input -- Diverse results were obtained, for instance it was found that due to the soil availability, the housing supply reaches the saturation point between 2040 and 2046 -- Finally this article could be considered as an example of how academic tools such as System Dynamics can be used by decisions makers in the government

Ação de anestésicos gerais em canais iônicos

Dissertação (mestrado)—Universidade de Brasília, Instituto de Ciências Biológicas, Departamento de Biologia Molecular, 2016.

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.