A new view of the bacterial cytosol environment

Nernst-Planck-Gaussian finite element modelling of Ca2+ electrodiffusion in amphibian striated muscle transverse tubule-sarcoplasmic reticular triadic junctional domains

Introduction

Intracellular Ca2+ signalling regulates membrane permeabilities, enzyme activity, and gene transcription amongst other functions. Large transmembrane Ca2+ electrochemical gradients and low diffusibility between cell compartments potentially generate short-lived, localised, high-[Ca2+] microdomains. The highest concentration domains likely form between closely apposed membranes, as at amphibian skeletal muscle transverse tubule-sarcoplasmic reticular (T-SR, triad) junctions.

Materials and methods

Finite element computational analysis characterised the formation and steady state and kinetic properties of the Ca2+ microdomains using established empirical physiological and anatomical values. It progressively incorporated Fick diffusion and Nernst-Planck electrodiffusion gradients, K+, Cl-, and Donnan protein, and calmodulin (CaM)-mediated Ca2+ buffering. It solved for temporal-spatial patterns of free and buffered Ca2+, Gaussian charge differences, and membrane potential changes, following Ca2+ release into the T-SR junction.

Results

Computational runs using established low and high Ca2+ diffusibility (D Ca2+) limits both showed that voltages arising from intracytosolic total [Ca2+] gradients and the counterions little affected microdomain formation, although elevated D Ca2+ reduced attained [Ca2+] and facilitated its kinetics. Contrastingly, adopting known cytosolic CaM concentrations and CaM-Ca2+ affinities markedly increased steady-state free ([Ca2+]free) and total ([Ca2+]), albeit slowing microdomain formation, all to extents reduced by high D Ca2+. However, both low and high D Ca2+ yielded predictions of similar, physiologically effective, [Ca2+-CaM]. This Ca2+ trapping by the relatively immobile CaM particularly increased [Ca2+] at the junction centre. [Ca2+]free, [Ca2+-CaM], [Ca2+], and microdomain kinetics all depended on both CaM-Ca2+ affinity and D Ca2+. These changes accompanied only small Gaussian (∼6 mV) and surface charge (∼1 mV) effects on tubular transmembrane potential at either D Ca2+.

Conclusion

These physical predictions of T-SR Ca2+ microdomain formation and properties are compatible with the microdomain roles in Ca2+ and Ca2+-CaM-mediated signalling but limited the effects on tubular transmembrane potentials. CaM emerges as a potential major regulator of both the kinetics and the extent of microdomain formation. These possible cellular Ca2+ signalling roles are discussed in relation to possible feedback modulation processes sensitive to the μM domain but not nM bulk cytosolic, [Ca2+]free, and [Ca2+-CaM], including ryanodine receptor-mediated SR Ca2+ release; Na+, K+, and Cl- channel-mediated membrane excitation and stabilisation; and Na+/Ca2+ exchange transport.

Use of high frequency electrorotation to identify cytoplasmic changes in cells non-disruptively

In this paper we demonstrate how the use of frequencies ranging from 50 kHz to 5 GHz in the analysis of cells by electrorotation can open the path to the identification of differences not detectable by conventional set-ups. Earlier works usually reported electrorotation devices operating below 20 MHz, limiting the response obtained to properties associated with the cell membrane. Those devices are thus unable to resolve the physiological properties in the cytoplasm. We used microwave-based technology to extend the frequency operation to 5 GHz. At high frequencies (from tens of MHz to GHz), the electromagnetic signal passes through the membrane and allows probing the cytoplasm. This enables several applications, such as cell classification, and viability analysis. Additionally, the use of conventional microfabrication techniques reduces the cost and complexity of analysis, compared to other non-invasive methods. We demonstrated the potential of this set-up by identifying two different populations of T-lymphocytes not distinguishable through visual assessment. We also assessed the effect of calcein on cell cytoplasmic properties and used it as a controlled experiment to demonstrate the possibility of this method to detect changes happening predominantly in the cytoplasm.

Mitochondrial signal transduction

