Dock Sem

Molecular Docking G. Schaftenaar

Docking Challenge Identification of the ligand’s correct binding geometry in the binding site ( Binding Mode ) Observation: Similar ligands can bind at quite different orientations in the active site.

Two main tasks of Docking Tools Sampling of conformational (Ligand) space Scoring protein-ligand complexes

Historically the first approaches. Protein and ligand fixed. Search for the relative orientation of the two molecules with lowest energy. FLOG (Flexible Ligands Oriented on Grid): each ligand represented by up to 25 low energy conformations. Rigid-body docking algorithms

Introducing flexibility: Whole molecule docking Monte Carlo methods (MC) Molecular Dynamics (MD) Simulated Annealing (SA) Genetic Algorithms (GA) Available in packages: AutoDock (MC,GA,SA) GOLD (GA) Sybyl (MD)

Monte Carlo Start with configuration A (energy E A ) Make random move to configuration B (energy E B ) Accept move when: E B < E A or if E B > E A except with probability P:

Molecular Dynamics force-field is used to calculate forces on each atom of the simulated system following Newton mechanics, calculate accelerations, velocities and new coordinates from the forces. (Force = mass times acceleration) The atoms are moved slightly with respect to a given time step

Simulated Annealing Finding a global minimium by lowering the temperature during the Monte Carlo/MD simulation

Genetic Algorithms Ligand translation, rotation and configuration variables constitute the genes Crossovers mixes ligand variables from parent configurations Mutations randomly change variables Natural selection of current generation based on fitness Energy scoring function determines fitness

Introducing flexibility: Fragment Based Methods build small molecules inside defined binding sites while maximizing favorable contacts. De Novo methods construct new molecules in the site. division into two major groups: Incremental construction (FlexX, Dock) Place & join.

Placing Fragments and Rigid Molecules All rigid-body docking methods have in common that superposition of point sets is a fundamental sub-problem that has to be solved efficiently : Geometric hashing Pose clustering Clique detection

Geometric hashing originates from computer vision Given a picture of a scene and a set of objects within the picture, both represented by points in 2d space, the goal is to recognize some of the models in the scene

Pose-Clustering For each triangle of receptor compute the transformation to each ligand matching triangle. Cluster transformations. Score the results.

Clique-Detection Nodes comprise of matches between protein and ligand Edges connect distance compatible pairs of nodes In a clique all pair of nodes are connected

Scoring Functions Shape & Chemical Complementary Scores Empirical Scoring Force Field Scoring Knowledge-based Scoring Consensus Scoring

Shape & Chemical Complementary Scores Divide accessible protein surface into zones: Hydrophobic Hydrogen-bond donating Hydrogen-bond accepting Do the same for the ligand surface Find ligand orientation with best complementarity score

Empirical Scoring Scoring parameters fit to reproduce Measured binding affinities (FlexX, LUDI, Hammerhead)

Empirical scoring Loss of entropy during binding Hydrogen-bonding Ionic interactions Aromatic interactions Hydrophobic interactions rot rot N G G G       0         bonds H neutral hb R f G . ,          . , int ionic io R f G         int arom arom R f G . ,         . . , cont lipo lipo R f G 

Force Field Scoring (Dock)             lig i prot j i ij ij ij ij nonbond r B r A E c 6 12 Nonbonding interactions (ligand-protein) : van der Waals - electrostatics Amber force field  ij j r q q

Knowledge-based Scoring Function Free energies of molecular interactions derived from structural information on Protein-ligand complexes contained in PDB Boltzmann-Like Statistics of Interatomic Contacts.       l p ref l p F P P      , exp ,  

Distribution of interatomic distances is converted into energy functions by inverting Boltzmann’s law. F P (N,O)

Potential of Mean Force (PMF)                 ij seg i corr Vol B ij r r f T k r F  _ ln   r ij seg  Number density of atom pairs of type ij at atom pair distance r Number density of atom pairs of type ij in reference sphere with radius R ij bulk  ij bulk 

Consensus Scoring Cscore: Integrate multiple scoring functions to produce a consensus score that is more accurate than any single function for predicting binding affinity.

Virtual screening by Docking Find weak binders in pool of non -binders Many false positives (96-100%) Consensus Scoring reduces rate of false positives

Concluding remarks Although the reliability of docking methods is not so high, they can provide new suggestions for protein-ligand interactions that otherwise may be overlooked Scoring functions are the Achilles’ heel of docking programs . False positives rates can be reduced using several scoring functions in a consensus-scoring strategy

Docking programs DOCK FlexX GOLD AutoDOCK Hammerhead FLOG

FLEXX Receptor is treated as rigid Incremental construction algorithm: Break Ligand up into rigid fragments Dock fragments into pocket of receptor Reassemble ligand from fragments in low energy conformations

How DOCK works Generate molecular surface of protein C avities in the receptor are used to define spheres (blue) ; the centr e s are potential locations for ligand atoms. thioketal in the HIV1-protease active site Sphere centr e s are m atched to ligand atoms, to determine possible orientations for the ligand. 10 4 orientations generated

GOLD (Genetic Optimisation for Ligand Docking) Performs automated docking with full acyclic ligand flexibility, partial cyclic ligand flexibility and partial protein flexibility in and around active site. Scoring: includes H-bonding term, pairwise dispersion potential (hydrophobic interactions), m olecular and mechanics term for internal energy. Analysis shows algorithm more likely to fail if ligand is large or highly flexible, and more likely to succeed if ligand is polar The GA is encoded to search for H-bonding networks first; Fitness function contains a term for dispersive interactions but takes no account of desolvation, thus underestimates The Hydrophobic Effect

Dock Sem

More Related Content

What's hot (19)

Viewers also liked (20)

Similar to Dock Sem (20)

Recently uploaded (20)

Dock Sem

Editor's Notes