Implicit solvent models

Implicit solvent models for biomolecular simulations are reviewed and their underlying statistical mechanical basis is discussed. The fundamental quantity that implicit models seek to approximate is the solute potential of mean force, which determines the statistical weight of solute conformations, and which is obtained by averaging over the solvent degrees of freedom. It is possible to express the total free energy as the reversible work performed in two successive steps. First, the solute is inserted in the solvent with zero atomic partial charges; second, the atomic partial charges of the solute are switched from zero to their full values. Consequently, the total solvation free energy corresponds to a sum of non-polar and electrostatic contributions. These two contributions are often approximated by simple geometrical models (such as solvent exposed area models) and by macroscopic continuum electrostatics, respectively. One powerful route is to approximate the average solvent density distribution around the solute, i.e. the solute-solvent density correlation functions, as in statistical mechanical integral equations. Recent progress with semi-analytical approximations makes continuum electrostatics treatments very efficient. Still more efficient are fully empirical, knowledge-based models, whose relation to explicit solvent treatments is not fully resolved, however. Continuum models that treat both solute and solvent as dielectric continua are also discussed, and the relation between the solute fluctuations and its macroscopic dielectric constant(s) clarified.     

Protein–solvent interaction

Protein folding and assembly can be manipulated in in vitro systems by co-solvents at high concentrations. A number of co-solvents that enhance protein stability and assembly have been shown to be excluded from the protein surface. Such co-solvent exclusion has been demonstrated by dialysis experiments and shown to be correlated with their effects on protein stability and assembly.     

Preferential and absolute interactions of solvent components with proteins in mixed solvent systems

An equation is derived relating total and preferential interactions of solvent components with macromolecules in a three-component system. Application of this equation to literature data shows that binding of non-polar solvents to proteins parallels the unfolding of the latter. This observation is discussed in terms of local inter-residue and residue-solvent interactions.     

Lability of chlorophylls in solvent

Journal of Applied Phycology - DMSO (dimethyl sulfoxide, (CH3)2SO) is an alternative solvent for the spectroscopic assay of chlorophylls (Chl) but has mainly been used on Chl a and b organisms, but...     

The thermodynamics of solvent exchange

A model for solvation in mixed solvents, which was developed for the free energy and preferential interaction [J. A. Schellman (1987), Biopolymers, Vol. 26, pp. 549–559; (1990), Biophysical Chemistry, Vol. 37, pp. 121–140; (1993), Biophysical Chemistry, Vol. 45, pp. 273–279], is extended in this paper to cover the thermal properties: enthalpy, entropy, and heat capacity. An important result is that the enthalpy of solvation H? responds directly to the fraction of site occupation. This differs from the free energy ? and preferential interaction Γ₃₂, which are measures of the excess binding above a random distribution of solvent molecules. In other words, the enthalpy is governed by K while ? and Γ₃₂ are governed by (K ? 1) where K is the equilibrium constant on a mole fraction scale [Schellman (1987)]. The solvation heat capacity C?p consists of two term: (1) the intrinsic heat capacity of species in solution with no change in composition, and (2) a term that accounts for the change in composition that accompanies solvent exchange. Binding to biological macromolecules is heterogeneous but experiementalists must use binding isotherms that assume the homogeneity of sites. Equations are developed for the interpretation of the experimental parameters (number of sites n_exp, equilibrium constant K_exp, and enthalpy, Δh_exp), when homogeneous formulas are applied to the heterogeneous case. It is shown that the experimental parameters for the occupation and enthalpy are simple functions of the moments of the distribution of equilibrium constants over the sites. In general, n_exp is greater than the true number of sites and K_exp is greater than the average of the equilibrium constants. The free energy and preferential interaction can be fit to a homogenious formula, but the parameters of the curve are not easily represented in terms of the moments of distributions over the sites. The strengths and deficiencies of this type of thermodynamic model are discussed. © 1994 John Wiley & Sons, Inc.     

Fifty years of solvent denaturation

The author has worked in the area of solvent denaturation and stabilization, off and on, for approximately 50 years. This paper is a personal view of the progress which has been made since 1950. The topic is limited to the development of thermodynamic molecular models for the interpretation of the unfolding and stabilization of protein structures. The story starts with the work in Kauzmann's laboratory shortly after World War II and proceeds through models for multisite binding, the linear denaturation curve, the special considerations required for understanding weak solvent exchange and a new model for the balance of solvent contact interactions and excluded volume.     

Computer-aided solvent screening for biocatalysis

A computer-aided solvent screening methodology is described and tested for biocatalytic systems composed of enzyme, essential water and substrates/products dissolved in a solvent medium, without cells. The methodology is computationally simple, using group contribution methods for calculating constrained properties related to chemical reaction equilibrium, substrate and product solubility, water solubility, boiling points, toxicity and others. Two examples are provided, covering the screening of solvents for lipase-catalyzed transesterification of octanol and inulin with vinyl laurate. Esterification of acrylic acid with octanol is also addressed. Solvents are screened and candidates identified, confirming existing experimental results. Although the examples involve lipases, the method is quite general, so there seems to be no preclusion against application to other biocatalysts.     

