SEARCH · PubMed Health

Results for “Solvent accessibility”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23Linked to original sources

PubMed1997

Hydrophobic folding units derived from dissimilar monomer structures and their interactions.

We have designed an automated procedure to cut a protein into compact hydrophobic folding units. The hydrophobic units are large enough to contain tertiary non-local interactions, reflecting potential nucleation sites during protein folding. The quality of a hydrophobic folding unit is evaluated by four criteria. The first two correspond to visual characterization of a structural domain, namely, compactness and extent of isolation. We use the definition of Zehfus and Rose (Zehfus MH, Rose GD, 1986, Biochemistry 25:35-340) to calculate the compactness of a cut protein unit. The isolation of a unit is based on the solvent accessible surface area (ASA) originally buried in the interior and exposed to the solvent after cutting. The third quantity is the hydrophobicity, equivalent to the fraction of the buried non-polar ASA with respect to the total non-polar ASA. The last criterion in the evaluation of a folding unit is the number of segments it includes. To conform with the rationale of obtaining hydrophobic units, which may relate to early folding events, the hydrophobic interactions are implicitly and explicitly applied in their generation and assessment. We follow Holm and Sander (Holm L, Sander C, 1994, Proteins 19:256-268) to reduce the multiple cutting-point problem to a one-dimensional search for all reasonable trial cuts. However, as here we focus on the hydrophobic cores, the contact matrix used to obtain the first non-trivial eigenvector contains only hydrophobic contracts, rather than all, hydrophobic and hydrophilic, interactions. This dataset of hydrophobic folding units, derived from structurally dissimilar single chain monomers, is particularly useful for investigations of the mechanism of protein folding. For cases where there are kinetic data, the one or more hydrophobic folding units generated for a protein correlate with the two or with the three-state folding process observed. We carry out extensive amino acid sequence order independent structural comparisons to generate a structurally non-redundant set of hydrophobic folding units for fold recognition and for statistical purposes.

Algorithms↗

PubMed2002

Anti-cooperativity and cooperativity in hydrophobic interactions: Three-body free energy landscapes and comparison with implicit-solvent potential functions for proteins.

Potentials of mean force (PMFs) of three-body hydrophobic association are investigated to gain insight into similar processes in protein folding. Free energy landscapes obtained from explicit simulations of three methanes in water are compared with that predicted by popular implicit-solvent effective potentials for the study of proteins. Explicit-water simulations show that for an extended range of three-methane configurations, hydrophobic association at 25 degrees C under atmospheric pressure is mostly anti-cooperative, that is, less favorable than if the interaction free energies were pairwise additive. Effects of free energy nonadditivity on the kinetic path of association and the temperature dependence of additivity are explored by using a three-methane system and simplified chain models. The prevalence of anti-cooperativity under ambient conditions suggests that driving forces other than hydrophobicity also play critical roles in protein thermodynamic cooperativity. We evaluate the effectiveness of several implicit-solvent potentials in mimicking explicit water simulated three-body PMFs. The favorability of the contact free energy minimum is found to be drastically overestimated by solvent accessible surface area (SASA). Both the SASA and a volume-based Gaussian solvent exclusion model fail to predict the desolvation barrier. However, this barrier is qualitatively captured by the molecular surface area model and a recent "hydrophobic force field." None of the implicit-solvent models tested are accurate for the entire range of three-methane configurations and several other thermodynamic signatures considered.

Animals↗

PubMed1998

Tertiary structures of amyloidogenic and non-amyloidogenic transthyretin variants: new model for amyloid fibril formation.

The most common form of hereditary systemic amyloidosis is familial amyloidotic polyneuropathy associated with single amino acid changes in the plasma protein transthyretin. So far, high resolution structures of only three amyloidogenic variants (Met30, Ser84, Ile122) and one non-amyloidogenic variant (Thr109) have been reported complemented by X-ray fiber diffraction studies and image reconstruction from electron micrographs of amyloid fibrils. To investigate the role of structural factors in this disease, we extended our studies to other transthyretin variants. We report crystallization and structural investigations of three amyloidogenic (Arg10, Ala60, Tyr77) and two non-amyloidogenic variants (Ser6, Met119). The similarity of these structures to normal transthyretin does not give direct clues to the fibril forming process. Since transthyretin amyloid fibrils contain a major fragment starting at position 49, besides the intact molecule, we calculated the solvent accessibility of residue 48. Indeed, all amyloidogenic variants show an increased main chain solvent exposure when compared to normal transthyretin and non-amyloidogenic variants, which can be postulated to result in increased susceptibility to proteolysis. After limited proteolysis, dimers are incapable of reassociation to native tetramers. We present a model for amyloid fibril formation based on formation of fibrils from N-terminal truncated dimers as building blocks.

