Author: rhoward

Published online 31 July 2023 in PLOS Computational Biology (doi 10.1371/journal.pcbi.1011255):

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

Christian Blau, Linnea Yvonnesdotter, Erik Lindahl

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.

Read the full publication here.

Published online 18 July 2023 in Journal of Chemical Information and Modeling (doi 10.1021/acs.jcim.3c00625):

Understanding drug skin permeation enhancers using molecular dynamics simulations

Christian Wennberg, Magnus Lundborg, Erik Lindahl, Lars Norlén

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.

Read the full publication here.

Published 11 July 2023 in Biophysical Journal (doi 10.1016/j.bpj.2023.05.033):

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

Linnea Yvonnesdotter, Urška Rovšnik, Christian Blau, Marie Lycksell, Rebecca J Howard, Erik Lindahl

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.

Read the full publication here.

Published 5 July 2023 in eLife (doi 10.7554/eLife.84808):

Determinants of sugar-induced influx in the mammalian fructose transporter GLUT5

Sarah E McComas, Tom Reichenbach, Darko Mitrovic, Claudia Alleva, Marta Bonaccorsi, Lucie Delemotte, David Drew

In mammals, glucose transporters (GLUT) control organism-wide blood-glucose homeostasis. In human, this is accomplished by 14 different GLUT isoforms, that transport glucose and other monosaccharides with varying substrate preferences and kinetics. Nevertheless, there is little difference between the sugar-coordinating residues in the GLUT proteins and even the malarial Plasmodium falciparum transporter PfHT1, which is uniquely able to transport a wide range of different sugars. PfHT1 was captured in an intermediate ‘occluded’ state, revealing how the extracellular gating helix TM7b has moved to break and occlude the sugar-binding site. Sequence difference and kinetics indicated that the TM7b gating helix dynamics and interactions likely evolved to enable substrate promiscuity in PfHT1, rather than the sugar-binding site itself. It was unclear, however, if the TM7b structural transitions observed in PfHT1 would be similar in the other GLUT proteins. Here, using enhanced sampling molecular dynamics simulations, we show that the fructose transporter GLUT5 spontaneously transitions through an occluded state that closely resembles PfHT1. The coordination of D-fructose lowers the energetic barriers between the outward- and inward-facing states, and the observed binding mode for D-fructose is consistent with biochemical analysis. Rather than a substrate-binding site that achieves strict specificity by having a high affinity for the substrate, we conclude GLUT proteins have allosterically coupled sugar binding with an extracellular gate that forms the high-affinity transition-state instead. This substrate-coupling pathway presumably enables the catalysis of fast sugar flux at physiological relevant blood-glucose concentrations.

Read the full publication here.

Published 5 July 2023 in eLife (doi 10.7554/eLife.84805):

Reconstructing the transport cycle in the sugar porter superfamily using coevolution-powered machine learning

Darko Mitrovic, Sarah E McComas, Claudia Alleva, Marta Bonaccorsi, David Drew, Lucie Delemotte

Sugar porters (SPs) represent the largest group of secondary-active transporters. Some members, such as the glucose transporters (GLUTs), are well known for their role in maintaining blood glucose homeostasis in mammals, with their expression upregulated in many types of cancers. Because only a few sugar porter structures have been determined, mechanistic models have been constructed by piecing together structural states of distantly related proteins. Current GLUT transport models are predominantly descriptive and oversimplified. Here, we have combined coevolution analysis and comparative modeling, to predict structures of the entire sugar porter superfamily in each state of the transport cycle. We have analyzed the state-specific contacts inferred from coevolving residue pairs and shown how this information can be used to rapidly generate free-energy landscapes consistent with experimental estimates, as illustrated here for the mammalian fructose transporter GLUT5. By comparing many different sugar porter models and scrutinizing their sequence, we have been able to define the molecular determinants of the transport cycle, which are conserved throughout the sugar porter superfamily. We have also been able to highlight differences leading to the emergence of proton-coupling, validating, and extending the previously proposed latch mechanism. Our computational approach is transferable to any transporter, and to other protein families in general.

Read the full publication here.

Published 28 June 2023 in the British Journal of Pharmacology (doi 10.1111/bph.16183):

