Computer-aided estimation of the hERG-mediated cardiotoxicity risk of potential drug components.

The hERG potassium channel is one of the most important anti-targets determining cardiotoxicity of potential drugs. Using fragmental descriptors and artificial neural networks, the predictive models of the relationship between the structure of organic compounds and their activity with respect to hERG were built, and the structural factors affecting it were analyzed. By their predictive ability and applicability domain, these models (N = 1000, Q 2 = 0.77, RMSE cv = 0.45 for affinity and N = 2886, Q 2 = 0.60, RMSE cv = 0.55 for channel inhibition) are superior to the previously published models and can be used to minimize the risk of cardiotoxicity during drug development.

Molecular design of proneurogenic and neuroprotective compounds-allosteric NMDA receptor modulators.

N-Methyl-D-aspartic acid (NMDA) receptor is a promising target for treatment of neurodegenerative diseases and other brain disorders as well as for designing proneurogenic compounds able to stimulate neurogenesis in adult brain. We analyzed the structure of the binding site of negative allosteric modulators in the amino-terminal domain of the NMDA receptor and identified possible modes of their binding as well as performed molecular design of new modulators that significantly differ from the known ones in structure and binding mode. In addition, we formed a focused library of chemical compounds with potential neuroprotective and proneurogenic properties, desirable set of pharmacokinetic properties, and low toxicity, which can be the basis for development of new-generation drugs.

Prediction of blood-brain barrier permeability of organic compounds.

Using fragmental descriptors and artificial neural networks, a predictive model of the relationship between the structure of organic compounds and their blood-brain barrier permeability was constructed and the structural factors affecting the readiness of this penetration were analyzed. This model (N = 529, Q 2 = 0.82, RMSE cv = 0.32) surpasses the previously published models in terms of the prediction accuracy and the applicability domain and can be used for the optimization of the pharmacokinetic parameters during drug development.

Novel bivalent positive allosteric modulators of AMPA receptor.

A positive allosteric modulator of AMPA receptors has been designed using computer-aided molecular modeling techniques. It possessed a record high experimentally confirmed potency in the picomolar concentration range and belongs to a new type of bivalent AMPA receptor ligands containing bicyclo[3.3.1]nonane scaffold. The suggested structure could serve as a basis for further optimization and development of drugs for the treatment of neurodegenerative diseases, cognition enhancement, and improvement of memory.

On mechanism of allosteric modulation of NMDA receptor via amino-terminal domains.

A possible mechanism of action of the allosteric modulators of NMDA (N-methyl-d-aspartate) receptors is proposed that involves the stabilization of the twisted closed-clamshell configuration of the amino-terminal domains of GluN1 and GluN2B subunits by negative modulators while positive modulators stabilize a roughly parallel tight arrangement of these domains. These respective motions may play an important role in the transition between the open-channel and closed-channel states of the receptor. In addition, some features of the negative modulator binding site found by means of the molecular dynamics study and pocket analysis can be used in the rational design of the allosteric NMDA receptor modulators.

Combined QSAR studies of inhibitor properties of O-phosphorylated oximes toward serine esterases involved in neurotoxicity, drug metabolism and Alzheimer's disease.

Oxime reactivation of serine esterases (EOHs) inhibited by organophosphorus (OP) compounds can produce O-phosphorylated oximes (POXs). Such oxime derivatives are of interest, because some of them can have greater anti-EOH potencies than the OP inhibitors from which they were derived. Accordingly, inhibitor properties of 58 POXs against four EOHs, along with pair-wise selectivities between them, have been analysed using different QSAR approaches. EOHs (with their abbreviations and consequences of inhibition in parentheses) comprised acetylcholinesterase (AChE: acute neurotoxicity; cognition enhancement), butyrylcholinesterase (BChE: inhibition of drug metabolism or stoichiometric scavenging of EOH inhibitors; cognition enhancement), carboxylesterase (CaE: inhibition of drug metabolism or stoichiometric scavenging of EOH inhibitors), and neuropathy target esterase (NTE: delayed neurotoxicity). QSAR techniques encompassed linear regression and backpropagation neural networks in conjunction with fragmental descriptors containing labelled atoms, Molecular Field Topology Analysis (MFTA), Comparative Molecular Similarity Index Analysis (CoMSIA), and molecular modelling. All methods provided mostly consistent and complementary information, and they revealed structural features controlling the 'esterase profiles', i.e. patterns of anti-EOH activities and selectivities of the compounds of interest. In addition, MFTA models were used to design a library of compounds having a cognition-enhancement esterase profile suitable for potential application to the treatment of Alzheimer's disease.

