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

Just passing through: Deploying aquaporins in microbial cell factories

As microbial membranes are naturally impermeable to even the smallest biomolecules, transporter proteins are physiologically essential for normal cell functioning. This makes transporters a key target area for engineering enhanced cell factories. As part of the wider cellular transportome, aquaporins (AQPs) are responsible for transporting small polar solutes, encompassing many compounds which are of great interest for industrial biotechnology, including cell feedstocks, numerous commercially relevant polyols and even weak organic acids. In this review, examples of cell factory engineering by targeting AQPs are presented. These AQP modifications aid in redirecting carbon fluxes and boosting bioconversions either by enhanced feedstock uptake, improved intermediate retention, increasing product export into the media or superior cell viability against stressors with applications in both bacterial and yeast production platforms. Additionally, the future potential for AQP deployment and targeting is discussed, showcasing hurdles and considerations of this strategy as well as recent advances and future directions in the field. By leveraging the natural diversity of AQPs and breakthroughs in channel protein engineering, these transporters are poised to be promising tools capable of enhancing a wide variety of biotechnological processes.

Fundamental limits on the rate of bacterial growth and their influence on proteomic composition

Despite abundant measurements of bacterial growth rate, cell size, and protein content, we lack a rigorous understanding of what sets the scale of these quantities and when protein abundances should (or should not) depend on growth rate. Here, we estimate the basic requirements and physical constraints on steady-state growth by considering key processes in cellular physiology across a collection of Escherichia coli proteomic data covering ≈4,000 proteins and 36 growth rates. Our analysis suggests that cells are predominantly tuned for the task of cell doubling across a continuum of growth rates; specific processes do not limit growth rate or dictate cell size. We present a model of proteomic regulation as a function of nutrient supply that reconciles observed interdependences between protein synthesis, cell size, and growth rate and propose that a theoretical inability to parallelize ribosomal synthesis places a firm limit on the achievable growth rate. A record of this paper's transparent peer review process is included in the supplemental information.

Genomic Investigation into the Virulome, Pathogenicity, Stress Response Factors, Clonal Lineages, and Phylogenetic Relationship of Escherichia coli Strains Isolated from Meat Sources in Ghana

Escherichia coli are among the most common foodborne pathogens associated with infections reported from meat sources. This study investigated the virulome, pathogenicity, stress response factors, clonal lineages, and the phylogenomic relationship of E. coli isolated from different meat sources in Ghana using whole-genome sequencing. Isolates were screened from five meat sources (beef, chevon, guinea fowl, local chicken, and mutton) and five areas (Aboabo, Central market, Nyorni, Victory cinema, and Tishegu) based in the Tamale Metropolis, Ghana. Following microbial identification, the E. coli strains were subjected to whole-genome sequencing. Comparative visualisation analyses showed different DNA synteny of the strains. The isolates consisted of diverse sequence types (STs) with the most common being ST155 (n = 3/14). Based Upon Related Sequence Types (eBURST) analyses of the study sequence types identified four similar clones, five single-locus variants, and two satellite clones (more distantly) with global curated E. coli STs. All the isolates possessed at least one restriction-modification (R-M) and CRISPR defence system. Further analysis revealed conserved stress response mechanisms (detoxification, osmotic, oxidative, and periplasmic stress) in the strains. Estimation of pathogenicity predicted a higher average probability score (P_score ≈ 0.937), supporting their pathogenic potential to humans. Diverse virulence genes that were clonal-specific were identified. Phylogenomic tree analyses coupled with metadata insights depicted the high genetic diversity of the E. coli isolates with no correlation with their meat sources and areas. The findings of this bioinformatic analyses further our understanding of E. coli in meat sources and are broadly relevant to the design of contamination control strategies in meat retail settings in Ghana.

The Ionophores CCCP and Gramicidin but Not Nigericin Inhibit Trypanosoma brucei Aquaglyceroporins at Neutral pH

Human African trypanosomiasis (HAT) is caused by Trypanosoma brucei parasites. The T. brucei aquaglyceroporin isoform 2, TbAQP2, has been linked to the uptake of pentamidine. Negative membrane potentials and transmembrane pH gradients were suggested to promote transport of the dicationic antitrypanosomal drug. Application of ionophores to trypanosomes further hinted at direct inhibition of TbAQP2 by carbonyl cyanide m-chlorophenyl hydrazone (CCCP). Here, we tested for direct effects of three classical ionophores (CCCP, nigericin, gramicidin) on the functionality of TbAQP2 and the related TbAQP3 at conditions that are independent from the membrane potential or a proton gradient. We expressed TbAQP2 and TbAQP3 in yeast, and determined permeability of uncharged glycerol at neutral pH using stopped-flow light scattering. The mobile proton carrier CCCP directly inhibited TbAQP2 glycerol permeability at an IC₅₀ of 2 µM, and TbAQP3 to a much lesser extent (IC₅₀ around 1 mM) likely due to different selectivity filter layouts. Nigericin, another mobile carrier, left both isoforms unaffected. The membrane-integral pore-forming gramicidin evenly inhibited TbAQP2 and TbAQP2 in the double-digit micromolar range. Our data exemplify the need for suitable controls to detect unwanted ionophore side effects even when used at concentrations that are typically recommended to disturb the transmembrane ion distribution.

