Water transporters: how so fast yet so selective?

High-Performance Biomimetic Water Channel: The Constructive Interplay of Interaction Parameters and Hydrophilic Doping Levels

In this study, we introduce a superfast biomimetic water channel mimicking the hydrophobicity scales of the Aquaporin (AQP) pore lining. Molecular dynamics simulation is used to scrutinize the impact of hydrophilic doping level in the nanotube and the water-wall interaction strength on water permeability. In the designed biomimetic channel, the constructive interplay of Lennard-Jones (LJ) ε parameters and hydrophilic doping levels increased the possibility of ultrafast water transport. Moreover, a unique set of LJ parameters is discovered for each biomimetic channel with different hydrophilic doping levels, enhancing water permeation. Inside high-performance biomimetic channels, water distribution surprisingly implies a varying pore geometry that narrows down in the middle, mimicking the pattern obtained from GplF pore analysis, evoking the narrow pore induced by the aromatic/arginine selectivity filter. This exciting accordance occurred as a result of tailoring specific hydrophilic arrays within the hydrophobic channel backbone by mimicking the AQP pore interior. The main takeaway of hydrophilic doping arrays implanted within the hydrophobic nanotube is to break the large barrier in the water-wall vdW energy profile into multiple reduced ones to increase water conduction. Consequently, the "water jumping" phenomenon in the middle of the biomimetic channel occurs under specific circumstances. The biomimetic channel with the highest value of water permeability of about 13.67 ± 0.66 × 10^-13 cm³·s^-1 exhibits the best mechanism for artificial water channels (AWCs), serving superfast water transport considering the low entrance barrier and weak water-wall interaction.

Topological analysis of Hedgehog acyltransferase, a multipalmitoylated transmembrane protein

Hedgehog proteins are secreted morphogens that play critical roles in development and disease. During maturation of the proteins through the secretory pathway, they are modified by the addition of N-terminal palmitic acid and C-terminal cholesterol moieties, both of which are critical for their correct function and localization. Hedgehog acyltransferase (HHAT) is the enzyme in the endoplasmic reticulum that palmitoylates Hedgehog proteins, is a member of a small subfamily of membrane-bound O-acyltransferase proteins that acylate secreted proteins, and is an important drug target in cancer. However, little is known about HHAT structure and mode of function. We show that HHAT is comprised of ten transmembrane domains and two reentrant loops with the critical His and Asp residues on opposite sides of the endoplasmic reticulum membrane. We further show that HHAT is palmitoylated on multiple cytosolic cysteines that maintain protein structure within the membrane. Finally, we provide evidence that mutation of the conserved His residue in the hypothesized catalytic domain results in a complete loss of HHAT palmitoylation, providing novel insights into how the protein may function in vivo.

Expression and functional analysis of the rice plasma-membrane intrinsic protein gene family

Plasma membrane intrinsic proteins (PIPs) are a subfamily of aquaporins that enable fast and controlled translocation of water across the membrane. In this study, we systematically identified and cloned ten PIP genes from rice. Based on the similarity of the amino acid sequences they encoded, these rice PIP genes were classified into two groups and designated as OsPIP1-1 to OsPIP1-3 and OsPIP2-1 to OsPIP2-7 following the nomenclature of PIP genes in maize. Quantitative RT-PCR analysis identified three root-specific and one leaf-specific OsPIP genes. Furthermore, the expression profile of each OsPIP gene in response to salt, drought and ABA treatment was examined in detail. Analysis on transgenic plants over-expressing of either OsPIP1 (OsPIP1-1) or OsPIP2 (OsPIP2-2) in wild-type Arabidopsis, showed enhanced tolerance to salt (100 mM of NaCl) and drought (200 mM of mannitol), but not to salt treatment of higher concentration (150 mM of NaCl). Taken together, these data suggest a distinct role of each OsPIP gene in response to different stresses, and should add a new layer to the understanding of the physiological function of rice PIP genes.

Acquired cadmium resistance in metallothionein-I/II(-/-) knockout cells: role of the T-type calcium channel Cacnalpha1G in cadmium uptake

