Protein-ligand binding site recognition using complementary binding-specific substructure comparison and sequence profile alignment

The Identification of Metal Ion Ligand-Binding Residues by Adding the Reclassified Relative Solvent Accessibility

Many proteins realize their special functions by binding with specific metal ion ligands during a cell's life cycle. The ability to correctly identify metal ion ligand-binding residues is valuable for the human health and the design of molecular drug. Precisely identifying these residues, however, remains challenging work. We have presented an improved computational approach for predicting the binding residues of 10 metal ion ligands (Zn2+, Cu2+, Fe2+, Fe3+, Co2+, Ca2+, Mg2+, Mn2+, Na+, and K+) by adding reclassified relative solvent accessibility (RSA). The best accuracy of fivefold cross-validation was higher than 77.9%, which was about 16% higher than the previous result on the same dataset. It was found that different reclassification of the RSA information can make different contributions to the identification of specific ligand binding residues. Our study has provided an additional understanding of the effect of the RSA on the identification of metal ion ligand binding residues.

Filamentation of the bacterial bi-functional alcohol/aldehyde dehydrogenase AdhE is essential for substrate channeling and enzymatic regulation

Acetaldehyde-alcohol dehydrogenase (AdhE) enzymes are a key metabolic enzyme in bacterial physiology and pathogenicity. They convert acetyl-CoA to ethanol via an acetaldehyde intermediate during ethanol fermentation in an anaerobic environment. This two-step reaction is associated to NAD⁺ regeneration, essential for glycolysis. The bifunctional AdhE enzyme is conserved in all bacterial kingdoms but also in more phylogenetically distant microorganisms such as green microalgae. It is found as an oligomeric form called spirosomes, for which the function remains elusive. Here, we use cryo-electron microscopy to obtain structures of Escherichia coli spirosomes in different conformational states. We show that spirosomes contain active AdhE monomers, and that AdhE filamentation is essential for its activity in vitro and function in vivo. The detailed analysis of these structures provides insight showing that AdhE filamentation is essential for substrate channeling within the filament and for the regulation of enzyme activity.

Evaluation of potential drugs against leishmaniasis targeting catalytic subunit of Leishmania donovani nuclear DNA primase using ligand based virtual screening, docking and molecular dynamics approaches

Leishmania donovani, causes leishmaniasis, a global health trouble with around 89 different countries and its population under its risk. Replication initiation events have been instrumental in regulating the DNA duplication and as the small subunit of L. donovani nuclear DNA primase (Ld-PriS) inherits the catalytic site, it plays a vital role in DNA replication. In this study we have aimed Ld-PriS for the first time as a prospective target for the application of drug against Leishmania parasite. 3-D structures of Ld-PriS were built and ligand-based virtual screening was performed using hybrid similarity recognition techniques. Ligands from the ZINC database were used for the screening purposes based on known DNA primase inhibitor Sphingosine as a query. Top 150 ligands were taken into consideration for molecular docking against the query protein (Ld-PriS) using PyRx and iGEMDOCK softwares. Top five compounds with the best docking score were selected for pharmacokinetic investigation and molecular dynamic simulation. These top five screened inhibitors showed very poor binding affinity toward the catalytic subunit of human primase indicating their safety toward the host normal replication mechanism. The top five compounds showed good pharmacokinetic profiles and ADMET predictions revealed good absorption, solubility, permeability, uniform distribution, proper metabolism, minimal toxicity and good bioavailability. Simulation studies upto 50 ns revealed the three leads ZINC000009219046, ZINC000025998119 and ZINC000004677901 bind with Ld-PriS throughout the simulation and there were no huge variations in their backbone suggesting that these three may play as potential lead compounds for developing new drug against leishmaniasis.Communicated by Ramaswamy H. Sarma.

PfCERLI1 is a conserved rhoptry associated protein essential for Plasmodium falciparum merozoite invasion of erythrocytes

The disease-causing blood-stage of the Plasmodium falciparum lifecycle begins with invasion of human erythrocytes by merozoites. Many vaccine candidates with key roles in binding to the erythrocyte surface and entry are secreted from the large bulb-like rhoptry organelles at the apical tip of the merozoite. Here we identify an essential role for the conserved protein P. falciparum Cytosolically Exposed Rhoptry Leaflet Interacting protein 1 (PfCERLI1) in rhoptry function. We show that PfCERLI1 localises to the cytosolic face of the rhoptry bulb membrane and knockdown of PfCERLI1 inhibits merozoite invasion. While schizogony and merozoite organelle biogenesis appear normal, biochemical techniques and semi-quantitative super-resolution microscopy show that PfCERLI1 knockdown prevents secretion of key rhoptry antigens that coordinate merozoite invasion. PfCERLI1 is a rhoptry associated protein identified to have a direct role in function of this essential merozoite invasion organelle, which has broader implications for understanding apicomplexan invasion biology.

Rhoptries are essential organelles for invasion of erythrocytes by Plasmodium. Here, the authors characterize the rhoptry-associated protein CERLI1 using quantitative super-resolution microscopy, showing that it is important for parasite invasion and secretion of rhoptry proteins including vaccine antigens.

In Silico Modeling of Crimean Congo Hemorrhagic Fever Virus Glycoprotein-N and Screening of Anti Viral Hits by Virtual Screening

Crimean-Congo hemorrhagic fever (CCHF) is a widespread zoonotic viral disease, caused by a tick-born virus Crimean-Congo hemorrhagic fever virus (CCHFV). This disease is endemic in Middle East, Asia, Africa and South-Eastern Europe with the mortality rate of 5–30%. CCHFV genome is composed of three segments: large, medium and small segments. M segment encodes a polyprotein (glycoprotein) so called glycoprotein N (Gn) which is considered as a potential druggable target for the effective therapy of CCHF. The complete structure of Gn is still not characterized. The aim of the current study is to predict the complete three-dimensional (3D-) structure of CCHFV Gn protein via threading-based modeling and investigate the residues crucial for binding with CCHFV envelop. The developed model displayed excellent stereo-chemical and geometrical properties. Subsequently structure based virtual screening (SBVS) was applied to discover novel inhibitors of Gn protein. A library of > 1300 anti-virals was selected from PubChem database and directed to the predicted binding site of Gn. The SBVS results led to the identification of thirty-seven compounds that inhibit the protein in computational analysis. Those 37 hits were subject to pharmacokinetic profiling which demonstrated that 30/37 compound possess safer pharmacokinetic properties. Thus, by specifically targeting Gn, less toxic and more potent inhibitors of CCHFV were identified in silico.

In silico analysis of non-synonymous single nucleotide polymorphisms (nsSNPs) in the human GJA3 gene associated with congenital cataract

BACKGROUND

Gap junction protein alpha 3 (GJA3), an important pathogenic gene of congenital cataracts, encodes the transmembrane protein connexin46, which functions as an intercellular channel for voltage and chemical gating by forming dodecamers. This study systematically collected nsSNP information for the GJA3 gene from SNP databases and literature and screened for nsSNPs with high risks of pathogenicity.

RESULTS

A total of 379 nsSNPs of GJA3 were identified. A total of 88 high-risk pathogenic GJA3 nsSNPs were found, including 31 published nsSNPs associated with congenital cataracts and 57 novel nsSNPs predicted by all eight online tools. The 88 high-risk pathogenic mutations, which are related to 67 amino acids in the wild-type sequences, cause a decrease in protein stability according to I-Mutant 3.0, MUpro and INPS. G2 and R33 were predicted to participate in post-translational modification and ligand binding by ModPred, RaptorX Binding and COACH. Additionally, high-risk mutations were likely to involve highly conserved sites, random coils, alpha helixes, and extracellular loops and were accompanied by changes in amino acid size, charge, hydrophobicity and spatial structure.

