On the role of metaheuristic optimization in bioinformatics

Metaheuristic algorithms are employed to solve complex and large-scale optimization problems in many different fields, from transportation and smart cities to finance. This paper discusses how metaheuristic algorithms are being applied to solve different optimization problems in the area of bioinformatics. While the text provides references to many optimization problems in the area, it focuses on those that have attracted more interest from the optimization community. Among the problems analyzed, the paper discusses in more detail the molecular docking problem, the protein structure prediction, phylogenetic inference, and different string problems. In addition, references to other relevant optimization problems are also given, including those related to medical imaging or gene selection for classification. From the previous analysis, the paper generates insights on research opportunities for the Operations Research and Computer Science communities in the field of bioinformatics

Performance Evaluation of Ingenious Crow Search Optimization Algorithm for Protein Structure Prediction

Protein structure prediction is one of the important aspects while dealing with critical diseases. An early prediction of protein folding helps in clinical diagnosis. In recent years, applications of metaheuristic algorithms have been substantially increased due to the fact that this problem is computationally complex and time-consuming. Metaheuristics are proven to be an adequate tool for dealing with complex problems with higher computational efficiency than conventional tools. The work presented in this paper is the development and testing of the Ingenious Crow Search Algorithm (ICSA). First, the algorithm is tested on standard mathematical functions with known properties. Then, the application of newly developed ICSA is explored on protein structure prediction. The efficacy of this algorithm is tested on a bench of artificial proteins and real proteins of medium length. The comparative analysis of the optimization performance is carried out with some of the leading variants of the crow search algorithm (CSA). The statistical comparison of the results shows the supremacy of the ICSA for almost all protein sequences

A Firefly-inspired method for protein structure prediction in lattice models

We introduce a Firefly-inspired algorithmic approach for protein structure prediction over two different lattice models in three-dimensional space. In particular, we consider three-dimensional cubic and three-dimensional face-centred-cubic (FCC) lattices. The underlying energy models are the Hydrophobic-Polar (H-P) model, the Miyazawa–Jernigan (M-J) model and a related matrix model. The implementation of our approach is tested on ten H-P benchmark problems of a length of 48 and ten M-J benchmark problems of a length ranging from 48 until 61. The key complexity parameter we investigate is the total number of objective function valuations required to achieve the optimum energy values for the H-P model or competitive results in comparison to published values for the M-J model. For H-P instances and cubic lattices, where data for comparison are available, we obtain an average speed-up over eight instances of 2.1, leaving out two extreme values (otherwise, 8.8). For six M-J instances, data for comparison are available for cubic lattices and runs with a population size of 100, where, a priori, the minimum free energy is a termination criterion. The average speed-up over four instances is 1.2 (leaving out two extreme values, otherwise 1.1), which is achieved for a population size of only eight instances. The present study is a test case with initial results for ad hoc parameter settings, with the aim of justifying future research on larger instances within lattice model settings, eventually leading to the ultimate goal of implementations for off-lattice models

A Firefly-inspired method for protein structure prediction in lattice models

We introduce a Firefly-inspired algorithmic approach for protein structure prediction over two different lattice models in three-dimensional space. In particular, we consider three-dimensional cubic and three-dimensional face-centred-cubic (FCC) lattices. The underlying energy models are the Hydrophobic-Polar (H-P) model, the Miyazawa–Jernigan (M-J) model and a related matrix model. The implementation of our approach is tested on ten H-P benchmark problems of a length of 48 and ten M-J benchmark problems of a length ranging from 48 until 61. The key complexity parameter we investigate is the total number of objective function valuations required to achieve the optimum energy values for the H-P model or competitive results in comparison to published values for the M-J model. For H-P instances and cubic lattices, where data for comparison are available, we obtain an average speed-up over eight instances of 2.1, leaving out two extreme values (otherwise, 8.8). For six M-J instances, data for comparison are available for cubic lattices and runs with a population size of 100, where, a priori, the minimum free energy is a termination criterion. The average speed-up over four instances is 1.2 (leaving out two extreme values, otherwise 1.1), which is achieved for a population size of only eight instances. The present study is a test case with initial results for ad hoc parameter settings, with the aim of justifying future research on larger instances within lattice model settings, eventually leading to the ultimate goal of implementations for off-lattice models

