Uncovering the basis for nonideal behavior of biological molecules

The molecular origin of the nonideal behavior for concentrated binary solutions of biochemical compounds is examined. The difference between activities expressed in the molar and molal conventions can be large. Considering the range from dilute to concentrated, we show that molar activity coefficients can be represented by simple but rigorous equations involving between one and three parameters only. We derive a universal relationship interconverting the scales of molarity and molality without requiring the density of the solution. The equations are developed from first principles using a statistical thermodynamic theory of molar activity coefficients. It is shown how to express activity coefficients in different concentration scales, and the advantages and disadvantages of using certain scales are discussed and compared with the experimental data. Several classes of biochemically relevant compounds, many of which are naturally occurring osmolytes, are discussed: six saccharides (glucose, xylose, maltose, mannose, raffinose, and sucrose), four polyols (glycerol, mannitol, erythritol, and sorbitol), five amino acids (glycine, alanine, sarcosine, glycine betaine, and proline), and urea. Of the 16 solutes, 10 could be described in terms of a single parameter that is due to pure first-order effects (packing, hydration, or space limitation). The remaining six exhibit significant second-order effects (solute-solute interactions) and require two additional parameters, one typically identified with the volume occupied per solute molecule in the pure solute (crystal or liquid) and the other with a self-association constant. The activity coefficients of the osmolytes roughly display the rank order found with respect to their ability to stabilize proteins. These findings are discussed in terms of the physical principles that give rise to the activity coefficients.     

Macromolecular components of the vitelline membrane of hen's egg. I. Membrane structure and its deterioration with age.

The vitelline membrane of hen's egg has been successfully solubilized in sodium dodecyl sulfate (SDS), guanidine hydrochloride and urea solutions, and its macromolecular components examined. SDS-gel electrophoresis of the membrane solution revealed the presence of three major components designated I, II, and III, all containing carbohydrate and protein. The approximate molecular weights of components I and II were 32,000 and 260,000, respectively, and the sedimentation coefficients were 2.2S and 4.3S. Component III was in an aggregated form which disintegrated into smaller components upon reduction with 2-mercaptoethanol. It was found that component II (4.3S component) deteriorated during storage of the egg with the concomitant formation of degraded components. The loss of this component was accompanied by a gradual decrease of the neutral sugar content of the vitelline membrane. On the basis of these data, the membrane structure and its deterioration during storage are discussed.     

A comparison of three models for the diffusion of oxygen in electrolyte solutions

