Molecular basis for the effect of urea and guanidinium chloride on the dynamics of unfolded polypeptide chains

Chemical denaturants are frequently used to unfold proteins and to characterize mechanisms and transition states of protein folding reactions. The molecular basis of the effect of urea and guanidinium chloride (GdmCl) on polypeptide chains is still not well understood. Models for denaturant--protein interaction include both direct binding and indirect changes in solvent properties. Here we report studies on the effect of urea and GdmCl on the rate constants (k(c)) of end-to-end diffusion in unstructured poly(glycine-serine) chains of different length. Urea and GdmCl both lead to a linear decrease of lnk(c) with denaturant concentration, as observed for the rate constants for protein folding. This suggests that the effect of denaturants on chain dynamics significantly contributes to the denaturant-dependence of folding rate constants for small proteins. We show that this linear dependency is the result of two additive non-linear effects, namely increased solvent viscosity and denaturant binding. The contribution from denaturant binding can be quantitatively described by Schellman's weak binding model with binding constants (K) of 0.62(+/-0.01)M(-1) for GdmCl and 0.26(+/-0.01)M(-1) for urea. In our model peptides the number of binding sites and the effect of a bound denaturant molecule on chain dynamics is identical for urea and GdmCl. The results further identify the polypeptide backbone as the major denaturant binding site and give an upper limit of a few nanoseconds for residence times of denaturant molecules on the polypeptide chain.     

Kinetics of Contact Formation and End-to-End Distance Distributions of Swollen Disordered Peptides

Unstructured polypeptide chains are subject to various degrees of swelling or compaction depending on the combination of solvent condition and amino acid sequence. Highly denatured proteins generally behave like random-coils with excluded volume repulsion, whereas in aqueous buffer more compact conformations have been observed for the low-populated unfolded state of globular proteins as well as for naturally disordered sequences. To quantitatively account for the different mechanisms inducing the swelling of polypeptides, we have examined three 14-residues peptides in aqueous buffer and in denaturant solutions, including the well characterized AGQ repeat as a reference and two variants, in which we have successively introduced charged side chains and removed the glycines. Quenching of the triplet state of tryptophan by close contact with cysteine has been used in conjunction with Förster resonance energy transfer to study the equilibrium and kinetic properties of the peptide chains. The experiments enable accessing end-to-end root mean-square distance, probability of end-to-end contact formation and intrachain diffusion coefficient. The data can be coherently interpreted on the basis of a simple chain model with backbone angles obtained from a library of coil segments of proteins and hard sphere repulsion at each C_α position. In buffered water, we find that introducing charges in a glycine-rich sequence induces a mild chain swelling and a significant speed-up of the intrachain dynamics, whereas the removal of the glycines results in almost a two-fold increase of the chain volume and a drastic slowing down. In denaturants we observe a pronounced swelling of all the chains, with significant differences between the effect of urea and guanidinium chloride.     

Individual cartilage aggrecan macromolecules and their constituent glycosaminoglycans visualized via atomic force microscopy

Atomic force microscopy was used in ambient conditions to directly image dense and sparse monolayers of bovine fetal epiphyseal and mature nasal cartilage aggrecan macromolecules adsorbed on mica substrates. Distinct resolution of the non-glycosylated N-terminal region from the glycosaminoglycan (GAG) brush of individual aggrecan monomers was achieved, as well as nanometer-scale resolution of individual GAG chain conformation and spacing. Fetal aggrecan core protein trace length (398+/-57 nm) and end-to-end length (257+/-87 nm) were both larger than that of mature aggrecan (352+/-88 and 226+/-81 nm, respectively). Similarly, fetal aggrecan GAG chain trace length (41+/-7 nm) and end-to-end (32+/-8 nm) length were both larger than that of mature aggrecan GAG (32+/-5 and 26+/-7 nm, respectively). GAG-GAG spacing along the core protein was significantly smaller in fetal compared to mature aggrecan (3.2+/-0.8 and 4.4+/-1.2nm, respectively). Together, these differences between the two aggrecan types were likely responsible for the greater persistence length of the fetal aggrecan (110 nm) compared to mature aggrecan (82 nm) calculated using the worm-like chain model. Measured dimensions and polymer statistical analyses were used in conjunction with the results of Western analyses, chromatographic, and carbohydrate electrophoresis measurements to better understand the dependence of aggrecan structure and properties on its constituent GAG chains.     

