An iris-like mechanism of pore dilation in the CorA magnesium transport system

Structural and biochemical analysis of the unique interactions of the Campylobacter jejuni CorA channel protein with divalent cations

CorA is a Mg2+ channel that plays a key role in the homeostasis of intracellular Mg2+ in bacteria and archaea. CorA consists of a cytoplasmic domain and a transmembrane domain and generates a Mg2+ pathway by forming a pentamer in the cell membrane. CorA gating is regulated via negative feedback by Mg2+, which is accommodated by the pentamerization interface of the CorA cytoplasmic domain (CorACD). The Mg2+-binding sites of CorACD differ depending on the species, suggesting that the Mg2+-binding modes and Mg2+-mediated gating mechanisms of CorA vary across prokaryotes. To define the Mg2+-binding mechanism of CorA in the Campylobacter jejuni pathogen, we structurally and biochemically characterized C. jejuni CorACD (cjCorACD). cjCorACD adopts a three-layered α/β/α structure as observed in other CorA orthologs. Interestingly, cjCorACD exhibited enhanced thermostability in the presence of Ca2+, Ni2+, Zn2+, or Mn2+ in addition to Mg2+, indicating that cjCorACD interacts with diverse divalent cations. This cjCorACD stabilization is mediated by divalent cation accommodation by negatively charged residues located at the bottom of the cjCorACD structure away from the pentamerization interface. Consistently, cjCorACD exists as a monomer irrespective of the presence of divalent cations. We concluded that cjCorACD binds divalent cations in a unique pentamerization-independent manner.

MRS2 missense variation at Asp216 abrogates inhibitory Mg2+ binding, potentiating cell migration and apoptosis resistance

Mitochondrial magnesium (Mg2+) is a crucial modulator of protein stability, enzymatic activity, ATP synthesis, and cell death. Mitochondrial RNA splicing protein 2 (MRS2) is the main Mg2+ channel in the inner mitochondrial membrane that mediates influx into the matrix. Recent cryo-electron microscopy (cryo-EM) human MRS2 structures exhibit minimal conformational changes at high and low Mg2+, yet the regulation of human MRS2 and orthologues by Mg2+ binding to analogous matrix domains has been well established. Further, a missense variation at D216 has been identified associated with malignant melanoma and MRS2 expression and activity is implicated in gastric cancer. Thus, to gain more mechanistic and functional insight into Mg2+ sensing by the human MRS2 matrix domain and the association with proliferative disease, we assessed the structural, biophysical, and functional effects of a D216Q mutant. We show that the D216Q mutation is sufficient to abrogate Mg2+-binding and associated conformational changes including increased α-helicity, stability, and monomerization. Further, we reveal that the MRS2 matrix domains interact with ~μM affinity, which is weakened by up to two orders of magnitude in the presence of Mg2+ for wild-type but unaffected for D216Q. Finally, we demonstrate the importance of Mg2+ sensing by MRS2 to prevent matrix Mg2+ overload as HeLa cells overexpressing MRS2 show enhanced Mg2+ uptake, cell migration, and resistance to apoptosis while MRS2 D216Q robustly potentiates these cancer phenotypes. Collectively, our findings further define the MRS2 matrix domain as a critical Mg2+ sensor that undergoes conformational and assembly changes upon Mg2+ interactions dependent on D216 to temper matrix Mg2+ overload.

Cryo-EM structures of human magnesium channel MRS2 reveal gating and regulatory mechanisms

Magnesium ions (Mg2+) play an essential role in cellular physiology. In mitochondria, protein and ATP synthesis and various metabolic pathways are directly regulated by Mg2+. MRS2, a magnesium channel located in the inner mitochondrial membrane, mediates the influx of Mg2+ into the mitochondrial matrix and regulates Mg2+ homeostasis. Knockdown of MRS2 in human cells leads to reduced uptake of Mg2+ into mitochondria and disruption of the mitochondrial metabolism. Despite the importance of MRS2, the Mg2+ translocation and regulation mechanisms of MRS2 are still unclear. Here, using cryo-EM we report the structures of human MRS2 in the presence and absence of Mg2+ at 2.8 Å and 3.3 Å, respectively. From the homo-pentameric structures, we identify R332 and M336 as major gating residues, which are then tested using mutagenesis and two cellular divalent ion uptake assays. A network of hydrogen bonds is found connecting the gating residue R332 to the soluble domain, potentially regulating the gate. Two Mg2+-binding sites are identified in the MRS2 soluble domain, distinct from the two sites previously reported in CorA, a homolog of MRS2 in prokaryotes. Altogether, this study provides the molecular basis for understanding the Mg2+ translocation and regulatory mechanisms of MRS2.

Cryo-EM structures of human magnesium channel MRS2 reveal gating and regulatory mechanisms

Magnesium ions (Mg2+) play an essential role in cellular physiology. In mitochondria, protein and ATP synthesis and various metabolic pathways are directly regulated by Mg2+. MRS2, a magnesium channel located in the inner mitochondrial membrane, mediates the influx of Mg2+ into the mitochondrial matrix and regulates Mg2+ homeostasis. Knockdown of MRS2 in human cells leads to reduced uptake of Mg2+ into mitochondria and disruption of the mitochondrial metabolism. Despite the importance of MRS2, the Mg2+ translocation and regulation mechanisms of MRS2 are still unclear. Here, using cryo-EM we determined the structure of human MRS2 in the presence and absence of Mg2+ at 2.8 Å and 3.3 Å, respectively. From the homo-pentameric structures, we identified R332 and M336 as major gating residues, which were then tested using mutagenesis and two cellular divalent ion uptake assays. A network of hydrogen bonds was found connecting the gating residue R332 to the soluble domain, potentially regulating the gate. Two Mg2+-binding sites were identified in the MRS2 soluble domain, distinct from the two sites previously reported in CorA, a homolog of MRS2 in prokaryotes. Altogether, this study provides the molecular basis for understanding the Mg2+ translocation and regulatory mechanisms of MRS2.

Molecular basis of Mg2+ permeation through the human mitochondrial Mrs2 channel

Mitochondrial RNA splicing 2 (Mrs2), a eukaryotic CorA ortholog, enables Mg2+ to permeate the inner mitochondrial membrane and plays an important role in mitochondrial metabolic function. However, the mechanism by which Mrs2 permeates Mg2+ remains unclear. Here, we report four cryo-electron microscopy (cryo-EM) reconstructions of Homo sapiens Mrs2 (hMrs2) under various conditions. All of these hMrs2 structures form symmetrical pentamers with very similar pentamer and protomer conformations. A special structural feature of Cl--bound R-ring, which consists of five Arg332 residues, was found in the hMrs2 structure. Molecular dynamics simulations and mitochondrial Mg2+ uptake assays show that the R-ring may function as a charge repulsion barrier, and Cl- may function as a ferry to jointly gate Mg2+ permeation in hMrs2. In addition, the membrane potential is likely to be the driving force for Mg2+ permeation. Our results provide insights into the channel assembly and Mg2+ permeation of hMrs2.

