A molecular model for lipid-protein interaction in membranes: the role of hydrophobic mismatch.

The interaction free energy between a hydrophobic, transmembrane, protein and the surrounding lipid environment is calculated based on a microscopic model for lipid organization. The protein is treated as a rigid hydrophobic solute of thickness dP, embedded in a lipid bilayer of unperturbed thickness doL. The lipid chains in the immediate vicinity of the protein are assumed to adjust their length to that of the protein (e.g., they are stretched when dP > doL) in order to bridge over the lipid-protein hydrophobic mismatch (dP-doL). The bilayer's hydrophobic thickness is assumed to decay exponentially to its asymptotic, unperturbed, value. The lipid deformation free energy is represented as a sum of chain (hydrophobic core) and interfacial (head-group region) contributions. The chain contribution is calculated using a detailed molecular theory of chain packing statistics, which allows the calculation of conformational properties and thermodynamic functions (in a mean-field approximation) of the lipid tails. The tails are treated as single chain amphiphiles, modeled using the rotational isometric state scheme. The interfacial free energy is represented by a phenomenological expression, accounting for the opposing effects of head-group repulsions and hydrocarbon-water surface tension. The lipid deformation free energy delta F is calculated as a function of dP-doL. Most calculations are for C14 amphiphiles which, in the absence of a protein, pack at an average area per head-group ao approximately equal to 32 A2 (doL approximately 24.5 A), corresponding to the fluid state of the membrane. When dP = doL, delta F > 0 and is due entirely to the loss of conformational entropy experienced by the chains around the protein. When dP > doL, the interaction free energy is further increased due to the enhanced stretching of the tails. When dP < doL, chain flexibility (entropy) increases, but this contribution to delta F is overcounted by the increase in the interfacial free energy. Thus, delta F obtains a minimum at dP-doL approximately 0. These qualitative interpretations are supported by detailed numerical calculations of the various contributions to the interaction free energy, and of chain conformational properties. The range of the perturbation of lipid order extends typically over few molecular diameters. A rather detailed comparison of our approach to other models is provided in the discussion.     

Thermodynamics of the temperature-induced unfolding of globular proteins.

The heat capacity, enthalpy, entropy, and Gibbs energy changes for the temperature-induced unfolding of 11 globular proteins of known three-dimensional structure have been obtained by microcalorimetric measurements. Their experimental values are compared to those we calculate from the change in solvent-accessible surface area between the native proteins and the extended polypeptide chain. We use proportionality coefficients for the transfer (hydration) of aliphatic, aromatic, and polar groups from gas phase to aqueous solution, we estimate vibrational effects, and we discuss the temperature dependence of each constituent of the thermodynamic functions. At 25 degrees C, stabilization of the native state of a globular protein is largely due to two favorable terms: the entropy of non-polar group hydration and the enthalpy of interactions within the protein. They compensate the unfavorable entropy change associated with these interactions (conformational entropy) and with vibrational effects. Due to the large heat capacity of nonpolar group hydration, its stabilizing contribution decreases quickly at higher temperatures, and the two unfavorable entropy terms take over, leading to temperature-induced unfolding.     

Theoretical study of aliphatic chain structure in mono- and bilayers

A model is proposed to explain the aliphatic chain structure in mono- and bilayers (monomolecular films, aqueous lipid dispersions and biologival membranes). The rotational isomerism around carbon-carbon bonds and chain interactions (dispersion forces) are taken into account. This model has been first used to calculate different conformational properties of normal hexan layers in terms of distance between chains (thickness, birefringence, melting entropy, order parameter S₃ introduced in spin label technique).     

Reorientational and Conformational Ordering Processes at Elevated Pressures in 1,2-Dioleoyl Phosphatidylcholine: A Raman and Infrared Spectroscopic Study