CONCLUSIONS

Eighty-eight high-risk pathogenic nsSNPs of GJA3 were screened out in the study, 57 of which were newly reported. The combination of multiple in silico tools is highly efficient for targeting pathogenic sites.

Genome-wide analysis of the HSP101/CLPB gene family for heat tolerance in hexaploid wheat

Heat Shock Protein 101 (HSP101), the homolog of Caseinolytic Protease B (CLPB) proteins, has functional conservation across species to play roles in heat acclimation and plant development. In wheat, several TaHSP101/CLPB genes were identified, but have not been comprehensively characterized. Given the complexity of a polyploid genome with its phenomena of homoeologous expression bias, detailed analysis on the whole TaCLPB family members is important to understand the genetic basis of heat tolerance in hexaploid wheat. In this study, a genome-wide analysis revealed thirteen members of TaCLPB gene family and their expression patterns in various tissues, developmental stages, and stress conditions. Detailed characterization of TaCLPB gene and protein structures suggested potential variations of the sub-cellular localization and their functional regulations. We revealed homoeologous specific variations among TaCLPB gene copies that have not been reported earlier. A study of the Chromosome 1 TaCLPB in four wheat genotypes demonstrated unique patterns of the homoeologous gene expression under moderate and extreme heat treatments. The results give insight into the strategies to improve heat tolerance by targeting one or some of the TaCLPB genes in wheat.

Exploring the computational methods for protein-ligand binding site prediction

Proteins participate in various essential processes in vivo via interactions with other molecules. Identifying the residues participating in these interactions not only provides biological insights for protein function studies but also has great significance for drug discoveries. Therefore, predicting protein-ligand binding sites has long been under intense research in the fields of bioinformatics and computer aided drug discovery. In this review, we first introduce the research background of predicting protein-ligand binding sites and then classify the methods into four categories, namely, 3D structure-based, template similarity-based, traditional machine learning-based and deep learning-based methods. We describe representative algorithms in each category and elaborate on machine learning and deep learning-based prediction methods in more detail. Finally, we discuss the trends and challenges of the current research such as molecular dynamics simulation based cryptic binding sites prediction, and highlight prospective directions for the near future.

In silico identification and evaluation of Bacillus subtilis cold shock protein B (cspB)-like plant RNA chaperones

Cold shock domain (CSD) proteins with nucleic acid binding properties are well conserved from bacteria to higher organisms. In bacteria, the cold shock proteins (CSPs) are single domain RNA chaperones, whereas in animals and plants, CSDs are accompanied by additional domains with roles in transcription regulation. Bacterial CSPs (Escherischia coli-cspA and Bacilus subtilis-cspB) have successfully imparted drought tolerance in transgenic plants; however, these cannot be deployed in food crops due to their low public acceptance of transgenics with bacterial genes. Therefore, this study aimed to identify CSPB-like proteins from plants that can be used for developing drought tolerant transgenic crops. Twelve single domain plant CSPs presenting >40% sequence identity with CSPB were identified. All 12 plant CSPs were modeled by homology modeling and refined by molecular dynamics simulation for 10 ns. Selected plant CSPs and CSPB exhibited high structural similarity (Tm-score: 0.63-0.86). Structure based phylogenetic analysis revealed that Triticum aestivum-csp1 and Aegilops tauschii-cspE are structurally closer to CSPB compared to their orthologs and paralogs. Molecular docking with three RNA molecules (5U, UC3U, and C2UC) indicates that Ricinus communis-csd1 and T. aestivum-csp1 have a binding pattern and docking scores similar to those of CSPB. Furthermore, MD simulations for 20 ns and analysis of RMSD, RMSF, Rg as well as the number of hydrogen bonds in all the three complexes revealed that plant CSP-RNA complexes behave in a similar manner to that of the CSPB-RNA complex, making them highly potential candidate genes for developing drought tolerance in transgenic plants. Communicated by Ramaswamy H. Sarma.

Investigation of the interaction of DAD1-LIKE LIPASE 3 (DALL3) with Selenium Binding Protein 1 (SBP1) in Arabidopsis thaliana

Phospholipase PLA₁-Iγ2 or otherwise DAD1-LIKE LIPASE 3 (DALL3) is a member of class I phospholipases and has a role in JA biosynthesis. AtDALL3 was previously identified in a yeast two-hybrid screening as an interacting protein of the Arabidopsis Selenium Binding Protein 1 (SBP1). In this work, we have studied AtDALL3 as an interacting partner of the Arabidopsis Selenium Binding Protein 1 (SBP1). Phylogenetic analysis showed that DALL3 appears in the PLA1-Igamma1, 2 group, paired with PLA1-Igammma1. The highest level of expression of AtDALL3 was observed in 10-day-old roots and in flowers, while constitutive levels were maintained in seedlings, cotyledons, shoots and leaves. In response to abiotic stress, DALL3 was shown to participate in the network of genes regulated by cadmium, selenite and selenate compounds. DALL3 promoter driven GUS assays revealed that the expression patterns defined were overlapping with the patterns reported for AtSBP1 gene, indicating that DALL3 and SBP1 transcripts co-localize. Furthermore, quantitative GUS assays showed that these compounds elicited changes in activity in specific cells files, indicating the differential response of DALL3 promoter. GFP::DALL3 studies by confocal microscopy demonstrated the localization of DALL3 in the plastids of the root apex, the plastids of the central root and the apex of emerging lateral root primordia. Additionally, we confirmed by yeast two hybrid assays the physical interaction of DALL3 with SBP1 and defined a minimal SBP1 fragment that DALL3 binds to. Finally, by employing bimolecular fluorescent complementation we demonstrated the in planta interaction of the two proteins.

In silico analysis of alternative splicing on drug-target gene interactions

Identifying and evaluating the right target are the most important factors in early drug discovery phase. Most studies focus on one protein ignoring the multiple splice-variant or protein-isoforms, which might contribute to unexpected therapeutic activity or adverse side effects. Here, we present computational analysis of cancer drug-target interactions affected by alternative splicing. By integrating information from publicly available databases, we curated 883 FDA approved or investigational stage small molecule cancer drugs that target 1,434 different genes, with an average of 5.22 protein isoforms per gene. Of these, 618 genes have ≥5 annotated protein-isoforms. By analyzing the interactions with binding pocket information, we found that 76% of drugs either miss a potential target isoform or target other isoforms with varied expression in multiple normal tissues. We present sequence and structure level alignments at isoform-level and make this information publicly available for all the curated drugs. Structure-level analysis showed ligand binding pocket architectures differences in size, shape and electrostatic parameters between isoforms. Our results emphasize how potentially important isoform-level interactions could be missed by solely focusing on the canonical isoform, and suggest that on- and off-target effects at isoform-level should be investigated to enhance the productivity of drug-discovery research.

Kaposi's Sarcoma-Associated Herpesvirus ORF66 Is Essential for Late Gene Expression and Virus Production via Interaction with ORF34

