Simulating the Skin Permeation Process of Ionizable Molecules.

It is commonly assumed that ionizable molecules, such as drugs, permeate through the skin barrier in their neutral form. By using molecular dynamics simulations of the charged and neutral states separately, we can study the dynamic protonation behavior during the permeation process. We have studied three weak acids and three weak bases and conclude that the acids are ionized to a larger extent than the bases, when passing through the headgroup region of the lipid barrier structure, at pH values close to their pKa. It can also be observed that even if these dynamic protonation simulations are informative, in the cases studied herein they are not necessary for the calculation of permeability coefficients. It is sufficient to base the calculations only on the neutral form, as is commonly done.

Allosteric Cholesterol Site in Glycine Receptors Characterized through Molecular Simulations.

Glycine receptors are pentameric ligand-gated ion channels that conduct chloride ions across postsynaptic membranes to facilitate fast inhibitory neurotransmission. In addition to gating by the glycine agonist, interactions with lipids and other compounds in the surrounding membrane environment modulate their function, but molecular details of these interactions remain unclear, in particular, for cholesterol. Here, we report coarse-grained simulations in a model neuronal membrane for three zebrafish glycine receptor structures representing apparent resting, open, and desensitized states. We then converted the systems to all-atom models to examine detailed lipid interactions. Cholesterol bound to the receptor at an outer-leaflet intersubunit site, with a preference for the open and desensitized versus resting states, indicating that it can bias receptor function. Finally, we used short atomistic simulations and iterative amino acid perturbations to identify residues that may mediate allosteric gating transitions. Frequent cholesterol contacts in atomistic simulations clustered with residues identified by perturbation analysis and overlapped with mutations influencing channel function and pathology. Cholesterol binding at this site was also observed in a recently reported pig heteromeric glycine receptor. These results indicate state-dependent lipid interactions relevant to allosteric transitions of glycine receptors, including specific amino acid contacts applicable to biophysical modeling and pharmaceutical design.

Functional and structural insights into activation of TRPV2 by weak acids.

Transient receptor potential (TRP) ion channels are involved in the surveillance or regulation of the acid-base balance. Here, we demonstrate that weak carbonic acids, including acetic acid, lactic acid, and CO2 activate and sensitize TRPV2 through a mechanism requiring permeation through the cell membrane. TRPV2 channels in cell-free inside-out patches maintain weak acid-sensitivity, but protons applied on either side of the membrane do not induce channel activation or sensitization. The involvement of proton modulation sites for weak acid-sensitivity was supported by the identification of titratable extracellular (Glu495, Glu561) and intracellular (His521) residues on a cryo-EM structure of rat TRPV2 (rTRPV2) treated with acetic acid. Molecular dynamics simulations as well as patch clamp experiments on mutant rTRPV2 constructs confirmed that these residues are critical for weak acid-sensitivity. We also demonstrate that the pore residue Glu609 dictates an inhibition of weak acid-induced currents by extracellular calcium. Finally, TRPV2-expression in HEK293 cells is associated with an increased weak acid-induced cytotoxicity. Together, our data provide new insights into weak acids as endogenous modulators of TRPV2.

Interactive computational and experimental approaches improve the sensitivity of periplasmic binding protein-based nicotine biosensors for measurements in biofluids.

We developed fluorescent protein sensors for nicotine with improved sensitivity. For iNicSnFR12 at pH 7.4, the proportionality constant for ∆F/F0vs [nicotine] (δ-slope, 2.7 µM-1) is 6.1-fold higher than the previously reported iNicSnFR3a. The activated state of iNicSnFR12 has a fluorescence quantum yield of at least 0.6. We measured similar dose-response relations for the nicotine-induced absorbance increase and fluorescence increase, suggesting that the absorbance increase leads to the fluorescence increase via the previously described nicotine-induced conformational change, the 'candle snuffer' mechanism. Molecular dynamics (MD) simulations identified a binding pose for nicotine, previously indeterminate from experimental data. MD simulations also showed that Helix 4 of the periplasmic binding protein (PBP) domain appears tilted in iNicSnFR12 relative to iNicSnFR3a, likely altering allosteric network(s) that link the ligand binding site to the fluorophore. In thermal melt experiments, nicotine stabilized the PBP of the tested iNicSnFR variants. iNicSnFR12 resolved nicotine in diluted mouse and human serum at 100 nM, the peak [nicotine] that occurs during smoking or vaping, and possibly at the decreasing levels during intervals between sessions. NicSnFR12 was also partially activated by unidentified endogenous ligand(s) in biofluids. Improved iNicSnFR12 variants could become the molecular sensors in continuous nicotine monitors for animal and human biofluids.

