Researchers in the Biophysics area at BIFI apply experimental and computational tools and methodologies in an interdisciplinary environment for understanding the behavior of biological systems, from molecules (proteins, nucleic acids, small molecules, etc.) to whole organisms and ecosystems, from a quantitative perspective. Proteins and organisms with biotechnological and/or biomedical relevance are of special interest. Research within this area has many applications in different fields:

Biotechnology (protein stabilization and modeling, protein function design, molecular machines, nanostructures, molecular modeling, electronic transport, catalysis, electromagnetic interaction with matter, and protein glycosylation and its role in signaling)

Biomedicine (drug discovery and design, pharmacological target identification and validation, protein-DNA interaction)

Biology (epidemiology, evolution, complex networks)

The Biophysics area at BIFI encompasses 10 lines of research:

  • Protein Folding and Molecular Design
  • Biomolecular Interactions
  • Protein glycosylation and its role in disease
  • Complex systems and networks
  • Physical modeling of biomolecules
  • Molecular dynamics and electronic structure
  • Flavoenzymes: action mechanisms and biotechnology
  • Protein Misfolding and Amyloid Aggregation
  • Clinical Diagnosis and Drug Delivery
  • Signal transduction and membrane protein therapies
Protein Folding and Molecular Design

Protein Folding and Molecular Design

Head of the Research Line:

Javier Sancho


Dr. Juan José Galano-Frutos
Dr. Alejandro Mahía
María Conde-Giménez
Sandra Salillas
Helena García
Ritwik Maity
Patricia Bruñen



figure1-angarica-et-al-2016The Protein Folding and Molecular Design group is devoted to the study and improvement of protein molecules. Our labs at BIFI and at the Faculty of Sciences are fully equipped for integrative research using techniques of Cellular and Molecular Biology, Biochemistry, Biophysics, and Computation. During the last 25 years we have studied the principles of protein stability, folding and binding, and the relationship between protein dynamics and function. Protein stability is a very important topic from both theoretical and practical perspectives, and we still lack enough quantitative understanding of it. We have used a variety of models, such as flavodoxin, to develop biotechnological and computational strategies to assist in protein stabilization using rational principles. Nowadays we develop statistical techniques to calculate accurately protein conformational stability from Molecular Dynamics simulations of x-ray structures and atomistic models of unfolded ensembles. We are also deeply involved in developing accurate tools for Genetic Interpretation,

We are quite active in the drug discovery field. In recent years, we have investigated protein conformational diseases in order to understand their molecular causes and to discover small molecules that can be developed into new drug therapies. In this endeavour, we combine classical Biochemistry and Biophysics with novel HTP screening techniques, Bioinformatics, Synthetic Chemistry and in vivo testing. Our projects include testing in animal models novel antimicrobials effective against the human pathogen Helicobacter pylori, and the discovery and rational improvement of pharmacological chaperones to rescue defective human enzymes (e.g. mutated PAH responsible for Phenylketonuria) and to inhibit the aggregation of a number of human peptides and proteins involved in amyloidogenesis.












Rather than in techniques, we are interested in problems so we usually combine experimental approaches with computational studies, as needed. We have developed atomistic models of the denatured ensemble (ProtSA), a genome scale Prion predictor (PrionScan), a predictor of locally unstable segments of proteins (ProteinLIPS), and Molecular Dynamics tools for the identification of SNPs involved in conformational diseases.

Our group is pleased to accept talented and motivated doctoral students or postdocs interested in our work. Contact us.


Relevant publications

1.-Molecular Dynamics Simulations for Genetic Interpretation in Protein Coding Regions: Where we Are, Where to Go and When. Juan José Galano-Frutos, Helena García-Cebollada, Javier Sancho. Briefings in Bioinformatics. 2020, In Press (doi: 10.1093/bib/bbz146).

2.-Design, Synthesis and Efficacy testing of nitroethylene- and 7-nitrobenzoxadiazol-based flavodoxin inhibitors against Helicobacter pylori drug-resistant clinical strains and in Helicobacter pylori-infected mice. Sandra Salillas, Miriam Alías, Valérie Michel, Alejandro Mahía, Ainhoa Lucía, Liliana Rodrigues, Jessica Bueno, Juan José Galano-Frutos, Hilde De Reuse, Adrián Velázquez-Campoy, José Alberto Carrodeguas, Carlos Sostres, Javier Castillo, José Antonio Aínsa, María Dolores Díaz-de-Villegas, Ángel Lanas, Eliette Touati, Javier Sancho. Journal of Medicinal Chemistry 2019, 62:6102-6115.

3.-Accurate Calculation of Barnase and SNase Folding Energetics using short MD simulations and an Atomistic Model of the Unfolded Ensemble. Evaluation of Force Fields and Water Models. Juan José Galano-Frutos, Javier Sancho. J. Chem. Inf. Model. 2019, 59:4350-4360.

4.- A novel small molecule inhibits α-synuclein aggregation, disrupts amyloid fibrils and prevents degeneration of dopaminergic neurons. J. Pujols, S. Peña-Díaz, D. F. Lázaro, F. Peccati, F. Pinheiro, D. González, A. Čarija, S. Navarro, M. Conde-Gimenez, J. García, S. Guardiola, E. Giralt, X. Salvatella, J. Sancho, M. Sodupe, T. F. Outeiro, E. Dalfó, and S. Ventura. Proc Natl Acad Sci USA (2018) 115:10481–10486.

5.- Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate. E. Lamazares, S. Vega, P. Ferreira, M. Medina, J. J. Galano, M. Martínez-Júlvez, A. Velázquez-Campoy & J. Sancho. Phys Chem Chem Phys. 19:19021-19031 (2017).

6.- Exploring the complete mutational space of the LDL receptor LA5 domain using molecular dynamics: Linking SNPs with disease phenotypes in familial hypercholesterolemia. V. E. Angarica, M. Orozco and J. Sancho. Human Molecular Genetics, 25:1233-1246 (2016).

7.- Rational stabilization of complex proteins: a divide and combine approach. E. Lamazares, I. Clemente, M. Bueno, A.Velázquez-Campoy, J. Sancho . Scientific Reports, 5, 9129 (2015).

8.- Improved flavodoxin inhibitors with potential therapeutic effects against Helicobacter pylori infection. J.J. Galano, M. Alías, R. Pérez, A. Velázquez-Campoy, P. S. Hoffman, J. Sancho. Journal of Medicinal Chemistry. 56-15, pp.6248-6258 (2013).

9.- Discovery of novel inhibitors of amyloid b-peptide 1-42 aggregation. L. C. López, S. Dos-Reis, A. Espargaró, J. A. Carrodeguas, M-L Maddelein, S Ventura, J Sancho. J. Med. Chem.55:9521-9530 (2012).

10.- ProtSA: A web application for calculating sequence specific protein solvent accessibilities in the unfolded ensemble. J. Estrada, P. Bernadó, M. Blackledge, J. Sancho. BMC Bioinformatics. 10:article104 (2009).


Main research projects

1.- PIREPRED: Cross-border Network for the interpretation of neonatal screening: from mutation to patient. EFA086/15 (ERDF-POCTEFA Interreg VA). 01/09/16-30/11/20. IP: Javier Sancho.

2.- Quantitative understanding of protein stability by modelling and simulation, and application to the kinetic and thermodynamic interpretation of variants in a single amino acid (QUPS)   PID2019-107293GB-I00: 01/06/20-30/05/23.MICIN. IP: Javier Sancho.

3.- FLAV4AMR: Flavodoxin inhibitors to kill resistant bacteria: EU (JPI-AMR) 2019/2022 IP: Javier Sancho.

4.Diseño y desarrollo de un nuevo test de fenilalanina para la gestión familiar de la fenilcetonuria. PID2019-107293GB-I00 (ISCIII). 2019-2021 IP: Javier Sancho.

5.- Diseño y desarrollo de nuevos fármacos para medicina personalizada y su aplicación a test diagnósticos. LMP30_18 (Gobierno de Aragón). 2019/2020 IP: Javier Sancho.



We maintain active collaborations with many national and international groups and networks: INPEC. See publications. See also project webs: RedMut, Pirepred.

Biomolecular Interactions

Biomolecular Interactions

Head of the Research Line:

Adrian Velazquez-Campoy


José Luis Neira (Universidad Miguel Hernández, Elche)
David Ortega Alarcón (doctorando FPI)
Ana Jiménez Alesanco (investigadora contratada)



All biological processes can be described as a sequence of interaction events between biomolecules and conformational changes coupled to those interactions. We work on three main aspects, from both basic and applied perspectives, related to biological molecules:

Biophysical characterization of proteins with biotechnological and/or biomedical relevance: conformational landscape, interactions with other biomolecules, conformational changes, phase diagrams, cooperative phenomena…

Development of experimental methodologies for the study of biomolecular interactions, cooperative phenomena and functional regulation in biological systems.

Development and implementation of high-throughput screening procedures for the identification of bioactive compounds able to modulate the function of pharmacological protein targets.


The biophysical characterization of a protein provides

– Information on the protein conformational landscape. The conformational landscape is the set of conformational states, structurally and energetically distinguishable, that are accessible and significantly populated. The distribution of the protein molecules among these conformational states depends on Gibbs free energy of each state, which in turn is a function of the environmental variables (temperature, pH, ionic strength,…), as well as the interaction with other biomolecules or mutations in the sequence. Because each conformational state may exhibit particular functional properties, the conformation and the function of a given protein are regulated through changes in the environmental variables and biochemical signaling (presence of interacting biomolecules) or mutations. This connection between conformation and function constitutes the basis for the allosteric behavior in proteins: modulation and control of protein conformation (and, therefore, protein function) through interactions with other biomolecules (e.g. ions, cofactors, inhibitors, etc.).

We apply biophysical tools and methods for studying the conformational landscape of proteins of biotechnological and and/or biomedical interest.

bmi_02 Information on the protein function, regulation and evolution. The physiological function of proteins, and its regulation, always relies on the interaction with other molecules. Therefore, the detailed structural and energetic characterization of its interactions is fundamental for understanding its function and regulation. The interaction with structural metal ions and cofactors is of special interest for protein function regulation. In addition, this information might shed light into evolutive aspects: mutation-induced and/or ligand-induced stabilization of alternative conformational and functional states may give rise to new functions.

We apply biophysical tools and methods to study binding interactions and its coupling with protein stability and the influence on the conformational landscape.

bmi_13 Information on the molecular and energetic basis of a disease. Very often defective proteins are causative agents for diseases. The mutation-induced and/or ligand-induced stabilization of non-native partially (or completely) unfolded states may lead to improper folding, inactivation and/or early degradation of a given protein, or to the formation of toxic aggregates. Likewise, mutations may distort or impair required interactions for proper function. The biophysical study of its conformational landscape provides information on conformational states that may be relevant regarding its function and regulation, as well as information on the impact of mutations on the interactions required for proper protein function.

We apply biophysical tools and methods to assess the impact of mutations and ligand binding on the stability and function of pharmacological protein targets.-

– Information for establishing strategies for the identification of bioactive compounds. In both infectious diseases (in which the inactivation of a key protein of the life cycle of the pathogen with an inhibitor is the usual therapeutic strategy) as well as in genetic or metabolic disorders (in which both the inactivation of an increased aberrant protein activity with an inhibitor or the rescue of a defective protein activity with a pharmacological chaperone represents alternative therapeutic strategies) it is necessary to modulate and control the function of a given protein. In order to achieve this goal we must identify low-molecular weight compounds: 1) able to interact with the protein; 2) able to modulate appropriately the protein function; 3) avoiding interactions with unwanted targets minimizing side-effects; and 4) exhibiting an optimal pharmacokinetic (ADMET) profile. The biophysical characterization of a protein target provides relevant information about points 1-3, representing a crucial element for designing, validating and optimizing molecular screening procedures.

We develop and implement high-throughput screening procedures for identifying bioactive compounds against pharmacological protein targets.


figurea-new        figureb-new

– Information for optimizing bioactive compounds (affinity, selectivity, induced-conformational changes in the target…). The biophysical study of protein interactions provides information on the relevant underlying intermolecular interactions (i.e. hydrogen bonds, van der Waals, electrostatic and hydrophobic interactions). The thermodynamic dissection of a given interaction allows estimating the enthalpic and entropic contributions to the Gibbs free energy of interaction. The partition of the Gibbs energy into its enthalpic and entropic terms (thermodynamic profile or signature of the binding) is fundamental for describing that interaction, because those two energetic terms reflect intermolecular interactions of different nature (regarding strength, specificity, susceptibility to mutations, and ability to redesign and engineering). It is known that ligands of similar affinity interacting with a protein will bind differently and will display different binding properties if they exhibit different thermodynamic signature; in particular, their interaction will be driven by different intermolecular forces, they will show different susceptibility to mutations in the protein target, and they will induce different conformational changes. Therefore, not only the binding affinity, but also the complete thermodynamic signature is relevant for understanding the mode of interaction and as a key decision criterion for the optimization of binding affinity and selectivity.
We develop and implement procedures for the thermodynamic characterization of binding interactions. A biophysical thermodynamic characterization of lead compound binding provides invaluable information for compound optimization.


Cooperative phenomena are inherent to the behavior of biomolecules and constitute the molecular and energetic basis for protein regulation and allostery. We may distinguish different cooperativity types or levels in biomolecules, which are intimately interrelated:

– Intrinsic protein structural stability. The conformational landscape and the intrinsic structural and energetic features of the protein conformation are key determinants of its behavior and function. Many proteins fold spontaneously and adopt a structured well-defined native state. Other proteins may populate non-native partially folded states with biological relevance. Moreover, many proteins do not fold into a well-defined structured conformational state unless they interact with another molecule (intrinsically disordered proteins), which may have implications regarding their function and regulation. Interestingly, conformational landscape and intrinsic structural stability are not conserved among structurally homologous proteins.
We study the structural stability of proteins, with particular interest on intrinsically (partially) unfolded proteins.


bmi_14 Coupling between binding and protein conformational changes. Very often protein interactions are accompanied by conformational changes. This is a reflection of the modulation of the conformational landscape of the protein through the interaction with a given ligand. The extent of the conformational change will depend on many factors, in particular, the intrinsic structural stability of the protein and the specific structural and functional features of the ligand. For example, it has been demonstrated the structural determinants in a ligand for binding affinity do not necessarily coincide with those determinants for allosteric regulation (which is mediated through conformational changes coupled to binding). Also, ligand binding induces protein stabilization giving rise to the distinction between strong or weak coupling between conformational changes and binding, as well as between induced-fit and conformational selection. In fact, induced-fit can be considered as a limit case of conformational selection when a very large Gibbs free energy gap exists between the two conformational states interconnected by ligand binding.
We study conformational changes coupled to ligand binding in proteins, as well as the interplay between folding and binding in protein function regulation and protein binding specificity.

bmi_10– Binding cooperativity. Homotropic cooperativity (binding cooperativity between molecules of the same ligand binding to several sites in a protein) and heterotropic cooperativity (binding cooperativity between molecules of different ligands binding to one or several sites in a protein) constitute the basis for protein function regulation and allostery. In a broad sense, allostery can be defined as the control and modulation of protein conformation (conformational landscape) through ligand binding. Therefore, any protein can be considered an allosterically regulated protein. The most common example is the pH dependency of the structural stability or the binding affinity of a ligand in a given protein; in this particular case, the binding or dissociation of protons taking place at ionizable functional groups in the protein (or in the ligand) alters the energetics of unfolding or ligand binding. Binding affinity and binding cooperativity are not conserved properties among structurally similar ligands.
We develop and implement methodologies for characterizing and assessing cooperative (homotropic and heterotropic) binding in proteins and its relation to the biological function.




Relevant publications

1.- J.L. Neira, J. Bintz, M. Arruebo, B. Rizzuti, T. Bonacci, S. Vega, A. Lanas, A. Velazquez-Campoy, J.L. Iovanna, O. Abian. Identification of a drug targeting an intrinsically disordered protein involved in pancreatic adenocarcinoma. Scientific Reports 2017 7:39732

2.- R. Claveria-Gimeno, P.M. Lanuza, I. Morales-Chueca, O.C. Jorge, S. Vega, O. Abian, M. Esteller, A. Velazquez-Campoy. The intervening domain from MeCP2 enhances the DNA affinity of the methyl binding domain and provides an independent DNA interaction site. Scientific Reports 2017 7:41635

3.- Biophysical screening for identifying pharmacological chaperones and inhibitors against conformational and infectious diseases. Velazquez-Campoy, J. Sancho, O. Abian, S. Vega. Current Drug Targets 2016 17:1492-1505.

4.- R. Sant’Anna, P. Gallego, L. Robinson, A. Pereira, Ne. Ferreira, F. Pinheiro, S. Esperante, I. Pallares, O. Huerta, R. Almeida, N. Reixach, R. Insa, A. Velazquez-Campoy, D. Reverter, N. Reig, S. Ventura. Repositioning Tolcapone as a potent inhibitor of transthyretin amyloidogenesis and associated cellular toxicity. Nature Communications 2016 7:10787.

5.- S. Vega, O. Abian, A. Velazquez-Campoy. A unified framework based on the binding polynomial for characterizing biological systems by isothermal titration calorimetry. Methods 2015, 76:99-115

6.- R.M. Buey, R. Ledesma-Amaro, A. Velazquez-Campoy, M. Balsera, M. Chagoyen, J.M. de Pereda, J.L. Revuelta. Guanine nucleotide binding to the Bateman domain mediates the allosteric inhibition of eukaryotic IMP dehydrogenases. R.M. Buey, R. Ledesma-Amaro, A. Velazquez-Campoy, M. Balsera, M. Chagoyen, J.M. de Pereda, J.L. Revuelta. Nature Communications 2015 6:8923.

