Human PIEZO1: Removing Inactivation

PIEZO1 is an inactivating eukaryotic cation-selective mechanosensitive ion channel. Two sites have been located in the channel that when individually mutated lead to xerocytotic anemia by slowing inactivation. By introducing mutations at two sites, one associated with xerocytosis and the other artificial, we were able to remove inactivation. The double mutant (DhPIEZO1) has a substitution of arginine for methionine (M2225R) and lysine for arginine (R2456K). The loss of inactivation was accompanied by ∼30-mmHg shift of the activation curve to lower pressures and slower rates of deactivation. The slope sensitivity of gating was the same for wild-type and mutants, indicating that the dimensional changes between the closed and open state are unaffected by the mutations. The unitary channel conductance was unchanged by mutations, so these sites are not associated with pore. DhPIEZO1 was reversibly inhibited by the peptide GsMTx4 that acted as a gating modifier. The channel kinetics were solved using complex stimulus waveforms and the data fit to a three-state loop in detailed balance. The reaction had two pressure-dependent rates, closed to open and inactivated to closed. Pressure sensitivity of the opening rate with no sensitivity of the closing rate means that the energy barrier between them is located near the open state. Mutant cycle analysis of inactivation showed that the two sites interacted strongly, even though they are postulated to be on opposite sides of the membrane.     

Human PIEZO1: Removing Inactivation

PIEZO1 is an inactivating eukaryotic cation-selective mechanosensitive ion channel. Two sites have been located in the channel that when individually mutated lead to xerocytotic anemia by slowing inactivation. By introducing mutations at two sites, one associated with xerocytosis and the other artificial, we were able to remove inactivation. The double mutant (DhPIEZO1) has a substitution of arginine for methionine (M2225R) and lysine for arginine (R2456K). The loss of inactivation was accompanied by ∼30-mmHg shift of the activation curve to lower pressures and slower rates of deactivation. The slope sensitivity of gating was the same for wild-type and mutants, indicating that the dimensional changes between the closed and open state are unaffected by the mutations. The unitary channel conductance was unchanged by mutations, so these sites are not associated with pore. DhPIEZO1 was reversibly inhibited by the peptide GsMTx4 that acted as a gating modifier. The channel kinetics were solved using complex stimulus waveforms and the data fit to a three-state loop in detailed balance. The reaction had two pressure-dependent rates, closed to open and inactivated to closed. Pressure sensitivity of the opening rate with no sensitivity of the closing rate means that the energy barrier between them is located near the open state. Mutant cycle analysis of inactivation showed that the two sites interacted strongly, even though they are postulated to be on opposite sides of the membrane.     

Glycerol inhibits water permeation through Plasmodium Falciparum aquaglyceroporin

Plasmodium falciparum aquaglyceroporin (PfAQP) is a multifunctional membrane protein in the plasma membrane of P. falciparum, the parasite that causes the most severe form of malaria. The current literature has established the science of PfAQP’s structure, functions, and hydrogen-bonding interactions but left unanswered the following fundamental question: does glycerol modulate water permeation through aquaglyceroporin that conducts both glycerol and water? This paper provides an affirmative answer to this question of essential importance to the protein’s functions. On the basis of the chemical-potential profile of glycerol from the extracellular bulk region, throughout PfAQP’s conducting channel, to the cytoplasmic bulk region, this study shows the existence of a bound state of glycerol inside aquaglyceroporin’s permeation pore, from which the dissociation constant is approximately 14 μM. A glycerol molecule occupying the bound state occludes the conducting pore through which permeating molecules line up in single file by hydrogen-bonding with one another and with the luminal residues of aquaglyceroporin. In this way, glycerol inhibits permeation of water and other permeants through aquaglyceroporin. The biological implications of this theory are discussed and shown to agree with the existent in vitro data. It turns out that the structure of aquaglyceroporin is perfect for the van der Waals interactions between the protein and glycerol to cause the existence of the bound state deep inside the conducting pore and, thus to play an unexpected but significant role in aquaglyceroporin’s functions.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Molecular mechanisms of estrogen recognition and 17-keto reduction by human 17β-hydroxysteroid dehydrogenase 1

