Structural and spectroscopic characterization of RufO indicates a new biological role in rufomycin biosynthesis.

Rufomycins constitute a class of cyclic heptapeptides isolated from actinomycetes. They are secondary metabolites that show promising treatment against Mycobacterium tuberculosis infections by inhibiting a novel drug target. Several nonproteinogenic amino acids are integrated into rufomycins, including a conserved 3-nitro-tyrosine. RufO, a cytochrome P450 (CYP)-like enzyme, was proposed to catalyze the formation of 3-nitro-tyrosine in the presence of O2 and NO. To define its biological function, the interaction between RufO and the proposed substrate tyrosine is investigated using various spectroscopic methods that are sensitive to the structural change of a heme center. However, a low- to high-spin state transition and a dramatic increase in the redox potential that are commonly found in CYPs upon ligand binding have not been observed. Furthermore, a 1.89-Å crystal structure of RufO shows that the enzyme has flexible surface regions, a wide-open substrate access tunnel, and the heme center is largely exposed to solvent. Comparison with a closely related nitrating CYP reveals a spacious and hydrophobic distal pocket in RufO, which is incapable of stabilizing a free amino acid. Molecular docking validates the experimental data and proposes a possible substrate. Collectively, our results disfavor tyrosine as the substrate of RufO and point to the possibility that the nitration occurs during or after the assembly of the peptides. This study indicates a new function of the unique nitrating enzyme and provides insights into the biosynthesis of nonribosomal peptides.

Inhibition of N-myristoyltransferase activity promotes androgen receptor degradation in prostate cancer.

BACKGROUND: Even though prostate cancer (PCa) patients initially respond to androgen deprivation therapy, some will eventually develop castration resistant prostate cancer (CRPC). Androgen receptor (AR) mediated cell signaling is a major driver in the progression of CRPC while only a fraction of PCa becomes AR negative. This study aimed to understand the regulation of AR levels by N-myristoyltransferase in PCa cells. METHODS: Two enantiomers, (1S,2S)- d-NMAPPD and (1R,2R)- d-NMAPPD (LCL4), were characterized by various methods (1 H and 13 C NMR, UHPLC, high-resolution mass spectra, circular dichroism) and evaluated for the ability to bind to N-myristoyltransferase 1 (NMT1) using computational docking analysis. structure-activity relationship analysis of these compounds led to the synthesis of (1R,2R)-LCL204 and evaluation as a potential NMT1 inhibitor utilizing the purified full length NMT1 enzyme. The NMT inhibitory activity wase determined by Click chemistry and immunoblotting. Regulation of NMT1 on tumor growth was evaluated in a xenograft tumor model. RESULTS: (1R,2R)- d-NMAPPD, but not its enantiomer (1S,2S)- d-NMAPPD, inhibited NMT1 activity and reduced AR protein levels. (1R,2R)-LCL204, a derivative of (1R,2R)- d-NMAPPD, inhibited global protein myristoylation. It also suppressed protein levels, nuclear translocation, and transcriptional activity of AR full-length or variants in PCa cells. This was due to enhanced ubiquitin and proteasome-mediated degradation of AR. Knockdown of NMT1 levels inhibited tumor growth and proliferation of cancer cells. CONCLUSION: Inhibitory efficacy on N-myristoyltransferase activity by d-NMAPPD is stereospecific. (1R,2R)-LCL204 reduced global N-myristoylation and androgen receptor protein levels at low micromolar concentrations in prostate cancer cells. pharmacological inhibition of NMT1 enhances ubiquitin-mediated proteasome degradation of AR. This study illustrates a novel function of N-myristoyltransferase and provides a potential strategy for treatment of CRPC.

Nicardipine is a putative EED inhibitor and has high selectivity and potency against chemoresistant prostate cancer in preclinical models.

