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 703 records · Page 39Linked to original sources

Ab initio prediction of the solution structures and populations of a cyclic pentapeptide in DMSO based on an implicit solvation model.

Using a recently developed statistical mechanics methodology, the solution structures and populations of the cyclic pentapeptide cyclo(D-Pro(1)-Ala(2)-Ala(3)-Ala(4)-Ala(5)) in DMSO are obtained ab initio, i.e., without using experimental restraints. An important ingredient of this methodology is a novel optimization of implicit solvation parameters, which in our previous publication [Baysal, C.; Meirovitch, H. J Am Chem Soc 1998, 120, 800-812] has been applied to a cyclic hexapeptide in DMSO. The molecule has been described by the simplified energy function E(tot) = E(GRO) + summation operator(k) sigma(k)A(k), where E(GRO) is the GROMOS force-field energy, sigma(k) and A(k) are the atomic solvation parameter (ASP) and the solvent accessible surface area of atom k. This methodology, which relies on an extensive conformational search, Monte Carlo simulations, and free energy calculations, is applied here with E(tot) based on the ASPs derived in our previous work, and for comparison also with E(GRO) alone. For both models, entropy effects are found to be significant. For E(tot), the theoretical values of proton-proton distances and (3)J coupling constants agree very well with the NMR results [Mierke, D. F.; Kurz, M.; Kessler, H. J Am Chem Soc 1994, 116, 1042-1049], while the results for E(GRO) are significantly worse. This suggests that our ASPs might be transferrable to other cyclic peptides in DMSO as well, making our methodology a reliable tool for an ab initio structure prediction; obviously, if necessary, parts of this methodology can also be incorporated in a best-fit analysis where experimental restraints are used.

Dimethyl Sulfoxide↗

SSALN: an alignment algorithm using structure-dependent substitution matrices and gap penalties learned from structurally aligned protein pairs.

In template-based modeling of protein structures, the generation of the alignment between the target and the template is a critical step that significantly affects the accuracy of the final model. This paper proposes an alignment algorithm SSALN that learns substitution matrices and position-specific gap penalties from a database of structurally aligned protein pairs. In addition to the amino acid sequence information, secondary structure and solvent accessibility information of a position are used to derive substitution scores and position-specific gap penalties. In a test set of CASP5 targets, SSALN outperforms sequence alignment methods such as a Smith-Waterman algorithm with BLOSUM50 and PSI_BLAST. SSALN also generates better alignments than PSI_BLAST in the CASP6 test set. LOOPP server prediction based on an SSALN alignment is ranked the best for target T0280_1 in CASP6. SSALN is also compared with several threading methods and sequence alignment methods on the ProSup benchmark. SSALN has the highest alignment accuracy among the methods compared. On the Fischer's benchmark, SSALN performs better than CLUSTALW and GenTHREADER, and generates more alignments with accuracy >50%, >60% or >70% than FUGUE, but fewer alignments with accuracy >80% than FUGUE. All the supplemental materials can be found at http://www.cs.cornell.edu/ approximately jianq/research.htm.

Algorithms↗

Fourier transform ion cyclotron resonance mass spectrometric detection of small Ca(2+)-induced conformational changes in the regulatory domain of human cardiac troponin C.

Troponin C (TnC), a calcium-binding protein of the thin filament of muscle, plays a regulatory role in skeletal and cardiac muscle contraction. NMR reveals a small conformational change in the cardiac regulatory N-terminal domain of TnC (cNTnC) on binding of Ca2+ such that the total exposed hydrophobic surface area increases very slightly from 3090 +/- 86 A2 for apo-cNTnC to 3108 +/- 71 A2 for Ca(2+)-cNTnC. Here, we show that measurement of solvent accessibility for backbone amide protons by means of solution-phase hydrogen/deuterium (H/D) exchange followed by pepsin digestion, high-performance liquid chromatography, and electrospray ionization high-field (9.4 T) Fourier transform Ion cyclotron resonance mass spectrometry is sufficiently sensitive to detect such small ligand binding-induced conformational changes of that protein. The extent of deuterium incorporation increases significantly on binding of Ca2+ for each of four proteolytic segments derived from pepsin digestion of the apo- and Ca(2+)-saturated forms of cNTnC. The present results demonstrate that H/D exchange monitored by mass spectrometry can be sufficiently sensitive to detect and identify even very small conformational changes in proteins, and should therefore be especially informative for proteins too large (or too insoluble or otherwise intractable) for NMR analysis.

