Efficient restraints for protein-protein docking by comparison of observed amino acid substitution patterns with those predicted from local environment

The discovery that the functions of most eukaryotic gene products are mediated through multi-protein complexes makes the prediction of protein interactions one of the most important current challenges in structural biology. Rigid-body docking methods can generate a large number of alternative candidates, but it is difficult to discriminate the near-native interactions from the large number of false positives. Many different scoring functions have been developed for this purpose, but in most cases, experimental and biological information is still required for accurate predictions. We explore here the use of evolutionary restraints in evaluating rigid-body docking geometries. In order to identify potential interface residues we identify functional residues based on the comparison of observed amino acid substitutions with those predicted from local environment. The interface residues identified by this method are correctly located in 85% of the cases. These predicted interface residues are used to define distance restraints that help to score rigid-body docking solutions. We have developed the pyDockRST software, which uses the percentage of satisfied distance restraints, together with the electrostatics and desolvation binding energy, to identify correct docking orientations. This methodology dramatically improves the docking results when compared to the use of energy criteria alone, and is able to find the correct orientation within the top 20 docking solutions in 80% of the cases.     

Kinship and sociality in coastal river otters: are they related?

Previous studies of coastal river otters (Lontra canadensis)in Prince William Sound, Alaska, USA, documented atypical socialorganization for mammals. Social groups were composed largelyof males, but some males remained solitary year-round and mostfemales were asocial. Because, in carnivores, groups are usuallycomposed of highly related individuals but group living alsoprovides advantages unrelated to kinship, we concurrently evaluatedthe role of relatedness and ecological benefits in socialityamong coastal river otters. By using DNA microsatellite analysisand radiotelemetry, we were able to reject the hypothesis thatsocial groups of otters were kin based. In addition, we foundno indication of kin avoidance, as would be expected from lowdispersal and high local competition. Sociality conferred noreproductive benefits or costs to otters; number of offspringand number of relatives in the population did not differ betweensocial and solitary animals. Solitary males were not older orlarger than were social males, and there was no relation betweenmale size and number of offspring, indicating that sexual selectiondid not mask a potential relation between sociality and reproductivesuccess. Among coastal river otters in this region, socialitycould be explained by the benefits obtained from cooperativeforaging on high-quality schooling pelagic fishes. Such benefitsdid not require association with kin, resulting in no selectionpressure for kin-based groups. The prediction that the degreeof sociality in the population will fluctuate relative to theabundance of schooling pelagic fishes merits further investigation.     

Prediction by graph theoretic measures of structural effects in proteins arising from non-synonymous single nucleotide polymorphisms

Recent analyses of human genome sequences have given rise to impressive advances in identifying non-synonymous single nucleotide polymorphisms (nsSNPs). By contrast, the annotation of nsSNPs and their links to diseases are progressing at a much slower pace. Many of the current approaches to analysing disease-associated nsSNPs use primarily sequence and evolutionary information, while structural information is relatively less exploited. In order to explore the potential of such information, we developed a structure-based approach, Bongo (Bonds ON Graph), to predict structural effects of nsSNPs. Bongo considers protein structures as residue-residue interaction networks and applies graph theoretical measures to identify the residues that are critical for maintaining structural stability by assessing the consequences on the interaction network of single point mutations. Our results show that Bongo is able to identify mutations that cause both local and global structural effects, with a remarkably low false positive rate. Application of the Bongo method to the prediction of 506 disease-associated nsSNPs resulted in a performance (positive predictive value, PPV, 78.5%) similar to that of PolyPhen (PPV, 77.2%) and PANTHER (PPV, 72.2%). As the Bongo method is solely structure-based, our results indicate that the structural changes resulting from nsSNPs are closely associated to their pathological consequences.     

Sequence analyses and comparative modeling of fly and worm fibroblast growth factor receptors indicate that the determinants for FGF and heparin binding are retained in evolution.