BACKGROUND: It is imperative to develop novel therapeutics to overcome chemoresistance, a significant obstacle in the clinical management of prostate cancer (PCa) and other cancers. METHODS: A phenotypic screen was performed to identify novel inhibitors of chemoresistant PCa cells. The mechanism of action of potential candidate(s) was investigated using in silico docking, and molecular and cellular assays in chemoresistant PCa cells. The in vivo efficacy was evaluated in mouse xenograft models of chemoresistant PCa. RESULTS: Nicardipine exhibited high selectivity and potency against chemoresistant PCa cells via inducing apoptosis and cell cycle arrest. Computational, molecular, and cellular studies identified nicardipine as a putative inhibitor of embryonic ectoderm development (EED) protein, and the results are consistent with a proposed mechanism of action that nicardipine destabilised enhancer of zeste homologue 2 (EZH2) and inhibited key components of noncanonical EZH2 signalling, including transducer and activator of transcription 3, S-phase kinase-associated protein 2, ATP binding cassette B1, and survivin. As a monotherapy, nicardipine effectively inhibited the skeletal growth of chemoresistant C4-2B-TaxR tumours. As a combination regimen, nicardipine synergistically enhanced the in vivo efficacy of docetaxel against C4-2 xenografts. CONCLUSION: Our findings provided the first preclinical evidence supporting nicardipine as a novel EED inhibitor that has the potential to be promptly tested in PCa patients to overcome chemoresistance and improve clinical outcomes.

Exploring and engineering PAM-diverse Streptococci Cas9 for PAM-directed bifunctional and titratable gene control in bacteria.

The RNA-guided Cas9s serve as powerful tools for programmable gene editing and regulation; their targeting scopes and efficacies, however, are always constrained by the PAM sequence stringency. Most Streptococci Cas9s, including the prototype SpCas9 from S. pyogenes, specifically recognize a canonical NGG PAM via a conserved RxR PAM-binding motif within the PAM-interaction (PI) domain. Here, SpCas9-based mining unveils three distinct and rarely presented PAM-binding motifs (QxxxR, QxQ and RxQ) among Streptococci Cas9 orthologs. With the catalytically-dead QxxxR-containing SedCas9 from S. equinus, we dissect its NAG PAM specificity and elucidate its underlying recognition mechanism via computational prediction and mutagenesis analysis. Replacing the SedCas9 PI domain with alternate PAM-binding motifs rewires its PAM specificity to NGG or NAA. Moreover, a semi-rational design with minimal mutation creates a SedCas9-NQ variant showing robust activity towards expanded NNG and NAA PAMs, based upon which we engineered a compact ω-SedCas9-NQ transcriptional regulator for PAM-directed bifunctional and titratable gene control. The ω-SedCas9-NQ mediated metabolic reprogramming of endogenous genes in Escherichia coli affords a 2.6-fold increase of 4-hydroxycoumarin production. This work reveals new Cas9 scaffolds with distinct PAM-binding motifs for PAM relaxation and creates a new PAM-diverse Cas9 variant for versatile gene control in bacteria.

Identifying Drug Candidates for COVID-19 with Large-Scale Drug Screening.

Papain-like protease (PLpro) is critical to COVID-19 infection. Therefore, it is a significant target protein for drug development. We virtually screened a 26,193 compound library against the PLpro of SARS-CoV-2 and identified several drug candidates with convincing binding affinities. The three best compounds all had better estimated binding energy than those of the drug candidates proposed in previous studies. By analyzing the docking results for the drug candidates identified in this and previous studies, we demonstrate that the critical interactions between the compounds and PLpro proposed by the computational approaches are consistent with those proposed by the biological experiments. In addition, the predicted binding energies of the compounds in the dataset showed a similar trend as their IC50 values. The predicted ADME and drug-likeness properties also suggested that these identified compounds can be used for COVID-19 treatment.

Morphological and physiological correlates of among- individual variation in basal metabolic rate in two passerine birds.