The reduction of inactive estrone (E1) to the active estrogen 17β-estradiol (E2) is catalyzed by type 1 17β-hydroxysteroid dehydrogenase (17HSD1). Crystallographic studies, modeling and activity measurement of mutants and chimeric enzymes have led to the understanding of its mechanism of action and the molecular basis for the estrogenic specificity. An electrophilic attack on the C17-keto oxygen by the Tyr 155 hydroxyl is proposed for initiation of the transition state. The active site is a hydrophobic pocket with catalytic residues at one end and the recognition machinery on the other. Tyr 155, Lys 159 and Ser 142 are essential for the activity. The presence of certain other amino acids near the substrate recognition end of the active site including His 152 and Pro 187 is critical to the shape complementarity of estrogenic ligands. His 221 and Glu 282 form hydrogen bonds with 3-hydroxyl of the aromatic A-ring of the ligand. This mechanism of recognition of E1 by 17HSD1 is similar to that of E2 by estrogen receptor α. In a ternary complex with NADP⁺ and equilin, an equine estrogen with C7=C8 double bond, the orientation of C17=O of equilin relative to the C4-hydride is more acute than the near normal approach of the hydride for the substrate. In the apo-enzyme structure, a substrate-entry loop (residues 186–201) is in an open conformation. The loop is closed in this complex and Phe 192 and Met 193 make contacts with the ligand. Residues of the entry loop could be partially responsible for the estrogenic specificity.     

Side-chain dynamics are critical for water permeation through aquaporin-1

Molecular dynamics simulations of aquaporin-1 embedded in a solvated lipid bilayer were carried out to investigate the mechanism of water permeation. The 2.2 Å resolution crystal structure of the bovine protein was used for five independent trajectories. During the equilibration and preparatory steps in which the protein was held fixed, water molecules inside the water channel adopted the same positions as observed in the crystal structure but they did not pass through the channel, suggesting that the dynamic motion of the protein is critical for water permeation. When the protein atoms were allowed to move, the side chains of the two asparagines in the two conserved Asn-Pro-Ala motifs near the center of the channel formed hydrogen bonds with water and helped water molecules move through the channel by actively aligning them for transport. The main-chain oxygen atoms, which were exposed to the pore surface in the crystal structure, also contributed to water transfer. Besides the constriction region observed in the crystal structure (Arg¹⁹⁷, Phe⁵⁸, His¹⁸², and Cys¹⁹¹), we found that His⁷⁶ and Val¹⁵⁵ act as a valve by dynamically blocking water permeation and helping control flow.     

Two Classes of Myosin Inhibitors,Para-nitroblebbistatin and Mavacamten,Stabilize β-Cardiac Myosin in Different Structural and Functional States

In addition to a conventional relaxed state, a fraction of myosins in the cardiac muscle exists in a low-energy consuming super-relaxed (SRX) state, which is kept as a reserve pool that may be engaged under sustained increased cardiac demand. The conventional relaxed and the super-relaxed states are widely assumed to correspond to a structure where myosin heads are in an open configuration, free to interact with actin, and a closed configuration, inhibiting binding to actin, respectively. Disruption of the myosin SRX population is an emerging model in different heart diseases, such as hypertrophic cardiomyopathy, which results in excessive muscle contraction, and stabilizing them using myosin inhibitors is budding as an attractive therapeutic strategy. Here we examined the structure–function relationships of two myosin ATPase inhibitors, mavacamten and para-nitroblebbistatin, and found that binding of mavacamten at a site different than para-nitroblebbistatin populates myosin into the SRX state. Para-nitroblebbistatin, binding to a distal pocket to the myosin lever arm near the nucleotide-binding site, does not affect the usual myosin SRX state but instead appears to render myosin into a new, perhaps much more inhibited, ‘ultra-relaxed’ state. X-ray scattering-based rigid body modeling shows that both mavacamten and para-nitroblebbistatin induce novel conformations in human β-cardiac heavy meromyosin that diverge significantly from the hypothetical open and closed states, and furthermore, mavacamten treatment causes greater compaction than para-nitroblebbistatin. Taken together, we conclude that mavacamten and para-nitroblebbistatin stabilize myosin in different structural states, and such states may give rise to different functional energy-sparing states.     

Reconstitution of water channel function and 2D-crystallization of human aquaporin 8

Among the thirteen human aquaporins (AQP0-12), the primary structure of AQP8 is unique. By sequence alignment it is evident that mammalian AQP8s form a separate subfamily distinct from the other mammalian aquaporins. The constriction region of the pore determining the solute specificity deviates in AQP8 making it permeable to both ammonia and H(2)O(2) in addition to water. To better understand the selectivity and gating mechanism of aquaporins, high-resolution structures are necessary. So far, the structure of three human aquaporins (HsAQP1, HsAQP4, and HsAQP5) have been solved at atomic resolution. For mammalian aquaporins in general, high-resolution structures are only available for those belonging to the water-specific subfamily (including HsAQP1, HsAQP4 and HsAQP5). Thus, it is of interest to solve structures of other aquaporin subfamily members with different solute specificities. To achieve this the aquaporins need to be overexpressed heterologously and purified. Here we use the methylotrophic yeast Pichia pastoris as a host for the overexpression. A wide screen of different detergents and detergent-lipid combinations resulted in the solubilization of functional human AQP8 protein and in well-ordered 2D crystals. It also became evident that removal of amino acids constituting affinity tags was crucial to achieve highly ordered 2D crystals diffracting to 3?.     

