4D-QSAR analysis of a set of propofol analogues: mapping binding sites for an anesthetic phenol on the GABA(A) receptor

CYCLOBUTANES IN SMALL MOLECULE DRUG CANDIDATES

Cyclobutanes are increasingly used in medicinal chemistry in the search for relevant biological properties. Important characteristics of the cyclobutane ring include its unique puckered structure, longer C−C bond lengths, increased C−C π‐character and relative chemical inertness for a highly strained carbocycle. This review will focus on contributions of cyclobutane rings in drug candidates to arrive at favorable properties. Cyclobutanes have been employed for improving multiple factors such as preventing cis/trans‐isomerization by replacing alkenes, replacing larger cyclic systems, increasing metabolic stability, directing key pharmacophore groups, inducing conformational restriction, reducing planarity, as aryl isostere and filling hydrophobic pockets.

Two Decades of 4D-QSAR: A Dying Art or Staging a Comeback?

A key question confronting computational chemists concerns the preferable ligand geometry that fits complementarily into the receptor pocket. Typically, the postulated ‘bioactive’ 3D ligand conformation is constructed as a ‘sophisticated guess’ (unnecessarily geometry-optimized) mirroring the pharmacophore hypothesis—sometimes based on an erroneous prerequisite. Hence, 4D-QSAR scheme and its ‘dialects’ have been practically implemented as higher level of model abstraction that allows the examination of the multiple molecular conformation, orientation and protonation representation, respectively. Nearly a quarter of a century has passed since the eminent work of Hopfinger appeared on the stage; therefore the natural question occurs whether 4D-QSAR approach is still appealing to the scientific community? With no intention to be comprehensive, a review of the current state of art in the field of receptor-independent (RI) and receptor-dependent (RD) 4D-QSAR methodology is provided with a brief examination of the ‘mainstream’ algorithms. In fact, a myriad of 4D-QSAR methods have been implemented and applied practically for a diverse range of molecules. It seems that, 4D-QSAR approach has been experiencing a promising renaissance of interests that might be fuelled by the rising power of the graphics processing unit (GPU) clusters applied to full-atom MD-based simulations of the protein-ligand complexes.

Syntheses of Analogues of Propofol: A Review

Propofol (2,6-diisopropylphenol) is an intravenous sedative/hypnotic agent that is used extensively for introduction and maintenance of general anaesthesia, sedation of critically ill patients and procedural sedation (e.g., endoscopy). Propofol has a rapid onset and offset of action and shows only minimal accumulation upon prolonged use. Propofol is only sparingly soluble in water and is currently marketed in 10% soybean oil-based lipid emulsion. Propofol’s anaesthetic properties were discovered over forty years ago, and it has been in clinical use for over thirty years. The main use of propofol remains as an anaesthetic but, over the years, analogues have been developed with varying properties from anticancer, anticonvulsant and antioxidant. In addition, large synthetic efforts have been made towards improving propofol’s water-solubility, its activity, and elucidating its structure–activity­ relationship and exact mechanism of action have been made. This review provides an overview of the research pertaining to propofol-like molecules and covers the efforts of synthetic chemists towards propofol analogues over the last 40 years.

1 Introduction

2 History

3 Early Work

4 Improving Water Solubility

5 The Importance of the Phenol

6 Exploring the Structure–Activity Relationship and Attempts to Improve Activity

7 Anticancer Activity

8 Anticonvulsant Properties

9 Antioxidant Activity

10 Photoactive Labelling to Elucidate Mechanism of Action

11 Photoregulation

12 Conclusion

Shared structural mechanisms of general anaesthetics and benzodiazepines

Most general anaesthetics and classical benzodiazepine drugs act through positive modulation of γ-aminobutyric acid type A (GABAA) receptors to dampen neuronal activity in the brain1-5. However, direct structural information on the mechanisms of general anaesthetics at their physiological receptor sites is lacking. Here we present cryo-electron microscopy structures of GABAA receptors bound to intravenous anaesthetics, benzodiazepines and inhibitory modulators. These structures were solved in a lipidic environment and are complemented by electrophysiology and molecular dynamics simulations. Structures of GABAA receptors in complex with the anaesthetics phenobarbital, etomidate and propofol reveal both distinct and common transmembrane binding sites, which are shared in part by the benzodiazepine drug diazepam. Structures in which GABAA receptors are bound by benzodiazepine-site ligands identify an additional membrane binding site for diazepam and suggest an allosteric mechanism for anaesthetic reversal by flumazenil. This study provides a foundation for understanding how pharmacologically diverse and clinically essential drugs act through overlapping and distinct mechanisms to potentiate inhibitory signalling in the brain.

A novel hybrid method named electron conformational genetic algorithm as a 4D QSAR investigation to calculate the biological activity of the tetrahydrodibenzazosines

To understand the structure-activity correlation of a group of tetrahydrodibenzazocines as inhibitors of 17β-hydroxysteroid dehydrogenase type 3, we have performed a combined genetic algorithm (GA) and four-dimensional quantitative structure-activity relationship (4D-QSAR) modeling study. The computed electronic and geometry structure descriptors were regulated as a matrix and named as electron-conformational matrix of contiguity (ECMC). A chemical property-based pharmacophore model was developed for series of tetrahydrodibenzazocines by EMRE software package. GA was employed to choose an optimal combination of parameters. A model has been developed for estimating anticancer activity quantitatively. All QSAR models were established with 40 compounds (training set), then they were considered for selective capability with additional nine compounds (test set). A statistically valid 4D-QSAR ( R training 2 = 0.856 , R test 2 = 0.851 and q² = 0.650) with good external set prediction was obtained.

Identifying Drugs that Bind Selectively to Intersubunit General Anesthetic Sites in the α1β3γ2 GABAAR Transmembrane Domain