Diffusion of oxygen through aqueous solutions is of great importance in biological systems. In this work, three models for the diffusion of oxygen through aqueous salt solutions are compared. One model uses mole fraction as the driving force (Fick's Law) and another uses chemical potential. The third model uses the gradient in oxygen activity as the driving force. This new model was chosen because of the availability of oxygen electrodes which directly measure oxygen activity in aqueous solution. These models have been used to reevaluate the technique of measuring O(2) diffusivities. We show that Pick's Law diffusion coefficients do not vary strongly with salt concentration as was erroneously reported in the literature. In addition, we compare the predicted O(2) fluxes of the three models over a wide range in O(2) concentrations. For oxygen concentrations of biological interest, the three models give identical predictions of the flux.     

Urea effect on Cu(2+)-induced DNA structural transitions in solution

Cu(2+) ion interaction with DNA in aqueous solutions containing urea (0-5 M) was studied by IR spectroscopy. It was shown that upon the Cu(2+) ion binding DNA transition into a compact form occurs. This transition is of positive cooperativity. We suppose that the mechanism of Cu(2+)-induced DNA compaction in solutions containing urea is not completely electrostatic. Urea addition to the DNA solution decreases the Cu(2+) ion concentration required to induce DNA compaction. As the urea content in solution rises, the binding constant of Cu(2+) ions interacting with DNA increases, going through the maximum in the case of 2 M solution; further increase of the urea content in solutions leads to decrease of the binding constant. DNA transition into the compact form under the Cu(2+) ion action is determined not only by the effects of the solution dielectric permeability but by the solvation effects; when changes of the dielectric permeability are small the solvation effects may prevail. Urea addition to the DNA solution also decreases cooperativity of the DNA compaction process. Perhaps, cooperativity of the DNA transition into the compact state depends on the ordered spatial structure of water adjacent to the macromolecule and decreases on the structure destruction.     

Preferential interactions in aqueous solutions of urea and KCl

A quantitative characterization of the thermodynamic effects due to interactions of salt ions and urea in aqueous solution is needed for rigorous analyses of the effects of changing urea concentration on biopolymer processes in solutions that also contain salt. Therefore, we investigate preferential interactions in aqueous solutions containing KCl and urea by using vapor pressure osmometry (VPO) to measure osmolality as a function of the molality of urea (component 3) over the range 0.09相似文献   

Inhibition of protein carbamylation in urea solution using ammonium-containing buffers

Urea solution is one of the most commonly employed protein denaturants for protease digestion in proteomic studies. However, it has long been recognized that urea solution can cause carbamylation at the N termini of proteins/peptides and at the side chain amino groups of lysine and arginine residues. Protein/peptide carbamylation blocks protease digestion and affects protein identification and quantification in mass spectrometry analysis by blocking peptide amino groups from isotopic/isobaric labeling and changing peptide charge states, retention times, and masses. In addition, protein carbamylation during sample preparation makes it difficult to study in vivo protein carbamylation. In this study, we compared the peptide carbamylation in urea solutions of different buffers and found that ammonium-containing buffers were the most effective buffers to inhibit protein carbamylation in urea solution. The possible mechanism of carbamylation inhibition by ammonium-containing buffers is discussed, and a revised procedure for the protease digestion of proteins in urea and ammonium-containing buffers was developed to facilitate its application in proteomic research.     

The effect of solution non-ideality on membrane transport in three-dimensional models of the renal concentrating mechanism

Previous models of the renal concentrating mechanism employ ideal approximations of solution thermodynamics for membrane transport calculation. In three-dimensional models of the renal medulla, predicted urine concentrations reach levels where there idealized approximations begin to break down. In this paper we derive equations that govern membrane transport for non-dilute solutions and use these equations in a three-dimensional model of the concentrating mechanism. New numerical methods were employed that are more stable than those employed previously. Compared to ideal solution models, the urea non-ideality tends to increase predicted osmolarities, whereas NaCl non-ideality decreases predictions.     

Multistationarity in mass action networks with applications to ERK activation

Ordinary Differential Equations (ODEs) are an important tool in many areas of Quantitative Biology. For many ODE systems multistationarity (i.e. the existence of at least two positive steady states) is a desired feature. In general establishing multistationarity is a difficult task as realistic biological models are large in terms of states and (unknown) parameters and in most cases poorly parameterized (because of noisy measurement data of few components, a very small number of data points and only a limited number of repetitions). For mass action networks establishing multistationarity hence is equivalent to establishing the existence of at least two positive solutions of a large polynomial system with unknown coefficients. For mass action networks with certain structural properties, expressed in terms of the stoichiometric matrix and the reaction rate-exponent matrix, we present necessary and sufficient conditions for multistationarity that take the form of linear inequality systems. Solutions of these inequality systems define pairs of steady states and parameter values. We also present a sufficient condition to identify networks where the aforementioned conditions hold. To show the applicability of our results we analyse an ODE system that is defined by the mass action network describing the extracellular signal-regulated kinase (ERK) cascade (i.e. ERK-activation).     

Reversible thermal denaturation of immobilized rhodanese

For the first time, the enzyme rhodanese had been refolded after thermal denaturation. This was previously not possible because of the strong tendency for the soluble enzyme to aggregate at temperatures above 37 degrees C. The present work used rhodanese that was covalently coupled to a solid support under conditions that were found to preserve enzyme activity. Rhodanese was immobilized using an N-hydroxymalonimidyl derivative of Sepharose containing a 6-carbon spacer. The number of immobilized competent active sites was measured by using [35S]SO3(2-) to form an active site persulfide that is the obligatory catalytic intermediate. Soluble enzyme was irreversibly inactivated in 10 min at 52 degrees C. The immobilized enzyme regained at least 30% of its original activity even after boiling for 20 min. The immobilized enzyme had a Km and Vmax that were each approximately 3 times higher than the corresponding values for the native enzyme. After preincubation at high temperatures, progress curves for the immobilized enzyme showed induction periods of up to 5 min before attaining apparently linear steady states. The pH dependence of the activity was the same for both the soluble and the immobilized enzyme. These results indicate significant stabilization of rhodanese after immobilization, and instabilities caused by adventitious solution components are not the sole reasons for irreversibility of thermal denaturation seen with the soluble enzyme. The results are consistent with models for rhodanese that invoke protein association as a major cause of inactivation of the enzyme. Furthermore, the induction period in the progress curves is consistent with studies which show that rhodanese refolding proceeds through intermediate states.     

Refolding of a high molecular weight protein: salt effect on collapse

Small-angle neutron scattering experiments were performed on dilute solutions of a high molecular weight protein (fibronectin, M = 580 kg/mol) in four cases: native conditions; unfolded state obtained by a denaturing agent (urea); and two badly refolded (or collapsed) states obtained by progressive elimination of the denaturing agent in salt-containing or salt-free solutions. Our main result is concerned by the conformation of the protein as the attempt for refolding is driven with or without salt. In salt-containing solution, we observe unambiguously that the protein chain collapses at large length scales but still obeys to a Gaussian statistics at short length scales. In other words, the globule embodies a large quantity of solvent compared to the compact situation. In salt-free solutions, the badly refolded protein is not globular but displays both a coil-like and an open conformation at large length scales and a local high density area. This behavior is discussed with respect to the scaling theories for polymers and polyampholytes.     

Calorimetric study of the thermal denaturation of beta-lactoglobulin in the presence of urea and phosphate ions]

