PubMed Health⌕ Search

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 505 records · Page 28Linked to original sources

Achieving stability and conformational specificity in designed proteins via binary patterning.

We have developed a method to determine the optimal binary pattern (arrangement of hydrophobic and polar amino acids) of a target protein fold prior to amino acid sequence selection in protein design studies. A solvent accessible surface is generated for a target fold using its backbone coordinates and "generic" side-chains, which are constructs whose size and shape are similar to an average amino acid. Each position is classified as hydrophobic or polar according to the solvent exposure of its generic side-chain. The method was tested by analyzing a set of proteins in the Protein Data Bank and by experimentally constructing and analyzing a set of engrailed homeodomain variants whose binary patterns were systematically varied. Selection of the optimal binary pattern results in a designed protein that is monomeric, well-folded, and hyperthermophilic. Homeodomain variants with fewer hydrophobic residues are destabilized, while additional hydrophobic residues induce aggregation. Binary patterning, in conjunction with a force field that models folded state energies, appears sufficient to satisfy two basic goals of protein design: stability and conformational specificity.

Animals↗

Overcoming solvent saturation-transfer artifacts in protein NMR at neutral pH. Application of pulsed field gradients in measurements of 1H-15N Overhauser effects.

Artifacts due to solvent saturation-transfer effects result in incorrect measurements of 1H-15N heteronuclear NOE (HNOE). These artifacts are commonly observed in aqueous protein solutions at neutral pH. We describe the application of pulsed field gradients (PFGs) together with long recycle delays in overcoming errors in HNOE measurements which arise from H2O solvent preirradiation and solvent saturation transfer. Even in the absence of explicit solvent irradiation, the HNOE pulse sequence itself results in a nonequilibrium spin-state distribution of solvent nuclei which can then be transferred by chemical exchange into amide-proton sites. This effect can be avoided by using PFGs for suppression of solvent H2O together with a recycle delay sufficiently long for the magnetization of water to relax back to equilibrium values during the preparation period. These effects were studied in 15N-enriched human type alpha transforming growth factor at pH 7.1. Comparisons of PFG-HNOE experiments with and without selective H2O irradiation and with different recycle times provide estimates of the effects of solvent irradiation on HNOE measurements, which are different for different amide nitrogen-15 nuclei. The amplitudes of these artificial HNOE enhancements are roughly correlated with solvent accessibilities of amide sites in the three-dimensional structure of hTGF alpha.

Artifacts↗

Using metal-catalyzed oxidation reactions and mass spectrometry to identify amino acid residues within 10 A of the metal in Cu-binding proteins.

Metal-catalyzed oxidation (MCO) reactions and mass spectrometry (MS) can be used together to determine the amino acids bound to Cu in a metalloprotein. Selective oxidation of only amino acids bound to Cu during the MCO/MS approach relies on proper choice of the types and concentrations of the reducing and oxidizing agents. We show here that if these MCO reagent concentrations are "detuned" or varied slightly from optimal conditions, nonmetal-bound amino acids close to Cu can also be oxidized in a controlled manner. Using Cu/Zn superoxide dismutase as a model system, we demonstrate that amino acids up to 10 A from Cu can be modified as long as they are readily accessible to diffusing reactive oxygen species. UV/VIS spectroscopy and measurements of oxidation kinetics provide evidence that the protein's structural integrity around Cu is maintained during the detuned MCO reactions. Because this method can identify amino acids around Cu that have low solvent accessibility, it should complement other radical-based protein surface-mapping techniques.

Amino Acids↗

Investigation of aromatic-backbone amide interactions in the model peptide acetyl-Phe-Gly-Gly-N-methyl amide using molecular dynamics simulations and protein database search.

