Size versus polarizability in protein-ligand interactions: binding of noble gases within engineered cavities in phage T4 lysozyme

Identifying the toxic mechanisms of emerging electronic contaminations liquid crystal monomers and the construction of a priority control list for graded control

Liquid crystal monomers (LCMs) are identified as emerging organic contaminations with largely unexplored health impacts. To elucidate their toxic mechanisms, support the establishment of environmental discharge and management standards, and promote effective LCMs control, this study constructs a database covering 20,545 potential targets of 1431 LCMs, highlighting 9 key toxic target proteins that disrupt the nervous system and metabolic functions. GO and KEGG pathway analysis suggests LCMs severely affect nervous system, linked to neurodegenerative diseases and mental health disorders, with toxicity variations driven by electronegativity and structural complexity of LCM terminal groups. To achieve tiered control of LCMs, construct toxicity risk control lists for 9 key toxic target proteins, suitable for the graded control of LCMs, management recommendations are provided based on toxicity levels. These lists were validated for reliability and offer reliable toxicity predictions for LCMs. SHAP analysis points to electronic properties, molecular shape, and structural characteristics of LCMs as primary health impact factors. As the first study integrating machine learning with computational toxicology to outline LCMs health impacts, it aims to enhance public understanding of LCM toxicity risks and support the development of environmental standards, effective management of LCM production and emissions, and reduction of public exposure risks.

The High-Pressure Freezing Laboratory for Macromolecular Crystallography (HPMX), an ancillary tool for the macromolecular crystallography beamlines at the ESRF