Basal metabolic rate (BMR) has been shown to be a highly phenotypic flexibility trait within species. A significant proportion of an individual's energy budget is accounted for by BMR, hence among-individual variation in this trait may affect other energetic processes, as well as fitness. In this study, we measured BMR, organ mass, mitochondrial respiration capacities and cytochrome c oxidase (COX) activities in muscle and liver and circulating levels of plasma triiodothyronine (T3) in Chinese bulbuls (Pycnonotus sinensis) and Eurasian tree sparrows (Passer montanus). Our results showed that heart and kidney mass was positively correlated with BMR in Chinese bulbuls, whereas liver and kidney mass was positively correlated with BMR in Eurasian tree sparrows. Regarding metabolic biochemical markers of tissues, state 4 respiration and COX activity in the muscles of the Chinese bulbuls was correlated with BMR, while state 4 respiration in the muscle and liver was correlated with BMR in Eurasian tree sparrows. T3 was significantly and positively correlated with BMR in Chinese bulbuls and Eurasian tree sparrows. Consistent with the above results, our findings suggest that T3 levels play an important role in modulating BMR in Chinese bulbuls and Eurasian tree sparrows. Moreover, individual variation in BMR can be explained partly by morphological and physiological mechanisms.

New therapeutics beyond amyloid-ß and tau for the treatment of Alzheimer's disease.

As the population ages, Alzheimer's disease (AD), the most common neurodegenerative disease in elderly people, will impose social and economic burdens to the world. Currently approved drugs for the treatment of AD including cholinesterase inhibitors (donepezil, rivastigmine, and galantamine) and an N-methyl-D-aspartic acid receptor antagonist (memantine) are symptomatic but poorly affect the progression of the disease. In recent decades, the concept of amyloid-ß (Aß) cascade and tau hyperphosphorylation leading to AD has dominated AD drug development. However, pharmacotherapies targeting Aß and tau have limited success. It is generally believed that AD is caused by multiple pathological processes resulting from Aß abnormality, tau phosphorylation, neuroinflammation, neurotransmitter dysregulation, and oxidative stress. In this review we updated the recent development of new therapeutics that regulate neurotransmitters, inflammation, lipid metabolism, autophagy, microbiota, circadian rhythm, and disease-modified genes for AD in preclinical research and clinical trials. It is to emphasize the importance of early diagnosis and multiple-target intervention, which may provide a promising outcome for AD treatment.

4-Methoxydalbergione is a potent inhibitor of human astroglioma U87 cells in vitro and in vivo.

Astroglioma is the most common primary tumor in the central nervous system without effective treatment strategies. Temozolomide (TMZ) is a chemotherapeutic drug to treat astroglioma but exhibits low potency and has side effects. Therefore, there is an urgent need to develop new compounds to treat astroglioma. Dalbergia sissoo Roxb was the source of Dalbergia odorifera in traditional Chinese medicine (TCM) and has been clinically used as an anti-tumor medicine. 4-Methoxydalbergione (4MOD) is purified from Dalbergia sissoo Roxb., and shows an inhibitory effect on osteosarcoma, but its effects on astroglioma have not been reported. Here, we evaluate its anti-astroglioma effects on both in vitro and in vivo models. In cultured astroglioma U87 cells, 4MOD inhibited cell proliferation and induced cell apoptosis in a time- and concentration-dependent manner. Compared with TMZ, 4MOD exhibited a tenfold greater potency of anti-astroglioma effects. 4MOD effectively stalled the cell cycle in G2 phase. Transcriptome sequencing (RNA-seq) showed that 4MOD upregulated 158 genes and downregulated 204 genes that are mainly enriched in cell membrane, cell division, cell cycle, p53, TNF, and MAPK signaling pathways, which may underlie its anti-tumor mechanisms. In a nude mouse xenograft model transplanted with U87 cells, 10 mg/kg 4MOD slowed down tumor growth rate, while at 30 mg/kg dose, it reduced tumor size. Collectively, this study demonstrates that 4MOD is a potent native compound that remarkably inhibits U87 astroglioma growth in both in vitro and in vivo models.

