Computational Biology

Theoretical methods and mathematical modeling can be applied to investigate complex physical systems that are found in biological organisms. Computer simulations are a common tool to analyze the interactions of biological structures at the molecular level, or the connections among different subsystems in large-scale networks, which are crucial to understand the basic mechanisms in the proteomics and metabolomics fields


Ligand binding. Receptor-ligand interactions are responsible of molecular recognition, binding, transport and release. These mechanisms are fundamental in many processes that are vital in any living organism, and can be exploited in several applications of nanotechnology. Molecular dynamics simulations and docking calculations are theoretical tools that allow us to determine a model of interaction between a receptor and a variety of ligands. These computational techniques are useful to predict structural and dynamic properties that determine the functionality of a molecular complex between transport proteins and small compounds, such as ligands of pharmaceutical interest.


Protein folding. Proteins are molecular ‘nanomachines’ that acquire their tridimensional structure through a spontaneous process of folding. This phenomenon is central for some of the basic mechanisms of life. Furthermore, their competitive processes (unfolding and misfolding) lead to a loss of function and are correlated to some serious pathologies. Protein folding is an intrinsically complex molecular process that needs to be investigated with a combined theoretical and experimental approach. Molecular dynamics simulations and other computational methods are able to reveal crucial details of the folding reaction.


Regulatory RNA. Overcoming the decades’ old view according to which they are mere intermediates between DNA and proteins in the central dogma of molecular biology, RNAs are increasingly recognized for the remarkable variety of functional roles they play in eukaryotic gene expression. Non-coding RNAs (ncRNAs), in particular, are deeply involved in developmental programs and disease, to the point of having become targets of novel therapeutic approaches that aim at either inhibiting or enhancing their functionality. Importantly, most ncRNAs act by protein-mediated binding to other nucleic acids, leading to complex inter-dependencies be- tween coding and ncRNAs, RNA-binding proteins and DNA. Their central role in regulation is perhaps best exemplified by microRNAs (miRNAs), small ncRNAs of 20–25 nucleotides (nt) that mediate post-transcriptional regulation (PTR) via gene expression silencing in animals. miRNAs are believed to affect the expression of about two-thirds of protein-coding genes in humans. Their apparent ubiquity, together with the strong topological heterogeneities that characterize the PTR network that maps out the known interactions between miRNAs and other RNA species (such as mRNAs), has lead to the idea that competition for shared miRNAs can cause an effective positive interaction between different transcripts, currently referred to as the ‘ceRNA effect’. In this context, we are interested in quantifying the role and effectiveness of miRNA-mediated RNA cross-talk, with the goal of clarifying (i) under which conditions it outperforms other regulatory mechanisms, for instance in processing gene expression noise; (b) whether it can carry a significant systemic role; and (c) its relevance in specific biological cases of differentiation and disease.


Cell growth and its biosynthetic costs. The coupling between the physiology of cell growth and cellular composition has been actively investigated since the 1940s. In exponentially growing bacteria, such interdependence is best expressed in a quantitative way by the bacterial ‘growth laws’ that directly relate the protein, DNA and RNA content of a cell to the growth rate. Many such laws have been experimentally characterized and many more are currently being probed at increasingly high resolution. The emerging scenario suggests that proteome organization in bacteria is actively regulated in response to the growth conditions. Several phenomenological models explain the origin of different growth laws at coarse-grained levels. By contrast, genome-scale approaches probing such relationships at the molecular level are far less developed. We have developed a mathematical modeling scheme called Constrained Allocation Flux Balance Analysis or CAFBA, in which the costs of gene expression are accounted for effectively through a single global constraint on metabolic fluxes, encodeing for the relative adjustment of proteome sectors at different growth rates. Using bacteria as the initial model organisms, we are interested in quantifying the trade-off between cell growth and its associated biosynthetic costs, generating testable predictions about the way in which the usage of metabolic pathways and protein expression levels are modulated by the growth conditions.


Cell-to-cell variability in exponentially growing bacteria. Current experimental techniques (see e.g. the ‘mother machine’) can probe physiological variability by characterizing e.g. growth rate distributions for bacterial populations at single cell resolution. These distributions reflect noise at various levels, from intracellular stochasticity in gene expression and metabolite levels to fluctuations in the extracellular medium. However, upon controlling the latter, they provide a window to analyze the role of noise in the genotype-phenotype relationship. Straightforward sam- pling of the feasible space predicted by mathematical models, however, does not explain the observed statistics. We are therefore interested in identifying a physical or biological principle that drives the selection of observed growth states and hence shed light on the origin of the observed phenotypic diversity.