Kaposi's sarcoma-associated herpesvirus (KSHV) is closely associated with B-cell and endothelial cell malignancies. After the initial infection, KSHV retains its viral genome in the nucleus of the host cell and establishes a lifelong latency. During lytic infection, KSHV-encoded lytic-related proteins are expressed in a sequential manner and are classified as immediate early, early, and late (L) gene transcripts. The transcriptional initiation of KSHV late genes is thought to require the complex formation of the viral preinitiation complex (vPIC), which may consist of at least 6 transcription factors (ORF18, -24, -30, -31, -34, and -66). However, the functional role of ORF66 in vPIC during KSHV replication remains largely unclear. Here, we generated ORF66-deficient KSHV using a bacterial artificial chromosome (BAC) system to evaluate its role during viral replication. While ORF66-deficient KSHV demonstrated mainly attenuated late gene expression and decreased virus production, viral DNA replication was unaffected. Chromatin immunoprecipitation analysis showed that ORF66 bound to the promoters of a late gene (K8.1) but did not bind to those of a latent gene (ORF72), an immediate early gene (ORF16), or an early gene (ORF46/47). Furthermore, we found that three highly conserved C-X-X-C sequences and a conserved leucine repeat in the C-terminal region of ORF66 were essential for the interaction with ORF34, the transcription of K8.1, and virus production. The interaction between ORF66 and ORF34 occurred in a zinc-dependent manner. Our data support a model in which ORF66 serves as a critical vPIC component to promote late viral gene expression and virus production.IMPORTANCE KSHV ORF66 is expressed during the early stages of lytic infection, and ORF66 and vPIC are thought to contribute significantly to late gene expression. However, the physiological importance of ORF66 in terms of vPIC formation remains poorly understood. Therefore, we generated an ORF66-deficient BAC clone and evaluated its viral replication. The results showed that ORF66 plays a critical role in virus production and the transcription of L genes. To our knowledge, this is the first report showing the function of ORF66 in virus replication using ORF66-deficient KSHV. We also clarified that ORF66 interacts with the transcription start site of the K8.1 gene, a late gene. Furthermore, we identified the ORF34-binding motifs in the ORF66 C terminus: three C-X-X-C sequences and a leucine-repeat sequence, which are highly conserved among beta- and gammaherpesviruses. Our study provides insights into the regulatory mechanisms of not only the late gene expression of KSHV but also those of other herpesviruses.

Computational prediction of active sites and ligands in different AHL quorum quenching lactonases and acylases

With the emergence of multidrug-resistant 'superbug', conventional treatments become obsolete. Quorum quenching (QQ), enzyme-dependent alteration of quorum sensing (QS), is now considered as a promising antimicrobial therapy because of its potentiality to impede virulence gene expression without resulting in growth inhibition and antibiotic resistance. In our study, we intended to compare between two major QQ enzyme groups (i.e., AHL lactonases and AHL acylases) in terms of their structural and functional aspects. The amino acid composition-based principal component analysis (PCA) suggested that probably there is no structural and functional overlapping between the two groups of enzymes as well as within the lactonase enzymes but the acylases may functionally be affected by one another. In subcellular localization analysis, we also found that most lactonases are cytoplasmic while acylases are periplasmic. Investigation on the secondary structural features showed random coil dominates over alpha-helix and beta-sheet in all evaluated enzymes. For structural comparison, the tertiary structures of the selected proteins were modelled and submitted to the PMDB database (Accession ID: PM0081007 to PM0081018). Interestingly, sequence alignment revealed the presence of several conserved domains important for functions in both protein groups. In addition, three amino acid residues, namely aspartic acid, histidine, and isoleucine, were common in the active sites of all protein models while most frequent ligands were found to be 3C7, FEO, and PAC. Importantly, binding interactions of predicted ligands were similar to that of native QS signal molecules. Furthermore, hydrogen bonds analysis suggested six proteins are more stable than others. We believe that the knowledge of this comparative study could be useful for further research in the development of QSbased universal antibacterial strategies.

Germline Mutation in KIF1Bβ Gene Associated with Loss of Heterozygosity: Usefulness of Next-Generation Sequencing in the Genetic Screening of Patients with Pheochromocytoma

The genetic approach of pheochromocytomas and paragangliomas has changed in the last two decades. Nowadays, we know that more than 40% of patients have a germline mutation in one of the susceptibility genes identified to date. Our aim is to underline how genetic diagnosis by next-generation sequencing (NGS) can improve the management of patients affected by pheochromocytomas and paragangliomas in our routine diagnostic screening. We reported a case presentation and next-generation sequencing analysis supported by in silico studies and evaluation of mitochondrial status in KIF1Bβ tissue. A 46-year-old male affected by a left secreting pheochromocytoma underwent surgery in 2017. After surgery, the normetanephrine levels decreased very slowly and a suspected abdominal lymph node was detected. We found a novel germline KIF1Bβ gene mutation, c.4052C > T, p. Pro1351Leu associated with tumor loss of heterozygosity, and resulted likely-pathogenetic by in silico studies. This mutation was also associated with an increased number of mitochondria through the electron microscopy compared with wild-type tissues as suggestive for mitochondria neoformation compensatory to the mitochondrial autophagic figures observed. Our results underline the usefulness of next-generation sequencing in the presence of multiple tumor predisposition genes and how, at the same time, its use may result challenging for the clinicians. To date, performing the genetic analysis according to the latest Consensus Statement is mandatory in patients affected by PHEO/PGL.

Oxidative stress induced by cadmium in lettuce (Lactuca sativa Linn.): Oxidative stress indicators and prediction of their genes

Environmental contamination with heavy metals is of concern as plants have the ability to absorb chemical toxicants facilitating the entry of toxic metals into the food chain. Lettuce (Lactuca sativa Linn.) was cultured in four nutrient solutions containing different concentrations of cadmium (0, 3, 6, and 9 mmol). The impact of heavy metal on the morphological features, antioxidant properties and antioxidant enzymes activity were investigated with primary focus on superoxide dismutase, ascorbate peroxidase, peroxidase and catalase enzymes. In silico methods were utilized in the study of the genes of these enzymes. Significant changes were observed in the morphological features of the plant with plants appearing stunted, more spherical and yellow in colour. A decrease in the dry mass of the plant was also detected. The Translocation factor (TF) for cadmium was significantly high in lettuce. Enhanced antioxidant enzymatic activity suggests that these enzymes are integrally involved in the defense mechanism of the plant to heavy metal stress. Also observed was an increase in total soluble protein, and total phenolic content. Total flavonoid content was not significantly affected. Fourteen genes encoding for ascorbate peroxidase and nineteen genes for superoxide dismutase were identified in lettuce. These enzymes varied from each other with regards to the number of exons and amino acids present, as well as their location within the cell. Plants exhibit various response mechanisms to combat heavy metal contamination.

Development of multi-epitope chimeric vaccine against Taenia solium by exploring its proteome: an in silico approach

Objective: Taenia solium is a neglected tropical disease; larvae of this parasite infect central nervous system i.e. Neurocysticercosis, and adults mature and survive into intestine i.e. Taeniasis. Globally more than 50 million people are at the risk of infection. This is one of the main etiological agents for onset of new early epilepsy in developing countries. However, there is no vaccine available to protect human from its infection. Hence, there is an urgent need for a good vaccine.Methods: We applied immune-informatics approach to design a multi-epitope chimeric vaccine consisting of both B and T-cell epitopes.Results: From the whole transcriptome of Taenia, we identified five suitable peptides present on cell membrane, epitope identification on these peptides were done by using various immunoinformatic software. Physiochemical properties were determined and the tertiary structure of vaccine was predicted, validated and refined, and to increase antigenicity we added linker to them. Best-modeled protein-complex was used for docking study with TLR1-2, TLR4, TLR3 and TLR7 and stability of molecular complex was determined by molecular dynamics simulation.Conclusions: Overall, we attempted to design an efficient subunit chimeric vaccine, which could stimulate humoral and cellular immune responses and could protect against both neurocysticercosis and taeniasis.