DeepDISE: DNA Binding Site Prediction Using a Deep Learning Method.

It is essential for future research to develop a new, reliable prediction method of DNA binding sites because DNA binding sites on DNA-binding proteins provide critical clues about protein function and drug discovery. However, the current prediction methods of DNA binding sites have relatively poor accuracy. Using 3D coordinates and the atom-type of surface protein atom as the input, we trained and tested a deep learning model to predict how likely a voxel on the protein surface is to be a DNA-binding site. Based on three different evaluation datasets, the results show that our model not only outperforms several previous methods on two commonly used datasets, but also demonstrates its robust performance to be consistent among the three datasets. The visualized prediction outcomes show that the binding sites are also mostly located in correct regions. We successfully built a deep learning model to predict the DNA binding sites on target proteins. It demonstrates that 3D protein structures plus atom-type information on protein surfaces can be used to predict the potential binding sites on a protein. This approach should be further extended to develop the binding sites of other important biological molecules.

Polyphenols as Potential Inhibitors of SARS-CoV-2 RNA Dependent RNA Polymerase (RdRp).

An increasing number of studies have demonstrated the antiviral nature of polyphenols, and many polyphenols have been proposed to inhibit SARS-CoV or SARS-CoV-2. Our previous study revealed the inhibitory mechanisms of polyphenols against DNA polymerase α and HIV reverse transcriptase to show that polyphenols can block DNA elongation by competing with the incoming NTPs. Here we applied computational approaches to examine if some polyphenols can also inhibit RNA polymerase (RdRp) in SARS-CoV-2, and we identified some better candidates than remdesivir, the FDA-approved drug against RdRp, in terms of estimated binding affinities. The proposed compounds will be further examined to develop new treatments for COVID-19.

Structural and Functional Enrichment Analyses for Antimicrobial Peptides.

Whether there is any inclination between structures and functions of antimicrobial peptides (AMPs) is a mystery yet to be unraveled. AMPs have various structures associated with many different antimicrobial functions, including antibacterial, anticancer, antifungal, antiparasitic and antiviral activities. However, none has yet reported any antimicrobial functional tendency within a specific category of protein/peptide structures nor any structural tendency of a specific antimicrobial function with respect to AMPs. Here, we examine the relationships between structures categorized by three structural classification methods (CATH, SCOP, and TM) and seven antimicrobial functions with respect to AMPs using an enrichment analysis. The results show that antifungal activities of AMPs were tightly related to the two-layer sandwich structure of CATH, the knottin fold of SCOP, and the first structural cluster of TM. The associations with knottin and TM Cluster 1 even sustained through the AMPs with a low sequence identity. Moreover, another significant mutual enrichment was observed between the third cluster of TM and anti-Gram-positive-bacterial/anti-Gram-negative-bacterial activities. The findings of the structure-function inclination further our understanding of AMPs and could help us design or discover new therapeutic potential AMPs.

Understand protein functions by comparing the similarity of local structural environments.

The three-dimensional structures of proteins play an essential role in regulating binding between proteins and their partners, offering a direct relationship between structures and functions of proteins. It is widely accepted that the function of a protein can be determined if its structure is similar to other proteins whose functions are known. However, it is also observed that proteins with similar global structures do not necessarily correspond to the same function, while proteins with very different folds can share similar functions. This indicates that function similarity is originated from the local structural information of proteins instead of their global shapes. We assume that proteins with similar local environments prefer binding to similar types of molecular targets. In order to testify this assumption, we designed a new structural indicator to define the similarity of local environment between residues in different proteins. This indicator was further used to calculate the probability that a given residue binds to a specific type of structural neighbors, including DNA, RNA, small molecules and proteins. After applying the method to a large-scale non-redundant database of proteins, we show that the positive signal of binding probability calculated from the local structural indicator is statistically meaningful. In summary, our studies suggested that the local environment of residues in a protein is a good indicator to recognize specific binding partners of the protein. The new method could be a potential addition to a suite of existing template-based approaches for protein function prediction.

