Delineating Protein-Protein Curvilinear Dissociation Pathways and Energetics with Naïve Multiple-Walker Umbrella Sampling Simulations.

The protein-protein interaction energetics can be obtained by calculating the potential of mean force (PMF) from umbrella sampling (US) simulations, in which samplings are often enhanced along a predefined vector as the reaction coordinate. However, any slight change in the vector may significantly vary the calculated PMF, and therefore the energetics using a random choice of vector may mislead. A non-predefined curve path-based sampling enhancement approach is a natural alternative, but was relatively less explored for protein-protein systems. In this work, dissociation of the barnase-barstar complex is simulated by implementing non-predefined curvilinear pathways in US simulations. A simple variational principle is applied to determine the lower bound PMF, which could be used to derive the standard free energy of binding. Two major dissociation pathways, which include interactions with the RNA-binding loop and the Val 36 to Gly 40 loop, are observed. Further, the proposed approach was used to discriminate the decoys from protein-protein docking studies. © 2019 Wiley Periodicals, Inc.

Recovery of the poisoned topoisomerase II for DNA religation: coordinated motion of the cleavage core revealed with the microsecond atomistic simulation.

Type II topoisomerases resolve topological problems of DNA double helices by passing one duplex through the reversible double-stranded break they generated on another duplex. Despite the wealth of information in the cleaving operation, molecular understanding of the enzymatic DNA ligation remains elusive. Topoisomerase poisons are widely used in anti-cancer and anti-bacterial therapy and have been employed to entrap the intermediates of topoisomerase IIß with religatable DNA substrate. We removed drug molecules from the structure and conducted molecular dynamics simulations to investigate the enzyme-mediated DNA religation. The drug-unbound intermediate displayed transitions toward the resealing-compliant configuration: closing distance between the cleaved DNA termini, B-to-A transformation of the double helix, and restoration of the metal-binding motif. By mapping the contact configurations and the correlated motions between enzyme and DNA, we identified the indispensable role of the linker preceding winged helix domain (WHD) in coordinating the movements of TOPRIM, the nucleotide-binding motifs, and the bound DNA substrate during gate closure. We observed a nearly vectorial transition in the recovery of the enzyme and identified the previously uncharacterized roles of Asn508 and Arg677 in DNA rejoining. Our findings delineate the dynamic mechanism of the DNA religation conducted by type II topoisomerases.

Review structure- and dynamics-based computational design of anticancer drugs.

Cancer is a class of highly complex diseases involving multiple genes and multiple cross-talks between signaling networks. Cancer cells may be developed from inherited defects or acquired damages of DNA. However, many cancers are resistant to treatment, and metastasis of cancers makes the disease even more intractable. Secondary malignancies are frequently observed after cancer chemotherapy. The call for more effective cancer therapy is obligatory. Using drug-cocktails that combine multiple anti-cancer agents working in different mechanisms has been a standard treatment of cancers to overcome the drug resistance problem. More recently, design of multiple ligands (may be more easily understood as "multiple target ligands"), i.e., single agents that target multiple biomolecules in a rational manner, receives increasing attention. For those who work on computational drug design, such tasks serve as new opportunities for achieving drugs with more effective pharmacological actions, in addition to designing compounds with better binding affinity, better selectivity, or to discovering compounds that can exert their actions allosterically. Some recent methodological developments on computational drug design are reviewed, and a few recent drug design efforts on a selected set of targets (topoisomerases, Ras proteins, protein kinases, and histone deacetylases) toward cancer treatment and cancer prevention are summarized.

Tetrahydropyran- and tetrahydrofuran-containing diarylheptanoids from Hedychium coronarium rhizomes.

Seven new diarylheptanoids (1-7) were isolated from the n-BuOH-soluble fraction of the rhizome of Hedychium coronarium. Hedycoropyrans A-C (1-3) contain a tetrahydropyran moiety, while hedycorofurans A-D (4-7) contain a tetrahydrofuran moiety, belonging to a rare structural class of diarylheptanoids. Their structures including stereochemistry were elucidated on the basis of 1D and 2D NMR and ECD spectroscopic analyses and HRESIMS data of the parent compounds and the isopropylidene derivatives of 4 and 7.

idTarget: a web server for identifying protein targets of small chemical molecules with robust scoring functions and a divide-and-conquer docking approach.