CB-Dock: a web server for cavity detection-guided protein-ligand blind docking

As the number of elucidated protein structures is rapidly increasing, the growing data call for methods to efficiently exploit the structural information for biological and pharmaceutical purposes. Given the three-dimensional (3D) structure of a protein and a ligand, predicting their binding sites and affinity are a key task for computer-aided drug discovery. To address this task, a variety of docking tools have been developed. Most of them focus on docking in the preset binding sites given by users. To automatically predict binding modes without information about binding sites, we developed a user-friendly blind docking web server, named CB-Dock, which predicts binding sites of a given protein and calculates the centers and sizes with a novel curvature-based cavity detection approach, and performs docking with a popular docking program, Autodock Vina. This method was carefully optimized and achieved ~70% success rate for the top-ranking poses whose root mean square deviation (RMSD) were within 2 Å from the X-ray pose, which outperformed the state-of-the-art blind docking tools in our benchmark tests. CB-Dock offers an interactive 3D visualization of results, and is freely available at http://cao.labshare.cn/cb-dock/.

Haplotype-Based Association and In Silico Studies of OPRM1 Gene Variants with Susceptibility to Opioid Dependence Among Addicted Iranians Undergoing Methadone Treatment

The associations of OPRM1 gene variants with opioid dependence have been demonstrated. This study investigated the association of rs495491, rs1799971 (A118G), rs589046, and rs10457090 variants of OPRM1 gene with opium dependence and their haplotypes among addicted individuals undergoing methadone treatment. Moreover, we investigated whether any of these variants were associated with libido dysfunction or insomnia among addicted people. A total of 404 individuals were genotyped by amplification refractory mutation system (ARMS) PCR. In silico studies were designed through homology modeling of A118G structures (N40 and D40) and docked with 41 FDA-approved drugs of OPRM1 protein by SWISS-MODEL, COACH, MolProbity, ProSA, Errat, Glide XP, and Autodock 4. Results revealed that rs495491, A118G, rs589046, and rs10457090 were significantly associated with opium dependence under recessive (P = 6.66E-10), dominant (P = 0.017), co-dominant (P = 0.001), and recessive (P = 9.28E-6) models of inheritance, respectively. Further analyses indicated three significant haplotypes including A-A-A-C (P-permutation < 1E-9), G-G-A-C (P-permutation = 0.04), and G-A-G-C (P-permutation = 8.69E-4). Genotype-phenotype associations of OPRM1 variants with insomnia and libido dysfunction showed no significant association. Docking showed the higher binding affinity of N40 rather than D40 model; however, methadone and morphine were bonded with D40 structure more powerful. Consequently, rs495491, A118G, rs589046, and rs10457090 were associated with opioid dependence among Iranians; also, A118G might be the most remarkable marker of OPRM1 owing to its vital structural roles.

Communication within East Antarctic Soil Bacteria

Antarctica, being the coldest, driest, and windiest continent on Earth, represents the most extreme environment in which a living organism can survive. Under constant exposure to harsh environmental threats, terrestrial Antarctica remains home to a great diversity of microorganisms, indicating that the soil bacteria must have adapted a range of survival strategies that require cell-to-cell communication. Survival strategies include secondary metabolite production, biofilm formation, bioluminescence, symbiosis, conjugation, sporulation, and motility, all of which are often regulated by quorum sensing (QS), a type of bacterial communication. Until now, such mechanisms have not been explored in terrestrial Antarctica. In this study, LuxI/LuxR-based quorum sensing (QS) activity was delineated in soil bacterial isolates recovered from Adams Flat, in the Vestfold Hills region of East Antarctica. Interestingly, we identified the production of potential homoserine lactones (HSLs) with chain lengths ranging from medium to long in 19 bacterial species using three biosensors, namely, Agrobacterium tumefaciens NTL4, Chromobacterium violaceum CV026, and Escherichia coli MT102, in conjunction with thin-layer chromatography (TLC). The majority of detectable HSLs were from Gram-positive species not previously known to produce HSLs. This discovery further expands our understanding of the microbial community capable of this type of communication, as well as provides insights into physiological adaptations of microorganisms that allow them to survive in the harsh Antarctic environment.IMPORTANCE Quorum sensing, a type of bacterial communication, is widely known to regulate many processes, including those that confer a survival advantage. However, little is known about communication by bacteria residing within Antarctic soils. Employing a combination of bacterial biosensors, analytical techniques, and genome mining, we found a variety of Antarctic soil bacteria speaking a common language, via LuxI/LuxR-based quorum sensing, thus potentially supporting survival in a mixed microbial community. This study reports potential quorum sensing activity in Antarctic soils and has provided a platform for studying physiological adaptations of microorganisms that allow them to survive in the harsh Antarctic environment.

Deleterious amino acid substitutions with a series of putative damaging effects on egg components are revealed in the ovalbumin gene family; an in silico approach

This study was conducted to identify the most deleterious nonsynonymous single nucleotide polymorphisms (nsSNPs) in the ovalbumin gene family, including OVALX, OVALY, and OVAL genes, which are involved in the synthesis of the most important components in the chickens’ eggs using a comprehensive in silico approach. Ten different computational servers were utilized to prioritize the possible deleterious effects of the retrieved nsSNPs in terms of structure, function, and stability. Results indicated entirely damaging effects of H365P in OVALX, I167T in OVALY, and V209G, L231P, F307C, and S317P in OVAL proteins. Further prediction tools showed that all of these deleterious nsSNPs were positioned in variable locations within several α-helix motifs in all studied ovalbumin proteins. Furthermore, all witnessed nsSNPs were predicted to be resided in the receptors binding sites, signifying remarkable involvement of such nsSNPs in damaging of the altered proteins. In conclusion, the present study provides the first inclusive data with regard to the most deleterious nsSNPs in OVALX, OVALY and OVAL genes in chickens. The present bioinformatics data may be useful for breeders who intend to raise chickens for egg production, in such a way the presence of any of these deleterious nsSNPs in any selected breed may possess several damaging effects on the egg components, which may impair egg production. Therefore, it can be stated that breeders have to confirm the absence of any of these deleterious nsSNPs before being proceeded further for large-scale egg-production purposes.

Moonlighting of Haemophilus influenzae heme acquisition systems contributes to the host airway-pathogen interplay in a coordinated manner

Nutrient iron sequestration is the most significant form of nutritional immunity and causes bacterial pathogens to evolve strategies of host iron scavenging. Cigarette smoking contains iron particulates altering lung and systemic iron homeostasis, which may enhance colonization in the lungs of patients suffering chronic obstructive pulmonary disease (COPD) by opportunistic pathogens such as nontypeable. NTHi is a heme auxotroph, and the NTHi genome contains multiple heme acquisition systems whose role in pulmonary infection requires a global understanding. In this study, we determined the relative contribution to NTHi airway infection of the four heme-acquisition systems HxuCBA, PE, SapABCDFZ, and HbpA-DppBCDF that are located at the bacterial outer membrane or the periplasm. Our computational studies provided plausible 3D models for HbpA, SapA, PE, and HxuA interactions with heme. Generation and characterization of single mutants in the hxuCBA, hpe, sapA, and hbpA genes provided evidence for participation in heme binding-storage and inter-bacterial donation. The hxuA, sapA, hbpA, and hpe genes showed differential expression and responded to heme. Moreover, HxuCBA, PE, SapABCDFZ, and HbpA-DppBCDF presented moonlighting properties related to resistance to antimicrobial peptides or glutathione import, together likely contributing to the NTHi-host airway interplay, as observed upon cultured airway epithelia and in vivo lung infection. The observed multi-functionality was shown to be system-specific, thus limiting redundancy. Together, we provide evidence for heme uptake systems as bacterial factors that act in a coordinated and multi-functional manner to subvert nutritional- and other sources of host innate immunity during NTHi airway infection.

