Crowding Effects on Association Reactions at Membranes

The effect of macromolecular crowding on the binding of ligands to a receptor near membranes is studied using Brownian dynamics simulations. The receptor is modeled as a reactive patch on a hard surface and the ligands and crowding agents are modeled as spheres that interact via a steep repulsive interaction potential. When a ligand collides with the patch, it reacts with probability p_rxn. The association rate constant (k_∞) can be decomposed into contributions from diffusion-limited (k_D) and reaction-limited (k_R) rates, i.e., 1/k_∞ = 1/k_D + 1/k_R. The simulations show that k_D is a nonmonotonic function of the volume fraction of crowding agents for receptors of small sizes. k_R is always an increasing function of the volume fraction of crowding agents, and the association rate constant k_∞ determined from both contributions has a qualitatively different dependence on the macromolecular crowding for high and low values of the reaction probability p_rxn. The simulation results are used to predict the velocity of the membrane protrusion driven by actin filament elongation. Based on the simple model where the protrusive force on the membrane is generated by the intercalation of actin monomers between the membrane and actin filament ends, we predict that crowding increases the local concentration of actin monomers near the filament ends and hence accelerates the membrane protrusion.     

Macromolecular Crowding Modulates Folding Mechanism of α/β Protein Apoflavodoxin

Protein dynamics in cells may be different from those in dilute solutions in vitro, because the environment in cells is highly concentrated with other macromolecules. This volume exclusion because of macromolecular crowding is predicted to affect both equilibrium and kinetic processes involving protein conformational changes. To quantify macromolecular crowding effects on protein folding mechanisms, we investigated the folding energy landscape of an α/β protein, apoflavodoxin, in the presence of inert macromolecular crowding agents, using in silico and in vitro approaches. By means of coarse-grained molecular simulations and topology-based potential interactions, we probed the effects of increased volume fractions of crowding agents (ϕ_c) as well as of crowding agent geometry (sphere or spherocylinder) at high ϕ_c. Parallel kinetic folding experiments with purified Desulfovibro desulfuricans apoflavodoxin in vitro were performed in the presence of Ficoll (sphere) and Dextran (spherocylinder) synthetic crowding agents. In conclusion, we identified the in silico crowding conditions that best enhance protein stability, and discovered that upon manipulation of the crowding conditions, folding routes experiencing topological frustrations can be either enhanced or relieved. Our test-tube experiments confirmed that apoflavodoxin''s time-resolved folding path is modulated by crowding agent geometry. Macromolecular crowding effects may be a tool for the manipulation of protein-folding and function in living cells.     

Coarse-grained molecular simulation of diffusion and reaction kinetics in a crowded virtual cytoplasm

We present a general-purpose model for biomolecular simulations at the molecular level that incorporates stochasticity, spatial dependence, and volume exclusion, using diffusing and reacting particles with physical dimensions. To validate the model, we first established the formal relationship between the microscopic model parameters (timestep, move length, and reaction probabilities) and the macroscopic coefficients for diffusion and reaction rate. We then compared simulation results with Smoluchowski theory for diffusion-limited irreversible reactions and the best available approximation for diffusion-influenced reversible reactions. To simulate the volumetric effects of a crowded intracellular environment, we created a virtual cytoplasm composed of a heterogeneous population of particles diffusing at rates appropriate to their size. The particle-size distribution was estimated from the relative abundance, mass, and stoichiometries of protein complexes using an experimentally derived proteome catalog from Escherichia coli K12. Simulated diffusion constants exhibited anomalous behavior as a function of time and crowding. Although significant, the volumetric impact of crowding on diffusion cannot fully account for retarded protein mobility in vivo, suggesting that other biophysical factors are at play. The simulated effect of crowding on barnase-barstar dimerization, an experimentally characterized example of a bimolecular association reaction, reveals a biphasic time course, indicating that crowding exerts different effects over different timescales. These observations illustrate that quantitative realism in biosimulation will depend to some extent on mesoscale phenomena that are not currently well understood.     

From microscopy data to in silico environments for in vivo-oriented simulations

