Search

The voltage-gated proton channel Hv1 mediates efflux of protons from the cell. Hv1 integrally contributes to various physiological processes including pH homeostasis and the respiratory burst of phagocytes. Inhibition of Hv1 may provide therapeutic avenues for the treatment of inflammatory diseases, breast cancer, and ischemic brain damage. In this work, we investigate two prototypical Hv1 inhibitors, 2-guanidinobenzimidazole (2GBI), and 5-chloro-2-guanidinobenzimidazole (GBIC), from an experimentally screened class of guanidine derivatives. Both compounds block proton conduction by binding the same site located on the intracellular side of the channel. However, when added to the extracellular medium, the compounds strongly differ in their ability to inhibit proton conduction, suggesting substantial differences in membrane permeability. Here, we compute the potential of mean force for each compound to permeate through the membrane using atomistic molecular dynamics simulations with the adaptive biasing force method. Our results rationalize the putative distinction between these two blockers with respect to their abilities to permeate the cellular membrane.

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Hv1 is a voltage-gated proton channel whose main function is to facilitate extrusion of protons from the cell. The development of effective channel blockers for Hv1 can lead to new therapeutics for the treatment of maladies related to Hv1 dysfunction. Although the mechanism of proton permeation in Hv1 remains to be elucidated, a series of small molecules have been discovered to inhibit Hv1. Here, we computed relative binding free energies of a prototypical Hv1 blocker on a model of human Hv1 in an open state. We used alchemical free energy perturbation techniques based on atomistic molecular dynamics simulations. The results support our proposed open state model and shed light on the preferred tautomeric state of the channel blocker. This work lays the groundwork for future studies on adapting the blocker molecule for more effective inhibition of Hv1.

The voltage-gated Hv1 proton channel is a ubiquitous membrane protein that has roles in a variety of cellular processes, including proton extrusion, pH regulation, production of reactive oxygen species, proliferation of cancer cells, and increased brain damage during ischemic stroke. A crystal structure of an Hv1 construct in a putative closed state has been reported, and structural models for the channel open state have been proposed, but a complete characterization of the Hv1 conformational dynamics under an applied membrane potential has been elusive. We report structural models of the Hv1 voltage-sensing domain (VSD), both in a hyperpolarized state and a depolarized state resulting from voltage-dependent conformational changes during a 10-μs-timescale atomistic molecular dynamics simulation in an explicit membrane environment. In response to a depolarizing membrane potential, the S4 helix undergoes an outward displacement, leading to changes in the VSD internal salt-bridge network, resulting in a reshaping of the permeation pathway and a significant increase in hydrogen bond connectivity throughout the channel. The total gating charge displacement associated with this transition is consistent with experimental estimates. Molecular docking calculations confirm the proposed mechanism for the inhibitory action of 2-guanidinobenzimidazole (2GBI) derived from electrophysiological measurements and mutagenesis. The depolarized structural model is also consistent with the formation of a metal bridge between residues located in the core of the VSD. Taken together, our results suggest that these structural models are representative of the closed and open states of the Hv1 channel.

BackgroundThe eye lens crystallins are highly soluble proteins that are required to last the lifespan of an organism due to low protein turnover in the lens. Crystallin aggregation leads to formation of light-scattering aggregates known as cataract. The G18V mutation of human γS-crystallin (γS-G18V), which is associated with childhood-onset cataract, causes structural changes throughout the N-terminal domain and increases aggregation propensity. The holdase chaperone protein αB-crystallin does not interact with wild-type γS-crystallin, but does bind its G18V variant. The specific molecular determinants of αB-crystallin binding to client proteins is incompletely charcterized. Here, a new variant of γS, γS-G18A, was created to test the limits of αB-crystallin selectivity.MethodsMolecular dynamics simulations were used to investigate the structure and dynamics of γS-G18A. The overall fold of γS-G18A was assessed by circular dichroism (CD) spectroscopy and intrinsic tryptophan fluorescence. Its thermal unfolding temperature and aggregation propensity were characterized by CD and DLS, respectively. Solution-state NMR was used to characterize interactions between αB-crystallin and γS-G18A.ResultsγS-G18A exhibits minimal structural changes, but has compromised thermal stability relative to γS-WT. The placement of alanine, rather than valine, at this highly conserved glycine position produces minor changes in hydrophobic surface exposure. However, human αB-crystallin does not bind the G18A variant, in contrast to previous observations for γS-G18V, which aggregates at physiological temperature.ConclusionsαB-crystallin is capable of distinguishing between aggregation-prone and function-preserving variants, and recognizing the transient unfolding or minor conformers that lead to aggregation in the disease-related variant.General significanceHuman αB-crystallin distinguishes between highly similar variants of a structural crystallin, binding the cataract-related γS-G18V variant, but not the function-preserving γS-G18A variant, which is monomeric at physiological temperature.

