Evaluating cepharanthine analogues as natural drugs against SARS-CoV-2.

Cepharanthine (CEP) is a natural biscoclaurine alkaloid of plant origin and was recently demonstrated to have anti-severe acute respiratory syndrome coronavirus 2 (anti-SARS-CoV-2) activity. In this study, we evaluated whether natural analogues of CEP may act as potential anti-coronavirus disease 2019 drugs. A total of 24 compounds resembling CEP were extracted from the KNApSAcK database, and their binding affinities to target proteins, including the spike protein and main protease of SARS-CoV-2, NPC1 and TPC2 in humans, were predicted via molecular docking simulations. Selected analogues were further evaluated by a cell-based SARS-CoV-2 infection assay. In addition, the efficacies of CEP and its analogue tetrandrine were assessed. A comparison of the docking conformations of these compounds suggested that the diphenyl ester moiety of the molecules was a putative pharmacophore of the CEP analogues.

Knowledge-based structural models of SARS-CoV-2 proteins and their complexes with potential drugs.

The World Health Organization (WHO) has declared the coronavirus disease 2019 (COVID-19) caused by the novel coronavirus SARS-CoV-2 a pandemic. There is, however, no confirmed anti-COVID-19 therapeutic currently. In order to assist structure-based discovery efforts for repurposing drugs against this disease, we constructed knowledge-based models of SARS-CoV-2 proteins and compared the ligand molecules in the template structures with approved/experimental drugs and components of natural medicines. Our theoretical models suggest several drugs, such as carfilzomib, sinefungin, tecadenoson, and trabodenoson, that could be further investigated for their potential for treating COVID-19.

Classification of ligand molecules in PDB with graph match-based structural superposition.

The fast heuristic graph match algorithm for small molecules, COMPLIG, was improved by adding a structural superposition process to verify the atom-atom matching. The modified method was used to classify the small molecule ligands in the Protein Data Bank (PDB) by their three-dimensional structures, and 16,660 types of ligands in the PDB were classified into 7561 clusters. In contrast, a classification by a previous method (without structure superposition) generated 3371 clusters from the same ligand set. The characteristic feature in the current classification system is the increased number of singleton clusters, which contained only one ligand molecule in a cluster. Inspections of the singletons in the current classification system but not in the previous one implied that the major factors for the isolation were differences in chirality, cyclic conformations, separation of substructures, and bond length. Comparisons between current and previous classification systems revealed that the superposition-based classification was effective in clustering functionally related ligands, such as drugs targeted to specific biological processes, owing to the strictness of the atom-atom matching.

Detecting structural similarity of ligand interactions in the lipid metabolic system including enzymes, lipid-binding proteins and nuclear receptors.

Nuclear receptors, intracellular lipid-binding proteins and metabolic enzymes are responsible for optimal metabolic homeostasis in higher organisms. Recent studies revealed the specific cooperation/competition among the subfamilies of these proteins. In this study, the nuclear receptor-lipid-binding protein-enzyme system, in which the interactions are mostly mediated by ligand molecules, was examined in terms of their ligand-binding structures to detect the similarity of interactions between functionally related subfamilies. The complex structures were dissected into single amino acid motifs for ligand fragment binding, and the presence and evolutionary origin of the motifs were compared among the protein families. As a result, functionally related nuclear receptor and enzyme pairs were found to share more motifs than expected, in agreement with the fact that the two families compete for the same ligand, and thus our study implies the possible co-evolution of the indirectly interacting protein system.

Structure based studies of the adaptive diversification process of congerins.

The isoforms of a fish galectin, congerins I and II, have several features that make them suitable for a study of accelerated process of molecular diversification based on 3D structures: They have been generated by a gene duplication, and still maintain 47% amino acid sequence identity to each other. Their genes show very high K A: /K S: ratio, and are though to be components of fish defense system. The crystal systems for a high-resolution analysis are known for both proteins. A series of works with biochemistry, molecular biology, and X-ray crystallography techniques have suggested that the two proteins might have evolved under differential selection pressures. Congerin I appeared to be a stabilized version of galectin-1. Congerin II was shown to be adapted to a new carbohydrate-ligand. The 3D structures of the wild type and mutant proteins have revealed the probable cause and consequence of the selection pressure responsible for the diversification of congerins.