Amyloid↗

PubMed1988

Modeling the three-dimensional structures of bacterial aminotransferases.

The refined crystallographic structure of the "closed" conformation of chicken mitochondrial aspartate aminotransferase has been used as a template for the construction of models of the two Escherichia coli aminotransferases encoded by the tyrB and aspC genes. The main results are as follows: (1) Only minor changes are required in the coordinates of the backbone atoms to accommodate the large number of substituted side chains. (2) All deletions and insertions required to allow maximum primary sequence alignment are on the solvent-accessible surface. (3) Charged residues are all located on the surface, in contact with solvent, except for certain conserved active site residues. (4) The close packing within the hydrophobic core is maintained. (5) The interactions between the subunits are maintained. (6) Modeling of tyrosine as an external aldimine into the active sites points to several residues that could be involved in determining the substrate specificities of these aminotransferases.

Binding Sites↗

PubMed1996

Hydrophobic distribution and spatial arrangement of amino acid residues in membrane proteins.

The analysis of known three-dimensional structures of membrane proteins provides an opportunity to understand their structure and stability. In this article we analyse the hydrophobic variation of amino acid residues at various ranges in membrane and aqueous parts of membrane proteins. The numerical indices for several properties of amino acid residues in membrane proteins, such as surrounding hydrophobicity, gain in surrounding hydrophobicity, hydrophobic gain ratio, accessible surface area, preference of amino acid residues in the interior and surface parts, solvent accessible reduction ratio and buriedness, were set up. The relative preference of amino acid residues at various positions of membrane proteins were obtained in a very realistic approach.

Amino Acids↗

PubMed2006

Transmembrane domain VII of the human apical sodium-dependent bile acid transporter ASBT (SLC10A2) lines the substrate translocation pathway.

Recent evidence implicating transmembrane (TM) segment 7 of the apical sodium-dependent bile acid transporter (ASBT) in substrate interaction warranted examination of its aqueous accessibility. Therefore, cysteine substitution of 22 consecutive amino acids was performed against a methanethiosulfonate (MTS)-resistant background (C270A). Activity and susceptibility to polar MTS derivatives [(2-aminoethyl)-methanethiosulfonate (MTSEA), [2-(trimethylammonium)ethyl]methanethiosulfonate (MTSET), and methanethiosulfonate ethylsulfonate (MTSES)] of mutants were evaluated in COS-1 cells. Thr289, Tyr293, Gln297, Ala301, Phe307, and Tyr308 represented loss-of-function mutants; furthermore, the measurable residual activities for T289C, Y293C, and A301C (<or=20% control) proved insensitive to MTS treatment. MTSES and MTSET inhibition was confined to residues lining the extracellular half of TM7; amino acids situated deeper within the membrane were unaffected. In contrast, the entire length of TM7 was susceptible to the relatively smaller MTSEA; moreover, MTSEA sensitivity was significantly amended by coapplication with substrates. This selective pattern of modification suggests that the highly conserved lower half of TM7 lies within a water-filled cavity easily accessible from the extracellular milieu, whereas residues approaching the cytosolic/membrane interface reside in pores for which accessibility is modulated by molecular volume. Functionally inactive and MTS-inaccessible residues (T289C, Y293C, Q297C, and A301C) within TM7 may play a structural role critical to transporter function; conversely, MTS-sensitive residues are spatially distinct and may demarcate a face of the TM involved in substrate translocation. In addition, computational analysis of solvent-accessible domains identified five key solvent pockets that predominantly line the hydrophilic face of TM7. Combined, our data suggest that TM7 plays a dominant role in the hASBT translocation process.

Amino Acid Sequence↗

PubMed2002

Characterization of the metal-binding site of bovine growth hormone through site-specific metal-catalyzed oxidation and high-performance liquid chromatography-tandem mass spectrometry.

