Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Molecular dynamics simulation insight into the mechanism of phenol adsorption at low coverages from aqueous solutions on microporous carbons

MD simulation studies showing the influence of porosity and carbon surface oxidation on phenol adsorption from aqueous solutions on carbons are reported. Based on a realistic model of activated carbon, three carbon structures with gradually changed microporosity were created. Next, a different number of surface oxygen groups was introduced. The pores with diameters around 0.6 nm are optimal for phenol adsorption and after the introduction of surface oxygen functionalities, adsorption of phenol decreases (in accordance with experimental data) for all studied models. This decrease is caused by a pore blocking effect due to the saturation of surface oxygen groups by highly hydrogen-bounded water molecules.

Scale bridging concepts in molecular simulation : coarse-graining and thermodynamic coupling

In dieser Arbeit wurden Simulation von Flüssigkeiten auf molekularer Ebene durchgeführt, wobei unterschiedliche Multi-Skalen Techniken verwendet wurden. Diese erlauben eine effektive Beschreibung der Flüssigkeit, die weniger Rechenzeit im Computer benötigt und somit Phänomene auf längeren Zeit- und Längenskalen beschreiben kann.rnrnEin wesentlicher Aspekt ist dabei ein vereinfachtes (“coarse-grained”) Modell, welches in einem systematischen Verfahren aus Simulationen des detaillierten Modells gewonnen wird. Dabei werden ausgewählte Eigenschaften des detaillierten Modells (z.B. Paar-Korrelationsfunktion, Druck, etc) reproduziert.rnrnEs wurden Algorithmen untersucht, die eine gleichzeitige Kopplung von detaillierten und vereinfachten Modell erlauben (“Adaptive Resolution Scheme”, AdResS). Dabei wird das detaillierte Modell in einem vordefinierten Teilvolumen der Flüssigkeit (z.B. nahe einer Oberfläche) verwendet, während der Rest mithilfe des vereinfachten Modells beschrieben wird.rnrnHierzu wurde eine Methode (“Thermodynamische Kraft”) entwickelt um die Kopplung auch dann zu ermöglichen, wenn die Modelle in verschiedenen thermodynamischen Zuständen befinden. Zudem wurde ein neuartiger Algorithmus der Kopplung beschrieben (H-AdResS) der die Kopplung mittels einer Hamilton-Funktion beschreibt. In diesem Algorithmus ist eine zur Thermodynamischen Kraft analoge Korrektur mit weniger Rechenaufwand möglich.rnrnAls Anwendung dieser grundlegenden Techniken wurden Pfadintegral Molekulardynamik (MD) Simulationen von Wasser untersucht. Mithilfe dieser Methode ist es möglich, quantenmechanische Effekte der Kerne (Delokalisation, Nullpunktsenergie) in die Simulation einzubeziehen. Hierbei wurde zuerst eine Multi-Skalen Technik (“Force-matching”) verwendet um eine effektive Wechselwirkung aus einer detaillierten Simulation auf Basis der Dichtefunktionaltheorie zu extrahieren. Die Pfadintegral MD Simulation verbessert die Beschreibung der intra-molekularen Struktur im Vergleich mit experimentellen Daten. Das Modell eignet sich auch zur gleichzeitigen Kopplung in einer Simulation, wobei ein Wassermolekül (beschrieben durch 48 Punktteilchen im Pfadintegral-MD Modell) mit einem vereinfachten Modell (ein Punktteilchen) gekoppelt wird. Auf diese Weise konnte eine Wasser-Vakuum Grenzfläche simuliert werden, wobei nur die Oberfläche im Pfadintegral Modell und der Rest im vereinfachten Modell beschrieben wird.

Molecular dynamics simulation of the structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays

The structural and dynamic properties of dioctadecyldimethylammoniums (DODDMA) intercalated into 2:1 layered clays are investigated using isothermal-isobaric (NPT) molecular dynamics (MD) simulation. The simulated results are in reasonably good agreement with the available experimental measurements, such as X-ray diffraction (XRD), atom force microscopy (AFM), Fourier transform infrared (FTIR), and nuclear magnetic resonance (NMR) spectroscopies. The nitrogen atoms are found to be located mainly within two layers close to the clay surface whereas methylene groups form a pseudoquadrilayer structure. The results of tilt angle and order parameter show that interior two-bond segments of alkyl chains prefer an arrangement parallel to the clay surface, whereas the segments toward end groups adopt a random orientation. In addition, the alkyl chains within the layer structure lie almost parallel to the clay surface whereas those out of the layer structure are essentially perpendicular to the surface. The trans conformations are predominant in all cases although extensive gauche conformations are observed, which is in agreement with previous simulations on n-butane. Moreover, an odd-even effect in conformation distributions is observed mainly along the chains close to the head and tail groups. The diffusion constants of both nitrogen atoms and methylene groups in these nanoconfined alkyl chains increase with the temperature and methelene position toward the tail groups.

The effect of pH on PAMAM dendrimer-siRNA complexation:endosomal considerations as determined by molecular dynamics simulation

Intracellular degradation of genes, most notably within the endo-lysosomal compartment is considered a significant barrier to (non-viral) gene delivery in vivo. Previous reports based on in vitro studies claim that carriers possessing a mixture of primary, secondary and tertiary amines are able to buffer the acidic environment within the endosome, allowing for timely release of their contents, leading to higher transfection rates. In this report, we adopt an atomistic molecular dynamics (MD) simulation approach, comparing the complexation of 21-bp siRNA with low-generation polyamidoamine (PAMAM) dendrimers (G0 and G1) at both neutral and acidic pHs, the latter of which mimics the degradative environment within maturing 'late-endosomes'. Our simulations reveal that the time taken for the dendrimer-gene complex (dendriplex) to reach equilibrium is appreciably longer at low pH and this is accompanied by more compact packaging of the dendriplex, as compared to simulations performed at neutral pH. We also note larger absolute values of calculated binding free energies of the dendriplex at low pH, indicating a higher dendrimer-nucleic acid affinity in comparison with neutral pH. These novel simulations provide a more detailed understanding of low molecular-weight polymer-siRNA behavior, mimicking the endosomal environment and provide input of direct relevance to the "proton sponge theory", thereby advancing the rational design of non-viral gene delivery systems.

Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

Biochemical and structural characterization of β-lactamases tem-180 and tem-201

With the constant development of new antibiotics, selective pressure is a force to reckon when investigating antibiotic resistance. Although advantageous for medical treatments, it leads to increasing resistance. It is essential to use more potent and toxic antibiotics. Enzymes capable of hydrolyzing antibiotics are among the most common ways of resistance and TEM variants have been detected in several resistant isolates. Due to the rapid evolution of these variants, complex phenotypes have emerged and the need to understand their biological activity becomes crucial. To investigate the biochemical properties of TEM-180 and TEM-201 several computational methodologies have been used, allowing the comprehension of their structure and catalytic activity, which translates into their biological phenotype. In this work we intent to characterize the interface between these proteins and the several antibiotics used as ligands. We performed explicit solvent molecular dynamics (MD) simulations of these complexes and studied a variety of structural and energetic features. The interfacial residues show a distinct behavior when in complex with different antibiotics. Nevertheless, it was possible to identify some common Hot Spots among several complexes – Lys73, Tyr105 and Glu166. The structural changes that occur during the Molecular Dynamic (MD) simulation lead to the conclusion that these variants have an inherent capacity of adapting to the various antibiotics. This capability might be the reason why they can hydrolyze antibiotics that have not been described until now to be degraded by TEM variants. The results obtained with computational and experimental methodologies for the complex with Imipenem have shown that in order to this type of enzymes be able to acylate the antibiotics, they need to be capable to protect the ligand from water molecules.

Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

High-affinity small molecule-phospholipid complex formation: Binding of siramesine to phosphatidic acid

Siramesine (SRM) is a sigma-2 receptor agonist which has been recently shown to inhibit growth of cancer cells. Fluorescence spectroscopy experiments revealed two distinct binding sites for this drug in phospholipid membranes. More specifically, acidic phospholipids retain siramesine on the bilayer surface due to a high-affinity interaction, reaching saturation at an apparent 1:1 drug-acidic phospholipid stoichiometry, where after the drug penetrates into the hydrocarbon core of the membrane. This behavior was confirmed using Langmuir films. Of the anionic phospholipids, the highest affinity, comparable to the affinities for the binding of small molecule ligands to proteins, was measured for phosphatidic acid (PA, mole fraction Of X-PA = 0.2 in phosphatidylcholine vesicles), yielding a molecular partition coefficient of 240 +/- 80 x 10(6). An MD simulation on the siramesine:PA interaction was in agreement with the above data. Taking into account the key role of PA as a signaling molecule promoting cell growth our results suggest a new paradigm for the development of anticancer drugs, viz. design of small molecules specifically scavenging phospholipids involved in the signaling cascades controlling cell behavior.

