Xuan-Yu Meng

Molecular Docking: A powerful approach for structure-based drug discovery

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.

Phosphoinositide Control of Membrane Protein Function: A Frontier Led by Studies on Ion Channels

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 [Ca(2+)]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.

The Molecular Mechanism by which PIP2 Opens the Intracellular G-Loop Gate of a Kir3.1 Channel

Inwardly rectifying potassium (Kir) channels are characterized by a long pore comprised of continuous transmembrane and cytosolic portions. A high-resolution structure of a Kir3.1 chimera revealed the presence of the cytosolic (G-loop) gate captured in the closed or open conformations. Here, we conducted molecular-dynamics simulations of these two channel states in the presence and absence of phosphatidylinositol bisphosphate (PIP(2)), a phospholipid that is known to gate Kir channels. Simulations of the closed state with PIP(2) revealed an intermediate state between the closed and open conformations involving direct transient interactions with PIP(2), as well as a network of transitional inter- and intrasubunit interactions. Key elements in the G-loop gating transition involved a PIP(2)-driven movement of the N-terminus and C-linker that removed constraining intermolecular interactions and led to CD-loop stabilization of the G-loop gate in the open state. To our knowledge, this is the first dynamic molecular view of PIP(2)-induced channel gating that is consistent with existing experimental data.

Molecular Mechanism of Species-Dependent Sweet Taste toward Artificial Sweeteners

The heterodimer of Tas1R2 and Tas1R3 is a broadly acting sweet taste receptor, which mediates mammalian sweet taste toward natural and artificial sweeteners and sweet-tasting proteins. Perception of sweet taste is a species-selective physiological process. For instance, artificial sweeteners aspartame and neotame taste sweet to humans, apes, and Old World monkeys but not to New World monkeys and rodents. Although specific regions determining the activation of the receptors by these sweeteners have been identified, the molecular mechanism of species-dependent sweet taste remains elusive. Using human/squirrel monkey chimeras, mutagenesis, and molecular modeling, we reveal that the different responses of mammalian species toward the artificial sweeteners aspartame and neotame are determined by the steric effect of a combination of a few residues in the ligand binding pocket. Residues S40 and D142 in the human Tas1R2, which correspond to residues T40 and E142 in the squirrel monkey Tas1R2, were found to be the critical residues for the species-dependent difference in sweet taste. In addition, human Tas1R2 residue I67, which corresponds to S67 in squirrel monkey receptor, modulates the higher affinity of neotame than of aspartame. Our studies not only shed light on the molecular mechanism of species-dependent sweet taste toward artificial sweeteners, but also provide guidance for designing novel effective artificial sweet compounds.

Selective Phosphorylation Modulates the PIP2 Sensitivity of the CaM-SK Channel Complex

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 a key cofactor for activation of small conductance Ca2+-activated potassium channels (SKs) by Ca(2+)-bound calmodulin (CaM). Removal of the endogenous PIP2 inhibits SKs. 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.

Building KCNQ1/KCNE1 Channel Models and Probing their Interactions by Molecular-Dynamics Simulations

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.

Structural Determinants of Phosphatidylinositol 4,5-Bisphosphate (PIP2) Regulation of BK Channel Activity through the RCK1 Ca2+ Coordination Site

Background: PIP2 has been reported to enhance Ca2+-driven gating, but the molecular determinants of this interplay are not known. Results: PIP2 interacts with specific basic residues and enhances Ca2+ gating through the αA-KDRDD-αB structural elements. Conclusion: The RCK1 Ca2+-binding site is coupled to PIP2. Significance: PIP2 is a key element in the regulation of BK channel activity. Big or high conductance potassium (BK) channels are activated by voltage and intracellular calcium (Ca2+). Phosphatidylinositol 4,5-bisphosphate (PIP2), a ubiquitous modulator of ion channel activity, has been reported to enhance Ca2+-driven gating of BK channels, but a molecular understanding of this interplay or even of the PIP2 regulation of this channels activity remains elusive. Here, we identify structural determinants in the KDRDD loop (which follows the αA helix in the RCK1 domain) to be responsible for the coupling between Ca2+ and PIP2 in regulating BK channel activity. In the absence of Ca2+, RCK1 structural elements limit channel activation through a decrease in the channels PIP2 apparent affinity. This inhibitory influence of BK channel activation can be relieved by mutation of residues that (a) connect either the RCK1 Ca2+ coordination site (Asp367 or its flanking basic residues in the KDRDD loop) to the PIP2-interacting residues (Lys392 and Arg393) found in the αB helix or (b) are involved in hydrophobic interactions between the αA and αB helix of the RCK1 domain. In the presence of Ca2+, the RCK1-inhibitory influence of channel-PIP2 interactions and channel activity is relieved by Ca2+ engaging Asp367. Our results demonstrate that, along with Ca2+ and voltage, PIP2 is a third factor critical to the integral control of BK channel activity.

