The Ca2+ permeation mechanism of the ryanodine receptor revealed by a multi-site ion model

The Effect of Lipopolysaccharides on the Electrostatic Properties of Gram-Negative General Porins from Enterobacteriaceae

We investigated, by using all-atom molecular dynamics simulations, the effect of the outer membrane of Gram-negative bacteria, composed in the outer leaflet by polar/charged lipopolysaccharides (LPS), on the electrostatic properties of general porins from the Enterobacteriaceae family. General porins constitute the main path for the facilitated diffusion of polar antibiotics through the outer membrane. As model system we selected OmpK36 from Klebsiella pneumoniae, the ortholog of OmpC from Escherichia coli. This species presents high variability of amino acid composition of porins, with the effect to increase its resistance to the penetration of antibiotics. The various properties we analyzed seem to indicate that LPS acts as an independent layer without affecting the internal electrostatic properties of OmpK36. The only apparent effect on the microsecond time scale we sampled is the appearance of calcium ions, when present at moderate concentration in solution, inside the pore. However, we noticed increased fluctuations of the polarization density and only minor changes on its average value.

Electrophysiological properties and structural prediction of the SARS-CoV-2 viroprotein E

COVID-19, the infectious disease caused by the most recently discovered coronavirus SARS- CoV-2, has caused millions of sick people and thousands of deaths all over the world. The viral positive-sense single-stranded RNA encodes 31 proteins among which the spike (S) is undoubtedly the best known. Recently, protein E has been reputed as a potential pharmacological target as well. It is essential for the assembly and release of the virions in the cell. Literature describes protein E as a voltage-dependent channel with preference towards monovalent cations whose intracellular expression, though, alters Ca2+ homeostasis and promotes the activation of the proinflammatory cascades. Due to the extremely high sequence identity of SARS-CoV-2 protein E (E-2) with the previously characterized E-1 (i.e., protein E from SARS-CoV) many data obtained for E-1 were simply adapted to the other. Recent solid state NMR structure revealed that the transmembrane domain (TMD) of E-2 self-assembles into a homo-pentamer, albeit the oligomeric status has not been validated with the full-length protein. Prompted by the lack of a common agreement on the proper structural and functional features of E-2, we investigated the specific mechanism/s of pore-gating and the detailed molecular structure of the most cryptic protein of SARS-CoV-2 by means of MD simulations of the E-2 structure and by expressing, refolding and analyzing the electrophysiological activity of the transmembrane moiety of the protein E-2, in its full length. Our results show a clear agreement between experimental and predictive studies and foresee a mechanism of activity based on Ca2+ affinity.

Structural characterization and functional insights into the type II secretion system of the poly-extremophile Deinococcus radiodurans

The extremophile bacterium D. radiodurans boasts a distinctive cell envelope characterized by the regular arrangement of three protein complexes. Among these, the Type II Secretion System (T2SS) stands out as a pivotal structural component. We used cryo-electron microscopy to reveal unique features, such as an unconventional protein belt (DR_1364) around the main secretin (GspD), and a cap (DR_0940) found to be a separated subunit rather than integrated with GspD. Furthermore, a novel region at the N-terminus of the GspD constitutes an additional second gate, supplementing the one typically found in the outer membrane region. This T2SS was found to contribute to envelope integrity, while also playing a role in nucleic acid and nutrient trafficking. Studies on intact cell envelopes show a consistent T2SS structure repetition, highlighting its significance within the cellular framework.

Calcium binding and permeation in TRPV channels: Insights from molecular dynamics simulations

