Computational Modelling Group

A composite likelihood approach to genome-wide data analyses.

Andrew Collins (Investigator), Jane Gibson, Ioannis Politopoulos

We describe composite likelihood-based analysis of a genome-wide breast cancer case-control sample by determining genome regions of fixed size on a linkage disequilibrium map which delimit comparable levels of linkage disequilibrium. Analysis of findings suggests further validation in more samples from other cohorts as well as the exploitation of novel computationally-intensive methods such as next-generation sequencing.

A habitat suitability model for predicting coral reef distributions in the Galápagos Islands

Terence Dawson (Investigator)

As part of a wider project developing a conservation strategy for the marine environment of the Galapagos Islands, this research used multi-variate modelling techniques to develop a habitat suitability prediction model for coral reefs.

Amorphous Computation, Random Graphs and Complex Biological Networks

Seth Bullock (Investigator)

This interdisciplinary research collaboration arose within the Simple Models of Complex Networks research cluster funded by the EPSRC www.epsrca.ac.uk through the Novel Computation Initiative. Here, leading groups from the Universities of Leeds, Sheffield, Nottingham, Southampton, Royal Holloway and King’s College and industrial partners BT are brought together for the first time to develop novel amorphous computation methods based on the theory of random graphs.

Antimicrobial Peptide and E. coli Membrane Interactions

Syma Khalid (Investigator), Thomas Piggot, Nils Berglund

Antimicrobial peptides (AMPs) are known to disrupt the membranes of bacterial cells such as E. coli. I work on investigating the nature of these interactions using molecular dynamics (MD) simulations.

Bioinformatic identification and physiological analysis of ethanol-related genes in C. elegans

Richard Edwards, Vincent O'Connor, Lindy Holden-Dye (Investigators), Ben Ient

Investigating the broad molecular, cellular and systems level impacts of acute and chronic ethanol in the nematode, Caenorhabditis elegans, as a model.

Cellular Automata Modelling of Membrane Formation and Protocell Evolution

Seth Bullock (Investigator), Stuart Bartlett

We simulate the meso-level behaviour of lipid-like particles in a range of chemical and physical environments. Self-organised protocellular structures can be shown to emerge spontaneously in systems with random, homogeneous initial conditions. We aim to take the current model further and explore more complex chemical scenarios in which a broad range of evolutionary behaviours will be exhibited.

Control and Prediction of the Organic Solid State

Richard Boardman

This project aims to produce a computer technology for the prediction of the crystal structure(s) of an organic molecule, that could be used even prior to the synthesis of the compound.

Such a computational study could be done relatively quickly to predict the dangers and opportunities of the solid phases of a molecule under development. Our project will develop the methods of experimental screening for polymorphs and their characterisation, and hence the combination will provide a major new technology for aiding industrial formulation.

DePuy Technology Partnership

Mark Taylor (Investigator), Adam Briscoe

This initiative concerns the transfer of knowledge between three key institutions (University of Southampton, University of Leeds and University of Hamburg) and DePuy International limited. The project is concerned with the ongoing advancement of technology used in orthopaedic devices.

Development of wide-ranging functionality in ONETEP

Chris-Kriton Skylaris (Investigator), Jacek Dziedzic

ONETEP is at the cutting edge of developments in first principles calculations. However, while the fundamental difficulties of performing accurate first-principles calculations with linear-scaling cost have been solved, only a small core of functionality is currently available in ONETEP which prevents its wide application. In this collaborative project between three Universities, the original developers of ONETEP will lead an ambitious workplan whereby the functionality of the code will be rapidly and significantly enriched.

Electrostatic embedded energy calculations of proteins, using the ONETEP DFT code

Chris-Kriton Skylaris (Investigator), Stephen Fox, Chris Pittock

Calculating the energy of a biomolecule in solvent, using quantum mechanics (QM) is possible, but extremely challenging, even with linear-scaling QM methods like ONETEP. Using electrostatic embedding, a novel twist on the existing QM/MM method is used to calculate the binding energy of a small ligand to a solvated protein, increasing the accuracy and realism of our general project work.

Flow and sedimentation processes in submarine meandering channels

Stephen Darby (Investigator)

