Beyond Rotatable Bond Counts: Capturing 3D Conformational Flexibility in a Single Descriptor

The hapten rigidity improves antibody performances in immunoassay for rifamycins: Immunovalidation and molecular mechanism

The immunogenicity of haptens determines the performance of the resultant antibody for small molecules. Rigidity is one of the basic physicochemical properties of haptens. However, few studies have investigated the effect of hapten rigidity on the strength of an immune response and overall antibody performance. Herein, we introduce three molecular descriptors that quantify hapten rigidity. By using of these descriptors, four rifamycin haptens with varied rigidity were designed. The structural and physicochemical feasibility of the designed haptens was then assessed by computational chemistry. Immunization demonstrated that the strength of induced immune responses, i.e., the titer and affinity of antiserum, was significantly increased with increased rigidity of haptens. Furthermore, molecular dynamic simulations demonstrated conformation constraint of rigid haptens contributed to the initial binding and activation of naïve B cells. Finally, a highly sensitive indirect competitive enzyme-linked immunosorbent assay was developed for detection of rifaximin, with an IC50 of 1.1 μg/L in buffer and a limit of detection of 0.2-11.3 μg/L in raw milk, river water, and soil samples. This work provides new insights into the effect of hapten rigidity on immunogenicity and offers new hapten design strategies for antibody discovery and vaccine development of small molecules.

Computer-aided drug discovery of c-Abl kinase inhibitors from plant compounds against chronic myeloid leukemia

Chronic myeloid leukemia (CML) is a hematological malignancy characterized by the neoplastic transformation of hematopoietic stem cells, driven by the Philadelphia (Ph) chromosome resulting from a translocation between chromosomes 9 and 22. This Ph chromosome harbors the breakpoint cluster region (BCR) and the Abelson (ABL) oncogene (BCR-ABL1) which have a constitutive tyrosine kinase activity. However, the tyrosine kinase activity of BCR-ABL1 have been identified as a key player in CML initiation and maintenance through c-Abl kinase. Despite advancements in tyrosine kinase inhibitors, challenges such as efficacy, safety concerns, and recurring drug resistance persist. This study aims to discover potential c-Abl kinase inhibitors from plant compounds with anti-leukemic properties, employing drug-likeness assessment, molecular docking, in silico pharmacokinetics (ADMET) screening, density function theory (DFT), and molecular dynamics simulations (MDS). Out of 58 screened compounds for drug-likeness, 44 were docked against c-Abl kinase. The top hit compound (isovitexin) and nilotinib (control drug) were subjected to rigorous analyses, including ADMET profiling, DFT evaluation, and MDS for 100 ns. Isovitexin demonstrated a notable binding affinity (-15.492 kcal/mol), closely comparable to nilotinib (-16.826 kcal/mol), showcasing a similar binding pose and superior structural stability and reactivity. While these findings suggest isovitexin as a potential c-Abl kinase inhibitor, further validation through urgent in vitro and in vivo experiments is imperative. This research holds promise for providing an alternative avenue to address existing CML treatment and management challenges.Communicated by Ramaswamy H. Sarma.

Exploring the impact of alignment media on RDC analysis of phosphorus-containing compounds: a molecular docking approach

Residual dipolar couplings (RDCs) are employed in NMR analysis when conventional methods, such as J-couplings and nuclear Overhauser effects (NOEs) fail. Low-energy (optimized) conformers are often used as input structures in RDC analysis programs. However, these low-energy structures do not necessarily resemble conformations found in anisotropic environments due to interactions with the alignment medium, especially if the analyte molecules are flexible. Considering interactions with alignment media in RDC analysis, we developed and evaluated a molecular docking-based approach to generate more accurate conformer ensembles for compounds in the presence of the poly-γ-benzyl-L-glutamate alignment medium. We designed chiral phosphorus-containing compounds that enabled us to utilize 31P NMR parameters for the stereochemical analysis. Using P3D/PALES software to evaluate diastereomer discrimination, we found that our conformer ensembles outperform moderately the standard, low-energy conformers in RDC analysis. To further improve our results, we (i) averaged the experimental values of the molecular docking-based conformers by applying the Boltzmann distribution and (ii) optimized the structures through normal mode relaxation, thereby enhancing the Pearson correlation factor R and even diastereomer discrimination in some cases. Nevertheless, we presume that significant differences between J-couplings in isotropic and in anisotropic environments may preclude RDC measurements for flexible molecules. Therefore, generating conformer ensembles based on molecular docking enhances RDC analysis for mildly flexible systems while flexible molecules may require applying more advanced approaches, in particular approaches including dynamical effects.