GABA_A receptors (GABA_ARs) are targets for important classes of clinical agents (e.g., anxiolytics, anticonvulsants, and general anesthetics) that act as positive allosteric modulators (PAMs). Previously, using photoreactive analogs of etomidate ([³H]azietomidate) and mephobarbital [[³H]1-methyl-5-allyl-5-(m-trifluoromethyl-diazirynylphenyl)barbituric acid ([³H]R-mTFD-MPAB)], we identified two homologous but pharmacologically distinct classes of general anesthetic binding sites in the α1β3γ2 GABA_AR transmembrane domain at β ⁺-α ^- (β ⁺ sites) and α ⁺-β ^-/γ ⁺-β ^- (β ^- sites) subunit interfaces. We now use competition photolabeling with [³H]azietomidate and [³H]R-mTFD-MPAB to identify para-substituted propofol analogs and other drugs that bind selectively to intersubunit anesthetic sites. Propofol and 4-chloro-propofol bind with 5-fold selectivity to β ⁺, while derivatives with bulkier lipophilic substitutions [4-(tert-butyl)-propofol and 4-(hydroxyl(phenyl)methyl)-propofol] bind with ∼10-fold higher affinity to β ^- sites. Similar to R-mTFD-MPAB and propofol, these drugs bind in the presence of GABA with similar affinity to the α ⁺-β ^- and γ ⁺-β ^- sites. However, we discovered four compounds that bind with different affinities to the two β ^- interface sites. Two of these bind with higher affinity to one of the β ^- sites than to the β ⁺ sites. We deduce that 4-benzoyl-propofol binds with >100-fold higher affinity to the γ ⁺-β ^- site than to the α ⁺-β ^- or β ⁺-α ^- sites, whereas loreclezole, an anticonvulsant, binds with 5- and 100-fold higher affinity to the α ⁺-β ^- site than to the β ⁺ and γ ⁺-β ^- sites. These studies provide a first identification of PAMs that bind selectively to a single intersubunit site in the GABA_AR transmembrane domain, a property that may facilitate the development of subtype selective GABA_AR PAMs.

Alkylphenol inverse agonists of HCN1 gating: H-bond propensity, ring saturation and adduct geometry differentially determine efficacy and potency

BACKGROUND AND PURPOSE

In models of neuropathic pain, inhibition of HCN1 is anti-hyperalgesic. 2,6-di-iso-propyl phenol (propofol) and its non-anesthetic congener, 2,6-di-tert-butyl phenol, inhibit HCN1 channels by stabilizing closed state(s).

EXPERIMENTAL APPROACH

Using in vitro electrophysiology and kinetic modeling, we systematically explore the contribution of ligand architecture to alkylphenol-channel coupling.

KEY RESULTS

When corrected for changes in hydrophobicity (and propensity for intra-membrane partitioning), the decrease in potency upon 1-position substitution (NCO∼OH >> SH >>> F) mirrors the ligands' H-bond acceptor (NCO > OH > SH >>> F) but not donor profile (OH > SH >>> NCO∼F). H-bond elimination (OH to F) corresponds to a ΔΔG of ∼4.5 kCal mol^-1 loss of potency with little or no disruption of efficacy. Substitution of compact alkyl groups (iso-propyl, tert-butyl) with shorter (ethyl, methyl) or more extended (sec-butyl) adducts disrupts both potency and efficacy. Ring saturation (with the obligate loss of both planarity and π electrons) primarily disrupts efficacy.

CONCLUSIONS AND IMPLICATIONS

A hydrophobicity-independent decrement in potency at higher volumes suggests the alkylbenzene site has a volume of ≥800 ų. Within this, a relatively static (with respect to ligand) H-bond donor contributes to initial binding with little involvement in generation of coupling energy. The influence of π electrons/ring planarity and alkyl adducts on efficacy reveals these aspects of the ligand present towards a face of the channel that undergoes structural changes during opening. The site's characteristics suggest it is "druggable"; introduction of other adducts on the ring may generate higher potency inverse agonists.

Interactions between model cell membranes and the neuroactive drug propofol

Vibrational sum frequency spectroscopy (VSFS) complemented by surface pressure isotherm and neutron reflectometry (NR) experiments were employed to investigate the interactions between propofol, a small amphiphilic molecule that currently is the most common general anaesthetic drug, and phospholipid monolayers. A series of biologically relevant saturated phospholipids of varying chain length from C₁₈ to C₁₄ were spread on either pure water or propofol (2,6-bis(1-methylethyl)phenol) solution in a Langmuir trough, and the change in the molecular structure of the film, induced by the interaction with propofol, was studied with respect to the surface pressure. The results from the surface pressure isotherm experiments revealed that propofol, as long as it remains at the interface, enhances the fluidity of the phospholipid monolayer. The VSF spectra demonstrate that for each phospholipid the amount of propofol in the monolayer region decreases with increasing surface pressure. Such squeeze out is in contrast to the enhanced interactions that can be exhibited by more complex amphiphilic molecules such as peptides. At surface pressures of 22-25 mN m^-1, which are relevant for biological cell membranes, most of the propofol has been expelled from the monolayer, especially in the case of the C₁₆ and C₁₈ phospholipids that adopt a liquid condensed phase packing of its alkyl tails. At lower surface pressures of 5 mN m^-1, the effect of propofol on the structure of the alkyl tails is enhanced when the phospholipids are present in a liquid expanded phase. Specifically, for the C₁₆ phospholipid, NR data reveal that propofol is located exclusively in the head group region, which is rationalized in the context of previous studies. The results imply a non-homogeneous distribution of propofol in the plane of real cell membranes, which is an inference that requires urgent testing and may help to explain why such low concentration of the drug are required to induce general anaesthesia.

A receptor dependent-4D QSAR approach to predict the activity of mutated enzymes

Screening and selection tools to obtain focused libraries play a key role in successfully engineering enzymes of desired qualities. The quality of screening depends on efficient assays; however, a focused library generated with a priori information plays a major role in effectively identifying the right enzyme. As a proof of concept, for the first time, receptor dependent - 4D Quantitative Structure Activity Relationship (RD-4D-QSAR) has been implemented to predict kinetic properties of an enzyme. The novelty of this study is that the mutated enzymes also form a part of the training data set. The mutations were modeled in a serine protease and molecular dynamics simulations were conducted to derive enzyme-substrate (E-S) conformations. The E-S conformations were enclosed in a high resolution grid consisting of 156,250 grid points that stores interaction energies to generate QSAR models to predict the enzyme activity. The QSAR predictions showed similar results as reported in the kinetic studies with >80% specificity and >50% sensitivity revealing that the top ranked models unambiguously differentiated enzymes with high and low activity. The interaction energy descriptors of the best QSAR model were used to identify residues responsible for enzymatic activity and substrate specificity.

The Role of the Hydroxyl Group in Propofol-Protein Target Recognition: Insights from ONIOM Studies