Mercury inhibits the L170C mutant of aquaporin Z by making waters clog the water channel

We conduct in silico experiments of the L170C mutant of the Escherichia coli aquaporin Z (AQPZ) with and without mercury bonded to residue Cys 170. We find that bonding mercury to Cys 170 does not induce consequential structural changes to the protein. We further find that mercury does not stick in the middle of the water channel to simply occlude water permeation, but resides on the wall of the water pore. However, we observe that the water permeation coefficient of L170C-Hg⁺ (with one mercury ion bonded to Cys 170) is approximately half of that of the mercury-free L170C. We examine the interactions between the mercury ion and the waters in its vicinity and find that five to six waters are strongly attracted by the mercury ion, occluding the space of the water channel. Therefore we conclude that mercury, at low concentration, inhibits AQPZ-L170C mutant by making water molecules clog the water channel.     

QM/MM computations reveal details of the acetyl-CoA synthase catalytic center

The “open” (A_open) and “closed” (A_closed) A-clusters of the acteyl-CoA synthase (ACS) enzyme from Moorella thermoacetica have been studied using a combined quantum mechanical (QM)/molecular mechanical (MM) approach. Geometry optimizations of the oxidized, one- and two-electron reduced A_open state have been carried out for the fully solvated ACS enzyme, and the CO ligand has been modeled in the reduced models. Using a combination of both α_open and α_closed protein scaffolds and the positions of metal atoms in these structures, we have been able to piece together critical parts of the catalytic cycle of ACS. We have replaced the unidentified exogenous ligand in the crystal structure with CO using both a square planar and tetrahedral proximal Ni atom. A one-electron reduced A-cluster that is characterized by a proximal Ni atom in a tetrahedral coordination pattern observed in both the A_open (lower occupancy proximal Ni) and A_closed (proximal Zn atom) geometries with three cysteine thiolates and a modeled CO ligand demonstrates excellent agreement with the crystal structure atomic positions, particularly with the displacement of the side chain ring of Phe512 which appears to serve as a structural gate for ligand binding. The QM/MM optimized geometry of the A-cluster of ACS with an uncoordinated, oxidized proximal nickel atom in a square planar geometry demonstrates poor agreement with the atomic coordinates taken from the crystal structure. Based on these calculations, we conclude that the square planar proximal nickel coordination that has been captured in the A_open structure does not correspond to the ligand-free, oxidized [Fe₄S₄]²⁺ − Ni_p²⁺ − Ni_d²⁺ state. Overall, these computations shed further light on the mechanistic details of protein conformational changes and electronic transitions involved in the ACS catalytic cycle.     

Crystal Structure of an Ammonia-Permeable Aquaporin

Aquaporins of the TIP subfamily (Tonoplast Intrinsic Proteins) have been suggested to facilitate permeation of water and ammonia across the vacuolar membrane of plants, allowing the vacuole to efficiently sequester ammonium ions and counteract cytosolic fluctuations of ammonia. Here, we report the structure determined at 1.18 Å resolution from twinned crystals of Arabidopsis thaliana aquaporin AtTIP2;1 and confirm water and ammonia permeability of the purified protein reconstituted in proteoliposomes as further substantiated by molecular dynamics simulations. The structure of AtTIP2;1 reveals an extended selectivity filter with the conserved arginine of the filter adopting a unique unpredicted position. The relatively wide pore and the polar nature of the selectivity filter clarify the ammonia permeability. By mutational studies, we show that the identified determinants in the extended selectivity filter region are sufficient to convert a strictly water-specific human aquaporin into an AtTIP2;1-like ammonia channel. A flexible histidine and a novel water-filled side pore are speculated to deprotonate ammonium ions, thereby possibly increasing permeation of ammonia. The molecular understanding of how aquaporins facilitate ammonia flux across membranes could potentially be used to modulate ammonia losses over the plasma membrane to the atmosphere, e.g., during photorespiration, and thereby to modify the nitrogen use efficiency of plants.     

Location of a common inhibitor binding site in the cytoplasmic vestibule of the cystic fibrosis transmembrane conductance regulator chloride channel pore