Metal-catalyzed oxidation was used to identify metal-binding His residues in bovine growth hormone (bGH), which has not been characterized well crystallographically due to a high propensity of bGH to aggregate. bGH was exposed to Cu(2+) and ascorbate (ascorbate/Cu(2+)/O(2)). 2-Oxo-His formation was identified by HPLC-tandem mass spectrometry (MS/MS) analysis of a tryptic digest. Two 2-oxo-His-containing fragments were detected, T2(O) (MH(2+)(2) = 748.8) and T20(O) (MH(+) = 528.3), both masses corresponding to the addition of only one oxygen atom (+16 amu) to the respective native fragments, T2 and T20. T2 contains (20)His and (22)His, and T20 contains (170)His. Quantitative HPLC-MS/MS analysis shows the following order of reactivity: (170)His >> (22)His > (20)His. Solvent-accessible surface area calculations determined (22)His and (170)His to be 26 and 35% solvent exposed, respectively, while (20)His is 65% solvent exposed. The presence of an analogous metal-binding site in human growth hormone, which is located in the hydrophobic core, and our experimental finding that oxidation was greatest for (22)His and (170)His in bGH suggests that (22)His and (170)His of bGH participate in metal binding. This result is supported by a previously predicted tertiary structure of bGH and compared with the location of metal-binding His residues of human growth hormone.

Amino Acid Sequence↗

PubMed2005

Finding gas diffusion pathways in proteins: application to O2 and H2 transport in CpI [FeFe]-hydrogenase and the role of packing defects.

We report on a computational investigation of the passive transport of H2 and O2 between the external solution and the hydrogen-producing active site of CpI [FeFe]-hydrogenase from Clostridium pasteurianum. Two distinct methodologies for studying gas access are discussed and applied: (1) temperature-controlled locally enhanced sampling, and (2) volumetric solvent accessibility maps, providing consistent results. Both methodologies confirm the existence and function of a previously hypothesized pathway and reveal a second major pathway that had not been detected by previous analyses of CpI's static crystal structure. Our results suggest that small hydrophobic molecules, such as H2 and O2, diffusing inside CpI, take advantage of well-defined preexisting packing defects that are not always apparent from the protein's static structure, but that can be predicted from the protein's dynamical motion. Finally, we describe two contrasting modes of intraprotein transport for H2 and O2, which in our model are differentiated only by their size.

Binding Sites↗

PubMed2002

Acetate binding at the photosystem II oxygen evolving complex: an S(2)-state multiline signal ESEEM study.

Previously, using acetate deuterated in the methyl hydrogen positions, we showed that acetate binds in close proximity to the Mn cluster/Y(.)(z) tyrosine dual spin complex in acetate-inhibited photosystem II (PSII) preparations exhibiting the "split" EPR signal arising from the S(2)-Y(.)(z) interaction [Force, D. A.; Randall, D. W.; Britt, R. D. Biochemistry 1997, 36, 12062-12070]. By using paramagnetic NO to quench the paramagnetism of Y(.)(z), we are able to observe the ESEEM spectrum of deuterated acetate interacting with only the Mn cluster. A good fit of the ESEEM data indicates two (2)H dipolar hyperfine couplings of 0.097 MHz and one of 0.190 MHz. Modeling of these dipolar interactions, using our "dangler" 3 + 1 model for the S(2)-state of the Mn cluster, reveals distances consistent with direct ligation of acetate to the Mn cluster. As acetate inhibition is competitive with the essential cofactor Cl(-), this suggests that Cl(-) ligates directly to the Mn cluster. The effect of acetate binding on the structure of the Mn cluster is investigated by comparing the Mn-histidine coupling in NO/acetate-treated PSII and untreated PSII using ESEEM. We find that the addition of acetate and NO does not affect the histidine ligation to the Mn cluster. We also investigate the ability of acetate to access Y(.)(z) in Mn-depleted PSII, a PSII preparation expected to be more solvent accessible than intact PSII. We detect no coupling between Y(.)(z) and acetate. We have previously shown that small alcohols such as methanol can ligate to the Mn cluster with ease, while larger alcohols such as 2-propanol, as well as DMSO, are excluded [Force, D. A.; Randall, D. W.; Lorigan, G. A.; Clemens, K. L.; Britt, R. D. J. Am. Chem. Soc. 1998, 120, 13321-13333]. We probe the effect of acetate binding on the ability of methanol and DMSO to bind to the Mn cluster. We find that methanol is able to bind to the Mn cluster in the presence of acetate. We detect no DMSO binding in the presence of acetate. Thus, acetate binding does not increase the affinity or accessibility for DMSO binding at the Mn cluster. We also explore the possibility that the acetate binding site is also a binding site for substrate water. By comparing the ratioed three-pulse ESEEM spectra of a control, untreated PSII sample in 50% D(2)O to an NO/acetate-treated PSII sample in 50% D(2)O, we find that the binding of acetate to the oxygen evolving complex of photosystem II displaces deuterons bound very closely to the Mn cluster.