The analogy of mitochondria as powerhouses has expired. Mitochondria are living, dynamic, maternally inherited, energy-transforming, biosynthetic, and signaling organelles that actively transduce biological information. We argue that mitochondria are the processor of the cell, and together with the nucleus and other organelles they constitute the mitochondrial information processing system (MIPS). In a three-step process, mitochondria (1) sense and respond to both endogenous and environmental inputs through morphological and functional remodeling; (2) integrate information through dynamic, network-based physical interactions and diffusion mechanisms; and (3) produce output signals that tune the functions of other organelles and systemically regulate physiology. This input-to-output transformation allows mitochondria to transduce metabolic, biochemical, neuroendocrine, and other local or systemic signals that enhance organismal adaptation. An explicit focus on mitochondrial signal transduction emphasizes the role of communication in mitochondrial biology. This framework also opens new avenues to understand how mitochondria mediate inter-organ processes underlying human health.

Building Structural Models of a Whole Mycoplasma Cell

Building structural models of entire cells has been a long-standing cross-discipline challenge for the research community, as it requires an unprecedented level of integration between multiple sources of biological data and enhanced methods for computational modeling and visualization. Here, we present the first 3D structural models of an entire Mycoplasma genitalium (MG) cell, built using the CellPACK suite of computational modeling tools. Our model recapitulates the data described in recent whole-cell system biology simulations and provides a structural representation for all MG proteins, DNA and RNA molecules, obtained by combining experimental and homology-modeled structures and lattice-based models of the genome. We establish a framework for gathering, curating and evaluating these structures, exposing current weaknesses of modeling methods and the boundaries of MG structural knowledge, and visualization methods to explore functional characteristics of the genome and proteome. We compare two approaches for data gathering, a manually-curated workflow and an automated workflow that uses homologous structures, both of which are appropriate for the analysis of mesoscale properties such as crowding and volume occupancy. Analysis of model quality provides estimates of the regularization that will be required when these models are used as starting points for atomic molecular dynamics simulations.

The Effect of the Osmotically Active Compound Concentration Difference on the Passive Water and Proton Fluxes across a Lipid Bilayer

The molecular details of the passive water flux across the hydrophobic membrane interior are still a matter of debate. One of the postulated mechanisms is the spontaneous, water-filled pore opening, which facilitates the hydrophilic connection between aqueous phases separated by the membrane. In the paper, we provide experimental evidence showing that the spontaneous lipid pore formation correlates with the membrane mechanics; hence, it depends on the composition of the lipid bilayer and the concentration of the osmotically active compound. Using liposomes as an experimental membrane model, osmotically induced water efflux was measured with the stopped-flow technique. Shapes of kinetic curves obtained at low osmotic pressure differences are interpreted in terms of two events: the lipid pore opening and water flow across the aqueous channel. The biological significance of the dependence of the lipid pore formation on the concentration difference of an osmotically active compound was illustrated by the demonstration that osmotically driven water flow can be accompanied by the dissipation of the pH gradient. The application of the Helfrich model to describe the probability of lipid pore opening was validated by demonstrating that the probability of pore opening correlates with the membrane bending rigidity. The correlation was determined by experimentally derived bending rigidity coefficients and probabilities of lipid pores opening.

In-cell protein landscapes: making the match between theory, simulation and experiment

Theory, computation and experiment have matched up for the folding of small proteins in vitro, a difficult feat because folding energy landscapes are fairly smooth and free energy differences between states are small. Smoothness means that protein structure and folding are susceptible to the local environment inside living cells. Theory, computation and experiment are now exploring cellular modulation of energy landscapes. Interesting concepts have emerged, such as co-evolution of protein surfaces with their cellular environment to reduce detrimental interactions. Here we look at very recent work beginning to bring together theory, simulations and experiments in the area of protein landscape modulation, to see what problems might be solved in the near future by combining these approaches.

The Similarities between Human Mitochondria and Bacteria in the Context of Structure, Genome, and Base Excision Repair System