Weakly polar interactions between the side-chain aromatic rings and hydrogens of backbone amides (Ar-HN) are found in unique conformational regions. To characterize these conformational regions and to elucidate factors that determine the conformation of the Ar-HN interactions, four 4-ns molecular dynamics simulations were performed using four different low-energy conformations obtained from simulated annealing and one extended conformation of the model tripeptide Ac-Phe-Gly-Gly-NH-CH(3) as starting structures. The Ar(i)-HN(i+1) interactions were 4 times more frequent than were Ar(i)-HN(i+2) interactions. Half of the conformations with Ar(i)-HN(i+2) interactions also contained an Ar(i)-HN(i+1) interaction. The solvent access surface area of the Phe side chain and of the amide groups of Phe1, Gly2, and Gly3 involved in Ar-HN interactions was significantly smaller than in residues not involved in such interactions. The number of hydrogen bonds between the solvent and Phe1, Gly2, and Gly3 amide groups was also lower in conformations with Ar-HN interactions. For each trajectory, structures that contained Ar(i)-HN(i), Ar(i)-HN(i+1), and Ar(i)-HN(i+2) interactions were clustered on the basis of similarity of selected torsion angles. Attraction energies between the aromatic ring and the backbone amide in representative conformations of the clusters ranged from -1.98 to -9.24 kJ mol(-1) when an Ar-HN interaction was present. The most representative conformations from the largest clusters matched well with the conformations from the Protein Data Bank of Phe-Gly-Gly protein fragments containing Ar-HN interactions.

Amides↗

The effect of vicinal polar and charged groups on hydrophobic hydration.

The use of a linear relationship between free energy of hydrophobic hydration and solvent-accessible apolar surface area has been helpful in interpreting the thermodynamics of biological macromolecules. However, a recent study (Y.-K. Cheng, P. J. Rossky, Nature 1998, Vol. 392, pp. 696-699) has established a substantial enthalpic dependence on biomolecular surface topography, originating from solvent hydrogen-bonding loss in a restrictive geometry. In this study, we use molecular dynamics simulations of 2-Zn insulin in water solvent to explore the further effect of vicinal polar or charged groups on hydrophobic hydration at a biomolecular surface. In contrast to the case for solvent more isolated from such polar solute influences, the binding energies of the water that is proximal to the hydrophobic dimeric interface of insulin and vicinal to polar and charged groups are comparable to the bulk solvent value, a result of compensating interaction primarily with the solute counterions. The results suggest a special importance for such polar/charged groups in biological processes involving hydrophobic surface regions of restricted geometry and also suggest a general route for tuning the hydrophobicity of interfaces.

Computer Simulation↗

Assessing the relative importance of the biophysical properties of amino acid substitutions associated with human genetic disease.

The inclusion of a mutation in a pathology-based database such as the Human Gene Mutation Database (HGMD) is a two-stage process: first, the mutation must occur at the DNA level, then it must cause a clinically detectable disease state. The likelihood of the latter step, termed the relative clinical observation likelihood (RCOL), can be regarded as a function of the structural/functional consequences of a mutation at the protein level. Following this paradigm, we modeled in silico all amino acid replacements that could potentially have arisen from an inherited single base pair substitution in five human genes encoding arylsulphatase A (ARSA), antithrombin III (SERPINC1), protein C (PROC), phenylalanine hydroxylase (PAH), and transthyretin (TTR). These proteins were chosen on the basis of 1) the availability of a crystallographic structure, and 2) a sufficiently large number of amino acid replacements being logged in HGMD. A total of 9,795 possible mutant structures were modeled and 20 different biophysical parameters assessed. Together with the HGMD-derived spectra of clinically detected mutations, these data allowed maximum likelihood estimation of RCOL profiles for the 20 parameters studied. Nine parameters (including energy difference between wild-type and mutant structures, accessibility of the mutated residue, and distance from the binding/active site) exhibited statistically significant variability in their RCOL profiles, indicating that mutation-associated changes affected protein function. As yet, however, a biological meaning could only be attributed to the RCOL profiles of solvent accessibility and, for three proteins, local energy change, disturbed geometry, and distance from the active center. The limited ability of the biophysical properties of mutations to explain clinical consequences is probably due to our current lack of understanding as to which amino acid residues are critical for protein folding. However, since the proteins examined here were unrelated, and our findings consistent, it may nevertheless prove possible to extrapolate to other proteins whose dysfunction underlies inherited disease.