Some calcium channels selectively permeate Ca2+, despite the high concentration of monovalent ions in the surrounding environment, which is essential for many physiological processes. Without atomistic and dynamical ion permeation details, the underlying mechanism of Ca2+ selectivity has long been an intensively studied, yet controversial, topic. This study takes advantage of the homologous Ca2+-selective TRPV6 and non-selective TRPV1 and utilizes the recently solved open-state structures and a newly developed multisite calcium model to investigate the ion binding and permeation features in TRPV channels by molecular dynamics simulations. Our results revealed that the open-state TRPV6 and TRPV1 show distinct ion binding patterns in the selectivity filter, which lead to different ion permeation features. Two Ca2+ ions simultaneously bind to the selectivity filter of TRPV6 compared with only one Ca2+ in the case of TRPV1. Multiple Ca2+ binding at the selectivity filter of TRPV6 permeated in a concerted manner, which could efficiently block the permeation of Na+. Cations of various valences differentiate between the binding sites at the entrance of the selectivity filter in TRPV6. Ca2+ preferentially binds to the central site with a higher probability of permeation, repelling Na+ to a peripheral site. Therefore, we believe that ion binding competition at the selectivity filter of calcium channels, including the binding strength and number of binding sites, determines Ca2+ selectivity under physiological conditions.

Structural basis for human Cav1.2 inhibition by multiple drugs and the neurotoxin calciseptine

Cav1.2 channels play crucial roles in various neuronal and physiological processes. Here, we present cryo-EM structures of human Cav1.2, both in its apo form and in complex with several drugs, as well as the peptide neurotoxin calciseptine. Most structures, apo or bound to calciseptine, amlodipine, or a combination of amiodarone and sofosbuvir, exhibit a consistent inactivated conformation with a sealed gate, three up voltage-sensing domains (VSDs), and a down VSDII. Calciseptine sits on the shoulder of the pore domain, away from the permeation path. In contrast, when pinaverium bromide, an antispasmodic drug, is inserted into a cavity reminiscent of the IFM-binding site in Nav channels, a series of structural changes occur, including upward movement of VSDII coupled with dilation of the selectivity filter and its surrounding segments in repeat III. Meanwhile, S4-5III merges with S5III to become a single helix, resulting in a widened but still non-conductive intracellular gate.

Computed Protein-Protein Enthalpy Signatures as a Tool for Identifying Conformation Sampling Problems

Understanding the thermodynamic signature of protein-peptide binding events is a major challenge in computational chemistry. The complexity generated by both components possessing many degrees of freedom poses a significant issue for methods that attempt to directly compute the enthalpic contribution to binding. Indeed, the prevailing assumption has been that the errors associated with such approaches would be too large for them to be meaningful. Nevertheless, we currently have no indication of how well the present methods would perform in terms of predicting the enthalpy of binding for protein-peptide complexes. To that end, we carefully assembled and curated a set of 11 protein-peptide complexes where there is structural and isothermal titration calorimetry data available and then computed the absolute enthalpy of binding. The initial "out of the box" calculations were, as expected, very modest in terms of agreement with the experiment. However, careful inspection of the outliers allows for the identification of key sampling problems such as distinct conformations of peptide termini not being sampled or suboptimal cofactor parameters. Additional simulations guided by these aspects can lead to a respectable correlation with isothermal titration calorimetry (ITC) experiments (R2 of 0.88 and an RMSE of 1.48 kcal/mol overall). Although one cannot know prospectively whether computed ITC values will be correct or not, this work shows that if experimental ITC data are available, then this in conjunction with computed ITC, can be used as a tool to know if the ensemble being simulated is representative of the true ensemble or not. That is important for allowing the correct interpretation of the detailed dynamics of the system with respect to the measured enthalpy. The results also suggest that computational calorimetry is becoming increasingly feasible. We provide the data set as a resource for the community, which could be used as a benchmark to help further progress in this area.

A systematic analysis of residue and risk of cyantraniliprole in the water-sediment system: Does metabolism reduce its environmental risk?

As a representative variety of diamide insecticides, cyantraniliprole has broad application prospects. In this study, the fate and risk of cyantraniliprole and its main metabolite J9Z38 in a water-sediment system were investigated. The present result showed that more J9Z38 was adsorbed in the sediment at the end of exposure. However, the bioaccumulation capacity of cyantraniliprole in zebrafish was higher than that of J9Z38. Cyantraniliprole had stronger influence on the antioxidant system and detoxification system of zebrafish than J9Z38. Moreover, cyantraniliprole induced more significant oxidative stress effect and more differentially expressed genes (DEGs) in zebrafish. Cyantraniliprole had significantly influence on the expression of RyR-receptor-related genes, which was confirmed by resolving their binding modes with key receptor proteins using AlphaFold2 and molecular docking techniques. In the sediment, both cyantraniliprole and J9Z38 had inhibitory effects on microbial community structure diversity and metabolic function, especially cyantraniliprole. The methane metabolism pathway, mediated by methanogens such as Methanolinea, Methanoregula, and Methanosaeta, may be the main pathway of degradation of cyantraniliprole and J9Z38 in sediments. The present results demonstrated that metabolism can reduce the environmental risk of cyantraniliprole in water-sediment system to a certain extent.

