Rank-based ant system method for non-linear QSPR analysis: QSPR studies of the solubility parameter.

The solubility parameter (δ) plays a unique role in the development of stable pharmaceutical formulations for assessing phase segregation during product synthesis. Understanding this parameter helps to determine how a drug substance will behave when processed or when dosed in vivo. The aim of this work was to develop a novel comprehensive yet rapid and accurate Quantitative Structure-Property Relationship (QSPR) method based on the rank-based ant system feature selection. The method was coupled with the multiple linear regression and support vector regression and applied to the assessment of solubility parameters for a diverse dataset of 1804 chemical compounds. The models were validated by solubility prediction of 360 test set compounds which were not used in building models. The developed models have high prediction power characterized by r (2) values 0.75 and 0.82, and RMSE values 1.96 and 1.65 (J/(cm(3)))(0.5) for the external test set. Various validation techniques and comparison results with the novel optimized support vector regression indicate that the developed models can be used to determine the solubility parameters for a diverse set of chemicals with an acceptable accuracy. The developed models can be beneficial for designing new chemical materials with desired solubility parameter values.

QSAR modeling of datasets with enantioselective compounds using chirality sensitive molecular descriptors.

Shape descriptors used in 3D QSAR studies naturally take into account chirality; however, for flexible and structurally diverse molecules such studies require extensive conformational searching and alignment. QSAR modeling studies of two datasets of fragrance compounds with complex stereochemistry using simple alignment-free chirality sensitive descriptors developed in our laboratories are presented. In the first investigation, 44 alpha-campholenic derivatives with sandalwood odor were represented as derivatives of several common structural templates with substituents numbered according to their relative spatial positions in the molecules. Both molecular and substituent descriptors were used as independent variables in MLR calculations, and the best model was characterized by the training set q2 of 0.79 and external test set r2 of 0.95. In the second study, several types of chirality descriptors were employed in combinatorial QSAR modeling of 98 ambergris fragrance compounds. Among 28 possible combinations of seven types of descriptors and four statistical modeling techniques, k nearest neighbor classification with CoMFA descriptors was initially found to generate the best models with the internal and external accuracies of 76 and 89%, respectively. The same dataset was then studied using novel atom pair chirality descriptors (cAP). The cAP are based on a modified definition of the atomic chirality, in which the seniority of the substituents is defined by their relative partial charge values: higher values correspond to higher seniorities. The resulting models were found to have higher predictive power than those developed with CoMFA descriptors; the best model was characterized by the internal and external accuracies of 82 and 94%, respectively. The success of modeling studies using simple alignment free chirality descriptors discussed in this paper suggests that they should be applied broadly to QSAR studies of many datasets when compound stereochemistry plays an important role in defining their activity.

Novel chirality descriptors derived from molecular topology.

Several series of novel chirality descriptors of chemical organic molecules have been introduced. The descriptors have been developed on the basis of conventional topological descriptors of molecular graphs. They include modified molecular connectivity indices, Zagreb group indices, extended connectivity, overall connectivity, and topological charge indices. These modified descriptors make use of an additional term called chirality correction, which is added to the vertex degrees of asymmetric atoms in a molecular graph. Chirality descriptors can be real or complex numbers. Advantages and drawbacks of different series of chirality descriptors are discussed. These descriptors circumvent the inability of conventional topological indices to distinguish chiral or enantiomeric isomers, which so far has been the major drawback of 2D descriptors as compared to true 3D descriptors (e.g., shape, molecular fields) of molecular structure. These novel chirality descriptors have been implemented in a quantitative structure-activity releationship (QSAR) study of a set of ecdysteroids with a high content of chiral and enantiomeric compounds using the k nearest neighbor QSAR method (kNN) recently developed in this laboratory. We show that the results of this study compare favorably with those obtained with the comparative molecular field analysis (CoMFA) applied to the same dataset. The novel chirality descriptors of molecular structure should find their applications in QSAR studies and related investigations of molecular sdatasets.

Molecular dataset diversity indices and their applications to comparison of chemical databases and QSAR analysis