Structural basis for the oligomerization-facilitated NLRP3 activation.

The NACHT-, leucine-rich-repeat-, and pyrin domain-containing protein 3 (NLRP3) is a critical intracellular inflammasome sensor and an important clinical target against inflammation-driven human diseases. Recent studies have elucidated its transition from a closed cage to an activated disk-like inflammasome, but the intermediate activation mechanism remains elusive. Here we report the cryo-electron microscopy structure of NLRP3, which forms an open octamer and undergoes a ~ 90° hinge rotation at the NACHT domain. Mutations on open octamer's interfaces reduce IL-1ß signaling, highlighting its essential role in NLRP3 activation/inflammasome assembly. The centrosomal NIMA-related kinase 7 (NEK7) disrupts large NLRP3 oligomers and forms NEK7/NLRP3 monomers/dimers which is a critical step preceding the assembly of the disk-like inflammasome. These data demonstrate an oligomeric cooperative activation of NLRP3 and provide insight into its inflammasome assembly mechanism.

Discovery of lipid binding sites in a ligand-gated ion channel by integrating simulations and cryo-EM.

Ligand-gated ion channels transduce electrochemical signals in neurons and other excitable cells. Aside from canonical ligands, phospholipids are thought to bind specifically to the transmembrane domain of several ion channels. However, structural details of such lipid contacts remain elusive, partly due to limited resolution of these regions in experimental structures. Here, we discovered multiple lipid interactions in the channel GLIC by integrating cryo-electron microscopy and large-scale molecular simulations. We identified 25 bound lipids in the GLIC closed state, a conformation where none, to our knowledge, were previously known. Three lipids were associated with each subunit in the inner leaflet, including a buried interaction disrupted in mutant simulations. In the outer leaflet, two intrasubunit sites were evident in both closed and open states, while a putative intersubunit site was preferred in open-state simulations. This work offers molecular details of GLIC-lipid contacts particularly in the ill-characterized closed state, testable hypotheses for state-dependent binding, and a multidisciplinary strategy for modeling protein-lipid interactions.

Interactive computational and experimental approaches improve the sensitivity of periplasmic binding protein-based nicotine biosensors for measurements in biofluids.

We developed fluorescent protein sensors for nicotine with improved sensitivity. For iNicSnFR12 at pH 7.4, the proportionality constant for ΔF/F0 vs [nicotine] (δ-slope, 2.7 µM-1) is 6.1-fold higher than the previously reported iNicSnFR3a. The activated state of iNicSnFR12 has a fluorescence quantum yield of at least 0.6. We measured similar dose-response relations for the nicotine-induced absorbance increase and fluorescence increase, suggesting that the absorbance increase leads to the fluorescence increase via the previously described nicotine-induced conformational change, the "candle snuffer" mechanism. Molecular dynamics (MD) simulations identified a binding pose for nicotine, previously indeterminate from experimental data. MD simulations also showed that Helix 4 of the periplasmic binding protein (PBP) domain appears tilted in iNicSnFR12 relative to iNicSnFR3a, likely altering allosteric network(s) that link the ligand binding site to the fluorophore. In thermal melt experiments, nicotine stabilized the PBP of the tested iNicSnFR variants. iNicSnFR12 resolved nicotine in diluted mouse and human serum at 100 nM, the peak [nicotine] that occurs during smoking or vaping, and possibly at the decreasing levels during intervals between sessions. NicSnFR12 was also partially activated by unidentified endogenous ligand(s) in biofluids. Improved iNicSnFR12 variants could become the molecular sensors in continuous nicotine monitors for animal and human biofluids.

Structure and dynamics of differential ligand binding in the human ρ-type GABAA receptor.

