Insights into protein flexibility: The relationship between normal modes and conformational change upon protein-protein docking.

Understanding protein interactions has broad implications for the mechanism of recognition, protein design, and assigning putative functions to uncharacterized proteins. Studying protein flexibility is a key component in the challenge of describing protein interactions. In this work, we characterize the observed conformational change for a set of 20 proteins that undergo large conformational change upon association (>2 A Calpha RMSD) and ask what features of the motion are successfully reproduced by the normal modes of the system. We demonstrate that normal modes can be used to identify mobile regions and, in some proteins, to reproduce the direction of conformational change. In 35% of the proteins studied, a single low-frequency normal mode was found that describes well the direction of the observed conformational change. Finally, we find that for a set of 134 proteins from a docking benchmark that the characteristic frequencies of normal modes can be used to predict reliably the extent of observed conformational change. We discuss the implications of the results for the mechanics of protein recognition.

Systems analysis of bacterial glycomes.

Bacteria produce an array of glycan-based structures including capsules, lipo-oligosaccharide and glycosylated proteins, which are invariably cell-surface-located. For pathogenic bacteria, such structures are involved in diverse roles in the life cycle of the bacterium, including adhesion, colonization, avoidance of predation and interactions with the immune system. Compared with eukaryotes, bacteria produce huge combinatorial variations of glycan structures, which, coupled to the lack of genetic data, has previously hampered studies on bacterial glycans and their role in survival and pathogenesis. The advent of genomics in tandem with rapid technological improvements in MS analysis has opened a new era in bacterial glycomics. This has resulted in a rich source of novel glycan structures and new possibilities for glycoprospecting and glycoengineering. However, assigning genetic information in predicted glycan biosynthetic pathways to the overall structural information is complex. Bioinformatic analysis is required, linked to systematic mutagenesis and functional analysis of individual genes, often from diverse biosynthetic pathways. This must then be related back to structural analysis from MS or NMR spectroscopy. To aid in this process, systems level analysis of the multiple datasets can be used to make predictions of gene function that can then be confirmed experimentally. The present paper exemplifies these advances with reference to the major gastrointestinal pathogen Campylobacter jejuni.

3D-Garden: a system for modelling protein-protein complexes based on conformational refinement of ensembles generated with the marching cubes algorithm.

MOTIVATION: Reliable structural modelling of protein-protein complexes has widespread application, from drug design to advancing our knowledge of protein interactions and function. This work addresses three important issues in protein-protein docking: implementing backbone flexibility, incorporating prior indications from experiment and bioinformatics, and providing public access via a server. 3D-Garden (Global And Restrained Docking Exploration Nexus), our benchmarked and server-ready flexible docking system, allows sophisticated programming of surface patches by the user via a facet representation of the interactors' molecular surfaces (generated with the marching cubes algorithm). Flexibility is implemented as a weighted exhaustive conformer search for each clashing pair of molecular branches in a set of 5000 models filtered from around approximately 340,000 initially. RESULTS: In a non-global assessment, carried out strictly according to the protocols for number of models considered and model quality of the Critical Assessment of Protein Interactions (CAPRI) experiment, over the widely-used Benchmark 2.0 of 84 complexes, 3D-Garden identifies a set of ten models containing an acceptable or better model in 29/45 test cases, including one with large conformational change. In 19/45 cases an acceptable or better model is ranked first or second out of 340,000 candidates. AVAILABILITY: http://www.sbg.bio.ic.ac.uk/3dgarden (server).

On the use of overlapping lattices for screening to find pairs of nearby points in two and three dimensions.

In designing an algorithm to find pairs of points that are within Euclidean distance d it is effective to use a screening procedure to reject most pairs of points that are far apart. A procedure based on multiple, overlapping lattices can efficiently identify close points and exclude distant ones.

Protein-protein docking using 3D-Dock in rounds 3, 4, and 5 of CAPRI.

In rounds 3-5 of CAPRI, the community-wide experiment on the comparative evaluation of protein-protein docking for structure prediction, we applied the 3D-Dock software package to predict the atomic structures of nine biophysical interactions. This approach starts with an initial grid-based shape complementarity search. The product of this is a large number of potential interacting conformations that are subsequently ranked by interface residue propensities and interaction energies. Refinement through detailed energetics and optimization of side-chain positions using a rotamer library is also performed. For rounds 3, 4, and 5 of the CAPRI evaluation, where possible, we clustered functional residues on the surfaces of the monomers as an indication of binding sites, using sequence based evolutionary conservations. In certain targets this provided a very useful tool for identifying the areas of interaction. During round 5, we also applied the techniques of side-chain trimming and geometrical clustering described in the literature. Of the nine target complexes in rounds 3-5, we predicted conformations that contained at least some correct contact residues for seven of these systems. For two of the targets, we submitted predictions that were considered as medium-quality. These were a nidogen-laminin complex for target 8 (T08) and a serine-threonine phosphatase bound to a targeting subunit (T14). For a further three target systems, we produced models that were rated as acceptable predictions.

Gene function hypotheses for the Campylobacter jejuni glycome generated by a logic-based approach.

Increasingly, experimental data on biological systems are obtained from several sources and computational approaches are required to integrate this information and derive models for the function of the system. Here, we demonstrate the power of a logic-based machine learning approach to propose hypotheses for gene function integrating information from two diverse experimental approaches. Specifically, we use inductive logic programming that automatically proposes hypotheses explaining the empirical data with respect to logically encoded background knowledge. We study the capsular polysaccharide biosynthetic pathway of the major human gastrointestinal pathogen Campylobacter jejuni. We consider several key steps in the formation of capsular polysaccharide consisting of 15 genes of which 8 have assigned function, and we explore the extent to which functions can be hypothesised for the remaining 7. Two sources of experimental data provide the information for learning-the results of knockout experiments on the genes involved in capsule formation and the absence/presence of capsule genes in a multitude of strains of different serotypes. The machine learning uses the pathway structure as background knowledge. We propose assignments of specific genes to five previously unassigned reaction steps. For four of these steps, there was an unambiguous optimal assignment of gene to reaction, and to the fifth, there were three candidate genes. Several of these assignments were consistent with additional experimental results. We therefore show that the logic-based methodology provides a robust strategy to integrate results from different experimental approaches and propose hypotheses for the behaviour of a biological system.