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.

Shared structural mechanisms of general anaesthetics and benzodiazepines.

Most general anaesthetics and classical benzodiazepine drugs act through positive modulation of γ-aminobutyric acid type A (GABAA) receptors to dampen neuronal activity in the brain1-5. However, direct structural information on the mechanisms of general anaesthetics at their physiological receptor sites is lacking. Here we present cryo-electron microscopy structures of GABAA receptors bound to intravenous anaesthetics, benzodiazepines and inhibitory modulators. These structures were solved in a lipidic environment and are complemented by electrophysiology and molecular dynamics simulations. Structures of GABAA receptors in complex with the anaesthetics phenobarbital, etomidate and propofol reveal both distinct and common transmembrane binding sites, which are shared in part by the benzodiazepine drug diazepam. Structures in which GABAA receptors are bound by benzodiazepine-site ligands identify an additional membrane binding site for diazepam and suggest an allosteric mechanism for anaesthetic reversal by flumazenil. This study provides a foundation for understanding how pharmacologically diverse and clinically essential drugs act through overlapping and distinct mechanisms to potentiate inhibitory signalling in the brain.

Bacterial protein domains with a novel Ig-like fold target human CEACAM receptors.

Streptococcus agalactiae, also known as group B Streptococcus (GBS), is the major cause of neonatal sepsis in humans. A critical step to infection is adhesion of bacteria to epithelial surfaces. GBS adhesins have been identified to bind extracellular matrix components and cellular receptors. However, several putative adhesins have no host binding partner characterised. We report here that surface-expressed ß protein of GBS binds to human CEACAM1 and CEACAM5 receptors. A crystal structure of the complex showed that an IgSF domain in ß represents a novel Ig-fold subtype called IgI3, in which unique features allow binding to CEACAM1. Bioinformatic assessment revealed that this newly identified IgI3 fold is not exclusively present in GBS but is predicted to be present in adhesins from other clinically important human pathogens. In agreement with this prediction, we found that CEACAM1 binds to an IgI3 domain found in an adhesin from a different streptococcal species. Overall, our results indicate that the IgI3 fold could provide a broadly applied mechanism for bacteria to target CEACAMs.

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.

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.

Exploring the Impact of Protein Chain Selection in Binding Energy Calculations with DFT.

Calculation of binding free energies between a protein and a ligand are highly desired for computer-aided drug design. Here we approximate the binding energies of ABL1, an enzyme which is the target for drugs used in the treatment of chronic myeloid leukaemia, with minimal models and density functional theory (DFT). Starting from the crystal structures of protein-drug complexes, we estimated the binding free energies having used all available individual molecules (protein chains) within each structure, not only a single one as commonly used, in order to see if the choice of the protein chain is important in such calculations. Differences were observed between chains in the same file. Energy decomposition analysis (EDA) revealed that the most important factors for binding were exchange, repulsion and electrostatics. The desolvation term varied dramatically between the inhibitors (between 4.2 and 92.3 kcal/mol). All functionals showed similar patterns in the EDA and in discriminating between the ligands. Non-covalent interactions (NCI) analysis was used to further explain the differences between protein chains and functionals. Overall, it is shown that small minimal models of a drug binding site can be useful to infer on the suitability of an initial crystal structure for further analysis such as EDA.

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.

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.

Trans vs. cis: a computational study of enasidenib resistance due to IDH2 mutations.

Isocitrate dehydrogenase 2 (IDH2) is a homodimeric enzyme that plays an important role in energy production. A mutation R140Q in one monomer makes the enzyme tumourigenic. Enasidenib is an effective inhibitor of IDH2/R140Q. A secondary mutation Q316E leads to enasidenib resistance. This mutation was hitherto only found in trans, i.e. where one monomer has the R140Q mutation and the other carries the Q316E mutation. It is not clear if the mutation only leads to resistance when in trans or if it has been discovered in trans only by chance, since it was only reported in two patients. Using molecular dynamics (MD) simulations we show that the binding of enasidenib to IDH2 is indeed much weaker when the Q316E mutation takes place in trans not in cis, which provides a molecular explanation for the clinical finding. This is corroborated by non-covalent interaction (NCI) analysis and DFT calculations. Whereas the MD simulations show a loss of one hydrogen bond upon the resistance mutation, NCI and energy decomposition analysis (EDA) reveal that a multitude of interactions are weakened.

Probing solution structure of the pentameric ligand-gated ion channel GLIC by small-angle neutron scattering.