A new mutual molecular dataset diversity index (MMDDI), individual molecular dataset diversity index (IMDDI), and volume ratio (VR) are proposed to assess molecular dataset diversity. MMDDI and IMDDI can serve as valuable instruments for selecting monomer pools for combinatorial synthesis and in decision making about acquiring new databases. MMDDI can also be used as one of the criteria to estimate the quality of quantitative structure-activity relationship (QSAR) models aimed at the prediction of biological activities. The indices can be calculated directly from molecular descriptor values. The procedures applied for MMDDI and IMDDI calculations allow one to automatically compile lists of compounds, which can simplify molecular diversity analyses and database searching. The information can also be used for forming training and test sets in QSAR analysis.

Validation of protein-based alignment in 3D quantitative structure-activity relationships with CoMFA models.

The predictive capabilities of protein-based alignment (PBA) and structure-based alignment (SBA) comparative molecular field analysis (CoMFA) models have been compared. 3D quantitative structure-activity relationship (3D QSAR) models have been derived for a series of N-benzylpiperidine derivatives which are potent acetylcholinesterase (AChE) inhibitors interesting for Alzheimer's disease. To establish a comparison with the classical SBA procedure, different assay models were derived by superposing ligand conformers that are docked to the AChE active site and by using the most active compound as the reference one. A Kohonen self organizing map (SOM) was applied to analyse the molecular diversity of the test set relative to that of the training set, in order to explain the influence of molecular diversity on the predictive power of the considered models. SBA 3D QSAR models have to be used to predict the inhibitory activity only for compounds belonging to subgroups included in the training set. The PBA 3D QSAR models appeared to have a higher predictability, even for compounds with a molecular diversity greater than that of the training set. This results from the fact that the protein helps to automatically select the active conformation which is fitting the 3D QSAR model.

A molecular mechanism for the cleavage of a disulfide bond as the primary function of agonist binding to G-protein-coupled receptors based on theoretical calculations supported by experiments.

A model of the binding site of delta-opioids in the extracellular region of the G-protein-coupled opioid receptor based on modelling studies is presented. The distance between Asp288 and the disulfide bridge (Cys121-Cys198) formed between the first and second extracellular loops was found to be short. This model is consistent with site-directed mutagenesis studies. The arrangement of the ligands found in the receptor led to the development of a reaction mechanism for the cleavage of the disulfide bond catalysed by the ligands. Semi-empirical quantum chemical PM3 and AM1 calculations as well as ab initio studies showed that the interaction between the carboxylic acid side chain of aspartic acid and the disulfide bond leads to the polarization of, and withdrawal of a proton from, the protonated nitrogen of the ligand to one of the sulfur atoms. A mixed sulfenic acid and carboxylic acid anhydrate is formed as an intermediate as well as a thiol. The accompanying cleavage of the disulfide bond may produce a conformational change in the extracellular loops such that the pore formed by the seven-helix bundle opens allowing entrance of the ligand, water and ions into the cell. Cleavage of the disulfide bond after opioid administration was demonstrated experimentally by flow-cytometric measurements employing CMTMR and monobromobimane-based analyses of membrane-located thiols. The suggested mechanism may explain, in a consistent way, the action of agonists and antagonists and is assumed to be common for many G-protein coupled receptors.

Probing of the receptor-binding sites of the H1 and H3 influenza A and influenza B virus hemagglutinins by synthetic and natural sialosides.