Increasement of O-acetylhomoserine production in Escherichia coli by modification of glycerol-oxidative pathway coupled with optimization of fermentation

OBJECTIVE

O-acetylhomoserine (OAH) is an important platform chemical to produce high-valuable chemicals. To improve the production of O-acetylhomoserine from glycerol, the glycerol-oxidative pathway was investigated and the optimization of fermentation with crude glycerol was carried out.

RESULTS

The glycerol-uptake system and glycerol-oxidative pathway were modified and O-acetyltransferase from Corynebacterium glutamicum was introduced into the engineered strain to produce O-acetylhomoserine. It was found that overexpression of glycerol 3-phosphate dehydrogenase improved the OAH production to 6.79 and 4.21 g/L from pure and crude glycerol, respectively. And the higher OAH production depending on higher level of transcription of glpD. Two-step statistical approach was employed to optimize the fermentation conditions. The significant effects of glycerol, ammonium chloride and yeast extract were screened applying Plackett-Burman design and were optimized further by employing the Response Surface Methodology. Under optimized conditions, the OAH production was up to 9.42 and 7.01 g/L when pure and crude glycerol were used in shake flask cultivations, respectively.

CONCLUSIONS

The enzymatic step catalyzing the oxidation of glycerol through GlpD was the key step for OAH production, which served the foundation for realization of a consistent OAH production from crude glycerol in the future.

Genome Mining and Comparative Pathogenomic Analysis of An Endemic Methicillin-Resistant Staphylococcus Aureus (MRSA) Clone, ST612-CC8-t1257-SCCmec_IVd(2B), Isolated in South Africa

This study undertook genome mining and comparative genomics to gain genetic insights into the dominance of the methicillin-resistant Staphylococcus aureus (MRSA) endemic clone ST612-CC8-t1257-SCCmec_IVd(2B), obtained from the poultry food chain in South Africa. Functional annotation of the genome revealed a vast array of similar central metabolic, cellular and biochemical networks within the endemic clone crucial for its survival in the microbial community. In-silico analysis of the clone revealed the possession of uniform defense systems, restriction-modification system (type I and IV), accessory gene regulator (type I), arginine catabolic mobile element (type II), and type 1 clustered, regularly interspaced, short palindromic repeat (CRISPR)Cas array (N = 7 ± 1), which offer protection against exogenous attacks. The estimated pathogenic potential predicted a higher probability (average Pscore ≈ 0.927) of the clone being pathogenic to its host. The clone carried a battery of putative virulence determinants whose expression are critical for establishing infection. However, there was a slight difference in their possession of adherence factors (biofilm operon system) and toxins (hemolysins and enterotoxins). Further analysis revealed a conserved environmental tolerance and persistence mechanisms related to stress (oxidative and osmotic), heat shock, sporulation, bacteriocins, and detoxification, which enable it to withstand lethal threats and contribute to its success in diverse ecological niches. Phylogenomic analysis with close sister lineages revealed that the clone was closely related to the MRSA isolate SHV713 from Australia. The results of this bioinformatic analysis provide valuable insights into the biology of this endemic clone.

Glycerol transport through the aquaglyceroporin GlpF: bridging dynamics and kinetics with atomic simulation

The aquaglyceroporin GlpF is a member of the aquaporin family. It selectively conducts small molecules, such as glycerol, across the cell membrane under a concentration gradient of the substrate. Atomistic molecular dynamics (MD) simulation would provide great insight into the substrate transport mechanism of GlpF and membrane channels alike. Ideally, non-equilibrium simulations under various concentration gradients of glycerol are desired to emulate the transportation in cells, but this kind of simulation is difficult due to a complicated system setup and high computational cost. Here, we present a new strategy to extract non-equilibrium kinetic information from equilibrium MD simulation. We first performed long-time (totally 22.5 μs) multi-copy equilibrium MD simulations of glycerol conduction through GlpF. Tens of times the spontaneous permeation of glycerol through GlpF was observed, allowing us to elucidate the detailed mechanism of the stereoselectivity for glycerol. Then we employed Markov state model (MSM) analysis of the MD trajectories to identify the intermediate states during glycerol transport and calculate the inter-state transition rate constants. Based on the results of MSM analysis, we built the kinetic models of glycerol transport and calculated the glycerol fluxes under various concentration gradients by solving the master equations. The results agree well with the experimental measurement at a certain glycerol concentration, and provide holistic information on the glycerol conduction capacity of GlpF. Our work demonstrates that long-time atomistic MD simulations can now bridge the microscopic dynamics and the kinetic description of substance transport through membrane channels, hopefully facilitating the engineering of new selective channels for various molecules.

Transcriptome and Proteome of Fish-Pathogenic Streptococcus agalactiae Are Modulated by Temperature

