Computer-aided design of polyketides with the required properties

An approach to rational design of new polyketides with the required spectrum of biological activity has been proposed. We have developed the BioGenPharm software, which generates combinatorial libraries of polyketides based on the user-defined input parameters, performs prediction of biological activity spectra for the generated structures and selection of molecuels with the required properties. PASS algorithm has been applied for prediction of polyketide activity spectra (http://www.ibmc.msk.ru/PASS). Validation of PASS prediction ability for polyketides was performed vs. the evaluation set containing 242 natural macrolides from the Dictionary of Natural Products. The mean prediction accuracy was 75.5%. The problem of choice of the cutting points for probability of the presence of activity (Pa), which provides optimal combination of such parameters as sensitivity, specificity, concordance has been considered. Applicability of the described method has been illustrated by generation of a virtual library of the erythromycin analogues and selection of substances with low probability of the hepatotoxic effect.     

The synthesis and hepatoprotective activity of esters of the lupane group triterpenoids

Hemisuccinates, hemiphthalates, acetylsalicylates, cinnamates, andp-methoxycinnamates of lupeol, betulin, and 3-O-acetylbetulin were synthesized via interaction with corresponding acid anhydrides or acid chlorides. A number of betulin esters in position 3 and 28 were shown to exhibit a pronounced hepatoprotective effect similar to that of betulin and silibor. These experimental data were in a good agreement with the computer prediction of their biological activity. Betulin 3,28-bishemiphthalate was more effective than carsil in models of experimental hepatitis caused by carbon tetrachloride, tetracycline, and ethanol.     

Capacities of computer evaluation of hidden potential of phytochemicals of medicinal plants of the traditional Indian Ayurvedic medicine

Applicability of our computer programs PASS and PharmaExpert for prediction of biological activity spectra of rather complex and structurally diverse phytocomponents of medicinal plants, both separately and in combinations has been evaluated. For this purpose we have created the web-resource containing known information about structural formulas and biological activity of 1906 phytocomponents of 50 Ayurvedic medicinal plants used in Traditional Indian Medicine (TIM) (http://ayurveda.pharmaexpert.ru). The PASS training set was updated by addition of information about structure and biological activity of 946 natural compounds; then the training procedure and validation were performed, to estimate the quality of PASS prediction. It was shown that the differences between the average accuracy of prediction obtained in leave-5%-out cross-validation (94.467%) and in leave-one-out cross-validation (94.605%) are very small thus demonstrating high predictive ability of the program. Results of biological activity spectra prediction for all phytocomponents included in our database coincided in 83.5% of cases with known experimental data. Additional types of biological activity predicted with high probability indicate further promising directions for further studies of certain phytocomponents of some medicinal plants. The analysis of prediction results of sets of phytocomponents in each of 50 medicinal plants was made by the PharmaExpert software. Based on this analysis, we found that the combination of phytocomponents from Passiflora incarnata may exhibit nootropic, anticonvulsant, and antidepressant effects. Experiments carried out in mice models confirmed the predicted effects of P. incarnata extracts.     

Computer Assessment of the Xenobiotic Metabolites Formation’s Probability in the Human Body

Biophysics - Abstract—Xenobiotics undergo biotransformation in the human body and the resulting metabolites may greatly differ in their physical, chemical, and biological properties from the...     

Role of the electrostatic interactions in pre-orientation of subunits in the formation of protein-protein complexes

Theoretical estimation of contribution of the electrostatic interactions to pre-orientation of ribonuclease subunits in process of complex formation was carried out. The subunit was considered as a multipole consisting of partial charges of all atoms of the molecule. The object of investigation was a system of two subunits with their centers of gravity fixed at some distance in vacuum. It was proposed that each subunit independently could rotate freely around its fixed center of gravity. The relative orientation states of the subunits in such system were searched at which the system has electrostatic energy minima (equilibrium states). In first approximation the equilibrium states were found using especially designed approximate method for electrostatic interaction energy calculation, which permitted to calculate and compare the energies of the system in 24(5) (approximately 8 10(6)) states with different mutual orientation of subunits. The angular coordinates of the found equilibrium states were further specified by calculation with gradient sliding method. Angular coordinates of the equilibrium states and the shapes of energy surface cuts along each coordinate angle were calculated also for the intersubunits distances diminished down to 50 angstroms. The dispersions of the angular coordinates of equilibrium states caused by heat movement (at T=300 degrees) and their changes with shortening the distance between centers of gravity of subunits were estimated. Mutual orientation of subunits in the equilibrium states of the system under consideration was found to be similar to their mutual orientations in complex. Also it was found that relaxation time of the system, caused by electrostatic interaction of subunits, after removing the system from an equilibrium state, is much less in vacuum than the mean time between their Brownian collisions at room temperature. It follows from these results that in the case of ribonuclease in vacuum the electrostatic interactions of its subunits must be strong enough to realize the effective pre-orientation of subunits during their Brownian approach from distances of the order 100 angstroms. Preliminary consideration taking into account the effect of surrounding water molecules on the electrostatic interactions of ribonuclease subunits showed that weakening of the interaction must be much less than in the case when one uses in its calculation the macroscopic dielectric permeability value equal to 80. So the results obtained for vacuum seem to be true for water solution also. More strict theoretical analysis of this problem will be carried out in the following publication.     

Design,Synthesis and Pharmacological Evaluation of Novel Vanadium-Containing Complexes as Antidiabetic Agents

Based on the data about structure and antidiabetic activity of twenty seven vanadium and zinc coordination complexes collected from literature we developed QSAR models using the GUSAR program. These QSAR models were applied to 10 novel vanadium coordination complexes designed in silico in order to predict their hypoglycemic action. The five most promising substances with predicted potent hypoglycemic action were selected for chemical synthesis and pharmacological evaluation. The selected coordination vanadium complexes were synthesized and tested in vitro and in vivo for their hypoglycemic activities and acute rat toxicity. Estimation of acute rat toxicity of these five vanadium complexes was performed using a freely available web-resource (http://way2drug.com/GUSAR/acutoxpredict.html). It has shown that the selected compounds belong to the class of moderate toxic pharmaceutical agents, according to the scale of Hodge and Sterner. Comparison with the predicted data has demonstrated a reasonable correspondence between the experimental and predicted values of hypoglycemic activity and toxicity. Bis{tert-butyl[amino(imino)methyl]carbamato}oxovanadium (IV) and sodium(2,2′-Bipyridyl)oxo-diperoxovanadate(V) octahydrate were identified as the most potent hypoglycemic agents among the synthesized compounds.     

Mapping of the active site of alcohol dehydrogenase with low-molecular ligands

In search of an active alcohol dehydrogenase inhibitor, the structure of which may serve as the basis for a potential drug design, the active site of alcohol dehydrogenase containing NAD and Zn²⁺ ions was mapped using the method of molecular mechanics. Molecular docking was performed using a number of ligands containing characteristic functional groups: formate ion, ammonia, ammonium ion, methanol, and methylamine. Sites of preferable binding were revealed for each ligand and arranged in order of decreasing energy of binding to the enzyme. A comparison of the predicted ligand-binding sites and the experimental data on the location of water and inhibitor binding sites in the known structures of corresponding alcohol dehydrogenase complexes indicated a coincidence of the complex formation sites, which confirms the validity of the method and provides the requirements for a highly effective inhibitor (the pharmacophore model).     

Molecular mechanisms of protein-protein recognition: whether the surface placed charged residues determine the recognition process?

We studied the structure and composition of contact areas in 812 different kind dimeric protein-protein complexes from Brookhaven data base (PDB ) in order to reveal their pecularities with regard to protein-protein recognition. We have found, that the large portion of complexes (approximately 70%) have oppositely charged residues in the contact areas (interfaces) on the subunits surfaces, which form electrostatic contacts - R:E, R:D, K:E, K:D, H:E, H:D. These results are consistent with the current view that high rate complex formation may be driven by the long-range electrostatic interaction between charged AA residues of subunits surfaces. However, there are many complexes among the studied ones (approximately 30%), which have no electrostatic contacts at all in their contact area. Thus a question arises: what forces account for high complex formation rates (i.e. for the distant orienting of subunits before encounter) by forming complexes where the surface contact areas lack electrostatic contacts? We believe that the long-range orienting electrostatic interaction of subunits may account for all cases of efficient complex formation if one drops the traditional view that protein subunits interact mainly through their surfaces. We suggest that the distant orienting being due to the electrostatic interaction between the whole aggregates of partial electric charges of atoms of each complex subunits. Our preliminary model calculations (unpublished) made for ribonuclease dimer (does not have electrostatic contacts) conform this suggestion.