Available experimental techniques cannot determine high-resolution three-dimensional structures of membrane proteins under a transmembrane voltage. Hence, the mechanism by which voltage-gated cation channels couple conformational changes within the four voltage sensor domains, in response to either depolarizing or polarizing transmembrane voltages, to opening or closing of the pore domain's ion channel remains unresolved. Single-membrane specimens, composed of a phospholipid bilayer containing a vectorially oriented voltage-gated K+ channel protein at high in-plane density tethered to the surface of an inorganic multilayer substrate, were developed to allow the application of transmembrane voltages in an electrochemical cell. Time-resolved neutron reflectivity experiments, enhanced by interferometry enabled by the multilayer substrate, were employed to provide directly the low-resolution profile structures of the membrane containing the vectorially oriented voltage-gated K+ channel for the activated, open and deactivated, closed states of the channel under depolarizing and hyperpolarizing transmembrane voltages applied cyclically. The profile structures of these single membranes were dominated by the voltage-gated K+ channel protein because of the high in-plane density. Importantly, the use of neutrons allowed the determination of the voltage-dependent changes in both the profile structure of the membrane and the distribution of water within the profile structure. These two key experimental results were then compared to those predicted by three computational modeling approaches for the activated, open and deactivated, closed states of three different voltage-gated K+ channels in hydrated phospholipid bilayer membrane environments. Of the three modeling approaches investigated, only one state-of-the-art molecular dynamics simulation that directly predicted the response of a voltage-gated K+ channel within a phospholipid bilayer membrane to applied transmembrane voltages by utilizing very long trajectories was found to be in agreement with the two key experimental results provided by the time-resolved neutron interferometry experiments.

Aquaporin 0 (AQP0) is essential for eye lens homeostasis as is regulation of its water permeability by Ca2+, which occurs through interactions with calmodulin (CaM), but the underlying molecular mechanisms are not well understood. Here, we use molecular dynamics (MD) simulations on the microsecond timescale under an osmotic gradient to explicitly model water permeation through the AQP0 channel. To identify any structural features that are specific to water permeation through AQP0, we also performed simulations of aquaporin 1 (AQP1) and a pure mixed lipid bilayer under the same conditions. The relative single-channel water osmotic permeability coefficients (pf) calculated from all of our simulations are in reasonable agreement with experiment. Our simulations allowed us to characterize the dynamics of the key structural elements that modulate the diffusion of water single-files through the AQP0 and AQP1 pores. We find that CaM binding influences the collective dynamics of the whole AQP0 tetramer, promoting the closing of both the extracellular and intracellular gates by inducing cooperativity between neighboring subunits.

