Acid-induced partly folded conformation resembling a molten globule state of xylanase from an alkalothermophilic Bacillus sp.

Nonnative protein structures having a compact secondary, but not rigid tertiary structure, have been increasingly observed as intermediate states in protein folding. We have shown for the first time during acid-induced unfolding of xylanase (Xyl II) the presence of a partially structured intermediate form resembling a molten globule state. The conformation and stability of Xyl II at acidic pH was investigated by equilibrium unfolding methods. Using intrinsic fluorescence and CD spectroscopic studies, we have established that Xyl II at pH 1.8 (A-state) retains the helical secondary structure of the native protein at pH 7.0, while the tertiary interactions are much weaker. At variance, from the native species (N-state), Xyl II in the A-state binds 1-anilino-8-sulfonic acid (ANS) indicating a considerable exposure of aromatic side chains. Lower concentration of Gdn HCl are required to unfold the A-state. For denaturation by Gdn HCl, the midpoint of the cooperative unfolding transition measured by fluorescence for the N-state is 3.5 +/- 0.1 M, which is higher than the value (2.2 +/- 0.1 M) observed for the A-state at pH 1.8. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its structural stability that is characteristic for native proteins.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Molecular motion of spin labeled side chains in the C-terminal domain of RGL2 protein: a SDSL-EPR and MD study

Five singly spin labeled side chains at surface sites in the C-terminal domain of RGL2 protein have been analyzed to investigate the general relationship between nitroxide side chain mobility and protein structure. At these sites, the structural perturbation produced by replacement of a native residue with a nitroxide side chain appears to be very slight at the level of the backbone fold. The primary determinants of the nitroxide side chain mobility are backbone dynamics and tertiary interactions. On the exposed surfaces of alpha-helices, the side chain mobility is not restricted by tertiary interactions but appears to be determined by backbone dynamics, while in loop sites, the side chain mobility is even higher. For a better understanding of the changes in the EPR spectral line shape, molecular dynamics simulations were performed and found in agreement with EPR spectral data.     

Molecular dynamics study of the dominant-negative E219K polymorphism in human prion protein

Human prion diseases are associated with misfolding or aggregation of the Human Prion Protein (HuPrP). Missense mutations in the HuPrP gene, contribute to conversion of HuPrP^C to HuPrP^Sc and amyloid formation. Based on our previous comprehensive study, three missense mutations, from two different functional groups, i.e. disease-related mutations, and protective mutations, were selected and extensive molecular dynamics simulations were performed on these three mutants to compare their dynamics and conformations with those of the wildtype HuPrP. In addition to simulations of monomeric forms of mutants, in order to study the dominant-negative effect of protective mutation (E219K), 30-ns simulations were performed on E219K-wildtype and wildtype-wildtype dimeric forms. Our results indicate that, although after 30-ns simulations the global three-dimensional structure of models remain fairly intact, the disease-related mutations (V210I and Q212P) introduce local structural changes, i.e. close contact changes and secondary structure changes, in addition to global flexibility changes. Furthermore, our results support the loss of hydrophobic interaction due to the mutations in hydrophobic core that has been reported by previous NMR and computational studies. On the other hand, this protective mutation (E219K) results in helix elongation, and significant increases of overall flexibility of E219K mutant during 30-ns simulation. In conclusion, the simulations of dimeric forms suggest that the dominant-negative effect of this protective mutation (E219K) is due to the incompatible structures and dynamics of allelic variants during conversion process.     

The equilibrium unfolding of Azotobacter vinelandii apoflavodoxin II occurs via a relatively stable folding intermediate.

A flavodoxin from Azotobacter vinelandii is chosen as a model system to study the folding of alpha/beta doubly wound proteins. The guanidinium hydrochloride induced unfolding of apoflavodoxin is demonstrated to be reversible. Apoflavodoxin thus can fold in the absence of the FMN cofactor. The unfolding curves obtained for wild-type, C69A and C69S apoflavodoxin as monitored by circular dichroism and fluorescence spectroscopy do not coincide. Apoflavodoxin unfolding occurs therefore not via a simple two-state mechanism. The experimental data can be described by a three-state mechanism of apoflavodoxin equilibrium unfolding in which a relatively stable intermediate is involved. The intermediate species lacks the characteristic tertiary structure of native apoflavodoxin as deduced from fluorescence spectroscopy, but has significant secondary structure as inferred from circular dichroism spectroscopy. Both spectroscopic techniques show that thermally-induced unfolding of apoflavodoxin also proceeds through formation of a similar molten globule-like species. Thermal unfolding of apoflavodoxin is accompanied by anomalous circular dichroism characteristics: the negative ellipticity at 222 nM increases in the transition zone of unfolding. This effect is most likely attributable to changes in tertiary interactions of aromatic side chains upon protein unfolding. From the presented results and hydrogen/deuterium exchange data, a model for the equilibrium unfolding of apoflavodoxin is presented.     

