Showing papers on "Cooperative binding published in 2011"

How does a drug molecule find its target binding site

Abstract: Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, detailed experimental characterization of the process by which such binding occurs has proven challenging. We conducted relatively long, unguided molecular dynamics simulations in which a ligand (the cancer drug dasatinib or the kinase inhibitor PP1) was initially placed at a random location within a box that also contained a protein (Src kinase) to which that ligand was known to bind. In several of these simulations, the ligand correctly identified its target binding site, forming a complex virtually identical to the crystallographically determined bound structure. The simulated trajectories provide a continuous, atomic-level view of the entire binding process, revealing persistent and noteworthy intermediate conformations and shedding light on the role of water molecules. The technique we employed, which does not assume any prior knowledge of the binding site's location, may prove particularly useful in the development of allosteric inhibitors that target previously undiscovered binding sites.

Pathway and mechanism of drug binding to G-protein-coupled receptors

Abstract: How drugs bind to their receptors--from initial association, through drug entry into the binding pocket, to adoption of the final bound conformation, or "pose"--has remained unknown, even for G-protein-coupled receptor modulators, which constitute one-third of all marketed drugs. We captured this pharmaceutically critical process in atomic detail using the first unbiased molecular dynamics simulations in which drug molecules spontaneously associate with G-protein-coupled receptors to achieve final poses matching those determined crystallographically. We found that several beta blockers and a beta agonist all traverse the same well-defined, dominant pathway as they bind to the β(1)- and β(2)-adrenergic receptors, initially making contact with a vestibule on each receptor's extracellular surface. Surprisingly, association with this vestibule, at a distance of 15 A from the binding pocket, often presents the largest energetic barrier to binding, despite the fact that subsequent entry into the binding pocket requires the receptor to deform and the drug to squeeze through a narrow passage. The early barrier appears to reflect the substantial dehydration that takes place as the drug associates with the vestibule. Our atomic-level description of the binding process suggests opportunities for allosteric modulation and provides a structural foundation for future optimization of drug-receptor binding and unbinding rates.

Binding leverage as a molecular basis for allosteric regulation.

Abstract: Allosteric regulation involves conformational transitions or fluctuations between a few closely related states, caused by the binding of effector molecules. We introduce a quantity called binding leverage that measures the ability of a binding site to couple to the intrinsic motions of a protein. We use Monte Carlo simulations to generate potential binding sites and either normal modes or pairs of crystal structures to describe relevant motions. We analyze single catalytic domains and multimeric allosteric enzymes with complex regulation. For the majority of the analyzed proteins, we find that both catalytic and allosteric sites have high binding leverage. Furthermore, our analysis of the catabolite activator protein, which is allosteric without conformational change, shows that its regulation involves other types of motion than those modulated at sites with high binding leverage. Our results point to the importance of incorporating dynamic information when predicting functional sites. Because it is possible to calculate binding leverage from a single crystal structure it can be used for characterizing proteins of unknown function and predicting latent allosteric sites in any protein, with implications for drug design.

Conversion of Molecular Information by Luminescent Nanointerface Self-Assembled from Amphiphilic Tb(III) Complexes

Abstract: A novel amphiphilic Tb3+ complex (TbL+) having anionic bis(pyridine) arms and a hydrophobic alkyl chain is developed. It spontaneously self-assembles in water and gives stable vesicles that show sensitized luminescence of Tb3+ ions at neutral pH. This TbL+ complex is designed to show coordinative unsaturation, i.e., water molecules occupy some of the first coordination spheres and are replaceable upon binding of phosphate ions. These features render TbL+ self-assembling receptor molecules which show increase in the luminescence intensity upon binding of nucleotides. Upon addition of adenosine triphosphate (ATP), significant amplification of luminescent intensity was observed. On the other hand, ADP showed moderately increased luminescence and almost no enhancement was observed for AMP. Very interestingly, the increase in luminescence intensity observed for ATP and ADP showed sigmoidal dependence on the concentration of added nucleotides. It indicates positive cooperative binding of these nucleotides to Tb...