Acetates↗

PubMed2003

A model for predicting likely sites of CYP3A4-mediated metabolism on drug-like molecules.

We have developed a rapid semiquantitative model for evaluating the relative susceptibilities of different sites on drug molecules to metabolism by cytochrome P450 3A4. The model is based on the energy necessary to remove a hydrogen radical from each site, plus the surface area exposure of the hydrogen atom. The energy of hydrogen radical abstraction is conventionally measured by AM1 semiempirical molecular orbital calculations. AM1 calculations show the following order of radical stabilities for the hydrogen atom abstractions: sp2 centers > heteroatom sp3 centers > carbon sp3 centers. Since AM1 calculations are too time intensive for routine work, we developed a statistical trend vector model, which is used to estimate the AM1 abstraction energy of a hydrogen atom from its local atomic environment. We carried out AM1 and trend vector calculations on 50 CYP3A4 substrates whose major sites of metabolism are known in the literature. A plot of the lowest hydrogen radical formation energy versus its sterically accessible surface area exposure for these 50 substrates shows that only those hydrogen atoms with solvent accessible surface area exposure > or = 8.0 A(2) are susceptible to CYP3A4-mediated metabolism. This approach forms the basis for our general model, which predicts sites on drugs that are susceptible to cytochrome P450 3A4-mediated hydrogen radical abstraction followed by a hydroxylation reaction. This model, in conjunction with specific enzyme site binding requirements, can aid in identifying possible sites of metabolism catalyzed by other cytochrome P450 enzymes.

Binding Sites↗

PubMed2002

A fast empirical approach to binding free energy calculations based on protein interface information.

Three useful variables from the interfaces of 20 protein-protein complexes were investigated. These variables are the side-chain accessible number (N(b)), the number of hydrophilic pairs (N(pair)) and buried a polar solvent accessible surface areas (DeltaDeltaASA(apol)). An empirical model based on the three variables was developed to describe the free energy of protein associations. As the results show, the side-chain accessible numbers characterize the loss of side-chain conformational entropy of protein interactions and the effective empirical function presented here has great capability for estimating the binding free energy. It was found that the variables of interface information capture most of the significant features of protein-protein association. Also, we applied the model based on the variables as a rescoring function to docking simulations and found that it has the potential to distinguish the 'true' binding mode. It is clear that the simple and empirical scale developed here is an attractive target function for calculating binding free energy for various biological processes to rational protein design.

Algorithms↗

PubMed1990

Conformations of disulfide bridges in proteins.

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of less than or equal to 2.0 A, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (phi, psi) values at Cys residues, bridge dihedral angles chi ss, chi i1, chi j1, chi i2, and chi j2, the distances C alpha i-C alpha j and C beta i-C beta j between the C alpha and C beta atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson [1981) Adv. Protein Chem. 34, 167-330). Disulfide bridging across antiparallel extended strands is observed in alpha-lytic protease, crambin, and beta-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.

Cystine↗

PubMed2003

RVP-net: online prediction of real valued accessible surface area of proteins from single sequences.

SUMMARY: RVP-net is an online program for the prediction of real valued solvent accessibility. All previous methods of accessible surface area (ASA) predictions classify amino acid residues into exposure states and named them buried or exposed based on different thresholds. Real values in some cases were generated by taking the mid points of these state thresholds. This is the first method, which provides a direct prediction of ASA without making exposure categories and achieves results better than 19% mean absolute error. To facilitate batch processing of several sequences, a standalone version of this tool is also provided. AVAILABILITY: Online predictions are available at http://www.netasa.org/rvp-net/. Standalone version of the program can be obtained from the corresponding author by E-mail request.

