Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain–side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein’s structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein–protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at . The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation. Figure Fuzzy-oil-drop model applied to identify the ligation site in lysozyme complexed with N-acetylglucosamine (PDB ID:1LMQ) in form of hydrophobicity deficiency (ΔH) profile and three-dimensional distribution of on protein surface