Hyphal Growth in Trichosporon asahii Is Accelerated by the Addition of Magnesium

Fungal dimorphism involves two morphologies: a unicellular yeast cell and a multicellular hyphal form. Invasion of hyphae into human cells causes severe opportunistic infections. The transition between yeast and hyphal forms is associated with the virulence of fungi; however, the mechanism is poorly understood. Therefore, we aimed to identify factors that induce hyphal growth of Trichosporon asahii, a dimorphic basidiomycete that causes trichosporonosis. T. asahii showed poor growth and formed small cells containing large lipid droplets and fragmented mitochondria when cultivated for 16 h in a nutrient-deficient liquid medium. However, these phenotypes were suppressed via the addition of yeast nitrogen base. When T. asahii cells were cultivated in the presence of different compounds present in the yeast nitrogen base, we found that magnesium sulfate was a key factor for inducing cell elongation, and its addition dramatically restored hyphal growth in T. asahii. In T. asahii hyphae, vacuoles were enlarged, the size of lipid droplets was decreased, and mitochondria were distributed throughout the cell cytoplasm and adjacent to the cell walls. Additionally, hyphal growth was disrupted due to treatment with an actin inhibitor. The actin inhibitor latrunculin A disrupted the mitochondrial distribution even in hyphal cells. Furthermore, magnesium sulfate treatment accelerated hyphal growth in T. asahii for 72 h when the cells were cultivated in a nutrient-deficient liquid medium. Collectively, our results suggest that an increase in magnesium levels triggers the transition from the yeast to hyphal form in T. asahii. These findings will support studies on the pathogenesis of fungi and aid in developing treatments. IMPORTANCE Understanding the mechanism underlying fungal dimorphism is crucial to discern its invasion into human cells. Invasion is caused by the hyphal form rather than the yeast form; therefore, it is important to understand the mechanism of transition from the yeast to hyphal form. To study the transition mechanism, we utilized Trichosporon asahii, a dimorphic basidiomycete that causes severe trichosporonosis since there are fewer studies on T. asahii than on ascomycetes. This study suggests that an increase in Mg2+, the most abundant mineral in living cells, triggers growth of filamentous hyphae and increases the distribution of mitochondria throughout the cell cytoplasm and adjacent to the cell walls in T. asahii. Understanding the mechanism of hyphal growth triggered by Mg2+ increase will provide a model system to explore fungal pathogenicity in the future.

The human MRS2 magnesium-binding domain is a regulatory feedback switch for channel activity

Mitochondrial RNA splicing 2 (MRS2) forms a magnesium (Mg2+) entry protein channel in mitochondria. Whereas MRS2 contains two transmembrane domains constituting a pore on the inner mitochondrial membrane, most of the protein resides within the matrix. Yet, the precise structural and functional role of this obtrusive amino terminal domain (NTD) in human MRS2 is unknown. Here, we show that the MRS2 NTD self-associates into a homodimer, contrasting the pentameric assembly of CorA, an orthologous bacterial channel. Mg2+ and calcium suppress lower and higher order oligomerization of MRS2 NTD, whereas cobalt has no effect on the NTD but disassembles full-length MRS2. Mutating-pinpointed residues-mediating Mg2+ binding to the NTD not only selectively decreases Mg2+-binding affinity ∼sevenfold but also abrogates Mg2+ binding-induced secondary, tertiary, and quaternary structure changes. Disruption of NTD Mg2+ binding strikingly potentiates mitochondrial Mg2+ uptake in WT and Mrs2 knockout cells. Our work exposes a mechanism for human MRS2 autoregulation by negative feedback from the NTD and identifies a novel gain of function mutant with broad applicability to future Mg2+ signaling research.

Slow conformational dynamics of the human A2A adenosine receptor are temporally ordered

A more complete depiction of protein energy landscapes includes the identification of different function-related conformational states and the determination of the pathways connecting them. We used total internal reflection fluorescence (TIRF) imaging to investigate the conformational dynamics of the human A2A adenosine receptor (A2AAR), a class A G protein-coupled receptor (GPCR), at the single-molecule level. Slow, reversible conformational exchange was observed among three different fluorescence emission states populated for agonist-bound A2AAR. Transitions among these states predominantly occurred in a specific order, and exchange between inactive and active-like conformations proceeded through an intermediate state. Models derived from molecular dynamics simulations with available A2AAR structures rationalized the relative fluorescence emission intensities for the highest and lowest emission states but not the transition state. This suggests that the functionally critical intermediate state required to achieve activation is not currently visualized among available A2AAR structures.

Mg2+-dependent conformational equilibria in CorA and an integrated view on transport regulation

The CorA family of proteins regulates the homeostasis of divalent metal ions in many bacteria, archaea, and eukaryotic mitochondria, making it an important target in the investigation of the mechanisms of transport and its functional regulation. Although numerous structures of open and closed channels are now available for the CorA family, the mechanism of the transport regulation remains elusive. Here, we investigated the conformational distribution and associated dynamic behaviour of the pentameric Mg²⁺ channel CorA at room temperature using small-angle neutron scattering (SANS) in combination with molecular dynamics (MD) simulations and solid-state nuclear magnetic resonance spectroscopy (NMR). We find that neither the Mg²⁺-bound closed structure nor the Mg²⁺-free open forms are sufficient to explain the average conformation of CorA. Our data support the presence of conformational equilibria between multiple states, and we further find a variation in the behaviour of the backbone dynamics with and without Mg²⁺. We propose that CorA must be in a dynamic equilibrium between different non-conducting states, both symmetric and asymmetric, regardless of bound Mg²⁺ but that conducting states become more populated in Mg²⁺-free conditions. These properties are regulated by backbone dynamics and are key to understanding the functional regulation of CorA.

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.

Asymmetric CorA Gating Mechanism as Observed by Molecular Dynamics Simulations

The CorA family of proteins plays a housekeeping role in the homeostasis of divalent metal ions in many bacteria and archaea as well as in mitochondria of eukaryotes, rendering it an important target to study the mechanisms of divalent transport and regulation across different life domains. Despite numerous studies, the mechanistic details of the channel gating and the transport of the metal ions are still not entirely understood. Here, we use all-atom and coarse-grained molecular dynamics simulations combined with in vitro experiments to investigate the influence of divalent cations on the function of CorA. Simulations reveal pronounced asymmetric movements of monomers that enable the rotation of the α7 helix and the cytoplasmic subdomain with the subsequent formation of new interactions and the opening of the channel. These computational results are functionally validated using site-directed mutagenesis of the intracellular cytoplasmic domain residues and biochemical assays. The obtained results infer a complex network of interactions altering the structure of CorA to allow gating. Furthermore, we attempt to reconcile the existing gating hypotheses for CorA to conclude the mechanism of transport of divalent cations via these proteins.

