Enhancing drug residence time by shielding of intra-protein hydrogen bonds: a case study on CCR2 antagonists. Supplementary Information
|
|
- Herbert Hubbard
- 5 years ago
- Views:
Transcription
1 Enhancing drug residence time by shielding of intra-protein hydrogen bonds: a case study on CCR2 antagonists Aniket Magarkar, Gisela Schnapp, Anna-Katharina Apel, Daniel Seeliger & Christofer S. Tautermann * Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Str. 65, D Biberach a.d. Riss, Germany KEYWORDS: CCR2, GPCR, water networks, drug residence time, buried hydrogen bonds Supplementary Information
2 Results: Section : Analysis of hydrogen bond in CCR2 and ligands MK-82 & 5a. A: Hydrogen bond : Glu29and MK82 Number of H-bonds Tme (ns) Number of H-bonds B: Hydrogen bond : Glu29 and 5a Tme (ns) Supplementary Information Figure : Interactions of ligands in the CCR2 binding pocket: The salt bridge between Glu29 and the ligands remains stable throughout the simulation. Section 2: Comparison of binding site occupied volume vs time for MK-82, 5a and 8 In order to understand the behavior of ligands MK-82, 5a and 8, we analyzed the volume occupied by ligands in the CCR2 binding site. First we calculated distance in the center of mass of protein and center of mass of ligand, as shown in Supporting Information Figure 2A, all of the three ligands are present in the binding site, throughout the simulation. (The error bars denote standard deviation calculated from all simulations). Then we plotted the volume of the ligands in CCR2 binding site vs time for three ligands. As can be seen in supporting Information Figure 2B, 5a occupies higher volume in binding site as compared to 8 and MK-82. Comparing volume occupied by 5a to 8, where the only difference in the structure is absence of Br in 8, it can be concluded that, higher volume of 5a is responsible for the broken water network in the CCR2 binding site. This disruption of water network in turn imparts rigidity to CCR2 binding site rigidity and slower off rate of for 5a.
3 B. MK-82 Ligand 5a Ligand 8 Volume (nm 3 ) Time (ns) Supporting Information Figure 2: (A) Distance between center of mass of protein and ligands (B) Occupied volume by ligands in the CCR2 binding site vs time. Section 3: Biased simulation In order to validate the hypothesis that extra water density is responsible for forming a dynamic H-bonding network and repeatedly disrupting the Tyr2-Glu29 H-bond, we set up a biased simulation, where the water molecules present near Tyr2, for MK-82 and 8 are treated in a different way than rest of the waters. The partial charges on these extra waters are set to zero (shown in Supporting Information Figure 3B in grayscale), to ensure that these water molecules will not be able to participate in H-bonding with other water molecules as well as Tyr2. The MK-82 and 8 systems are then simulated for 5 ns with 5 replicates, with position restraints in z-direction (perpendicular to the membrane) on selected water molecule s oxygen atoms as shown in Supporting Information Figure 3. The position restraints applied (see method section) ensure that the selected water molecules can t be replaced by other normal water molecules, which have H-bonding abilities. Analysis of these trajectories showed that the Tyr2-Glu29 H-bond remained stable throughout the simulation time. This finding allows the interpretation that the Tyr2-Glu29 interaction is usually part of a highly dynamic H-bond network, which loses all dynamics when taking out the additional water molecule.
4 A Tyr2 Glu29 Charged waters B Tyr2 Glu29 Supporting Information Figure 3: Effect of water on the Tyr2-Glu29 H-bond (A) For MK-82 and 8, water molecules form a network of H-bonds, which disrupts the Tyr2-Glu29 interaction (B) When the charges on these molecules were turned off, they are no longer able to influence Tyr2-Glu29 interaction and Tyr2-Glu29 H-bonds are stable throughout 5 ns of simulations for MK-82 (C) and 8 (D) Section 4: Choice of compound 8 over compound 5b for hypothesis validation In order to test the our hypothesis that, extra water density near Tyr2 as shown in Figure 3A (main text) (highlighted in the green circle), is responsible for short RT of MK-82, by disrupting the rigidity of the binding pocket and keeping it flexible. In order to test this hypothesis with further simulations, we pick an additional ligand from the study by Vilums et al. with following criteria: ) chemical structure should be very similar to 5a, 2) CCR2 affinity should be very similar to 5a and 3) RT should be much shorter than for 5a, if possible closer to the RT of MK-82. Compound 8 (Figure ) from their dataset is chosen as it fulfills all of the above-mentioned criteria as the only change in structure is the deletion of the Br in 5a. Although 5b, a diastereomer of 5a has much shorter residence time as compared to 5a, fulfills the criteria no., the difference in affinity is much larger as compared to that of 5a. For this reason, ligand 8 was chosen over 5b. Section 5: Methods Number of H-bonds Initial system preparation: C MK-82 Tyr2-Glu29 H-bond Time (ns) Uncharged waters In total we have performed MD simulations for 5 systems as shown in Supporting Information Table. To start with, the recently solved structure of CCR2/MK-82, (PDB id: 6GPS) was used as a starting point with minor modifications, to resemble the wild type structure of CCR2. Specifically the rubredoxin fusion on intracellular loop 3 was manually removed, still keeping residues R23 and R24. The missing loop between these residues was modeled and energy minimized employing standard homology modeling tools as implemented in MOE (Molecular Operating Environment (MOE), 28.; Chemical Computing Group). Also the Asn75 mutant introduced for thermostabilisation was back-mutated to Gly D 8 Tyr2-Glu29 H-bond Time (ns)
5 to match the wild-type form of CCR2. The final structure was protonated by the Protonate3D procedure as implemented in MOE. Based on the structural similarity of ligands, 5a and 8 with MK-82, they were manually placed in the binding pocket to obtain initial configurations for the MD systems. Tethered energy minimization protocols were applied (first minimization: only H atoms relaxed, second minimization: tethering heavy atoms to.5, third minimization: tethering the protein backbone to.5 all with the Amber:EHT force field as implemented in MOE) to relax ligand poses and remove minor clashes. These structures were then embedded in a membrane bilayer consisting of 86 POPC molecules with 82 water molecules, and Cl - ions were added to neutralize the charges. The system was then equilibrated in following 2 stages. ) Position restraints on protein heavy atoms (fc= KJ mol - nm - ) and ligand (fc= KJ mol - nm - ) and lipid headgroup (fc= KJ mol - nm - ) for ns 2) Position restraints on protein backbone (fc= KJ mol - nm - ) for ns Following this procedure, as shown in supporting table ns or 5 ns of production run was performed for the respective systems. Simulation parameters: All simulations were performed with the GROMACS 5..4 package 2. The protein was parameterized using the AMBER ff99sb-ildn 3 force field, water molecules were described with the SPC/E 4 watermodel, and for lipids S-lipids 5 parameters were used. Newton s equations of motion were integrated by employing the leap-frog algorithm with a time step of 4 fs with virtual sites 6,7. A cutoff of. nm was applied to short-range electrostatic interactions while long-range electrostatics was calculated using the particle mesh Ewald method 8. Van der Waals interactions were truncated at. nm. Covalent bonds containing hydrogen atoms were constrained by the LINCS algorithm 9, and water molecules were held rigid by the SETTLE algorithm. The temperature of the system was maintained at 3 K by the velocity rescaling thermostat and the Parrinello Rahman barostat was used to keep the system pressure at bar. The water network analysis was performed by in-house developed protocol as described in previous studies 2. All visualizations were performed in VMD 3. For the biased system setup for systems 4 and 5, charges on the selected waters were set to by defining a variant of SPC/E water molecules in a separate topology file. A force constant of 5 KJ mol - nm - were applied in z-direction in order to waters to maintain the position and not to be replaced by normal SPC/E water molecules. For H-bond analysis Gromacs analysis tools were used with default parameters, that is, distance cutoff of.35nm. System Ligand Replicas Simulation Length MK-82 5 Unbiased ns 2 5a 5 Unbiased ns Unbiased ns 4 MK-82 5 Biased (uncharged waters) 5 ns Biased (uncharged waters) 5 ns Supporting Information Table : Summary of simulated systems
6 Section 6: Calculation of free energies ΔG based on K i and residence time Supporting Information Figure 4: Energy profiles of two similar compounds with different K i and residence time values. Assuming Boltzmann equilibrium, the change in the free energy of the bound states of two ligands ( G BS) is directly related to the ratio of their K i-values via: G #$ = R T ln, K.[ligand ] K. [ligand 2] 8 For the change in the activation barrier ( G TS), the deviation of the residence time (RT) from G BS-driven behavior is calculated by: G 9$ = R T ln : K.[ligand ] ] ; R T ln :RT[ligand K. [ligand 2] RT[ligand 2] ; with R being the molar gas constant and T stands for the temperature.
7 References: () Vilums, M.; Zweemer, A. J. M.; Barmare, F.; van der Gracht, A. M. F.; Bleeker, D. C. T.; Yu, Z.; de Vries, H.; Gross, R.; Clemens, J.; Krenitsky, P.; Brussee, J.; Stamos, D.; Saunders, J.; Heitman, L. H.; IJzerman, A. P. When Structure Affinity Relationships Meet Structure Kinetics Relationships: 3- ((Inden--Yl)Amino)--Isopropyl-Cyclopentane--Carboxamides as CCR2 Antagonists. Eur. J. Med. Chem. 25, 93, (2) Pronk, S.; Páll, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M. R.; Smith, J. C.; Kasson, P. M.; van der Spoel, D.; Hess, B.; Lindahl, E. GROMACS 4.5: a High-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit. Bioinformatics 23, 29 (7), (3) Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J. L.; Dror, R. O.; Shaw, D. E. Improved Side-Chain Torsion Potentials for the Amber ff99sb Protein Force Field. Proteins 2, 78 (8), (4) Berendsen, H. J. C.; Grigera, J. R.; Straatsma, T. P. The Missing Term in Effective Pair Potentials. J. Phys Chem. 987, 9 (24), (5) Jämbeck, J. P. M.; Lyubartsev, A. P. Another Piece of the Membrane Puzzle: Extending Slipids Further. J Chem. Theory Comput. 22, 9 (), (6) Feenstra, K. A.; Hess, B.; Berendsen, H. J. C. Improving Efficiency of Large Time-Scale Molecular Dynamics Simulations of Hydrogen-Rich Systems. J. Comp. Chem. 999, 2 (8), (7) Melcr, J.; Bonhenry, D.; Timr, Š.; Jungwirth, P. Transmembrane Potential Modeling: Comparison Between Methods of Constant Electric Field and Ion Imbalance. J Chem. Theory Comput. 26, 2 (5), (8) Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 995, 3 (9), (9) Hess, B.; Bekker, H.; Berendsen, H. LINCS: a Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 997, 8 (2), () Jämbeck, J. P. M.; Lyubartsev, A. P. An Extension and Further Validation of an All-Atomistic Force Field for Biological Membranes. J. Chem. Theory Comput. 22, 8 (8), () Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: a New Molecular Dynamics Method. J. Appl. Phy. 98, 52 (2), (2) Tautermann, C. S.; Seeliger, D.; Kriegl, J. M. What Can We Learn From Molecular Dynamics Simulations for GPCR Drug Design? Comput. Struct. Biotechnol. J. 25, 3, 2. (3) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 996, 4 (),
MARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
More informationJournal of Pharmacology and Experimental Therapy-JPET#172536
A NEW NON-PEPTIDIC INHIBITOR OF THE 14-3-3 DOCKING SITE INDUCES APOPTOTIC CELL DEATH IN CHRONIC MYELOID LEUKEMIA SENSITIVE OR RESISTANT TO IMATINIB Manuela Mancini, Valentina Corradi, Sara Petta, Enza
More informationSupporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor
Supporting Material for Microscopic origin of gating current fluctuations in a potassium channel voltage sensor J. Alfredo Freites, * Eric V. Schow, * Stephen H. White, and Douglas J. Tobias * * Department
More informationSupplemntary Infomation: The nanostructure of. a lithium glyme solvate ionic liquid at electrified. interfaces
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 207 Supplemntary Infomation: The nanostructure of a lithium glyme solvate ionic liquid
More informationRoutine access to millisecond timescale events with accelerated molecular dynamics
Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING
More informationSupplementary Information. Surface Microstructure Engenders Unusual Hydrophobicity in. Phyllosilicates
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supplementary Information Surface Microstructure Engenders Unusual Hydrophobicity in Phyllosilicates
More informationTailoring the Properties of Quadruplex Nucleobases for Biological and Nanomaterial Applications
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Supporting Information for: Tailoring the Properties of Quadruplex Nucleobases
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 28 Electronic Supplementary Information ph-dependent Cooperativity and Existence of
More informationUnit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement
Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement Xin Yin, Yeqi Shi, Yanbing Wei, Yongho Joo, Padma Gopalan, Izabela Szlufarska,
More informationStructure Investigation of Fam20C, a Golgi Casein Kinase
Structure Investigation of Fam20C, a Golgi Casein Kinase Sharon Grubner National Taiwan University, Dr. Jung-Hsin Lin University of California San Diego, Dr. Rommie Amaro Abstract This research project
More informationSupporting Information Quantitative characterization of the binding and unbinding of millimolar drug fragments with molecular dynamics simulations
Supporting Information Quantitative characterization of the binding and unbinding of millimolar drug fragments with molecular dynamics simulations Simulation Details All binding simulations were run on
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationStructural heterogeneity of the μ-opioid receptor s conformational ensemble in the. apo state
Structural heterogeneity of the μ-opioid receptor s conformational ensemble in the apo state Diniz M. Sena Jr 1,2,3,&, Xiaojing Cong 1,2,&,*, Alejandro Giorgetti 2,4, Achim Kless 5, Paolo Carloni 1,2,*
More informationState Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academic of Sciences, Dalian , P. R. China.
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2015 Supplementary Information for: Theoretical exploration of MgH 2 and graphene nano-flake in cyclohexane:
More informationSupporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of
Supporting Information: Improved Parametrization of Lithium, Sodium, Potassium, and Magnesium ions for All-Atom Molecular Dynamics Simulations of Nucleic Acid Systems. Jejoong Yoo and Aleksei Aksimentiev,,
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationComparing crystal structure of M.HhaI with and without DNA1, 2 (PDBID:1hmy and PDBID:2hmy),
Supporting Information 1. Constructing the starting structure Comparing crystal structure of M.HhaI with and without DNA1, 2 (PDBID:1hmy and PDBID:2hmy), we find that: the RMSD of overall structure and
More informationSupporting information:
Supporting information: A simple method to measure protein side-chain mobility using NMR chemical shifts. Mark V. Berjanskii, David S. Wishart* Department of Computing Science, University of Alberta, Edmonton,
More informationSupporting Information for. Hydrogen Bonding Structure at Zwitterionic. Lipid/Water Interface
Supporting Information for Hydrogen Bonding Structure at Zwitterionic Lipid/Water Interface Tatsuya Ishiyama,, Daichi Terada, and Akihiro Morita,, Department of Applied Chemistry, Graduate School of Science
More informationSupporting Information
Supporting Information Structure and Dynamics of Uranyl(VI) and Plutonyl(VI) Cations in Ionic Liquid/Water Mixtures via Molecular Dynamics Simulations Katie A. Maerzke, George S. Goff, Wolfgang H. Runde,
More informationCarlo Camilloni, Andrea Cavalli,, and Michele Vendruscolo INTRODUCTION
pubs.acs.org/jpcb Assessment of the Use of NMR Chemical Shifts as Replica-Averaged Structural Restraints in Molecular Dynamics Simulations to Characterize the Dynamics of Proteins Carlo Camilloni, Andrea
More informationTowards fast and accurate binding affinity. prediction with pmemdgti: an efficient. implementation of GPU-accelerated. Thermodynamic Integration
Towards fast and accurate binding affinity prediction with pmemdgti: an efficient implementation of GPU-accelerated Thermodynamic Integration Tai-Sung Lee,, Yuan Hu, Brad Sherborne, Zhuyan Guo, and Darrin
More informationSUPPLEMENTARY MATERIAL. Supplementary material and methods:
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2015 SUPPLEMENTARY MATERIAL Supplementary material and methods: - Computational
More informationMolecular modeling. A fragment sequence of 24 residues encompassing the region of interest of WT-
SUPPLEMENTARY DATA Molecular dynamics Molecular modeling. A fragment sequence of 24 residues encompassing the region of interest of WT- KISS1R, i.e. the last intracellular domain (Figure S1a), has been
More informationSupporting Information
Supporting Information The Molecular Mechanism of P2Y 1 Receptor Activation Shuguang Yuan,* H. C. Stephen Chan, Horst Vogel, Slawomir Filipek, Raymond C. Stevens,* and Krzysztof Palczewski* anie_201605147_sm_miscellaneous_information.pdf
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationSAM Teacher s Guide Protein Partnering and Function
SAM Teacher s Guide Protein Partnering and Function Overview Students explore protein molecules physical and chemical characteristics and learn that these unique characteristics enable other molecules
More informationUniversal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information
Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information B. Shadrack Jabes, Dusan Bratko, and Alenka Luzar Department of Chemistry,
More informationMedical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium.
Supporting Information Towards peptide vaccines against Zika virus: Immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins Muhammad Usman Mirza
More informationSupporting Information
pk a values of titrable amino acids at the water/membrane interface Vitor H. Teixeira, Diogo Vila-Viçosa, Pedro B. P. S. Reis, and Miguel Machuqueiro Centro de Química e Bioquímica, Departamento de Química
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationSupporting Information
Supporting Information ph-responsive self-assembly of polysaccharide through a rugged energy landscape Brian H. Morrow, Gregory F. Payne, and Jana Shen Department of Pharmaceutical Sciences, School of
More informationMulti-scale approaches in description and design of enzymes
Multi-scale approaches in description and design of enzymes Anastassia Alexandrova and Manuel Sparta UCLA & CNSI Catalysis: it is all about the barrier The inside-out protocol: Big Aim: development of
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2016 Electronic Supplementary Information Active Site Structure and Absorption Spectrum of Channelrhodopsin-2
More informationMolecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains
Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains Fatemeh Khalili-Araghi, Emad Tajkhorshid, Benoît Roux, and Klaus Schulten Department of Physics, Department of Biochemistry,
More informationProceedings of the ASME 2009 International Mechanical Engineering Congress & Exposition
Proceedings of the ASME 9 International Mechanical Engineering Congress & Exposition IMECE9 November 3-9, Lake Buena Vista, Florida, USA Proceedings of the ASME International Mechanical Engineering Congress
More informationMolecular Dynamics. Molecules in motion
Molecular Dynamics Molecules in motion 1 Molecules in mo1on Molecules are not sta1c, but move all the 1me Source: h9p://en.wikipedia.org/wiki/kine1c_theory 2 Gasses, liquids and solids Gasses, liquids
More informationT H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y. jgp
S u p p l e m e n ta l m at e r i a l jgp Lee et al., http://www.jgp.org/cgi/content/full/jgp.201411219/dc1 T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y S u p p l e m e n ta l D I S C U S
More informationSupplemental Material for Global Langevin model of multidimensional biomolecular dynamics
Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics Norbert Schaudinnus, Benjamin Lickert, Mithun Biswas and Gerhard Stock Biomolecular Dynamics, Institute of Physics,
More informationProbing α/β Balances in Modified Amber Force Fields from a Molecular Dynamics Study on a ββα Model Protein (1FSD)
Benchmark Study of Modified Amber Force Fields on ββα Protein Bull. Korean Chem. Soc. 214, Vol. 35, No. 6 1713 http://dx.doi.org/1.512/bkcs.214.35.6.1713 Probing α/β Balances in Modified Amber Force Fields
More informationSupplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached
Supplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached to HpUreI. Urea hydrolysis products 2NH 3 and 1CO 2
More informationEnhancing Specificity in the Janus Kinases: A Study on the Thienopyridine. JAK2 Selective Mechanism Combined Molecular Dynamics Simulation
Electronic Supplementary Material (ESI) for Molecular BioSystems. This journal is The Royal Society of Chemistry 2015 Supporting Information Enhancing Specificity in the Janus Kinases: A Study on the Thienopyridine
More informationThermodynamics. Entropy and its Applications. Lecture 11. NC State University
Thermodynamics Entropy and its Applications Lecture 11 NC State University System and surroundings Up to this point we have considered the system, but we have not concerned ourselves with the relationship
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationQuantum Interference and Selectivity through Biological Ion Channels
1 SUPPLEMENTARY INFORMATION Quantum Interference and Selectivity through Biological Ion Channels Vahid Salari 1,3,*, Hamidreza Naeij 2, Afshin Shafiee 2,3 1 Department of Physics, Isfahan University of
More informationChoosing weights for simulated tempering
PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural
More informationAmino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information
Amino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information Grzegorz Nawrocki and Marek Cieplak Institute of Physics, Polish Academy of Sciences,
More informationUser Guide for LeDock
User Guide for LeDock Hongtao Zhao, PhD Email: htzhao@lephar.com Website: www.lephar.com Copyright 2017 Hongtao Zhao. All rights reserved. Introduction LeDock is flexible small-molecule docking software,
More informationDISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING:
DISCRETE TUTORIAL Agustí Emperador Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: STRUCTURAL REFINEMENT OF DOCKING CONFORMATIONS Emperador
More informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationDrug binding and subtype selectivity in G-protein-coupled receptors
Drug binding and subtype selectivity in G-protein-coupled receptors Albert C. Pan The Aspect of Time in Drug Design Schloss Rauischholzhausen Marburg, Germany Thursday, March 27 th, 2014 D. E. Shaw Research
More informationDiffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study
Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron
More informationChapter 4. Glutamic Acid in Solution - Correlations
Chapter 4 Glutamic Acid in Solution - Correlations 4. Introduction Glutamic acid crystallises from aqueous solution, therefore the study of these molecules in an aqueous environment is necessary to understand
More informationHomology modeling. Dinesh Gupta ICGEB, New Delhi 1/27/2010 5:59 PM
Homology modeling Dinesh Gupta ICGEB, New Delhi Protein structure prediction Methods: Homology (comparative) modelling Threading Ab-initio Protein Homology modeling Homology modeling is an extrapolation
More informationNature Structural & Molecular Biology: doi: /nsmb Supplementary Figure 1
Supplementary Figure 1 Crystallization. a, Crystallization constructs of the ET B receptor are shown, with all of the modifications to the human wild-type the ET B receptor indicated. Residues interacting
More informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
More informationCalcium Binding to Calmodulin by Molecular. Dynamics with Effective Polarization
Calcium Binding to Calmodulin by Molecular Dynamics with Effective Polarization Miriam Kohagen, Martin Lepsik, and Pavel Jungwirth Institute of Organic Chemistry and Biochemistry, Academy of Sciences of
More informationtconcoord-gui: Visually Supported Conformational Sampling of Bioactive Molecules
Software News and Updates tconcoord-gui: Visually Supported Conformational Sampling of Bioactive Molecules DANIEL SEELIGER, BERT L. DE GROOT Computational Biomolecular Dynamics Group, Max-Planck-Institute
More informationNMR study of complexes between low molecular mass inhibitors and the West Nile virus NS2B-NS3 protease
University of Wollongong Research Online Faculty of Science - Papers (Archive) Faculty of Science, Medicine and Health 2009 NMR study of complexes between low molecular mass inhibitors and the West Nile
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.1821 Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways Kai J. Kohlhoff 1,4 *, Diwakar Shukla 1,2, Morgan Lawrenz 2, Gregory R. Bowman 2,
More informationSupporting Information
Supporting Information Constant ph molecular dynamics reveals ph-modulated binding of two small-molecule BACE1 inhibitors Christopher R. Ellis 1,, Cheng-Chieh Tsai 1,, Xinjun Hou 2, and Jana Shen 1, 1
More informationEquilibrated atomic models of outward-facing P-glycoprotein and effect of ATP binding on structural dynamics (Supplementary Information)
Equilibrated atomic models of outward-facing P-glycoprotein and effect of ATP binding on structural dynamics (Supplementary Information) Lurong Pan 1 and Stephen G. Aller 2 * 1,2 Department of Pharmacology
More informationPosition of Transmembrane Helix 6 Determines Receptor G Protein Coupling Specificity
Position of Transmembrane Helix 6 Determines Receptor G Protein Coupling Specificity Alexander S. Rose,, Matthias Elgeti, Ulrich Zachariae, Helmut Grubmüller, Klaus Peter Hofmann,, Patrick Scheerer,,#
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationAccurate Structure Prediction and Conformational Analysis of Cyclic Peptides with Residue-Specific Force Fields
Accurate Structure Prediction and Conformational Analysis of Cyclic Peptides with Residue-Specific Force Fields Hao Geng, Fan Jiang, *, and Yun-Dong Wu *,, School of Chemical Biology and Biotechnology,
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationExploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials
Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research
More information7.012 Problem Set 1. i) What are two main differences between prokaryotic cells and eukaryotic cells?
ame 7.01 Problem Set 1 Section Question 1 a) What are the four major types of biological molecules discussed in lecture? Give one important function of each type of biological molecule in the cell? b)
More informationRanjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India. 1 st November, 2013
Hydration of protein-rna recognition sites Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India 1 st November, 2013 Central Dogma of life DNA
More informationMicroscopic analysis of protein oxidative damage: effect of. carbonylation on structure, dynamics and aggregability of.
Microscopic analysis of protein oxidative damage: effect of carbonylation on structure, dynamics and aggregability of villin headpiece Drazen etrov 1,2,3 and Bojan Zagrovic 1,2,3,* Supplementary Information
More informationOxygen Binding in Hemocyanin
Supporting Information for Quantum Mechanics/Molecular Mechanics Study of Oxygen Binding in Hemocyanin Toru Saito and Walter Thiel* Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470
More informationWhat binds to Hb in addition to O 2?
Reading: Ch5; 158-169, 162-166, 169-174 Problems: Ch5 (text); 3,7,8,10 Ch5 (study guide-facts); 1,2,3,4,5,8 Ch5 (study guide-apply); 2,3 Remember Today at 5:30 in CAS-522 is the second chance for the MB
More informationSupplementary Information
Supplementary Information Resveratrol Serves as a Protein-Substrate Interaction Stabilizer in Human SIRT1 Activation Xuben Hou,, David Rooklin, Hao Fang *,,, Yingkai Zhang Department of Medicinal Chemistry
More informationKd = koff/kon = [R][L]/[RL]
Taller de docking y cribado virtual: Uso de herramientas computacionales en el diseño de fármacos Docking program GLIDE El programa de docking GLIDE Sonsoles Martín-Santamaría Shrödinger is a scientific
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/330/6002/341/dc1 Supporting Online Material for Atomic-Level Characterization of the Structural Dynamics of Proteins David E. Shaw,* Paul Maragakis, Kresten Lindorff-Larsen,
More informationpmx: Automated Protein Structure and Topology Generation for Alchemical Perturbations
pmx: Automated Protein Structure and Topology Generation for Alchemical Perturbations Vytautas Gapsys,* [a] Servaas Michielssens, [a] Daniel Seeliger,* [b] and Bert L. de Groot* [a] Computational protein
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationChapters 11 and 12: Intermolecular Forces of Liquids and Solids
1 Chapters 11 and 12: Intermolecular Forces of Liquids and Solids 11.1 A Molecular Comparison of Liquids and Solids The state of matter (Gas, liquid or solid) at a particular temperature and pressure depends
More informationComputational modeling of G-Protein Coupled Receptors (GPCRs) has recently become
Homology Modeling and Docking of Melatonin Receptors Andrew Kohlway, UMBC Jeffry D. Madura, Duquesne University 6/18/04 INTRODUCTION Computational modeling of G-Protein Coupled Receptors (GPCRs) has recently
More informationUnfolding CspB by means of biased molecular dynamics
Chapter 4 Unfolding CspB by means of biased molecular dynamics 4.1 Introduction Understanding the mechanism of protein folding has been a major challenge for the last twenty years, as pointed out in the
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature11085 Supplementary Tables: Supplementary Table 1. Summary of crystallographic and structure refinement data Structure BRIL-NOP receptor Data collection Number of crystals 23 Space group
More informationSupporting information
Electronic Supplementary Material ESI for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Supporting information Understanding three-body contributions to coarse-grained force
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature11524 Supplementary discussion Functional analysis of the sugar porter family (SP) signature motifs. As seen in Fig. 5c, single point mutation of the conserved
More informationComplex Structure Modeling
Molecular Modeling and Simulation On July 5, 2017 Complex Structure Modeling Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More informationIntroduction to Structure Preparation and Visualization
Introduction to Structure Preparation and Visualization Created with: Release 2018-4 Prerequisites: Release 2018-2 or higher Access to the internet Categories: Molecular Visualization, Structure-Based
More informationProteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability
Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Dr. Andrew Lee UNC School of Pharmacy (Div. Chemical Biology and Medicinal Chemistry) UNC Med
More informationSUPPLEMENTARY INFORMATION
Supplementary Results DNA binding property of the SRA domain was examined by an electrophoresis mobility shift assay (EMSA) using synthesized 12-bp oligonucleotide duplexes containing unmodified, hemi-methylated,
More informationMolecular Dynamics in practice with GROMACS
Molecular Dynamics in practice with GROMACS GROMACS is one of the wold s fastest software package for molecular dynamics simulations. One can find many helpful materials, manual as well as to download
More informationHow is molecular dynamics being used in life sciences? Davide Branduardi
How is molecular dynamics being used in life sciences? Davide Branduardi davide.branduardi@schrodinger.com Exploring molecular processes with MD Drug discovery and design Protein-protein interactions Protein-DNA
More informationarxiv: v1 [physics.chem-ph] 30 Jul 2012
Is ion channel selectivity mediated by confined water? Diego Prada-Gracia and Francesco Rao School of Soft Matter, Freiburg Institute for Advanced Studies, Freiburg, Germany arxiv:27.6953v [physics.chem-ph]
More informationWhat makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces: Electronic Supplementary
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry B. This journal is The Royal Society of Chemistry 21 What makes a good graphene-binding peptide? Adsorption of amino acids and
More informationDon t forget to bring your MD tutorial. Potential Energy (hyper)surface
Don t forget to bring your MD tutorial Lab session starts at 1pm You will have to finish an MD/SMD exercise on α-conotoxin in oxidized and reduced forms Potential Energy (hyper)surface What is Force? Energy
More informationToward an Understanding of GPCR-ligand Interactions. Alexander Heifetz
Toward an Understanding of GPCR-ligand Interactions Alexander Heifetz UK QSAR and ChemoInformatics Group Conference, Cambridge, UK October 6 th, 2015 Agenda Fragment Molecular Orbitals (FMO) for GPCR exploration
More informationIdentifying Interaction Hot Spots with SuperStar
Identifying Interaction Hot Spots with SuperStar Version 1.0 November 2017 Table of Contents Identifying Interaction Hot Spots with SuperStar... 2 Case Study... 3 Introduction... 3 Generate SuperStar Maps
More informationCHEM 4170 Problem Set #1
CHEM 4170 Problem Set #1 0. Work problems 1-7 at the end of Chapter ne and problems 1, 3, 4, 5, 8, 10, 12, 17, 18, 19, 22, 24, and 25 at the end of Chapter Two and problem 1 at the end of Chapter Three
More information2. In regards to the fluid mosaic model, which of the following is TRUE?
General Biology: Exam I Sample Questions 1. How many electrons are required to fill the valence shell of a neutral atom with an atomic number of 24? a. 0 the atom is inert b. 1 c. 2 d. 4 e. 6 2. In regards
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationIntroduction FULL PAPER. DOI: /jcc Journal of Computational Chemistry 2017, DOI: /jcc
Efficient Estimation of Binding Free Energies between Peptides and an MHC Class II Molecule Using Coarse- Grained Molecular Dynamics Simulations with a Weighted Histogram Analysis Method Ming Huang, Wenjun
More information