Control of ion selectivity in LeuT: two Na+ binding sites with two different mechanisms

Crown-Ether Crystal Channel Membranes with Subnanometer Pores for Selective Na+ Transport

Emulating biological sodium ion channels to achieve high selectivity and rapid Na+ transport is important for water desalination, energy conversion, and separation processes. However, the development of artificial ion channels, especially multichannels, to achieve high ion selectivity, remains a challenge. In this work, we demonstrate the fabrication of ion channel membranes utilizing crown-ether crystals (DA18C6-nitrate crystals), which feature extremely consistent subnanometer pores. The polyethylene terephthalate (PET) membranes were initially subjected to amination, followed by the in situ growth of DA18C6-nitrate crystals to establish ordered multichannels aimed at facilitating selective Na+ conductance. These channels allow rapid Na+ transport while inhibiting the migration of other ions (K+ and Ca2+). The Na+ transport rate was 2.15 mol m-2 h-1, resulting in the Na+/K+ and Na+/Ca2+ selectivity ratios of 6.53 and 12.56, respectively. Due to the immobilization of the crown-ether ring, when the size of the transmembrane ion exceeded that of the crown-ether ring's cavity, the ions had to undergo a dehydration process to pass through the channel. This resulted in the ions encountering a higher energy barrier upon entering the channel, making it more difficult for them to permeate. However, the size of Na+ was compatible with the cavity of the crown-ether ring and was able to displace the hydrated layer effectively, facilitating selective Na+ translocation. In summary, this research offers a promising approach for the future development of functionalized ion channels and efficient membrane materials tailored for high-performance Na+ separation.

Flow-based bioconjugation of coumarin phosphatidylethanolamine probes: Optimised synthesis and membrane molecular dynamics studies

A series of phosphatidylethanolamine fluorescent probes head-labelled with 3-carboxycoumarin was prepared by an improved bioconjugation approach through continuous flow synthesis. The established procedure, supported by a design of experiment (DoE) set-up, resulted in a significant reduction in the reaction time compared to the conventional batch method, in addition to a minor yield increase. The characterization of these probes was enhanced by an in-depth molecular dynamics (MD) study of the behaviour of a representative probe of this family, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine labelled with 3-carboxycoumarin (POPE-COUM), in bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/1-stearoyl-2-linoleoyl-sn-glycero-3-phosphocholine (SLPC) 2:1, mimicking the composition of the egg yolk lecithin membranes recently used experimentally by our group to study POPE-COUM as a biomarker of the oxidation state and integrity of large unilamellar vesicles (LUVs). The MD simulations revealed that the coumarin group is oriented towards the bilayer interior, leading to a relatively internal location, in agreement with what is observed in the nitrobenzoxadiazole fluorophore of commercial head-labelled NBD-PE probes. This behaviour is consistent with the previously stated hypothesis that POPE-COUM is entirely located within the LUVs structure. Hence, the delay on the oxidation of the probe in the oxygen radical absorbance capacity (ORAC) assays performed is related with the inaccessibility of the probe until alteration of the LUV structure occurs. Furthermore, our simulations show that POPE-COUM exerts very little global and local perturbation on the host bilayer, as evaluated by key properties of the unlabelled lipids. Together, our findings establish PE-COUM as suitable fluorescent lipid analogue probes.

Mechanistic understanding of KCNQ1 activating polyunsaturated fatty acid analogs

The KCNQ1 channel is important for the repolarization phase of the cardiac action potential. Loss of function mutations in KCNQ1 can cause long QT syndrome (LQTS), which can lead to cardiac arrythmia and even sudden cardiac death. We have previously shown that polyunsaturated fatty acids (PUFAs) and PUFA analogs can activate the cardiac KCNQ1 channel, making them potential therapeutics for the treatment of LQTS. PUFAs bind to KCNQ1 at two different binding sites: one at the voltage sensor (Site I) and one at the pore (Site II). PUFA interaction at Site I shifts the voltage dependence of the channel to the left, while interaction at Site II increases maximal conductance. The PUFA analogs, linoleic-glycine and linoleic-tyrosine, are more effective than linoleic acid at Site I, but less effective at Site II. Using both simulations and experiments, we find that the larger head groups of linoleic-glycine and linoleic-tyrosine interact with more residues than the smaller linoleic acid at Site I. We propose that this will stabilize the negatively charged PUFA head group in a position to better interact electrostatically with the positively charges in the voltage sensor. In contrast, the larger head groups of linoleic-glycine and linoleic-tyrosine compared with linoleic acid prevent a close fit of these PUFA analogs in Site II, which is more confined. In addition, we identify several KCNQ1 residues as critical PUFA-analog binding residues, thereby providing molecular models of specific interactions between PUFA analogs and KCNQ1. These interactions will aid in future drug development based on PUFA-KCNQ1 channel interactions.

Interaction of Near-Infrared (NIR)-Light Responsive Probes with Lipid Membranes: A Combined Simulation and Experimental Study

Cancer is considered a major societal challenge for the next decade worldwide. Developing strategies for simultaneous diagnosis and treatment has been considered a promising tool for fighting cancer. For this, the development of nanomaterials incorporating prototypic near-infrared (NIR)-light responsive probes, such as heptamethine cyanines, has been showing very promising results. The heptamethine cyanine-incorporating nanomaterials can be used for a tumor's visualization and, upon interaction with NIR light, can also produce a photothermal/photodynamic effect with a high spatio-temporal resolution and minimal side effects, leading to an improved therapeutic outcome. In this work, we studied the interaction of 12 NIR-light responsive probes with lipid membrane models by molecular dynamics simulations. We performed a detailed characterization of the location, orientation, and local perturbation effects of these molecules on the lipid bilayer. Based on this information, the probes were divided into two groups, predicting a lower and higher perturbation of the lipid bilayer. From each group, one molecule was selected for testing in a membrane leakage assay. The experimental data validate the hypothesis that molecules with charged substituents, which function as two polar anchors for the aqueous phase while spanning the membrane thickness, are more likely to disturb the membrane by the formation of defects and pores, increasing the membrane leakage. The obtained results are expected to contribute to the selection of the most suitable molecules for the desired application or eventually guiding the design of probe modifications for achieving an optimal interaction with tumor cell membranes.

Modeling Zeta Potential for Nanoparticles in Solution: Water Flexibility Matters

