An Atomistic Model for Simulations of the General Anesthetic Isoflurane (original) (raw)
Related papers
Atomistic models of general anesthetics for use in in silico biological studies
The journal of physical chemistry. B, 2014
While small molecules have been used to induce anesthesia in a clinical setting for well over a century, a detailed understanding of the molecular mechanism remains elusive. In this study, we utilize ab initio calculations to develop a novel set of CHARMM-compatible parameters for the ubiquitous modern anesthetics desflurane, isoflurane, sevoflurane, and propofol for use in molecular dynamics (MD) simulations. The parameters generated were rigorously tested against known experimental physicochemical properties including dipole moment, density, enthalpy of vaporization, and free energy of solvation. In all cases, the anesthetic parameters were able to reproduce experimental measurements, signifying the robustness and accuracy of the atomistic models developed. The models were then used to study the interaction of anesthetics with the membrane. Calculation of the potential of mean force for inserting the molecules into a POPC bilayer revealed a distinct energetic minimum of 4-5 kcal/m...
Conformational Exploration of Enflurane in Solution and in a Biological Environment
The journal of physical chemistry. A, 2015
Enflurane is a fluorinated volatile anesthetic, whose bioactive conformation is not known. Actually, a few studies have reported on the conformations of enflurane in nonpolar solution and gas phase. The present computational and spectroscopic (infrared and NMR) work shows that three pairs of isoenergetic conformers take place in the gas phase, neat liquid, polar, and nonpolar solutions. According to docking studies, a single conformation is largely preferred over its isoenergetic isomers to complex with the active site of Integrin LFA-1 enzyme (PDB code: 3F78 ), where the widely used anesthetic isoflurane (a constitutional isomer of enflurane) is known to bind. Weak hydrogen bonding from an electrostatic interaction between the CHF2 hydrogen and the central CF2 fluorines was not found to rule the conformational isomerism of enflurane. Moreover, intramolecular interactions based on steric, electrostatic, and hyperconjugative effects usually invoked to describe the anomeric effect are...
Interactions of Anesthetics with the Water−Hexane Interface. A Molecular Dynamics Study
The Journal of Physical Chemistry B, 1997
The free energy profiles characterizing the transfer of nine solutes across the liquid-vapor interfaces of water and hexane and across the water-hexane interface were calculated from molecular dynamics simulations. Among the solutes were n-butane and three of its halogenated derivatives, as well as three halogenated cyclobutanes. The two remaining molecules, dichlorodifluoromethane and 1,2-dichloroperfluoroethane, belong to series of halo-substituted methanes and ethanes, described in previous studies (J. Chem. Phys. 1996, 104, 3760; Chem. Phys. 1996, 204, 337). Each series of molecules contains structurally similar compounds that differ greatly in anesthetic potency. The accuracy of the simulations was tested by comparing the calculated and the experimental free energies of solvation of all nine compounds in water and in hexane. In addition, the calculated and the measured surface excess concentrations of n-butane at the water liquid-vapor interface were compared. In all cases, good agreement with experimental results was found. At the water-hexane interface, the free energy profiles for polar molecules exhibited significant interfacial minima, whereas the profiles for nonpolar molecules did not. The existence of these minima was interpreted in terms of a balance between the free energy contribution arising from solute-solvent interactions and the work to form a cavity that accommodates the solute. These two contributions change monotonically, but oppositely, across the interface. The interfacial solubilities of the solutes, obtained from the free energy profiles, correlate very well with their anesthetic potencies. This is the case even when the Meyer-Overton hypothesis, which predicts a correlation between anesthetic potency and solubility in oil, fails.
Water solvent and local anesthetics: A computational study
International Journal of Quantum Chemistry, 2007
There are various experimental studies regarding the toxicity and the time of action of local anesthetics, which contain general insights about their pharmacological and physicochemical properties. Although a detailed microscopic analysis of the local anesthetics would contribute to understanding these properties, there are relatively few theoretical studies about these molecules. In this article, we present the results from calculations performed for three local anesthetics: tetracaine, procaine, and lidocaine, both in their charged and uncharged forms, in aqueous environment. We have used the density functional theory and molecular dynamics simulations to study the structural characteristics of these compounds. The radial distribution function g(r) was used to examine the structure of water molecules surrounding different regions of the local anesthetics. We demonstrated the nonhomogeneous character of the anesthetics with respect to their affinity to water solvent molecules as well as the modifications in their affinity to water caused by changes in their charge state. We also observed that the biological potency of the anesthetics is more related to the behavior of specific groups within the molecule, which are responsible for the interaction with the lipid phase of membranes, rather than the general properties of the molecule as a whole.
Interactions of anesthetics with the membrane-water interface
Chemical Physics, 1996
Although the potency of conventional anesthetics correlates with lipophilicity, an affinity to water also is essential. It was recently found that compounds with very low affinities to water do not produce anesthesia regardless of their lipophilicity. This finding implies that clinical anesthesia might arise because of interactions at molecular sites near the interface of neuronal membranes with the aqueous environment and, therefore, might require increased concentrations of anesthetic molecules at membrane interfaces. As an initial test of this hypothesis, we calculated in molecular dynamics simulations the free energy profiles for the transfer of anesthetic l,l,Ztrifluoroethane and nonanesthetic perlluoroethane across water-membrane and water-hexane interfaces. Consistent with the hypothesis, it was found that trifluoroethane, but not perhuoroethane, exhibits a free energy minimum and, therefore, increased concentrations at both interfaces. The transfer of trifluoroethane from water to the nonpolar hexane or interior of the membrane is accompanied by a considerable, solvent-induced shift in the conformational equilibrium around the C-C bond. a single solute molecule and a lamella of 495 water However, at levels of the agent lower than 15% the measured MAC becomes unreliable.
Atomistic Model for Simulations of the Sedative Hypnotic Drug 2,2,2-Trichloroethanol
ACS Omega, 2018
2,2,2-Trichloroethanol (TCE) is the active form of the sedative hypnotic drug chloral hydrate, one of the oldest sleep medications in the market. Understanding of TCE's action mechanisms to its many targets, particularly within the ion channel family, could benefit from the state-ofthe-art computational molecular studies. In this direction, we employed de novo modeling aided by the force field toolkit to develop CHARMM36-compatible TCE parameters. The classical potential energy function was calibrated targeting molecular conformations, local interactions with water molecules, and liquid bulk properties. Reference data comes from both tabulated thermodynamic properties and ab initio calculations at the MP2 level. TCE solvation free energy calculations in water and oil reproduce a lipophilic, yet nonhydrophobic, behavior. Indeed, the potential mean force profile for TCE partition through the phospholipid bilayer reveals the sedative's preference for the interfacial region. The calculated partition coefficient also matches experimental measures. Further validation of the proposed parameters is supported by the model's ability to recapitulate quenching experiments demonstrating TCE binding to bovine serum albumin.
Isoflurane Alters the Structure and Dynamics of GLIC
Biophysical Journal, 2011
Pentameric ligand-gated ion channels are targets of general anesthetics. Although the search for discrete anesthetic binding sites has achieved some degree of success, little is known regarding how anesthetics work after the events of binding. Using the crystal structures of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), which is sensitive to a variety of general anesthetics, we performed multiple molecular dynamics simulations in the presence and absence of the general anesthetic isoflurane. Isoflurane bound to several locations within GLIC, including the transmembrane pocket identified crystallographically, the extracellular (EC) domain, and the interface of the EC and transmembrane domains. Isoflurane also entered the channel after the pore was dehydrated in one of the simulations. Isoflurane disrupted the quaternary structure of GLIC, as evidenced in a striking association between the binding and breakage of intersubunit salt bridges in the EC domain. The pore-lining helix experienced lateral and inward radial tilting motion that contributed to the channel closure. Isoflurane binding introduced strong anticorrelated motions between different subunits of GLIC. The demonstrated structural and dynamical modulations by isoflurane aid in the understanding of the underlying mechanism of anesthetic inhibition of GLIC and possibly other homologous pentameric ligand-gated ion channels.
Physical Chemistry Chemical Physics, 2012
Propofol is a general anaesthetic that exerts its action by interaction with the GABA A receptor. Crystallographic studies suggest that there is a direct interaction between propofol and the phenolic residue of a tyrosine in the channel. In this study we create propofolÁ Á Áphenol clusters by their co-expansion in jets. The complex is probed using a set of mass-resolved spectroscopic strategies: 2-color REMPI, UV/UV hole-burning, IR/UV double resonance and the novel technique IR/IR/UV triple resonance. The existence of at least six different isomers in the expansion is demonstrated. All the isomers are stabilized by interactions between their aromatic rings. Additionally, in some conformers the OH moieties form hydrogen bonds in some of the isomers, with propofol and phenol alternating their donor-acceptor roles, while in others the -OHÁ Á ÁOH angle points to a dipole-dipole interaction. Interpretation of the data in the light of dispersion-corrected DFT calculations shows that shallow barriers separate all the isomers, both in the ground and excited electronic states. Comparison of the structures of the complex with the X-ray diffraction data is also offered.
Computational Studies on the Interactions of Inhalational Anesthetics with Proteins
Accounts of Chemical Research, 2010
D espite the widespread clinical use of anesthetics since the 19th century, a clear understanding of the mechanism of anesthetic action has yet to emerge. On the basis of early experiments by Meyer, Overton, and subsequent researchers, the cell's lipid membrane was generally concluded to be the primary site of action of anesthetics. However, later experiments with lipid-free globular proteins, such as luciferase and apoferritin, shifted the focus of anesthetic action to proteins. Recent experimental studies, such as photoaffinity labeling and mutagenesis on membrane proteins, have suggested specific binding sites for anesthetic molecules, further strengthening the proteocentric view of anesthetic mechanism. With the increased availability of high-resolution crystal structures of ion channels and other integral membrane proteins, as well as the availability of powerful computers, the structure-function relationship of anesthetic-protein interactions can now be investigated in atomic detail.