Allosteric control of ligand-binding affinity using engineered conformation-specific effector proteins

Abstract: We describe a phage display methodology for engineering synthetic antigen binders (sABs) that recognize either the apo or the ligand-bound conformation of maltose-binding protein (MBP). sABs that preferentially recognize the maltose-bound form of MBP act as positive allosteric effectors by substantially increasing the affinity for maltose. A crystal structure of a sAB bound to the closed form of MBP reveals the basis for this allosteric effect. We show that sABs that recognize the bound form of MBP can rescue the function of a binding-deficient mutant by restoring its natural affinity for maltose. Furthermore, the sABs can enhance maltose binding in vivo, as they provide a growth advantage to bacteria under low-maltose conditions. The results demonstrate that structure-specific sABs can be engineered to dynamically control ligand-binding affinities by modulating the transition between different conformations.

Modulation of the stability of the Salmonella fourU-type RNA thermometer

Abstract: RNA thermometers are translational control elements that regulate the expression of bacterial heat shock and virulence genes. They fold into complex secondary structures that block translation at low temperatures. A temperature increase releases the ribosome binding site and thus permits translation initiation. In fourU-type RNA thermometers, the AGGA sequence of the SD region is paired with four consecutive uridines. We investigated the melting points of the wild-type and mutant sequences. It was decreased by 5°C when a stabilizing GC basepair was exchanged by an AU pair or increased by 11°C when an internal AG mismatch was converted to a GC pair, respectively. Stabilized or destabilized RNA structures are directly correlated with decreased or increased in vivo gene expression, respectively. Mg2+ also affected the melting point of the fourU thermometer. Variations of the Mg2+ concentration in the physiological range between 1 and 2 mM translated into a 2.8°C shift of the melting point. Thus, Mg2+ binding to the hairpin RNA is regulatory relevant. Applying three different NMR techniques, two Mg2+ binding sites were found in the hairpin structure. One of these binding sites could be identified as outer sphere binding site that is located within the fourU motif. Binding of the two Mg2+ ions exhibits a positive cooperativity with a Hill coefficient of 1.47. Free energy values delta G for Mg2+ binding determined by NMR are in agreement with data determined from CD measurements. Physiological Mg2+ concentrations reduce enthalpy and entropy values of uncorrelated base pair opening processes for almost all nucleobases.

Sequence-dependent cooperative binding of p53 to DNA targets and its relationship to the structural properties of the DNA targets

Abstract: The prime mechanism by which p53 acts as a tumor suppressor is as a transcription factor regulating the expression of diverse downstream genes. The DNA-binding domain of p53 (p53DBD) interacts with defined DNA sites and is the main target for mutations in human primary tumors. Here, we show that the CWWG motif, found in the center of each consensus p53 half-site, is a key player in p53/DNA interactions. Gel-mobility-shift assays provide a unique opportunity to directly observe the various oligomeric complexes formed between p53DBD and its target sites. We demonstrate that p53DBD binds to p53 consensus sites containing CATG with relatively low cooperativity, as both dimers and tetramers, and with even lower cooperativity to such sites containing spacer sequences. p53DBD binds to sites containing CAAG and CTAG with measurable affinity only when imbedded in two contiguous p53 half-sites and only as tetramers (with very high cooperativity). There are three orders-of-magnitude difference in the cooperativity of interaction between sites differing in their non-contacted step, and further two orders-of-magnitude difference as a function of spacer sequences. By experimentally measuring the global structural properties of these sites, by cyclization kinetics of DNA minicircles, we correlate these differences with the torsional flexibility of the binding sites.

Characterization of calmodulin binding domains in TRPV2 and TRPV5 C-tails.