Streptococcus agalactiae is one of the most important pathogens associated with streptococcosis outbreaks in Nile tilapia farms worldwide. High water temperature (above 27°C) has been described as a predisposing factor for the disease in fish. At low temperatures (below 25°C), fish mortalities are not usually observed in farms. Temperature variation can modulate the expression of genes and proteins involved in metabolism, adaptation, and bacterial pathogenicity, thus increasing or decreasing the ability to infect the host. This study aimed to evaluate the transcriptome and proteome of a fish-pathogenic S. agalactiae strain SA53 subjected to in vitro growth at different temperatures using a microarray and label-free shotgun LC-HDMS^E approach. Biological triplicates of isolates were cultured in BHIT broth at 22 or 32°C for RNA and protein isolation and submitted for transcriptomic and proteomic analyses. In total, 1,730 transcripts were identified in SA53, with 107 genes being differentially expressed between the temperatures evaluated. A higher number of genes related to metabolism, mainly from the phosphotransferase system (PTS) and ATP-binding cassette (ABC) transport system, were upregulated at 32°C. In the proteome analysis, 1,046 proteins were identified in SA53, of which 81 were differentially regulated between 22 and 32°C. Proteins involved in defense mechanisms, lipid transport and metabolism, and nucleotide transport and metabolism were upregulated at 32°C. A higher number of interactions were observed in proteins involved in nucleotide transport and metabolism. We observed a low correlation between the transcriptome and proteome datasets. Our study indicates that the transcriptome and proteome of a fish-adapted S. agalactiae strain are modulated by temperature, particularly showing differential expression of genes/proteins involved in metabolism, virulence factors, and adaptation.

The GlpF residue Trp219 is part of an amino-acid cluster crucial for aquaglyceroporin oligomerization and function

The vestibule loop regions of aquaglyceroporins are involved in accumulation of glycerol inside the channel pore. Even though most loop regions do not show high sequence similarity among aquaglyceroporins, loop E is highly conserved in aquaglyceroporins, but not in members of the homologous aquaporins. Specifically, a tryptophan residue is extremely conserved within this loop. We have investigated the role of this residue (Trp219) that deeply protrudes into the protein and potentially interacts with adjacent loops, using the E. coli aqualgyeroporin GlpF as a model. Replacement of Trp219 affects the activity of GlpF and impairs the stability of the tetrameric protein. Furthermore, we have identified an amino acid cluster involving Trp219 that stabilizes the GlpF tetramer. Based on our results we propose that Trp219 is key for formation of a defined vestibule structure, which is crucial for glycerol accumulation as well as for the stability of the active GlpF tetramer.

Computing osmotic permeabilities of aquaporins AQP4, AQP5, and GlpF from near-equilibrium simulations

Measuring or computing the single-channel permeability of aquaporins/aquaglyceroporins (AQPs) has long been a challenge. The measured values scatter over an order of magnitude but the corresponding Arrhenius activation energies converge in the current literature. Osmotic flux through an AQP was simulated as water current forced through the channel by kilobar hydraulic pressure or theoretically approximated as single-file diffusion. In this paper, we report large scale simulations of osmotic current under sub M gradient through three AQPs (water channels AQP4 and AQP5 and glycerol-water channel GlpF) using the mature particle mesh Ewald technique (PME) for which the established force fields have been optimized with known accuracy. These simulations were implemented with hybrid periodic boundary conditions devised to avoid the artifactitious mixing across the membrane in a regular PME simulation. The computed single-channel permeabilities at 5°C and 25°C are in agreement with recently refined experiments on GlpF. The Arrhenius activation energies extracted from our simulations for all the three AQPs agree with the in vitro measurements. The single-file diffusion approximations from our large-scale simulations are consistent with the current literature on smaller systems. From these unambiguous agreements among the in vitro and in silico studies, we observe the quantitative accuracy of the all-atom force fields of the current literature for water-channel biology. We also observe that AQP4, that is particularly rich in the central nervous system, is more efficient in water conduction and more temperature-sensitive than other water-only channels (excluding glycerol channels that also conduct water when not inhibited by glycerol).

Simultaneous conversion of glucose and xylose to 3-hydroxypropionic acid in engineered Escherichia coli by modulation of sugar transport and glycerol synthesis

Escherichia coli expressing the Lactobacillus brevis dhaB1B2B3 and dhaR1R2 clusters and Pseudomonas aeruginosa aldhH was engineered to produce 3-HP from glucose and xylose via the glycerol biosynthetic pathway. Glycerol, a key precursor for 3-HP biosynthesis was produced by overexpression of the GPD1 and GPP2 genes from Saccharomyces cerevisiae. For relief of carbon catabolite repression, deletion of the chromosomal ptsG gene and overexpression of the endogenous xylR gene rendered engineered E. coli JHS01300/pCPaGGRm to utilize glucose and xylose simultaneously and to produce glycerol at 0.48 g/g yield and 0.35 g/L-h productivity. Finally, engineered E. coli JHS01300/pELDRR+pCPaGGRm produced 29.4 g/L of 3-HP with 0.54 g/L-h productivity and 0.36 g/g yield in a sugar-limited fed-batch fermentation. It was concluded that dual modulation of sugar transport and glycerol biosynthesis is a promising strategy for efficient conversion of glucose and xylose to 3-HP.

Glycerol modulates water permeation through Escherichia coli aquaglyceroporin GlpF