In our previous study, we introduced a combination methodology of Fluorescence Correlation Spectroscopy (FCS) and Transmission Electron Microscopy (TEM), which is powerful to investigate the effect of intracellular environment to biochemical reaction processes. Now, we developed a reconstruction method of realistic simulation spaces based on our TEM images. Interactive raytracing visualization of this space allows the perception of the overall 3D structure, which is not directly accessible from 2D TEM images. Simulation results show that the diffusion in such generated structures strongly depends on image post-processing. Frayed structures corresponding to noisy images hinder the diffusion much stronger than smooth surfaces from denoised images. This means that the correct identification of noise or structure is significant to reconstruct appropriate reaction environment in silico in order to estimate realistic behaviors of reactants in vivo. Static structures lead to anomalous diffusion due to the partial confinement. In contrast, mobile crowding agents do not lead to anomalous diffusion at moderate crowding levels. By varying the mobility of these non-reactive obstacles (NRO), we estimated the relationship between NRO diffusion coefficient (D_nro) and the anomaly in the tracer diffusion (α). For D_nro=21.96 to 44.49 μm²/s, the simulation results match the anomaly obtained from FCS measurements. This range of the diffusion coefficient from simulations is compatible with the range of the diffusion coefficient of structural proteins in the cytoplasm. In addition, we investigated the relationship between the radius of NRO and anomalous diffusion coefficient of tracers by the comparison between different simulations. The radius of NRO has to be 58 nm when the polymer moves with the same diffusion speed as a reactant, which is close to the radius of functional protein complexes in a cell.     

Atomistic Modeling of Macromolecular Crowding Predicts Modest Increases in Protein Folding and Binding Stability

Theoretical models predict that macromolecular crowding can increase protein folding stability, but depending on details of the models (e.g., how the denatured state is represented), the level of stabilization predicted can be very different. In this study, we represented the native and denatured states atomistically, with conformations sampled from explicit-solvent molecular dynamics simulations at room temperature and high temperature, respectively. We then designed an efficient algorithm to calculate the allowed fraction, f, when the protein molecule is placed inside a box of crowders. That a fraction of placements of the protein molecule is disallowed because of volume exclusion by the crowders leads to an increase in chemical potential, given by Δμ = −k_BT lnf. The difference in Δμ between the native and denatured states predicts the effect of crowding on the folding free energy. Even when the crowders occupied 35% of the solution volume, the stabilization reached only 1.5 kcal/mol for cytochrome b₅₆₂. The modest stabilization predicted is consistent with experimental studies. Interestingly, a mixture of different sized crowders was found to exert a greater effect than the sum of the individual species of crowders. The stabilization of crowding on the binding stability of barnase and barstar, based on atomistic modeling of the proteins, was similarly modest. These findings have profound implications for macromolecular crowding inside cells.     

The Folding Transition-State Ensemble of a Four-Helix Bundle Protein: Helix Propensity as a Determinant and Macromolecular Crowding as a Probe

The four-helix bundle protein Rd-apocyt b₅₆₂, a redesigned stable variant of apocytochrome b₅₆₂, exhibits two-state folding kinetics. Its transition-state ensemble has been characterized by Φ-value analysis. To elucidate the molecular basis of the transition-state ensemble, we have carried out high-temperature molecular dynamics simulations of the unfolding process. In six parallel simulations, unfolding started with the melting of helix I and the C-terminal half of helix IV, and followed by helix III, the N-terminal half of helix IV and helix II. This ordered melting of the helices is consistent with the conclusion from native-state hydrogen exchange, and can be rationalized by differences in intrinsic helix propensity. Guided by experimental Φ-values, a putative transition-state ensemble was extracted from the simulations. The residue helical probabilities of this transition-state ensemble show good correlation with the Φ-values. To further validate the putative transition-state ensemble, the effect of macromolecular crowding on the relative stability between the unfolded ensemble and the transition-state ensemble was calculated. The resulting effect of crowding on the folding kinetics agrees well with experimental observations. This study shows that molecular dynamics simulations combined with calculation of crowding effects provide an avenue for characterize the transition-state ensemble in atomic details.     

Non-linear effects of macromolecular crowding on enzymatic activity of multi-copper oxidase