Estimation of the diffusion-limited rate of microtubule assembly.

Microtubule assembly is a complex process with individual microtubules alternating stochastically between extended periods of assembly and disassembly, a phenomenon known as dynamic instability. Since the discovery of dynamic instability, molecular models of assembly have generally assumed that tubulin incorporation into the microtubule lattice is primarily reaction-limited. Recently this assumption has been challenged and the importance of diffusion in microtubule assembly dynamics asserted on the basis of scaling arguments, with tubulin gradients predicted to extend over length scales exceeding a cell diameter, approximately 50 microns. To assess whether individual microtubules in vivo assemble at diffusion-limited rates and to predict the theoretical upper limit on the assembly rate, a steady-state mean-field model for the concentration of tubulin about a growing microtubule tip was developed. Using published parameter values for microtubule assembly in vivo (growth rate = 7 microns/min, diffusivity = 6 x 10(-12) m2/s, tubulin concentration = 10 microM), the model predicted that the tubulin concentration at the microtubule tip was approximately 89% of the concentration far from the tip, indicating that microtubule self-assembly is not diffusion-limited. Furthermore, the gradients extended less than approximately 50 nm (the equivalent of about two microtubule diameters) from the microtubule tip, a distance much less than a cell diameter. In addition, a general relation was developed to predict the diffusion-limited assembly rate from the diffusivity and bulk tubulin concentration. Using this relation, it was estimated that the maximum theoretical assembly rate is approximately 65 microns/min, above which tubulin can no longer diffuse rapidly enough to support faster growth.     

Characterization of beta-connectin (titin 2) from striated muscle by dynamic light scattering.

Connectin (titin) is a large filamentous protein (single peptide) with a molecular mass of approximately 3 MDa, contour length approximately 900 nm, and diameter approximately 4 nm, and resides in striated muscle. Connectin links the thick filaments to the Z-lines in a sarcomere and produces a passive elastic force when muscle fiber is stretched. The aim of this study is to elucidate some aspects of physical properties of isolated beta-connectin (titin 2), a proteolytic fragment of connectin, by means of dynamic light-scattering (DLS) spectroscopy. The analysis of DLS spectra for beta-connectin gave the translational diffusion coefficient of 3.60 x 10(-8) cm2/s at 10 degrees C (or the hydrodynamic radius of 44.1 nm), molecular mass little smaller than 3.0 MDa (for a literature value of sedimentation coefficient), the root-mean-square end-to-end distance of 163 nm (or the radius of gyration of 66.6 nm), and the Kuhn segment number of 30 and segment length of 30 nm (or the persistence length of 15 nm). These results permitted to estimate the flexural rigidity of 6.0 x 10(-20) dyn x cm2 for filament bending, and the elastic constant of 7 dyn/cm for extension of one persistence length. Based on a simple model, implications of the present results in muscle physiology are discussed.     

Dynamic force spectroscopy of E. coli P pili

Surface organelles (so-called pili) expressed on the bacterial membrane mediate the adhesion of Escherichia coli causing urinary tract infection. These pili possess some extraordinary elongation properties that are assumed to allow a close bacterium-to-host contact even in the presence of shear forces caused by urine flow. The elongation properties of P pili have therefore been assessed for low elongation speeds (steady-state conditions). This work reports on the behavior of P pili probed by dynamic force spectroscopy. A kinetic model for the unfolding of a helixlike chain structure is derived and verified. It is shown that the unfolding of the quaternary structure of the PapA rod takes place at a constant force that is almost independent of elongation speed for slow elongations (up to approximately 0.4 mum/s), whereas it shows a dynamic response with a logarithmic dependence for fast elongations. The results provide information about the energy landscape and reaction rates. The bond length and thermal bond opening and closure rates for the layer-to-layer bond have been assessed to approximately 0.76 nm, approximately 0.8 Hz, and approximately 8 GHz, respectively. The results also support a previously constructed sticky-chain model for elongation of the PapA rod that until now had been experimentally verified only under steady-state conditions.     