Mitochondria emerged from bacterial ancestors during endosymbiosis and are crucial for cellular processes such as energy production and homeostasis, stress responses, cell survival, and more. They are the site of aerobic respiration and adenosine triphosphate (ATP) production in eukaryotes. However, oxidative phosphorylation (OXPHOS) is also the source of reactive oxygen species (ROS), which are both important and dangerous for the cell. Human mitochondria contain mitochondrial DNA (mtDNA), and its integrity may be endangered by the action of ROS. Fortunately, human mitochondria have repair mechanisms that allow protecting mtDNA and repairing lesions that may contribute to the occurrence of mutations. Mutagenesis of the mitochondrial genome may manifest in the form of pathological states such as mitochondrial, neurodegenerative, and/or cardiovascular diseases, premature aging, and cancer. The review describes the mitochondrial structure, genome, and the main mitochondrial repair mechanism (base excision repair (BER)) of oxidative lesions in the context of common features between human mitochondria and bacteria. The authors present a holistic view of the similarities of mitochondria and bacteria to show that bacteria may be an interesting experimental model for studying mitochondrial diseases, especially those where the mechanism of DNA repair is impaired.

The ambivalent role of water at the origins of life

Life as we know it would not exist without water. However, water molecules not only serve as a solvent and reactant but can also promote hydrolysis, which counteracts the formation of essential organic molecules. This conundrum constitutes one of the central issues in origin of life. Hydrolysis is an important part of energy metabolism for all living organisms but only because, inside cells, it is a controlled reaction. How could hydrolysis have been regulated under prebiotic settings? Lower water activities possibly provide an answer: geochemical sites with less free and more bound water can supply the necessary conditions for protometabolic reactions. Such conditions occur in serpentinising systems, hydrothermal sites that synthesise hydrogen gas via rock-water interactions. Here, we summarise the parallels between biotic and abiotic means of controlling hydrolysis in order to narrow the gap between biochemical and geochemical reactions and briefly outline how hydrolysis could even have played a constructive role at the origin of molecular self-organisation.

In-Cell Protein-Protein Contacts: Transient Interactions in the Crowd

Proteins in vivo are immersed in a crowded environment of water, ions, metabolites, and macromolecules. In-cell experiments highlight how transient weak protein-protein interactions promote (via functional "quinary structure") or hinder (via competitive binding or "sticking") complex formation. Computational models of the cytoplasm are expensive. We tackle this challenge with an all-atom model of a small volume of the E. coli cytoplasm to simulate protein-protein contacts up to the 5 μs time scale on the special-purpose supercomputer Anton 2. We use three CHARMM-derived force fields: C22*, C36m, and C36mCU (with CUFIX corrections). We find that both C36m and C36mCU form smaller contact surfaces than C22*. Although CUFIX was developed to reduce protein-protein sticking, larger contacts are observed with C36mCU than C36m. We show that the lifespan Δt of protein-protein contacts obeys a power law distribution between 0.03 and 3 μs, with ∼90% of all contacts lasting <1 μs (similar to the time scale for downhill folding).

In Vivo Water Dynamics in Shewanella oneidensis Bacteria at High Pressure

Following observations of survival of microbes and other life forms in deep subsurface environments it is necessary to understand their biological functioning under high pressure conditions. Key aspects of biochemical reactions and transport processes within cells are determined by the intracellular water dynamics. We studied water diffusion and rotational relaxation in live Shewanella oneidensis bacteria at pressures up to 500 MPa using quasi-elastic neutron scattering (QENS). The intracellular diffusion exhibits a significantly greater slowdown (by −10–30%) and an increase in rotational relaxation times (+10–40%) compared with water dynamics in the aqueous solutions used to resuspend the bacterial samples. Those results indicate both a pressure-induced viscosity increase and slowdown in ionic/macromolecular transport properties within the cells affecting the rates of metabolic and other biological processes. Our new data support emerging models for intracellular organisation with nanoscale water channels threading between macromolecular regions within a dynamically organized structure rather than a homogenous gel-like cytoplasm.

Dynamic cluster formation determines viscosity and diffusion in dense protein solutions