Enzymes catalyze biochemical reactions in highly crowded environments where the amount of macromolecules may occupy up to 40% of the volume. Here we report how cell-like conditions tune catalytic parameters for the monomeric multi-copper oxidase, Saccharomyces cerevisiae Fet3p, in vitro. At low amounts of crowding agent, we detect increases in both of K_M (weaker substrate binding) and k_cat (improved catalytic efficiency), whereas at higher crowding levels, both parameters were reduced. Presence of crowding agents does not affect Fet3p structural content but increases thermal resistance. The observations are compatible with ordering of a non-optimal substrate-binding site and restricted internal dynamics as a result of excluded volume effects making the protein less structurally ‘strained’.     

Combined Effects of Agitation,Macromolecular Crowding,and Interfaces on Amyloidogenesis

Amyloid formation and accumulation is a hallmark of protein misfolding diseases and is associated with diverse pathologies including type II diabetes and Alzheimer''s disease (AD). In vitro, amyloidogenesis is widely studied in conditions that do not simulate the crowded and viscous in vivo environment. A high volume fraction of most biological fluids is occupied by various macromolecules, a phenomenon known as macromolecular crowding. For some amyloid systems (e.g. α-synuclein) and under shaking condition, the excluded volume effect of macromolecular crowding favors aggregation, whereas increased viscosity reduces the kinetics of these reactions. Amyloidogenesis can also be catalyzed by hydrophobic-hydrophilic interfaces, represented by the air-water interface in vitro and diverse heterogeneous interfaces in vivo (e.g. membranes). In this study, we investigated the effects of two different crowding polymers (dextran and Ficoll) and two different experimental conditions (with and without shaking) on the fibrilization of amyloid-β peptide, a major player in AD pathogenesis. Specifically, we demonstrate that, during macromolecular crowding, viscosity dominates over the excluded volume effect only when the system is spatially non homogeneous (i.e. an air-water interface is present). We also show that the surfactant activity of the crowding agents can critically influence the outcome of macromolecular crowding and that the structure of the amyloid species formed may depend on the polymer used. This suggests that, in vivo, the outcome of amyloidogenesis may be affected by both macromolecular crowding and spatial heterogeneity (e.g. membrane turn-over). More generally, our work suggests that any factors causing changes in crowding may be susceptibility factors in AD.     

Quantification of Excluded Volume Effects on the Folding Landscape of Pseudomonas aeruginosa Apoazurin In Vitro

Proteins fold and function inside cells that are crowded with macromolecules. Here, we address the role of the resulting excluded volume effects by in vitro spectroscopic studies of Pseudomonas aeruginosa apoazurin stability (thermal and chemical perturbations) and folding kinetics (chemical perturbation) as a function of increasing levels of crowding agents dextran (sizes 20, 40, and 70 kDa) and Ficoll 70. We find that excluded volume theory derived by Minton quantitatively captures the experimental effects when crowding agents are modeled as arrays of rods. This finding demonstrates that synthetic crowding agents are useful for studies of excluded volume effects. Moreover, thermal and chemical perturbations result in free energy effects by the presence of crowding agents that are identical, which shows that the unfolded state is energetically the same regardless of method of unfolding. This also underscores the two-state approximation for apoazurin’s unfolding reaction and suggests that thermal and chemical unfolding experiments can be used in an interchangeable way. Finally, we observe increased folding speed and invariant unfolding speed for apoazurin in the presence of macromolecular crowding agents, a result that points to unfolded-state perturbations. Although the absolute magnitude of excluded volume effects on apoazurin is only on the order of 1–3 kJ/mol, differences of this scale may be biologically significant.     

Quantification of Excluded Volume Effects on the Folding Landscape of Pseudomonas aeruginosa Apoazurin In?Vitro