Subtype-specific modulation of human Kv7 channels by the anticonvulsant cannabidiol through a lipid-exposed pore-domain site

Michael Pökl, Akshay Sridhar, Damon JA Frampton, Veronika A Linhart, Lucie Delemotte, Sara I Liin

Background and Purpose
Cannabidiol (CBD) is used clinically as an anticonvulsant. However, its precise mechanism of action has remained elusive. CBD was recently demonstrated to enhance the activity of the neuronal KV7.2/7.3 channel, which may be one important contributor to CBD’s anticonvulsant effect. Curiously, CBD inhibits the closely related cardiac KV7.1/KCNE1 channel. However, whether and how CBD affects other KV7 subtypes remains unstudied, and the CBD interaction sites mediating these diverse effects remain unknown.

Experimental approach
Here, we used electrophysiology, molecular dynamics simulations, molecular docking, and site-directed mutagenesis to address these open questions.

Key results
We found that CBD modulates the activity of all human KV7 subtypes and that the effects are subtype-dependent. CBD enhanced the activity of KV7.2-7.5 subtypes, seen as a V50 shift towards more negative voltages or increased maximum conductance. In contrast, CBD inhibited the KV7.1 and KV7.1/KCNE1 channels, seen as a V50 shift towards more positive voltages and reduced conductance. In KV7.2 and KV7.4, we propose a CBD interaction site at the subunit interface in the pore domain which overlaps with the interaction site of other compounds, notably the anticonvulsant retigabine. However, CBD relies on other residues for its effects than the conserved tryptophan that is critical for retigabine effects. We propose a similar, though not identical CBD site in KV7.1, with a non-conserved phenylalanine being important.

Conclusions and implications
We identify novel targets of CBD, contributing to a better understanding of CBD’s clinical effects, and provide mechanistic insights into how CBD modulates different KV7 subtypes.

Read the full publication here.

Published 21 June 2023 in Journal of Physical Chemistry Letters (doi 10.1021/acs.jpclett.3c00803):

Modulation of pore opening of eukaryotic sodium channels by π-helices in S6

Koushik Choudhury, Lucie Delemotte

Voltage-gated sodium channels are heterotetrameric sodium selective ion channels that play a central role in electrical signaling in excitable cells. With recent advances in structural biology, structures of eukaryotic sodium channels have been captured in several distinct conformations corresponding to different functional states. The secondary structure of the pore lining S6 helices of subunits DI, DII, and DIV has been captured with both short π-helix stretches and in fully α-helical conformations. The relevance of these secondary structure elements for pore gating is not yet understood. Here, we propose that a π-helix in at least DI-S6, DIII-S6, and DIV-S6 results in a fully conductive state. On the other hand, the absence of π-helix in either DI-S6 or DIV-S6 yields a subconductance state, and its absence from both DI-S6 and DIV-S6 yields a nonconducting state. This work highlights the impact of the presence of a π-helix in the different S6 helices of an expanded pore on pore conductance, thus opening new doors toward reconstructing the entire conformational landscape along the functional cycle of Nav Channels and paving the way to the design of state-dependent modulators.

Read the full publication here.

Published 21 June 2023 in eLife (doi 10.7554/eLife.80303):

Interplay between VSD, pore, and membrane lipids in electromechanical coupling in HCN channels

Ahmad Elbahnsi, John Cowgill, Verena Burtscher, Linda Wedemann, Luise Zeckey, Baron Chanda, Lucie Delemotte

