Extension of the fragment method to calculate amino acid zwitterion and side chain partition coefficients

The fragment method of calculating partition coefficients (P) has been extended to include the common amino acids (AAs). The results indicate that polar and charged side chains influence the hydrophobicity of atoms in the side chain in a predictable manner. Field effects, as evidenced through polar proximity factors and bond factors, need to be considered for accurate estimation of transfer phenomena. The calculated log P and delta G degree ' values of the 20 AAs agree well with the observed values. Pro calculates to be more hydrophilic than the observed log P. Hydrophobicity scales for peptide side chain residues are compared and evaluated in terms of suitability. Calculated pi values for nonpolar side chain residues agree well with the observed values; calculated values for uncharged polar side chain residues deviate by about 0.6 log units except for Gln and Cys; and polar side chain residues with charged side chains calculate as too hydrophilic. Reasons for the differences are explored. We also suggest that tightly bound water to polar moieties in amino acids and peptides may be transferred into the octanol phase during partitioning experiments. A quantitative methodology is presented which characterizes the thermodynamic partitioning of groups and individual atoms in amino acids and proteins.     

Development of an isotope dilution mass spectrometry assay for HbA1c based on enzyme-cleaved peptide analysis

The hydrophobic contributions of 17 individual peptides, fused to the N-terminal of Bacillus stearothermophilus lactate dehydrogenase (LDH) were studied by hydrophobic interaction chromatography (HIC) and aqueous two-phase system (ATPS). The constructs were sequenced from a protein library designed with a five-amino acid randomised region in the N-terminal of an LDH protein. The 17 LDH variants and an LDH control lacking the randomised region were expressed in Escherichia coli. HIC and ATPS behaviour of the proteins indicated significant differences in protein hydrophobicity, even though the modifications caused only 1% increase in protein molecular weight and 2% variation in isoelectric points. HIC and ATPS results correlated well (R(2) = 0.89). Protein expression was clearly affected by N-terminal modification, but there was no evidence that the modification affected protein activity. A GluAsnAlaAspVal modification resulted in increased protein expression. In most cases, HIC and ATPS results compared favourably with those predicted on the basis of 34 amino acid residue hydrophobicity scales; assuming exposure of tag residues to solution. Exceptions included LeuAlaGlyValIle and LeuTyrGlyCysIle modifications, which were predicted, assuming full solution exposure, to be more hydrophobic than observed.     

Immobilized lipase Candida sp. 99-125 catalyzed methanolysis of glycerol trioleate: solvent effect

The immobilized lipase Candida sp. 99-125 catalyzed methanolysis of glycerol trioleate was studied in twelve different solvents in order to deduce the solvent effect through an attempt to correlate the highest yield with such solvent properties as hydrophobicity (log P), dielectric constant (epsilon), and Hildebrand solubility parameter (delta). The results showed that the conversion of glycerol trioleate and yield of oleic acid methyl ester were quite dependent on the solvent. The catalyst lipase in various solvents also needed different optimum amount of water to keep its maximum activity, and generally this lipase in more hydrophobic solvents required more water. The correlation between the highest yield and log P value was found to be reasonable except deviation of data points of certain solvents, while no obvious correlation existed between the other two parameters, dielectric constant (epsilon) and Hildebrand solubility parameter (delta), and the enzyme activity. The study revealed that more hydrophobic solvents such as n-hexane or cyclohexane were more suitable solvents for Candida sp. 99-125 catalyzed transesterification of glycerol trioleate to oleic acid methyl ester.     

Statistical distribution of hydrophobic residues along the length of protein chains. Implications for protein folding and evolution.

We consider in this paper the statistical distribution of hydrophobic residues along the length of protein chains. For this purpose we used a binary hydrophobicity scale which assigns hydrophobic residues a value of one and non-hydrophobes a value of zero. The resulting binary sequences are tested for randomness using the standard run test. For the majority of the 5,247 proteins examined, the distribution of hydrophobic residues along a sequence cannot be distinguished from that expected for a random distribution. This suggests that (a) functional proteins may have originated from random sequences, (b) the folding of proteins into compact structures may be much more permissive with less sequence specificity than previously thought, and (c) the clusters of hydrophobic residues along chains which are revealed by hydrophobicity plots are a natural consequence of a random distribution and can be conveniently described by binomial statistics.     