Abstract: The transient receptor potential channels TRPV2 and TRPV5 belong to the vanilloid TRP subfamily. TRPV2 is highly similar to TRPV1 and shares many common properties with it. TRPV5 (and also its homolog TRPV6) is a rather distinct member of the TRPV subfamily. It is distant for being strictly Ca(2+)-selective and features quite different properties from the rest of the TRPV subfamily. It is known that TRP channels are regulated by calmodulin in a calcium-dependent manner. In our study we identified a calmodulin binding site on the C-termini of TRPV2 (654-683) and TRPV5 (587-616) corresponding to the consensus CaM binding motif 1-5-10. The R679 and K681 single mutants of TRPV2 caused a 50% decrease in binding affinity and a double mutation of K661/K664 of the same peptide lowered the binding affinity by up to 75%. A double mutation of R606/K607 and triple mutation of R594/R606/R610 in TRPV5 C-terminal peptide resulted in the total loss of binding affinity to calmodulin. These results demonstrate that the TRPV2 C-tail and TRPV5 C-tail contain calmodulin binding sites and that the basic residues are strongly involved in TRP channel binding to calmodulin.

Thermodynamic properties of the specific binding between Ag+ ions and C:C mismatched base pairs in duplex DNA.

Abstract: Metal-mediated base pairs formed by the interaction between metal ions and artificial bases in oligonucleotides have been developed for potential applications in nanotechnology. We recently found that a natural C:C mismatched base pair bound to an Ag(+) ion to generate a novel metal-mediated base pair in duplex DNA. Preparation of the novel C-Ag-C base pair involving natural bases is more convenient than that of metal-mediated base pairs involving artificial bases because time-consuming base synthesis is not required. Here, we examined the thermodynamic properties of the binding between the Ag(+) ion and each of single and double C:C mismatched base pair in duplex DNA by isothermal titration calorimetry. The Ag(+) ion specifically bound to the C:C mismatched base pair at a 1:1 molar ratio with 10(6) M(-1) binding constant, which was significantly larger than those for nonspecific metal ion-DNA interactions. The specific binding between the Ag(+) ion and the single C:C mismatched base pair was mainly driven by the positive dehydration entropy change and the negative binding enthalpy change. In the interaction between the Ag(+) ion and each of the consecutive and interrupted double C:C mismatched base pairs, stoichiometric binding at a 1:1 molar ratio was achieved in each step of the first and second Ag(+) binding. The binding affinity for the second Ag(+) binding was similar to that for the first Ag(+) binding. Stoichiometric binding without interference and negative cooperativity may be favorable for aligning multiple Ag(+) ions in duplex DNA for applications of the metal-mediated base pairs in nanotechnology.

Steric and Allosteric Effects of Fatty Acids on the Binding of Warfarin to Human Serum Albumin Revealed by Molecular Dynamics and Free Energy Calculations

Abstract: Human serum albumin (HSA) binds with drugs and fatty acids (FAs). This study was initiated to elucidate the relationship between the warfarin binding affinity of HSA and the positions of bound FA molecules. Molecular dynamics simulations of 11 HSA-warfarin-myristate complexes were performed. HSA-warfarin binding free energy was then calculated for each of the complexes by the molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) method. The results indicated that the magnitude of the binding free energy was smaller in HSA-warfarin complexes that had 4 or more myristate molecules than in complexes with no myristate molecules. The unfavorable effect on the HSA-warfarin binding affinity was caused sterically by the binding of a myristate molecule to the FA binding site closest to the warfarin binding site. On the other hand, the magnitude of HSA-warfarin binding free energy was largest when 3 myristate molecules were bound to the high-affinity sites. The strongest HSA-warfarin binding was attributable to favorable entropic contribution related to larger atomic fluctuations of the amino acid residues at the warfarin binding site. In the binding of 2 myristate molecules to the sites with the highest and second-highest affinities, allosteric modulation that enhanced electrostatic interactions between warfarin and some of the amino acid residues around the warfarin binding site was observed. This study clarified the structural and energetic properties of steric/allosteric effects of FAs on the HSA-warfarin binding affinity and illustrated the approach to analyze protein–ligand interactions in situations such that multiple ligands bind to the other sites of the protein.