7.- K. Gonzalez-Arzola, I. Diaz-Moreno, A. Cano-Gonzalez, A. Diaz-Quintana, A. Velazquez-Campoy, B. Moreno-Beltrán, A. Lopez-Rivas, M.A. De la Rosa. Structural basis for inhibition of the histone chaperone activity of SET/TAF-Iβ by cytochrome c.  Proceedings of the National Academy of Sciences USA 2015 112:9908-9913.

8.- J.J. Galano, M. Alias, R. Perez, A. Velazquez-Campoy, P.S. Hoffman, J. Sancho. Improved flavodoxin inhibitors with potential therapeutic effects against Helicobacter pylori infection. Journal of Medicinal Chemistry 2013 56:6248-6258.

9.- O. Abian, S. Vega, J. Sancho, A. Velazquez-Campoy. Allosteric inhibitors of the NS3 protease from the hepatitis C virus. PLoS ONE 2013 8:e69773.

10.- S. Chopra, A. Palencia, C. Virus, A. Tripathy, B.R. Temple, A. Velazquez-Campoy, S. Cusack, J.S. Reader. Plant tumour biocontrol agent employs a tRNA-dependent mechanism to inhibit leucyl-tRNA synthetase. Nature Communications 2013 4:1417.

11.- P.T. Martins, A. Velazquez-Campoy, W.L. Vaz, R.M. Cardoso, J. Valerio, M.J. Moreno. Kinetics and thermodynamics of chlorpromazine interaction with lipid bilayers: Effect of charge and cholesterol. Journal of the American Chemical Society 2012 134:4184-4195.

12.- O. Abian, S. Vega, J.L. Neira, A. Velazquez-Campoy. Conformational stability of hepatitis C virus NS3 protease. Biophysical Journal 2010 99:3811-3820.

13. N. Cremades, A. Velazquez-Campoy, M. Martinez-Julvez, J.L. Neira, I. Perez-Dorado, J. Hermoso-Dominguez, P. Jimenez, A. Lanas, P.S. Hoffman, J. Sancho. Discovery of specific flavodoxin inhibitors as potential therapeutic agents against Helicobacter pylori infection. ACS Chemical Biology 2009 4:928-938.

14.- O. Abian, J.L. Neira, A. Velazquez-Campoy. Thermodynamics of zinc binding to hepatitis C virus NS3 protease: A folding by binding event. Proteins: Structure, Function and Bioinformatics 2009 77:624-636.

15.- A.L. Pey, M. Ying, N. Cremades, A. Velazquez-Campoy, T. Scherer, B. Thöny, J. Sancho, A. Martinez. Identification of pharmacological chaperones as new therapeutic agents to treat phenylketonuria. Journal of Clinical Investigation 2008 118:2858-2867.



Main research projects

1.- Structure, energetics, and simulation of (partially) unfolded conformations in proteins. Towards quantitative atomic models for the protein stability. Ministerio de Economía y Competitividad (BFU2016-78232-P). 2016-2019. PI: Adrián Velázquez Campoy / Javier Sancho (co-IPs) (Universidad de Zaragoza – BIFI).

2.- New infrastructure for the Laboratorio Avanzado de Cribado e Interacciones Moleculares de Aragon (LACRIMA). Diputación General de Aragón y Ministerio de Economía y Competitividad (Proyecto de Infraestructuras Científico-Tecnológicas UNZA15-EE-3250). Universidad de Zaragoza. 2016-2018. PI: Javier Sancho (Universidad de Zaragoza – BIFI).

3.- Between atom and cell: Integrating molecular biophysics approaches for biology and healthcare (MOBIEU). European Cooperation in Science and Technology (eCOST, COST Action CA15126). ARBRE (Association of Resources for Biophysical Research in Europe). 2016-2020. PI: Patrick England (Institute Pasteur). Adrián Velázquez-Campoy (Management Committee).

4.- Validation of a new, quick, non-invasive diagnostic method in serum for early detection of pancreatic cancer (PANCal). Asociación Española de Gastroenterología – Club Español del Páncreas. Instituto Investigación Sanitaria de Aragón, Hospital Clínico Universitario Lozano Blesa, Hospital Universitario Miguel Servet, Hospital de Donosti, Hospital de Barbastro, Universidad de Zaragoza – Instituto BIFI. 2015. PI: Olga Abián (Instituto de Investigación Sanitaria de Aragón (IIS) – Universidad de Zaragoza – BIFI).

5.- Protein stability: Basic principles of (partially) unfolded states and molecular studies on conformational disorders. Ministerio de Economía y Competitividad (BFU2013-47064-P). Universidad de Zaragoza. 2014-2017. PI: Adrián Velázquez Campoy / Javier Sancho (co-IPs) (Universidad de Zaragoza – BIFI).

6.- NEUROMED – Diagnosis and treatment of three neurodegenerative diseases (Parkinson, Phenylketonuria and TTR amyloidosis). ERDF – SUDOE Interreg IVB (SOE4/P1/E831 – Neuromed). Universidad de Zaragoza, Universidad Autónoma de Barcelona, Instituto de Biologia Molecular e Celular (Porto), Universidade de Coimbra, CNRS, INSERM (Languedoc-Roussillon). 2014-2015. PI: Javier Sancho (Universidad de Zaragoza – BIFI).

7.- Bioavailability of amphiphilic ligands – Drugs and metabolites. Ministerio de Ciencia e Innovación (Proyectos de Movilidad – Acciones Integradas, PRI-AIBPT-2011-1025). Universidad de Zaragoza, Universidad de Coimbra. 2012-2013. PI: Adrián Velázquez Campoy (Universidad de Zaragoza – BIFI).

8.- NS3 protease from hepatitis C virus: Identification of competitive and allosteric inhibitors. Ministerio de Ciencia e Innovación (BFU2010-19451). Universidad de Zaragoza. 2010-2013. PI: Adrián Velázquez Campoy (Universidad de Zaragoza – BIFI).

9.- NS3 protease from hepatitis C virus: Competitive and allosteric inhibitors, and molecular basis of drug resistance. Universidad de Zaragoza (UZ2009-BIO-05). Universidad de Zaragoza. 2010. PI: Adrián Velázquez Campoy (Universidad de Zaragoza – BIFI).

10.- Conformational equilibrium of the NS3 protease from hepatitis C virus: Non-native states and identification of a new type of allosteric inhibitors. Diputación General de Aragón (PI044/09). Universidad de Zaragoza. 2009-2011. PI: Adrián Velázquez Campoy (Universidad de Zaragoza – BIFI).




  • Olga Abian (Instituto Aragones de Ciencias de la Salud and BIFI, Spain)
  • Juan Ausio (University of Victoria, Canada)
  • Rui Brito (Universidade de Coimbra, Portugal)
  • Pierpaolo Bruscolini (Universidad de Zaragoza and BIFI, Spain)
  • Jose A. Carrodeguas (Universidad de Zaragoza and BIFI, Spain)
  • Irene Diaz-Moreno (Instituto de Bioquimica Vegetal y Fotosintesis – CSIC, Spain)
  • Maria Fillat (Universidad de Zaragoza and BIFI, Spain)
  • Marcos R. Fontes (São Paulo State University, Brazil)
  • Ernesto Freire (The Johns Hopkins University, USA)
  • Enrique Garcia-Hernandez (Universidad Nacional Autonoma de Mexico, Mexico)
  • Ruben Martinez-Buey (Universidad de Salamanca, Spain)
  • Milagros Medina (Universidad de Zaragoza and BIFI, Spain)
  • Maria João Moreno (Universidade de Coimbra, Portugal)
  • Arturo Muga (Universidad del Pais Vasco, Spain)
  • Jose A. Navarro (Instituto de Bioquimica Vegetal y Fotosintesis – CSIC, Spain)
  • Julian Pardo (Universidad de Zaragoza and IIS-Aragon, Spain)
  • Santiago Ramon-Maiques (CNIO, Spain)
  • David Reverter (Universidad Autonoma de Barcelona, Spain)
  • Javier Sancho (Universidad de Zaragoza and BIFI, Spain)
  • Jayaraman Sivaraman (National University of Singapore, Singapore)
  • Maria A. Urbaneja (Universidad del Pais Vasco, Spain)
  • Salvador Ventura (Universidad Autonoma de Barcelona, Spain)
Protein glycosylation and its role in disease

Protein glycosylation and its role in disease

Head of the Research Line:

Ramón Hurtado Guerrero


Ana García García
Javier Macías León
Andrés Manuel González Ramírez
Víctor Taleb
Erandi Lira Navarrete



Our group is interested in the study of glycosyltransferases, glycosylhydrolases and carbohydrate-binding proteins/modules involved in human diseases. We use protein X-ray crystallography as the main tool complemented with enzymology, inhibition studies, etc, in order to study the molecular mechanisms of enzymes that are involved in the synthesis, modification and degradation of glycoconjugates, oligo and polysaccharides (see more relevant publications below).

We majorly work in Protein O-glycosylation and glycosyltransferases responsible of this post-translational modification. In particular we are currently working in several glycosyltransferases such as Protein O-fucosyltransferases 1 and 2 (POFUT1 and 2), FUT8 and the large family of GalNAc-Ts. While POFUT1 and 2 fucosylate folded domains such as EFG and TSR domains, respectively, GalNAc-Ts mainly glycosylate unstructured regions present in a large number of proteins such as mucins . FUT8 is responsible of the so-called core fucosylation (Figure 1).

Figure 1. Scheme summarizing the structure of FUT8 complexed to GDP and a biantennary N-glycan.


Furthermore we are interested in the elucidation of the reaction coordinates and the molecular mechanism by using transition state analogues or Michaelis complexes. Finally these studies will be important for the design of new future drugs with potential therapeutic applications.

Although not related to the Glycobiology field, we have a close collaboration with Prof. Guadix and Dr Conejo-García from University of Granada in the development of new inhibitors against human choline kinase α1.



Relevant publications

1.- Structural insights into mechanism and specificity of O-GlcNAc transferase. Clarke AJ*, Hurtado-Guerrero R*, Pathak S*, Schüttelkopf AW*, et al. EMBO J. 2008, 27(20): 2780-8. *equal contributation as first author.

2.- Molecular mechanisms of O-GlcNAcylation. Hurtado-Guerrero R, Dorfmueller HC, van Aalten DM. Curr Opin Struct Biol. 2008, 18(5): 551-7.

3.- Recent structural and mechanistic insights into post-translational enzymatic glycosylation. Ramon Hurtado-Guerrero* and Gideon J Davies. Current Opinion in Chemical Biology. 16(5-6):479-87, 2012. *corresponding author

4.- The Vibrio cholerae colonization factor GbpA possesses a modular structure that governs binding to different host surfaces. Wong E, Vaaje-Kolstad G, Ghosh A, Ramon Hurtado-Guerrero, et al. PLoS Pathogen, 2012.

5.- The mechanism of allosteric coupling in choline kinase α1 revealed by a rationally designed inhibitor. María Sahún-Roncero, Belén Rubio-Ruiz, Giorgio Saladino, Ana Conejo-García, Antonio Espinosa, Adrián Velázquez-Campoy, Francesco Luigi Gervasio, Antonio Entrena and Ramon Hurtado-Guerrero*. Angewandte Chemie (selected as VIP), 52(17):4582-6, 2013. *corresponding author

6.- Combined structural snapshots and metadynamics reveal a substrate-guided front-face reaction for polypeptide GalNAc-transferase T2. Erandi Lira-Navarrete, Javier Iglesias-Fernández, Wesley F. Zandberg, Ismael Compañón, Yun Kong, Francisco Corzana, B. Mario Pinto, Henrik Clausen, Jesús M. Peregrina, David Vocadlo, Carme Rovira* and Ramon Hurtado-Guerrero*. Angewandte Chemie International Edition, 53(31):8206-10, 2014. *corresponding author

7.- Dynamic interplay between catalytic and lectin domains of GalNAc transferases modulates protein O-glycosylation. Erandi Lira-Navarrete, Matilde de las Rivas, Ismael Compañón, María Carmen Pallarés, Yun Kong, Javier Iglesias-Fernández, Gonçalo J. L. Bernardes, Jesús M. Peregrina, Carme Rovira, Pau Bernadó, Pierpaolo Bruscolini, Henrik Clausen, Anabel Lostao, Francisco Corzana, and Ramon Hurtado-Guerrero*. Nature Communications, 5;6:6937, 2015. *corresponding author

8.- X-ray Structures decipher the Non-equivalence of Serine and Threonine O-glycosylation points: Implications for the Molecular Recognition of the Tn Antigen by an anti-MUC1 Antibody. Nuria Martínez-Sáez,‡ Jorge Castro-López,‡ Jessika Valero-González,‡ David Madariaga, Ismael Compañón, Víctor J. Somovilla, Míriam Salvadó, Juan L. Asensio, Jesús Jiménez-Barbero, Alberto Avenoza, Jesús H. Busto, Gonçalo J. L. Bernardes, Jesús M. Peregrina,* Ramón Hurtado-Guerrero,* and Francisco Corzana*. Angewandte Chemie International, 54(34):9830-9834, 2015. *corresponding author

9.- Mucin Architecture behind the Immune Response: Design, Evaluation and Conformational Analysis of an Antitumor Vaccine Derived from an Unnatural MUC1 Fragment. Martínez-Sáez N, Supekar NT, Wolfert MA, Bermejo IA, Hurtado-Guerrero R, Asensio JL, Jiménez-Barbero J, Busto JH, Avenoza A, Boon G-J, Peregrina JM, and Corzana F. Chemical Science, 2016. DOI: 10.1039/C5SC04039F

10.- A proactive role of water molecules in acceptor recognition by Protein-O-fucosyltransferase 2. Jessika Valero-González, Christina Leonhard-Melief , Erandi Lira-Navarrete, Gonzalo Jiménez-Osés, Cristina Hernández-Ruiz, María Carmen Pallarés, Inmaculada Yruela, Deepika Vasudevan, Anabel Lostao, Francisco Corzana, Hideyuki Takeuchi, Robert S. Haltiwanger, and Ramon Hurtado-Guerrero*. Nature Chemical Biology, accepted manuscript, 2016. DOI:10.1038/nchembio.2019. * corresponding author

11.- A trapped covalent intermediate of a glycoside hydrolase on the pathway to transglycosylation. Insights from experiments and QM/MM simulations. Lluís Raich, Vladimir Borodkin, Wenxia Fang, Jorge Castro-López, Daan van Aalten*, Ramon Hurtado-Guerrero* and Carme Rovira*. Journal of the American Chemical Society (JACS), 2016. DOI: 10.1021/jacs.5b10092. * corresponding authors

12.- The interdomain flexible linker of the polypeptide GalNAc transferases (GalNAc-Ts) dictates their long-range glycosylation preferences. Matilde de las Rivas, Erandi Lira-Navarrete, et al., and Ramon Hurtado-Guerrero*.Nature Communications, 2017. DOI: 10.1038/s41467-017-02006-0. *corresponding author.

13.- The use of fluoroproline in MUC1 antigen enables efficient detection of antibodies in patients with prostate cancer. Somovilla Víctor, et al., Hurtado-Guerrero Ramon, et al., Corzana Francisco. JACS, 2017. DOI:10.1021/jacs.7b09447.

14.– Water Sculpts the Distinctive Shapes and Dynamics of the Tn Antigens: Implications for their Molecular Recognition. Iris A. Bermejo, et al., Ramón Hurtado-Guerrero*. JACS, 2018 DOI: 10.1021/jacs.8b04801. *corresponding authors.

15.- Dynamic acetylation of cytosolic phosphoenolpyruvate carboxykinase toggles enzyme activity between gluconeogenic and anaplerotic reactions. Pedro Latorre-Muro, Josue Baeza, Eric A. Armstrong, Ramón Hurtado-Guerrero, et al., José A. Carrodeguas* and John M. Denu*. Molecular Cell, 2018, just accepted. Authors.

16.- Structural basis for arginine glycosylation of host substrates by bacterial effector proteins. Jun Bae Park, Young Hun Kim, Youngki Yoo, Juyeon Kim, Sung-Hoon Jun, Jin Won Cho, Samir El Qaidi, Samuel Walpole, Serena Monaco, Ana A. García-García, Miaomiao Wu, Michael P. Hays, Ramon Hurtado-Guerrero*, Jesus Angulo*, Philip R. Hardwidge, Jeon-Soo Shin* and Hyun-Soo Cho*. Nature Communications, 2018, just accepted. *joint corresponding authorship.

17.- Structural and mechanistic insights into the catalytic domain-mediated short-range glycosylation preferences of GalNAc-T4. de las Rivas Matilde, Paul Daniel Earnest James, Coelho Helena, Lira-Navarrete Erandi, Raich Lluís, Companon Ismael, Diniz Ana, Lagartera Laura, Jiménez-Barbero Jesús, Clausen Henrik, Rovira Carme, Marcelo Filipa, Corzana Francisco, Gerken Thomas* and Hurtado-Guerrero, Ramón*. ACS Central Science, 2018, just accepted. *joint corresponding authorship.

18. Polypeptide GalNAc-Ts: from redundancy to specificit.  de las Rivas Matilde, Lira-Navarrete, Erandi, Gerken Thomas* and Hurtado-Guerrero, Ramon*. . Current Opinion in Structural Biology, 2019. doi: 10.1016/ *joint corresponding authorship.