The neurotransmitter γ-aminobutyric acid (GABA) drives critical inhibitory processes in and beyond the nervous system, partly via ionotropic type-A receptors (GABAARs). Pharmacological properties of ρ-type GABAARs are particularly distinctive, yet the structural basis for their specialization remains unclear. Here, we present cryo-EM structures of a lipid-embedded human ρ1 GABAAR, including a partial intracellular domain, under apo, inhibited, and desensitized conditions. An apparent resting state, determined first in the absence of modulators, was recapitulated with the specific inhibitor (1,2,5,6-tetrahydropyridin-4-yl)methylphosphinic acid and blocker picrotoxin and provided a rationale for bicuculline insensitivity. Comparative structures, mutant recordings, and molecular simulations with and without GABA further explained the sensitized but slower activation of ρ1 relative to canonical subtypes. Combining GABA with picrotoxin also captured an apparent uncoupled intermediate state. This work reveals structural mechanisms of gating and modulation with applications to ρ-specific pharmaceutical design and to our biophysical understanding of ligand-gated ion channels.

Structural insights into opposing actions of neurosteroids on GABAA receptors.

γ-Aminobutyric acid type A (GABAA) receptors mediate fast inhibitory signaling in the brain and are targets of numerous drugs and endogenous neurosteroids. A subset of neurosteroids are GABAA receptor positive allosteric modulators; one of these, allopregnanolone, is the only drug approved specifically for treating postpartum depression. There is a consensus emerging from structural, physiological and photolabeling studies as to where positive modulators bind, but how they potentiate GABA activation remains unclear. Other neurosteroids are negative modulators of GABAA receptors, but their binding sites remain debated. Here we present structures of a synaptic GABAA receptor bound to allopregnanolone and two inhibitory sulfated neurosteroids. Allopregnanolone binds at the receptor-bilayer interface, in the consensus potentiator site. In contrast, inhibitory neurosteroids bind in the pore. MD simulations and electrophysiology support a mechanism by which allopregnanolone potentiates channel activity and suggest the dominant mechanism for sulfated neurosteroid inhibition is through pore block.

Understanding Drug Skin Permeation Enhancers Using Molecular Dynamics Simulations.

Our skin constitutes an effective permeability barrier that protects the body from exogenous substances but concomitantly severely limits the number of pharmaceutical drugs that can be delivered transdermally. In topical formulation design, chemical permeation enhancers (PEs) are used to increase drug skin permeability. In vitro skin permeability experiments can measure net effects of PEs on transdermal drug transport, but they cannot explain the molecular mechanisms of interactions between drugs, permeation enhancers, and skin structure, which limits the possibility to rationally design better new drug formulations. Here we investigate the effect of the PEs water, lauric acid, geraniol, stearic acid, thymol, ethanol, oleic acid, and eucalyptol on the transdermal transport of metronidazole, caffeine, and naproxen. We use atomistic molecular dynamics (MD) simulations in combination with developed molecular models to calculate the free energy difference between 11 PE-containing formulations and the skin's barrier structure. We then utilize the results to calculate the final concentration of PEs in skin. We obtain an RMSE of 0.58 log units for calculated partition coefficients from water into the barrier structure. We then use the modified PE-containing barrier structure to calculate the PEs' permeability enhancement ratios (ERs) on transdermal metronidazole, caffeine, and naproxen transport and compare with the results obtained from in vitro experiments. We show that MD simulations are able to reproduce rankings based on ERs. However, strict quantitative correlation with experimental data needs further refinement, which is complicated by significant deviations between different measurements. Finally, we propose a model for how to use calculations of the potential of mean force of drugs across the skin's barrier structure in a topical formulation design.

Gentle and fast all-atom model refinement to cryo-EM densities via a maximum likelihood approach.

Better detectors and automated data collection have generated a flood of high-resolution cryo-EM maps, which in turn has renewed interest in improving methods for determining structure models corresponding to these maps. However, automatically fitting atoms to densities becomes difficult as their resolution increases and the refinement potential has a vast number of local minima. In practice, the problem becomes even more complex when one also wants to achieve a balance between a good fit of atom positions to the map, while also establishing good stereochemistry or allowing protein secondary structure to change during fitting. Here, we present a solution to this challenge using a maximum likelihood approach by formulating the problem as identifying the structure most likely to have produced the observed density map. This allows us to derive new types of smooth refinement potential-based on relative entropy-in combination with a novel adaptive force scaling algorithm to allow balancing of force-field and density-based potentials. In a low-noise scenario, as expected from modern cryo-EM data, the relative-entropy based refinement potential outperforms alternatives, and the adaptive force scaling appears to aid all existing refinement potentials. The method is available as a component in the GROMACS molecular simulation toolkit.