Among aquaglyceroporins that transport both water and glycerol across the cell membrane, Escherichia coli glycerol uptake facilitator (GlpF) is the most thoroughly studied. However, one question remains: Does glycerol modulate water permeation? This study answers this fundamental question by determining the three-dimensional potential of mean force of glycerol along the permeation path through GlpF's conducting pore. There is a deep well near the Asn-Pro-Ala (NPA) motifs (6.5kcal/mol below the bulk level) and a barrier near the selectivity filter (10.1kcal/mol above the well bottom). This profile owes its existence to GlpF's perfect steric arrangement: The glycerol-protein van der Waals interactions are attractive near the NPA but repulsive elsewhere in the conducting pore. In light of the single-file nature of waters and glycerols lining up in GlpF's amphipathic pore, it leads to the following conclusion: Glycerol modulates water permeation in the μM range. At mM concentrations, GlpF is glycerol-saturated and a glycerol residing in the well occludes the conducting pore. Therefore, water permeation is fully correlated to glycerol dissociation that has an Arrhenius activation barrier of 6.5kcal/mol. Validation of this theory is based on the existent in vitro data, some of which have not been given the proper attention they deserved: The Arrhenius activation barriers were found to be 7kcal/mol for water permeation and 9.6kcal/mol for glycerol permeation; The presence of up to 100mM glycerol did not affect the kinetics of water transport with very low permeability, in apparent contradiction with the existent theories that predicted high permeability (0M glycerol).

Design of a recombinant Escherichia coli for producing L-phenylalanine from glycerol

A recombinant Escherichia coli was engineered to produce the commercially important amino acid L-phenylalanine (L-Phe) using glycerol as the carbon source. Compared to the conventionally used glucose and sucrose, glycerol is a less expensive carbon source. As phenylalanine dehydrogenase (PheDH) activity is involved in the last step of L-Phe synthesis in E. coli, a phenylalanine dehydrogenase gene (phedh) from the thermotolerant Bacillus lentus was cloned into pRSFDuet-1 (pPheDH) and expressed in E. coli BL21(DE3). The resulting clone had a limited ability to produce L-Phe from glycerol, possibly because of a poor glycerol uptake by the cell, or an inability to excrete L-Phe, or both. Therefore, yddG gene encoding an aromatic amino acid exporter and glpF gene encoding a glycerol transport facilitator were coexpressed with the phedh in a reengineered E. coli. In a glycerol medium, the maximum L-Phe production rates of the clones pPY (phedh and yddG genes) and pPYF (phedh, yddG and glpF genes) were 1.4- and 1.8-fold higher than the maximum production rate of the pPheDH clone. The better producing pPYF clone was further evaluated in a 5 l stirred-tank fermenter (37 °C, an aeration rate of 1 vvm, an agitation speed of 400 rpm). In the fermenter, the maximum concentration of L-Phe (366 mg/l) was achieved in a much shorter period compared to in the shake flasks. In the latter, the highest titer of L-Phe was only 76 % of the maximum value attained in the fermenter.

Discrete-state representation of ion permeation coupled to fast gating in a model of CLC-chloride channels: analytic estimation of the state-to-state rate constants

Analytical estimation of state-to-state rate constants is carried out for a recently developed discrete state model of chloride ion motion in a CLC chloride channel (Coalson and Cheng, J. Phys. Chem. B 2010, 114, 1424). In the original presentation of this model, the same rate constants were evaluated via three-dimensional Brownian dynamics simulations. The underlying dynamical theory is an appropriate single- or multiparticle three-dimensional Smoluchowski equation. Taking advantage of approximate geometric symmetries (based on the details of the model channel geometry), well-known formulas for state-to-state transition rates are appealed to herein and adapted as necessary to the problem at hand. Rates of ionic influx from a bulk electrolyte reservoir to the nearest binding site within the channel pore are particularly challenging to compute analytically because they reflect multi-ion interactions (as opposed to single-ion dynamics). A simple empirical correction factor is added to the single-ion rate constant formula in this case to account for the saturation of influx rate constants with increasing bulk Cl(-) concentration. Overall, the agreement between all analytically estimated rate constants is within a factor of 2 of those computed via three-dimensional Brownian dynamics simulations, and often better than this. Current-concentration curves obtained using rate constants derived from these two different computational approaches agree to within 25%.

Crowding effects in non-equilibrium transport through nano-channels

Transport through nano-channels plays an important role in many biological processes and industrial applications. Gaining insights into the functioning of biological transport processes and the design of man-made nano-devices requires an understanding of the basic physics of such transport. A simple exclusion process has proven to be very useful in explaining the properties of several artificial and biological nano-channels. It is particularly useful for modeling the influence of inter-particle interactions on transport characteristics. In this paper, we explore several models of the exclusion process using a mean field approach and computer simulations. We examine the effects of crowding inside the channel and in its immediate vicinity on the mean flux and the transport times of single molecules. Finally, we discuss the robustness of the theory's predictions with respect to the crucial characteristics of the hindered diffusion in nano-channels that need to be included in the model.

Fluxes of non-interacting and strongly repelling particles through a single conical channel: Analytical results and their numerical tests