Molecular dynamics simulations of human alpha-lactalbumin: changes to the structural and dynamical properties of the protein at low pH.

Two 700-ps molecular dynamics simulations of human alpha-lactalbumin have been compared. Both were initiated from an X-ray structure determined at pH 6.5. One simulation was designed to represent native conditions and the other the protein in solution at pH 2.0 without a bound calcium ion. The low pH conditions were modelled by protonating the aspartate, glutamate, and histidine side chains and the protein C-terminus. Significant changes were observed for the C-terminal region of the sequence in the simulation at low pH. Most notably an alpha-helix, helix D, and the C-terminal 3(10) helix were substantially disrupted relative to the simulation at high pH. These perturbations to the native fold are similar to those observed in an X-ray structure of alpha-lactalbumin at pH 4.2. In addition, larger fluctuations about side chain torsion angles were observed in the low pH simulation than in that corresponding to the higher pH. These structural and dynamical changes might be representative of the early stages of the transition to the molten-globule state of the protein known to be formed under low pH conditions in solution.     

The Effect of Electrostatics on the Marginal Cooperativity of an Ultrafast Folding Protein

Proteins fold up by coordinating the different segments of their polypeptide chain through a network of weak cooperative interactions. Such cooperativity results in unfolding curves that are typically sigmoidal. However, we still do not know what factors modulate folding cooperativity or the minimal amount that ensures folding into specific three-dimensional structures. Here, we address these issues on BBL, a small helical protein that folds in microseconds via a marginally cooperative downhill process (Li, P., Oliva, F. Y., Naganathan, A. N., and Muñoz, V. (2009) Proc. Natl. Acad. Sci. USA. 106, 103–108). Particularly, we explore the effects of salt-induced screening of the electrostatic interactions in BBL at neutral pH and in acid-denatured BBL. Our results show that electrostatic screening stabilizes the native state of the neutral and protonated forms, inducing complete refolding of acid-denatured BBL. Furthermore, without net electrostatic interactions, the unfolding process becomes much less cooperative, as judged by the broadness of the equilibrium unfolding curve and the relaxation rate. Our experiments show that the marginally cooperative unfolding of BBL can still be made twice as broad while the protein retains its ability to fold into the native three-dimensional structure in microseconds. This result demonstrates experimentally that efficient folding does not require cooperativity, confirming predictions from theory and computer simulations and challenging the conventional biochemical paradigm. Furthermore, we conclude that electrostatic interactions are an important factor in determining folding cooperativity. Thus, electrostatic modulation by pH-salt and/or mutagenesis of charged residues emerges as an attractive tool for tuning folding cooperativity.     

Folding and unfolding of gammaTIM monomers and dimers

Kinetic simulations of the folding and unfolding of triosephosphate isomerase (TIM) from yeast were conducted using a single monomer gammaTIM polypeptide chain that folds as a monomer and two gammaTIM chains that fold to the native dimer structure. The basic protein model used was a minimalist Gō model using the native structure to determine attractive energies in the protein chain. For each simulation type--monomer unfolding, monomer refolding, dimer unfolding, and dimer refolding--thirty simulations were conducted, successfully capturing each reaction in full. Analysis of the simulations demonstrates four main conclusions. First, all four simulation types have a similar "folding order", i.e., they have similar structures in intermediate stages of folding between the unfolded and folded state. Second, despite this similarity, different intermediate stages are more or less populated in the four different simulations, with 1), no intermediates populated in monomer unfolding; 2), two intermediates populated with beta(2)-beta(4) and beta(1)-beta(5) regions folded in monomer refolding; 3), two intermediates populated with beta(2)-beta(3) and beta(2)-beta(4) regions folded in dimer unfolding; and 4), two intermediates populated with beta(1)-beta(5) and beta(1)-beta(5) + beta(6) + beta(7) + beta(8) regions folded in dimer refolding. Third, simulations demonstrate that dimer binding and unbinding can occur early in the folding process before complete monomer-chain folding. Fourth, excellent agreement is found between the simulations and MPAX (misincorporation proton alkyl exchange) experiments. In total, this agreement demonstrates that the computational Gō model is accurate for gammaTIM and that the energy landscape of gammaTIM appears funneled to the native state.     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.     

Structural details,pathways, and energetics of unfolding apomyoglobin

