Ribonucleotide Reductase Requires Subunit Switching in Hypoxia to Maintain DNA Replication.

Cells exposed to hypoxia experience replication stress but do not accumulate DNA damage, suggesting sustained DNA replication. Ribonucleotide reductase (RNR) is the only enzyme capable of de novo synthesis of deoxyribonucleotide triphosphates (dNTPs). However, oxygen is an essential cofactor for mammalian RNR (RRM1/RRM2 and RRM1/RRM2B), leading us to question the source of dNTPs in hypoxia. Here, we show that the RRM1/RRM2B enzyme is capable of retaining activity in hypoxia and therefore is favored over RRM1/RRM2 in order to preserve ongoing replication and avoid the accumulation of DNA damage. We found two distinct mechanisms by which RRM2B maintains hypoxic activity and identified responsible residues in RRM2B. The importance of RRM2B in the response to tumor hypoxia is further illustrated by correlation of its expression with a hypoxic signature in patient samples and its roles in tumor growth and radioresistance. Our data provide mechanistic insight into RNR biology, highlighting RRM2B as a hypoxic-specific, anti-cancer therapeutic target.

Building predictive Markov models of ion channel permeation from molecular dynamics simulations.

Molecular dynamics (MD) simulation of biological processes has always been a challenging task due to the long timescales of the processes involved and the large amount of output data to handle. Markov state models (MSMs) have been introduced as a powerful tool in this area of research, as they provide a mechanistically comprehensible synthesis of the large amount of MD data and, at the same time, can be used to rapidly estimate experimental properties of biological processes. Herein, we propose a method for building MSMs of ion channel permeation from MD trajectories, which directly evaluates the current flowing through the channel from the model's transition matrix (T), which is crucial for comparing simulations and experimental data. This is achieved by including in the model a flux matrix that summarizes information on the charge moving across the channel between each pair of states of the MSM and can be used in conjunction with T to predict the ion current. A procedure to drastically reduce the number of states in the MSM while preserving the estimated ion current is also proposed. Applying the method to the KcsA channel returned an MSM with five states with significant equilibrium occupancy, capable of accurately reproducing the single-channel ion current from microsecond MD trajectories.

Early Steps in C-Type Inactivation of the hERG Potassium Channel.

Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K+-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.

A Potential Route of Capsaicin to Its Binding Site in the TRPV1 Ion Channel.

Transient receptor potential (TRP) ion channels are important pharmacological targets because of their role in the perception of pain, and so, understanding their chemical regulation is essential for the development of analgesic drugs. Among the currently known TRP channel chemical agonists, capsaicin, the active compound of chili pepper, is probably the most exhaustively studied. The availability of the three-dimensional structure of the vanilloid receptor 1 (TRPV1) has fueled computational studies revealing the molecular details of capsaicin binding modes. Although this is a significant step, a comprehensible binding mechanism or pathway is invaluable for targeting TRP channels in modern pharmacology. In the present work, free-energy and enhanced sampling techniques have been used to explore a possible membrane-mediated pathway for capsaicin to enter the TRPV1 binding pocket where capsaicin accesses the protein starting at the extracellular milieu through the outer leaflet and into its binding site in the protein. The main states visited along this route have been characterized and include (i) a bound state in agreement with the binding mode "head-down, tail-up" and (ii) an alternative state corresponding to a "head-up, tail-down" binding mode. In agreement with previous reports, binding is mediated by both hydrogen bonds and van der Waals interactions, and residue Y511 is crucial for stabilizing the bound state and during the binding process. Together, these results provide a foundation to further understand TRPV channels, and they could be used to guide therapeutic design of selective inhibitors potentially leading to novel avenues for pharmacological applications targeting the TRPV1 channel.

Antibiotic resistance and host immune evasion in Staphylococcus aureus mediated by a metabolic adaptation.