Binding site matching in rational drug design: algorithms and applications

Interactions between proteins and small molecules are critical for biological functions. These interactions often occur in small cavities within protein structures, known as ligand-binding pockets. Understanding the physicochemical qualities of binding pockets is essential to improve not only our basic knowledge of biological systems, but also drug development procedures. In order to quantify similarities among pockets in terms of their geometries and chemical properties, either bound ligands can be compared to one another or binding sites can be matched directly. Both perspectives routinely take advantage of computational methods including various techniques to represent and compare small molecules as well as local protein structures. In this review, we survey 12 tools widely used to match pockets. These methods are divided into five categories based on the algorithm implemented to construct binding-site alignments. In addition to the comprehensive analysis of their algorithms, test sets and the performance of each method are described. We also discuss general pharmacological applications of computational pocket matching in drug repurposing, polypharmacology and side effects. Reflecting on the importance of these techniques in drug discovery, in the end, we elaborate on the development of more accurate meta-predictors, the incorporation of protein flexibility and the integration of powerful artificial intelligence technologies such as deep learning.

Identification of a metallothionein gene in honey bee Apis mellifera and its expression profile in response to Cd, Cu and Pb exposure

Metallothioneins are ubiquitous proteins important in metal homeostasis and detoxification. However, they have not previously been identified in honey bees or other Hymenoptera, where metallothioneins could be of ecophysiological and ecotoxicological significance. Better understanding of the molecular responses to stress induced by toxic metals could contribute to honey bee conservation. In addition, honey bee metallothionein could represent a biomarker for monitoring environmental quality. Here we identify and characterize a metallothionein gene in Apis mellifera (AmMT). AmMT is 1,680 bp long and encodes a 48 amino acids protein with 15 cysteines and no aromatic residues. A metal response element upstream of the start codon, coupled with numerous cis-regulatory elements indicate the functional context of AmMT. Molecular modelling predicts several transition metal binding sites, and comparative phylogenetic analysis revealed five putative metallothionein proteins in three other hymenoptera species. AmMT was characterized by cloning the full-length coding sequence of the putative metallothionein. Recombinant AmMT was found to increase metal tolerance upon overexpression in Escherichia coli supplemented with Cd, Cu or Pb. Finally, in laboratory tests on honey bees, gene expression profiles showed a dose-dependant relationship between Cd, Cu and Pb concentrations present in food and AmMT expression, while field experiments showed induction of AmMT in bees from an industrial site compared to those from an urban area. These studies suggest that AmMT has metal binding properties in agreement with a possible role in metal homeostasis. Further functional and structural characterization of metallothionein in honey bees and other Hymenoptera are necessary.

SXGBsite: Prediction of Protein-Ligand Binding Sites Using Sequence Information and Extreme Gradient Boosting

The prediction of protein-ligand binding sites is important in drug discovery and drug design. Protein-ligand binding site prediction computational methods are inexpensive and fast compared with experimental methods. This paper proposes a new computational method, SXGBsite, which includes the synthetic minority over-sampling technique (SMOTE) and the Extreme Gradient Boosting (XGBoost). SXGBsite uses the position-specific scoring matrix discrete cosine transform (PSSM-DCT) and predicted solvent accessibility (PSA) to extract features containing sequence information. A new balanced dataset was generated by SMOTE to improve classifier performance, and a prediction model was constructed using XGBoost. The parallel computing and regularization techniques enabled high-quality and fast predictions and mitigated overfitting caused by SMOTE. An evaluation using 12 different types of ligand binding site independent test sets showed that SXGBsite performs similarly to the existing methods on eight of the independent test sets with a faster computation time. SXGBsite may be applied as a complement to biological experiments.

Intelligent Design of 14-3-3 Docking Proteins Utilizing Synthetic Evolution Artificial Intelligence (SYN-AI)

The ability to write DNA code from scratch will allow for the discovery of new and interesting chemistries as well as allowing the rewiring of cell signal pathways. Herein, we have utilized synthetic evolution artificial intelligence (SYN-AI) to intelligently design a set of 14-3-3 docking genes. SYN-AI engineers synthetic genes utilizing a parental gene as an evolution template. Wherein, evolution is fast-forwarded by transforming template gene sequences to DNA secondary and tertiary codes based upon gene hierarchical structural levels. The DNA secondary code allows identification of genomic building blocks across an orthologous sequence space comprising multiple genomes. Where, the DNA tertiary code allows engineering of supersecondary structures. SYN-AI constructed a library of 10 million genes that was reduced to three structurally functional 14-3-3 docking genes by applying natural selection protocols. Synthetic protein identity was verified utilizing Clustal Omega sequence alignments and Phylogeny.fr phylogenetic analysis. Wherein, we were able to confirm the three-dimensional structure utilizing I-TASSER and protein-ligand interactions utilizing COACH and Cofactor. The conservation of allosteric communications was confirmed utilizing elastic and anisotropic network models. Whereby, we utilized elNemo and ANM2.1 to confirm conservation of the 14-3-3 ζ amphipathic groove. Notably, to the best of our knowledge, we report the first 14-3-3 docking genes to be written from scratch.

Understanding the Effect of Structural Diversity in WRKY Transcription Factors on DNA Binding Efficiency through Molecular Dynamics Simulation

A precise understanding of the molecular mechanism involved in stress conditions has great importance for crop improvement. Biomolecules, such as WRKY proteins, which are the largest transcription factor family that is widely distributed in higher plants, plays a significant role in plant defense response against various biotic and abiotic stressors. In the present study, an extensive homology-based three-dimensional model construction and subsequent interaction study of WRKY DNA-binding domain (DBD) in CcWRKY1 (Type I), CcWRKY51 (Type II), and CcWRKY70 (Type III) belonging to pigeonpea, a highly tolerant crop species, was performed. Evaluation of the generated protein models was done to check their reliability and accuracy based on the quantitative and qualitative parameters. The final model was subjected to investigate the comparative binding analysis of different types of WRKY–DBD with DNA-W-box (a cis-acting element) by protein–DNA docking and molecular dynamics (MD) simulation. The DNA binding specificity with WRKY variants was scrutinized through protein–DNA interaction using the HADDOCK server. The stability, as well as conformational changes of protein–DNA complex, was investigated through molecular dynamics (MD) simulations for 100 ns using GROMACS. Additionally, the comparative stability and dynamic behavior of each residue of the WRKY–DBD type were analyzed in terms of root mean square deviation (RMSD), root mean square fluctuation (RMSF)values of the backbone atoms for each frame taking the minimized structure as a reference. The details of DNA binding activity of three different types of WRKY–DBD provided here will be helpful to better understand the regulation of WRKY gene family members in plants.

Structure-function relationships of the 5-oxoprolinase subunit A: Guiding biological sciences students down the path less traveled

