Biomedical subjects

Jonathan W Essex

Publications and source records attributed to Jonathan W Essex.

11 recordsLinked to original sources

PubMed2006

Protein-ligand binding affinity predictions by implicit solvent simulations: a tool for lead optimization?

Continuum electrostatics is combined with rigorous free-energy calculations in an effort to deliver a reliable and efficient method for in silico lead optimization. The methodology is tested by calculation of the relative binding free energies of a set of inhibitors of neuraminidase, cyclooxygenase2, and cyclin-dependent kinase 2. The calculated free energies are compared to the results obtained with explicit solvent simulations and empirical scoring functions. For cyclooxygenase2, deficiencies in the continuum electrostatics theory are identified and corrected with a modified simulation protocol. For neuraminidase, it is shown that a continuum representation of the solvent leads to markedly different protein-ligand interactions compared to the explicit solvent simulations, and a reconciliation of the two protocols is problematic. Cyclin-dependent kinase 2 proves more challenging, and none of the methods employed in this study yield high quality predictions. Despite the differences observed, for these systems, the use of an implicit solvent framework to predict the ranking of congeneric inhibitors to a protein is shown to be faster, as accurate or more accurate than the explicit solvent protocol, and superior to empirical scoring schemes.

Cyclin-Dependent Kinase 2↗

PubMed2006

Three hydrolases and a transferase: comparative analysis of active-site dynamics via the BioSimGrid database.

Comparative molecular dynamics (MD) simulations enable us to explore the conformational dynamics of the active sites of distantly related enzymes. We have used the BioSimGrid (http://www.biosimgrid.org) database to facilitate such a comparison. Simulations of four enzymes were analyzed. These included three hydrolases and a transferase, namely acetylcholinesterase, outer-membrane phospholipase A, outer-membrane protease T, and PagP (an outer-membrane enzyme which transfers a palmitate chain from a phospholipid to lipid A). A set of 17 simulations were analyzed corresponding to a total of approximately 0.1 micros simulation time. A simple metric for active-site integrity was used to demonstrate the existence of clusters of dynamic conformational behaviour of the active sites. Small (i.e. within a cluster) fluctuations appear to be related to the function of an enzymatically active site. Larger fluctuations (i.e. between clusters) correlate with transitions between catalytically active and inactive states. Overall, these results demonstrate the potential of a comparative MD approach to analysis of enzyme function. This approach could be extended to a wider range of enzymes using current high throughput MD simulation and database methods.

Acetylcholinesterase↗

PubMed2005

An analysis of the feasibility of short read sequencing.

Several methods for ultra high-throughput DNA sequencing are currently under investigation. Many of these methods yield very short blocks of sequence information (reads). Here we report on an analysis showing the level of genome sequencing possible as a function of read length. It is shown that re-sequencing and de novo sequencing of the majority of a bacterial genome is possible with read lengths of 20-30 nt, and that reads of 50 nt can provide reconstructed contigs (a contiguous fragment of sequence data) of 1000 nt and greater that cover 80% of human chromosome 1.

Chromosomes, Human, Pair 1↗

PubMed2005

Grid computing and biomolecular simulation.

Biomolecular computer simulations are now widely used not only in an academic setting to understand the fundamental role of molecular dynamics on biological function, but also in the industrial context to assist in drug design. In this paper, two applications of Grid computing to this area will be outlined. The first, involving the coupling of distributed computing resources to dedicated Beowulf clusters, is targeted at simulating protein conformational change using the Replica Exchange methodology. In the second, the rationale and design of a database of biomolecular simulation trajectories is described. Both applications illustrate the increasingly important role modern computational methods are playing in the life sciences.

Biopolymers↗

PubMed2005

Grid-based dynamic electronic publication: a case study using combined experiment and simulation studies of crown ethers at the air/water interface.

Conventional paper-based scientific publications are limited in the amount of data and interaction they can provide. Simply placing an electronic copy of such a publication on the web makes for easier distribution, but more can be achieved by making use of the technology to navigate through the paper, to show the links between calculated parameters and the data, and provide access to the chain of calculated results all the way back to the raw data and ultimately the laboratory notebook. While the paper version of this publication is in a relatively conventional form, an interactive demonstrator version (available at http://epaper.combe.chem.soton.ac.uk and as an installable package) illustrates the concepts of publication@source, whereby all the figures and data presented in the paper are linked back to the original raw data together with a description of the processes by which the raw data were analysed. This level of interactivity is achieved using semantic technologies, which have the additional advantage of making the final document subsequently available and navigable by automated techniques. We present the combined information from experimental studies of surface tension and second harmonic generation (SHG) on the behaviour of benzo-15-crown-5 at the solution/air interface, together with a molecular dynamics computer simulation, to demonstrate how the simulation aids the interpretation of the SHG experiment. The adsorption isotherm was determined using SHG and fitted to a Langmuir isotherm giving Delta ads G0=26 kJ mol-1.

Air↗

PubMed2005

Mechanism and structure-activity relationships of norspermidine-based peptidic inhibitors of trypanothione reductase.

A library of polyamine-peptide conjugates based around some previously identified inhibitors of trypanothione reductase was synthesised by parallel solid-phase chemistry and screened. Kinetic analysis of library members established that subtle structural changes altered their mechanism of action, switching between competitive and non-competitive inhibition. The mode of action of the non-competitive inhibitors was investigated in detail by a variety of techniques including enzyme kinetic analysis (looking at both NADPH and trypanothione disulfide substrates), gel filtration chromatography and analytical ultracentrifugation, leading to the identification of an allosteric mode of inhibition.

Chromatography, Gel↗

PubMed2004

The parameterization and validation of generalized born models using the pairwise descreening approximation.

Generalized Born Surface Area (GBSA) models for water using the Pairwise Descreening Approximation (PDA) have been parameterized by two different methods. The first method, similar to that used in previously reported parameterizations, optimizes all parameters against the experimental free energies of hydration of organic molecules. The second method optimizes the PDA parameters to compensate only for systematic errors of the PDA. The best models are compared to Poisson-Boltzmann calculations and applied to the computation of potentials of mean force (PMFs) for the association of various molecules. PMFs present a more rigorous test of the ability of a solvation model to correctly reproduce the screening of intermolecular interactions by the solvent, than its accuracy at predicting free energies of hydration of small molecules. Models derived with the first method are sometimes shown to fail to compute accurate potentials of mean force because of large errors in the computation of Born radii, while no such difficulties are observed with the second method. Furthermore, accurate computation of the Born radii appears to be more important than good agreement with experimental free energies of solvation. We discuss the source of errors in the potentials of mean force and suggest means to reduce them. Our findings suggest that Generalized Born models that use the Pairwise Descreening Approximation and that are derived solely by unconstrained optimization of parameters against free energies of hydration should be applied to the modeling of intermolecular interactions with caution.

Journal Article↗

PubMed2004

BioSimGrid: towards a worldwide repository for biomolecular simulations.

BioSimGrid is a database for biomolecular simulations, or, a "Protein Data Bank extended in time" for molecular dynamics trajectories. We describe the implementation details: architecture, data schema, deposition, and analysis modules. We encourage the simulation community to explore BioSimGrid and work towards a common trajectory exchange format.

Acetylcholinesterase↗

PubMed2003

FDS: flexible ligand and receptor docking with a continuum solvent model and soft-core energy function.

The docking of flexible small molecule ligands to large flexible protein targets is addressed in this article using a two-stage simulation-based method. The methodology presented is a hybrid approach where the first component is a dock of the ligand to the protein binding site, based on deriving sets of simultaneously satisfied intermolecular hydrogen bonds using graph theory and a recursive distance geometry algorithm. The output structures are reduced in number by cluster analysis based on distance similarities. These structures are submitted to a modified Monte Carlo algorithm using the AMBER-AA molecular mechanics force field with the Generalized Born/Surface Area (GB/SA) continuum model. This solvent model is not only less expensive than an explicit representation, but also yields increased sampling. Sampling is also increased using a rotamer library to direct some of the protein side-chain movements along with large dihedral moves. Finally, a softening function for the nonbonded force field terms is used, enabling the potential energy function to be slowly turned on throughout the course of the simulation. The docking procedure is optimized, and the results are presented for a single complex of the arabinose binding protein. It was found that for a rigid receptor model, the X-ray binding geometry was reproduced and uniquely identified based on the associated potential energy. However, when side-chain flexibility was included, although the X-ray structure was identified, it was one of three possible binding geometries that were energetically indistinguishable. These results suggest that on relaxing the constraint on receptor flexibility, the docking energy hypersurface changes from being funnel-like to rugged. A further 14 complexes were then examined using the optimized protocol. For each complex the docking methodology was tested for a fully flexible ligand, both with and without protein side-chain flexibility. For the rigid protein docking, 13 out of the 15 test cases were able to find the experimental binding mode; this number was reduced to 11 for the flexible protein docking. However, of these 11, in the majority of cases the experimental binding mode was not uniquely identified, but was present in a cluster of low energy structures that were energetically indistinguishable. These results not only support the presence of a rugged docking energy hypersurface, but also suggest that it may be necessary to consider the possibility of more than one binding conformation during ligand optimization.

Algorithms↗

PubMed2002

Fluoride-selective binding in a new deep cavity calix[4]pyrrole: experiment and theory.

A new "super-extended cavity" tetraacetylcalix[4]pyrrole derivative was synthesized and characterized, and X-ray crystal structures of complexes bound to fluoride and acetonitrile were obtained. The binding behavior of this receptor was investigated by NMR titration, and the complex was found to exclusively bind fluoride ions in DMSO-d(6). This unusual binding behavior was investigated by Monte Carlo free energy perturbation simulations and Poisson calculations, and the ion specificity was seen to result from the favorable electrostatic interactions that the fluoride gains by sitting lower in the phenolic cavity of the receptor. The effect of water present in the DMSO on the calculated free energies of binding was also investigated. Owing to the use of a saturated ion solution, the effect of contaminating water is small in this case; however, it has the potential to be very significant at lower ion concentrations. Finally, the adaptive umbrella WHAM protocol was investigated and optimized for use in binding free energy calculations, and its efficiency was compared to that of the free energy perturbation calculations; adaptive umbrella WHAM was found to be approximately two times more efficient. In addition, structural evidence demonstrates that the protocol explores a wider conformational range than free energy perturbation and should therefore be the method of choice. This paper represents the first complete application of this methodology to "alchemical" changes.

Journal Article↗

PubMedDate not supplied

A computer-aided drug discovery system for chemistry teaching.

The Schools Malaria Project (http://emalaria.soton.ac.uk/) brings together school students with university researchers in the hunt for a new antimalaria drug. The design challenge being offered to students is to use a distributed drug search and selection system to design potential antimalaria drugs. The system is accessed via a Web interface. This e-science project displays the results of the trials in an accessible manner, giving students an opportunity for discussion and debate both with peers and with the university contacts. The project has been implemented by using distributed computing techniques, spreading computer load over a network of machines that cross institutional boundaries, forming a grid. This provides access to greater computing power and allows a much more complex and detailed formulation of the drug design problem to be tackled for research, teaching, and learning.

Animals↗