Proteins fold and function inside cells that are crowded with macromolecules. Here, we address the role of the resulting excluded volume effects by in vitro spectroscopic studies of Pseudomonas aeruginosa apoazurin stability (thermal and chemical perturbations) and folding kinetics (chemical perturbation) as a function of increasing levels of crowding agents dextran (sizes 20, 40, and 70 kDa) and Ficoll 70. We find that excluded volume theory derived by Minton quantitatively captures the experimental effects when crowding agents are modeled as arrays of rods. This finding demonstrates that synthetic crowding agents are useful for studies of excluded volume effects. Moreover, thermal and chemical perturbations result in free energy effects by the presence of crowding agents that are identical, which shows that the unfolded state is energetically the same regardless of method of unfolding. This also underscores the two-state approximation for apoazurin’s unfolding reaction and suggests that thermal and chemical unfolding experiments can be used in an interchangeable way. Finally, we observe increased folding speed and invariant unfolding speed for apoazurin in the presence of macromolecular crowding agents, a result that points to unfolded-state perturbations. Although the absolute magnitude of excluded volume effects on apoazurin is only on the order of 1–3 kJ/mol, differences of this scale may be biologically significant.     

Macromolecular crowding-induced molten globule states of the alkali pH-denatured proteins

Structural and molecular properties extracted from circular dichroism (CD), tryptophan fluorescence and 1-anilino-8-napthalene sulfonate (ANS) binding experiments suggest that the high concentration of synthetic crowding agents (dextran 40, dextran 70 and ficoll 70) stabilizes and refolds the base-denatured ferricytochrome c (Ferricyt c) and lysozyme (Lyz) at pH 12.9 (±0.1) to molten globule (MG) states (C_B-states). These results further revealed that the C_B-states resemble the generic properties of MG-states. Thermodynamic analysis of thermal denaturation curves of base-denatured Ferricyt c and Lyz at pH 12.9 (±0.1) under variable concentrations of crowding agents (dextran 40, dextran 70 and ficoll 70) revealed that the crowder presence increases the thermal stability of base-denatured proteins and also prevents the cold denaturation of Ferricyt c. The results further showed that the nature, size and shape of crowder influence the crowding-mediated increase in secondary structure stabilization and thermal stability of base-denatured Ferricyt c and Lyz. Analysis of kinetic and thermodynamic parameters measured for CO association reaction of alkaline ferrocytochrome c (Ferrocyt c) at pH 12.9 (±0.1) under variable concentrations of crowding agents (dextran 40, dextran 70 and ficoll 70) revealed that the crowder presence reduces the level of structural fluctuation of M80-containing Ω-loop that control CO association to alkaline Ferrocyt c.     

Encounter probabilities between male and female rotifers

The relationship between population density and sexual reproduction is explored for planktonic rotifer populations. It has been argued that selection favors sex at population maxima because of favorable probabilities of encounter. Factors determining the probability of encounter are examined. Experimentally determined parameters for males and females (swimming speed, duration of fertility, net fecundity, and encounter distance) are used to investigate density effects on the probability of encounter. As population density increases, the probability of encounter also increases sigmoidally. The probability of a female encountering a male during her fertile period begins to increase rapidly at ≈ 50 males/1. The increase continues until it becomes asymptotic at ≈ 1000 males/1, where the probability of a female encountering, at least one male during her fertile period approaches one. The duration of female fertility is less than that of males, causing the female curve to be displaced toward higher population densities. Critical density, the population density where R_o′ is exactly equal to one, for males and females is 11.7 and 141 individuals/1, respectively. The precision of the model was evaluated by observing the frequency of male-female encounters during 212 5-min observation periods. The observed encounter rate was 1.49/min, compared to a predicted rate of 1.53/ min. These values are not significantly different as determined by t-test. The constraints placed on the timing of sex by population density are discussed.     

Simulation of molecular crowding effects on an Alzheimer's α-amyloid peptide

Fibril formation by the Alzheimer's β-amyloid (Aβ) peptide in brain tissue is integral to the Alzheimer's disease pathology. Understanding the conformational properties and the mechanisms triggering aggregation of the Aβ peptides, at an atomic level of detail, is of crucial importance for the design of effective therapeutic agents against this disease. In this work, the conformational transitions and dynamic properties of an amyloidogenic peptide fragment (Aβ_10-35) were studied by molecular dynamics simulations in systems modeling infinite dilution and the presence of macromolecular crowding agents (CA). The model system consists of the peptide described with an atomistic force field, the CA represented by inert, quasi-hard spheres and a continuum solvent model. This combined model allowed the simulations to be extended to 100 ns each. Simulations were carried out starting from a completely extended structure, a β-strand structure, and four nuclear magnetic resonance structures in dilute aqueous solution. For all structures, two additional simulations were performed that included the inert CA in the solution and occupied approx 30 and 40% of the volume, respectively. For two of the nuclear magnetic resonance structures, additional simulations were carried out with 35% volume fraction of CA to further examine the diffusive behavior of the peptide. The peptide adopted a collapsed coil conformation in all simulations. The results of the simulations in dilute solution showed reasonable qualitative agreement with experimental and other simulation results, whereas the presence of volume excluding agents resulted in some distinct changes in properties (e.g., an increase in the appearance of transient β-structure or decreases in diffusivity with increasing CA concentration). At the same time, internal motion such as order parameter or atomic root mean square fluctuations showed less systematic responses to volume exclusion.     

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