The overall aim of this project is to generate a step-change in our understanding of the interactions between flow,
morphology & sedimentology within an active submarine channel fed by saline density currents. This central aim will be addressed through a combination of field measurements and innovative numerical modelling of gravity current morphodynamics

How far can we stretch the MARTINI?

Syma Khalid (Investigator), Ric Gillams

To date, coarse-grained lipid models have generally been parameterised to ensure the correct prediction of structural properties of membranes, such as the area per lipid and the bilayer thickness. The work described here explores the extent to which coarse-grained models are able to predict correctly bulk properties of lipids (phase behaviour) as well as the mechanical properties, such as lateral pressure profiles and stored elastic stress in bilayers. Such an evaluation is crucial for understanding the predictive capabilities of coarse-grained models.

Hybrid quantum and classical free energy methods in computational drug optimisation

Jonathan Essex, Chris-Kriton Skylaris (Investigators), Christopher Cave-Ayland

This work is based around the application of thermodynamics and quantum mechanics to the field of computational drug design and optimisation. Through the application of these theories the calculation of the physical properties of drug-like molecules is possible and hence some predictive power for their pharmaceutical activity in vivo can be obtained.

Identification of novel Crustacean Pathogen Receptor Proteins

Richard Edwards, Chris Hauton, Timothy Elliott (Investigators), Oyindamola Lawal, Lloyd Mushambadzi

We are mining EST libraries (sequence fragments of expressed genes) for novel proteins that might play a role in the immune response of crustaceans.

Identification of phage DNA, common insertion sites and their effect on genes within S.pneumoniae

Richard Edwards, Amy Dean

This study seeks to find if there are any common insertion sites across different strains of S.pneumoniae and discover genes that undergo frequent mutation due to phages and if these mutations can be linked to virulence of the strains.

Identifying factors required for DNA methylation using the imprinting control protein ZFP57

Deborah Mackay (Investigator)

Mutation of ZFP57 in humans is associated with widespread loss of DNA methylation at imprinted genes, and clinical features including congenital anomalies and developmental delay (Mackay et al, 08). This indicates that ZFP57 is required for DNA methylation of imprinted genes necessary for normal development.
We propose to identify the DNA sequences targeted by ZFP57, and its protein cofactors. This work will give insight into the biology of imprinting, indicate mechanisms of disease, and identify novel imprinted genes.

Immunotherapy Research: Modelling MHC Class I Complex Assembly

Timothy Elliott, Jorn Werner (Investigators), Alistair Bailey

This project uses mathematical modelling and simulation to investigate mechanisms by which our cells process and present biological information that is used by our immune system to distinguish between healthy and diseased cells.

Impacts of Climate and Sea-Level Change on Coastal Gullies

Stephen Darby (Investigator), Chris Hackney, Julian Leyland

Content to follow

Imprinting Disorders Finding Out Why

Karen Temple (Investigator)

We are conducting a research project to determine the cause and clinical impact of widespread imprinting aberrations in human development. We are recruiting patients with possible or definite imprinting disorders (due to methylation loss or gain at an imprinted loci)

including Silver Russell syndrome, Transient Neonatal diabetes, Beckwith Wiedemann syndrome, Angelman syndrome Prader Willi syndrome, UPD 14 syndromes and Pseudohypoparathyroidism.

Integrated in silico prediction of protein-protein interaction motifs

Richard Edwards (Investigator), Kieren Lythgow

Many vital protein-protein interactions are mediated by Short Linear Motifs (SLiMs) which are short proteins typically 5-15 amino acids long containing only a few positions crucial to function. This project integrates a number of leading computational techniques to predict novel SLiMs and add crucial detail to protein-protein interaction networks.

Interactome-wide prediction of short linear protein interaction motifs in humans

Richard Edwards (Investigator)

Short Linear Motifs (SLiMs) are important in many protein-protein interactions. In previous work, we have developed a computational tool, SLiMFinder, which places the interpretation of evidence for motifs within a statistical framework with high specificity, and subsequently enhanced sensitivity through application of conservation-based sequence masking. We are now applying these tools to a comprehensive set of human protein-protein interactions in order to predict novel human SLiMs in silico.