Nonequilibrium molecular dynamics simulations were performed to study the electrokinetic properties of five mainstream TIPxP water models (namely, TIP3P-FB, TIP3Pm, TIP4P-FB, TIP4P-Ew, and TIP4P/2005) in NaCl aqueous solutions in the presence of a negatively charged TiO₂ surface. The impact of solvent flexibility and system geometry on the electro-osmotic (EO) mobility and flow direction was systematically assessed and compared. We found that lack of water flexibility decelerates the forward EO flow of aqueous solutions at moderate (0.15 M) or high (0.30 M) NaCl concentrations, in some special cases to such an extent that EO flow reversal occurs. Zeta potential (ZP) values were then determined from the bulk EO mobilities using the Helmholtz-Smoluchowski formula. The straight comparison against available experimental data strongly suggests that water flexibility improves the ZP determination of NaCl solutions adjacent to a realistic TiO₂ surface under neutral pH conditions.

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

The pH Response of a Peptoid Oligomer

Polypeptoids are N-substituted glycine polymers, which differ from peptides in the placement of the side chain on the amide nitrogen rather than the C_α carbon. A peptoid with a chiral side chain containing both an aromatic group and carboxylic acid has a structure that responds to pH changes. All-atom molecular dynamics simulations using a force field specifically tuned for peptoids were carried out with an advanced sampling method for the peptoid (S)-N-(1-carboxy-2-phenylethyl)glycine in the high and low pH limits. The simulations show that the structure changes from mostly cis amide bonds at low pH to mostly trans bonds at high pH. The structural changes are driven by side chain-backbone hydrogen bonds at low pH and side chain repulsions and increased water contact at high pH.

Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study

Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.

The impact of pathogenic and artificial mutations on Claudin-5 selectivity from molecular dynamics simulations

Tight-junctions (TJs) are multi-protein complexes between adjacent endothelial or epithelial cells. In the blood-brain-barrier (BBB), they seal the paracellular space and the Claudin-5 (Cldn5) protein forms their backbone. Despite the fundamental role in brain homeostasis, little is known on Cldn5-based TJ assemblies. Different structural models were suggested, with Cldn5 protomers generating paracellular pores that restrict the passage of ions and small molecules. Recently, the first Cldn5 pathogenic mutation, G60R, was identified and shown to induce Cl^--selective channels and Na⁺ barriers in BBB TJs, providing an excellent opportunity to validate the structural models. Here, we used molecular dynamics to study the permeation of ions and water through two distinct G60R-Cldn5 paracellular architectures. Only the so-called Pore I reproduces the functional modification observed in experiments, displaying a free energy (FE) minimum for Cl^- and a barrier for Na⁺ consistent with anionic selectivity. We also studied the artificial Q57D and Q63D mutations in the constriction region, Q57 being conserved in Cldns except for cation permeable homologs. In both cases, we obtain FE profiles consistent with facilitated passage of cations. Our calculations provide the first in-silico description of a Cldn5 pathogenic mutation, further assessing the TJ Pore I model and yielding new insight on BBB's paracellular selectivity.

How Accurate Is the Egg-Box Model in Describing the Binding of Calcium to Polygalacturonate? A Molecular Dynamics Simulation Study

We performed molecular dynamics (MD) simulations of octameric galacturonate, GalA₈, chains in the presence of Ca²⁺ in a ratio of R = [Ca²⁺]/[GalA] = 0.25 in order to determine to which extent the popular "egg-box model" (EBM) is able to describe the association between Ca²⁺ cations and polygalacturonate (polyGalA) chains. To this aim, we slightly revised the empirical parameters for the interaction between Ca²⁺ and the carboxylate oxygen atoms of GalA units so as to reproduce the experimental Ca²⁺-GalA association constant. We also defined an ad hoc order parameter, referred to as the egg-box score (EBS), that quantifies any deviation of the local coordination geometry of calcium cations with respect to an "ideal" EBM coordination geometry. The results reveal that the local coordination geometry of Ca²⁺ cations bound to polyGalA chains differs from that of the EBM. Moreover, polyGalA chains exhibit significant conformational disorder, and the cross-link angles formed between polyGalA chain axes are broadly distributed. Overall, the present study suggests that the EBM fails to describe accurately the association modes between calcium and polyGalA chains at a molar ratio R of 0.25.

The voltage-sensing domain of a hERG1 mutant is a cation-selective channel

A cationic leak current known as an "omega current" may arise from mutations of the first charged residue in the S4 of the voltage sensor domains of sodium and potassium voltage-gated channels. The voltage-sensing domains (VSDs) in these mutated channels act as pores allowing nonspecific passage of cations, such as Li+, K+, Cs+, and guanidinium. Interestingly, no omega currents have been previously detected in the nonswapped voltage-gated potassium channels such as the human-ether-a-go-go-related (hERG1), hyperpolarization-activated cyclic nucleotide-gated, and ether-a-go-go channels. In this work, we discovered a novel omega current by mutating the first charged residue of the S4 of the hERG1, K525 to serine. To characterize this omega current, we used various probes, including the hERG1 pore domain blocker, dofetilide, to show that the omega current does not require cation flux via the canonical pore domain. In addition, the omega flux does not cross the conventional selectivity filter. We also show that the mutated channel (K525S hERG1) conducts guanidinium. These data are indicative of the formation of an omega current channel within the VSD. Using molecular dynamics simulations with replica-exchange umbrella sampling simulations of the wild-type hERG1 and the K525S hERG1, we explored the molecular underpinnings governing the cation flow in the VSD of the mutant. We also show that the wild-type hERG1 may form water crevices supported by the biophysical surface accessibility data. Overall, our multidisciplinary study demonstrates that the VSD of hERG1 may act as a cation-selective channel wherein a mutation of the first charged residue in the S4 generates an omega current. Our simulation uncovers the atomistic underpinning of this mechanism.

Atomistic Modeling of Jet Formation in Charged Droplets

The first atomistic simulations that reveal the mechanism of Rayleigh fission are presented. It is demonstrated that simple ion or macroion ejection takes place through droplet deformation from a spherical into a distinct "tear" shape that contains a conical protrusion. We assert that the latter state is a free-energy minimum along an order parameter that measures the degree of droplet asphericity. The charged droplet's long-time evolution proceeds by alternating between the two minima above and below the critical value that are reached through solvent evaporation and ion ejection, respectively. For the first time, this mechanism allows one to explain the nature of the progeny droplets and the percentage of charge lost during fission. The cone half angle is estimated and found to be in good agreement with the value predicted from the solution of the electrostatic equation for the dielectric liquid. It is found that the conical deformation is independent of the effect of electrohydrodynamic forces reported in experiments. Contrary to the experimental observations of two diametrically opposite jets for droplets suspended in the electric field, we find that a single jet is formed at the Rayleigh limit. The study provides insight into the mechanism of capture of a macroion in jets appearing in electrospray ionization mass spectrometry (ESI-MS) experiments and may explain the tolerance of the ESI-MS spectrum to salt contamination of the sample.