Electrochemical measurement of lateral diffusion coefficients of ubiquinones and plastoquinones of various isoprenoid chain lengths incorporated in model bilayers.

The long-range diffusion coefficients of isoprenoid quinones in a model of lipid bilayer were determined by a method avoiding fluorescent probe labeling of the molecules. The quinone electron carriers were incorporated in supported dimyristoylphosphatidylcholine layers at physiological molar fractions (<3 mol%). The elaborate bilayer template contained a built-in gold electrode at which the redox molecules solubilized in the bilayer were reduced or oxidized. The lateral diffusion coefficient of a natural quinone like UQ10 or PQ9 was 2.0 +/- 0.4 x 10(-8) cm2 s(-1) at 30 degrees C, two to three times smaller than the diffusion coefficient of a lipid analog in the same artificial bilayer. The lateral mobilities of the oxidized or reduced forms could be determined separately and were found to be identical in the 4-13 pH range. For a series of isoprenoid quinones, UQ2 or PQ2 to UQ10, the diffusion coefficient exhibited a marked dependence on the length of the isoprenoid chain. The data fit very well the quantitative behavior predicted by a continuum fluid model in which the isoprenoid chains are taken as rigid particles moving in the less viscous part of the bilayer and rubbing against the more viscous layers of lipid heads. The present study supports the concept of a homogeneous pool of quinone located in the less viscous region of the bilayer.     

Rate of intrachain contact formation in an unfolded protein: temperature and denaturant effects

We have measured the effect of temperature and denaturant concentration on the rate of intrachain diffusion in an unfolded protein. After photodissociating a ligand from the heme iron of unfolded horse cytochrome c, we use transient optical absorption spectroscopy to measure the time scale of the diffusive motions that bring the heme, located at His18, into contact with its native ligand, Met80. Measuring the rate at which this 62 residue intrachain loop forms under both folding and unfolding conditions, we find a significant effect of denaturant on the chain dynamics. The diffusion of the chain accelerates as denaturant concentration decreases, with the contact formation rate approaching a value near approximately 6x10(5) s(-1) in the absence of denaturant. This result agrees well with an extrapolation from recent loop formation measurements in short synthetic peptides. The temperature dependence of the rate of contact formation indicates an Arrhenius activation barrier, Ea approximately 20 kJ/mol, at high denaturant concentrations, comparable to what is expected from solvent viscosity effects alone. Although Ea increases by several kBT as denaturant concentration decreases, the overall rate of diffusion nevertheless increases. These results indicate that inter-residue energetic interactions do not control conformational diffusion in unfolded states, even under folding conditions.     

Mathematical and experimental analyses of antibody transport in hollow-fiber-based specific antibody filters