Effects of solution crowding on actin polymerization reveal the energetic basis for nucleotide-dependent filament stability

Actin polymerization is a fundamental cellular process involved in cell structure maintenance, force generation, and motility. Phosphate release from filament subunits following ATP hydrolysis destabilizes the filament lattice and increases the critical concentration (C_c) for assembly. The structural differences between ATP- and ADP-actin are still debated, as well as the energetic factors that underlie nucleotide-dependent filament stability, particularly under crowded intracellular conditions. Here, we investigate the effect of crowding agents on ATP- and ADP-actin polymerization and find that ATP-actin polymerization is largely unaffected by solution crowding, while crowding agents lower the C_c of ADP-actin in a concentration-dependent manner. The stabilities of ATP- and ADP-actin filaments are comparable in the presence of physiological amounts (∼ 30% w/v) and types (sorbitol) of low molecular weight crowding agents. Crowding agents act to stabilize ADP-F-actin by slowing subunit dissociation. These observations suggest that nucleotide hydrolysis and phosphate release per se do not introduce intrinsic differences in the in vivo filament stability. Rather, the preferential disassembly of ADP-actin filaments in cells is driven through interactions with regulatory proteins. Interpretation of the experimental data according to osmotic stress theory implicates water as an allosteric regulator of actin activity and hydration as the molecular basis for nucleotide-dependent filament stability.     

Quantitative measures of femoral fracture repair in rats derived by micro-computed tomography

Although fracture healing is frequently studied in pre-clinical models of long bone fractures using rodents, there is a dearth of objective quantitative techniques to assess successful healing. Biomechanical testing is possibly the most quantitative and relevant to a successful clinical outcome, but it is a destructive technique providing little insight into the cellular mechanisms associated with healing. The advent of X-ray computed tomography (CT) has provided the opportunity to quantitatively and non-destructively assess bone structure and density, but it is unknown how measurements derived using this technology relate to successful healing. To examine possible relationships, we used a pre-clinical model to test for statistically significant correlations between quantitative characteristics of the callus by micro-CT (μCT) and the bending strength, stiffness, and energy-to-failure of the callus as assessed by three-point bending of excised bones. A closed, transverse fracture was generated in the mid-shaft of rat femurs by impact loading. Shortly thereafter, the rats received a one-time, local injection of either the vehicle or one of four doses of lovastatin. Following sacrifice after 4 weeks of healing, fractured femurs were extracted for μCT analysis and then three-point bending. Setting the region of interest to be 3.2 mm above and below the fracture line, we acquired standard and new μCT-derived measurements. The mineralized callus volume and the mineral density of the callus correlated positively with callus strength (r_xy=?0.315, p=0.016 and r_xy=0.444, p<0.0005, respectively) and stiffness (r_xy=?0.271, p=0.040 and r_xy=0.325, p=0.013, respectively), but the fraction of the callus that mineralized and the moment of inertia of the callus did not. This fraction did correlate with energy-to-failure (r_xy=?0.343, p=0.0085). Of the μCT-derived measurements, quantifying defects within the outer bridging cortices of the callus produced the strongest correlation with both callus strength (r_xy=0.557, p<0.0001) and stiffness (r_xy=0.468, p=0.0002). By both reducing structural defects and increasing mineralization, lovastatin appears to increase the callus strength.     

Inheritance of Traits Controlled by Odd Chromosomes Using Data on Transmission of Monosomic Addition S<Superscript>t</Superscript> Chromosome of the <Emphasis Type="Italic">Elymus Sibiricus</Emphasis> Genome