Protein folding is often difficult to characterize experimentally because of the transience of intermediate states, and the complexity of the protein-solvent system. Atomistic simulations, which could provide more detailed information, have had to employ highly simplified models or high temperatures, to cope with the long time scales of unfolding; direct simulation of folding is even more problematic. We report a fully atomistic simulation of the acid-induced unfolding of apomyoglobin in which the protonation of acidic side-chains to simulate low pH is sufficient to induce unfolding at room temperature with no added biasing forces or other unusual conditions; and the trajectory is validated by comparison to experimental characterization of intermediate states. Novel insights provided by their analysis include: characterization of a dry swollen globule state forming a barrier to initial unfolding or final folding; observation of cooperativity in secondary and tertiary structure formation and its explanation in terms of dielectric environments; and structural details of the intermediate and the completely unfolded states. These insights involve time scales and levels of structural detail that are presently beyond the range of experiment, but come within reach through the simulation methods described here. An implicit solvation model is used to analyze the energetics of protein folding at various pH and ionic strength values, and a reasonable estimate of folding free energy is obtained. Electrostatic interactions are found to disfavor folding.     

About the structural role of disulfide bridges in serum albumins: evidence from protein simulated unfolding

The role of the 17 disulfide (S-S) bridges in preserving the native conformation of human serum albumin (HSA) is investigated by performing classical molecular dynamics (MD) simulations on protein structures with intact and, respectively, reduced S-S bridges. The thermal unfolding simulations predict a clear destabilization of the protein secondary structure upon reduction of the S-S bridges as well as a significant distortion of the tertiary structure that is revealed by the changes in the protein native contacts fraction. The effect of the S-S bridges reduction on the protein compactness was tested by calculating Gibbs free energy profiles with respect to the protein gyration radius. The theoretical results obtained using the OPLS-AA and the AMBER ff03 force fields are in agreement with the available experimental data. Beyond the validation of the simulation method, the results here reported provide new insights into the mechanism of the protein reductive/oxidative unfolding/folding processes. It is predicted that in the native conformation of the protein, the thiol (-SH) groups belonging to the same reduced S-S bridge are located in potential wells that maintain them in contact. The -SH pairs can be dispatched by specific conformational transitions of the peptide chain located in the neighborhood of the cysteine residues.     

Identifying the site of initial tertiary structure disruption during apomyoglobin unfolding.