A Computational Model for Kinetic Studies of Cadherin Binding and Clustering.

Cadherin is a cell-surface transmembrane receptor that mediates calcium-dependent cell-cell adhesion and is a major component of adhesive junctions. The formation of intercellular adhesive junctions is initiated by trans binding between cadherins on adjacent cells, which is followed by the clustering of cadherins via the formation of cis interactions between cadherins on the same cell membranes. Moreover, classical cadherins have multiple glycosylation sites along their extracellular regions. It was found that aberrant glycosylation affects the adhesive function of cadherins and correlates with metastatic phenotypes of several cancers. However, a mechanistic understanding of cadherin clustering during cell adhesion and the role of glycosylation in this process is still lacking. Here, we designed a kinetic model that includes multistep reaction pathways for cadherin clustering. We further applied a diffusion-reaction algorithm to numerically simulate the clustering process using a recently developed coarse-grained model. Using experimentally measured rates of trans binding between soluble E-cadherin extracellular domains, we conducted simulations of cadherin-mediated cell-cell binding kinetics, and the results are quantitatively comparable to experimental data from micropipette experiments. In addition, we show that incorporating cadherin clustering via cis interactions further increases intercellular binding. Interestingly, a two-phase kinetic profile was derived under the assumption that glycosylation regulates the kinetic rates of cis interactions. This two-phase profile is qualitatively consistent with experimental results from micropipette measurements. Therefore, our computational studies provide new, to our knowledge, insights into the molecular mechanism of cadherin-based cell adhesion.

Experimental and Theoretical Investigation of Effects of Ethanol and Acetic Acid on Carcinogenic NDMA Formation in Simulated Gastric Fluid.

N-nitrosodimethylamine (NDMA), as a representative of endogenously formed N-nitroso compounds (NOCs), has become the focus of considerable research interest due to its unusually high carcinogenicity. In this study, effects of ethanol and acetic acid on the formation of NDMA from dimethylamine (DMA) and nitrite in simulated gastric fluid (SGF) were investigated. Experimental results showed that ethanol in the concentrations of 1-8% (v/v) and acetic acid in the concentrations of 0.01-8% (v/v) exhibit inhibitory and promotion effects on the formation of NDMA, respectively. Moreover, they are both in a dose-dependent manner with the largest inhibition/promotion rate reaching ∼70%. Further experimental investigations indicate that ethanol and acetic acid are both able to scavenge nitrite in SGF. It implies that there are interactions of ethanol and acetic acid with nitrite or nitrite-related nitrosating agents rather than DMA. Theoretical calculations confirm the above experimental results and demonstrate that ethanol and acetic acid can both react with nitrite-related nitrosating agents to produce ethyl nitrite (EtONO) and acetyl nitrite (AcONO), respectively. Furthermore, the reactivities of ethyl nitrite, acetyl nitrite, and dinitrogen trioxide reacting with DMA were found in the order of AcONO > N2O3 ≫ EtONO. This is probably the main reason why there are completely different effects of ethanol and acetic acid on NDMA formation. On the basis of the above results, two requirements for a potential inhibitor of NOCs formation in SGF were provided. The results obtained in this study will be helpful in better understanding the inhibition/promotion mechanisms of compounds on NDMA formation in SGF and searching for protective substances to prevent carcinogenic NOCs formation.

Decomposing the space of protein quaternary structures with the interface fragment pair library.