Bioinformatics was recently introduced as a module for both undergraduate and postgraduate biological sciences students at our institution. Our experience shows that inquiry-based hands-on exercises provide the most efficient approach to bioinformatic straining. In this article, we report a structural bioinformatics project carried out by Master degree students to determine structure-function relationships of the uncharacterized prokaryotic 5-oxoprolinase subunit A (PxpA). PxpA associates with the PxpBC complex to form a functional 5-oxoprolinase enzyme for conversion of 5-oxoproline to L-glutamate. Although the exact role of PxpA is yet to be determined, it has been demonstrated that PxpBC catalyses the first step of the reaction, which is phosphorylation of 5-oxoproline. Here, we provide evidence that PxpA is involved in the last two steps of the reaction:decyclization of the labile phosphorylated 5-oxoproline to the equally labile γ-glutamylphosphate, and subsequent dephosphorylation to L-glutamate. Structural bioinformatics analysis of four putative PxpA structures revealed that PxpA adopts a non-canonical TIM barrel fold with well-characterized TIM barrel enzyme features. These include a C-terminal groove comprising potentially essential conserved amino acid residues organized into putative motifs. Phylogenetic analysis suggests a relationship between taxonomic grouping and PxpA oligomerization. PxpA forms a tunnel upon ligand binding, thus suggesting that the PxpABC complex employs the mechanism of substrate channeling to protect labile intermediates. Ultimately, students were able to form a testable hypothesis on the function of PxpA, an achievement we consider encouraging other students to emulate. © 2019 International Union of Biochemistry and Molecular Biology, 47(6):620-631, 2019.

Molecular docking study of lamellarin analogues and identification of potential inhibitors of HIV-1 integrase strand transfer complex by virtual screening

Molecular docking has been applied to elucidate the binding of lamellarin analogues with HIV-1 integrase strand transfer complex (PDB ID: 5U1C). The results suggest hydrogen bond interaction with residue Glu92 is key, and stabilisation by π-π stacking interactions with DNA base is chiefly influential to strand transfer activity. Other residues involved in hydrogen bonding are Cys65, His67, Asp64, Asp116 and chelation with Mg2+ ion was seen for certain analogues. Furthermore, hydrophobic interactions can be accounted for several amino acids including Asp64, Cys65, Asp116, His67, Glu92, Tyr143, Phe121, Gly118, Pro142 and Val72, as well as the DNA base. The molecular docking results are in line with the reported literatures of other inhibitors and strand transfer activity observed previously by Faulkner. We further employed molecular docking simulation to virtually screen and identified 4 novel potential inhibitors of HIV-1 integrase strand transfer complex from a Chembridge diversity collection of 25,132 small molecule compounds; Chembridge ID compound codes: 22850303, 27553460, 24578440 and 27591056. The candidates clearly formed hydrogen bonding interactions with important residues: His67 and Glu92. In addition, hydrophobic interactions were seen with residues similar to interactions with lamellarin analogues. The calculated drug-like scores are suggestive of these compounds to have clinical potential and ADMET predictions implied of their acceptable pharmacokinetic and toxicity profiles.

Evolutionary significance and functional characterization of streptomycin adenylyltransferase from Serratia marcescens

Complete functional annotations of proteins are essential to understand the role and mechanisms in pathogenesis. Aminoglycoside nucleotidyltransferases are the subclasses of aminoglycosides modifying enzymes conferring resistance to organisms. Insight into the structural and functional understanding of nucleotidyltransferase family protein provides vital information to combat pathogenesis. Phylogenetic analysis is employed to identify the evolutionary significance and common motif's present among the homologs of nucleotidyltransferase family protein. Structure, sequence based approaches and molecular docking were implemented to predict the exact function of the protein. Wide distribution of the nucleotidyltransferase family protein in gram-positive and gram-negative organisms are evidenced from phylogenetic analysis. Five common motifs were present in all the homolog's of nucleotidyltransferase family protein. Sequence-structure based functional annotations predicts that the targeted protein function as ATP-Mg dependent streptomycin adenylyltransferase. Structural comparisons and docking studies correlate well with the identified function. The complete function of nucleotidyltransferase family protein was identified as Streptomycin adenylyltransferase and it could be targeted as a potential therapeutic target to overcome antibiotic resistance.Communicated by Ramaswamy H. SarmaAbbreviationsAACaminoglycoside acetyltransferasesAMEaminoglycoside modifying enzymeANTaminoglycoside nucleotidyltransferasesAPHaminoglycoside phosphotransferasesATPadenosine triphosphateCASTpcomputer atlas and surface topography of proteinsDUFdomains of unknown functionGlidegrid-based ligand docking with energeticHMMhidden Markov modelMASTmotif alignment and search toolMEGAmolecular evolutionary genetics analysisMEMEmultiple Em for motif elicitationMSAmultiple sequence alignmentNMPnucleoside monophosphateNTPnucleoside triphosphateNTnucleotidyltransferaseOPLSoptimized potential for liquid simulationXPextra precision.

Uncovering the structure and function of Pseudomonas aeruginosa periplasmic proteins by an in silico approach

Pseudomonas aeruginosa is an opportunistic human pathogen highly relevant from a biomedical viewpoint. It is one of the main causes of infection in hospitalized patients and a major cause of mortality of cystic fibrosis patients. This is also due to its ability to develop resistance to antibiotics by various mechanisms. Therefore, it is urgent and desirable to identify novel targets for the development of new antibacterial drugs against Pseudomonas aeruginosa. In this work this problem was tackled by an in silico approach aimed at providing a reliable structural model and functional annotation for the Pseudomonas aeruginosa periplasmic proteins for which these data are not available yet. A total of 83 protein sequences were analyzed, and the corresponding structural models were built, leading to the identification of 32 periplasmic 'substrate-binding proteins', 14 enzymes and 4 proteins with different functions, including lipids and metals binding. The most interesting cases were found within the 'enzymes' group with the identification of a lipase, which can be regarded as a virulence factor, a protease involved in the assembly of β-barrel membrane proteins and a l,d-transpeptidase, which could contribute to confer resistance to β-lactam antibiotics to the bacterium.Communicated by Ramaswamy H. Sarma.

Computational methods and tools for binding site recognition between proteins and small molecules: from classical geometrical approaches to modern machine learning strategies

In the current "genomic era" the number of identified genes is growing exponentially. However, the biological function of a large number of the corresponding proteins is still unknown. Recognition of small molecule ligands (e.g., substrates, inhibitors, allosteric regulators, etc.) is pivotal for protein functions in the vast majority of the cases and knowledge of the region where these processes take place is essential for protein function prediction and drug design. In this regard, computational methods represent essential tools to tackle this problem. A significant number of software tools have been developed in the last few years which exploit either protein sequence information, structure information or both. This review describes the most recent developments in protein function recognition and binding site prediction, in terms of both freely-available and commercial solutions and tools, detailing the main characteristics of the considered tools and providing a comparative analysis of their performance.

Targeting of copper-trafficking chaperones causes gene-specific systemic pathology in Drosophila melanogaster: prospective expansion of mutational landscapes that regulate tumor resistance to cisplatin