Staphylococcus aureus is a notorious human bacterial pathogen with considerable capacity to develop antibiotic resistance. We have observed that human infections caused by highly drug-resistant S. aureus are more prolonged, complicated, and difficult to eradicate. Here we describe a metabolic adaptation strategy used by clinical S. aureus strains that leads to resistance to the last-line antibiotic, daptomycin, and simultaneously affects host innate immunity. This response was characterized by a change in anionic membrane phospholipid composition induced by point mutations in the phospholipid biosynthesis gene, cls2, encoding cardiolipin synthase. Single cls2 point mutations were sufficient for daptomycin resistance, antibiotic treatment failure, and persistent infection. These phenotypes were mediated by enhanced cardiolipin biosynthesis, leading to increased bacterial membrane cardiolipin and reduced phosphatidylglycerol. The changes in membrane phospholipid profile led to modifications in membrane structure that impaired daptomycin penetration and membrane disruption. The cls2 point mutations also allowed S. aureus to evade neutrophil chemotaxis, mediated by the reduction in bacterial membrane phosphatidylglycerol, a previously undescribed bacterial-driven chemoattractant. Together, these data illustrate a metabolic strategy used by S. aureus to circumvent antibiotic and immune attack and provide crucial insights into membrane-based therapeutic targeting of this troublesome pathogen.

The ligand-bound state of a G protein-coupled receptor stabilizes the interaction of functional cholesterol molecules.

Cholesterol is a major component of mammalian plasma membranes that not only affects the physical properties of the lipid bilayer but also is the function of many membrane proteins including G protein-coupled receptors. The oxytocin receptor (OXTR) is involved in parturition and lactation of mammals and in their emotional and social behaviors. Cholesterol acts on OXTR as an allosteric modulator inducing a high-affinity state for orthosteric ligands through a molecular mechanism that has yet to be determined. Using the ion channel-coupled receptor technology, we developed a functional assay of cholesterol modulation of G protein-coupled receptors that is independent of intracellular signaling pathways and operational in living cells. Using this assay, we discovered a stable binding of cholesterol molecules to the receptor when it adopts an orthosteric ligand-bound state. This stable interaction preserves the cholesterol-dependent activity of the receptor in cholesterol-depleted membranes. This mechanism was confirmed using time-resolved FRET experiments on WT OXTR expressed in CHO cells. Consequently, a positive cross-regulation sequentially occurs in OXTR between cholesterol and orthosteric ligands.

Ion Conduction Mechanism as a Fingerprint of Potassium Channels.

K+-channels are membrane proteins that regulate the selective conduction of potassium ions across cell membranes. Although the atomic mechanisms of K+ permeation have been extensively investigated, previous work focused on characterizing the selectivity and occupancy of the binding sites, the role of water molecules in the conduction process, or the identification of the minimum energy pathways enabling permeation. Here, we exploit molecular dynamics simulations and the analytical power of Markov state models to perform a comparative study of ion conduction in three distinct channel models. Significant differences emerged in terms of permeation mechanisms and binding site occupancy by potassium ions and/or water molecules from 100 µs cumulative trajectories. We found that, at odds with the current paradigm, each system displays a characteristic permeation mechanism, and thus, there is not a unique way by which potassium ions move through K+-channels. The high functional diversity of K+-channels can be attributed in part to the differences in conduction features that have emerged from this work. This study provides crucial information and further inspiration for wet-lab chemists designing new synthetic strategies to produce versatile artificial ion channels that emulate membrane transport for their applications in diagnosis, sensors, the next generation of water treatment technologies, etc., as the ability of synthetic channels to transport molecular ions across a bilayer in a controlled way is usually governed through the choice of metal ions, their oxidation states, or their coordination geometries.

An Integrated Mass Spectrometry and Molecular Dynamics Simulations Approach Reveals the Spatial Organization Impact of Metal-Binding Sites on the Stability of Metal-Depleted Metallothionein-2 Species.