Fluoride permeation mechanism of the Fluc channel in liposomes revealed by solid-state NMR

Solid-state nuclear magnetic resonance (ssNMR) methods can probe the motions of membrane proteins in liposomes at the atomic level and propel the understanding of biomolecular processes for which static structures cannot provide a satisfactory description. In this work, we report our study on the fluoride channel Fluc-Ec1 in phospholipid bilayers based on ssNMR and molecular dynamics simulations. Previously unidentified fluoride binding sites in the aqueous vestibules were experimentally verified by 19F-detected ssNMR. One of the two fluoride binding sites in the polar track was identified as a water molecule by 1H-detected ssNMR. Meanwhile, a dynamic hotspot at loop 1 was observed by comparing the spectra of wild-type Fluc-Ec1 in variant buffer conditions or with its mutants. Therefore, we propose that fluoride conduction in the Fluc channel occurs via a "water-mediated knock-on" permeation mechanism and that loop 1 is a key molecular determinant for channel gating.

Curvature Matters: Modeling Calcium Binding to Neutral and Anionic Phospholipid Bilayers

In this work, the influence of membrane curvature on the Ca²⁺ binding to phospholipid bilayers is investigated by means of molecular dynamics simulations. In particular, we compared Ca²⁺ binding to flat, elastically buckled, or uniformly bent zwitterionic and anionic phospholipid bilayers. We demonstrate that Ca²⁺ ions bind preferably to the concave membrane surfaces in both types of bilayers. We also show that the membrane curvature leads to pronounced changes in Ca²⁺ binding including differences in free ion concentrations, lipid coordination distributions, and the patterns of ion binding to different chemical groups of lipids. Moreover, these effects differ substantially for the concave and convex membrane monolayers. Comparison between force fields with either full or scaled charges indicates that charge scaling results in reduction of the Ca²⁺ binding to curved phosphatidylserine bilayers, while for phosphatidylcholine membranes, calcium binds only weakly for both force fields.

A cooperative knock-on mechanism underpins Ca2+-selective cation permeation in TRPV channels

The selective exchange of ions across cellular membranes is a vital biological process. Ca2+-mediated signaling is implicated in a broad array of physiological processes in cells, while elevated intracellular concentrations of Ca2+ are cytotoxic. Due to the significance of this cation, strict Ca2+ concentration gradients are maintained across the plasma and organelle membranes. Therefore, Ca2+ signaling relies on permeation through selective ion channels that control the flux of Ca2+ ions. A key family of Ca2+-permeable membrane channels is the polymodal signal-detecting transient receptor potential (TRP) ion channels. TRP channels are activated by a wide variety of cues including temperature, small molecules, transmembrane voltage, and mechanical stimuli. While most members of this family permeate a broad range of cations non-selectively, TRPV5 and TRPV6 are unique due to their strong Ca2+ selectivity. Here, we address the question of how some members of the TRPV subfamily show a high degree of Ca2+ selectivity while others conduct a wider spectrum of cations. We present results from all-atom molecular dynamics simulations of ion permeation through two Ca2+-selective and two non-selective TRPV channels. Using a new method to quantify permeation cooperativity based on mutual information, we show that Ca2+-selective TRPV channel permeation occurs by a three-binding site knock-on mechanism, whereas a two-binding site knock-on mechanism is observed in non-selective TRPV channels. Each of the ion binding sites involved displayed greater affinity for Ca2+ over Na+. As such, our results suggest that coupling to an extra binding site in the Ca2+-selective TRPV channels underpins their increased selectivity for Ca2+ over Na+ ions. Furthermore, analysis of all available TRPV channel structures shows that the selectivity filter entrance region is wider for the non-selective TRPV channels, slightly destabilizing ion binding at this site, which is likely to underlie mechanistic decoupling.

