Prediction of protein binding sites using physical and chemical descriptors and the support vector machine regression method

In this paper a new continuous variable called core-ratio is defined to describe the probability for a residue to be in a binding site, thereby replacing the previous binary description of the interface residue using 0 and 1. So we can use the support vector machine regression method to fit the core-ratio value and predict the protein binding sites. We also design a new group of physical and chemical descriptors to characterize the binding sites. The new descriptors are more effective, with an averaging procedure used. Our test shows that much better prediction results can be obtained by the support vector regression （SVR） method than by the support vector classification method.

An Approximated Voxel Approach for the Identification and Modelling of Ligand-Binding Sites

Most protein-ligand interactions take place on surfaces and include but not limited to factors such as chemical composition, hydrophobicity, electronegavitiy and shape complementarity. Past studies showed that protein-protein interactions occur on comparatively fiat regions whereas protein-ligand bindings involve crevices. In the search for such sites various approaches have been designed and developed each of which is algorithmically unique. The use of grid units or voxels has been demonstrated in early studies with relatively good results obtained. We present here an approximated approach comprising of the use of voxels and computer vision methods in the search for ligand-binding areas. Each test protein is modelled and analysed in 2D with all corresponding residues graphically presented for successfully identified sites. The study was carried out on 2 sets of proteins： FK506-bound proteins and heme-bound proteins with promising results obtained for all test cases.

Increased monoamine oxidase activity and imidazoline binding sites in insulin-resistant adipocytes from obese Zucker rats

BACKGROUND Despite overt insulin resistance,adipocytes of genetically obese Zucker rats accumulate the excess of calorie intake in the form of lipids.AIM To investigate whether factors can replace or reinforce insulin lipogenic action by exploring glucose uptake activation by hydrogen peroxide,since it is produced by monoamine oxidase(MAO)and semicarbazide-sensitive amine oxidase(SSAO)in adipocytes.METHODS 3H-2-deoxyglucose uptake(2-DG)was determined in adipocytes from obese and lean rats in response to insulin or MAO and SSAO substrates such as tyramine and benzylamine.14C-tyramine oxidation and binding of imidazolinic radioligands[3H-Idazoxan,3H-(2-benzofuranyl)-2-imidazoline]were studied in adipocytes,the liver,and muscle.The influence of in vivo administration of tyramine+vanadium on glucose handling was assessed in lean and obese rats.RESULTS 2-DG uptake and lipogenesis stimulation by insulin were dampened in adipocytes from obese rats,when compared to their lean littermates.Tyramine and benzylamine activation of hexose uptake was vanadate-dependent and was also limited,while MAO was increased and SSAO decreased.These changes were adipocyte-specific and accompanied by a greater number of imidazoline I2 binding sites in the obese rat,when compared to the lean.In vitro,tyramine precluded the binding to I2 sites,while in vivo,its administration together with vanadium lowered fasting plasma levels of glucose and triacylglycerols in obese CONCLUSION The adipocytes from obese Zucker rats exhibit increased MAO activity and imidazoline binding site number.However,probably as a consequence of SSAO down-regulation,the glucose transport stimulation by tyramine is decreased as much as that of insulin in these insulin-resistant adipocytes.The adipocyte amine oxidases deserve more studies with respect to their putative contribution to the management of glucose and lipid handling.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat(Triticum aestivum.L.,BBAADD)is hypothesized to increase its adaptability and/or plasticity.However,the molecular basis of expression divergence remains unclear.Squamosa promoter-binding protein-like(SPL)transcription factors are critical for a wide array of biological processes.In this study,we constructed expression regulatory networks by combining DAP-seq for 40 SPLs,ATACseq,and RNA-seq.Our findings indicate that a group of low-affinity SPL binding regions(SBRs)were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif.The SBRs including the low-affinity ones are evolutionarily conserved,enriched GWAS signals related to important agricultural traits.However,those SBRs are highly diversified among the cis-regulatory regions(CREs)of syntenic genes,with less than 8%SBRs coexisting in triad genes,suggesting that CRE variations are critical for subgenome differentiations.Knocking out of Ta SPL7A/B/D and Ta SPL15A/B/D subfamily further proved that both high-and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes.Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Addition of heparin binding sites strongly increases the bone forming capabilities of BMP9 in vivo

