In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.

Nuclear magnetic resonance structure shows that the severe acute respiratory syndrome coronavirus-unique domain contains a macrodomain fold

The nuclear magnetic resonance (NMR) structure of a central segment of the previously annotated severe acute respiratory syndrome (SARS)-unique domain (SUD-M, for "middle of the SARS-unique domain") in SARS coronavirus (SARS-CoV) nonstructural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3 residues 528 to 648, and there is a flexibly extended N-terminal tail with the residues 513 to 527 and a C-terminal flexible tail of residues 649 to 651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527 to 651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly(A) and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In a further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1 ''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows three-dimensional structure homology with several helicases and nucleoside triphosphate-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

Morphology and myofiber composition of skeletal musculature of the forelimb in young and aged wild type and myostatin null mice

Most current research into therapeutic approaches to muscle diseases involves the use of the mouse as an experimental model. Furthermore, a major strategy to alleviate myopathic symptoms through enhancing muscle growth and regeneration is to inhibit the action of myostatin (Mstn), a transforming growth factor-beta (TGF-beta) family member that inhibits muscle growth. Presently, however, no study has expanded the morphological analysis of mouse skeletal muscle beyond a few individual muscles of the distal hindlimb, through which broad conclusions have been based. Therefore, we have initially undertaken an expansive analysis of the skeletal musculature of the mouse forelimb and highlighted the species-specific differences between equivalent muscles of the rat, another prominently used experimental model. Subsequently, we examined the musculature of the forelimb in both young and old adult wild-type (mstn(+/+)) and myostatin null (mstn(-/-)) mice and assessed the potential beneficial and detrimental effects of myostatin deletion on muscle morphology and composition during the aging process. We showed that: (1) the forelimb muscles of the mouse display a more glycolytic phenotype than those of the rat; (2) in the absence of myostatin, the induced myofiber hyperplasia, hypertrophy, and glycolytic conversion all occur in a muscle-specific manner; and, importantly, (3) the loss of myostatin significantly alters the dynamics of postnatal muscle growth and impairs age-related oxidative myofiber conversion.

NMR structure shows that the SARS-unique domain contains a macrodomain fold

The NMR structure of a central segment of the previously annotated "SARS-unique domain" (SUD-M; "middle of the SARS-unique domain") in the SARS coronavirus (SARS-CoV) non-structural protein 3 (nsp3) has been determined. SUD-M(513-651) exhibits a macrodomain fold containing the nsp3-residues 528-648, and there is a flexibly extended N-terminal tail with the residues 513-527 and a C-terminal flexible tail of residues 649-651. As a follow-up to this initial result, we also solved the structure of a construct representing only the globular domain of residues 527-651 [SUD-M(527-651)]. NMR chemical shift perturbation experiments showed that SUD-M(527-651) binds single-stranded poly-A and identified the contact area with this RNA on the protein surface, and electrophoretic mobility shift assays then confirmed that SUD-M has higher affinity for purine bases than for pyrimidine bases. In further search for clues to the function, we found that SUD-M(527-651) has the closest three-dimensional structure homology with another domain of nsp3, the ADP-ribose-1''-phosphatase nsp3b, although the two proteins share only 5% sequence identity in the homologous sequence regions. SUD-M(527-651) also shows 3D structure homology with several helicases and NTP-binding proteins, but it does not contain the motifs of catalytic residues found in these structural homologues. The combined results from NMR screening of potential substrates and the structure-based homology studies now form a basis for more focused investigations on the role of the SARS-unique domain in viral infection.

Differential expression of mRNAs encoding co-receptor (betaglycan) and activin type-I preovulatory and prehierarchical follicles of the putative inhibin and type-II receptors in the laying hen ovary