Cryo-EM structures of human SPCA1a reveal the mechanism of Ca2+/Mn2+ transport into the Golgi apparatus

Secretory pathway Ca²⁺/Mn²⁺ ATPase 1 (SPCA1) actively transports cytosolic Ca²⁺ and Mn²⁺ into the Golgi lumen, playing a crucial role in cellular calcium and manganese homeostasis. Detrimental mutations of the ATP2C1 gene encoding SPCA1 cause Hailey-Hailey disease. Here, using nanobody/megabody technologies, we determined cryo-electron microscopy structures of human SPCA1a in the ATP and Ca²⁺/Mn²⁺-bound (E1-ATP) state and the metal-free phosphorylated (E2P) state at 3.1- to 3.3-Å resolutions. The structures revealed that Ca²⁺ and Mn²⁺ share the same metal ion-binding pocket with similar but notably different coordination geometries in the transmembrane domain, corresponding to the second Ca²⁺-binding site in sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA). In the E1-ATP to E2P transition, SPCA1a undergoes similar domain rearrangements to those of SERCA. Meanwhile, SPCA1a shows larger conformational and positional flexibility of the second and sixth transmembrane helices, possibly explaining its wider metal ion specificity. These structural findings illuminate the unique mechanisms of SPCA1a-mediated Ca²⁺/Mn²⁺ transport.

Mechanism of Calcium Permeation in a Glutamate Receptor Ion Channel

The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) are neurotransmitter-activated cation channels ubiquitously expressed in vertebrate brains. The regulation of calcium flux through the channel pore by RNA-editing is linked to synaptic plasticity while excessive calcium influx poses a risk for neurodegeneration. Unfortunately, the molecular mechanisms underlying this key process are mostly unknown. Here, we investigated calcium conduction in calcium-permeable AMPAR using Molecular Dynamics (MD) simulations with recently introduced multisite force-field parameters for Ca²⁺. Our calculations are consistent with experiment and explain the distinct calcium permeability in different RNA-edited forms of GluA2. For one of the identified metal binding sites, multiscale Quantum Mechanics/Molecular Mechanics (QM/MM) simulations further validated the results from MD and revealed small but reproducible charge transfer between the metal ion and its first solvation shell. In addition, the ion occupancy derived from MD simulations independently reproduced the Ca²⁺ binding profile in an X-ray structure of an NaK channel mimicking the AMPAR selectivity filter. This integrated study comprising X-ray crystallography, multisite MD, and multiscale QM/MM simulations provides unprecedented insights into Ca²⁺ permeation mechanisms in AMPARs, and paves the way for studying other biological processes in which Ca²⁺ plays a pivotal role.

Conformational dynamics and putative substrate extrusion pathways of the N-glycosylated outer membrane factor CmeC from Campylobacter jejuni

The outer membrane factor CmeC of the efflux machinery CmeABC plays an important role in conferring antibiotic and bile resistance to Campylobacter jejuni. Curiously, the protein is N-glycosylated, with the glycans playing a key role in the effective function of this system. In this work we have employed atomistic equilibrium molecular dynamics simulations of CmeC in a representative model of the C. jejuni outer membrane to characterise the dynamics of the protein and its associated glycans. We show that the glycans are more conformationally labile than had previously been thought. The extracellular loops of CmeC visit the open and closed states freely suggesting the absence of a gating mechanism on this side, while the narrow periplasmic entrance remains tightly closed, regulated via coordination to solvated cations. We identify several cation binding sites on the interior surface of the protein. Additionally, we used steered molecular dynamics simulations to elucidate translocation pathways for a bile acid and a macrolide antibiotic. These, and additional equilibrium simulations suggest that the anionic bile acid utilises multivalent cations to climb the ladder of acidic residues that line the interior surface of the protein.

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

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

Calcium Ion Binding to the Mutants of Calmodulin: A Structure-Based Computational Predictive Model of Binding Affinity Using a Charge Scaling Approach in Molecular Dynamics Simulation