Pentameric ligand-gated ion channels undergo subtle conformational cycling to control electrochemical signal transduction in many kingdoms of life. Several crystal structures have now been reported in this family, but the functional relevance of such models remains unclear. Here, we used small-angle neutron scattering (SANS) to probe ambient solution-phase properties of the pH-gated bacterial ion channel GLIC under resting and activating conditions. Data collection was optimized by inline paused-flow size-exclusion chromatography, and exchanging into deuterated detergent to hide the micelle contribution. Resting-state GLIC was the best-fit crystal structure to SANS curves, with no evidence for divergent mechanisms. Moreover, enhanced-sampling molecular-dynamics simulations enabled differential modeling in resting versus activating conditions, with the latter corresponding to an intermediate ensemble of both the extracellular and transmembrane domains. This work demonstrates state-dependent changes in a pentameric ion channel by SANS, an increasingly accessible method for macromolecular characterization with the coming generation of neutron sources.

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.

Allosteric enhancement of the BCR-Abl1 kinase inhibition activity of nilotinib by cobinding of asciminib.

Inhibitors that bind competitively to the ATP binding pocket in the kinase domain of the oncogenic fusion protein BCR-Abl1 are used successfully in targeted therapy of chronic myeloid leukemia (CML). Such inhibitors provided the first proof of concept that kinase inhibition can succeed in a clinical setting. However, emergence of drug resistance and dose-dependent toxicities limit the effectiveness of these drugs. Therefore, treatment with a combination of drugs without overlapping resistance mechanisms appears to be an appropriate strategy. In the present work, we explore the effectiveness of combination therapies of the recently developed allosteric inhibitor asciminib with the ATP-competitive inhibitors nilotinib and dasatinib in inhibiting the BCR-Abl1 kinase activity in CML cell lines. Through these experiments, we demonstrate that asciminib significantly enhances the inhibition activity of nilotinib, but not of dasatinib. Exploring molecular mechanisms for such allosteric enhancement via systematic computational investigation incorporating molecular dynamics, metadynamics simulations, and density functional theory calculations, we found two distinct contributions. First, binding of asciminib triggers conformational changes in the inactive state of the protein, thereby making the activation process less favorable by ∼4 kcal/mol. Second, the binding of asciminib decreases the binding free energies of nilotinib by ∼3 and ∼7 kcal/mol for the wildtype and T315I-mutated protein, respectively, suggesting the possibility of reducing nilotinib dosage and lowering risk of developing resistance in the treatment of CML.

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.

Structural basis for allosteric transitions of a multidomain pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels (pLGICs) are allosteric receptors that mediate rapid electrochemical signal transduction in the animal nervous system through the opening of an ion pore upon binding of neurotransmitters. Orthologs have been found and characterized in prokaryotes and they display highly similar structure-function relationships to eukaryotic pLGICs; however, they often encode greater architectural diversity involving additional amino-terminal domains (NTDs). Here we report structural, functional, and normal-mode analysis of two conformational states of a multidomain pLGIC, called DeCLIC, from a Desulfofustis deltaproteobacterium, including a periplasmic NTD fused to the conventional ligand-binding domain (LBD). X-ray structure determination revealed an NTD consisting of two jelly-roll domains interacting across each subunit interface. Binding of Ca2+ at the LBD subunit interface was associated with a closed transmembrane pore, with resolved monovalent cations intracellular to the hydrophobic gate. Accordingly, DeCLIC-injected oocytes conducted currents only upon depletion of extracellular Ca2+; these were insensitive to quaternary ammonium block. Furthermore, DeCLIC crystallized in the absence of Ca2+ with a wide-open pore and remodeled periplasmic domains, including increased contacts between the NTD and classic LBD agonist-binding sites. Functional, structural, and dynamical properties of DeCLIC paralleled those of sTeLIC, a pLGIC from another symbiotic prokaryote. Based on these DeCLIC structures, we would reclassify the previous structure of bacterial ELIC (the first high-resolution structure of a pLGIC) as a "locally closed" conformation. Taken together, structures of DeCLIC in multiple conformations illustrate dramatic conformational state transitions and diverse regulatory mechanisms available to ion channels in pLGICs, particularly involving Ca2+ modulation and periplasmic NTDs.

An angular motion of a conserved four-helix bundle facilitates alternating access transport in the TtNapA and EcNhaA transporters.