Chloride transport by the cystic fibrosis transmembrane conductance regulator (CFTR) Cl- channel is inhibited by a broad range of organic anions that enter the channel pore from its cytoplasmic end, physically occluding the Cl- permeation pathway. These open channel blocker molecules are presumed to bind within a relatively wide pore inner vestibule that shows little discrimination between different large anions. The present study uses patch clamp recording to identify a pore-lining lysine residue, Lys-95, that acts to attract large blocker molecules into this inner vestibule. Mutations that remove the fixed positive charge associated with this amino acid residue dramatically weaken the blocking effects of five structurally unrelated open channel blockers (glibenclamide, 4,4'-dinitrostilbene-2,2'-disulfonic acid, lonidamine, 5-nitro-2-(3-phenylpropylamino)benzoic acid, and taurolithocholate-3-sulfate) when applied to the cytoplasmic face of the membrane. Mutagenesis of Lys-95 also induced amino acid side chain charge-dependent rectification of the macroscopic current-voltage relationship, consistent with the fixed positive charge on this residue normally acting to attract Cl- ions from the intracellular solution into the pore. These results identify Lys-95 as playing an important role in attracting permeant anions into the channel pore inner vestibule, probably by an electrostatic mechanism. This same electrostatic attraction mechanism also acts to attract larger anionic molecules into the relatively wide inner vestibule, where these substances bind to block Cl- permeation. Thus, structurally diverse open channel blockers of CFTR appear to share a common molecular mechanism of action that involves interaction with a positively charged amino acid side chain located in the inner vestibule of the pore.     

Model of multifractal gating of single ionic channels in biological membranes

A physico-mathematical model of the gating machinery of single ionic channels in biological membranes has been developed. In the paradigm of this model, gating particles are subjected to: (i) deterministic friction force responsible for interactions of gating particles with the surrounding solution; (ii) deterministic potential force depending on the structure and conformational state of the channel pore (the latter is controlled by the transmembrane voltage V and regulates the motion of particles overcoming potential barriers on going from the closed (open) to the open (closed) state of the channel); (iii) deterministic force responsible for interactions of water molecules with hydrophobic sites in the channel pore, and, finally, (iv) stochastic thermal fluctuation force. The model affords adequate approximation of experimental data.     

Dynamics and energetics of water permeation through the aquaporin channel

Structural properties of water inside bovine aquaporin-1 are investigated by molecular simulation. The calculations, which are based on the recently determined X-ray structure at 2.2 A resolution (Sui et al., Nature 2001;414:872-878), are carried out on one monomeric subunit immersed in a water-n-octane-water bilayer. Molecular dynamics (MD) simulations suggest that His182, a fully conserved residue in the channel pore, is protonated in the delta position. Furthermore, they reveal a highly ordered water structure in the channel, induced by the electrostatic properties of the protein. Multiple-steering MD simulations are used to calculate the free-energy of water diffusion. To the best of our knowledge, this represents the first free-energy calculation based on the new, high-resolution structure of the pore. The calculated barrier is 2.5 kcal/mol, and it is associated to water permeation through the Asn-Pro-Ala (NPA) region of the pore, where water molecules are only hydrogen-bonded with themselves. These findings are fully consistent with those based on the previous MD studies on the human protein (de Groot and Grubmüller, Science 2001;294:2353-2357).     

Interaction between the Pore and a Fast Gate of the Cardiac Sodium Channel

Permeant ions affect a fast gating process observed in human cardiac sodium channels (Townsend, C., H.A. Hartmann, and R. Horn. 1997. J. Gen. Physiol. 110:11–21). Removal of extracellular permeant ions causes a reduction of open probability at positive membrane potentials. These results suggest an intimate relationship between the ion-conducting pore and the gates of the channel. We tested this hypothesis by three sets of manipulations designed to affect the binding of cations within the pore: application of intracellular pore blockers, mutagenesis of residues known to contribute to permeation, and chemical modification of a native cysteine residue (C373) near the extracellular mouth of the pore. The coupling between extracellular permeant ions and this fast gating process is abolished both by pore blockers and by a mutation that severely affects selectivity. A more superficial pore mutation or chemical modification of C373 reduces single channel conductance while preserving both selectivity of the pore and the modulatory effects of extracellular cations. Our results demonstrate a modulatory gating role for a region deep within the pore and suggest that the structure of the permeation pathway is largely preserved when a channel is closed.     

Molecular Mechanisms of How Mercury Inhibits Water Permeation through Aquaporin-1: Understanding by Molecular Dynamics Simulation