Bone Morphogenetic proteins(BMPs)like BMP2 and BMP7 have shown great potential in the treatment of severe bone defects.In recent in vitro studies,BMP9 revealed the highest osteogenic potential compared to other BMPs,possibly due to its unique signaling pathways that differs from other osteogenic BMPs.However,in vivo the bone forming capacity of BMP9-adsorbed scaffolds is not superior to BMP2 or BMP7.In silico analysis of the BMP9 protein sequence revealed that BMP9,in contrast to other osteogenic BMPs such as BMP2,completely lacks so-called heparin binding motifs that enable extracellular matrix(ECM)interactions which in general might be essential for the BMPs’osteogenic function.Therefore,we genetically engineered a new BMP9 variant by adding BMP2-derived heparin binding motifs to the N-terminal segment of BMP9′s mature part.The resulting protein(BMP9 HB)showed higher heparin binding affinity than BMP2,similar osteogenic activity in vitro and comparable binding affinities to BMPR-II and ALK1 compared to BMP9.However,remarkable differences were observed when BMP9 HB was adsorbed to collagen scaffolds and implanted subcutaneously in the dorsum of rats,showing a consistent and significant increase in bone volume and density compared to BMP2 and BMP9.Even at 10-fold lower BMP9 HB doses bone tissue formation was observed.This innovative approach of significantly enhancing the osteogenic properties of BMP9 simply by addition of ECM binding motifs,could constitute a valuable replacement to the commonly used BMPs.The possibility to use lower protein doses demonstrates BMP9 HB’s high translational potential.

Evolutionary trace analysis of eukaryotic DNA topoisomerase I superfamily:Identification of novel antitumor drug binding site

The studies of novel inhibitors of DNA topoisomerase I(Topo I)have already be-come very promising in cancer chemotherapy.Identifying the new drug-binding residues is playing an important role in the design and optimization of Topo I inhibitors.The designed com-pounds may have novel scaffolds,thus will be helpful to overcome the toxicities of current camptothecin(CPT)drugs and may provide a solution to cross resistance with these drugs.Mul-tiple sequence alignments were performed on eukaryotic DNA topoisomerase I superfamily and thus the evolutionary tree was constructed.The Evolutionary Trace method was applied to iden-tify functionally important residues of human Topo I.It has been demonstrated that class-specific hydrophobic residues Ala351,Met428,Pro431 are located around the 7,9-position of CPT,indi-cating suitable substitution of hydrophobic group on CPT will increase antitumor activity.The conservative residue Lys436 in the superfamily is of particular interest and new CPT derivatives designed based on this residue may greatly increase water solubility of such drugs.It has also been demonstrated that the residues Asn352 and Arg364 were conservative in the superfamily,whose mutation will render CPT resistance.As our molecular docking studies demonstrated they did not make any direct interaction with CPT,they are important drug-binding site residues for future design of novel non-camptothecin lead compounds.This work provided a strong basis for the design and synthesis of novel highly potent CPT derivatives and virtual screening for novel lead compounds.

Arabidopsis AtADF1 is Functionally Affected by Mutations on Actin Binding Sites