Lagrangian modelling of ecosystem dynamics at the Bermuda Atlantic Time-series Study station

Tom Anderson, Seth Bullock (Investigators), Melissa Saeland

Focus in the marine ecosystem modelling community is starting to shift towards the use of Lagrangian, agent-based models as these are believed to produce more realistic results. The basic assumptions behind these models have not been thoroughly tested, and this project aims to undertake a detailed study of Lagrangian marine ecosystem models, before creating one to investigate the dynamics at the Bermuda Atlantic Time-series Study station (BATS).

Lyotropic phase transitions of lipids studied by CG MD simulation and experimental techniques

Syma Khalid (Investigator), Josephine Corsi

A study of the phase behaviour of cationic lipid - DNA complexes such as those used for transfection by coarse grained molecular dynamics simulation. Lipid systems studied include DOPE, DOPE/DNA and DOPE/DOTAP/DNA. Structural parameters and phase behaviour observed computationally have been compared with those gained using Small Angle X-ray Scattering (SAXS) and polarising light microscopy techniques.

Mass Spec identification of proteins utilising EST libraries

Richard Edwards, Maria Debora Iglesias-Rodriguez (Investigators), Bethan Jones

Expressed Sequence Tag (EST) data presents a particular challenge for the identification of proteins using mass spectrometry (MS): it is often redundant (multiple copies of the same gene), consists primarily of short fragments of coding sequence, contains many sequencing errors and is generally poorly annotated. We are developing computational pipelines to maximise robust protein identifications from EST data despite these challenges.

Mathematical modelling of plant nutrient uptake

Tiina Roose (Investigator)

In this project I will describe a model of plant water and nutrient uptake and how to translate this model and experimental data from the single root scale to the root branching structure scale.

Measuring biomolecules - improvements to the spectroscopic ruler

Pavlos Lagoudakis, Tom Brown (Investigators), Jan Junis Rindermann, James Richardson

The spectroscopic ruler is a technique to measure the geometry of biomolecules on the nm scale by labeling them with pairs of fluorescent markers and measuring distance dependent non-radiative energy transfer between them. The remaining uncertainty in the application of the technique originates from the unknown orientation between the optical dipole moments of the fluorescent markers, especially when the molecule undergoes thermal fluctuations in physiological conditions. Recently we introduced a simulation based method for the interpretation of the fluorescence decay dynamics of the markers that allows us to retrieve both the average orientation and the extent of directional fluctuations of the involved dipole moments.

Modelling Macro-Nutrient Release & Fate Resulting from Sediment Resuspension in Shelf Seas

Chris Wood

This study involves adapting a previously-published model to take into account the effect resuspension events (both natural and anthropogenic) may have on nutrient dynamics at the sediment-water interface, and hence produce better estimates for the total nutrient budgets for shelf seas.

Molecular Fragments in Inhibitor Design

Jonathan Essex (Investigator), Michael Bodnarchuk

Fragment-Based Drug Discovery (FBDD) has emerged as an important tool in the drug discovery process. Instead of screening entire drug molecules, FBDD screens molecular fragments; constituents which make up drug molecules. A computational approach to identifying fragment binding is currently being sought which also yield binding free energy estimation.

Multi-objective design optimisation of coronary stents

Neil Bressloff, Georges Limbert (Investigators), Sanjay Pant

Stents are tubular type scaffolds that are deployed (using an inflatable balloon on a catheter), most commonly to recover the shape of narrowed (diseased) arterial segments. Despite the widespread clinical use of stents in cardiovascular intervention, the presence of such devices can cause adverse responses leading to fatality or to the need for further treatment. The most common unwanted responses of inflammation are in-stent restenosis and thrombosis. Such adverse biological responses in a stented artery are influenced by many factors, including the design of the stent. This project aims at using multi-objective optimisation techniques to find an optimum family of coronary stents which are more resistant to the processes of in-stent restenosis (IR) and stent thrombosis (ST).

Multi-scale simulations of bacterial outer-membrane proteins

Syma Khalid (Investigator), Jamie Parkin