Hyperpolarized-activated and cyclic nucleotide-gated (HCN) channels are the only members of the voltage-gated ion channel superfamily in mammals that open upon hyperpolarization, conferring them pacemaker properties that are instrumental for rhythmic firing of cardiac and neuronal cells. Activation of their voltage-sensor domains (VSD) upon hyperpolarization occurs through a downward movement of the S4 helix bearing the gating charges, which triggers a break in the alpha-helical hydrogen bonding pattern at the level of a conserved Serine residue. Previous structural and molecular simulation studies had however failed to capture pore opening that should be triggered by VSD activation, presumably because of a low VSD/pore electromechanical coupling efficiency and the limited timescales accessible to such techniques. Here, we have used advanced modeling strategies, including enhanced sampling molecular dynamics simulations exploiting comparisons between non-domain swapped voltage-gated ion channel structures trapped in closed and open states to trigger pore gating and characterize electromechanical coupling in HCN1. We propose that the coupling mechanism involves the reorganization of the interfaces between the VSD helices, in particular S4, and the pore-forming helices S5 and S6, subtly shifting the balance between hydrophobic and hydrophilic interactions in a ‘domino effect’ during activation and gating in this region. Remarkably, our simulations reveal state-dependent occupancy of lipid molecules at this emergent coupling interface, suggesting a key role of lipids in hyperpolarization-dependent gating. Our model provides a rationale for previous observations and a possible mechanism for regulation of HCN channels by the lipidic components of the membrane.

Read the full publication here.

Published 23 February 2023 in ACS Chemical Neuroscience (doi 10.1021/acschemneuro.2c00783):

Structural elucidation of ivermectin binding to α7nAChR and the induced channel desensitization

Vasyl Bondarenko, Qiang Chen, Kevin Singewald, Nandan Haloi, Tommy S. Tillman, Rebecca J Howard, Erik Lindahl, Yan Xu, Pei Tang

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.

Read the full publication here.

Published 14 December 2022 in Journal of General Physiology (doi 10.1085/jgp.202213214):

An α–π transition in S6 shapes the conformational cycle of the bacterial sodium channel NavAb

Koushik Choudhury, Rebecca J Howard, Lucie Delemotte

Voltage-gated sodium channels play an important role in electrical signaling in excitable cells. In response to changes in membrane potential, they cycle between nonconducting and conducting conformations. With recent advances in structural biology, structures of sodium channels have been captured in several distinct conformations, which are thought to represent different functional states. However, it has been difficult to capture the intrinsically transient open state. We recently showed that a proposed open state of the bacterial sodium channel NavMs was not conductive and that a conformational change involving a transition to a π-helix in the pore-lining S6 helix converted this structure into a conducting state. However, the relevance of this structural feature in other sodium channels, and its implications for the broader gating cycle, remained unclear. Here, we propose a comparable open state of another class of bacterial channel from Aliarcobacter butzleri (NavAb) with characteristic pore hydration, ion permeation, and drug binding properties. Furthermore, we show that a π-helix transition can lead to pore opening and that such a conformational change blocks fenestrations in the inner helix bundle. We also discover that a region in the C-terminal domain can undergo a disordering transition proposed to be important for pore opening. These results support a role for a π-helix transition in the opening of NavAb, enabling new proposals for the structural annotation and drug modulation mechanisms in this important sodium channel model.

Read the full publication here.

Published 8 December 2022 in Proceedings of the National Academy of Sciences of the USA (doi 10.1073/pnas.2210669119):

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

Marie Lycksell, Urška Rovšnik, Anton Hanke, Anne Martel, Rebecca J Howard, Erik Lindahl

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.

Read the full publication here.

Published 18 October 2022 in Biophysical Journal (doi 10.1016/j.bpj.2022.09.009):

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

Magnus Lundborg, Christian Wennberg, Jack Lidmar, Berk Hess, Erik Lindahl, Lars Norlén

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 K_Pexp-calc mean square difference was 0.9 cm² h⁻².

Read the full publication here.

Published 17 October 2022 in Proceedings of the National Academy of Sciences of the USA (doi 10.1073/pnas.2208081119):

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

Yuxuan Zhuang, Colleen M Noviello, Ryan E Hibbs, Rebecca J Howard, Erik Lindahl

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that modulates neuronal excitability, largely by allowing Ca²⁺ 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 Ca²⁺. 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.

Read the full publication here.

Published 7 October 2022 in PLoS Computational Biology (doi 10.1371/journal.pcbi.1010583):

Markov state modelling reveals heterogeneous drug-inhibition mechanism of Calmodulin

Annie M Westerlund, Akshay Sridhar, Leo Dahl, Alma Andersson, Anna-Yaroslava Bodnar, Lucie Delemotte