Computational study of ion permeation through claudin-4 paracellular channels

Claudins (Cldns) form a large family of protein homologs that are essential for the assembly of paracellular tight junctions (TJs), where they form channels or barriers with tissue-specific selectivity for permeants. In contrast to several family members whose physiological role has been identified, the function of claudin 4 (Cldn4) remains elusive, despite experimental evidence suggesting that it can form anion-selective TJ channels in the renal epithelium. Computational approaches have recently been employed to elucidate the molecular basis of Cldns' function, and hence could help in clarifying the role of Cldn4. In this work, we use structural modeling and all-atom molecular dynamics simulations to transfer two previously introduced structural models of Cldn-based paracellular complexes to Cldn4 to reproduce a paracellular anion channel. Free energy calculations for ionic transport through the pores allow us to establish the thermodynamic properties driving the ion-selectivity of the structures. While one model shows a cavity permeable to chloride and repulsive to cations, the other forms barrier to the passage of all the major physiological ions. Furthermore, our results confirm the charge selectivity role of the residue Lys65 in the first extracellular loop of the protein, rationalizing Cldn4 control of paracellular permeability.

Density Functional Theory Study on the Selective Reductive Amination of Aldehydes and Ketones over Their Reductions to Alcohols Using Sodium Triacetoxyborohydride

Reductive amination is one of the most important methods to synthesize amines, having a wide application in the pharmaceutical, fine chemicals, and materials industries. In general, the reaction begins with dehydration between a carbonyl compound and an amine compound, forming an imine, which is then reduced to an alkylated amine product. Sodium triacetoxyborohydride (STAB) is a popular choice for the reducing agent as it shows selectivity for imines over aldehydes and ketones, which is particularly important in direct reductive amination where the imine and carbonyl compounds are present concurrently. Here, we analyze the reaction pathways of acid-catalyzed direct reductive amination in 1,2-dichloroethane (DCE) with acetaldehyde and methylamine. We find that the transition states for the formation and subsequent reduction of Z-methylethylideneimine (resultant aldimine from acetaldehyde and methylamine) have lower energies than the reduction of acetaldehyde. Transition state structures for the hydride transfers are organized by the Lewis-acidic sodium ion. Additionally, reduction reactions with formaldehyde and acetone and their imine derivatives (with methylamine) are investigated, and again, the hydride transfer to the resultant aldimine or ketimine is lower in energy than that of their parent carbonyl compound.

The Single Residue K12 Governs the Exceptional Voltage Sensitivity of Mitochondrial Voltage-Dependent Anion Channel Gating

The voltage-dependent anion channel (VDAC) is a β-barrel channel of the mitochondrial outer membrane (MOM) that passively transports ions, metabolites, polypeptides, and single-stranded DNA. VDAC responds to a transmembrane potential by "gating," i.e. transitioning to one of a variety of low-conducting states of unknown structure. The gated state results in nearly complete suppression of multivalent mitochondrial metabolite (such as ATP and ADP) transport, while enhancing calcium transport. Voltage gating is a universal property of β-barrel channels, but VDAC gating is anomalously sensitive to transmembrane potential. Here, we show that a single residue in the pore interior, K12, is responsible for most of VDAC's voltage sensitivity. Using the analysis of over 40 μs of atomistic molecular dynamics (MD) simulations, we explore correlations between motions of charged residues inside the VDAC pore and geometric deformations of the β-barrel. Residue K12 is bistable; its motions between two widely separated positions along the pore axis enhance the fluctuations of the β-barrel and augment the likelihood of gating. Single channel electrophysiology of various K12 mutants reveals a dramatic reduction of the voltage-induced gating transitions. The crystal structure of the K12E mutant at a resolution of 2.6 Å indicates a similar architecture of the K12E mutant to the wild type; however, 60 μs of atomistic MD simulations using the K12E mutant show restricted motion of residue 12, due to enhanced connectivity with neighboring residues, and diminished amplitude of barrel motions. We conclude that β-barrel fluctuations, governed particularly by residue K12, drive VDAC gating transitions.

Impacts of targeting different hydration free energy references on the development of ion potentials

Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol^-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H⁺). The choice of ΔG°(H⁺) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K⁺, Mg²⁺, Al³⁺, and Cl^- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H⁺) for the cation and an increase in ΔG°(H⁺) for the anion. The increase in ΔG°(H⁺) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H⁺) value of -259.8 kcal mol^-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H⁺) for Mg²⁺ model is better for a reasonable description of the metal binding network.

Binding of polar and hydrophobic molecules at the LiCoO2 (001)-water interface: force field development and molecular dynamics simulations

A classical model in the framework of the INTERFACE force field has been developed for treating the LiCoO₂ (LCO) (001)/water interface. In comparison to ab initio molecular dynamics (MD) simulations based on density functional theory, MD simulations using the classical model lead to generally reliable descriptions of interfacial properties, such as the density distribution of water molecules. Water molecules in close contact with the LCO surface form a strongly adsorbed layer, which leads to a free energy barrier for the adsorption of polar or charged molecules to the LCO surface. Moreover, due to the strong hydrogen bonding interactions with the LCO surface, the first water layer forms an interface that exhibits hydrophobic characters, leading to favorable adsorption of non-polar molecules to the interface. Therefore, despite its highly polar nature, the LCO (001) surface binds not only polar/charged but also non-polar solutes. As an application, the model is used to analyze the adsorption of reduced nicotinamide adenine dinucleotide (NADH) and its molecular components to the LCO (001) surface in water. The results suggest that recently observed redox activity of NADH at the LCO/water interface was due to the co-operativity between the ribose component, which drives binding to the LCO surface, and the nicotinamide moiety, which undergoes oxidation.

A molecular switch controls the impact of cholesterol on a Kir channel

Cholesterol is one of the main components found in plasma membranes and is involved in lipid-dependent signaling enabled by integral membrane proteins such as inwardly rectifying potassium (Kir) channels. Similar to other ion channels, most of the Kir channels are down-regulated by cholesterol. One of the very few notable exceptions is Kir3.4, which is up-regulated by this important lipid. Here, we discovered and characterized a molecular switch that controls the impact (up-regulation vs. down-regulation) of cholesterol on Kir3.4. Our results provide a detailed molecular mechanism of tunable cholesterol regulation of a potassium channel.