The plant actin depolymerizing factor （ADF） binds to both monomeric and filamentous actin, and is directly involved in the depolymerization of actin filaments. To better understand the actin binding sites of the Arabidopsis thaliana L. AtADF1, we generated mutants of AtADF1 and investigated their functions in vitro and in vivo. Analysis of mutants harboring amino acid substitutions revealed that charged residues （Arg98 and Lys100） located at the α-helix 3 and forming an actin binding site together with the N-terminus are essential for both G- and F-actin binding. The basic residues on the β-strand 5 （K82/A） and the α-helix 4 （R135/A, R137/A） form another actin binding site that is important for F-actin binding. Using transient expression of CFP-tagged AtADF1 mutant proteins in onion （Allium cepa） peel epidermal cells and transgenic Arabidopsis thaliana L. plants overexpressing these mutants, we analyzed how these mutant proteins regulate actin organization and affect seedling growth. Our results show that the ADF mutants with a lower affinity for actin filament binding can still be functional, unless the affinity for actin monomers is also affected. The G-actin binding activity of the ADF plays an essential role in actin binding, depolymerization of actin polymers, and therefore in the control of actin organization.

Functional characterization of the NF-kB binding site in the human NOD2 promoter

Nucleotide-binding and oligomerization domain 2(NOD2),a member of the NOD protein family,plays an important role in innate immunity.In response to pathogen attack,NOD2 stimulates cytokine and defensin production by activating nuclear factor(NF)-kB,a key transcription factor responsible for mediating downstream reactions.However,the mechanism linking NOD2 regulation and NF-kB activation is poorly understood.Using bioinformatics,we found a completely preserved canonical NF-kB binding site in the NOD2 core promoter(216 to 225 bp)in both humans and chimpanzees.The functional role of this NF-kB binding site was investigated using the enhanced green fluorescent protein(EGFP)reporter system,site-directed mutagenesis,the NF-kB activation inhibitor(JSH-23)and the chromatin immunoprecipitation(ChIP)assay.The results show that the NF-kB binding site is critical for regulation of the NOD2 gene.Either deletion of the NF-kB binding elements within the NOD2 promoter or treatment with an NF-kB activation inhibitor could lead to a significant loss of NOD2 promoter activity as detected by reporter gene assay.The canonical NF-kB binding site was bound by NF-kB as determined by the ChIP method.Based on these results,we suggest a positive feedback regulation between NF-kB and NOD2,which may represent an efficient mechanism in response to pathogen invasion.

Systematic identification and annotation of multiple-variant compound effects at transcription factor binding sites in human genome

Understanding the functional effects of genetic variants is crucial in modern genomics and genetics. Transcription factor binding sites （TFBSs） are one of the most important cis-regulatory elements. While multiple tools have been developed to assess functional effects of genetic variants at TFBSs, they usually assume that each variant works in isolation and neglect the potential ＂interference＂ among multiple variants within the same TFBS. In this study, we presented COPE-TFBS （Context-Oriented Predictor for variant Effect on Transcription Factor Binding Site）, a novel method that considers sequence context to accurately predict variant effects on TFBSs. We systematically re-analyzed the sequencing data from both the 1000 Genomes Project and the Genotype-Tissue Expression （GTEx） Project via COPE-TFBS, and identified numbers of novel TFBSs, transformed TFBSs and discordantly annotated TFBSs resulting from multiple variants, further highlighting the necessity of sequence context in accurately annotating genetic variants.

Structure-Based Prediction of Transcription Factor Binding Sites

Transcription Factors（TFs） are a very diverse family of DNA-binding proteins that play essential roles in the regulation of gene expression through binding to specific DNA sequences. They are considered as one of the prime drug targets since mutations and aberrant TF-DNA interactions are implicated in many diseases.Identification of TF-binding sites on a genomic scale represents a critical step in delineating transcription regulatory networks and remains a major goal in genomic annotations. Recent development of experimental high-throughput technologies has provided valuable information about TF-binding sites at genome scale under various physiological and developmental conditions. Computational approaches can provide a cost-effective alternative and complement the experimental methods by using the vast quantities of available sequence or structural information. In this review we focus on structure-based prediction of transcription factor binding sites. In addition to its potential in genomescale predictions, structure-based approaches can help us better understand the TF-DNA interaction mechanisms and the evolution of transcription factors and their target binding sites. The success of structure-based methods also bears a translational impact on targeted drug design in medicine and biotechnology.

Studies of Site Specific DNA Binding of Small Peptides by Competitive Assays with Hoechst 33258