Automated simulation-based membrane protein refinement into cryo-EM data.

The resolution revolution has increasingly enabled single-particle cryogenic electron microscopy (cryo-EM) reconstructions of previously inaccessible systems, including membrane proteins-a category that constitutes a disproportionate share of drug targets. We present a protocol for using density-guided molecular dynamics simulations to automatically refine atomistic models into membrane protein cryo-EM maps. Using adaptive force density-guided simulations as implemented in the GROMACS molecular dynamics package, we show how automated model refinement of a membrane protein is achieved without the need to manually tune the fitting force ad hoc. We also present selection criteria to choose the best-fit model that balances stereochemistry and goodness of fit. The proposed protocol was used to refine models into a new cryo-EM density of the membrane protein maltoporin, either in a lipid bilayer or detergent micelle, and we found that results do not substantially differ from fitting in solution. Fitted structures satisfied classical model-quality metrics and improved the quality and the model-to-map correlation of the x-ray starting structure. Additionally, the density-guided fitting in combination with generalized orientation-dependent all-atom potential was used to correct the pixel-size estimation of the experimental cryo-EM density map. This work demonstrates the applicability of a straightforward automated approach to fitting membrane protein cryo-EM densities. Such computational approaches promise to facilitate rapid refinement of proteins under different conditions or with various ligands present, including targets in the highly relevant superfamily of membrane proteins.

MDverse: Shedding Light on the Dark Matter of Molecular Dynamics Simulations.

The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2,000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation, such as temperature and simulation length, and identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore collected MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and increase populating and standardizing metadata to reuse this valuable matter.

Structural Elucidation of Ivermectin Binding to α7nAChR and the Induced Channel Desensitization.

The α7 nicotinic acetylcholine receptor (α7nAChR) mediates signaling in the central nervous system and cholinergic anti-inflammatory pathways. Ivermectin is a positive allosteric modulator of a full-length α7nAChR and an agonist of the α7nAChR construct containing transmembrane (TMD) and intracellular (ICD) domains, but structural insights of the binding have not previously been determined. Here, combining nuclear magnetic resonance as a primary experimental tool with Rosetta comparative modeling and molecular dynamics simulations, we have revealed details of ivermectin binding to the α7nAChR TMD + ICD and corresponding structural changes in an ivermectin-induced desensitized state. Ivermectin binding was stabilized predominantly by hydrophobic interactions from interfacial residues between adjacent subunits near the extracellular end of the TMD, where the inter-subunit gap was substantially expanded in comparison to the apo structure. The ion-permeation pathway showed a profile distinctly different from the resting-state profile but similar to profiles of desensitized α7nAChR. The ICD also exhibited structural changes, including reorientation of the MX and h3 helices relative to the channel axis. The resulting structures of the α7nAChR TMD + ICD in complex with ivermectin provide opportunities for discovering new modulators of therapeutic potential and exploring the structural basis of cytoplasmic signaling under different α7nAChR functional states.

Biophysical characterization of calcium-binding and modulatory-domain dynamics in a pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels (pLGICs) perform electrochemical signal transduction in organisms ranging from bacteria to humans. Among the prokaryotic pLGICs, there is architectural diversity involving N-terminal domains (NTDs) not found in eukaryotic relatives, exemplified by the calcium-sensitive channel (DeCLIC) from a Desulfofustis deltaproteobacterium, which has an NTD in addition to the canonical pLGIC structure. Here, we have characterized the structure and dynamics of DeCLIC through cryoelectron microscopy (cryo-EM), small-angle neutron scattering (SANS), and molecular dynamics (MD) simulations. In the presence and absence of calcium, cryo-EM yielded structures with alternative conformations of the calcium-binding site. SANS profiles further revealed conformational diversity at room temperature beyond that observed in static structures, shown through MD to be largely attributable to rigid-body motions of the NTD relative to the protein core, with expanded and asymmetric conformations improving the fit of the SANS data. This work reveals the range of motion available to the DeCLIC NTD and calcium-binding site, expanding the conformational landscape of the pLGIC family. Further, these findings demonstrate the power of combining low-resolution scattering, high-resolution structural, and MD simulation data to elucidate interfacial interactions that are highly conserved in the pLGIC family.