Folding patterns of porin and bacteriorhodopsin.

Porin spans the outer membrane of Escherichia coli with most of the protein embedded within the membrane. It lacks pronounced hydrophobic domains and consists predominantly of beta-pleated sheet. These observations require the accommodation of polar and ionizable residues in an environment that has a low dielectric constant. Owing to a currently limited understanding of the constraints governing membrane protein structure, a minimal approach to structure prediction is proposed that identifies segments causing polypeptides to reverse their direction (turn identification). The application of this procedure avoids hydrophobicity parameters and yields a model of porin which is in good agreement with all experimental data available. The presence of polar and ionizable residues within membrane boundaries implies a dense (saturating) network of hydrogen bond donor and acceptor groups. Application to a paradigm of hydrophobic membrane proteins, bacteriorhodopsin, reveals a pattern consistent with its alpha-helical folding. The postulated structure includes significantly more polar residues in the membrane domain than have been assumed previously, suggesting that there are also hydrogen bonding networks in bacteriorhodopsin. Extensive networks permeating protein interior and surfaces would explain the extraordinary stability and the tight interactions between functional units in the formation of crystalline arrays of both proteins.     

Use of chemically modified proteins to study the effect of a single protein property on partitioning in aqueous two-phase systems: Effect of surface hydrophobicity

Two different series of hydrophobically modified proteins were partitioned in a number of aqueous two-phase systems (ATPS) to investigate the effect of hydrophobicity as a single property on partitioning. The modified proteins were derived from beta-lactoglobulin and bovine serum albumin (BSA). Measurement of the surface hydrophobicity of the proteins is important; hydrophobic interaction chromatography (HIC) was used for this purpose. The resolution of the systems (R) in terms of protein surface hydrophobicity and the intrinsic hydrophobicity (log P(0)) of the systems was established. The effect of the addition of NaCl to PEG/phosphate and PEG/dextran systems was analyzed in terms of the hydrophobicity difference between the phases and their ability to promote hydrophobic interactions between the protein surface and the PEG molecules. The values for R and log P(0) differed somewhat depending on which group of modified proteins was used for partitioning. The addition of NaCl to PEG/phosphate systems promoted an increase in the values of R, showing an important effect on the resolution of the systems for protein surface hydrophobicity (twice as high when compared with systems without NaCl). For PEG/dextran systems, the addition of 9% NaCl (w/w) promoted an improvement in the resolution toward surface hydrophobicity with an increase of 60% on the value of R. (c) 1996 John Wiley & Sons, Inc.     

Regeneration and Fusion of Mycelial Protoplasts of Tricholoma matsutake

The Hill-reaction inhibitory activity of 30 derivatives of 2-difluoromethylthio-1,3,5-triazine was determined using chlorella. The structure-activity relationships were analyzed using the hydrophobic parameter (log P), the electronic substituent constant (c) and the steric substituent constant (Es). The Hill-reaction inhibitory activity showed a parabolic relationship to log P and the steric substituent parameters at 4,6-dialkylamino- positions in the triazine ring. The equation revealed that optimum hydrophobicity and size of the substituents were necessary for high activity.     

Characteristic features of amino acid residues in coiled-coil protein structures

Detailed analyses of protein structures provide an opportunity to understand conformation and function in terms of amino acid sequence and composition. In this work, we have systematically analyzed the characteristic features of the amino acid residues found in alpha-helical coiled-coils and, in so doing, have developed indices for their properties, conformational parameters, surrounding hydrophobicity and flexibility. As expected, there is preference for hydrophobic (Ala, Leu), positive (Lys, Arg) and negatively (Glu) charged residues in coiled-coil domains. However, the surrounding hydrophobicity of residues in coiled-coil domains is significantly less than that for residues in other regions of coiled-coil proteins. The analysis of temperature factors in coiled-coil proteins shows that the residues in these domains are more stable than those in other regions. Further, we have delineated the medium- and long-range contacts in coiled-coil domains and compared the results with those obtained for other (non-coiled-coil) parts of the same proteins and non-coiled-coil helical segments of globular proteins. The residues in coiled-coil domains are largely influenced by medium-range contacts, whereas long-range interactions play a dominant role in other regions of these same proteins as well as in non-coiled-coil helices. We have also revealed the preference of amino acid residues to form cation-pi interactions and we found that Arg is more likely to form such interactions than Lys. The parameters developed in this work can be used to understand the folding and stability of coiled-coil proteins in general.     

