Modeling approaches for ligand-based 3D similarity.

3D ligand-based similarity approaches are widely used in the early phases of drug discovery for tasks such as hit finding by virtual screening or compound design with quantitative structure-activity relationships. Here in we review widely used software for performing such tasks. Some techniques are based on relatively mature technology, shape-based similarity for instance. Typically, these methods remained in the realm of the expert user, the experienced modeler. However, advances in implementation and speed have improved usability and allow these methods to be applied to databases comprising millions of compounds. There are now many reports of such methods impacting drug-discovery projects. As such, the medicinal chemistry community has become the intended market for some of these new tools, yet they may consider the wide array and choice of approaches somewhat disconcerting. Each method has subtle differences and is better suited to certain tasks than others. In this article we review some of the widely used computational methods via application, provide straightforward background on the underlying theory and provide examples for the interested reader to pursue in more detail. In the new era of preclinical drug discovery there will be ever more pressure to move faster and more efficiently, and computational approaches based on 3D ligand similarity will play an increasing role in in this process.

A comparison of ligand based virtual screening methods and application to corticotropin releasing factor 1 receptor.

Ligand based virtual screening approaches were applied to the CRF1 receptor. We compared ECFP6 fingerprints, FTrees, Topomers, Cresset FieldScreen, ROCS OpenEye shape Tanimoto, OpenEye combo-score and OpenEye electrostatics. The 3D methods OpenEye Shape Tanimoto, combo-score and Topomers performed the best at separating actives from inactives in retrospective experiments. By virtue of their higher enrichment the same methods identified more active scaffolds. However, amongst a given number of active compounds the Cresset and OpenEye electrostatic methods contained more scaffolds and returned ranked compounds with greater diversity. A selection of the methods were employed to recommend compounds for screening in a prospective experiment. New CRF1 actives antagonists were found. The new actives contained different underlying chemical architecture to the query molecules, results indicative of successful scaffold-hopping.

Conformational dynamics of a lipid-interacting protein: MD simulations of saposin B.

Saposin B is a water soluble alpha-helical protein which can bind to membranes and extract selected lipids, especially cerebroside sulfates. The X-ray structure of saposin B is homodimeric. There are two conformations of the dimer in the crystal-one with a closed central cavity (the AB dimer) and one (the CD dimer) with a more open cavity. We have conducted a series of short (5 ns) molecular dynamics simulations of saposin B, starting from both the AB and CD conformations and with/without bound lipid and/or water molecules within the central hydrophobic cavity. The more open (CD) dimer showed greater conformational drift than the AB dimer. The conformational drift was also somewhat higher in the absence of bound lipid. Two more extended (30 ns) simulations of AB and CD dimers were performed and analyzed in terms of changes in intersubunit packing within the dimers. The AB dimer remained largely unchanged in conformation over the duration of the extended simulation. In contrast, the CD dimer underwent a substantial conformational change corresponding to a 'scissor' motion of the two monomers so as to compress the central cavity to a more closed conformation than that seen in the AB dimer structure. A H-bond between the Q53 and Y54 side chains of the alpha3 helices of the two opposing monomers seems to hold the dimer in this 'scissor-closed' conformation. We suggest that a cycle of conformational changes, expanding and compressing the central cavity of the saposin B dimer, may play a key role in facilitating lipid extraction from bilayers.

SGTx1, a Kv channel gating-modifier toxin, binds to the interfacial region of lipid bilayers.

SGTx1 is a gating-modifier toxin that has been shown to inhibit the voltage-gated potassium channel Kv2.1. SGTx1 is thought to bind to the S3b-S4a region of the voltage-sensor, and is believed to alter the energetics of gating. Gating-modifier toxins such as SGTx1 are of interest as they can be used to probe the structure and dynamics of their target channels. Although there are experimental data for SGTx1, its interaction with lipid bilayer membranes remains to be characterized. We performed atomistic and coarse-grained molecular dynamics simulations to study the interaction of SGTx1 with a POPC and a 3:1 POPE/POPG lipid bilayer membrane. We reveal the preferential partitioning of SGTx1 into the water/membrane interface of the bilayer. We also show that electrostatic interactions between the charged residues of SGTx1 and the lipid headgroups play an important role in stabilizing SGTx1 in a bilayer environment.

Vstx1, a modifier of Kv channel gating, localizes to the interfacial region of lipid bilayers.

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.