Algorithms↗

PubMed2003

Free energy landscape of protein folding in water: explicit vs. implicit solvent.

The Generalized Born (GB) continuum solvent model is arguably the most widely used implicit solvent model in protein folding and protein structure prediction simulations; however, it still remains an open question on how well the model behaves in these large-scale simulations. The current study uses the beta-hairpin from C-terminus of protein G as an example to explore the folding free energy landscape with various GB models, and the results are compared to the explicit solvent simulations and experiments. All free energy landscapes are obtained from extensive conformation space sampling with a highly parallel replica exchange method. Because solvation model parameters are strongly coupled with force fields, five different force field/solvation model combinations are examined and compared in this study, namely the explicit solvent model: OPLSAA/SPC model, and the implicit solvent models: OPLSAA/SGB (Surface GB), AMBER94/GBSA (GB with Solvent Accessible Surface Area), AMBER96/GBSA, and AMBER99/GBSA. Surprisingly, we find that the free energy landscapes from implicit solvent models are quite different from that of the explicit solvent model. Except for AMBER96/GBSA, all other implicit solvent models find the lowest free energy state not the native state. All implicit solvent models show erroneous salt-bridge effects between charged residues, particularly in OPLSAA/SGB model, where the overly strong salt-bridge effect results in an overweighting of a non-native structure with one hydrophobic residue F52 expelled from the hydrophobic core in order to make better salt bridges. On the other hand, both AMBER94/GBSA and AMBER99/GBSA models turn the beta-hairpin in to an alpha-helix, and the alpha-helical content is much higher than the previously reported alpha-helices in an explicit solvent simulation with AMBER94 (AMBER94/TIP3P). Only AMBER96/GBSA shows a reasonable free energy landscape with the lowest free energy structure the native one despite an erroneous salt-bridge between D47 and K50. Detailed results on free energy contour maps, lowest free energy structures, distribution of native contacts, alpha-helical content during the folding process, NOE comparison with NMR, and temperature dependences are reported and discussed for all five models.

Amino Acids↗

PubMed1984

An analysis of incorrectly folded protein models. Implications for structure predictions.

Proteins with homologous amino acid sequences have similar folds and it has been assumed that an unknown three-dimensional structure can be obtained from a known homologous structure by substituting new side-chains into the polypeptide chain backbone, followed by relatively small adjustment of the model. To examine this approach of structure prediction and, more generally, to isolate the characteristics of native proteins, we constructed two incorrectly folded protein models. Sea-worm hemerythrin and the variable domain of mouse immunoglobulin K-chain, two proteins with no sequence homology, were chosen for study; the former is composed of a bundle of four alpha-helices and the latter consists of two 4-stranded beta-sheets. Using an automatic computer procedure, hemerythrin side-chains were substituted into the immunoglobulin domain and vice versa. The structures were energy-minimized with the program CHARMM and the resulting structures compared with the correctly folded forms. It was found that the incorrect side-chains can be incorporated readily into both types of structures (alpha-helices, beta-sheets) with only small structural adjustments. After constrained energy-minimization, which led to an average atomic co-ordinate shift of no more than 0.7 to 0.9 A, the incorrectly folded models arrived at potential energy values comparable to those of the correct structures. Detailed analysis of the energy results shows that the incorrect structures have less stabilizing electrostatic, van der Waals' and hydrogen-bonding interactions. The difference is particularly pronounced when the electrostatic and van der Waals' energy terms are calculated by modified equations that include an approximate representation of solvent effects. The incorrectly folded structures also have a significantly larger solvent-accessible surface and a greater fraction of non-polar side-chain atoms exposed to solvent. Examination of their interior shows that the packing of side-chains at the secondary structure interfaces, although corresponding to sterically allowed conformations, deviates from the characteristics found in normal proteins. The analysis of incorrectly folded structures has made it clear that the absence of bad non-bonded contacts, though necessary, is not sufficient to demonstrate the validity of model-built structures and that modeling of homologous structures has to be accompanied by a thorough quantitative evaluation of the results. Further, certain features that characterize native proteins are made evident by their absence in misfolded models.

Amino Acid Sequence↗

PubMed1997