Amino Acid Substitution↗

Purification and cation binding properties of the recombinant human S100 calcium-binding protein A3, an EF-hand motif protein with high affinity for zinc.

The calcium-binding protein S100A3 is an unusual member of the S100 family, characterized by its very high content of Cys. In order to study the biochemical, cation-binding, and conformational properties, we produced and purified the recombinant human protein in Escherichia coli. The recombinant protein forms noncovalent homodimers, tetramers, and polymers in vitro with a subunit molecular weight of 11,712. The Zn(2+)-binding parameters of S100A3 were studied by equilibrium gel filtration and yielded a stoichiometry of four Zn2+ per monomer with a [Zn2+]0.5 of 11 microM and a Hill coefficient of 1.4 at physiological ionic strength. The affinity for Ca2+ is too low to be determined by direct methods (KCa > 10 mM). Ca(2+)- and Zn(2+)-binding can be followed by optical methods: the Trp-45 fluorescence is high in the metal-free form and addition of Zn2+ and Ca2+, but not of Mg2+, leads to a 4-fold quenching. Ca2+ and Zn2+ promote also quite similar conformational changes in the Tyr and Trp environment as monitored by difference spectrophotometry. Fluorescence titrations with Zn2+ confirmed that there is one set of high affinity binding sites with a [Zn2+]0.5 of 8 microM and a Hill coefficient of 1.3. Binding of Zn2+ to a second set of low affinity sites induces protein precipitation. Fluorescence titrations with Ca2+ confirmed the very low affinity of S100A3 for this ion with a [Ca2+]0.5 of 30 mM and slight negative cooperativity. Mg2+ has no effect on this binding curve. Of the 10 Cys residues in S100A3, 5 only are free thiols, and accessible to 5,5'-dithiobis(2-nitro-benzoic acid); they display a high reactivity in the metal-free and Ca2+ form, but a 20-fold lowered reactivity in the Zn2+ form of S100A3. Ca(2+)-binding promotes the formation of a solvent-accessible hydrophobic surface as monitored by the 60-fold fluorescence increase of 2-p-toluidinylnaphthalene-6-sulfonate, whereas Zn2+ has no noticeable influence. Our data indicate that Ca2+ and Zn2+ do not bind to the same sites and that under physiological conditions S100A3 is a Zn(2+)-binding rather than a Ca(2+)-binding protein; nevertheless, very specific conformational changes are introduced by either Ca2+ or Zn2+. Since no Zn(2+)-binding motif of known structure was identified in the primary sequence of S100A3, the results are suggestive for a novel Zn(2+)-binding motif.

Base Sequence↗

Calcium-dependent stabilization of the central sequence between Met(76) and Ser(81) in vertebrate calmodulin.

