Biomedical subjects

Peter V Coveney

Publications and source records attributed to Peter V Coveney.

At least 19 recordsLinked to original sources

PubMed2006

A computational protocol for the integration of the monotopic protein prostaglandin H2 synthase into a phospholipid bilayer.

Prostaglandin H2 synthase (PGHS) synthesizes PGH2, a prostaglandin precursor, from arachidonic acid and was the first monotopic enzyme to have its structure experimentally determined. Both isozymes of PGHS are inhibited by nonsteroidal antiinflammatory drugs, an important class of drugs that are the primary means of relieving pain and inflammation. Selectively inhibiting the second isozyme, PGHS-2, minimizes the gastrointestinal side-effects. This had been achieved by the new PGHS-2 selective NSAIDs (i.e., COX-2 inhibitors) but it has been recently suggested that they suffer from additional side-effects. The design of these drugs only made use of static structures from x-ray crystallographic experiments. Investigating the dynamics of both PGHS-1 and PGHS-2 using classical molecular dynamics is expected to generate new insight into the differences in behavior between the isozymes, and therefore may allow improved PGHS-2 selective inhibitors to be designed. We describe a molecular dynamics protocol that integrates PGHS monomers into phospholipid bilayers, thereby producing in silico atomistic models of the PGHS system. Our protocol exploits the vacuum created beneath the protein when several lipids are removed from the top leaflet of the bilayer. The protein integrates into the bilayer during the first 5 ns in a repeatable process. The integrated PGHS monomer is stable and forms multiple hydrogen bonds between the phosphate groups of the lipids and conserved basic residues (Arg, Lys) on the protein. These interactions stabilize the system and are similar to interactions observed for transmembrane proteins.

Algorithms↗

PubMed2006

Stress response and structural transitions in sheared gyroidal and lamellar amphiphilic mesophases: Lattice-Boltzmann simulations.

We report on the stress response of gyroidal and lamellar amphiphilic mesophases to steady shear simulated using a bottom-up lattice-Boltzmann model for amphiphilic fluids and sliding periodic (Lees-Edwards) boundary conditions. We study the gyroid per se (above the sponge-gyroid transition, of high crystallinity) and the molten gyroid (within such a transition, of shorter-range order). We find that both mesophases exhibit shear thinning, more pronounced and at lower strain rates for the molten gyroid. At late times after the onset of shear, the skeleton of the crystalline gyroid becomes a structure of interconnected irregular tubes and toroidal rings, mostly oriented along the velocity ramp imposed by the shear, in contradistinction with free-energy Langevin-diffusion studies which yield a much simpler structure of disentangled tubes. We also compare the shear stress and deformation of lamellar mesophases with and without amphiphile when subjected to the same shear flow applied normal to the lamellae. We find that the presence of amphiphile allows (a) the shear stress at late times to be higher than in the case without amphiphile, and (b) the formation of rich patterns on the sheared interface, characterized by alternating regions of high and low curvature.

Journal Article↗

PubMed2006

Layered double hydroxide minerals as possible prebiotic information storage and transfer compounds.

One of the fundamental difficulties when considering the origin of life on Earth is the identification of an emergent system that not only replicated, but also had the capacity to undergo discrete mutation in such a way that following generations might inherit and pass on the mutation. We speculate that the layered double hydroxide (LDH) minerals are plausible candidates for a proto-RNA molecule. We describe a hypothetical LDH-like system which, when intercalated with certain anions, forms crystals with a high degree of internal order giving rise to novel information storage structures in which replication fidelity is maintained, a concept we use to propose an explanation for interstratification in terephthalate LDHs. The external surfaces of these hypothetical crystals provide active sites whose structure and chemistry is dictated by the internal information content of the LDH. Depending on the LDH polytype, the opposing external surfaces of a crystal may give rise to reactive sites that are either complementary or mirror images of each other, and so may be chiral. We also examine similarities between these proposed "proto-RNA" structures and the DNA that encodes the hereditary information in life today, concluding with a hypothetical scenario wherein these proto-RNA molecules predated the putative RNA-world.

Hydroxides↗

PubMed2005

Force field validation for nucleic acid simulations: comparing energies and dynamics of a DNA dodecamer.