Using a diffusion model of particle dynamics in the channel, we study entropic effects in channel-facilitated transport. We derive general expressions for the fluxes of non-interacting particles and particles that strongly repel each other through the channel of varying cross section area, assuming that the transport is driven by the difference in particle concentrations on the two sides of the membrane. For a special case of a right truncated cone expanding in the left-to-right direction, we show how the fluxes depend on the geometric parameters of the channel and on the particle concentrations. For non-interacting particles the flux is direction-independent in the sense that inversion of the concentration difference leads to the inversion of the direction of the flux without changing its magnitude. This symmetry is broken for repelling particles: The flux in the left-to-right direction exceeds its right-to-left counterpart. Our theoretical predictions are supported by three-dimensional Brownian dynamics simulations.

Arsenic transport in prokaryotes and eukaryotic microbes

Aquaporins (AQPs) and aquaglyceroporins facilitate transport of a broad spectrum of substrates such as water, glycerol and other small uncharged solutes. More recently, AQPs ave also been shown to facilitate diffusion of metalloids such as arsenic (As) and antimony (Sb). At neutral pH, the trivalent forms of these metalloids are structurally similar to glycerol and hence they can enter cells through AQPs. As- and Sb-containing compounds are toxic to cells, yet both metalloids are used as chemotherapeutic agents for treating acute promyelocytic leukemia and diseases caused by protozoan parasites. In this chapter, we will review the role of AQPs and other proteins in metalloid transport in prokaryotes and eukaryotic microbes.

Enhancement of transport selectivity through nano-channels by non-specific competition

The functioning of living cells requires efficient and selective transport of materials into and out of the cell, and between different cellular compartments. Much of this transport occurs through nano-scale channels that do not require large scale molecular re-arrangements (such as transition from a 'closed' to an 'open' state) and do not require a direct input of metabolic energy during transport. Nevertheless, these 'always open' channels are highly selective and pass only their cognate molecules, while efficiently excluding all others; indeed, these channels can efficiently transport specific molecules even in the presence of a vast excess of non-specific molecules. Such biological transporters have inspired the creation of artificial nano-channels. These channels can be used as nano-molecular sorters, and can also serve as testbeds for examining modes of biological transport. In this paper, we propose a simple kinetic mechanism that explains how the selectivity of such 'always open' channels can be based on the exclusion of non-specific molecules by specific ones, due to the competition for limited space inside the channel. The predictions of the theory account for the behavior of the nuclear pore complex and of artificial nanopores that mimic its function. This theory provides the basis for future work aimed at understanding the selectivity of various biological transport phenomena.

Discrete-state representation of ion permeation coupled to fast gating in a model of ClC chloride channels: comparison to multi-ion continuous space Brownian dynamics simulations

A discrete-state model of chloride ion motion in a ClC chloride channel is constructed, following a previously developed multi-ion continuous space model of the same system (Cheng, M. H.; Mamonov, A. B.; Dukes, J. W.; Coalson, R. D. J. Phys. Chem. B 2007, 111, 5956) that included a simplistic representation of the fast gate in this channel. The reducibility of the many-body continuous space to the eight discrete-state model considered in the present work is examined in detail by performing three-dimensional Brownian dynamics simulations of each allowed state-to-state transition in order to extract the appropriate rate constant for this process, and then inserting the pairwise rate constants thereby obtained into an appropriate set of kinetic master equations. Experimental properties of interest, including the rate of Cl(-) ion permeation through the open channel and the average rate of closing of the fast gate as a function of bulk Cl(-) ion concentrations in the intracellular and extracellular electrolyte reservoirs are computed. Good agreement is found between the results obtained via the eight discrete-state model versus the multi-ion continuous space model, thereby encouraging continued development of the discrete-state model to include more complex behaviors observed experimentally in these channels.

Effects of jamming on nonequilibrium transport times in nanochannels

Many biological channels perform highly selective transport without direct input of metabolic energy and without transitions from a "closed" to an "open" state during transport. Mechanisms of selectivity of such channels serve as an inspiration for creation of artificial nanomolecular sorting devices and biosensors. To elucidate the transport mechanisms, it is important to understand the transport on the single molecule level in the experimentally relevant regime when multiple particles are crowded in the channel. In this Letter we analyze the effects of interparticle crowding on the nonequilibrium transport times through a finite-length channel by means of analytical theory and computer simulations.

Cargo surface hydrophobicity is sufficient to overcome the nuclear pore complex selectivity barrier

To fulfil their function, nuclear pore complexes (NPCs) must discriminate between inert proteins and nuclear transport receptors (NTRs), admitting only the latter. This specific permeation is thought to depend on interactions between hydrophobic patches on NTRs and phenylalanine-glycine (FG) or related repeats that line the NPC. Here, we tested this premise directly by conjugating different hydrophobic amino-acid analogues to the surface of an inert protein and examining its ability to cross NPCs unassisted by NTRs. Conjugation of as few as four hydrophobic moieties was sufficient to enable passage of the protein through NPCs. Transport of the modified protein proceeded with rates comparable to those measured for the innate protein when bound to an NTR and was relatively insensitive both to the nature and density of the amino acids used to confer hydrophobicity. The latter observation suggests a non-specific, small, and plant interaction network between cargo and FG repeats.

Selectivity and cooperativity of modulatory ions in a neurotransmitter receptor

