DOCSLIB.ORG

MODELLER 9V7 Manual

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
MODELLER 9V7 Manual MODELLER A Program for Protein Structure Modeling Release 9v7, r6923 Andrej Saliˇ with help from Ben Webb, M.S. Madhusudhan, Min-Yi Shen, Marc A. Martı-Renom, Narayanan Eswar, Frank Alber, Maya Topf, Baldomero Oliva, Andr´as Fiser, Roberto S´anchez, Bozidar Yerkovich, Azat Badretdinov, Francisco Melo, John P. Overington, and Eric Feyfant email: modeller-care AT salilab.org URL http://salilab.org/modeller/ 2009/06/12 ii Contents Copyright notice xix Acknowledgments xxi 1 Introduction 1 1.1 What is Modeller?............................................. 1 1.2 Modeller bibliography.......................................... 2 1.3 Obtaining and installing the program . ........... 3 1.4 Bugreports........................................ .......... 4 1.5 Method for comparative protein structure modeling by Modeller ................... 5 1.6 Using Modeller forcomparativemodeling................................ 8 1.6.1 Preparing input files . ..... 8 1.6.2 Running Modeller ......................................... 9 2 Automated comparative modeling with automodel 11 2.1 Simpleusage ........................................ ......... 11 2.2 More advanced usage . ......... 12 2.2.1 Including water molecules, HETATM residues, and hydrogen atoms . ...... 12 2.2.2 Changing the default optimization and refinement protocol . ...... 14 2.2.3 Getting a very fast and approximate model . ......... 14 2.2.4 Building a model from multiple templates . ........ 15 2.2.5 Building an all hydrogen model . ........ 16 2.2.6 Refining only part of the model . ..... 16 2.2.7 Including disulfide bridges . 18 2.2.8 Providing your own restraints file . ............ 19 2.2.9 Using your own initial model . ......... 19 2.2.10 Adding additional restraints to the defaults . ........... 20 2.2.11 Building multi-chain models with symmetry . ........... 21 2.2.12 Accessing output data after modeling is complete . ....... 22 2.2.13 Fully automated alignment and modeling . ......... 23 2.3 Loopoptimization .................................. ............ 24 2.3.1 Automatic loop refinement after model building . ....... 24 2.3.2 Defining loop regions for refinement . ..... 25 2.3.3 Refining an existing PDB file . ...... 26 iii iv CONTENTS 3 Frequently asked questions and history 27 3.1 Frequently asked questions (FAQ) and examples . ......... 27 3.2 Modeller updates ............................................. 34 3.2.1 Changes since release 9v6 . ....... 34 3.2.2 Changes since release 9v5 . ....... 35 4 Comparative modeling class reference 37 4.1 automodel reference . 37 4.1.1 automodel() — prepare to build one or more comparative models . ......... 37 4.1.2 automodel.library schedule — select optimization schedule . 38 4.1.3 automodel.md level — control the model refinement level . 38 4.1.4 automodel.outputs — all output data for generated models . ...... 38 4.1.5 automodel.rand method — control initial model randomization . 38 4.1.6 automodel.generate method — control initial model generation . 38 4.1.7 automodel.max var iterations — select length of optimizations . 38 4.1.8 automodel.repeat optimization — number of times to repeat optimization . 38 4.1.9 automodel.max molpdf — objective function cutoff . 39 4.1.10 automodel.initial malign3d — initial template alignment . 39 4.1.11 automodel.starting model — first model to build . 39 4.1.12 automodel.ending model — last model to build . 39 4.1.13 automodel.final malign3d — final template-model alignment . 39 4.1.14 automodel.write intermediates — write intermediate files during optimization . 39 4.1.15 automodel.trace output — control optimization output . 39 4.1.16 automodel.get optimize actions() — get actions to carry out during the initial optimization . 39 4.1.17 automodel.get refine actions() — get actions to carry out during the refinement . 40 4.1.18 automodel.select atoms() — select region for optimization . 40 4.1.19 automodel.auto align() — generate an automatic initial alignment . 40 4.1.20 automodel.very fast() — request rapid optimization . 40 4.1.21 automodel.make() — build all models . ......... 40 4.1.22 automodel.cluster() — cluster all built models . ......... 40 4.1.23 automodel.special restraints() — add additional restraints . 41 4.1.24 automodel.nonstd restraints() — add restraints on ligands . 41 4.1.25 automodel.special patches() — add additional patches to the topology . 41 4.1.26 automodel.user after single model() — analyze or refine each model . 41 4.1.27 automodel.get model filename() — get the model PDB name . 41 4.1.28 automodel.use parallel job() — parallelize model building . 42 4.2 allhmodel reference . 42 4.2.1 allhmodel() — prepare to build all-hydrogen models . ........ 42 4.3 loopmodel reference . 42 4.3.1 loopmodel() — prepare to build models with loop refinement . ...... 42 4.3.2 loopmodel.loop.md level — control the loop model refinement level . 42 4.3.3 loopmodel.loop.max var iterations — select length of optimizations . 43 4.3.4 loopmodel.loop.library schedule — select optimization schedule . 43 4.3.5 loopmodel.loop.starting model — first loop model to build . 43 CONTENTS v 4.3.6 loopmodel.loop.ending model — last loop model to build . 43 4.3.7 loopmodel.loop.outputs — all output data for generated loop models . ......... 43 4.3.8 loopmodel.select loop atoms() — select region for loop optimization . 43 4.3.9 loopmodel.get loop model filename() — get the model PDB name . 43 4.3.10 loopmodel.user after single loop model() — analyze or refine each loop model . 43 4.3.11 loopmodel.read potential() — read in the loop modeling potential . 44 4.3.12 loopmodel.build ini loop() — create the initial conformation of the loop . 44 4.4 dope loopmodel reference.......................................... 44 4.4.1 dope loopmodel() — prepare to build models with DOPE loop refinement . 44 4.5 dopehr loopmodel reference......................................... 44 5 Modeller general reference 45 5.1 Miscellaneous rules and features of Modeller .............................. 45 5.1.1 Modeller system.......................................... 45 5.1.2 Controlling breakpoints and the amount of output . ........... 45 5.1.3 Filenaming...................................... ........ 45 5.1.4 Filetypes ...................................... ......... 46 5.2 Stereochemical parameters and molecular topology . ........... 46 5.2.1 Modeling residues with non-existing or incomplete entries in the topology and parameter libraries 47 5.3 Spatialrestraints .................................. ............. 47 5.3.1 Specification of restraints . .......... 48 5.3.2 Specification of pseudo atoms . ...... 50 5.3.3 Excluded pairs . ...... 52 5.3.4 Rigidbodies ..................................... ........ 52 5.3.5 Symmetryrestraints ............................... .......... 53 6 Modeller command reference 57 6.1 Key for command descriptions . ....... 57 6.2 The environ class: Modeller environment ............................... 57 6.2.1 environ() — create a new Modeller environment........................ 57 6.2.2 environ.io — default input parameters . ....... 58 6.2.3 environ.edat — default objective function parameters . ........ 58 6.2.4 environ.libs — Modeller libraries................................. 58 6.2.5 environ.schedule scale — energy function scaling factors . 58 6.2.6 environ.dendrogram() — clustering . ........ 58 6.2.7 environ.principal components() — clustering . 58 6.2.8 environ.system() — execute system command . ........ 59 6.2.9 environ.make pssmdb() — Create a database of PSSMs given a list of profiles . 59 6.3 The energy data class: objective function parameters . 61 6.3.1 energy data() — create a new set of objective function parameters . 61 6.3.2 energy data.contact shell — nonbond distance cutoff . 61 6.3.3 energy data.update dynamic — nonbond recalculation threshold . 61 6.3.4 energy data.sphere stdv — soft-sphere standard deviation . 62 6.3.5 energy data.dynamic sphere — calculate soft-sphere overlap restraints . 62 vi CONTENTS 6.3.6 energy data.dynamic lennard — calculate Lennard-Jones restraints . 62 6.3.7 energy data.dynamic coulomb — calculate Coulomb restraints . 62 6.3.8 energy data.dynamic modeller — calculate non-bonded spline restraints . 62 6.3.9 energy data.excl local — exclude certain local pairs of atoms . 62 6.3.10 energy data.radii factor—scaleatomicradii ........................... 62 6.3.11 energy data.lennard jones switch — Lennard-Jones switching parameters . 62 6.3.12 energy data.coulomb switch — Coulomb switching parameters . 63 6.3.13 energy data.relative dielectric — relative dielectric . 63 6.3.14 energy data.covalent cys — use disulfide bridges in residue distance . 63 6.3.15 energy data.nonbonded sel atoms — control interaction with picked atoms . 63 6.3.16 energy data.nlogn use — select non-bond list generation algorithm . 63 6.3.17 energy data.energy terms — user-defined global energy terms . 63 6.4 The io data class: coordinate file input parameters . 64 6.4.1 io data() — create a new input parameters object . 64 6.4.2 io data.hetatm — whether to read HETATM records . 64 6.4.3 io data.hydrogen — whether to read hydrogen atoms . 64 6.4.4 io data.water — whether to read water molecules . 64 6.4.5 io data.atom files directory — search path for coordinate files . 64 6.5 The Libraries class: stereochemical parameters and molecular topology . 65 6.5.1 Libraries.topology — topology library information . ................... 65 6.5.2 Libraries.parameters — parameter library information . .............. 65 6.5.3 Topology.append() — append residue topology library . ........... 65 6.5.4 Topology.clear()
Load more