19. Structure-based design of potent tumor-associated antigens: modulation of peptide presentation by single atom O/S or O/Se substitutions at the glycosidic linkage. Compañón Ismael, Guerreiro Ana, Mangini Vincenzo Castro-López Jorge, Escudero-Casao Margarita, Avenoza Alberto, Busto Jesús, Castillón Sergio, Jiménez-Barbero Jesús, Asensio Juan, Jiménez-Osés Gonzalo, Boutureira Omar, Peregrina Jesús, Hurtado-Guerrero Ramón, Fiammengo Roberto, Bernardes Gonçalo and Corzana Francisco. . JACS, 2019. doi: 10.1021/jacs.8b13503

20. Mechanisms of redundancy and specificity of the Aspergillus fumigatus Crh transglycosylases. Fang Wenxia, Sanz Ana Belén, Galan Bartual Sergio, Wang Bin, Ferenbach Andrew, Farkaš Vladimír, Hurtado-Guerrero Ramon, Arroyo Javier and van Aalten DMF.  Nature Communications, 2019. doi: 10.1038/s41467-019-09674-0.

21. Molecular basis for Fibroblast growth factor 23 (FGF23) O-glycosylation and regulation by polypeptide GalNAc-T3. Matilde de las Rivas, Earnest Daniel, Yoshiki Narimatsu, Ismael Compañón, Kentaro Kato, Pablo Hermosilla, Aurélien Thureau, Laura Ceballos-Laita, Helena Coelho, Pau Bernado, Filipa Marcelo, Lars Hansen, Ryota Maeda, Anabel Lostao, Francisco Corzana, Henrik Clausen, Thomas Gerken and Ramón Hurtado-Guerrero*.  Nature Chemical Biology, 2020. doi: 10.1038/s41589-019-0444-x * corresponding autor.

22.  Structural basis for substrate specificity and catalysis of α1,6-fucosyltransferase. Ana García-García, Laura Ceballos-Laita, Sonia Serna, Raik Artschwager, Niels C. Reichardt, Francisco Corzana and Ramon Hurtado-Guerrero*. Just accepted in Nature Communications, 2020. doi: 10.1038/s41467-020-14794-z. * corresponding autor.

23. Recent advances in the design of Choline kinase α inhibitors and the molecular basis of their inhibition. Belén Rubio-Ruiz, Lucía Serán-Aguilera, Ramón Hurtado-Guerrero*, and Ana Conejo-García*.  Medicinal Research Reviews. doi: 10.1002/med.21746. *Joint corresponding authorship

Main research projects

1.- Study of glycosyltransferases involved in the Notch signalling pathway
MICINN, 2011-2014. 120,000 euros. PI: Ramón Hurtado-Guerrero

2.- Study of protein-carbohydrate interactions involved in human diseases. MEC, 2014-2016. 76,000 euros. PI: Ramón Hurtado-Guerrero

3.- BFU2016-75633-P (MICINN), 2017-2019. 220,000 euros plus extra funding for one PhD student. PI: Ramón Hurtado-Guerrero

4.- FEDER proyects for infrastructure and in particular for LACRIMA (2016-2017). UNZA15-EE-3250. PI: Javier Sancho. coPIs: Ramón Hurtado-Guerrero and others. 386,872 euros

5.- FEDER proyects for infrastructure and in particular for the adquistion of a FRET machine (2016-2017). UNZA15-EE-2922. PI: Nunilo Cremades. coPIs: Ramón Hurtado-Guerrero and others. 323,049 euros

6.- Crowdfunding for development of antifungal drugs. PI: Ramón Hurtado-Guerrero. 6600 euros, 2018-2019.

7.- Collaborative research network in Glycobiology “GLYCOBIOCHEM”. CTQ2017-90719-REDT. PI: Niels Reichardt, coPI: Ramón Hurtado-Guerrero and others. 17,000 euros, 01072018-30062020.

8.- DGA Group. Chemical Biology and Computation. E34_17R. PI: Pedro Merino, coPI: Ramón Hurtado-Guerrero. 45,971 euros, 2018-2019.

9.- Bases moleculares de la O-glicosilación tipo mucina y su aplicación en el tratamiento de tumores. LMP58_18 (DGA). 85,500 euros, 2019-2020. PI: Ramón Hurtado-Guerrero

10.- FEDER proyect for infrastructure (EQC2019-005847-P). Implementation of complex cellular models for drug discovery: Advanced equipment for the production, purification and analysis of proteins. PI: Milagros Medina. coPIs: Ramón Hurtado-Guerrero and others. 463,086.87 euros (2017-2020).

11.- Desarrollo de nanobodies y “adhirons” frente al dominio de unión de la glicoproteína “spike” del virus SARS-CoV-2 como tratamiento para la enfermedad CoVid19. DGA, 275,000 euros, 2020-2021. PI: Ramón Hurtado-Guerrero, coPI: Julián Pardo

12.- Desarrollo y validación de un test de inmunocromatografía para su aplicación en el diagnóstico serológico frente a SARSCoV2 en animales de compañía. Proyecto Santander-Unizar. 8500 euros, 2020-2021. PI: Pérez Cabrejas, M. Dolores, coPI among others: Ramón Hurtado-Guerrero.

13.- PID2019-105451GB-I00 (MICINN), 2020-2022. 266,200 euros plus extra funding for one PhD student. PI: Ramón Hurtado-Guerrero

14. Collaborative research network in Glycobiology “GLYCOBIOCHEM”. CTQ2017-90719-REDT. PI: Niels Reichardt, coPI: Ramón Hurtado-Guerrero and others. 17,000 euros, 01072018-30062020.

15.- COST “Innogly”. PI: Luigi Lay, coPIs: Ramón Hurtado-Guerrero and others. 2019-2022



Robert Haltiwanger, The University of Georgia
Henrik Clausen, University of Copenhagen
Daan van Aalten, University of Dundee
Philip Hardwidge, Kansas State University
Tom Gerken, Case Western Reserver University
Pedro Merino, Universidad de Zaragoza
Francisco Corzana, Universidad de La Rioja
Carme Rovira, Universidad de Barcelona
Filipa Marcelo, New University of Lisbon
Julián Pardo, Universidad de Zaragoza
Anabel Lostao, Universidad de Zaragoza

Complex Systems
and Networks (COSNET Lab)

Complex Systems and Networks (COSNET Lab)

Head of the Research Line:

Professor Yamir Moreno


Yamir Moreno (Group Leader)
Sandro Meloni (Researcher)
Carlos Gracia-Lázaro (Researcher)
Emanuele Cozzo (Postdoctoral researcher)
Sergio Arregui (PhD Student)
Alberto Aleta (PhD Student)
Felipe M. Cardoso (PhD Student)
Pablo Piedrahita (PhD Student)

BIFI members linked to the research line who work at the Department of Condensed Matter Physics, Faculty of Sciences, University of Zaragoza:

Professor Luis Mario Floría Peralta
Jesús Gómez-Gardeñes



The Complex Systems and Networks Lab (COSNET) was founded in 2003 by Professor Yamir Moreno and has an extensive scientific background in Statistical Physics and in the Physics of Complex Systems. During the last decade, the group has consolidated as a leading team at international level in the study of several topics related to the structure and dynamics of Complex Networks, Epidemiology, Systems Biology, Multiplex Networks, synchronization phenomena, technological and social networks, specially, in the emergence and evolution of online social movements, Evolutionary Game Theory, as well as in the analysis of human collective behavior and in the development of large-scale experiments to study cooperation in humans.

Complex Network Theory

Networked systems are all around us. The accumulated evidence that complex systems cannot be fully understood by studying only their isolated constituents, has given rise to the birth of a new movement of interest and research in the study of complex networks. The expectancy is that understanding and modeling the structure of a complex network would lead to a better cottoning on its dynamical and functional behavior. Though modern network theory has produced a number of relevant results in the last few years, it is still in its infancy, particularly, when it comes to applications in real systems and to the comprehension of the relation between structure and function (dynamics). The main purpose of this research line is the study of complex networks and the collective behavior of dynamical agents that interact among them following the couplings given by the topology of these complex networks. As a result, the fundamental objectives are:

  • The simultaneous characterization of the interactions and dynamics at a local scale and the study of their integration into a global and coherent dynamics at a system-wide scale.
  • The study of how global dynamics affects local interactions.
  • The statistical characterization of real networks.
  • The design of realistic models.
  • The study of other dynamical processes on complex networks and the emergence of collective behavior.
  • The development and application of analytical tools to study complex networks.
  • Fostering of a community of multidisciplinary scientists, who master the discipline of complex systems and use it for their research.
  • Identify the best course of action to transfer the acquired knowledge from basic sciences to the application level for a proper characterization and exploitation of real systems.

Complex network theory is particularly advantageous to explore several aspects of complexity. The fundamental idea is to discover the structure of interactions between the components of the system and the emergent behavior of many-body systems coupled to the underlying structure. This would improve our understanding and modeling capabilities so that we may control or predict the dynamics and function of complex networked systems. In addition, this approach does not rely on a detailed knowledge of the system’s constituents, but in the analysis of the relationship between them and therefore, allows obtaining universal results that can be generalized with relative ease (the study of epidemic spreading processes is equivalent to the spread of computer viruses). For example, biological networks like protein interaction networks, share many structural (scale-freeness) and dynamical (functional modules) features with other seemingly different systems such as the Internet and interaction patterns in social systems. Thus, systems as diverse as peer-to-peer networks, neural systems, socio-technical phenomena or complex biological networks can be studied with a unique theoretical and computational tool.

On the other hand, there are still many unknown systems and processes in which the new discoveries and techniques developed in the last years can shed light and provide novel results. For example, why scale-free networks are ubiquitous in Nature? Are there universal principles that govern the growth and evolution of these networks? How the dynamics of local interactions are spread and integrated into the global scale? And, in turn, how the macroscopic behavior of the system modulates interactions at lower scales? Finally, are there common patterns that can be identified not only structurally, but also in the functional organization of these systems? The latter and other questions should be studied by adopting new perspectives and approaches based on multidisciplinary research. The research results have important applications in problems such as: the modeling of biological processes at the molecular and cellular levels (metabolism, gene expression, etc.), the study of epidemic dynamics, the characterization of transport and diffusion processes in networks and communication technologies, several synchronization phenomena and the emergence of collective effects with applications in neuroscience and social systems. The ultimate goal is to understand the general principles governing these biological, social and technological systems in order to be able to predict, design and control the behavior of a wide range of real systems.

“The Science of Complexity”

The Science of Complexity has become extremely important in the study of many real systems. Complexity theory is based on the holistic principle antireductionist that considers the whole greater than the sum of its parts. This perspective allows an accurate analysis of different phenomena so as to predict the evolution of them over time. Complex behavior occurs when many interactions at the local scale collectively lead to unpredictable larger-scale outcomes. However, only very recently scientists have started to reconsider the traditional reductionism viewpoint that has frequently driven science. The accumulated evidence that systems as complex as a group of social animals, or the cells of a living system, cannot be fully understood by simply reducing them to a sum of their fundamental parts, has produced an increasingly large interest in the study of complex systems. These studies are revealing and explaining a range of emergent system behaviors and providing a deeper understanding of entire systems and their responses, with often surprising and unexpected results.

Multiplex Networks

Until very recently, most of the studies on the structure and dynamics of networks considered them as single layer graphs, in which all the interactions, no matter whether they were obtained in different conditions, appeared aggregated. However, as more data about real systems became available, it was evident that many systems are indeed made up by different interaction layers that are interdependent. Think of, for instance, a transportation system, where different transportation modes coexist in space. An accurate representation of such a system would represent each mode of transportation as one layer, with connections among them representing the possibility of commuting from one to another: from a bus line to a metro line, to a tram, etc. The same applies to online social systems, where an individual could be in several platforms (Twitter, Facebook, Google Plus, etc.) concurrently, thus communicating and being exposed to info through many information channels. A final example in the biological domain is given by cellular processes, where there are many biochemical pathways operating at the same time.

This new structural paradigm poses important challenges. First, we need to figure out when the multilayer structure is important and then, develop new metrics, algorithms and representations of such systems. In addition, dynamical processes such as diffusion and spreading dynamics should be reassessed, as it is expected that new phenomenology arises. In our group, we have been studying these networks since a few years and have already produced a significant number of contributions that have had a deep impact on the general subject of network science. This includes the study of diffusion processes, the spreading of multistrain and interacting diseases, new metrics and formalisms to characterize the structure of multilayer networks and the dynamics of multichannel information dissemination. We continue to explore this line and plan to tackle several interesting problems in the biological domain.




Epidemic spreading is a central issue in a variety of fields. In this context, both epidemic modeling and data collection about contact networks are intensively contributing to development of an accurate computational and theoretical framework to simulate how diseases spread and evolve as well as to the search of efficient immunization and vaccination policies. However, the lack of a complete description of connectivity maps and the singularities of the transmission mechanisms make the analysis extremely difficult. Much effort should still be invested in the design of efficient and reliable epidemic spreading models that benefits from data of connectivity patterns given by social maps of contacts (directionality, age, gender, etc.) and also to understand the dynamics of multi-strain or interacting diseases. Our research is mainly focus on the study of situations in which multiple pathogens coexist within the same host population, including systems of competing pathogens (e.g., seasonal influenza) or the so-called syndemic systems (e.g. HIV and Tuberculosis), i.e., the convergence of two or more diseases that act synergistically in a population to magnify the burden of disease. For these latter scenarios, considering multiple networks of contacts or layers is paramount.


Epidemiological Framework & Multilevel Approach

The study of epidemic spreading processes from a quantitative point of view is a vast field under intense investigation for nearly a century. Mathematical and computational modeling of infectious diseases is a collective endeavor of scientists coming from many scientific fields, from applied mathematicians and epidemiologists, to computer scientists and physicists. Physicist’s approaches to problems in Epidemiology invoke statistical physics, the theory of phase transitions and critical phenomena, to grasp the macroscopic behavior of epidemic outbreaks. It is not adventurous to claim that one of the main successful frameworks is the Mean-Field (MF) approximation, where homogeneity and isotropy are hypothesized to reduce the complexity of the system under study. This approach is mainly used to study the spreading of a disease in a system in which individuals are identified with the nodes of a (complex) network of contacts. Another important alternative approach aimed at describing the large-scale spreading of infectious diseases mathematically is the so-called meta-population approach. This framework describes a set of spatially structured interacting subpopulations as a network whose links denote the mobility of individuals across subpopulations. Each subpopulation consists of a number of individuals that are divided into several classes according to their dynamical state with respect to the modeled disease, for instance: susceptible, infected, removed, etc. The internal compartmental dynamics models the contagion dynamics by considering that people in the same subpopulation are in contact and may change their state according to their interactions and the disease dynamics. Finally, subpopulations also interact and exchange individuals due to mobility from one subpopulation to another.

Being able to understand analytically how the various model components influence the dynamics and interact is not always possible. More often, this is achieved numerically. We therefore pay particular attention to the development of efficient numerical methods that complement, when possible, analytical insights. The constraint of modeling a real process in which myriads of sources of complexity are involved makes the combination of massive data bases with relatively simple models the only sensible approach to model fitting, projection and understanding. Finally, it is worth mentioning that epidemiological research nowadays faces problems related to the lack of appropriate, disease-specific theoretical and computational models to understand the transmission mechanisms behind global public health threats. This is amplified by our current limited knowledge about the interplay among the various scales involved in the transmission of infectious diseases at the global scale.

Most epidemic models are developed assuming that the spreading process takes place on a single level (be it a single population, a meta-population system or a network of contacts). Therefore, pressing problems rooted at the interdependency of multi-scales call for the development of a whole new set of theoretical and simulation approaches. Our research is related the emergence and evolution of global diseases from the new paradigm of multi-scale, interdependent complex systems. Taking into consideration not only a single interaction level but also the interdependency between many levels and scales is a radical change, a change that makes a substantial difference as it allows to address problems such as comorbidity and the spreading of persistent infections and multi-strain diseases, among others. The outcomes may be rewarding and the resulting framework would be invaluable to better comprehend global health threats.

As a result, our goal is to develop a contemporary epidemiological framework that integrates the many aspects involved in the spreading of global diseases: from single networks of contacts to metapopulation systems, to the interaction between different diseases and strains and to the influence of human behavioral changes and mobility patterns.

In this context, our main objectives are:

  • To integrate individual-level approaches into meta-population schemes with the aim of adding 
further realism to the structure of subpopulations.
  • To generate new mathematical and computational modeling schemes to accurately describe the 
impact of human behavioral changes on the course of an epidemic.
  • To develop models that will allow addressing the effects of the interaction between persistent 
infections (and specifically, between TB and AIDS) in structured populations.
  • To significantly advance in a brand new theoretical and modeling framework to better grasp the role 
of competition between cross-reacting strains of infection in multi-scale diseases.
  • To compare the structural and dynamical properties of epidemic multi-scale models looking for similarities, organizational principles and universal dynamical patterns with other multi-level, interdependent complex systems.


Structure and Dynamics of Online Social Systems

The collection and mining of large amounts of digital data is allowing gaining deep insights into the structure and dynamics of large-scale socio-technical systems. In order to advance more in our current understanding of such systems, we need to develop a proper and efficient way to handle these datasets as well as to implement new algorithms and modeling tools that will ultimately allow to extract new knowledge from them. In this respect, it has become apparent that online social systems are playing a role in the development of different collective phenomena that manifest in the off-line world. Therefore, the analysis of the emergence and evolution of influential social movements in social platforms like Twitter is extremely useful if we aim at making further predictions about human collective behavior.