Ions play a modulatory role in many proteins. Kainate receptors, members of the ionotropic glutamate receptor family, require both monovalent anions and cations in the extracellular milieu for normal channel activity. Molecular dynamics simulations and extensive relative binding free energy calculations using thermodynamic integration were performed to elucidate the rank order of binding of monovalent cations, using x-ray crystal structures of the GluR5 kainate receptor dimers with bound cations from the alkali metal family. The simulations show good agreement with experiments and reveal that the underlying backbone structure of the binding site is one of the most rigid regions of the protein. A simplified model where the partial charge of coordinating oxygens was varied suggests that selectivity arises from the presence of two carboxylate groups. Furthermore, using a potential of mean force derived from umbrella sampling, we show that the presence of cations lower the energy barrier for anion approach and binding in the buried anion binding cavity.

Effects of multiple occupancy and interparticle interactions on selective transport through narrow channels: theory versus experiment

Many biological and artificial transport channels function without direct input of metabolic energy during a transport event and without structural rearrangements involving transitions from a closed to an open state. Nevertheless, such channels are able to maintain efficient and selective transport. It has been proposed that attractive interactions between the transported molecules and the channel can increase the transport efficiency and that the selectivity of such channels can be based on the strength of the interaction of the specifically transported molecules with the channel. Herein, we study the transport through narrow channels in a framework of a general kinetic theory, which naturally incorporates multiparticle occupancy of the channel and non-single-file transport. We study how the transport efficiency and the probability of translocation through the channel are affected by interparticle interactions in the confined space inside the channel, and establish conditions for selective transport. We compare the predictions of the model with the available experimental data and find good semiquantitative agreement. Finally, we discuss applications of the theory to the design of artificial nanomolecular sieves.

Efficiency, selectivity, and robustness of nucleocytoplasmic transport

All materials enter or exit the cell nucleus through nuclear pore complexes (NPCs), efficient transport devices that combine high selectivity and throughput. NPC-associated proteins containing phenylalanine–glycine repeats (FG nups) have large, flexible, unstructured proteinaceous regions, and line the NPC. A central feature of NPC-mediated transport is the binding of cargo-carrying soluble transport factors to the unstructured regions of FG nups. Here, we model the dynamics of nucleocytoplasmic transport as diffusion in an effective potential resulting from the interaction of the transport factors with the flexible FG nups, using a minimal number of assumptions consistent with the most well-established structural and functional properties of NPC transport. We discuss how specific binding of transport factors to the FG nups facilitates transport, and how this binding and competition between transport factors and other macromolecules for binding sites and space inside the NPC accounts for the high selectivity of transport. We also account for why transport is relatively insensitive to changes in the number and distribution of FG nups in the NPC, providing an explanation for recent experiments where up to half the total mass of the FG nups has been deleted without abolishing transport. Our results suggest strategies for the creation of artificial nanomolecular sorting devices.

The DNA at the heart of our cells is contained in the nucleus. This nucleus is surrounded by a barrier in which are buried gatekeepers, termed nuclear pore complexes (NPCs), which allow the quick and efficient passage of certain materials while excluding all others. It has long been known that materials must bind to the NPC to be transported across it, but how this binding translates into selective passage through the NPC has remained a mystery. Here we describe a theory to explain how the NPC works. Our theory accounts for the observed characteristics of NPC–mediated transport, and even suggests strategies for the creation of artificial nanomolecular sorting devices.

On achieving experimental accuracy from molecular dynamics simulations of flexible molecules: aqueous glycerol

The rotational isomeric states (RIS) of glycerol at infinite dilution have been characterized in the aqueous phase via a 1 micros conventional molecular dynamics (MD) simulation, a 40 ns enhanced sampling replica exchange molecular dynamics (REMD) simulation, and a reevaluation of the experimental NMR data. The MD and REMD simulations employed the GLYCAM06/AMBER force field with explicit treatment of solvation. The shorter time scale of the REMD sampling method gave rise to RIS and theoretical scalar 3J(HH) coupling constants that were comparable to those from the much longer traditional MD simulation. The 3J(HH) coupling constants computed from the MD methods were in excellent agreement with those observed experimentally. Despite the agreement between the computed and the experimental J-values, there were variations between the rotamer populations computed directly from the MD data and those derived from the experimental NMR data. The experimentally derived populations were determined utilizing limiting J-values from an analysis of NMR data from substituted ethane molecules and may not be completely appropriate for application in more complex molecules, such as glycerol. Here, new limiting J-values have been derived via a combined MD and quantum mechanical approach and were used to decompose the experimental 3J(HH) coupling constants into population distributions for the glycerol RIS.

Diffusion model of solute dynamics in a membrane channel: mapping onto the two-site model and optimizing the flux