Ovarian follicle development is primarily regulated by an interplay between the pituitary gonadotrophins, LH and FSH, and ovary-derived steroids. Increasing evidence implicates regulatory roles of transforming growth factor-beta (TGF beta) superfamily members, including inhibins and activins. The aim of this study was to identify the expression of mRNAs encoding key receptors of the inhibin/activin system in ovarian follicles ranging from 4 mm in diameter to the dominant F1 follicle (similar to 40 turn). Ovaries were collected (n=16) from inid-sequence hens maintained on a long-day photoschedule (16h of light:8 h of darkness). All follicles removed were dissected into individual granulosa and thecal layers. RNA was extracted and cDNA synthesized. Real-time quantitative PCR was used to quantify the expression of niRNA encoding betaglycan, activin receptor (ActR) subtypes (type-I, -IIA and -IIB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH); receptor expression data were normalized to GAPDH expression. Detectable levels of ActRI, -IIA and -IIB and the inhibin co-receptor (betaglycan) expression were found in all granulosa and thecal layers analysed. Granulosa ActRI mRNA peaked (P < 0(.)05) in 8-9(.)9 mm follicles, whereas ActRIIA rose significantly from 6-7(.)9 mm to 8-9(.)9 nun, before filling to F3/2; levels then rose sharply (3-fold) to F1 levels. Granulosa betaglycan niRNA expression rose 3-fold from 4-5(.)9 min to 8-9(.)9 mm, before falling 4-fold to F3/2; levels then rose sharply (4-fold) to F1 levels. ActRIIB levels did not vary significantly during follicular development. Thecal ActRI mRNA expression was similar from 4-7(.)9 mm then decreased significantly to a nadir at the F4 position, before increasing 2-fold to the F1 (P < 0(.)05). Although thecal ActRIIB and -IIA expression did not vary significantly from 4 nim to F3, ActRIIB expression increased significantly (2-fold) from F3 to F1 and ActIIA, increased 22-fold from F2 to F1 (P < 0(.)05). Thecal betaglycan fell to a nadir at F6 after follicle selection; levels then increased significantly to F2, before filling similar to 50% in the F I. In all follicles studied expression of betaglycan and ActRI (granulosa: 1-0(.)65, P < 0-001, n=144/group; theca: r=0(.)49, P < 0-001, n=144/group) was well correlated. No significant correlations were identified between betaglycan and ActRIIA or -IIB. Considering all follicles analysed, granulosa mRNA expression of betaglycan, ActRI ActRIIA and ActRIIB were all significantly lower than in corresponding thecal tissue (betaglycan, 11(.)4-fold; ActRIIB, 5(.)1-fold; ActR(.) 3-8-fold: ActRIIA, 2(.)8-fold). The co-localization of type-I and -II activin receptors and betaglycan on granulosa and thecal cells are consistent with a local auto/paracrine role of inhibins and activins in modulating ovarian follicle development, selection and progression in the domestic fowl.

Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes

Homeobox genes encode DNA-binding proteins, many of which are implicated in the control of embryonic development. Evolutionarily, most homeobox genes fall into two related clades: the ANTP and the PRD classes. Some genes in ANTP class, notably Hox, ParaHox, and NK genes, have an intriguing arrangement into physical clusters. To investigate the evolutionary history of these gene clusters, we examined homeobox gene chromosomal locations in the cephalochordate amphioxus, Branchiostoma floridae. We deduce that 22 amphioxus ANTP class homeobox genes localize in just three chromosomes. One contains the Hox cluster plus AmphiEn, AmphiMnx, and AmphiDll. The ParaHox cluster resides in another chromosome, whereas a third chromosome contains the NK type homeobox genes, including AmphiMsx and ArnphiTlx. By comparative analysis we infer that clustering of ANTP class homeobox genes evolved just once, during a series of extensive cis-duplication events of genes early in animal evolution. A trans-duplication event occurred later to yield the Hox and ParaHox gene clusters on different chromosomes. The results obtained have implications for understanding the origin of homeobox gene clustering, the diversification of the ANTP class of homeobox genes, and the evolution of animal genomes.

Differential expression of mRNAs encoding co-receptor (betaglycan) and activin type-I preovulatory and prehierarchical follicles of the putative inhibin and type-II receptors in the laying hen ovary