BACKGROUND: The physical interactions between proteins constitute the basis of protein quaternary structures. They dominate many biological processes in living cells. Deciphering the structural features of interacting proteins is essential to understand their cellular functions. Similar to the space of protein tertiary structures in which discrete patterns are clearly observed on fold or sub-fold motif levels, it has been found that the space of protein quaternary structures is highly degenerate due to the packing of compact secondary structure elements at interfaces. Therefore, it is necessary to further decompose the protein quaternary structural space into a more local representation. RESULTS: Here we constructed an interface fragment pair library from the current structure database of protein complexes. After structural-based clustering, we found that more than 90% of these interface fragment pairs can be represented by a limited number of highly abundant motifs. These motifs were further used to guide complex assembly. A large-scale benchmark test shows that the native-like binding is highly likely in the structural ensemble of modeled protein complexes that were built through the library. CONCLUSIONS: Our study therefore presents supportive evidences that the space of protein quaternary structures can be represented by the combination of a small set of secondary-structure-based packing at binding interfaces. Finally, after future improvements such as adding sequence profiles, we expect this new library will be useful to predict structures of unknown protein-protein interactions.

LISE: a server using ligand-interacting and site-enriched protein triangles for prediction of ligand-binding sites.

LISE is a web server for a novel method for predicting small molecule binding sites on proteins. It differs from a number of servers currently available for such predictions in two aspects. First, rather than relying on knowledge of similar protein structures, identification of surface cavities or estimation of binding energy, LISE computes a score by counting geometric motifs extracted from sub-structures of interaction networks connecting protein and ligand atoms. These network motifs take into account spatial and physicochemical properties of ligand-interacting protein surface atoms. Second, LISE has now been more thoroughly tested, as, in addition to the evaluation we previously reported using two commonly used small benchmark test sets and targets of two community-based experiments on ligand-binding site predictions, we now report an evaluation using a large non-redundant data set containing >2000 protein-ligand complexes. This unprecedented test, the largest ever reported to our knowledge, demonstrates LISE's overall accuracy and robustness. Furthermore, we have identified some hard to predict protein classes and provided an estimate of the performance that can be expected from a state-of-the-art binding site prediction server, such as LISE, on a proteome scale. The server is freely available at http://lise.ibms.sinica.edu.tw.

A multiscale model for simulating binding kinetics of proteins with flexible linkers.

The kinetics of protein interactions are essential determinants in many cellular processes such as signal transduction and transcriptional regulation. Many proteins involved in these functions contain intrinsic disordered regions. This makes conformational flexibility become an unneglectable factor when studying the binding kinetic of these proteins. Compared with the binding of rigid proteins that is limited by diffusions, the binding mechanisms of proteins with internal flexibility are much more complicated. Using a small protein that contains two domains and a connecting loop as a testing system, we developed a multiscale simulation framework to study the role of flexible linkers in regulating kinetics of protein binding. The association and dissociation processes were implemented by a coarse-grained Monte-Carlo algorithm, while the conformational changes of the flexible linker were captured from all-atom molecular dynamic simulations. Our simulations illustrated that the presence of the extended domain linker can enhance the rate of protein association. On the other hand, the full-length flexible molecule is more difficult to dissociate than its two rigid domains but much easier than the molecule with a rigid linker. Overall, our studies demonstrated that both kinetics and thermodynamics of protein binding are closely modulated by the dynamic features of linker regions.

A coarse-grained model for the simulations of biomolecular interactions in cellular environments.

The interactions of bio-molecules constitute the key steps of cellular functions. However, in vivo binding properties differ significantly from their in vitro measurements due to the heterogeneity of cellular environments. Here we introduce a coarse-grained model based on rigid-body representation to study how factors such as cellular crowding and membrane confinement affect molecular binding. The macroscopic parameters such as the equilibrium constant and the kinetic rate constant are calibrated by adjusting the microscopic coefficients used in the numerical simulations. By changing these model parameters that are experimentally approachable, we are able to study the kinetic and thermodynamic properties of molecular binding, as well as the effects caused by specific cellular environments. We investigate the volumetric effects of crowded intracellular space on bio-molecular diffusion and diffusion-limited reactions. Furthermore, the binding constants of membrane proteins are currently difficult to measure. We provide quantitative estimations about how the binding of membrane proteins deviates from soluble proteins under different degrees of membrane confinements. The simulation results provide biological insights to the functions of membrane receptors on cell surfaces. Overall, our studies establish a connection between the details of molecular interactions and the heterogeneity of cellular environments.