An Approach to comparing protein structures and origami models - Part 2. Multi-domain proteins

Protein structure is an important field of research, with particular significance in its potential applications in biomedicine and nanotechnology. In a recent study, we presented a general approach for comparing protein structures and origami models and demonstrated it with single-domain proteins. For example, the analysis of the α-helical barrel of the outer membrane protein A (OmpA) suggests that there are similar patterns between its structure and the Kresling origami model, providing insight into structure-activity relationships. Here we demonstrate that our approach can be expanded beyond single-domain proteins to also include multi-domain proteins, and to study dynamic processes of biomolecules. Two examples are given: (1) The eukaryotic chaperonin (TRiC) protein is compared with a newly generated origami model, and with an origami model that is constructed from two copies of the Flasher origami model, and (2) the CorA Magnesium transport system is compared with a newly generated origami model and with an origami model that combines the Kresling and Flasher origami models. Based on the analysis of the analog origami models, it is indicated that it is possible to identify building blocks for constructing assembled origami models that are analogous to protein structures. In addition, it is identified that the expansion/collapse mechanisms of the TRiC and CorA are auxetic. Namely, these proteins require a single motion for synchronized folding along two or three axes.

Biophysical Insight on the Membrane Insertion of an Arginine-Rich Cell-Penetrating Peptide

Cell-penetrating peptides (CPPs) are short peptides that can translocate and transport cargoes into the intracellular milieu by crossing biological membranes. The mode of interaction and internalization of cell-penetrating peptides has long been controversial. While their interaction with anionic membranes is quite well understood, the insertion and behavior of CPPs in zwitterionic membranes, a major lipid component of eukaryotic cell membranes, is poorly studied. Herein, we investigated the membrane insertion of RW16 into zwitterionic membranes, a versatile CPP that also presents antibacterial and antitumor activities. Using complementary approaches, including NMR spectroscopy, fluorescence spectroscopy, circular dichroism, and molecular dynamic simulations, we determined the high-resolution structure of RW16 and measured its membrane insertion and orientation properties into zwitterionic membranes. Altogether, these results contribute to explaining the versatile properties of this peptide toward zwitterionic lipids.

Automated cryo-EM structure refinement using correlation-driven molecular dynamics

We present a correlation-driven molecular dynamics (CDMD) method for automated refinement of atomistic models into cryo-electron microscopy (cryo-EM) maps at resolutions ranging from near-atomic to subnanometer. It utilizes a chemically accurate force field and thermodynamic sampling to improve the real-space correlation between the modeled structure and the cryo-EM map. Our framework employs a gradual increase in resolution and map-model agreement as well as simulated annealing, and allows fully automated refinement without manual intervention or any additional rotamer- and backbone-specific restraints. Using multiple challenging systems covering a wide range of map resolutions, system sizes, starting model geometries and distances from the target state, we assess the quality of generated models in terms of both model accuracy and potential of overfitting. To provide an objective comparison, we apply several well-established methods across all examples and demonstrate that CDMD performs best in most cases.

The hydrophobic nature of a novel membrane interface regulates the enzyme activity of a voltage-sensing phosphatase

Voltage-sensing phosphatases (VSP) contain a voltage sensor domain (VSD) similar to that of voltage-gated ion channels but lack a pore-gate domain. A VSD in a VSP regulates the cytoplasmic catalytic region (CCR). However, the mechanisms by which the VSD couples to the CCR remain elusive. Here we report a membrane interface (named ‘the hydrophobic spine’), which is essential for the coupling of the VSD and CCR. Our molecular dynamics simulations suggest that the hydrophobic spine of Ciona intestinalis VSP (Ci-VSP) provides a hinge-like motion for the CCR through the loose membrane association of the phosphatase domain. Electrophysiological experiments indicate that the voltage-dependent phosphatase activity of Ci-VSP depends on the hydrophobicity and presence of an aromatic ring in the hydrophobic spine. Analysis of conformational changes in the VSD and CCR suggests that the VSP has two states with distinct enzyme activities and that the second transition depends on the hydrophobic spine.

Loss of cytosolic Mg2+ binding sites in the Thermotoga maritima CorA Mg2+ channel is not sufficient for channel opening

The CorA Mg²⁺ channel is a homopentamer with five-fold symmetry. Each monomer consists of a large cytoplasmic domain and two transmembrane helices connected via a short periplasmic loop. In the Thermotoga maritima CorA crystal structure, a Mg²⁺ is bound between D89 of one monomer and D253 of the adjacent monomer (M1 binding site). Release of Mg²⁺ from these sites has been hypothesized to cause opening of the channel. We generated mutants to disrupt Mg²⁺ interaction with the M1 site. Crystal structures of the D89K/D253K and D89R/D253R mutants, determined to 3.05 and 3.3 Å, respectively, showed no significant structural differences with the wild type structure despite absence of Mg²⁺ at the M1 sites. Both mutants still appear to be in the closed state. All three mutant CorA proteins exhibited transport of ⁶³Ni²⁺, indicating functionality. Thus, absence of Mg²⁺ from the M1 sites neither causes channel opening nor prevents function. We also provide evidence that the T. maritima CorA is a Mg²⁺ channel and not a Co²⁺ channel.

Mapping the Interface of a GPCR Dimer: A Structural Model of the A2A Adenosine and D2 Dopamine Receptor Heteromer

The A_2A adenosine (A_2AR) and D₂ dopamine (D₂R) receptors form oligomers in the cell membrane and allosteric interactions across the A_2AR–D₂R heteromer represent a target for development of drugs against central nervous system disorders. However, understanding of the molecular determinants of A_2AR–D₂R heteromerization and the allosteric antagonistic interactions between the receptor protomers is still limited. In this work, a structural model of the A_2AR–D₂R heterodimer was generated using a combined experimental and computational approach. Regions involved in the heteromer interface were modeled based on the effects of peptides derived from the transmembrane (TM) helices on A_2AR–D₂R receptor–receptor interactions in bioluminescence resonance energy transfer (BRET) and proximity ligation assays. Peptides corresponding to TM-IV and TM-V of the A_2AR blocked heterodimer interactions and disrupted the allosteric effect of A_2AR activation on D₂R agonist binding. Protein–protein docking was used to construct a model of the A_2AR–D₂R heterodimer with a TM-IV/V interface, which was refined using molecular dynamics simulations. Mutations in the predicted interface reduced A_2AR–D₂R interactions in BRET experiments and altered the allosteric modulation. The heterodimer model provided insights into the structural basis of allosteric modulation and the technique developed to characterize the A_2AR–D₂R interface can be extended to study the many other G protein-coupled receptors that engage in heteroreceptor complexes.