Raman and infrared spectra of fully hydrated bilayers of 1,2-dioleoyl phosphatidylcholine (DOPC) were measured at increasing hydrostatic pressures up to -37 kbar. Under ambient conditions aqueous dispersions of DOPC are in the liquid crystalline state. The application of an external hydrostatic pressure induces conformational and dynamic ordering processes in DOPC, which trigger a first-order structural phase transition at 5 kbar from a disordered liquid crystalline state to a highly ordered gel state. In the gel phase the methylene chains of each molecule are fully extended and the two all-trans chain segments on both sides of the rigid cis double bond form a bent structure. The bent oleoyl chains in each molecule, as well as in neighboring molecules are packed parallel to each other. To achieve this parallel interchain packing, the double bonds of the sn-1 and sn-2 chains of each molecule must be aligned at the same position with respect to the bilayer interface which is achieved by a rotation of the C—C bonds in the glycerol moiety in the head group. The extremely strong interchain interactions in the gel phase of DOPC are unique for this lipid with cis dimono-unsaturated acyl chains. Our experimental results suggest that in the pressure-induced gel phase of DOPC the olefinic CH bonds are rotated out of the phase of the bent oleoyl chains and that the oleoyl chains of opposing bilayers bend towards opposite directions.     

Structural changes and fluctuations of proteins: I. A statistical thermodynamic model

A general theory of the structural changes and fluctuations of proteins has been proposed based on statistical thermodynanic considerations at the chain level.The “structure” of protein was assumed to be characterized by the state of secondary bonds between unique pairs of specific sites on peptide chains. Every secondary bond changes between the bonded and unboned states by thermal agitation and the “structure” is continuously fluctuating. The free energy of the “structural state” that is defined by the fraction of secondary bonds in the bonded state has been expressed by the bond energy, the cooperative interaction between bonds, the mixing entropy of bonds, and the entropy of polypeptide chains. The most probable “structural state” can be simply determined by graphical analysis and the effect of temperature or solvent composition on it is discussed. The temperature dependence of the free energy, the probability distribution of structural states and the specific heat have been calculated for two examples of structural change.The theory predicts two different types of structural changes from the ordered to disordered state, a “structural transition” and a “gradual structural change” with rising temperature, In the “structural transition”, the probability distribution has two maxima in the temperature range of transition. In the “gradual structural change”, the probability distribution has only one maximum during the change.A considerable fraction of secondary bonds is in the unbonded state and is always fluctuating even in the ordered state at room temperature. Such structural fluctuations in a single protein molecule have been discussed quantitatively.The theory is extended to include small molecules which bind to the protein molecule and affect the structural state. The changes of structural state caused by specific and non-specific binding and allosteric effects are explained in a unified manner.     

Electrostatically induced change of the conformational order of charged lipid membranes

Raman spectroscopy and X-ray diffraction are used to investigate the influence of surface charges on the structure of ionizable lipid membranes of dimyristoylmethylphosphatidic acid. The membrane surface charge density is regulated by varying the pH of the aqueous phase. Changes of the conformational order of the lipid chains are determined from the intensity of the CC stretch chain vibrations around 1100 cm^?1 in a lipid Raman spectrum. In going from an electrical neutral to a negatively charged membrane, the conformational order is reduced by 5% in the ordered and by 9% in the fluid membrane phase, corresponding to 0.6 and 0.8 CC bonds, respectively, which change from a trans to a gauche conformation. The electrostatically induced conformational change is mainly concentrated at the lipid chain ends as indicated by the spectral variations of the 890 cm^?1 CH₃ rocking band of the chain termini. The X-ray diffraction experiments show that increasing the surface charge density in the ordered membrane phase leads to a lateral expansion of the packing of the lipid polar groups, whereas the packing of the lipid chains in a plane perpendicular to the chain axes remains constant, indicating an increase of the tilt of the lipid chains from δ = 10° (pH 3) to δ = 27° (pH 9).     

Raman spectroscopic study of polycrystalline mono- and polyunsaturated 1-eicosanoyl-d(39)-2-eicosenoyl-sn-glycero-3-phosphocholines: bilayer lipid clustering and acyl chain order and disorder characteristics