Cholesterol decreases the activity of the majority of ion channels while increasing the activity of only a few, yet it remains unclear how. Here, we used the inwardly rectifying potassium channel Kir3.4, which is up-regulated by cholesterol, as a tool to address this question. Employing mutagenesis and electrophysiology, we discovered a molecular switch that controls the impact of cholesterol on the channel. Through a single point mutation at position 182 in the transmembrane domain of Kir3.4, we converted the cholesterol-driven up-regulation of the channel into down-regulation. Microseconds-long coarse-grained and atomistic molecular dynamics simulations revealed that the effect of the point mutation propagated toward the selectivity filter of the channel whose conformation controls the conductance of the channel. Planar lipid bilayer experiments validated these results, showing that although cholesterol up-regulated Kir3.4 by increasing its open probability, cholesterol down-regulated the mutant by decreasing its conductance. Further studies underscored the role of mutation-specific alterations of cholesterol distribution in proximity to the channel in cholesterol’s impact on channel activity, highlighting the role of subtle molecular differences in determining how cholesterol distributes around proteins and affects their function.

Development of CHARMM Additive Potential Energy Parameters for α-Methyl Amino Acids

Potential energy parameters for α-methyl amino acids were generated with ab initio calculations on α-methyl-N-acetylalanyl-N'-methylamide (the α-methyl "alanine dipeptide") which served as an input to a grid-based correction to the backbone torsional potential (known as CMAP) consistent with the CHARMM36m additive protein force field. The new parameters were validated by comparison with experimentally determined helicities of the 22 residue C-terminal peptide (H10) from apolipoprotein A1 and five α-methylated variants in water and 0.3:0.7 trifluoroethanol (TFE)/water. Conventional molecular dynamics simulation totaling 30 μs for each peptide is in overall good agreement with the experiment, including the increased helicity in 30% TFE. An additional 500 ns of simulation using two-dimensional dihedral biasing (bpCMAP) replica exchange reduced left-handed conformations, increased right-handed helices, and thereby mostly decreased agreement with the experiment. Analysis of side chain-side chain salt bridges suggests that the overestimation of the helical content may be, in part, due to such interactions. The increased helicity of the peptides in 30% TFE arises from decreased hydrogen bonding of the backbone atoms to water and a concomitant increase in intramolecular backbone hydrogen bonds.

Open-state structure and pore gating mechanism of the cardiac sodium channel

The heartbeat is initiated by voltage-gated sodium channel Na_V1.5, which opens rapidly and triggers the cardiac action potential; however, the structural basis for pore opening remains unknown. Here, we blocked fast inactivation with a mutation and captured the elusive open-state structure. The fast inactivation gate moves away from its receptor, allowing asymmetric opening of pore-lining S6 segments, which bend and rotate at their intracellular ends to dilate the activation gate to ∼10 Å diameter. Molecular dynamics analyses predict physiological rates of Na⁺ conductance. The open-state pore blocker propafenone binds in a high-affinity pose, and drug-access pathways are revealed through the open activation gate and fenestrations. Comparison with mutagenesis results provides a structural map of arrhythmia mutations that target the activation and fast inactivation gates. These results give atomic-level insights into molecular events that underlie generation of the action potential, open-state drug block, and fast inactivation of cardiac sodium channels, which initiate the heartbeat.

A general mechanism of KCNE1 modulation of KCNQ1 channels involving non-canonical VSD-PD coupling

Voltage-gated KCNQ1 channels contain four separate voltage-sensing domains (VSDs) and a pore domain (PD). KCNQ1 expressed alone opens when the VSDs are in an intermediate state. In cardiomyocytes, KCNQ1 co-expressed with KCNE1 opens mainly when the VSDs are in a fully activated state. KCNE1 also drastically slows the opening of KCNQ1 channels and shifts the voltage dependence of opening by >40 mV. We here show that mutations of conserved residues at the VSD-PD interface alter the VSD-PD coupling so that the mutant KCNQ1/KCNE1 channels open in the intermediate VSD state. Using recent structures of KCNQ1 and KCNE beta subunits in different states, we present a mechanism by which KCNE1 rotates the VSD relative to the PD and affects the VSD-PD coupling of KCNQ1 channels in a non-canonical way, forcing KCNQ1/KCNE1 channels to open in the fully-activated VSD state. This would explain many of the KCNE1-induced effects on KCNQ1 channels.

Determinants of ion selectivity in ASIC1a- and ASIC2a-containing acid-sensing ion channels

Trimeric acid-sensing ion channels (ASICs) contribute to neuronal signaling by converting extracellular acidification into excitatory sodium currents. Previous work with homomeric ASIC1a implicates conserved leucine (L7') and consecutive glycine-alanine-serine (GAS belt) residues near the middle, and conserved negatively charged (E18') residues at the bottom of the pore in ion permeation and/or selectivity. However, a conserved mechanism of ion selectivity throughout the ASIC family has not been established. We therefore explored the molecular determinants of ion selectivity in heteromeric ASIC1a/ASIC2a and homomeric ASIC2a channels using site-directed mutagenesis, electrophysiology, and molecular dynamics free energy simulations. Similar to ASIC1a, E18' residues create an energetic preference for sodium ions at the lower end of the pore in ASIC2a-containing channels. However, and in contrast to ASIC1a homomers, ion permeation through ASIC2a-containing channels is not determined by L7' side chains in the upper part of the channel. This may be, in part, due to ASIC2a-specific negatively charged residues (E59 and E62) that lower the energy of ions in the upper pore, thus making the GAS belt more important for selectivity. This is confirmed by experiments showing that the L7'A mutation has no effect in ASIC2a, in contrast to ASIC1a, where it eliminated selectivity. ASIC2a triple mutants eliminating both L7' and upper charges did not lead to large changes in selectivity, suggesting a different role for L7' in ASIC2a compared with ASIC1a channels. In contrast, we observed measurable changes in ion selectivity in ASIC2a-containing channels with GAS belt mutations. Our results suggest that ion conduction and selectivity in the upper part of the ASIC pore may differ between subtypes, whereas the essential role of E18' in ion selectivity is conserved. Furthermore, we demonstrate that heteromeric channels containing mutations in only one of two ASIC subtypes provide a means of functionally testing mutations that render homomeric channels nonfunctional.

Investigating the Mechanism of Sodium Binding to SERT Using Direct Simulations

The serotonin transporter (SERT) terminates neurotransmission by transporting serotonin from the synapse into the pre-synaptic nerve terminal. Altered SERT function leads to several neurological diseases including depression, anxiety, mood disorders, and attention deficit hyperactivity disorders (ADHD). Accordingly SERT is the target for their pharmacological treatments, but also targeted by multiple drugs of abuse. Transport of serotonin by SERT is energized by the transmembrane electrochemical gradient of sodium. We used extensive molecular dynamics simulations to investigate the process of sodium binding to SERT, which is the first step in the transport cycle that leads to serotonin uptake. Comparing data from 51 independent simulations, we find a remarkably well-defined path for sodium entry and could identify two transient binding sites, while observing binding kinetics that are comparable to experimental data. Importantly, the structure and dynamics of the sodium binding sites indicate that sodium binding is accompanied by an induced-fit mechanism that leads to new conformations and reduces local dynamics.