The binding of calcium ions (Ca²⁺) to the calcium-binding proteins (CBPs) controls a plethora of regulatory processes. Among the roles played by CBPs in several diseases, the onset and progress of some cardiovascular diseases are caused by mutations in calmodulin (CaM), an important member of CBPs. Rationalization and prediction of the binding affinity of Ca²⁺ ions to the CaM can play important roles in understanding the origin of cardiovascular diseases. However, there is no robust structure-based computational method for predicting the binding affinity of Ca²⁺ ions to the different forms of CBPs in general and CaM in particular. In the current work, we have devised a fast yet accurate computational technique to accurately calculate the binding affinity of Ca²⁺ to the different forms of CaM. This method combines the well-known molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method and a charge scaling approach developed by previous authors that takes care of the polarization of CaM and Ca²⁺ ions. Our detailed analysis of the different components of binding free energy shows that subtle changes in electrostatics and van der Waals contribute to the difference in the binding affinity of mutants from that of the wild type (WT), and the charge scaling approach is superior in calculating these subtle changes in electrostatics as compared to the nonpolarizable force field used in this work. A statistically significant regression model made from our binding free energy calculations gives a correlation coefficient close to 0.8 to the experimental results. This structure-based predictive model can open up a new strategy to understand and predict the binding of Ca²⁺ to the mutants of CBPs, in general.

Binding Energy and Free Energy of Calcium Ion to Calmodulin EF-Hands with the Drude Polarizable Force Field

Calcium ions are important messenger molecules in cells, which bind calcium-binding proteins to trigger many biochemical processes. We constructed four model systems, each containing one EF-hand loop of calmodulin with one calcium ion bound, and investigated the binding energy and free energy of Ca²⁺ by the quantum mechanics symmetry-adapted perturbation theory (SAPT) method and the molecular mechanics with the additive CHARMM36m (C36m) and the polarizable Drude force fields (FFs). Our results show that the explicit introduction of polarizability in the Drude not only yields considerably improved agreement with the binding energy calculated from the SAPT method but is also able to capture each component of the binding energies including electrostatic, induction, exchange, and dispersion terms. However, binding free energies computed with the Drude and the C36m FFs both deviated significantly from the experimental measurements. Detailed analysis indicated that one of main reasons might be that the strong interactions between Ca²⁺ and the side chain nitrogen of Asn/Gln in the Drude FF caused the distorted coordination geometries of calcium. Our work illustrated the importance of polarization in modeling ion-protein interactions and the difficulty in generating accurate and balanced FF models to represent the polarization effects.

The Impact of Mutation L138F/L210F on the Orai Channel: A Molecular Dynamics Simulation Study

The calcium release-activated calcium channel, composed of the Orai channel and the STIM protein, plays a crucial role in maintaining the Ca²⁺ concentration in cells. Previous studies showed that the L138F mutation in the human Orai1 creates a constitutively open channel independent of STIM, causing severe myopathy, but how the L138F mutation activates Orai1 is still unclear. Here, based on the crystal structure of Drosophila melanogaster Orai (dOrai), molecular dynamics simulations for the wild-type (WT) and the L210F (corresponding to L138F in the human Orai1) mutant were conducted to investigate their structural and dynamical properties. The results showed that the L210F dOrai mutant tends to have a more hydrated hydrophobic region (V174 to F171), as well as more dilated basic region (K163 to R155) and selectivity filter (E178). Sodium ions were located deeper in the mutant than in the wild-type. Further analysis revealed two local but essential conformational changes that may be the key to the activation. A rotation of F210, a previously unobserved feature, was found to result in the opening of the K163 gate through hydrophobic interactions. At the same time, a counter-clockwise rotation of F171 occurred more frequently in the mutant, resulting in a wider hydrophobic gate with more hydration. Ultimately, the opening of the two gates may facilitate the opening of the Orai channel independent of STIM.

Structural and energetic analysis of metastable intermediate states in the E1P-E2P transition of Ca2+-ATPase