De Martino

CNR Researcher



Full Professor



CNR Researcher

Facilities and Labs

S.Li.M. Lab @ Roma


  1. M Mori et al. Constrained Allocation Flux Balance Analysis, PLOS Comp Biol 12:e1004913 (2016). DOI:10.1371/journal.pcbi.1004913
  2. D De Martino et al, Growth against entropy in bacterial metabolism: the phenotypic trade-off behind empirical growth rate distributions in E. coli, Phys Biol 13:036005 (2016). DOI:1088/1478-3975/13/3/036005
  3. S Grigolon et al, Noise Processing by MicroRNA-Mediated Circuits: the Incoherent Feed-Forward Loop, Revisited, Heliyon 2:e00095 (2016). DOI:  1016/j.heliyon.2016.e00095
  4. Martirosyan A et al, Probing the Limits to MicroRNA-Mediated Control of Gene Expression, PLOS Comp Biol 12(1): e1004715 (2016). DOI: 10.1371/journal.pcbi.1004715
  5. Evoli, L. Mobley, R. Guzzi, B. Rizzuti, Multiple   binding modes of ibuprofen in human serum albumin identified by absolute binding free energy calculations,    bioRxiv, 8, 1-27, (2016) doi:10.1101/068502
  6. Neira, B. Rizzuti, J. L. Iovanna, Determinants of the pKa values of ionizable residues in an intrinsically disordered protein, Archives of Biochemistry and Biophysics, 595, 1-16, (2016) doi: 10.1016/
  7. Capuani F et al, Quantitative constraint based computational model of tumor-to-stroma coupling via lactate shuttle, Sci Rep 5:11880 (2015). DOI: 10.1038/srep11880
  8. Rizzuti, R. Bartucci, L. Sportelli, R. Guzzi, Fatty acid binding into the highest affinity site of human serum albumin observed in molecular dynamics simulation, Archives of Biochemistry and Biophysics, 579, 18-25, (2015) doi: 1016/
  9. Evoli, R. Guzzi, B. Rizzuti, Molecular simulations of ß-lactoglobulin complexed with fatty acids reveal the structural basis of ligand affinity to internal and possible external binding sites, Proteins: Structure, Function, and Bioinformatics, 82, 2609-2619, (2014)     doi: 1002/prot.24625
  10. Pantusa, R. Bartucci, B. Rizzuti, Stability of trans-resveratrol associated with transport proteins, Journal of Agricultural and Food Chemistry, 62, 4384-4391, (2014) doi: 1021/jf405584a
  11. D De Martino et al. Inferring metabolic phenotypes from the exometabolome through a thermodynamic variational principle. New J Phys 16: 115018 (2014). DOI:  1088/1367-2630/16/11/115018
  12. M Figliuzzi et al, RNA based regulation: dynamics and response to perturbations of competing RNAs. Biophys J 107:1011 (2014). DOI: 1016/j.bpj.2014.06.035
  13. A De Martino et al, Identifying all moiety conservation laws in genome-scale metabolic networks. PLOS ONE 9:e100750 (2014). DOI: 10.1371/journal.pone.0100750
  14. A Seganti et al. Searching for feasible stationary states in reaction networks by solving a Boolean constraint satisfaction problem. Phys Rev E 89:022139 (2014). DOI: 1103/PhysRevE.89.022139

