|
|
||||||||
Department of Biomolecular Sciences, UMIST, Manchester M60 1QD, UK
Reprint requests to: Jim Warwicker, Department of Biomolecular Sciences, UMIST, PO Box 88, Manchester M60 1QD, UK; e-mail: jim.warwicker{at}umist.ac.uk; fax: 44 (0)161 236 0409.
(RECEIVED July 3, 2002; FINAL REVISION September 13, 2002; ACCEPTED September 23, 2002)
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.0222702.
| Abstract |
|---|
|
|
|---|
Keywords: Conformational entropy; rotamers; dimerization; protein-protein interactions
| Introduction |
|---|
|
|
|---|
An understanding of the molecular details of protein recognition sites would allow for the automated modeling of protein complexes from known monomer structures. This is pursued in the computational field of protein docking, whereby relative conformational space is searched and the resulting conformers ranked according to some force field or scoring function (Smith and Sternberg 2002). Surface shape complementarity is the most accurate predictor of proteinprotein complexes at present (Ritchie and Kemp 2000). However, this is most effective with protein components extracted from a known complex, rather than component structures that have been experimentally determined outside of the complex.
It is intriguing that promising complexes can be docked with the relatively simple potentials that describe shape complementarity, and that surface plasticity is a critical factor beyond this (Brady and Sharp 1997; Kimura et al. 2001). Finding the most discriminating potential, including shape variability, is an important goal if accurate prediction of protein complexes is to be achieved. One approach to identifying elements of an effective potential for ranking docked configurations is to characterize the properties of known interfaces. This can be achieved by examining the chemical characteristics and residue propensities of the interfacial region in the context of the rest of the protein surface (Lo Conte et al. 1999; Valdar and Thornton 2001a, 2001b), or by carrying out different ranking methodologies on arbitrary surface regions (Jones and Thornton 1997a, 1997b). In the latter method key features of proteinprotein interactions are determined using "surface patches," whereby equal-sized regions of the protein surface are characterized in terms of a variety of factors, such as planarity, hydrophobicity, and residue propensity. Another useful approach is to use evolutionary data to identify conserved residues on the protein surface, as these tend to be involved in binding and/or recognition (Lichtarge et al. 1997; Hu et al. 2000; Elcock and McCammon 2001). This method is less useful for structures without an adequate number of close homologs.
The current study extends analysis of interfacial properties (Ponstingl et al. 2000; Valdar and Thornton 2001a, 2001b) with the addition of a quantity that is the ratio of the estimated change in side-chain conformational entropy (S) upon complexation and the solvent accessible surface area (A). This quantity was calculated for total surface area (S/A), and was examined both for individual residues and over surface patches (Jones and Thornton 1997a). The entropic and accessible area components are likely to be major contributors to binding/docking. While precise two-component docking will be determined by complementarity, of interest in this initial study is the intrinsic propensity of a surface, in terms of minimal side-chain entropy loss for a maximal surface area burial.
The loss of side-chain conformational entropy upon complexation has been discussed in the context of protein folding (Lee et al. 1994; Creamer 2000) and proteinprotein complexation (Doig and Sternberg 1995; Brady and Sharp 1997). These discussions have given rise to values of about 1.5R per residue equivalent to 3.74 kJ mol-1 at 300 K, upon folding. Estimates of side-chain conformational entropy were made with a modification of the mean field algorithm of Koehl and Delarue (Koehl and Delarue 1994). In scanning residues or areas at an interface, it is assumed that monomer side-chain conformational entropy and solvent accessible surface of each interfacial residue would be lost in the complex. It may be expected that more favorable binding surfaces correlate with a smaller loss of side-chain conformational entropy and a larger burial of (nonpolar) surface area. The hypothesis that the magnitude of S/A may be smaller for interfacial regions is tested for a set of 25 homodimers and 14 heterodimers (Jones and Thornton 1997a).
| Results |
|---|
|
|
|---|
S) were found to be either interfacial or in direct contact with interfacial residues. Interfacial residues were assessed as those with changes in solvent accessible area on dimerization. Residues with only one rotamer (Pro, Gly and Ala) as well as disulphide bridges (Cys) were fixed in the mean field calculations.
Some interfacial residues were found to have zero or very small
S upon complexation. These were mainly residues with relatively few allowed rotamers (e.g., Thr, Ser, Val), for which a small
S would be expected. One or two more flexible residues (e.g., Glu, Gln, His) in each interface were also found to have a relatively small
S. This lack of change upon complexation was in conjunction with a low monomeric S value, suggesting that they are in a conformation that is favorable for dimer formation.
To investigate differences between interfacial and noninterfacial surfaces, separate
S/
A values were calculated, where these quantities are given in units of R per 100 Å2. As seen in Table 1
, the side chains of 21 out of the 25 homodimers were less flexible (up to 34.9%) at the interface than elsewhere. The heterodimer data (Table 2
) shows a similar trend for the large protomer set, but the small protomer set has interfaces that are only slightly less flexible than the overall surface. Overall, surfaces in the small protomer set exhibit less flexibility (per unit area) than those in the large protomer or homodimer sets.
|
|
S/
A provides information beyond that found in the crystallographic B-factor. In terms of a flexibility measure,
S/
A will be less dependent on crystal contacts than the B-factor.
|
|
S/
A. This result suggests that although simple occurrence of residue type contributes to conformational flexibility at the homodimer interface, a large part of the decreased flexibility generally observed at homodimer interfaces (Table 1
|
S/
A were determined for each patch.
Initial calculations showed that 68% of the interfacial patches over the 25 homodimers, 66% over the nine larger protomer heterodimers, but only 8% over the 13 smaller protomer heterodimers appeared in the first 10 percentile (i.e., ranked first) for patch analysis (Table 1
, Fig. 4
). As previously demonstrated, the difference in
S/
A at the interface and the rest of the protein surface is not mainly due to any intrinsic inflexibility of the side chains at the interface, but is largely due to side-chain restriction with respect to the rest of the protein surface. Therefore, patch analysis was carried out over a range of patch sizes in an attempt to analyze the interface and, in particular, any specific residues which dominate the interface.
|
S/
A or by percentage overlap with the true interface. It should be noted that in this method not one of the randomly generated patches completely overlaps with the true interface; the range for the most overlapping patch in an individual protein is 51.9%95.2%. This is indicative of the shape difference between the automated patches and the true interface, with the automated patches being approximately circular and contiguous over the surface, which is rarely the case for the true interface. Plotting the top hit of each against patch size for a structure yields different information (Fig. 5
S/
A hit (most inflexible patch) line shows that after an initial decrease in
S/
A over the first six patches there is a steady increase in
S/
A as the patch size increases, indicating that beyond a threshold patch size (e.g., six for 3sdh)
S/
A will tend to increase with the number of residues in each patch.
|
S/
A as the patch size increases. However, the individual patches are generally more variable with distinct peaks and troughs in the data. The troughs identify regions of low
S/
A in the interfacial surface region. This allows us to determine whether the interface is uniformly inflexible or is dominated by subregions of inflexibility. For example, in the plot for 3sdh (Fig. 5
Additionally, the interface can be analyzed on a per residue basis to see if there are residues that dominate (Fig. 6
). Generally, only two or three residues per interface stand out. They tend to be either fixed residues which are relatively solvent exposed (e.g., Val or Pro), thereby contributing low side-chain flexibility, or large residues (e.g., Lys), which confer higher flexibility to the interface. However, these stand-out residues still only alter the average
S/
A for the interfacial patch by at most ±
8%.
|
Ssc). Buried nonpolar area can also be determined for the structures, and an associated free energy estimated with an effective surface tension of 0.1 kJ mol-1 per Å2 of buried nonpolar area, that is within the generally used range (Raschke et al. 2001). At 300 K these contributions to protein dimerization (
Gnp,
Gsc) are estimated in Table 3
Gsc can also make a significant (unfavorable) contribution to binding, consistent with the view that evolved modulation of
Gsc will impact on affinity.
|
| Discussion |
|---|
|
|
|---|
S/
A gives insight beyond that from crystallographic B-factor or intrinsic side-chain flexibility at the interface.
Of the four homodimers that are more flexible at the interface than elsewhere, two (3sdh and 3ssi) have very small differences of -0.03 and -0.01, respectively. Of all the homodimers, 3ssi has the lowest overall
S/
A (Table 1
), and
Ssc (dimerization) is small (Table 3
). These features may relate to the extensive exposed ß-sheet structure, which is partially used to form the dimer interface. Several water molecules are involved in a hydrogen-bonding network across the interface for 3sdh. Such interactions are not included in our analysis, but they may be indicative of a degree of functional flexibility in this dimeric haemoglobin interface (Royer 1994). Indeed, this flexibility could underlie our observation of relatively high
S/
A for the interface in a monomer context.
Structures 2wrp and 2ccy have relatively large negative differences between the interfacial and noninterfacial regions (Table 1
). Tryptophan repressor (2wrp) is an intertwined dimer, so that reference to an extracted monomer state is probably inappropriate in this case. The very large negative
S/
A difference value for the cytochrome c 2ccy (Finzel et al. 1985) is due to the highest conformational entropy at the interface for the whole dataset and a relatively low conformational entropy at the noninterface (Table 1
). In an attempt to determine whether this was an isolated example, homologous cytochrome c proteins from four sources were also examined Chromatium vinosum: 1bbh (Ren et al. 1993), Alcaligenes denitrificans: 1cgo (Dobbs et al. 1996), Alcaligenes xylosoxidans: 1e83 (Lawson et al. 2000), and Rhodocyclus gelatinosus: 1jaf (Archer et al. 1997) (Table 4
). All the structures are four-helical bundle homodimers with a monomeric C
RMSD of
1.9 Å from 2ccy using the combinatorial extension method (Shindyalov and Bourne 1998). Despite the structural similarity, these structures yield a
S/
A difference between -10.7% and 14.5% (Table 4
), indicating that the underlying helical framework does not determine the original result with 2ccy. The range is most likely due to substantial sequence variation on top of a well-conserved dimeric structural framework.
|
Ssc, the loss of conformational entropy upon complexation, will make a significant contribution to calculations of binding affinity, which is borne out by our observation that the interface is generally less flexible than the rest of the protein surface suggesting that an interface can be preconditioned through restriction of side-chain rotamers. One might expect such an effect to be highlighted for systems with particularly tight binding. Table 5
S/
A differences for three endonuclease colicin structures and their associated immunity proteins: 1emv and 7cei are both DNases and 1e44 an RNase. DNase immunity proteins have dissociation constants in the femtomolar region (Wallis et al. 1995), and have significant structural and sequence similarities (Kühlmann et al. 2000). The RNases are structurally dissimilar to the DNases, and their binding is less well characterized.
|
S/
A % difference of all the dimers here analyzed (38.1 and 43.3%, respectively). In addition, they have the lowest interfacial
S/
A, suggesting that the large % difference is due to the interface being extremely inflexible. The endonuclease partners to the immunity proteins for 1emv and 7cei do not show this level of interfacial inflexibility, reflected also in relatively low rankings in patch analysis (Table 5
In an attempt to explore further the relationship between the strength of dimer association and side-chain conformational entropy a set of known structures with experimentally determined association free energies (Horton and Lewis 1992) were studied (Table 6
). While the interfacial regions generally exhibit lower
S/
A than noninterfacial regions, there is no correlation apparent between this property and the association free energy. This result is consistent with a view of interfacial energetics in which binding energy results from a complex combination of multiple factors.
|
S/
A at the interface than elsewhere. For the heterodimers, this is mostly the larger protomer. Whereas
S/
A is about equal for interfacial regions of the larger and smaller protomers, it tends to be larger for noninterfacial regions in the larger protomer (Table 2
S/
A, thereby reducing surface conformational entropy relative to the larger protomer set (Table 2The fact that computational protein docking is so much easier with separated components from the complex, than with structures solved outside of the complex, points to a degree of induced fit mechanism for interface formation. Generally, a mean field algorithm (Koehl and Delarue 1994) for side-chain placement can play a role in refining computed dockings, but the lower resolution question that we have asked is whether there is a detectable interfacial signature in terms of side-chain flexibility per unit surface area. Our results indicate that this is the case, although the size of the signal is not constant over the systems studied, and will also be coupled to properties such as individual monomer (folded) stability and other binding energy components. Clearly, proteins that can adapt their binding sites to bind several different ligands will probably not have a low conformational entropy at the interface, and their binding energy is likely to be dominated by other components (DeLano et al. 2000).
This work lends itself to progression in at least two directions. First, addition of the
S/
A quantity to the patch analysis approach to prediction of potential interface regions (Jones and Thornton 1997b) is promising, particularly if combined with clustering techniques to improve overlap between computationally generated patches and true interfaces. Second, the scale of calculated
Ssc suggests that such a term should be included in empirical estimates of binding affinities. Indeed, the range of
Ssc in Table 3
(about 80 kJ/mole) is equivalent to a change of
1014 in association constant, although these mean field values are probably overestimates considering that side-chain/side-chain correlations and other potential terms could lead to a narrowing of the rotamer probability distribution (Koehl and Delarue 1994). The first suggested direction relates to a monomer-based method for estimating interfacial propensity, while the second (
Ssc) method is relevant for analysis of complexes (either experimental or computed).
| Materials and methods |
|---|
|
|
|---|
Definition of the interface
Surface residues were classed as those with an accessible surface area (A) in the monomer of
0.1 Å2. The interface was defined as the set of residues for which A decreased by
0.1 Å2 upon dimerization. The probe radius was 1.4 Å. All HETATM records were retained except for nonfunctional ligands (e.g., glycol) and water oxygen atoms.
Side-chain conformational entropy
A mean field method for calculating side-chain/side-chain interactions in the context of a rotamer set (Koehl and Delarue 1994) has previously been adapted to look at the packing and maximum possible solvent accessibility of ionisable side-chains in proteins (J. Warwicker, unpubl.). Koehl and Delarue (1994) calculated the effective potential (E) on rotamer k of side-chain i as:
![]() | ((1)) |
In looking at side-chain packing, the Lennard-Jones VdW interaction parameters used to describe the potential in Equation 1
(Koehl and Delarue 1994) were replaced with hard sphere collisions, such that VdW overlap is disallowed subject to an overall relaxation of VdW radii that is incremented until a packing solution is found. In this adaptation of the method, Equation 1
reduces to:
![]() | ((2)) |
The backbone-independent rotamer set of Tuffery (Tuffery et al. 1997) was used. In this adaptation of the algorithm, which was designed to survey reasonable alternative packing solutions for side-chains in known structures, the experimental rotamers were also included, thereby minimizing the required VdW relaxation. Where studying configurations for which experimental side-chain rotamers are not available, the Lennard-Jones potential formalism (Koehl and Delarue 1994) is preferable because the current algorithm would impose a uniform VdW relaxation across all interactions, leading to wide-scale clashes.
The conformational matrix was used to estimate the conformational entropy of side chains (Hill 1956; Koehl and Delarue 1994):
![]() | ((3)) |
Surface and interface properties
In the context of proteinprotein interactions, the ratio of conformational entropy to solvent accessible area (A) was calculated for each residue (i), Si/Ai, and this measure for a patch of residues (constituting either the true interface or a computed test patch) as
Si/
Ai. Such sums over residues were calculated for various surface regions, and the i subscript dropped for simplicity. Anp is used to denote the nonpolar contribution to the total solvent accessible area, A. Percentage difference between the interfacial and the noninterfacial surface regions was calculated as:
![]() | ((4)) |
Patch analysis
As previously described (Jones and Thornton 1997a), each patch was defined with a central surface-accessible residue and built using a defined number of nearest-neighbor residues (based on C
C
atom distances). Variable patch sizes were created by starting with n = 1 and continuing up to the size of the experimental interface. Solvent vectors (Jones and Thornton 1997a) were not applied to these patches. The number of patches generated per structure is a function of the patch size and the proteins surface area, but generally did not exceed 400.
Percentage overlap of the patches with the experimental interface was determined as:
![]() |
| Acknowledgments |
|---|
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| References |
|---|
|
|
|---|
Berman, H.M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T.N., Weissig, H., Shindyalov, I.N., and Bourne, P.E. 2000. The protein data bank. Nucleic Acids Res. 28: 235242.
Brady, G.P. and Sharp, K.A. 1997. Entropy in protein folding and in proteinprotein interactions. Curr. Opin. Struct. Biol. 7: 215221.[CrossRef][Medline]
Creamer, T.P. 2000. Side-chain conformational entropy in protein unfolded states. Proteins Struct.Funct. Genet. 40: 443450.
DeLano, W.L., Ultsch, M.H., de Vos, A.M., and Wells, J.A. 2000. Convergent solutions to binding at a proteinprotein interface. Science 287: 12791283.
Dobbs, A.J., Anderson, B.F., Faber, H.R., and Baker, E.N. 1996. Three-dimensional structure of cytochrome c from two Alcaligenes species and the implications for four-helix bundle structures. Acta Crystallogr. D D52: 356368.[CrossRef]
Doig, A.J. and Sternberg, M.J.E. 1995. Side-chain conformational entropy in protein-folding. Protein Sci. 4: 22472251.[Abstract]
Elcock, A.H. and McCammon, J.A. 2001. Identification of protein oligomerization states by analysis of interface conservation. Proc. Natl. Acad. Sci. 98: 29902994.
Finzel, B.C., Weber, P.C., Hardman, K.D., and Salemme, F.R. 1985. Structure of ferricytochrome c from Rhodospirillum molischianum at 1.67Å resolution. J. Mol. Biol. 186: 627643.[CrossRef][Medline]
Glaser, F., Steinberg, D.M., Vakser, I.A., and Ben-Tal, N. 2001. Residue frequencies and pairing preferences at proteinprotein interfaces. Proteins Struct. Funct. Genet. 43: 89102.[CrossRef][Medline]
Hill, T.L. 1956. Statistical mechanics. McGraw Hill, New York.
Horton, N. and Lewis, M. 1992. Calculation of the free energy of association for protein complexes. Protein Sci. 1: 169181.[Abstract]
Hu, Z., Ma, B., Wolfson, H., and Nussinov, R. 2000. Conservation of polar residues as hot spots at protein interfaces. Proteins Struct. Funct. Genet. 39: 331342.[CrossRef][Medline]
Ito, T., Chiba, T., and Yoshida, M. 2001. Exploring the protein interactome using comprehensive two-hybrid projects. Trends Biotechnol. 19: S23S27.[CrossRef][Medline]
Jones, S. and Thornton, J.M. 1996. Principles of proteinprotein interactions. Proc. Natl. Acad. Sci. 93: 1320.
. 1997a. Analysis of proteinprotein interaction sites using surface patches. J. Mol. Biol. 272: 121132.[CrossRef][Medline]
. 1997b. Prediction of proteinprotein interaction sites using patch analysis. J. Mol. Biol. 272: 133143.[CrossRef][Medline]
Kimura, S.R., Brower, R.C., Vajda, S., and Camacho, C.J. 2001. Dynamical view of the positions of key side chains in proteinprotein recognition. Biophys. J. 80: 635642.
Koehl, P. and Delarue, M. 1994. Application of a self-consistent mean field theory to predict protein side-chains conformation and estimate their conformational entropy. J. Mol. Biol. 239: 249275.[CrossRef][Medline]
Kühlmann, U.C., Pommer, A.J., Moore, G.R., James, R., and Kleanthous, C. 2000. Specificity in proteinprotein interactions: The structural basis for dual recognition in endonuclease colicinimmunity protein complexes. J. Mol. Biol. 301: 11631178.[CrossRef][Medline]
Lawson, D.M., Stevenson, C.E.M., Andrew, C.R., and Eady, R.R. 2000. Unprecedented proximal binding of nitric oxide to heme: Implications for guanylate cyclase. EMBO J. 19: 56615671.[CrossRef][Medline]
Lee, K.H., Xie, D., Freire, E., and Amzel, L.M. 1994. Estimation of changes in side-chain configurational entropy in binding and foldingGeneral methods and application to helix formation. Proteins Struct. Funct. Genet. 20: 6884.[CrossRef][Medline]
Lichtarge, O., Yamamoto, K.R., and Cohen, F.E. 1997. Identification of functional surface of the zinc binding domains of intracellular receptors. J. Mol. Biol. 274: 325337.[CrossRef][Medline]
Lo Conte, L., Chothia, C., and Janin, J. 1999. The atomic structure of proteinprotein recognition sites. J. Mol. Biol. 285: 21772198.[CrossRef][Medline]
Luker, G.D., Sharma, V., Pica, C.M., Dahlheimer, J.L., Li, W., Ochesky, J., Ryan, C.E., Piwnica-Worms, H., and Piwnica-Worms, D. 2002. Noninvasive imaging of proteinprotein interactions in living animals. Proc. Natl. Acad. Sci. 99: 69616966.
Ponstingl, H., Henrick, K., and Thornton, J. M. 2000. Discriminating between homodimeric and monomeric proteins in the crystalline state. Proteins Struct. Funct. Genet. 41: 4757.[CrossRef][Medline]
Raschke, T.M., Tsai, J., and Levitt, M. 2001. Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water. Proc. Natl. Acad. Sci. 98: 59655969.
Ray, P., Pimenta, H., Paulmurugan, R., Berger, F., Phelps, M.E., Iyer, M., and Gambhir, S.S. 2002. Noninvasive quantitative imaging of proteinprotein interactions in living subjects. Proc. Natl. Acad. Sci. 99: 31053110.
Ren, Z., Meyer, T., and McRee, D.E. 1993. Atomic structure of a cytochrome c with an unusual ligand-controlled dimer dissociation at 1.8Å resolution. J. Mol. Biol. 234: 433445.[CrossRef][Medline]
Ritchie, D.W. and Kemp, G.J.L. 2000. Protein docking using spherical polar Fourier correlations. Proteins Struct. Funct. Genet. 39: 178194.[CrossRef][Medline]
Royer, W.E. 1994. High-resolution crystallographic analysis of a co-operative dimer hemoglobin. J. Mol. Biol. 235: 657681.[CrossRef][Medline]
Shindyalov, I.N. and Bourne, P.E. 1998. Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng. 11: 739747.
Smith, G.R. and Sternberg, M.J.E. 2002. Prediction of proteinprotein interactions by docking methods. Curr. Opin. Struct. Biol. 12: 2835.[CrossRef][Medline]
Tuffery, P., Etchebest, C., and Hazout, S. 1997. Prediction of protein side chain conformations: A study on the influence of backbone accuracy on conformation stability in the rotamer space. Protein Eng. 10: 361372.
Valdar, W.S.J. and Thornton, J.M. 2001a. Conservation helps to identify biologically relevant crystal contacts. J. Mol. Biol. 313: 399416.[CrossRef][Medline]
. 2001b. Proteinprotein interfaces: Analysis of amino acid conservation in homodimers. Proteins Struct. Funct. Genet. 42: 108124.[CrossRef][Medline]
Wallis, R., Moore, G.R., James, R., and Kleanthous, C. 1995. Proteinprotein interactions in colicin E9 DNase-immunity protein complexes. Diffusion-controlled association and femtomolar binding for the cognate complex. Biochemistry 34: 1374313750.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:
![]() |
Y. C. Chen and C. Lim Common physical basis of macromolecule-binding sites in proteins Nucleic Acids Res., November 6, 2008; (2008) gkn868v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
C.-e. A. Chang, W. A. McLaughlin, R. Baron, W. Wang, and J. A. McCammon Entropic contributions and the influence of the hydrophobic environment in promiscuous protein-protein association PNAS, May 27, 2008; 105(21): 7456 - 7461. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Grueninger, N. Treiber, M. O. P. Ziegler, J. W. A. Koetter, M.-S. Schulze, and G. E. Schulz Designed Protein-Protein Association Science, January 11, 2008; 319(5860): 206 - 209. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.-X. Zhou and S. Qin Interaction-site prediction for protein complexes: a critical assessment Bioinformatics, September 1, 2007; 23(17): 2203 - 2209. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Wollert, D. W. Heinz, and W.-D. Schubert Thermodynamically reengineering the listerial invasion complex InlA/E-cadherin PNAS, August 28, 2007; 104(35): 13960 - 13965. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Liang, C. Zhang, S. Liu, and Y. Zhou Protein binding site prediction using an empirical scoring function Nucleic Acids Res., August 7, 2006; 34(13): 3698 - 3707. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Warwicker Improved pKa calculations through flexibility based sampling of a water-dominated interaction scheme Protein Sci., October 22, 2004; 13(10): 2793 - 2805. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Moutevelis and J. Warwicker Prediction of pKa and redox properties in the thioredoxin superfamily Protein Sci., October 22, 2004; 13(10): 2744 - 2752. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. R. Caffrey, S. Somaroo, J. D. Hughes, J. Mintseris, and E. S. Huang Are protein-protein interfaces more conserved in sequence than the rest of the protein surface? Protein Sci., January 1, 2004; 13(1): 190 - 202. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Ma, T. Elkayam, H. Wolfson, and R. Nussinov Protein-protein interactions: Structurally conserved residues distinguish between binding sites and exposed protein surfaces PNAS, May 13, 2003; 100(10): 5772 - 5777. [Abstract] [Full Text] [PDF] |
||||
| ||||||||