Propofol (PFL, 1-hydroxyl-2,6-diisopropylbenzene) is currently used widely as one of the most well-known intravenous anesthetics to relieve surgical suffering, but its mechanism of action is not yet clear. Previous experimental studies have demonstrated that the hydroxyl group of PFL plays a dominant role in the molecular recognition of PFL with receptors that lead to hypnosis. To further explore the mechanism of anesthesia induced by PFL in the present work, the exact binding features and interaction details of PFL with three important proteins, human serum albumin (HSA), the pH-gated ion channel from Gloeobacter violaceus (GLIC), and horse spleen apoferritin (HSAF), were investigated systematically by using a rigorous three-layer ONIOM (M06-2X/6-31+G*:PM6:AMBER) method. Additionally, to further characterize the possible importance of such hydroxyl interactions, a similar set of calculations was carried out on the anesthetically inactive fropofol (FFL, 1-fluoro-2,6-diisopropylbenzene) in which the fluorine was substituted for the hydroxyl. According to the ONIOM calculations, atoms in molecules (AIM) analyses, and electrostatic potential (ESP) analyses, the significance of hydrogen bond, halogen bond, and hydrophobic interactions in promoting proper molecular recognition was revealed. The binding interaction energies of PFL with different proteins were generally larger than FFL and are a significant determinant of their differential anesthetic efficacies. Interestingly, although the hydrogen-bonding effect of the hydroxyl moiety was prominent in propofol, the substitution of the 1-hydroxyl by a fluorine atom did not prevent FFL from binding to the protein via a halogen-bonding interaction. It therefore became clear that multiple specific interactions rather than just hydrogen or halogen bonds must be taken into account to explain the different anesthesia endpoints caused by PFL and FFL. The contributions of key residues in ligand-receptor binding were also quantified, and the calculated results agreed with many available experimental observations. This work will provide complementary insights into the molecular mechanisms of anesthetic action for PFL from a robust theoretical point of view. This will not only assist in interpreting experimental observations but will also help to develop working hypotheses for further experiments and future drug design.

Curcumol allosterically modulates GABA(A) receptors in a manner distinct from benzodiazepines

Inhibitory A type γ-aminobutyric acid receptors (GABA_ARs) play a pivotal role in orchestrating various brain functions and represent an important molecular target in neurological and psychiatric diseases, necessitating the need for the discovery and development of novel modulators. Here, we show that a natural compound curcumol, acts as an allosteric enhancer of GABA_ARs in a manner distinct from benzodiazepines. Curcumol markedly facilitated GABA-activated currents and shifted the GABA concentration-response curve to the left in cultured hippocampal neurons. When co-applied with the classical benzodiazepine diazepam, curcumol further potentiated GABA-induced currents. In contrast, in the presence of a saturating concentration of menthol, a positive modulator for GABA_AR, curcumol failed to further enhance GABA-induced currents, suggesting shared mechanisms underlying these two agents on GABA_ARs. Moreover, the benzodiazepine antagonist flumazenil did not alter the enhancement of GABA response by curcumol and menthol, but abolished that by DZP. Finally, mutations at the β2 or γ2 subunit predominantly eliminated modulation of recombinant GABA_ARs by curcumol and menthol, or diazepam, respectively. Curcumol may therefore exert its actions on GABA_ARs at sites distinct from benzodiazepine sites. These findings shed light on the future development of new therapeutics drugs targeting GABA_ARs.

A Novel Bifunctional Alkylphenol Anesthetic Allows Characterization of γ-Aminobutyric Acid, Type A (GABAA), Receptor Subunit Binding Selectivity in Synaptosomes

Propofol, an intravenous anesthetic, is a positive modulator of the GABAA receptor, but the mechanistic details, including the relevant binding sites and alternative targets, remain disputed. Here we undertook an in-depth study of alkylphenol-based anesthetic binding to synaptic membranes. We designed, synthesized, and characterized a chemically active alkylphenol anesthetic (2-((prop-2-yn-1-yloxy)methyl)-5-(3-(trifluoromethyl)-3H-diazirin-3-yl)phenol, AziPm-click (1)), for affinity-based protein profiling (ABPP) of propofol-binding proteins in their native state within mouse synaptosomes. The ABPP strategy captured ∼4% of the synaptosomal proteome, including the unbiased capture of five α or β GABAA receptor subunits. Lack of γ2 subunit capture was not due to low abundance. Consistent with this, independent molecular dynamics simulations with alchemical free energy perturbation calculations predicted selective propofol binding to interfacial sites, with higher affinities for α/β than γ-containing interfaces. The simulations indicated hydrogen bonding is a key component leading to propofol-selective binding within GABAA receptor subunit interfaces, with stable hydrogen bonds observed between propofol and α/β cavity residues but not γ cavity residues. We confirmed this by introducing a hydrogen bond-null propofol analogue as a protecting ligand for targeted-ABPP and observed a lack of GABAA receptor subunit protection. This investigation demonstrates striking interfacial GABAA receptor subunit selectivity in the native milieu, suggesting that asymmetric occupancy of heteropentameric ion channels by alkylphenol-based anesthetics is sufficient to induce modulation of activity.

A Cysteine Substitution Probes β3H267 Interactions with Propofol and Other Potent Anesthetics in α1β3γ2L γ-Aminobutyric Acid Type A Receptors

BACKGROUND

Anesthetic contact residues in γ-aminobutyric acid type A (GABAA) receptors have been identified using photolabels, including two propofol derivatives. O-propofol diazirine labels H267 in β3 and α1β3 receptors, whereas m-azi-propofol labels other residues in intersubunit clefts of α1β3. Neither label has been studied in αβγ receptors, the most common isoform in mammalian brain. In αβγ receptors, other anesthetic derivatives photolabel m-azi-propofol-labeled residues, but not βH267. The authors' structural homology model of α1β3γ2L receptors suggests that β3H267 may abut some of these sites.

METHODS

Substituted cysteine modification-protection was used to test β3H267C interactions with four potent anesthetics: propofol, etomidate, alphaxalone, and R-5-allyl-1-methyl-5-(m-trifluoromethyl-diazirinylphenyl) barbituric acid (mTFD-MPAB). The authors expressed α1β3γ2L or α1β3H267Cγ2L GABAA receptors in Xenopus oocytes. The authors used voltage clamp electrophysiology to assess receptor sensitivity to γ-aminobutyric acid (GABA) and anesthetics and to compare p-chloromercuribenzenesulfonate modification rates with GABA versus GABA plus anesthetics.

RESULTS