We are developing hollow fiber-based specific antibody filters (SAFs) that selectively remove antibodies of a given specificity directly from whole blood, without separation of the plasma and cellular blood components and with minimal removal of plasma proteins other than the targeted pathogenic antibodies. A principal goal of our research is to identify the primary mechanisms that control antibody transport within the SAF and to use this information to guide the choice of design and operational parameters that maximize the SAF-based antibody removal rate. In this study, we formulated a simple mathematical model of SAF-based antibody removal and performed in vitro antibody removal experiments to test key predictions of the model. Our model revealed three antibody transport regimes, defined by the magnitude of the Damk?hler number Da (characteristic antibody-binding rate/characteristic antibody diffusion rate): reaction-limited (Da /= 10). For a given SAF geometry, blood flow rate, and antibody diffusivity, the highest antibody removal rate was predicted for diffusion-limited antibody transport. Additionally, for diffusion-limited antibody transport the predicted antibody removal rate was independent of the antibody-binding rate and hence was the same for any antibody-antigen system and for any patient within one antibody-antigen system. Using SAF prototypes containing immobilized bovine serum albumin (BSA), we measured anti-BSA removal rates consistent with transport in the intermediate regime (Da approximately 3). We concluded that initial SAF development work should focus on achieving diffusion-limited antibody transport by maximizing the SAF antibody-binding capacity (hence maximizing the characteristic antibody-binding rate). If diffusion-limited antibody transport is achieved, the antibody removal rate may be raised further by increasing the number and length of the SAF fibers and by increasing the blood flow rate through the SAF.     

Segmental conformational disorder and dynamics in the intrinsically disordered protein α-synuclein and its chain length dependence

Conformational ensembles of fully disordered natural polypeptides represent the starting point of protein refolding initiated by transfer to folding conditions. Thus, understanding the transient properties and dimensions of such peptides under folding conditions is a necessary step in the understanding of their subsequent folding behavior. Such ensembles can also undergo alternative folding and form amyloid structures, which are involved in many neurological degenerative diseases. Here, we performed a structural study of this initial state using time-resolved fluorescence resonance energy transfer analysis of a series of eight partially overlapping double-labeled chain segments of the N-terminal and NAC domains of the α-synuclein molecule. The distributions of end-to-end distance and segmental intramolecular diffusion coefficients were simultaneously determined for eight labeled chain segments. We used the coefficient of variation, C_v, as a measure of the conformational heterogeneity (i.e., structural disorder). With the exception of two segments, the C_vs were characteristic of a fully disordered state of the chain. Subtle deviations from this behavior at the segment labeled in the NAC domain and the segment at the N termini reflected subtle conformational bias that might be related to the initiation of transition to amyloid aggregates. The chain length dependence of the mean segmental end-to-end distance followed a power law as predicted by Flory, but the dependence was steeper than previously predicted, probably due to the contribution of the excluded volume effect, which is more dominant for shorter-chain segments. The observed intramolecular diffusion coefficients (< 10 to ∼ 25 ?²/ns) are only an order of magnitude lower than the common diffusion coefficients of low molecular weight probes. This diffusion coefficient increased with chain length, probably due to the cumulative contributions of minor bond rotations along the chain. These results gave us a reference both for characteristics of a natural unfolded polypeptide at the moment of initiation of folding and for detection of possible initiation sites of the amyloid transition.     

Effects of urea and guanidine hydrochloride on peptide and nonpolar groups

The free energy transfer of several N-acetyl(glycine)n ethyl esters (n = 1-3) and side chain derivatives (Ala, Val, Nva, Leu, Nle, and Phe) from water to urea and guanidine hydrochloride solutions has been determined from the solubility and distribution coefficients of these compounds between aqueous and nonaqueous phases. These uncharged model peptides, unlike the amino acids used for a similar study, avoid complication due to charge effects for the transfer process. The compounds with an increase in the number of glycyl groups show additivity of the group free energy toward the transfer from water to urea solution but not to guanidine hydrochloride solution. The derivatives with a side chain show that the principle of group additivity does not hold true for the aliphatic side chains for the transfer to either urea or guanidine hydrochloride solutions. In fact, the free energy of transfer of the side chains, viz., aliphatic ones, is found to be energetically unfavorable in moderately high denaturant concentration. Phenylalanyl, the only aromatic side chain studied here, showed a favorable free energy of transfer to the denaturant solutions. In addition, the values of the favorable free energy obtained in this study are much smaller than the values obtained from the study of the amino acids. The transfer of the glycyl group to the denaturant solutions is exothermic whereas the transfer of the side chains is endothermic in nature.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hydrogen-Bond Driven Loop-Closure Kinetics in Unfolded Polypeptide Chains