Polycrystalline lipid samples of a series of mono- and polyunsaturated, double bond positional isomers of 1-eicosanoyl-d(39)-2-eicosenoyl-sn-glycero-3-phosphocholines [C(20-d(39)):C(20:1 Delta(j))PC, with j = 5, 8, 11, or 13; C(20-d(39)):C(20:2 Delta(11,14))PC; and C(20-d(39)):C(20:3 Delta(11, 14,17))PC] were investigated using vibrational Raman spectroscopy to assess the acyl chain packing order-disorder characteristics and putative bilayer cluster formation of the isotopically differentiated acyl chains. Perdeuteration of specifically the saturated sn-1 acyl chains for these bilayer systems enables each chain's intra- and intermolecular conformational and organizational properties to be evaluated separately. Various saturated chain methylene CD(2) and carbon-carbon (C&bond;C) stretching mode peak height intensity ratios and line width parameters for the polycrystalline samples demonstrate a high degree of sn-1 chain order that is unaffected by either the double bond placement or number of unsaturated bonds within the sn-2 chain. In contrast, the unsaturated sn-2 chain spectral signatures reflect increasing acyl chain conformational disorder as either the cis double bond is generally repositioned toward the chain terminus or the number of double bonds increases from one to three. The lipid bilayer chain packing differences observed between the sn-1 and sn-2 chains of this series of monounsaturated and polyunsaturated 20 carbon chain lipids suggest the existence of laterally distributed microdomains predicated on the formation of highly ordered, saturated sn-1 chain clusters.     

Structural determinants of miscibility in surface films of galactosylceramide and phosphatidylcholine: effect of unsaturation in the galactosylceramide acyl chain

The Langmuir film balance technique has been used to define the surface structure and determine the mixing behavior of galactosylceramide (GalCer) and phosphatidylcholines in surface phases. To determine the effect of unsaturation on surface behavior, chain-pure GalCer species containing either oleoyl (18:1 delta 9), eicosenoyl (20:1 delta 11), or eicosadienoyl (20:2 delta 11,14) fatty acyl chains were synthesized. Using bovine brain GalCer as a reference, surface pressure versus molecular area (phi-A) isotherms of the pure lipids were measured and analyzed by determining their compressibilities and by using an equation of state for lipid monolayers. This information, when coupled with surface potential versus molecular area (delta V-A) analyses, provides insights into GalCer surface structure in terms of molecular packing and orientation. Lipid mixing behavior was determined by classical approaches which involve analyzing the average molecular area, the average surface dipole moment, and surface pressure as a function of film composition. The results indicate that, in contrast to the complex mixing behavior displayed by bovine brain GalCer and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), chain-pure GalCer species containing either oleoyl, eicosenoyl, or eicosadienoyl fatty acyl chains are miscible with POPC over the entire composition range. Moreover, increasing amounts of GalCer containing eicosenoyl acyl chains systematically elevate dipalmitoylphosphatidylcholine's (DPPC) liquid-expanded-to-liquid-condensed transition pressure. Such behavior is consistent with GalCer being miscible with the liquid-expanded phase of DPPC. Thus, fatty acyl unsaturation is a critical parameter governing the mixing behavior of GalCer with phosphatidylcholine.     

How membrane chain-melting phase-transition temperature is affected by the lipid chain asymmetry and degree of unsaturation: an effective chain-length model.

Hydrocarbon effects on the lipid chain-melting phase-transition temperature are analyzed. The membrane fluidization temperature is shown to increase with the effective chain length, which is proportional to the thickness of the well-packed hydrocarbon region. The latter, as a rule, increases with the length of the longest ordered and aligned segment on each chain. This conclusion is independent of the cause for the reduced chain packing in membrane interior: chain unsaturation (which effectively decouples the two hydrocarbon segments disjoined by a double bond) or chain asymmetry (which causes the terminal hydrocarbon segments to lose close contact) both affect the bilayer chain-melting phase-transition temperature comparably on the effective chain-length scale. Thermodynamic consequences of the trans unsaturation are approximately 50% smaller than the effects of the double bonds in the cis conformation, owing to the smaller membrane perturbation by the former double bonds. A simple quantitative model is introduced for the analysis of the phospholipid chain-melting phase behavior. This new model permits quantitative predictions of the chain-melting transition temperature solely on the basis of the known lipid chemical composition. It also explains lipid sensitivity to the hydrocarbon type and attachment. The model agreement with the experimental data is usually better than to within 99% and thus comparable to experimental scatter, even when only a few or no adjustable parameters are used. The membrane fluidization temperature is calculated for a number of potentially interesting, also as yet unexplored, phospholipids, and the biological significance of the effective chain-length concept is discussed.     

