在本报告中，开发了两种算法，算法1和算法2.开发算法1以搜索具有短的活性区域的全长蛋白的多肽链的相互作用。开发了算法2，以确定在从N末端到C末端的方向上形成二聚体时全长蛋白之间最活跃的相互作用位点。使用蛋白质MDM2，NAP1，P53进行数值计算。对于现代蛋白质组学，蛋白质相互作用的研究和预测是非常重要的任务，因为它们确定蛋白质从细胞水平到整个生物体的功能。对于已知结构所知的蛋白质，根据其三级结构的构象的已知数据搜索分子间相互作用减少了搜索两个相互作用分子表面的几何互补性的问题，并将其触点建模，所谓的分子对接。分子对接的任务是构象搜索算法的任务，其由于蛋白质分子的扭转角度的变化而减少了所形成的生物复合物的构象空间。在大多数情况下，现代化构象搜索算法在相对较短的时间内找到一般接近实验发现的结构的构象。然而，存在对对接的成功产生重大影响的因素，这通常不会在标准算法中考虑。一种这种因素是靶蛋白的构象迁移率。 The mobility range can be different beginning with a small adjustment of the side chains and ending with scale domain movements. These movements play an important role. At first glance, the most logical solution to this problem is to take into account the mobility of the protein in a docking program. Unfortunately, modern computational tools do not allow such modeling to be performed in an acceptable time frame since a protein molecule is very large, and allowing for mobility over all degrees of freedom can lead to a so-called combinatorial explosion (an astronomical increase in the number of possible variants). Only in some programs is there a limited mobility of protein binding sites (usually at the level of a small adaptation of conformations of the side chains of the active center residues). Another approach to this problem consists in docking the same protein in several different conformations and then selecting the best solutions from each docking run. The third approach is to find a universal structure of the target protein in which docking would produce fairly good results for different classes of ligands. In this case, the number of missed (but correct) solutions decreases, but the number of incorrect options also increases significantly. It should also be noted that most programs for the theoretical docking of proteins work according to the following principle: one protein is fixed in space, and the second is rotated around it in a variety of ways. At the same time, for each rotation configuration, estimates are made for the evaluation function. The evaluation function is based on surface complementarity (the mutual correspondence of complementary structures (macromolecules, radicals), determined by their chemical properties), electrostatic interactions, van der Waals repulsion and so on. The problem with this approach is that calculations throughout the configuration space require a lot of time, rarely leading to a single solution, which in turn does not allow us to speak of the uniqueness of the target protein and ligand interaction variant. So in the work while modeling by the methods of molecular dynamics, from 200 to 10 000 possible combinations of the formation of a protein complex with a ligand were found. Such a large number of modifications, along with the lack of a criterion for selecting the most probable variants of the bound structures of biological complexes (which would allow a radical reduction in their number) makes it very difficult to interpret the theoretical results obtained for practical use, namely, the finding of catalytic centers and a qualitative assessment of the dissociation constant of interacting substances. In contrast to the above computer simulation algorithms, mathematical algorithms have been developed in this chapter that allow determining the detection of proteins active regions and detecting the stability of different regions of protein complexes (linear docking) by analyzing the potential energy matrix of pairwise electrostatic interaction between different sites of the biological complex, such as the homodimer of the histone chaperone Nap1-Nap1, the heterodimer of the p53 Mdm2 proteins, and the homodimer Mdm2 Mdm2, which are responsible for the entry of a whole protein molecule into biochemical reactions.
T V Koshlan