Rodolfo Briones

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. We have combined molecular dynamics simulations with fluorescence spectroscopy of the prokaryotic homolog GltPh and patch-clamp recordings of mammalian EAATs to determine how these transporters conduct anions. Whereas outward- and inward-facing GltPh conformations are nonconductive, lateral movement of the glutamate transport domain from intermediate transporter conformations results in formation of an anion-selective conduction pathway. Fluorescence quenching of inserted tryptophan residues indicated the entry of anions into this pathway, and mutations of homologous pore-forming residues had analogous effects on GltPh simulations and EAAT2/EAAT4 measurements of single-channel currents and anion/cation selectivities. These findings provide a mechanistic framework of how neurotransmitter transporters can operate as anion-selective and ligand-gated ion channels.

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. In the current study, we present a collection of methods for an efficient local membrane property calculation, comprising bilayer thickness, area per lipid, deuterium order parameters, Gaussian and mean curvature. The local membrane property calculation allows for a direct mapping of the membrane features, which subsequently can be used for further analysis and visualization of the processes of interest. The main features of the described methods are highlighted in a number of membrane systems, namely: a pure dimyristoyl-phosphatidyl-choline (DMPC) bilayer, a fusion peptide interacting with a membrane, voltage-dependent anion channel protein embedded in a DMPC bilayer, cholesterol enriched bilayer and a coarse grained simulation of a curved palmitoyl-oleoyl-phosphatidyl-choline lipid membrane. The local membrane property analysis proves to provide an intuitive and detailed view on the observables that are otherwise interpreted as averaged bilayer properties.

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. We addressed these questions by using molecular dynamics simulations and crystallographic refinement, and calculated time-averaged densities of dimyristoyl-phosphatidylcholine lipids around AQP0. Our simulations demonstrate that, although the experimentally determined crystallographic lipid positions are constrained by the crystal packing, they appropriately describe the behavior of unconstrained lipids around an individual AQP0 tetramer, and thus likely represent physiologically relevant lipid positions.While the acyl chains were well localized, the lipid head groups were not. Furthermore, in silico mutations showed that electrostatic interactions do not play a major role attracting these phospholipids towards AQP0. Instead, the mobility of the protein crucially modulates the lipid localization and explains the difference in lipid density between extracellular and cytoplasmic leaflets. Moreover, our simulations support a general mechanism in which membrane proteins laterally diffuse accompanied by several layers of localized lipids, with the positions of the annular lipids being influenced the most by the protein surface. We conclude that the acyl chains rather than the head groups define the positions of dimyristoyl-phosphatidylcholine lipids around AQP0. Lipid localization is largely determined by the mobility of the protein surface, whereas hydrogen bonds play an important but secondary role.

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. Several examples illustrate that the partial least-squares-based functional mode analysis successfully reveals the collective dynamics underlying the fluctuations in selected functional order parameters. Applications to T4 lysozyme, the Trp-cage, the aquaporin channels Aqy1 and hAQP1, and the CLC-ec1 chloride antiporter are presented in which the active site geometry, the hydrophobic solvent-accessible surface, channel gating dynamics, water permeability (p(f)), and a dihedral angle are defined as functional order parameters. The Aqy1 case reveals a gating mechanism that connects the inner channel gating residues with the protein surface, thereby providing an explanation of how the membrane may affect the channel. hAQP1 shows how the p(f) correlates with structural changes around the aromatic/arginine region of the pore. The CLC-ec1 application shows how local motions of the gating Glu(148) couple to a collective motion that affects ion affinity in the pore.

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. However, only 20 CLCNKB mutations have been functionally analyzed, due to technical difficulties regarding ClC-Kb functional expression in heterologous systems. This review provides an overview of recent progress in the functional consequences of CLCNKB mutations on ClC-Kb chloride channel activity. It has been observed that 1) all ClC-Kb mutants have an impaired expression at the membrane; and 2) a minority of the mutants combines reduced membrane expression with altered pH-dependent channel gating. Although further investigation is needed to fully characterize disease pathogenesis, Bartter syndrome type 3 probably belongs to the large family of conformational diseases, in which the mutations destabilize channel structure, inducing ClC-Kb retention in the endoplasmic reticulum and accelerated channel degradation.

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. Experiments with VDAC-1 mutants identified amino acids that regulate the gating process. However, the mechanisms for how these factors regulate VDAC-1, and which changes they trigger in the channel, are still unknown. In this study, molecular dynamics simulations and single-channel experiments of VDAC-1 show agreement for the current-voltage relationships of an “open” channel and they also show several subconducting transient states that are more cation selective in the simulations. We observed voltage-dependent asymmetric distortions of the VDAC-1 barrel and the displacement of particular charged amino acids. We constructed conformational models of the protein voltage response and the pore changes that consistently explain the protein conformations observed at opposite voltage polarities, either in phosphatidylethanolamine or phosphatidylcholine membranes. The submicrosecond VDAC-1 voltage response shows intrinsic structural changes that explain the role of key gating amino acids and support some of the current gating hypotheses. These voltage-dependent protein changes include asymmetric barrel distortion, its interaction with the membrane, and significant displacement of N-α amino acids.

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. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.

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