Enhancement of low GABA (eliciting 5% of maximum) responses by equihypnotic concentrations of all four anesthetics was similar in α1β3γ2L and α1β3H267Cγ2L receptors (n > 3). Direct activation of α1β3H267Cγ2L receptors, but not α1β3γ2L, by mTFD-MPAB and propofol was significantly greater than the other anesthetics. Modification of β3H267C by p-chloromercuribenzenesulfonate (n > 4) was rapid and accelerated by GABA. Only mTFD-MPAB slowed β3H267C modification (approximately twofold; P = 0.011).

CONCLUSIONS

β3H267 in α1β3γ2L GABAA receptors contacts mTFD-MPAB, but not propofol. The study results suggest that β3H267 is near the periphery of one or both transmembrane intersubunit (α+/β- and γ+/β-) pockets where both mTFD-MPAB and propofol bind.

Long-term action of propofol on cognitive function and hippocampal neuroapoptosis in neonatal rats

Propofol is a short-acting anesthetic and generally is utilized for the induction and maintenance of anesthesia in pediatrics and adults. However, whether repeated use of propofol affects long-term cognitive function remains unclear. This study investigated the effects of propofol on cognitive function and hippocampal neuroapoptosis in neonatal rat. A total of 112 male newborn 7-day old Sprague-Dawley rats were randomly divided into 8 groups (n=14 rats per group) and intraperitoneally injected either with saline or propofol at 50, 100, and 150 mg/kg/day for 5 consecutive days. Four non-surgical groups were assigned as Con1, P50, P100, and P150. Four surgical groups were received an appendicectomy under propofol anesthesia and assigned as Con2, SP50, SP100, SP150. After 2 months raising, cognitive function, hippocampal neuroapoptosis, and intracephalic inflammatory cytokines were evaluated. There was no obvious effect on the cognitive function and neuroapoptosis after repeated use of propofol at a low dose for 5 days, whereas repeated use of propofol at a middle/high dose significantly increase the expression of apoptotic factors (caspase-3 and Bax), pro-inflammatory cytokines (IL-1β, IL-6, and TNF-α), and impair the cognitive function. Thus, our data suggest that repeated use of propofol at a low dose may be safe during the period of brain growth spurt. Using propofol at a recommended or higher dose for anaesthesia may lead to the cognitive defects, attributed to hippocampal neuroapoptosis and the overexpression of pro-inflammatory cytokines in the brain.

Role for the propofol hydroxyl in anesthetic protein target molecular recognition

Propofol is a widely used intravenous general anesthetic. We synthesized 2-fluoro-1,3-diisopropylbenzene, a compound that we call "fropofol", to directly assess the significance of the propofol 1-hydroxyl for pharmacologically relevant molecular recognition in vitro and for anesthetic efficacy in vivo. Compared to propofol, fropofol had a similar molecular volume and only a small increase in hydrophobicity. Isothermal titration calorimetry and competition assays revealed that fropofol had higher affinity for a protein site governed largely by van der Waals interactions. Within another protein model containing hydrogen bond interactions, propofol demonstrated higher affinity. In vivo, fropofol demonstrated no anesthetic efficacy, but at high concentrations produced excitatory activity in tadpoles and mice; fropofol also antagonized propofol-induced hypnosis. In a propofol protein target that contributes to hypnosis, α1β2γ2L GABAA receptors, fropofol demonstrated no significant effect alone or on propofol positive allosteric modulation of the ion channel, suggesting an additional requirement for the 1-hydroxyl within synaptic GABAA receptor site(s). However, fropofol caused similar adverse cardiovascular effects as propofol by a dose-dependent depression of myocardial contractility. Our results directly implicate the propofol 1-hydroxyl as contributing to molecular recognition within protein targets leading to hypnosis, but not necessarily within protein targets leading to side effects of the drug.

Impact of malnutrition on propofol consumption and recovery time among patients undergoing laparoscopic gastrointestinal surgery

BACKGROUND

Malnutrition is a major health problem, especially in hospitalized patients as it can be closely related to many post-operative complications. However, research on malnutrition and its effect on the outcome of general anesthesia have been largely neglected. Here we investigated malnutrition status on propofol consumption and recovery time among patients undergoing laparoscopic gastrointestinal surgery under general anesthesia.

METHODS

One hundred and one patients were recruited between January and June 2012 at Tongji Hospital and assigned into three groups according to Nutritional Risk Screening Tool 2002 score. A standard combined general anesthesia procedure was performed under regular monitoring. The dosage of propofol needed for induction, consumption during maintenance and recovery time were recorded.

RESULTS

When compared with normal nutritional status individuals, the propofol dosage at induction was significantly decreased about 4.3% in moderate malnutritional status patients (P < 0.01) and about 16.8% in severely malnutritional status patients (P < 0.01). The average consumption of propofol was also significantly lower in malnourished individuals; for moderate malnutritional, the decrease was about 20% (P < 0.01) while for the severely malnutritional, it was 30% (P < 0.01) when compared with normal nutritional status individuals. For the recovery time of propofol anesthesia, the patients with severe malnutritional status awoke average 6.8 min later than those normally nourished (P < 0.01), but those patients with moderate malnutrition status did not (P = 0.885).

CONCLUSION

The present results indicate that the dosage and recovery time of propofol does change in malnourished individuals. Therefore, malnutrition may somehow affect the outcome of general anesthesia.

Assessment of homology templates and an anesthetic binding site within the γ-aminobutyric acid receptor

BACKGROUND

Anesthetics mediate portions of their activity via modulation of the γ-aminobutyric acid receptor (GABAaR). Although its molecular structure remains unknown, significant progress has been made toward understanding its interactions with anesthetics via molecular modeling.

METHODS

The structure of the torpedo acetylcholine receptor (nAChRα), the structures of the α4 and β2 subunits of the human nAChR, the structures of the eukaryotic glutamate-gated chloride channel (GluCl), and the prokaryotic pH-sensing channels, from Gloeobacter violaceus and Erwinia chrysanthemi, were aligned with the SAlign and 3DMA algorithms. A multiple sequence alignment from these structures and those of the GABAaR was performed with ClustalW. The Modeler and Rosetta algorithms independently created three-dimensional constructs of the GABAaR from the GluCl template. The CDocker algorithm docked a congeneric series of propofol derivatives into the binding pocket and scored calculated binding affinities for correlation with known GABAaR potentiation EC50s.

RESULTS