Characterization of the length dependence of end-to-end loop-closure kinetics in unfolded polypeptide chains provides an understanding of early steps in protein folding. Here, loop-closure in poly-glycine-serine peptides is investigated by combining single-molecule fluorescence spectroscopy with molecular dynamics simulation. For chains containing more than 10 peptide bonds loop-closing rate constants on the 20–100 nanosecond time range exhibit a power-law length dependence. However, this scaling breaks down for shorter peptides, which exhibit slower kinetics arising from a perturbation induced by the dye reporter system used in the experimental setup. The loop-closure kinetics in the longer peptides is found to be determined by the formation of intra-peptide hydrogen bonds and transient β-sheet structure, that accelerate the search for contacts among residues distant in sequence relative to the case of a polypeptide chain in which hydrogen bonds cannot form. Hydrogen-bond-driven polypeptide-chain collapse in unfolded peptides under physiological conditions found here is not only consistent with hierarchical models of protein folding, that highlights the importance of secondary structure formation early in the folding process, but is also shown to speed up the search for productive folding events.     

Kinetics of aggregation of synthetic beta-amyloid peptide.

beta-Amyloid peptide is the major protein component of senile plaques and cerebrovascular amyloid deposits in patients with Alzheimer's disease. The peptide deposits extracellularly in the form of amyloid fibrils, in a cross-beta conformation. beta-amyloid peptide is a 39- to 43-residue segment of a normal membrane precursor protein. In this work, a peptide homologous to the first 40 amino acids of beta-amyloid peptide, beta(1-40), was synthesized and characterized. beta(1-40) exhibited a sharp change in solubility near physiological pH and gel formation at concentrations of 3 mg/ml or greater. Circular dichroism indicated that beta(1-40) contained approximately two-thirds beta-structure, but no alpha-helical character. Quasi-elastic and classical light scattering measurements showed that beta(1-40) aggregated end-to-end in solution, reaching average molecular weights greater than 4 x 10(6) after 13 days. The aggregates were best modeled as rigid rods of 5 nm diameter, similar to the diameter of amyloid fibrils purified from plaques. A mathematical model based on diffusion-limited aggregation was developed to describe the kinetics of aggregation.     

Evidence for a Partially Structured State of the Amylin Monomer

Islet amyloid polypeptide (amylin) is the main component in amyloid deposits formed in type II diabetes. We used triplet quenching to probe the dynamics of contact formation between the N-terminal disulfide loop and a C-terminal tryptophan in monomeric amylins from human and rat. Quenching rates measured in the absence of denaturant are four times larger than those in 6 M guanidinium chloride, indicating a decrease in the average end-to-end distance (collapse) at low denaturant concentrations. We were surprised to find an even greater (sevenfold) increase in quenching rates on removal of denaturant for a hydrophilic control peptide containing the disulfide loop compared to the same peptide without the loop (twofold change). These results suggest that collapse is driven by backbone-backbone and backbone-side chain interactions involving the disulfide loop portion of the chain rather than by the formation of side-chain hydrophobic contacts. Molecular dynamics simulations of the control peptide show that the collapse results from hydrogen-bonding interactions between the central residues of the chain and the disulfide loop. The quenching experiments also indicate that the monomer of the human, amyloidogenic form of amylin is more compact than the rat form, which does not form amyloid. We discuss these newly observed differences between human and rat amylin in solution and their possible relation to aggregation and to the physiological function of amylin binding to the calcitonin receptor.     

Computer simulation of hydrodynamic properties of semiflexible macromolecules: randomly broken chains, wormlike chains, and analysis of properties of DNA