Agonist-induced dimer dissociation as a macromolecular step in G protein-coupled receptor signaling

G protein-coupled receptors (GPCRs) constitute the largest family of cell surface receptors. They can exist and act as dimers, but the requirement of dimers for agonist-induced signal initiation and structural dynamics remains largely unknown. Frizzled 6 (FZD₆) is a member of Class F GPCRs, which bind WNT proteins to initiate signaling. Here, we show that FZD₆ dimerizes and that the dimer interface of FZD₆ is formed by the transmembrane α-helices four and five. Most importantly, we present the agonist-induced dissociation/re-association of a GPCR dimer through the use of live cell imaging techniques. Further analysis of a dimerization-impaired FZD₆ mutant indicates that dimer dissociation is an integral part of FZD₆ signaling to extracellular signal-regulated kinases1/2. The discovery of agonist-dependent dynamics of dimers as an intrinsic process of receptor activation extends our understanding of Class F and other dimerizing GPCRs, offering novel targets for dimer-interfering small molecules.

Frizzled 6 (FZD₆) is a G protein-coupled receptor (GPCR) involved in several cellular processes. Here, the authors use live cell imaging and spectroscopy to show that FZD₆ forms dimers, whose association is regulated by WNT proteins and that dimer dissociation is crucial for FZD₆ signaling.

Asymmetry in structural response of inner and outer transmembrane segments of CorA protein by a coarse-grain model

Structure of CorA protein and its inner (i.corA) and outer (o.corA) transmembrane (TM) components are investigated as a function of temperature by a coarse-grained Monte Carlo simulation. Thermal response of i.corA is found to differ considerably from that of the outer component, o.corA. Analysis of the radius of gyration reveals that the inner TM component undergoes a continuous transition from a globular conformation to a random coil structure on raising the temperature. In contrast, the outer transmembrane component exhibits an abrupt (nearly discontinuous) thermal response in a narrow range of temperature. Scaling of the structure factor shows a globular structure of i.corA at a low temperature with an effective dimension D ∼ 3 and a random coil at a high temperature with D ∼ 2. The residue distribution in o.corA is slightly sparser than that of i.corA in a narrow thermos-responsive regime. The difference in thermos-response characteristics of these components (i.corA and o.corA) may reflect their unique transmembrane functions.

Cryo-EM Structures of the Magnesium Channel CorA Reveal Symmetry Break upon Gating

CorA, the major Mg(2+) uptake system in prokaryotes, is gated by intracellular Mg(2+) (KD ∼ 1-2 mM). X-ray crystallographic studies of CorA show similar conformations under Mg(2+)-bound and Mg(2+)-free conditions, but EPR spectroscopic studies reveal large Mg(2+)-driven quaternary conformational changes. Here, we determined cryo-EM structures of CorA in the Mg(2+)-bound closed conformation and in two open Mg(2+)-free states at resolutions of 3.8, 7.1, and 7.1 Å, respectively. In the absence of bound Mg(2+), four of the five subunits are displaced to variable extents (∼ 10-25 Å) by hinge-like motions as large as ∼ 35° at the stalk helix. The transition between a single 5-fold symmetric closed state and an ensemble of low Mg(2+), open, asymmetric conformational states is, thus, the key structural signature of CorA gating. This mechanism is likely to apply to other structurally similar divalent ion channels.

Interactions of Polyethylenimines with Zwitterionic and Anionic Lipid Membranes

Interactions between polyethylenimines (PEIs) and phospholipid membranes are of fundamental importance for various biophysical applications of these polymers such as gene delivery. Despite investigations into the nature of these interactions, their molecular basis remains poorly understood. In this article, we combined experimental methods and atomistic molecular dynamics (MD) simulations to obtain comprehensive insight into the effect of linear and branched PEIs on zwitterionic and anionic bilayers used as simple models of mammalian cellular membranes. Our results show that PEIs adsorb only partially on the surface of zwitterionic membranes by forming hydrogen bonds to the lipid headgroups, whereas a large part of the polymer chains dangles freely in the aqueous phase. In contrast, PEIs readily adhere to and insert into the anionic membrane. The attraction of the polymer chains to the membrane is due to electrostatic interactions as well as hydrogen bonding between the amine groups of PEI and the phosphate groups of lipids. These interactions were found to induce a substantial reorganization of the bilayer in the polymer vicinity due to the reorientation of lipid molecules. The lipid headgroups were pulled toward the center of the membrane, which can facilitate transmembrane translocations of anionic lipids. Furthermore, the PEI-lipid interactions affect the stability of liposomal dispersions, but we did not see any evidence of disruption of the vesicular structures into small fragments at polymer concentrations typically used in gene therapy. Our results provide a detailed molecular-level description of the lipid organization in the membrane in the presence of polycations that can be useful in understanding their mechanisms of in vitro and in vivo cytotoxicity.

Inner and Outer Coordination Shells of Mg(2+) in CorA Selectivity Filter from Molecular Dynamics Simulations

Structural data of CorA Mg(2+) channels show that the five Gly-Met-Asn (GMN) motifs at the periplasmic loop of the pentamer structure form a molecular scaffold serving as a selectivity filter. Unfortunately, knowledge about the cation selectivity of Mg(2+) channels remains limited. Since Mg(2+) in aqueous solution has a strong first hydration shell and apparent second hydration sphere, the coordination structure of Mg(2+) in a CorA selectivity filter is expected to be different from that in bulk water. Hence, this study investigated the hydration structure and ligand coordination of Mg(2+) in a selectivity filter of CorA using molecular dynamics (MD) simulations. The simulations reveal that the inner-shell structure of Mg(2+) in the filter is not significantly different from that in aqueous solution. The major difference is the characteristic structural features of the outer shell. The GMN residues engage indirectly in the interactions with the metal ion as ligands in the second shell of Mg(2+). Loss of hydrogen bonds between inner- and outer-shell waters observed from Mg(2+) in bulk water is mostly compensated by interactions between waters in the first solvation shell and the GMN motif. Some water molecules in the second shell remain in the selectivity filter and become less mobile to support the metal binding. Removal of Mg(2+) from the divalent cation sensor sites of the protein had an impact on the structure and metal binding of the filter. From the results, it can be concluded that the GMN motif enhances the affinity of the metal binding site in the CorA selectivity filter by acting as an outer coordination ligand.

Potential of mean force analysis of the self-association of leucine-rich transmembrane α-helices: difference between atomistic and coarse-grained simulations