We develop a detailed description of protein translational and rotational diffusion in concentrated solution on the basis of all-atom molecular dynamics simulations in explicit solvent. Our systems contain up to 540 fully flexible proteins with 3.6 million atoms. In concentrated protein solutions (100 mg/mL and higher), the proteins ubiquitin and lysozyme, as well as the protein domains third IgG-binding domain of protein G and villin headpiece, diffuse not as isolated particles, but as members of transient clusters between which they constantly exchange. A dynamic cluster model nearly quantitatively explains the increase in viscosity and the decrease in protein diffusivity with protein volume fraction, which both exceed the predictions from widely used colloid models. The Stokes-Einstein relations for translational and rotational diffusion remain valid, but the effective hydrodynamic radius grows linearly with protein volume fraction. This increase follows the observed increase in cluster size and explains the more dramatic slowdown of protein rotation compared with translation. Baxter's sticky-sphere model of colloidal suspensions captures the concentration dependence of cluster size, viscosity, and rotational and translational diffusion. The consistency between simulations and experiments for a diverse set of soluble globular proteins indicates that the cluster model applies broadly to concentrated protein solutions, with equilibrium dissociation constants for nonspecific protein-protein binding in the Kd ≈ 10-mM regime.

The combination of K+ deficiency with other environmental stresses: What is the outcome?

Potassium (K⁺ ) is a macronutrient known for its high mobility and positive charge, which allows efficient and fast control of the electrical balance and osmotic potential in plant cells. Such features allow K⁺ to remarkably contribute to plant stress adaptation. Some agricultural lands are deficient in K⁺ , imposing a stress that reduces crop yield and makes fertilization a common practice. However, individual stress conditions in the field are rare, and crops usually face a combination of different stresses. As plant response to a stress combination cannot always be deduced from individual stress action, it is necessary to gain insights into the specific mechanisms that connect K⁺ homeostasis with other stress effects to improve plant performance in the context of climate change. Surprisingly, plant responses to environmental stresses under a K⁺ -limiting scenario are poorly understood. In the present review, we summarize current knowledge and find substantial gaps regarding specific outcomes of K⁺ deficiency in addition to other environmental stresses. In this regard, combined nutrient deficiencies of K⁺ and other macronutrients are covered in the first part of the review and interactions arising from K⁺ deficiency with salinity, drought and biotic factors in the second part. Information available so far suggests a prominent role of potassium and nitrate transport systems and their regulatory proteins in the response of plants to several stress combinations. Thus, such molecular pathways, which are located at the crossroad between K⁺ homeostasis and environmental stresses, could be considered biotechnological targets in future studies.

Towards capturing cellular complexity: combining encapsulation and macromolecular crowding in a reverse micelle

This paper studies the orientational structure and dynamics of multi-protein systems under confinement and discusses the implications on biological cells.

Macromolecular crowding and the importance of proper hydration for the structure and dynamics of protein solutions

Extensive computational studies of ubiquitin crowding with a special focus on protein hydration directly visible in dielectric spectra.

Crowding in Cellular Environments at an Atomistic Level from Computer Simulations

The effects of crowding in biological environments on biomolecular structure, dynamics, and function remain not well understood. Computer simulations of atomistic models of concentrated peptide and protein systems at different levels of complexity are beginning to provide new insights. Crowding, weak interactions with other macromolecules and metabolites, and altered solvent properties within cellular environments appear to remodel the energy landscape of peptides and proteins in significant ways including the possibility of native state destabilization. Crowding is also seen to affect dynamic properties, both conformational dynamics and diffusional properties of macromolecules. Recent simulations that address these questions are reviewed here and discussed in the context of relevant experiments.

Biomolecular interactions modulate macromolecular structure and dynamics in atomistic model of a bacterial cytoplasm