Identification of possible protein targets of small chemical molecules is an important step for unravelling their underlying causes of actions at the molecular level. To this end, we construct a web server, idTarget, which can predict possible binding targets of a small chemical molecule via a divide-and-conquer docking approach, in combination with our recently developed scoring functions based on robust regression analysis and quantum chemical charge models. Affinity profiles of the protein targets are used to provide the confidence levels of prediction. The divide-and-conquer docking approach uses adaptively constructed small overlapping grids to constrain the searching space, thereby achieving better docking efficiency. Unlike previous approaches that screen against a specific class of targets or a limited number of targets, idTarget screen against nearly all protein structures deposited in the Protein Data Bank (PDB). We show that idTarget is able to reproduce known off-targets of drugs or drug-like compounds, and the suggested new targets could be prioritized for further investigation. idTarget is freely available as a web-based server at http://idtarget.rcas.sinica.edu.tw.

Drug-induced conformational population shifts in topoisomerase-DNA ternary complexes.

Type II topoisomerases (TOP2) are enzymes that resolve the topological problems during DNA replication and transcription by transiently cleaving both strands and forming a cleavage complex with the DNA. Several prominent anti-cancer agents inhibit TOP2 by stabilizing the cleavage complex and engendering permanent DNA breakage. To discriminate drug binding modes in TOP2-α and TOP2-ß, we applied our newly developed scoring function, dubbed AutoDock4RAP, to evaluate the binding modes of VP-16, m-AMSA, and mitoxantrone to the cleavage complexes. Docking reproduced crystallographic binding mode of VP-16 in a ternary complex of TOP2-ß with root-mean-square deviation of 0.65 Å. Molecular dynamics simulation of the complex confirmed the crystallographic binding mode of VP-16 and the conformation of the residue R503. Drug-related conformational changes in R503 have been observed in ternary complexes with m-AMSA and mitoxantrone. However, the R503 rotamers in these two simulations deviate from their crystallographic conformations, indicating a relaxation dynamics from the conformations determined with the drug replacement procedure. The binding mode of VP-16 in the cleavage complex of TOP2-α was determined by the conjoint use of docking and molecular dynamics simulations, which fell within a similar binding pocket of TOP2-ß cleavage complex. Our findings may facilitate more efficient design efforts targeting TOP2-α specific drugs.

AutoBind: automatic extraction of protein-ligand-binding affinity data from biological literature.

MOTIVATION: Determination of the binding affinity of a protein-ligand complex is important to quantitatively specify whether a particular small molecule will bind to the target protein. Besides, collection of comprehensive datasets for protein-ligand complexes and their corresponding binding affinities is crucial in developing accurate scoring functions for the prediction of the binding affinities of previously unknown protein-ligand complexes. In the past decades, several databases of protein-ligand-binding affinities have been created via visual extraction from literature. However, such approaches are time-consuming and most of these databases are updated only a few times per year. Hence, there is an immediate demand for an automatic extraction method with high precision for binding affinity collection. RESULT: We have created a new database of protein-ligand-binding affinity data, AutoBind, based on automatic information retrieval. We first compiled a collection of 1586 articles where the binding affinities have been marked manually. Based on this annotated collection, we designed four sentence patterns that are used to scan full-text articles as well as a scoring function to rank the sentences that match our patterns. The proposed sentence patterns can effectively identify the binding affinities in full-text articles. Our assessment shows that AutoBind achieved 84.22% precision and 79.07% recall on the testing corpus. Currently, 13 616 protein-ligand complexes and the corresponding binding affinities have been deposited in AutoBind from 17 221 articles. AVAILABILITY: AutoBind is automatically updated on a monthly basis, and it is freely available at http://autobind.csie.ncku.edu.tw/ and http://autobind.mc.ntu.edu.tw/. All of the deposited binding affinities have been refined and approved manually before being released.

Biomimetic affinity sensor for the ultrasensitive detection of neonicotinoids.