Other Selected Publications:

  1. FA Massucci et al, Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective. BMC Sys Biol 7:103 (2013). DOI:   1186/1752-0509-7-103
  2. M Figliuzzi et al, MicroRNAs as a selective channel of communication between competing RNAs, Biophys J 104:1203 (2013). DOI: 1016/j.bpj.2013.01.012
  3. A Seganti et al. Boolean constraint satisfaction problems for reaction networks, J Stat Mech P09009 (2013). DOI: 10.1088/1742-5468/2013/09/P09009
  4. D De Martino et al. Counting and correcting thermodynamically infeasible flux cycles in genome-scale metabolic networks, Metabolites 3:946 (2013). DOI: 3390/metabo3040946
  5. FA Massucci et al. A novel methodology to estimate metabolic flux distributions in constraint-based models, Metabolites 3:838 (2013). DOI: 3390/metabo3030838
  6. Evoli, Guzzi, B. Rizzuti, Dynamics and unfolding pathway of chimeric azurin variants: insights from molecular dynamics simulation, Journal of Biological Inorganic Chemistry, 18, 739-749, (2013) doi: 10.1007/s00775-013-1017-1
  7. Rizzuti, V. Daggett, Using simulations to provide the framework for experimental protein folding studies, Archives of Biochemistry and Biophysics 531, 128-135, (2013)   doi: 1016/
  8. Guzzi, Rizzuti, R. Bartucci, Dynamics and binding affinity of spin-labelled stearic acids in ß-lactoglobulin: evidences from EPR spectroscopy and molecular dynamics simulation, Journal of Physical Chemistry B, 116, 11608-11615 (2012) doi: 10.1021/jp3074392
  9. A De Martino, D De Martino, R Mulet and G Uguzzoni. Reaction networks as systems for resource allocation: a variational principle for non-equilibrium steady states. PLoS ONE 7:e39849 (2012). DOI: 1371/journal.pone.0039849
  10. D De Martino, M Figliuzzi, A De Martino and E Marinari. A scalable algorithm to explore the Gibbs energy landscape of genome-scale metabolic networks. PLoS Comp Biol 8:e1002562 (2012). DOI: 1371/journal.pcbi.1002562

Latest News

ICONS - A Symposium on Colloidal Nanocrystals


October 10-11, 2019


Cnr Nanotec Lecce


The symposium focuses on colloidal semiconductor nanocrystals (also known as quantum dots), which are a central topic in materials science and nanotechnology nowadays. The event will bring together renowned scientists in this research field discussing on fundamentals and future directions of this promising class of materials.  The program consists of two half-day sessions with four talks each covering synthesis, characterization, properties, and applications of colloidal semiconductor nanocrystals. The detailed program will be soon available.

The symposium is intended for anyone that is interested on this topic and will be free of charge. Attendees may present their own work in the dedicated poster session. Registration, with an eventual abstract submission for poster presentation, is required (please, send an e-mail to: The event is co-organized by Consiglio Nazionale delle Ricerche - Istituto di Nanotecnologia and Istituto Italiano di Tecnologia.

The final program is now available for download: ICoNS_Program[/vc_column_text][/vc_column][/vc_row]


 01 luglio 2019 - ore 14:15


Cnr Nanotec Lecce


Realizzato nell'ambito delle attività del progetto "TecnoMed Puglia - Tecnopolo per la medicina di precisione", il meeting è dedicato allo studio delle malattie neurodegenerative: dai nuovi biomarcatori alle piu recenti modellizzazioni, per una migliore comprensione dei meccanismi di base e quindi per lo sviluppo di terapie sempre più ritagliate sul singolo paziente.

EIT RawMaterials Roadshow

21 giugno 2019 ore 09:00 – 15:00


Lecce, Aula Fermi Edificio Aldo Romano, Campus Ekotecne, Via Lecce-Monteroni


Farà tappa a Lecce il prossimo 21 giugno, presso l’Aula Fermi dell’edificio IBIL all’interno del Campus Ecotekne, l’EIT RawMaterials, la piattaforma per il sostegno all’innovazione finanziata dall’Istituto Europeo di Innovazione e Tecnologia (EIT).


L’EIT ha creato le cosiddette KIC – Knowledge Innovation Community, comunità che mirano alla promozione dell’innovazione e della formazione in Europa in settori cruciali, sostenendo l’imprenditorialità e favorendo il passaggio di nuove idee dalla fase di incubazione al mercato.


La EIT RawMaterials si impegna ad affrontare la sfida globale dell’approvvigionamento delle materie prime in Europa attraverso programmi e progetti che mirano allo sviluppo di tecnologia nell’intera catena di valore delle materie prime: dall’esplorazione delle risorse, all’industria mineraria, dai processi metallurgici alla sostituzione delle materie prime critiche o tossiche, dal riciclo dei materiali dei prodotti a fine vita sino alla progettazione di prodotti per l’economia circolare. Nell’ambito dei programmi di sviluppo a livello regionale, la EIT RawMaterials ha creato un Hub nella Regione Puglia coordinato da ENEA, al fine di incrementare il coinvolgimento degli ecosistemi locali nelle attività della KIC e del suo partenariato.


Il MedinHub avrà inoltre l’obiettivo di raggiungere nuove organizzazioni e promuovere la partecipazione delle industrie e delle PMI più innovative, nonché il coinvolgimento delle prestigiose università e centri di ricerca dell’area.


Link per la registrazione: