- 46 Views
- Uploaded on
- Presentation posted in: General

The Basic Technology Research Programme

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Proof of Concept Studies & Consortia Building Networks

- Cross research council endeavour
- administered by EPSRC

- Funding for research to create a new technology
- Change the way we do science
- Underpin the future industrial base

- 15 research projects funded up to April 2003
- Total funding for this period - £41M
- To support large, long term, high risk, high impact research consortia
- Encourage investigation of speculative ideas

- Two levels of funding
- One year start up
- Full grant up to five years

- Two types of start up funding
- Proof of concept
- Consortia building networking

- One year funding up to £100K
- Research to investigate feasibility of developing the new technology
- Output – a business case for the next step of investigation to be submitted in May 2004
- Basic Technology Programme
- Existing Research Council initiatives
- DTI programmes

- Involvement of the users of the new technology at a very early stage
- Funding to form networks & hold workshops

- Basic Technology Funding for October 2003 to September 2004
- Proof of concept
- Consortia building networking

- Academic partners
- University of Portsmouth
- University of Erlangen
- University of Southampton
- University of Oxford
- University of Aberdeen

- Development of techniques to describe irregular solids & surfaces
- Development of projection & pattern recognition techniques for non-planar colour-coded surfaces
- spherical harmonics, molecular topology

- Conformational analysis
- Rigid body dynamics incorporating surface features
- rigid parts of molecule treated as anisotropic solids linked by rotatable bonds

- Investigate how best to generate prediction models using surface properties that define a low dimensional chemical space
- QSAR, pattern recognition, artificial intelligence, analysis of surfaces

- Bench marking using Grid computing

- High throughput virtual docking
- Physical property mapping
- ADMET prediction
- Long time-period simulation techniques
- Crystallisation and solubility
- Prediction of tautomers
- Chemical reactivity and metabolism

Letchworth, 16th March 2004

- Molecular Surfaces and Property Calculation
- RGB Encoding & Pattern Recognition
- Conformational Analysis
- Rigid Body Molecular Dynamics
- Analysis of Variables & QSAR models
- Grid Computing
- Consortium Building

Small

Consensus Set of 74 Drug Molecules (diverse)

QSAR set (31 CoMFA steroids)

Medium

WDI subset (2,400 comps)

Harvard Chembank dataset (2,000 comps)

Large

WDI (50,000)

Maybridge (50,000)

Allopurinol

3D co-ordinates from CORINA

QM calculations with VAMP

Local Properties and surfaces from ParaSurf

Surfaces

Isodensity Surfaces

Shrink Wrap

Marching Cube

Surfaces fit to Spherical Harmonics

Properties

MEP, LIE, LEA and LP

Encoded at points on the surface

Encoded as Spherical Harmonic Expansions

Each Local Property encoded as a colour

LIE encoded on Red channel

LEA encoded on Green Channel

LP encoded on Blue Channel

Alternative Encoding

LIE

LEA

Absolute value of MEP

Efficient All Atom MD analysis (DASH)

Treated as time series (not Cluster Analysis)

Scales linearly with simulation length

No need for arbitrary choice of number of clusters

Can be analysed using Markov Chain methodology

Well founded methodology e.g. CNS / XPLOR (Axel T. Brunger, Stanford University)

Idea is to use rigid groups to model flexibility:

In the ligand

and the protein binding site.

Allows time-steps of 10fs to 20fs.

LIE

LEA

LP

MEP

LIE

1

0.44

0.26

0.39

LEA

0.44

1

0.58

0.47

LP

0.26

0.58

1

-0.1

MEP

0.39

0.47

-0.1

1

34 descriptors based on Normal Distribution

Principal Components

Spherical Harmonic Co-efficients

Maximum value of the local ionization energy

Minimum value of the local ionization energy

Mean value of the local ionization energy

Range of the local ionization energy

Variance in the local ionization energy

Moments

Order 1 – Mean

Order 2 – Variance

Order 3 – Skewness

Order 4 – Kurtosis

Overlapping Gaussians

Derived from previous work on MD analysis

Models derived from Local Properties

Surface Integral Model for Solvation Energy

RMS Error ~ 0.75 Kcal

Drug Likeness

SOMs trained on WDI (drugs) & Maybridge (general)

Parameters from PC of Local Property Descriptors

Medium sized datasets superimposed on SOMs

ParaSurf compiled on

SGI IRIX

Windows

Linux (SUSE)

IBM AIX

Future Platforms

SUN Solaris

GRID enabling at Portsmouth (Mark Baker), Southampton and Oxford.

SGI R10k, 256MB

VAMP ~ 30s/compound

ParaSurf ~ 10s/compound

Intel 1.8 Xeon/ AMD Athlon XP-2000+

ParaSurf ~ 2s/compound

SGI FUEL Workstation R14K

ParaSurf ~ 2s/compound

- Properties can be calculated
- Properties can be RGB encoded
- Properties are local
- Properties can be used for QSAR models

- Comparing and recognising 3D objects is an active research area in robotics and AI.
- Fast methods have been developed for database indexing.
- Rotationally invariant descriptors of 3D objects are possible.

- Can we recognise similar surfaces?
- Can we recognise similar surfacefragments?
- Can we identify the most similar surface to our target?
- How do we compare field descriptors on the molecular surface?

- Internal coordinates e.g. a distance matrix.
- Energy distributions based on the spherical harmonics.
- The spherical harmonic coefficients.
- Radial integration, radial scanning, and invariant moments.

Two different approaches:

- Using spherical harmonic molecular surfaces [J. Comp. Chem. 20(4) 383-395; Ritchie and Kemp 2000; University of Aberdeen].
- Partial molecular alignment via local structure analysis [J. Chem. Inf. Comput. Sci. 40(2) 503-512 ; Robinson, Lyne and Richards 1999; University of Oxford].

A grid is placed on a ParaSurf surface in order to reduce

the number of surface points from 4038 to 55.

- We do not know which points on the two surfaces need to be aligned with each other.
- The essential approach is:
all surface points on one surface are compared with all points on the other.

- For two surfaces, with M and N points, MN possible alignments are possible:
- we want to reduce this large search space!

The voting pairs can have a critical effect

on the quality of the surface alignment.

- A voting table may list all matching pairs of surface points (i.e. all possible alignments).
- A smart editing of votes within the voting table can enable speed and accuracy.
- We want to only consider alignments between similar local features on the surfaces.
- The more false votes we have in the voting table the harder it is to find the optimum alignment.

P1

The internal distance matrix

can be used to distinguish

between surface points.

P3

By comparing rows and

columns from distance matrices

of different surfaces we can

detect similar surface features.

P2

Similar local features, or interest points, on the molecular surface can be identified using a distance matrix.

For a point on each surface:

- Arrays of internal surface point distances are calculated for both points i.e. dist1[], dist2[].
- After a crude alignment, the absolute difference of dist1[] and dist2[] indicates the similarity of this pair of points.

The optimum alignment is composed of a rotation R and a translation T.

- Apply the current rotation r:
- Score the translation vectors t =p – q of all voting pairs (p,q) using a gravitational potential:
- High potentials identify clusters of similar translation vectors.
- The vector with the highest potential is the optimum translation T.

- Scoring all r gives R and T.

Translation vectors

(x,y coordinates plotted)

Some voting pairs for example rotations

4

3

2

1

A

B

A

B

A

B

A

B

- Surfaces are aligned (using a quick search method; e.g. 45º rotations).
- Best N alignments are selected.
- Each alignment is gently perturbed and optimised using the field descriptors.

- Align using the field descriptors’ values to identify suitable voting pairs:
- only match on similar field descriptors.

- Filtering can be achieved by aligning the fields separately.
- More accurate alignments can be generated by combining field values.

- Voting pairs:
- The distance between points in surface grid.
- The number of voting pairs.
- Identifying and selecting local features.
- How to represent the fields at interest points.

- Scoring:
- Scoring function to identify the correct rotation and translation (e.g. gravitational potential).
- Target function to compare different surface alignments (e.g. RMSD).

- Optimising the alignments.

Characterize the behaviour of a property

f : S

on amolecular surface S, in terms of a directed graph G on S derived from the gradient vector field

x =grad f(x)

Vertices (G) =fixed points of grad f (= critical points of f ).

Edges (G) = stable and unstable manifolds of the saddle points.

- minima
- saddles
- maxima

- Similarity
- Pattern recognition methods
- Maximal common subgraphs

- Complementarity
- Compare ligand graph with graph induced on ligand
- by receptor

- QSAR
- Topological indices

- S = Connolly Surface
- f(x) = Electrostatic Potential = ∑ q(i)/ d(x,i)
- Method
- Locate critical points of f (Newton-Raphson).
- Linearize at saddles, find eigenvectors of Hessian( f ).
- Integrate gradient vector field forward in time from 2 points on
- unstable eigenvector, backward in time from 2 points on stable
- eigenvector (Runge-Kutta).
- Integrate to boundary of Connolly surface patch, then continue
- on adjacent patch until reaching another critical point.

8 maxima

7 minima

13 saddles

#maxima – #saddles + #minima = (S) = 2

- Implementation for
S = spherical harmonic surface

f = MEP, LIE, LEA and LP

- Use images of triangulation points as starting points for Newton-Raphson search for critical points.
- Automatic differentiation.

Critical features

Dave Whitley

Pattern matching

on surfaces

Martin Swain

Molecular

surfaces

QM properties

presented on

surface

Data reduction

and QSAR

Brian Hudson

Compound

screening

Spherical harmonic

representation

Dave Ritchie

- High-throughput ligand docking
- Superimposition of ligand and a “negative” of the receptor

- Use of the fields to drive simulation
- Use of the fields to derive intermolecular forces
- Rigid-body motions – long time-step MD
- Free energy calculations

- Rapid screening using computationally fast approaches
- 3D fields – Andy Vinter

- On reduced set:
- Semi-empirical property calculations and alignments

- On most interesting molecules:
- Density-functional or ab-initio calculations and alignment

- More accurate molecular representations are used as appropriate, as resources allow