Identification of PUFA interaction sites on the cardiac potassium channel KCNQ1

Polyunsaturated fatty acids (PUFAs), but not saturated fatty acids, modulate ion channels such as the cardiac KCNQ1 channel, although the mechanism is not completely understood. Using both simulations and experiments, we find that PUFAs interact directly with the KCNQ1 channel via two different binding sites: one at the voltage sensor and one at the pore. These two amphiphilic binding pockets stabilize the negatively charged PUFA head group by electrostatic interactions with R218, R221, and K316, while the hydrophobic PUFA tail is selectively stabilized by cassettes of hydrophobic residues. The rigid saturated tail of stearic acid prevents close contacts with KCNQ1. By contrast, the mobile tail of PUFA linoleic acid can be accommodated in the crevice of the hydrophobic cassette, a defining feature of PUFA selectivity in KCNQ1. In addition, we identify Y268 as a critical PUFA anchor point underlying fatty acid selectivity. Combined, this study provides molecular models of direct interactions between PUFAs and KCNQ1 and identifies selectivity mechanisms. Long term, this understanding may open new avenues for drug development based on PUFA mechanisms.

Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations

Activation of voltage-gated ion channels is regulated by conformational changes of the voltage sensor domains (VSDs), four water- and ion-impermeable modules peripheral to the central, permeable pore domain. Anomalous currents, defined as ω-currents, have been recorded in response to mutations of residues on the VSD S4 helix and associated with ion fluxes through the VSDs. In humans, gene defects in the potassium channel Kv7.2 result in a broad range of epileptic disorders, from benign neonatal seizures to severe epileptic encephalopathies. Experimental evidence suggests that the R207Q mutation in S4, associated with peripheral nerve hyperexcitability, induces ω-currents at depolarized potentials, but the fine structural details are still elusive. In this work, we use atom-detailed molecular dynamics simulations and a refined model structure of the Kv7.2 VSD in the active conformation in a membrane/water environment to study the effect of R207Q and four additional mutations of proven clinical importance. Our results demonstrate that the R207Q mutant shows the most pronounced increase of hydration in the internal VSD cavity, a feature favoring the occurrence of ω-currents. Free energy and kinetics calculations of sodium permeation through the native and mutated VSD indicate as more favorable the formation of a cationic current in the latter. Overall, our simulations establish a mechanistic linkage between genetic variations and their physiological outcome, by providing a computational description that includes both thermodynamic and kinetic features of ion permeation associated with ω-currents.

Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel

Despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ∼50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water-channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.

Refinement of a cryo-EM structure of hERG: Bridging structure and function

The human-ether-a-go-go-related gene (hERG) encodes the voltage-gated potassium channel (KCNH2 or Kv11.1, commonly known as hERG). This channel plays a pivotal role in the stability of phase 3 repolarization of the cardiac action potential. Although a high-resolution cryo-EM structure is available for its depolarized (open) state, the structure surprisingly did not feature many functionally important interactions established by previous biochemical and electrophysiology experiments. Using molecular dynamics flexible fitting (MDFF), we refined the structure and recovered the missing functionally relevant salt bridges in hERG in its depolarized state. We also performed electrophysiology experiments to confirm the functional relevance of a novel salt bridge predicted by our refinement protocol. Our work shows how refinement of a high-resolution cryo-EM structure helps to bridge the existing gap between the structure and function in the voltage-sensing domain (VSD) of hERG.

Molecular Characterization of the Surface Excess Charge Layer in Droplets

The surface excess charge layer (SECL) in droplets has often been associated with distinct chemistry. We examine the effect of the nature of ions in the composition and structure of SECL by using molecular dynamics. We find that in the presence of simple ions the thickness of SECL is invariant not only with respect to droplet size but also with respect to the nature of the ions. In the presence of simple ions, this layer has a thickness of ∼1.5-1.7 nm but in the presence of macroions it may extend to ∼2.0 nm. The proportion of ions contained in SECL depends on the nature of the ions and the droplet size. For the same droplet size, I^- and model H₃O⁺ ions show considerably higher concentration than Na⁺ and Cl^- ions. We identify the maximum ion concentration region, which, in nanodrops, may partially overlap with SECL. As the relative shape fluctuations decrease when microdrop size is approached, the overlap between SECL and maximum ion concentration region increases. We suggest the extension of the bilayer droplet structure assumed in the equilibrium partitioning model of Enke to include the maximum ion concentration region that may not coincide with SECL in nanodrops. We compute the ion concentrations in SECL, which are those that should enter the kinetic equation in the ion-evaporation mechanism, instead of the overall drop ion concentration that has been used.

Structural basis for voltage-sensor trapping of the cardiac sodium channel by a deathstalker scorpion toxin

Voltage-gated sodium (Na_V) channels initiate action potentials in excitable cells, and their function is altered by potent gating-modifier toxins. The α-toxin LqhIII from the deathstalker scorpion inhibits fast inactivation of cardiac Na_V1.5 channels with IC₅₀ = 11.4 nM. Here we reveal the structure of LqhIII bound to Na_V1.5 at 3.3 Å resolution by cryo-EM. LqhIII anchors on top of voltage-sensing domain IV, wedged between the S1-S2 and S3-S4 linkers, which traps the gating charges of the S4 segment in a unique intermediate-activated state stabilized by four ion-pairs. This conformational change is propagated inward to weaken binding of the fast inactivation gate and favor opening the activation gate. However, these changes do not permit Na⁺ permeation, revealing why LqhIII slows inactivation of Na_V channels but does not open them. Our results provide important insights into the structural basis for gating-modifier toxin binding, voltage-sensor trapping, and fast inactivation of Na_V channels.

Origin of the Phosphoprotein Phosphatase (PPP) sequence family in Bacteria: Critical ancestral sequence changes, radiation patterns and substrate binding features

Background

Phosphoprotein phosphatases (PPP) belong to the PPP Sequence family, which in turn belongs to the broader metallophosphoesterase (MPE) superfamily. The relationship between the PPP Sequence family and other members of the MPE superfamily remains unresolved, in particular what transitions took place in an ancestral MPE to ultimately produce the phosphoprotein specific phosphatases (PPPs).

Methods