Mammalian metallothioneins (MTs) are a group of cysteine-rich proteins that bind metal ions in two α- and ß-domains and represent a major cellular Zn(II)/Cu(I) buffering system in the cell. At cellular free Zn(II) concentrations (10-11-10-9 M), MTs do not exist in fully loaded forms with seven Zn(II)-bound ions (Zn7MTs). Instead, MTs exist as partially metal-depleted species (Zn4-6MT) because their Zn(II) binding affinities are on the nano- to picomolar range comparable to the concentrations of cellular Zn(II). The mode of action of MTs remains poorly understood, and thus, the aim of this study is to characterize the mechanism of Zn(II) (un)binding to MTs, the thermodynamic properties of the Zn1-6MT2 species, and their mechanostability properties. To this end, native mass spectrometry (MS) and label-free quantitative bottom-up and top-down MS in combination with steered molecular dynamics simulations, well-tempered metadynamics (WT-MetaD), and parallel-bias WT-MetaD (amounting to 3.5 µs) were integrated to unravel the chemical coordination of Zn(II) in all Zn1-6MT2 species and to explain the differences in binding affinities of Zn(II) ions to MTs. Differences are found to be the result of the degree of water participation in MT (un)folding and the hyper-reactive character of Cys21 and Cys29 residues. The thermodynamics properties of Zn(II) (un)binding to MT2 are found to differ from those of Cd(II), justifying their distinctive roles. The potential of this integrated strategy in the investigation of numerous unexplored metalloproteins is attested by the results highlighted in the present study.

Multiple Interactions of Glucose with the Extra-Membranous Loops of GLUT1 Aid Transport.

Molecular dynamics simulations amounting to ≈8 µs demonstrate that the glucose transporter GLUT1 undergoes structural fluctuations mediated by the fluidity of the lipid bilayer and the proximity to glucose. The fluctuations of GLUT1 increase as the glucose concentration is raised. These fluctuations are more pronounced when the lipid bilayer is in the fluid compared to the gel phase. Glucose interactions are confined to the extra-membranous residues when the lipid is in the gel phase but diffuses into the transmembrane regions in the fluid phase. Proximity of glucose to GLUT1 causes asynchronous expansions of key bottlenecks at the internal and external openings of the central pore. This is accomplished only by small conformational changes at the single residue level that lower the resistance to glucose movements, thereby permitting unsteered glucose and water movements along the entire length of the pore. When glucose is near salt bridges located at the external and internal openings of the central pore, the distance separating the polar amino acid residues guarding these apertures tends to increase in both fluid and gel phases. It is evident that the multiplicity of glucose interactions, obtained with high concentrations, amplifies the structural fluctuations in GLUT1. The findings that most of the salt bridges and the bottlenecks appear to be operated by glucose proximity suggest that the main triggers to activation of transport are located within the solvent accessible linker regions in the extramembranous zones.

Capturing the Molecular Mechanism of Anesthetic Action by Simulation Methods.

Significant computational efforts have been focused toward exposing the molecular mechanisms of anesthesia in recent years. In the past decade, this has been aided considerably by a momentous increase in the number of high-resolution structures of ion channels, which are putative targets for the anesthetic agents, as well as advancements in high-performance computing technologies. In this review, typical simulation methods to investigate the behavior of model membranes and membrane-protein systems are briefly reviewed, and related computational studies are surveyed. Both lipid- and protein-mediated mechanisms of anesthetic action are scrutinized, focusing on the behavior of ion channels in the latter case.

Ion-triggered selectivity in bacterial sodium channels.

Since the availability of the first crystal structure of a bacterial Na+ channel in 2011, understanding selectivity across this family of membrane proteins has been the subject of intense research efforts. Initially, free energy calculations based on molecular dynamics simulations revealed that although sodium ions can easily permeate the channel with their first hydration shell almost intact, the selectivity filter is too narrow for efficient conduction of hydrated potassium ions. This steric view of selectivity was subsequently questioned by microsecond atomic trajectories, which proved that the selectivity filter appears to the permeating ions as a highly degenerate, liquid-like environment. Although this liquid-like environment looks optimal for rapid conduction of Na+, it seems incompatible with efficient discrimination between similar ion species, such as Na+ and K+, through steric effects. Here extensive molecular dynamics simulations, combined with Markov state model analyses, reveal that at positive membrane potentials, potassium ions trigger a conformational change of the selectivity toward a nonconductive metastable state. It is this transition of the selectivity filter, and not steric effects, that prevents the outward flux of K+ at positive membrane potentials. This description of selectivity, triggered by the nature of the permeating ions, might have implications on the current understanding of how ion channels, and in particular bacterial Na+ channels, operate at the atomic scale.