Interaction of transmembrane (TM) proteins is important in many biological processes. Large-scale computational studies using coarse-grained (CG) simulations are becoming popular. However, most CG model parameters have not fully been calibrated with respect to lateral interactions of TM peptide segments. Here, we compare the potential of mean forces (PMFs) of dimerization of TM helices obtained using a MARTINI CG model and an atomistic (AT) Berger lipids-OPLS/AA model (AT(OPLS)). For helical, tryptophan-flanked, leucine-rich peptides (WL15 and WALP15) embedded in a parallel configuration in an octane slab, the AT(OPLS) PMF profiles showed a shallow minimum (with a depth of approximately 3 kJ/mol; i.e., a weak tendency to dimerize). A similar analysis using the CHARMM36 all-atom model (AT(CHARMM)) showed comparable results. In contrast, the CG analysis generally showed steep PMF curves with depths of approximately 16-22 kJ/mol, suggesting a stronger tendency to dimerize compared to the AT model. This CG > AT discrepancy in the propensity for dimerization was also seen for dilauroylphosphatidylcholine (DLPC)-embedded peptides. For a WL15 (and WALP15)/DLPC bilayer system, AT(OPLS) PMF showed a repulsive mean force for a wide range of interhelical distances, in contrast to the attractive forces observed in the octane system. The change from the octane slab to the DLPC bilayer also mitigated the dimerization propensity in the CG system. The dimerization energies of CG (AALALAA)3 peptides in DLPC and dioleoylphosphatidylcholine bilayers were in good agreement with previous experimental data. The lipid headgroup, but not the length of the lipid tails, was a key causative factor contributing to the differences between octane and DLPC. Furthermore, the CG model, but not the AT model, showed high sensitivity to changes in amino acid residues located near the lipid-water interface and hydrophobic mismatch between the peptides and membrane. These findings may help interpret CG and AT simulation results on membrane proteins.

Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA

Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg²⁺ into cells, we examine early events in the relaxation of Mg²⁺ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or “wetting” events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg²⁺ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction.

This study shows how rapid wetting/dewetting transitions in the pores of ion channels participate in the control of biological ion permeation. Ion channels catalyze ionic permeation across non-polar membranes via water-filled pores. However, non-polar stretches or hydrophobic bottlenecks are present in the pores of many ion channels. To clarify the relationship between channel regulation, pore hydration, and ion permeation, we examine how the slow relaxation of magnesium channel CorA from its closed state towards its open state modulates wetting of its hydrophobic bottleneck. Results provide a quantitative description of wetting and dewetting probabilities and kinetics and a quantitative relationship between the extent of pore hydration and the energetics of ion permeation, consistent with a mechanism of hydrophobic gating.

Identification of the antiepileptic racetam binding site in the synaptic vesicle protein 2A by molecular dynamics and docking simulations

Synaptic vesicle protein 2A (SV2A) is an integral membrane protein necessary for the proper function of the central nervous system and is associated to the physiopathology of epilepsy. SV2A is the molecular target of the anti-epileptic drug levetiracetam and its racetam analogs. The racetam binding site in SV2A and the non-covalent interactions between racetams and SV2A are currently unknown; therefore, an in silico study was performed to explore these issues. Since SV2A has not been structurally characterized with X-ray crystallography or nuclear magnetic resonance, a three-dimensional (3D) model was built. The model was refined by performing a molecular dynamics simulation (MDS) and the interactions of SV2A with the racetams were determined by docking studies. A reliable 3D model of SV2A was obtained; it reached structural equilibrium during the last 15 ns of the MDS (50 ns) with remaining structural motions in the N-terminus and long cytoplasmic loop. The docking studies revealed that hydrophobic interactions and hydrogen bonds participate importantly in ligand recognition within the binding site. Residues T456, S665, W666, D670 and L689 were important for racetam binding within the trans-membrane hydrophilic core of SV2A. Identifying the racetam binding site within SV2A should facilitate the synthesis of suitable radio-ligands to study treatment response and possibly epilepsy progression.

A cleanroom sleeping environment's impact on markers of oxidative stress, immune dysregulation, and behavior in children with autism spectrum disorders

BACKGROUND

An emerging paradigm suggests children with autism display a unique pattern of environmental, genetic, and epigenetic triggers that make them susceptible to developing dysfunctional heavy metal and chemical detoxification systems. These abnormalities could be caused by alterations in the methylation, sulfation, and metalloprotein pathways. This study sought to evaluate the physiological and behavioral effects of children with autism sleeping in an International Organization for Standardization Class 5 cleanroom.

METHODS

Ten children with autism, ages 3-12, slept in a cleanroom for two weeks to evaluate changes in toxin levels, oxidative stress, immune dysregulation, and behavior. Before and after the children slept in the cleanroom, samples of blood and hair and rating scale scores were obtained to assess these changes.

RESULTS

Five children significantly lowered their concentration of oxidized glutathione, a biomarker of oxidative stress. The younger cohort, age 5 and under, showed significantly greater mean decreases in two markers of immune dysregulation, CD3% and CD4%, than the older cohort. Changes in serum magnesium, influencing neuronal regulation, correlated negatively while changes in serum iron, affecting oxygenation of tissues, correlated positively with age. Changes in serum benzene and PCB 28 concentrations showed significant negative correlations with age. The younger children demonstrated significant improvements on behavioral rating scales compared to the older children. In a younger pair of identical twins, one twin showed significantly greater improvements in 4 out of 5 markers of oxidative stress, which corresponded with better overall behavioral rating scale scores than the other twin.

CONCLUSIONS

Younger children who slept in the cleanroom altered elemental levels, decreased immune dysregulation, and improved behavioral rating scales, suggesting that their detoxification metabolism was briefly enhanced. The older children displayed a worsening in behavioral rating scale performance, which may have been caused by the mobilization of toxins from their tissues. The interpretation of this exploratory study is limited by lack of a control group and small sample size. The changes in physiology and behavior noted suggest that performance of larger, prospective controlled studies of exposure to nighttime or 24 hour cleanroom conditions for longer time periods may be useful for understanding detoxification in children with autism.

TRIAL REGISTRATION

Clinical Trial Registration Number NCT02195401 (Obtained July 18, 2014).

Sodium ion binding pocket mutations and adenosine A2A receptor function