Spin-label electron paramagnetic resonance (EPR) provides optimal resolution of dynamic and conformational heterogeneity on the nanosecond time-scale and was used to assess the structure of the sequence between Met(76) and Ser(81) in vertebrate calmodulin (CaM). Previous fluorescence resonance energy transfer and anisotropy measurements indicate that the opposing domains of CaM are structurally coupled and the interconnecting central sequence adopts conformationally distinct structures in the apo-form and following calcium activation. In contrast, NMR data suggest that the opposing domains of CaM undergo independent rotational dynamics and that the sequence between Met(76) and Ser(81) in the central sequence functions as a flexible linker that connects two structurally independent domains. However, these latter measurements also resolve weak internuclear interactions that suggest the formation of transient helical structures that are stable on the nanosecond time-scale within the sequence between Met(76) and Asp(80) in apo-CaM (H. Kuboniwa, N. Tjandra, S. Grzekiek, H. Ren, C. B. Klee, and A. Bax, 1995, Nat. Struct. Biol. 2:768-776). This reported conformational heterogeneity was resolved using site-directed mutagenesis and spin-label EPR, which detects two component spectra for 1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl)-methanethiosulfonate spin labels (MTSSL) bound to CaM mutants T79C and S81C that include a motionally restricted component. In comparison to MTSSL bound within stable helical regions, the fractional contribution of the immobilized component at these positions is enhanced upon the addition of small amounts of the helicogenic solvent trifluoroethanol (TFE). These results suggest that the immobilized component reflects the formation of stable secondary structures. Similar spectral changes are observed upon calcium activation, suggesting a calcium-dependent stabilization of the secondary structure. No corresponding changes are observed in either the solvent accessibility to molecular oxygen or the maximal hyperfine splitting. In contrast, more complex spectral changes in the line-shape and maximal hyperfine splitting are observed for spin labels bound to sites that undergo tertiary contact interactions. These results suggest that spin labels at solvent-exposed positions within the central sequence are primarily sensitive to backbone fluctuations and that either TFE or calcium binding stabilizes the secondary structure of the sequence between Met(76) and Ser(81) and modulates the structural coupling between the opposing domains of CaM.

Amino Acid Motifs↗

Effects of anion proximity in peptide primary sequence on the rate and mechanism of leucine oxidation.

Hydroxyl radical surface mapping is a useful tool for investigating protein structure and folding. The rate of protein side-chain oxidation by the hydroxyl radical is known to be affected primarily by the chemical reactivity of the side chain and the accessibility of the reactive site to the radical. Efforts have been made to determine the inherent rate of stable product formation of each amino acid side chain, so that the rate of oxidation of an amino acid can be used to accurately estimate the average solvent accessibility of the amino acid side chain in the folded protein. However, the effects of nearby primary sequence on peptide oxidation have not been studied. Here, we examine the amounts of various oxidation products of a small peptide consisting of one leucine and one aspartic acid separated by zero to five glycine residues, as well as with modification of the N- and C-terminus. We find that the relative amounts of certain oxidation products can be heavily influenced by the primary structure of the surrounding peptide. The formation of many products, including hydroxylation, is inhibited by proximity to negative charges, while the formation of other products showed more complicated responses to changing primary sequence.

Amino Acid Sequence↗

Side-chain torsional potentials: effect of dipeptide, protein, and solvent environment.

Side-chain torsional potentials in the bovine pancreatic trypsin inhibitor are calculated from empirical energy functions by use of the known X-ray structure of the protein and the rigid-geometry mapping technique. The potentials are analyzed to determine the roles and relative importance of contributions from the dipeptide backbone, the protein, and the crystalline environment of solvent and other protein molecules. The structural characteristics of the side chains determine two major patterns of energy surfaces, E(X1,X2): a gamma-branched pattern and a pattern for longer, straight side chains (Arg, Lys, Glu, and Met). Most of the dipeptide potential curves and surfaces have a local minimum corresponding to the side-chain torsional angles in the X-ray structure. Addition of the protein forces sharpens and/or selects from these minima, providing very good agreement with the experimental conformation for most side chains at the surface or in the core of the protein. Inclusion of the crystalline environment produces still better results, especially for the side chains extending away from the protein. The results are discussed in terms of the details of the interactions due to the surrounding, calculated solvent-accessibility figures and the temperature factors derived from the crystallographic refinement of the pancreatic trypsin inhibitor.

Animals↗

Reduction of furin-nicked Pseudomonas exotoxin A: an unfolding story.