This article describes the High-Pressure Freezing Laboratory for Macromolecular Crystallography (HPMX) at the ESRF, and highlights new and complementary research opportunities that can be explored using this facility. The laboratory is dedicated to investigating interactions between macromolecules and gases in crystallo, and finds applications in many fields of research, including fundamental biology, biochemistry, and environmental and medical science. At present, the HPMX laboratory offers the use of different high-pressure cells adapted for helium, argon, krypton, xenon, nitrogen, oxygen, carbon dioxide and methane. Important scientific applications of high pressure to macromolecules at the HPMX include noble-gas derivatization of crystals to detect and map the internal architecture of proteins (pockets, tunnels and channels) that allows the storage and diffusion of ligands or substrates/products, the investigation of the catalytic mechanisms of gas-employing enzymes (using oxygen, carbon dioxide or methane as substrates) to possibly decipher intermediates, and studies of the conformational fluctuations or structure modifications that are necessary for proteins to function. Additionally, cryo-cooling protein crystals under high pressure (helium or argon at 2000 bar) enables the addition of cryo-protectant to be avoided and noble gases can be employed to produce derivatives for structure resolution. The high-pressure systems are designed to process crystals along a well defined pathway in the phase diagram (pressure-temperature) of the gas to cryo-cool the samples according to the three-step `soak-and-freeze method'. Firstly, crystals are soaked in a pressurized pure gas atmosphere (at 294 K) to introduce the gas and facilitate its interactions within the macromolecules. Samples are then flash-cooled (at 100 K) while still under pressure to cryo-trap macromolecule-gas complexation states or pressure-induced protein modifications. Finally, the samples are recovered after depressurization at cryo-temperatures. The final section of this publication presents a selection of different typical high-pressure experiments carried out at the HPMX, showing that this technique has already answered a wide range of scientific questions. It is shown that the use of different gases and pressure conditions can be used to probe various effects, such as mapping the functional internal architectures of enzymes (tunnels in the haloalkane dehalogenase DhaA) and allosteric sites on membrane-protein surfaces, the interaction of non-inert gases with proteins (oxygen in the hydrogenase ReMBH) and pressure-induced structural changes of proteins (tetramer dissociation in urate oxidase). The technique is versatile and the provision of pressure cells and their application at the HPMX is gradually being extended to address new scientific questions.

Local Xenon-Protein Interaction Produces Global Conformational Change and Allosteric Inhibition in Lysozyme

Noble gases have well-established biological effects, yet their molecular mechanisms remain poorly understood. Here, we investigated, both experimentally and computationally, the molecular modes of xenon (Xe) action in bacteriophage T4 lysozyme (T4L). By combining indirect gassing methods with a colorimetric lysozyme activity assay, a reversible, Xe-specific (20 ± 3)% inhibition effect was observed. Accelerated molecular dynamic simulations revealed that Xe exerts allosteric inhibition on the protein by expanding a C-terminal hydrophobic cavity. Xe-induced cavity expansion results in global conformational changes, with long-range transduction distorting the active site where peptidoglycan binds. Interestingly, the peptide substrate binding site that enables lysozyme specificity does not change conformation. Two T4L mutants designed to reshape the C-terminal Xe cavity established a correlation between cavity expansion and enzyme inhibition. This work also highlights the use of Xe flooding simulations to identify new cryptic binding pockets. These results enrich our understanding of Xe-protein interactions at the molecular level and inspire further biochemical investigations with noble gases.

Evolutionary Role of Water-Accessible Cavities in Src Homology 2 (SH2) Domains

Protein excited states are fundamental in the understanding of biological function, despite the fact they are hardly observed using traditional biophysical methodologies. Pressure perturbation coupled with nuclear magnetic resonance (NMR) spectroscopy is a powerful physicochemical tool to glance at these low-populated high-energy states on a residue-by-residue basis and underpin mechanistic insights into protein functionalities. Here we performed pressure titrations using NMR spectroscopy and relaxation dispersion experiments to identify the low-lying energetic states of the c-Abl SH2 domain. By showing that the SH2 excited state contains a hydrated hydrophobic cavity, fast-exchange motions, and highly conserved residues facing the water-accessible hole, we discuss the implications of water-protein interactions in SH2 modules achieving high-affinity binding and promiscuous phospho-Tyr peptide recognition.

Filling of a water-free void explains the allosteric regulation of the β1-adrenergic receptor by cholesterol

Recent high-pressure NMR results indicate that the preactive conformation of the β₁-adrenergic receptor (β₁AR) harbours completely empty cavities of ~100 ų volume, which disappear in the active conformation of the receptor. Here we have localized these cavities using X-ray crystallography of xenon-derivatized β₁AR crystals. One of the cavities is in direct contact with the cholesterol-binding pocket. Solution NMR shows that addition of the cholesterol analogue cholesteryl hemisuccinate impedes the formation of the active conformation of detergent-solubilized β₁AR by blocking conserved G protein-coupled receptor microswitches, concomitant with an affinity reduction of both isoprenaline and G protein-mimicking nanobody Nb80 for β₁AR detected by isothermal titration calorimetry. This wedge-like action explains the function of cholesterol as a negative allosteric modulator of β₁AR. A detailed understanding of G protein-coupled receptor regulation by cholesterol by filling of a dry void and the easy scouting for such voids by xenon may provide new routes for the development of allosteric drugs.

High-pressure crystallography shows noble gas intervention into protein-lipid interaction and suggests a model for anaesthetic action

In this work we examine how small hydrophobic molecules such as inert gases interact with membrane proteins (MPs) at a molecular level. High pressure atmospheres of argon and krypton were used to produce noble gas derivatives of crystals of three well studied MPs (two different proton pumps and a sodium light-driven ion pump). The structures obtained using X-ray crystallography showed that the vast majority of argon and krypton binding sites were located on the outer hydrophobic surface of the MPs – a surface usually accommodating hydrophobic chains of annular lipids (which are known structural and functional determinants for MPs). In conformity with these results, supplementary in silico molecular dynamics (MD) analysis predicted even greater numbers of argon and krypton binding positions on MP surface within the bilayer. These results indicate a potential importance of such interactions, particularly as related to the phenomenon of noble gas-induced anaesthesia.

Noble gases are known to interact with proteins and can be good anaesthetics in hyperbaric conditions. This study identifies argon and krypton binding sites on membrane proteins and proposes as a hypothesis that noble gases, by altering protein/lipid contacts, may affect protein function.

Discovery of the Xenon-Protein Interactome Using Large-Scale Measurements of Protein Folding and Stability

The intermolecular interactions of noble gases in biological systems are associated with numerous biochemical responses, including apoptosis, inflammation, anesthesia, analgesia, and neuroprotection. The molecular modes of action underlying these responses are largely unknown. This is in large part due to the limited experimental techniques to study protein-gas interactions. The few techniques that are amenable to such studies are relatively low-throughput and require large amounts of purified proteins. Thus, they do not enable the large-scale analyses that are useful for protein target discovery. Here, we report the application of stability of proteins from rates of oxidation (SPROX) and limited proteolysis (LiP) methodologies to detect protein-xenon interactions on the proteomic scale using protein folding stability measurements. Over 5000 methionine-containing peptides and over 5000 semi-tryptic peptides, mapping to ∼1500 and ∼950 proteins, respectively, in the yeast proteome, were assayed for Xe-interacting activity using the SPROX and LiP techniques. The SPROX and LiP analyses identified 31 and 60 Xe-interacting proteins, respectively, none of which were previously known to bind Xe. A bioinformatics analysis of the proteomic results revealed that these Xe-interacting proteins were enriched in those involved in ATP-driven processes. A fraction of the protein targets that were identified are tied to previously established modes of action related to xenon's anesthetic and organoprotective properties. These results enrich our knowledge and understanding of biologically relevant xenon interactions. The sample preparation protocols and analytical methodologies developed here for xenon are also generally applicable to the discovery of a wide range of other protein-gas interactions in complex biological mixtures, such as cell lysates.

Method for the Identification of Potentially Bioactive Argon Binding Sites in Protein Families

Argon belongs to the group of chemically inert noble gases, which display a remarkable spectrum of clinically useful biological properties. In an attempt to better understand noble gases, notably argon's mechanism of action, we mined a massive noble gas modeling database which lists all possible noble gas binding sites in the proteins from the Protein Data Bank. We developed a method of analysis to identify among all predicted noble gas binding sites the potentially relevant ones within protein families which are likely to be modulated by Ar. Our method consists in determining within structurally aligned proteins the conserved binding sites whose shape, localization, hydrophobicity, and binding energies are to be further examined. This method was applied to the analysis of two protein families where crystallographic noble gas binding sites have been experimentally determined. Our findings indicate that among the most conserved binding sites, either the most hydrophobic one and/or the site which has the best binding energy corresponds to the crystallographic noble gas binding sites with the best occupancies, therefore the best affinity for the gas. This method will allow us to predict relevant noble gas binding sites that have potential pharmacological interest and thus potential Ar targets that will be prioritized for further studies including in vitro validation.

Comparative study of the effects of high hydrostatic pressure per se and high argon pressure on urate oxidase ligand stabilization

The stability of the tetrameric enzyme urate oxidase in complex with excess of 8-azaxanthine was investigated either under high hydrostatic pressure per se or under a high pressure of argon. The active site is located at the interface of two subunits, and the catalytic activity is directly related to the integrity of the tetramer. This study demonstrates that applying pressure to a protein–ligand complex drives the thermodynamic equilibrium towards ligand saturation of the complex, revealing a new binding site. A transient dimeric intermediate that occurs during the pressure-induced dissociation process was characterized under argon pressure and excited substates of the enzyme that occur during the catalytic cycle can be trapped by pressure. Comparison of the different structures under pressure infers an allosteric role of the internal hydrophobic cavity in which argon is bound, since this cavity provides the necessary flexibility for the active site to function.

Neuroprotective Effects of the Inert Gas Argon on Experimental Traumatic Brain Injury In Vivo with the Controlled Cortical Impact Model in Mice

Simple Summary

Traumatic brain injuries remain one of the leading causes of death in the western world and developing countries. There is an urgent need for causal therapies for such injuries. The noble gas argon has already shown promising results in in-vitro models. The influence of argon on the extent of damage after a craniocerebral trauma will be investigated in this study, in vivo, in mice. After the trauma, the animals were examined for neurological impairments and their brains were removed to detect brain edema and microscopically detectable alterations.

Abstract

Argon has shown neuroprotective effects after traumatic brain injury (TBI) and cerebral ischemia in vitro and in focal cerebral ischemia in vivo. The purpose of this study is to show whether argon beneficially impacts brain contusion volume (BCV) as the primary outcome parameter, as well as secondary outcome parameters, such as brain edema, intracranial pressure (ICP), neurological outcome, and cerebral blood flow (CBF) in an in-vivo model. Subjects were randomly assigned to either argon treatment or room air. After applying controlled cortical impact (CCI) onto the dura with 8 m/s (displacement 1 mm, impact duration 150 ms), treatment was administered by a recovery chamber with 25%, 50%, or 75% argon and the rest being oxygen for 4 h after trauma. Two control groups received room air for 15 min and 24 h, respectively. Neurological testing and ICP measurements were performed 24 h after trauma, and brains were removed to measure secondary brain damage. The primary outcome parameter, BCV, and the secondary outcome parameter, brain edema, were not significantly reduced by argon treatment at any concentration. There was a highly significant decrease in ICP at 50% argon (p = 0.001), and significant neurological improvement (beamwalk missteps) at 25% and 50% argon (p = 0.01; p = 0.049 respectively) compared to control.

Molecular Sensing with Host Systems for Hyperpolarized 129Xe

Hyperpolarized noble gases have been used early on in applications for sensitivity enhanced NMR. 129Xe has been explored for various applications because it can be used beyond the gas-driven examination of void spaces. Its solubility in aqueous solutions and its affinity for hydrophobic binding pockets allows "functionalization" through combination with host structures that bind one or multiple gas atoms. Moreover, the transient nature of gas binding in such hosts allows the combination with another signal enhancement technique, namely chemical exchange saturation transfer (CEST). Different systems have been investigated for implementing various types of so-called Xe biosensors where the gas binds to a targeted host to address molecular markers or to sense biophysical parameters. This review summarizes developments in biosensor design and synthesis for achieving molecular sensing with NMR at unprecedented sensitivity. Aspects regarding Xe exchange kinetics and chemical engineering of various classes of hosts for an efficient build-up of the CEST effect will also be discussed as well as the cavity design of host molecules to identify a pool of bound Xe. The concept is presented in the broader context of reporter design with insights from other modalities that are helpful for advancing the field of Xe biosensors.

Evolution of the androgen receptor: Perspectives from human health to dancing birds

Androgenic hormones orchestrate the development and activation of diverse reproductive phenotypes across vertebrates. Although extensive work investigates how selection for these traits modifies individual elements of this signaling system (e.g., hormone or androgen receptor [AR] levels), we know less about natural variation in the AR sequence across vertebrates. Our knowledge of AR sequence mutations is largely limited to work in human patients or cell-lines, providing a framework to contextualize single mutations at the expense of evolutionary timescale. Here we unite both perspectives in a review that explores the functional significance of AR on a domain-by-domain basis, using existing knowledge to highlight how and why each region might evolve. We then examine AR sequence variation on different timescales by examining sequence variation in clades originating in the Cambrian (vertebrates; >500 mya) and Cretaceous (birds; >65 mya). In each case, we characterize how the receptor has changed over time and discuss which regions are most likely to evolve in response to selection. Overall, domains that are required for androgenic signaling to function (e.g., DNA- and ligand-binding) tend to be conserved. Meanwhile, areas that interface with co-regulatory molecules can exhibit notable variation even between closely related species. We propose that accumulating mutations in regulatory regions is one way that AR structure might act as a substrate for selection to guide the evolution of reproductive traits. By synthesizing literature across disciplines and highlighting the evolutionary potential of specific AR regions, we hope to inspire new avenues of integrative research into endocrine system evolution.

Massive in Silico Study of Noble Gas Binding to the Structural Proteome

Noble gases are chemically inert, and it was therefore thought they would have little effect on biology. Paradoxically, it was found that they do exhibit a wide range of biological effects, many of which are target-specific and potentially useful and some of which have been demonstrated in vivo. The underlying mechanisms by which useful pharmacology, such as tissue and neuroprotection, anti-addiction effects, and analgesia, is elicited are relatively unexplored. Experiments to probe the interactions of noble gases with specific proteins are more difficult with gases than those with other chemicals. It is clearly impractical to conduct the large number of gas-protein experiments required to gain a complete picture of noble gas biology. Given the simplicity of atoms as ligands, in silico methods provide an opportunity to gain insight into which noble gas-protein interactions are worthy of further experimental or advanced computational investigation. Our previous validation studies showed that in silico methods can accurately predict experimentally determined noble gas binding sites in X-ray structures of proteins. Here, we summarize the largest reported in silico reverse docking study involving 127 854 protein structures and the five nonradioactive noble gases. We describe how these computational screening methods are implemented, summarize the main types of interactions that occur between noble gases and target proteins, describe how the massive data set that this study generated can be analyzed (freely available at group18.csiro.au), and provide the NDMA receptor as an example of how these data can be used to understand the molecular pharmacology underlying the biology of the noble gases. We encourage chemical biologists to access the data and use them to expand the knowledge base of noble gas pharmacology, and to use this information, together with more efficient delivery systems, to develop "atomic drugs" that can fully exploit their considerable and relatively unexplored potential in medicine.

On identifying collective displacements in apo-proteins that reveal eventual binding pathways

Binding of small molecules to proteins often involves large conformational changes in the latter, which open up pathways to the binding site. Observing and pinpointing these rare events in large scale, all-atom, computations of specific protein-ligand complexes, is expensive and to a great extent serendipitous. Further, relevant collective variables which characterise specific binding or un-binding scenarios are still difficult to identify despite the large body of work on the subject. Here, we show that possible primary and secondary binding pathways can be discovered from short simulations of the apo-protein without waiting for an actual binding event to occur. We use a projection formalism, introduced earlier to study deformation in solids, to analyse local atomic displacements into two mutually orthogonal subspaces—those which are “affine” i.e. expressible as a homogeneous deformation of the native structure, and those which are not. The susceptibility to non-affine displacements among the various residues in the apo- protein is then shown to correlate with typical binding pathways and sites crucial for allosteric modifications. We validate our observation with all-atom computations of three proteins, T4-Lysozyme, Src kinase and Cytochrome P450.

Designing drugs which target specific proteins involved in diseases consumes a lot of time and effort in the pharmaceutical industry. In recent times, in silico design of drugs using all-atom molecular modelling has started to provide crucial inputs. Even so, discovery of binding pathways of small molecules both at the primary binding site, as well as sites for allosteric control, is time consuming and often fortuitous. We provide here a framework within which critical conformational changes likely to occur during binding are quantified from statistical analysis of configurations of proteins in their apo, or inactive form, greatly simplifying identification of target residues. We illustrate this idea by analysing ligand binding pathways for three proteins T4- Lysozyme, P450 and Src kinase, which are active respectively in the immune system, metabolism and cancer.

Decoding the Rich Biological Properties of Noble Gases: How Well Can We Predict Noble Gas Binding to Diverse Proteins?

The chemically inert noble gases display a surprisingly rich spectrum of useful biological properties. Relatively little is known about the molecular mechanisms behind these effects. It is clearly not feasible to conduct large numbers of pharmacological experiments on noble gases to identify activity. Computational studies of the binding of noble gases and proteins can address this paucity of information and provide insight into mechanisms of action. We used bespoke computational grid calculations to predict the positions of energy minima in the interactions of noble gases with diverse proteins. The method was validated by quantifying how well simulations could predict binding positions in 131 diverse protein X-ray structures containing 399 Xe and Kr atoms. We found excellent agreement between calculated and experimental binding positions of noble gases. 94 % of all crystallographic xenon atoms were within 1 Xe van der Waals (vdW) diameter of a predicted binding site, and 97 % lay within 2 vdW diameters. 100 % of crystallographic krypton atoms were within 1 Kr vdW diameter of a predicted binding site. We showed the feasibility of large-scale computational screening of all ≈60 000 unique structures in the Protein Data Bank. This will elucidate biochemical mechanisms by which these novel 'atomic drugs' elicit their valuable biochemical properties and identify new medical uses.

Atomic resolution mechanism of ligand binding to a solvent inaccessible cavity in T4 lysozyme

Ligand binding sites in proteins are often localized to deeply buried cavities, inaccessible to bulk solvent. Yet, in many cases binding of cognate ligands occurs rapidly. An intriguing system is presented by the L99A cavity mutant of T4 Lysozyme (T4L L99A) that rapidly binds benzene (~106 M-1s-1). Although the protein has long served as a model system for protein thermodynamics and crystal structures of both free and benzene-bound T4L L99A are available, the kinetic pathways by which benzene reaches its solvent-inaccessible binding cavity remain elusive. The current work, using extensive molecular dynamics simulation, achieves this by capturing the complete process of spontaneous recognition of benzene by T4L L99A at atomistic resolution. A series of multi-microsecond unbiased molecular dynamics simulation trajectories unequivocally reveal how benzene, starting in bulk solvent, diffuses to the protein and spontaneously reaches the solvent inaccessible cavity of T4L L99A. The simulated and high-resolution X-ray derived bound structures are in excellent agreement. A robust four-state Markov model, developed using cumulative 60 μs trajectories, identifies and quantifies multiple ligand binding pathways with low activation barriers. Interestingly, none of these identified binding pathways required large conformational changes for ligand access to the buried cavity. Rather, these involve transient but crucial opening of a channel to the cavity via subtle displacements in the positions of key helices (helix4/helix6, helix7/helix9) leading to rapid binding. Free energy simulations further elucidate that these channel-opening events would have been unfavorable in wild type T4L. Taken together and via integrating with results from experiments, these simulations provide unprecedented mechanistic insights into the complete ligand recognition process in a buried cavity. By illustrating the power of subtle helix movements in opening up multiple pathways for ligand access, this work offers an alternate view of ligand recognition in a solvent-inaccessible cavity, contrary to the common perception of a single dominant pathway for ligand binding.

Nuclear magnetic resonance-based determination of dioxygen binding sites in protein cavities

The location and ligand accessibility of internal cavities in cysteine-free wild-type T4 lysozyme was investigated using O₂ gas-pressure NMR spectroscopy and molecular dynamics (MD) simulation. Upon increasing the concentration of dissolved O₂ in solvent to 8.9 mM, O₂ -induced paramagnetic relaxation enhancements (PREs) to the backbone amide and side chain methyl protons were observed, specifically around two cavities in the C-terminal domain. To determine the number of O₂ binding sites and their atomic coordinates from the 1/r⁶ distance dependence of the PREs, we established an analytical procedure using Akaike's Information Criterion, in combination with a grid-search. Two O₂ -accessible sites were identified in internal cavities: One site was consistent with the xenon-binding site in the protein in crystal, and the other site was established to be a novel ligand-binding site. MD simulations performed at 10 and 100 mM O₂ revealed dioxygen ingress and egress as well as rotational and translational motions of O₂ in the cavities. It is therefore suggested that conformational fluctuations within the ground-state ensemble transiently develop channels for O₂ association with the internal protein cavities.

Mapping Hydrophobic Tunnels and Cavities in Neuroglobin with Noble Gas under Pressure

Internal cavities are crucial for conformational flexibility of proteins and can be mapped through noble gas diffusion and docking. Here we investigate the hydrophobic cavities and tunnel network in neuroglobin (Ngb), a hexacoordinated heme protein likely to be involved in neuroprotection, using crystallography under noble gas pressure, mostly at room temperature. In murine Ngb, a large internal cavity is involved in the heme sliding mechanism to achieve binding of gaseous ligands through coordination to the heme iron. In this study, we report that noble gases are hosted by two major sites within the internal cavity. We propose that these cavities could store oxygen and allow its relay in the heme proximity, which could correspond to NO location in the nitrite-reductase function of Ngb. Thanks to a recently designed pressurization cell using krypton at high pressure, a new gas binding site has been characterized that reveals an alternate pathway for gaseous ligands. A new gas binding site on the proximal side of the heme has also been characterized, using xenon pressure on a Ngb mutant (V140W) that binds CO with a similar rate and affinity to the wild-type, despite a reshaping of the internal cavity. Moreover, this study, to our knowledge, provides new insights into the determinants of the heme sliding mechanism, suggesting that the shift at the beginning of helix G precedes and drives this process.

Xenon for tunnelling analysis of the efflux pump component OprN

Tripartite efflux pumps are among the main actors responsible for antibiotics resistance in Gram-negative bacteria. In the last two decades, structural studies gave crucial information about the assembly interfaces and the mechanistic motions. Thus rigidifying the assembly seems to be an interesting way to hamper the drug efflux. In this context, xenon is a suitable probe for checking whether small ligands could act as conformational lockers by targeting hydrophobic cavities. Here we focus on OprN, the outer membrane channel of the MexEF efflux pump from Pseudomonas aeruginosa. After exposing OprN crystals to xenon gas pressure, 14 binding sites were observed using X-ray crystallography. These binding sites were unambiguously characterized in hydrophobic cavities of OprN. The major site is observed in the sensitive iris-like region gating the channel at the periplasmic side, built by the three key-residues Leu 405, Asp 109, and Arg 412. This arrangement defines along the tunnel axis a strong hydrophobic/polar gradient able to enhance the passive efflux mechanism of OprN. The other xenon atoms reveal strategic hydrophobic regions of the channel scaffold to target, with the aim to freeze the dynamic movements responsible of the open/close conformational equilibrium in OprN.

Argon Delays Initiation of Liver Regeneration after Partial Hepatectomy in Rats

BACKGROUND

The liver can heal up to restitutio ad integrum following damage resulting from various causes. Different studies have demonstrated the protective effect of argon on various cells and organs. To the best of our knowledge, the organ-protective effects of the noble gas argon on the liver have not yet been investigated, although argon appears to influence signal paths that are well-known mediators of liver regeneration. We hypothesized that argon inhalation prior to partial hepatectomy (70%) has a positive effect on the initiation of liver regeneration in rats.

METHODS

Partial hepatectomy (70%) with or without inhaled argon (50 vol%) was performed for 1 h. Liver tissue was harvested after 3, 36, and 96 h to analyze the mRNA and protein expression of hepatocyte growth factor (HGF), interleukin-6 (IL-6), tumor necrosis factor-α, and extracellular signal-regulated kinase 1/2. Histological tissue samples were prepared for immunohistochemistry (bromodeoxyuridine [BrdU], Ki-67, and TUNEL) and blood was analyzed regarding the effects of argon on liver function. Statistical analyses were performed using 1-way ANOVA followed by the post hoc Tukey-Kramer test.

RESULTS

After 3 h, the primary outcome parameter of hepatocyte proliferation was significantly reduced with argon 50 vol% inhalation in comparison to nitrogen inhalation (BrdU: 15.7 ± 9.7 vs. 7.7 ± 3.1 positive cells/1,000 hepatocytes, p = 0.013; Ki-67: 17.6 ± 13.3 vs. 4.7 ± 5.4 positive cells/1,000 hepatocytes, p = 0.006). This was most likely mediated by significant downregulation of HGF (after 3 h: 5.2 ± 3.2 vs. 2.3 ± 1.0 fold, p = 0.032; after 96 h: 2.1 ± 0.5 vs. 1.3 ± 0.3 fold, p = 0.029) and IL-6 (after 3 h: 43.7 ± 39.6 vs. 8.5 ± 9.2 fold, p = 0.032). Nevertheless, we could detect no significant effect on the weight of the residual liver, liver-body weight ratio, or liver blood test results after argon inhalation.

CONCLUSION

Impairment of liver regeneration was apparent after argon 50 vol% inhalation that was most probably mediated by downregulation of HGF and IL-6 in the initial phase. However, the present study was not adequately powered to prove that argon has detrimental effects on the liver. Further studies are needed to evaluate the effects of argon on livers with preexisting conditions as well as on ischemia-reperfusion models.

The Molecular Pathway of Argon-Mediated Neuroprotection

The noble gas argon has attracted increasing attention in recent years, especially because of its neuroprotective properties. In a variety of models, ranging from oxygen-glucose deprivation in cell culture to complex models of mid-cerebral artery occlusion, subarachnoid hemorrhage or retinal ischemia-reperfusion injury in animals, argon administration after individual injury demonstrated favorable effects, particularly increased cell survival and even improved neuronal function. As an inert molecule, argon did not show signs of adverse effects in the in vitro and in vivo model used, while being comparably cheap and easy to apply. However, the molecular mechanism by which argon is able to exert its protective and beneficial characteristics remains unclear. Although there are many pieces missing to complete the signaling pathway throughout the cell, it is the aim of this review to summarize the known parts of the molecular pathways and to combine them to provide a clear insight into the cellular pathway, starting with the receptors that may be involved in mediating argons effects and ending with the translational response.

Structural Basis for Xenon Inhibition in a Cationic Pentameric Ligand-Gated Ion Channel

GLIC receptor is a bacterial pentameric ligand-gated ion channel whose action is inhibited by xenon. Xenon has been used in clinical practice as a potent gaseous anaesthetic for decades, but the molecular mechanism of interactions with its integral membrane receptor targets remains poorly understood. Here we characterize by X-ray crystallography the xenon-binding sites within both the open and "locally-closed" (inactive) conformations of GLIC. Major binding sites of xenon, which differ between the two conformations, were identified in three distinct regions that all belong to the trans-membrane domain of GLIC: 1) in an intra-subunit cavity, 2) at the interface between adjacent subunits, and 3) in the pore. The pore site is unique to the locally-closed form where the binding of xenon effectively seals the channel. A putative mechanism of the inhibition of GLIC by xenon is proposed, which might be extended to other pentameric cationic ligand-gated ion channels.

Pred-binding: large-scale protein-ligand binding affinity prediction

Drug target interactions (DTIs) are crucial in pharmacology and drug discovery. Presently, experimental determination of compound-protein interactions remains challenging because of funding investment and difficulties of purifying proteins. In this study, we proposed two in silico models based on support vector machine (SVM) and random forest (RF), using 1589 molecular descriptors and 1080 protein descriptors in 9948 ligand-protein pairs to predict DTIs that were quantified by Ki values. The cross-validation coefficient of determination of 0.6079 for SVM and 0.6267 for RF were obtained, respectively. In addition, the two-dimensional (2D) autocorrelation, topological charge indices and three-dimensional (3D)-MoRSE descriptors of compounds, the autocorrelation descriptors and the amphiphilic pseudo-amino acid composition of protein are found most important for Ki predictions. These models provide a new opportunity for the prediction of ligand-receptor interactions that will facilitate the target discovery and toxicity evaluation in drug development.

Detecting O2 binding sites in protein cavities

Internal cavities are important elements in protein structure, dynamics, stability and function. Here we use NMR spectroscopy to investigate the binding of molecular oxygen (O₂) to cavities in a well-studied model for ligand binding, the L99A mutant of T4 lysozyme. On increasing the O₂ concentration to 8.9 mM, changes in ¹H, ¹⁵N, and ¹³C chemical shifts and signal broadening were observed specifically for backbone amide and side chain methyl groups located around the two hydrophobic cavities of the protein. O₂-induced longitudinal relaxation enhancements for amide and methyl protons could be adequately accounted for by paramagnetic dipolar relaxation. These data provide the first experimental demonstration that O₂ binds specifically to the hydrophobic, and not the hydrophilic cavities, in a protein. Molecular dynamics simulations visualized the rotational and translational motions of O₂ in the cavities, as well as the binding and egress of O₂, suggesting that the channel consisting of helices D, E, G, H, and J could be the potential gateway for ligand binding to the protein. Due to strong paramagnetic relaxation effects, O₂ gas-pressure NMR measurements can detect hydrophobic cavities when populated to as little as 1%, and thereby provide a general and highly sensitive method for detecting oxygen binding in proteins.

Polyelectrolyte conformational transition in aqueous solvent mixture influenced by hydrophobic interactions and hydrogen bonding effects: PAA-water-ethanol

Molecular dynamics simulations of poly(acrylic acid) PAA chain in water-ethanol mixture were performed for un-ionized and ionized cases at different degree-of-ionization 0%, 80% and 100% of PAA chain by Na(+) counter-ions and co-solvent (ethanol) concentration in the range 0-90vol% ethanol. Aspects of structure and dynamics were investigated via atom pair correlation functions, number and relaxation of hydrogen bonds, nearest-neighbor coordination numbers, and dihedral angle distribution function for back-bone and side-groups of the chain. With increase in ethanol concentration, chain swelling is observed for un-ionized chain (f=0) and on the contrary chain shrinkage is observed for partially and fully ionized cases (i.e., f=0.8 and 1). For un-ionized PAA, with increase in ethanol fraction ϕeth the number of PAA-ethanol hydrogen bonds increases while PAA-water decreases. Increase in ϕeth leads to PAA chain expansion for un-ionized case and chain shrinkage for ionized case, in agreement with experimental observations on this system. For ionized-PAA case, chain shrinkage is found to be influenced by intermolecular hydrogen bonding with water as well as ethanol. The localization of ethanol molecules near the un-ionized PAA backbone at higher levels of ethanol is facilitated by a displacement of water molecules indicating presence of specific ethanol hydration shell, as confirmed by results of the RDF curves and coordination number calculations. This behavior, controlled by hydrogen bonding provides a significant contribution to such a conformational transition behavior of the polyelectrolyte chain. The interactions between counter-ions and charges on the PAA chain also influence chain collapse. The underlying origins of polyelectrolyte chain collapse in water-alcohol mixtures are brought out for the first time via explicit MD simulations by this study.

Cavity as a source of conformational fluctuation and high-energy state: high-pressure NMR study of a cavity-enlarged mutant of T4 lysozyme

Although the structure, function, conformational dynamics, and controlled thermodynamics of proteins are manifested by their corresponding amino acid sequences, the natural rules for molecular design and their corresponding interplay remain obscure. In this study, we focused on the role of internal cavities of proteins in conformational dynamics. We investigated the pressure-induced responses from the cavity-enlarged L99A mutant of T4 lysozyme, using high-pressure NMR spectroscopy. The signal intensities of the methyl groups in the (1)H/(13)C heteronuclear single quantum correlation spectra, particularly those around the enlarged cavity, decreased with the increasing pressure, and disappeared at 200 MPa, without the appearance of new resonances, thus indicating the presence of heterogeneous conformations around the cavity within the ground state ensemble. Above 200 MPa, the signal intensities of >20 methyl groups gradually decreased with the increasing pressure, without the appearance of new resonances. Interestingly, these residues closely matched those sensing a large conformational change between the ground- and high-energy states, at atmospheric pressure. (13)C and (1)H NMR line-shape simulations showed that the pressure-induced loss in the peak intensity could be explained by the increase in the high-energy state population. In this high-energy state, the aromatic side chain of F114 gets flipped into the enlarged cavity. The accommodation of the phenylalanine ring into the efficiently packed cavity may decrease the partial molar volume of the high-energy state, relative to the ground state. We suggest that the enlarged cavity is involved in the conformational transition to high-energy states and in the volume fluctuation of the ground state.

Volatile anesthetics inhibit sodium channels without altering bulk lipid bilayer properties

Although general anesthetics are clinically important and widely used, their molecular mechanisms of action remain poorly understood. Volatile anesthetics such as isoflurane (ISO) are thought to alter neuronal function by depressing excitatory and facilitating inhibitory neurotransmission through direct interactions with specific protein targets, including voltage-gated sodium channels (Na(v)). Many anesthetics alter lipid bilayer properties, suggesting that ion channel function might also be altered indirectly through effects on the lipid bilayer. We compared the effects of ISO and of a series of fluorobenzene (FB) model volatile anesthetics on Na(v) function and lipid bilayer properties. We examined the effects of these agents on Na(v) in neuronal cells using whole-cell electrophysiology, and on lipid bilayer properties using a gramicidin-based fluorescence assay, which is a functional assay for detecting changes in lipid bilayer properties sensed by a bilayer-spanning ion channel. At clinically relevant concentrations (defined by the minimum alveolar concentration), both the FBs and ISO produced prepulse-dependent inhibition of Na(v) and shifted the voltage dependence of inactivation toward more hyperpolarized potentials without affecting lipid bilayer properties, as sensed by gramicidin channels. Only at supra-anesthetic (toxic) concentrations did ISO alter lipid bilayer properties. These results suggest that clinically relevant concentrations of volatile anesthetics alter Na(v) function through direct interactions with the channel protein with little, if any, contribution from changes in bulk lipid bilayer properties. Our findings further suggest that changes in lipid bilayer properties are not involved in clinical anesthesia.

Crystallographic Studies with Xenon and Nitrous Oxide Provide Evidence for Protein-dependent Processes in the Mechanisms of General Anesthesia

Background:

The mechanisms by which general anesthetics, including xenon and nitrous oxide, act are only beginning to be discovered. However, structural approaches revealed weak but specific protein–gas interactions.

Methods:

To improve knowledge, we performed x-ray crystallography studies under xenon and nitrous oxide pressure in a series of 10 binding sites within four proteins.

Results:

Whatever the pressure, we show (1) hydrophobicity of the gas binding sites has a screening effect on xenon and nitrous oxide binding, with a threshold value of 83% beyond which and below which xenon and nitrous oxide, respectively, binds to their sites preferentially compared to each other; (2) xenon and nitrous oxide occupancies are significantly correlated respectively to the product and the ratio of hydrophobicity by volume, indicating that hydrophobicity and volume are binding parameters that complement and oppose each other’s effects; and (3) the ratio of occupancy of xenon to nitrous oxide is significantly correlated to hydrophobicity of their binding sites.

Conclusions:

These data demonstrate that xenon and nitrous oxide obey different binding mechanisms, a finding that argues against all unitary hypotheses of narcosis and anesthesia, and indicate that the Meyer–Overton rule of a high correlation between anesthetic potency and solubility in lipids of general anesthetics is often overinterpreted. This study provides evidence that the mechanisms of gas binding to proteins and therefore of general anesthesia should be considered as the result of a fully reversible interaction between a drug ligand and a receptor as this occurs in classical pharmacology.

Argon: systematic review on neuro- and organoprotective properties of an "inert" gas

Argon belongs to the group of noble gases, which are regarded as chemically inert. Astonishingly some of these gases exert biological properties and during the last decades more and more reports demonstrated neuroprotective and organoprotective effects. Recent studies predominately use in vivo or in vitro models for ischemic pathologies to investigate the effect of argon treatment. Promising data has been published concerning pathologies like cerebral ischemia, traumatic brain injury and hypoxic ischemic encephalopathy. However, models applied and administration of the therapeutic gas vary. Here we provide a systematic review to summarize the available data on argon's neuro- and organoprotective effects and discuss its possible mechanism of action. We aim to provide a summary to allow further studies with a more homogeneous setting to investigate possible clinical applications of argon.

Predicting when biliary excretion of parent drug is a major route of elimination in humans

Biliary excretion is an important route of elimination for many drugs, yet measuring the extent of biliary elimination is difficult, invasive, and variable. Biliary elimination has been quantified for few drugs with a limited number of subjects, who are often diseased patients. An accurate prediction of which drugs or new molecular entities are significantly eliminated in the bile may predict potential drug-drug interactions, pharmacokinetics, and toxicities. The Biopharmaceutics Drug Disposition Classification System (BDDCS) characterizes significant routes of drug elimination, identifies potential transporter effects, and is useful in understanding drug-drug interactions. Class 1 and 2 drugs are primarily eliminated in humans via metabolism and will not exhibit significant biliary excretion of parent compound. In contrast, class 3 and 4 drugs are primarily excreted unchanged in the urine or bile. Here, we characterize the significant elimination route of 105 orally administered class 3 and 4 drugs. We introduce and validate a novel model, predicting significant biliary elimination using a simple classification scheme. The model is accurate for 83% of 30 drugs collected after model development. The model corroborates the observation that biliarily eliminated drugs have high molecular weights, while demonstrating the necessity of considering route of administration and extent of metabolism when predicting biliary excretion. Interestingly, a predictor of potential metabolism significantly improves predictions of major elimination routes of poorly metabolized drugs. This model successfully predicts the major elimination route for poorly permeable/poorly metabolized drugs and may be applied prior to human dosing.

Functional relevance of the internal hydrophobic cavity of urate oxidase

Urate oxidase from Aspergillus flavus is a 135 kDa homo-tetramer which has a hydrophobic cavity buried within each monomer and located close to its active site. Crystallographic studies under moderate gas pressure and high hydrostatic pressure have shown that both gas presence and high pressure would rigidify the cavity leading to an inhibition of the catalytic activity. Analysis of the cavity volume variations and functional modifications suggest that the flexibility of the cavity would be an essential parameter for the active site efficiency. This cavity would act as a connecting vessel to give flexibility to the neighboring active site, and its expansion under pure oxygen pressure reveals that it might serve as a transient reservoir on its pathway to the active site.

Argon gas: a potential neuroprotectant and promising medical therapy

Argon is a noble gas element that has demonstrated narcotic and protective abilities that may prove useful in the medical field. The earliest records of argon gas have exposed its ability to exhibit narcotic symptoms at hyperbaric pressures greater than 10 atmospheres with more recent evidence seeking to display argon as a potential neuroprotective agent. The high availability and low cost of argon provide a distinct advantage over using similarly acting treatments such as xenon gas. Argon gas treatments in models of brain injury such as in vitro Oxygen-Glucose-Deprivation (OGD) and Traumatic Brain Injury (TBI), as well as in vivo Middle Cerebral Artery Occlusion (MCAO) have largely demonstrated positive neuroprotective behavior. On the other hand, some warning has been made to potential negative effects of argon treatments in cases of ischemic brain injury, where increases of damage in the sub-cortical region of the brain have been uncovered. Further support for argon use in the medical field has been demonstrated in its use in combination with tPA, its ability as an organoprotectant, and its surgical applications. This review seeks to summarize the history and development of argon gas use in medical research as mainly a neuroprotective agent, to summarize the mechanisms associated with its biological effects, and to elucidate its future potential.

Structure and dynamics of an imidazoline nitroxide side chain with strongly hindered internal motion in proteins

A disulfide-linked imidazoline nitroxide side chain (V1) has a similar and highly constrained internal motion at diverse topological sites in a protein, unlike that for the disulfide-linked pyrroline nitroxide side chain (R1) widely used in site directed spin labeling EPR. Crystal structures of V1 at two positions in a helix of T4 Lysozyme and quantum mechanical calculations suggest the source of the constraints as intra-side chain interactions of the disulfide sulfur atoms with both the protein backbone and the 3-nitrogen in the imidazoline ring. These interactions apparently limit the conformation of the side chain to one of only three possible rotamers, two of which are observed in the crystal structure. An inter-spin distance measurement in frozen solution using double electron-electron resonance (DEER) gives a value essentially identical to that determined from the crystal structure of the protein containing two copies of V1, indicating that lattice forces do not dictate the rotamers observed. Collectively, the results suggest the possibility of predetermining a unique rotamer of V1 in helical structures. In general, the reduced rotameric space of V1 compared to R1 should simplify interpretation of inter-spin distance information in terms of protein structure, while the highly constrained internal motion is expected to extend the dynamic range for characterizing large amplitude nanosecond backbone fluctuations.

Temperature and pressure dependence of azurin stability as monitored by tryptophan fluorescence and phosphorescence. The case of F29A mutant

The effects of a single-point, F29A, cavity-forming mutation on the unfolding thermodynamic parameters of azurin from Pseudomonas aeruginosa and on the internal dynamics of the protein fold under pressure were probed by the fluorescence and phosphorescence emission of Trp48, deeply buried in the compact hydrophobic core of the macromolecule. Pressure-induced unfolding, monitored by the shift in the fluorescence spectrum, led to a volume change of 70-90mlmol(-1). The difference in the unfolding volume between F29A and wild type azurin was smaller than the volume of the space theoretically created in the mutant, indicating that the cavity is, at least partially, filled with water molecules. The complex temperature dependence of the unfolding volume, for temperatures up to 20°C, suggests the formation of an expanded form of the protein and highlights how the packing efficiency of azurin appears to contribute to the magnitude of internal void volume at any given temperature. Changes in flexibility of the protein matrix around the chromophore were monitored by the intrinsic phosphorescence lifetime. At 40°C the application of pressure in the predenaturation range initially decreases the internal flexibility of azurin, the trend eventually reverting on approaching unfolding. The main difference between wild type and the cavity mutant is the inversion point which happens at 300MPa for wild type and at 150MPa for F29A. This suggests that, for the cavity mutant, pressure-induced internal hydration is more dominant than any compaction of the globular fold at relatively low pressures.

Tunnels modulate ligand flux in a heme nitric oxide/oxygen binding (H-NOX) domain

Interior topological features, such as pockets and channels, have evolved in proteins to regulate biological functions by facilitating the diffusion of biomolecules. Decades of research using the globins as model heme proteins have clearly highlighted the importance of gas pockets around the heme in controlling the capture and release of O(2). However, much less is known about how ligand migration contributes to the diverse functions of other heme protein scaffolds. Heme nitric oxide/oxygen binding (H-NOX) domains are a conserved family of gas-sensing heme proteins with a divergent fold that are critical to numerous signaling pathways. Utilizing X-ray crystallography with xenon, a tunnel network has been shown to serve as a molecular pathway for ligand diffusion. Structure-guided mutagenesis results show that the tunnels have unexpected effects on gas-sensing properties in H-NOX domains. The findings provide insights on how the flux of biomolecules through protein scaffolds modulates protein chemistry.

Structure and dynamics of a conformationally constrained nitroxide side chain and applications in EPR spectroscopy

A disulfide-linked nitroxide side chain (R1) is the most widely used spin label for determining protein topology, mapping structural changes, and characterizing nanosecond backbone motions by site-directed spin labeling. Although the internal motion of R1 and the number of preferred rotamers are limited, translating interspin distance measurements and spatial orientation information into structural constraints is challenging. Here, we introduce a highly constrained nitroxide side chain designated RX as an alternative to R1 for these applications. RX is formed by a facile cross-linking reaction of a bifunctional methanethiosulfonate reagent with pairs of cysteine residues at i and i + 3 or i and i + 4 in an α-helix, at i and i + 2 in a β-strand, or with cysteine residues in adjacent strands in a β-sheet. Analysis of EPR spectra, a crystal structure of RX in T4 lysozyme, and pulsed electron-electron double resonance (ELDOR) spectroscopy on an immobilized protein containing RX all reveal a highly constrained internal motion of the side chain. Consistent with the constrained geometry, interspin distance distributions between pairs of RX side chains are narrower than those from analogous R1 pairs. As an important consequence of the constrained internal motion of RX, spectral diffusion detected with ELDOR reveals microsecond internal motions of the protein. Collectively, the data suggest that the RX side chain will be useful for distance mapping by EPR spectroscopy, determining spatial orientation of helical segments in oriented specimens, and measuring structural fluctuations on the microsecond time scale.

Measurement of radon and xenon binding to a cryptophane molecular host

Xenon and radon have many similar properties, a difference being that all 35 isotopes of radon ((195)Rn-(229)Rn) are radioactive. Radon is a pervasive indoor air pollutant believed to cause significant incidence of lung cancer in many geographic regions, yet radon affinity for a discrete molecular species has never been determined. By comparison, the chemistry of xenon has been widely studied and applied in science and technology. Here, both noble gases were found to bind with exceptional affinity to tris-(triazole ethylamine) cryptophane, a previously unsynthesized water-soluble organic host molecule. The cryptophane-xenon association constant, K(a)=42,000 ± 2,000 M(-1) at 293 K, was determined by isothermal titration calorimetry. This value represents the highest measured xenon affinity for a host molecule. The partitioning of radon between air and aqueous cryptophane solutions of varying concentration was determined radiometrically to give the cryptophane-radon association constant K(a)=49,000 ± 12,000 M(-1) at 293 K.

Protein dynamics and pressure: what can high pressure tell us about protein structural flexibility?

After a brief overview of NMR and X-ray crystallography studies on protein flexibility under pressure, we summarize the effects of hydrostatic pressure on the native fold of azurin from Pseudomonas aeruginosa as inferred from the variation of the intrinsic phosphorescence lifetime and the acrylamide bimolecular quenching rate constants of the buried Trp residue. The pressure/temperature response of the globular fold and modulation of its dynamical structure is analyzed both in terms of a reduction of internal cavities and of the hydration of the polypeptide. The study of the effect of two single point cavity forming mutations, F110S and I7S, on the unfolding volume change (ΔV(0)) of azurin and on the internal dynamics of the protein fold under pressure demonstrate that the creation of an internal cavity will enhance the plasticity and lower the stability of the globular structure. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Molecular targets and mechanisms for ethanol action in glycine receptors

Glycine receptors (GlyRs) are recognized as the primary mediators of neuronal inhibition in the spinal cord, brain stem and higher brain regions known to be sensitive to ethanol. Building evidence supports the notion that ethanol acting on GlyRs causes at least a subset of its behavioral effects and may be involved in modulating ethanol intake. For over two decades, GlyRs have been studied at the molecular level as targets for ethanol action. Despite the advances in understanding the effects of ethanol in vivo and in vitro, the precise molecular sites and mechanisms of action for ethanol in ligand-gated ion channels in general, and in GlyRs specifically, are just now starting to become understood. The present review focuses on advances in our knowledge produced by using molecular biology, pressure antagonism, electrophysiology and molecular modeling strategies over the last two decades to probe, identify and model the initial molecular sites and mechanisms of ethanol action in GlyRs. The molecular targets on the GlyR are covered on a global perspective, which includes the intracellular, transmembrane and extracellular domains. The latter has received increasing attention in recent years. Recent molecular models of the sites of ethanol action in GlyRs and their implications to our understanding of possible mechanism of ethanol action and novel targets for drug development in GlyRs are discussed.

Trypsin-ligand binding free energy calculation with AMOEBA

The binding free energies of several benzamidine -like inhibitors to trypsin were examined using a polarizable molecular mechanics potential. All the computed binding free energies are in good agreement with the experimental data. From free energy decomposition, electrostatic interaction was indicated to be the driving force for the binding. MD simulations show that the ligands form hydrogen bonds with trypsin and water molecules nearby in a competitive fashion. While the binding free energy is independent of the ligand dipole moment, it shows a strong correlation with the ligand molecular polarizability.