Dynamic disorder of the lipid bilayer presents a challenge for establishing structure-function relationships in membranous systems. The resulting structural heterogeneity is especially evident for peripheral and spontaneously inserting membrane proteins, which are not constrained by the well-defined transmembrane topology and exert their action in the context of intimate interaction with lipids. Here, we propose a concerted approach combining depth-dependent fluorescence quenching with Molecular Dynamics simulation to decipher dynamic interactions of membrane proteins with the lipid bilayers. We apply this approach to characterize membrane-mediated action of the diphtheria toxin translocation domain. First, we use a combination of the steady-state and time-resolved fluorescence spectroscopy to characterize bilayer penetration of the NBD probe selectively attached to different sites of the protein into membranes containing lipid-attached nitroxyl quenching groups. The constructed quenching profiles are analyzed with the Distribution Analysis methodology allowing for accurate determination of transverse distribution of the probe. The results obtained for 12 NBD-labeled single-Cys mutants are consistent with the so-called Open-Channel topology model. The experimentally determined quenching profiles for labeling sites corresponding to L350, N373, and P378 were used as initial constraints for positioning TH8-9 hairpin into the lipid bilayer for Molecular Dynamics simulation. Finally, we used alchemical free energy calculations to characterize protonation of E362 in soluble translocation domain and membrane-inserted conformation of its TH8-9 fragment. Our results indicate that membrane partitioning of the neutral E362 is more favorable energetically (by ~ 6 kcal/mol), but causes stronger perturbation of the bilayer, than the charged E362.

It is now well established by numerous experimental and computational studies that the adsorption propensities of inorganic anions conform to the Hofmeister series. The adsorption propensities of inorganic cations, such as the alkali metal cations, have received relatively little attention. Here we use a combination of liquid-jet X-ray photoelectron experiments and molecular dynamics simulations to investigate the behavior of K+ and Li+ ions near the interfaces of their aqueous solutions with halide ions. Both the experiments and the simulations show that Li+ adsorbs to the aqueous solution-vapor interface, while K+ does not. Thus, we provide experimental validation of the "surfactant-like" behavior of Li+ predicted by previous simulation studies. Furthermore, we use our simulations to trace the difference in the adsorption of K+ and Li+ ions to a difference in the resilience of their hydration shells.

Hydrophobic amino acids are abundant in transmembrane (TM) helices of membrane proteins. Charged residues are sparse, apparently due to the unfavorable energetic cost of partitioning charges into nonpolar phases. Nevertheless, conserved arginine residues within TM helices regulate vital functions, such as ion channel voltage gating and integrin receptor inactivation. The energetic cost of arginine in various positions along hydrophobic helices has been controversial. Potential of mean force (PMF) calculations from atomistic molecular dynamics simulations predict very large energetic penalties, while in vitro experiments with Sec61 translocons indicate much smaller penalties, even for arginine in the center of hydrophobic TM helices. Resolution of this conflict has proved difficult, because the in vitro assay utilizes the complex Sec61 translocon, while the PMF calculations rely on the choice of simulation system and reaction coordinate. Here we present the results of computational and experimental studies that permit direct comparison with the Sec61 translocon results. We find that the Sec61 translocon mediates less efficient membrane insertion of Arg-containing TM helices compared with our computational and experimental bilayer-insertion results. In the simulations, a combination of arginine snorkeling, bilayer deformation, and peptide tilting is sufficient to lower the penalty of Arg insertion to an extent such that a hydrophobic TM helix with a central Arg residue readily inserts into a model membrane. Less favorable insertion by the translocon may be due to the decreased fluidity of the endoplasmic reticulum (ER) membrane compared with pure palmitoyloleoyl-phosphocholine (POPC). Nevertheless, our results provide an explanation for the differences between PMF- and experiment-based penalties for Arg burial.

The YidC/Oxa1/Alb3 family of membrane proteins function to insert proteins into membranes in bacteria, mitochondria, and chloroplasts. Recent X-ray structures of YidC from Bacillus halodurans and Escherichia coli revealed a hydrophilic groove that is accessible from the lipid bilayer and the cytoplasm. Here, we explore the water accessibility within the conserved core region of the E. coli YidC using in vivo cysteine alkylation scanning and molecular dynamics (MD) simulations of YidC in POPE/POPG membranes. As expected from the structure, YidC possesses an aqueous membrane cavity localized to the membrane inner leaflet. Both the scanning data and the MD simulations show that the lipid-exposed transmembrane helices 3, 4, and 5 are short, leading to membrane thinning around YidC. Close examination of the MD data reveals previously unrecognized structural features that are likely important for protein stability and function.