Recently we identified a sodium ion binding pocket in a high-resolution structure of the human adenosine A2A receptor. In the present study we explored this binding site through site-directed mutagenesis and molecular dynamics simulations. Amino acids in the pocket were mutated to alanine, and their influence on agonist and antagonist affinity, allosterism by sodium ions and amilorides, and receptor functionality was explored. Mutation of the polar residues in the Na(+) pocket were shown to either abrogate (D52A(2.50) and N284A(7.49)) or reduce (S91A(3.39), W246A(6.48), and N280A(7.45)) the negative allosteric effect of sodium ions on agonist binding. Mutations D52A(2.50) and N284A(7.49) completely abolished receptor signaling, whereas mutations S91A(3.39) and N280A(7.45) elevated basal activity and mutations S91A(3.39), W246A(6.48), and N280A(7.45) decreased agonist-stimulated receptor signaling. In molecular dynamics simulations D52A(2.50) directly affected the mobility of sodium ions, which readily migrated to another pocket formed by Glu13(1.39) and His278(7.43). The D52A(2.50) mutation also decreased the potency of amiloride with respect to ligand displacement but did not change orthosteric ligand affinity. In contrast, W246A(6.48) increased some of the allosteric effects of sodium ions and amiloride, whereas orthosteric ligand binding was decreased. These new findings suggest that the sodium ion in the allosteric binding pocket not only impacts ligand affinity but also plays a vital role in receptor signaling. Because the sodium ion binding pocket is highly conserved in other class A G protein-coupled receptors, our findings may have a general relevance for these receptors and may guide the design of novel synthetic allosteric modulators or bitopic ligands.

Solvent cavitation under solvophobic confinement

The stability of liquids under solvophobic confinement can tip in favor of the vapor phase, nucleating a liquid-to-vapor phase transition that induces attractive forces between confining surfaces. In the case of water adjacent to hydrophobic surfaces, experimental and theoretical evidence support confinement-mediated evaporation stabilization of biomolecular and colloidal assemblies. The macroscopic thermodynamic theory of cavitation under confinement establishes the connection between the size of the confining surfaces, interfacial free energies, and bulk solvent pressure with the critical evaporation separation and interfacial forces. While molecular simulations have confirmed the broad theoretical trends, a quantitative comparison based on independent measurements of the interfacial free energies and liquid-vapor coexistence properties has, to the best of our knowledge, not yet been performed. To overcome the challenges of simulating a large number of systems to validate scaling predictions for a three-dimensional fluid, we simulate both the forces and liquid-vapor coexistence properties of a two-dimensional Lennard-Jones fluid confined between solvophobic plates over a range of plate sizes and reservoir pressures. Our simulations quantitatively agree with theoretical predictions for solvent-mediated forces and critical evaporation separations once the length dependence of the solvation free energy of an individual confining plate is taken into account. The effective solid-liquid line tension length dependence results from molecular scale correlations for solvating microscopic plates and asymptotically decays to the macroscopic value for plates longer than 150 solvent diameters. The success of the macroscopic thermodynamic theory at describing two-dimensional liquids suggests application to surfactant monolayers to experimentally confirm confinement-mediated cavitation.

General rules for the arrangements and gating motions of pore-lining helices in homomeric ion channels

The pore-lining helix (PLH) bundles are central to the function of all ion channels, as their conformational rearrangements dictate channel gating. Here we explore all plausible oligomeric arrangements of the PLH bundles within homomeric ion channels by building models using generic restraints. In particular, the distance between two neighbouring PLHs was bounded both below and above in order to avoid steric clash and allow proper packing. The resulting models provide a theoretical representation of the accessible space for oligomeric arrangements. While the represented space is confined, it encompasses nearly all the ion channel PLH bundles for which the structures are currently known. For a multitude of channels, gating models suggested by paths within the confined accessible space are in qualitative agreement with those established in previous structural and computational studies.

Rearrangements of the pore-lining helix (PLH) bundles of ion channels are central to their gating mechanisms. Here, Dai et al. use a modelling approach to define the general rules that govern the arrangements and gating motions of the PLHs in homomeric ion channels.

Molecular mechanism of Mg2+-dependent gating in CorA

CorA is the major transport system responsible for Mg²⁺ uptake in bacteria and can functionally substitute for its homologue Mrs2p in the yeast inner mitochondrial membrane. Although several CorA crystal structures are available, the molecular mechanism of Mg²⁺ uptake remains to be established. Here we use EPR spectroscopy, electrophysiology and molecular dynamic simulations to show that CorA is regulated by cytoplasmic Mg²⁺ acting as a ligand and elucidate the basic conformational rearrangements responsible for Mg²⁺-dependent gating. Mg²⁺ unbinding at the divalent cation sensor triggers a conformational change that leads to the inward motion of the stalk helix, which propagates to the pore forming transmembrane helix TM1. Helical tilting and rotation in TM1 generates an iris-like motion that increases the diameter of the permeation pathway, triggering ion conduction. This work establishes the molecular basis of a Mg²⁺-driven negative feedback loop in CorA as the key physiological event controlling Mg²⁺ uptake and homeostasis in prokaryotes.

Peregrination of the selectivity filter delineates the pore of the human voltage-gated proton channel hHV1

Extraordinary selectivity is crucial to all proton-conducting molecules, including the human voltage-gated proton channel (hH_V1), because the proton concentration is >10⁶ times lower than that of other cations. Here we use “selectivity filter scanning” to elucidate the molecular requirements for proton-specific conduction in hH_V1. Asp¹¹², in the middle of the S1 transmembrane helix, is an essential part of the selectivity filter in wild-type (WT) channels. After neutralizing Asp¹¹² by mutating it to Ala (D112A), we introduced Asp at each position along S1 from 108 to 118, searching for “second site suppressor” activity. Surprisingly, most mutants lacked even the anion conduction exhibited by D112A. Proton-specific conduction was restored only with Asp or Glu at position 116. The D112V/V116D channel strikingly resembled WT in selectivity, kinetics, and ΔpH-dependent gating. The S4 segment of this mutant has similar accessibility to WT in open channels, because R211H/D112V/V116D was inhibited by internally applied Zn²⁺. Asp at position 109 allowed anion permeation in combination with D112A but did not rescue function in the nonconducting D112V mutant, indicating that selectivity is established externally to the constriction at F150. The three positions that permitted conduction all line the pore in our homology model, clearly delineating the conduction pathway. Evidently, a carboxyl group must face the pore directly to enable conduction. Molecular dynamics simulations indicate reorganization of hydrogen bond networks in the external vestibule in D112V/V116D. At both positions where it produces proton selectivity, Asp frequently engages in salt linkage with one or more Arg residues from S4. Surprisingly, mean hydration profiles were similar in proton-selective, anion-permeable, and nonconducting constructs. That the selectivity filter functions in a new location helps to define local environmental features required to produce proton-selective conduction.

The role of a sodium ion binding site in the allosteric modulation of the A(2A) adenosine G protein-coupled receptor

The function of G protein-coupled receptors (GPCRs) can be modulated by a number of endogenous allosteric molecules. In this study, we used molecular dynamics, radioligand binding, and thermostability experiments to elucidate the role of the recently discovered sodium ion binding site in the allosteric modulation of the human A(2A) adenosine receptor, conserved among class A GPCRs. While the binding of antagonists and sodium ions to the receptor was noncompetitive in nature, the binding of agonists and sodium ions appears to require mutually exclusive conformational states of the receptor. Amiloride analogs can also bind to the sodium binding pocket, showing distinct patterns of agonist and antagonist modulation. These findings suggest that physiological concentrations of sodium ions affect functionally relevant conformational states of GPCRs and can help to design novel synthetic allosteric modulators or bitopic ligands exploiting the sodium ion binding pocket.