Ion pumps (or P-type ATPases) are membrane proteins, which transport ions through biological membranes against a concentration gradient, a function essential for many biological processes, such as muscle contraction, neurotransmission, and metabolism. Molecular mechanisms underlying active ion transport by ion pumps have been investigated by biochemical experiments and high-resolution structure analyses. Here, the transition of sarcoplasmic reticulum Ca²⁺-ATPase upon dissociation of Ca²⁺ is investigated using atomistic molecular dynamics simulations. We find intermediate structures along the pathway are stabilized by transient interactions between A- and P-domains as well as lipid molecules in the transmembrane helices.

Sarcoplasmic reticulum (SR) Ca²⁺-ATPase transports two Ca²⁺ ions from the cytoplasm to the SR lumen against a large concentration gradient. X-ray crystallography has revealed the atomic structures of the protein before and after the dissociation of Ca²⁺, while biochemical studies have suggested the existence of intermediate states in the transition between E1P⋅ADP⋅2Ca²⁺ and E2P. Here, we explore the pathway and free energy profile of the transition using atomistic molecular dynamics simulations with the mean-force string method and umbrella sampling. The simulations suggest that a series of structural changes accompany the ordered dissociation of ADP, the A-domain rotation, and the rearrangement of the transmembrane (TM) helices. The luminal gate then opens to release Ca²⁺ ions toward the SR lumen. Intermediate structures on the pathway are stabilized by transient sidechain interactions between the A- and P-domains. Lipid molecules between TM helices play a key role in the stabilization. Free energy profiles of the transition assuming different protonation states suggest rapid exchanges between Ca²⁺ ions and protons when the Ca²⁺ ions are released toward the SR lumen.

Mechanism of calcium potentiation of the α7 nicotinic acetylcholine receptor

The α7 nicotinic acetylcholine receptor (nAChR) is among the most abundant types of nAChR in the brain, yet the ability of nerve-released ACh to activate α7 remains enigmatic. In particular, a major population of α7 resides in extra-synaptic regions where the ACh concentration is reduced, owing to dilution and enzymatic hydrolysis, yet ACh shows low potency in activating α7. Using high-resolution single-channel recording techniques, we show that extracellular calcium is a powerful potentiator of α7 activated by low concentrations of ACh. Potentiation manifests as robust increases in the frequency of channel opening and the average duration of the openings. Molecular dynamics simulations reveal that calcium binds to the periphery of the five ligand binding sites and is framed by a pair of anionic residues from the principal and complementary faces of each site. Mutation of residues identified by simulation prevents calcium from potentiating ACh-elicited channel opening. An anionic residue is conserved at each of the identified positions in all vertebrate species of α7. Thus, calcium associates with a novel structural motif on α7 and is an obligate cofactor in regions of limited ACh concentration.

Nonselective cation permeation in an AMPA-type glutamate receptor

Fast excitatory synaptic transmission in the central nervous system relies on the AMPA-type glutamate receptor (AMPAR). This receptor incorporates a nonselective cation channel, which is opened by the binding of glutamate. Although the open pore structure has recently became available from cryo-electron microscopy (Cryo-EM), the molecular mechanisms governing cation permeability in AMPA receptors are not understood. Here, we combined microsecond molecular dynamic (MD) simulations on a putative open-state structure of GluA2 with electrophysiology on cloned channels to elucidate ion permeation mechanisms. Na⁺, K⁺, and Cs⁺ permeated at physiological rates, consistent with a structure that represents a true open state. A single major ion binding site for Na⁺ and K⁺ in the pore represents the simplest selectivity filter (SF) structure for any tetrameric cation channel of known structure. The minimal SF comprised only Q586 and Q587, and other residues on the cytoplasmic side formed a water-filled cavity with a cone shape that lacked major interactions with ions. We observed that Cl^- readily enters the upper pore, explaining anion permeation in the RNA-edited (Q586R) form of GluA2. A permissive architecture of the SF accommodated different alkali metals in distinct solvation states to allow rapid, nonselective cation permeation and copermeation by water. Simulations suggested Cs⁺ uses two equally populated ion binding sites in the filter, and we confirmed with electrophysiology of GluA2 that Cs⁺ is slightly more permeant than Na⁺, consistent with serial binding sites preferentially driving selectivity.