Ovarian follicle development is primarily regulated by an interplay between the pituitary gonadotrophins, LH and FSH, and ovary-derived steroids. Increasing evidence implicates regulatory roles of transforming growth factor-beta (TGF beta) superfamily members, including inhibins and activins. The aim of this study was to identify the expression of mRNAs encoding key receptors of the inhibin/activin system in ovarian follicles ranging from 4 mm in diameter to the dominant F1 follicle (similar to 40 turn). Ovaries were collected (n=16) from inid-sequence hens maintained on a long-day photoschedule (16h of light:8 h of darkness). All follicles removed were dissected into individual granulosa and thecal layers. RNA was extracted and cDNA synthesized. Real-time quantitative PCR was used to quantify the expression of niRNA encoding betaglycan, activin receptor (ActR) subtypes (type-I, -IIA and -IIB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH); receptor expression data were normalized to GAPDH expression. Detectable levels of ActRI, -IIA and -IIB and the inhibin co-receptor (betaglycan) expression were found in all granulosa and thecal layers analysed. Granulosa ActRI mRNA peaked (P < 0(.)05) in 8-9(.)9 mm follicles, whereas ActRIIA rose significantly from 6-7(.)9 mm to 8-9(.)9 nun, before filling to F3/2; levels then rose sharply (3-fold) to F1 levels. Granulosa betaglycan niRNA expression rose 3-fold from 4-5(.)9 min to 8-9(.)9 mm, before falling 4-fold to F3/2; levels then rose sharply (4-fold) to F1 levels. ActRIIB levels did not vary significantly during follicular development. Thecal ActRI mRNA expression was similar from 4-7(.)9 mm then decreased significantly to a nadir at the F4 position, before increasing 2-fold to the F1 (P < 0(.)05). Although thecal ActRIIB and -IIA expression did not vary significantly from 4 nim to F3, ActRIIB expression increased significantly (2-fold) from F3 to F1 and ActIIA, increased 22-fold from F2 to F1 (P < 0(.)05). Thecal betaglycan fell to a nadir at F6 after follicle selection; levels then increased significantly to F2, before filling similar to 50% in the F I. In all follicles studied expression of betaglycan and ActRI (granulosa: 1-0(.)65, P < 0-001, n=144/group; theca: r=0(.)49, P < 0-001, n=144/group) was well correlated. No significant correlations were identified between betaglycan and ActRIIA or -IIB. Considering all follicles analysed, granulosa mRNA expression of betaglycan, ActRI ActRIIA and ActRIIB were all significantly lower than in corresponding thecal tissue (betaglycan, 11(.)4-fold; ActRIIB, 5(.)1-fold; ActR(.) 3-8-fold: ActRIIA, 2(.)8-fold). The co-localization of type-I and -II activin receptors and betaglycan on granulosa and thecal cells are consistent with a local auto/paracrine role of inhibins and activins in modulating ovarian follicle development, selection and progression in the domestic fowl.

Neoglycolipid probes prepared via oxime ligation for microarray analysis of oligosaccharide-protein interactions

Neoglycolipid technology is the basis of a microarray platform for assigning oligosaccharide ligands for carbohydrate-binding proteins. The strategy for generating the neoglycolipid probes by reductive amination results in ring opening of the core monosaccharides. This often limits applicability to short-chain saccharides, although the majority of recognition motifs are satisfactorily presented with neoglycolipids of longer oligosaccharides. Here, we describe neoglycolipids prepared by oxime ligation. We provide evidence from NMR studies that a significant proportion of the oxime-linked core monosaccharide is in the ring-closed form, and this form selectively interacts with a carbohydrate-binding protein. By microarray analyses we demonstrate the effective presentation with oxime-linked neoglycolipids of (1) Lewis(x) trisaccharide to antibodies to Lewisx, (2) sialyllactose analogs to the sialic acid-binding receptors, siglecs, and (3) N-glycans to a plant lectin that requires an intact N-acetylglucosamine core.

Gene transcription in Daphnia magna: effects of acute exposure to a carbamate insecticide and an acetanilide herbicide.