Structural characterization of protein unfolding intermediates [Kiefhaber et al. (1995) Nature 375, 513; Hoeltzli et al.(1995) Proc. Natl. Acad. Sci. U.S.A. 92, 9318], which until recently were thought to be nonexistent, is beginning to give information on the mechanism of unfolding. To test for apomyoglobin unfolding intermediates, we monitored kinetics of urea-induced denaturation by stop-flow tryptophan fluorescence and quench-flow amide hydrogen exchange. Both measurements yield a single, measurable kinetic phase of identical rate, indicating that the reaction is highly cooperative. A burst phase in fluorescence, however, suggests that an intermediate is rapidly formed. To structurally characterize it, we carried out stop-flow thiol-disulfide exchange studies of 10 single cysteine-containing mutants. Cysteine probes buried at major sites of helix-helix pairing revealed that side chains throughout the protein unpack and become accessible to the labeling reagent [5, 5'-dithiobis (2-nitrobenzoic acid)] with one of two rates. Probes located at all helical-packing interfaces-except for one-become exposed at the rate of global unfolding as determined by fluorescence and hydrogen exchange measurements. In contrast, probes located at the A-E helical interface undergo complete thiol-disulfide exchange within the mixing dead time of 6 ms. These results point to the existence of a burst-phase unfolding intermediate that contains globally intact hydrogen bonds but locally disrupted side-chain packing interactions. Dissolution of secondary and tertiary structure are therefore not tightly coupled processes. We suggest that disruption of tertiary structure may be a stepwise process that begins at the weakest point of the native fold, as determined by native-state hydrogen-exchange parameters.     

More stable structure of wheat germ lipase at low pH than its native state

Wheat germ lipase is a cereal lipase which is a monomeric protein. In the present study we sought to structurally characterize this protein along with equilibrium unfolding in solution. Conformational changes occurring in the protein with varying pH, were monitored by circular dichroism (CD) spectroscopy, fluorescence emission spectroscopy, binding of hydrophobic dye, 1-anilino 8-naphthalenesulfonic acid (ANS) and dynamic light scattering (DLS). Our study showed that acid denaturation of lipase lead to characterization of multiple monomeric intermediates. Native protein at pH 7.0 showed far-UV spectrum indicating mixed structure with both alpha and beta-type of characteristics. Activity of lipase was found to fall on either sides of pH 7.0–8.0. Acid-unfolded state was characterized at pH 4.0 with residual secondary structure, disrupted tertiary spectrum and red-shifted fluorescence spectrum with decreased intensity. Further decrease in pH lead to formation of secondary structure and acid-induced molten globule state was found to be stabilized at pH 1.4, with exposed tryptophan residues and hydrophobic patches. Notably, interesting finding of this study was characterization of acid-induced state at pH 0.8 with higher secondary structure content than native lipase, regain in tertiary spectrum and induction of compact conformation. Although enzymatically inactive, acid-induced state at pH 0.8 was found to be structurally more stable than native lipase, as shown by chemical and thermal denaturation profiles.     

Geofold: topology-based protein unfolding pathways capture the effects of engineered disulfides on kinetic stability

Protein unfolding is modeled as an ensemble of pathways, where each step in each pathway is the addition of one topologically possible conformational degree of freedom. Starting with a known protein structure, GeoFold hierarchically partitions (cuts) the native structure into substructures using revolute joints and translations. The energy of each cut and its activation barrier are calculated using buried solvent accessible surface area, side chain entropy, hydrogen bonding, buried cavities, and backbone degrees of freedom. A directed acyclic graph is constructed from the cuts, representing a network of simultaneous equilibria. Finite difference simulations on this graph simulate native unfolding pathways. Experimentally observed changes in the unfolding rates for disulfide mutants of barnase, T4 lysozyme, dihydrofolate reductase, and factor for inversion stimulation were qualitatively reproduced in these simulations. Detailed unfolding pathways for each case explain the effects of changes in the chain topology on the folding energy landscape. GeoFold is a useful tool for the inference of the effects of disulfide engineering on the energy landscape of protein unfolding.     

Alternatively folded states of an immunoglobulin

Well-defined, non-native protein structures of low stability have been increasingly observed as intermediates in protein folding or as equilibrium structures populated under specific solvent conditions. These intermediate structures, frequently referred to as molten globule states, are characterized by the presence of secondary structure, a lack of significant tertiary contacts, increased hydrophobicity and partial specific volume as compared to native structures, and low cooperativity in thermal unfolding. The present study demonstrates that under acidic conditions (pH less than 3) the antibody MAK33 can assume a folded stable conformation. This A-state is characterized by a high degree of secondary structure, increased hydrophobicity, a native-like maximum wavelength of fluorescence emission, and a tendency toward slow aggregation. A prominent feature of this low-pH conformation is the stability against denaturant and thermal unfolding that is manifested in highly cooperative reversible phase transitions indicative of the existence of well-defined tertiary contacts. These thermodynamic results are corroborated by the kinetics of folding from the completely unfolded chain to the alternatively folded state at pH 2. The given data suggest that MAK33 at pH 2 adopts a cooperative structure that differs from the native immunoglobulin fold at pH 7. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its extraordinary structural stability that is characteristic for native protein structures.     

The acid-induced folded state of Sac7d is the native state

Sac7d unfolds at low pH in the absence of salt, with the greatest extent of unfolding obtained at pH 2. We have previously shown that the acid unfolded protein is induced to refold by decreasing the pH to 0 or by addition of salt (McCrary BS, Bedell J. Edmondson SP, Shriver JW, 1998, J Mol Biol 276:203-224). Both near-ultraviolet circular dichroism spectra and ANS fluorescence enhancements indicate that the acid- and salt-induced folded states have a native fold and are not molten globular. 1H,15N heteronuclear single quantum coherence NMR spectra confirm that the native, acid-, and salt-induced folded states are essentially identical. The most significant differences in amide 1H and 15N chemical shifts are attributed to hydrogen bonding to titrating carboxyl side chains and through-bond inductive effects. The 1H NMR chemical shifts of protons affected by ring currents in the hydrophobic core of the acid- and salt-induced folded states are identical to those observed in the native. The radius of gyration of the acid-induced folded state at pH 0 is shown to be identical to that of the native state at pH 7 by small angle X-ray scattering. We conclude that acid-induced collapse of Sac7d does not lead to a molten globule but proceeds directly to the native state. The folding of Sac7d as a function of pH and anion concentration is summarized with a phase diagram that is similar to those observed for other proteins that undergo acid-induced folding except that the A-state is encompassed by the native state. These results demonstrate that formation of a molten globule is not a general property of proteins that are refolded by acid.     

Stability of the glycerol facilitator in detergent solutions

Understanding membrane protein folding and stability is required for a molecular explanation of function and for the development of interventions in membrane protein folding diseases. Stable aqueous detergent solutions of the Escherichia coli glycerol facilitator in its native oligomeric state have been difficult to prepare as the protein readily unfolds and forms nonspecific aggregates. Here, we report a study of the structure and stability of the glycerol facilitator in several detergent solutions by Blue Native and sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD), and fluorescence. Protein tetramers were prepared in neutral dodecyl maltoside (DDM) and in zwitterionic lysomyristoylphosphatidylcholine (LMPC) detergent solutions that are stable during SDS-PAGE. Thermal unfolding experiments show that the protein is more stable in LMPC than in DDM. Tertiary structure unfolds before quaternary and some secondary structure in LMPC, whereas unfolding is more cooperative in DDM. The high stability of the protein in DDM is evident from the unfolding half-life of 8 days in 8 M urea, suggesting that hydrophobic interactions contribute to the stability. The protein unfolds readily in LMPC below pH 6, whereas the tetramer remains intact at pH 4 in DDM. At pH 4 in DDM, the protein is more sensitive than at neutral pH to unfolding by SDS and the effect is reversible. At pH 3 in DDM, the tetramer unfolds, losing its tertiary structure but retaining native helical structure which melts at significantly lower temperatures than in the native tetramer. The glycerol facilitator prepared in SDS is mainly monomeric and has ~10% less alpha-helix than the native protein. CD suggests that it forms a condensed structure with non-native tertiary contacts highly similar to the state observed in LMPC at low pH. The implications of the results for in vitro and in vivo folding of the protein are discussed.     