The Bilayer Collective Properties Govern the Interaction of an HIV-1 Antibody with the Viral Membrane.

Efficient engagement with the envelope glycoprotein membrane-proximal external region (MPER) results in robust blocking of viral infection by a class of broadly neutralizing antibodies (bnAbs) against human immunodeficiency virus (HIV). Developing an accommodation surface that engages with the viral lipid envelope appears to correlate with the neutralizing potency displayed by these bnAbs. The nature of the interactions established between the antibody and the lipid is nonetheless a matter of debate, with some authors arguing that anti-MPER specificity arises only under pathological conditions in autoantibodies endowed with stereospecific binding sites for phospholipids. However, bnAb-lipid interactions are often studied in systems that do not fully preserve the biophysical properties of lipid bilayers, and therefore, questions on binding specificity and the effect of collective membrane properties on the interaction are still open. Here, to evaluate the specificity of lipid interactions of an anti-MPER bnAb (4E10) in an intact membrane context, we determine quantitatively its association with lipid bilayers by means of scanning fluorescence correlation spectroscopy and all-atom molecular dynamic simulations. Our data support that 4E10 establishes electrostatic and hydrophobic interactions with the viral membrane surface and that the collective physical properties of the lipid bilayer influence 4E10 dynamics therein. We conclude that establishment of peripheral, nonspecific electrostatic interactions with the viral membrane through accommodation surfaces may assist high-affinity binding of HIV-1 MPER epitope at membrane interfaces. These findings highlight the importance of considering antibody-lipid interactions in the design of antibody-based anti-HIV strategies.

Mechanism of Molecular Oxygen Diffusion in a Hypoxia-Sensing Prolyl Hydroxylase Using Multiscale Simulation.

The chronic response of animals to hypoxia is mediated by the αß-heterodimeric hypoxia-inducible transcription factors (α,ß-HIFs) which upregulate the expression of sets of genes that work to ameliorate the effects of limiting dioxygen. The HIF prolyl hydroxylase domain enzymes (PHDs) are Fe(II)- and 2-oxoglutarate-dependent oxygenases that act as hypoxia-sensing components of the HIF system: prolyl-hydroxylation signals for dioxgen availability-dependent HIF-α degradation via the ubiquitin proteasome system. The unusual kinetic properties of the PHDs, in particular a high Km for dioxygen and slow reaction with dioxygen, are proposed to enable their hypoxia-sensing role. An understanding of how dioxygen is delivered to, and binds at, the active site of the PHDs is important for the development of a chemical understanding of the hypoxic response. We employed a combined multiscale approach involving classical atomistic equilibrium and nonequilibrium MD simulations combined with QM/MM trajectories to investigate dioxygen diffusion to, and binding at, the active site in the PHD2.Fe(II).2OG.HIF substrate complex; PHD2 is the most important of the three human PHDs. The transport of dioxygen to the active site is described; dioxygen transport follows a single well-defined hydrophobic tunnel, formed from both enzyme and substrate elements to reach the PHD2 active site. The results provide estimates for rate constants that define a diffusion-reaction model for dioxygen:PHD2 interactions; in combination with reported biophysical analyses they provide chemical insight into the basis of the slow reaction of PHD2 with dioxygen. They imply that the reversible binding of dioxygen is central to the hypoxia-sensing capacity of the PHDs and that different PHD HIF-α substrate combinations might have different dioxygen sensitivity profiles. The extent of HIF-α substrate prolyl hydroxylation, which signals for subsequent HIF-α degradation, may thus be a manifestation of the equilibrium between dioxygen in bulk solution and dioxygen bound to the PHD2.Fe.2OG.HIF-α substrate complex.

Conduction and Gating Properties of the TRAAK Channel from Molecular Dynamics Simulations with Different Force Fields.