Copper, a transition metal, is an essential component for normal growth and development. It acts as a critical co-factor of many enzymes that play key roles in diverse cellular processes. The present study attempts to investigate the regulatory functions decisively controlling copper trafficking during development and aging of the Drosophila model system. Hence, through engagement of the GAL4/UAS genetic platform and RNAi technology, we herein examined the in vivo significance of Atox1 and CCS genes, products of which pivotally govern cellular copper trafficking in fly tissue pathophysiology. Specifically, we analyzed the systemic effects of their targeted downregulation on the eye, wing, neuronal cell populations and whole-body tissues of the fly. Our results reveal that, in contrast to the eye, suppression of their expression in the wing leads to a notable increase in the percentage of malformed organs observed. Furthermore, we show that Atox1 or CCS gene silencing in either neuronal or whole-body tissues can critically affect the viability and climbing capacity of transgenic flies, while their double-genetic targeting suggests a rather synergistic mode of action of the cognate protein products. Interestingly, pharmacological intervention with the anti-cancer drug cisplatin indicates the major contribution of CCS copper chaperone to cisplatin's cellular trafficking, and presumably to tumor resistance often acquired during chemotherapy. Altogether, it seems that Atox1 and CCS proteins serve as tissue/organ-specific principal regulators of physiological Drosophila development and aging, while their tissue-dependent downregulation can provide important insights for Atox1 and CCS potential exploitation as predictive gene biomarkers of cancer-cell chemotherapy responses.

Summary: We demonstrate the essential roles of Atox1 and CCS copper-trafficking chaperones in Drosophila development and aging. We also provide insights for their therapeutic exploitation as cisplatin regulators during cancer chemotherapy.

I-TASSER gateway: A protein structure and function prediction server powered by XSEDE

There is an increasing gap between the number of known protein sequences and the number of proteins with experimentally characterized structure and function. To alleviate this issue, we have developed the I-TASSER gateway, an online server for automated and reliable protein structure and function prediction. For a given sequence, I-TASSER starts with template recognition from a known structure library, followed by full-length atomic model construction by iterative assembly simulations of the continuous structural fragments excised from the template alignments. Functional insights are then derived from comparative matching of the predicted model with a library of proteins with known function. The I-TASSER pipeline has been recently integrated with the XSEDE Gateway system to accommodate pressing demand from the user community and increasing computing costs. This report summarizes the configuration of the I-TASSER Gateway with the XSEDE-Comet supercomputer cluster, together with an overview of the I-TASSER method and milestones of its development.

Positive and Negative Regulation of the Virulence-Associated Coronafacoyl Phytotoxin in the Potato Common Scab Pathogen Streptomyces scabies

The potato common scab pathogen Streptomyces scabies produces N-coronafacoyl-l-isoleucine (CFA-Ile), which is a member of the coronafacoyl family of phytotoxins that are synthesized by multiple plant pathogenic bacteria. The CFA-Ile biosynthetic gene cluster contains a regulatory gene, cfaR, which directly controls the expression of the phytotoxin structural genes. In addition, a gene designated orf1 encodes a predicted ThiF family protein and is cotranscribed with cfaR, suggesting that it also plays a role in the regulation of CFA-Ile production. In this study, we demonstrated that CfaR is an essential activator of coronafacoyl phytotoxin production, while ORF1 is dispensable for phytotoxin production and may function as a helper protein for CfaR. We also showed that CFA-Ile inhibits the ability of CfaR to bind to the promoter region driving expression of the phytotoxin biosynthetic genes and that elevated CFA-Ile production by overexpression of both cfaR and orf1 in S. scabies increases the severity of disease symptoms induced by the pathogen during colonization of potato tuber tissue. Overall, our study reveals novel insights into the regulatory mechanisms controlling CFA-Ile production in S. scabies and it provides further evidence that CFA-Ile is an important virulence factor for this organism.

Mice expressing the variant rs1143679 allele of ITGAM (CD11b) show impaired DC-mediated T cell proliferation

Genome-wide association studies (GWAS) and functional genomic analyses have implicated several ITGAM (CD11b) single-nucleotide polymorphisms (SNPs) in the development of SLE and other disorders. ITGAM encodes the αM chain of the β2 integrin Mac-1, a receptor that plays important roles in myeloid cell functions. The ITGAM SNP rs1143679, which results in an arginine to histidine change at amino acid position 77 of the CD11b protein, has been shown to reduce binding to several ligands and to alter Mac-1-mediated cellular response in vitro. Importantly, however, the potential contribution of this SNP variant to the initiation and/or progression of immune and inflammatory processes in vivo remains unexplored. Herein, we describe for the first time the generation and characterization of a mouse line expressing the 77His variant of CD11b. Surprisingly, we found that 77His did not significantly affect Mac-1-mediated leukocyte migration and activation as assessed using thioglycollate-induced peritonitis and LPS/TNF-α-induced dermal inflammation models. In contrast, expression of this variant did alter T cell immunity, as evidenced by significantly reduced proliferation of ovalbumin (OVA)-specific transgenic T cells in 77His mice immunized with OVA. Reduced antigen-specific T cell proliferation was also observed when either 77His splenic dendritic cells (DCs) or bone marrow-derived DCs were used as antigen-presenting cells (APCs). Although more work is necessary to determine how this alteration might influence the development of SLE or other diseases, these in vivo findings suggest that the 77His variant of CD11b can compromise the ability of DCs to induce antigen-driven T cell proliferation.

Study of the contamination rate and change in growth features of lettuce (Lactuca sativa Linn.) in response to cadmium and a survey of its phytochelatin synthase gene

Crops can become contaminated when grown in soils containing heavy metals. Cadmium is a heavy metal that poses a significant health risk to humans. The purpose of this study was to evaluate the effect of cadmium on lettuce (Lactuca sativa Linn) and the contamination risk of lettuce grown in cadmium environments. The results showed that photosynthesis and growth parameters were significantly affected by cadmium. Lettuce has the ability to absorb large amounts of cadmium from the contaminated environment and so is a cadmium hyperaccumulator plant. The study showed that approximately 35% of the total absorbed cadmium is transmitted to aerial and edible parts of lettuce. This study was undertaken as lettuce has the ability to absorb and accumulate high levels of cadmium. There are however are no reports on the PCS gene and the potential for high cadmium accumulation in lettuce. The bioinformatics study revealed that lettuce has two phytochelatin synthase genes that produce 6 PCSs through splicing leading to the ability of lettuce to store high levels of cadmium. These six sequences although different in length have high similarity. Sequence structure, cellular location, three-dimensional structure, phylogeny and a comparison of their catalytic power were evaluated. The high accumulation of cadmium in lettuce and the presence of several PCSs contribute to the accumulation of cadmium in aerial tissues. The cultivation of lettuce in contaminated environments led us to evaluate suspected farms for the presence of cadmium in produce. Lettuce grown in industrial environments contaminated with cadmium can pose a serious threat to human health.

A novel protein descriptor for the prediction of drug binding sites

Background

Binding sites are the pockets of proteins that can bind drugs; the discovery of these pockets is a critical step in drug design. With the help of computers, protein pockets prediction can save manpower and financial resources.

Results

In this paper, a novel protein descriptor for the prediction of binding sites is proposed. Information on non-bonded interactions in the three-dimensional structure of a protein is captured by a combination of geometry-based and energy-based methods. Moreover, due to the rapid development of deep learning, all binding features are extracted to generate three-dimensional grids that are fed into a convolution neural network. Two datasets were introduced into the experiment. The sc-PDB dataset was used for descriptor extraction and binding site prediction, and the PDBbind dataset was used only for testing and verification of the generalization of the method. The comparison with previous methods shows that the proposed descriptor is effective in predicting the binding sites.

Conclusions