Systematic Comparison of Experimental Crystallographic Geometries and Gas-Phase Computed Conformers for Torsion Preferences

We performed exhaustive torsion sampling on more than 3 million compounds using the GFN2-xTB method and performed a comparison of experimental crystallographic and gas-phase conformers. Many conformer sampling methods derive torsional angle distributions from experimental crystallographic data, limiting the torsion preferences to molecules that must be stable, synthetically accessible, and able to be crystallized. In this work, we evaluate the differences in torsional preferences of experimental crystallographic geometries and gas-phase computed conformers from a broad selection of compounds to determine whether torsional angle distributions obtained from semiempirical methods are suitable priors for conformer sampling. We find that differences in torsion preferences can be mostly attributed to a lack of available experimental crystallographic data with small deviations derived from gas-phase geometry differences. GFN2 demonstrates the ability to provide accurate and reliable torsional preferences that can provide a basis for new methods free from the limitations of experimental data collection. We provide Gaussian-based fits and sampling distributions suitable for torsion sampling and propose an alternative to the widely used "experimental torsion and knowledge distance geometry" (ETKDG) method using quantum torsion-derived distance geometry (QTDG) methods.

An Ecosystem for Digital Reticular Chemistry

The vastness of the materials design space makes it impractical to explore using traditional brute-force methods, particularly in reticular chemistry. However, machine learning has shown promise in expediting and guiding materials design. Despite numerous successful applications of machine learning to reticular materials, progress in the field has stagnated, possibly because digital chemistry is more an art than a science and its limited accessibility to inexperienced researchers. To address this issue, we present mofdscribe, a software ecosystem tailored to novice and seasoned digital chemists that streamlines the ideation, modeling, and publication process. Though optimized for reticular chemistry, our tools are versatile and can be used in nonreticular materials research. We believe that mofdscribe will enable a more reliable, efficient, and comparable field of digital chemistry.

In Silico Study Approach on a Series of 50 Polyphenolic Compounds in Plants; A Comparison on the Bioavailability and Bioactivity Data

Fifty (50) phytocompounds from several subclasses of polyphenols, chosen based on their abundance in the plant world, were analyzed through density functional methods, using computational tools to evaluate their oral availability and particular bioactivity on several cell modulators; key descriptors and molecular features related to the electron density and electrostatic potential for the lowest energy conformers of the investigated molecules were computed. An analysis of the bioactivity scores towards six cell modulators (GPCR ligand, ion channel modulator, kinase inhibitor, nuclear receptor ligand, protease inhibitor and enzyme inhibitor) was also achieved, in the context of investigating their potential side effects on the human digestive processes. Summarizing, computational results confirmed in vivo and in vitro data regarding the high bioavailability of soy isoflavones and better bioavailability of free aglycones in comparison with their esterified and glycosylated forms. However, by a computational approach analyzing Lipinski’s rule, apigenin and apigenin-7-O-rhamnoside, naringenin, hesperetin, genistein, daidzin, biochanin A and formonetin in the flavonoid series and all hydroxycinnamic acids and all hydroxybenzoic acids excepting ellagic acid were proved to have the best bioavailability data; rhamnoside derivatives, the predominant glycosides in green plants, which were reported to have the lowest bioavailability values by in vivo studies, were revealed to have the best bioavailability data among the studied flavonoids in the computational approach. Results of in silico screening on the phenolic derivatives series also revealed their real inhibitory potency on the six parameters studied, showing a remarkable similitude between the flavonoid series, while flavonoids were more powerful natural cell modulators than the phenyl carboxylic acids tested. Thus, it can be concluded that there is a need for supplementation with digestive enzymes, mainly in the case of individuals with low digestive efficiency, to obtain the best health benefits of polyphenols in humans.

Discovering Relationships between OSDAs and Zeolites through Data Mining and Generative Neural Networks