Differential interactions of resting, activated, and desensitized states of the α7 nicotinic acetylcholine receptor with lipidic modulators.

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that modulates neuronal excitability, largely by allowing Ca2+ permeation. Agonist binding promotes transition from a resting state to an activated state, and then rapidly to a desensitized state. Recently, cryogenic electron microscopy (cryo-EM) structures of the human α7 receptor in nanodiscs were reported in multiple conformations. These were selectively stabilized by inhibitory, activating, or potentiating compounds. However, the functional annotation of these structures and their differential interactions with unresolved lipids and ligands remain incomplete. Here, we characterized their ion permeation, membrane interactions, and ligand binding using computational electrophysiology, free-energy calculations, and coarse-grained molecular dynamics. In contrast to nonconductive structures in apparent resting and desensitized states, the structure determined in the presence of the potentiator PNU-120596 was consistent with an activated state permeable to Ca2+. Transition to this state was associated with compression and rearrangement of the membrane, particularly in the vicinity of the peripheral MX helix. An intersubunit transmembrane site was implicated in selective binding of either PNU-120596 in the activated state or cholesterol in the desensitized state. This substantiates functional assignment of all three lipid-embedded α7-receptor structures with ion-permeation simulations. It also proposes testable models of their state-dependent interactions with lipophilic ligands, including a mechanism for allosteric modulation at the transmembrane subunit interface.

Skin permeability prediction with MD simulation sampling spatial and alchemical reaction coordinates.

A molecular-level understanding of skin permeation may rationalize and streamline product development, and improve quality and control, of transdermal and topical drug delivery systems. It may also facilitate toxicity and safety assessment of cosmetics and skin care products. Here, we present new molecular dynamics simulation approaches that make it possible to efficiently sample the free energy and local diffusion coefficient across the skin's barrier structure to predict skin permeability and the effects of chemical penetration enhancers. In particular, we introduce a new approach to use two-dimensional reaction coordinates in the accelerated weight histogram method, where we combine sampling along spatial coordinates with an alchemical perturbation virtual coordinate. We present predicted properties for 20 permeants, and demonstrate how our approach improves correlation with ex vivo/in vitro skin permeation data. For the compounds included in this study, the obtained log KPexp-calc mean square difference was 0.9 cm2 h-2.

Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case.

Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.

Structural and dynamic mechanisms of GABAA receptor modulators with opposing activities.

γ-Aminobutyric acid type A (GABAA) receptors are pentameric ligand-gated ion channels abundant in the central nervous system and are prolific drug targets for treating anxiety, sleep disorders and epilepsy. Diverse small molecules exert a spectrum of effects on γ-aminobutyric acid type A (GABAA) receptors by acting at the classical benzodiazepine site. They can potentiate the response to GABA, attenuate channel activity, or counteract modulation by other ligands. Structural mechanisms underlying the actions of these drugs are not fully understood. Here we present two high-resolution structures of GABAA receptors in complex with zolpidem, a positive allosteric modulator and heavily prescribed hypnotic, and DMCM, a negative allosteric modulator with convulsant and anxiogenic properties. These two drugs share the extracellular benzodiazepine site at the α/γ subunit interface and two transmembrane sites at ß/α interfaces. Structural analyses reveal a basis for the subtype selectivity of zolpidem that underlies its clinical success. Molecular dynamics simulations provide insight into how DMCM switches from a negative to a positive modulator as a function of binding site occupancy. Together, these findings expand our understanding of how GABAA receptor allosteric modulators acting through a common site can have diverging activities.

A general tight-binding based energy decomposition analysis scheme for intermolecular interactions in large molecules.

In this work, a general tight-binding based energy decomposition analysis (EDA) scheme for intermolecular interactions is proposed. Different from the earlier version [Xu et al., J. Chem. Phys. 154, 194106 (2021)], the current tight-binding based density functional theory (DFTB)-EDA is capable of performing interaction analysis with all the self-consistent charge (SCC) type DFTB methods, including SCC-DFTB2/3 and GFN1/2-xTB, despite their different formulas and parameterization schemes. In DFTB-EDA, the total interaction energy is divided into frozen, polarization, and dispersion terms. The performance of DFTB-EDA with SCC-DFTB2/3 and GFN1/2-xTB for various interaction systems is discussed and assessed.