Data-driven modelling of biological multi-scale processes

Biological processes involve a variety of spatial and temporal scales. A holistic understanding of many biological processes therefore requires multi-scale models which capture the relevant properties on all these scales. In this manuscript we review mathematical modelling approaches used to describe the individual spatial scales and how they are integrated into holistic models. We discuss the relation between spatial and temporal scales and the implication of that on multi-scale modelling. Based upon this overview over state-of-the-art modelling approaches, we formulate key challenges in mathematical and computational modelling of biological multi-scale and multi-physics processes. In particular, we considered the availability of analysis tools for multi-scale models and model-based multi-scale data integration. We provide a compact review of methods for model-based data integration and model-based hypothesis testing. Furthermore, novel approaches and recent trends are discussed, including computation time reduction using reduced order and surrogate models, which contribute to the solution of inference problems. We conclude the manuscript by providing a few ideas for the development of tailored multi-scale inference methods.Comment: This manuscript will appear in the Journal of Coupled Systems and Multiscale Dynamics (American Scientific Publishers

Engineering protein assemblies with allosteric control via monomer fold-switching

The macromolecular machines of life use allosteric control to self-assemble, dissociate and change shape in response to signals. Despite enormous interest, the design of nanoscale allosteric assemblies has proven tremendously challenging. Here we present a proof of concept of allosteric assembly in which an engineered fold switch on the protein monomer triggers or blocks assembly. Our design is based on the hyper-stable, naturally monomeric protein CI2, a paradigm of simple two-state folding, and the toroidal arrangement with 6-fold symmetry that it only adopts in crystalline form. We engineer CI2 to enable a switch between the native and an alternate, latent fold that self-assembles onto hexagonal toroidal particles by exposing a favorable inter-monomer interface. The assembly is controlled on demand via the competing effects of temperature and a designed short peptide. These findings unveil a remarkable potential for structural metamorphosis in proteins and demonstrate key principles for engineering protein-based nanomachinery.This work was supported by the European Research Council (grant ERC-2012-ADG- 323059 to V.M.) and by the PRODESTECH network funded through the CONSOLIDER program from the Spanish Government (grant CSD2009-00088). L.A.C. acknowledges support from Ministry of Economy and Competitiveness through grants BIO2016- 78768-P and RYC-2013-13197. V.M. acknowledges additional support from the W.M. Keck Foundation and from the CREST Center for Cellular and Biomomolecular Machines (grant NSF-CREST-1547848). J.M.V. acknowledges additional support from Ministry of Economy and Competitiveness through grant BFU2016-75984. F.M.R. and A.R. thank the staff from the ALBA synchrotron (Spain) for assistance with the XALOC beamline. Structural data are deposited in the Protein Data Bank with accession codes 6QIY (X-ray CI2 classical geometry) and 6QIZ (X-ray CI2 domain swapped) and EMD- 4568 (cryo-EM CI2eng assembly)

The emergence of biofilms:Computational and experimental studies

The response of biofilms to any external stimuli is a cumulative response aggregated from individual bacteria residing within the biofilm. The organizational complexity of biofilm can be studied effectively by understanding bacterial interactions at cell level. The overall aim of the thesis is to explore the complex evolutionary behaviour of bacterial biofilms. This thesis is divided into three major studies based on the type of perturbation analysed in the study. The first study analyses the physics behind the development of mushroom-shaped structures from the influence of nutrient cues in biofilms. Glazier-Graner-Hogeweg model is used to simulate the cell characteristics. From the study, it is observed that chemotaxis of bacterial cells towards nutrient source is one of the major precursors for formation of mushroom-shaped structures. The objective of the second study is to analyse the impact of environmental conditions on the inter-biofilm quorum sensing (QS) signalling. Using a hybrid convection-diffusion-reaction model, the simulations predict the diffusivity of QS molecules, the spatiotemporal variations of QS signal concentrations and the competition outcome between QS and quorum quenching mutant bacterial communities. The mechanical effects associated with the fluid-biofilm interaction is addressed in the third study. A novel fluid-structure interaction model based on fluid dynamics and structural energy minimization is developed in the study. Model simulations are used to analyse the detachment and surface effects of the fluid stresses on the biofilm. In addition to the mechanistic models described, a separate study is carried out to estimate the computational efficiency of the biofilm simulation models