Daphnia magna is a key invertebrate in the freshwater environment and is used widely as a model in ecotoxicological measurements and risk assessment. Understanding the genomic responses of D. magna to chemical challenges will be of value to regulatory authorities worldwide. Here we exposed D. magna to the insecticide methomyl and the herbicide propanil to compare phenotypic effects with changes in mRNA expression levels. Both pesticides are found in drainage ditches and surface water bodies standing adjacent to crops. Methomyl, a carbamate insecticide widely used in agriculture, inhibits acetylcholinesterase, a key enzyme in nerve transmission. Propanil, an acetanilide herbicide, is used to control grass and broad-leaf weeds. The phenotypic effects of single doses of each chemical were evaluated using a standard immobilisation assay. Immobilisation was linked to global mRNA expression levels using the previously estimated 48h-EC(1)s, followed by hybridization to a cDNA microarray with more than 13,000 redundant cDNA clones representing >5000 unique genes. Following exposure to methomyl and propanil, differential expression was found for 624 and 551 cDNAs, respectively (one-way ANOVA with Bonferroni correction, Pproteins, ATP synthesis-related proteins), moulting (e.g., chitin-binding proteins, cuticular proteins) and protein biosynthesis (e.g., ribosomal proteins, transcription factors). Methomyl induced the transcription of genes involved in specific processes such as ion homeostasis and xenobiotic metabolism. Propanil highly promoted haemoglobin synthesis and up-regulated genes specifically related to defence mechanisms (e.g., innate immunity response systems) and neuronal pathways. Pesticide-specific toxic responses were found but there is little evidence for transcriptional responses purely restricted to genes associated with the pesticide target site or mechanism of toxicity.

Analysis of aluminium content and iron homeostasis in nipple aspirate fluids from healthy women and breast cancer-affected patients

Aluminium is not a physiological component of the breast but has been measured recently in human breast tissues and breast cyst fluids at levels above those found in blood serum or milk. Since the presence of aluminium can lead to iron dyshomeostasis, levels of aluminium and iron-binding proteins (ferritin, transferrin) were measured in nipple aspirate fluid (NAF), a fluid present in the breast duct tree and mirroring the breast microenvironment. NAFs were collected noninvasively from healthy women (NoCancer; n = 16) and breast cancer-affected women (Cancer; n = 19), and compared with levels in serum (n = 15) and milk (n = 45) from healthy subjects. The mean level of aluminium, measured by ICP-mass spectrometry, was significantly higher in Cancer NAF (268.4 ± 28.1 μg l−1; n = 19) than in NoCancer NAF (131.3 ± 9.6 μg l−1; n = 16; P < 0.0001). The mean level of ferritin, measured through immunoassay, was also found to be higher in Cancer NAF (280.0 ± 32.3 μg l−1) than in NoCancer NAF (55.5 ± 7.2 μg l−1), and furthermore, a positive correlation was found between levels of aluminium and ferritin in the Cancer NAF (correlation coefficient R = 0.94, P < 0.001). These results may suggest a role for raised levels of aluminium and modulation of proteins that regulate iron homeostasis as biomarkers for identification of women at higher risk of developing breast cancer. The reasons for the high levels of aluminium in NAF remain unknown but possibilities include either exposure to aluminium-based antiperspirant salts in the adjacent underarm area and/or preferential accumulation of aluminium by breast tissues.

Aluminium and human breast diseases

The human breast is exposed to aluminium from many sources including diet and personal care products, but dermal application of aluminium-based antiperspirant salts provides a local long-term source of exposure. Recent measurements have shown that aluminium is present in both tissue and fat of the human breast but at levels which vary both between breasts and between tissue samples from the same breast. We have recently found increased levels of aluminium in noninvasively collected nipple aspirate fluids taken from breast cancer patients (mean 268±28 g/l) compared with control healthy subjects (mean 131±10 g/l) providing evidence of raised aluminium levels in the breast microenvironment when cancer is present. The measurement of higher levels of aluminium in type I human breast cyst fluids (median 150g/l) compared with human serum (median 6g/l) or human milk (median 25g/l) warrants further investigation into any possible role of aluminium in development of this benign breast disease. Emerging evidence for aluminium in several breast structures now requires biomarkers of aluminium action in order to ascertain whether the presence of aluminium has any biological impact. To this end, we report raised levels of proteins that modulate iron homeostasis (ferritin, transferrin) in parallel with raised aluminium in nipple aspirate fluids in vivo, and we report overexpression of mRNA for several S100 calcium binding proteins following long-term exposure of MCF-7 human breast cancer cells in vitro to aluminium chlorhydrate.

