Free-energy landscape of glycerol permeation through aquaglyceroporin GlpF determined from steered molecular dynamics simulations

The free-energy landscape of glycerol permeation through the aquaglyceroporin GlpF has been estimated in the literature by the nonequilibrium method of steered molecular dynamics (SMD) simulations and by the equilibrium method of adaptive biasing force (ABF) simulations. However, the ABF results qualitatively disagree with the SMD results that were based on the Jarzynski equality (JE) relating the equilibrium free-energy difference to the nonequilibrium work of the irreversible pulling experiments. In this paper, I present a new SMD study of the glycerol permeation through GlpF to explore the free-energy profile of glycerol along the permeation channel. Instead of the JE in terms of thermodynamic work, I use the fluctuation-dissipation theorem (FDT) of Brownian dynamics (BD), in terms of mechanical work, for extracting the free-energy difference from the nonequilibrium work of irreversible pulling experiments. The results of this new SMD-BD-FDT study are in agreement with the experimental data and with the ABF results.     

The structural pathway for water permeation through sodium-glucose cotransporters

Although water permeation across cell membranes occurs through several types of membrane proteins, the only permeation mechanism resolved at atomic scale is that through aquaporins. Crystallization of the Vibrio parahaemolyticus sodium-galactose transporter (vSGLT) allows investigation of putative water permeation pathways through both vSGLT and the homologous human Na-glucose cotransporter (hSGLT1) using computational methods. Grand canonical Monte Carlo and molecular dynamics simulations were used to stably insert water molecules in both proteins, showing the presence of a water-filled pathway composed of ∼100 water molecules. This provides a structural basis for passive water permeation that is difficult to reconcile with the water cotransport hypothesis. Potential-of-mean-force calculations of water going through the crystal structure of vSGLT shows a single barrier of 7.7 kCal mol⁻¹, in agreement with previously published experimental data for cotransporters of the SGLT family. Electrophysiological and volumetric experiments performed on hSGLT1-expressing Xenopus oocytes showed that the passive permeation pathway exists in different conformational states. In particular, experimental conditions that aim to mimic the conformation of the crystal structure displayed passive water permeability. These results provide groundwork for understanding the structural basis of cotransporter water permeability.     

Diffusion of glycerol through Escherichia coli aquaglyceroporin GlpF

The glycerol uptake facilitator, GlpF, a major intrinsic protein found in Escherichia coli, selectively conducts water and glycerol across the inner membrane. The free energy landscape characterizing the assisted transport of glycerol by this homotetrameric aquaglyceroporin has been explored by means of equilibrium molecular dynamics over a timescale spanning 0.12 μs. To overcome the free energy barriers of the conduction pathway, an adaptive biasing force is applied to the glycerol molecule confined in each of the four channels. The results illuminate the critical role played by intramolecular relaxation on the diffusion properties of the permeant. These free energy calculations reveal that glycerol tumbles and isomerizes on a timescale comparable to that spanned by its adaptive-biasing-force-assisted conduction in GlpF. As a result, reorientation and conformational equilibrium of glycerol in GlpF constitute a bottleneck in the molecular simulations of the permeation event. A profile characterizing the position-dependent diffusion of the permeant has been determined, allowing reaction rate theory to be applied for investigating conduction kinetics based on the measured free energy landscape.     

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.     

Molecular dynamics simulation of water permeation through the alpha-hemolysin channel

The alpha-hemolysin (AHL) nanochannel is a non-selective channel that allows for uncontrolled transport of small molecules across membranes leading to cell death. Although it is a bacterial toxin, it has promising applications, ranging from drug delivery systems to nano-sensing devices. This study focuses on the transport of water molecules through an AHL nanochannel using molecular dynamics (MD) simulations. Our results show that AHL can quickly transport water across membranes. The first-passage time approach was used to estimate the diffusion coefficient and the mean exit time. To study the energetics of transport, the potential of mean force (PMF) of a water molecule along the AHL nanochannel was calculated. The results show that the energy barriers of water permeation across a nanopore are always positive along the channel and the values are close to thermal energy (k_BT). These findings suggest that the observed quick permeation of water is due to small energy barriers and a hydrophobic inner channel surface resulting in smaller friction. We speculate that these physical mechanisms are important in how AHL causes cell death.     

Dihydroxyacetone and methylglyoxal as permeants of the Plasmodium aquaglyceroporin inhibit parasite proliferation

The aquaglyceroporin of Plasmodium falciparum (PfAQP) is a bi-functional channel with permeability for water and solutes. Its functions supposedly are in osmotic protection of parasites and in facilitation of glycerol permeation for glycerolipid biosynthesis. Here, we show PfAQP permeability for the glycolysis-related metabolites methylglyoxal, a cytotoxic byproduct, and dihydroxyacetone, a ketotriose. AQP3, the red cell aquaglyceroporin, also passed dihydroxacetone but excluded methylglyoxal. Proliferation of malaria parasites was inhibited by methylglyoxal with an IC50 around 200 microM. Surprisingly, also dihydroxyacetone, which is an energy source in human cells, was antiproliferative in chloroquine-sensitive and resistant strains with an IC50 around 3 mM. We expressed P. falciparum glyceraldehyde 3-phosphate dehydrogenase (PfGAPDH) to examine whether it is inhibited by either carbonyl compound. Methylglyoxal did not affect PfGAPDH on incubation with 2.5 mM for 20 h. Treatment with 2.5 mM dihydroxyacetone, however, abolished PfGAPDH activity within 6 h. Aquaglyceroporin permeability for glycolytic metabolites may thus be of physiological significance.     