Important questions exist regarding the quality of force fields used in molecular dynamics (MD) simulations and their interoperable use with other available MD implementations. NAMD is one of the most efficient and scalable parallel molecular dynamics codes for large-scale biomolecular simulations in the open source domain. It is the aim of this article to analyze and compare the dynamics of a benchmark DNA dodecamer d(CTTTTGCAAAAG)2 system, including its binding to a specific drug molecule arising from the use of various simulation protocols in NAMD using Amber98, with the dynamics arising from simulations of the same dodecamer using Amber98 in the AMBER package, one of the most well-established simulation codes for nucleic acids. Based upon a set of validation benchmarks, the details of which are discussed, we find that nucleic acid simulations using NAMD give meaningful results and that the essential features of the resulting dynamics are similar to those arising from the AMBER package. This sets the stage for reliable large-scale simulations of nucleic acids using NAMD.

Computer Simulation↗

PubMed2005

Modelling biological complexity: a physical scientist's perspective.

We discuss the modern approaches of complexity and self-organization to understanding dynamical systems and how these concepts can inform current interest in systems biology. From the perspective of a physical scientist, it is especially interesting to examine how the differing weights given to philosophies of science in the physical and biological sciences impact the application of the study of complexity. We briefly describe how the dynamics of the heart and circadian rhythms, canonical examples of systems biology, are modelled by sets of nonlinear coupled differential equations, which have to be solved numerically. A major difficulty with this approach is that all the parameters within these equations are not usually known. Coupled models that include biomolecular detail could help solve this problem. Coupling models across large ranges of length- and time-scales is central to describing complex systems and therefore to biology. Such coupling may be performed in at least two different ways, which we refer to as hierarchical and hybrid multiscale modelling. While limited progress has been made in the former case, the latter is only beginning to be addressed systematically. These modelling methods are expected to bring numerous benefits to biology, for example, the properties of a system could be studied over a wider range of length- and time-scales, a key aim of systems biology. Multiscale models couple behaviour at the molecular biological level to that at the cellular level, thereby providing a route for calculating many unknown parameters as well as investigating the effects at, for example, the cellular level, of small changes at the biomolecular level, such as a genetic mutation or the presence of a drug. The modelling and simulation of biomolecular systems is itself very computationally intensive; we describe a recently developed hybrid continuum-molecular model, HybridMD, and its associated molecular insertion algorithm, which point the way towards the integration of molecular and more coarse-grained representations of matter. The scope of such integrative approaches to complex systems research is circumscribed by the computational resources available. Computational grids should provide a step jump in the scale of these resources; we describe the tools that RealityGrid, a major UK e-Science project, has developed together with our experience of deploying complex models on nascent grids. We also discuss the prospects for mathematical approaches to reducing the dimensionality of complex networks in the search for universal systems-level properties, illustrating our approach with a description of the origin of life according to the RNA world view.

Animals↗

PubMed2005

Flux boundary conditions in particle simulations.

Flux boundary conditions are interesting in a number of contexts ranging from multiscale simulations to simulations of molecular hydrodynamics in nanoscale systems. Here we introduce, analyze, and test a general scheme to impose boundary conditions that simultaneously control the momentum and energy flux into open particle systems The scheme is shown to handle far from equilibrium simulations. It acquires its main characteristics from the requirement that it fulfills the second law of thermodynamics and thus minimizes the entropy production, when it is applied to reversible processes. It is shown both theoretically and through simulations that the scheme emulates the effect of an extended particle system as far as particle number fluctuations, temperature, and density profiles are concerned. The numerical scheme is further shown to be accurate and stable in both equilibrium and far from equilibrium contexts.

Journal Article↗

PubMed2005

WEDS: a Web services-based environment for distributed simulation.

Web services have the potential to radically enhance the ability of researchers to make use of distributed computing resources, but jargon and a plethora of standards make their use almost impossible for the scientist without prior experience of the necessary technologies. A powerful and simple WSRF-based middleware scheme is presented, designed to let scientists remotely deploy single or multiple instances of a pre-existing code across multiple resources, and giving steering, visualization and workflow functionality with only simple modifications to program code. It is hoped that the development and implementation of such a toolkit will be relevant not only to the problem of deploying workstation-class codes in real time, but also the move towards more tractable alternatives to the Globus toolkit for deployment of processes in a high-performance computing environment.

Computer Simulation↗

PubMed2005

Large-scale lattice Boltzmann simulations of complex fluids: advances through the advent of computational Grids.