The structural basis of bacterial manganese import

Metal ions are essential for all forms of life. In prokaryotes, ATP-binding cassette (ABC) permeases serve as the primary import pathway for many micronutrients including the first-row transition metal manganese. However, the structural features of ionic metal transporting ABC permeases have remained undefined. Here, we present the crystal structure of the manganese transporter PsaBC from Streptococcus pneumoniae in an open-inward conformation. The type II transporter has a tightly closed transmembrane channel due to "extracellular gating" residues that prevent water permeation or ion reflux. Below these residues, the channel contains a hitherto unreported metal coordination site, which is essential for manganese translocation. Mutagenesis of the extracellular gate perturbs manganese uptake, while coordination site mutagenesis abolishes import. These structural features are highly conserved in metal-specific ABC transporters and are represented throughout the kingdoms of life. Collectively, our results define the structure of PsaBC and reveal the features required for divalent cation transport.

Probing the transport of Ni(II) ions through the internal tunnels of the Helicobacter pylori UreDFG multimeric protein complex

The survival of several pathogenic bacteria, such as Helicobacter pylori (Hp), relies on the activity of the nickel-dependent enzyme urease. Nickel insertion into urease is mediated by a multimeric chaperone complex (HpUreDFG) that is responsible for the transport of Ni(II) from a conserved metal binding motif located in the UreG dimer (CPH motif) to the catalytic site of the enzyme. The X-ray structure of HpUreDFG revealed the presence of water-filled tunnels that were proposed as a route for Ni(II) translocation. Here, we probe the transport of Ni(II) through the internal tunnels of HpUreDFG, from the CPH motif to the external surface of the complex, using microsecond-long enhanced molecular dynamics simulations. The results suggest a "bucket-brigade" mechanism whereby Ni(II) can be transported through a series of stations found along these internal pathways.

Calcium Release Channels: Structure and Function of IP3 Receptors and Ryanodine Receptors

Ca²⁺-release channels are giant membrane proteins that control the release of Ca²⁺ from the endoplasmic and sarcoplasmic reticulum. The two members, ryanodine receptors (RyRs) and inositol-1,4,5-trisphosphate Receptors (IP₃Rs), are evolutionarily related and are both activated by cytosolic Ca²⁺. They share a common architecture, but RyRs have evolved additional modules in the cytosolic region. Their massive size allows for the regulation by tens of proteins and small molecules, which can affect the opening and closing of the channels. In addition to Ca²⁺, other major triggers include IP₃ for the IP₃Rs, and depolarization of the plasma membrane for a particular RyR subtype. Their size has made them popular targets for study via electron microscopic methods, with current structures culminating near 3Å. The available structures have provided many new mechanistic insights int the binding of auxiliary proteins and small molecules, how these can regulate channel opening, and the mechanisms of disease-associated mutations. They also help scrutinize previously proposed binding sites, as some of these are now incompatible with the structures. Many questions remain around the structural effects of post-translational modifications, additional binding partners, and the higher-order complexes these channels can make in situ. This review summarizes our current knowledge about the structures of Ca²⁺-release channels and how this informs on their function.

Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations

Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.

Atomistic Details of Charge/Space Competition in the Ca2+ Selectivity of Ryanodine Receptors

Ryanodine receptors (RyRs) are ion channels responsible for the fast release of Ca²⁺ from the sarco/endoplasmic reticulum to the cytosol and show a selectivity of Ca²⁺ over monovalent cations. By utilizing a recently developed multisite Ca²⁺ model in molecular dynamic simulations, we show that multiple cations accumulate in the upper selectivity filter of RyRs, and the small size and high valence of Ca²⁺ make it preferable to K⁺ in competition for space in this confined region of negative electrostatic potential. The presence of Ca²⁺ in the upper selectivity filter significantly increases the energy barrier of K⁺ permeation, while the presence of K⁺ has little impact on the Ca²⁺ permeation. Our results provide the atomistic details of the charge/space competition mechanism for the ion selectivity of RyRs, which ensures the robustness of their Ca²⁺ release function. The mechanism could be utilized in protein- and nanoengineering for valence selectivity of ion species.