Amino Acid Sequence↗

Combining evolutionary and structural information for local protein structure prediction.

We study the effects of various factors in representing and combining evolutionary and structural information for local protein structural prediction based on fragment selection. We prepare databases of fragments from a set of non-redundant protein domains. For each fragment, evolutionary information is derived from homologous sequences and represented as estimated effective counts and frequencies of amino acids (evolutionary frequencies) at each position. Position-specific amino acid preferences called structural frequencies are derived from statistical analysis of discrete local structural environments in database structures. Our method for local structure prediction is based on ranking and selecting database fragments that are most similar to a target fragment. Using secondary structure type as a local structural property, we test our method in a number of settings. The major findings are: (1) the COMPASS-type scoring function for fragment similarity comparison gives better prediction accuracy than three other tested scoring functions for profile-profile comparison. We show that the COMPASS-type scoring function can be derived both in the probabilistic framework and in the framework of statistical potentials. (2) Using the evolutionary frequencies of database fragments gives better prediction accuracy than using structural frequencies. (3) Finer definition of local environments, such as including more side-chain solvent accessibility classes and considering the backbone conformations of neighboring residues, gives increasingly better prediction accuracy using structural frequencies. (4) Combining evolutionary and structural frequencies of database fragments, either in a linear fashion or using a pseudocount mixture formula, results in improvement of prediction accuracy. Combination at the log-odds score level is not as effective as combination at the frequency level. This suggests that there might be better ways of combining sequence and structural information than the commonly used linear combination of log-odds scores. Our method of fragment selection and frequency combination gives reasonable results of secondary structure prediction tested on 56 CASP5 targets (average SOV score 0.77), suggesting that it is a valid method for local protein structure prediction. Mixture of predicted structural frequencies and evolutionary frequencies improve the quality of local profile-to-profile alignment by COMPASS.

Algorithms↗

Analysis of transmembrane segment 7 of the dipeptide transporter hPepT1 by cysteine-scanning mutagenesis.

To investigate the involvement of transmembrane segment 7 (TMS7) of hPepT1 in forming the putative central aqueous channel through which the substrate traverses, we individually mutated each of the 21 amino acids in TMS7 to a cysteine and analyzed the mutated transporters using the scanning cysteine accessibility method. Y287C- and M292C-hPepT1 did not express at the plasma membrane. Out of the remaining 19 transporters, three (F293C-, L296C-, and F297C-hPepT1) showed negligible glycyl-sarcosine (gly-sar) uptake activity and may play an important role in defining the overall hPepT1 structure. K278C-hPepT1 showed approximately 40% activity and the 15 other transporters exhibited more than 50% gly-sar uptake when compared with wild type (WT)-hPepT1. Gly-sar uptake for the 16 active transporters containing cysteine mutations was then measured in the presence of 2.5 mM 2-aminoethyl methanethiosulfonate hydrobromide (MTSEA) or 1 mM [2-(trimethylammonium) ethyl] methanethiosulfonate bromide (MTSET). Gly-sar uptake was significantly inhibited for each of the 16 single cysteine mutants in the presence of 2.5 mM MTSEA. In contrast, significant inhibition of uptake was only observed for K278C-, M279C-, V280C-, T281C-, M284C-, L286C-, P291C-, and D298C-hPepT1 in the presence of 1 mM MTSET. MTSET modification of R282C-hPepT1 resulted in a significant increase in gly-sar uptake. To investigate this further, we mutated WT-hPepT1 to R282A-, R282E-, and R282K-hPepT1. R282E-hPepT1 showed a 43% reduction in uptake activity, whereas R282A- and R282K-hPepT1 had activities comparable with WT-hPepT1, suggesting a role for the Arg-282 positive charge in substrate translocation. Most of the amino acids that were MTSET-sensitive upon cysteine mutation, including R282C, are located toward the intracellular end of TMS7. Hence, our results suggest that TMS7 of hPepT1 is relatively solvent-accessible along most of its length but that the intracellular end of the transmembrane domain is particularly so. From a structure-function perspective, we speculate that the extracellular end of TMS7 may shift following substrate binding, providing the basis for channel opening and substrate translocation.

Amino Acid Substitution↗

Estimating the number of bound waters in Gd(III) complexes revisited. Improved methods for the prediction of q-values.