Significance Sterol regulatory element-binding protein (SREBP) signaling is responsible for transcriptional regulation of cellular lipid homeostasis. Aberrant regulation of the feedback loop that switches on transcription of related genes upon depletion of fatty acids and cholesterol abrogates cellular integrity and is a hallmark of many cancers and other diseases. Regulated intramembrane proteolysis is a conserved mechanism also involved in other signaling processes like Wnt and Notch and defines cellular regulation from homeostasis to proliferation, compartmentalization, and differentiation. An understanding of basic features, however, such as substrate selectivity and specificity, has remained elusive to date. Understanding the proteolysis of the SREBP precursor is an important goal both for understanding and medically intervening in SREBP signaling and related diseases as well as for understanding other regulated intramembrane proteolysis-dependent processes. 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. In particular, this includes substrate recognition, what properties render a given substrate amenable for RIP, and how the lipid environment affects the substrate cleavage. Members of the sterol regulatory element-binding protein (SREBP) family of transcription factors are critical regulators of genes involved in cholesterol/lipid homeostasis. After site-1 protease cleavage of the inactive SREBP transmembrane precursor protein, RIP of the anchor intermediate by site-2 protease generates the mature transcription factor. In this work, we have investigated the labile anchor intermediate of SREBP-1 using NMR spectroscopy. Surprisingly, NMR chemical shifts, site-resolved solvent exposure, and relaxation studies show that the cleavage site of the lipid-signaling protein intermediate bears rigid α-helical topology. An evolutionary conserved motif, by contrast, interrupts the secondary structure ∼9–10 residues C-terminal of the scissile bond and acts as an inducer of conformational flexibility within the carboxyl-terminal transmembrane region. These results are consistent with molecular dynamics simulations. Topology, stability, and site-resolved dynamics data suggest that the cleavage of the α-helical substrate in the case of RIP may be associated with a hinge motion triggered by the molecular environment.

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. We illustrate the approach by calculating the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K+ channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against experimental data. We then use the gating charge analysis to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na+ ions in Na+-coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with experiments, and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions.

Lipid Protein Interactions and Dynamical Properties of VDAC 1 channel

VDAC-1 (Voltage Dependent Anion Channel) is one of the main components of the outer mitochondrial membrane. It is responsible for the transport of ATP and other anions, and it is involved in apoptosis and cancer [1]. The X-ray and NMR structures [2-4] showed VDAC as a 19 beta-barrel structure with an N-terminal alpha-helix bound to its interior. Mutations of E73 facing the membrane could be associated with local perturbations of the surrounding lipids. Remarkably, these changes are coupled to the intrinsic VDAC dynamics. To address this we measured the the local average thickness of the membrane and several structural properties of several E73 mutants of VDAC-1. Solution NMR in LDAO and molecular dynamics of VDAC-1 inserted in DPMC phospholipid patches were carried out. Mutation or chemical modification of E73 strongly reduces the micro- to millisecond dynamics in solution. The results show the main distortions of the membrane are located around E73. The major fluctuations of VDAC barrel were found to be correlated with the charge of E73. The motion amplitude described as PCA eigenvectors show structural deformations of VDAC mainly around E73X. The motions correspond to changes of the whole beta-barrel structure and align with the position of E73X. These results help to understand the intrinsic dynamic of VDAC and its possible interaction mechanism with lipid membranes.[1] Zaid et al., Cell Death Differ. 12, 751 (2005)[2] Hiller et al., Science 321, 1208 (2008)[3] Ujwal et al., PNAS 105, 17742 (2008)[4] Bayrhuber et al., PNAS 105, 15370 (2008)