The steady-state flux through a singly occupied membrane channel is found for both discrete and continuum models of the solute dynamics in the channel. The former describes the dynamics as nearest-neighbor jumps between N sites, while the latter assumes that the molecule diffuses in a one-dimensional potential of mean force. For both models it is shown that the flux is the same as that for a simple two-site model with appropriately chosen rate constants, which contain all the relevant information about the more detailed dynamics. An interesting consequence of single occupancy is that the flux has a maximum as a function of the channel-solute interaction. If this interaction is too attractive, the molecule will never leave the channel, thus blocking it for the passage of other molecules. If it is too repulsive, the solute molecule will never enter the channel. Thus the flux vanishes in the two limits and, hence, has a maximum somewhere in-between. In the framework of the diffusion model, we find the optimal intrachannel potential of mean force that maximizes the flux using the calculus of variations. For a symmetric channel this potential is flat and occupies the entire channel. In the general case of an asymmetric channel, the optimal potential is obtained by tilting the optimal flat potential for the corresponding symmetric channel around the channel center, so that the solute is driven towards the reservoir with the lower solute concentration by a constant force. This implies that the flux is higher when the solute binding near the channel exit is stronger than that near the entrance.

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 micros. 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.

Comparative functional analysis of aquaporins/glyceroporins in mammals and anurans

Maintenance of fluid homeostasis is critical to establishing and maintaining normal physiology. The landmark discovery of membrane water channels (aquaporins; AQPs) ushered in a new area in osmoregulatory biology that has drawn from and contributed to diverse branches of biology, from molecular biology and genomics to systems biology and evolution, and from microbial and plant biology to animal and translational physiology. As a result, the study of AQPs provides a unique and integrated backdrop for exploring the relationships between genes and genome systems, the regulation of gene expression, and the physiologic consequences of genetic variation. The wide species distribution of AQP family members and the evolutionary conservation of the family indicate that the control of membrane water flux is a critical biological process. AQP function and regulation is proving to be central to many of the pathways involved in individual physiologic systems in both mammals and anurans. In mammals, AQPs are essential to normal secretory and absorptive functions of the eye, lung, salivary gland, sweat glands, gastrointestinal tract, and kidney. In urinary, respiratory, and gastrointestinal systems, AQPs are required for proper urine concentration, fluid reabsorption, and glandular secretions. In anurans, AQPs are important in mediating physiologic responses to changes in the external environment, including those that occur during metamorphosis and adaptation from an aquatic to terrestrial environment and thermal acclimation in anticipation of freezing. Therefore, an understanding of AQP function and regulation is an important aspect of an integrated approach to basic biological research.

Major intrinsic proteins (MIPs) in plants: a complex gene family with major impacts on plant phenotype

The ubiquitous cell membrane proteins called aquaporins are now firmly established as channel proteins that control the specific transport of water molecules across cell membranes in all living organisms. The aquaporins are thus likely to be of fundamental significance to all facets of plant growth and development affected by plant-water relations. A majority of plant aquaporins have been found to share essential structural features with the human aquaporin and exhibit water-transporting ability in various functional assays, and some have been shown experimentally to be of critical importance to plant survival. Furthermore, substantial evidence is now available from a number of plant species that shows differential gene expression of aquaporins in response to abiotic stresses such as salinity, drought, or cold and clearly establishes the aquaporins as major players in the response of plants to conditions that affect water availability. This review summarizes the function and regulation of these genes to develop a greater understanding of the response of plants to water insufficiency, and particularly, to identify tolerant genotypes of major crop species including wheat and rice and plants that are important in agroforestry.

Protein mobility and diffusive barriers inEscherichia coli: consequences of osmotic stress

The effect of osmotic stress on the intracellular diffusion of proteins in Escherichia coli was studied, using a pulsed version of fluorescence recovery after photo-bleaching, pulsed-FRAP. This method employs sequences of laser pulses which only partly bleach the fluorophores in a cell. Because the cell size and geometry are taken into account, pulsed-FRAP enables to measure diffusion in very small cells of different shapes. We found that upon an osmotic upshock from 0.15 to 0.6 Osm, imposed by NaCl or sorbitol, the apparent intracellular diffusion (D) of mobile green fluorescent protein (GFP) decreased from 3.2 to 0.4 microm(2) s(-1), whereas the membrane permeable glycerol had no effect. Exposing E. coli cells to higher osmolalities (> 0.6 Osm) led to compartmentalization of the GFP into discrete pools, from where the GFP could not escape. Although free diffusion through the cell was hindered, the mobility of GFP in these pools was still relatively high (D approximately 0.4 microm(2) s(-1)). The presence of osmoprotectants restored the effect of osmotic stress on the protein mobility and apparent compartmentalization. Also, lowering the osmolality from 0.6 Osm back to 0.15 Osm restored the mobility of GFP. The implications of these findings in terms of heterogeneities and diffusive barriers inside the cell are discussed.

Comparative high-density microarray analysis of gene expression during growth of Lactobacillus helveticus in milk versus rich culture medium