Two literature computational methods for the prediction of the number of inner-sphere aqua ligands, q, have been applied to a test set of seven Gd(aminocarboxylate) complexes. The first method is based on the hypothesis that q should be proportional to the solvent-accessible surface area of the ligand-complexed Gd ion (Castonguay, L. A., Treasurywala, A. M., Caulfield, T. J., Jaeger, E. P., and Kellar, K. E. (1999) Prediction of q-Values and Conformations of Gadolinium Chelates for Magnetic Resonance Imaging. Bioconjugate Chem. 10, 958). The second method is based on the hypothesis that q-values can be deduced by examining series of steric energy versus ionic radii plots as a function of coordination number (Reichert, D. E., Hancock, R. D., and Welch, M. J. (1996) Molecular Mechanics Investigation of Gadolinium(III) Complexes. Inorg. Chem. 35, 7013). This study identifies deficiencies in these methods and, with respect to the first method, describes some apparent errors. Although neither method was reliable at predicting q, two alternate approaches based on either molecular mechanics strain thresholds or exposed Gd surface area thresholds are shown to predict observed q-values for all Gd aminocarboxylate complexes in the test set.

Computer Simulation↗

Ion-pair partition of quarternary ammonium drugs: the influence of counter ions of different lipophilicity, size, and flexibility.