Multiple pesticides are often used in combination to protect crops from pests. This makes rapid on-site detection of pesticide contamination challenging. Herein, we describe a method for simultaneous detection of diverse neonicotinoid pesticides using a sensor that combines neonicotinoid-specific odorant-binding protein 2 (OBP2), which was cloned from an insect chemical sensing protein and modified gold nanoparticles with local surface plasmon resonance (LSPR)-based digital nanoplasmonometry (DiNM). When neonicotinoid pesticides bind to OBP2 on gold nanoparticles, the induced LSPR shift peak wavelength is too small to be measured using conventional LSPR immunoassays. DiNM records and compares the scattered image intensity in two adjacent wavelength bands, A and B, centered on the LSPR peak. It considers both the peak shift and the relative intensity change in these two bands, resulting in a significant LSPR signal enhancement. Then the spectral-image contrast was computed as the signal response. Using this approach, we obtained excellent limits of detection (LODs) of 1.4, 1.5, and 4.5 ppb for the neonicotinoids imidacloprid, acetamiprid, and dinotefuran, respectively. Blind tests demonstrated high positive and negative rates for teas, approximately 85 and 100%, respectively. Recombinant OBP2 produced in E. coli offers several advantages over antibodies, including high yield, time savings, and cost effectiveness. Moreover, this method is highly selective and sensitive to neonicotinoids, making it practical for field use.

A Curvilinear-Path Umbrella Sampling Approach to Characterizing the Interactions Between Rapamycin and Three FKBP12 Variants.

Rapamycin is an immunosuppressant macrolide that exhibits anti-proliferative properties through inhibiting the mTOR kinase. In fact, the drug first associates with the FKBP12 enzyme before interacting with the FRB domain of its target. Despite the availability of structural and thermodynamic information on the interaction of FKBP12 with rapamycin, the energetic and mechanistic understanding of this process is still incomplete. We recently reported a multiple-walker umbrella sampling simulation approach to characterizing the protein-protein interaction energetics along curvilinear paths. In the present paper, we extend our investigations to a protein-small molecule duo, the FKBP12•rapamycin complex. We estimate the binding free energies of rapamycin with wild-type FKBP12 and two mutants in which a hydrogen bond has been removed, D37V and Y82F. Furthermore, the underlying mechanistic details are analyzed. The calculated standard free energies of binding agree well with the experimental data, and the roles of the hydrogen bonds are shown to be quite different for each of these two mutated residues. On one hand, removing the carboxylate group of D37 strongly destabilizes the association; on the other hand, the hydroxyl group of Y82 is nearly unnecessary for the stability of the complex because some nonconventional, cryptic, indirect interaction mechanisms seem to be at work.

Low-resolution description of the conformational space for intrinsically disordered proteins.

Intrinsically disordered proteins (IDP) are at the center of numerous biological processes, and attract consequently extreme interest in structural biology. Numerous approaches have been developed for generating sets of IDP conformations verifying a given set of experimental measurements. We propose here to perform a systematic enumeration of protein conformations, carried out using the TAiBP approach based on distance geometry. This enumeration was performed on two proteins, Sic1 and pSic1, corresponding to unphosphorylated and phosphorylated states of an IDP. The relative populations of the obtained conformations were then obtained by fitting SAXS curves as well as Ramachandran probability maps, the original finite mixture approach RamaMix being developed for this second task. The similarity between profiles of local gyration radii provides to a certain extent a converged view of the Sic1 and pSic1 conformational space. Profiles and populations are thus proposed for describing IDP conformations. Different variations of the resulting gyration radius between phosphorylated and unphosphorylated states are observed, depending on the set of enumerated conformations as well as on the methods used for obtaining the populations.

Robust scoring functions for protein-ligand interactions with quantum chemical charge models.

Ordinary least-squares (OLS) regression has been used widely for constructing the scoring functions for protein-ligand interactions. However, OLS is very sensitive to the existence of outliers, and models constructed using it are easily affected by the outliers or even the choice of the data set. On the other hand, determination of atomic charges is regarded as of central importance, because the electrostatic interaction is known to be a key contributing factor for biomolecular association. In the development of the AutoDock4 scoring function, only OLS was conducted, and the simple Gasteiger method was adopted. It is therefore of considerable interest to see whether more rigorous charge models could improve the statistical performance of the AutoDock4 scoring function. In this study, we have employed two well-established quantum chemical approaches, namely the restrained electrostatic potential (RESP) and the Austin-model 1-bond charge correction (AM1-BCC) methods, to obtain atomic partial charges, and we have compared how different charge models affect the performance of AutoDock4 scoring functions. In combination with robust regression analysis and outlier exclusion, our new protein-ligand free energy regression model with AM1-BCC charges for ligands and Amber99SB charges for proteins achieve lowest root-mean-squared error of 1.637 kcal/mol for the training set of 147 complexes and 2.176 kcal/mol for the external test set of 1427 complexes. The assessment for binding pose prediction with the 100 external decoy sets indicates very high success rate of 87% with the criteria of predicted root-mean-squared deviation of less than 2 Å. The success rates and statistical performance of our robust scoring functions are only weakly class-dependent (hydrophobic, hydrophilic, or mixed).

SLITHER: a web server for generating contiguous conformations of substrate molecules entering into deep active sites of proteins or migrating through channels in membrane transporters.

Many proteins use a long channel to guide the substrate or ligand molecules into the well-defined active sites for catalytic reactions or for switching molecular states. In addition, substrates of membrane transporters can migrate to another side of cellular compartment by means of certain selective mechanisms. SLITHER (http://bioinfo.mc.ntu.edu.tw/slither/or http://slither.rcas.sinica.edu.tw/) is a web server that can generate contiguous conformations of a molecule along a curved tunnel inside a protein, and the binding free energy profile along the predicted channel pathway. SLITHER adopts an iterative docking scheme, which combines with a puddle-skimming procedure, i.e. repeatedly elevating the potential energies of the identified global minima, thereby determines the contiguous binding modes of substrates inside the protein. In contrast to some programs that are widely used to determine the geometric dimensions in the ion channels, SLITHER can be applied to predict whether a substrate molecule can crawl through an inner channel or a half-channel of proteins across surmountable energy barriers. Besides, SLITHER also provides the list of the pore-facing residues, which can be directly compared with many genetic diseases. Finally, the adjacent binding poses determined by SLITHER can also be used for fragment-based drug design.

The effect of different colored light emitting diode illumination on egg laying performance, egg qualities, blood hormone levels and behavior patterns in Brown Tsaiya duck.

OBJECTIVE: The objective of this experiment was to investigate the effects of different colors produced by light emitting diode (LED) on Brown Tsaiya ducks. METHODS: A total of 144 female Brown Tsaiya ducks were randomly allocated into three individual cage rearing chambers with different LED illumination colors as treatments. Three different treatments were: i) white color, ii) blue color, and iii) red color. The experiment periods were from ducks 21 to 49 weeks of age, determined traits included i) egg laying performance, ii) feed intake, iii) egg shell breaking strength, iv) egg shell thickness, v) egg Haugh unit, vi) egg weight, vii) serum Estradiol and Progesterone concentration, and viii) behavior pattern. RESULTS: The results indicated that when compared with white and blue color, red color could stimulate ducks sexual maturation and raised the egg laying performance. The red light group was also observed to have the highest feed intake among three treatments. The blue treatment had the lowest egg shell breaking strength and the highest egg weight among three treatments, nevertheless, no significant difference was observed among three treatments on egg shell thickness and egg Haugh unit. The red light group had higher serum estradiol concentration than the white and blue groups, but no significant difference among treatments on the serum Progesterone concentration was found. The results of behavior pattern indicated that red light group showed more feeding and less resting behavior compared to the blue light group. CONCLUSION: We found a potential of applying red light illumination in the indoor laying duck raising system with positive results on egg laying performance and acceptable egg weight, equivalent egg qualities compared to white and blue light.

Molecular dynamics simulations of the rotary motor F(0) under external electric fields across the membrane.

The membrane-bound component F(0), which is a major component of the F(0)F(1)-ATP synthase, works as a rotary motor and plays a central role in driving the F(1) component to transform chemiosmotic energy into ATP synthesis. We conducted molecular dynamics simulations of b(2)-free F(0) in a 1-palmitoyl-2-oleoyl-phosphatidylcholine lipid bilayer for tens of nanoseconds with two different protonation states of the cAsp-61 residue at the interface of the a-c complex in the absence of electric fields and under electric fields of +/-0.03 V/nm across the membrane. To our surprise, we observed that the upper half of the N-terminal helix of the c(1) subunit rotated about its axis clockwise by 30 degrees . An energetic analysis revealed that the electrostatic repulsion between this N-terminal helix and subunit c(12) was a major contributor to the observed rotation. A correlation map analysis indicated that the correlated motions of residues in the interface of the a-c complex were significantly reduced by external electric fields. The deuterium order parameter (S(CD)) profile calculated by averaging all the lipids in the F(0)-bound bilayer was not very different from that of the pure bilayer system, in agreement with recent (2)H solid-state NMR experiments. However, by delineating the lipid properties according to their vicinity to F(0), we found that the S(CD) profiles of different lipid shells were prominently different. Lipids close to F(0) formed a more ordered structure. Similarly, the lateral diffusion of lipids on the membrane surface also followed a shell-dependent behavior. The lipids in the proximity of F(0) exhibited very significantly reduced diffusional motion. The numerical value of S(CD) was anticorrelated with that of the diffusion coefficient, i.e., the more ordered lipid structures led to slower lipid diffusion. Our findings will help elucidate the dynamics of F(0) depending on the protonation state and electric field, and may also shed some light on the interactions between the motor F(0) and its surrounding lipids under physiological conditions, which could help to rationalize its extraordinary energy conversion efficiency.

