Stability estimation through multivariate approach among solasodine-rich lines of Solanum khasianum (C.B. Clarke): an important industrial plant.

Solanum khasianum is a medicinally important plant that is a source of steroidal alkaloids 'solasodine.' It has various industrial applications, including oral contraceptives and other pharmaceutical uses. The present study was based on 186 germplasm of S. khasianum, which were analyzed for the stability of economically important traits like solasodine content and fruit yield. The collected germplasm was planted during Kharif 2018, 2019, and 2020 in RCBD with three replications at the experimental farm of CSIR-NEIST, Jorhat, Assam, India. A multivariate approach for stability analysis was applied to identify stable germplasm of S. khasianum for economically important traits. The germplasm was analyzed for additive main effects and multiplicative interaction (AMMI), GGE biplot, multi-trait stability index, and Shukla's variance which were evaluated for three environments. The AMMI ANOVA revealed significant GE interaction for all the studied traits. The stable and high-yielding germplasm was identified from the AMMI biplot, GGE biplot, Shukla's variance value, and MTSI plot analysis. Lines no. 90, 85, 70, 107, and 62 were identified as highly stable fruit yielders while, lines no. 1, 146, and 68 were identified as stable high solasodine lines. However, considering both traits, i.e., high fruit yield and solasodine content, MTSI analysis was performed which showed that lines 1, 85, 70,155, 71, 114, 65, 86, 62, 116, 32, and 182 could be used in a breeding program. Thus, this identified germplasm can be considered for further varietal development and could be used in a breeding program. The findings of the present study would be beneficial for the S. khasianum breeding program.

Molecular property diagnostic suite (MPDS): Development of disease-specific open source web portals for drug discovery.

Molecular property diagnostic suite (MPDS) is a Galaxy-based open source drug discovery and development platform. MPDS web portals are designed for several diseases, such as tuberculosis, diabetes mellitus, and other metabolic disorders, specifically aimed to evaluate and estimate the drug-likeness of a given molecule. MPDS consists of three modules, namely data libraries, data processing, and data analysis tools which are configured and interconnected to assist drug discovery for specific diseases. The data library module encompasses vast information on chemical space, wherein the MPDS compound library comprises 110.31 million unique molecules generated from public domain databases. Every molecule is assigned with a unique ID and card, which provides complete information for the molecule. Some of the modules in the MPDS are specific to the diseases, while others are non-specific. Importantly, a suitably altered protocol can be effectively generated for another disease-specific MPDS web portal by modifying some of the modules. Thus, the MPDS suite of web portals shows great promise to emerge as disease-specific portals of great value, integrating chemoinformatics, bioinformatics, molecular modelling, and structure- and analogue-based drug discovery approaches.

How good are publicly available web services that predict bioactivity profiles for drug repurposing?

Drug repurposing provides a non-laborious and less expensive way for finding new human medicines. Computational assessment of bioactivity profiles shed light on the hidden pharmacological potential of the launched drugs. Currently, several freely available computational tools are available via the Internet, which predict multitarget profiles of drug-like compounds. They are based on chemical similarity assessment (ChemProt, SuperPred, SEA, SwissTargetPrediction and TargetHunter) or machine learning methods (ChemProt and PASS). To compare their performance, this study has created two evaluation sets, consisting of (1) 50 well-known repositioned drugs and (2) 12 drugs recently patented for new indications. In the first set, sensitivity values varied from 0.64 (TarPred) to 1.00 (PASS Online) for the initial indications and from 0.64 (TarPred) to 0.98 (PASS Online) for the repurposed indications. In the second set, sensitivity values varied from 0.08 (SuperPred) to 1.00 (PASS Online) for the initial indications and from 0.00 (SuperPred) to 1.00 (PASS Online) for the repurposed indications. Thus, this analysis demonstrated that the performance of machine learning methods surpassed those of chemical similarity assessments, particularly in the case of novel repurposed indications.

Differentiation of isomeric 2-aryldimethyltetrahydro-5-quinolinones by electron ionization and electrospray ionization mass spectrometry.

A series of isomeric 2-aryl-6,6-dimethyltetrahydro-5-quinolinones (set I) and 2-aryl-7,7-dimethyltetrahydro-5-quinolinones (set II) were studied under positive ion electron ionization (EI) and electrospray ionization (ESI) techniques. Under EI conditions, the molecular ions were found to be less stable in set I isomers, and they resulted in abundant fragment ions, i.e., [M-CH(3)](+), [M-CO](+.), [M-HCO](+), [M-(CH(3),CO)](+), and [M-(CH(3),CH(2)O)](+), when compared with set II isomers. In addition, the set I isomers showed specific fragment ions corresponding to [M-OH](+) and [M-OCH(3)](+). The retro-Diels-Alder (RDA) product ion was always higher in set II isomers. The ESI mass spectra produced [M + H](+) ions, and their decomposition showed favorable loss of CH(3) radical, CH(4) and C(2)H(6) molecules in set I isomers. The set II isomers, however, showed predominant RDA product ions, and specific loss of H(2)O. The selectivity in EI and ESI was attributed to the instability of set I isomers by the presence of a gem-dimethyl group at the α-position, and it was supported by the data from model compounds without a gem-dimethyl group. Density functional theory (DFT) calculations successfully corroborated the fragmentation pathways for diagnostic ions. This study revealed the effect of a gem-dimethyl group located at the α-position to the carbonyl having aromatic/unsaturated carbon on the other side of the carbonyl group.

The performance of computational techniques in locating the charge separated intermediates in organocatalytic transformations.