The denaturation of beta-lactoglobulin in solution with different content of urea and phosphates has been studied calorimetrically. It has been shown that the increase of phosphate ion concentration in solution leads to an increase of beta-lactoglobulin stability, while increase of urea concentration leads to an opposite effect. The variation of these components in solution practically does not influence the value of the heat capacity increment of beta-lactoglobulin in the considered temperature region. Accordingly the denaturation enthalpy is a linear function of temperature whose slope does not differ for solution with urea concentration less than 4.4 M. However, the absolute value of denaturation enthalpy in these solutions at corresponding temperatures differs significantly due to the heat effect of additional urea solvation during transition to the denatured state. The latter leads to a decrease of the overall denaturation enthalpy and, as a result, a shift of the enthalpy plot to higher temperatures providing conditions for studying the thermodynamic and structural characteristics of the molecule in the cold denatured-state.     

Effects of cell volume regulating osmolytes on glycerol 3-phosphate binding to triosephosphate isomerase

During cell volume regulation, intracellular concentration changes occur in both inorganic and organic osmolytes in order to balance the extracellular osmotic stress and maintain cell volume homeostasis. Generally, salt and urea increase the Km's of enzymes and trimethylamine N-oxide (TMAO) counteracts these effects by decreasing Km's. The hypothesis to account for these effects is that urea and salt shift the native state ensemble of the enzyme toward conformers that are substrate-binding incompetent (BI), while TMAO shifts the ensemble toward binding competent (BC) species. Km's are often complex assemblies of rate constants involving several elementary steps in catalysis, so to better understand osmolyte effects we have focused on a single elementary event, substrate binding. We test the conformational shift hypothesis by evaluating the effects of salt, urea, and TMAO on the mechanism of binding glycerol 3-phosphate, a substrate analogue, to yeast triosephosphate isomerase. Temperature-jump kinetic measurements promote a mechanism consistent with osmolyte-induced shifts in the [BI]/[BC] ratio of enzyme conformers. Importantly, salt significantly affects the binding constant through its effect on the activity coefficients of substrate, enzyme, and enzyme-substrate complex, and it is likely that TMAO and urea affect activity coefficients as well. Results indicate that the conformational shift hypothesis alone does not account for the effects of osmolytes on Km's.     