Biological macromolecules function in highly crowded cellular environments. The structure and dynamics of proteins and nucleic acids are well characterized in vitro, but in vivo crowding effects remain unclear. Using molecular dynamics simulations of a comprehensive atomistic model cytoplasm we found that protein-protein interactions may destabilize native protein structures, whereas metabolite interactions may induce more compact states due to electrostatic screening. Protein-protein interactions also resulted in significant variations in reduced macromolecular diffusion under crowded conditions, while metabolites exhibited significant two-dimensional surface diffusion and altered protein-ligand binding that may reduce the effective concentration of metabolites and ligands in vivo. Metabolic enzymes showed weak non-specific association in cellular environments attributed to solvation and entropic effects. These effects are expected to have broad implications for the in vivo functioning of biomolecules. This work is a first step towards physically realistic in silico whole-cell models that connect molecular with cellular biology.

Roles and Transport of Sodium and Potassium in Plants

The two alkali cations Na(+) and K(+) have similar relative abundances in the earth crust but display very different distributions in the biosphere. In all living organisms, K(+) is the major inorganic cation in the cytoplasm, where its concentration (ca. 0.1 M) is usually several times higher than that of Na(+). Accumulation of Na(+) at high concentrations in the cytoplasm results in deleterious effects on cell metabolism, e.g., on photosynthetic activity in plants. Thus, Na(+) is compartmentalized outside the cytoplasm. In plants, it can be accumulated at high concentrations in vacuoles, where it is used as osmoticum. Na(+) is not an essential element in most plants, except in some halophytes. On the other hand, it can be a beneficial element, by replacing K(+) as vacuolar osmoticum for instance. In contrast, K(+) is an essential element. It is involved in electrical neutralization of inorganic and organic anions and macromolecules, pH homeostasis, control of membrane electrical potential, and the regulation of cell osmotic pressure. Through the latter function in plants, it plays a role in turgor-driven cell and organ movements. It is also involved in the activation of enzymes, protein synthesis, cell metabolism, and photosynthesis. Thus, plant growth requires large quantities of K(+) ions that are taken up by roots from the soil solution, and then distributed throughout the plant. The availability of K(+) ions in the soil solution, slowly released by soil particles and clays, is often limiting for optimal growth in most natural ecosystems. In contrast, due to natural salinity or irrigation with poor quality water, detrimental Na(+) concentrations, toxic for all crop species, are present in many soils, representing 6 % to 10 % of the earth's land area. Three families of ion channels (Shaker, TPK/KCO, and TPC) and 3 families of transporters (HAK, HKT, and CPA) have been identified so far as contributing to K(+) and Na(+) transport across the plasmalemma and internal membranes, with high or low ionic selectivity. In the model plant Arabidopsis thaliana, these families gather at least 70 members. Coordination of the activities of these systems, at the cell and whole plant levels, ensures plant K(+) nutrition, use of Na(+) as a beneficial element, and adaptation to saline conditions.

Peptide folding in the presence of interacting protein crowders

Using Monte Carlo methods, we explore and compare the effects of two protein crowders, BPTI and GB1, on the folding thermodynamics of two peptides, the compact helical trp-cage and the β-hairpin-forming GB1m3. The thermally highly stable crowder proteins are modeled using a fixed backbone and rotatable side-chains, whereas the peptides are free to fold and unfold. In the simulations, the crowder proteins tend to distort the trp-cage fold, while having a stabilizing effect on GB1m3. The extent of the effects on a given peptide depends on the crowder type. Due to a sticky patch on its surface, BPTI causes larger changes than GB1 in the melting properties of the peptides. The observed effects on the peptides stem largely from attractive and specific interactions with the crowder surfaces, and differ from those seen in reference simulations with purely steric crowder particles.

Assessing the potential of atomistic molecular dynamics simulations to probe reversible protein-protein recognition and binding

Protein-protein recognition and binding are governed by diffusion, noncovalent forces and conformational flexibility, entangled in a way that only molecular dynamics simulations can dissect at high resolution. Here we exploited ubiquitin's noncovalent dimerization equilibrium to assess the potential of atomistic simulations to reproduce reversible protein-protein binding, by running submicrosecond simulations of systems with multiple copies of the protein at millimolar concentrations. The simulations essentially fail because they lead to aggregates, yet they reproduce some specificity in the binding interfaces as observed in known covalent and noncovalent ubiquitin dimers. Following similar observations in literature we hint at electrostatics and water descriptions as the main liable force field elements, and propose that their optimization should consider observables relevant to multi-protein systems and unfolded proteins. Within limitations, analysis of binding events suggests salient features of protein-protein recognition and binding, to be retested with improved force fields. Among them, that specific configurations of relative direction and orientation seem to trigger fast binding of two molecules, even over 50 Å distances; that conformational selection can take place within surface-to-surface distances of 10 to 40 Å i.e. well before actual intermolecular contact; and that establishment of contacts between molecules further locks their conformations and relative orientations.