Stabilization of proteins by enhancement of inter-residue hydrophobic contacts: lessons of T4 lysozyme and barnase

Although the hydrophobic interactions are considered as the main contributors to the protein stability, not much examples of protein stabilization by rational increasing of this type of interactions still can be found in literature. This is partly due to the lack of proper theoretical "measure" of hydrophobic interactions and their changes upon mutations. In the present paper the molecular hydrophobicity potential approach is used to assess how the changes in type and the strength of inter-residue contacts upon single amino acid mutations are correlated with the changes in thermodynamic stability of T4 lysozyme and barnase mutants, and which factors affect these correlations. Mutations changing unfavorable hydrophilic-to-hydrophobic contacts into favorable hydrophobic were found to enhance the thermodynamic stability in more than 81 % of cases, if these mutations do not create steric bumps and do not involve proline residues and hydrogen-bonded side-chains. Mutations increasing hydrophobic contributions (according to molecular hydrophobicity potential formalism) lead to increase of thermodynamic stability in more than 94% of cases for certain type of mutations (i.e., mutations not involving charged residues, Pro and residues with side-chain hydrogen bonds, when these mutations do not introduce steric bumps and do not involve strongly exposed residues and residues situated at helix N- and C-cap positions). For this type of mutations the correlation was found between the change in hydrophobic contributions of mutated residues deltaCphob and thermodynamic parameters deltaTm (change in melting temperature) and deltadeltaG (change in free energy of unfolding). Although the correlation coefficients were larger if the experimental structures of mutants were used for the calculations (correlation coefficients r(exp) deltaC,deltaT = .85 and r(exp) deltaC,deltadeltaG = 0.87) than if the modeled structures were used instead (r(mod) deltaC,deltaT = 0.74 and r(mod)deltaC,deltadeltaG = 0.76), the modelled structures of mutants in the vast majority of cases can be used for qualitative predictition of the protein stabilization. Basing on the analysis of mutations increasing hydrophobic contributions in T4 lysozyme the substitution matrix was derived, which can be used to decide which new residue should be put instead the old one to increase the stability of protein. The estimation shows that the number of potential mutation sites for enhancement of hydrophobic interactions in T4 lysozyme is quite large, and only approximately 10 per cent of them were studied thus far. Basing on the current analysis of T4 lysozyme and barnase mutations the algorithm for increasing of protein stability via increasing of hydrophobic interactions for the proteins with known spatial structure is proposed.     

Clustering of large hydrophobes in the hydrophobic core of two-stranded alpha-helical coiled-coils controls protein folding and stability

The de novo design and biophysical characterization of two 60-residue peptides that dimerize to fold as parallel coiled-coils with different hydrophobic core clustering is described. Our goal was to investigate whether designing coiled-coils with identical hydrophobicity but with different hydrophobic clustering of non-polar core residues (each contained 6 Leu, 3 Ile, and 7 Ala residues in the hydrophobic core) would affect helical content and protein stability. The disulfide-bridged P3 and P2 differed dramatically in alpha-helical structure in benign conditions. P3 with three hydrophobic clusters was 98% alpha-helical, whereas P2 was only 65% alpha-helical. The stability profiles of these two analogs were compared, and the enthalpy and heat capacity changes upon denaturation were determined by measuring the temperature dependence by circular dichroism spectroscopy and confirmed by differential scanning calorimetry. The results showed that P3 assembled into a stable alpha-helical two-stranded coiled-coil and exhibited a native protein-like cooperative two-state transition in thermal melting, chemical denaturation, and calorimetry experiments. Although both peptides have identical inherent hydrophobicity (the hydrophobic burial of identical non-polar residues in equivalent heptad coiled-coil positions), we found that the context dependence of an additional hydrophobic cluster dramatically increased stability of P3 (Delta Tm approximately equal to 18 degrees C and Delta[urea](1/2) approximately equal to 1.5 M) as compared with P2. These results suggested that hydrophobic clustering significantly stabilized the coiled-coil structure and may explain how long fibrous proteins like tropomyosin maintain chain integrity while accommodating polar or charged residues in regions of the protein hydrophobic core.     

Cellular response mechanisms in Pseudomonas aeruginosa PseA during growth in organic solvents