The "hockey sticks" effect revisited: the conformational and electronic properties of 3,7-dithia-1,5-diazabicyclo[3.3.1]nonane from the QTAIM perspective.

The conformational effects in bicyclo[3.3.1]nonanes, while thoroughly studied, have not yet received the full theoretical explanation. R. F. W. Bader's quantum theory of atoms in molecules presents unique opportunities for studying the stereoelectronic interactions (SEI) and weak intramolecular bonding leading to these effects. Here, we report the study of 3,7-dithia-1,5-diazabicyclo[3.3.1]nonane by means of the topological analysis of the calculated (MP2(full)/6-311++G**) and experimental (X-ray derived) charge density to reveal the origins of the so-called "hockey sticks" effect observed in similar compounds. A new explanation of the relative stability of bicyclo[3.3.1]nonane conformers based on the analysis of the QTAIM atomic energies is given. The H···H and S···S interactions in bicyclo[3.3.1]nonane and its dithia derivatives are shown to be significant factors contributing to the differences in the relative stability of the conformers.

Antioxidative activity of ferrocenes bearing 2,6-di-tert-butylphenol moieties.

The antioxidative activity of ferrocenes bearing either 2,6-di-tert-butylphenol or phenyl groups has been compared using DPPH (1,1-diphenyl-2-picrylhydrazyl) test and in the study of the in vitro impact on lipid peroxidation in rat brain homogenate and on some characteristics of rat liver mitochondria. The results of DPPH test at 20 degrees C show that the activity depends strongly upon the presence of phenolic group but is improved by the influence of ferrocenyl fragment. The activity of N-(3,5-di-tert-butyl-4-hydroxyphenyl)iminomethylferrocene (1), for instance, was 88.4%, which was higher than the activity of a known antioxidant 2,6-di-tert-butyl-4-methylphenol (BHT) (48.5%), whereas the activity of N-phenyl-iminomethylferrocene 2 was almost negligible -2.9%. The data obtained demonstrate that the compounds with 2,6-di-tert-butylphenol moiety are significantly more active than the corresponding phenyl analogues in the in vitro study of lipid peroxidation in rat brain homogenate. Ferrocene 1 performs a promising behavior as an antioxidant and inhibits the calcium-dependent swelling of mitochondria. These results allow us to propose the potential cytoprotective (neuroprotective) effect of ditopic compounds containing antioxidant 2,6-di-tert-butylphenol group and redox active ferrocene fragment.

New Test System for Serine/Threonine Protein Kinase Inhibitors Screening: E. coli APHVIII/Pk25 design.

An efficient test system for serine/threonine protein kinase inhibitors screening has been developed based on theE. coliprotein system APHVIII/Pk25. Phosphorylation of aminoglycoside phosphotransferase VIII (APHVIII) by protein kinases enhances resistance of the bacterial cell to aminoglycoside antibiotics, e.g. kanamycin. Addition of protein kinase inhibitors prevents phosphorylation and increases cell sensitivity to kanamycin. We have obtained modifications of APHVIII in which phosphorylatable Ser146 was encompassed into the canonical autophosphorylation sequence ofStreptomyces coelicolorPk25 protein kinase. Mutant and wild-typeaphVIII were cloned intoE. coliwith the catalytic domain ofpk25. As a result of the expression of these genes, accumulation of corresponding proteins was clearly observed. Extracted from bacterial lysates, Pk25 demonstrated its ability to autophosphorylate. It was shown that variants ofE. colicontaining bothaphVIIIand рк25were more resistant to kanamycin than those carrying onlyaphVIII. Protein kinase inhibitors of the indolylmaleimide class actively inhibited Pk25 and reduced cell resistance to kanamycin. Modeling of APHVIII and Pk25 3D structures showed that pSer146 is an analog of phosphoserine in the ribose pocket of protein kinase A. Pk25 conformation was similar to that of РknB ofMycobacterium tuberculosis. Potential indolylmaleimide inhibitors were docked into the ATP-binding pocket of Pk25. The designed test system can be used for the primary selection of ATP-competitive small molecule protein kinase inhibitors.

