Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations

Insights into the Nature of Anesthetic-Protein Interactions: An ONIOM Study

Anesthetics have been employed widely to relieve surgical suffering, but their mechanism of action is not yet clear. For over a century, the mechanism of anesthesia was previously thought to be via lipid bilayer interactions. In the present work, a rigorous three-layer ONIOM(M06-2X/6-31+G*:PM6:AMBER) method was utilized to investigate the nature of interactions between several anesthetics and actual protein binding sites. According to the calculated structural features, interaction energies, atomic charges, and electrostatic potential surfaces, the amphiphilic nature of anesthetic-protein interactions was demonstrated for both inhalational and injectable anesthetics. The existence of hydrogen and halogen bonding interactions between anesthetics and proteins was clearly identified, and these interactions served to assist ligand recognition and binding by the protein. Within all complexes of inhalational or injectable anesthetics, the polarization effects play a dominant role over the steric effects and induce a significant asymmetry in the otherwise symmetric atomic charge distributions of the free ligands in vacuo. This study provides new insight into the mechanism of action of general anesthetics in a more rigorous way than previously described. Future rational design of safer anesthetics for an aging and more physiologically vulnerable population will be predicated on this greater understanding of such specific interactions.

Sevoflurane binds and allosterically blocks integrin lymphocyte function-associated antigen-1

BACKGROUND

Volatile anesthetics have been shown to modify immune cell functions via several mechanisms, some of which have been only partially elucidated. We demonstrated that isoflurane inhibits primary leukocyte integrin lymphocyte function-associated antigen-1 (LFA-1) by binding to the allosteric cavity critical for conformational activation to its high-affinity form. It remains to be determined whether the allosteric inhibition of LFA-1 by isoflurane can be generalized to other anesthetics such as sevoflurane.

METHODS

The effects of sevoflurane on the ability of LFA-1 to bind to its counter-ligand, intercellular adhesion molecule-1, was studied in leukocytes by flow cytometry. To examine whether sevoflurane acts directly on LFA-1, we measured ligand-binding using beads coated with purified LFA-1 protein. To distinguish between competitive versus allosteric inhibition, we analyzed the effects of sevoflurane on both wild-type and mutant-locked high-affinity LFA-1. One-way analysis of variance was employed for statistical analysis of the data. Nuclear magnetic resonance spectroscopy was used to identify sevoflurane binding site(s).

RESULTS

Sevoflurane at clinically relevant concentrations inhibited the ligand-binding function of LFA-1 in leukocytes as well as in cell-free assays (P<0.05). Sevoflurane blocked wild-type but not locked high-affinity LFA-1, thereby demonstrating an allosteric mode of inhibition. Nuclear magnetic resonance spectroscopy revealed that sevoflurane bound to the allosteric cavity, to which LFA-1 allosteric antagonists and isoflurane also bind.

CONCLUSIONS

This study suggests that sevoflurane also blocks the activation-dependent conformational changes of LFA-1 to the high-affinity form. The allosteric mode of action exemplified by sevoflurane and isoflurane via LFA-1 might represent one of the underlying mechanisms of anesthetic-mediated immunomodulation.

The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR

Water is an essential component for many biological processes. Pauling proposed that water might play a critical role in general anesthesia by forming water clathrates around anesthetic molecules. To examine potential involvement of water in general anesthesia, we analyzed water within alpha4beta2 nAChR, a putative protein target hypersensitive to volatile anesthetics. Experimental structure-derived closed- and open-channel nAChR systems in a fully hydrated lipid bilayer were examined using all-atom molecular dynamics simulations. At the majority of binding sites in alpha4beta2 nAChR, halothane replaced the slow-exchanging water molecules and caused a regional water population decrease. Only two binding sites had an increased quantity of water in the presence of halothane, where water arrangements resemble clathrate-like structures. The small number of such clathrate-like water clusters suggests that the formation of water clathrates is unlikely to be a primary cause for anesthesia. Despite the decrease in water population at most of the halothane binding sites, the number of sites that can be occupied transiently by water is actually increased in the presence of halothane. Many of these water sites were located between two subunits or in regions containing agonist binding sites or critical structural elements for transducing agonist binding to channel gating. Changes in water sites in the presence of halothane affected water-mediated protein-protein interactions and the protein dynamics, which can have direct impact on protein function. Collectively, water contributes to the action of anesthetics in proteins by mediating interactions between protein subunits and altering protein dynamics, instead of forming water clathrates around anesthetics.