Listeria monocytogenes is a solvent tolerant organism secreting a solvent stable lipase: potential biotechnological applications

Biotechnology Letters - The emerging biobased economy will require robust, adaptable, organisms for the production and processing of biomaterials as well as for bioremediation. Recently, the search...     

Organic solvent tolerance of halophilic archaea

Organic solvent tolerance was tested in type strains of type species of the sixteen genera of Halobacteriaceae, the halophilic archaea. Most of the strains were observed to grow in the presence of hexylether (log Pow=5.1), but none grew in the presence of n-octane (log Pow=4.9) except Halogeometricum borinquense JCM 10706T and Halorubrum saccharovorum JCM 8865T. On the other hand, two strains, Haloarcula spp. OHF-1 and 2 isolated from a French solar salt were found to show stronger tolerance even to isooctane (log Pow=4.8). Growth of some strains was retarded by the presence of n-decane but reached to the same cell densities at late stationary phase. Final cell densities of some strains were greatly repressed by the presence of the solvent.     

On protein solubility in organic solvent

Solubility of a model protein, hen egg-white lysozyme, was investigated in a wide range of neat nonaqueous solvents and binary mixtures thereof. All solvents that are protic, very hydrophilic, and polar readily dissolve more than 10 mg/mL of lysozyme (lyophilized from aqueous solution of pH 6.0). Only a marginal correlation was found between the lysozyme solubility in a non-aqueous solvent and the letter's dielectric constant or Hildebrand solubility parameter, and no correlation was observed with the dipole moment. Lysozyme dissolved in dimethyl sulfoxide (DMSO) could be precipitated by adding protein nondissolving co-solvents, although the enzyme had a tendency to form supersaturated solutions in such mixtures. The solubility of lysozyme, both in an individual solvent (1,5-pentanediol) and in binary solvent mixtures (DMSO/acetonitrile), markedly increased when the pH of the enzyme aqueous solution prior to lyophilization was moved away from the proteins's isoelectric point. (c) 1994 John Wiley & Sons, Inc.     

PredAcc: prediction of solvent accessibility

SUMMARY: PredAcc is a tool for predicting the solvent accessibility of protein residues from the sequence at different relative accessibility levels (0-55%). The prediction rate varies between 70. 7% (for 25% relative accessibility) and 85.7% (for 0% relative accessibility). Amino acids are predicted in four categories: almost certainly hidden and almost certainly exposed with a given a posteriori prediction error, probably hidden and probably exposed otherwise. AVAILABILITY: http://condor.urbb.jussieu.fr/PredAccCfg.html CONTACT: tuffery@urbb.jussieu.fr     

Microbial ecology of chlorinated solvent biodegradation

This study focused on the microbial ecology of tetrachloroethene (PCE) degradation to trichloroethene, cis‐1,2‐dichloroethene and vinyl chloride to evaluate the relationship between the microbial community and the potential accumulation or degradation of these toxic metabolites. Multiple soil microcosms supplied with different organic substrates were artificially contaminated with PCE. A thymidine analogue, bromodeoxyuridine (BrdU), was added to the microcosms and incorporated into the DNA of actively replicating cells. We compared the total and active bacterial communities during the 50‐day incubations by using phylogenic microarrays and 454 pyrosequencing to identify microorganisms and functional genes associated with PCE degradation to ethene. By use of this integrative approach, both the key community members and the ecological functions concomitant with complete PCE degradation could be determined, including the presence and activity of microbial community members responsible for producing hydrogen and acetate, which are critical for Dehalococcoides‐mediated PCE degradation. In addition, by correlation of chemical data and phylogenic microarray data, we identified several bacteria that could potentially oxidize hydrogen. These results demonstrate that PCE degradation is dependent on some microbial community members for production of appropriate metabolites, while other members of the community compete for hydrogen in soil at low redox potentials.     

Implicit solvent simulation models for biomembranes

Fully atomic simulation strategies are infeasible for the study of many processes of interest to membrane biology, biophysics and biochemistry. We review various coarse-grained simulation methodologies with special emphasis on methods and models that do not require the explicit simulation of water. Examples from our own research demonstrate that such models have potential for simulating a variety of biologically relevant phenomena at the membrane surface.     

Selective solvent delignification for fermentation enhancement

Cellulose and hemicellulose in renewable biomass resources such as cornstover and wheat straw have been examined as substrates for the production of ethanol. A mixed culture of selected strains of Clostridium thermocellum and Clostridium thermosaccharolyticum are used to accomplish both the hydrolysis and fermentation of these carbohydrates in a single step. However, lignin and related phenolic materials are shown to diminish the rate, extent, and yield at which these carbohydrates can be utilized for ethanol production. In order to overcome this problem, a selective solvent pretreatment with alkaline-ethanol-water mixtures was examined for the delignification of cellulosic biomass under conditions where very little loss of fermentable carbohyrates results. Under optimal conditions, up to 67% of the initial lignin in cornstover can be extracted while 95% of the alpha-cellulose and pentosan carbohydrates remain insoluble. Subsequent mixed culture fermentation of the treated material has shown a 400% increase in the rate of degradation and greater than 85% utilization of the substrate. The effects of various extraction parameters on delignification kinetics and subsequent fermentation performance are discussed.