Simplified models, including implicit-solvent and coarse-grained models, are useful tools to investigate the physical properties of biological macromolecules of large size, like protein complexes, large DNA/RNA strands and chromatin fibres. While advanced Monte Carlo techniques are quite efficient in sampling the conformational space of such models, the availability of realistic potentials is still a limitation to their general applicability. The recent development of a computational scheme capable of designing potentials to reproduce any kind of experimental data that can be expressed as thermal averages of conformational properties of the system has partially alleviated the problem. Here we present a program that implements the optimization of the potential with respect to the experimental data through an iterative Monte Carlo algorithm and a rescaling of the probability of the sampled conformations. The Monte Carlo sampling includes several types of moves, suitable for different kinds of system, and various sampling schemes, such as fixed-temperature, replica-exchange and adaptive simulated tempering. The conformational properties whose thermal averages are used as inputs currently include contact functions, distances and functions of distances, but can be easily extended to any function of the coordinates of the system. Program summary: Program title: MonteGrappa Catalogue identifier: AEUO_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEUO_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 139,987 No. of bytes in distributed program, including test data, etc.: 1,889,541 Distribution format: tar.gz Programming language: C Computer: Any computer with C compilers Operating system: Linux, Unix, OSX RAM: Bytes depend on the size of the system, typically 4 GB Classification: 3, 16.1 External routines: gsl, MPI (optional) Nature of problem: Optimize an interaction potential for coarse-grained models of biopolymers based on experimental data expressed as averages of conformational properties Solution method: Iterative Monte Carlo sampling coupled with minimization of the chi2 between experimental and back-calculated data making use of a reweighting algorithm Running time: Hours to days, depending on the complexity of the problem.

MonteGrappa : an iterative Monte Carlo program to optimize biomolecular potentials in simplified models / G. Tiana, F. Villa, Y. Zhan, R. Capelli, C. Paissoni, P. Sormanni, E. Heard, L. Giorgetti, R. Meloni. - In: COMPUTER PHYSICS COMMUNICATIONS. - ISSN 0010-4655. - 186(2015), pp. 93-104. [10.1016/j.cpc.2014.09.012]

MonteGrappa : an iterative Monte Carlo program to optimize biomolecular potentials in simplified models

G. Tiana;R. Capelli;C. Paissoni;R. Meloni
2015

Abstract

Simplified models, including implicit-solvent and coarse-grained models, are useful tools to investigate the physical properties of biological macromolecules of large size, like protein complexes, large DNA/RNA strands and chromatin fibres. While advanced Monte Carlo techniques are quite efficient in sampling the conformational space of such models, the availability of realistic potentials is still a limitation to their general applicability. The recent development of a computational scheme capable of designing potentials to reproduce any kind of experimental data that can be expressed as thermal averages of conformational properties of the system has partially alleviated the problem. Here we present a program that implements the optimization of the potential with respect to the experimental data through an iterative Monte Carlo algorithm and a rescaling of the probability of the sampled conformations. The Monte Carlo sampling includes several types of moves, suitable for different kinds of system, and various sampling schemes, such as fixed-temperature, replica-exchange and adaptive simulated tempering. The conformational properties whose thermal averages are used as inputs currently include contact functions, distances and functions of distances, but can be easily extended to any function of the coordinates of the system. Program summary: Program title: MonteGrappa Catalogue identifier: AEUO_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEUO_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 139,987 No. of bytes in distributed program, including test data, etc.: 1,889,541 Distribution format: tar.gz Programming language: C Computer: Any computer with C compilers Operating system: Linux, Unix, OSX RAM: Bytes depend on the size of the system, typically 4 GB Classification: 3, 16.1 External routines: gsl, MPI (optional) Nature of problem: Optimize an interaction potential for coarse-grained models of biopolymers based on experimental data expressed as averages of conformational properties Solution method: Iterative Monte Carlo sampling coupled with minimization of the chi2 between experimental and back-calculated data making use of a reweighting algorithm Running time: Hours to days, depending on the complexity of the problem.
Coarse-grained models; Force fields; Monte Carlo methods
Settore FIS/03 - Fisica della Materia
2015
Article (author)
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/243597
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 11
  • ???jsp.display-item.citation.isi??? 10
social impact