Unresponsive correlated motion in alpha7 nAChR to halothane binding explains its functional insensitivity to volatile anesthetics

Neuronal nicotinic acetylcholine receptors (nAChRs) have been implicated as targets for general anesthetics, but the functional responses to anesthetic modulation vary considerably among different subtypes of nAChRs. Inhaled general anesthetics, such as halothane, could effectively inhibit the channel activity of the alpha4beta2 nAChR but not the homologous alpha7 nAChR. To understand why alpha7 is insensitive to inhaled general anesthetics, we performed multiple sets of 20 ns molecular dynamics (MD) simulations on the closed- and open-channel alpha7 in the absence and presence of halothane and critically compared the results with those from our studies on the alpha4beta2 nAChR (Liu et al. J. Phys. Chem. B 2009, 113, 12581 and Liu et al. J. Phys. Chem. B 2010, 114, 626). Several halothane binding sites with fairly high binding affinities were identified in both closed- and open-channel alpha7, suggesting that a lack of sensitive functional responses of the alpha7 nAChR to halothane in the previous experiments was unlikely due to a lack of halothane interaction with alpha7. The binding affinities of halothane in alpha7 seemed to be protein conformation-dependent. Overall, halothane affinity was higher in the closed-channel alpha7. Halothane binding to alpha7 did not induce profound changes in alpha7 structure and dynamics that could be related to the channel function. In contrast, correlated motion of the open-channel alpha4beta2 was reduced substantially in the presence of halothane, primarily due to the more susceptible nature of beta2 to anesthetic modulation. The amphiphilic extracellular and transmembrane domain interface of the beta2 subunit is attractive to halothane and is susceptible to halothane perturbation, which may be why alpha4beta2 is functionally more sensitive to halothane than alpha7.

Computational studies on the interactions of inhalational anesthetics with proteins

Despite the widespread clinical use of anesthetics since the 19th century, a clear understanding of the mechanism of anesthetic action has yet to emerge. On the basis of early experiments by Meyer, Overton, and subsequent researchers, the cell's lipid membrane was generally concluded to be the primary site of action of anesthetics. However, later experiments with lipid-free globular proteins, such as luciferase and apoferritin, shifted the focus of anesthetic action to proteins. Recent experimental studies, such as photoaffinity labeling and mutagenesis on membrane proteins, have suggested specific binding sites for anesthetic molecules, further strengthening the proteocentric view of anesthetic mechanism. With the increased availability of high-resolution crystal structures of ion channels and other integral membrane proteins, as well as the availability of powerful computers, the structure-function relationship of anesthetic-protein interactions can now be investigated in atomic detail. In this Account, we review recent experiments and related computer simulation studies involving interactions of inhalational anesthetics and proteins, with a particular focus on membrane proteins. Globular proteins have long been used as models for understanding the role of protein-anesthetic interactions and are accordingly examined in this Account. Using selected examples of membrane proteins, such as nicotinic acetyl choline receptor (nAChR) and potassium channels, we address the issues of anesthetic binding pockets in proteins, the role of conformation in anesthetic effects, and the modulation of local as well as global dynamics of proteins by inhaled anesthetics. In the case of nicotinic receptors, inhalational anesthetic halothane binds to the hydrophobic cavity close to the M2-M3 loop. This binding modulates the dynamics of the M2-M3 loop, which is implicated in allosterically transmitting the effects to the channel gate, thus altering the function of the protein. In potassium channels, anesthetic molecules preferentially potentiate the open conformation by quenching the motion of the aromatic residues implicated in the gating of the channel. These simulations suggest that low-affinity drugs (such as inhalational anesthetics) modulate the protein function by influencing local as well as global dynamics of proteins. Because of intrinsic experimental limitations, computational approaches represent an important avenue for exploring the mode of action of anesthetics. Molecular dynamics simulations-a computational technique frequently used in the general study of proteins-offer particular insight in the study of the interaction of inhalational anesthetics with membrane proteins.