Multiple structure alignments of templates revealed a clear consensus of residue locations relevant to anesthetic effects except for torpedo nAChR. Within the GABAaR models generated from GluCl, the residues notable for modulating anesthetic action within transmembrane segments 1, 2, and 3 converged on the intersubunit interface between α and β subunits. Docking scores of a propofol derivative series into this binding site showed strong linear correlation with GABAaR potentiation EC50.

CONCLUSION

Consensus structural alignment based on homologous templates revealed an intersubunit anesthetic binding cavity within the transmembrane domain of the GABAaR, which showed a correlation of ligand docking scores with experimentally measured GABAaR potentiation.

The comparison of preemptive effects of propofol, remifentanil and ketamine on post-operative pain scores and analgesic requirements in elective lower abdominal surgery under general anesthesia: A randomized, double-blinded study

BACKGROUND

In this randomized, double-blinded study, we investigated the preemptive effects of propofol, remifentanil or ketamine on post-operative pain scores and analgesic requirements in elective lower abdominal surgeries under general anesthesia during the first 24 h of post-operative period.

MATERIALS AND METHODS

Seventy five patients, American Society of Anesthesiologists physical status I or II candidate for elective lower abdominal surgery under general anesthesia were randomized to three groups (25 each). According to their allocated group, patients received either propofol 0.25 mg/kg, remifentanil 0.25 mic/kg or ketamine 0.3 mg/kg as preemptive analgesia immediately after the induction of general anesthesia. Post-operative pain scores with a numerical rating scale (visual analogue scale 0-10) were assessed and analgesic requirements and side-effects were compared through analysis using the SPSS version 18 in the post-operative period; post-anesthesia care unit 2, 6, 12 and 24 h.

RESULTS

Patients' demographics were similar in all groups. The pain scores were significantly lower in remifentanil group immediately after recovery and also at 2 and 6 h post-operatively, but it reversed at 12 and 24 h after recovery comparing with propofol and ketamine. However, the mean of administered morphine in the first 24 h was significantly lower in propofol group (18.97 ± 6.6) comparing with remifentanil group (21.96 ± 6.55) and ketamine group (24.26 ± 5.84) (P value, 0.01).

CONCLUSION

Prophylactic preemptive single dose of intravenous (IV) 0.25 mg/kg propofol significantly decreased post-operative analgesia requirements comparing with IV 0.3 mg/kg ketamine or 0.25 μg/kg remifentanil.

Robust photoregulation of GABA(A) receptors by allosteric modulation with a propofol analogue

Photochemical switches represent a powerful method for improving pharmacological therapies and controlling cellular physiology. Here we report the photo-regulation of GABA_A receptors (GABA_ARs) by a derivative of propofol (2,6-diisopropylphenol), a GABA_AR allosteric modulator, that we have modified to contain photo-isomerizable azobenzene. Using α₁β₂γ₂ GABA_ARs expressed in Xenopus laevis oocytes and native GABA_ARs of isolated retinal ganglion cells, we show that the trans-azobenzene isomer of the new compound (trans-MPC088), generated by visible light (wavelengths ~440 nm), potentiates the GABA-elicited response and at higher concentrations directly activates the receptors. cis-MPC088, generated from trans-MPC088 by UV light (~365 nm), produces little if any receptor potentiation/activation. In cerebellar slices, MPC088 co-applied with GABA affords bidirectional photo-modulation of Purkinje cell membrane current and spike-firing rate. The findings demonstrate photo-control of GABA_ARs by an allosteric ligand and open new avenues for fundamental and clinically oriented research on GABA_ARs, a major class of neurotransmitter receptors in the central nervous system.

Quantitative structure-activity relationships of monoterpenoid binding activities to the housefly GABA receptor

BACKGROUND

Monoterpenoids are a large group of plant secondary metabolites. Many of these naturally occurring compounds have shown good insecticidal potency on pest insects. Previous studies in this laboratory have indicated that some monoterpenoids have positive modulatory effects on insect GABA receptors. In this study, the key properties of monoterpenoids involved in monoterpenoid binding activity at the housefly GABA receptor were determined by developing quantitative structure-activity relationship (QSAR) models, and the relationship between the toxicities of these monoterpenoids and their GABA receptor binding activities was evaluated.

RESULTS

Two QSAR models were determined for nine monoterpenoids showing significant effects on [³H]-TBOB binding and for nine p-menthane analogs with at least one oxygen atom attached to the ring. The Mulliken charges on certain carbon atoms, the log P value and the total energy showed significant relationships with binding activities to the housefly GABA receptor in these two QSAR models.

CONCLUSIONS

From the QSAR models, some chemical and structural parameters, including the electronic properties, hydrophobicity and stability of monoterpenoid molecules, were suggested to be strongly involved in binding activities to the housefly GABA receptor. These findings will help to understand the mode of action of these natural insecticides, and provide guidance to predict more monoterpenoid insecticides.

Recognition of anesthetic barbiturates by a protein binding site: a high resolution structural analysis

Barbiturates potentiate GABA actions at the GABA_A receptor and act as central nervous system depressants that can induce effects ranging from sedation to general anesthesia. No structural information has been available about how barbiturates are recognized by their protein targets. For this reason, we tested whether these drugs were able to bind specifically to horse spleen apoferritin, a model protein that has previously been shown to bind many anesthetic agents with affinities that are closely correlated with anesthetic potency. Thiopental, pentobarbital, and phenobarbital were all found to bind to apoferritin with affinities ranging from 10–500 µM, approximately matching the concentrations required to produce anesthetic and GABAergic responses. X-ray crystal structures were determined for the complexes of apoferritin with thiopental and pentobarbital at resolutions of 1.9 and 2.0 Å, respectively. These structures reveal that the barbiturates bind to a cavity in the apoferritin shell that also binds haloalkanes, halogenated ethers, and propofol. Unlike these other general anesthetics, however, which rely entirely upon van der Waals interactions and the hydrophobic effect for recognition, the barbiturates are recognized in the apoferritin site using a mixture of both polar and nonpolar interactions. These results suggest that any protein binding site that is able to recognize and respond to the chemically and structurally diverse set of compounds used as general anesthetics is likely to include a versatile mixture of both polar and hydrophobic elements.

p-(4-Azipentyl)propofol: a potent photoreactive general anesthetic derivative of propofol