Metallothioneins (MTs) are cytoplasmic proteins that sequester certain divalent cations and are considered a primary cellular defense against the toxic transition metal cadmium (Cd(2+)). MT-I/II(-/-) knockout [MT(-/-)] cells are available and serve as an excellent tool to study non-MT-related mechanisms in metal tolerance. In the current study, Cd(2+)-resistant MT(-/-) (CdR) and CdR revertant (CdR-rev) cell lines were developed and characterized to investigate non-MT-mediated cellular protection mechanisms. Resistance to Cd(2+) was approximately 70-fold higher in CdR than the parental MT(-/-) cell line (IC(50) = 20 versus 0.3 microM, respectively) and was stable in the absence of Cd(2+) for 35 days. Accumulation of Cd(2+) by the CdR cell line was reduced by approximately 95% compared with parental cells, primarily because of a decreased Cd(2+) uptake. Cd(2+) uptake by the MT(-/-) parental cell line was independent of sodium, energy, and electrogenic potential. Uptake was saturable (K(m) = 65 nM; V(max) = 4.9 pmol/mg/min) and pH-dependent (maximal at pH 6.5-7). Potent inhibitors of Cd(2+) uptake included Zn(2+) (IC(50) = 7 microM), Mn(2+) (IC(50) = 0.4 microM), and the T-type Ca(2+) channel antagonist mibefradil (IC(50) = 5 microM), whereas other metals (including Fe(2+)) and L-type Ca(2+) channel antagonists had little effect. Immunoblot and real-time reverse transcription-polymerase chain reaction analysis indicated that the Cacnalpha(1G) T-type Ca(2+) channel was expressed at a reduced level in CdR compared with the parental MT(-/-) cell line, suggesting it is important for Cd(2+) uptake. The CdR1-rev cell line was found to have a Cd(2+) uptake and sensitivity level in between that of the CdR1 and MT(-/-) cell lines. Consistent with this was an intermediate expression of Cacnalpha(1G) in the CdR-rev cell line. These data suggest that decreased expression of Cacnalpha(1G) protects cells from Cd(2+) exposure by limiting Cd(2+) uptake.

Water Cluster Growth in Hydrophobic Solid Nanospaces

The growth mechanism of water clusters in carbon nanopores is clearly elucidated by in situ small-angle X-ray scattering (SAXS) studies and grand canonical Monte Carlo (GCMC) simulations at 293-313 K. Water molecules are isolated from each other in hydrophobic nanopores below relative pressures (P/P(0)) of 0.5. Water molecules associate with each other to form clusters of about 0.6 nm in size at P/P(0)=0.6, accompanied by a remarkable aggregation of these clusters. The complete filling of carbon nanopores finishes at about P/P(0)=0.8. The correlation length analysis of SAXS profiles leads to the proposal of a growth mechanism for these water clusters and the presence of the critical cluster size of 0.6 nm leads to extremely stable clusters of water molecules in hydrophobic nanopores. Once a cluster of the critical size is formed in hydrophobic nanopores, the predominant water adsorption begins to fill carbon nanopores.

Studies of proton translocations in biological systems: simulating proton transport in carbonic anhydrase by EVB-based models

Proton transport (PTR) processes play a major role in bioenergetics and thus it is important to gain a molecular understanding of these processes. At present the detailed description of PTR in proteins is somewhat unclear and it is important to examine different models by using well-defined experimental systems. One of the best benchmarks is provided by carbonic anhydrase III (CA III), because this is one of the few systems where we have a clear molecular knowledge of the rate constant of the PTR process and its variation upon mutations. Furthermore, this system transfers a proton between several water molecules, thus making it highly relevant to a careful examination of the "proton wire" concept. Obtaining a correlation between the structure of this protein and the rate of the PTR process should help to discriminate between alternative models and to give useful clues about PTR processes in other systems. Obviously, obtaining such a correlation requires a correct representation of the "chemistry" of PTR between different donors and acceptors, as well as the ability to evaluate the free energy barriers of charge transfer in proteins, and to simulate long-time kinetic processes. The microscopic empirical valence bond (Warshel, A., and R. M. Weiss. 1980. J. Am. Chem. Soc. 102:6218-6226; and Aqvist, J., and A. Warshel. 1993. Chem. Rev. 93:2523-2544) provides a powerful way for representing the chemistry and evaluating the free energy barriers, but it cannot be used with the currently available computer times in direct simulation of PTR with significant activation barriers. Alternatively, one can reduce the empirical valence bond (EVB) to the modified Marcus' relationship and use semimacroscopic electrostatic calculations plus a master equation to determine the PTR kinetics (Sham, Y., I. Muegge, and A. Warshel. 1999. Proteins. 36:484-500). However, such an approximation does not provide a rigorous multisite kinetic treatment. Here we combine the useful ingredients of both approaches and develop a simplified EVB effective potential that treats explicitly the chain of donors and acceptors while considering implicitly the rest of the protein/solvent system. This approach can be used in Langevin dynamics simulations of long-time PTR processes. The validity of our new simplified approach is demonstrated first by comparing its Langevin dynamics results for a PTR along a chain of water molecules in water to the corresponding molecular dynamics simulations of the fully microscopic EVB model. This study examines dynamics of both models in cases of low activation barriers and the dependence of the rate on the energetics for cases with moderate barriers. The study of the dependence on the activation barrier is next extended to the range of higher barriers, demonstrating a clear correlation between the barrier height and the rate constant. The simplified EVB model is then examined in studies of the PTR in carbonic anhydrase III, where it reproduces the relevant experimental results without the use of any parameter that is specifically adjusted to fit the energetics or dynamics of the reaction in the protein. We also validate the conclusions obtained previously from the EVB-based modified Marcus' relationship. It is concluded that this approach and the EVB-based model provide a reliable, effective, and general tool for studies of PTR in proteins. Finally in view of the behavior of the simulated result, in both water and the CA III, we conclude that the rate of PTR in proteins is determined by the electrostatic energy of the transferred proton as long as this energy is higher than a few kcal/mol.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Realistic Simulations of Proton Transport along the Gramicidin Channel: Demonstrating the Importance of Solvation Effects