The presence of a large number of fibroblast growth factors (FGFs) and multiple splice forms of their receptors (FGFRs) in higher vertebrates makes the three-dimensional (3D) analysis of FGF interactions with their receptors a formidable task. The situation differs in Caenorhabditis elegans (worm) and Drosophila melanogaster (fruit fly), where only one or two FGF and FGFR sequences have been identified. Structural studies of the FGF-FGFR complexes in such primitive organisms should reveal the basic features of the ligand-receptor interactions as they first emerged through evolution. We have analysed the sequences of worm and fly FGFs and FGFRs and used the recently determined crystal structure of the human FGF1-FGFR2-heparin ternary complex [Pellegrini, L., Burke, D.F., von Delft, F., Mulloy, B. and Blundell, T.L. (2000) Nature 407, 1029-34] to construct 3D models of the homologous complexes. In spite of a low sequence similarity with their human counterparts, key structural features required for ligand-receptor and protein-heparin binding in humans are conserved in the fly and worm FGF-FGFR-heparin complexes. Analyses of the models show that tertiary interactions that are not conserved in sequence are maintained through novel interactions or complementary mutations in the fly and worm sequences. The overall charge distributions observed in the human FGF-FGFR-heparin complex are retained in the fly and worm models. The arginine residue at position 253 in the linker region between the Ig-like domains D2 and D3 in the wild type fly and worm sequences is particularly striking, as the Pro253Arg mutation in humans is responsible for Apert syndrome. This change may enhance the affinity of receptors for their FGF molecules as observed in Apert mutants.     

A six-stranded double-psi beta barrel is shared by several protein superfamilies.

BACKGROUND: Six-stranded beta barrels with a pseudo-twofold axis are found in several proteins. One group comprises a Greek-key structure with all strands antiparallel; an example is the N-terminal domain of ferredoxin reductase. Others involve parallel strands forming two psi structures (the double-psi beta barrel). A recently discovered example of the latter class is aspartate-alpha-decarboxylase (ADC) from Escherichia coli, a pyruvoyl-dependent tetrameric enzyme involved in the synthesis of pantothenate. RESULTS: Visual inspection and automated database searches identified the six-stranded double-psi beta barrel in ADC, Rhodobacter sphaeroides dimethylsulfoxide (DMSO) reductase, E. coli formate dehydrogenase H (FDHH), the plant defense protein barwin, Humicola insolens endoglucanase V (EGV) and, with a circular permutation, in the aspartic proteinases. Structure-based sequence alignments revealed several interactions including hydrophobic contacts or sidechain-mainchain hydrogen bonds that position the middle beta strand under a psi loop, which may significantly contribute to stabilizing the fold. The identification of key interactions allowed the filtering of weak sequence similarities to some of these proteins, which had been detected by sequence database searches. This led to the prediction of the double-psi beta-barrel domain in several families of proteins in eukaryotes and archaea. CONCLUSIONS: The structure comparison and clustering study of double-psi beta barrels suggests that there could be a common homodimeric ancestor to ADC, FDHH and DMSO reductase, and also to barwin and EGV. There are other protein families with unknown structure that are likely to adopt the same fold. In the known structures, the protein active sites cluster around the psi loop, indicating that its rigidity, protrusion and free mainchain functional groups may be well suited to providing a framework for catalysis.     

Modeling alpha-helical transmembrane domains: the calculation and use of substitution tables for lipid-facing residues.

Amino acid substitution tables are calculated for residues in membrane proteins where the side chain is accessible to the lipid. The analysis is based upon the knowledge of the three-dimensional structures of two homologous bacterial photosynthetic reaction centers and alignments of their sequences with the sequences of related proteins. The patterns of residue substitutions show that the lipid-accessible residues are less conserved and have distinctly different substitution patterns from the inaccessible residues in water-soluble proteins. The observed substitutions obtained from sequence alignments of transmembrane regions (identified from, e.g., hydrophobicity analysis) can be compared with the patterns derived from the substitution tables to predict the accessibility of residues to the lipid. A Fourier transform method, similar to that used for the calculation of a hydrophobic moment, is used to detect periodicity in the predicted accessibility that is compatible with the presence of an alpha-helix. If the putative transmembrane region is identified as helical, then the buried and exposed faces can be discriminated. The presence of charged residues on the lipid-exposed face can help to identify the regions that are in contact with the polar environment on the borders of the bilayer, and the construction of a meaningful three-dimensional model is then possible. This method is tested on an alignment of bacteriorhodopsin and two related sequences for which there are structural data at near atomic resolution.     

The 3-D structure of HIV-1 proteinase and the design of antiviral agents for the treatment of AIDS

A proteinase is essential for replication of HIV. Cloning and chemical synthesis have provided a sufficient supply of HIV-1 proteinase for the determination of its three-dimensional structure. Analogies between the structures of HIV-1 proteinase and the mammalian enzyme renin, which is involved in the control of blood pressure, have given important clues concerning the design of specific inhibitors that have antiviral activity.