Calmodulin (CaM) is a calcium sensor which binds and regulates a wide range of target-proteins. This implicitly enables the concentration of calcium to influence many downstream physiological responses, including muscle contraction, learning and depression. The antipsychotic drug trifluoperazine (TFP) is a known CaM inhibitor. By binding to various sites, TFP prevents CaM from associating to target-proteins. However, the molecular and state-dependent mechanisms behind CaM inhibition by drugs such as TFP are largely unknown. Here, we build a Markov state model (MSM) from adaptively sampled molecular dynamics simulations and reveal the structural and dynamical features behind the inhibitory mechanism of TFP-binding to the C-terminal domain of CaM. We specifically identify three major TFP binding-modes from the MSM macrostates, and distinguish their effect on CaM conformation by using a systematic analysis protocol based on biophysical descriptors and tools from machine learning. The results show that depending on the binding orientation, TFP effectively stabilizes features of the calcium-unbound CaM, either affecting the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the bound domain. The conclusions drawn from this work may in the future serve to formulate a complete model of pharmacological modulation of CaM, which furthers our understanding of how these drugs affect signaling pathways as well as associated diseases.

Read the full publication here.

Available online 27 September 2022 in Physics of Fluids (doi 10.1063/5.0121144):

Asymmetry of wetting and de-wetting on high-friction surfaces originates from the same molecular physics

Michele Pellegrino, Berk Hess

The motion of three-phase contact lines is one of the most relevant research topics of micro- and nano-fluidics. According to many hydrodynamic and molecular models, the dynamics of contact lines is assumed overdamped and dominated by localised liquid-solid friction, entailing the existence of a mobility relation between contact line speed and microscopic contact angle. We present and discuss a set of non-equilibrium atomistic Molecular Dynamics simulations of water nanodroplets spreading on or confined between silica-like walls, showing the existence of the aforementioned relation and its invariance under wetting modes (‘spontaneous’ or ‘forced’). Upon changing the wettability of the walls, it has been noticed that more hydrophilic substrates are easier to wet rather than de-wet; we show how this asymmetry can be automatically captured by a contact line friction model that accounts for the molecular transport between liquid layers. A simple examination of the order and orientation of near-contact-line water molecules corroborates the physical foundation of the model. Lastly, we propose an approach to discriminate between contact line friction models which overcomes the limitations of experimental resolution. This work constitutes a stepping stone towards demystifying wetting dynamics on high-friction hydrophilic substrates and underlines the relevance of contact line friction in modelling the motion of three-phase contact lines.

Read the full publication here.

Published 21 September 2022 in Journal of Chemical Theory and Computation (doi 10.1021/acs.jctc.2c00516):

Scalable Constant pH Molecular Dynamics in GROMACS

Noora Aho*, Pavel Buslaev*, Anton Jansen, Paul Bauer, Gerrit Groenhof*, Berk Hess*

Molecular dynamics (MD) computer simulations are used routinely to compute atomistic trajectories of complex systems. Systems are simulated in various ensembles, depending on the experimental conditions one aims to mimic. While constant energy, temperature, volume, and pressure are rather straightforward to model, pH, which is an equally important parameter in experiments, is more difficult to account for in simulations. Although a constant pH algorithm based on the λ-dynamics approach by Brooks and co-workers was implemented in a fork of the GROMACS molecular dynamics program, uptake has been rather limited, presumably due to the poor scaling of that code with respect to the number of titratable sites. To overcome this limitation, we implemented an alternative scheme for interpolating the Hamiltonians of the protonation states that makes the constant pH molecular dynamics simulations almost as fast as a normal MD simulation with GROMACS. In addition, we implemented a simpler scheme, called multisite representation, for modeling side chains with multiple titratable sites, such as imidazole rings. This scheme, which is based on constraining the sum of the λ-coordinates, not only reduces the complexity associated with parametrizing the intramolecular interactions between the sites but also is easily extendable to other molecules with multiple titratable sites. With the combination of a more efficient interpolation scheme and multisite representation of titratable groups, we anticipate a rapid uptake of constant pH molecular dynamics simulations within the GROMACS user community.

Read the full publication here, and accompanying manuscript here.

*corresponding authors