GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations.

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations.

Gating Charge Calculations by Computational Electrophysiology Simulations.

Electrical cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the molecular details of voltage sensing remain insufficiently understood. Here, we present a molecular dynamics simulation-based method to determine the underlying charge movement across the membrane-the gating charge-by measuring electrical capacitor properties of membrane-embedded proteins.

Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects.

Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution.

Voltage Dependence of Conformational Dynamics and Subconducting States of VDAC-1.

The voltage-dependent anion channel 1 (VDAC-1) is an important protein of the outer mitochondrial membrane that transports energy metabolites and is involved in apoptosis. The available structures of VDAC proteins show a wide β-stranded barrel pore, with its N-terminal α-helix (N-α) bound to its interior. Electrophysiology experiments revealed that voltage, its polarity, and membrane composition modulate VDAC currents.

The membrane anchor of the transcriptional activator SREBP is characterized by intrinsic conformational flexibility.

Regulated intramembrane proteolysis (RIP) is a conserved mechanism crucial for numerous cellular processes, including signaling, transcriptional regulation, axon guidance, cell adhesion, cellular stress responses, and transmembrane protein fragment degradation. Importantly, it is relevant in various diseases including Alzheimer's disease, cardiovascular diseases, and cancers. Even though a number of structures of different intramembrane proteases have been solved recently, fundamental questions concerning mechanistic underpinnings of RIP and therapeutic interventions remain.

ClC-K chloride channels: emerging pathophysiology of Bartter syndrome type 3.

The mutations in the CLCNKB gene encoding the ClC-Kb chloride channel are responsible for Bartter syndrome type 3, one of the four variants of Bartter syndrome in the genetically based nomenclature. All forms of Bartter syndrome are characterized by hypokalemia, metabolic alkalosis, and secondary hyperaldosteronism, but Bartter syndrome type 3 has the most heterogeneous presentation, extending from severe to very mild. A relatively large number of CLCNKB mutations have been reported, including gene deletions and nonsense or missense mutations.

Mechanisms of anion conduction by coupled glutamate transporters.

Excitatory amino acid transporters (EAATs) are essential for terminating glutamatergic synaptic transmission. They are not only coupled glutamate/Na(+)/H(+)/K(+) transporters but also function as anion-selective channels. EAAT anion channels regulate neuronal excitability, and gain-of-function mutations in these proteins result in ataxia and epilepsy.

CLCNKB mutations causing mild Bartter syndrome profoundly alter the pH and Ca2+ dependence of ClC-Kb channels.

ClC-Kb, a member of the ClC family of Cl(-) channels/transporters, plays a major role in the absorption of NaCl in the distal nephron. CLCNKB mutations cause Bartter syndrome type 3, a hereditary renal salt-wasting tubulopathy. Here, we investigate the functional consequences of a Val to Met substitution at position 170 (V170M, α helix F), which was detected in eight patients displaying a mild phenotype.

Computational analysis of local membrane properties.

In the field of biomolecular simulations, dynamics of phospholipid membranes is of special interest. A number of proteins, including channels, transporters, receptors and short peptides are embedded in lipid bilayers and tightly interact with phospholipids. While the experimental measurements report on the spatial and/or temporal average membrane properties, simulation results are not restricted to the average properties.

Partial least-squares functional mode analysis: application to the membrane proteins AQP1, Aqy1, and CLC-ec1.

We introduce an approach based on the recently introduced functional mode analysis to identify collective modes of internal dynamics that maximally correlate to an external order parameter of functional interest. Input structural data can be either experimentally determined structure ensembles or simulated ensembles, such as molecular dynamics trajectories. Partial least-squares regression is shown to yield a robust solution to the multidimensional optimization problem, with a minimal and controllable risk of overfitting, as shown by extensive cross-validation.

β-Barrel mobility underlies closure of the voltage-dependent anion channel.

The voltage-dependent anion channel (VDAC) is the major protein in the outer mitochondrial membrane, where it mediates transport of ATP and ADP. Changes in its permeability, induced by voltage or apoptosis-related proteins, have been implicated in apoptotic pathways. The three-dimensional structure of VDAC has recently been determined as a 19-stranded β-barrel with an in-lying N-terminal helix.

Molecular driving forces defining lipid positions around aquaporin-0.

Lipid-protein interactions play pivotal roles in biological membranes. Electron crystallographic studies of the lens-specific water channel aquaporin-0 (AQP0) revealed atomistic views of such interactions, by providing high-resolution structures of annular lipids surrounding AQP0. It remained unclear, however, whether these lipid structures are representative of the positions of unconstrained lipids surrounding an individual protein, and what molecular determinants define the lipid positions around AQP0.

Functional dynamics in the voltage-dependent anion channel.

The voltage-dependent anion channel (VDAC), located in the outer mitochondrial membrane, acts as a gatekeeper for the entry and exit of mitochondrial metabolites. Here we reveal functional dynamics of isoform one of VDAC (VDAC1) by a combination of solution NMR spectroscopy, Gaussian network model analysis, and molecular dynamics simulation. Micro- to millisecond dynamics are significantly increased for the N-terminal six β-strands of VDAC1 in micellar solution, in agreement with increased B-factors observed in the same region in the bicellar crystal structure of VDAC1.

Voltage-dependent and -independent titration of specific residues accounts for complex gating of a ClC chloride channel by extracellular protons.

The ClC transport protein family comprises both Cl(-) ion channel and H(+)/Cl(-) and H(+)/NO(3)(-) exchanger members. Structural studies on a bacterial ClC transporter reveal a pore obstructed at its external opening by a glutamate side-chain which acts as a gate for Cl(-) passage and in addition serves as a staging post for H(+) exchange. This same conserved glutamate acts as a gate to regulate Cl(-) flow in ClC channels.