Water activity coefficients in aqueous amino acid solutions by molecular dynamics simulation: 1. Force field development

New force fields for molecular dynamics (MD) simulation of aqueous zwitterionic amino acid simulations were developed. These were especially designed to calculate activity coefficient of water in amino acid solutions with high accuracy. For example, aqueous solutions of the following amino acids were considered: glycine, alanine, α-aminobutyric acid, α-aminovalerianic acid, valine and leucine. The force fields were obtained by quantum chemical calculations using B3LYP/6-31G and MP2/6-311(d,p) model theories in combination with the Merz–Kollmann–Singh scheme. To further increase the accuracy of the force field, a polarised continuum was considered in all quantum chemical calculations. Water activity coefficients obtained from MD using different all-purpose literature force fields, namely, OPLS, AMBER ff03 and GROMOS 53A6 as well as experimental data are compared with the results utilising the new force field. The new force field is shown to give better results compared with experimental data than existing force fields.     

Effects of polyethylene glycol and dextrans on immunoprecipitations: a two-cross immunodiffusion study

Precipitating titers and immunochemical titers obtained in a wide range of antigen-to-antibody concentration ratios by the two-cross immunodiffusion technique are compared with the corresponding laser light scatter precipitin curves. The two-cross immunodiffusion technique has also been applied to investigate whether polyethylene glycol of molecular mass 6000 and dextrans of molecular masses from 10,000 to 2,000,000 enhance the immunoprecipitation processes of the system human serum IgG-rabbit immune serum at pH 5.5 and 8.1 at 20 degrees C. It was found that the significant increase of precipitating titers of both precipitating components in the presence of polyethylene glycol is a consequence of a strong decrease of solubility of the primary antigen-antibody complex. The decrease of solubility does not affect the immunochemical titer of the immune serum, indicating stoichiometrical invariance of the precipitate at the equivalence. The apparent strong decrease of diffusion coefficients of both antigen and antibody in 20- and 40-g/liter polyethylene glycol solution is attributed to increase of viscosity of the solutions and to a partial self-association of protein molecules due to steric exclusion. In 40-g/liter polyethylene glycol solutions at pH 5.5 every fourth molecular entity of antigen and every third molecular entity of antibody are present in the form of a two-molecular self-associate, whereas in 20-g/liter polyethylene glycol solutions only 1% of antigen molecules and 8% of antibody molecules are associated. With the increase of pH to 8.1 the self-association of protein molecules is strongly further enhanced. Dextrans in 20-g/liter solutions, without regard to their relative molecular masses, do not influence precipitating titers and solubility of the antigen-antibody system at equivalence and do not enhance self-association of protein molecules. The strong decrease of diffusion coefficients of immunoglobulin G antigen and antibodies in dextran solutions is solely attributed to the increase of viscosity of the dextran solutions; hence there was no evidence of interaction of dextrans with serum IgG proteins.     

Osmoregulation in red blood cells of Bufo viridis

The effect of hyperosmotic solution of NaCl, urea and mannitol on Bufo viridis red blood cells were studied. The percentage of water content in B. viridis red blood cells decreased significantly in NaCl and mannitol hypertonic solutions compared to urea hypertonic solution. The urea concentration found in red blood cells in a urea hypertonic solution was significantly higher than in red blood cells acclimated to NaCl and mannitol hypertonic solutions. The Na+ concentration was significantly lower in red blood cells immersed in urea hypertonic solution than in red blood cells immersed in hypertonic NaCl and mannitol solutions. However, the K+ concentration increased at a similar rate in three different hypertonic solutions.     