Aims:  Solvent-tolerant bacteria have emerged as a new class of micro-organisms able to grow at high concentrations of toxic solvents. Such bacteria and their solvent-stable enzymes are perceived to be useful for biotransformations in nonaqueous media. In the present study, the solvent-responsive features of a lipase–producing, solvent-tolerant strain Pseudomonas aeruginosa PseA have been investigated to understand the cellular mechanisms followed under solvent-rich conditions.
Methods and Results:  The solvents, cyclohexane and tetradecane with differing log P -values (3·2 and 7·6 respectively), have been used as model systems. Effect of solvents on (i) the cell morphology and structure (ii) surface hydrophobicity and (iii) permeability of cell membrane have been examined using transmission electron microscopy, atomic force microscopy and other biochemical techniques. The results show that (i) less hydrophobic (low log P -value) solvent cyclohexane alters the cell membrane integrity and (ii) cells adapt to organic solvents by changing morphology, size, permeability and surface hydrophobicity. However, no such changes were observed in the cells grown in tetradecane.
Conclusions:  It may be concluded that P. aeruginosa PseA responds differently to solvents of different hydrophobicities. Bacterial cell membrane is more permeable to less hydrophobic solvents that eventually accumulate in the cytoplasm, while highly hydrophobic solvents have lesser tendency to access the membrane.
Significance and Impact of the Study:  To the best of our knowledge, these are first time observations that show that way of bacterial solvent adaptability depends on nature of solvent. Difference in cellular responses towards solvents of varying log P -values (hydrophobicity) might prove useful to search for a suitable solvent for carrying out whole-cell biocatalysis.     

Fine structure analysis of a protein folding transition state; distinguishing between hydrophobic stabilization and specific packing

Developing a detailed understanding of the structure and energetics of protein folding transition states is a key step in describing the folding process. The phi-value analysis approach allows the energetic contribution of side-chains to be mapped out by comparing wild-type with individual mutants where conservative changes are introduced. Studies where multiple substitutions are made at individual sites are much rarer but are potentially very useful for understanding the contribution of each element of a side-chain to transition state formation, and for distinguishing the relative importance of specific packing versus hydrophobic interactions. We have made a series of conservative mutations at multiple buried sites in the N-terminal domain of L9 in order to assess the relative importance of specific side-chain packing versus less specific hydrophobic stabilization of the transition state. A total of 28 variants were prepared using both naturally occurring and non-naturally occurring amino acids at six sites. Analysis of the mutants by NMR and CD showed no perturbation of the structure. There is no correlation between changes in hydrophobicity and changes in stability. In contrast, there is excellent linear correlation between the hydrophobicity of a side-chain and the log of the folding rate, ln(k(f)). The correlation between ln(k(f)) and the change in hydrophobicity holds even for substitutions that change the shape and/or size of a side-chain significantly. For most sites, the correlation with the logarithm of the unfolding rate, ln(k(u)), is much worse. Mutants with more hydrophobic amino acid substitutions fold faster, and those with less hydrophobic amino acid substitutions fold slower. The results show that hydrophobic interactions amongst core residues are an important driving force for forming the transition state, and are more important than specific tight packing interactions. Finally, a number of substitutions lead to negative phi-values and the origin of these effects are described.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.     

QSAR study on phenolic activity: need of positive hydrophobic term (log P) in QSAR

The phenolic activity (log 1/C) of a large series of phenols against L1210 leukaemia cells was modelled using physicochemical parameters other than conventional electronic and steric parameters. Attempts have also been made to examine need or otherwise of the hydrophobic parameter, log P, in such studies. The results have shown that contribution of log P in modelling log 1/C is favourable.     

Magnitude of the hydrophobic effect at central versus peripheral sites in protein-protein interfaces

Hydrophobic interactions are essential for stabilizing protein-protein complexes, whose interfaces generally consist of a central cluster of hot spot residues surrounded by less important peripheral residues. According to the O-ring hypothesis, a condition for high affinity binding is solvent exclusion from interacting residues. This hypothesis predicts that the hydrophobicity at the center is significantly greater than at the periphery, which we estimated at 21 cal mol(-1) A(-2). To measure the hydrophobicity at the center, structures of an antigen-antibody complex where a buried phenylalanine was replaced by smaller hydrophobic residues were determined. By correlating structural changes with binding free energies, we estimate the hydrophobicity at this central site to be 46 cal mol(-1) A(-2), twice that at the periphery. This context dependence of the hydrophobic effect explains the clustering of hot spots at interface centers and has implications for hot spot prediction and the design of small molecule inhibitors.     