Bindings of hMRP1 transmembrane peptides with dodecylphosphocholine and dodecyl-β-d-maltoside micelles: a molecular dynamics simulation study

In this paper, we describe molecular dynamics simulation results of the interactions between four peptides (mTM10, mTM16, TM17 and KTM17) with micelles of dodecylphosphocholine (DPC) and dodecyl-β-d-maltoside (DDM). These peptides represent three transmembrane fragments (TM10, 16 and 17) from the MSD1 and MSD2 membrane-spanning domains of an ABC membrane protein (hMRP1), which play roles in the protein functions. The peptide-micelle complex structures, including the tryptophan accessibility and dynamics were compared to circular dichroism and fluorescence studies obtained in water, trifluoroethanol and with micelles. Our work provides additional results not directly accessible by experiments that give further support to the fact that these peptides adopt an interfacial conformation within the micelles. We also show that the peptides are more buried in DDM than in DPC, and consequently, that they have a larger surface exposure to water in DPC than in DDM. As noted previously by simulations and experiments we have also observed formation of cation-π bonds between the phosphocholine DPC headgroup and Trp peptide residue. Concerning the peptide secondary structures (SS), we find that in TFE their initial helical conformations are maintained during the simulation, whereas in water their initial SS are lost after few nanoseconds of simulation. An intermediate situation is observed with micelles, where the peptides remain partially folded and more structured in DDM than in DPC. Finally, our results show no sign of β-strand structure formation as invoked by far-UV CD experiments even when three identical peptides are simulated either in water or with micelles.

Construction and validation of a homology model of the human voltage-gated proton channel hHV1