Using Iridis to run multiple simulations, I aim to simulate the outer membrane proteins of Pseudomonas aeruginosa, using X-ray crystal structures of proteins only recently resolved by Bert van den Berg, University of Massachusetts. By modelling the proteins in a realistic P. aeruginosa outer membrane, I am to gain insight into the binding of these proteins to specific substrates and their function.

Multiscale modelling of biological membranes

Jonathan Essex (Investigator), Mario Orsi

Biological membranes are complex and fascinating systems, characterised by proteins floating in a sea of lipids. Biomembranes, besides being the fundamental structures employed by nature to encapsulate cells, play crucial roles in many phenomena indispensable for life, such as growth, energy storage, and in general information transduction via neural activity. In this project, we develop and apply multiscale computational models to simulate biological membranes and obtain molecular-level insights into fundamental structures and phenomena.

Multiscale Modelling of Electrochemical Processes in Neurons

John Chad (Investigator), Stuart George

Using asymptotic expansions to determine how the signalling behaviour of neurons is related to their microstructure.

Multiscale Simulation of Cellular Calcium Signalling

Hans Fangohr, Jonathan Essex (Investigators), Dan Mason

Calcium ions play a vitally important role in signal transduction and are key to many cellular processes including muscle contraction and cell apoptosis (cell death). This importance has made calcium an active area in biomedical science and mathematical modelling.

MXL Project

Mark Taylor, Junfen Shi (Investigators)

‘MXL’ is short for “Enhanced patient safety by computational Modelling from clinically available X-rays to minimise the risk of overload and instability for optimised function and Longevity”. This is an international EU-funded project which the Bioengineering Sciences Research Group at Southampton is involved in. For more information, visit http://www.m-x-l.eu

OMSys Towards a system model of a bacterial outer membrane

Syma Khalid (Investigator)

Many bacteria have an outer membrane which is the interface between the cell and its environment. The components of this membrane are well studied at an individual level, but there is a need to model and understand the outer membrane as a whole. In this project we aim to develop such a model of a bacterial outer membrane, linking computer simulations of the component molecules through to a more "systems biology" approach to modelling the outer membrane as a whole. Such an approach to modelling an OM must be multi-scale i.e. it must embrace a number of levels ranging from atomistic level modelling of e.g. the component proteins through to higher level "agent-based" modelling of the interplay of multiple components within the outer membrane as a whole. The different levels of description will be integrated to enable predictive modelling in order to explore the roles of outer membrane changes in e.g. antibiotic resistance.

Probing the oligomeric state and interaction surface of Fukutin Transmembrane Domain in lipid bilayer via Molecular Dynamics simulations

Nils Berglund

Fukutin Transmembrane Domain (FK1TMD) is localised to the endoplasmic reticulum or Golgi Apparatus within the cell where it is believed to function as a glycosyltransferase. Its localisation within the cell is thought to be mediated by the interaction of its N-terminal transmembrane domain with the lipid bilayers surrounding these compartments, each of which possess a distinctive lipid composition. Studies have revealed that the N-terminal transmembrane domain of FK1TMD exists as dimer within dilauroylphosphatidylcholine bilayers and this interaction is driven by interactions between a characteristic TXXSS motif. Furthermore residues close to the N-terminus that have previously been shown to play a key role in the clustering of lipids are shown to play a key role in anchoring the protein in the membrane.

Reconstructing past lake conditions using sediment cores

Lake sediments can be analysed for the reconstruction of past environmental conditions, and past abundances of different species. These data are the first step in the creation of a simulation model which will investigate the dramatic fluctuations in environmental conditions in the East African Rift Valley soda lakes.

Selection pressure for language and theory-of-mind in monkeys

Jason Noble (Investigator)

To what extent are the alarm calls of putty-nosed monkeys likely to be a good model for human language evolution? Simulation is used to classify evolutionary trajectories as either plausible or implausible, and to put lower bounds on the cognitive complexity required to perform particular behaviours.

Simulation modelling of habitat permeability for mammalian wildlife

Patrick Doncaster, Jason Noble (Investigators), Angela Watkins

Using and integrating least-cost models and agent-based simulations to explore the way in which mammals interact with, and hence move, through fragmented landscapes.

Spatially Embedded Complex Systems Engineering

Seth Bullock (Investigator)