A single, bi-functional aquaglyceroporin in blood-stage Plasmodium falciparum malaria parasites.

The malaria parasite Plasmodium falciparum faces drastic osmotic changes during kidney passages and is engaged in the massive biosynthesis of glycerolipids during its development in the blood-stage. We identified a single aquaglyceroporin (PfAQP) in the nearly finished genome of P. falciparum with highest similarity to the Escherichia coli glycerol facilitator (50.4%), but both canonical Asn-Pro-Ala (NPA) motifs in the pore region are changed to Asn-Leu-Ala (NLA) and Asn-Pro-Ser (NPS), respectively. Expression in Xenopus oocytes renders them highly permeable for both water and glycerol. Sugar alcohols up to five carbons and urea pass the pore. Mutation analyses of the NLA/NPS motifs showed their structural importance, but the symmetrical pore properties were maintained. PfAQP is expressed in blood-stage parasites throughout the development from rings via trophozoites to schizonts and is localized to the parasite but not to the erythrocyte cytoplasm or membrane. Its unique bi-functionality indicates functions in the protection from osmotic stress and efficiently provides access to the serum glycerol pool for the use in ATP generation and primarily in the phospholipid synthesis.     

Rapid evaluation of oxygen and water permeation through microplate sealing tapes

Eight commercially available microplate sealing tapes and 10 other suitable materials (transparent wound dressings) are compared qualitatively in terms of their ability to minimize water evaporation from a multiwell plate while maintaining the oxygen supply as high as possible, which is necessary for applications like aerobic growth. The transparency and sterility of the products are considered as well. All evaluated commercially available sealing tapes fall into one of the following two classes: (1) O(2) transfer is comparable to that of an unsealed plate, but water vapor retention is relatively low, or (2) O(2) transfer via the sealing is slower, but the water retention capability is comparably high. All but one of the evaluated wound dressings fall under the second class. That dressing, however, constitutes a compromise by showing both moderate O(2) permeability and medium water retention. But the estimated mass transport in a microtiter plate sealed with this dressing is about 5 times slower than that of an unsealed 96 well plate. The aim of this publication is to enable the reader to choose a microtiter plate sealing from the materials evaluated within this work and to use the rapid methods described herein to easily perform tests of additional sealing materials.     

Novel Structure In the Pellicular Complex of Plasmodium Falciparum Gametocytes1

ABSTRACT. Using transmission electron microscopy, transverse dense bands were found to be associated with subpellicular microtubules and inner membraner in Plasmodium falciparum gametocytes. These dense bands may act as supportive structures to maintain the parallel arrangement of the microtubules, and/or to connect them to the inner membranes.     

Bilirubin inhibits Plasmodium falciparum growth through the generation of reactive oxygen species

Free heme is very toxic because it generates highly reactive hydroxyl radicals ((.)OH) to cause oxidative damage. Detoxification of free heme by the heme oxygenase (HO) system is a very common phenomenon by which free heme is catabolized to form bilirubin as an end product. Interestingly, the malaria parasite, Plasmodium falciparum, lacks an HO system, but it forms hemozoin, mainly to detoxify free heme. Here, we report that bilirubin significantly induces oxidative stress in the parasite as evident from the increased formation of lipid peroxide, decrease in glutathione content, and increased formation of H(2)O(2) and (.)OH. Bilirubin can effectively inhibit hemozoin formation also. Furthermore, results indicate that bilirubin inhibits parasite growth and induces caspase-like protease activity, up-regulates the expression of apoptosis-related protein (Gene ID PFI0450c), and reduces the mitochondrial membrane potential. (.)OH scavengers such as mannitol, as well as the spin trap alpha-phenyl-n-tert-butylnitrone, effectively protect the parasite from bilirubin-induced oxidative stress and growth inhibition. These findings suggest that bilirubin, through the development of oxidative stress, induces P. falciparum cell death and that the malaria parasite lacks an HO system probably to protect itself from bilirubin-induced cell death as a second line of defense.     

Flow cytometry-assisted rapid isolation of recombinant Plasmodium berghei parasites exemplified by functional analysis of aquaglyceroporin