We synthesized 2,6-diisopropyl-4-[3-(3-methyl-3H-diazirin-3-yl)propyl]phenol (p-(4-azipentyl)propofol), or p-4-AziC5-Pro, a novel photoactivable derivative of the general anesthetic propofol. p-4-AziC5-Pro has an anesthetic potency similar to that of propofol. Like propofol, the compound potentiates inhibitory GABA(A) receptor current responses and allosterically modulates binding to both agonist and benzodiazepine sites, assayed on heterologously expressed GABA(A) receptors. p-4-AziC5-Pro inhibits excitatory current responses of nACh receptors expressed in Xenopus oocytes and photoincorporates into native nACh receptor-enriched Torpedo membranes. Thus, p-4-AziC5-Pro is a functional general anesthetic that both modulates and photoincorporates into Cys-loop ligand-gated ion channels, making it an excellent candidate for use in identifying propofol binding sites.

Cyclohexanol analogues are positive modulators of GABA(A) receptor currents and act as general anaesthetics in vivo

GABA(A) receptors meet all the pharmacological criteria required to be considered important general anaesthetic targets. In the following study, the modulatory effects of various commercially available and novel cyclohexanols were investigated on recombinant human γ-aminobutyric acid (GABA(A), α(1)β(2)γ(2s)) receptors expressed in Xenopus oocytes, and compared to the modulatory effects on GABA currents observed with exposures to the intravenous anaesthetic agent, propofol. Submaximal EC(20) GABA currents were typically enhanced by co-applications of 3-300 μM cyclohexanols. For instance, at 30 μM 2,6-diisopropylcyclohexanol (a novel compound) GABA responses were increased ~3-fold (although similar enhancements were achieved at 3 μM propofol). As regards rank order for modulation by the cyclohexanol analogues at 30 μM, the % enhancements for 2,6-dimethylcyclohexanol~2,6-diethylcyclohexanol~2,6-diisopropylcyclohexanol~2,6-di-sec-butylcyclohexanol ≫2,6-di-tert-butylcyclohexanol~4-tert-butylcyclohexanol>cyclohexanol~cyclopentanol~2-methylcyclohexanol. We further tested the potencies of the cyclohexanol analogues as general anaesthetics using a tadpole in vivo assay. Both 2,6-diisopropylcyclohexanol and 2,6-dimethylcyclohexanol were effective as anaesthetics with EC(50)s of 14.0 μM and 13.1 μM respectively, while other cyclohexanols with bulkier side chains were less potent. In conclusion, our data indicate that cyclohexanols are both positive modulators of GABA(A) receptors currents and anaesthetics. The positioning and size of the alkyl groups at the 2 and 6 positions on the cyclohexanol ring were critical determinants of activity.

Quantitative bioactivity prediction and pharmacophore identification for benzotriazine derivatives using the electron conformational-genetic algorithm in QSAR

The electron conformational-genetic algorithm (EC-GA), a sophisticated hybrid approach combining the GA and EC methods, has been employed for a 4D-QSAR procedure to identify the pharmacophore for benzotriazines as sarcoma inhibitors and for quantitative prediction of activity. The calculated geometry and electronic structure parameters of every atom and bond of each molecule are arranged in a matrix described as the electron-conformational matrix of contiguity (ECMC). By comparing the ECMC of one of the most active compounds with other ECMCs we were able to obtain the features of the pharmacophore responsible for the activity, as submatrices of the template known as electron conformational submatrices of activity. The GA was used to select the most important descriptors and to predict the theoretical activity of training and test sets. The predictivity of the model was internally validated. The best QSAR model was selected, having r² = 0.9008, standard error = 0.0510 and cross-validated squared correlation coefficient, q² = 0.8192.

Potentiating action of propofol at GABAA receptors of retinal bipolar cells

PURPOSE

Propofol (2,6-diisopropyl phenol), a widely used systemic anesthetic, is known to potentiate GABA(A) receptor activity in a number of CNS neurons and to produce changes in electroretinographically recorded responses of the retina. However, little is known about propofol's effects on specific retinal neurons. The authors investigated the action of propofol on GABA-elicited membrane current responses of retinal bipolar cells, which have both GABA(A) and GABA(C) receptors.

METHODS

Single, enzymatically dissociated bipolar cells obtained from rat retina were treated with propofol delivered by brief application in combination with GABA or other pharmacologic agents or as a component of the superfusing medium.

RESULTS

When applied with GABA at subsaturating concentrations and with TPMPA (a known GABA(C) antagonist), propofol markedly increased the peak amplitude and altered the kinetics of the response. Propofol increased the response elicited by THIP (a GABA(A)-selective agonist), and the response was reduced by bicuculline (a GABA(A) antagonist). The response to 5-methyl I4AA, a GABA(C)-selective agonist, was not enhanced by propofol. Serial brief applications of (GABA + TPMPA + propofol) led to a progressive increase in peak response amplitude and, at higher propofol concentrations, additional changes that included a prolonged time course of response recovery. Pre-exposure of the cell to perfusing propofol typically enhanced the rate of development of potentiation produced by (GABA + TPMPA + propofol) applications.

CONCLUSIONS

Propofol exerts a marked and selective potentiation on GABA(A) receptors of retinal bipolar cells. The data encourage the use of propofol in future studies of bipolar cell function.

Spike avalanches exhibit universal dynamics across the sleep-wake cycle

BACKGROUND

Scale-invariant neuronal avalanches have been observed in cell cultures and slices as well as anesthetized and awake brains, suggesting that the brain operates near criticality, i.e. within a narrow margin between avalanche propagation and extinction. In theory, criticality provides many desirable features for the behaving brain, optimizing computational capabilities, information transmission, sensitivity to sensory stimuli and size of memory repertoires. However, a thorough characterization of neuronal avalanches in freely-behaving (FB) animals is still missing, thus raising doubts about their relevance for brain function.

METHODOLOGY/PRINCIPAL FINDINGS

To address this issue, we employed chronically implanted multielectrode arrays (MEA) to record avalanches of action potentials (spikes) from the cerebral cortex and hippocampus of 14 rats, as they spontaneously traversed the wake-sleep cycle, explored novel objects or were subjected to anesthesia (AN). We then modeled spike avalanches to evaluate the impact of sparse MEA sampling on their statistics. We found that the size distribution of spike avalanches are well fit by lognormal distributions in FB animals, and by truncated power laws in the AN group. FB data surrogation markedly decreases the tail of the distribution, i.e. spike shuffling destroys the largest avalanches. The FB data are also characterized by multiple key features compatible with criticality in the temporal domain, such as 1/f spectra and long-term correlations as measured by detrended fluctuation analysis. These signatures are very stable across waking, slow-wave sleep and rapid-eye-movement sleep, but collapse during anesthesia. Likewise, waiting time distributions obey a single scaling function during all natural behavioral states, but not during anesthesia. Results are equivalent for neuronal ensembles recorded from visual and tactile areas of the cerebral cortex, as well as the hippocampus.