Complete atomistic model of a bacterial cytoplasm for integrating physics, biochemistry, and systems biology

A model for the cytoplasm of Mycoplasma genitalium is presented that integrates data from a variety of sources into a physically and biochemically consistent model. Based on gene annotations, core genes expected to be present in the cytoplasm were determined and a metabolic reaction network was reconstructed. The set of cytoplasmic genes and metabolites from the predicted reactions were assembled into a comprehensive atomistic model consisting of proteins with predicted structures, RNA, protein/RNA complexes, metabolites, ions, and solvent. The resulting model bridges between atomistic and cellular scales, between physical and biochemical aspects, and between structural and systems views of cellular systems and is meant as a starting point for a variety of simulation studies.

An evaluation of genetically encoded FRET-based biosensors for quantitative metabolite analyses in vivo

A broad range of genetically-encoded fluorescence biosensors has been developed, allowing the detection of signaling intermediates and metabolites in real time. Many of these biosensors are based on Foerster Resonance Energy Transfer (FRET). The two biosensors of the well-known "Venus-flytrap" type exemplarily studied in this work are composed of a central sugar binding protein flanked by two green fluorescent protein derivatives, namely ECFP as well as Citrine and EYFP, respectively. In order to evaluate FRET-based biosensors as an in vivo tool for quantitative metabolite analyses, we have thoroughly studied the effects of pH, buffer salts, ionic strength, temperature and several intracellular metabolites on the signal intensity of both biosensors and both fluorescence proteins. Almost all micro-environmental variations led to considerably different FRET signals, because either the fluorescent proteins or the metabolite binding domains were affected by the tested parameters. This resulted not only in altered FRET ratios between the apo state and the saturated state but also in significant shifts of the apparent binding constant. This underlines the necessity of careful controls in order to allow reliable quantitative measurements in vivo.

Dissecting the effects of concentrated carbohydrate solutions on protein diffusion, hydration, and internal dynamics

We present herein a thorough description of the effects of high glucose concentrations on the diffusion, hydration and internal dynamics of ubiquitin, as predicted from extensive molecular dynamics simulations on several systems described at fully atomistic level. We observe that the protein acts as a seed that speeds up the natural propensity of glucose to cluster at high concentration; the sugar molecules thus aggregate around the protein trapping it inside a dynamic cage. This process extensively dehydrates the protein surface, restricts the motions of the remaining water molecules, and drags the large-scale, collective motions of protein atoms slowing down the rate of exploration of the conformational space despite only a slight dampening of fast, local dynamics. We discuss how these effects could be relevant to the function of sugars as preservation agents in biological materials, and how crowding by small sticky molecules could modulate proteins across different reaction coordinates inside the cellular cytosol.

A mathematical model of metabolism and regulation provides a systems-level view of how Escherichia coli responds to oxygen

The efficient redesign of bacteria for biotechnological purposes, such as biofuel production, waste disposal or specific biocatalytic functions, requires a quantitative systems-level understanding of energy supply, carbon, and redox metabolism. The measurement of transcript levels, metabolite concentrations and metabolic fluxes per se gives an incomplete picture. An appreciation of the interdependencies between the different measurement values is essential for systems-level understanding. Mathematical modeling has the potential to provide a coherent and quantitative description of the interplay between gene expression, metabolite concentrations, and metabolic fluxes. Escherichia coli undergoes major adaptations in central metabolism when the availability of oxygen changes. Thus, an integrated description of the oxygen response provides a benchmark of our understanding of carbon, energy, and redox metabolism. We present the first comprehensive model of the central metabolism of E. coli that describes steady-state metabolism at different levels of oxygen availability. Variables of the model are metabolite concentrations, gene expression levels, transcription factor activities, metabolic fluxes, and biomass concentration. We analyze the model with respect to the production capabilities of central metabolism of E. coli. In particular, we predict how precursor and biomass concentration are affected by product formation.