Today, techno-social systems pervade our world. ICT systems are leading to profound transformations in the way our society self-organize and generate and use information as well as on how humans interact, travel, behave, etc. Such transformations are not always traduced in an improved human society, as they can also lead to new forms of instabilities or make social and economic systems more fragile. Moreover, the ongoing technological revolution is taking place at a never-saw pace, which makes it even more urgent the development of tools that allow understanding the emerging properties of complex techno-social systems and to anticipate the consequences of new regulations, actions, or systems’ failures.

Although outstanding results have recently been obtained in modeling collective behavior in techno- social systems, we have not yet progressed enough in basic theoretical aspects and in the application of the generated knowledge to the characterization and understanding of real social phenomena. Moreover, with the unprecedented amount of data at our disposal nowadays, new challenges arise.

Present-day tools simply fail to keep up with the shifting challenges that the interaction between humans and complex ICT systems poses. For instance, think of a techno-social system like online social networks in which individuals engage in a multitude of categorical (politics, science, sports, technology, etc.) layers, giving to each of them a specific weight. Could we predict or simply understand how likely it is that a given rumor, idea or belief reaches system-wide proportions? Is there a general mechanism behind such kind of social phenomenon or are there different mechanisms unique to each categorical level or social niche? Are influential individuals globally defined or instead are they system’s dependent? Can we understand what gives rise to social uproar as recently witnessed? Or even more, are we in a position to answer those questions that we would ideally like to ask?



Evolutionary Game Theory and The Socio-Technical Man

Evolutionary Game Theory is a branch of mathematics that analyzes the strategic interaction between rational or irrational agents that can be individuals, groups, corporations, etc. The game consists of a set of players and strategies. Each player earns a reward that depends not only on his/her own strategy, but also on the strategy carried out by his/her competitors. This theory is very useful to reveal and predict human behavior and provides a framework and analytical tools for understanding a wide range of phenomena that occur in real life and are linked to decision-making by individuals or groups of individuals who interact with each other (see contributions of John von Neumann, Oskar Morgenstern and John Nash). The Prisoner’s Dilemma is the most studied model in Evolutionary Game Theory and is a symmetric, bi-personal, finite, static and non-zero-sum game in which players must choose between two strategies: Cooperate or defect. From an individual perspective, defect is always the best strategy. Nevertheless, cooperation gives the highest payoff collectively. This is the dilemma.

If we want to describe the Socio-technical man, we need to understand first, some basic problems such as: how humans interact with the environment and other individuals, how cooperative behavior emerges and survives, and how social networks evolve and shape the way we communicate and interact with each other. To this end, we should develop new ways to analyze existing data, perform controlled experiments with groups –of different sizes- of humans facing social dilemmas and hypothetical scenarios and come out with new theoretical and computational methods and algorithms. By doing this, we will be in the position to better understand what are the factors determining human behavior in a plethora of situations, and hopefully, provide hints to, e.g., policy makers, with the aim of creating a better and more sustainable Society for the future. On its turn, given the universality of many concepts and methods of Complexity Science, we are also sure that the new methods will contribute to the development of other areas (for instance, the Physics of non-equilibrium systems, Economics and even Ecology).

Finally, Evolutionary Game Theory is also very useful tool to study different phenomena and biological processes such as the ecology of bacterial population and the evolution of virus and species, in general.


The aim of the Laboratory of Experimental Economics “NECTUNT LAB” is the study of cooperative phenomena in humans as well as the applications of Evolutionary Game Theory in different fields of science. Methodologically, we make use of this theory and computer simulations to build more realistic models based on our findings in controlled experiments, both onsite and online. Currently, we are involved in a European Research Project (IBSEN). The main objective of the project is to develop the largest simulator of human behavior up to date, with the ultimate goal of answering questions such as: What are the mechanisms and real motivations that promote the emergence and evolution of cooperation in humans? How individuals behave in different contexts? How financial bubbles are formed?

On December 20, 2011, the largest experiment developed up to now with thousands of subjects playing a Prisoner’s Dilemma Game took place in Zaragoza, Spain. Almost 1300 students coming from 42 schools and institutes of the Autonomous Community of Aragon were involved in this unprecedented real-time experiment. The design, software and visualization platform were developed by BIFI researchers.


Relevant publications

1.- Complex Networks: Structure and Dynamics. S. Boccaletti, V. Latora, Y. Moreno, M. Chávez and D.-U. Hwang. Physics Reports 424, 175-308 (2006).

2.- Heterogeneous networks do not promote cooperation when humans play a Prisoner’s Dilemma. C. Gracia-Lázaro, A. Ferrer, G. Ruíz, A. Tarancón, J. A. Cuesta, A. Sánchez, and Y. Moreno. Proceedings of the National Academy of Sciences USA 109, 12922-12926 (2012).

3.- Evolutionary dynamics of group interactions on structured populations – A review. M. Perc, J. Gómez-Gardeñes, A. Szolnoki, L. M. Floría and Y. Moreno. Journal of the Royal Society Interface 10, 20120997 (2013).

4.- Host mobility drives pathogen competition in spatially structured populations. C. Poletto, S. Meloni, V. Colizza, Y. Moreno and A. Vespignani. PLoS Computational Biology 9 (8): e1003169 (2013).

5.- The role of hidden influentials in the diffusion of online information cascades. R. A Baños, J. Borge-Holthoefer and Y. Moreno. EPJ Data Science 2:6 (2013).

6.- Multilayer Networks. M. Kivela, A. Arenas, M. Barthelemy, J. P. Gleeson, Y. Moreno, and M. A. Porter. Journal of Complex Networks 2, 203-271 (2014).

7.- Behavioral transition with age in social dilemmas: From reciprocal youth to persistent response in adulthood. M. Gutiérrez-Roig, C. Gracia-Lázaro, J. Perelló, Y. Moreno, and A. Sánchez. Nature Communications 5:4362, doi: 10.1038/ncomms5362 (2014).

8.- Dynamics of interacting diseases. J. Sanz, C. -Y. Xia, S. Meloni and Y. Moreno. Physical Review X 4, 041005 (2014).

9.- Reputation drives cooperative behavior and network formation in human groups. J. A. Cuesta, C. Gracia-Lázaro, A. Ferrer, Y. Moreno, and A. Sánchez. Scientific Reports 5:7843 (2015).

10.- Characterizing two-pathogen competition in spatially structured environments. C. Poletto, S. Meloni, A. Van Metre, V. Colizza, Y. Moreno and A. Vespignani. Scientific Reports 5:7895 (2015).

11.- Dynamic instability of cooperation due to diverse activity patterns in evolutionary social dilemmas. C. -Y. Xia, S. Meloni, M. Perc and Y. Moreno, “”, Europhysics Letters 109, 58002 (2015).

12.- Spatiotemporal characterization of information-driven collective phenomena through transfer entropy. J. Borge-Holthoefer, N. Perra, B. Gonçalves, S. Gonzalez-Bailón, A. Arenas, Y. Moreno, and A. Vespignani. Science Advances 2, e1501158 (2016).

13.- Modeling the effects of network structure, competition and memory time on social spreading phenomena. J. P. Gleeson, K. P. O’Sullivan, R. A. Baños, Y. Moreno. Physical Review X 6, 021019 (2016).

14.- Humans conform to a reduced set of behavioral phenotypes when facing social dilemmas. J. Poncela-Casasnovas, M. Gutiérrez-Roig, C. Gracia-Lázaro, J. Vicens, J. Gomez-Gardeñes, J. Perelló, Y. Moreno, J. Duch, and A. Sánchez. Science Advances 2, e1600451 (2016).


Main research projects

1.- BRIDGING THE GAP: FROM INDIVIDUAL BEHAVIOR TO THE SOCIO-TECHNICAL MAN (IBSEN), European Commission. H2020 FET Open, Project number 662725, 2015-2018.

2.- DISTRIBUTED GLOBAL FINANCIAL SYSTEMS FOR SOCIETY (DOLFINS), European Commission. H2020 FET Proactive GSS, Project number 640772, 2015-2017.

3.- FOUNDATIONAL RESEARCH ON MULTILEVEL COMPLEX NETWORKS AND SYSTEMS (MULTIPLEX), European Commission. FET Proactive IP Project number 317532, 2012-2016.

4.- MATHEMATICAL FRAMEWORK FOR MULTIPLEX NETWORKS (PLEXMATH), European Commission. FET Proactive STREP Project number 317614, 2012-2015.



  • Alessandro Vespignani (Sternberg Distinguished University Professor College of Computer and Information Science, College of Science, Bouvé College of Health Sciences, Northeastern University, Boston, US)
  • Vittoria Colizza (Inserm and UPMC Université Paris 06, Faculté de Médecine, Paris & ISI Foundation, Turin, Italy)
  • Chiara Poletto (Researcher (Chargé de Recherche 2ème classe) EPICX-Lab, iPLESP, INSERM & UPMC UMR-S 1136)
  • Carlos Martín (Group of Mycobacterial Genetics, Faculty of Medicine, University of Zaragoza, Spain)
  • James Gleeson (Department of Mathematics and Statistics, University of Limerick, Ireland)
  • Javier Borge-Holthoefer (Internet Interdisciplinary Institute, IN3, Group: CoSIN3, Barcelona, Catalunya)-External BIFI Member.
  • Sandra González-Bailón (University of Pennsylvania’s Annenberg School for Communication, US)
  • Francisco Rodriguez (Department of Applied Mathematics and Statistics, Institute of Mathematics and Computer Science, University of São Paulo, Brazil)
  • Angel Sánchez (Universidad Carlos III, Madrid, Spain. External BIFI Member)
  • José Cuesta (Universidad Carlos III, Madrid, Spain. External BIFI Member)
  • Josep Perelló (OpenSystems-UB, Departament de Física Fonamental, Universitat de Barcelona)
  • Matjaz Perc (Faculty of Natural Sciences and Mathematics, University of Maribor, Slovenia)
  • Zhen Wang (Research Associate, School of Computer Science and Engineering, Nanyang Technological University)


The content of this Website is protected by Copyright © law
*The reproduction of these texts is strictly prohibited. They may only be used for scientific dissemination and citing always the original source.

Physical modeling of biomolecules

Physical modeling of biomolecules

Head of the Research Line:

Pierpaolo Bruscolini


Pierpaolo Bruscolini, Researcher UZ
Antonio Rey Gayo, Full Professor UCM
Ana Ma Rubio Caparrós, Associate professor UCM
Fernando Falo Forniés, Full professor UZ
Alessandro Fiasconaro, Assistant Professor UZ



The research line of “Physical modeling of biomolecules” is articulated around three research groups which share the same approach to the physics of biological molecules as well as the analysis techniques from statistical mechanics.

Bruscolini’s Group (Statistical-physics modeling of biomolecules)

We apply coarse-grained models and statistical mechanics methods to protein folding, protein function, protein design and sequencing, seeking the best balance between the quantitative accuracy of the predictions (increasing with the complexity of the model) and the viability in terms of computational costs and time (increasing with the simplicity of the approach). We are also interested in statistical inference and its applications, and we are using it to investigate protein signaling and redesign. Our view is that the big-data explosion that we are witnessing in biology can be best tackled with simple models and approaches that help rationalizing and understanding the biological processes.

In more detail, our recent activity encompasses:

a) Protein Folding: we use simple models, and especially the Wako-Saito-Muñoz-Eaton model, for which an exact solution can be given (PRL 88, 258101, 2002). In particular, we have developed the latter model to provide a quantitative description of the experimental data from different experimental techniques.

b) Protein design and “antibody humanization”: we have applied statistical-physics methods to protein design at fixed backbone. More recently, we have worked on a statistical inference approach to the “humanization” of antibodies; that is, the design of a human antibody sequence, to be used in pre-clinical studies, starting from a murine antibody, by optimizing a “humanness score” derived by statistical inference from public databases of human sequences. Anyone interested in applications/collaborations please contact: pier “at”

c) Protein Sequencing: we are interested in de-novo sequencing by Tandem Mass Spectrometry, and have proposed a novel algorithm, T-NovoMS, which relies on the mapping of the sequencing of a parent peptide on the thermodynamics of a one-dimensional system, with the experimental MSMS spectrum acting as an external field. We have developed a web server implementing the method:

d) Protein function: we have applied a simple worm-like-chain model to explain the interplay between the catalytic and lectine domain of GalNAC transferase 2 in determining the structure and function of the enzyme.

Rey’s Group (Computer modeling of protein folding and aggregation)

We use molecular modeling and coarse-grained models to explore the stability of the native structure in proteins and the protein folding process.

We design ourselves the simulation models in order to properly analyze structural, thermodynamic and kinetic aspects of the folding process. The most important part in this design is the interaction potential, which has to permit to reach the native conformation from the unfolded state. In the last decade, we have been mainly devoted to structure-based (Gō-type) models, also known as native-centric models. They consider the contacts between amino acid pairs which can be found in the native structure, in order to define the attractive interactions which rank the different conformations sampled along the simulation. We have designed simulation models which combine these potentials with mean field interactions which take into account the chemical sequence of the considered protein sequence. In addition, we have introduced new interactions to describe the hydrogen bonds present in the protein backbone, which play a fundamental role in the study of the aggregation processes which frequently appear when proteins are kept at moderate or high concentrations.

Our main sampling technique to study the thermal transition between the native and unfolded states of the desired proteins has been an in-house designed Monte Carlo method, together with replica exchange (parallel tempering). Our research work involves a strong methodological part, by dissecting the different physicochemical contributions to the behavior of the simulation models. By using this approach to the problem, we have been able to build models which can be widely used is different systems of interest in relevant scientific and applied problems: two state proteins vs downhill folding, proteins with thermodynamic intermediates, knotted proteins, unfolding transitions under pressure, proteins in highly confined environments or the competition between folding and aggregation.

     2       3


Falo´s Group (Modeling and Simulation of biological systems)

We use mesoscopic models at different scales to address different biological problems:

a) Transport by molecular motors. We investigate the transport of cargos inside cell by processive directional motors (kinesin and dynein). Our work is focused now on the study of the influence of microtubule network structure in cargo transport and as a mean of effective motor-motor interaction.

b) Translocation of polymers. The polymers passage through membranes is an important phenomenon both from the point of view of biological and technological processes. We are developing models of translocation driven by time-dependent forces (deterministic or stochastic) to simulate a molecular motor.

c) Free energy landscape of biomolecules: we apply Markov Network Models to the process of mechanical and thermal unfolding of proteins. Besides, we are studying the modelling of the mechanical unbinding of molecular complexes in order to understand the underlying free energy profile.

d) Models of non-canonical DNA structures: We focus on the behavior of guanine quadruplex of DNA (G-quaduplex) under external forces as those generated by AFM or optical tweezers. We are developing mesoscopic models of such molecules as well as performing “all atom” molecular dynamics simulations. Our goal is to bridge the gap between the theory and simulation regimes and the realistic experimental conditions.

e) Other fields of study comprise biological water behaviour, system biology of cellular differentiation, DNA models of denaturation and unzipping.


Relevant publications

1.- Active translocation of a semifexible polymer assisted by an ATP-based molecular motor. A. Fiasconaro, J.J. Mazo and F. Falo. Scientific Reports 7, 4188 (2017).

2.-Mesoscopic model for DNA G-quadruplex unfolding A. E. Bergues-Pupo, I. Gutiérrez, J. R. Arias-Gonzalez F. Falo, and A. Fiasconaro. Scientific Reports 7, 11756 (2017).

3.-Humanization of Antibodies using a Statistical Inference Approach. A. Clavero- Álvarez , T. Di Mambro , S. Pérez-Gaviro , M. Magnani, P. Bruscolini. Scientific Reports, accepted.

4.-Steric confinement and enhanced local flexibility assist knotting in simple models of protein folding. Miguel Soler, Antonio Rey, and Patrícia FN Faísca. Phys. Chem. Chem. Phys. DOI: 10.1039/c6cp05086g (2016).

5.- Dynamic interplay between catalytic and lectin domains of GalNAc-transferases modulates protein O-glycosylation. Lira-Navarrete, E.; De Las Rivas; Compañón, I.; Pallarés, M. C.; Kong, Y.; Iglesias-Fernández, J.; Bernardes, G. J. L.; Peregrina, J. M.; Rovira, C.; Bernadó, P.; Bruscolini, P.; Clausen, H.; Lostao, A.; Corzana, F.; Hurtado-Guerrero, R. Nature Communications 6: 6937, 2015.

6.- Mapping the Topography of a Protein Energy Landscape. Hutton, R. D.; Wilkinson, J.; Faccin, M.; Sivertsson, E. M.; Pelizzola, A.; Lowe, A. R.; Bruscolini, P.; Itzhaki, L. S. J. Am. Chem Soc. 137 – 46, pp. 14610 – 14625. 2015.

7.- Active polymer translocation in the 3d domain. A. Fiasconaro, J. J. Mazo, F. Falo. Physical Review E 91, 022113 (2015).

8.- An integrative approach for modeling and simulation of Heterocyst pattern formation in Cyanobacteria filaments. A. Torres-Sánchez, J. Gómez-Gardeñes, F. Falo. PLoS Comput Biol 11(3) e1004129. (2015).

9.- Role of the central cations in the mechanical unfolding of DNA and RNA G-quadruplexes. A. E. Bergues-Pupo, J. R. Arias-Gonzalez, M. C. Morón, A. Fiasconaro and F. Falo. Nucleic Acids Research 43(15) 7638-7647 (2015).

10.- How determinant is N-terminal to C-terminal coupling for protein folding?. Heinrich Krobath, Antonio Rey, and Patrícia FN Faísca. Phys. Chem. Chem. Phys. 17, 3512 – 3524 (2015).

11.- Intermediates in the folding equilibrium of repeat proteins from the TPR family. Vicente González and Antonio Rey. Eur. Biophys. J. 43, 433 – 443 (2014).