CONCLUSIONS/SIGNIFICANCE

Altogether, the data provide a comprehensive link between behavior and brain criticality, revealing a unique scale-invariant regime of spike avalanches across all major behaviors.

4D-QSAR: perspectives in drug design

Drug design is a process driven by innovation and technological breakthroughs involving a combination of advanced experimental and computational methods. A broad variety of medicinal chemistry approaches can be used for the identification of hits, generation of leads, as well as to accelerate the optimization of leads into drug candidates. The quantitative structure–activity relationship (QSAR) formalisms are among the most important strategies that can be applied for the successful design new molecules. This review provides a comprehensive review on the evolution and current status of 4D-QSAR, highlighting present challenges and new opportunities in drug design.

3D pharmacophore mapping using 4D QSAR analysis for the cytotoxicity of lamellarins against human hormone-dependent T47D breast cancer cells

4D quantitative structure-activity relationship (QSAR) and 3D pharmacophore models were built and investigated for cytotoxicity using a training set of 25 lamellarins against human hormone dependent T47D breast cancer cells. Receptor-independent (RI) 4D QSAR models were first constructed from the exploration of eight possible receptor-binding alignments for the entire training set. Since the training set is small (25 compounds), the generality of the 4D QSAR paradigm was then exploited to devise a strategy to maximize the extraction of binding information from the training set and to also permit virtual screening of diverse lamellarin chemistry. 4D QSAR models were sought for only six of the most potent lamellarins of the training set as well as another subset composed of lamellarins with constrained ranges in molecular weight and lipophilicity. This overall modeling strategy has permitted maximizing 3D pharmacophore information from this small set of structurally complex lamellarins that can be used to drive future analog synthesis and the selection of alternate scaffolds. Overall, it was found that the formation of an intermolecular hydrogen bond and the hydrophobic interactions for substituents on the E ring most modulate the cytotoxicity against T47D breast cancer cells. Hydrophobic substitutions on the F-ring can also enhance cytotoxic potency. A complementary high-throughput virtual screen to the 3D pharmacophore models, a 4D fingerprint QSAR model, was constructed using absolute molecular similarity. This 4D fingerprint virtual high-throughput screen permits a larger range of chemistry diversity to be assayed than with the 4D QSAR models. The optimized 4D QSAR 3D pharmacophore model has a leave-one-out cross-correlation value of xv-r2 = 0.947, while the optimized 4D fingerprint virtual screening model has a value of xv-r2 = 0.719. This work reveals that it is possible to develop significant QSAR, 3D pharmacophore, and virtual screening models for a small set of lamellarins showing cytotoxic behavior in breast cancer screens that can guide future drug development based upon lamellarin chemistry.

Comparative molecular field analysis and synthetic validation of a hydroxyamide-propofol binding and functional block of neuronal voltage-dependent sodium channels

Voltage gated sodium channels represent an important therapeutic target for a number of neurological disorders including epilepsy. Unfortunately, medicinal chemistry strategies for discovering new classes of antagonist for trans-membrane ion channels have been limited to mostly broad screening compound arrays. We have developed new sodium channel antagonist based on a propofol scaffold using the ligand based strategy of comparative molecular field analysis (CoMFA). The resulting CoMFA model was correlated and validated to provide insights into the design of new antagonists and to prioritize synthesis of these new structural analogs (compounds 4 and 5) that satisfied the steric and electrostatic model. Compounds 4 and 5 were evaluated for [(3)H]-batrachotoxinin-A-20-alpha-benzoate ([(3)H]-BTX-B) displacement yielding IC(50)'s of 22 and 5.7 microM, respectively. We further examined the potency of these two compounds to inhibit neuronal sodium currents recorded from cultured hippocampal neurons. At a concentration of 50 microM, compounds 4 and 5 tonically inhibited sodium channels currents by 59+/-7.8% (n=5) and 70+/-7.5% (n=7), respectively. This clearly demonstrates that these compounds functionally antagonize native neuronal sodium channel currents. In summary, we have shown that CoMFA can be effectively used to discover new classes of sodium channel antagonists.

A unitary anesthetic binding site at high resolution

Propofol is the most widely used injectable general anesthetic. Its targets include ligand-gated ion channels such as the GABA(A) receptor, but such receptor-channel complexes remain challenging to study at atomic resolution. Until structural biology methods advance to the point of being able to deal with systems such as the GABA(A) receptor, it will be necessary to use more tractable surrogates to probe the molecular details of anesthetic recognition. We have previously shown that recognition of inhalational general anesthetics by the model protein apoferritin closely mirrors recognition by more complex and clinically relevant protein targets; here we show that apoferritin also binds propofol and related GABAergic anesthetics, and that the same binding site mediates recognition of both inhalational and injectable anesthetics. Apoferritin binding affinities for a series of propofol analogs were found to be strongly correlated with the ability to potentiate GABA responses at GABA(A) receptors, validating this model system for injectable anesthetics. High resolution x-ray crystal structures reveal that, despite the presence of hydrogen bond donors and acceptors, anesthetic recognition is mediated largely by van der Waals forces and the hydrophobic effect. Molecular dynamics simulations indicate that the ligands undergo considerable fluctuations about their equilibrium positions. Finally, apoferritin displays both structural and dynamic responses to anesthetic binding, which may mimic changes elicited by anesthetics in physiologic targets like ion channels.

Modeling the LPS neutralization activity of anti-endotoxins

Bacterial lipopolysaccharides (LPS), also known as endotoxins, are major structural components of the outer membrane of Gram-negative bacteria that serve as a barrier and protective shield between them and their surrounding environment. LPS is considered to be a major virulence factor as it strongly stimulates the secretion of pro-inflammatory cytokines which mediate the host immune response and culminating in septic shock. Quantitative structure-activity relationship studies of the LPS neutralization activities of anti-endotoxins were performed using charge and quantum chemical descriptors. Artificial neural network implementing the back-propagation algorithm was selected for the multivariate analysis. The predicted activities from leave-one-out cross-validation were well correlated with the experimental values as observed from the correlation coefficient and root mean square error of 0.930 and 0.162, respectively. Similarly, the external testing set also yielded good predictivity with correlation coefficient and root mean square error of 0.983 and 0.130. The model holds great potential for the rational design of novel and robust compounds with enhanced neutralization activity.