QSAR study on permeability of hydrophobic compounds with artificial membranes

We previously reported a classical quantitative structure-activity relationship (QSAR) equation for permeability coefficients (P(app-pampa)) by parallel artificial membrane permeation assay (PAMPA) of structurally diverse compounds with simple physicochemical parameters, hydrophobicity at a particular pH (logP(oct) and |pK(a)-pH|), hydrogen-accepting ability (SA(HA)), and hydrogen-donating ability (SA(HD)); however, desipramine, imipramine, and testosterone, which have high logP(oct) values, were excluded from the derived QSAR equation because their measured P(app-pampa) values were lower than calculated. In this study, for further investigation of PAMPA permeability of hydrophobic compounds, we experimentally measured the P(app-pampa) of more compounds with high hydrophobicity, including several pesticides, and compared the measured P(app-pampa) values with those calculated from the QSAR equation. As a result, compounds having a calculated logP(app-pampa)>-4.5 showed lower measured logP(app-pampa) than calculated because of the barrier of the unstirred water layer and the membrane retention of hydrophobic compounds. The bilinear QSAR model explained the PAMPA permeability of the whole dataset of compounds, whether hydrophilic or hydrophobic, with the same parameters as the equation in the previous study. In addition, PAMPA permeability coefficients correlated well with Caco-2 cell permeability coefficients. Since Caco-2 cell permeability is effective for the evaluation of human oral absorption of compounds, the proposed bilinear model for PAMPA permeability could be useful for not only effective screening for several drug candidates but also the risk assessment of chemicals and agrochemicals absorbed by humans.     

Determination of the hydrophobicity of local anesthetic agents

Hydrophobicity, a term used to describe a fundamental physicochemical property of local anesthetics, was in the past obtained by octanol/buffer partitioning. It has been suggested that the octanol method, despite its obvious advantages, also has some drawbacks. HPLC has become an attractive alternative for the measurement of hydrophobicity and has been applied to local anesthetics recently. However, the methods in current use for measuring the hydrophobicity of local anesthetics suffer from a number of limitations and remain obscure. This study introduces a new HPLC method for measuring the hydrophobicity of eight local anesthetics in current clinical use. Using a C(18) derivatized polystyrene-divinylbenzene stationary phase HPLC column, the log k'(w) values of local anesthetics were determined by measuring the capacity factor k'(i) in the process of chromatographic separation using a hydrophobic stationary phase and a hydrophilic mobile phase. A rapid reversed-phase HPLC method was developed to directly measure log k'(w) of eight local anesthetics. A high correlation between log k'(w) and hydrophobicity (log P(oct)) from the traditional shake-flask method was obtained for the local anesthetics, demonstrating the reliability of the method. The results reveal an improved method for measuring the hydrophobicity of the local anesthetic agents in the unionized form. This simple, sensitive and reproducible approach may serve as a valuable tool for describing the physicochemical properties of novel local anesthetics.     

Immiscible organic solvent inactivation of urease, chymotrypsin, lipase, and ribonuclease: separation of dissolved solvent and interfacial effects

A new technique with controlled interface generation allows separation and quantitation of enzyme inactivation by both solvent/aqueous interface and dissolved solvent. This has now been used in n-butanol, isopropylether, 2-octanone, n-hexane, n-butylbenzene, and n-tridecane. Ribonuclease was stable with all the solvent/aqueous interfaces studied. Chymotrypsin was mainly inactivated by the more hydrophobic solvent/aqueous interfaces, whereas lipase was only inactivated by the less hydrophobic solvent/aqueous interfaces. Urease was inactivated by some interfaces, but not all, without an obvious trend. Thus, the commonly expected simple relationship with solvent polarity (e.g., log P) does not apply when interfacial inactivation is determined specifically. Greater dissolved solvent inactivation occurred with the more polar solvents, though only a general trend was apparent with log P. A better correlation was noted with the Hilde-brand solubility parameter. Interfacial effects are discussed with reference to enzyme molecular weight, denaturation temperature, hydrophobicity, and adiabatic compressibility. (c) 1994 John Wiley & Sons, Inc.