Mapping the suramin-binding sites of human neutrophil elastase: investigation by fluorescence resonance energy transfer and molecular modeling.

Neutrophil elastase (NE), a mediator of inflammation, binds with high affinity numerous anionic molecules including suramin, a polysulfated naphthylurea, which inhibits it with a Ki of 0.2 microM and a 4:1 suramin:NE stoichiometry and thus constitutes a potential therapeutic agent. In an attempt to locate the suramin molecules on NE, we investigated the NE-suramin interaction using steady-state and time-resolved fluorescence spectroscopy. The time-resolved intensity decay of NE, a protein with three Trp residues, in positions 27, 141, and 237 (chymotrypsin numbering system) was best described by a three-exponential function with lifetimes ranging from 0.22 to 2.28 ns. Comparison of the accessibility of the three lifetime classes to the fluorescence quenchers acrylamide and iodide with the computed solvent accessibility of the three Trp residues in the crystal structure of NE indicates that the main, if not the sole, contribution to the 2.28 ns lifetime class is brought about by the fully buried Trp 141 residue. The addition of suramin to NE induces a sharp decrease in NE fluorescence and a corresponding increase in suramin fluorescence due to an efficient fluorescence resonance energy transfer (FRET) between the Trp residues of NE, acting as donors, and the naphthalene rings of suramin, behaving as acceptors. From the fate of the longest lifetime class in the presence of variable suramin concentrations, we deduce that two suramins are bound at less than 17 A from Trp 141, whereas the two others are located at least 29 A from Trp 141. Moreover, neither the binding of suramin to NE nor the FRET process was modified when NE was complexed with a peptide chloromethylketone inhibitor, suggesting that suramin does not directly interfere with the substrate binding site of NE. These data were used as constraints to model the NE-suramin complex.

Acrylamide↗

PubMed1997

Role of macromolecular hydration in the binding of the Escherichia coli cyclic AMP receptor to DNA.

The osmotic stress technique was used to measure the changes in macromolecular hydration that accompany binding of the Escherichia coli CAP protein to its transcription-regulatory site (C1) in the lactose promoter and that accompany the transfer of CAP from site C1 to nonspecific genomic DNA. Formation of the C1 complex is accompanied by the net release of 79 +/- 11 water molecules. If all water molecules were released from macromolecular surfaces, this result would be consistent with a net reduction of solvent-accessible surface area of 711 +/- 189 A2. This area is only slightly smaller than the solvent-inaccessible macromolecular interface in crystalline CAP-DNA complexes. The transfer of CAP from site C1 to nonspecific sites is accompanied by the net uptake of 56 +/- 10 water molecules. Taken with the water stoichiometry of sequence-specific binding, this value implies that formation of a nonspecific complex is accompanied by the net release of 2-44 water molecules. The enhanced stabilities of CAP-DNA complexes with increased osmolality (decreased water activity) may contribute to the ability of E.coli cells to tolerate dehydration and/or high external salt concentrations.

DNA, Bacterial↗

PubMed2000

Localization and quantification of hydrophobicity: the molecular free energy density (MolFESD) concept and its application to sweetness recognition.

A method for the localization, the quantification, and the analysis of hydrophobicity of a molecule or a molecular fragment is presented. It is shown that the free energy of solvation for a molecule or the transfer free energy from one solvent to another can be represented by a surface integral of a scalar quantity, the molecular free energy surface density (MolFESD), over the solvent accessible surface of that molecule. This MolFESD concept is based on a model approach where the solvent molecules are considered to be small in comparison to the solute molecule, and the solvent can be represented by a continuous medium with a given dielectric constant. The transfer energy surface density for a 1-octanol/water system is empirically determined employing a set of atomic increment contributions and distance dependent membership functions measuring the contribution of the increments to the surface value of the MolFESD. The MolFESD concept can be well used for the quantification of the purely hydrophobic contribution to the binding constants of molecule-receptor complexes. This is demonstrated with the sweeteners sucrose and sucralose and various halogen derivatives. Therein the relative sweetness, which is assumed to be proportional to the binding constant, nicely correlates to the surface integral over the positive, hydrophobic part of the MolFESD, indicating that the sweetness receptor can be characterized by a highly flexible hydrophobic pocket instead of a localized binding site.

Models, Molecular↗