What makes a molecule an anaesthetic? Studies on the mechanisms of anaesthesia using a physicochemical approach

Recent studies of mechanisms of anaesthesia have been mainly 'target orientated', investigating the activity of both volatile and i.v. agents at putative sites of action. An alternative approach is one that is 'ligand orientated', focusing on the properties of molecules that define their immobilizing ability and secondly define their potency. The use of conventional descriptors (such as non-polar solubility or the octanol-water partition coefficient [Log P]) are limited in their utility as predictors of potency as they represent three-dimensional molecular properties as a one-dimensional parameter. Using different computer-based molecular modelling methods (molecular similarity studies and comparative molecular field analysis [CoMFA]), we have identified the molecular bases of the activity of structurally diverse anaesthetics, such that they can be described as a single model based on the spatial distribution of molecular bulk and electrostatic potential. The same approach can also be used to model other properties of anaesthetic agents, such as cardiovascular depression. The present data suggest that, for the i.v. agents, it may be difficult to separate immobilizing (anaesthetic) activity and cardiovascular depression within a single molecule.

Tryptophan mutations at azi-etomidate photo-incorporation sites on alpha1 or beta2 subunits enhance GABAA receptor gating and reduce etomidate modulation

The potent general anesthetic etomidate produces its effects by enhancing GABA(A) receptor activation. Its photolabel analog [(3)H]azi-etomidate labels residues within transmembrane domains on alpha and beta subunits: alphaMet236 and betaMet286. We hypothesized that these methionines contribute to etomidate sites formed at alpha-beta subunit interfaces and that increasing side-chain bulk and hydrophobicity at either locus would mimic etomidate binding and block etomidate effects. Channel activity was electrophysiologically quantified in alpha(1)beta(2)gamma(2L) receptors with alpha(1)M236W or beta(2)M286W mutations, in both the absence and the presence of etomidate. Measurements included spontaneous activation, GABA EC(50), etomidate agonist potentiation, etomidate direct activation, and rapid macrocurrent kinetics. Both alpha(1)M236W and beta(2)M286W mutations induced spontaneous channel opening, lowered GABA EC(50), increased maximal GABA efficacy, and slowed current deactivation, mimicking effects of etomidate on alpha(1)beta(2)gamma(2L) channels. These changes were larger with alpha(1)M236W than with beta(2)M286W. Etomidate (3.2 muM) reduced GABA EC(50) much less in alpha(1)M236Wbeta(2)gamma(2L) receptors (2-fold) than in wild type (23-fold). However, etomidate was more potent and efficacious in directly activating alpha(1)M236Wbeta(2)gamma(2L) compared with wild type. In alpha(1)beta(2)M286Wgamma(2L) receptors, etomidate induced neither agonist-potentiation nor direct channel activation. These results support the hypothesis that alpha(1)Met236 and beta(2)Met286 are within etomidate sites that allosterically link to channel gating. Although alpha(1)M236W produced the larger impact on channel gating, beta(2)M286W produced more profound changes in etomidate sensitivity, suggesting a dominant role in drug binding. Furthermore, quantitative mechanistic analysis demonstrated that wild-type and mutant results are consistent with the presence of only one class of etomidate sites mediating both agonist potentiation and direct activation.

Construction of 4D-QSAR models for use in the design of novel p38-MAPK inhibitors

The p38-mitogen-activated protein kinase (p38-MAPK) plays a key role in lipopolysaccharide-induced tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1) release during the inflammatory process, emerging as an attractive target for new anti-inflammatory agents. Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis [Hopfinger et al., J. Am. Chem. Soc., 119 (1997) 10509] was applied to a series of 33 (a training set of 28 and a test set of 5) pyridinyl-imidazole and pyrimidinyl-imidazole inhibitors of p38-MAPK, with IC50 ranging from 0.11 to 2100 nM [Liverton et al., J. Med. Chem., 42 (1999) 2180]. Five thousand conformations of each analogue were sampled from a molecular dynamics simulation (MDS) during 50 ps at a constant temperature of 303 K. Each conformation was placed in a 2 angstroms grid cell lattice for each of three trial alignments. 4D-QSAR models were constructed by genetic algorithm (GA) optimization and partial least squares (PLS) fitting, and evaluated by leave-one-out cross-validation technique. In the best models, with three to six terms, the adjusted cross-validated squared correlation coefficients, Q2adj, ranged from 0.67 to 0.85. Model D (Q2adj = 0.84) was identified as the most robust model from alignment 1, and it is representative of the other best models. This model encompasses new molecular regions as containing pharmacophore sites, such as the amino-benzyl moiety of pyrimidine analogs and the N1-substituent in the imidazole ring. These regions of the ligands should be further explored to identify better anti-inflammatory inhibitors of p38-MAPK.

Computational chemistry approaches to drug discovery in signal transduction

The advent of therapeutic strategies aimed at targeting specific macromolecular components of deregulated signaling pathways associated with particular disease states has given rise to the idea that it should be possible to design ligands as drug candidates to these targets from first principles. This concept has been beckoning for a long time but structure-based ligand design only became feasible once it was possible to determine the 3-D structures of molecular targets at atomic resolution. However, structure-based design turned out to be difficult, chiefly because under physiological conditions both receptors and ligands are not static but they behave dynamically. While it is possible to design ligands with high steric and electronic complementarity to a receptor site, it is always uncertain how biologically relevant the assumed conformations of both ligand and receptor actually are. The fact that it remains beyond our current abilities to predict with sufficient accuracy the affinity between hypothetical ligand and receptor poses is in part connected with this problem and continues to confound the reliable prediction of drug-like ligands for therapeutic targets. Nevertheless, significant progress has been made and so-called virtual screening methods that use computational methods to dock candidate ligands into receptor sites and to score the resulting complexes are now used routinely as one of the components in drug discovery screening campaigns. Here an overview is given of the underlying principles, implementations, and applications of structure-guided computational design technologies. Although the emphasis is on receptor-based strategies, mention will also be made of some of the more established ligand-based approaches, such as similarity analyses and quantitative structure-activity relationship methods.