The formation of zwitterionic adducts between neutral nucleophiles such as NMe(3) and PMe(3) with neutral electrophiles such as methyl vinyl ketone (MVK) has been studied with a wide variety of theoretical methods. It has been found that hybrid density functional methods such as B3LYP are not capable of describing these zwitterionic structures as minima on the potential energy surface. This is also true for combinations of MP2 theory with basis sets lacking diffuse basis functions. The mPW1K hybrid functional, in contrast, correctly describes zwitterionic adducts as true intermediates on the PES. On the basis of this insight, a new version of the G3 compound energy scheme has been developed for the accurate description of zwitterionic structures. It has also been verified that modifications of the B2-PLYP double-hybrid functional are equally capable of the proper description of zwitterionic adducts. The applicability of this latter class of methods to a larger dataset involving combinations of different nucleophiles and electrophiles has been documented.

Molecular modeling studies of phenoxypyrimidinyl imidazoles as p38 kinase inhibitors using QSAR and docking.

p38 Kinase plays a vital role in inflammation mediated by tumor necrosis factor-alpha (TNFalpha) and interleukin-1beta (IL-1beta) pathways and inhibitors of p38 kinase provide effective approach for the treatment of inflammatory diseases. Pyridinyl and pyrimidinyl imidazoles, selectively inhibit p38alpha MAP kinase, are useful in the treatment of inflammatory diseases like rheumatoid arthritis. Three dimensional quantitative structure-activity relationship studies (3D-QSAR) involving comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA) and molecular docking were performed on 44 phenoxypyrimidinyl imidazole p38 kinase inhibitors to find out the structural relationship with the activity. The best predictive CoMFA model with atom fit alignment resulted in cross-validated r(2) value of 0.553, noncross-validated r(2) value of 0.908 and standard error of estimate 0.187. Similarly the best predictive CoMSIA model was derived with q(2) of 0.508, noncross-validated r(2) of 0.894 and standard error of estimate of 0.197, comprising steric, electrostatic, hydrophobic and hydrogen bond donor fields. These models were able to predict the activity of test set molecules efficiently within an acceptable error range. GOLD and FlexX were employed to dock the inhibitors into the active site of the p38 kinase and these docking studies revealed the vital interactions and binding conformation of the inhibitors. The information rendered by 3D-QSAR models and the docking interactions may afford valuable clues to optimize the lead and design new potential inhibitors.

Assessment of theoretical methods for the calculation of methyl cation affinities.

The methyl cation affinity (MCA; 298 K) of a variety of neutral and anionic bases has been examined computationally with a wide variety of theoretical methods. These include high-level composite procedures such as W1, G3, G3B3, and G2, conventional ab initio methods such as CCSD(T) and MP2, as well as a selection of density functional theory (DFT) methods. Experimental results for a variety of small model systems are well reproduced with practically all these methods, and the performance of DFT based methods are far superior in comparison to their MP2 analogs for these small models. For larger model, systems including motifs frequently encountered in organocatalysts, the performance deteriorates somewhat for DFT methods, while it improves significantly for MP2, rendering the former methods unreliable for common organic bases. Thus, MP2 calculations performed in combination with basis sets such as 6-31+G(2d, p) or larger, appear to offer a practical and reliable approach to compute MCAs of organic bases.

Heterobuckybowls: a theoretical study on the structure, bowl-to-bowl inversion barrier, bond length alternation, structure-inversion barrier relationship, stability, and synthetic feasibility.

Hybrid density functional theory (DFT) calculations at the B3LYP/cc-pVDZ level have been performed on a series of heterobuckybowls, 3X, C(18)X(3)H(6) (X = O, NH, CH(2), BH, S, PH, PH(3), Si, SiH(2), and AlH). The minimum energy conformations and the transition states for bowl-to-bowl inversion, where the geometry is bowl shaped, are computed and characterized by frequency calculations. The geometries of heterotrindenes, 2X, C(12)X(3)H(6) (X = O, NH, CH(2), BH, S, PH, PH(3), Si, SiH(2), and AlH), were obtained, and the bond length alternation (Delta) in the central benzenoid ring shows remarkable sensitivity as a function of substituent with a wide range of fluctuations (-0.014 to +0.092 A). The Delta computed in 2BH was found to be comparable with the highest bond alternation reported to date in benzenoid frameworks. The inversion dynamics of these heterobowls and their bowl depths were fit to a mixed quartic/quadratic function. The size of the heteroatom seems to exclusively control the bowl depth and rigidity as well as the synthetic feasibility. In contrast, the bond length alternation seems to be controlled by electronic factors and not by the size of the substituted atom either in trindenes or in heterosumanenes. The thermodynamic stability of this class of compounds is very much comparable with trithiasumanene (3S), which has been synthesized recently. The chemical hardness (eta) was measured to assess the stability of the heterosumanenes. The strain energy buildup in a sequential ring closure strategy along two synthetic routes, namely a triphenylene route and a trindene route, were explored, and the trindene route was found to be highly favorable for making such compounds compared to the triphenylene route. However, in both routes the ease of the synthetic feasibility increases as the size of the heteroatom increases.

Tailoring the curvature, bowl rigidity and stability of heterobuckybowls: theoretical design of synthetic strategies towards heterosumanenes.

Quantum mechanical calculations predict that larger heteroatom substituents on the periphery increase the feasibility of the crucial third ring closure in sumanene and are responsible for the accompanying modulations in the curvature, rigidity, stability and some of the physicochemical properties of the resulting heterosumanenes. Systematic application of semiempirical, ab initio, and DFT methods reveal that the qualitative trends obtained and our principal conclusions are independent of level of theory, albeit with minor quantitative differences.