430 likes | 446 Views
BioSimGRID is a worldwide repository that integrates biomolecular simulation data, enabling comparative analysis and promoting collaboration in computational systems biology. It provides a centralized platform for deposition, analysis, and dissemination of biomolecular simulation data, enhancing the efficiency and accessibility of research in structural biology and drug discovery.
E N D
BioSimGRID and BioSimGRID ’lite’ • Towards a worldwide repository for biomolecular simulation • www.biosimgrid.org Philip C Biggin http://indigo1.biop.ox.ac.uk phil@biop.ox.ac.uk
Overview • Introduction • - Motivation • - Consortium • - Case studies – added value from comparisons • Design • - Architecture • - Data schema • How to use • - Deposition • - Analysis • - Worldwide application • The Future • Towards computational systems biology
Current Paradigm for MD Simulations • Target selection: literature based; interesting protein/problem • System preparation: highly interactive; slow; idiosyncratic • Simulation: diversity of protocols • Analysis: highly interactive; slow; idiosyncratic • Dissemination: traditional – papers, posters, talks • Archival: ‘archive’ data … and then mislay the tape! • No third party involvement
Integrating Simulations and Structural Biology of Proteins Sequence alignment Biomedically relevant homologue(s) Novel structure (RCSB) bacterial K channel bioinformatics & structural biology Homology model(s) mammalian K channel Biomolecular simulation database MD simulations dynamics in membrane BioSimGRID Comparative analysis Interaction site dynamics Evaluation/refinement of model Biological and pharmacological simulation & modelling e.g. drug discovery drug docking calculations drug discovery
York Nottingham Oxford RAL London Bristol Southampton Consortium • Oxford:Mark Sansom, Paul Jeffreys, Bing Wu, Kaihsu Tai • Southampton:Jon Essex, Simon Cox, Stuart Murdock, Muan Hong Ng, Hans Fogohr, Steven Johnston • London:David Moss • Nottingham:Charlie Laughton • York:Leo Caves • Bristol: Adrian Mulholland
Comparative Simulations: Drug Receptors • Why? – increase significance of results • Sampling – long simulations and multiple simulations • Sampling via biology – exploiting evolution • Biology emerges from comparisons… • e.g. mammalian receptor vs. bacterial binding protein glutamate • Rat GluR2 EC fragment • Major receptor in mammalian brains – drug target • MD simulations with/without bound ligands • Analyse inter-domain motions D1 D2
GluR2 – Flexibility & Gating… Kainate empty Glutamate > >> “ON” “OFF” 4 empty • Flexibility depends on ligand occupancy & species • Gating mechanism – decrease in flexibility on channel activation • But … incomplete sampling • Need: longer simulations & comparative simulations 3 +Glu RMSD (Å) 2 1 +Kai 0 2.0 0 0.5 1.0 1.5 time (ns)
+ Gln GlnBP – A Bacterial Binding Protein MD Simulation X-ray structures empty Gln bound empty Gln bound • GlnBP – bacterial 2-domain periplasmic binding protein • Similar fold to mammalian GluR2 • X-ray shows ligand binding induces domain closure • MD shows ligand binding reduces inter-domain motions - cf. GluR2 simulations
Case Study 2.. OMPLA AChE Acetylcholinesterase Outer-membrane phospholipase
So how do compare… • Similar active sites or similar motions • Different structures • Simulated with different MD packages (analysis difficult if not visualization) • On different hard drives/tapes/CDs/DVDs. • Under different graduate students’ desks • Under different postdocs’ beds • In different rubbish bins!
Answer… Create a wordwide repository of molecular simulations…. BioSimGrid = BioSimDB + Toolkits + Integration
Apache / Tomcat / SSL / Python Authentication Authorisation Accounting Trajectory Query Tool SQL Editor Data Retrieval Tool Data Deposition Tool Video/Img Engine HTML Generator Analysis Tool Database BioSimGrid Architecture… GUI Web Application Python Application HTTP(S) SSH Service TCP/IP Middle- ware BioSim Data Engine / Storage Resource Broker TCP/IP DB/Data Flat Files
Cross-software Analysis… • BioSimDB = PDB (or NDB) for MD • enable discovery of new science (cf. genomics/proteomic initiatives) CHARMM NAMD BioSimDB LAMMPS AMBER GROMACS TINKER
It’s a Distributed Database • Nobody has enough disk space in one place anyway • Distributed and duplicate • Any piece of information is stored in at least two sites • …for resilience
Database Database Current Architecture soton.biosimgrid.org oxford.biosimgrid.org BioSim Data Engine Services BioSim Data Engine Services IDA MCAT SRB Server SRB Server MCAT SRB Agent IDA SRB Agent F/F Interface DB Interface F/F Interface DB Interface F/F Engine DB Engine F/F Engine DB Engine Flat Files Cache Flat Files Cache
Data Schema • The hierachy is like that in the PDB: • Chain residue atom coordinate • …but also extended in the time dimension: frames
Metadata.. • …is the data about data • MD setup, parameters, instantaneous properties, etc. • People currently write this in papers • People forget something • The disciplined way:- • …structured schema
Deposition… Unified deposition for trajectories from any packages.
BioSimDB Toolkit • Analysis tools • Radius of Gyration • Surface and Volume • RMSD/RMSF • Centre of Mass • Inter-atomic distances • Distance matrix • Internal angles • Principal Component Analysis • Average structure
New workflow with BioSimGrid • Target selection: literature based; interesting protein/problem • Perform simulation (or use someone else’s) • Protocals more systematically recorded/checked/confirmed • Archive data to BioSimGrid • Analyse shared data (either locally or distributed) • Dissemination: traditional – papers, posters, talks • Store results in BioSimGrid • Third parties can analyse data you deposit
That’s dandy - but who is this aimed at? • Novice and Expert.. • Novice (web/GUI) • Makes selections • Guided through the options • Can only do specific things • Difficult to make mistakes • Expert (employ scripting) • Python interpreter • Much available • Reasonably unrestricted
Example sessions Even in script mode the syntax is quite informative:- FC = FrameCollection(`2, 100-200`) myRMSD = RMSD(FC) myRMSD.createPNG() Provide biochemists with little computational experience a means of analysing computational data and obtain meaningful results.
Example sessions Viewlet of a session; Demo4.html
BioSimGrid ‘Lite’ • Light version before final rollout • Provides equilibrated lipid bilayer boxes • Also provides ontogeny: How the box came about… • …metadata • …equilibration process (all the frames)
Deliverables to Date… • Database schema • Sample database (with test trajectories) • Prototype shared between 2 sites • Analysis tools – preliminary versions (about 14 tools) • Interface to database for data retrieval • Python hosting environment
Roadmap • Dec 2002 – project started • July 2003 – (internal) prototype • September 2003 – working prototype (All Hands meeting) • November 2003 – test ‘real world’ applications • December 2003 – multi-site prototype • 2004 – multi-site deposition of data • 2005 – open up to additional groups for deposition/testing
If you are interested… The team would like to hear from interested parties especially with new ideas etc • Benefits to you • New directions are implemented • Toolkit suits your needs • Shared development of code • Faster and more thorough development • BioSimGrid Benefits • Larger user community • More work gets done • Code is efficient. • BioSimGrid and community is successful
Future Directions in the GRID context • HTMD – simulations coupled to structural genomics • Diamond light source • Computational system biology – virtual outer membrane • HPCx • Multiscale biomolecular simulations – from QM/MM to meso-scale modelling • GRID-enabled simulations BioSimGrid
Structural Genomics & HTMD synchrotron compute GRID MD database novel biology… • Overall vision – simulation as an integral component of structural genomics • Needs capacity computation – GRID? • MD database (distributed) – BioSimGRID
d- d- d- d- TolC OMPLA Pi d+ OpcA OmpT OmpX FhuA LamB OmpA OmpF PhoE TonB MalE FhuD Pi PiBP Towards a Virtual Outer Membrane (vOM) • First step towards computational systems biology – a suitable system • Bacterial OMs – 5 or 6 proteins = 90% of protein content • Structures or good homology models of proteins are available • Complex lipid – outer leaflet is lipopolysaccharide (LPS) • Minimum system size ca. 2.5x106 atoms; simulation times ca. 50 ns • cf. current FhuA – 80,000 atoms & 10 ns – need HPCx
Multiscale Biomolecular Simulations QM (Bristol) Drug-binding (Southampton) Protein Motions (Oxford) Drug Diffusion (London) • Membrane bound enzymes – major drug targets (cf. ibruprofen, anti-depressants, endocannabinoids) • Complex multi-scale problem: QM/MM; ligand binding; membrane/protein fluctuations; diffusive motion of substrates/drugs in multiple phases • Need for GRID-based integrated simulations
References… • K. Tai, S. Murdock, B.Wu, MH Ng, S. Johnston, H. Fangohr, S. Cox, P Jeffreys, J. Essex, M.S.P. Sansom. Org. Biomol. Chem :: Under review • MH Ng, S. Johnston, S. Murdock, B. Wu, K. Tai, H. fangohr, S. Cox, J. Essex, M.S.P. Sansom, P.Jeffrey. • UK E-Science Programme All Hands Meeting 2004 :: Accepted. • 3. Python Website – www.python.org • 4. BioSimGrid – www.biosimgrid.org
Acknowledgements Oxford e-Science Center Professor Paul Jeffreys Dr Bing Wu (database management) Matthew Dovey Ivaylo Kostadinov Oxford Professor Mark Sansom Dr Carmen Domene Dr Alessandro Grottesi Dr Andrew Hung Dr Daniele Bemporad Dr Shozeb Haider Dr Kaihsu Tai (curation and integration) Dr George Patargias Oliver Beckstein Jennifer Johnston Syma Khalid Jorge Pikunic Pete Bond Zara Sands Jonathan Cuthbertson Sundeep Deol Jeff Campbell Yalini Pathy Loredana Vaccaro Shiva Amiri Katherine Cox Robert d’Rozario John Holyoake Samantha Kaye Anthony Ivetac Sylvanna Ho Southampton Dr Stuart Murdock (generic analysis tools) Dr Muan Hong Ng (data retrieval) Dr Hans Fangohr Steven Johnston Prof Simon Cox Dr Jon Essex Elsewhere Leo Caves (York) Charles Laughton (Nottingham) David Moss (Birkbeck) Oliver Smart (Birmingham) Adrian Mulholland (Bristol) Marc Baaden (Paris) BBSRC DTI The Wellcome Trust GSK EC (TMR) OeSC (EPSRC & DTI)EPSRC OSC(JIF) MRC
More information… team@biosimgrid.org www.biosimgrid.org