Molecular dynamics simulations of Hydrogenobacter thermophilus cytochrome c552: comparisons of the wild-type protein, a b-type variant, and the apo state

Molecular dynamic simulations have been performed for wild-type Hydrogenobacter thermophilus cytochrome c(552), a b-type variant of the protein, and the apo state with the heme prosthetic group removed. In the b-type variant, Cys 10 and Cys 13 were mutated to alanine residues, and so the heme group was no longer covalently bound to the protein. Two 8-ns simulations have been performed for each system at 298 and 360 K. The simulations of the wild-type protein at 298 K show a very close agreement with experimental NMR data. A fluxional process involving the side chain of Met 59, which coordinates to the heme iron, is observed in accord with proposals from NMR studies. Overall, the structure and dynamical behavior of the protein during the simulations of the b-type variant is closely similar to that of the wild-type protein. However, side chains in the heme-binding site show larger fluctuations in the b-type variant simulation at 360 K. In addition, structural changes are seen for a number of residues close to the heme group, particularly Gly 22 and Ser 51. The simulations of the apo state show significant conformational changes for residues 50-59. These residues form a loop region, which packs over the heme group in the wild-type protein and hydrogen bonds to the heme propionate groups. In the absence of heme, in the apo state simulations, these residues form short but persistent regions of beta-sheet secondary structure. These could provide nucleation sites for the conversion to amyloid fibrils.     

Mapping the Structure of an Integral Membrane Protein under Semi-Denaturing Conditions by Laser-Induced Oxidative Labeling and Mass Spectrometry

Little is known about the structural properties of semi-denatured membrane proteins. The current study employs laser-induced oxidative labeling of methionine side chains in combination with electrospray mass spectrometry and optical spectroscopy for gaining insights into the conformation of bacteriorhodopsin (BR) under partially denaturing conditions. The native protein shows extensive oxidation at M32, M68, and M163, which are located in solvent-accessible loops. In contrast, M20 (helix A), M56/60 (helix B), M118 (helix D), M145 (helix E), and M209 (helix G) are strongly protected, consistent with the known protein structure. Exposure of the protein to acidic conditions leads to a labeling pattern very similar to that of the native state. The absence of large-scale conformational changes at low pH is in agreement with recent crystallography data. Solubilization of BR in SDS induces loss of the retinal chromophore concomitant with collapse of the binding pocket, thereby precluding solvent access to the protein interior. Tryptophan fluorescence data confirm the presence of a large protein core that remains protected from water. However, oxidative labeling indicates partial unfolding of helices A and D in SDS. Irreversible thermal denaturation of the protein at 100 °C induces a labeling pattern quite similar to that seen upon SDS exposure. Labeling experiments on refolded bacterioopsin reveal a native-like structure, but with partial unfolding of helix D. Our data suggest that noncovalent contacts with the retinal chromophore in native BR play an important role for the stability of this particular helix. Overall, the present work illustrates the viability of using laser-induced oxidative labeling as a novel tool for characterizing structural changes of membrane proteins in response to alterations of their solvent environment.