12.- Design of a rotamer library for coarse-grained models in protein folding simulations. María Larriva and Antonio Rey. J. Chem. Inform. Model. 54, 302 – 313 (2014).

13.- Binary interactions between dendrimer molecules. A simulation study. Ana M. Rubio, Carl C. McBride and Juan J. Freire. Macromolecules, 47, 5379 – 5387 (2014).

14.- MS/MS spectra interpretation as a statistical-mechanics problem. Faccin, M.; Bruscolini, P. Analytical Chemistry. 85 – 10, pp. 4884 – 4892. 2013.

15.- Quantitative prediction of protein folding behaviors from a simple statistical model. Bruscolini, P.; Naganathan, A. N. J. Am. Chem. Soc. 133 – 14, pp. 5372 – 5379. 2011.

16.- Influence of direct motor-motor interaction in models for cargo transport by a single team of motors. S. Bouzat, F. Falo. Physical Biology 7, 046009 (2010).

17.- Computational Protein Design with Side-Chain Conformational Entropy. Sciretti, D.; Bruscolini,P.; Pelizzola,A.; Pretti,M.; Jaramillo,A. Proteins. Structure Function and Bioinformatics 74 – 1, pp. 176 – 191. 2009.

18.- Exploring the free energy landscape: From dynamics to networks and back. D. Prada-Gracia, J. Gómez-Gardeñes, P. Echenique, F. Falo. PLoS Computational Biology 5(6): e1000415. (2009) 9 pages.


Main research projects

1.-FIS2017-87519-P: Abordando la complejidad de Sistemas Sociotécnicos, Biológicos y Naturales National Project. Founding Agency: AGENCIA ESTATAL DE INVESTIGACIÓN PIs: Yamir Moreno Vega and Fernando Falo. Star-end date: 1-1-18/31-12-20.


3.- E36_17R : Grupo de Investigación: Física Estadística y No Lineal. Regional Project. Founding Agency: DIPUTACIÓN GENERAL DE ARAGÓN. PI. Juan José Mazo. 2018-2020. Start-end date: 1-1-17/ 31-12-19

4.- E30-17R : Grupo de Investigación: Supercomputación y Física de Sistemas Complejos y Biológicos (COMPHYS). Founding Agency: DIPUTACIÓN GENERAL DE ARAGÓN. PI. David Iñiguez Dieste. 2018-2020. Start-end date: 1-1-17/ 31-12-19

National Project. PI: Juan José Mazo Torres. Funding agency: MINECO. MINISTERIO DE ECONOMIA Y COMPETITIVIDAD. Start-end date: 01/01/2015 – 31/12/2017.

6.- FIS2011-25167: Redes, Biofísica y Ciencia No Lineal. National Project. PI: Juan José Mazo Torres. Funding agency: Ministerio de Ciencia e Innovación. Start-end date: 1-1-2012 / 31-7-2015.

7.- FIS2009-13364-C02-01. “ACERCAMIENTO COMPUTACIONAL A LA COMPLEJIDAD EN REDES,PROTEINAS, Y SISTEMAS DE MUCHOS AGENTES”. Ámbito geográfico: Nacional. PI: Pierpaolo Bruscolini. Funding agency: MINISTERIO DE CIENCIA E INNOVACION. Start-end date: 01/01/2010 – 31/12/2012.

8.- (FIS2008-01240).Dinámica y Estructura de Sistemas Complejos. Ámbito geográfico: Nacional. PI: Juan José Mazo Torres. Funding agency: MINISTERIO DE CIENCIA E INNOVACION. Start-end date: 2009-2011.

9.- FIS2006-12781-C02-01 “COMPLEJIDAD EN PROTEÍNAS, REDES Y SISTEMAS DE MUCHOS AGENTES”. Ámbito geográfico: Nacional. PI: Pierpaolo Bruscolini. Funding agency: MINISTERIO DE EDUCACION Y CIENCIA. Start-end date: 01/10/2006 – 30/09/2009.

10.- FIS2009-13364-C02-02 “Modelos físicos para la simulación de tránsitos conformacionales en proteínas”. National Project. PI: Antonio Rey Gayo. Funding agency:Ministerio de Ciencia e Innovación. Proyecto Start-end date: 2010- 2013.

11.- S2009/PPQ-1551: “Química a alta presión, QUIMAPRES”. Regional Project. PI: Antonio Rey Gayo (at the UCM-SIMPOL group). Valentín García Baonza (coordinator). Funding agency: Comunidad de Madrid. Start-end date: 2010 – 2013.



Bruscolini’s group:

  • Alessandro Pelizzola, Dept of Physics, Politecnico di Torino (Italy).
  • Antonio Rey, Dept. of Chemistry, UCM, Madrid (Spain).
  • Laura Itzhaki, Dept. of Pharmacy, Cambridge University (UK).
  • Javier Sancho Sanz, BIFI, Universidad de Zaragoza (Spain).
  • Ramón Hurtado, BIFI, Universidad de Zaragoza (Spain).
  • Milagros Medina, BIFI, Universidad de Zaragoza (Spain).
  • Fernando Falo Fornies, BIFI, Universidad de Zaragoza (Spain).
  • Sergio Perez Gaviro, BIFI, Centro Universitario de la Defensa, Zaragoza, Spain.
  • Mauro Magnani, University of Urbino, (Italy)

Rey’s Group

  • Patrícia F.N. Faísca (University of Lisbon, Portugal).
  • Valentín García Baonza (UCM, Dept. Química Física I).
  • Pierpaolo Bruscolini (UZ y BIFI).
  • Manuel Ángel Ramos (UCM, Dept. Matemática Aplicada).

Falo’s Group

  • Sebastián Bouzat, CEA Bariloche, Argentina).
  • Juan José Mazo (ICMA and Universidad de Zaragoza).
  • Jesús Gómez-Gardeñes (BIFI and University of Zaragoza).
  • María del Carmen Morón (ICMA and University of Zaragoza).
  • Jesús Bergues (Universidad San Jorge, Zaragoza).
  • Pierpaolo Bruscolini (UZ y BIFI).
  • Ricardo Arías González. Instituto de Nanociencia de Madrid (IMDEA).
  • Anabel Lostao. Instituto de Nanociencia de Aragón.
Molecular dynamics and electronic structure

Molecular dynamics and electronic structure

Head of the Research Line:

Jesús Clemente-Gallardo

Researchers: Permanent members

José Luis Alonso Buj (Departamento de Física Teórica, Facultad de Ciencias)
Alberto Castro Barrigón (ARAID, Edificio I+D)
Jesús Clemente Gallardo (Departamento de Física Teórica, Edificio I+D)
Fernando Falceto Blecua (Departamento de Física Teórica, Facultad de Ciencias)
Víctor Gopar (Departamento de Física Teórica, Facultad de Ciencias)
Víctor Polo Ortiz (Departamento de Química Física, Facultad de Ciencias)

Researchers: PhD students

Adrián Gómez Pueyo (Departamento de Física Teórica, Edificio I+D)
Jorge Alberto Jover Galtier (Departamento de Física Teórica, Facultad de Ciencias)



In our group, we are mainly concerned with the application of theoretical and computational tools to the study of the behavior of biological and solid state systems. Most of our methods are based on quantum mechanics and in the tools required to combine it efficiently with classical mechanical methods. We deal with many different aspects, from the most theoretical to the most applied. The following are our main lines of work:

A.- Hybrid quantum-classical models: non-adiabatic dynamics of molecular systems

Quantum Mechanics admits a description in terms of tensorial objects defined on the space of physical states. It has been proved very helpful to describe important magnitudes such as the entanglement or the purity of a given system, while being formally analogous to the tensorial description of a classical mechanical system. Both formalisms admit Hamiltonian descriptions for the most common dynamical situations, and this fact allows a straightforward generalization to statistical systems with that type of dynamics for the microstates.

When describing molecular systems, it is common to approximate some degrees of freedom (usually those corresponding to the nuclei) and describe them as classical objects, but there are some degrees of freedom (usually part or the electronic system) which must be described as quantum objects. The tensorial description allows us to combine them in a hybrid quantum-classical description which can be made to maintain several mathematical properties and this provides us with a mathematically rigorous description of the corresponding statistical system. With these tools a firm theoretical background can be constructed to enlarge traditional adiabatic methods such as those based on the Ehrenfest formalism in order to include properties such as decoherence of the electronic dynamics.

Nonetheless, several properties of the resulting models are not fully understood yet. We are studying the main properties of the equilibrium distributions of hybrid systems, in particular in the thermodynamic limit. We are also studying their dynamical properties, in particular what concerns the decoherence of the electronic states. And finally, we are trying to understand ab initio how to model the interaction of a hybrid system with an environment by incorporating stochastic effects in the hybrid dynamics.

B.- Foundations and applications of time-dependent density-functional theory

Density-functional methods have become the most successful techniques for electronic structure calculations; the time-dependent variant is the choice for spectroscopic studies that involve excited electronic states. We develop algorithms and code for the application of these theoretical tools and we belong to the team of developers of the code Octopus (

From the code we also consider relevant applications in Physics, mostly centered on molecules and nanostructures, with focus on the irradiation of the systems with large intensity lasers, and the non-linear phenomena triggered by these usually ultra-short pulses.

C.-Analysis of chemical phenomena using DFT

In this research line, theoretical studies at the DFT level are carried out for molecular systems of interest in the fields of material science and catalysis. Working together with experimental groups of inorganic and physical chemistry, theoretical calculations are nowadays a fundamental tool for the understanding of chemical concepts. Hence, important questions regarding reaction mechanisms, spectroscopic properties or molecular orbital analysis are elucidated from first principles. The theoretical results not only provide a detailed explanation of the obtained results but serve as a guide for the rational design of new molecular species with desired properties.

D.- Quantum optimal control theory

The properties and the microscopic structure of matter can be modified at the electronic level with external knobs: most notably with ultra-fast custom-shaped intense laser pulses, but also through the modification of variables such as temperature, pressure, doping, presence of external static electric or magnetic fields, etc. We propose the use of first-principles electronic structure techniques (for example, time-dependent density-functional theory, TDDFT) in connection with control theories, to design materials, or prepare states of matter, with optimal values of selected properties. The rapidly evolving capabilities of advanced laser sources (pulse durations, intensities, shaping freedom) has triggered intense research in quantum optimal control theory and other control schemes – yet few studies have directly addressed the fast, attosecond scale movement of electrons with an ab initio technique such as TDDFT, that also promises the capacity of dealing with large systems. TDDFT offers the possibility of dealing with highly non-linear response of atomic and molecular systems, even though its accuracy and reliability are permanently improving, and still require theoretical and methodological advances.

E.- Wave transport phenomena in complex systems.

Transport of quantum and classical waves through complex systems has attracted a lot of interest from both fundamental and practical interests. For instance, transport of quantum waves, or matter waves, such as electrons and photons are at the forefront of research in condensed matter. This has been motivated by the creation of new materials such as graphene and topological insulators.

In this research line we investigate the effects of the presence of disorder such as impurities on different quantities related to the quantum transport of electrons through small samples made of graphene, topological insulator materials, and normal metals. The presence of disorder gives a random character to the transport electronic properties of the materials. Thus, in this research line we are particularly interested in studying the statistical properties of quantities, such as the conductance, shot noise, concurrence, etc.

Regarding classical waves, it turns out that many phenomena that one can observe in quantum waves (electrons) can be also seen in classical (macroscopic) systems such as microwave waveguides. This is a signature of the universality of the wave phenomena that we study. Therefore, our theoretical framework has been applied to studied the transmission of electromagnetic waves, and other quantities, in experimental microwave setups.


Relevant publications

1. Conductance of 1D quantum wires with anomalous electron-wavefunction localization. Ilias Amanatidis, Ioannis Kleftogiannis, Fernando Falceto, Victor A. Gopar. Phys. Rev. B 85, 235450 (2012).

2. Ehrenfest dynamics is purity non-preserving: A necessary ingredient for decoherence. J. L. Alonso, J. Clemente-Gallardo,J. C. Cuchí, P. Echenique, F. Falceto. Journal of Chemical Physics 137, 054106, 2012.

3. Non-adiabatic effects within a single thermally averaged potential energy surface: Thermal expansion and reaction rates of small molecules. J. L. Alonso, A. Castro, J. Clemente-Gallardo, P. Echenique, J. J. Mazo, V. Polo, A. Rubio and D. Zueco. Journal of Chemical Physics 137, 22A533, 2012.

4. Conductance through disordered graphene nanoribbons: Standard and anomalous electron localization. Ioannis Kleftogiannis, Ilias Amanatidis, Victor A. Gopar. Phys. Rev. B 88, 205414 (2013).

5. Comment on “Correlated electron-nuclear dynamics: Exact factorization of the molecular wavefunction” [J. Chem. Phys. 137, 22A530 (2012)]. J. L. Alonso J. Clemente-Gallardo P. Echenique-Robba and J. A. Jover-Galtier. J. Chem. Phys. 139, 087101, 2013.

6. Beyond Anderson Localization in 1D: Anomalous Localization of Microwaves in Random Waveguides. A. A. Fernandez-Marin, J. A. Mendez-Bermudez, J. Carbonell, F. Cervera, J. Sanchez-Dehesa, and Victor A. Gopar. Phys. Rev. Lett., 113, 233901 (2014).

7. Optimal control of high-harmonic generation by intense few-cycle pulses. Solanpaa, J.; Budagosky, J. A.; Shvetsov-Shilovski, N. I.; et al. Physical Review A 90, 053402. 2014.

8. Nonextensive thermodynamic functions in the Schrödinger-Gibbs ensemble. J. L. Alonso, A. Castro, J. Clemente-Gallardo, J. C. Cuchí, P. Echenique-Robba, J. G. Esteve and F. Falceto. Phys. Rev. E 91 022137, 2015 .

9. NH Activation of Ammonia by [{M(-OMe)(cod)}2] (M = Ir, Rh) Complexes: A DFT Study. Vélez, E.; Betoré M. P.; Casado, M. A.; Polo, V. Organometallics, 2015, 34, 3959–3966.

10. Solvent-Free Iridium-Catalyzed Reactivity of CO2 with Secondary Amines and Hydrosilanes. Julian, A.; Polo, V.; Jaseer, E. A.; Fernandez-Alvarez, F. J.; Oro, L. A. ChemCatChem 2015, 7, 3895–3902.

11. Enhancing and controlling single-atom high-harmonic generation spectra: a time-dependent density-functional scheme. A. Castro, A. Rubio, and E. K. U. Gross. Eur. Phys. J. B 88, 191 (2015).

12. Oxidative Addition of the N–H Bond of Ammonia to Iridium Bis(phosphane) Complexes: A Combined Experimental and Theoretical Study. Betoré, M. P.; Casado, M. A.; García-Orduña, P.; Lahoz, F. J.; Polo, V.; Oro, L. A. Organometallic, 2016, 35, 720-731.

13. Alkoxycarbonylation of α,β-Unsaturated Amides Catalyzed by Palladium (II) Complexes: A DFT Study of the Mechanism. Suleiman, R. K..; Polo, V.; El Ali, B. RSC Advances, 2016, 6, 8440-8448.

14. Theoretical shaping of femtosecond laser pulses for molecular photo-dissociation with control techniques based on Ehrenfest’s dynamics and time-dependent density-functional theory. Alberto Castro. ChemPhysChem 17, 1439 (2016).

15. Tailored pump-probe transient spectroscopy with time-dependent density-functional theory: controlling absorption spectra. Jessica Walkenhorst, Umberto De Giovannini,1, Alberto Castro, and Angel Rubio. Eur. Phys. J. B 89, 128 (2016).


Main research projects

1.- FP7-NMP-2011-SMALL-5: Dynamics and control in nanostructures for magnetic recording and energy applications, 7th framework programme, EU

2.- FIS2013-46159-C3-2-P. TEORÍA DE SISTEMAS HÍBRIDOS CLÁSICO-CUÁNTICOS: EQUILIBRIO, DINÁMICA Y CONTROL., 01/07/2014-31/12/2017, 45980€, IP: Alberto Castro

3.- FIS2014-61301-EXP: “Una ruta nueva en la bsqueda del funcional exacto de la teorı́a de funcionales de la densidad” (. Entidad financiadora: Ministerio de Economiı́a y Competitividad. Investigador principal: Alberto Castro. Financiación: 35000€ Duración del proyecto: 09/2015 – 08/2017





  • Juan Carlos Cuchí, Universidad de Lérida, Spain
  • Ángel Rubio, Universidad del Pais Vasco, Spain
Flavoenzymes: action mechanisms and biotechnology

Flavoenzymes: action mechanisms and biotechnology

Head of the Research Line:

Dr. Milagros Medina


Dr. Marta Martínez Júlvez
Dr. Patricia Ferreira Neila
Silvia Romero Tamayo
Ernesto Anoz Carbonel
Martha Minjárez Sáenz
Nerea Novo Huerta
Sergio Boneta Martinez


Flavoenzymes are versatile and diverse biomolecules that contribute to maintain and build cellular structures, to move molecules among subcellular compartments, to eliminate toxins, to facilitate chemical transformations in biosynthetic pathways, while are particularly in charge of cellular bioenergetics and signaling. They have as cofactors the riboflavin (RF, vitamin B2) derivatives flavin mononucleotide (FMN) and/or flavin adenine dinucleotide (FAD), which confer them unique and versatile properties. All organisms contain key flavoproteins and flavoenzymes, and many of them are becoming interesting as therapeutic targets or biotechnological tools. Understanding of their mechanisms for their biotechnological applications is the main goal of the “Flavoenzymes: mechanisms of action and biotechnology” group.

Our hypothesis is that we need to understand the molecular mechanisms of flavoenzymes with key metabolic functions in bacteria, parasites, plants or animals, in order to exploit their potential in several areas of biotechnology. In our group, we apply an interdisciplinary combination of tools of biochemistry, biophysics, cell biology and computational biology in order to investigate the mechanisms and factors that provide versatility to several flavoenzyme-dependent systems, and to use this information in new biotechnological and therapeutic applications. The group’s publications reflect our contributions in key areas in the understanding of flavoenzymes. Given the wide range of interest and methodologies in our research, we maintain close collaborations with specialists from other disciplines, both within BIFI and UNIZAR, as well as with others from outside institutions.


The human body uses more than 80 flavoenzymes that perform very different functions. If an error, mutation, occurs in one of them, it might aberrantly affect its activity, stability and/or interplay with other biomolecules, altering cellular functions. Two thirds of flavin-dependent human proteins are associated with disorders caused by allelic variants, making of interest the development of molecules able to rescue native conformation and function. To treat the cause of these diseases, it is advisable the use of personalized molecular therapies able to rescue each variant functionality. Implementation of such actions requires the in advance understanding of relationships between the pathological and molecular mechanisms in which the protein is involved. The human AIF (hAIF) family consists of three flavoproteins that have in common: i) mitochondrial localization, ii) FAD- and NADH-dependent oxidoreductase domains involved in an -albeit yet not well understood- oxidoreductase function and iii) an apoptosis-inducing activityAIF1, also known as AIF, if the family prototype.

