Computer-aided research on the relationship between molecular structures of natural compounds (NC) and their biological activities have been carried out extensively because the molecular structures of new drug candidates are usually analogous to or derived from the molecular structures of NC. In order to express the relationship physically realistically using a computer, it is essential to have a molecular descriptor set that can adequately represent the characteristics of the molecular structures belonging to the NC's chemical space. Although several topological descriptors have been developed to describe the physical, chemical, and biological properties of organic molecules, especially synthetic compounds, and have been widely used for drug discovery researches, these descriptors have limitations in expressing NC-specific molecular structures. To overcome this, we developed a novel molecular fingerprint, called Natural Compound Molecular Fingerprints (NC-MFP), for explaining NC structures related to biological activities and for applying the same for the natural product (NP)-based drug development. NC-MFP was developed to reflect the structural characteristics of NCs and the commonly used NP classification system. NC-MFP is a scaffold-based molecular fingerprint method comprising scaffolds, scaffold-fragment connection points (SFCP), and fragments. The scaffolds of the NC-MFP have a hierarchical structure. In this study, we introduce 16 structural classes of NPs in the Dictionary of Natural Product database (DNP), and the hierarchical scaffolds of each class were calculated using the Bemis and Murko (BM) method. The scaffold library in NC-MFP comprises 676 scaffolds. To compare how well the NC-MFP represents the structural features of NCs compared to the molecular fingerprints that have been widely used for organic molecular representation, two kinds of binary classification tasks were performed. Task I is a binary classification of the NCs in commercially available library DB into a NC or synthetic compound. Task II is classifying whether NCs with inhibitory activity in seven biological target proteins are active or inactive. Two tasks were developed with some molecular fingerprints, including NC-MFP, using the 1-nearest neighbor (1-NN) method. The performance of task I showed that NC-MFP is a practical molecular fingerprint to classify NC structures from the data set compared with other molecular fingerprints. Performance of task II with NC-MFP outperformed compared with other molecular fingerprints, suggesting that the NC-MFP is useful to explain NC structures related to biological activities. In conclusion, NC-MFP is a robust molecular fingerprint in classifying NC structures and explaining the biological activities of NC

In the past 2 years, since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), multiple SARS-CoV-2 variants have emerged. Whenever a new variant emerges, considerable time is required to analyze the binding affinity of the viral surface proteins to human angiotensin-converting enzyme 2 (hACE2) and monoclonal antibodies. To efficiently predict the binding affinities associated with hACE2 and monoclonal antibodies in a short time, herein, we propose a method applying statistical analysis to simulations performed using molecular and quantum mechanics. This method efficiently predicted the trend of binding affinity for the binding of the spike protein of each variant of SARS-CoV-2 to hACE2 and individually to eight commercial monoclonal antibodies. Additionally, this method accurately predicted interaction energy changes in the crystal structure for 10 of 13 mutated residues in the Omicron variant, showing a significant change in the interaction energy of hACE2. S375F was found to be a mutation that majorly changed the binding affinity of the spike protein to hACE2 and the eight monoclonal antibodies. Our proposed analysis method enables the prediction of the binding affinity of new variants to hACE2 or to monoclonal antibodies in a shorter time compared to that utilized by the experimental method.

No abstract

X-linked inhibitor of apoptosis protein (XIAP) is an important regulator of cancer cell survival whose BIR3 domain (XIAP-BIR3) recognizes the Smac N-terminal tetrapeptide sequence (AVPI), making it an attractive protein-protein interaction (PPI) target for cancer therapies. We used the fragment molecular orbital (FMO) method to study the binding modes and affinities between XIAP-BIR3 and a series of its inhibitors (1–8) that mimic the AVPI binding motif; the inhibitors had common interactions with key residues in a hot spot region of XIAP-BIR3 (P1–P4 subpockets) with increased binding affinity mainly attributed to specific interactions with the P1 and P4 subpockets. Based on the structural information from FMO results, we proposed a novel XIAP natural product inhibitor, neoeriocitrin 10, which was derived from our preciously reported XIAP-BIR3 inhibitor 9, can be used as a highly potent candidate for XIAP-BIR3 inhibition. We also performed pair interaction energy decomposition analysis to investigate the binding energies between specific binding residues and individual ligands, showing that the novel natural product neoeriocitrin 10 had a higher binding affinity than epicatechin gallate 9. Molecular docking and dynamics simulations were performed to explore the mode of binding between 10 and XIAP-BIR3, demonstrating that 10 binds more strongly to the P1 and P4 pockets than 9. Overall, we present a novel natural product, neoeriocitrin 10, and demonstrate that the FMO method can be used to identify hot spots in PPIs and design new compounds for XIAP inhibition.