Organic structure directing agents (OSDAs) play a crucial role in the synthesis of micro- and mesoporous materials especially in the case of zeolites. Despite the wide use of OSDAs, their interaction with zeolite frameworks is poorly understood, with researchers relying on synthesis heuristics or computationally expensive techniques to predict whether an organic molecule can act as an OSDA for a certain zeolite. In this paper, we undertake a data-driven approach to unearth generalized OSDA-zeolite relationships using a comprehensive database comprising of 5,663 synthesis routes for porous materials. To generate this comprehensive database, we use natural language processing and text mining techniques to extract OSDAs, zeolite phases, and gel chemistry from the scientific literature published between 1966 and 2020. Through structural featurization of the OSDAs using weighted holistic invariant molecular (WHIM) descriptors, we relate OSDAs described in the literature to different types of cage-based, small-pore zeolites. Lastly, we adapt a generative neural network capable of suggesting new molecules as potential OSDAs for a given zeolite structure and gel chemistry. We apply this model to CHA and SFW zeolites generating several alternative OSDA candidates to those currently used in practice. These molecules are further vetted with molecular mechanics simulations to show the model generates physically meaningful predictions. Our model can automatically explore the OSDA space, reducing the amount of simulation or experimentation needed to find new OSDA candidates.

Novel Pyran-Linked Phthalazinone-Pyrazole Hybrids: Synthesis, Cytotoxicity Evaluation, Molecular Modeling, and Descriptor Studies

A series of novel pyran-linked phthalazinone-pyrazole hybrids were designed and synthesized by a facile one-pot three-component reaction employing substituted phthalazinone, 1H-pyrazole-5-carbaldehyde, and active methylene compounds. Optimization studies led to the identification of L-proline and ethanol as efficient catalyst and solvent, respectively. This was followed by evaluation of anticancer activity against solid tumor cell lines of lung and cervical carcinoma that displayed IC₅₀ values in the range of 9.8–41.6 µM. Molecular modeling studies were performed, and crucial interactions with the target protein were identified. The drug likeliness nature of the compounds and molecular descriptors such as molecular flexibility, complexity, and shape index were also calculated to understand the potential of the synthesized molecules to act as lead-like molecule upon further detailed biological investigations as well as 3D-QSAR studies.

One class classification as a practical approach for accelerating π-π co-crystal discovery

The implementation of machine learning models has brought major changes in the decision-making process for materials design. One matter of concern for the data-driven approaches is the lack of negative data from unsuccessful synthetic attempts, which might generate inherently imbalanced datasets. We propose the application of the one-class classification methodology as an effective tool for tackling these limitations on the materials design problems. This is a concept of learning based only on a well-defined class without counter examples. An extensive study on the different one-class classification algorithms is performed until the most appropriate workflow is identified for guiding the discovery of emerging materials belonging to a relatively small class, that being the weakly bound polyaromatic hydrocarbon co-crystals. The two-step approach presented in this study first trains the model using all the known molecular combinations that form this class of co-crystals extracted from the Cambridge Structural Database (1722 molecular combinations), followed by scoring possible yet unknown pairs from the ZINC15 database (21 736 possible molecular combinations). Focusing on the highest-ranking pairs predicted to have higher probability of forming co-crystals, materials discovery can be accelerated by reducing the vast molecular space and directing the synthetic efforts of chemists. Further on, using interpretability techniques a more detailed understanding of the molecular properties causing co-crystallization is sought after. The applicability of the current methodology is demonstrated with the discovery of two novel co-crystals, namely pyrene-6H-benzo[c]chromen-6-one (1) and pyrene-9,10-dicyanoanthracene (2).

Machine learning using one class classification on a database of existing co-crystals enables the identification of co-formers which are likely to form stable co-crystals, resulting in the synthesis of two co-crystals of polyaromatic hydrocarbons.

Flexibility in early drug discovery: focus on the beyond-Rule-of-5 chemical space

Large and flexible compounds are of interest in pharmaceutical programs aimed at challenging protein targets that cannot be modulated by Rule of Five (Ro5)-compliant small molecules. Given their particular structural features, early drug discovery is now in charge of identifying which molecular descriptors should be used in the often called beyond-Rule-of-5 (bRo5) chemical space. Here, we focus on flexibility descriptors. First, we discuss the concept of flexibility and then focus on the number of rotatable bonds (NRot), the most common in silico descriptor. After identifying the pros and cons of NRot, we discuss how Kier's index Φ can replace NRot, and the limits of 3D descriptors. Finally, we show how a misuse of NRot and Φ can result in incorrect interpretations of the impact of flexibility in the bRo5 space and how flexibility has potential in the prospective design of orally bioavailable bRo5 drug candidates.

