Type II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G2019S.
|
|
- Sophie Hancock
- 5 years ago
- Views:
Transcription
1 Type II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G219S. Min Liu@&*, Samantha A. Bender%*, Gregory D Cuny@, Woody Sherman, Marcie Glicksman@ Soumya S. Harvard NeuroDiscovery Center, Harvard University, 65 Landsdowne St., #452, Cambridge, MA 2139 % Department of Neurology, Brigham and Women's Hospital # Center for Neurologic Diseases, Brigham and Women's Hospital Schrodinger, 12 W. 45th Street, New York, NY
2 Materials and Methods: Homology Modeling and docking: The LRRK2 kinase domain, between residues , was modeled using Modeller. Briefly, the main criteria in homology modeling were template selection and sequence alignment between the target and the template. The top three hits based on sequence identity were the ack1 kinase (31% identity) and B-raf (33% identity). In this case, the template of B-raf kinase which had 33% sequence identity was used for homology modeling since this enzyme had higher sequence conservation around the active site region compared to the other kinases. The Cα RMSD and the backbone RMS deviations for the model and the template crystal structure were < Å and < 1.2 Å respectively. The best model was subjected to geometric evaluations using PROCHECK with an overall G-value of -4. Ramachandran plots indicated that >9% of the residues are in the allowed region of the map. Standard bond lengths and bond angles of the model were determined using WHAT IF with an RMS-Z score of.889 and.91 suggesting that the model is of high quality. In addition, we used Prime, Rosetta and SwissModel to generate models for LRRK2 kinase domain for the same set of residues. The The Cα RMSD and the backbone RMS deviations between any pair of structures was < 1.2Å. A second set of models were generated using the recently determined structure of Roco 4 (pdb code: 4FG) by Wittinghoffer and co-workers. Structural comparison between this model and the model of LRRK2 generated using B-Raf template were virtually identical. Key residue positions including the catalytic and hydrophobic spines described by Taylor and coworkers were in the expected positions. Key catalytic residues including H1994, D1996 were in
3 the expected positions. Comparison of the hinge region revealed that the B-raf hinge was closer in sequence compared to the Roco 4 hinge and was therefore chosen to be model for further studies. Docking of ATP molecule was performed using GLIDE v2.2 (Schrödinger inc.). A 1Å search grid was used using the center of mass of glycine rich loop and the activation loop as the point of origin. Maestro v9.2 (Schrödinger inc.) was used to interactively place the Mg2+ atom near the γ-phosphate of ATP within reasonable bonding distance and geometry based on structures of other kinases in the protein database. Superpositioning wherever necessary was carried out using the modeling software Coot v6.2. Loop modeling and preliminary MD thermalization of the DYG-in and DYG-out conformations: The conformation switch between the active and inactive conformations requires two correlated motions at a protein level. The first motion involves repositioning of the N-terminal β-sheet region and C-helix (Figure 1c) breaking a salt-bridge between K196 and E192 when the kinase switches from active to inactive form. The second motion involves a flip of the DYG-motif, where D217 and Y218 exchange positions (DYG-in to DYG-out). A two step modeling procedure was carried out using Prime and HingeProt to built a DYG-out and open kinase conformation. Briefly, loop conformation were sampled with the Prime extended sampling settings. Side chains within 7.5Å of all residues remained flexible and all structures within a 5 kcal/mol energy window were saved. Prime uses a hierarchal approach where the two major
4 steps are dihedral angle sampling of the backbone to generate possible loop conformations (no steric clashes with other protein parts or inhibitors is ensured) followed by side-chain optimization of loop structures further chosen for optimization. The algorithm for side-chain optimization uses sampling from a highly detailed (1 resolution) rotamer library constructed by Xiang and Honig from a database of 297 proteins. The final models were subject to MD thermalization to determine if they were stable structure. In the thermalization scheme, MD simulations of the two conformational states of the kinase (WT, G219S and I22T) described above were carried out using the Desmond 3. (D.E Shaw Research and Schrödinger inc.) and the OPLS-AA 25 force field. Water solvation was described by immersion of the DYG-in (closed) and DYG-out (open) structures in 26, TIP3P waters. Electroneutrality was ensured by adding three Na+ ions to the system. Periodic boundary conditions and the particle mesh Ewald (PME) method (to account for long-range electrostatic interactions) were used throughout. Bonds involving hydrogens were constrained using the SHAKE algorithm, and a time step integration of 2fs were used for all simulations. A steepest-descent minimization and thermalization scheme was applied to all of the initial structures. The systems were heated from to 3 K in 2ps, keeping the Cα atoms fixed in their original positions. In the next step, all of the constraints were lifted, and the equilibration was continued in the isobaric isothermal ensemble with Nose Hoover thermostats for 5. ns.
5 Metadynamic simulations: After MD thermalization, the DYG-in (closed) and DYG-out (open) conformations of LRRK2 were submitted to the Yale MolMov morphing server ( Starting from our initial and final structures, we obtained ten interpolated intermediates. Each of the ten morphing conformations was then hydrated and thermalized at 3 K using the scheme described above. Notably, some morphing conformations were quite stable during the thermalization, while others showed a rapid increase of the rmsd of the Cα atoms. Finally total of seven conformations (two initial and five intermediates) starting from the DYG-in (closed state) to the DYG-out (open state), and the five intermediate morphing conformers were used to generate the guess path for metadynamicsbased computations. The parameters of the Gaussians were chosen on the basis of pilot calculations as a compromise between accuracy of the simulation coinciding with the morphing path (see supporting information for details) and the computational cost. After several pilot calculations, we defined two variables, s and z, which are able to describe the position of a point in configurational space closely resembling the morphing path generated using the Prime modeled DYG-in and DYG-out states. In case of distance measurements, two sets of atoms whose center of masses were positioned on either K196 or E192 were used as one collective variable (z). The center of mass of the activation loop was calculated based on the position of the Cα atoms of residues Y218, A221, C224, R226 and M227 (s). For metadynamics simulations based on torsion angles to measure change in conformational freedom of the loop during a G219 S219 and I22 T22 mutations, the φ and ψ angles of G219 /S219 and were used
6 as collective variables to measure rotational conformational barriers. The weighted histogram analysis method was applied to a series of simulations in order to obtain a free-energy profile using s and z.
7 (a) (b) [s] [z] Y218 Inactive A221 C224 K196 E192 open R226 M227 active closed (c) [φ]/[ψ] Y218 (DYG-in) Y218 (DYG-out) Figure S1, Supporting information
8 (a) WT 1.2 A D v/[e], min-1 v/[e], min-1 G219S 5 um 25 um 12.5 um 6.25 um um um [ATP], µ M 1.2 (kcat)atp min-1 (kcat)atp, min-1 [ATP], µ M B E C (kcat/km)atp min-1 (kcat/km)atp, min-1µ M [Imatinib], µ M [Imatinib], µ M [Imatinib], µ M Figure S2, Supporting information F [Imatinib], µ M
9 (b) WT G219S 3. kcat, min-1 kcat, min A D 1 kcat/km, min-1 µ M-1 kcat/km, min-1µ M-1 B C.8 8 Ki,app, µ M Ki,app, µ M 5 [Sorafenib], µ M E [Sorafenib], um 1 3 [Sorafenib], µ M [Sorafenib], um F [ATP], um Figure S2 (continued), Supporting information [ATP], µ M
10 (c) G219S A um 5 um 2.5 um 1.25 um.6 um um D 5 um 2.5 um 1.25 um.6 um.3 um um v/[e], min-1 v/[e], min WT kcat, min-1 kcat, min [Bosutinib], um 3 kcat/km, min-1µ M-1 4 kcat/km 1 12 E [Bosutinib], um C [ATP], um [ATP], um B [Bosutinib], um Figure S2 (continued), Supporting information F [Bosutinib], um 1 12
11 (a) N-terminal β-sheet domain Helix 1 Hinge loop ATP Activation loop C-lobe hydrophobic helix Figure S3, Supporting information
12 (b) K196 E192 ATP binding hinge ATP D217 G219S 2+ Mg C224 H1998 D1994 K1996 Figure S3 (continued), supporting information C225
13 (c) ATP binding pocket (DYG-out) ATP binding pocket (DYG-in) DFG in Figure S3 (continued), supporting information DFG-out
14 LRRK2G219 LRRK2S kcal/mol Figure S4, supporting information
15 Bosutinib DYG-in Imatinib DYG-in Figure S5, supporting information Sorafenib DYG-in Ponatinib DYG-in
16 Bosutinib DYG-out Imatinib DYG-out Figure S5 (continued), supporting information Sorafenib DYG-out Ponatinib DYG-out
17 Table S1: DFG-in (Ramachandran angles) Phi Psi LRRK2 A:ASP 217 A:TYR 218 A:GLY 219 ckit A:ASP 81 A:PHE 811 A:GLY Aurora A:ASP 274 A:PHE 275 A:GLY Src A:ASP 44 A:PHE 45 A:GLY cabl A:ASP 381 A:PHE 382 A:GLY B-raf A:ASP 594 A:PHE 595 A:GLY EPHA A:ASP 764 A:PHE 765 A:GLY LCK A:ASP 382 A:PHE 383 A:GLY MK14 A:ASP 168 A:PHE 169 A:GLY
18 Table S1 (continued) DFG-out (Ramachandran angles) LRRK2 A:ASP 217 A:TYR 218 A:GLY 219 Phi Psi ckit A:ASP 81 A:PHE 811 A:GLY Aurora A:ASP 274 A:PHE 275 A:GLY Src A:ASP 44 A:PHE 45 A:GLY cabl A:ASP 381 A:PHE 382 A:GLY B-raf A:ASP 593 A:PHE 594 A:GLY EPHA A:ASP 764 A:PHE 765 A:GLY LCK A:ASP 382 A:PHE 383 A:GLY MK14 A:ASP 168 A:PHE 169 A:GLY
19 Figure legends: Figure S1: Animated views of the collective variables used in the metadynamics simulation (a) This shows the movement of the N-terminal ATP binding β-sheet domain bearing K196 and a slight rotation of C-helix bearing E192. In the closed form the two residues are with 2.Å of each other leading to the formation of a salt-bridge while in the open form they move away >4.Å from each other. This center of masses for this motion represents the first collective variable (z) in our metadynamic simulations. (b) The position of the activation loop shown in active (green) and inactive (red) form of LRRK2. The CPK spheres shown in green represent the atoms used to define the center of mass for the loop and the collective variable (z) when it switches between active and inactive form. (c) Flip of the DYG motif of LRRK2 mediated by rotation of φ/ψ angles of the activation loop. The figure here is generated using the model coordinates for WT LRRK2. Figure S2: (a) Inhibition of LRRK2-catalyzed LRRKtide phosphorylation by Imatinib. Initial velocities was measured for the mutant G219S (A) and WT (D) as a function of [ATP] at different [Imatinib] as shown in the panel and at a fixed LRRKtide concentration of 5 µm. B & C (for G219S) or E & F (for WT): Imatinib concentration dependencies of (kcat)atp and (kcat/km)atp apparent values derived from analysis of the data of panel A & D. (b) Inhibition of LRRK2-catalyzed LRRKtide phosphorylation by Sorafenib. A & B (for G219S) or D & E (for WT): Sorafenib concentration dependencies of (kcat)atp and (kcat/km)atp apparent values derived from analysis of the initial velocity data. C (for the mutant) & F (for WT): Ki,app values were plotted against [ATP]. Due to the apparent partial inhibition of Sorafenib at higher ATP concentration, it is difficult to determine the inhibition mechanism using the replot methods for WT. Instead, the shape of the plot of Ki,app vs. [ATP] was used. The Ki,app values of Sorafenib linearly increases with [ATP] for the mutant G219S, suggesting an ATP competitive
20 inhibition. The Ki,app values of Sorafenib is independent of [ATP] for WT, suggesting an ATP noncompetitive inhibition. (c) Inhibition of LRRK2-catalyzed LRRKtide phosphorylation by Bosutinib. Initial velocities was measured for the mutant G219S (A) and WT (D) as a function of [ATP] at different [Bosutinit] as shown in the panel and at a fixed LRRKtide concentration of 5 µm. B & C (for G219S) or E & F (for WT): Bosutinib concentration dependencies of (kcat)atp and (kcat/km)atp apparent values derived from analysis of the data of panel A & D. Figure S3: (a) Overall architecture of LRRK2 kinase domain generated using Modeller. The structure shows an archetypal kinase fold with two major domains. The N-terminal β-sheet rich domain and the largely α-helical C-terminal domain with the ATP docked into the binding cleft. The positions of the Parkinson s disease linked mutations G219S and I22T are shown in CPK. (b) Structural details of the ATP binding site of LRRK2 with ATP docked in. ATP phosphate groups make hydrogen bonds with the backbone atoms of the glycine-rich loop. The adenine substructure makes hydrogen bonds with the hinge region. In addition, the DYG motif on the activation loop stabilizes a Mg 2+ ion, which in turns makes salt-bridge like interactions with the phosphates of ATP. Residues from the catalytic loop K1996 and H1998 provide additional stabilizing interactions to the bound ATP. (c) Surface representation of the DYG-in and DYG-out structures of LRRK2. Green denotes hydrophobic residues, red denotes negatively charged and blue denotes positively charged residues. In the DYG-out conformation, the hydrophobic Tyr occupies part of the ATP binding pocket Figure S4: Free energy surface (FES) generated from metadynamic simulation of LRRK WT and G219S (top panel) generated using the φ/ψ angles of G219 and S219. There is a significant lowering of accessible conformation is observed for G219S relative to WT.
21 Figure S5: Ligand interaction diagrams for G219S in the DYG-in state and WT in DYG-out state with IFD docked DFG-out inhibitors, bosutinib, imatinib, sorafenib and poantinib. This figure accompanies figure 5 in the main text showing all of the residues that interact with the ligand. The hydrophobic residues are shown in green, negatively charged in red and positively charged in blue. A grey hue around the ligand atoms indicate solvent exposure and pink arrows indicate hydrogen from donor to acceptor.
Protein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron.
Protein Dynamics The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Below is myoglobin hydrated with 350 water molecules. Only a small
More informationIntroduction to Comparative Protein Modeling. Chapter 4 Part I
Introduction to Comparative Protein Modeling Chapter 4 Part I 1 Information on Proteins Each modeling study depends on the quality of the known experimental data. Basis of the model Search in the literature
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 informationDesign of a Novel Globular Protein Fold with Atomic-Level Accuracy
Design of a Novel Globular Protein Fold with Atomic-Level Accuracy Brian Kuhlman, Gautam Dantas, Gregory C. Ireton, Gabriele Varani, Barry L. Stoddard, David Baker Presented by Kate Stafford 4 May 05 Protein
More informationFlexibility and Constraints in GOLD
Flexibility and Constraints in GOLD Version 2.1 August 2018 GOLD v5.6.3 Table of Contents Purpose of Docking... 3 GOLD s Evolutionary Algorithm... 4 GOLD and Hermes... 4 Handling Flexibility and Constraints
More informationBioinformatics. Macromolecular structure
Bioinformatics Macromolecular structure Contents Determination of protein structure Structure databases Secondary structure elements (SSE) Tertiary structure Structure analysis Structure alignment Domain
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 informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/3/4/e1600663/dc1 Supplementary Materials for A dynamic hydrophobic core orchestrates allostery in protein kinases Jonggul Kim, Lalima G. Ahuja, Fa-An Chao, Youlin
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 informationComputational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data
Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data ABSTRACT Keyword Lei Shi 1 Advisor: Gianluigi Veglia 1,2 Department of Chemistry 1, & Biochemistry, Molecular
More informationExamples of Protein Modeling. Protein Modeling. Primary Structure. Protein Structure Description. Protein Sequence Sources. Importing Sequences to MOE
Examples of Protein Modeling Protein Modeling Visualization Examination of an experimental structure to gain insight about a research question Dynamics To examine the dynamics of protein structures To
More informationHOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target.
HOMOLOGY MODELING Homology modeling, also known as comparative modeling of protein refers to constructing an atomic-resolution model of the "target" protein from its amino acid sequence and an experimental
More informationDocking. GBCB 5874: Problem Solving in GBCB
Docking Benzamidine Docking to Trypsin Relationship to Drug Design Ligand-based design QSAR Pharmacophore modeling Can be done without 3-D structure of protein Receptor/Structure-based design Molecular
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 informationBuilding 3D models of proteins
Building 3D models of proteins Why make a structural model for your protein? The structure can provide clues to the function through structural similarity with other proteins With a structure it is easier
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 informationSUPPLEMENTARY INFORMATION
Figure S1. Secondary structure of CAP (in the camp 2 -bound state) 10. α-helices are shown as cylinders and β- strands as arrows. Labeling of secondary structure is indicated. CDB, DBD and the hinge are
More informationStructure to Function. Molecular Bioinformatics, X3, 2006
Structure to Function Molecular Bioinformatics, X3, 2006 Structural GeNOMICS Structural Genomics project aims at determination of 3D structures of all proteins: - organize known proteins into families
More informationThe Schrödinger KNIME extensions
The Schrödinger KNIME extensions Computational Chemistry and Cheminformatics in a workflow environment Jean-Christophe Mozziconacci Volker Eyrich Topics What are the Schrödinger extensions? Workflow application
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 informationComputational Biology 1
Computational Biology 1 Protein Function & nzyme inetics Guna Rajagopal, Bioinformatics Institute, guna@bii.a-star.edu.sg References : Molecular Biology of the Cell, 4 th d. Alberts et. al. Pg. 129 190
More informationModeling for 3D structure prediction
Modeling for 3D structure prediction What is a predicted structure? A structure that is constructed using as the sole source of information data obtained from computer based data-mining. However, mixing
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 informationCMPS 3110: Bioinformatics. Tertiary Structure Prediction
CMPS 3110: Bioinformatics Tertiary Structure Prediction Tertiary Structure Prediction Why Should Tertiary Structure Prediction Be Possible? Molecules obey the laws of physics! Conformation space is finite
More informationCMPS 6630: Introduction to Computational Biology and Bioinformatics. Tertiary Structure Prediction
CMPS 6630: Introduction to Computational Biology and Bioinformatics Tertiary Structure Prediction Tertiary Structure Prediction Why Should Tertiary Structure Prediction Be Possible? Molecules obey the
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 informationSUPPLEMENTARY FIGURES
SUPPLEMENTARY FIGURES Supplementary Figure 1 Protein sequence alignment of Vibrionaceae with either a 40-residue insertion or a 44-residue insertion. Identical residues are indicated by red background.
More informationTHE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION
THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION AND CALIBRATION Calculation of turn and beta intrinsic propensities. A statistical analysis of a protein structure
More informationComputational chemical biology to address non-traditional drug targets. John Karanicolas
Computational chemical biology to address non-traditional drug targets John Karanicolas Our computational toolbox Structure-based approaches Ligand-based approaches Detailed MD simulations 2D fingerprints
More informationSupplementary Information. The protease GtgE from Salmonella exclusively targets. inactive Rab GTPases
Supplementary Information The protease GtgE from Salmonella exclusively targets inactive Rab GTPases Table of Contents Supplementary Figures... 2 Supplementary Figure 1... 2 Supplementary Figure 2... 3
More informationExperimental and Computational Mutagenesis to Investigate the. Positioning of a General Base within an Enzyme Active Site
Experimental and Computational Mutagenesis to Investigate the Positioning of a General Base within an Enzyme Active Site Jason P. Schwans, Philip Hanoian, Benjamin J. Lengerich, Fanny Sunden, Ana Gonzalez
More informationIntroduction to" Protein Structure
Introduction to" Protein Structure Function, evolution & experimental methods Thomas Blicher, Center for Biological Sequence Analysis Learning Objectives Outline the basic levels of protein structure.
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 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 informationBSc and MSc Degree Examinations
Examination Candidate Number: Desk Number: BSc and MSc Degree Examinations 2018-9 Department : BIOLOGY Title of Exam: Molecular Biology and Biochemistry Part I Time Allowed: 1 hour and 30 minutes Marking
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 informationProgress of Compound Library Design Using In-silico Approach for Collaborative Drug Discovery
21 th /June/2018@CUGM Progress of Compound Library Design Using In-silico Approach for Collaborative Drug Discovery Kaz Ikeda, Ph.D. Keio University Self Introduction Keio University, Tokyo, Japan (Established
More informationSignal Transduction Phosphorylation Protein kinases. Misfolding diseases. Protein Engineering Lysozyme variants
Signal Transduction Phosphorylation Protein kinases Misfolding diseases Protein Engineering Lysozyme variants Cells and Signals Regulation The cell must be able to respond to stimuli Cellular activities
More informationPhysiochemical Properties of Residues
Physiochemical Properties of Residues Various Sources C N Cα R Slide 1 Conformational Propensities Conformational Propensity is the frequency in which a residue adopts a given conformation (in a polypeptide)
More informationBioengineering & Bioinformatics Summer Institute, Dept. Computational Biology, University of Pittsburgh, PGH, PA
Pharmacophore Model Development for the Identification of Novel Acetylcholinesterase Inhibitors Edwin Kamau Dept Chem & Biochem Kennesa State Uni ersit Kennesa GA 30144 Dept. Chem. & Biochem. Kennesaw
More informationICM-Chemist-Pro How-To Guide. Version 3.6-1h Last Updated 12/29/2009
ICM-Chemist-Pro How-To Guide Version 3.6-1h Last Updated 12/29/2009 ICM-Chemist-Pro ICM 3D LIGAND EDITOR: SETUP 1. Read in a ligand molecule or PDB file. How to setup the ligand in the ICM 3D Ligand Editor.
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 informationProtein structure prediction. CS/CME/BioE/Biophys/BMI 279 Oct. 10 and 12, 2017 Ron Dror
Protein structure prediction CS/CME/BioE/Biophys/BMI 279 Oct. 10 and 12, 2017 Ron Dror 1 Outline Why predict protein structure? Can we use (pure) physics-based methods? Knowledge-based methods Two major
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 informationMD Simulation in Pose Refinement and Scoring Using AMBER Workflows
MD Simulation in Pose Refinement and Scoring Using AMBER Workflows Yuan Hu (On behalf of Merck D3R Team) D3R Grand Challenge 2 Webinar Department of Chemistry, Modeling & Informatics Merck Research Laboratories,
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 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 informationBiomolecules: lecture 10
Biomolecules: lecture 10 - understanding in detail how protein 3D structures form - realize that protein molecules are not static wire models but instead dynamic, where in principle every atom moves (yet
More informationONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole
ONETEP PB/SA: Application to G-Quadruplex DNA Stability Danny Cole Introduction Historical overview of structure and free energy calculation of complex molecules using molecular mechanics and continuum
More informationA conserved P-loop anchor limits the structural dynamics that mediate. nucleotide dissociation in EF-Tu.
Supplemental Material for A conserved P-loop anchor limits the structural dynamics that mediate nucleotide dissociation in EF-Tu. Evan Mercier 1,2, Dylan Girodat 1, and Hans-Joachim Wieden 1 * 1 Alberta
More informationSUPPLEMENTARY INFORMATION. doi: /nature07461
Figure S1 Electrophysiology. a ph-activation of. Two-electrode voltage clamp recordings of Xenopus oocytes expressing in comparison to waterinjected oocytes. Currents were recorded at 40 mv. The ph of
More informationDOCKING TUTORIAL. A. The docking Workflow
2 nd Strasbourg Summer School on Chemoinformatics VVF Obernai, France, 20-24 June 2010 E. Kellenberger DOCKING TUTORIAL A. The docking Workflow 1. Ligand preparation It consists in the standardization
More informationComputational Study of the DFG-Flip Conformational Transition in c Abl and c Src Tyrosine Kinases
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/jpcb Computational
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 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 information5.1. Hardwares, Softwares and Web server used in Molecular modeling
5. EXPERIMENTAL The tools, techniques and procedures/methods used for carrying out research work reported in this thesis have been described as follows: 5.1. Hardwares, Softwares and Web server used in
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 informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/309/5742/1868/dc1 Supporting Online Material for Toward High-Resolution de Novo Structure Prediction for Small Proteins Philip Bradley, Kira M. S. Misura, David Baker*
More informationOrientational degeneracy in the presence of one alignment tensor.
Orientational degeneracy in the presence of one alignment tensor. Rotation about the x, y and z axes can be performed in the aligned mode of the program to examine the four degenerate orientations of two
More informationPolypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation
Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation Jakob P. Ulmschneider and William L. Jorgensen J.A.C.S. 2004, 126, 1849-1857 Presented by Laura L. Thomas and
More informationIntroduction to Computational Structural Biology
Introduction to Computational Structural Biology Part I 1. Introduction The disciplinary character of Computational Structural Biology The mathematical background required and the topics covered Bibliography
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 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 informationALL LECTURES IN SB Introduction
1. Introduction 2. Molecular Architecture I 3. Molecular Architecture II 4. Molecular Simulation I 5. Molecular Simulation II 6. Bioinformatics I 7. Bioinformatics II 8. Prediction I 9. Prediction II ALL
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 informationEuropean Life Science Bootcamp: Case Studies in Enhanced Sampling Methods.
European Life Science Bootcamp: Case Studies in Enhanced Sampling Methods davide.branduardi@schrodinger.com Outline Limitations of molecular dynamics Overcoming them with Enhanced sampling Pragmatic use
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 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 informationProtein Structure Prediction II Lecturer: Serafim Batzoglou Scribe: Samy Hamdouche
Protein Structure Prediction II Lecturer: Serafim Batzoglou Scribe: Samy Hamdouche The molecular structure of a protein can be broken down hierarchically. The primary structure of a protein is simply its
More informationSupplementary Figures
1 Supplementary Figures Supplementary Figure 1 Type I FGFR1 inhibitors (a) Chemical structures of a pyrazolylaminopyrimidine inhibitor (henceforth referred to as PAPI; PDB-code of the FGFR1-PAPI complex:
More informationSupplementary Materials for
www.sciencesignaling.org/cgi/content/full/5/243/ra68/dc1 Supplementary Materials for Superbinder SH2 Domains Act as Antagonists of Cell Signaling Tomonori Kaneko, Haiming Huang, Xuan Cao, Xing Li, Chengjun
More informationAn Introduction to Metabolism
An Introduction to Metabolism Chapter 8 Objectives Distinguish between the following pairs of terms: catabolic and anabolic pathways; kinetic and potential energy; open and closed systems; exergonic and
More informationProtein Structure Prediction
Page 1 Protein Structure Prediction Russ B. Altman BMI 214 CS 274 Protein Folding is different from structure prediction --Folding is concerned with the process of taking the 3D shape, usually based on
More informationConformational Searching using MacroModel and ConfGen. John Shelley Schrödinger Fellow
Conformational Searching using MacroModel and ConfGen John Shelley Schrödinger Fellow Overview Types of conformational searching applications MacroModel s conformation generation procedure General features
More informationUseful background reading
Overview of lecture * General comment on peptide bond * Discussion of backbone dihedral angles * Discussion of Ramachandran plots * Description of helix types. * Description of structures * NMR patterns
More informationComputational Study of Gleevec and G6G Reveals Molecular Determinants of Kinase Inhibitor Selectivity
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/jacs Computational
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 informationL718Q mutant EGFR escapes covalent inhibition by stabilizing. a non-reactive conformation of the lung cancer drug. osimertinib
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2018 Electronic Supplementary Information (ESI) for L718Q mutant EGFR escapes covalent inhibition
More information4 Examples of enzymes
Catalysis 1 4 Examples of enzymes Adding water to a substrate: Serine proteases. Carbonic anhydrase. Restrictions Endonuclease. Transfer of a Phosphoryl group from ATP to a nucleotide. Nucleoside monophosphate
More informationNMR Structure of the HWE Kinase Associated Response Regulator Sma0114 in Its Activated State
pubs.acs.org/biochemistry NMR Structure of the HWE Kinase Associated Response Regulator Sma0114 in Its Activated State Sarah R. Sheftic, Emma White, Daniel J. Gage, and Andrei T. Alexandrescu* Department
More informationBuild_model v User Guide
Build_model v.2.0.1 User Guide MolTech Build_model User Guide 2008-2011 Molecular Technologies Ltd. www.moltech.ru Please send your comments and suggestions to contact@moltech.ru. Table of Contents Input
More informationBiophysics 490M Project
Biophysics 490M Project Dan Han Department of Biochemistry Structure Exploration of aa 3 -type Cytochrome c Oxidase from Rhodobacter sphaeroides I. Introduction: All organisms need energy to live. They
More informationSchrodinger ebootcamp #3, Summer EXPLORING METHODS FOR CONFORMER SEARCHING Jas Bhachoo, Senior Applications Scientist
Schrodinger ebootcamp #3, Summer 2016 EXPLORING METHODS FOR CONFORMER SEARCHING Jas Bhachoo, Senior Applications Scientist Numerous applications Generating conformations MM Agenda http://www.schrodinger.com/macromodel
More informationschematic diagram; EGF binding, dimerization, phosphorylation, Grb2 binding, etc.
Lecture 1: Noncovalent Biomolecular Interactions Bioengineering and Modeling of biological processes -e.g. tissue engineering, cancer, autoimmune disease Example: RTK signaling, e.g. EGFR Growth responses
More informationβ1 Structure Prediction and Validation
13 Chapter 2 β1 Structure Prediction and Validation 2.1 Overview Over several years, GPCR prediction methods in the Goddard lab have evolved to keep pace with the changing field of GPCR structure. Despite
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Structure of human carbamoyl phosphate synthetase: deciphering the on/off switch of human ureagenesis Sergio de Cima, Luis M. Polo, Carmen Díez-Fernández, Ana I. Martínez, Javier
More informationCAP 5510 Lecture 3 Protein Structures
CAP 5510 Lecture 3 Protein Structures Su-Shing Chen Bioinformatics CISE 8/19/2005 Su-Shing Chen, CISE 1 Protein Conformation 8/19/2005 Su-Shing Chen, CISE 2 Protein Conformational Structures Hydrophobicity
More informationProtein Structure. W. M. Grogan, Ph.D. OBJECTIVES
Protein Structure W. M. Grogan, Ph.D. OBJECTIVES 1. Describe the structure and characteristic properties of typical proteins. 2. List and describe the four levels of structure found in proteins. 3. Relate
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature11054 Supplementary Fig. 1 Sequence alignment of Na v Rh with NaChBac, Na v Ab, and eukaryotic Na v and Ca v homologs. Secondary structural elements of Na v Rh are indicated above the
More informationReceptor Based Drug Design (1)
Induced Fit Model For more than 100 years, the behaviour of enzymes had been explained by the "lock-and-key" mechanism developed by pioneering German chemist Emil Fischer. Fischer thought that the chemicals
More informationMetabolism and Enzymes
Energy Basics Metabolism and Enzymes Chapter 5 Pgs. 77 86 Chapter 8 Pgs. 142 162 Energy is the capacity to cause change, and is required to do work. Very difficult to define quantity. Two types of energy:
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationLecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions
More informationSUPPLEMENTARY INFORMATION
doi:1.138/nature1737 Supplementary Table 1 variant Description FSEC - 2B12 a FSEC - 6A1 a K d (leucine) c Leucine uptake e K (wild-type like) K (Y18F) K (TS) K (TSY) K288A mutant, lipid facing side chain
More informationProgramme Last week s quiz results + Summary Fold recognition Break Exercise: Modelling remote homologues
Programme 8.00-8.20 Last week s quiz results + Summary 8.20-9.00 Fold recognition 9.00-9.15 Break 9.15-11.20 Exercise: Modelling remote homologues 11.20-11.40 Summary & discussion 11.40-12.00 Quiz 1 Feedback
More informationProtein Structure Prediction and Protein-Ligand Docking
Protein Structure Prediction and Protein-Ligand Docking Björn Wallner bjornw@ifm.liu.se Jan. 24, 2014 Todays topics Protein Folding Intro Protein structure prediction How can we predict the structure of
More informationThe Structure and Functions of Proteins
Wright State University CORE Scholar Computer Science and Engineering Faculty Publications Computer Science and Engineering 2003 The Structure and Functions of Proteins Dan E. Krane Wright State University
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature11744 Supplementary Table 1. Crystallographic data collection and refinement statistics. Wild-type Se-Met-BcsA-B SmCl 3 -soaked EMTS-soaked Data collection Space
More informationStructural insight into the binding of C60-derivatives with enoyl-pyruvate transferase from Helicobacter pylori
www.bioinformation.net Volume 13(6) Hypothesis Structural insight into the binding of C60-derivatives with enoyl-pyruvate transferase from Helicobacter pylori Mohammad Teimouri 1,*, Muhammad Junaid 2,
More informationTable 1. Crystallographic data collection, phasing and refinement statistics. Native Hg soaked Mn soaked 1 Mn soaked 2
Table 1. Crystallographic data collection, phasing and refinement statistics Native Hg soaked Mn soaked 1 Mn soaked 2 Data collection Space group P2 1 2 1 2 1 P2 1 2 1 2 1 P2 1 2 1 2 1 P2 1 2 1 2 1 Cell
More information