Inhibition of histone deacetylase activity is a novel function of the antifolate drug methotrexate.

Methotrexate (MTX) is a dihydrofolate reductase (DHFR) inhibitor widely used for treating human cancers, and overexpression of histone deacetylase (HDAC) is usually found in tumors. HDAC inhibitors (HDACi) can reactivate tumor suppressor genes and serve as potential anti-cancer drugs. In this study, we found that MTX shared structural similarity with some HDACi and molecular modeling showed that MTX indeed docks into the active site of HDLP, a bacterial homologue of HDAC. Subsequent in vitro assay demonstrated MTX's inhibition on HDAC activity in human cancer cells. The global acetylation of histone H3 was also induced by MTX. Moreover, MTX inhibited immunoprecipitated HDAC1/2 activity but not their protein levels. This study provides evidence that MTX inhibits HDAC activity.

Quantitative phosphoproteomic analysis identifies the potential therapeutic target EphA2 for overcoming sorafenib resistance in hepatocellular carcinoma cells.

Limited therapeutic options are available for advanced-stage hepatocellular carcinoma owing to its poor diagnosis. Drug resistance to sorafenib, the only available targeted agent, is commonly reported. The comprehensive elucidation of the mechanisms underlying sorafenib resistance may thus aid in the development of more efficacious therapeutic agents. To clarify the signaling changes contributing to resistance, we applied quantitative phosphoproteomics to analyze the differential phosphorylation changes between parental and sorafenib-resistant HuH-7 cells. Consequently, an average of ~1500 differential phosphoproteins were identified and quantified, among which 533 were significantly upregulated in resistant cells. Further bioinformatic integration via functional categorization annotation, pathway enrichment and interaction linkage analysis led to the discovery of alterations in pathways associated with cell adhesion and motility, cell survival and cell growth and the identification of a novel target, EphA2, in resistant HuH-7R cells. In vitro functional analysis indicated that the suppression of EphA2 function impairs cell proliferation and motility and, most importantly, overcomes sorafenib resistance. The attenuation of sorafenib resistance may be achieved prior to its development through the modulation of EphA2 and the subsequent inhibition of Akt activity. Binding analyses and in silico modeling revealed a ligand mimic lead compound, prazosin, that could abate the ligand-independent oncogenic activity of EphA2. Finally, data obtained from in vivo animal models verified that the simultaneous inhibition of EphA2 with sorafenib treatment can effectively overcome sorafenib resistance and extend the projected survival of resistant tumor-bearing mice. Thus our findings regarding the targeting of EphA2 may provide an effective approach for overcoming sorafenib resistance and may contribute to the management of advanced hepatocellular carcinoma.

Protemot: prediction of protein binding sites with automatically extracted geometrical templates.

UNLABELLED: Geometrical analysis of protein tertiary substructures has been an effective approach employed to predict protein binding sites. This article presents the Protemot web server that carries out prediction of protein binding sites based on the structural templates automatically extracted from the crystal structures of protein-ligand complexes in the PDB (Protein Data Bank). The automatic extraction mechanism is essential for creating and maintaining a comprehensive template library that timely accommodates to the new release of PDB as the number of entries continues to grow rapidly. The design of Protemot is also distinctive by the mechanism employed to expedite the analysis process that matches the tertiary substructures on the contour of the query protein with the templates in the library. This expediting mechanism is essential for providing reasonable response time to the user as the number of entries in the template library continues to grow rapidly due to rapid growth of the number of entries in PDB. This article also reports the experiments conducted to evaluate the prediction power delivered by the Protemot web server. Experimental results show that Protemot can deliver a superior prediction power than a web server based on a manually curated template library with insufficient quantity of entries. AVAILABILITY: http://protemot.csie.ntu.edu.tw/step1.cgi http://bioinfo.mc.ntu.edu.tw/protemot/step1.cgi.