Solution structure of the Zβ domain of human DNA-dependent activator of IFN-regulatory factors and its binding modes to B- and Z-DNAs

Abstract: The DNA-dependent activator of IFN-regulatory factors (DAI), also known as DLM-1/ZBP1, initiates an innate immune response by binding to foreign DNAs in the cytosol. For full activation of the immune response, three DNA binding domains at the N terminus are required: two Z-DNA binding domains (ZBDs), Zα and Zβ, and an adjacent putative B-DNA binding domain. The crystal structure of the Zβ domain of human DAI (hZβDAI) in complex with Z-DNA revealed structural features distinct from other known Z-DNA binding proteins, and it was classified as a group II ZBD. To gain structural insights into the DNA binding mechanism of hZβDAI, the solution structure of the free hZβDAI was solved, and its bindings to B- and Z-DNAs were analyzed by NMR spectroscopy. Compared to the Z-DNA–bound structure, the conformation of free hZβDAI has notable alterations in the α3 recognition helix, the “wing,” and Y145, which are critical in Z-DNA recognition. Unlike some other Zα domains, hZβDAI appears to have conformational flexibility, and structural adaptation is required for Z-DNA binding. Chemical-shift perturbation experiments revealed that hZβDAI also binds weakly to B-DNA via a different binding mode. The C-terminal domain of DAI is reported to undergo a conformational change on B-DNA binding; thus, it is possible that these changes are correlated. During the innate immune response, hZβDAI is likely to play an active role in binding to DNAs in both B and Z conformations in the recognition of foreign DNAs.

A negative cooperativity mechanism of human CYP2E1 inferred from molecular dynamics simulations and free energy calculations.

Abstract: Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. Thi...

The Role of Transmembrane Domain 3 in the Actions of Orthosteric, Allosteric, and Atypical Agonists of the M4 Muscarinic Acetylcholine Receptor

Abstract: Despite the discovery of a diverse range of novel agonists and allosteric modulators of the M4 muscarinic acetylcholine (ACh) receptor (mAChR), little is known about how such ligands activate the receptor. We used site-directed mutagenesis of conserved residues in transmembrane 3 (TMIII), a key region involved in G protein-coupled receptor activation, to probe the binding and function of prototypical orthosteric mAChR agonists, allosteric modulators, and “atypical” agonists. We found that most mutations did not affect the binding of the allosteric modulators, with the exception of W108 3.28 A and L109 3.29 A (which may contribute directly to the interface between allosteric and orthosteric sites) and mutation D112 3.32 N (which may cause a global disruption of a hydrogen bond network). Although numerous mutations affected signaling, we did not identify amino acids that were important for the functional activity of any one class of agonist (orthosteric, allosteric, or atypical) to the exclusion of any others, suggesting that TMIII is key for the transmission of stimulus irrespective of the agonist. We also identified two key residues, Trp108 3.28 and Asp112 3.32 , that are essential for the transmission of binding cooperativity between 3-amino-5-chloro-6-methoxy-4-methyl-thieno[2,3-b]pyridine2-carboxylic acid cyclopropylamide (LY2033298) and ACh. Finally, we found that LY2033298 was able to rescue functionally impaired signaling of ACh at the majority of mutants tested in a manner that was inversely correlated with the ACh signaling efficacy, indicating that a key part of the mechanism of the positive cooperativity mediated by LY2033298 on the endogenous agonist involves a global drive of the receptor toward an active conformation.

Strong binding affinity of a zinc-porphyrin-based receptor for halides through the cooperative effects of quadruple C-H hydrogen bonds and axial ligation.