With a view to finding out precisely how small peptides recognize a particular binding site of DNA, we have accomplished DNA binding studies of two peptides, H-Tyr-Arg-OH (YR) and H-Gly-Gly-His-OH (GGH) by using measurements in comparison with the binding between DNA and Hoechst 33258. The inhibition mode by YR and GGH to DNA binding of Hoechst 33258 was analyzed by Lineweaver-Burk plot which shows the plot of typical competitive inhibition at concentration of Hoechst 33258 from 3.66 ( 10-9 mol / L to 1.09 ( 10-8 mol / L. And it is concluded that YR binds to DNA in its minor groove (AT rich regions) with a binding constant K = 1.02 ( 108 (mol / L)-1. The GGH(s specificity is reduced at high concentration because it can also bind GC base pair.

TrFAST: A Tool to Predict Signaling Pathway-specific Transcription Factor Binding Sites

Recent advances in the development of high-throughput tools have significantly revolutionized our understanding of molecular mech- anisms underlying normal and dysfunctional biological processes. Here we present a novel computational tool, transcription factor search and analysis tool （TrFAST）, which was developed for the in silico analysis of transcription factor binding sites （TFBSs） of sig- naling pathway-specific TFs. TrFAST facilitates searching as well as comparative analysis of regulatory motifs through an exact pattern matching algorithm followed by the graphical representation of matched binding sites in multiple sequences up to 50 kb in length. TrFAST is proficient in reducing the number of comparisons by the exact pattern matching strategy. In contrast to the pre-existing tools that find TFBS in a single sequence, TrFAST seeks out the desired pattern in multiple sequences simultaneously. It counts the GC con- tent within the given multiple sequence data set and assembles the combinational details of consensus sequence（s） located at these regions, thereby generating a visual display based on the abundance of unique pattern. Comparative regulatory region analysis of multi- ple orthologous sequences simultaneously enhances the features of TrFAST and provides a significant insight into study of conservation of non-coding cis-regulatory elements. TrFAST is freely available at http：//www.fi-pk.com/trfast.html.

Nucleosomal Context of Binding Sites Influences Transcription Factor Binding Affinity and Gene Regulation

Transcription factor （TF） binding to its DNA target site plays an essential role in gene regulation. The location, orientation and spacing of transcription factor binding sites （TFBSs） also affect regulatory function of the TF. However, how nucleosomal context of TFBSs influences TF binding and subsequent gene regulation remains to be elucidated. Using genome-wide nucleosome positioning and TF binding data in budding yeast, we found that binding affinities of TFs to DNA tend to decrease with increasing nucleosome occupancy of the associated binding sites. We further demonstrated that nucleosomal context of binding sites is correlated with gene regulation of the corresponding TF. Nucleosome-depleted TFBSs are linked to high gene activity and low expression noise, whereas nucleosome-covered TFBSs are associated with low gene activity and high expression noise. Moreover, nucleosome-covered TFBSs tend to disrupt coexpression of the corresponding TF target genes. We conclude that nucleosomal context of binding sites influences TF binding affinity, subsequently affecting the regulation of TFs on their target genes. This emphasizes the need to include nucleosomal context of TFBSs in modeling gene regulation.

Computational prediction and validation of specific EmbR binding site on PknH

Tuberculosis drug resistance continues to threaten global health but the underline molecular mechanisms are not clear.Ethambutol(EMB),one of the well-known first-line drugs in tuberculosis treatment is,unfortunately,not free from drug resistance problems.Genomic studies have shown that some genetic mutations in Mycobacterium tuberculosis(Mtb)EmbR,and EmbC/A/B genes cause EMB resistance.EmbR-PknH pair controls embC/A/B operon,which encodes EmbC/A/B genes,and EMB interacts with EmbA/B proteins.However,the EmbR binding site on PknH was unknown.We conducted molecular simulation on the EmbR-peptides binding structures and discovered phosphorylated PknH 273-280(N′-HEALS^(P)DPD-C′)makesβstrand with the EmbR FHA domain,asβ-MoRF(MoRF;molecular recognition feature)does at its binding site.Hydrogen bond number analysis also supported the peptides’β-MoRF forming activity at the EmbR FHA domain.Also,we discovered that previously known phosphorylation residues might have their chronological order according to the phosphorylation status.The discovery validated that Mtb PknH 273-280(N′-HEALSDPD-C′)has reliable EmbR binding affinity.This approach is revolutionary in the computer-aided drug discovery field,because it is the first trial to discover the protein-protein interaction site,and find binding partner in nature from this site.

Structure effect of nucleotides in terbium(Ⅲ)-nucleotide fluorescent reaction—-New evidence for the binding sites of terbium(Ⅲ) on nucleotides

Upon addition of Tb^(3+) to 16 nucleotides and homopolynucleotides, all of them showed a characteristic green emission from Tb^(3+), but with much different intensity, upon excitation in the aromatic region of bases. The result suggested that nucleotides with at least one carbonyl group in nucleotide bases are better enhancers to the fluorescence of Tb3+. The complexes of ATP, GDP and GTP with T5^(3+) are synthesized as two types of models. Guanine tpye nucleotides with one carbonyl group in the bases are the best enhancers, while adenine type nucleotides with no carbonyl group in the bases are poorest enhancers to the fluorescence of Tb^(3+). Comparing the IR spectra of ATP, GTP, GDP and their Tb^(3+) complexes suggested that C-6 carbonyl group in GTP and GDP may be involved in complex formation, which may be responsible for the effective energy transfer. This is further supported by comparing the UV spectra of ATP, Poly(A), GTP, and Poly(G) with their Tb^(3+) complexes in water solution.

3-Arylisothiazoloquinols as potent ligands for the benzodiazepine site of GABA_(A) receptors

3-Arylisothiazolo[5,4-b]quinolin-4(9H)-ones and 3-arylisoxazolo[5,4-b]quinolin-4(9H)-ones were synthesized and assayed for affinity for the benzodiazepine binding site of the GABAA receptors. While the 3-arylisothiazoloquinolin-4-ones were found to be potent ligands, with affinities (expressed as the affinity Ki value) down to 1 nM, the 3-arylisoxazoloquinolin-4-ones are less potent. This is suggested to depend on sterical repulsive interaction of the 3-arylisoxazoloquinolin-4-ones with the receptor essential volume of the binding site, and a higher electron density at the nitrogen in the azole ring (N-2) as well as the carbonyl oxygen in the isothiazoloquinolin-4-ones enabling them to interact stronger with hydrogen bond donor sites at the binding site.

基于3A组装的多拷贝基因策略优化重组水蛭素变体Ⅲ的生产

水蛭素变体Ⅲ(Hirudin varidantⅢ,Hv3)是一种从水蛭中提取的活性成分,在预防和治疗白内障方面具有潜在作用。为了制备足量Hv3用于进一步的临床前应用研究,该研究开发了一种在谷氨酸棒杆菌(Corynebacterium glutamicum)中高效、稳定表达重组水蛭素变体Ⅲ(recombinant Hirudin varidantⅢ,Rhv3)的办法。首先,在C.glutamicum中构建含Ptac启动子和CspA信号肽的Rhv3分泌表达菌株,比较其在不同培养基中的表达情况。结果表明在BHI培养基中Rhv3活性最高,达到3.02×10^(3) ATU/L。为进一步提升Rhv3产量,通过核糖体结合位点(ribosome binding site,RBS)筛选,选用RBS1应用于Rhv3多拷贝菌株构建和表达。然后利用3A组装技术构建含不同信号肽和不同拷贝数的Rhv3分泌表达菌株进一步优化其产量。通过放大发酵培养,Rhv3产量达到1.89 g/L,活性达到10.91×10^(3) ATU/L。最后利用镍柱对Rhv3进行简单纯化,其纯度达到90%。该异源表达策略有效提高Rhv3表达量,为Rhv3的重组生产提供了一种质量可靠、低成本的表达体系。

Cloning and Analysis of a Disease Resistance Gene Homolog from Soybean

Conserved domains e.g. nucleotide binding site (NBS) were found in several cloned plant disease resistance genes. Based on the NBS domain, resistance gene analogs (RGAs) have been isolated previously and were used as probes to screen a soybean (Glycine max L. Merr.) cDNA library. A full-length cDNA, KR3, was obtained by screening the library and rapid amplification of cDNA ends (RACE) method. Sequence analysis revealed that the cDNA is 2 353 bp in length and the open reading frame (ORF) codes for a polypeptide of 636 amino acids with a Toll-Interleukin-1 receptor (TIR) and a NBS domain. Sequence alignment showed that it was similar to N gene of tobacco. The phylogenetic tree analysis of R proteins with NBS from higher plants was performed. The KR3 gene has low copies in soybean genome and its expression was induced by exogenous salicylic acid (SA).

Interactions Between Anticancer Pt(II) complexes and Human Erythrocyte Spectrin *

The interactions between human erythrocyte spectrin(SP) and Pt(II) complexes with different composition and configuration were studied by fluorescence and circular dichroism spectra. The results showed that there are 4.7×10 2 binding sites of cisplatin(CDDP) in a spectrin tetramer(SPT). Among them, about 70 sites with apparent binding constant K 1】3.47×10 6 were of highest affinity, 1.8×10 2 sites with K 2 = 3.47×10 6 were of high affinity, and other 2.2×10 2 sites with K 3 = 8.77×10 5 were of low affinity. The conformation change of spectrin, depending on the concentration of Pt(II) complex and molar ratio(R) of Pt(II) complex to spectrin, was induced by the binding of Pt(II) complexes. It indicated that the interaction of both CDDP and cis diaquodiamine platinum(DADP) with SP followed a two step first order kinetic process in the first stage (1 h), and the kinetic constants were determined. In the second stage, the induced conformation change, polymerization and depolymerization of SP were probably involved. It was noticed that in the reaction of SP and Pt(II) complexes with 1,2 cyclohexanediammine isomers as chiral carrier ligand, stereo matching played a more important role than the affinity of Pt(II) to thiol groups of SP.

Analysis of three types of resistance gene analogs in Pm U region from Triticum urartu

Resistance gene analog(RGA) screening of mapped disease-resistant genes not only helps to clone these genes but also helps to develop efficient molecular markers for resistance breeding. The present study focused on the PmU region located on chromosome 7 Au L of Triticum urartu, and recently, a nucleotide binding site(NBS)-encoding gene, Pm60, was cloned from the same chromosome arm. In this research, NBS, protein kinase(PK), and ATP-binding cassette(ABC), the three disease resistance-related gene families, were analyzed within PmU region by using informatics tools, and an expression experiment was conducted to verify their functions in vivo. Comparative genomic analysis revealed that 126 RGAs were included on chromosome 7 Au L, and 30 of the RGAs as well as Pm60 were found in the Pm U region. Transcriptome database analysis of T. urartu revealed 14 PmU-RGAs with expression data, and three PmU-NBSs exhibited significant changes in expression after inoculation with Blumeria graminis f. sp. tritici(Bgt); TRIUR314879 was up-regulated, while TRIUR300450 and TRIUR306270 were down-regulated. Cluster analysis showed that these three PmU-NBSs were clustered far from the cloned wheat resistance genes. Then, qRT-PCR was performed to investigate the expression of 14 PmU-RGAs and Pm60 after inoculation with Bgt race E09; the results showed that Pm60 was specifically expressed in UR206 which carrying PmU, but not in susceptible UR203; while TRIUR314879 was significantly up-regulated and TRIUR300450 was downregulated in UR206 after inoculation. These results indicated that PmU is Pm60, and TRIUR314879 and TRIUR300450 may also be involved in the defense against Bgt.