Asymmetric CorA Gating Mechanism as Observed by Molecular Dynamics Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Nemchinova M
PROVIDER: S-EPMC8154316 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:CorA is the major transport system responsible for Mg(2+) 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(2+) uptake remains to be established. Here we use electron paramagnetic resonance spectroscopy, electrophysiology and molecular dynamic simulations to show that CorA is regulated by cytoplasmic Mg(2+) acting as a ligand and elucidate the basic conformational rearrangements responsible for Mg(2+)-dependent gating. Mg(2+) 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(2+)-driven negative feedback loop in CorA as the key physiological event controlling Mg(2+) uptake and homeostasis in prokaryotes.
| S-EPMC4066822 | biostudies-literature
Project description:Aquaporins are protein channels located across the cell membrane with the role of conducting water or other small sugar alcohol molecules (aquaglyceroporins). The high-resolution X-ray structure of the human aquaporin 5 (HsAQP5) shows that HsAQP5, as all the other known aquaporins, exhibits tetrameric structure. By means of molecular dynamics simulations we analyzed the role of spontaneous fluctuations on the structural behavior of the human AQP5. We found that different conformations within the tetramer lead to a distribution of monomeric channel structures, which can be characterized as open or closed. The switch between the two states of a channel is a tap-like mechanism at the cytoplasmic end which regulates the water passage through the pore. The channel is closed by a translation of the His67 residue inside the pore. Moreover, water permeation rate calculations revealed that the selectivity filter, located at the other end of the channel, regulates the flow rate of water molecules when the channel is open, by locally modifying the orientation of His173. Furthermore, the calculated permeation rates of a fully open channel are in good agreement with the reported experimental value.
| S-EPMC3614956 | biostudies-literature
Project description:Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to fatal cardiac arrhythmias such as catecholaminergic polymorphic ventricular tachycardia (CPVT). While many CPVT mutations are associated with an increase in Ca2+ leak from the sarcoplasmic reticulum, the mechanistic details of RyR2 channel gating are not well understood, and this poses a barrier in the development of new pharmacological treatments. To address this, we explore the gating mechanism of the RyR2 using molecular dynamics (MD) simulations. We test the effect of changing the conformation of certain structural elements by constructing chimera RyR2 structures that are derived from the currently available closed and open cryo-electron microscopy (cryo-EM) structures, and we then use MD simulations to relax the system. Our key finding is that the position of the S4-S5 linker (S4S5L) on a single subunit can determine whether the channel as a whole is open or closed. Our analysis reveals that the position of the S4S5L is regulated by interactions with the U-motif on the same subunit and with the S6 helix on an adjacent subunit. We find that, in general, channel gating is crucially dependent on high percent occupancy interactions between adjacent subunits. We compare our interaction analysis to 49 CPVT1 mutations in the literature and find that 73% appear near a high percent occupancy interaction between adjacent subunits. This suggests that disruption of cooperative, high percent occupancy interactions between adjacent subunits is a primary cause of channel leak and CPVT in mutant RyR2 channels.
| S-EPMC9720719 | biostudies-literature
Project description:Gastric H,K-ATPase is an electroneutral transmembrane pump that moves protons from the cytoplasm of the parietal cell into the gastric lumen in exchange for potassium ions. The mechanism of transport against the established electrochemical gradients includes intermediate conformations in which the transferred ions are trapped (occluded) within the membrane domain of the pump. The pump cycle involves switching between the E1 and E2P states. Molecular dynamics simulations on homology models of the E2P and E1 states were performed to investigate the mechanism of K(+) movement in this enzyme. We performed separate E2P simulations with one K(+) in the luminal channel, one K(+) ion in the occlusion site, two K(+) ions in the occlusion site, and targeted molecular dynamics from E2P to E1 with two K(+) ions in the occlusion site. The models were inserted into a lipid bilayer system and were stable over the time course of the simulations, and K(+) ions in the channel moved to a consistent location near the center of the membrane domain, thus defining the occlusion site. The backbone carbonyl oxygen from residues 337 through 342 on the nonhelical turn of M4, as well as side-chain oxygen from E343, E795, and E820, participated in the ion occlusion. A single water molecule was stably bound between the two K(+) ions in the occlusion site, providing an additional ligand and partial shielding the positive charges from one another. Targeted molecular dynamics was used to transform the protein from the E2P to the E1 state (two K(+) ions to the cytoplasm). This simulation identified the separation of the water column in the entry channel as the likely gating mechanism on the luminal side. A hydrated exit channel also formed on the cytoplasmic side of the occlusion site during this simulation. Hence, water molecules became available to hydrate the ions. The movement of the M1M2 transmembrane segments, and the displacement of residues Q159, E160, Q110, and T152 during the conformational change, as well as the motions of E343 and L346, acted as the cytoplasmic-side gate.
| S-EPMC2527272 | biostudies-literature
Project description:P2X receptors are trimeric ATP-gated cation channels participating in diverse physiological processes. How ATP binding triggers channel opening remains unclear. Here the gating mechanism of a P2X receptor was studied by normal mode analysis and molecular dynamics (MD) simulations. Based on the resting-state crystal structure, a normal mode involving coupled motions of three β-strands (β1, β13, and β14) at the trimeric interface of the ligand-binding ectodomain and the pore-lining helix (TM2) in the transmembrane domain (TMD) was identified. The resulting widening of the fenestrations above the TMD and opening of the transmembrane pore produce known signatures of channel activation. In MD simulations, ATP was initially placed in the putative binding pocket (defined by four charged residues located in β1, β13 and β14) in two opposite orientations, with the adenine either proximal or distal to the TMD. In the proximal orientation, the triphosphate group extends outward to draw in the four charged residues, leading to closure of β13/β14 toward β1. The adenine ring, wedged between β1 and β13, acts as a fulcrum for the β14 lever, turning a modest closure around the triphosphate group into significant opening of the pre-TM2 loop. The motions of these β-strands are similar to those in the putative channel-activation normal mode. In the distal orientation, the ATP stabilizes the trimeric interface and the closure of the pre-TM2 loop, possibly representing desensitization. Our computational studies produced the first complete model, supported by experimental data, for how ATP binding triggers activation of a P2X receptor.
| S-EPMC3306669 | biostudies-literature
Project description:CorA, a divalent-selective channel in the metal ion transport superfamily, is the major Mg2+-influx pathway in prokaryotes. CorA structures in closed (Mg2+-bound), and open (Mg2+-free) states, together with functional data showed that Mg2+-influx inhibits further Mg2+-uptake completing a regulatory feedback loop. While the closed state structure is a symmetric pentamer, the open state displayed unexpected asymmetric architectures. Using high-speed atomic force microscopy (HS-AFM), we explored the Mg2+-dependent gating transition of single CorA channels: HS-AFM movies during Mg2+-depletion experiments revealed the channel's transition from a stable Mg2+-bound state over a highly mobile and dynamic state with fluctuating subunits to asymmetric structures with varying degree of protrusion heights from the membrane. Our data shows that at Mg2+-concentration below Kd, CorA adopts a dynamic (putatively open) state of multiple conformations that imply structural rearrangements through hinge-bending in TM1. We discuss how these structural dynamics define the functional behavior of this ligand-dependent channel.
| S-EPMC6927688 | biostudies-literature
Project description:Serine proteases are involved in many fundamental physiological processes, and control of their activity mainly results from the fact that they are synthetized in an inactive form that becomes active upon cleavage. Three decades ago Martin Karplus's group performed the first molecular dynamics simulations of trypsin, the most studied member of the serine protease family, to address the transition from the zymogen to its active form. Based on the computational power available at the time, only high frequency fluctuations, but not the transition steps, could be observed. By performing accelerated molecular dynamics (aMD) simulations, an interesting approach that increases the configurational sampling of atomistic simulations, we were able to observe the N-terminal tail insertion, a crucial step of the transition mechanism. Our results also support the hypothesis that the hydrophobic effect is the main force guiding the insertion step, although substantial enthalpic contributions are important in the activation mechanism. As the N-terminal tail insertion is a conserved step in the activation of serine proteases, these results afford new perspective on the underlying thermodynamics of the transition from the zymogen to the active enzyme.
| S-EPMC4241106 | biostudies-literature
Project description:A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific ?-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26-27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K(+) channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K(+) and Na(+) concentration gradients and the electrostatic potential difference of living cells.
| S-EPMC3787774 | biostudies-literature
Project description:One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.
| S-EPMC3508910 | biostudies-literature
Project description:Loop-gating is one of two voltage-dependent mechanisms that regulate the open probability of connexin channels. The loop-gate permeability barrier is formed by a segment of the first extracellular loop (E1) (the parahelix) and appears to be accompanied by straightening of the bend angle between E1 and the first transmembrane domain (TM1). Here, all-atom molecular dynamics simulations are used to identify and characterize interacting van der Waals and electrostatic networks that stabilize the parahelices and TM1/E1 bend angles of the open Cx26 hemichannel. Dynamic fluctuations in an electrostatic network in each subunit are directly linked to the stability of parahelix structure and TM1/E1 bend angle in adjacent subunits. The electrostatic network includes charged residues that are pore-lining and thus positioned to be voltage sensors. We propose that the transition to the closed state is initiated by voltage-driven disruption of the networks that stabilize the open-state parahelix configuration, allowing the parahelix to protrude into the channel pore to form the loop-gate barrier. Straightening of the TM1/E1 bend appears to be a consequence of the reorganization of the interacting networks that accompany the conformational change of the parahelix. The electrostatic network extends across subunit boundaries, suggesting a concerted gating mechanism.
| S-EPMC3309406 | biostudies-literature