To compare features of the receptor-binding sites (RBSs) of different influenza virus hemagglutinins (HA), binding of a number of synthetic sialic acid (SA) analogs and natural sialosides by a panel of about 30 human influenza A and B virus strains was studied in a competitive ligand binding assay. For all the viruses tested, the N-acetyl group of Neu5Ac, as well as the natural orientation of the carboxylic group at C2 and the hydroxylic group at C4, was essential for binding. Significant type- and subtype-specific differences were observed in virus recognition of asialic parts of sialosides. H1 strains, unlike H3 and type B viruses, were found to bind alpha 2-6-sialyl-N-acetyllactosamine with about an order of magnitude higher affinity than alpha 2-6-sialyllactose (6'SL). The H1 viruses and the H3 strains with Gln in position 226 of HA, but not the H3 strains with Leu-226, bound 6'SL with a lower affinity than alpha 2-3-sialyllactose; this effect correlated clearly with the preferential binding by the former viruses of unsubstituted alpha Neu5Ac compared to methyl alpha-glycoside of Neu5Ac. Thus, differentiation between the types of the SA-Gal linkage by the A viruses appeared to depend, at least partially, upon the recognition by the HA of the first hydrocarbon group of the aglycon. Type B virus strains were distinct in having a lower affinity for the Neu5Ac moiety and in providing a higher contribution of the asialic portions of sialosides to the HA-ligand interactions. The last effects are presumably due to the amino acid insertions in the type B HA surrounding the RBS, which makes the receptor-binding pocket deeper. The results obtained in the present investigation indicate that while the functional groups of Neu5Ac studied are recognized by the RBSs of all influenza viruses, the magnitude of their contribution to the binding energy, as well as the contribution of the asialic portion of the receptor, may vary in dependence upon the virus type, subtype, and strain.

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. III. Probable binding conformations of mu-agonists with phenylalanine in position 3.

Theoretical conformational analysis was carried out for several tetrapeptide analogues of beta-casomorphin and dermorphin containing a Phe residue in position 3. Sets of low-energy backbone structures of the mu-selective peptides [N-Me-Phe3, D-Pro4]-morphiceptin and Tyr-D-Orn-Phe-Asp-NH2 were obtained. These sets of structures were compared for geometrical similarity between themselves and with the low-energy conformations found for the delta-selective peptide Tyr-D-Cys-Phe-D-Pen-OH and nonactive peptide Tyr-Orn-Phe-Asp-NH2. Two pairs of geometrically similar conformations of mu-selective peptides, sharing no similarity with the conformations of peptides showing low affinity to the mu-receptor, were selected as two alternative models of probable mu-receptor-bound backbone conformations. Both models share geometrical similarity with the low-energy structures of the linear mu-selective peptide Tyr-D-Ala-Phe-Phe-NH2. Putative binding conformations of Tyr1 and Phe3 side chains are also discussed.

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. II. Models of mu- and delta-receptor-bound structures for analogues containing Phe4.

Models of mu- and delta-receptor-bound backbone conformations of enkephalin cyclic analogues containing Phe4 were determined by comparing geometrical similarity among the previously found low-energy backbone structures of [D-Cys2,Cys5]-enkephalinamide, [D-Cys2,D-Cys5]-enkephalinamide, [D-Pen2,L-Pen5]-enkephalin and [D-Pen2,D-Pen5]-enkephalin. The present mu-receptor-bound conformation resembles a beta-I bend in the peptide backbone centred on the Gly3-Phe4 region. Two slightly different models were found for the delta-receptor-bound conformation; both of them are more extended than the mu-receptor-bound conformation and include a gamma-turn (or a gamma-like turn) on the Gly3 residue. Energetically favourable rotamers of Tyr and Phe side chains were also determined for the mu- and delta-conformations. The present models of mu- and delta-conformations share geometrical similarity with the low-energy structures of Leu-enkephalin and the Tyr-D-Lys-Gly-Phe-analogue.

Conformational features responsible for binding of cyclic analogues of enkephalin to opioid receptors. I. Low-energy peptide backbone conformers of analogues containing Phe4.

Low-energy peptide backbone conformers were found by means of energy calculation for several cyclic analogues of enkephalin in an attempt to assess models for receptor-bound conformations for opioid receptors of the mu- and delta-types. They included [D-Cys2, L-Cys5]- and [D-Cys2, D-Cys5]-enkephalinamides showing moderate preference for mu-receptors, the delta-selective compounds [D-Pen2, L-Pen5] and [D-Pen2, D-Pen5]-enkephalins and Tyr-D-Lys-Gly-Phe- analogue possessing very high affinity to receptors of the mu-type. The low-energy conformers obtained for these analogues were in good agreement with the results of calculations by other authors and with experimental evidence. All of the analogues contain a Phe residue in position 4 of the peptide chain which facilitates the eventual search for geometrical similarity between the low-energy backbone conformers of different analogues in question.

Conformational aspects of antiviral activity of deoxyguanosine acyclic analogues.

Conformational possibilities of a series of deoxyguanosine analogues possessing or lacking antiviral activity were evaluated using methods of the molecular mechanics. Comparison of the spatial structures of acyclic analogues with one another and with the spatial structures of deoxyguanosine demonstrates restricted conformational mobility for compounds devoid of activity. The level of sterically allowed superposition of functional groups from the acyclic moieties of analogues and the corresponding atomic centres of deoxyribose could serve as a criterion of activity. The superposition could be performed in two different ways through either of the nonhydrogen substituents at the C1' atom in the five-membered ring.