During the last 2.5 years, the RealityGrid project has allowed us to be one of the few scientific groups involved in the development of computational Grids. Since smoothly working production Grids are not yet available, we have been able to substantially influence the direction of software and Grid deployment within the project. In this paper, we review our results from large-scale three-dimensional lattice Boltzmann simulations performed over the last 2.5 years. We describe how the proactive use of computational steering, and advanced job migration and visualization techniques enabled us to do our scientific work more efficiently. The projects reported on in this paper are studies of complex fluid flows under shear or in porous media, as well as large-scale parameter searches, and studies of the self-organization of liquid cubic mesophases.

Computer Graphics↗

PubMed2005

Hybrid molecular-continuum fluid models: implementation within a general coupling framework.

Over the past three years we have been developing a new approach for the modelling and simulation of complex fluids. This approach is based on a multiscale hybrid scheme, in which two or more contiguous subdomains are dynamically coupled together. One subdomain is described by molecular dynamics while the other is described by continuum fluid dynamics; such coupled models are of considerable importance for the study of fluid dynamics problems in which only a restricted aspect requires a fully molecular representation. Our model is representative of the generic set of coupled models whose algorithmic structure presents interesting opportunities for deployment on a range of architectures including computational grids. Here we describe the implementation of our HybridMD code within a coupling framework that facilitates flexible deployment on such architectures.

Computer Simulation↗

PubMed2005

Grid-based steered thermodynamic integration accelerates the calculation of binding free energies.

The calculation of binding free energies is important in many condensed matter problems. Although formally exact computational methods have the potential to complement, add to, and even compete with experimental approaches, they are difficult to use and extremely time consuming. We describe a Grid-based approach for the calculation of relative binding free energies, which we call Steered Thermodynamic Integration calculations using Molecular Dynamics (STIMD), and its application to Src homology 2 (SH2) protein cell signalling domains. We show that the time taken to compute free energy differences using thermodynamic integration can be significantly reduced: potentially from weeks or months to days of wall-clock time. To be able to perform such accelerated calculations requires the ability to both run concurrently and control in realtime several parallel simulations on a computational Grid. We describe how the RealityGrid computational steering system, in conjunction with a scalable classical MD code, can be used to dramatically reduce the time to achieve a result. This is necessary to improve the adoption of this technique and further allows more detailed investigations into the accuracy and precision of thermodynamic integration. Initial results for the Src SH2 system are presented and compared to a reported experimental value. Finally, we discuss the significance of our approach.

Computer Graphics↗

PubMed2005

Peptide recognition by the T cell receptor: comparison of binding free energies from thermodynamic integration, Poisson-Boltzmann and linear interaction energy approximations.

The binding to the T cell receptor of wild-type and variant HTLV-1 Tax peptide complexed to the major histocompatibility complex has been investigated by means of molecular dynamics simulations. The binding free energy difference is calculated using the molecular mechanics Poisson-Boltzmann surface area and linear interaction energy methods. These methods extract useful information on the binding energetics from simulations of the physical states of the ligands, which are more computationally expedient than the commonly used thermodynamic integration method. The successful reproduction of the relative binding free energies shows that these methods can be useful for free energy calculations and the rational design of drugs and vaccines.

Animals↗

PubMed2005

Large scale molecular dynamics simulation of native and mutant dihydropteroate synthase-sulphanilamide complexes suggests the molecular basis for dihydropteroate synthase drug resistance.

Antibiotic resistance is hampering the efficacy of drugs in the treatment of several pathological infections. Dihydropteroate synthase (DHPS) has been targeted by sulphonamide inhibitors for the past 60 years and has developed different amino acid mutations to survive sulpha drug action. We couple homology modelling techniques and massively parallel molecular dynamics simulations to study both the drug-bound and apo forms of native and mutant DHPS. Simulations of the complex between sulphanilamide and Streptomyces pneumoniae, DHPS shows how sulphanilamide is able to position itself close to 6-hydroxymethyl-7, 8-dihydropteridine-phosphate in a suitable position for the enzymatic transformation whereas in the mutant complex the sulpha drug is expelled from the catalytic site. Our simulations, therefore, provide insight into the molecular basis for drug resistance with S. pneumoniae DHPS.

Anti-Bacterial Agents↗

PubMed2005

Scientific Grid computing.