General anesthetic binding to neuronal alpha4beta2 nicotinic acetylcholine receptor and its effects on global dynamics

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a target for general anesthetics. Currently available experimental structural information is inadequate to understand where anesthetics bind and how they modulate the receptor motions essential to function. Using our published open-channel structure model of alpha4beta2 nAChR, we identified and evaluated six amphiphilic interaction sites for the volatile anesthetic halothane via flexible ligand docking and subsequent 20-ns molecular dynamics simulations. Halothane binding energies ranged from -6.8 to -2.4 kcal/mol. The primary binding sites were located at the interface of extracellular and transmembrane domains, where halothane perturbed conformations of, and widened the gap among, the Cys loop, the beta1-beta2 loop, and the TM2-TM3 linker. The halothane with the highest binding affinity at the interface between the alpha4 and beta2 subunits altered interactions between the protein and nearby lipids by competing for hydrogen bonds. Gaussian network model analyses of the alpha4beta2 nAChR structures at the end of 20-ns simulations in the absence or presence of halothane revealed profound changes in protein residue mobility. The concerted motions critical to protein function were also perturbed considerably. Halothane's effect on protein dynamics was not confined to the residues adjacent to the binding sites, indicating an action on a more global scale.

Volatile anesthetic binding to proteins is influenced by solvent and aliphatic residues

The main objective of this work was to characterize VA binding sites in multiple anesthetic target proteins. A computational algorithm was used to quantify the solvent exclusion and aliphatic character of amphiphilic pockets in the structures of VA binding proteins. VA binding sites in the protein structures were defined as the pockets with solvent exclusion and aliphatic character that exceeded minimum values observed in the VA binding sites of serum albumin, firefly luciferase, and apoferritin. We found that the structures of VA binding proteins are enriched in these pockets and that the predicted binding sites were consistent with experimental determined binding locations in several proteins. Autodock3 was used to dock the simulated molecules of 1,1,1,2,2-pentafluoroethane, difluoromethyl 1,1,1,2-tetrafluoroethyl ether, and sevoflurane and the isomers of halothane and isoflurane into these potential binding sites. We found that the binding of the various VA molecules to the amphiphilic pockets is driven primarily by VDW interactions and to a lesser extent by weak hydrogen bonding and electrostatic interactions. In addition, the trend in Delta G binding values follows the Meyer-Overton rule. These results suggest that VA potencies are related to the VDW interactions between the VA ligand and protein target. It is likely that VA bind to sites with a high degree of solvent exclusion and aliphatic character because aliphatic residues provide favorable VDW contacts and weak hydrogen bond donors. Water molecules occupying these sites maintain pocket integrity, associate with the VA ligand, and diminish the unfavorable solvation enthalpy of the VA. Water molecules displaced into the bulk by the VA ligand may provide an additional favorable enthalpic contribution to VA binding. Anesthesia is a component of many health related procedures, the outcomes of which could be improved with a better understanding of the molecular targets and mechanisms of anesthetic action.

The volatile anesthetic isoflurane perturbs conformational activation of integrin LFA-1 by binding to the allosteric regulatory cavity