It is a phylogenetically preserved mitochondrial flavoenzyme that shares homology with different reductase families. In intact mitochondria, AIF exhibits NADH-oxide-reductase activity through its flavin cofactor, FAD, which has been linked to the stability and biogenesis of respiratory complexes I and III through interaction with other proteins, such as CHCHD4 (homolog in humans to MIA40). However, neither their role as reductase nor the molecular mechanisms underlying the biogenesis of mitochondrial complexes are currently understood. In addition, when mitochondria receive an apoptotic stimulus, AIF is translocated to the nucleus, where it interacts with other nuclear proteins forming the degradosome that induces chromatinolysis. Likewise, the involvement of mutations in the sequence of AIF as direct responsible for various neurodegenerative disorders that lead to a decrease in the efficiency of the mitochondrial oxidative phosphorylation system has been described.
Thus, both the interest in the design of new therapies to modulate the independent caspase apoptosis and the understanding of the molecular basis of the pathologies related to single mutations in AIF has increased in recent years. This has turned AIF into a potential target for the treatment of pathological alterations (cancer or degenerative diseases) in which this protein either causes a defect or an excess of apoptosis, or alters the mitochondrial bioenergetics. The redox state of AIF influences its conformation and its monomer-dimer equilibrium and, by extension, its pro-apoptotic function because these factors modulate the interaction with other proteins.
AIF1, also known as AIF, if the family prototype. It is a phylogenetically preserved mitochondrial flavoenzyme that shares homology with different reductase families. In intact mitochondria, AIF exhibits NADH-oxide-reductase activity through its flavin cofactor, FAD, which has been linked to the stability and biogenesis of respiratory complexes I and III through interaction with other proteins, such as CHCHD4 (homolog in humans to MIA40). However, neither their role as reductase nor the molecular mechanisms underlying the biogenesis of mitochondrial complexes are currently understood. In addition, when mitochondria receive an apoptotic stimulus, AIF is translocated to the nucleus, where it interacts with other nuclear proteins forming the degradosome that induces chromatinolysis. Likewise, the involvement of mutations in the sequence of AIF as direct responsible for various neurodegenerative disorders that lead to a decrease in the efficiency of the mitochondrial oxidative phosphorylation system has been described.

Since one of the possibilities of modulating the pro-apoptotic activity of AIF is the regulation of its redox activity, it is imperative to answer multiple questions related to the mechanism and cellular function of the mentioned activity.

Our research on AIF aims to provide some of these responses by understanding the molecular mechanisms linked to the processes of oxide-reduction of AIF and the consequences of its redox activity in the interaction with its protein partners in the different cellular compartments where it’s found


Pathogenic bacteria usually survive in intracellular niches deprived of nutrients when infecting mammalian cells. Among those nutrients is RF, a vitamin precursor of FMN and FAD.

Thus, flavins and flavoproteins are relevantcomponents in regulating bacterial physiology against stresses and in the biology of host cells infection and colonization, closely relating flavin metabolism of pathogens to survive in situations of oxidative stress, microanaerobic conditions or lack of nutrients. Some flavoproteins are involved in the pathogen:host cell interaction, or as primary photoreceptors in signalling response of bacterium against stress. Cell free FMN levels also control expression of particular genes through FMN riboswitches, while many pathogens are unable to acquire FMN and FAD from their environment. Finally, most flavoprotein traits and mechanisms are species-specific, and the flavoprotein content appears widely divergent. These facts make flavin homeostasis and bacterial flavoproteomes interesting targets for the development of chemotherapies to fight infectious diseases. In this context, it is important to compare flavoproteins from different species with respect to ligand binding and reaction mechanism, because, even when substrates are invariant, significant differences may prevail. Such comparative bacterial flavoenzymology will help in the development of species-specific antimicrobials. We are making our contributions by identifying potential inhibitors of a phytopathogenic bacteria reductase, as well as of the FAD synthases (FADS) from Streptococcus pneumoniae and Mycobacterium tuberculosis.


The biosynthesis of FMN and FAD, their cell homeostasis and their distribution to client apo-proteins are critical and coordinated processes for cell energy efficiency, being also of interest for their sanitary and/or biotechnological (providing FMN and FAD to flavoproteins for heterogeneous and multi-functional biocatalysts) exploitation. Bacterial FADSs have both ATP:riboflavin kinase (RFK, FMN production) and ATP:FMN:adenylyltransferase (FMNAT, FAD production) activities, which, as we have shown, are performed by almost independent protein modules. These two activities are performed by two independent enzymes in eukaryotes, known as HsRFK and HsFADS in humans. The bacterial RFK module shows sequence and structural homology with HsRFK, while the FMNAT module neither presents sequence nor structural similarity with the proteins that produce FAD in mammals. The monofunctional mammalian FADS/FMNAT has been suggested to operate as a chaperone of the flavin, providing physical interaction with its client proteins. Such a possibility has not been explored in the bifunctional bacterial FADSs, and nothing is known regarding HsRFK. Since the understanding of the molecular bases of such canalization might be fundamental to achieve the identification of molecules inhibiting the FMN and FAD production, as well as positively or negatively modulating their canalization to the client protein, we consider further efforts must be done to clarify these facts. We are using different human and bacterial proteins to address this topic.


FAD dependent reductases of bacteria: In collaboration with the groups of Dr. E. Ceccarelli and Dr. E. Orellano of the National University of Rosario, Argentina, we are comparatively studying bacterial and plastidic enzymes. Since in many bacteria, these enzymes are essential or are involved in the response to oxidative stress, they can be treated as interesting pharmacological targets for the treatment of infections caused by pathogens.

In collaboration with Dr. M. Balsera of the Instituto de Recursos Naturales y Agrobiología de Salamanca we are kinetically evaluating the catalytic mechanisms of a variety of bacterial thioredoxin reductases that control the redox state of thioredoxins—widely occurring proteins that regulate a spectrum of enzymes by dithiol-disulfide exchange.

Oxidoreductases in human metabolism: In collaboration with Dr. A. Pey of the University of Granada, we are studying the catalytic mechanism of the antioxidant and cancer-associated NAD(H):quinone oxidoreductase 1 enzyme. Kinetic studies are being used to evaluate conformational dynamics and functional cooperativity, as well as to better understand the effects of cancer-associated single amino acid variants and post-translational modifications in this protein of high relevance in cancer progression and treatment.

Dehydrogenases and oxidases: In collaboration with Dr. Martínez, from Centro de Investigaciones Biológicas, CSIC, Madrid, we have described the catalytic mechanism of the aryl-alcohol oxidase fungal flavoenzyme and identified its key residues. This enzyme is involved in the degradation of lignin and the use of plant biomass for biotechnological purposes. The structure-function studies carried out have allowed to describe in molecular detail its stereoselective mechanism of oxidation of aromatic alcohols and aldehydes and to explore its potential for the production of enantiomers or the synthesis of biopolymers of interest. In collaboration with Dr. Nonato of the Universidade de São Paulo, Brazil, we are characterizing other flavoenzymes of biomedical interest such as the dihydroorotate dehydrogenase of Leshmania major and the L-amino acid oxidase of snake venom.

In addition, we also maintain specific collaborations with various groups with activity in the field of flavoenzymes for the kinetic and structural characterization of their flavoproteins or systems.

Keywords: flavoenzymes as biocatalysts • flavoenzyme mechanisms • rapid kinetic methods • enzymatic structure and dynamics • chemical biology


  • Production of native and mutant proteins through the use of protein engineering techniques.
  • Homologous and heterologous protein expression in different microorganisms.
  • Purification of proteins (electrophoresis, chromatographic methods, HPLC, FPLC, …).
  • Absorption spectrometry: kinetic studies in steady-state, differential spectroscopy, binding affinity, midpoint reduction potentials determination.
  • Use of transient kinetic techniques such as stopped-flow and laser flash photolysis.
  • Interaction and kinetic assays under anaerobic conditions.
  • Fluorescence spectroscopy and circular dichroism.
  • Protein crystallization and x-ray diffraction for the determination of 3D protein and protein complexes structures.
  • Advance Electron Paramagnetic Resonance related techniques (ESEEM, HYSCORE, ENDOR).
  • Isothermal Titration and Differential Scanning Calorimetries.
  • Computational Biology Methods: Docking, Molecular Dynamics and QM/MM simulations.


1.- The catalytic cycle of the antioxidant and cancer-associated human NQO1 enzyme: hydride transfer, conformational dynamics and functional cooperativity. Ernesto Anoz-Carbonell, David J. Timson, Angel Pey and Milagros Medina. Antioxidants 9, 972. DOI: 10.3390/antiox9090772 (2020)

2.- The apoptosis-inducing factor family: Moonlighting proteins in the crosstalk between mitochondria and nuclei. Nerea Novo, Patricia Ferreira, Milagros Medina M. IUBMB Life 2020,1–14. DOI: 10.1002/iub.2390 (2020)

3.- In silico discovery and biological validation of ligands of FAD synthase, a promising new antimicrobial target. Isaias Lans, Ernesto Anoz-Carbonell, Karen Palacio-Rodríguez, José Antonio Aínsa, Milagros Medina and Pilar Cossio. PLOS Computational Biology, 19. DOI: 10.1371/journal.pcbi.1007898 (2020)

4.- Human riboflavin kinase: Speciesspecific traits in the biosynthesis of the FMN cofactor. Ernesto Anoz‐Carbonell, Maribel Rivero, Victor Polo, Adrián Velázquez‐Campoy, Milagros Medina. The FASEB Journal, 34, 8. DOI: 10.1096/fj.202000566R (2020)

5.- Redox and ligand binding dependent conformational ensembles in the human apoptosis-inducing factor regulate its pro-life and cell death functions. Raquel Villanueva, Silvia Romero-Tamayo, Ruben Laplaza, Juan Martínez-Olivan, Adrián Velázquez-Campoy, Javier Sancho, Patricia Ferreira and Milagros Medina. Antioxidants and Redox Signaling. 20;30(18):2013-2029. DOI: 10.1089/ars.2018.7658. (2019)

6.- Towards the competent conformation for catalysis in the ferredoxin-NADP+ reductase from the Brucella ovis pathogen. Daniel Pérez-Amigot, Victor Taleb, Sergio Boneta, Ernesto Anoz-Carbonell, Maria Sebastián, Adrián Velázquez-Campoy, Victor Polo, Marta Martínez-Júlvez, and Milagros Medina. Biochim Biophys Acta Bioenerg, 1860:148058. DOI: 10.1016/j.bbabio.2019.148058. (2019)

7.- Self-sustained enzymatic cascade for the production of 2, 5-furandicarboxylic acid from 5-methoxymethylfurfural. Juan Carro, Elena Fernández-Fueyo, Carmen Fernández-Alonso, Javier Cañada, René Ullrich, Martin Hofrichter, Miguel Alcalde, Patricia Ferreira, Ángel T. Martínez. Biothechnology for Biofuels 11 – 1, pp. 86 (2018)

8.- Discovery of antimicrobial compounds targeting bacterial FAD synthetase. María Sebastián, Ernesto Anoz-Carbonel, Begoña Gracia, Pilar Cossio, José Antonio Aínsa, Isaías Lans and Milagros Medina. J. Enz. Inhib. Med. Chem. 33, 241-254 (2018)

9.- Identification of inhibitors targeting the Ferredoxin-NADP+ reductase from the Xanthomonas citri subsp citri phytopathogenic bacteria. Marta Martínez-Júlvez, Guillermina Goñi, Daniel Pérez- Amigot, Rubén Laplaza, Irina Ionescu, Silvana Petrocelli, María Laura Tondo, Javier Sancho, Elena Orellano and Milagros Medina. Molecules 23. E29. (2018)

10.- Multiple implications of an active site phenylalanine in the catalysis of aryl-alcohol oxidase. Juan Carro, Pep Amengual-Rigo, Ferran Sancho, Milagros Medina, Víctor Guallar, Patricia Ferreira, Ángel T. Martínez. Scientific Reports 8(1):8121 (2018)

11.- Proline dehydrogenase from Thermus thermophilus does not discriminate between FAD and FMN as cofactor. Mieke M.E.  Huijbers, Marta Martínez-Júlvez, Adrie H. Westphal, Estela Delgado-Arciniega, Milagros Medina and Willem J.H. van Berkel. Scientific Reports 7:43880 (2017)

12.- The trimer interface in the quaternary structure of the bifunctional prokaryotic FAD synthetase from Corynebacterium ammoniagenes. Ana Serrano, María Sebastián, Sonia Arilla-Luna, Silvia Baquedano, Beatriz Herguedas, Adrián Velázquez-Campoy, Marta, Martínez-Júlvez and Milagros Medina. Scientific Reports 7:404 (2017)

13.- Protein dynamics promote hydride tunnelling in substrate oxidation by aryl- alcohol oxidase. Juan Carro, Marta Martínez-Júlvez, Milagros Medina, Angel T. Martínez and Patricia Ferreira. Phys.Chem.Chem.Phys. 19, 28666 (2017)     

14.- Redox proteins of hydroxylating bacterial dioxygenases establish a regulatory cascade that prevents gratuitous induction of tetralin biodegradation genes. Laura Ledesma-García, Ana Sánchez-Azqueta, Milagros Medina, Francisca Reyes-Ramírez and Eduardo Santero. Scientific Reports 6:23848 (2016)

15.- Key residues regulating the reductase activity of the human mitochondrial apoptosis inducing factor. Raquel Villanueva, Carlos Marcuello, Alejandro Usón, M. Dolores Miramar, Maria Luisa Peleato, Ana Lostao, Santos A. Susin, Patricia Ferreira, and Milagros Medina. Biochemistry 54, 5175-5184 (2015)

16.- Dynamics of the active site architecture in plant-type Ferredoxin-NADP+ reductases catalytic complexes. Ana Sánchez-Azqueta, Daniela L. Catalano-Dupuy, Arleth López-Rivero, María Laura Tondo, Elena G. Orellano, Eduardo A. Ceccarelli, and Milagros Medina. BBA-Bioenergetics 1837, 1730-1738 (2014)

17.- Structural insights into the coenzyme mediated monomer-dimer transition of the pro-apoptotic Apoptosis Inducing Factor. Patricia Ferreira, Raquel Villanueva, Marta Martínez-Júlvez, Beatriz Herguedas, Carlos Marcuello, Patricio Fernandez-Silva, Lauriane Cabon, Juan A. Hermoso, Anabel Lostao, Santos A. Susin and Milagros Medina. Biochemistry 53, 4204-4215 (2014)

18.- Theoretical study of the mechanism of the hydride transfer between Ferredoxin NADPreductase and NADP+. The role of Tyrosine 303. Isaías Lans, Milagros Medina, Edina Rosta, Gerhard Hummer, Mireia García-Viloca, José M. Lluch, and Àngels González-Lafont. J. Am. Chem. Soc. 134, 20544-20553 (2012)



1.- Flavoenzymes in Health, Disease and Drug Discovery (FIHDADD). PID2019-103901GB-I00. Spanish Ministry of Science, Innovation and Universities. June 2020-June 2023. University of Zaragoza. Research Leaders: Milagros Medina and Patricia Ferreira.

2.- Paramagnetic Species in Catalysis Research. A Unified Approach Towards Heterogeneous, Homogeneous and Enzyme Catalysis (PARACAT). 813209. Program Marie Sklodovswka Curie, EU. Innovative Training Networks (ITN). January 2019-December 2022. Universities of Zaragoza, Turín, Amberes, Cardiff y Leipzig. Research Leader: Inés García Rubio from UNIZAR.

3.- Structural Biology Reference Group. E35_20R. Diputación General de Aragón (DGA). January 2020-December 2022. University of Zaragoza. Research Leader: Milagros Medina.

4.- Flavoenzymes: molecular mechanisms and targets, pathologies and Biotechnological applications. BIO2016-75183-P. Spanish Ministry of Economy and Competitiviness (MINECO). January 2017-December 2019. University of Zaragoza. Research Leader: Milagros Medina.

5.- Flavoenzyme dependent systems: from action mechanisms to biotechnological and sanitary applications. BIO2013-42978-P. Spanish Ministry of Economy and Competitiviness (MINECO). January2014-December 2016. University of Zaragoza. Research Leader: Milagros Medina.