We use structural and sequence alignment data, phylogenetic tree analysis, sequence signature (Weblogo) analysis, in silico protein-peptide modeling data, and in silico mutagenesis to trace a likely route of evolution from MPEs to the PPP Sequence family. Hidden Markov Model (HMM) based iterative database search strategies were utilized to identify PPP Sequence Family members from numerous bacterial groups.

Results

Using Mre11 as proxy for an ancestral nuclease-like MPE we trace a possible evolutionary route that alters a single active site substrate binding His-residue to yield a new substrate binding accessory, the "2-Arg-Clamp". The 2-Arg-Clamp is not found in MPEs, but is present in all PPP Sequence family members, where the phosphomonesterase reaction predominates. Variation in position of the clamp arginines and a supplemental sequence loop likely provide substrate specificity for each PPP Sequence family group.

Conclusions

Loss of a key substrate binding His-in MPEs opened the path to bind novel substrates and evolution of the 2-Arg-Clamp, a sequence change seen in both bacterial and eukaryotic phosphoprotein phosphatases.General significance: We establish a likely evolutionary route from nuclease-like MPE to PPP Sequence family enzymes, that includes the phosphoprotein phosphatases.

Allosteric Coupling Between Drug Binding and the Aromatic Cassette in the Pore Domain of the hERG1 Channel: Implications for a State-Dependent Blockade

Human-ether-a-go-go-related channel (hERG1) is the pore-forming domain of the delayed rectifier K⁺ channel in the heart which underlies the I_Kr current. The channel has been extensively studied due to its propensity to bind chemically diverse group of drugs. The subsequent hERG1 block can lead to a prolongation of the QT interval potentially leading to an abnormal cardiac electrical activity. The recently solved cryo-EM structure featured a striking non-swapped topology of the Voltage-Sensor Domain (VSD) which is packed against the pore-domain as well as a small and hydrophobic intra-cavity space. The small size and hydrophobicity of the cavity was unexpected and challenges the already-established hypothesis of drugs binding to the wide cavity. Recently, we showed that an amphipathic drug, ivabradine, may favorably bind the channel from the lipid-facing surface and we discovered a mutant (M651T) on the lipid facing domain between the VSD and the PD which inhibited the blocking capacity of the drug. Using multi-microseconds Molecular Dynamics (MD) simulations of wild-type and M651T mutant hERG1, we suggested the block of the channel through the lipid mediated pathway, the opening of which is facilitated by the flexible phenylalanine ring (F656). In this study, we characterize the dynamic interaction of the methionine-aromatic cassette in the S5-S6 helices by combining data from electrophysiological experiments with MD simulations and molecular docking to elucidate the complex allosteric coupling between drug binding to lipid-facing and intra-cavity sites and aromatic cassette dynamics. We investigated two well-established hERG1 blockers (ivabradine and dofetilide) for M651 sensitivity through electrophysiology and mutagenesis techniques. Our electrophysiology data reveal insensitivity of dofetilide to the mutations at site M651 on the lipid facing side of the channel, mirroring our results obtained from docking experiments. Moreover, we show that the dofetilide-induced block of hERG1 occurs through the intracellular space, whereas little to no block of ivabradine is observed during the intracellular application of the drug. The dynamic conformational rearrangement of the F656 appears to regulate the translocation of ivabradine into the central cavity. M651T mutation appears to disrupt this entry pathway by altering the molecular conformation of F656.

Mechanism of homodimeric cytokine receptor activation and dysregulation by oncogenic mutations

Homodimeric class I cytokine receptors are assumed to exist as preformed dimers that are activated by ligand-induced conformational changes. We quantified the dimerization of three prototypic class I cytokine receptors in the plasma membrane of living cells by single-molecule fluorescence microscopy. Spatial and spatiotemporal correlation of individual receptor subunits showed ligand-induced dimerization and revealed that the associated Janus kinase 2 (JAK2) dimerizes through its pseudokinase domain. Oncogenic receptor and hyperactive JAK2 mutants promoted ligand-independent dimerization, highlighting the formation of receptor dimers as the switch responsible for signal activation. Atomistic modeling and molecular dynamics simulations based on a detailed energetic analysis of the interactions involved in dimerization yielded a mechanistic blueprint for homodimeric class I cytokine receptor activation and its dysregulation by individual mutations.

Ion Binding and Selectivity of the Na+/H+ Antiporter MjNhaP1 from Experiment and Simulation

Cells employ membrane-embedded antiporter proteins to control their pH, salt concentration, and volume. The large family of cation/proton antiporters is dominated by Na⁺/H⁺ antiporters that exchange sodium ions against protons, but homologous K⁺/H⁺ exchangers have recently been characterized. We show experimentally that the electroneutral antiporter NhaP1 of Methanocaldococcus jannaschii (MjNhaP1) is highly selective for Na⁺ ions. We then characterize the ion selectivity in both the inward-open and outward-open states of MjNhaP1 using classical molecular dynamics simulations, free energy calculations, and hybrid quantum/classical (QM/MM) simulations. We show that MjNhaP1 is highly selective for binding of Na⁺ over K⁺ in the inward-open state, yet it is only weakly selective in the outward-open state. These findings are consistent with the function of MjNhaP1 as a sodium-driven deacidifier of the cytosol that maintains a high cytosolic K⁺ concentration in environments of high salinity. By combining experiment and computation, we gain mechanistic insight into the Na⁺/H⁺ transport mechanism and help elucidate the molecular basis for ion selectivity in cation/proton exchangers.

The Pore-Lipid Interface: Role of Amino-Acid Determinants of Lipophilic Access by Ivabradine to the hERG1 Pore Domain

Abnormal cardiac electrical activity is a common side effect caused by unintended block of the promiscuous drug target human ether-à-go-go-related gene (hERG1), the pore-forming domain of the delayed rectifier K⁺ channel in the heart. hERG1 block leads to a prolongation of the QT interval, a phase of the cardiac cycle that underlies myocyte repolarization detectable on the electrocardiogram. Even newly released drugs such as heart-rate lowering agent ivabradine block the rapid delayed rectifier current I_Kr, prolong action potential duration, and induce potentially lethal arrhythmia known as torsades de pointes. In this study, we describe a critical drug-binding pocket located at the lateral pore surface facing the cellular membrane. Mutations of the conserved M651 residue alter ivabradine-induced block but not by the common hERG1 blocker dofetilide. As revealed by molecular dynamics simulations, binding of ivabradine to a lipophilic pore access site is coupled to a state-dependent reorientation of aromatic residues F557 and F656 in the S5 and S6 helices. We show that the M651 mutation impedes state-dependent dynamics of F557 and F656 aromatic cassettes at the protein-lipid interface, which has a potential to disrupt drug-induced block of the channel. This fundamentally new mechanism coupling the channel dynamics and small-molecule access from the membrane into the hERG1 intracavitary site provides a simple rationale for the well established state-dependence of drug blockade. SIGNIFICANCE STATEMENT: The drug interference with the function of the cardiac hERG channels represents one of the major sources of drug-induced heart disturbances. We found a novel and a critical drug-binding pocket adjacent to a lipid-facing surface of the hERG1 channel, which furthers our molecular understanding of drug-induced QT syndrome.