Background Dog-associated infections are related to more than 70 human diseases. Given that the health diagnosis of a dog requires expertise of the veterinarian, an artificial intelligence model for detecting dog diseases could significantly reduce time and cost required for a diagnosis and efficiently maintain animal health. Objective We collected normal and multispectral images to develop classification model of each three dog skin diseases (bacterial dermatosis, fungal infection, and hypersensitivity allergic dermatosis). The single models (normal image- and multispectral image-based) and consensus models were developed used to four CNN model architecture (InceptionNet, ResNet, DenseNet, MobileNet) and select well-performed model. Results For single models, such as normal image- or multispectral image-based model, the best accuracies and Matthew’s correlation coefficients (MCCs) for validation data set were 0.80 and 0.64 for bacterial dermatosis, 0.70 and 0.36 for fungal infection, and 0.82 and 0.47 for hypersensitivity allergic dermatosis. For the consensus models, the best accuracies and MCCs for the validation set were 0.89 and 0.76 for the bacterial dermatosis data set, 0.87 and 0.63 for the fungal infection data set, and 0.87 and 0.63 for the hypersensitivity allergic dermatosis data set, respectively, which supported that the consensus models of each disease were more balanced and well-performed. Conclusions We developed consensus models for each skin disease for dogs by combining each best model developed with the normal and multispectral images, respectively. Since the normal images could be used to determine areas suspected of lesion of skin disease and additionally the multispectral images could help confirming skin redness of the area, the models achieved higher prediction accuracy with balanced performance between sensitivity and specificity.

Six kuwanon derivatives (A/B/C/E/H/J) extracted from the roots of Morus alba L. were evaluated to determine their cyclooxygenase (COX)-1 and 2 inhibitory effects. Cyclooxygenase (COX) is known as the target enzyme of nonsteroidal anti-inflammatory drugs (NSAIDs), which are the most widely used therapeutic agents for pain and inflammation. Among six kuwanon derivatives, kuwanon A showed selective COX-2 inhibitory activity, almost equivalent to that of celecoxib, a known COX inhibitor. Kuwanon A showed high COX-2 inhibitory activity (IC50 = 14 μM) and a selectivity index (SI) range of >7.1, comparable to celecoxib (SI > 6.3). To understand the mechanisms underlying this effect, we performed docking simulations, fragment molecular orbital (FMO) calculations, and pair interaction energy decomposition analysis (PIEDA) at the quantum-mechanical level. As a result, kuwanon A had the strongest interaction with Arg120 and Tyr355 at the gate of the COX active site (−7.044 kcal/mol) and with Val89 in the membrane-binding domain (−6.599 kcal/mol). In addition, kuwanon A closely bound to Val89, His90, and Ser119, which are residues at the entrance and exit routes of the COX active site (4.329 Å). FMO calculations and PIEDA well supported the COX-2 selective inhibitory action of kuwanon A. It showed that the simulation and modeling results and experimental evidence were consistent.

Using 126 endogenous lichen fungus (ELF) extracts, inhibitory activities against monoamine oxidases (MAOs) and cholinesterases (ChEs) were evaluated. Among them, extract ELF29 of the endogenous fungus Diaporthe mahothocarpus of the lichen Cladonia symphycarpia showed the highest inhibitory activity against hMAO-A. Compounds alternariol (AT), 5′-hydroxy-alternariol (HAT), and mycoepoxydiene (MED), isolated from the extract, had potent inhibitory activities against hMAO-A with IC50 values of 0.020, 0.31, and 8.68 µM, respectively. AT, HAT, and MED are reversible competitive inhibitors of hMAO-A with Ki values of 0.0075, 0.116, and 3.76 µM, respectively. The molecular docking studies suggested that AT, HAT, and MED had higher binding affinities for hMAO-A (−9.1, −6.9, and −5.6 kcal/mol, respectively) than for hMAO-B (−6.3, −5.2, and −3.7 kcal/mol, respectively). The relative tight binding might result from a hydrogen bond interaction of the three compounds with a Tyr444 residue in hMAO-A, whereas no hydrogen bond interaction was proposed in hMAO-B. In silico pharmacokinetics, the three compounds showed high gastrointestinal absorption without violating Lipinski’s five rules, but only MED showed high probability to cross the blood–brain barrier. These results suggest that AT, HAT, and MED are candidates for treating neuropsychiatric disorders, such as depression and cardiovascular disease.

To describe the physically realistic solvation free energy surface of a molecule in a solvent, a generalized version of the solvation free energy density (G-SFED) calculation method has been developed. In the G-SFED model, the contribution from the hydrogen bond (HB) between a solute and a solvent to the solvation free energy was calculated as the product of the acidity of the donor and the basicity of the acceptor of an HB pair. The acidity and basicity parameters of a solute were derived using the summation of acidities and basicities of the respective acidic and basic functional groups of the solute, and that of the solvent was experimentally determined. Although the contribution of HBs to the solvation free energy could be evenly distributed to grid points on the surface of a molecule, the G-SFED model was still inadequate to describe the angle dependency of the HB of a solute with a polarizable continuum solvent. To overcome this shortcoming of the G-SFED model, the contribution of HBs was formulated using the geometric parameters of the grid points described in the HB coordinate system of the solute. We propose an HB angle dependency incorporated into the G-SFED model, i.e., the G-SFED-HB model, where the angular-dependent acidity and basicity densities are defined and parametrized with experimental data. The G-SFED-HB model was then applied to calculate the solvation free energies of organic molecules in water, various alcohols and ethers, and the log P values of diverse organic molecules, including peptides and a protein. Both the G-SFED model and the G-SFED-HB model reproduced the experimental solvation free energies with similar accuracy, whereas the distributions of the SFED on the molecular surface calculated by the G-SFED and G-SFED-HB models were quite different, especially for molecules having HB donors or acceptors. Since the angle dependency of HBs was included in the G-SFED-HB model, the SFED distribution of the G-SFED-HB model is well described as compared to that of the G-SFED model.