The dimeric 44-residue E5 protein of bovine papillomavirus is the smallest known naturally occurring oncoprotein. This transmembrane protein binds to the transmembrane domain (TMD) of the platelet-derived growth factor β receptor (PDGFβR), causing dimerization and activation of the receptor. Here, we use Rosetta membrane modeling and all-atom molecular dynamics simulations in a membrane environment to develop a chemically detailed model of the E5 protein/PDGFβR complex. In this model, an active dimer of the PDGFβR TMD is sandwiched between two dimers of the E5 protein. Biochemical experiments showed that the major PDGFβR TMD complex in mouse cells contains two E5 dimers and that binding the PDGFβR TMD to the E5 protein is necessary and sufficient to recruit both E5 dimers into the complex. These results demonstrate how E5 binding induces receptor dimerization and define a molecular mechanism of receptor activation based on specific interactions between TMDs.

Voltage-sensing domains (VSDs) sense changes in the membrane electrostatic potential and, through conformational changes, regulate a specific function. The VSDs of wild-type voltage-dependent K+, Na+, and Ca2+ channels do not conduct ions, but they can become ion-permeable through pathological mutations in the VSD. Relatively little is known about the underlying mechanisms of conduction through VSDs. The most detailed studies have been performed on Shaker K+ channel variants in which ion conduction through the VSD is manifested in electrophysiology experiments as a voltage-dependent inward current, the so-called omega current, which appears when the VSDs are in their resting state conformation. Only monovalent cations appear to permeate the Shaker VSD via a pathway that is believed to be, at least in part, the same as that followed by the S4 basic side chains during voltage-dependent activation. We performed μs-time scale atomistic molecular dynamics simulations of a cation-conducting variant of the Shaker VSD under applied electric fields in an experimentally validated resting-state conformation, embedded in a lipid bilayer surrounded by solutions containing guanidinium chloride or potassium chloride. Our simulations provide insights into the Shaker VSD permeation pathway, the protein-ion interactions that control permeation kinetics, and the mechanism of voltage-dependent activation of voltage-gated ion channels.

Human voltage-gated proton (hHv1) channels are crucial for regulating essential biological processes such as immune cell respiratory burst, sperm capacitation, and cancer cell migration. Despite the significant concentration difference between protons and other ions in physiological conditions, hHv1 demonstrates remarkable proton selectivity. Our calculations of single-proton, cation, and anion permeation free energy profiles quantitatively demonstrate that the proton selectivity of the wild-type channel originates from its strong proton affinity via the titration of the key residues D112 and D174, although the channel imposes similar kinetic blocking effects for protons compared to other ions. A two-proton knock-on model is proposed to mathematically explain the electrophysiological measurements of the pH-dependent proton conductance in the conductive state. Moreover, it is shown that the anion selectivity of the D112N mutant channel is tied to impaired proton transport and substantial anion leakage. [ABSTRACT FROM AUTHOR]

Aquaporin 0 (AQP0), the major intrinsic protein of the eye lens, plays a vital role in maintaining lens clarity by facilitating the transport of water across lens fiber cell membranes. AQP0 reduces its osmotic water permeability constant (Pf) in response to increases in the external calcium concentration, an effect that is mediated by an interaction with the calcium-binding messenger protein, calmodulin (CaM), and phosphorylation of the CaM-binding site abolishes calcium sensitivity. Despite recent structural characterization of the AQP0-CaM complex, the mechanism by which CaM modulates AQP0 remains poorly understood. By combining atomistic molecular dynamics simulations and oocyte permeability assays, we conclude that serine phosphorylation of AQP0 does not inhibit CaM binding to the whole AQP0 protein. Instead, AQP0 phosphorylation alters calcium sensitivity by modifying the AQP0-CaM interaction interface, particularly at an arginine-rich loop that connects the fourth and fifth transmembrane helices. This previously unexplored loop, which sits outside of the canonical CaM-binding site on the AQP0 cytosolic face, mechanically couples CaM to the pore-gating residues of the second constriction site. We show that this allosteric loop is vital for CaM regulation of the channels, facilitating cooperativity between adjacent subunits and regulating factors such as serine phosphorylation. Similar allosteric interactions may also mediate CaM modulation of the properties of other CaM-regulated proteins.