Tetramers are the activation-competent species of the HOCl-specific transcription factor HypT

Hypochlorous acid (HOCl) is an important component of the immune system and is produced by neutrophils to kill invading microorganisms. The transcription factor HypT is specifically activated by HOCl by methionine oxidation and protects Escherichia coli cells from the detrimental effects of HOCl. HypT forms dodecameric ring-like oligomers. Binding of HypT to DNA induces dissociation of the dodecamers into dimers and tetramers, thus forming the DNA-binding species. To dissect HypT dissociation, binding to DNA, and activation, we aimed to dissociate the dodecamers independently of DNA and to analyze HOCl-dependent activation in vitro. We found that HypT dodecamers dissociated into tetramers in the presence of l-arginine and NaCl, which was reversible upon dilution of the additive. Making use of the reversible dissociation, we generated mixed assemblies consisting of wild-type and mutant HypT subunits and determined that mutant subunits with reduced thermal stability were stabilized by wild-type HypT in the mixed assembly. HypT tetramers, as present at high NaCl concentrations, were stabilized against thermal unfolding and aggregation triggered by high HOCl concentrations. Importantly, in vitro activation by HOCl of HypT tetramers was completed within 1 min, whereas activation of dodecamers required 1 h for completion. Furthermore, activation of HypT tetramers required stoichiometric amounts of HOCl instead of an excess of HOCl, as observed for dodecamers. This supports the idea that small HypT oligomers are the activation-competent species, whereas the dodecamers are a storage form. Our study reveals the importance of the dynamic oligomeric state for HypT activation by HOCl.

Molecular dynamics simulations of highly crowded amino acid solutions: comparisons of eight different force field combinations with experiment and with each other

Although it is now commonly accepted that the highly crowded conditions encountered inside biological cells have the potential to significantly alter the thermodynamic properties of biomolecules, it is not known to what extent the thermodynamics of fundamental types of interactions such as salt bridges and hydrophobic interactions are strengthened or weakened by high biomolecular concentrations. As one way of addressing this question we have performed a series of all-atom explicit solvent molecular dynamics (MD) simulations to investigate the effect of increasing solute concentration on the behavior of four types of zwitterionic amino acids in aqueous solution. We have simulated systems containing glycine, valine, phenylalanine or asparagine at concentrations of 50, 100, 200 and 300 mg/ml. Each molecular system has been simulated for 1 μs in order to obtain statistically converged estimates of thermodynamic parameters, and each has been conducted with 8 different force fields and water models; the combined simulation time is 128 μs. The density, viscosity, and dielectric increments of the four amino acids calculated from the simulations have been compared to corresponding experimental measurements. While all of the force fields perform well at reproducing the density increments, discrepancies for the viscosity and dielectric increments raise questions both about the accuracy of the simulation force fields and, in certain cases, the experimental data. We also observe large differences between the various force fields' descriptions of the interaction thermodynamics of salt bridges and, surprisingly, these differences also lead to qualitatively different predictions of their dependences on solute concentration. For the aliphatic interactions of valine sidechains, fewer differences are observed between the force fields, but significant differences are again observed for aromatic interactions of phenylalanine sidechains. Taken together, the results highlight the potential power of using explicit-solvent simulation methods to understand behavior in concentrated systems but also hint at potential difficulties in using these methods to obtain consistent views of behavior in intracellular environments.

Reaching new levels of realism in modeling biological macromolecules in cellular environments