Effects of various numbers and positions of cis double bonds in the sn-2 acyl chain of phosphatidylethanolamine on the chain-melting temperature

In an attempt to investigate systematically the effects of various single and multiple cis carbon-carbon double bonds in the sn-2 acyl chains of natural phospholipids on membrane properties, we have de novo synthesized unsaturated C20 fatty acids comprised of single or multiple methylene-interrupted cis double bonds. Subsequently, 15 molecular species of phosphatidylethanolamine (PE) with sn-1 C20-saturated and sn-2 C20-unsaturated acyl chains were semi-synthesized by acylation of C20-lysophosphatidylcholine with unsaturated C20 fatty acids followed by phospholipase D-catalyzed base-exchange reaction in the presence of excess ethanolamine. The gel-to-liquid crystalline phase transitions of these 15 mixed-chain PE, in excess H2O, were investigated by high resolution differential scanning calorimetry. In addition, the energy-minimized structures of these sn-1 C20-saturated/sn-2 C20-unsaturated PE were simulated by molecular mechanics calculations. It is shown that the successive introduction of cis double bonds into the sn-2 acyl chain of C(20):C(20)PE can affect the gel-to-liquid crystalline phase transition temperature, Tm, of the lipid bilayer in some characteristic ways; moreover, the effect depends critically on the position of cis double bonds in the sn-2 acyl chain. Specifically, we have constructed a novel Tm diagram for the 15 species of unsaturated PE, from which the effects of the number and the position of cis double bonds on Tm can be examined simultaneously in a simple, direct, and unifying manner. Interestingly, the characteristic Tm profiles exhibited by different series of mixed-chain PE with increasing degree of unsaturation can be interpreted in terms of structural changes associated with acyl chain unsaturation.     

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.     

Mean-field calculations of chain packing and conformational statistics in lipid bilayers: comparison with experiments and molecular dynamics studies.

A molecular, mean-field theory of chain packing statistics in aggregates of amphiphilic molecules is applied to calculate the conformational properties of the lipid chains comprising the hydrophobic cores of dipalmitoyl-phosphatidylcholine (DPPC), dioleoyl-phosphatidylcholine (DOPC), and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers in their fluid state. The central quantity in this theory, the probability distribution of chain conformations, is evaluated by minimizing the free energy of the bilayer assuming only that the segment density within the hydrophobic region is uniform (liquidlike). Using this distribution we calculate chain conformational properties such as bond orientational order parameters and spatial distributions of the various chain segments. The lipid chains, both the saturated palmitoyl (-(CH2)14-CH3) and the unsaturated oleoyl (-(CH2)7-CH = CH-(CH2)7-CH3) chains are modeled using rotational isomeric state schemes. All possible chain conformations are enumerated and their statistical weights are determined by the self-consistency equations expressing the condition of uniform density. The hydrophobic core of the DPPC bilayer is treated as composed of single (palmitoyl) chain amphiphiles, i.e., the interactions between chains originating from the same lipid headgroup are assumed to be the same as those between chains belonging to different molecules. Similarly, the DOPC system is treated as a bilayer of oleoyl chains. The POPC bilayer is modeled as an equimolar mixture of palmitoyl and oleoyl chains. Bond orientational order parameter profiles, and segment spatial distributions are calculated for the three systems above, for several values of the bilayer thickness (or, equivalently, average area/headgroup) chosen, where possible, so as to allow for comparisons with available experimental data and/or molecular dynamics simulations. In most cases the agreement between the mean-field calculations, which are relatively easy to perform, and the experimental and simulation data is very good, supporting their use as an efficient tool for analyzing a variety of systems subject to varying conditions (e.g., bilayers of different compositions or thicknesses at different temperatures).     

A measure of conformational entropy change during thermal protein unfolding using neutron spectroscopy