We introduce a definition of Grid computing which is adhered to throughout this Theme Issue. We compare the evolution of the World Wide Web with current aspirations for Grid computing and indicate areas that need further research and development before a generally usable Grid infrastructure becomes available. We discuss work that has been done in order to make scientific Grid computing a viable proposition, including the building of Grids, middleware developments, computational steering and visualization. We review science that has been enabled by contemporary computational Grids, and associated progress made through the widening availability of high performance computing.

Computer Simulation↗

PubMed2005

Molecular basis of peptide recognition by the TCR: affinity differences calculated using large scale computing.

Free energy calculations of the wild-type and the variant human T cell lymphotropic virus type 1 Tax peptide presented by the MHC to the TCR have been performed using large scale massively parallel molecular dynamics simulations. The computed free energy difference (-1.86 +/- 0.44 kcal/mol) using alchemical mutation-based thermodynamic integration agrees well with experimental data (-2.9 +/- 0.2 kcal/mol). Our simulations exploit state-of-the-art hardware and codes whose algorithms have been optimized for supercomputing platforms. This enables us to simulate larger, more realistic biological systems for longer durations without the imposition of artificial constraints.

Binding Sites↗

PubMed2005

Symmetry-breaking in chiral polymerisation.

We propose a model for chiral polymerisation and investigate its symmetric and asymmetric solutions. The model has a source species which decays into left- and right-handed types of monomer, each of which can polymerise to form homochiral chains; these chains are susceptible to 'poisoning' by the opposite-handed monomer. Homochiral polymers are assumed to influence the proportion of each type of monomer formed from the precursor. We show that for certain parameter values a positive feedback mechanism makes the symmetric steady-state solution unstable. The kinetics of polymer formation are then analysed in the case where the system starts from zero concentrations of monomers and chains. We show that following a long induction time, extremely large concentrations of polymers are formed for a short time, during this time an asymmetry introduced into the system by a random external perturbation may be massively amplified. The system then approaches one of the steady-state solutions described above.

Models, Chemical↗

PubMed2004

Energy controlled insertion of polar molecules in dense fluids.

We present a method to search low energy configurations of polar molecules in the complex potential energy surfaces associated with dense fluids. The search is done in the configurational space of the translational and rotational degrees of freedom of the molecule, combining steepest-descent and Newton-Raphson steps which embed information on the average sizes of the potential energy wells obtained from prior inspection of the liquid structure. We perform a molecular dynamics simulation of a liquid water shell which demonstrates that the method enables fast and energy-controlled water molecule insertion in aqueous environments. The algorithm finds low energy configurations of incoming water molecules around three orders of magnitude faster than direct random insertion. This method represents an important step towards dynamic simulations of open systems and it may also prove useful for energy-biased ensemble average calculations of the chemical potential.

Journal Article↗

PubMed2004

Hybrid molecular-continuum fluid dynamics.

We describe recent developments in the hybrid atomistic/continuum modelling of dense fluids. We discuss the general implementation of mass, momentum and energy transfers between a region described by molecular dynamics and a neighbouring domain described by the Navier-Stokes equations for unsteady continuum fluid flow.

Journal Article↗

PubMed2004

Large-scale grid-enabled lattice Boltzmann simulations of complex fluid flow in porous media and under shear.

Well-designed lattice Boltzmann codes exploit the essentially embarrassingly parallel features of the algorithm and so can be run with considerable efficiency on modern supercomputers. Such scalable codes permit us to simulate the behaviour of increasingly large quantities of complex condensed matter systems. In the present paper, we present some preliminary results on the large-scale three-dimensional lattice Boltzmann simulation of binary immiscible fluid flows through a porous medium, derived from digitized X-ray micro-tomographic data of Bentheimer sandstone, and from the study of the same fluids under shear. Simulations on such scales can benefit considerably from the use of computational steering, and we describe our implementation of steering within the lattice Boltzmann code, called LB3D, making use of the RealityGrid steering library. Our large-scale simulations benefit from the new concept of capability computing, designed to prioritize the execution of big jobs on major supercomputing resources. The advent of persistent computational grids promises to provide an optimal environment in which to deploy these mesoscale simulation methods, which can exploit the distributed nature of computer, visualization and storage resources to reach scientific results rapidly; we discuss our work on the grid-enablement of lattice Boltzmann methods in this context.

Journal Article↗