The molecular and structural basis of anesthetic interactions with conformations and functionalities of cell surface receptors remains to be elucidated. We have demonstrated that the widely used volatile anesthetic isoflurane blocks the activation-dependent conformational conversion of integrin lymphocyte function associated antigen-1 (LFA-1), the major leukocyte cell adhesion molecule, to a high-affinity configuration. Perturbation of LFA-1 activation by isoflurane at clinically relevant concentrations leads to the inhibition of T-cell interactions with target cells as well as ligand-triggered intracellular signaling. Nuclear magnetic resonance spectroscopy reveals that isoflurane binds within a cavity in the LFA-1 ligand-binding domain, which is a previously identified drug-binding site for allosteric small-molecule antagonists that stabilize LFA-1 in a low-affinity conformation. These results provide a potential mechanism for the immunomodulatory properties of isoflurane.

Four-alpha-helix bundle with designed anesthetic binding pockets. Part I: structural and dynamical analyses

The four-alpha-helix bundle mimics the transmembrane domain of the Cys-loop receptor family believed to be the protein target for general anesthetics. Using high resolution NMR, we solved the structure (Protein Data Bank ID: 2I7U) of a prototypical dimeric four-alpha-helix bundle, (Aalpha(2)-L1M/L38M)(2,) with designed specific binding pockets for volatile anesthetics. Two monomers of the helix-turn-helix motif form an antiparallel dimer as originally designed, but the high-resolution structure exhibits an asymmetric quaternary arrangement of the four helices. The two helices from the N-terminus to the linker (helices 1 and 1') are associated with each other in the dimer by the side-chain ring stacking of F12 and W15 along the long hydrophobic core and by a nearly perfect stretch of hydrophobic interactions between the complementary pairs of L4, L11, L18, and L25, all of which are located at the heptad e position along the helix-helix dimer interface. In comparison, the axes of the two helices from the linker to the C-terminus (helices 2 and 2') are wider apart from each other, creating a lateral access pathway around K47 from the aqueous phase to the center of the designed hydrophobic core. The site of the L38M mutation, which was previously shown to increase the halothane binding affinity by approximately 3.5-fold, is not part of the hydrophobic core presumably involved in the anesthetic binding but shows an elevated transverse relaxation (R(2)) rate. Qualitative analysis of the protein dynamics by reduced spectral density mapping revealed exchange contributions to the relaxation at many residues in the helices. This observation was confirmed by the quantitative analysis using the Modelfree approach and by the NMR relaxation dispersion measurements. The NMR structures and Autodock analysis suggest that the pocket with the most favorable amphipathic property for anesthetic binding is located between the W15 side chains at the center of the dimeric hydrophobic core, with the possibility of two additional minor binding sites between the F12 and F52 ring stacks of each monomer. The high-resolution structure of the designed anesthetic-binding protein offers unprecedented atomistic details about possible sites for anesthetic-protein interactions that are essential to the understanding of molecular mechanisms of general anesthesia.

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.

Four-alpha-helix bundle with designed anesthetic binding pockets. Part II: halothane effects on structure and dynamics

As a model of the protein targets for volatile anesthetics, the dimeric four-alpha-helix bundle, (Aalpha(2)-L1M/L38M)(2), was designed to contain a long hydrophobic core, enclosed by four amphipathic alpha-helices, for specific anesthetic binding. The structural and dynamical analyses of (Aalpha(2)-L1M/L38M)(2) in the absence of anesthetics (another study) showed a highly dynamic antiparallel dimer with an asymmetric arrangement of the four helices and a lateral accessing pathway from the aqueous phase to the hydrophobic core. In this study, we determined the high-resolution NMR structure of (Aalpha(2)-L1M/L38M)(2) in the presence of halothane, a clinically used volatile anesthetic. The high-solution NMR structure, with a backbone root mean-square deviation of 1.72 A (2JST), and the NMR binding measurements revealed that the primary halothane binding site is located between two side-chains of W15 from each monomer, different from the initially designed anesthetic binding sites. Hydrophobic interactions with residues A44 and L18 also contribute to stabilizing the bound halothane. Whereas halothane produces minor changes in the monomer structure, the quaternary arrangement of the dimer is shifted by about half a helical turn and twists relative to each other, which leads to the closure of the lateral access pathway to the hydrophobic core. Quantitative dynamics analyses, including Modelfree analysis of the relaxation data and the Carr-Purcell-Meiboom-Gill transverse relaxation dispersion measurements, suggest that the most profound anesthetic effect is the suppression of the conformational exchange both near and remote from the binding site. Our results revealed a novel mechanism of an induced fit between anesthetic molecule and its protein target, with the direct consequence of protein dynamics changing on a global rather than a local scale. This mechanism may be universal to anesthetic action on neuronal proteins.