Thermal unfolding of proteins at high temperatures is caused by a strong increase of the entropy change which lowers Gibbs free energy change of the unfolding transition (DeltaG(unf) = DeltaH - TDeltaS). The main contributions to entropy are the conformational entropy of the polypeptide chain itself and ordering of water molecules around hydrophobic side chains of the protein. To elucidate the role of conformational entropy upon thermal unfolding in more detail, conformational dynamics in the time regime of picoseconds was investigated with neutron spectroscopy. Confined internal structural fluctuations were analyzed for alpha-amylase in the folded and the unfolded state as a function of temperature. A strong difference in structural fluctuations between the folded and the unfolded state was observed at 30 degrees C, which increased even more with rising temperatures. A simple analytical model was used to quantify the differences of the conformational space explored by the observed protein dynamics for the folded and unfolded state. Conformational entropy changes, calculated on the basis of the applied model, show a significant increase upon heating. In contrast to indirect estimates, which proposed a temperature independent conformational entropy change, the measurements presented here, demonstrated that the conformational entropy change increases with rising temperature and therefore contributes to thermal unfolding.     

Effects of hydrated water on protein unfolding

The conformational stability of a protein in aqueous solution is described in terms of the thermodynamic properties such as unfolding Gibbs free energy, which is the difference in the free energy (Gibbs function) between the native and random conformations in solution. The properties are composed of two contributions, one from enthalpy due to intramolecular interactions among constituent atoms and chain entropy of the backbone and side chains, and the other from the hydrated water around a protein molecule. The hydration free energy and enthalpy at a given temperature for a protein of known three-dimensional structure can be calculated from the accessible surface areas of constituent atoms according to a method developed recently. Since the hydration free energy and enthalpy for random conformations are computed from those for an extended conformation, the thermodynamic properties of unfolding are evaluated quantitatively. The evaluated hydration properties for proteins of known transition temperature (Tm) and unfolding enthalpy (delta Hm) show an approximately linear dependence on the number of constituent heavy atoms. Since the unfolding free energy is zero at Tm, the enthalpy originating from interatomic interactions of a polypeptide chain and the chain entropy are evaluated from an experimental value of delta Hm and computed properties due to the hydrated water around the molecule at Tm. The chain enthalpy and entropy thus estimated are largely compensated by the hydration enthalpy and entropy, respectively, making the unfolding free energy and enthalpy relatively small. The computed temperature dependences of the unfolding free energy and enthalpy for RNase A, T4 lysozyme, and myoglobin showed a good agreement with the experimental ones.(ABSTRACT TRUNCATED AT 250 WORDS)     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

153 Wonderful roles of the entropy in protein dynamics,binding and folding