Recommended publications
  • MODELLER 10.1 Manual
    MODELLER A Program for Protein Structure Modeling Release 10.1, r12156 Andrej Saliˇ with help from Ben Webb, M.S. Madhusudhan, Min-Yi Shen, Guangqiang Dong, Marc A. Martı-Renom, Narayanan Eswar, Frank Alber, Maya Topf, Baldomero Oliva, Andr´as Fiser, Roberto S´anchez, Bozidar Yerkovich, Azat Badretdinov, Francisco Melo, John P. Overington, and Eric Feyfant email: modeller-care AT salilab.org URL https://salilab.org/modeller/ 2021/03/12 ii Contents Copyright notice xxi Acknowledgments xxv 1 Introduction 1 1.1 What is Modeller?............................................. 1 1.2 Modeller bibliography....................................... .... 2 1.3 Obtainingandinstallingtheprogram. .................... 3 1.4 Bugreports...................................... ............ 4 1.5 Method for comparative protein structure modeling by Modeller ................... 5 1.6 Using Modeller forcomparativemodeling. ... 8 1.6.1 Preparinginputfiles . ............. 8 1.6.2 Running Modeller ......................................... 9 2 Automated comparative modeling with AutoModel 11 2.1 Simpleusage ..................................... ............ 11 2.2 Moreadvancedusage............................... .............. 12 2.2.1 Including water molecules, HETATM residues, and hydrogenatoms .............. 12 2.2.2 Changing the default optimization and refinement protocol ................... 14 2.2.3 Getting a very fast and approximate model . ................. 14 2.2.4 Building a model from multiple templates . .................. 15 2.2.5 Buildinganallhydrogenmodel
    [Show full text]
  • MODELLER 10.0 Manual
    MODELLER A Program for Protein Structure Modeling Release 10.0, r12005 Andrej Saliˇ with help from Ben Webb, M.S. Madhusudhan, Min-Yi Shen, Guangqiang Dong, Marc A. Martı-Renom, Narayanan Eswar, Frank Alber, Maya Topf, Baldomero Oliva, Andr´as Fiser, Roberto S´anchez, Bozidar Yerkovich, Azat Badretdinov, Francisco Melo, John P. Overington, and Eric Feyfant email: modeller-care AT salilab.org URL https://salilab.org/modeller/ 2021/01/25 ii Contents Copyright notice xxi Acknowledgments xxv 1 Introduction 1 1.1 What is Modeller?............................................. 1 1.2 Modeller bibliography....................................... .... 2 1.3 Obtainingandinstallingtheprogram. .................... 3 1.4 Bugreports...................................... ............ 4 1.5 Method for comparative protein structure modeling by Modeller ................... 5 1.6 Using Modeller forcomparativemodeling. ... 8 1.6.1 Preparinginputfiles . ............. 8 1.6.2 Running Modeller ......................................... 9 2 Automated comparative modeling with AutoModel 11 2.1 Simpleusage ..................................... ............ 11 2.2 Moreadvancedusage............................... .............. 12 2.2.1 Including water molecules, HETATM residues, and hydrogenatoms .............. 12 2.2.2 Changing the default optimization and refinement protocol ................... 14 2.2.3 Getting a very fast and approximate model . ................. 14 2.2.4 Building a model from multiple templates . .................. 15 2.2.5 Buildinganallhydrogenmodel
    [Show full text]
  • Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 11 June 2019 doi:10.20944/preprints201906.0091.v1 Peer-reviewed version available at Biomolecules 2019, 9; doi:10.3390/biom9040123 Molecular Docking Studies of a Cyclic Octapeptide-Cyclosaplin from Sandalwood Abheepsa Mishra1, 2,* and Satyahari Dey1 1Plant Biotechnology Laboratory, Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur-721302, West Bengal, India; [email protected] (A.M.); [email protected] (S.D.) 2Department of Internal Medicine, The University of Texas Southwestern Medical Center, 5323 Harry Hines Blvd, Dallas, TX 75390, USA *Correspondence: [email protected]; Tel.:+1(518) 881-9196 Abstract Natural products from plants such as, chemopreventive agents attract huge attention because of their low toxicity and high specificity. The rational drug design in combination with structure based modeling and rapid screening methods offer significant potential for identifying and developing lead anticancer molecules. Thus, the molecular docking method plays an important role in screening a large set of molecules based on their free binding energies and proposes structural hypotheses of how the molecules can inhibit the target. Several peptide based therapeutics have been developed to combat several health disorders including cancers, metabolic disorders, heart-related, and infectious diseases. Despite the discovery of hundreds of such therapeutic peptides however, only few peptide-based drugs have made it to the market. Moreover, until date the activities of cyclic peptides towards molecular targets such as protein kinases, proteases, and apoptosis related proteins have never been explored. In this study we explore the in silico kinase and protease inhibitor potentials of cyclosaplin as well as study the interactions of cyclosaplin with other cancer-related proteins.
    [Show full text]
  • New Computational Approaches for Investigating the Impact of Mutations on the Transglucosylation Activity of Sucrose Phosphorylase Enzyme
    Thesis submitted to the Université de la Réunion for award of Doctor of Philosophy in Sciences speciality in bioinformatics New computational approaches for investigating the impact of mutations on the transglucosylation activity of sucrose phosphorylase enzyme Mahesh VELUSAMY Thesis to be presented on 18th December 2018 in front of the jury composed of Mme Isabelle ANDRE, Directeur de Recherche CNRS, INSA de TOULOUSE, Rapporteur M. Manuel DAUCHEZ, Professeur, Université de Reims Champagne Ardennes, Rapporteur M. Richard DANIELLOU, Professeur, Université d'Orléans, Examinateur M. Yves-Henri SANEJOUAND, Directeur de Recherche CNRS, Université de Nantes, Examinateur Mme Irène MAFFUCCI, Maître de Conférences, Université de Technologie de Compiègne, Examinateur Mme Corinne MIRAL, Maître de Conférences HDR, Université de Nantes, Examinateur M. Frédéric CADET, Professeur, Université de La Réunion, Co-directeur de thèse M. Bernard OFFMANN, Professeur, Université de Nantes, Directeur de thèse ெ்்ந்ி ிைு்ூுத் ெ்யாம் ெ்த உதி்ு ையகு் ானகு் ஆ்ற் அிு. -ிு்ு் ுதி் அ்ப் ுுக், எனு த்ைத ேுாி, தா் க்ூி, அ்ண் ு்ு்ுமா், த்ைக பா்பா, ஆ்தா ுு்மா், ீனா, அ்ண் ்ுக், ெபிய்பா, ெபிய்மா ம்ு் உுுைணயா் இு்த அைண்ு ந்ப்கு்ு் எனு மனமா்்த ந்ி. இ்த ஆ்ி்ைக ுுைம அை3த்ு ுுுத்்காரண், எனு ஆ்ி்ைக இய்ுன் ேபராிிய் ெப்னா்் ஆஃ்ேம். ப்ேு துண்கி் நா் மனதாு், ெபாுளாதார அளிு் க்3்ி் இு்தேபாு, என்ு இ்ெனாு த்ைதயாகே இு்ு எ்ைன பா்்ு்ெகா்3ா். ுி்பாக, எனு ூ்றா் ஆ்ு இுிி், அ்்ு எ்ளோ தி்ப்3 க3ைமக் ம்ு் ிர்ிைனக் இு்தாு், அைத்ெபாு்பு்தாு, அ் என்ு ெ்த ெபாுளாதார உதி, ப்கைB்கழக பிு ம்ு் இதர ி்ாக ்ப்த்ப்3 உதிகு்ு எ்ன ைகமா்ு ெகாு்தாு் ஈ3ாகாு.
    [Show full text]
  • University of Huddersfield Repository
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by University of Huddersfield Repository University of Huddersfield Repository Liang, Shuo, Holmes, Violeta and Kureshi, Ibad Hybrid Computer Cluster with High Flexibility Original Citation Liang, Shuo, Holmes, Violeta and Kureshi, Ibad (2012) Hybrid Computer Cluster with High Flexibility. In: IEEE Cluster 2012, 24-28 September 2012, Beijing, China. This version is available at http://eprints.hud.ac.uk/15840/ The University Repository is a digital collection of the research output of the University, available on Open Access. Copyright and Moral Rights for the items on this site are retained by the individual author and/or other copyright owners. Users may access full items free of charge; copies of full text items generally can be reproduced, displayed or performed and given to third parties in any format or medium for personal research or study, educational or not-for-profit purposes without prior permission or charge, provided: • The authors, title and full bibliographic details is credited in any copy; • A hyperlink and/or URL is included for the original metadata page; and • The content is not changed in any way. For more information, including our policy and submission procedure, please contact the Repository Team at: [email protected]. http://eprints.hud.ac.uk/ Hybrid Computer Cluster with High Flexibility Shuo Liang Dr Violeta Holmes Ibad Kureshi School of Computing and Engineering School of Computing and Engineering School of Computing and Engineering University of Huddersfield University of Huddersfield University of Huddersfield Huddersfield, United Kingdom Huddersfield, United Kingdom Huddersfield, United Kingdom Email: [email protected] Email: [email protected] Email: [email protected] Abstract—In this paper we present a cluster middleware, the others support multi-platform, supplying different user designed to implement a Linux-Windows Hybrid HPC Cluster, interfaces and document support.
    [Show full text]
  • Graduate Brochure 2014 Draft.Pub
    THETHE UNIVERSITYUNIVERSITY OFOF MISSOURIMISSOURI——ST.ST. LOUISLOUIS 315 Benton Hall (MC 27) One University Blvd. St. Louis, MO 63121-4400 http://www.umsl.edu/chemistry/ DepartmentDepartment ofof ChemistryChemistry && BiochemistryBiochemistry Phone: (314) 516-5311 Fax: (314) 516-5342 ST.LOUIS The St. Louis metropolitan region, located around the ence Center has many educational interactive exhibits confluence of the Missouri and Mississippi Rivers, has including an Omnimax theatre. And finally the St. Louis evolved through the centuries from the homes of ancient Zoo, home to more than 18,000 exotic animals, many of Native American civilizations to major fur trading cen- them rare and endangered, and set in the rolling hills, ter to the "Gateway to the West," marked by lasting lakes and glades of Forest Park, is a great place to visit. French, Spanish and African influences. Today the area is a town with hospitable people rooted in a myriad of Activities and attractions are many, and St. Louis has regional, ethnic, and cultural traditions reflective of our a great variety of night life choices. Blues clubs and complex world, well supported by easy access to parks, restaurants are tucked away in the red brick buildings of educational centers, sports venues, museums, and histor- the historic Soulard neighborhood, while Dixieland and ic sites. blues dinner cruises are available through the port of St. Louis. Clubs and restaurants are alive until early morn- Great ethnic and classic neighborhoods characterize ing hours in the converted warehouses of the Landing, the region. A cross-section of the area can provide ex- north of the Arch, and along the newly-developed "loft" amples of wonderful Victorian architecture, museum, neighborhood of Washington Avenue.
    [Show full text]
  • Structural Bioinformatics
    Structural Bioinformatics Sanne Abeln K. Anton Feenstra Centre for Integrative Bioinformatics (IBIVU), and Department of Computer Science, Vrije Universiteit, De Boelelaan 1081A, 1081 HV Amsterdam, Netherlands December 4, 2017 arXiv:1712.00425v1 [q-bio.BM] 1 Dec 2017 2 Abstract This chapter deals with approaches for protein three-dimensional structure prediction, starting out from a single input sequence with unknown struc- ture, the `query' or `target' sequence. Both template based and template free modelling techniques are treated, and how resulting structural models may be selected and refined. We give a concrete flowchart for how to de- cide which modelling strategy is best suited in particular circumstances, and which steps need to be taken in each strategy. Notably, the ability to locate a suitable structural template by homology or fold recognition is crucial; without this models will be of low quality at best. With a template avail- able, the quality of the query-template alignment crucially determines the model quality. We also discuss how other, courser, experimental data may be incorporated in the modelling process to alleviate the problem of missing template structures. Finally, we discuss measures to predict the quality of models generated. Structural Bioinformatics c Abeln & Feenstra, 2014-2017 Contents 7 Strategies for protein structure model generation5 7.1 Template based protein structure modelling . .7 7.1.1 Homology based Template Finding . .7 7.1.2 Fold recognition . .8 7.1.3 Generating the target-template alignment . 10 7.1.4 Generating a model . 11 7.1.5 Loop or missing substructure modelling . 12 7.2 Template-free protein structure modelling .
    [Show full text]
  • Biodynamo: a New Platform for Large-Scale Biological Simulation
    Lukas Johannes Breitwieser, BSc BioDynaMo: A New Platform for Large-Scale Biological Simulation MASTER’S THESIS to achieve the university degree of Diplom-Ingenieur Master’s degree programme Software Engineering and Management submitted to Graz University of Technology Supervisor Ass.Prof. Dipl.-Ing. Dr.techn. Christian Steger Institute for Technical Informatics Advisor Ass.Prof. Dipl.-Ing. Dr.techn. Christian Steger Dr. Fons Rademakers (CERN) November 23, 2016 AFFIDAVIT I declare that I have authored this thesis independently, that I have not used other than the declared sources/resources, and that I have explicitly indicated all ma- terial which has been quoted either literally or by content from the sources used. The text document uploaded to TUGRAZonline is identical to the present master‘s thesis. Date Signature This thesis is dedicated to my parents Monika and Fritz Breitwieser Without their endless love, support and encouragement I would not be where I am today. iii Acknowledgment This thesis would not have been possible without the support of a number of people. I want to thank my supervisor Christian Steger for his guidance, support, and input throughout the creation of this document. Furthermore, I would like to express my deep gratitude towards the entire CERN openlab team, Alberto Di Meglio, Marco Manca and Fons Rademakers for giving me the opportunity to work on this intriguing project. Furthermore, I would like to thank them for vivid discussions, their valuable feedback, scientific freedom and possibility to grow. My gratitude is extended to Roman Bauer from Newcastle University (UK) for his helpful insights about computational biology detailed explanations and valuable feedback.
    [Show full text]
  • CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations
    molecules Article CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations 1,2,3, 2,3, 1, 1 Alessandro Marchetto y , Zeineb Si Chaib y , Carlo Alberto Rossi y, Rui Ribeiro , Sergio Pantano 4, Giulia Rossetti 2,5,6,* and Alejandro Giorgetti 1,2,* 1 Department of Biotechnology, University of Verona, 37134 Verona, Italy; [email protected] (A.M.); [email protected] (C.A.R.); [email protected] (R.R.) 2 Computational Biomedicine, Institute for Neuroscience and Medicine (INM-9) and Institute for Advanced Simulations (IAS-5), Forschungszentrum Jülich, 52425 Jülich, Germany; [email protected] 3 Faculty of Mathematics, Computer Science and Natural Sciences, RWTH Aachen University,52062 Aachen, Germany 4 Institut Pasteur de Montevideo, Mataojo 2020, 11400 Montevideo, Uruguay; [email protected] 5 Jülich Supercomputing Center (JSC), Forschungszentrum Jülich, 52428 Jülich, Germany 6 Department of Hematology,Oncology,Hemostaseology,and Stem Cell Transplantation, RWTH Aachen University, 52062 Aachen, Germany * Correspondence: [email protected] (G.R.); [email protected] (A.G.) These authors contributed equally to this work. y Academic Editors: Carmen Domene and Mercedes Alfonso-Prieto Received: 30 October 2020; Accepted: 28 November 2020; Published: 15 December 2020 Abstract: Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations.
    [Show full text]
  • Methods for the Refinement of Protein Structure 3D Models
    International Journal of Molecular Sciences Review Methods for the Refinement of Protein Structure 3D Models Recep Adiyaman and Liam James McGuffin * School of Biological Sciences, University of Reading, Reading RG6 6AS, UK; [email protected] * Correspondence: l.j.mcguffi[email protected]; Tel.: +44-0-118-378-6332 Received: 2 April 2019; Accepted: 7 May 2019; Published: 1 May 2019 Abstract: The refinement of predicted 3D protein models is crucial in bringing them closer towards experimental accuracy for further computational studies. Refinement approaches can be divided into two main stages: The sampling and scoring stages. Sampling strategies, such as the popular Molecular Dynamics (MD)-based protocols, aim to generate improved 3D models. However, generating 3D models that are closer to the native structure than the initial model remains challenging, as structural deviations from the native basin can be encountered due to force-field inaccuracies. Therefore, different restraint strategies have been applied in order to avoid deviations away from the native structure. For example, the accurate prediction of local errors and/or contacts in the initial models can be used to guide restraints. MD-based protocols, using physics-based force fields and smart restraints, have made significant progress towards a more consistent refinement of 3D models. The scoring stage, including energy functions and Model Quality Assessment Programs (MQAPs) are also used to discriminate near-native conformations from non-native conformations. Nevertheless, there are often very small differences among generated 3D models in refinement pipelines, which makes model discrimination and selection problematic. For this reason, the identification of the most native-like conformations remains a major challenge.
    [Show full text]
  • Exercise 6: Homology Modeling
    Workshop in Computational Structural Biology (81813), Spring 2019 Exercise 6: Homology Modeling Homology modeling, sometimes referred to as comparative modeling, is a commonly-used alternative to ab-initio modeling (which requires zero knowledge about the native structure, but also overlooks much of our understanding of structural diversity). Homology modeling leverages the massive amount of structural knowledge available today to guide our search of the native structure. Practically, it assumes that sequence similarity implies structural similarity, i.e. assumes a protein homologous to our query protein cannot be very different from it. Therefore, we can roughly "thread" the query sequence onto the template (homologue) structure, and exploit existing knowledge to guide our modeling. Contents Homology modeling with Rosetta 1 Overview 1 Template Selection 1 Loop Modeling 2 Running the full simulation 4 Bonus section! Homology modeling with automated/pre-calculated servers 6 MODBASE 6 I-Tasser 6 Homology modeling with Rosetta Overview In this section, we will model the structure of the Erbin PDZ domain, using homologue template structures, and analyze the quality of the models. The structure of this protein has been solved (PDB ID 1MFG), so we will be able to compare prediction with experiment. As taught in class, classical homology modeling in Rosetta is composed of four phases: 1. Template selection 2. Template-query alignment 3. Rebuilding loop regions 4. Full-atom refinement of the model Our model protein is the Erbin PDZ domain: Erbin contains a class-I PDZ domain that binds to the C-terminal region of the receptor tyrosine kinase ErbB2. PDB ID 1MFG is the crystal structure of the human Erbin PDZ bound to the peptide EYLGLDVPV corresponding to the C-terminal residues of human ErbB2.
    [Show full text]
  • Advances in Rosetta Protein Structure Prediction on Massively Parallel Systems
    UC San Diego UC San Diego Previously Published Works Title Advances in Rosetta protein structure prediction on massively parallel systems Permalink https://escholarship.org/uc/item/87g6q6bw Journal IBM Journal of Research and Development, 52(1) ISSN 0018-8646 Authors Raman, S. Baker, D. Qian, B. et al. Publication Date 2008 Peer reviewed eScholarship.org Powered by the California Digital Library University of California Advances in Rosetta protein S. Raman B. Qian structure prediction on D. Baker massively parallel systems R. C. Walker One of the key challenges in computational biology is prediction of three-dimensional protein structures from amino-acid sequences. For most proteins, the ‘‘native state’’ lies at the bottom of a free- energy landscape. Protein structure prediction involves varying the degrees of freedom of the protein in a constrained manner until it approaches its native state. In the Rosetta protein structure prediction protocols, a large number of independent folding trajectories are simulated, and several lowest-energy results are likely to be close to the native state. The availability of hundred-teraflop, and shortly, petaflop, computing resources is revolutionizing the approaches available for protein structure prediction. Here, we discuss issues involved in utilizing such machines efficiently with the Rosetta code, including an overview of recent results of the Critical Assessment of Techniques for Protein Structure Prediction 7 (CASP7) in which the computationally demanding structure-refinement process was run on 16 racks of the IBM Blue Gene/Le system at the IBM T. J. Watson Research Center. We highlight recent advances in high-performance computing and discuss future development paths that make use of the next-generation petascale (.1012 floating-point operations per second) machines.
    [Show full text]