Abstract: A new type of host compound (1), tetraphenyl zinc-porphyrin (ZnTPP) that contains four triazole groups at the ortho-position of each phenyl group, has been synthesized and characterized by using (1)H, (13)C NMR, and MALDI-TOF-MS analyses. The host-guest complex formation between 1 and halides was investigated by using (1)H NMR spectroscopy in [D(6)]DMSO. The triazole, benzyl, and phenylene proton signals were shifted upfield by the addition of halides in the form of tetrabutylammonium salts, which implies that the triazole protons in 1 are allocated very closely to the porphyrin ring and are directed toward the binding pocket over the porphyrin ring during the formation of hydrogen bonds. The UV/Vis absorption spectra showed that both the Soret and Q band absorptions of 1 underwent a strong redshift due to the addition of halides. Compound 1 exhibited surprisingly strong binding affinities for the halides, where the association constants for Cl(-), Br(-), and I(-) binding were estimated to be larger than 10(8), 1.79×10(7), and 1.84×10(5) M(-1), respectively. The UV/Vis absorption changes and the result of competitive titration using 4-tert-butylpyridine indicated that the cooperative effects of axial coordination and C-H···X hydrogen bond interactions resulted in the strong binding affinity of 1 to halides.

A β-to-β 2,5-thienylene-bridged cyclic porphyrin tetramer: its rational synthesis and 1 : 2 binding mode with C60

Abstract: A β-to-β2,5-thienylene-bridged cyclic porphyrin tetramer was rationally synthesized via a concise synthetic route. The tetraporphyrin exhibits a positive cooperative binding ability to C60 and demonstrates a new potential of the nonplanar, distorted cyclic porphyrin arrays.

The SRSF1 linker induces semi-conservative ESE binding by cooperating with the RRMs

Abstract: SR proteins promote spliceosome formation by recognizing exonic splicing enhancers (ESEs) during pre-mRNA splicing. Each SR protein binds diverse ESEs using strategies that are yet to be elucidated. Here, we show that the RNA-binding domain (RBD) of SRSF1 optimally binds to decameric purine rich ESE sequences although locations of purines are not stringently specified. The presence of uracils either within or outside of the recognition site is detrimental for binding with SRSF1. The entire RBD, comprised of two RRMs and a glycine-rich linker, is essential for ESE binding. Mutation within each segment reduced or nearly abolished binding, suggesting that these segments mediate cooperative binding. The linker plays a decisive role in organizing ESE binding. The flanking basic regions of the linker appear to communicate with each other in bringing the two RRMs close together to form the complex with RNA. Our study thus suggests semi-conservative adaptable interaction between ESE and SRSF1, and such binding mode is not only essential for the recognition of plethora of physiological ESE sequences but may also be essential for the interaction with various factors during the spliceosome assembly.

Combinatorial Regulation by MrpC2 and FruA Involves Three Sites in the fmgE Promoter Region during Myxococcus xanthus Development

Abstract: Starvation causes cells in a dense population of Myxococcus xanthus to change their gliding movements and construct mounds. Short-range C-signaling between rod-shaped cells within mounds induces gene expression that promotes differentiation into spherical spores. Several C-signal-dependent genes have been shown to be regulated by cooperative binding of two transcription factors to the promoter region. These FruA- and MrpC2-regulated genes (designated fmg) each exhibit a different arrangement of binding sites. Here, we describe fmgE, which appears to be regulated by three sites for cooperative binding of FruA and MrpC2. Chromatin immunoprecipitation analysis showed that association of MrpC2 and/or its longer form, MrpC with the fmgE promoter region, depends on FruA, consistent with cooperative binding of the two proteins in vivo. Electrophoretic mobility shift assays with purified His10-MrpC2 and FruA-His6 indicated cooperative binding in vitro to three sites in the fmgE promoter region. The effects of mutations on binding in vitro and on expression of fmgE-lacZ fusions correlated site 1 (at about position −100 relative to the transcriptional start site) with negative regulation and site 2 (just upstream of the promoter) and site 3 (at about position +100) with positive regulation. Site 3 was bound by His10-MrpC2 alone, or the combination of His10-MrpC2 and FruA-His6, with the highest affinity, followed by site 1 and then site 2, supporting a model in which site 3 recruits MrpC2 and FruA to the fmgE promoter region, site 1 competes with site 2 for transcription factor binding, and site 2 occupancy is required to activate the promoter but only occurs when C-signaling produces a high concentration of active FruA.