6.- Seguimiento de los cambios conformacionales del dominio apoptótico del Factor de Inducción de Apoptosis Humano (hAIF) con marcaje selectivo de espín y espectroscopia de EPR. University of Zaragoza. January 2015-December 2015. Research Leader: Patricia Ferreira.

7.- Retos enzimáticos, químicos y de ingeniería para la utilización de los recursos agroforestales no alimentarios (lignocelulosa) en una bio-economía más sostenible y menos contaminante. AC2015-00008-00-00. Spanish Ministry of Economy and Competitiviness (MINECO). July 2015-March 2017. CIB-CSIC/ University of Zaragoza and others. Research Leader: Susana Camarero Fernández.

8.- Structure-based drug design for diagnosis and treatment of neurological diseases: dissecting and modulating complex function in the monoaminergic systems of the brain. COST ACTION CM1103. EU. November 2011-November 2015. University of St. Andrews+18. Research Leader: Rona Ramsay.

9.- Catalytic mechanisms in flavoenzymes: keys for their biotechnological and therapeutic use. BIO2010-14983. Spanish Ministry of Science and Innovation. January 2011- September 2014. University of Zaragoza. Research Leader: Milagros Medina.

10.- Propiedades pro-apoptótica y oxido-reductasa de la proteína mitocondrial AIF (factor de inducción de apoptosis): ¿Funciones biológicas independientes o complementarias? Ref. 212363 (JIUZ-2012-BIO-01). University of Zaragoza. January 2013-December 2013. Research Leader: Patricia Ferreira.


Enzymatic Composition and Enzymatic process for the production of 2,5-furandicarboxylic acid from 5-methoxymethylfurfural using said enzymatic composition. J. Carro, E. Fernández-Fueyo, M. Alcalde, P. Ferreira, R. Ullrich, M. Hofrichter, AT. Martinez. CSIC, University of Zaragoza and Technical University of Dresden. 2017.




Dr. Adrián Velázquez-Campoy
Dr. Pier Paolo Bruscolini
Dr. Ramón Hurtado-Guerrero
Dr. José Alberto Carrodeguas
Dr. José Antonio Aínsa
Dr. Patricio Fernández
Dr. Raquel Moreno
Dr. Víctor Polo
Dr. Javier Sancho

From other Institutions

Instituto de Nanociencia de Aragón (INA), Zaragoza
– Dr. Ana Isabel Gracia Lostao
Estación Experimental de Aula Dei, Zaragoza
– Dr. Inmaculada Yruela
Universidad Nacional de Rosario, Rosario, Argentina
– Prof. Néstor Carrillo
– Dr. Eduardo Ceccarelli
– Dr. Elena Orellano
Centro de Investigaciones Biológicas, CSIC, Madrid
– Dr. Ángel Martínez
Instituto de Bioquímica Vegetal y Fotosíntesis, Universidad de Sevilla-CSIC. Sevilla
– Dr. Manuel Hervás
– Dr. José A. Navarro
Centre de Recherche des Cordeliers, Paris, France 
– Dr. Santos Susín
Universita degli Studi di Bari
– Dr. Maria Barile
Universidad de Sao Paolo, Brasil
– Dr. M Cristina Nonato
Clark University, Biology Department, USA
– Dr. David Scott Hibbett
Universidad de Granada
– Dr. Angel L. Pey
Centro de Investigación y Tecnología Agroalimentaria de Aragón
– Dr. Pilar María Muñoz Álvaro
Instituto de Recursos Naturales y Agrobiología de Salamanca  
– Dr. Mónica Balsera Diéguez
University of Antioquia UdeA, Medellin, Colombia.
– Dr. Pilar Cossio
Universitat Jaume I, Castellón
– Dr. Vicent Moliner

Protein Misfolding and Amyloid Aggregation

Protein Misfolding and Amyloid Aggregation

Head of the Research Line:

Dr. Nunilo Cremades


José Daniel Camino, PhD student
Pablo Gracia, PhD student
Diego De la Fuente, PhD student
David Polanco Irisarri, PhD student



The phenomenon of protein misfolding and amyloid aggregation has emerged in recent years as a subject of fundamental importance in a wide range of scientific disciplines such as physics, chemistry, biology and medicine. This explosion of interest on the process of protein misfolding and amyloid aggregation has primarily resulted from the recognition that approximately 50 human diseases and disorders are associated with this process, some of them among the most common and debilitating medical conditions in the modern world, including Alzheimer’s, Parkinson’s and type II diabetes.

Despite the social and economical impact of some of these diseases, little is known about the molecular origins and mechanisms of protein amyloid aggregation and its associated toxicity. The research conducted in the group led by Dr. Nunilo Cremades aims to address these fundamental questions by combining a wide range of biophysical techniques, including state-of-the-art single-molecule fluorescence, with cell biology experiments. In addition, we aim to identify new protein targets for the development of early diagnostic tools as well as more effective therapeutics for this type of diseases.


Figure 1. The development and application of single-molecule fluorescence techniques have allowed us to investigate in unprecedented detail amyloid aggregation and to discover new possible therapeutic targets. Cremades N. et al. Cell 2012.


Figure 2. We have recently been able to purify highly stable amyloid oligomeric species of alpha-synuclein, the protein whose aggregation and deposition is linked with the development of Parkinson’s disease, and show that these species are highly cytotoxic and have general properties common to other amyloid oligomers. By combining a wide range of biophysical methods with cryo-EM image reconstruction techniques we were able to obtain three-dimensional structural models and reveal the quaternary structural architectures of toxic alpha-synuclein amyloid oligomers. Chen SW. et al. PNAS USA 2015.

Figure 3. We have studied the role of water in the modulation of the energetic barrier of alpha-synuclein amyloid nucleation, and we have found that under conditions of low water activity homogeneous nucleation is orders of magnitude faster than homogeneous or heterogeneous nucleation under high protein hydration conditions. We have also demonstrated that alpha-synuclein homogeneous nucleation is highly favoured in the interior of biomolecular condensates, also referred to as membrane-less organelles, generated by liquid-liquid phase separation. Camino et al. Chem. Sci. 2020.


The group is also interested in understanding the role of intrinsic structural disorder in protein function and disease.


Figure 4. Typical phase diagram for a natively ordered and intrinsically disordered protein. We recently characterised the energy landscape of a naturally occurring intrinsically disordered enzyme. This protein acquires a wide range of molten globule-like, pre-molten globule-like and random coil-like conformations depending on the solution conditions. Zambelli and Cremades et al. Mol. Biosyst. 2012.


Figure 5. The presence of a specific type of structural disorder (molten globule-like conformations) correlates with the ability of human lysozyme to form amyloid fibrils (I56T mutant responsible for a hereditary form of systemic amyloidosis). Dhulesia and Cremades et al. J. Am. Chem. Soc. 2010


Relevant publications

1.The extent of protein hydration dictates the preference for heterogeneous or homogeneous nucleation generating either parallel or antiparallel b-sheet alpha-synuclein aggregates. Camino J.D., Gracia P., Chen S.W., Sot J., de la Arada I., Sebastián V., Arrondo J.L.R., Goñi F.M., Dobson C.M. and Cremades N. Chem. Sci. 2020, DOI: 10.1039/D0SC05297C.

2. Defining alpha-synuclein species responsible for Parkinson´s disease phenotypes in mice. Froula J.M., Castellana-Cruz M., Anabtawi N.M., Camino J.D., Chen S.W., Thrasher D.R., Freire J., Yazdi AA., Fleming S., Dobson C.M., Kumita J.R., Cremades N. (co-corresponding author) and Volpicelli-Daley L.A. J. Biol. Chem. 2019, 294(27):10392-10406. Selected paper of the year in JBC in the area “Molecular basis of disease”.

3. The contribution of biophysical and structural studies of protein self-assembly to the design of therapeutic strategies for amyloid diseases. Cremades N. and Dobson C.M. Neurobiol. Dis. 2018, 109:178-190.

4. Structural basis of membrane disruption and cellular toxicity by alpha-synuclein oligomers. Fusco G., Chen S.W., Williamson P.T.F., Cascella R., Perni M., Jarvis J.A., Cecchi C., Vendruscolo M., Chiti F., Cremades N., Ying L., Dobson C.M. and De Simone A. Science 2017, 358(6369): 1440-1443.

5. Inhibition of alpha-synuclein fibril elongation by Hsp70 is governed by a kinetic binding competition between alpha-synuclein species. Aprile F.A., Arosio P., Fusco G., Chen S.W., Kumita J.R., Dhulesia A., Tortora P., Knowles T.P., Vendruscolo M., Dobson C.M. and Cremades N. Biochemistry 2017, 56(9): 1177-1180.

6. Amyloid-b and a-synuclein decreases the level of metal-catalyzed reactive oxygen species by radical scavenging and redox silencing. Pedersen J.T., Chen S.W., Borg C.B., Ness S., Bahl J.M., Heegaard N.H., Dobson C.M., Hemmingsen L., Cremades N. (co-corresponding author) and Teilum K. J. Am. Chem. Soc. 2016, 138(12):3966-9.

7. Alpha-synuclein oligomers interact with metal ions to induce oxidative stress and neuronal death in Parkinson’s disease. Deas E., Cremades N., Angelova P.R., Ludtmann M.H., Yao Z., Chen S.W., Horrocks M.H., Banushi B., Little D., Devine M.J., Gissen P., Klenerman D., Dobson C.M., Wood N.W., Gandhi S. & Abramov AY. Antioxid. Redox Signal. 2016, 24(7):376-91.

8. Structural characterization of toxic oligomers that are kinetically trapped during α-synuclein fibril formation. Chen S.W., Drakulic S., Deas E., Ouberai M., Aprile F.A., Arranz R., Ness S., Roodveldt C., Guilliams T., De-Genst E.J., Klenerman D., Wood N.W., Knowles T.P.J., Alfonso C., Rivas G., Abramov A.Y., Valpuesta J.M., Dobson C.M. and Cremades N. Proc. Natl. Acad. Sci U.S.A. 2015, 112(16):E1994-2003.

9. Direct observation of the interconversion of normal and toxic forms of alpha-synuclein. Cremades N., Cohen S.I., Deas E., Abramov A.Y., Chen A.Y., Orte A., Sandal M., Clarke R.W., Dunne P., Aprile F.A., Bertoncini C.W., Wood N.W., Knowles T.P.J., Dobson C.M. & Klenerman D. Cell 2012, 149(5): 1048-59.


Main research projects

Main current research projects

  • PGC2018-096335-B-100: “Understanding Parkinson´s disease alpha-synuclein amyloid formation in the cell: From single-molecule biophysics in cell-like test tubes to organs-on-chips”. Sponsor: Ministerio Español de Ciencia, Innovación y Universidades, Agencia Estatal de Investigación / EU (Fondos FEDER). IP: Dra. Nunilo Cremades y Dr. José A. Carrodeguas.1/01/2019 – 31/12/2022
  • UZCUD2020-BIO-01: “Caracterización de las primeras interacciones entre moléculas de alfa-sinucleína que llevan a la formación de agregados tóxicos involucrados en la enfermedad de Párkinson”. Sponsor: Universidad de Zaragoza. IP: Dra. Nunilo Cremades y Dra. Inés García. 1/10/2020 – 30/09/2021



Laura Volpicelli-Daley (University of Alabama, USA)
Fabrizio Chiti (University of Florence, Italy)
Janet Kumita (University of Cambridge, UK)
Douglas V. Laurents (IQFR-CSIC, Spain)
Felix Goñi (Basque Centre for Biophysics, Spain)
Arturo Muga (Basque Centre for Biophysics, Spain)
Paola Picotti (ETH Zurich, Switzerland)



Dra. Nunilo Cremades

Instituto de Biocomputación y Física de Sistemas Complejos
Universidad de Zaragoza
Mariano Esquillor S/N. Edificio I+D+i
Zaragoza 50018, Spain
Teléfono: +34-876555417

Clinical Diagnosis and Drug Delivery

Clinical Diagnosis and Drug Delivery

Head of the Research Line:

Olga Abian Franco


Rafael Clavería, Predoctoral Student
María Arruebo, Predoctoral Student
Alberto Rodrigo, Predoctoral Student
Arturo Vinuesa, Predoctoral Student



Clinical Diagnosis

Differential Scanning Calorimetry (DSC) has recently emerged as a promising technique that provides useful information about serum/plasma interactomics (composition in proteins and metabolites, as well as their interactions). As a result of the illness, serum/plasma composition is altered and it is possible to discriminate between healthy individuals and patients with certain diseases.

A previous study performed in our group (25 healthy subjects and 60 gastric adenocarcinoma patients) showed that there were significant differences in the variables obtained from the calorimetric profiles of healthy and gastric adenocarcinoma patients and also between patients with different disease stages.


The validation and implementation of a methodology (DIGCAL) as a useful tool in diagnosing and monitoring relevant tumor pathologies (pancreatic ductal adenocarcinoma, preneoplasic pancreatic cystic lesions and stomach cancer), as well as in the isolation and identification of potential tumor biomarkers, are proposed in this project.

Calorimetric profiles of serum from healthy subjects and patients with the three types of cancer at different stages will be obtained before and after 6 months of clinical treatment. The multiparametric analysis of the thermal profiles will allow us to establish a clinical protocol for: 1/ screening a certain tumor process; 2/ classifying patients according to their tumour stage; 3/ monitoring the progression or control of the illness; 4/ identifying specific biomarkers for each disease type.
DIGCAL could reach the clinical level not only as a diagnosis tool, but also as a monitoring and tracking system of patients during the therapeutic treatment, adding value in prognostic or pharmacologic decisions. At the same time, identified potential biomarkers could be part of a user-friendly diagnosis kit for primary points of care.

Drug Delivery

Drug delivery is the method or process of administering a pharmaceutical compound to achieve a therapeutic effect in humans or animals. Drug delivery technologies modify drug release profile, absorption, distribution and elimination for the benefit of improving product efficacy and safety, as well as patient convenience and compliance. Drug release is from: diffusion, degradation, swelling, and affinity-based mechanisms. Most common routes of administration include the preferred non-invasive peroral (through the mouth), topical (skin), transmucosal (nasal, buccal/sublingual, vaginal, ocular and rectal) and inhalation routes. Many medications such as peptide and protein, antibody, vaccine and gene based drugs, in general may not be delivered using these routes because they might be susceptible to enzymatic degradation or cannot be absorbed into the systemic circulation efficiently due to molecular size and charge issues to be therapeutically effective. For this reason many protein and peptide drugs have to be delivered by injection or a nanoneedle array. For example, many immunizations are based on the delivery of protein drugs and are often done by injection.

5       6

Current efforts in the area of drug delivery include the development of targeted delivery in which the drug is only active in the target area of the body (for example, in cancerous tissues) and sustained release formulations in which the drug is released over a period of time in a controlled manner from a formulation. In order to achieve efficient targeted delivery, the designed system must avoid the host’s defense mechanisms and circulate to its intended site of action. Types of sustained release formulations include liposomes, drug loaded biodegradable microspheres and drug polymer conjugates. In this sense, Nanoparticles (NP) in Biomedicine represent a promising technology for drug transport and release. There are many possibilities for NP´s surface functionalization and so, many strategies for including drugs in them (allowing NP going mainly to their acting site) can be developed.

This research line represents a new strategy for including some antiviral compounds active against hepatitis C virus (HCV), that were previously developed in this group.

Several materials have been used:

1/ Cyclodextrins:

The chemical structure of CDs, cyclic oligosaccharides composed of α-1,4-glycosidic-linked glycosyl residues, provides them structural and physico-chemical properties that allow their use as molecular carriers.

In their hydrophobic cavity a wide range of compounds ranging from ions to proteins can be trapped. In addition, CDs exhibit low toxicity and low immunogenicity, and they have been used in the pharmaceutical field by promoting CD-drug complexes in order to improve the absorption, distribution, metabolism, excretion, and toxicity (ADMET)-related properties of a drug (eg, solubility, stability, delivery and release, membrane permeability and absorption, toxicity). Currently, more than 30 products can be found in the market based on CD complexes.

2/ Shell Cross-Linked Polymeric Micelles

Cross-linked polymeric micelles (CLPM), formed by amphiphilic block copolymers, have been successful for biomedical applications. irst, they can fulfill those general requirements for drug delivery systems: water solubility, low toxicity, to increase the stability of the drug inside the living organisms, to facilitate cellular uptake compared to free drug, and to produce its controlled release at a specific location. Second, the amphiphilic nature of the constituent polymer results in a hydrophobic core and a hydrophilic shell that allows encapsulation of both types of drug. Third, these nanoparticles offer further stability under high dilution conditions, below the critical micellar concentration, compared to other polymeric micelles.


Indeed, cross-linking avoids its disintegration in the bloodstream and the release of the drug before reaching the target cell. Particularly, fixation of the micelle structure by light-induced covalent cross-linking, mostly employing acrylate reactive groups, represents a clean and effective procedure to prepare stable polymer micelles that can hold either water-soluble and non-water soluble molecules and transport them through the bloodstream.

3/ Block Copolymers Micelles

Polymeric drug carriers are one of the current challenges of nanomedicine. Since the concept of physical drug encapsulation within polymeric aggregates was introduced, a significant number of polymer assemblies have been identified. In particular, the construction of amphiphilic block copolymer-based drug carriers is a subject of great interest and a stimulating topic of interdisciplinary research in chemistry, biology and materials science. In aqueous media, self-assembly of amphiphilic block copolymers (BCs) to minimize energetically unfavorable hydrophobic water interactions can lead to a variety of polymeric nanostructures including especially appealing spherical micelles and vesicles.