Upon entering mammalian cells, Pseudomonas exotoxin A (PE) is proteolytically processed by furin to produce an N-terminal fragment of 28 kDa and a C-terminal fragment of 37 kDa. Cleavage is followed by the reduction of a key disulfide bond (cysteines 265-287). This combination of proteolysis and reduction releases the 37 kDa C-terminal fragment, which then translocates to the cytosol where it ADP-ribosylates elongation factor 2 and inhibits protein synthesis. To investigate toxin reduction, furin-nicked PE or a hypercleavable mutant, PEW281A, was subjected to various treatments and then analyzed for fragment production. Reduction was evident only when unfolding conditions and a reducing agent were applied. Thermal unfolding of PE, as evidenced by changes in alpha-helical content and increased sensitivity to trypsin, rendered nicked toxin susceptible to protein disulfide isomerase- (PDI-) mediated reduction. When subcellular fractions from toxin-sensitive cells were incubated with nicked PE, toxin unfolding and reducing activities were present in the membrane fraction but not the soluble fraction. These data indicate that PE reduction is a two-step process: unfolding that allows access to the Cys265-287 disulfide bond, followed by reduction of the sulfur-sulfur bond by PDI or a PDI-like enzyme. With regard to cellular processing, we propose that the toxin's three-dimensional structure retains a "closed" conformation that restricts solvent access to the Cys265-287 disulfide bond until after a cell-mediated unfolding event.

ADP Ribose Transferases↗

CUPSAT: prediction of protein stability upon point mutations.

CUPSAT (Cologne University Protein Stability Analysis Tool) is a web tool to analyse and predict protein stability changes upon point mutations (single amino acid mutations). This program uses structural environment specific atom potentials and torsion angle potentials to predict DeltaDeltaG, the difference in free energy of unfolding between wild-type and mutant proteins. It requires the protein structure in Protein Data Bank format and the location of the residue to be mutated. The output consists information about mutation site, its structural features (solvent accessibility, secondary structure and torsion angles), and comprehensive information about changes in protein stability for 19 possible substitutions of a specific amino acid mutation. Additionally, it also analyses the ability of the mutated amino acids to adapt the observed torsion angles. Results were tested on 1538 mutations from thermal denaturation and 1603 mutations from chemical denaturation experiments. Several validation tests (split-sample, jack-knife and k-fold) were carried out to ensure the reliability, accuracy and transferability of the prediction method that gives >80% prediction accuracy for most of these validation tests. Thus, the program serves as a valuable tool for the analysis of protein design and stability. The tool is accessible from the link http://cupsat.uni-koeln.de.

Databases, Protein↗

The accessibility of iron at the active site of recombinant human phenylalanine hydroxylase to water as studied by 1H NMR paramagnetic relaxation. Effect of L-Phe and comparison with the rat enzyme.

The high-spin (S = 5/2) Fe(III) ion at the active site of recombinant human phenylalanine hydroxylase (PAH) has a paramagnetic effect on the longitudinal relaxation rate of water protons. This effect is proportional to the concentration of enzyme, with a paramagnetic molar-relaxivity value at 400 MHz and 25 degrees C of 1. 3 (+/- 0.03) x 10(3) s-1 M-1. The value of the Arrhenius activation energy (Ea) for the relaxation rate was -14.4 +/- 1.1 kJ/mol for the resting enzyme, indicating a fast exchange of water protons in the paramagnetic environment. The frequency dependence of the relaxation rate also supported this hypothesis. Thus, the recombinant human PAH appears to have a more solvent-accessible catalytic iron than the rat enzyme, in which the water coordinated to the metal is slowly exchanging with the solvent. These findings may be related to the level of basal activity before activation for these enzymes, which is higher for human than for rat PAH. In the presence of saturating (5 mM) concentrations of the substrate L-Phe, the paramagnetic molar relaxivity for human PAH decreased to 0.72 (+/- 0.05) x 10(3) s-1 M-1 with no significant change in the Ea. Effective correlation times (tauC) of 1.8 (+/- 0.3) x 10(-10) and 1.25 (+/- 0.2) x 10(-10) s-1 were calculated for the enzyme and the enzyme-substrate complex, respectively, and most likely represent the electron spin relaxation rate (tauS) for Fe(III) in each case. Together with the paramagnetic molar-relaxivity values, the tauC values were used to estimate Fe(III)-water distances. It seems that at least one of the three water molecules coordinated to the iron in the resting rat and human enzymes is displaced from coordination on the binding of L-Phe at the active site.