The nature of proton transduction (PTR) through a file of water molecules, along the gramicidin A (gA) channel, has long been considered as being highly relevant to PTR in biological systems. Previous attempts to model this process implied that the so-called Grotthuss mechanism and the corresponding orientation of the water file plays a major role. The present work reexamines the PTR in gA by combining a fully microscopic empirical valence bond (EVB) model and a recently developed simplified EVB-based model with Langevin dynamics (LD) simulations. The full model is used first to evaluate the free energy profile for a stepwise PTR process. The corresponding results are then used to construct the effective potential of the simplified EVB. This later model is then used in Langevin dynamics simulations, taking into account the correct physics of possible concerted motions and the effect of the solvent reorganization. The simulations reproduce the observed experimental trend and lead to a picture that is quite different from that assumed previously. It is found that the PTR in gA is controlled by the change in solvation energy of the transferred proton along the channel axis. Although the time dependent electrostatic fluctuations of the channel and water dipoles play their usual role in modulating the proton-transfer process (Proc. Natl. Acad. Sci. U.S.A. 1984, 81, 444), the PTR rate is mainly determined by the free energy profile. Furthermore, the energetics of the reorientation of the unprotonated water file do not appear to provide a consistent way of assessing the activation barrier for the PTR process. It seems to us that in the case of gA, and probably other systems with significant electrostatic barriers for the transfer of the proton charge, the PTR rate is controlled by the electrostatic barrier. This finding has clear consequences with regards to PTR processes in biological systems.

What really prevents proton transport through aquaporin? Charge self-energy versus proton wire proposals

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.

Glycerol conductance and physical asymmetry of the Escherichia coli glycerol facilitator GlpF

The aquaglyceroporin GlpF is a transmembrane channel of Escherichia coli that facilitates the uptake of glycerol by the cell. Its high glycerol uptake rate is crucial for the cell to survive in very low glycerol concentrations. Although GlpF allows both influx and outflux of glycerol, its structure, similar to the structure of maltoporin, exhibits a significant degree of asymmetry. The potential of mean force characterizing glycerol in the channel shows a corresponding asymmetry with an attractive vestibule only at the periplasmic side. In this study, we analyze the potential of mean force, showing that a simplified six-step model captures the kinetics and yields a glycerol conduction rate that agrees well with observation. The vestibule improves the conduction rate by 40% and 75% at 10- micro M and 10-mM periplasmic glycerol concentrations, respectively. In addition, neither the conduction rate nor the conduction probability for a single glycerol (efficiency) depends on the orientation of GlpF. GlpF appears to conduct equally well in both directions under physiological conditions.

Mechanisms of selectivity in channels and enzymes studied with interactive molecular dynamics

Interactive molecular dynamics, a new modeling tool for rapid investigation of the physical mechanisms of biological processes at the atomic level, is applied to study selectivity and regulation of the membrane channel protein GlpF and the enzyme glycerol kinase. These proteins facilitate the first two steps of Escherichia coli glycerol metabolism. Despite their different function and architecture the proteins are found to employ common mechanisms for substrate selectivity: an induced geometrical fit by structurally homologous binding sites and an induced rapid dipole moment reversal. Competition for hydrogen bonding sites with water in both proteins is critical for substrate motion. In glycerol kinase, it is shown that the proposed domain motion prevents competition with water, in turn regulating the binding of glycerol.