Identification of a Novel Allosteric Binding Site in the CXCR2 Chemokine Receptor

Abstract: We have shown previously that different chemical classes of small-molecule antagonists of the human chemokine CXCR2 receptor interact with distinct binding sites of the receptor. Although an intracellular binding site for diarylurea CXCR2 antagonists, such as N-(2-bromophenyl)-N′-(7-cyano-1H-benzotriazol-4-yl)urea (SB265610), and thiazolopyrimidine compounds was recently mapped by mutagenesis studies, we now report on an imidazolylpyrimidine antagonist binding pocket in the transmembrane domain of CXCR2. Using different CXCR2 orthologs, chimeric proteins, site-directed mutagenesis, and in silico modeling, we have elucidated the binding mode of this antagonist. Our in silico-guided mutagenesis studies indicate that the ligand binding cavity for imidazolylpyrimidine compounds in CXCR2 is located between transmembrane (TM) helices 3 (Phe1303.36), 5 (Ser2175.44, Phe2205.47), and 6 (Asn2686.52, Leu2716.55) and suggest that these antagonists enter CXCR2 via the TM5-TM6 interface. It is noteworthy that the same interface is postulated as the ligand entry channel in the opsin receptor and is occupied by lipid molecules in the recently solved crystal structure of the CXCR4 chemokine receptor, suggesting a general ligand entrance mechanism for nonpolar ligands to G protein-coupled receptors. The identification of a novel allosteric binding cavity in the TM domain of CXCR2, in addition to the previously identified intracellular binding site, shows the diversity in ligand recognition mechanisms by this receptor and offers new opportunities for the structure-based design of small allosteric modulators of CXCR2 in the future.

Negatively cooperative binding of high-density lipoprotein to the HDL receptor SR-BI.

Abstract: Scavenger receptor class B, type I (SR-BI), is a high-density lipoprotein (HDL) receptor, which also binds low-density lipoprotein (LDL), and mediates the cellular selective uptake of cholesteryl esters from lipoproteins. SR-BI also is a coreceptor for hepatitis C virus and a signaling receptor that regulates cell metabolism. Many investigators have reported that lipoproteins bind to SR-BI via a single class of independent (not interacting), high-affinity binding sites (one site model). We have reinvestigated the ligand concentration dependence of 125I-HDL binding to SR-BI and SR-BI-mediated specific uptake of [3H]CE from [3H]CE-HDL using an expanded range of ligand concentrations (<1 μg of protein/mL, lower than previously reported). Scatchard and nonlinear least-squares model fitting analyses of the binding and uptake data were both inconsistent with a single class of independent binding sites binding univalent lipoprotein ligands. The data are best fit by models in which SR-BI has either two independen...

Glycine hinges with opposing actions at the acetylcholine receptor-channel transmitter binding site.

Abstract: The extent to which agonists activate synaptic receptor-channels depends on both the intrinsic tendency of the unliganded receptor to open and the amount of agonist binding energy realized in the channel-opening process. We examined mutations of the nicotinic acetylcholine receptor transmitter binding site (α subunit loop B) with regard to both of these parameters. αGly147 is an “activation” hinge where backbone flexibility maintains high values for intrinsic gating, the affinity of the resting conformation for agonists and net ligand binding energy. αGly153 is a “deactivation” hinge that maintains low values for these parameters. αTrp149 (between these two glycines) serves mainly to provide ligand binding energy for gating. We propose that a concerted motion of the two glycine hinges (plus other structural elements at the binding site) positions αTrp149 so that it provides physiologically optimal binding and gating function at the nerve-muscle synapse.

Molecular Dynamics and Docking Studies on Cardiac Troponin C

Abstract: Cardiac troponin C (cTnC) is the Ca²⁺ dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 A, 5.5 A, 6 A, 6.5 A, 7.0 A and 7.5 A were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.