In vitro evolutionary thermostabilization of congerin II: a limited reproduction of natural protein evolution by artificial selection pressure.

The thermostability of the conger eel galectin, congerin II, was improved by in vitro evolutionary protein engineering. Two rounds of random PCR mutagenesis and selection experiments increased the congerin II thermostability to a level comparative to its naturally thermostable isoform, congerin I. The crystal structures of the most thermostable double mutant, Y16S/T88I, and the related single mutants, Y16S and T88I, were determined at 2.0 angstroms, 1.8 angstroms, and 1.6 angstroms resolution, respectively. The exclusion of two interior water molecules by the Thr88Ile mutation, and the relief of adjacent conformational stress by the Tyr16Ser mutation were the major contributions to the thermostability. These features in the congerin II mutants are similar to those observed in congerin I. The natural evolution of congerin genes, with the K(A)/K(S) ratio of 2.6, was accelerated under natural selection pressures. The thermostabilizing selection pressure artificially applied to congerin II mimicked the implied natural pressure on congerin I. The results showed that the artificial pressure made congerin II partially reproduce the natural evolution of congerin I.

An empirical approach for structure-based prediction of carbohydrate-binding sites on proteins.

A computer program system was developed to predict carbohydrate-binding sites on three-dimensional (3D) protein structures. The programs search for binding sites by referring to the empirical rules derived from the known 3D structures of carbohydrate-protein complexes. A total of 80 non-redundant carbohydrate-protein complex structures were selected from the Protein Data Bank for the empirical rule construction. The performance of the prediction system was tested on 50 known complex structures to determine whether the system could detect the known binding sites. The known monosaccharide-binding sites were detected among the best three predictions in 59% of the cases, which covered 69% of the polysaccharide-binding sites in the target proteins, when the performance was evaluated by the overlap between residue patches of predicted and known binding sites.

Crystal structure of a conger eel galectin (congerin II) at 1.45A resolution: implication for the accelerated evolution of a new ligand-binding site following gene duplication.

The crystal structure of congerin II, a galectin family lectin from conger eel, was determined at 1.45A resolution. The previously determined structure of its isoform, congerin I, had revealed a fold evolution via strand swap; however, the structure of congerin II described here resembles other prototype galectins. A comparison of the two congerin genes with that of several other galectins suggests acceralated evolution of both congerin genes following gene duplication. The presence of a Mes (2-[N-morpholino]ethanesulfonic acid) molecule near the carbohydrate-binding site in the crystal structure points to the possibility of an additional binding site in congerin II. The binding site consists of a group of residues that had been replaced following gene duplication suggesting that the binding site was built under selective pressure. Congerin II may be a protein specialized for biological defense with an affinity for target carbohydrates on parasites' cell surface.

The speciation of conger eel galectins by rapid adaptive evolution.

Many cases of accelerated evolution driven by positive Darwinian selection are identified in the genes of venomous and reproductive proteins. This evolutional phenomenon might have important consequences in their gene-products' functions, such as multiple specific toxins for quick immobilization of the prey and the establishment of barriers to fertilization that might lead to speciation, and in the molecular evolution of novel genes. Recently, we analyzed the molecular evolution of two galectins isolated from the skin mucus of conger eel (Conger myriaster), named congerins I and II, by cDNA cloning and X-ray structural analysis, and we found that they have evolved in the rapid adaptive manner to emergence of a new structure including strand-swapping and a unique new ligand-binding site. In this review article we summarize and discuss the molecular evolution, especially the rapid adaptive evolution, and the structure-function relationships of conger eel galectins.