Animals↗

Mutations in alpha-helical solvent-exposed sites of eglin c have long-range effects: evidence from molecular dynamics simulations.

Eglin c is a small protease inhibitor whose structural and thermodynamic properties have been well studied. Previous thermodynamic measurements on mutants at solvent-accessible positions in the protein's helix have shown the unexpected result that the data could be best fit by the inclusion of residue- and position-specific parameters to the model. To explore the origins of this surprising result, long molecular dynamics simulations in explicit solvent have been performed. These simulations indicate specific long-range interactions between the solvent-exposed residues in the eglin c alpha-helix and binding loop, an unexpected observation for such a small protein. The residues involved in the interaction are on opposite sides of the protein, about 25 A apart. Simulations of alanine substitutions at the solvent-exposed helix positions, arginine 22, glutamic acid 23, threonine 26, and leucine 27, show both small and large perturbations of eglin c dynamics. Two mutations exhibit large impacts on the long-range helix-loop interactions. Previous stability measurements (Yi et al., Biochemistry 2003;42:7594-7603) had indicated that an alanine substitution at position 27 was less stabilizing than at other solvent-exposed positions in the helix. The L27A mutation effects observed in these simulations suggest that the position-dependent loss of stability measured in wet bench experiments is derived from changes in dynamics that involve long-range interactions; thus, these simulations support the hypothesis that solvent-exposed positions in helices are not always equivalent.

Computational Biology↗

A molecular dynamics analysis of protein structural elements.

The relation between protein secondary structure and internal motions was examined by using molecular dynamics to calculate positional fluctuations of individual helix, beta-sheet, and loop structural elements in free and substrate-bound hen egg-white lysozyme. The time development of the fluctuations revealed a general correspondence between structure and dynamics; the fluctuations of the helices and beta-sheets converged within the 101 psec period of the simulation and were lower than average in magnitude, while the fluctuations of the loop regions were not converged and were mostly larger than average in magnitude. Notable exceptions to this pattern occurred in the substrate-bound simulation. A loop region (residues 101-107) of the active site cleft had significantly reduced motion due to interactions with the substrate. Moreover, part of a loop and a 3(10) helix (residues of 67-88) not in contact with the substrate showed a marked increase in fluctuations. That these differences in dynamics of free and substrate-bound lysozyme did not result simply from sampling errors was established by an analysis of the variations in the fluctuations of the two halves of the 101 psec simulation of free lysozyme. Concerted transitions of four to five mainchain phi and psi angles between dihedral wells were shown to be responsible for large coordinate shifts in the loops. These transitions displaced six or fewer residues and took place either abruptly, in 1 psec or less, or with a diffusive character over 5-10 psec. Displacements of rigid secondary structures involved longer timescale motions in bound lysozyme; a 0.5 A rms change in the position of a helix occurred over the 55 psec simulation period. This helix reorientation within the protein appears to be a response to substrate binding. There was little correlation between the solvent accessible surface area and the dynamics of the different structural elements.

Binding Sites↗

Solvent-induced distortions and the curvature of alpha-helices.

