Simulation Environment for Life Sciences @ BSC

Objective: Overview of simulation technologies used in Life Sciences and their specific adaptation to HPC environment.

Place: UPC Campus Nord premises, Vertex building, room VS218.


12 March 13 March
09.00 - 10.30  Welcome & Introduction (JLG) 09.00 - 09.30  Simulation DBs (MoDEL) (PA)
10.30 - 11.00  Break 09:30 - 10.30 Simulation Data Management (PA)
11.00 - 11.45  Atomistic MD Algorithm (JLG) 10.30 - 11.00 Break
11.45 - 12.30  Algorithm improvements & HPC (JLG) 11.00 - 12.00  Coarse Grained MD (AE)
12.30 - 14.00   Lunch break 12.00 - 13.00  Application Examples. Simulation of Nucleic Acids (PD)
14.00 - 15.00  Simulation Setup (AH) 13.00 - 14.00  Lunch break
15.00 - 16.00   Setup and Analysis Hands On (AH) 14.00 - 16.00  Coarse Grained MD manual slides Hands on (AE)
16.00 - 16.30  Break 16.00 - 16.30  Break&
16.30 - 18.00  Setup and Analysis Hands On (AH) 16.30 - 18.00  Free Hands on

AE: Agustí Emperador (IRB) , AH: Adam Hospital (IRB), JLG: Josep Ll. Gelpi (BSC-IRB- UB), PD: Pablo Dans (IRB), PA: Pau Andrio (BSC), RG: Ramón Goñi (BSC),

Software to be installed locally

Alternatively you may install a VirtualBox VM (2Gb)

Web tools

Selected references & URLs

General Review

Orozco M, Orellana L, Hospital A, Naganathan AN, Emperador A, Carrillo O, Gelpi JL. Coarse-grained representation of protein flexibility. Foundations, successes, and shortcomings. Adv Protein Chem Struct Biol 2011, 85:183-215.

Orozco M, Luque FJ. Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. Chem Rev 2000, 100:4187-4226.

Larsson P, Hess B, Lindahl E. Algorithm improvements for molecular dynamics simulations. Wiley Interdisciplinary Reviews-Computational Molecular Science 2010, 1:93-108.

Buch I, Harvey MJ, Giorgino T, Anderson DP, De Fabritiis G. High-throughput all-atom molecular dynamics simulations using distributed computing. J Chem Inf Model 2010, 50:397-403.

Hospital, A, Gelpi, J.L. High-throughput molecular dynamics simulations. Toward a dynamic PDB. WIRE 2013 (Early view) DOI: 10.1002/wcms.1142


Mackerell AD, Wiorkiewiczkuczera J, Karplus M. An all-atom empirical energy function for the simulation of nucleic-acids. Journal of the American Chemical Society 1995, 117:11946-11975.

MacKerell AD, Bashford D, Bellott M, Dunbrack RL, Evanseck JD, Field MJ, Fischer S, Gao J, Guo H, Ha S, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. Journal of Physical Chemistry B 1998, 102:3586-3616.

Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM, Ferguson DM, Spellmeyer DC, Fox T, Caldwell JW, Kollman PA. A 2nd generation force-field for the simulation of proteins, nucleic-acids, and organic-molecules. Journal of the American Chemical Society 1995, 117:5179-5197

Kaminski GA, Friesner RA, Tirado-Rives J, Jorgensen WL. Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides. Journal of Physical Chemistry B 2001, 105:6474-6487

MD Codes and helper applications

Harvey M, Giupponi G, De Fabritiis G. ACEMD: Accelerated molecular dynamics simulations in the microseconds timescale. J. Chem. Theory and Comput 2009, 5.            

Case DA, Darden TA, Cheatham I, T.E., Simmerling CL, Wang J, Duke RE, Luo R, Walker RC, Zhang W, Merz KM, et al. AMBER 12. University of California, San Francisco. 2012

Brooks BR, Brooks CL, 3rd, Mackerell AD, Jr., Nilsson L, Petrella RJ, Roux B, Won Y, Archontis G, Bartels C, Boresch S, et al. CHARMM: the biomolecular simulation program. J Comput Chem 2009, 30:1545-1614.

Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. Journal of Chemical Theory and Computation 2008, 4:435-447.

Hospital A, Andrio P, Fenollosa C, Cicin-Sain D, Orozco M, Gelpi JL. MDWeb and MDMoby: an integrated web-based platform for molecular dynamics simulations. Bioinformatics 2012, 28:1278-1279.

Nelson MT, Humphrey W, Gursoy A, Dalke A, Kale LV, Skeel RD, Schulten K. NAMD: A parallel, object oriented molecular dynamics program. International Journal of Supercomputer Applications and High Performance Computing 1996, 10:251-268.

Trajectory Databases

Rueda M, Ferrer-Costa C, Meyer T, Perez A, Camps J, Hospital A, Gelpi JL, Orozco M. A consensus view of protein dynamics. Proc Natl Acad Sci U S A 2007, 104:796-801.

Simms AM, Toofanny RD, Kehl C, Benson NC, Daggett V. Dynameomics: design of a computational lab workflow and scientific data repository for protein simulations. Protein Eng Des Sel 2008, 21:369-377.

Meyer T, D'Abramo M, Hospital A, Rueda M, Ferrer-Costa C, Perez A, Carrillo O, Camps J, Fenollosa C, Repchevsky D, et al. MoDEL (Molecular Dynamics Extended Library): a database of atomistic molecular dynamics trajectories. Structure 2010, 18:1399-1409.

Analysis tools

Camps J, Carrillo O, Emperador A, Orellana L, Hospital A, Rueda M, Cicin-Sain D, D'Abramo M, Gelpi JL, Orozco M. FlexServ: an integrated tool for the analysis of protein flexibility. Bioinformatics 2009, 25:1709-1710

PCAsuite. Compression based on Essential dynamics
Meyer T, Ferrer-Costa C, Perez A, Rueda M, Bidon-Chanal A, Luque FJ, Laughton CA, Orozco M. Essential dynamics: A tool for efficient trajectory compression and management. Journal of Chemical Theory and Computation 2006, 2:251-258