In recent years, the K2P family of potassium channels has been the subject of intense research activity. Owing to the complex function and regulation of this family of ion channels, it is common practice to complement experimental findings with the atomistic description provided by computational approaches such as molecular dynamics (MD) simulations, especially, in light of the unprecedented timescales accessible at present. However, despite recent substantial improvements, the accuracy of MD simulations is still undermined by the intrinsic limitations of force fields. Here, we systematically assessed the performance of the most popular force fields employed to study ion channels at timescales that are orders of magnitude greater than the ones accessible when these energy functions were first developed. Using 32 µs of trajectories, we investigated the dynamics of a member of the K2P ion channel family, the TRAAK channel, using two established force fields in simulations of biological systems: AMBER and CHARMM. We found that while results are comparable on the nanosecond timescales, significant inconsistencies arise at microsecond timescales.

Differential Stability of Aurein 1.2 Pores in Model Membranes of Two Probiotic Strains.

Aurein 1.2 is an antimicrobial peptide from the skin secretion of an Australian frog. In the previous experimental work, we reported a differential action of aurein 1.2 on two probiotic strains Lactobacillus delbrueckii subsp. bulgaricus (CIDCA 331) and Lactobacillus delbrueckii subsp. lactis (CIDCA 133). The differences found were attributed to the bilayer compositions. Cell cultures and CIDCA 331-derived liposomes showed higher susceptibility than the ones derived from the CIDCA 133 strain, leading to content leakage and structural disruption. Here, we used molecular dynamics simulations to explore these systems at the atomistic level. We hypothesize that if the antimicrobial peptides organized themselves to form a pore, it will be more stable in membranes that emulate the CIDCA 331 strain than in those of the CIDCA 133 strain. To test this hypothesis, we simulated preassembled aurein 1.2 pores embedded into bilayer models that emulate the two probiotic strains. It was found that the general behavior of the systems depends on the composition of the membrane rather than the preassemble system characteristics. Overall, it was observed that aurein 1.2 pores are more stable in the CIDCA 331 model membranes. This fact coincides with the high susceptibility of this strain against antimicrobial peptide. In contrast, in the case of the CIDCA 133 model membranes, peptides migrate to the water-lipid interphase, the pore shrinks, and the transport of water through the pore is reduced. The tendency of glycolipids to make hydrogen bonds with peptides destabilizes the pore structures. This feature is observed to a lesser extent in CIDCA 331 due to the presence of anionic lipids. Glycolipid transverse diffusion (flip-flop) between monolayers occurs in the pore surface region in all the cases considered. These findings expand our understanding of the antimicrobial peptide resistance properties of probiotic strains.

Role of Zinc and Magnesium Ions in the Modulation of Phosphoryl Transfer in Protein Tyrosine Phosphatase 1B.

While the majority of phosphatases are metalloenzymes, the prevailing model for the reactions catalyzed by protein tyrosine phosphatases does not involve any metal ion, yet both metal cations and oxoanions affect their enzymatic activity. Mg2+ and Zn2+ activate and inhibit, respectively, protein tyrosine phosphatase 1B (PTP1B). Molecular dynamics simulations, metadynamics, and quantum chemical calculations in combination with experimental investigations demonstrate that Mg2+ and Zn2+ compete for the same binding site in the active site only in the closed conformation of the enzyme in its phosphorylated state. The two cations have different effects on the arrangements and activities of water molecules that are necessary for the hydrolysis of the phosphocysteine intermediate in the second catalytic step of the reaction. Remarkable differences between the established structural enzymology of PTP1B investigated ex vivo and the function of PTP1B in vivo become evident. Different reaction pathways are viable when the presence of metal ions and their cellular concentrations are considered. The findings suggest that the substrate delivers the inhibitory Zn2+ ion to the active site. The inhibition and activation can be ascribed to the different coordination chemistries of Zn2+ and Mg2+ ions and the orientation of the metal-coordinated water molecules. Metallochemistry adds an additional dimension to the regulation of PTP1B and presumably other members of this enzyme family.

Exposure of the HIV-1 broadly neutralizing antibody 10E8 MPER epitope on the membrane surface by gp41 transmembrane domain scaffolds.