Mitochondria-targeting peptide SS-31 attenuates ferroptosis via inhibition of the p38 MAPK signaling pathway in the hippocampus of epileptic rats.

Ferroptosis is a newly identified form of non-apoptotic regulated cell death (RCD) andplaysanimportantrole in epileptogenesis. The p38 mitogen-activated protein kinase (p38 MAPK) pathway has been confirmed to be involved in ferroptosis. The mitochondria-targeting antioxidant Elamipretide (SS-31) can reduce the generation of lipid peroxidation and the buildup of reactive oxygen species (ROS). Collectively, our present study was to decipher whether SS-31 inhibits ferroptosis via the p38 MAPK signaling pathway in the rat epilepsy model induced by pilocarpine (PILO).Adult male Wistar rats were randomly divided into four groups: control group (CON group), epilepsy group (EP group), SS-31 treatment group (SS group), and p38 MAPK inhibitor (SB203580) treatment group (SB group). Our results demonstrated that the rat hippocampal neurons after epilepsy were followed by accumulated iron and malondialdehyde (MDA) content, upregulated phosphorylated p38 MAPK protein (P-p38) and nuclear factor erythroid 2-related factor 2 (Nrf2) levels, reduced glutathione peroxidase 4 (Gpx4) content, and depleted glutathione (GSH) activity. Morphologically, mitochondrial ultrastructural damage under electron microscopy was manifested by a partial increase in outer membrane density, disappearance of mitochondrial cristae, and mitochondrial shrinkage. SS-31 and SB203580 treatment blocked the initiation and progression of ferroptosis in the hippocampus of epileptic rats via reducing the severity of epileptic seizures, reversing the expression of Gpx4, P-p38 , decreasing the levels of iron and MDA, as well as increasing the activity of GSH and Nrf2. To summarize, our findings proved that ferroptosis was coupled with the pathology of epilepsy, and SS-31 can inhibit PILO-induced seizures by preventing ferroptosis, which may be connected to the inhibition of p38 MAPK phosphorylation, highlighting the potential therapeutic value for targeting ferroptosis process in individuals with seizure-related diseases.

Ligand-binding site prediction using ligand-interacting and binding site-enriched protein triangles.

MOTIVATION: Knowledge about the site at which a ligand binds provides an important clue for predicting the function of a protein and is also often a prerequisite for performing docking computations in virtual drug design and screening. We have previously shown that certain ligand-interacting triangles of protein atoms, called protein triangles, tend to occur more frequently at ligand-binding sites than at other parts of the protein. RESULTS: In this work, we describe a new ligand-binding site prediction method that was developed based on binding site-enriched protein triangles. The new method was tested on 2 benchmark datasets and on 19 targets from two recent community-based studies of such predictions, and excellent results were obtained. Where comparisons were made, the success rates for the new method for the first predicted site were significantly better than methods that are not a meta-predictor. Further examination showed that, for most of the unsuccessful predictions, the pocket of the ligand-binding site was identified, but not the site itself, whereas for some others, the failure was not due to the method itself but due to the use of an incorrect biological unit in the structure examined, although using correct biological units would not necessarily improve the prediction success rates. These results suggest that the new method is a valuable new addition to a suite of existing structure-based bioinformatics tools for studies of molecular recognition and related functions of proteins in post-genomics research. AVAILABILITY: The executable binaries and a web server for our method are available from http://sourceforge.net/projects/msdock/ and http://lise.ibms.sinica.edu.tw, respectively, free for academic users.