How C-terminal carboxyamidation alters the biological activity of peptides from the venom of the eumenine solitary wasp

Inflammatory peptides display different types of post-transcriptional modifications, such as C-terminal amidation, that alter their biological activity. Here we describe the structural and molecular dynamics features of the mast cell degranulating peptide, eumenine mastoparan-AF (EMP-AF-NH2), found in the venom of the solitary wasp, and of its carboxyl-free C-terminal form (EMP-AF-COO-) characterized by a reduced activity. Circular dichroism indicates that both peptides switch from a random coil conformation in water to a helical structure in TFE and SDS micelles. NMR data, in 30% TFE, reveal that the two peptides fold into an alpha-helix spanning most of their length, while they differ in terms of molecular rigidity. To understand the origins of the conformational flexibility observed in the case of EMP-AF-COO-, a 5 ns MD simulation was carried out for each peptide, in an explicit water/TFE environment. The results show that the two peptides differ in an H-bond between Leu14 NH2 and the backbone carbonyl of Ile11. The loss of that H-bond in EMP-AF-COO- leads to a significant modification of its structural dynamics. In fact, as evidenced by essential dynamics analysis, while EMP-AF-NH2 exists mainly as a rigid structure, EMP-AF-COO- presents two helical stretches that fluctuate in some sort of independent fashion. We conclude that the diverse biological activity of the two peptides is not simply due to the reduction of the net positive charge, as generally suggested, but also to a structural perturbation of the amphipathic alpha-helix that affects their ability to perturb the cell membrane.

Why silver deposition is so fast: Solving the enigma of metal deposition

A perfect match: Silver deposition is one of the fastest electrochemical reactions, even though the Ag+ ion loses more than 5 eV solvation energy in the process. This phenomenon, an example of the enigma of metal deposition, was investigated by a combination of MD simulations, DFT, and specially developed theory. At the surface, the Ag+ ion experiences a strong interaction with the sp band of silver, which catalyzes the reaction. Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Interaction of cyclic and linear labaditin peptides with anionic and zwitterionic micelles

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Interaction of cyclic and linear Labaditin peptides with anionic and zwitterionic micelles

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Energielandschaft eines kleinen Modellglasbildners

In dieser Arbeit wurde ein Glasbildner vom Lennard-Jones-Typ mittels einer detaillierten Analyse seiner potentiellen Energielandschaft (PEL) untersucht. Bei hinreichend tiefen Temperaturen ist allein die Kenntnis der effektiven Dichte der Energie der Minima notwendig, um das Konfigurationsintegral und damit thermodynamische Größen des Systems zu berechnen. Aus den Zeitreihen der Energie der Minima zu verschiedenen Temperaturen konnte die effektive Dichte der Energie der Minima bestimmt werden. Durch einen Vergleich der unmittelbar in der MD-Simulation gemessenen thermodynamischen Größen mit Berechnungen in harmonischer Approximation mittels der effektiven Dichte kann direkt auf anharmonische Beiträge in diesen Größen geschlossen werden. Die Minima bestimmen die lokale Struktur der PEL. Doch neben den Minima finden sich auch nichtlokale topologische Elemente. Vornehmlich niederenergetische Minima bilden Täler in der PEL. In dem Temperaturbereich der unterkühlten Flüssigkeit ist das System für lange Zeiten in diesen Tälern gefangen. Die Existenz der Täler verlangsamt die Einteilchen-Dynamik dramatisch und wird letztlich zum zeitlimitierenden Schritt bei der vollständigen Relaxation des Systems. Dies zeigt die Bedeutung der Täler, der langlebigen metastabilen Zustände, bei tiefen Temperaturen für die komplexe Ortsraumdynamik des Systems, die sich zum Beispiel in der nichtexponentiellen Relaxation von Korrelationsfunktionen oder den dynamischen Heterogenitäten widerspiegeln.