Lactobacillus helveticus CNRZ32 is used by the dairy industry to modulate cheese flavor. The compilation of a draft genome sequence for this strain allowed us to identify and completely sequence 168 genes potentially important for the growth of this organism in milk or for cheese flavor development. The primary aim of this study was to investigate the expression of these genes during growth in milk and MRS medium by using microarrays. Oligonucleotide probes against each of the completely sequenced genes were compiled on maskless photolithography-based DNA microarrays. Additionally, the entire draft genome sequence was used to produce tiled microarrays in which noninterrupted sequence contigs were covered by consecutive 24-mer probes and associated mismatch probe sets. Total RNA isolated from cells grown in skim milk or in MRS to mid-log phase was used as a template to synthesize cDNA, followed by Cy3 labeling and hybridization. An analysis of data from annotated gene probes identified 42 genes that were upregulated during the growth of CNRZ32 in milk (P < 0.05), and 25 of these genes showed upregulation after applying Bonferroni's adjustment. The tiled microarrays identified numerous additional genes that were upregulated in milk versus MRS. Collectively, array data showed the growth of CNRZ32 in milk-induced genes encoding cell-envelope proteinases, oligopeptide transporters, and endopeptidases as well as enzymes for lactose and cysteine pathways, de novo synthesis, and/or salvage pathways for purines and pyrimidines and other functions. Genes for a hypothetical phosphoserine utilization pathway were also differentially expressed. Preliminary experiments indicate that cheese-derived, phosphoserine-containing peptides increase growth rates of CNRZ32 in a chemically defined medium. These results suggest that phosphoserine is used as an energy source during the growth of L. helveticus CNRZ32.

Identity of distributions of direct uphill and downhill translocation times for particles traversing membrane channels

We study the distribution of direct translocation times for particles passing through membrane channels connecting two reservoirs. The direct translocation time is a conditional first-passage time defined as the residence time of the particle in the channel while passing to the other side of the membrane directly, i.e., without returning to the reservoir from which it entered. We show that the distributions of direct translocation times are identical for translocation in both directions, independent of any asymmetry in the potential across the channel and, hence, the translocation probabilities.

What Makes an Aquaporin a Glycerol Channel? A Comparative Study of AqpZ and GlpF

The recent availability of high-resolution structures of two structurally highly homologous, but functionally distinct aquaporins from the same species, namely Escherichia coli AqpZ, a pure water channel, and GlpF, a glycerol channel, presents a unique opportunity to understand the mechanism of substrate selectivity in these channels. Comparison of the free energy profile of glycerol conduction through AqpZ and GlpF reveals a much larger barrier in AqpZ (22.8 kcal/mol) than in GlpF (7.3 kcal/mol). In either channel, the highest barrier is located at the selectivity filter. Analysis of substrate-protein interactions suggests that steric restriction of AqpZ is the main contribution to this large barrier. Another important difference is the presence of a deep energy well at the periplasmic vestibule of GlpF, which was not found in AqpZ. The latter difference can be attributed to the more pronounced structural asymmetry of GlpF, which may play a role in attracting glycerol.

Computational studies of membrane channels

The determination of the structure of several members of the K+ channel and aquaporin family represents a unique opportunity to explain the mechanism of these biomolecular systems. With their ability to go beyond static structures, molecular dynamics simulations offer a unique route for relating functional properties to membrane channel structure. The recent progress in this area is reviewed.

Lipids do influence protein function-the hydrophobic matching hypothesis revisited

A topical review of the current state of lipid-protein interactions is given with focus on the physical interactions between lipids and integral proteins in lipid-bilayer membranes. The concepts of hydrophobic matching and curvature stress are revisited in light of recent data obtained from experimental and theoretical studies which demonstrate that not only do integral proteins perturb the lipids, but the physical state of the lipids does also actively influence protein function. The case of the trans-membrane water-channel protein aquaporin GlpF from E. coli imbedded in lipid-bilayer membranes is discussed in some detail. Numerical data obtained from Molecular Dynamics simulations show on the one side that the lipid bilayer adapts to the channel by a hydrophobic matching condition which reflects the propensity of the lipid molecules for forming curved structures. On the other side, it is demonstrated that the transport function of the channel is modulated by the matching condition and/or the curvature stress in a lipid-specific manner.

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.

Fluctuation-driven molecular transport through an asymmetric membrane channel

Channel proteins that selectively conduct molecules across cell membranes often exhibit an asymmetric structure. By means of a stochastic model, we argue that channel asymmetry in the presence of nonequilibrium fluctuations, fueled by the cell's metabolism as observed recently, can dramatically influence the transport through such channels by a ratchetlike mechanism. For an aquaglyceroporin that conducts water and glycerol, we show that a previously determined asymmetric glycerol potential leads to enhanced inward transport of glycerol, but for unfavorably high glycerol concentrations also to enhanced outward transport that protects a cell against poisoning.

Theory and simulation of water permeation in aquaporin-1

We discuss the difference between osmotic permeability pf and diffusion permeability pd of single-file water channels and demonstrate that the pf/pd ratio corresponds to the number of effective steps a water molecule needs to take to permeate a channel. While pd can be directly obtained from equilibrium molecular dynamics simulations, pf can be best determined from simulations in which a chemical potential difference of water has been established on the two sides of the channel. In light of this, we suggest a method to induce in molecular dynamics simulations a hydrostatic pressure difference across the membrane, from which pf can be measured. Simulations using this method are performed on aquaporin-1 channels in a lipid bilayer, resulting in a calculated pf of 7.1 x 10(-14) cm(3)/s, which is in close agreement with observation. Using a previously determined pd value, we conclude that pf/pd for aquaporin-1 measures approximately 12. This number is explained in terms of channel architecture and conduction mechanism.