The translational and rotational diffusion coefficients and the intrinsic viscosity of semiflexible, randomly broken, and wormlike chains have been obtained by Monte Carlo simulation in the context of the rigid-body treatment. Both approximate and rigorous rigid-body hydrodynamics are used, so that the error introduced by the approximate methods can be evaluated. A randomly broken chain and a wormlike chain having the same contour length and persistence length have the same radius of gyration but different values for any of the hydrodynamic properties. The two types of chains are compared in this regard. Considering that the cross section of the chain is represented by a cylinder better than by a string of spheres, we devise a cylindrical correction to be applied to the results simulated for chains of beads. Application is made to the analysis of experimental data for the translational and rotational coefficients of DNA fragments with up to 10(3) base pairs, obtaining the persistence length for each model. The values for the wormlike chain agree well with model-independent values obtained from radii of gyration and with other literature data at varying ionic strength. The randomly broken chain is equally able to reproduce the experimental length dependence of the properties, but the resulting persistence length may be too high.     

Measuring dynamic flexibility of the coil state of a helix-forming peptide

To investigate the dynamic flexibility of the coil state of a helix-forming peptide the end-to-end contact rate was determined. Nanosecond optical excitation of tryptophan at one end of a 22 residue, alanine peptide populates a long-lived triplet state which is quenched upon close contact with a cyclic disulfide attached to the opposite end. Analysis of the decay of the triplet population using a two-state model for helix formation yields the diffusion-limited end-to-end contact rate of the coil state of the peptide as well as the helix-->coil and coil-->helix rates. The helix-coil rates are very similar to those previously measured in laser temperature-jump experiments. The end-to-end contact rate of 1.1 x 10(7) s(-1) in the coil state is tenfold faster than the rate for a disordered peptide with threonine substituted for alanine and, somewhat surprisingly, is about twice the rate for a disordered glycine-containing peptide. These differences are discussed in terms of the theory of Szabo, Schulten and Schulten. The rates should provide important new benchmarks for testing the accuracy of atomistic molecular dynamics simulations.     

Protein stabilization by urea and guanidine hydrochloride

The urea, guanidine hydrochloride, salt, and temperature dependence of the rate of dissociation of CO from a nonequilibrium state of CO-bound native ferrocytochrome c has been studied at pH 7. The heme iron of ferrocytochrome c in the presence of denaturing concentrations of guanidine hydrochloride (GdnHCl) and urea prepared in 0.1 M phosphate, pH 7, binds CO. When the unfolded protein solution is diluted 101-fold into CO-free folding buffer, the protein chain refolds completely, leaving the CO molecule bonded to the heme iron. Subsequently, slow thermal dissociation of the CO molecule yields to the heme coordination of the native M80 ligand. Thus, the reaction monitors the rate of thermal conversion of the CO-liganded native ferrocytochrome c to the M80-liganded native protein. The rate of this reaction, k(diss), shows a characteristic dependence on the presence of nondenaturing concentrations of the denaturants in the reaction medium. The rate decreases by approximately 1.9-3-fold as the concentration of GdnHCl in the refolding medium increases from nearly 0 to approximately 2.1 M. Similarly, the rate decreases by 1.8-fold as the urea concentration is raised from 0.l to approximately 5 M. At still higher concentrations of the denaturants the denaturing effect sets in, the protein is destabilized, and hence the CO dissociation rate increases sharply. The activation energy of the reaction, E(a), increases when the denaturant concentration in the reaction medium is raised: from 24.1 to 28.3 kcal mol(-1) for a 0.05-2.1 M rise in GdnHCl and from 25.2 to 26.9 kcal mol(-1) for a 0.1-26.9 M increase in urea. Corresponding to these increases in denaturant concentrations are also increases in the activation entropy, S(diss)/R, where R is the gas constant of the reaction. The denaturant dependence of these kinetic and thermodynamic parameters of the CO dissociation reaction suggests that binding interactions with GdnHCl and urea can increase the structural and energetic stability of ferrocytochrome c up to the limit of the subdenaturing concentrations of the additives. NaCl and Na(2)SO(4), which stabilize proteins through their salting-in effect, also decrease the rate with a corresponding increase in activation entropy of CO dissociation from CO-bound native ferrocytochrome c, lending support to the view that low concentrations of GdnHCl and urea stabilize proteins. These results have direct relevance to the understanding and interpretation of the free energy-denaturant relationship and protein folding chevrons.     