The Cytosolic GH Loop Regulates the Phosphatidylinositol 4,5-Bisphosphate-induced Gating Kinetics of Kir2 Channels

Background: The detailed mechanism of PIP2-induced Kir channel gating remains elusive. Results: Specific mutations increase the flexibility of the cytosolic GH loop and accelerate the PIP2-induced gating kinetics of Kir2 channels. Conclusion: Interactions of the GH loop with the N terminus regulate the PIP2-induced gating kinetics of Kir2 channels. Significance: We identify a novel region in Kir channels involved in the control of PIP2-induced gating. Inwardly rectifying K+ (Kir) channels set the resting membrane potential and regulate cellular excitability. The activity of Kir channels depends critically on the phospholipid PIP2. The molecular mechanism by which PIP2 regulates Kir channel gating is poorly understood. Here, we utilized a combination of computational and electrophysiological approaches to discern structural elements involved in regulating the PIP2-induced gating kinetics of Kir2 channels. We identify a novel role for the cytosolic GH loop. Mutations that directly or indirectly affect GH loop flexibility (e.g. V223L, E272G, D292G) increase both the on- and especially the off-gating kinetics. These effects are consistent with a model in which competing interactions between the CD and GH loops for the N terminus regulate the gating of the intracellular G loop gate.

Computational Approaches for Modeling GPCR Dimerization

Growing experimental evidences suggest that dimerization and oligomerization are important for G Protein- Coupled Receptors (GPCRs) function. The detailed structural information of dimeric/oligomeric GPCRs would be very important to understand their function. Although it is encouraging that recently several experimental GPCR structures in oligomeric forms have appeared, experimental determination of GPCR structures in oligomeric forms is still a big challenge, especially in mimicking the membrane environment. Therefore, development of computational approaches to predict dimerization of GPCRs will be highly valuable. In this review, we summarize computational approaches that have been developed and used for modeling of GPCR dimerization. In addition, we introduce a novel two-dimensional Brownian Dynamics based protein docking approach, which we have recently adapted, for GPCR dimer prediction.

Unifying Mechanism of Controlling Kir3 Channel Activity by G Proteins and Phosphoinositides

The question that started with the pioneering work of Otto Loewi in the 1920s, to identify how stimulation of the vagus nerve decreased heart rate, is approaching its 100th year anniversary. In the meantime, we have learned that the neurotransmitter acetylcholine acting through muscarinic M2 receptors activates cardiac potassium (Kir3) channels via the βγ subunits of G proteins, an important effect that contributes to slowing atrial pacemaker activity. Concurrent stimulation of M1 or M3 receptors hydrolyzes PIP2, a signaling phospholipid essential to maintaining Kir3 channel activity, thus causing desensitization of channel activity and protecting the heart from overinhibition of pacemaker activity. Four mammalian members of the Kir3 subfamily, expressed in heart, brain, endocrine organs, etc., are modulated by a plethora of stimuli to regulate cellular excitability. With the recent great advances in ion channel structural biology, three-dimensional structures of Kir3 channels with PIP2 and the Gβγ subunits are now available. Mechanistic insights have emerged that explain how modulatory control of activity feeds into a core mechanism of channel-PIP2 interactions to regulate the conformation of channel gates. This complex but beautiful system continues to surprise us for almost 100 years with an apparent wisdom in its intricate design.