PURPOSE: The ion-pair partition of quaternary ammonium (QA) pharmacons with organic counter ions of different lipophilicity, size, shape and flexibility was studied to elucidate relationships between ion-pair formation and chemical structure. METHODS: The apparent partition coefficient (P') of 4 QAs was measured in octanol/pH 7.4 phosphate buffer system by the shake-flask method as a function of molar excess of ten counter ions (Y), namely: mesylate (MES), acetate (AC), pyruvate (PYRU), nicotinate (NIC), hydrogenfumarate (HFUM), hydrogenmaleate (HMAL), p-toluenesulfonate (PTS), caproate (CPR), deoxycholate (DOC) and prostaglandin E1 anion (PGE1). RESULTS: Based on 118 of highly precise logP' values (SD< 0.05), the intrinsic lipophilicity (without external counter ions) and the ion-pair partition of QAs (with different counter ions) were characterized. Linear correlation was found between the logP' of ion-pairs and the size of the counter ions described by the solvent accessible surface area (SASA). The lipophilicity increasing effect of the counter ions were quantified and the following order was established: DOC approximate to PGE1 >> CPR approximate to PTS >> NIC approximate to HMAL >> PYRU approximate to AC approximate to MES approximate to HFUM. Analyzing the lipophilicity/molar ratio (QA:Y) profile, the differences in the ion-pair formation were shown and attributed to the differences in the flexibility/rigidity and size both of QA and Y. CONCLUSIONS: Since the largest (in average, 300 X) lipophilicity enhancement was found by the influence of DOC and PGE1 and considerable (on average 40 X) increase was observed by CPR and PTS, it was concluded that bile acids and prostaglandin anions may play a significant role in the ion-pair transport of quaternary ammonium drugs and caproic acid and p-toluenesulfonic acid may be useful salt forming agents to improve the pharmacokinetics of hydrophilic drugs.

Absorption↗

Engineered antibody-drug conjugates with defined sites and stoichiometries of drug attachment.

The chimeric anti-CD30 IgG1, cAC10, conjugated to eight equivalents of monomethyl auristatin E (MMAE) was previously shown to have potent antitumor activity against CD30-expressing tumors xenografts in mice. Moreover, the therapeutic index was increased by lowering the stoichiometry from 8 drugs/antibody down to 2 or 4. Limitations of such 'partially-loaded' conjugates are low yield (10-30%) as they are purified from mixtures with variable stoichiometry (0-8 drugs/antibody), and heterogeneity as the 2 or 4 drugs are distributed over eight possible cysteine conjugation sites. Here, the solvent-accessible cysteines that form the interchain disulfide bonds in cAC10 were replaced with serine, to reduce the eight potential conjugation sites down to 4 or 2. These Cys-->Ser antibody variants were conjugated to MMAE in near quantitative yield (89-96%) with defined stoichiometries (2 or 4 drugs/antibody) and sites of drug attachment. The engineered antibody-drug conjugates have comparable antigen-binding affinities and in vitro cytotoxic activities with corresponding purified parental antibody-drug conjugates. Additionally, the engineered and parental antibody-drug conjugates have similar in vivo properties including antitumor activity, pharmacokinetics and maximum tolerated dose. Our strategy for generating antibody-drug conjugates with defined sites and stoichiometries of drug loading is potentially broadly applicable to other antibodies as it involves engineering of constant domains.

Animals↗

Evaluation of a fast implicit solvent model for molecular dynamics simulations.

A solvation term based on the solvent accessible surface area (SASA) is combined with the CHARMM polar hydrogen force field for the efficient simulation of peptides and small proteins in aqueous solution. Only two atomic solvation parameters are used: one is negative for favoring the direct solvation of polar groups and the other positive for taking into account the hydrophobic effect on apolar groups. To approximate the water screening effects on the intrasolute electrostatic interactions, a distance-dependent dielectric function is used and ionic side chains are neutralized. The use of an analytical approximation of the SASA renders the model extremely efficient (i.e., only about 50% slower than in vacuo simulations). The limitations and range of applicability of the SASA model are assessed by simulations of proteins and structured peptides. For the latter, the present study and results reported elsewhere show that with the SASA model it is possible to sample a significant amount of folding/unfolding transitions, which permit the study of the thermodynamics and kinetics of folding at an atomic level of detail.

Computer Simulation↗

Fluorescence and 13C NMR determination of side-chain and backbone dynamics of synthetic melittin and melittin analogues in isotropic solvents.

The dynamics in isotopic solvents of selectively 13C labeled synthetic melittin and three analogues have been investigated by using NMR and fluorescence techniques both separately and in combination. In conjunction with the "model-free" approach to interpretation of NMR relaxation data [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4570], the availability of steady-state fluorescence anisotropy and lifetime data augment T1, T2, and NOE data to provide quantitative information about fluorophore dynamics in these peptides. A method is presented for using combined fluorescence and NMR data to obtain technique- and model-independent values for parameters describing local motion of 13C-labeled fluorophores in peptides and proteins. The dynamics of melittin and melittin analogues are found to be consistent with structural characteristics inferred from CD, fluorescence, and NMR spectral information presented in the preceding paper (Weaver et al., 1989). In particular, the mobility of the random coil peptide monomers is shown to be quite similar, while side-chain as well as peptide backbone motion in the aggregated or oligomeric species differs markedly among the analogues. For melittin itself, experimentally determined overall rotational correlation times for the monomer and tetramer agree very well with values predicted on the basis of solvent-accessible protein surface area. The local dynamics of selectively 13C-labeled Trp-19 and Gly-12 residues of melittin are also found to be consistent with peptide structure. In random coil melittin monomer, a specific model for the motion indicates that the Trp side chain moves through an approximate angle of +/- 71 degrees about the beta-gamma bond with a correlation time of 159 +/- 24 ps. In melittin tetramer, the indole moiety is spatially more confined with a flip angle of +/- 37 degrees, yet demonstrates an increased rate of motion with a correlation time of 56 +/- 8 ps. The constrained mobility of the Trp-19 side chain is consistent with motional constraints inferred from the X-ray structure of melittin tetramer. These results show that protein side-chain motion, even of moieties as large as indole, can occur on the picosecond time scale and that these motions are reasonably similar to those inferred from molecular dynamics simulations.

Bee Venoms↗

Mapping the G-actin binding surface of cofilin using synchrotron protein footprinting.

Cofilin is an actin regulatory protein that binds to both monomeric and filamentous actin, and has filament severing activity. Although crystal structures for the monomeric forms of both G-actin and cofilin have been described, the structure of the binary cofilin-G-actin complex is not available. Synchrotron protein footprinting is used to identify specific side chain residues on the cofilin surface that are buried in the formation of the cofilin-G-actin binary complex. Exposure to synchrotron X-rays results in stable oxidative modifications of aromatic, aliphatic, and sulfur-containing side chains, with the rate of modification for a particular residue being dependent on its intrinsic reactivity and solvent accessibility. The rates of modification were monitored for a number of peptides generated by digestion of oxidized cofilin, both in isolation and in its binary complex with G-actin. After binding to G-actin takes place, a significant decrease in modification rates, indicating protection of side chain groups, is seen for cofilin peptides corresponding to residues 4-20, 10-17, 83-96, 91-105, and 106-117. A number of other peptides show no change in reactivity, and are presumed to represent regions distal to the binding site. Tandem mass spectrometry demonstrates that residues Leu 13, Pro 94, Met 99, and Leu 108 and 112 directly participate in the binding interface. These results are generally consistent with, and complementary to, the results of previous site-directed mutagenesis studies and extend our understanding of the G-actin binding surface of cofilin.

Actin Depolymerizing Factors↗

Disentangling the perturbational effects of amino acid substitutions in the DNA-binding domain of p53.

The spectrum of somatic cancer-associated missense mutations in the human TP53 gene was studied in order to assess the potential structural and functional importance of various intra-molecular properties associated with these substitutions. Relating the observed frequency of particular amino acid substitutions in the p53 DNA-binding domain to their expected frequency, as calculated from DNA sequence-dependent mutation rates, yielded estimates of their relative clinical observation likelihood (RCOL). Several biophysical properties were found to display significant covariation with RCOL values. Thus RCOL values were observed to decrease with increasing solvent accessibility of the substituted residue and with increasing distance from the p53 DNA-binding and Zn2+ -binding sites. The number of adverse steric interactions introduced by an amino acid replacement was found to be positively correlated with its RCOL value, irrespective of the magnitude of the interactions. A gain in hydrogen bond number was found to be only half as likely to come to clinical attention as mutations involving either a reduction or no change in hydrogen bond number. When the difference in potential energy between the wild-type and mutant DNA-binding domains was considered, RCOL values exhibited a minimum around changes of zero. Finally, classification of mutated residues in terms of their protein/solvent environment yielded, for somatic p53 mutations, RCOL values that resembled those previously determined for inherited mutations of human factor IX causing haemophilia B, suggesting that similar mechanisms may be responsible for the mutation-related perturbation of biological function in different protein folds.

Amino Acid Substitution↗

Populations of hydrophobic amino acids within protein globular domains: identification of conserved "topohydrophobic" positions.

The 3D structural comparison of families of divergent homologous domains revealed two main populations of hydrophobic amino acids, one with a low and the other with a significantly higher mean solvent accessibility, allowing two regions of the core of protein globular domains to be distinguished. The side chains of hydrophobic amino acids in topologically conserved positions (positions in the structural alignment where only hydrophobic amino acids are found), which we call topohydrophobic positions, are considerably less dispersed than those of the other amino acids (hydrophobic or not). Mean distances between gravity centers of amino acids in topohydrophobic positions are significantly shorter than those for non-topohydrophobic positions and show that the corresponding amino acids are almost all in direct contact in the inner core of globular domains. This study also showed that the small number of topohydrophobic positions is a characteristic of the structural differences between proteins of a family. This criterion is independent of the sequence identity between the sequences and of the root-mean-square distance between their corresponding structures. Using sensitive sequence alignment processes it will be possible, for many protein families, to identify topohydrophobic positions from sequences only.

Amino Acid Sequence↗

Comparison of protein surfaces using a genetic algorithm.

A genetic algorithm (GA) is described which is used to compare the solvent-accessible surfaces of two proteins or fragments of proteins, represented by a dot surface calculated using the Connolly algorithm. The GA is used to move one surface relative to the other to locate the most similar surface region between the two. The matching process is enhanced by the use of the surface normals and shape terms provided by the Connolly program and also by a simple hydrogen-bonding descriptor and an additional shape descriptor. The algorithm has been tested in applications ranging from the comparison of small surface patches to the comparison of whole protein surfaces, and it has performed correctly in all cases. Examples of the matches are given and a quantitative analysis of the quality of the matches is performed. A number of possible future enhancements to the program are described which would allow the GA to be used for more complex surface comparisons.

Alcohol Dehydrogenase↗

Statistically significant dependence of the Xaa-Pro peptide bond conformation on secondary structure and amino acid sequence.

BACKGROUND: A reliable prediction of the Xaa-Pro peptide bond conformation would be a useful tool for many protein structure calculation methods. We have analyzed the Protein Data Bank and show that the combined use of sequential and structural information has a predictive value for the assessment of the cis versus trans peptide bond conformation of Xaa-Pro within proteins. For the analysis of the data sets different statistical methods such as the calculation of the Chou-Fasman parameters and occurrence matrices were used. Furthermore we analyzed the relationship between the relative solvent accessibility and the relative occurrence of prolines in the cis and in the trans conformation. RESULTS: One of the main results of the statistical investigations is the ranking of the secondary structure and sequence information with respect to the prediction of the Xaa-Pro peptide bond conformation. We observed a significant impact of secondary structure information on the occurrence of the Xaa-Pro peptide bond conformation, while the sequence information of amino acids neighboring proline is of little predictive value for the conformation of this bond. CONCLUSION: In this work, we present an extensive analysis of the occurrence of the cis and trans proline conformation in proteins. Based on the data set, we derived patterns and rules for a possible prediction of the proline conformation. Upon adoption of the Chou-Fasman parameters, we are able to derive statistically relevant correlations between the secondary structure of amino acid fragments and the Xaa-Pro peptide bond conformation.

Amino Acid Sequence↗

Molecular diversity sample generation on the basis of quantum-mechanical computations and principal component analysis.

The present study introduces a new strategy of selection of a maximum diversity sample of n compounds from N available in a molecular database. This strategy can be useful in pharmacological screening, combinatorial chemistry or parallel synthesis planning. It consists of first describing the compounds by means of parameters derived from quantum mechanical computations (water solvation deltaG, benzene solvation deltaG, octanol solvation deltaG, dipolar moment), as well as standard molecular parameters such as solvent-accessible surface area and molecular weight. Solvation parameters are used because of the importance of this phenomenon in the pharmacological behaviour. Redundant information in the description of the compounds is eliminated by using principal components (PC) instead of the original descriptors. Based on the similarity between the N compounds in the PC space, they are classified into n groups by k-means cluster analysis. The compounds that are nearest to the centroid of each cluster constituted the maximum diversity sample. When practical difficulties exist for the use of one of the proposed compounds, another also close to the cluster centroid can substitute for it. This strategy has been tested in the selection of a sample of 50 amines from the 923 available in the Aldrich catalogue. The results have been contrasted with those obtained from an optimal, distance-based experimental design, resulting in an 86% of agreement between both approaches. An R(2)-like diversity coefficient has been used to assess the quality of the proposed solutions.

Amines↗

Sequence homology and structural analysis of the clostridial neurotoxins.

The clostridial neurotoxins (CNTs), comprised of tetanus neurotoxin (TeNT) and the seven serotypes of botulinum neurotoxin (BoNT A-G), specifically bind to neuronal cells and disrupt neurotransmitter release by cleaving proteins involved in synaptic vesicle membrane fusion. In this study, multiple CNT sequences were analyzed within the context of the 1277 residue BoNT/A crystal structure to gain insight into the events of binding, pore formation, translocation, and catalysis that are required for toxicity. A comparison of the TeNT-binding domain structure to that of BoNT/A reveals striking differences in their surface properties. Further, the solvent accessibility of a key tryptophan in the C terminus of the BoNT/A-binding domain refines the location of the ganglioside-binding site. Data collected from a single frozen crystal of BoNT/A are included in this study, revealing slight differences in the binding domain orientation as well as density for a previously unobserved translocation domain loop. This loop and the conservation of charged residues with structural proximity to putative pore-forming sequences lend insight into the CNT mechanism of pore formation and translocation. The sequence analysis of the catalytic domain revealed an area near the active-site likely to account for specificity differences between the CNTs. It revealed also a tertiary structure, highly conserved in primary sequence, which seems critical to catalysis but is 30 A from the active-site zinc ion. This observation, along with an analysis of the 54 residue "belt" from the translocation domain are discussed with respect to the mechanism of catalysis.

Amino Acid Sequence↗

Atom depth as a descriptor of the protein interior.

Atom depth, defined as the distance (dpx, A) of a nonhydrogen atom from its closest solvent-accessible protein neighbor, provides a simple but precise description of the protein interior. Mean residue depths can be easily computed and are very sensitive to structural features. From the analysis of the average and maximum atom depths of a set of 136 protein structures, we derive a limit of approximately 200 residues for protein and protein domain size. The average and maximum atom depths in a protein are related to its size but not to the fold type. From the same set of structures, we calculated the mean residue depths for the 20 amino acid types, and show that they correlate well with hydrophobicity scales. We show that dpx values can be used to partition atoms in discrete layers according to their depth and to identify atoms that, although buried, are potential targets for posttranslational modifications like phosphorylation. Finally, we find a correlation between highly conserved residues in structural neighbors of the same fold type, and their mean residue depth in the reference structure.

Algorithms↗