Arabidopsis annexin1 mediates the radical-activated plasma membrane Ca2+ - and K+ -permeable conductance in root cells

Plant cell growth and stress signaling require Ca2+ influx through plasma membrane transport proteins that are regulated by reactive oxygen species. In root cell growth, adaptation to salinity stress, and stomatal closure, such proteins operate downstream of the plasma membrane NADPH oxidases that produce extracellular superoxide anion, a reactive oxygen species that is readily converted to extracellular hydrogen peroxide and hydroxyl radicals, OH_. In root cells, extracellular OH_ activates a plasma membrane Ca2+-permeable conductance that permits Ca2+ influx. In Arabidopsis thaliana, distribution of this conductance resembles that of annexin1 (ANN1). Annexins are membrane binding proteins that can form Ca2+-permeable conductances in vitro. Here, the Arabidopsis loss-of-function mutant for annexin1 (Atann1) was found to lack the root hair and epidermal OH_-activated Ca2+- and K+-permeable conductance. This manifests in both impaired root cell growth and ability to elevate root cell cytosolic free Ca2+ in response to OH_. An OH_-activated Ca2+ conductance is reconstituted by recombinant ANN1 in planar lipid bilayers. ANN1 therefore presents as a novel Ca2+-permeable transporter providing a molecular link between reactive oxygen species and cytosolic Ca2+ in plants.

Regulation of mitogen-activated protein kinases in cardiac myocytes through the small G protein Rac1.

Small guanine nucleotide-binding proteins of the Ras and Rho (Rac, Cdc42, and Rho) families have been implicated in cardiac myocyte hypertrophy, and this may involve the extracellular signal-related kinase (ERK), c-Jun N-terminal kinase (JNK), and/or p38 mitogen-activated protein kinase (MAPK) cascades. In other systems, Rac and Cdc42 have been particularly implicated in the activation of JNKs and p38-MAPKs. We examined the activation of Rho family small G proteins and the regulation of MAPKs through Rac1 in cardiac myocytes. Endothelin 1 and phenylephrine (both hypertrophic agonists) induced rapid activation of endogenous Rac1, and endothelin 1 also promoted significant activation of RhoA. Toxin B (which inactivates Rho family proteins) attenuated the activation of JNKs by hyperosmotic shock or endothelin 1 but had no effect on p38-MAPK activation. Toxin B also inhibited the activation of the ERK cascade by these stimuli. In transfection experiments, dominant-negative N17Rac1 inhibited activation of ERK by endothelin 1, whereas activated V12Rac1 cooperated with c-Raf to activate ERK. Rac1 may stimulate the ERK cascade either by promoting the phosphorylation of c-Raf or by increasing MEK1 and/or -2 association with c-Raf to facilitate MEK1 and/or -2 activation. In cardiac myocytes, toxin B attenuated c-Raf(Ser-338) phosphorylation (50 to 70% inhibition), but this had no effect on c-Raf activity. However, toxin B decreased both the association of MEK1 and/or -2 with c-Raf and c-Raf-associated ERK-activating activity. V12Rac1 cooperated with c-Raf to increase expression of atrial natriuretic factor (ANF), whereas N17Rac1 inhibited endothelin 1-stimulated ANF expression, indicating that the synergy between Rac1 and c-Raf is potentially physiologically important. We conclude that activation of Rac1 by hypertrophic stimuli contributes to the hypertrophic response by modulating the ERK and/or possibly the JNK (but not the p38-MAPK) cascades.