The most critical bottleneck in the generation of recombinant Plasmodium berghei parasites is the mandatory in vivo cloning step following successful genetic manipulation. This study describes a new technique for rapid selection of recombinant P. berghei parasites. The method is based on flow cytometry to isolate isogenic parasite lines and represents a major advance for the field, in that it will speed the generation of recombinant parasites as well as cut down on animal use significantly. High expression of GFP during blood infection, a prerequisite for robust separation of transgenic lines by flow cytometry, was achieved. Isogenic recombinant parasite populations were isolated even in the presence of a 100-fold excess of wild-type (WT) parasites. Aquaglyceroporin (AQP) loss-of-function mutants and parasites expressing a tagged AQP were generated to validate this approach. aqp^? parasites grow normally within the WT phenotypic range during blood infection of NMRI mice. Similarly, colonization of the insect vector and establishment of an infection after mosquito transmission were unaffected, indicating that AQP is dispensable for life cycle progression in vivo under physiological conditions, refuting its use as a suitable drug target. Tagged AQP localized to perinuclear structures and not the parasite plasma membrane. We suggest that flow-cytometric isolation of isogenic parasites overcomes the major roadblock towards a genome-scale repository of mutant and transgenic malaria parasite lines.     

Induction of Oxidant Stress by Iron Available in Advanced Forms of Plasmodium Falciparum

Oxidative stress has been incriminated as a deleterious factor in the development of malaria parasites. Various chemical reductones which can undergo cyclic oxidation and reduction, such as ascorbate have been shown to cause oxidative stress to red blood cells. This, naturally-occurring and redox-active compound, can induce the formation of active oxygen derived species, such as superoxide radicals (.O^?₂), hydrogen peroxide (H₂O₂) and hydroxyl radical (OH.), The formation of the hydroxyl radical, the ultimate deleterious species, is mediated by the redox-active and available transition metals iron and copper in the Haber-Weiss reaction.

During the development of the parasite, hemoglobin is progressively digested and a concurrent release of high levels of iron-containing breakdown products takes place within the red blood cell. Indications for the progressive increase in redox-active iron during the growth of P. falciparum have been recently found in our lab: a) adventitious ascorbatc proved highly detrimental to the parasite when added to the mature forms. In contrast, if the parasitized erythrocytes were in the early phase following invasion, and only low levels of iron-containing structures had been liberated. then the observed effect was a small promotion of parasite development. b) erythrocytes containing mature parasites were more potent than erythrocytes containing ring forms as a source for redox-active iron in the acerbate-driven metal-mediated degradation of DNA. The addition of extracts from parasitized erythrocytes and ascorbate to DNA causcd a dose and time dependent DNA degradation. Non-infected erythrocytes had no effect. These findings could also propose that the parasite-dependent accumulation of redox-active forms of iron within the erythrocytes serve as a biological clock triggering the rupture of the red blood cell membrane at the right moment, when the parasite reaches its maturity.     

Nucleoside permeation in mouse erythrocytes infected with Plasmodium yoelii

In normal mouse erythrocytes, nucleoside permeation was almost completely blocked in the presence of binding site-saturating concentrations of nitrobenzylthioinosine, whereas permeation in erythrocytes infected with the malarial parasite, Plasmodium yoelii, was substantial under these conditions, suggesting the presence of a permeation mechanism of low sensitivity to nitrobenzylthioinosine in the infected cells. Binding sites for nitrobenzylthioinosine were more numerous on infected erythrocytes than on uninfected cells. When mice infected with P. yoelii were treated with combinations of tubercidin and nitrobenzylthioinosine 5'-monophosphate, progression of parasitemia was delayed and survival times were increased.     

Glycerol 3-phosphate inhibits swarming and aggregation of Myxococcus xanthus

We have cloned a gene of Myxococcus xanthus with similarities to the permease for glycerol 3-phosphate (G3P) of other bacteria. Expression of the gene increased significantly during the first hours of starvation. Swarming of the wild-type strain was inhibited and aggregation was delayed by G3P. Conversely, a DeltaglpT strain aggregated even on rich medium. These results indicate that G3P may function to regulate the timing of aggregation in M. xanthus.     

Monocarboxylic acid permeation through lipid bilayer membranes

Summary The membrane permeability coefficients for the homologous monocarboxylic acids, formic through hexanoic, as well as benzoic and salicylic, were determined for egg phosphatidylcholine-decane planar bilayer membranes. The permeabilities of formic, acetic and propionic acid were also determined for solvent-free phosphatidylethanolamine bilayers. Permeability coefficients were calculated from tracer fluxes measured under otherwise symmetrical conditions, and precautions were taken to ensure that the values were not underestimated due to unstirred layer effects. The relation between the nonionic (HA) permeability (P ^m ) and the hexadecane/water partition coefficient (K _p ) was: log ^m =0.90 log K_p+0.87 (correlation coefficient=0.996). Formic acid was excluded from the analysis because its permeability was sixfold higher than predicted by the other acids. The permeabilities for solvent-free membranes were similar to those for decanecontaining membranes. The exceptionally high permeability of formic acid and the high correlation of the other permeabilities to the hexadecane/water partition coefficient is a pattern that conforms with other nonelectrolyte permeabilities through bilayers. Similarly, the mean incremental free energy change per methylene group (G-CH₂-) was –764 cal mol^–1, similar to other homologous solutes in other membrane systems. However, much less negative G values (–120, to –400 cal mol^–1) were previously reported for fatty acids permeating bilayers and biological membranes. These values are due primarily to unstirred layer effects, metabolism and binding to membranes and other cell components.