On the basis of the winter bread wheat cultivar Obryi, two independent disomic addition lines BC₁₂F_∞ with the chromosome of the E. sibiricus S^t genome are created. A practical algorithm for determining the probabilities of transmission of the odd chromosome separately through male and female gametes in selfpollination of hemizygous hybrids from the equation p²–(1 + f₁–f₄) × p + f₁ = 0 is proposed, where p is the probability of the formation of viable gametes with the considered chromosome and f₁ and f₄ are the empirical frequencies of the corresponding homozygotes with and without the trait. The probability of transmission of an alien univalent chromosome through pollen (p_♂) is associated with the frequency of its transmission through the egg cell (p_♀) in backcrosses and in self-pollination (1–f₄) by the equation p_♂ = 1–f₄/(1–p_♀). The calculated empirically dependent estimates of the probabilities of transmission of the added chromosome through the egg cell p_♀ = 18.7% and through pollen p_♂ = 4.3% correspond to the empirical frequencies obtained for backcrosses. The coefficients of the gamete selection V_♀ = 0.748 and V_♂ = 0.172 are calculated, and the expected segregation for the alien trait controlled by a dominant gene located in the added chromosome is determined—with the trait: without the trait is 0.222: 0.778 in F₂; 0.187: 0.813 in equational and 0.043: 0.957 in certational backcrosses.     

Variations of Dynamic Contrast-Enhanced Magnetic Resonance Imaging in Evaluation of Breast Cancer Therapy Response: A Multicenter Data Analysis Challenge

Pharmacokinetic analysis of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) time-course data allows estimation of quantitative parameters such as K^trans (rate constant for plasma/interstitium contrast agent transfer), v_e (extravascular extracellular volume fraction), and v_p (plasma volume fraction). A plethora of factors in DCE-MRI data acquisition and analysis can affect accuracy and precision of these parameters and, consequently, the utility of quantitative DCE-MRI for assessing therapy response. In this multicenter data analysis challenge, DCE-MRI data acquired at one center from 10 patients with breast cancer before and after the first cycle of neoadjuvant chemotherapy were shared and processed with 12 software tools based on the Tofts model (TM), extended TM, and Shutter-Speed model. Inputs of tumor region of interest definition, pre-contrast T₁, and arterial input function were controlled to focus on the variations in parameter value and response prediction capability caused by differences in models and associated algorithms. Considerable parameter variations were observed with the within-subject coefficient of variation (wCV) values for K^trans and v_p being as high as 0.59 and 0.82, respectively. Parameter agreement improved when only algorithms based on the same model were compared, e.g., the K^trans intraclass correlation coefficient increased to as high as 0.84. Agreement in parameter percentage change was much better than that in absolute parameter value, e.g., the pairwise concordance correlation coefficient improved from 0.047 (for K^trans) to 0.92 (for K^trans percentage change) in comparing two TM algorithms. Nearly all algorithms provided good to excellent (univariate logistic regression c-statistic value ranging from 0.8 to 1.0) early prediction of therapy response using the metrics of mean tumor K^trans and k_ep (= K^trans/v_e, intravasation rate constant) after the first therapy cycle and the corresponding percentage changes. The results suggest that the interalgorithm parameter variations are largely systematic, which are not likely to significantly affect the utility of DCE-MRI for assessment of therapy response.     

The Influence of Crowding Conditions on the Thermodynamic Feasibility of Metabolic Pathways

Intracellular reactions are carried out in a crowded medium where the macromolecules occupy ∼40% of the total volume. This decrease in the available volume affects the activity of the reactants. Scaled particle theory is used for the estimation of the activity coefficients of the metabolites, and thereby for the assessment of the impact of the presence of background molecules, on the estimation of the Gibbs free energy change (Δ_rG) of the reactions. The lactic acid pathway and the central carbon metabolism of Actinobacillus succinogenes for the production of succinic acid from glycerol have been used as illustrative case studies. Results suggest the importance of maintaining intracellular crowded regions to favor the feasibility of a pathway that in other circumstances would be infeasible. Moreover, the crowding conditions may change the directionality of reactions and can modify the feasible range of fluxes estimated for a metabolic system compared with those obtained at standard biological conditions.