Organized as bilayers, phospholipids are the fundamental building blocks of cellular membranes and determine many of their biological functions. Interactions between the two leaflets of the bilayer (interleaflet coupling) have been implicated in the passage of information through membranes. However, physically, the meaning of interleaflet coupling is ill defined and lacks a structural basis. Using all-atom molecular dynamics simulations of fluid phospholipid bilayers of five different lipids with differing degrees of acyl-chain asymmetry, we have examined interleaflet mixing to gain insights into coupling. Reasoning that the transbilayer distribution of terminal methyl groups is an appropriate measure of interleaflet mixing, we calculated the transbilayer distributions of the acyl chain terminal methyl groups for five lipids: dioleoylphosphatidylcholine (DOPC), palmitoyloleoylphosphatidylcholine (POPC), stearoyloleoylphosphatidylcholine (SOPC), oleoylmyristoylphosphatidylcholine (OMPC), and dimyristoylphosphatidylcholine (DMPC). We observed in all cases very strong mixing across the bilayer midplane that diminished somewhat with increasing acyl-chain ordering defined by methylene order parameters. A hallmark of the interleaflet coupling idea is complementarity, which postulates that lipids with short alkyl chains in one leaflet will preferentially associate with lipids with long alkyl chains in the other leaflet. Our results suggest a much more complicated picture for thermally disordered bilayers that we call distributed complementarity, as measured by the difference in the peak positions of the sn-1 and sn-2 methyl distributions in the same leaflet.

Carnivorous plants represent a so far underexploited reservoir of novel proteases with potentially useful activities. Here we investigate 44 cysteine proteases from the Cape sundew, Drosera capensis, predicted from genomic DNA sequences. D. capensis has a large number of cysteine protease genes; analysis of their sequences reveals homologs of known plant proteases, some of which are predicted to have novel properties. Many functionally significant sequence and structural features are observed, including targeting signals and occluding loops. Several of the proteases contain a new type of granulin domain. Although active site residues are conserved, the sequence identity of these proteases to known proteins is moderate to low; therefore, comparative modeling with all-atom refinement and subsequent atomistic MD-simulation is used to predict their 3D structures. The structure prediction data, as well as analysis of protein structure networks, suggest multifarious variations on the papain-like cysteine protease structural theme. This in silico methodology provides a general framework for investigating a large pool of sequences that are potentially useful for biotechnology applications, enabling informed choices about which proteins to investigate in the laboratory.

Upon endoplasmic reticulum Ca(2+) store depletion, Orai channels in the plasma membrane are activated directly by endoplasmic reticulum-resident STIM proteins to generate the Ca(2+)-selective, Ca(2+) release-activated Ca(2+) (CRAC) current. After the molecular identification of Orai, a plethora of functional and biochemical studies sought to compare Orai homologs, determine their stoichiometry, identify structural domains responsible for the biophysical fingerprint of the CRAC current, identify the physiological functions, and investigate Orai homologs as potential therapeutic targets. Subsequently, the solved crystal structure of Drosophila Orai (dOrai) substantiated many findings from structure-function studies, but also revealed an unexpected hexameric structure. In this review, we explore Orai channels as elucidated by functional and biochemical studies, analyze the dOrai crystal structure and its implications for Orai channel function, and present newly available information from molecular dynamics simulations that shed light on Orai channel gating and permeation.

The voltage sensor domain (VSD) of voltage-gated cation (e.g., Na(+), K(+)) channels central to neurological signal transmission can function as a distinct module. When linked to an otherwise voltage-insensitive, ion-selective membrane pore, the VSD imparts voltage sensitivity to the channel. Proteins homologous with the VSD have recently been found to function themselves as voltage-gated proton channels or to impart voltage sensitivity to enzymes. Determining the conformational changes associated with voltage gating in the VSD itself in the absence of a pore domain thereby gains importance. We report the direct measurement of changes in the scattering-length density (SLD) profile of the VSD protein, vectorially oriented within a reconstituted phospholipid bilayer membrane, as a function of the transmembrane electric potential by time-resolved X-ray and neutron interferometry. The changes in the experimental SLD profiles for both polarizing and depolarizing potentials with respect to zero potential were found to extend over the entire length of the isolated VSD's profile structure. The characteristics of the changes observed were in qualitative agreement with molecular dynamics simulations of a related membrane system, suggesting an initial interpretation of these changes in terms of the VSD's atomic-level 3-D structure.