The alpha-helix defined in 1951 by Pauling et al. on the basis of model building and X-ray fibre diffraction data has 3.65 residues per turn (n) achieved with planar peptides, torsion angles of phi = -48 degrees and psi = -57 degrees and hydrogen bonds which are close to linear. Although X-ray analyses of proteins have confirmed the general correctness of the model for the helix, recent high resolution (1.7-1.0 A) diffraction studies have shown that the parameters described by Pauling et al. and later by Perutz and Arnott and Wonacott are not a good description of the alpha-helices in globular proteins, where the mean values of phi, psi are usually close to -63 degrees, -42 degrees. Here we show that these values arise as a mean of two significantly different classes in amphipathic helices depending on whether the peptide carbonyl oxygen is hydrogen bonded to a solvent or polar side-chain atom. The hydrogen bonds made by the hydrophilic carbonyls to the NH groups within helices are longer and less linear than those involving hydrophobic carbonyls. We also show that these effects are associated with a significant curvature of helices in globular proteins. For example, the alpha-helix in avian pancreatic peptide (aPP) has a radius of curvature of approximately 70 A. These results are of significance in the packing of helices in fibrous and globular proteins, in the calculation of their dipole moments, solvent accessibilities and internal energies, and in the theoretical estimation of spectroscopic properties such as circular dichroism and Raman scattering.

Pancreatic Polypeptide↗

The cytoplasmic fragment of the aspartate receptor displays globally dynamic behavior.

A number of cloned soluble fragments if the bacterial chemotaxis transmembrane receptors retain partial function. Prior studies of a fragment corresponding to the cytoplasmic domain (c-fragment) of the Escherichia coli aspartate receptor have correlated the signaling state of mutant receptors with the oligomerization state of the c-fragments: equilibria of smooth-swimming mutants are shifted toward oligomeric states; tumble mutants are shifted toward monomeric states [Long, D. G., & Weis, R. M. (1992) Biochemistry 31, 9904-9911]. We have applied several experimental probes of local and global structural flexibility to two signaling states, the wild-type (monomeric) and S461L smooth mutant (predominantly dimeric) c-fragments. Featureless near-UV CD spectra are observed, which indicate that the single Trp residue is in a symmetric environment (most likely averaged by fluctuations) and suggest that the C-termini of both proteins are highly mobile. Both proteins undergo extremely rapid proteolysis and enhance ANS fluorescence, which indicates that many sites are accessible to trypsin cleavage and hydrophobic sites are accessible to ANS binding. The global nature of the flexibility is demonstrated by 1H NMR studies. Lack of chemical shift dispersion suggests that fluctuations average the environments of side chains and backbone protons. Rapid exchange of 99% of the observable amide protons suggests that these fluctuations give high solvent accessibility to nearly the entire backbone. This evidence indicates that both monomeric and dimeric c-fragments are globally flexible proteins, with properties similar to "molten-globule" states. The significance of this flexibility depends on whether it is retained in functioning receptors: the c-fragment structure may lack important tertiary contacts, protein-protein interactions, or topological constraints needed to stabilize a nondynamic native structure, or the cytoplasmic domain of the native receptor may retain flexibility which may be modulated in the mechanism of transmembrane signaling.

Anilino Naphthalenesulfonates↗

Coupling nonpolar and polar solvation free energies in implicit solvent models.

Recent studies on the solvation of atomistic and nanoscale solutes indicate that a strong coupling exists between the hydrophobic, dispersion, and electrostatic contributions to the solvation free energy, a facet not considered in current implicit solvent models. We suggest a theoretical formalism which accounts for coupling by minimizing the Gibbs free energy of the solvent with respect to a solvent volume exclusion function. The resulting differential equation is similar to the Laplace-Young equation for the geometrical description of capillary interfaces but is extended to microscopic scales by explicitly considering curvature corrections as well as dispersion and electrostatic contributions. Unlike existing implicit solvent approaches, the solvent accessible surface is an output of our model. The presented formalism is illustrated on spherically or cylindrically symmetrical systems of neutral or charged solutes on different length scales. The results are in agreement with computer simulations and, most importantly, demonstrate that our method captures the strong sensitivity of solvent expulsion and dewetting to the particular form of the solvent-solute interactions.

Computer Simulation↗