The 10E8 antibody achieves near-pan neutralization of HIV-1 by targeting the remarkably conserved gp41 membrane-proximal external region (MPER) and the connected transmembrane domain (TMD) of the HIV-1 envelope glycoprotein (Env). Thus, recreating the structure that generates 10E8-like antibodies is a major goal of the rational design of anti-HIV vaccines. Unfortunately, high-resolution information of this segment in the native Env is lacking, limiting our understanding of the behavior of the crucial 10E8 epitope residues. In this report, two sequences, namely, MPER-TMD1 (gp41 residues 671-700) and MPER-TMD2 (gp41 residues 671-709) were compared both experimentally and computationally, to assess the TMD as a potential membrane integral scaffold for the 10E8 epitope. These sequences were selected to represent a minimal (MPER-TMD1) or full-length (MPER-TMD2) TMD membrane anchor according to mutagenesis results reported by Yue et al. (2009) J. Virol. 83, 11,588. Immunochemical assays revealed that MPER-TMD1, but not MPER-TMD2, effectively exposed the MPER C-terminal stretch, harboring the 10E8 epitope on the surface of phospholipid bilayers containing a cholesterol concentration equivalent to that of the viral envelope. Molecular dynamics simulations, using the recently resolved TMD trimer structure combined with the MPER in a cholesterol-enriched model membrane confirmed these results and provided an atomistic mechanism of epitope exposure which revealed that TMD truncation at position A700 combined with N-terminal addition of lysine residues positively impacts epitope exposure. Overall, these results provide crucial insights into the design of effective MPER-TMD derived immunogens.

Stereospecific Interactions of Cholesterol in a Model Cell Membrane: Implications for the Membrane Dipole Potential.

Cholesterol is a major constituent of the plasma membrane in higher order eukaryotes. The effect of cholesterol on the structure and organisation of cell membranes has been studied extensively by both experimental and computational means. In recent years, a wealth of data has been accumulated illustrating how subtle differences in the structure of cholesterol equate to considerable changes in the physical properties of the membrane. The effect of cholesterol stereoisomers, in particular, has been established, identifying a direct link with the activity of specific membrane proteins. In this study, we perform extensive molecular dynamics simulations of phospholipid bilayers containing three isomers of cholesterol, the native form (nat-cholesterol), the enantiomer of the native form (ent-cholesterol), and an epimer of cholesterol that differs by the orientation of the polar hydroxyl group (epi-cholesterol). Based on these simulations, an atomic-level description of the stereospecific cholesterol-phospholipid interactions is provided, establishing a potential mechanism for the perturbation of membrane properties, specifically the membrane dipole potential.

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology.

This is a perspective article entitled "Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology" which is following a CECAM meeting with the same name.

Location and Character of Volatile General Anesthetics Binding Sites in the Transmembrane Domain of TRPV1.

It has been proposed that general anesthesia results from direct multisite interactions with multiple and diverse ion channels in the brain. An understanding of the mechanisms by which general anesthetics modulate ion channels is essential to clarify their underlying behavior and their role in reversible immobilization and amnesia. Despite the fact that volatile general anesthetics are drugs that primarily induce insensitivity to pain, they have been reported to sensitize and active the vanilloid-1 receptor, TRPV1, which is known to mediate the response of the nervous system to certain harmful stimuli and which plays a crucial role in the pain pathway. Currently, the mechanism of action of anesthetics is unknown and the precise molecular sites of interaction have not been identified. Here, using ∼2.5 µs of classical molecular dynamics simulations and metadynamics, we explore these enigmas. Binding sites are identified and the strength of the association is further characterized using alchemical free-energy calculations. Anesthetic binding/unbinding proceeds primarily through a membrane-embedded pathway, and subsequently, a complex scenario is established involving multiple binding sites featuring single or multiple occupancy states of two small volatile drugs. One of the five anesthetic binding sites reported was previously identified experimentally, and another one, importantly, is identical to that of capsaicin, one of the chemical stimuli that activate TRPV1. However, in contrast to capsaicin, isoflurane and chloroform binding free-energies render modest to no association compared to capsaicin, suggesting a different activation mechanism. Uncovering chloroform and isoflurane modulatory sites will further our understanding of the TRPV1 molecular machinery and open the possibility of developing site-specific drugs.