Relevant publications

1.- Polymeric micelles from block copolymers containing 2,6-diacylaminopyridine units for encapsulation of hydrophobic drugs.
Concellón A, Clavería-Gimeno E, Velázquez-Campoy A, Abian O, Piñol M, Oriol L.RSC Advances, 2016, 6, 24066–24075.

2.- Biophysical Screening for Identifying Pharmacological Chaperones and Inhibitors against Conformational and Infectious Diseases. Velazquez-Campoy A, Sancho J, Abian O, Vega S. Curr Drug Targets. 2016 Jan 31. [Epub ahead of print]

3.- Cysteine Mutational Studies Provide Insight into a Thiol-Based Redox Switch Mechanism of Metal and DNA Binding in FurA from Anabaena sp. PCC 7120. Botello-Morte L, Pellicer S, Sein-Echaluce VC, Contreras LM, Neira JL, Abian O, Velázquez-Campoy A, Peleato ML, Fillat MF, Bes MT. Antioxid Redox Signal. 2016, 24, 173-185.

4.- On the link between conformational changes, ligand binding and heat capacity. Vega S, Abian O, Velazquez-Campoy A. Biochim Biophys Acta. 2015 Oct 14. pii: S0304-4165(15)00274-3.

5.- Shell Cross-Linked Polymeric Micelles as Camptothecin Nanocarriers for anti-HCV Therapy. Jiménez-Pardo I, González-Pastor R, Lancelot A, Claveria-Gimeno R, Velázquez-Campoy A, Abian O, Ros MB, and Sierra T. Macromol. Biosci. 2015, 15, 1381–1391.

6.- Su1990 A New Technology for the Classification of Patients With Gastric Adenocarcinoma Based on Differential Scanning Calorimetry Serum Thermograms. Vega S, Garcia-Gonzalez MA, Lanas A, Velazquez-Campoy A, Abian O.*. Gastroenterology 148(4): S-569 · April 2015

7.- Rescuing compound bioactivity in a secondary cell-based screening by using γ-cyclodextrin as a molecular carrier. Clavería-Gimeno R, Vega S, Grazu V, De la Fuente JM, Lanas A, Velazquez-Campoy A, Abian O.*. International Journal of Nanomedicine 2015, 10: 2249-2259.

8.- Deconvolution Analysis for Classifying Gastric Adenocarcinoma Patients Based on Differential Scanning Calorimetry Serum Thermograms. Vega S, Garcia-Gonzalez MA, Lanas A, Velazquez-Campoy A, Abian O.*. Scientific reports, 5:7988 (2015).

9.- Ionic liquids in water: a green and simple approach to improve activity and selectivity of lipases.
Filice M, Romero O,  Abian O, De las Rivas B and Palomo J.M. RSC Advances 4: 49115- 49122 (2014).

10.- A unified framework based on the binding polynomial for characterizing biological systems by isothermal titration calorimetry. Vega S., Abian O.*and Velazquez-Campoy A. Methods. 2014 pp: S1046-2023(14) 00316-8.

11.- Allosteric Inhibitors of the NS3 Protease From the Hepatitis C Virus. Abian O.*, Vega S., Sancho J., Velazquez-Campoy A. PLOSone 2013, 8 (7): 69773.

12.- NS3 protease from hepatitis C virus: Biophysical studies on an intrinsically disordered protein domain. Vega S., Neira J.L., Marcuello C.,  Lostao A., Abian O. * and Velazquez-Campoy A. International Journal of Molecular Sciences 2013, 14: 13282-13306.

13.- Experimental Validation of In Silico Target Predictions on Synergistic Protein Targets. Cortes-Ciriano I, Koutsoukas A, Abian O, Velazquez-Campoy A and Bender A. MedChemComm 2013, 4, 278–288.

14.- Altering the interfacial activation mechanism of a lipase by solid-phase selective chemical modification. López-Gallego F, Abian O, Guisán JM. Biochemistry. 2012 Sep 4;51(35):7028-36.

15.- Semisynthetic peptide-lipase conjugates for improved biotransformations. Romero O, Filice M, de las Rivas B, Carrasco-Lopez C, Klett J, Morreale A, Hermoso JA, Guisan JM, Abian O, Palomo JM. Chem Commun (Camb). 2012 Sep 18;48(72):9053-5.


Main research projects

Ongoing Projects:

1.- Analysis of protein/metabolites interactions in plasma serum using calorimetry: application as a quick and noninvasive diagnostic method for early detection and monitoring of tumoral digestive diseases (DIGCAL). Funding Institution: Health Institute Carlos III. From: January 2016 To: December 2018. Principal Investigator (PI): Olga Abian Franco.

2.- Validación de un nuevo método diagnóstico en suero, rápido no invasivo para detección precoz de cáncer de páncreas (PANCal). Funding Institution: Asociación Española de Gastroenterología (AEG). From: 2015 To: 2016. Principal Investigator (PI): Olga Abian Franco.

Former Projects:

3.- Adapted Nanoparticles for transport and specific release of drugs against hepatitis C virus (VHC).
Funding Institution: Health Institute Carlos III. From: 2011 To: 2013. Principal Investigator (PI): Olga Abian Franco.

4.- Implementation of in vitro e in vivo studies of anti-infectious compounds effective against Helicobacter pylori y el HCV (hepatitis C virus). Funding Institution: Health Institute Carlos III. From: 2008 To: 2011. Principal Investigator (PI): Olga Abian Franco.



Collaborations from BIFI

Dr. Adrián Velázquez-Campoy
Prof. Javier Sancho
Dr. Jose Luis Neira

Collaborators from other Institutions

James Graham Brown Cancer Center, University of Louisville, Louisville, KY, EEUU
PhD. Nichola Garbett

Instituto Aragonés de Ciencias de la Salud (I+CS), Zaragoza
Prof. Angel Lanas
Dra. Trinidad Serrano
Dra. Estela Solanas

Instituto de Ciencia de Materiales de Aragón (ICMA). Química Orgánica. Facultad de Ciencias.
Dra. Teresa Sierra
Prof. Luis Oriol
Prof. Milagros Piñol

Instituto de Nanociencia (INA), Universidad de Zaragoza, Zaragoza
Dr. Jesús Martínez de Lafuente
Dra. Valeria Grazú
Dra. Berta Saez

Universidad de San Jorge (USJ)
Prof. Victor López
Prof. Elisa Langa

Instituto de Catálisis y Petroleoquímica, CSIC, Madrid.
Dr. Jose Miguel Palomo
Dr. Fernando López Gallego
Dr. Jose Manuel Guisan

Universidad de Zaragoza
Dr. Jose Antonio Ainsa

Unidad de Investigación Traslacional, Hospital Universitario Miguel Servet, Zaragoza
Dra. Pilar Alfonso
Dra. Pilar Giraldo
Dr. Miguel Pocovi

Servicio de Microbiología-INIBIC. Complejo Hospitalario Universitario A Coruña, La Coruña
Dr. Francisco José Pérez-Llarena

Signal transduction and membrane protein therapies

Signal transduction and membrane protein therapies

Head of the Research Line:

Javier García Nafría

Personal group webpage:


Sandra Arroyo Urea, PhD student
Sofía Nadal Delgado, TFG student
Iris del Val García, TFG student



Brain membrane protein complexes

Membrane proteins represent 30% of the human genome, mediate intercellular communication and are the target for 60% of all the drugs currently used to treat disease. Membrane proteins play a especial role in the brain, where cellular communication is the basis for brain function including consciousness and memory formation. Our goal is to understand how neuronal receptors work, understanding signal detection and transduction to the intracellular millieu. We pose special emphasis on the integration of signals at the membrane by membrane protein complexes. For this purpose we use an integrated approach of structural biology techniques (Cryo-electron microscopy and X-ray crystallography), biophysics and biochemical assays. High-resolution cryo-electron microscopy currently plays a central role since we can now determine structures of membrane protein complexes that were previously intractable by X-ray crystallography.

Protein engineering tools

Membrane proteins are intrinsically unstable outside the membrane environment and to overcome this problem we develop our own protein engineering tools. We have develop a molecular cloning system that simplifies all cloning protocols (García-Nafría, Sci. Rep, 2016 and Watson, JBC, 2019) , allowing to create in a cheap, quick and simple manner an array of membrane protein constructs that can be tested for expression and stability. However, this system can be applied in any biomedical research laboratory in any type of project.

Methodologies used in our research

  • Molecular cloning.
  • Protein engineering.
  • Protein production in bacteria, insect and mammalian cells.
  • Protein purification.
  • Surface plasmon resonance (Biacore).
  • Microscale thermophoresis.
  • Macromolecular X-ray crystallography.
  • Single-particle Cryo-electron microscopy.
  • Biochemical cell assays (Western Blot, fluorescence, functional assays).


Relevant publications

1.Cryo-Electron Microscopy: Moving Beyond X-ray Crystal Structures for Drug receptors and Drug development.García-Nafría J.* and Tate CG*.  Annual Reviews in Pharmacology and Toxicology. 2020. In press. (*corresponding author)

2. Architecture of the heteromeric GluA1/2 AMPA receptor in complex with the auxiliary subunit TARP γ8. Herguedas B, Watson JF, Ho H, Cais O, García-Nafría J, Greger IH. Science. 2019 Mar 14. pii: eaav9011

3. Cryo-EM structure of the serotonin 5-HT1B receptor coupled to heterotrimeric Go. García-Nafría J., Nehme R., Edwards P., Tate CG. Nature. 2018, 558 (7711). 620-62.

4. Cryo-EM structure of the adenosine A2A receptor coupled to an engineered heterotrimeric G protein.García-Nafría J., Lee Y., Bai X., Carpenter B. & Tate CG. Elife. 2018, May 4;7. pii: e35946.

5. Structure and organization of heteromeric AMPA-type glutamate receptors. Herguedas B*, García-Nafria J*, Cais O, Fernandez-Leiro R, Ho H, Krieger J, Greger IH. Science. 2016, Apr 29;352(6285):aad3873.

6. IVA cloning: A single-tube universal cloning system exploiting bacterial In Vivo Assembly. García-Nafria J*, Watson JF, Greger IH. Scientific Reports.  2016, Jun 6;6:27459. (*corresponding author).


Main research projects

1.- BASE ESTRUCTURAL DE LA ACTIVACION Y UNION DE FARMACOS A RECEPTORES DE DOPAMINA HUMANOS. Ministerio de Ciencia, Innovación y Universidades. Convocatoria 2018 Retos I+D+I. RTI2018-095629-J-I00. 189,970 euros.



Ramón Hurtado Guerrero (BIFI, Zaragoza, Spain).
Nicole C. Kaneider (Cambridge University, Cambridge, UK).
Rafael Fernandez Leiro (CNIO, Madrid, Spain).
Silvia Hernandez Ainsa (INA, Zaragoza, Spain)
Jesus del Barrio Lasheras (INA, Zaragoza, Spain).
Jake Watson (MRC Laboratory of Molecular Biology, Cambridge, UK).
Maria José Sanchez Barrena (CSIC, Madrid, Spain).
Beatriz Herguedas Francés (BIFI, Zaragoza, Spain).

Enzyme modulation & Reaction Mechanisms

Enzyme modulation & Reaction Mechanisms

Head of the Research Line:

Pedro Merino


Rosa Matute
Marisol Matute
Mattia Ghirardello
Giulia Monciino



Our group focus its activity on pursuing new insights in to our understanding of health related biological processes. With a wide experience on asymmetric organic synthesis and well-equipped laboratories in the Faculty of Sciences, the group, working in the field of Chemical Biology is interested in the design and synthesis of small molecules ‑including glycomimetics and nitrogenated compounds- acting as key modulators and/or inhibitors of target enzymes associated to specific biological functions. We are particularly interested in glycosyltransferases and transglycosylases among other enzymes. This activity is carried out by combining a series of multidisciplinary tools and techniques including in-house developed synthetic methodologies using metal- and organic catalysis, biocomputational approaches (docking and molecular dynamics) and advanced spectroscopic techniques (i.e. STD-NMR).

Figure 1. Glycosyltransferase GalNAc-T2: Molecular dynamics study and interactions with a designed ligand of the

Figure 2. Glycosyltransferase POFUT1: SN1-like catalytic mechanism with formation of an intimate ion pair in the transition state

In addition to biocomputational methodologies the group is also interested in studying reaction mechanisms by using QM calculations and modern topological approaches such as ELF (electron localization function) and NCI (non-covalent interactions) analyses.

Figure 3. QM and topological studies of reaction mechanisms


1.- Racemic Hemiacetals as Oxygen-Centered Pronucleophiles Triggering Cascade 1,4-Addition/Michael Reaction through Dynamic Kinetic Resolution under Iminium catalysis. Reaction Development and Mechanistic Insights, Orue, A.; Uria, U.; Roca-López, D.; Delso, I.; Reyes, E.; Carrillo, L.; Merino, P.; Vicario, J. L., Chem Sci. 2017, 8, 2904-2913.

2.- Regioselectivity Change in the Organocatalytic Enantioselective (3+2) Cycloaddition with Nitrones Through Cooperative H-Bonding Catalysis/Iminium Activation, Prieto, L.; Juste-Navarro, V.; Uria, U.; Delso, I.; Reyes, E.; Tejero, T.; Carrillo, L.; Merino, P.; Vicario, J. L., Chem Eur. J. 2017, 23, 2764-2768.

3.- Chemical Approaches to Inhibitors of Isoprenoid Biosynthesis: Targeting Farnesyl and Geranylgeranyl Pyrophosphate Synthases, Merino, P.; Maiuolo, L.; Delso, I.; Algieri, V.; De Nino, A.; Tejero, T., RSC Adv. 2017, 7, 10947-10967.

4.- Glycomimetics Targeting Glycosyltransferases: Synthetic, Computational and Structural Studies of Less-Polar Conjugates, Ghirardello, M.; de las Rivas, M.; Lacetera, A.; Delso, I.; Lira-Navarrete, E.; Tejero, T.; Martin-Santamaria, S.; Hurtado-Guerrero, R.; Merino, P., Chem Eur. J. 2016, 22, 7215-7224.

5.- Azomethine Ylides from Nitrones: Using Catalytic n-BuLi for the Totally Stereoselective Synthesis of trans-2-Alkyl-3-Oxazolines, Juste-Navarro, V.; Delso, I.; Tejero, T.; Merino, P., Chem Eur. J. 2016, 22, 11527-11532.

6.- Rational design of glycomimetic compounds targeting the Saccharomyces cerevisae transglycosylase Gas2, Delso, I.; Valero-González, J.; Marca, E.; Tejero, T.; Hurtado-Guerrero, R.; Merino, P., Chem Biol. Drug Design 2016, 87, 163-170.

7.- Nucleoside Diphosphate Sugar Analogs Targeting Glycosyltransferases, Merino, P.; Delso, I.; Tejero, T.; Ghirardello, M.; Juste-Navarro, V., Asian J. Org. Chem. 2016, 1413-1427.

8.- Design and Synthesis of Dopaminergic Agonists, Matute, M. S.; Matute, R.; Merino, P., Curr Med. Chem. 2016, 23, 2790-2825.

9.- A Natural Ternary Complex Trapped in Crystal Reveals the Catalytic Mechanism of a Retaining Glycosyltransferase, Albesa-Jové, D.; Mendoza, M. F.; Rodrigo-Unzueta, A.; Gomollón-Bel, F.; Cifuente, J.; Urresti, S.; Comino, N.; Gómez, H.; Romero-García, J.; Lluch, J. M.; Sancho-Vaello, E.; Biarnés, X.; Planas, A.; Merino, P.; Masgrau, L.; Guerin, M. E., Angew. Chem. Int. Ed. 2015, 54, 9898-9902.

10.- Theoretical Elucidation of the Mechanism of the Cycloaddition between Nitrone Ylides and Electron-deficient Alkenes, Merino, P.; Tejero, T.; Diez Martinez, A., J. Org. Chem. 2014, 79, 2189-2202.


1.- Rational Design and Stereoselective Synthesis of Glycomimetics. MINECO, 2014-2016 CTQ2013-44367-C2-1-P. 151.200 EURO. PI: Pedro Merino

2.- Rational Design of Glycomimetics Inhibitors of Glycosyltransferases. MINECO, 2017-2019 CTQ2016-76155-R. 172.800 EURO. PI: Pedro Merino


Prof. Jose L. Vicario, Universidad del Pais Vasco
Prof. P. Caramella, University of Pavia, Italy
Prof. Eric Oldfield, University of Illinois, USA
Prof. Ramon Rios, University of Southampton, United Kingdom
Prof. Daan van Aalten, University of Dundee, United Kingdom
Prof. Ugo Chiacchio, University of Catania, Italy
Prof. Loredana Maiuolo, University of Calabria, Italy
Prof. Luis Domingo, Universidad de Valencia
Prof. Santos Fustero, Universidad de Valencia
Prof. Jose C. Menendez, Universidad Complutense de Madrid
Dr. Marcelo Guerin, CIC Biogune, Bilbao
Dr. Laura Masgrau, Universidad Autonoma de Barcelona
Dr. Sonsoles Martin-Santamaría, CIB, CSIC, Madrid
Dr. Ramón Hurtado-Guerrero, BIFI, Universidad de Zaragoza
Dr. Victor Polo, BIFI, Universidad de Zaragoza
Dr. Francisco Corzana, Universidad de La Rioja

Contact Us

We're not around right now. But you can send us an email and we'll get back to you, asap.

Not readable? Change text. captcha txt

Start typing and press Enter to search