Prediction of protein-protein binding sites

A new algorithm to predict protein-protein binding sites using conservation of both protein surface structure and physical-chemical properties in structurally similar proteins is developed. Binding-site residues in proteins are known to be more conserved than the rest of the surface, and finding local surface similarities by comparing a protein to its structural neighbors can potentially reveal the location of binding sites on this protein. This approach, which has previously been used to predict binding sites for small ligands, is now extended to predict protein-protein binding sites. Examples of binding-site predictions for a set of proteins, which have previously been studied for sequence conservation in protein-protein interfaces, are given. The predicted binding sites and the actual binding sites are in good agreement. Our algorithm for finding conserved surface structures in a set of similar proteins is a useful tool for the prediction of protein-protein binding sites.



Graph theory

A new algorithm for finding a maximum clique in an undirected graph is described. An approximate coloring algorithm has been improved and used to provide bounds to the size of the maximum clique in a basic algorithm which finds a maximum clique. This basic algorithm was then extended to include dynamically varying bounds. The resulting algorithm is significantly faster than the comparable algorithm.



Virtual screening

We have screened the NCI diversity set against an enzyme involved in bacterial cell wall biosynthesis. Contrary to currently available antibiotics which target extracelular enzymes, we targeted an intracelular step of the biosynthesis. We obtained three active compounds that exhibited a very good antibacterial activity on cultured bacteria.



Biomolecular simulations

Molecular dynamics simulation has been performed to investigate the structural properties of perifosine and its synthetic spin-labeled alkylphospholipid analogues. The conformations adopted by these compounds in water and in a dipalmitoylphosphatidylcholine bilayer as a function of the presence and position of the N-oxyl- 4',4'-dimethyloxazolidine ring (doxyl group) have been investigated by all-atom molecular dynamics. No predominant conformation was observed in water, but the molecules adopt specific orientations and conformations in the lipid bilayer. As is expected, alkyl chains tend to insert into the hydrophobic core, while charged groups stay at the lipid-water interface. A doxyl group in the middle of the alkyl chain moves up to the interface region, thus preventing adoption of the extended conformation. Compounds with a doxyl group close to the polar head group adopt conformations similar to that of unlabeled perifosine within the first nanoseconds of simulation. When the doxyl group is at the end of alkyl chain, the spin-labeled molecule needs more time to reach equilibrium. These results indicate a considerable effect of the doxyl position within the alkyl chain on its localization in the lipid bilayer and can be extended further to other similar spin probes used in the electron paramagnetic resonance spectroscopy of biological membranes.