Kinetics and mechanism of inhibition of serine esterases by fluorinated aminophosphonates.

This paper reviews previously published data and presents new results to address the hypothesis that fluorinated aminophosphonates (FAPs), (RO)(2)P(O)C(CF(3))(2)NHS(O)(2)C(6)H(5), R=alkyl, inhibit serine esterases by scission of the P-C bond. Kinetics studies demonstrated that FAPs are progressive irreversible inhibitors of acetylcholinesterase (AChE, EC 3.1.1.7.), butyrylcholinesterase (BChE, EC 3.1.1.8.), carboxylesterase (CaE, EC 3.1.1.1.), and neuropathy target esterase (NTE, EC 3.1.1.5.), consistent with P-C bond breakage. Chemical reactivity experiments showed that diMe-FAP and diEt-FAP react with water to yield the corresponding dialkylphosphates and (CF(3))(2)CHNHS(O)(2)C(6)H(5), indicating lability of the P-C bond. X-ray crystallography of diEt-FAP revealed an elongated (and therefore weaker) P-C bond (1.8797 (13)A) compared to P-C bonds in dialkylphosphonates lacking alpha-CF(3) groups (1.805-1.822A). Semi-empirical and non-empirical molecular modeling of diEt-FAP and (EtO)(2)P(O)C(CH(3))(2)NHS(O)(2)C(6)H(5) (diEt-AP), which lacks CF(3) groups, indicated lengthening and destabilization of the P-C bond in diEt-FAP compared to diEt-AP. Active site peptide adducts formed by reacting diEt-FAP with BChE and diBu-FAP with NTE catalytic domain (NEST) were identified using peptide mass mapping with mass spectrometry (MS). Mass shifts (mean+/-SE, average mass) for peaks corresponding to active site peptides with diethylphosphoryl and monoethylphosphoryl adducts on BChE were 136.1+/-0.1 and 108.0+/-0.1Da, respectively. Corresponding mass shifts for dibutylphosphoryl and monobutylphosphoryl adducts on NEST were 191.8+/-0.2 and 135.5+/-0.1Da, respectively. Each of these values was statistically identical to the theoretical mass shift for each dialkylphosphoryl and monoalkylphosphoryl species. The MS results demonstrate that inhibition of BChE and NEST by FAPs yields dialkylphosphoryl and monoalkylphosphoryl adducts, consistent with phosphorylation via P-C bond cleavage and aging by net dealkylation. Taken together, predictions from enzyme kinetics, chemical reactivity, X-ray crystallography, and molecular modeling were confirmed by MS and support the hypothesis that FAPs inhibit serine esterases via scission of the P-C bond.

Design of topological indices: computer-oriented approach.

A novel method is suggested for constructing topological indices (TIs) of molecular graphs which models human logic. This method is described in terms of a block scheme, consisting of the mutually connected elementary blocks. In each block the simple transformations of a molecular graph are fulfilled. A variant of the transformation is selected from the list of possible variants. Every TI is obtained as a result of the sequential execution of a number of operations, corresponding to some 'walk' on the block scheme. This walk can be selected both randomly and by the investigator. The suggested method can serve as a basis for the development of the respective computer program which may be used for the automatic construction of any number of TIs of differing nature. By this process one can also obtain the TIs that are unlikely to be constructed manually, due to their complexity. The set of obtained TIs may be used for building the structure-property models. In the case of an unsatisfactory result the obtained set of TIs may be changed using the described generator of TIs. A number of examples of application of the suggested approach for the building QSAR/QSPR models is given.