The entropy, which is central to the second law of thermodynamics, determines that the thermal energy always flows spontaneously from regions of higher temperature to regions of lower temperature. In the protein–solvent thermodynamic system, the entropy is defined as a measure of how evenly the thermal energy would distribute over the entire system (Liu et al., 2012). Such tendency to distribute energy as evenly as possible will reduce the state of order of the initial system, and hence, the entropy can be regarded as an expression of the disorder, or randomness of the system (Yang et al., 2012). For a protein–solvent system under a constant solvent condition, the origin of entropy is the thermal energy stored in atoms, which makes atoms jostle around and bump onto one another, thus leading to vibrations of the covalent bonds connecting two atoms (occurring on the fs timescale) and the rotational and translational motions of amino acid side chain groups (occurring on ps timescale) and water molecules. These motions break the noncovalent bonds around structural regions that are weakly constrained thereby triggering the competitive interactions among residues or between residues and water molecules leading ultimately to the loop motions (occurring on ns timescale) around the protein surface. The loop motions can further transmit either through the water network around the protein surface or via specific structural components (such as the hinge-bending regions) over the entire protein molecule leading to large concerted motions (occurring on μs to s timescales) that are most relevant to protein functions (Amadei, Linssen & Berendsen, 1993; Tao, Rao & Liu, 2010). Thus, the multiple hierarchies of the protein dynamics on distinct timescales (Henzler-Wildman & Kern, 2007) are a consequence of the cascade amplification of the microscopic motions of atoms and groups for which the entropy originating from atomic thermal energy is most fundamental. In the case of protein–ligand binding, the importance of the entropy is embodied in the following aspects. (i) The release of the water molecule kinetic energy (which is a process of the solvent entropy maximization) will cause Brownian motions of individual water molecules which result in strong Brownian bombardments to solute molecules causing molecule wanders/diffusions and subsequent accident contacts/collisions between proteins and ligands. (ii) Such collisions will inevitably cause water molecule displacement and, if the contact interfaces are properly complementary, the requirement to increase the solvent entropy would further displace the water network around the binding interfaces thus leading to the formation the initial protein-ligand complex. (iii) In the initial complex, the loose association of the two partners provide the opportunity for protein to increase conformational entropy, thus triggering the conformational adjustments through competitive interaction between protein residues and ligand, leading ultimately to the formation of tightly associated complex (Liu et al., 2012). In the protein folding process, the first stage, i.e. the rapid hydrophobic collapse (Agashe, Shastry & Udgaonkar, 1995; Dill, 1985), is in fact driven by the effect of the solvent entropy maximization. Specifically, the requirement to maintain as many as possible the dynamic hydrogen bonds among the water molecules will squeeze/sequestrate the hydrophobic amino acid side chains into the interior of the folding intermediates and expose the polar/charged side chains onto the intermediate surface. This will minimize the solvent accessible surface area of the folding intermediates and as thus maximize the entropy of the solvent. The resulting molten globule states (Ohgushi & Wada, 1983) may contain a few secondary structural components and native tertiary contacts, while many native contacts, or close residue–residue interactions present in the native state have not yet formed. However, the nature to increase the protein conformational entropy can trigger a further conformational adjustment process, i.e. the conformational entropy increase breaks the transient secondary or tertiary contacts and triggers the competitive interactions among protein residues and between residues and water. This process may repeat many rounds until the negative enthalpy change resulting from the noncovalent formations can overcompensate for protein conformational entropy loss. In summary, we consider that the tendency to maximize the entropy of the protein–solvent system, which originates from the atomic thermal energy, is the most fundamental driving factor for protein folding, binding, and dynamics, whereas the enthalpy reduction, an opposing factor that tends to make the system become ordered, can compensate for the effect of entropy loss to ultimately allow the system to reach equilibrium at the free energy minima, either global or local.     

Conformational stability and activity of ribonuclease T1 with zero, one, and two intact disulfide bonds

Ribonuclease T1 has two disulfide bonds linking cysteine residues 2-10 and 6-103. We have prepared a derivative of ribonuclease T1 in which one disulfide bond is broken and the cysteine residues carboxymethylated, (2-10)-RCM-T1, and three derivatives in which both disulfides are broken and the cysteine residues reduced, R-T1, carboxamidomethylated, RCAM-T1, or carboxymethylated, RCM-T1. The RNA hydrolyzing activity of these proteins has been measured, and urea and thermal denaturation studies have been used to determine conformational stability. The activity, melting temperature, and conformational stability of the proteins are: ribonuclease T1 (100%, 59.3 degrees C, 10.2 kcal/mol), (2-10)-RCM-T1 (86%, 53.3 degrees C, 6.8 kcal/mol), R-T1 (53%, 27.2 degrees C, 3.0 kcal/mol), RCAM-T1 (43%, 21.2 degrees C, 1.5 kcal/mol), and RCM-T1 (35%, 16.6 degrees C, 0.9 kcal/mol). Thus, the conformational stability is decreased by 3.4 kcal/mol when one disulfide bond is broken and by 7.2-9.3 kcal/mol when both disulfide bonds are broken. It is quite remarkable that RNase T1 can fold and function with both disulfide bonds broken and the cysteine residues carboxymethylated. The large decrease in the stability is due mainly to an increase in the conformational entropy of the unfolded protein which results when the constraints of the disulfide bonds on the flexibility are removed. We propose a new equation for predicting the effect of a cross-link on the conformational entropy of a protein: delta Sconf = -2.1 - (3/2)R 1n n, where n is the number of residues between the side chains which are cross-linked. This equation gives much better agreement with experimental results than other forms of this equation which have been used previously.     

On the conformational, physical properties and functions of polyunsaturated acyl chains