There is ongoing debate regarding the mechanism through which cation/proton antiporters (CPAs), like Thermus thermophilus NapA (TtNapA) and Escherichia coli NapA (EcNhaA), alternate between their outward- and inward-facing conformations in the membrane. CPAs comprise two domains, and it is unclear whether the transition is driven by their rocking-bundle or elevator motion with respect to each other. Here we address this question using metadynamics simulations of TtNapA, where we bias conformational sampling along two axes characterizing the two proposed mechanisms: angular and translational motions, respectively. By applying the bias potential for the two axes simultaneously, as well as to the angular, but not the translational, axis alone, we manage to reproduce each of the two known states of TtNapA when starting from the opposite state, in support of the rocking-bundle mechanism as the driver of conformational change. Next, starting from the inward-facing conformation of EcNhaA, we sample what could be its long-sought-after outward-facing conformation and verify it using cross-linking experiments.

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.

Regulation of a pentameric ligand-gated ion channel by a semiconserved cationic lipid-binding site.

Pentameric ligand-gated ion channels (pLGICs) are crucial mediators of electrochemical signal transduction in various organisms from bacteria to humans. Lipids play an important role in regulating pLGIC function, yet the structural bases for specific pLGIC-lipid interactions remain poorly understood. The bacterial channel ELIC recapitulates several properties of eukaryotic pLGICs, including activation by the neurotransmitter GABA and binding and modulation by lipids, offering a simplified model system for structure-function relationship studies. In this study, functional effects of noncanonical amino acid substitution of a potential lipid-interacting residue (W206) at the top of the M1-helix, combined with detergent interactions observed in recent X-ray structures, are consistent with this region being the location of a lipid-binding site on the outward face of the ELIC transmembrane domain. Coarse-grained and atomistic molecular dynamics simulations revealed preferential binding of lipids containing a positive charge, particularly involving interactions with residue W206, consistent with cation-π binding. Polar contacts from other regions of the protein, particularly M3 residue Q264, further support lipid binding via headgroup ester linkages. Aromatic residues were identified at analogous sites in a handful of eukaryotic family members, including the human GABAA receptor ε subunit, suggesting conservation of relevant interactions in other evolutionary branches. Further mutagenesis experiments indicated that mutations at this site in ε-containing GABAA receptors can change the apparent affinity of the agonist response to GABA, suggesting a potential role of this site in channel gating. In conclusion, this work details type-specific lipid interactions, which adds to our growing understanding of how lipids modulate pLGICs.

Structure of the DP1-DP2 PolD complex bound with DNA and its implications for the evolutionary history of DNA and RNA polymerases.

PolD is an archaeal replicative DNA polymerase (DNAP) made of a proofreading exonuclease subunit (DP1) and a larger polymerase catalytic subunit (DP2). Recently, we reported the individual crystal structures of the DP1 and DP2 catalytic cores, thereby revealing that PolD is an atypical DNAP that has all functional properties of a replicative DNAP but with the catalytic core of an RNA polymerase (RNAP). We now report the DNA-bound cryo-electron microscopy (cryo-EM) structure of the heterodimeric DP1-DP2 PolD complex from Pyrococcus abyssi, revealing a unique DNA-binding site. Comparison of PolD and RNAPs extends their structural similarities and brings to light the minimal catalytic core shared by all cellular transcriptases. Finally, elucidating the structure of the PolD DP1-DP2 interface, which is conserved in all eukaryotic replicative DNAPs, clarifies their evolutionary relationships with PolD and sheds light on the domain acquisition and exchange mechanism that occurred during the evolution of the eukaryotic replisome.

Structure of the human ClC-1 chloride channel.

ClC-1 protein channels facilitate rapid passage of chloride ions across cellular membranes, thereby orchestrating skeletal muscle excitability. Malfunction of ClC-1 is associated with myotonia congenita, a disease impairing muscle relaxation. Here, we present the cryo-electron microscopy (cryo-EM) structure of human ClC-1, uncovering an architecture reminiscent of that of bovine ClC-K and CLC transporters. The chloride conducting pathway exhibits distinct features, including a central glutamate residue ("fast gate") known to confer voltage-dependence (a mechanistic feature not present in ClC-K), linked to a somewhat rearranged central tyrosine and a narrower aperture of the pore toward the extracellular vestibule. These characteristics agree with the lower chloride flux of ClC-1 compared with ClC-K and enable us to propose a model for chloride passage in voltage-dependent CLC channels. Comparison of structures derived from protein studied in different experimental conditions supports the notion that pH and adenine nucleotides regulate ClC-1 through interactions between the so-called cystathionine-ß-synthase (CBS) domains and the intracellular vestibule ("slow gating"). The structure also provides a framework for analysis of mutations causing myotonia congenita and reveals a striking correlation between mutated residues and the phenotypic effect on voltage gating, opening avenues for rational design of therapies against ClC-1-related diseases.

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.