The topological similarity of voltage-gated proton channels (H_V1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H_V1s have a similar structure. We present the construction and validation of a homology model of the human H_V1 (hH_V1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH_V1 by exploiting the homology of hH_V1 with VSDs of K⁺ and Na⁺ channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH_V1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp¹¹². Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH_V1. Each Arg residue in the S4 helix of hH_V1 was replaced by His to test accessibility using Zn²⁺ as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg²¹¹, the third Arg residue in the S4 helix in hH_V1, remains accessible to the internal solution and is located near the charge transfer center, Phe¹⁵⁰.

The structure and regulation of magnesium selective ion channels

The magnesium ion (Mg(2+)) is the most abundant divalent cation within cells. In man, Mg(2+)-deficiency is associated with diseases affecting the heart, muscle, bone, immune, and nervous systems. Despite its impact on human health, little is known about the molecular mechanisms that regulate magnesium transport and storage. Complete structural information on eukaryotic Mg(2+)-transport proteins is currently lacking due to associated technical challenges. The prokaryotic MgtE and CorA magnesium transport systems have recently succumbed to structure determination by X-ray crystallography, providing first views of these ubiquitous and essential Mg(2+)-channels. MgtE and CorA are unique among known membrane protein structures, each revealing a novel protein fold containing distinct arrangements of ten transmembrane-spanning α-helices. Structural and functional analyses have established that Mg(2+)-selectivity in MgtE and CorA occurs through distinct mechanisms. Conserved acidic side-chains appear to form the selectivity filter in MgtE, whereas conserved asparagines coordinate hydrated Mg(2+)-ions within the selectivity filter of CorA. Common structural themes have also emerged whereby MgtE and CorA sense and respond to physiologically relevant, intracellular Mg(2+)-levels through dedicated regulatory domains. Within these domains, multiple primary and secondary Mg(2+)-binding sites serve to staple these ion channels into their respective closed conformations, implying that Mg(2+)-transport is well guarded and very tightly regulated. The MgtE and CorA proteins represent valuable structural templates to better understand the related eukaryotic SLC41 and Mrs2-Alr1 magnesium channels. Herein, we review the structure, function and regulation of MgtE and CorA and consider these unique proteins within the expanding universe of ion channel and transporter structural biology.

Conical lipids in flat bilayers induce packing defects similar to that induced by positive curvature

In biological membranes, changes in lipid composition or mechanical deformations produce defects in the geometrical arrangement of lipids, thus allowing the adsorption of certain peripheral proteins. Here, we perform molecular dynamics simulations on bilayers containing a cylindrical lipid (PC) and a conical lipid (DOG). Profiles of atomic density and lateral pressure across the bilayer show differences in the acyl chain region due to deeper partitioning of DOG compared to PC. However, such analyses are less informative for the interfacial region where peripheral proteins adsorb. To circumvent this limitation, we develop, to our knowledge, a new method of membrane surface analysis. This method allows the identification of chemical defects, where hydrocarbon chains are accessible to the solvent, and geometrical defects, i.e., voids deeper than the glycerol backbone. The size and number of both types of defects increase with the number of monounsaturated acyl chains in PC and with the introduction of DOG, although the defects do not colocalize with the conical lipid. Interestingly, the size and probability of the defects promoted by DOG resemble those induced by positive curvature, thus explaining why conical lipids and positive curvature can both drive the adsorption of peripheral proteins that use hydrophobic residues as membrane anchors.

Amphipathic lipid packing sensor motifs: probing bilayer defects with hydrophobic residues

Sensing membrane curvature allows fine-tuning of complex reactions that occur at the surface of membrane-bound organelles. One of the most sensitive membrane curvature sensors, the Amphipathic Lipid Packing Sensor (ALPS) motif, does not seem to recognize the curved surface geometry of membranes per se; rather, it recognizes defects in lipid packing that arise from membrane bending. In a companion paper, we show that these defects can be mimicked by introducing conical lipids in a flat lipid bilayer, in agreement with experimental observations. Here, we use molecular-dynamics (MD) simulations to characterize ALPS binding to such lipid bilayers. The ALPS motif recognizes lipid-packing defects by a conserved mechanism: peptide partitioning is driven by the insertion of hydrophobic residues into large packing defects that are preformed in the bilayer. This insertion induces only minor modifications in the statistical distribution of the free packing defects. ALPS insertion is severely hampered when monounsaturated lipids are replaced by saturated lipids, leading to a decrease in packing defects. We propose that the hypersensitivity of ALPS motifs to lipid packing defects results from the repetitive use of hydrophobic insertions along the monotonous ALPS sequence.

Magnesium deficiency phenotypes upon multiple knockout of Arabidopsis thaliana MRS2 clade B genes can be ameliorated by concomitantly reduced calcium supply

Plant MRS2 membrane protein family members have been shown to play important roles in magnesium uptake and homeostasis. Single and double knockouts for two Arabidopsis thaliana genes, AtMRS2-1 and AtMRS2-5, have previously not shown significant phenotypes even under limiting Mg(2+) supply although both are strongly expressed already in early seedlings. Together with AtMRS2-10, these genes form clade B of the AtMRS2 gene family. We now succeeded in obtaining homozygous AtMRS2-1/10 double and AtMRS2-1/5/10 triple knockout lines after selection under increased magnesium supply. Although wilting early, both new mutant lines develop fully and are also fertile under standard magnesium supply, but show severe developmental retardation under limiting Mg(2+) concentrations. To investigate nutrient dependency of germination and seedling development under various conditions, including variable supplies of Mg(2+), Ca(2+), Zn(2+), Mn(2+), Co(2+), Cd(2+) and Cu(2+), in a reproducible and economical way, we employed a small-scale liquid culturing system in 24-well plate set-ups. This allowed the growth and monitoring of individual plantlets of different mutant lines under several nutritional conditions in parallel, and the scoring and statistical evaluation of developmental stages and biomass accumulation. Detrimental effects of higher concentrations of these elements were similar in mutants and the wild type. However, growth retardation phenotypes seen upon hydroponic cultivation under low Mg(2+) could be ameliorated when Ca(2+) concentrations were concomitantly lowered, supporting indications for an important interplay of these two most abundant divalent cations in the nutrient homeostasis of plants.

Catalysis of Na+ permeation in the bacterial sodium channel Na(V)Ab

Determination of a high-resolution 3D structure of voltage-gated sodium channel Na(V)Ab opens the way to elucidating the mechanism of ion conductance and selectivity. To examine permeation of Na(+) through the selectivity filter of the channel, we performed large-scale molecular dynamics simulations of Na(V)Ab in an explicit, hydrated lipid bilayer at 0 mV in 150 mM NaCl, for a total simulation time of 21.6 μs. Although the cytoplasmic end of the pore is closed, reversible influx and efflux of Na(+) through the selectivity filter occurred spontaneously during simulations, leading to equilibrium movement of Na(+) between the extracellular medium and the central cavity of the channel. Analysis of Na(+) dynamics reveals a knock-on mechanism of ion permeation characterized by alternating occupancy of the channel by 2 and 3 Na(+) ions, with a computed rate of translocation of (6 ± 1) × 10(6) ions⋅s(-1) that is consistent with expectations from electrophysiological studies. The binding of Na(+) is intimately coupled to conformational isomerization of the four E177 side chains lining the extracellular end of the selectivity filter. The reciprocal coordination of variable numbers of Na(+) ions and carboxylate groups leads to their condensation into ionic clusters of variable charge and spatial arrangement. Structural fluctuations of these ionic clusters result in a myriad of ion binding modes and foster a highly degenerate, liquid-like energy landscape propitious to Na(+) diffusion. By stabilizing multiple ionic occupancy states while helping Na(+) ions diffuse within the selectivity filter, the conformational flexibility of E177 side chains underpins the knock-on mechanism of Na(+) permeation.

Exploring the structure and function of Thermotoga maritima CorA reveals the mechanism of gating and ion selectivity in Co2+/Mg2+ transport

The CorA family of divalent cation transporters utilizes Mg²⁺ and Co²⁺ as primary substrates. The molecular mechanism of its function, including ion selectivity and gating, has not been fully characterized. Recently we reported a new structure of a CorA homologue from Methanocaldococcus jannaschii, which provided novel structural details that offered the conception of a unique gating mechanism involving conversion of an open hydrophilic gate into a closed hydrophobic one. In the present study we report functional evidence for this novel gating mechanism in the Thermotoga maritima CorA together with an improved crystal structure of this CorA to 2.7 Å (1 Å=0.1 nm) resolution. The latter reveals the organization of the selectivity filter to be similar to that of M. jannaschii CorA and also the previously unknown organization of the second signature motif of the CorA family. The proposed gating is achieved by a helical rotation upon the binding of a metal ion substrate to the regulatory binding sites. Additionally, our data suggest that the preference of this CorA for Co²⁺ over Mg²⁺ is controlled by the presence of threonine side chains in the channel. Finally, the roles of the intracellular metal-binding sites have been assigned to increased thermostability and regulation of the gating. These mechanisms most likely apply to the entire CorA family as they are regulated by the highly conserved amino acids.

Molecular dynamics simulations of membrane-sugar interactions

It is well documented that disaccharides in general and trehalose (TRH) in particular strongly affect physical properties and functionality of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH by means of molecular dynamics (MD) computer simulations. Ten different TRH concentrations were studied in the range wTRH = 0-0.20 (w/w). The potential of mean force (PMF) for DMPC bilayer-TRH interactions was determined using two different force fields, and was subsequently used in a simple analytical model for description of sugar binding at the membrane interface. The MD results were in good agreement with the predictions of the model. The net affinities of TRH for the DMPC bilayer derived from the model and MD simulations were compared with experimental results. The area per lipid increases and the membrane becomes thinner with increased TRH concentration, which is interpreted as an intercalation effect of the TRH molecules into the polar part of the lipids, resulting in conformational changes in the chains. These results are consistent with recent experimental observations. The compressibility modulus related to the fluctuations of the membrane increases dramatically with increased TRH concentration, which indicates higher order and rigidity of the bilayer. This is also reflected in a decrease (by a factor of 15) of the lateral diffusion of the lipids. We interpret these observations as a formation of a glassy state at the interface of the membrane, which has been suggested in the literature as a hypothesis for the membrane-sugar interactions.

Detergent-mediated protein aggregation

Because detergents are commonly used to solvate membrane proteins for structural evaluation, much attention has been devoted to assessing the conformational bias imparted by detergent micelles in comparison to the native environment of the lipid bilayer. Here, we conduct six 500-ns simulations of a system with >600,000 atoms to investigate the spontaneous self assembly of dodecylphosphocholine detergent around multiple molecules of the integral membrane protein PagP. This detergent formed equatorial micelles in which acyl chains surround the protein's hydrophobic belt, confirming existing models of the detergent solvation of membrane proteins. In addition, unexpectedly, the extracellular and periplasmic apical surfaces of PagP interacted with the headgroups of detergents in other micelles 85 and 60% of the time, respectively, forming complexes that were stable for hundreds of nanoseconds. In some cases, an apical surface of one molecule of PagP interacted with an equatorial micelle surrounding another molecule of PagP. In other cases, the apical surfaces of two molecules of PagP simultaneously bound a neat detergent micelle. In these ways, detergents mediated the non-specific aggregation of folded PagP. These simulation results are consistent with dynamic light scattering experiments, which show that, at detergent concentrations ≥600 mM, PagP induces the formation of large scattering species that are likely to contain many copies of the PagP protein. Together, these simulation and experimental results point to a potentially generic mechanism of detergent-mediated protein aggregation.