Effect of statin use on the risk of medically attended acute respiratory illness among influenza vaccinated elderly.

OBJECTIVES: The immunomodulatory effects of statins may reduce the immune response induced by influenza vaccines. However, evidence regarding the effect of statin use on the effectiveness of seasonal influenza vaccines against medically attended acute respiratory illness (MAARI) in the elderly remains scarce. METHODS: We conducted a retrospective cohort study using data from Taiwan's National Health Insurance Research Database. Elderly adults aged ≧ 66 years who were vaccinated with seasonal influenza vaccines during the 2007-2008 to 2012-2013 influenza seasons were enrolled for this analysis. We compared the risk of MAARI between statin and non-statin users. Propensity score matching and conditional logistic regression models were used to analyze the data. RESULTS: A total of 440,180 elderly were included in this study. In general, the risk of MAARI was higher in statin users than non-statin users (odds ratio [OR]: 1.03, 95% confidence interval [CI]: 1.02-1.05). Statin exposure after vaccination was associated with a higher risk of MAARI (OR: 1.05, 95% CI: 1.02-1.07). Among different statin agents, simvastatin and lovastatin use was associated with a significant increase in the risk of MAARI (ORsimvastatin: 1.14, 95% CI: 1.10-1.18; ORlovastatin: 1.18, 95% CI: 1.12-1.25). CONCLUSIONS: Statin exposure, especially simvastatin and lovastatin, was associated with a higher risk of MAARI in the seasonal influenza vaccinated elderly. Future studies exploring the differences between individual statins and mechanisms of their immunomodulatory effects are necessary.

Structural insights into the gating of DNA passage by the topoisomerase II DNA-gate.

Type IIA topoisomerases (Top2s) manipulate the handedness of DNA crossovers by introducing a transient and protein-linked double-strand break in one DNA duplex, termed the DNA-gate, whose opening allows another DNA segment to be transported through to change the DNA topology. Despite the central importance of this gate-opening event to Top2 function, the DNA-gate in all reported structures of Top2-DNA complexes is in the closed state. Here we present the crystal structure of a human Top2 DNA-gate in an open conformation, which not only reveals structural characteristics of its DNA-conducting path, but also uncovers unexpected yet functionally significant conformational changes associated with gate-opening. This structure further implicates Top2's preference for a left-handed DNA braid and allows the construction of a model representing the initial entry of another DNA duplex into the DNA-gate. Steered molecular dynamics calculations suggests the Top2-catalyzed DNA passage may be achieved by a rocker-switch-type movement of the DNA-gate.

MEDock: a web server for efficient prediction of ligand binding sites based on a novel optimization algorithm.

The prediction of ligand binding sites is an essential part of the drug discovery process. Knowing the location of binding sites greatly facilitates the search for hits, the lead optimization process, the design of site-directed mutagenesis experiments and the hunt for structural features that influence the selectivity of binding in order to minimize the drug's adverse effects. However, docking is still the rate-limiting step for such predictions; consequently, much more efficient algorithms are required. In this article, the design of the MEDock web server is described. The goal of this sever is to provide an efficient utility for predicting ligand binding sites. The MEDock web server incorporates a global search strategy that exploits the maximum entropy property of the Gaussian probability distribution in the context of information theory. As a result of the global search strategy, the optimization algorithm incorporated in MEDock is significantly superior when dealing with very rugged energy landscapes, which usually have insurmountable barriers. This article describes four different benchmark cases that span a diverse set of different types of ligand binding interactions. These benchmarks were compared with the use of the Lamarckian genetic algorithm (LGA), which is the major workhorse of the well-known AutoDock program. These results demonstrate that MEDock consistently converged to the correct binding modes with significantly smaller numbers of energy evaluations than the LGA required. When judged by a threshold of the number of energy evaluations consumed in the docking simulation, MEDock also greatly elevates the rate of accurate predictions for all benchmark cases. MEDock is available at http://medock.csie.ntu.edu.tw/ and http://bioinfo.mc.ntu.edu.tw/medock/.