Review and Modeling of Crystal Growth of Atropisomers from Solutions

In this paper, theories on anisotropic crystal growth and crystallization of atropisomers are reviewed and a model for anisotropic crystal growth from solution containing slow inter-converting conformers is presented. The model applies to systems with growth-dominated crystallization from solutions and assumes that only one conformation participates in the solute integration step and is present in the crystal lattice. Other conformers, defined as the wrong conformers, must convert to the right conformer before they can assemble to the crystal lattice. The model presents a simple implicit method for evaluating the growth inhibition effect by the wrong conformers. The crystal growth model applies to anisotropic growth in two main directions, namely a slow-growing face and a fast-growing face and requires the knowledge of solute crystal face integration coefficients in both directions. A parameter estimation algorithm was derived to extract those coefficients from data about temporal concentration and crystal size during crystallization and was designed to have a short run time, while providing a high-resolution estimation. The model predicts a size-dependent growth rate and simulations indicated that for a given seed size and solvent system and for an isothermal anti-solvent addition crystallization, the seed loading and the supersaturation at seeding are the main factors impacting the final aspect ratio. The model predicts a decrease of the growth inhibition effect by the wrong conformer with increasing temperature, likely due to faster equilibration between conformers and/or a decrease of the population of the wrong conformer, if of low energy, at elevated temperatures. Finally, the model predicts that solute surface integration becomes the rate-limiting mechanism for high solute integration activation energies, resulting in no impact of the WC on the overall crystal growth process.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

Privileged substructures for anti-sickling activity via cheminformatic analysis

Sickle cell disease (SCD), an autosomal recessive genetic disorder, has been recognized by the World Health Organization (WHO) as a major public health problem as it affects 300 000 individuals worldwide. Complications arising from SCD include anemia, microvascular occlusion, severe pain, stokes, renal dysfunction and infections. A lucrative therapeutic strategy is to employ anti-sickling agents that can disrupt the formation of the HbS polymer. This study therefore employed cheminformatic approaches, encompassing classification structure–activity relationship (CSAR) modeling, to deduce the privileged substructures giving rise to the anti-sickling activity of an investigated set of 115 compounds, followed by substructure analysis. Briefly, the compiled compounds were described by fingerprint descriptors and used in the construction of CSAR models via several machine learning algorithms. The modelability of the data set, as exemplified by the MODI index, was determined to be in the range of 0.70–0.84. The predictive performance was deduced by the accuracy, sensitivity, specificity and Matthews correlation coefficient, which was found to be statistically robust, whereby the former three parameters afforded values in excess of 0.7 while the latter statistical parameter provided a value greater than 0.5. An analysis of the top 20 important substructure descriptors for anti-sickling activity revealed that 10 important features were significant in the differentiation of actives from inactives, as illustrated by aromaticity/conjugation (e.g. SubFPC287, SubFPC171 and SubFPC5), carbonyl groups (e.g. SubFPC137, SubFPC139, SubFPC49 and SubFPC135) and miscellaneous groups (e.g. SubFPC303, SubFPC302 and SubFPC275). Furthermore, an analysis of the structure–activity relationship revealed that the length of alkyl chains, choice of functional moiety and position of substitution on the benzene ring may affect the anti-sickling activity of these compounds. Thus, this knowledge is anticipated to be useful for guiding the design of robust compounds against the gelling activity of HbS, as preliminarily demonstrated in the data-driven compound design presented herein.

Cheminformatic approaches (classification structure–activity relationship models based on 12 fingerprint classes) were employed for deducing privileged substructures giving rise to the anti-sickling activity of an investigated set of 115 compounds.

Application of two-dimensional binary fingerprinting methods for the design of selective Tankyrase I inhibitors

In the present study, five important binary fingerprinting techniques were used to model novel flavones for the selective inhibition of Tankyrase I. From the fingerprints used: the fingerprint atom pairs resulted in a statistically significant 2D QSAR model using a kernel-based partial least square regression method. This model indicates that the presence of electron-donating groups positively contributes to activity, whereas the presence of electron withdrawing groups negatively contributes to activity. This model could be used to develop more potent as well as selective analogues for the inhibition of Tankyrase I. Schematic representation of 2D QSAR work flow.