Molecular docking has become an increasingly important tool for drug discovery. In this review, we present a brief introduction of the available molecular docking methods, and their development and applications in drug discovery. The relevant basic theories, including sampling algorithms and scoring functions, are summarized. The differences in and performance of available docking software are also discussed. Flexible receptor molecular docking approaches, especially those including backbone flexibility in receptors, are a challenge for available docking methods. A recently developed Local Move Monte Carlo (LMMC) based approach is introduced as a potential solution to flexible receptor docking problems. Three application examples of molecular docking approaches for drug discovery are provided.

Anionic phospholipids are critical constituents of the inner leaflet of the plasma membrane, ensuring appropriate membrane topology of transmembrane proteins. Additionally, in eukaryotes, the negatively charged phosphoinositides serve as key signals not only through their hydrolysis products but also through direct control of transmembrane protein function. Direct phosphoinositide control of the activity of ion channels and transporters has been the most convincing case of the critical importance of phospholipid-protein interactions in the functional control of membrane proteins. Furthermore, second messengers, such as [Ca2+]i, or posttranslational modifications, such as phosphorylation, can directly or allosterically fine-tune phospholipid-protein interactions and modulate activity. Recent advances in structure determination of membrane proteins have allowed investigators to obtain complexes of ion channels with phosphoinositides and to use computational and experimental approaches to probe the dynamic mechanisms by which lipid-protein interactions control active and inactive protein states.

Phosphatidylinositol bisphosphate (PIP2) regulates the activities of many membrane proteins including ion channels through direct interactions. However, the affinity of PIP2 is so high for some channel proteins that its physiological role as a modulator has been questioned. Here we show that PIP2 is an important cofactor for activation of small conductance Ca2+-activated potassium channels (SK) by Ca2+-bound calmodulin (CaM). Removal of the endogenous PIP2 inhibits SK channels. The PIP2-binding site resides at the interface of CaM and the SK C-terminus. We further demonstrate that the affinity of PIP2 for its target proteins can be regulated by cellular signaling. Phosphorylation of CaM T79, located adjacent to the PIP2-binding site, by Casein Kinase 2 reduces the affinity of PIP2 for the CaM-SK channel complex by altering the dynamic interactions among amino acid residues surrounding the PIP2-binding site. This effect of CaM phosphorylation promotes greater channel inhibition by G-protein-mediated hydrolysis of PIP2.

The slow delayed rectifier (I(KS)) channel is composed of KCNQ1 (pore-forming) and KCNE1 (auxiliary) subunits, and functions as a repolarization reserve in the human heart. Design of I(KS)-targeting anti-arrhythmic drugs requires detailed three-dimensional structures of the KCNQ1/KCNE1 complex, a task made possible by Kv channel crystal structures (templates for KCNQ1 homology-modeling) and KCNE1 NMR structures. Our goal was to build KCNQ1/KCNE1 models and extract mechanistic information about their interactions by molecular-dynamics simulations in an explicit lipid/solvent environment. We validated our models by confirming two sets of model-generated predictions that were independent from the spatial restraints used in model-building. Detailed analysis of the molecular-dynamics trajectories revealed previously unrecognized KCNQ1/KCNE1 interactions, whose relevance in I(KS) channel function was confirmed by voltage-clamp experiments. Our models and analyses suggest three mechanisms by which KCNE1 slows KCNQ1 activation: by promoting S6 bending at the Pro hinge that closes the activation gate; by promoting a downward movement of gating charge on S4; and by establishing a network of electrostatic interactions with KCNQ1 on the extracellular surface that stabilizes the channel in a pre-open activated state. Our data also suggest how KCNE1 may affect the KCNQ1 pore conductance.