SECSE brought together an interdisciplinary team of scientists working on an ambitious three-and-a-half year project titled. The research cluster spanned neuroscience, artificial intelligence, geography, and complex systems in an attempt to understand the role of spatial organization and spatial processes in complex networks within the domains of neural control, geo-information systems and distributed IT systems such as those implicated in air-traffic control.

Statistical model of the knee

Mark Taylor (Investigator), Francis Galloway, Prasanth Nair

Development of methods for large scale computational testing of a tibial tray incorporating inter-patient variability.

Sustainable domain-specific software generation tools for extremely parallel particle-based simulations

Chris-Kriton Skylaris (Investigator)

A range of particle based methods (PBM) are currently used to simulate materials in chemistry, engineering, physics and biophysics. The 4 types of PBM considered directly in the proposed are molecular dynamics (MD), the ONETEP quantum mechanics-based program, discrete element modelling (DEM), and smoothed particle hydrodynamics (SPH).
The overall research objective is to develop a sustainable tool that will deliver, in the future, cutting edge research applicable to applications ranging from dam engineering to atomistic drug design.

Tag based transcriptome analysis of gene expression in a promising green algae

Richard Edwards, Andreas Johansson

We use SuperSAGE in combination with next-generation sequencing to compare differences in gene expression between selected mutants and the wild type of a green algae. The data in the form of millions of 26 bp tags representing short stretches of expressed genes, will be analysed to find patterns of variation in gene expression under different conditions.

The autotransporter β domain: insights into structure and function through multi scale molecular dynamics simulations

Syma Khalid (Investigator), Daniel Holdbrook, Thomas Piggot

We are performing a series of molecular dynamics simulations involving all autotransporters with known structure. We aim to identify key structural and dynamic properties in this family of proteins.

The ONETEP project

Chris-Kriton Skylaris (Investigator), Stephen Fox, Chris Pittock, Alvaro Ruiz-Serrano, Jacek Dziedzic

Program for large-scale quantum mechanical simulations of matter from first principles quantum mechanics. Based on theory and algorithms we have developed for linear-scaling density functional theory calculations on parallel computers.

The Origins of Communication Revisited

Jason Noble (Investigator), Jordi Arranz

Quinn (2001) sought to demonstrate that communication be- tween simulated agents could be evolved without pre-defined communication channels. Quinn’s work was exciting because it showed the potential for ALife models to look at the real origin of communication; however, the work has never been replicated. In order to test the generality of Quinn’s result we use a similar task but a completely different agent architecture. We find that qualitatively similar behaviours emerge, but it is not clear whether they are genuinely communicative. We extend Quinn’s work by adding perceptual noise and internal state to the agents in order to promote ritualization of the nascent signal. Results were inconclusive; philosophical implications are discussed.

Tissue Engineering

Tiina Roose (Investigator)

This project deals with applying mathematical and computational modelling techniques to answer questions that are useful for tissue engineering applications.

Transgenerational inheritance of allergy in a multi generational cohort

John Holloway (Investigator)

The aim of this project is to determine the vertical transmission of DNA methylation by identification of CpG sites by microarray analysis of 450,000 CpG sites in 252 women of the IoW cohort that are associated with allergic sensitization and testing whether the identified methylation patterns are vertically transmitted to their offspring and whether modifiable environmental conditions during gestation affect DNA methylation.

Water Molecules in Protein Binding Sites

Jonathan Essex (Investigator), Michael Bodnarchuk

Water molecules are commonplace in protein binding sites, although the true location of them can often be hard to predict from crystallographic methods. We are developing tools which enable the location and affinity of water molecules to be found.

µ-VIS Computed Tomography Centre

Ian Sinclair, Richard Boardman, Dmitry Grinev, Philipp Thurner, Simon Cox, Jeremy Frey, Mark Spearing, Kenji Takeda (Investigators)

A dedicated centre for computed tomography (CT) at Southampton, providing complete support for 3D imaging science, serving Engineering, Biomedical, Environmental and Archaeological Sciences. The centre encompasses five complementary scanning systems supporting resolutions down to 200nm and imaging volumes in excess of one metre: from a matchstick to a tree trunk, from an ant's wing to a gas turbine blade.