Validation of a handheld display device for an expandable labeled magnitude scale (LMS)

Palmtop computers provide a possible avenue for the convenient collection of subjective ratings from individuals outside of a fixed laboratory setting. One disadvantage of these computers is the small size of their display screens, which may reduce the resolution of responses available as compared with standard display screens. One plausible solution to this problem is to use a scale that expands contingent on an initial response made by the subject, so that the final response is made from a scale with finer resolution. To validate this approach, we compared taste intensity judgments of six sucrose solutions (0.03, 0.06, 0.12, 0.24, 0.48, and 0.96 M), using a labeled magnitude scale either presented in expandable form on a palmtop computer (Palm scale) or in conventional (nonexpandable) form on a standard 17' PC monitor (PC scale). Twenty-four subjects rated all six sucrose solutions thrice, using both scale types, the different scales being used on different days of testing. The scales led to very similar taste intensity ratings at all but the lowest concentration, which was rated less intense on the Palm scale. The Palm scale was used with slightly less precision than the PC scale for the weakest solution concentrations. In summary, the responses of the two scales were similar enough to validate the use of the expandable scale on the palmtop computer outside the laboratory setting.     

A kinetic description of antifreeze glycoprotein activity

The antifreeze glycoproteins (AFGP) of polar fish have the ability to depress the freezing temperature of water approximately 500 times the amount expected based on the number of AFGP molecules in solution; yet AFGP solutions have a purely colligative melting point depression. The difference of solution melting and freezing temperatures is the antifreeze activity of AFGP. One characteristic of AFGP activity that requires further examination is the effect of concentration on antifreeze activity, especially whether the activity saturates at high concentrations or the measured activity increases ad infinitum. This study first surveys the activity of the various antifreeze components from both Pagothenia borchgrevinki and the Arg-containing antifreeze glycoprotein from Eleginus gracilis (EgAF). It was found that all AFGP components examined have a plateau in activity at high concentration, but the actual value of the plateau activity differs between the different length AFGP components and between AFGP and EgAF. While the low molecular weight components of both AFGP and EgAF lose activity at deep supercooling, at high concentration activity is restored. The activity data is then shown to fit a reversible kinetic model of AFGP activity, and the coefficients obtained are used to compare the activity differences between AFGP components and between AFGP and EgAF. The model is also shown to describe the activity of the antifreeze protein of the fish Pseudopleuronectes americanus and the thermal hysteresis protein of the insect, Tenebrio molitor.     

Comparison of Guanidine Hydrochloride (GdnHCl) and Urea Denaturation on Inactivation and Unfolding of Human Placental Cystatin (HPC)

The activity and conformational change of human placental cystatin (HPC), a low molecular weight thiol proteinase inhibitor (12,500) has been investigated in presence of guanidine hydrochloride (GdnHCl) and urea. The denaturation of HPC was followed by activity measurements, fluorescence spectroscopy and Circular Dichroism (CD) studies. Increasing the denaturant concentration significantly enhanced the inactivation and unfolding of HPC. The enzyme was 50% inactivated at 1.5 M GdnHCl or 3 M urea. Up to 1.5 M GdnHCl concentration there was quenching of fluorescence intensity compared to native form however at 2 M concentration intensity increased and emission maxima had 5 nm red shift with complete unfolding in 4–6 M range. The mid point of transition was in the region of 1.5–2 M. In case of urea denaturation, the fluorescence intensity increased gradually with increase in the concentration of denaturant. The protein unfolded completely in 6–8 M concentration of urea with a mid-point of transition at 3 M. CD spectroscopy shows that the ellipticity of HPC has increased compared to that of native up to 1.5 M GdnHCl and then there is gradual decrease in ellipticity from 2 to 5 M concentration. At 6 M GdnHCl the protein had random coil conformation. For urea the ellipticity decreases with increase in concentration showing a sigmoidal shaped transition curve with little change up to 1 M urea. The protein greatly loses its structure at 6 M urea and at 8 M it is a random coil. The urea induced denaturation follows two-state rule in which Native→Denatured state transition occurs in a single step whereas in case of GdnHCl, intermediates or non-native states are observed at lower concentrations of denaturant. These intermediate states are possibly due to stabilizing properties of guanidine cation (Gdn⁺) at lower concentrations, whereas at higher concentrations it acts as a classical denaturant.     