The conformational properties of the acyls of biological membranes--hydrocarbon chains with isolated cis double bonds--were studied by computer simulation. The Monte Carlo method was used, with continuous variation of bond rotation angles within the (0, 360 degree) range considered. It has been shown, that if all double bonds of molecules are separated only by one methylene group, and their number in the chain is maximum, the molecule is characterized by the highest equilibrium flexibility (at temperatures only encountered by biological systems) as compared to any similar molecules. It is such a structure which is inherent to docosahexaenoic acid. The above molecule coefficient that characterizes the temperature sensitivity of the molecule sizes is 10-times lower than that of a saturated chain. The polyunsaturated chain segment with high probability assumes the extended (in perfect crystal structures the 'angle iron-shaped') conformation when all the molecules are efficiently packed below the phase-transition temperatures. The annular lipid layer of embedded enzymes is assumed to be enriched with polyunsaturated fatty acid acyls. The above physical properties of polyunsaturated chains are bound to favour the maintenance of the proper conformational mobility of biomembrane enzymes, to relax the negative influence of environmental temperature changes on their activity. When freezing biological membranes they are bound to provide the molecule packing which is free of high tensions.     

High state of order of isolated bacterial lipopolysaccharide and its possible contribution to the permeation barrier property of the outer membrane.

The conformational properties of the isolated S form of Salmonella sp. lipopolysaccharide (LPS), of Re mutant LPS, and of free lipid A were investigated by using X-ray diffraction and conformational energy calculations. The data obtained showed that LPS in a dried, in a hydrated, and probably also in an aqueous dispersion state is capable of forming bilayered lamellar arrangements similar to phospholipids. From the bilayer packing periodicities, a geometrical model of the extensions of the LPS regions lipid A, 2-keto-3-deoxyoctulosonic acid, and O-specific chain along the membrane normal could be calculated. Furthermore, the lipid A component was found to assume a remarkably high ordered conformation: its fatty acid chains were tightly packed in a dense hexagonal lattice with a center-to-center distance of 0.49 nm. The hydrophilic backbone of lipid A showed a strong tendency to form domains in the membrane, resulting in a more or less parallel arrangement of lipid A units. According to model calculations, the hydrophilic backbone of lipid A appears to be oriented approximately 45 degrees to the membrane surface, which would lead to a shed roof-like appearance of the surface structure in the indentations of which the 2-keto-3-deoxyoctulosonic acid moiety would fit. In contrast, the O-specific chains assume a low ordered, heavily coiled conformation. Comparison of these structural properties with those known for natural phospholipids in biological membranes indicates that the high state of order of the lipid A portion of LPS might be an important factor in the structural role and permeation barrier functions of LPS in the outer membrane of gram-negative bacteria.     

Affinity of phosphatidylcholine molecular species for the bovine phosphatidylcholine and phosphatidylinositol transfer proteins. Properties of the sn-1 and sn-2 acyl binding sites

Both the phosphatidylcholine transfer protein (PC-TP) and the phosphatidylinositol transfer protein (PI-TP) act as carriers of phosphatidylcholine (PC) molecules between membranes. To study the structure of the acyl binding sites of these proteins, the affinity of 32 distinct natural and related PC molecular species was determined by using a previously developed fluorometric competition assay. Marked differences in affinity between species were observed with both proteins. Affinity vs lipid hydrophobicity (determined by reverse-phase HPLC) plots displayed a well-defined maximum indicating that the acyl chain hydrophobicity is an important determinant of binding of a phospholipid molecule by these transfer proteins. However, besides the overall lipid hydrophobicity, steric properties of the individual acyl chains contribute considerably to the affinity, and PC-TP and PI-TP respond differently to modifications of the acyl chain structure. The affinity of PC-TP increased steadily with increasing unsaturation of the sn-2 acyl moiety, resulting in high affinity for species containing four and six double bonds in the sn-2 chain, whereas the affinity of PI-TP first increased up to two to three double bonds and then declined. These data, as well as the distinct effects of sn-2 chain double bond position and bromination, indicate that the sn-2 acyl chain binding sites of the two proteins are structurally quite different. The sn-1 acyl binding sites are dissimilar as well, since variation of the length of saturated sn-1 chain affected the affinity differently. The data are discussed in terms of the structural organization of the sn-1 and sn-2 acyl binding sites of PC-TP and PI-TP.(ABSTRACT TRUNCATED AT 250 WORDS)