An increasing number of studies are aimed at modeling cellular environments in a comprehensive and realistic fashion. A major challenge in these efforts is how to bridge spatial and temporal scales over many orders of magnitude. Furthermore, there are additional challenges in integrating different aspects ranging from questions about biomolecular stability in crowded environments to the description of reactive processes on cellular scales. In this review, recent studies with models of biomolecules in cellular environments at different levels of detail are discussed in terms of their strengths and weaknesses. In particular, atomistic models, implicit representations of cellular environments, coarse-grained and spheroidal models of biomolecules, as well as the inclusion of reactive processes via reaction-diffusion models are described. Furthermore, strategies for integrating the different models into a comprehensive description of cellular environments are discussed.

Computer Simulations of the Bacterial Cytoplasm

Ever since the pioneering work of Minton, it has been recognized that the highly crowded interior of biological cells has the potential to cause dramatic changes to both the kinetics and thermodynamics of protein folding and association events relative to behavior that might be observed in dilute solution conditions. One very productive way to explore the effects of crowding on protein behavior has been to use macromolecular crowding agents that exclude volume without otherwise strongly interacting with the protein under study. An alternative, complementary approach to understanding the potential differences between behavior in vivo and in vitro is to develop simulation models that explicitly attempt to model intracellular environments at the molecular scale, and that thereby can be used to directly monitor biophysical behavior in conditions that accurately mimic those encountered in vivo. It is with studies of this type that the present review will be concerned. We review in detail four published studies that have attempted to simulate the structure and dynamics of the bacterial cytoplasm and that have each explored different biophysical aspects of the cellular interior. While each of these studies has yielded important new insights, there are important questions that remain to be resolved in terms of determining the relative contributions made by energetic and hydrodynamic interactions to the diffusive behavior of macromolecules and to the thermodynamics of protein folding and associations in vivo. Some possible new directions for future generation simulation models of the cytoplasm are outlined.

Mechanisms determining the fate of dispersed bacterial communities in new environments

Recent work has shown that dispersal has an important role in shaping microbial communities. However, little is known about how dispersed bacteria cope with new environmental conditions and how they compete with local resident communities. To test this, we implemented two full-factorial transplant experiments with bacterial communities originating from two sources (freshwater or saline water), which were incubated, separately or in mixes, under both environmental conditions. Thus, we were able to separately test for the effects of the new environment with and without interactions with local communities. We determined community composition using 454-pyrosequencing of bacterial 16S rRNA to specifically target the active fraction of the communities, and measured several functional parameters. In absence of a local resident community, the net functional response was mainly affected by the environmental conditions, suggesting successful functional adaptation to the new environmental conditions. Community composition was influenced both by the source and the incubation environment, suggesting simultaneous effects of species sorting and functional plasticity. In presence of a local resident community, functional parameters were higher compared with those expected from proportional mixes of the unmixed communities in three out of four cases. This was accompanied by an increase in the relative abundance of generalists, suggesting that competitive interactions among local and immigrant taxa could explain the observed 'functional overachievement'. In summary, our results suggest that environmental filtering, functional plasticity and competition are all important mechanisms influencing the fate of dispersed communities.

Standardized assay medium to measure Lactococcus lactis enzyme activities while mimicking intracellular conditions

Knowledge of how the activity of enzymes is affected under in vivo conditions is essential for analyzing their regulation and constructing models that yield an integrated understanding of cell behavior. Current kinetic parameters for Lactococcus lactis are scattered through different studies and performed under different assay conditions. Furthermore, assay conditions often diverge from conditions prevailing in the intracellular environment. To establish uniform assay conditions that resemble intracellular conditions, we analyzed the intracellular composition of anaerobic glucose-limited chemostat cultures of L. lactis subsp. cremoris MG 1363. Based on this, we designed a new assay medium for enzyme activity measurements of growing cells of L. lactis, mimicking as closely as practically possible its intracellular environment. Procedures were optimized to be carried out in 96-well plates, and the reproducibility and dynamic range were checked for all enzyme activity measurements. The effects of freezing and the carryover of ammonium sulfate from the addition of coupling enzymes were also established. Activities of all 10 glycolytic and 4 fermentative enzymes were measured. Remarkably, most in vivo-like activities were lower than previously published data. Yet, the ratios of V(max) over measured in vivo fluxes were above 1. With this work, we have developed and extensively validated standard protocols for enzyme activity measurements for L. lactis.