Nonlinear sorting, curvature generation, and crowding of endophilin N-BAR on tubular membranes

The curvature of biological membranes is controlled by membrane-bound proteins. For example, during endocytosis, the sorting of membrane components, vesicle budding, and fission from the plasma membrane are mediated by adaptor and accessory proteins. Endophilin is a peripherally binding membrane protein that functions as an endocytic accessory protein. Endophilin's membrane tubulation capacity is well known. However, to understand the thermodynamic and mechanical aspects of endophilin function, experimental measurements need to be compared to quantitative theoretical models. We present measurements of curvature sorting and curvature generation of the endophilin A1 N-BAR domain on tubular membranes pulled from giant unilamellar vesicles. At low concentration, endophilin functions primarily as a membrane curvature sensor; at high concentrations, it also generates curvature. We determine the spontaneous curvature induced by endophilin and observe sigmoidal curvature/composition coupling isotherms that saturate at high membrane tensions and protein solution concentrations. The observation of saturation is supported by a strong dependence of lateral diffusion coefficients on protein density on the tether membrane. We develop a nonlinear curvature/composition coupling model that captures our experimental observations. Our model predicts a curvature-induced phase transition among two states with varying protein density and membrane curvature. This transition could act as a switch during endocytosis.     

A simple model for solvation in mixed solvents. Applications to the stabilization and destabilization of macromolecular structures

The properties of a simple model for solvation in mixed solvents are explored in this paper. The model is based on the supposition that solvent replacement is a simple one-for-one substitution reaction at macromolecular sites which are independent of one another. This leads to a new form for the binding polynomial in which all terms are associated with ligand interchange rather than ligand addition. The principal solvent acts as one of the ligands. Thermodynamic analysis then shows that thermodynamic binding (i.e., selective interaction) depends on the properties of K'-1, whereas stoichiometric binding (site occupation) depends on K'. K' is a 'practical' interchange equilibrium constant given by (f3/f1)K, where K is the true equilibrium constant for the interchange of components 3 and 1 on the site and f3 and f4 denote their respective activity coefficients on the mole fraction scale. Values of K' less than unity lead to negative selective interaction. It is selective interaction and not occupation number which determines the thermodynamic effects of solvation. When K' greater than 100 on the mole fraction scale or K' greater than 2 on the molality scale (in water), the differences between stoichiometric binding and selective interaction become less than 1%. The theory of this paper is therefore necessary only for very weak binding constants. When K'-1 is small, large concentrations of the added solvent component are required to produce a thermodynamic effect. Under these circumstances the isotherms for the selective interaction and for the excess (or transfer) free energy are strongly dependent on the behavior of the activity coefficients of both solvent components. Two classes of behavior are described depending on whether the components display positive or negative deviations from Raoult's law. Examples which are discussed are aqueous solutions of urea and guanidinium chloride for positive deviations and of sucrose and glucose for negative deviations. Examination of the few studies which have been reported in the literature shows that most of the qualitative features of the stabilization of proteins by sugars and their destabilization by urea and guanidinium chloride are faithfully represented with the model. This includes maxima in the free energy of stabilization and destabilization, decreased and zero selective interaction at high concentrations, etc. These phenomena had no prior explanation. Deficiencies in the model as a representation of solvation in aqueous solution are discussed in the appendix.