Dynamics of unfolded polypeptide chains in crowded environment studied by fluorescence correlation spectroscopy

Proteins have evolved to fold and function within a cellular environment that is characterized by high macromolecular content. The earliest step of protein folding represents intrachain contact formation of amino acid residues within an unfolded polypeptide chain. It has been proposed that macromolecular crowding can have significant effects on rates and equilibria of biomolecular processes. However, the kinetic consequences on intrachain diffusion of polypeptides have not been tested experimentally, yet. Here, we demonstrate that selective fluorescence quenching of the oxazine fluorophore MR121 by the amino acid tryptophan (Trp) in combination with fast fluorescence correlation spectroscopy (FCS) can be used to monitor end-to-end contact formation rates of unfolded polypeptide chains. MR121 and Trp were incorporated at the terminal ends of polypeptides consisting of repetitive units of glycine (G) and serine (S) residues. End-to-end contact formation and dissociation result in "off" and "on" switching of MR121 fluorescence and underlying kinetics can be revealed in FCS experiments with nanosecond time resolution. We revisit previous experimental studies concerning the dependence of end-to-end contact formation rates on polypeptide chain length, showing that kinetics can be described by Gaussian chain theory. We further investigate effects of solvent viscosity and temperature on contact formation rates demonstrating that intrachain diffusion represents a purely diffusive, entropy-controlled process. Finally, we study the influence of macromolecular crowding on polypeptide chain dynamics. The data presented demonstrate that intrachain diffusion is fast in spite of hindered diffusion caused by repulsive interactions with macromolecules. Findings can be explained by effects of excluded volume reducing chain entropy and therefore accelerating the loop search process. Our results suggest that within a cellular environment the early formation of structural elements in unfolded proteins can still proceed quite efficiently in spite of hindered diffusion caused by high macromolecular content.     

Thermodynamic analysis of chain-melting transition temperatures for monounsaturated phospholipid membranes: dependence on cis-monoenoic double bond position.

Unsaturated phospholipid is the membrane component that is essential to the dynamic environment needed for biomembrane function. The dependence of the chain-melting transition temperature, T(t), of phospholipid bilayer membranes on the position, n(u), of the cis double bond in the glycerophospholipid sn-2 chain can be described by an expression of the form T(t) = T(t)(c)(1 + h'(c)|n(u) - n(c)|)/(1 + s'(c)|n(u) - n(c)|), where n(c) is the chain position of the double bond corresponding to the minimum transition temperature, T(t)(c), for constant diacyl lipid chain lengths. This implies that the incremental transition enthalpy (and entropy) contributed by the sn-2 chain is greater for whichever of the chain segments, above or below the double-bond position, is the longer. The critical position, n(c), of the double bond is offset from the center of the sn-2 chain by an approximately constant amount, deltan(c) approximately 1. 5 C-atom units. The dependence of the parameters T(t)(c), h'(c), and s'(c) on sn-1 and sn-2 chain lengths can be interpreted consistently when allowance is made for the chain packing mismatch between the sn-1 and sn-2 chains. The length of the sn-2 chain is reduced by approximately 0.8 C-atom units by the cis double bond, in addition to a shortening by approximately 1.3 C-atom units by the bent configuration at the C-2 position. Based on this analysis, a general thermodynamic expression is proposed for the dependence of the chain-melting transition temperature on the position of the cis double bond and on the sn-1 and sn-2 chain lengths. The above treatment is restricted mostly to double-bond positions close to the center of the sn-2 chain. For double bonds positioned closer to the carboxyl or terminal methyl ends of the sn-2 chain, the effects on transition enthalpy can be considerably larger. They may be interpreted by the same formalism, but with different characteristic parameters, h'(c) and s'(c), such that the shorter of the chain segments makes a considerably smaller contribution to the calorimetric properties of the chain-melting transition.