Isothermal Titration Calorimetry of Be2+ with Phosphatidylserine Models Guides All-Atom Force-Field Development for Lipid-Ion Interactions

Beryllium has multiple industrial applications but exposure to its dust during manufacturing is associated with developing chronic inflammation in lungs known as berylliosis. Besides binding to specific alleles of MHC-II, Be²⁺ was recently found to compete with Ca²⁺ for binding sites on phosphatidylserine-containing membranes and inhibit recognition of this lipid by phagocytes. Computational studies of possible molecular targets for this small toxic dication are impeded by the absence of a reliable force field. This study introduces parameters for Be²⁺ for the CHARMM36 additive force field that represent interactions with water, including free energy of hydration and ion-monohydrate interaction energy and separation distance; and interaction parameters describing Be²⁺ affinity for divalent ion binding sites on lipids, namely phosphoryl and carboxylate oxygens. Results from isothermal titration calorimetry experiments for the binding affinities of Be²⁺ to dimethyl phosphate and acetate ions reveal that Be²⁺ strongly binds to phosphoryl groups. Revised interaction parameters for Be²⁺ with these types of oxygens reproduce experimental affinities in solution simulations. Surface tensions calculated from simulations of DOPS monolayers with varied concentrations of Be²⁺ are compared with prior results from Langmuir monolayer experiments, verifying the compacting effect that produces greater surface tensions (lower pressures) for Be²⁺-bound monolayers at the same surface area in comparison with K⁺. The new parameters will enable simulations that should reveal the mechanism of Be²⁺ interference with molecular recognition and signaling processes.

Structure of voltage-dependent anion channel-tethered bilayer lipid membranes determined using neutron reflectivity

Neutron reflectivity (NR) has emerged as a powerful technique to study the structure and behavior of membrane proteins at planar lipid interfaces. Integral membrane proteins (IMPs) remain a significant challenge for NR owing to the difficulty of forming complete bilayers with sufficient protein density for scattering techniques. One strategy to achieve high protein density on a solid substrate is the capture of detergent-stabilized, affinity-tagged IMPs on a nitrilotriacetic acid (NTA)-functionalized self-assembled monolayer (SAM), followed by reconstitution into the lipids of interest. Such protein-tethered bilayer lipid membranes (ptBLMs) have the notable advantage of a uniform IMP orientation on the substrate. Here, NR is used to provide a structural characterization of the ptBLM process from formation of the SAM to capture of the detergent-stabilized IMP and lipid reconstitution. The mitochondrial outer-membrane voltage-dependent anion channel (VDAC), which controls the exchange of bioenergetic metabolites between mitochondria and the cytosol, was used as a model β-barrel IMP. Molecular dynamics simulations were used for comparison with the experimental results and to inform the parameters of the physical models describing the NR data. The detailed structure of the SAM is shown to depend on the density of the NTA chelating groups. The relative content of detergent and protein in surface-immobilized, detergent-stabilized VDAC is measured, while the reconstituted lipid bilayer is shown to be complete to within a few percent, using the known atomic structure of VDAC. Finally, excess lipid above the reconstituted bilayer, which is of consequence for more indirect structural and functional studies, is shown to be present.

Molecular Dynamics of Ion Conduction through the Selectivity Filter of the NaVAb Sodium Channel

The determination of the atomic structures of voltage-gated bacterial sodium channels using X-ray crystallography has provided a first view of this family of membrane proteins. Molecular dynamics simulations offer one approach to clarify the underlying mechanism of permeation and selectivity in these channels. However, it appears that the intracellular gate of the pore domain is either closed or only open partially in the available X-ray structures. The lack of structure with a fully open intracellular gate poses a special challenge to computational studies aimed at simulating ion conduction. To circumvent this problem, we simulated a model of the Na_VAb channel in which the transmembrane S5 and S6 helices of the pore domain have been truncated to provide direct open access to the intracellular entryway to the pore. Molecular dynamics simulations were carried out over a range of membrane potential and ion concentration of sodium and potassium. The simulations show that the Na_VAb selectivity filter is essentially a cationic pore supporting the conduction of ions at a rate comparable to aqueous diffusion with no significant selectivity for sodium. Conductance and selectivity vary as a function of ion concentration for both cations. Permeation occurs primarily via a knock-on mechanism for both sodium and potassium, although the ion ordering in single file along the pore is not strictly maintained. The character of the outward current appears quite different from the inward current, with a buildup on ions in the selectivity filter prior to escape toward the extracellular side, indicating the presence of a rectification effect that is overcome by nonphysiological applied voltages.

Structural basis for gating pore current in periodic paralysis

Potassium-sensitive Hypokalemic and Normokalemic Periodic Paralysis (HypoPP, NormoPP) are inherited skeletal muscle diseases characterized by episodes of flaccid muscle weakness^1,2. They are caused by mutations in one gating charge in an S4 transmembrane segment in the voltage sensor (VS) of voltage-gated sodium channel Na_v1.4 or calcium channel Ca_v1.1^1,2. Mutations of the outermost arginine gating charges (R1 and R2) cause HypoPP^1,2 by creating a pathogenic gating pore in the VS through which cations leak in the resting state^3,4. Mutations of the third arginine gating charge (R3) cause NormoPP⁵ owing to cationic leak in activated/inactivated states⁶. Here we present high-resolution structures of these pathogenic gating pores in the model bacterial sodium channel Na_VAb^7,8. Mutation of R2 in Na_VAb gives gating pore current in resting states, whereas mutation of R3 gives gating pore current in activated/inactivated states. Mutations R2G and R3G have no effect on backbone structures of VS, but create aqueous space near the hydrophobic constriction site (HCS) that controls gating charge movement through VS. The R3G mutation extends the extracellular aqueous cleft completely through the activated VS. Although the R2G mutation does not create a continuous aqueous pathway in the activated state, molecular modeling of the resting state reveals a complete water-accessible pathway. Crystal structures of Na_VAb/R2G in complex with guanidinium define a potential drug target site. Molecular dynamics simulations illustrate the mechanism of Na⁺ permeation through the mutant gating pore in concert with conformational fluctuations of gating charge R4. Our results reveal pathogenic mechanisms of periodic paralysis at the atomic level and suggest designs of drugs that may prevent ionic leak and provide symptomatic relief from these episodic diseases.