A new protein descriptor is proposed for the prediction of the drug binding sites of proteins. This method combines the three-dimensional structure of a protein and non-bonded interactions with small molecules to involve important factors influencing the formation of binding site. Analysis of the experiments indicates that the descriptor is robust for site prediction.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3058-0) contains supplementary material, which is available to authorized users.

Frontline Science: Antagonism between regular and atypical Cxcr3 receptors regulates macrophage migration during infection and injury in zebrafish

The CXCR3‐CXCL11 chemokine‐signaling axis plays an essential role in infection and inflammation by orchestrating leukocyte trafficking in human and animal models, including zebrafish. Atypical chemokine receptors (ACKRs) play a fundamental regulatory function in signaling networks by shaping chemokine gradients through their ligand scavenging function, while being unable to signal in the classic G‐protein‐dependent manner. Two copies of the CXCR3 gene in zebrafish, cxcr3.2 and cxcr3.3, are expressed on macrophages and share a highly conserved ligand‐binding site. However, Cxcr3.3 has structural characteristics of ACKRs indicative of a ligand‐scavenging role. In contrast, we previously showed that Cxcr3.2 is an active CXCR3 receptor because it is required for macrophage motility and recruitment to sites of mycobacterial infection. In this study, we generated a cxcr3.3 CRISPR‐mutant to functionally dissect the antagonistic interplay among the cxcr3 paralogs in the immune response. We observed that cxcr3.3 mutants are more susceptible to mycobacterial infection, whereas cxcr3.2 mutants are more resistant. Furthermore, macrophages in the cxcr3.3 mutant are more motile, show higher activation status, and are recruited more efficiently to sites of infection or injury. Our results suggest that Cxcr3.3 is an ACKR that regulates the activity of Cxcr3.2 by scavenging common ligands and that silencing the scavenging function of Cxcr3.3 results in an exacerbated Cxcr3.2 signaling. In human, splice variants of CXCR3 have antagonistic functions and CXCR3 ligands also interact with ACKRs. Therefore, in zebrafish, an analogous regulatory mechanism appears to have evolved after the cxcr3 gene duplication event, through diversification of conventional and atypical receptor variants.

CoABind: a novel algorithm for Coenzyme A (CoA)- and CoA derivatives-binding residues prediction

Motivation

Coenzyme A (CoA)-protein binding plays an important role in various cellular functions and metabolic pathways. However, no computational methods can be employed for CoA-binding residues prediction.

Results

We developed three methods for the prediction of CoA- and CoA derivatives-binding residues, including an ab initio method SVMpred, a template-based method TemPred and a consensus-based method CoABind. In SVMpred, a comprehensive set of features are designed from two complementary sequence profiles and the predicted secondary structure and solvent accessibility. The engine for classification in SVMpred is selected as the support vector machine. For TemPred, the prediction is transferred from homologous templates in the training set, which are detected by the program HHsearch. The assessment on an independent test set consisting of 73 proteins shows that SVMpred and TemPred achieve Matthews correlation coefficient (MCC) of 0.438 and 0.481, respectively. Analysis on the predictions by SVMpred and TemPred shows that these two methods are complementary to each other. Therefore, we combined them together, forming the third method CoABind, which further improves the MCC to 0.489 on the same set. Experiments demonstrate that the proposed methods significantly outperform the state-of-the-art general-purpose ligand-binding residues prediction algorithm COACH. As the first-of-its-kind method, we anticipate CoABind to be helpful for studying CoA-protein interaction.

Availability and implementation

http://yanglab.nankai.edu.cn/CoABind.

Supplementary information

Supplementary data are available at Bioinformatics online.

Involvement of peroxisome proliferator-activated receptor γ in anticonvulsant activity of α-asaronol against pentylenetetrazole-induced seizures in zebrafish

In mammals, peroxisome proliferators activated receptors (PPARs), the nuclear hormone receptors, have been reported to be involved in seizure control. Selective agonists and antagonists of PPARs raise seizure thresholds and suppress seizures, respectively. In this study, we evaluated the anticonvulsant effects of α-asaronol, a metabolic product of α-asarone, on pentylenetetrazole (PTZ)-induced seizures in zebrafish and investigated the underlying mechanisms. As a result, α-asaronol ameliorated seizures with increase of seizure latency, as well as decrease of seizure-like behavior, c-fos expression, and abnormal neuronal discharge in a concentration dependent manner. By comparing gene expression profiles of zebrafish undergoing seizures and α-asaronol pretreated zebrafish, we found that α-asaronol attenuate seizures through increase of PPAR γ expression, while PPAR γ antagonist GW9662 inhibit the anti-seizures actions of α-asaronol. Moreover, molecular docking simulation implied the physical interaction between α-asaronol and PPAR γ. The overall results indicated that the anticonvulsant effects of α-asaronol are regulated through PPAR γ-mediated pathway, which shed light on development of α-asaronol as a potential antiepileptic drug. In addition, it is for first time to report that PPAR γ is associated with seizures in zebrafish, supporting previous evidence that zebrafish is a suitable alternative for studying seizures.

A novel CDKN2A variant (p16L117P ) in a patient with familial and multiple primary melanomas

Germline mutations in CDKN2A (p16) are commonly found in patients with family history of melanoma or personal history of multiple primary melanomas. The p16 tumor suppressor gene regulates cell cycle progression and senescence through binding of cyclin-dependent kinases (CDK) and also regulates cellular oxidative stress independently of cell cycle control. We identified a germline missense (c.350T>C, p.Leu117Pro) CDKN2A mutation in a patient who had history of four primary melanomas, numerous nevi, and self-reported family history of melanoma. This particular CDKN2A mutation has not been previously reported in prior large studies of melanoma kindreds or patients with multiple primary melanomas. Compared with wild-type p16, the p16L117P mutant largely retained binding capacity for CDK4 and CDK6 but exhibited impaired capacity for repressing cell cycle progression and inducing senescence, while retaining its ability to reduce mitochondrial reactive oxygen species. Structural modeling predicted that the Leu117Pro mutation disrupts a putative adenosine monophosphate (AMP) binding pocket involving residue 117 in the fourth ankyrin domain. Identification of this new likely pathogenic variant extends our understanding of CDKN2A in melanoma susceptibility and implicates AMP as a potential regulator of p16.

Insight into the binding of a synthetic nitro-flavone derivative with human poly (ADP-ribose) polymerase 1

Flavones are important bioactive compounds, many of which are effective in cancer therapy for their ability to target enzymes related to DNA repair and cell proliferation. In this report, the interaction of a synthetic nitroflavone, 2,4-nitrophenylchromen-4-one (4NCO) with human poly (ADP-ribose) polymerase 1 (hPARP1) was investigated to explore its inhibitory action. Its interaction with hPARP1 was compared with that of other inhibitors through molecular docking studies. Further insight into the 4NCO-hPARP1 interaction was obtained from competitive docking and molecular dynamic simulation studies. In silico mutagenesis studies and per-residue interaction energy calculations were carried out. Quantitative Structure Activity Relationship analysis was also performed to calculate its predictive percent inhibitory activity. Our results indicated that 4NCO exhibited competitive mode of binding to hPARP1. It formed a stable interaction with the protein thereby hindering any further molecular interaction to render it inactive with a predictive inhibition of 96%. It also had good ADMET properties and showed best Autodock binding free energy values compared to other known inhibitors. 4NCO showed good hPARP1 inhibitory properties with higher bioavailability and lower probability of getting effluxed. Development of inhibitors against hPARP1 is important for cell proliferative disorders, where 4NCO can be predicted as a potential new drug.