Aquaporin (AQP) functions as a water-conducting pore. Mercury inhibits the water permeation through AQP. Although site-directed mutagenesis has shown that mercury binds to Cys¹⁸⁹ during the inhibition process, it is not fully understood how this inhibits the water permeation through AQP1. We carried out 40 ns molecular dynamics simulations of bovine AQP1 tetramer with mercury (Hg-AQP1) or without mercury (Free AQP1). In Hg-AQP1, Cys¹⁹¹ (Cys¹⁸⁹ in human AQP1) is converted to Cys-SHg⁺ in each monomer. During each last 10 ns, we observed water permeation events occurred 23 times in Free AQP1 and never in Hg-AQP1. Mercury binding did not influence the whole structure, but did induce a collapse in the orientation of several residues at the ar/R region. In Free AQP1, backbone oxygen atoms of Gly¹⁹⁰, Cys¹⁹¹, and Gly¹⁹² lined, and were oriented to, the surface of the water pore channel. In Hg-AQP1, however, the SHg⁺ of Cys191-SHg⁺ was oriented toward the outside of the water pore. As a result, the backbone oxygen atoms of Gly¹⁹⁰, Cys¹⁹¹, and Gly¹⁹² became disorganized and the ar/R region collapsed, thereby obstructing the permeation of water. We suggest that mercury disrupts the water pore of AQP1 through local conformational changes in the ar/R region.     

Determinants of Water Permeability through Nanoscopic Hydrophilic Channels

Naturally occurring pores show a variety of polarities and sizes that are presumably directly linked to their biological function. Many biological channels are selective toward permeants similar or smaller in size than water molecules, and therefore their pores operate in the regime of single-file water pores. Intrinsic factors affecting water permeability through such pores include the channel-membrane match, the structural stability of the channel, the channel geometry and channel-water affinity. We present an extensive molecular dynamics study on the role of the channel geometry and polarity on the water osmotic and diffusive permeability coefficients. We show that the polarity of the naturally occurring peptidic channels is close to optimal for water permeation, and that the water mobility for a wide range of channel polarities is essentially length independent. By systematically varying the geometry and polarity of model hydrophilic pores, based on the fold of gramicidin A, the water density, occupancy, and permeability are studied. Our focus is on the characterization of the transition between different permeation regimes in terms of the structure of water in the pores, the average pore occupancy and the dynamics of the permeating water molecules. We show that a general relationship between osmotic and diffusive water permeability coefficients in the single-file regime accounts for the time averaged pore occupancy, and that the dynamics of the permeating water molecules through narrow non single file channels effectively behaves like independent single-file columns.     

Computer simulation of water in cytochrome c oxidase

Statistical mechanics and molecular dynamics simulations have been carried out to study the distribution and dynamics of internal water molecules in bovine heart cytochrome c oxidase (CcO). CcO is found to be capable of holding plenty of water, which in subunit I alone amounts to about 165 molecules. The dynamic characterization of these water molecules is carried out. The nascent water molecules produced in the redox reaction at the heme a(3)-CuB binuclear site form an intriguing chain structure. The chain begins at the position of Glu242 at the end of the D channel, and has a fork structure, one branch of which leads to the binuclear center, and the other to the propionate d of heme a(3). The branch that leads to the binuclear center has dynamic access both to the site where the formation of water occurs, and to delta-nitrogen of His291. From the binuclear center, the chain continues to run into the K channel. The stability of this hydrogen bond network is examined dynamically. The catalytic site is located at the hydrophobic region, and the nascent water molecules are produced at the top of the energy hill. The energy gradient is utilized as the mechanism of water removal from the protein. The water exit channels are explored using high-temperature dynamics simulations. Two putative channels for water exit from the catalytic site have been identified. One is leading directly toward Mg(2+) site. However, this channel is only open when His291 is dissociated from CuB. If His291 is bound to CuB, the only channel for water exit is the one that originates at E242 and leads toward the middle of the membrane. This is the same channel that is presumably used for oxygen supply.     

The Outer Pore and Selectivity Filter of TRPA1

TRPA1 (transient-receptor-potential-related ion channel with ankyrin domains) is a direct receptor or indirect effector for a wide variety of nociceptive signals, and thus is a compelling target for development of analgesic pharmaceuticals such as channel blockers. Recently, the structure of TRPA1 was reported, providing insights into channel assembly and pore architecture. Here we report whole-cell and single-channel current recordings of wild-type human TRPA1 as well as TRPA1 bearing point mutations of key charged residues in the outer pore. These measurements demonstrate that the glutamate at position 920 plays an important role in collecting cations into the mouth of the pore, by changing the effective surface potential by ~16 mV, while acidic residues further out have little effect on permeation. Electrophysiology experiments also confirm that the aspartate residue at position 915 represents a constriction site of the TRPA1 pore and is critical in controlling ion permeation.