Calmodulin (CaM) is a universal regulatory protein that communicates the presence of calcium to its molecular targets and correspondingly modulates their function. This key signaling protein is important for controlling the activity of hundreds of membrane channels and transporters. However, understanding of the structural mechanisms driving CaM regulation of full-length membrane proteins has remained elusive. In this study, we determined the pseudoatomic structure of full-length mammalian aquaporin-0 (AQP0, Bos taurus) in complex with CaM, using EM to elucidate how this signaling protein modulates water-channel function. Molecular dynamics and functional mutation studies reveal how CaM binding inhibits AQP0 water permeability by allosterically closing the cytoplasmic gate of AQP0. Our mechanistic model provides new insight, only possible in the context of the fully assembled channel, into how CaM regulates multimeric channels by facilitating cooperativity between adjacent subunits.

Several laboratories have carried out molecular dynamics (MD) simulations of arginine interactions with lipid bilayers and found that the energetic cost of placing arginine in lipid bilayers is an order of magnitude greater than observed in molecular biology experiments in which Arg-containing transmembrane helices are inserted across the endoplasmic reticulum membrane by the Sec61 translocon. We attempt here to reconcile the results of the two approaches. We first present MD simulations of guanidinium groups alone in lipid bilayers, and then, to mimic the molecular biology experiments, we present simulations of hydrophobic helices containing single Arg residues at different positions along the helix. We discuss the simulation results in the context of molecular biology results and show that the energetic discrepancy is reduced, but not eliminated, by considering free energy differences between Arg at the interface and at the center of the model helices. The reduction occurs because Arg snorkeling to the interface prevents Arg from residing in the bilayer center where the energetic cost of desolvation is highest. We then show that the problem with MD simulations is that they measure water-to-bilayer free energies, whereas the molecular biology experiments measure the energetics of partitioning from translocon to bilayer, which raises the fundamental question of the relationship between water-to-bilayer and water-to-translocon partitioning. We present two thermodynamic scenarios as a foundation for reconciliation of the simulation and molecular biology results. The simplest scenario is that translocon-to-bilayer partitioning is independent of water-to-bilayer partitioning; there is no thermodynamic cycle connecting the two paths.

Despite the growing number of atomic-resolution membrane protein structures, direct structural information about proteins in their native membrane environment is scarce. This problem is particularly relevant in the case of the highly charged S1-S4 voltage-sensing domains responsible for nerve impulses, where interactions with the lipid bilayer are critical for the function of voltage-activated ion channels. Here we use neutron diffraction, solid-state nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics simulations to investigate the structure and hydration of bilayer membranes containing S1-S4 voltage-sensing domains. Our results show that voltage sensors adopt transmembrane orientations and cause a modest reshaping of the surrounding lipid bilayer, and that water molecules intimately interact with the protein within the membrane. These structural findings indicate that voltage sensors have evolved to interact with the lipid membrane while keeping energetic and structural perturbations to a minimum, and that water penetrates the membrane, to hydrate charged residues and shape the transmembrane electric field.

The molecular basis underlying the rich phase behavior of globular proteins remains poorly understood. We use atomistic multiscale molecular simulations to model the solution‐state conformational dynamics and interprotein interactions of γD‐crystallin and its P23T‐R36S mutant, which drastically limits the protein solubility, at both infinite dilution and at a concentration where the mutant fluid phase and crystalline phase coexist. We find that while the mutant conserves the protein fold, changes to the surface exposure of residues in the neighborhood of residue‐36 enhance protein–protein interactions and develop specific protein–protein contacts found in the protein crystal lattice. [ABSTRACT FROM AUTHOR]