Dynamics of firefly luciferase inhibition by general anesthetics: Gaussian and anisotropic network analyses

The new crystal structures of the product-bound firefly luciferase combined with the previously determined substrate-free structures allow for a detailed analysis of the dynamics basis for the luciferase enzymatic activities. Using the Gaussian network model and the anisotropic network model, we show here that the superposition of the three slowest anisotropic network model modes, consisting of the bending, rotating, and rocking motions of the C-domain, accounts for large rearrangement of domains from the substrate-free (open) to product-bound (closed) conformation and thus constitutes a critical component of the enzyme's functions. The analysis also offers a unique platform to reexamine the molecular mechanism of the anesthetic inhibition of the firefly luciferase. Through perturbing the protein backbone network by introducing additional nodes to represent anesthetics, we found that the presence of two representative anesthetics, halothane and n-decanol, in different regions of luciferase had distinctively different effects on the protein's global motion. Only at the interface of the C- and N-domains did the anesthetics cause the most profound reduction in the overall flexibility of the C-domain and the concomitant increase in the flexibility of the loop, where the substitution of a conserved lysine residue was found experimentally to lead to >2-3 orders of magnitude reduction in activity. These anesthetic-induced dynamics changes can alter the normal function of the protein, appearing as an epiphenomenon of an "inhibition". The implication of the study is that a leading element for general anesthetic action on proteins is to disrupt the modes of motion essential to protein functions.

Prediction of volatile anesthetic binding sites in proteins

Computational methods designed to predict and visualize ligand protein binding interactions were used to characterize volatile anesthetic (VA) binding sites and unoccupied pockets within the known structures of VAs bound to serum albumin, luciferase, and apoferritin. We found that both the number of protein atoms and methyl hydrogen, which are within approximately 8 A of a potential ligand binding site, are significantly greater in protein pockets where VAs bind. This computational approach was applied to structures of calmodulin (CaM), which have not been determined in complex with a VA. It predicted that VAs bind to [Ca(2+)](4)-CaM, but not to apo-CaM, which we confirmed with isothermal titration calorimetry. The VA binding sites predicted for the structures of [Ca(2+)](4)-CaM are located in hydrophobic pockets that form when the Ca(2+) binding sites in CaM are saturated. The binding of VAs to these hydrophobic pockets is supported by evidence that halothane predominantly makes contact with aliphatic resonances in [Ca(2+)](4)-CaM (nuclear Overhauser effect) and increases the Ca(2+) affinity of CaM (fluorescence spectroscopy). Our computational analysis and experiments indicate that binding of VA to proteins is consistent with the hydrophobic effect and the Meyer-Overton rule.

Blind docking of drug-sized compounds to proteins with up to a thousand residues

Blind docking was introduced for the detection of possible binding sites and modes of peptide ligands by scanning the entire surface of protein targets. In the present study, the method is tested on a group of drug-sized compounds and proteins with up to a thousand amino acid residues. Both proteins from complex structures and ligand-free proteins were used as targets. Robustness, limitations and future perspectives of the method are discussed. It is concluded that blind docking can be used for unbiased mapping of the binding patterns of drug candidates.