Molecular dynamics simulations of lipid nanodiscs

A lipid nanodisc is a discoidal lipid bilayer stabilized by proteins, peptides, or polymers on its edge. Nanodiscs have two important connections to structural biology. The first is associated with high-density lipoprotein (HDL), a particle with a variety of functionalities including lipid transport. Nascent HDL (nHDL) is a nanodisc stabilized by Apolipoprotein A-I (APOA1). Determining the structure of APOA1 and its mimetic peptides in nanodiscs is crucial to understanding pathologies related to HDL maturation and designing effective therapies. Secondly, nanodiscs offer non-detergent membrane-mimicking environments and greatly facilitate structural studies of membrane proteins. Although seemingly similar, natural and synthetic nanodiscs are different in that nHDL is heterogeneous in size, due to APOA1 elasticity, and gradually matures to become spherical. Synthetic nanodiscs, in contrast, should be homogenous, stable, and size-tunable. This report reviews previous molecular dynamics (MD) simulation studies of nanodiscs and illustrates convergence and accuracy issues using results from new multi-microsecond atomistic MD simulations. These new simulations reveal that APOA1 helices take 10-20 μs to rearrange on the nanodisc, while peptides take 2 μs to migrate from the disc surfaces to the edge. These systems can also become kinetically trapped depending on the initial conditions. For example, APOA1 was trapped in a biologically irrelevant conformation for the duration of a 10 μs trajectory; the peptides were similarly trapped for 5 μs. It therefore remains essential to validate MD simulations of these systems with experiments due to convergence and accuracy issues. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Competition between Li+ and Na+ in sodium transporters and receptors: Which Na+-Binding sites are "therapeutic" Li+ targets?

Li⁺ (turquoise), the better charge acceptor, can displace Na⁺ (purple) bound by only one or two aa residues in buried sites. Thus, Li⁺ can displace Na⁺ bound by Asp^– and Ser in the A_2AAR/β₁AR receptor and enhance the metal site's stability, thus prohibiting structural distortions induced by agonist binding, leading to lower cytosolic levels of activated G-proteins, which are hyperactive in bipolar disorder patients.

Sodium (Na⁺) acts as an indispensable allosteric regulator of the activities of biologically important neurotransmitter transporters and G-protein coupled receptors (GPCRs), which comprise well-known drug targets for psychiatric disorders and addictive behavior. How selective these allosteric Na⁺-binding sites are for the cognate cation over abiogenic Li⁺, a first-line drug to treat bipolar disorder, is unclear. Here, we reveal how properties of the host protein and its binding cavity affect the outcome of the competition between Li⁺ and Na⁺ for allosteric binding sites in sodium transporters and receptors. We show that rigid Na⁺-sites that are crowded with multiple protein ligands are well-protected against Li⁺ attack, but their flexible counterparts or buried Na⁺-sites containing only one or two protein ligands are vulnerable to Li⁺ substitution. These findings suggest a novel possible mode of Li⁺ therapeutic action: By displacing Na⁺ bound by ≤2 protein ligands in buried GPCR sites and stabilizing the receptor's inactive state, Li⁺ could prohibit conformational changes to an active state, leading to lower cytosolic levels of activated guanine nucleotide-binding proteins, which are hyperactive/overexpressed in bipolar disorder patients.

Structural properties of apolipoprotein A-I mimetic peptides that promote ABCA1-dependent cholesterol efflux

Peptides mimicking the major protein of highdensity lipoprotein (HDL), apolipoprotein A-I (apoA-I), are promising therapeutics for cardiovascular diseases. Similar to apoA-I, their atheroprotective property is attributed to their ability to form discoidal HDL-like particles by extracting cellular cholesterol and phospholipids from lipid microdomains created by the ABCA1 transporter in a process called cholesterol efflux. The structural features of peptides that enable cholesterol efflux are not well understood. Herein, four synthetic amphipathic peptides denoted ELK, which only contain Glu, Leu, Lys, and sometimes Ala, and which have a wide range of net charges and hydrophobicities, were examined for cholesterol efflux. Experiments show that ELKs with a net neutral charge and a hydrophobic face that subtends an angle of at least 140° are optimal for cholesterol efflux. All-atom molecular dynamics simulations show that peptides that are effective in promoting cholesterol efflux stabilize HDL nanodiscs formed by these peptides by the orderly covering of the hydrophobic acyl chains on the edge of the disc. In contrast to apoA-I, which forms an anti-parallel double belt around the HDL, active peptides assemble in a mostly anti-parallel “picket fence” arrangement. These results shed light on the efflux ability of apoA-I mimetics and inform the future design of such therapeutics.

A kinetic account for amphetamine-induced monoamine release

This study on serotonin transporters shows that amphetamine-induced monoamine release requires cooperative substrate and cosubstrate binding. A kinetic model is presented that can account qualitatively and quantitatively for the releasing action of amphetamines.

The plasmalemmal monoamine transporters for dopamine, norepinephrine, and serotonin (SERT) are targets for amphetamines. In vivo, amphetamines elicit most, if not all, of their actions by triggering monoamine efflux. This is thought to be accomplished by an amphetamine-induced switch from the forward-transport to the substrate-exchange mode. The mechanism underlying this switch has remained elusive; available kinetic models posit that substrates and cosubstrate Na⁺ ions bind either in a random or in a sequential order. Neither can account for all reported experimental observations. We used electrophysiological recordings to interrogate crucial conformational transitions associated with the binding of five different substrates (serotonin, para-chloroamphetamine, and the high-affinity naphthyl-propan-amines PAL-287, PAL-1045, and PAL-1046) to human SERT expressed in HEK293 cells; specifically, we determined the relaxation kinetics of SERT from a substrate-loaded to a substrate-free state at various intracellular and extracellular Na⁺ concentrations. These rates and their dependence on intracellular and extracellular Na⁺ concentrations differed considerably between substrates. We also examined the effect of K⁺ on substrate affinity and found that K⁺ enhanced substrate dissociation. A kinetic model was developed, which allowed for random, but cooperative, binding of substrate and Na⁺ (or K⁺). The synthetic data generated by this model recapitulated the experimental observations. More importantly, the cooperative binding model accounted for the releasing action of amphetamines without any digression from alternating access. To the best of our knowledge, this model is the first to provide a mechanistic framework for amphetamine-induced monoamine release and to account for the findings that some substrates are less efficacious than others in promoting the substrate-exchange mode.

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.