Effect of photoinhibition on algal photosynthesis: a dynamic model

Recent evidence from algal physiology and molecular biologyconfirms that photoinhibition is directly related to D1 proteindamage and recovery, and D1 protein damage leads to a decreasein electron transfer or an increase in turnover time of theelectron transfer chain. In this study, the turnover time ofthe electron transfer chain is defined as a function of therelative concentration of D1 protein in reaction centre II andthe photoinhibition processes due to D1 protein degradationare incorporated into a model of photosynthesis, initiated byDubinsky et al. (Plant Cell Physiol., 27, 1335–1349, 1986)and developed by Sakshaug et al. (Limnol. Oceanogr., 34, 198–205,1989). D1 protein damage is assumed to be both light and D1protein concentration dependent, and to be proportional to thecross-section of PSII (     

From pull-down data to protein interaction networks and complexes with biological relevance

Motivation: Recent improvements in high-throughput Mass Spectrometry(MS) technology have expedited genome-wide discovery of protein–proteininteractions by providing a capability of detecting proteincomplexes in a physiological setting. Computational inferenceof protein interaction networks and protein complexes from MSdata are challenging. Advances are required in developing robustand seamlessly integrated procedures for assessment of protein–proteininteraction affinities, mathematical representation of proteininteraction networks, discovery of protein complexes and evaluationof their biological relevance. Results: A multi-step but easy-to-follow framework for identifyingprotein complexes from MS pull-down data is introduced. It assessesinteraction affinity between two proteins based on similarityof their co-purification patterns derived from MS data. It constructsa protein interaction network by adopting a knowledge-guidedthreshold selection method. Based on the network, it identifiesprotein complexes and infers their core components using a graph-theoreticalapproach. It deploys a statistical evaluation procedure to assessbiological relevance of each found complex. On Saccharomycescerevisiae pull-down data, the framework outperformed othermore complicated schemes by at least 10% in F₁-measure and identified610 protein complexes with high-functional homogeneity basedon the enrichment in Gene Ontology (GO) annotation. Manual examinationof the complexes brought forward the hypotheses on cause offalse identifications. Namely, co-purification of differentprotein complexes as mediated by a common non-protein molecule,such as DNA, might be a source of false positives. Protein identificationbias in pull-down technology, such as the hydrophilic bias couldresult in false negatives. Contact: samatovan{at}ornl.gov Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Jonathan Wren Present address: Department of Biomedical Informatics, VanderbiltUniversity, Nashville, TN 37232. The authors wish it to be known that, in their opinion, thefirst two authors should be regarded as joint First Authors.     

Celestial3D: a novel method for 3D visualization of familial data

Summary: Traditional two-dimensional (2D) software programsfor drawing pedigrees are limited when dealing with extendedpedigrees. In successive generations, the number of individualsgrows exponentially, leading to an unworkable amount of spacerequired in the horizontal direction for 2D displays. In addition,it is not always possible to place closely related individualsnear each other due to the lack of space in 2Ds. To addressthese issues we have developed three-dimensional (3D) pedigreedrawing techniques to enable clearer visualization of extendedpedigrees. Currently no other methods are available for displayingextended pedigrees in 3Ds. We have made freely available a softwaretool—‘Celestial3D’—that implements thesenovel techniques. Availability: Freely available to non-commercial users Contact: celestial3d{at}genepi.org.au Supplementary information: www.genepi.org.au/celestial3d Associate Editor: Martin Bishop ¹A more extensive list of software tools appears in the SupplementaryMaterial.     

A combined molecular modeling study on a series of pyrazole/isoxazole based human Hsp90α inhibitors

Inhibition of the protein chaperone Hsp90α is a promising approach for cancer therapy. In this work, a molecular modeling study combining pharmacophore model, molecular docking and three-dimensional quantitative structure-activity relationships (3D-QSAR) was performed to investigate a series of pyrazole/isoxazole scaffold inhibitors of human Hsp90α. The pharmacophore model can provide the essential features required for the biological activities of the inhibitors. The molecular docking study can give insight into the binding mode between Hsp90α and its inhibitors. 3D-QSAR based on CoMFA and CoMSIA models were performed from three different strategies for conformational selection and alignment. The receptor-based models gave the most statistically significant results with cross-validated q ² values of 0.782 and 0.829 and r ² values of 0.909 and 0.968, for CoMFA and CoMSIA respectively. Furthermore, the 3D contour maps superimposed within the binding site of Hsp90α could help to understand the pivotal interaction and the structural requirements for potent Hsp90α inhibitors. The results show 4-position of pyrazole/isoxazole ring requires bulky and hydrophobic groups, and bulky and electron repulsion substituent of 5-amides is favorable for enhancing activity. This study will be helpful for the rational design of new potent Hsp90α inhibitors.     

3D-Garden: a system for modelling protein-protein complexes based on conformational refinement of ensembles generated with the marching cubes algorithm

Motivation: Reliable structural modelling of protein–proteincomplexes has widespread application, from drug design to advancingour knowledge of protein interactions and function. This workaddresses three important issues in protein–protein docking:implementing backbone flexibility, incorporating prior indicationsfrom experiment and bioinformatics, and providing public accessvia a server. 3D-Garden (Global And Restrained Docking ExplorationNexus), our benchmarked and server-ready flexible docking system,allows sophisticated programming of surface patches by the uservia a facet representation of the interactors’ molecularsurfaces (generated with the marching cubes algorithm). Flexibilityis implemented as a weighted exhaustive conformer search foreach clashing pair of molecular branches in a set of 5000 modelsfiltered from around 340 000 initially. Results: In a non-global assessment, carried out strictly accordingto the protocols for number of models considered and model qualityof the Critical Assessment of Protein Interactions (CAPRI) experiment,over the widely-used Benchmark 2.0 of 84 complexes, 3D-Gardenidentifies a set of ten models containing an acceptable or bettermodel in 29/45 test cases, including one with large conformationalchange. In 19/45 cases an acceptable or better model is rankedfirst or second out of 340 000 candidates. Availability: http://www.sbg.bio.ic.ac.uk/3dgarden (server) Contact: v.lesk{at}ic.ac.uk Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Burkhard Rost     

Segmentation of cDNA microarray spots using markov random field modeling

Motivation: Spot segmentation is a critical step in microarraygene expression data analysis. Therefore, the performance ofsegmentation may substantially affect the results of subsequentstages of the analysis, such as the detection of differentiallyexpressed genes. Several methods have been developed to segmentmicroarray spots from the surrounding background. In this study,we have proposed a new approach based on Markov random field(MRF) modeling and tested its performance on simulated and realmicroarray images against a widely used segmentation methodbased on Mann–Whitney test adopted by QuantArray software(Boston, MA). Spot addressing was performed using QuantArray.We have also devised a simulation method to generate microarrayimages with realistic features. Such images can be used as goldstandards for the purposes of testing and comparing differentsegmentation methods, and optimizing segmentation parameters. Results: Experiments on simulated and 14 actual microarray imagesets show that the proposed MRF-based segmentation method candetect spot areas and estimate spot intensities with higheraccuracy. Availability: The algorithms were implemented in Matlab^TM (TheMathworks, Inc., Natick, MA) environment. The codes for MRF-basedsegmentation and image simulation methods are available uponrequest. Contact: demirkaya{at}ieee.org     

Effects of Blue Light Pretreatments on Internode Extension Growth in Mustard Seedlings after the Transition to Darkness: Analysis of the Interaction with Phytochrome

The effects of blue light (B) pretreatments on internode extensiongrowth and their possible interaction with phytochrome mediatedresponses were examined in Sinapis alba seedlings grown for11 d under 280 µmol m^–2 s^–1 of continuousblue-deficient light from low pressure sodium lamps (SOX). SupplementaryB (16 µmol m^–2 s^–1) caused no detectable inhibitionof the first internode growth rate under continuous SOX, butgrowth rate was inhibited after transfer to darkness. This effect,and the growth promotion caused by far-red bend-of-day' lightpulses were additive. The addition of B at 16 µmol m^–2s^–1 during 11 d, or only during the first 9 or 10 d orthe latest 0.75, 1 or 2 d of the SOX pretreatment caused approximatelythe same extent of inhibition after the transition to darkness.A single hour of supplementary B before darkness caused morethan 50% of the maximum inhibition. However, 24 h of lower fluencerates of B (4 or 7 µmol m^–2 s^–1) were ineffective.Covering the internode during the supplementary B period didnot prevent the response to B after the transition to darkness.Far-red light given simultaneously with B (instead of the SOXbackground) reduced the inhibitory effect of B. Above a given threshold fluence rate, B perceived mainly inthe leaves inhibits extension growth in subsequent darkness,provided that high phytochrome photo-equilibria are presentduring the irradiation with B. Once triggered, this effect doesnot interact significantly with the ‘end-of-day’phytochrome effect. Key words: Blue light, extension growth, phytochrome     

Pharmacophore based 3D-QSAR modeling and free energy analysis of VEGFR-2 inhibitors

VEGFR-2, a transmembrane tyrosine kinase receptor is responsible for angiogenesis and has been an attractive target in treating cancers. The inhibition mechanism of structurally diverse urea derivatives, reported as VEGFR-2 inhibitors, was explored by pharmacophore modeling, QSAR, and molecular dynamics based free energy analysis.The pharmacophore hypothesis AADRR, resulted in a highly significant atom based 3D-QSAR model (r² = 0.94 and q² = 0.84). Binding free energy analysis of the docked complexes of highly active and inactive compounds, after 7 ns MD simulation, revealed the importance of van der Waals interaction in VEGFR-2 inhibition. The decomposition of binding free energy on a per residue basis disclosed that the residues in hinge region and hydrophobic pocket play a role in discriminating the active and inactive inhibitors. Thus, the present study proposes a pharmacophore hypothesis representing the identified interactions pattern and its further application as a template in screening databases to identify novel VEGFR-2 inhibitor scaffolds.     

Automated image alignment for 2D gel electrophoresis in a high-throughput proteomics pipeline

Motivation: The quest for high-throughput proteomics has revealeda number of challenges in recent years. Whilst substantial improvementsin automated protein separation with liquid chromatography andmass spectrometry (LC/MS), aka ‘shotgun’ proteomics,have been achieved, large-scale open initiatives such as theHuman Proteome Organization (HUPO) Brain Proteome Project haveshown that maximal proteome coverage is only possible when LC/MSis complemented by 2D gel electrophoresis (2-DE) studies. Moreover,both separation methods require automated alignment and differentialanalysis to relieve the bioinformatics bottleneck and so makehigh-throughput protein biomarker discovery a reality. The purposeof this article is to describe a fully automatic image alignmentframework for the integration of 2-DE into a high-throughputdifferential expression proteomics pipeline. Results: The proposed method is based on robust automated imagenormalization (RAIN) to circumvent the drawbacks of traditionalapproaches. These use symbolic representation at the very earlystages of the analysis, which introduces persistent errors dueto inaccuracies in modelling and alignment. In RAIN, a third-ordervolume-invariant B-spline model is incorporated into a multi-resolutionschema to correct for geometric and expression inhomogeneityat multiple scales. The normalized images can then be compareddirectly in the image domain for quantitative differential analysis.Through evaluation against an existing state-of-the-art methodon real and synthetically warped 2D gels, the proposed analysisframework demonstrates substantial improvements in matchingaccuracy and differential sensitivity. High-throughput analysisis established through an accelerated GPGPU (general purposecomputation on graphics cards) implementation. Availability: Supplementary material, software and images usedin the validation are available at http://www.proteomegrid.org/rain/ Contact: g.z.yang{at}imperial.ac.uk Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: David Rocke     

The Hydrostatic and Hydrodynamic Volumes of Polyols in Aqueous Solutions and Their Sweet Taste

The tastes and solution properties of sugar alcohols were studiedin an attempt to illuminate the mechanism of sweet taste chemoreception.The SMURF method was used to measure taste time-intensity ofaqueous solutions of sugar alcohols and the results were interpretedusing the Stevens power function and kinetic parameters. Theapparent molar volumes, apparent specific volumes, partial molarvolumes, partial specific volumes and intrinsic viscositiesof the solutions were studied. Apparent molar volume reflectsthe size of the molecule in a hydrostatic state whereas intrinsicviscosity gives a measure of the size of the molecules in ahydrodynamic state. Generally the apparent molar volumes ofthe polyols are 6–13% greater than those of the parentsugars, indicating less interaction with the water structure.Apparent specific volume values can predict taste quality, andthe average apparent specific volume for the sugar alcoholsstudied fits within the central part of the sweet range, i.e.0.5–0.68 cm³/g, which accords with their ability to elicita pure sweet taste response. Intensities and persistences ofsweetness in the polyols followed the same trend as intrinsicviscosities. Chem. Senses 22: 149–161, 1997.     

Polycyclic dinitriles: a novel class of potent GABAergic insecticides provides a new radioligand, [3H]BIDN

The polycyclic dinitriles are a potent class of insecticides which are non-competitive GABA (γ-aminobutyric a acid) antagonists acting at the convulsant site. Comparison with other classes of GABA convulsant site ligands using molecular modelling has shown significant structural similarities. We have developed a pharmacophore model which unifies this class and some previous classes of GABA convulsants. Key pharmacophore elements are a polarizable functionality separated by a fixed distance from two H-bond accepting elements. This model is based on information from X-ray crystal structures and Sybyl using the Tripos force field. Using this pharmacophore model, numerous structural modifications were explored to enhance understanding of structure-activity relationships at the GABA receptor convulsant site of insects and mammals. A radiolabelled bicyclic dinitrile, [³H]BIDN ([³H]3,3-bis-trifluoromethyl-bicyclo[2,2,1] heptane-2,2-dicarbonitrile), was prepared from this area of chemistry and was used as a probe for the interaction of polycyclic dinitriles at the target site.     

Practical FDR-based sample size calculations in microarray experiments

Motivation: Owing to the experimental cost and difficulty inobtaining biological materials, it is essential to considerappropriate sample sizes in microarray studies. With the growinguse of the False Discovery Rate (FDR) in microarray analysis,an FDR-based sample size calculation is essential. Method: We describe an approach to explicitly connect the samplesize to the FDR and the number of differentially expressed genesto be detected. The method fits parametric models for degreeof differential expression using the Expectation–Maximizationalgorithm. Results: The applicability of the method is illustrated withsimulations and studies of a lung microarray dataset. We proposeto use a small training set or published data from relevantbiological settings to calculate the sample size of an experiment. Availability: Code to implement the method in the statisticalpackage R is available from the authors. Contact: jhu{at}mdanderson.org     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Influence of the 3'-UTR-length of mKIAA cDNAs and their Sequence Features to the mRNA Expression Level in the Brain

We have previously described the sequence features of 1500 mouseKIAA (mKIAA) genes in comparison with those of human KIAA genes(Okazaki, N., Kikuno, R., Inamoto, S., Hara, Y., Nagase, T.,Ohara, O., and Koga, H. 2002, DNA Res., 9, 179–188; Okazaki,N., Kikuno, R., Ohara, R., Inamoto, S., Aizawa, H., Yuasa, S.,Nakajima, D., Nagase, T., Ohara, O., and Koga, H. 2003, DNARes., 10, 35–48; Okazaki, N., Kikuno, R., Ohara, R., Inamoto,S., Koseki, H., Hiraoka, S., Saga, Y., Nagase, T., Ohara, O.,and Koga, H. 2003, DNA Res., 10, 167–180; and Okazaki,N., F-Kikuno, R., Ohara, R., Inamoto, S., Koseki, H., Hiraoka,S., Saga, Y., Seino, S., Nishimura, M., Kaisho, T., Hoshino,K., Kitamura, H., Nagase, T., Ohara, O., and Koga, H. 2004,DNA Res., 11, 205–218). To validate the orthologous relationshipbetween mKIAA and KIAA genes in detail, we examined their chromosomalpositions and evolutionary rate of synonymous substitutionsand confirmed that >93% of the mKIAA/KIAA gene pairs areorthologous. During the sequence analysis of mKIAA genes, wefound that 3'-untranslated region (3'-UTR) lengths of mKIAAand KIAA genes are extremely long. In the meanwhile, we havealso examined the tissue-specific expression of 1700 mKIAA genesusing cDNA microarray and verified predominantly their expressionin adult brain (Koga, H., Yuasa, S., Nagase, T., Shimada, K.,Nagano, M., Imai, K., Ohara, R., Nakajima, D., Murakami, M.,Kawai, M., Miki, F., Magae, J., Inamoto, S., Okazaki, N., Ohara,O. 2004, DNA Res., 11, 293–304). To connect these twoevidences, we statistically analysed the relationship betweenthem by using the mKIAA genes. Consequently, a positive correlationwas observed between the 3'-UTR lengths and the relative expressionintensities in adult brain. Furthermore, we searched sequenceelements in the 3'-UTR possibly related with their expressionand found some candidates regarding the brain-specific expression.     

Genetic Sensitivity to 6-n-Propylthiouracil (PROP) and Hedonic Responses to Bitter and Sweet Tastes

Genetically mediated sensitivity to the bitter taste of 6-n-propylthiouracil(PROP) has been associated with greater acuity for bitter andfor some sweet tastes. Thus far, few studies have explored therelationship between PROP taste sensitivity and hedonic responsesto bitter and sweet. In this study, 87 normal-weight young womenwere divided into PROP non-tasters (n = 18), regular tasters(n = 49), and supertasters (n = 20), based on their PROP detectionthresholds and the scaling of five suprathreshold solutionsof PROP and NaCl. Non-tasters had thresholds >1.8 x 10^–4mol/l PROP. Supertasters had thresholds <3.2 x 10^–5mol/l PROP and PROP/NaCl ratios >1.70. As expected, dislikeof the bitter taste of PROP was determined by its perceivedintensity, which was greater among supertasters than among regulartasters or non-tasters. Significant correlations were observedbetween PROP taste thresholds and the sum of intensity ratings(r = –0.61) and between summed intensity and summed hedonicratings (r = –0.80). PROP taste sensitivity was weaklylinked to enhanced perception of sweet taste, but did not predicthedonic responses to sucrose or to saccharin solutions. Giventhat the dislike of PROP solutions is determined by their perceivedintensity, hedonic responses to PROP solutions may provide arapid way of screening for PROP taster status. Chem. Senses22: 27–37, 1997.     

In silico approaches to identify novel myeloid cell leukemia-1 (Mcl-1) inhibitors for treatment of cancer

Myeloid cell leukemia-1 (Mcl-1) has been a validated and attractive target for cancer therapy. Over-expression of Mcl-1 in many cancers allows cancer cells to evade apoptosis and contributes to the resistance to current chemotherapeutics. Here, we identified new Mcl-1 inhibitors using a multi-step virtual screening approach. First, based on two different ligand-receptor complexes, 20 pharmacophore models were established by simultaneously using ‘Receptor-Ligand Pharmacophore Generation’ method and manual build feature method, and then carefully validated by a test database. Then, pharmacophore-based virtual screening (PB-VS) could be performed by using the 20 pharmacophore models. In addition, docking study was used to predict the possible binding poses of compounds, and the docking parameters were optimized before performing docking-based virtual screening (DB-VS). Moreover, a 3D QSAR model was established by applying the 55 aligned Mcl-1 inhibitors. The 55 inhibitors sharing the same scaffold were docked into the Mcl-1 active site before alignment, then the inhibitors with possible binding conformations were aligned. For the training set, the 3D QSAR model gave a correlation coefficient r² of 0.996; for the test set, the correlation coefficient r² was 0.812. Therefore, the developed 3D QSAR model was a good model, which could be applied for carrying out 3D QSAR-based virtual screening (QSARD-VS). After the above three virtual screening methods orderly filtering, 23 potential inhibitors with novel scaffolds were identified. Furthermore, we have discussed in detail the mapping results of two potent compounds onto pharmacophore models, 3D QSAR model, and the interactions between the compounds and active site residues.     

Expanding Southwest Pacific mitochondrial haplogroups P and Q

doi:10.1093/molbev/msi142 Mol Biol Evol. 22:1506–1517. 2005. In the report of mitochondrial DNA haplogroup variation in NearOceania, the section on the calculation     

Pharmacophore mapping of arylamino-substituted benzo[b]thiophenes as free radical scavengers

Predictive pharmacophore models have been developed for a series of arylamino-substituted benzo[b]thiophenes exhibiting free radical scavenging activity. 3D pharmacophore models were generated using a set of 20 training set compounds and subsequently validated by mapping 6 test set compounds using Discovery Studio 2.1 software. Further model validation was performed by randomizing the data using Fischer’s validation technique at the 95% confidence level. The most predictive pharmacophore model developed using the conformers obtained from the BEST method showed a correlation coefficient (r) of 0.942 and consisted of three features: hydrogen bond donor, hydrogen bond acceptor and aromatic ring. Acceptable values of external validation parameters, like R_pred² R_{{\rm{pred}}}^2 (0.853) and r_{m( test )}² r_{m\left( {test} \right)}^2 (0.844), also implied that the external predictivity of the model was significant. The development of further pharmacophore models using conformers obtained from the FAST method yielded a few models with good predictivity, with the best one (r = 0.904) consisting of two features: hydrogen bond donor and hydrogen bond acceptor. Significant values of external validation parameters, R_pred² R_{{\rm{pred}}}^2 (0.913) and r_{m( test )}² r_{m\left( {test} \right)}^2 (0.821), also reflect the high predictive ability of the model. Again, Fischer validation results implied that the models developed were robust enough and their good results were not based on mere chance. These validation approaches indicate the reliability of the predictive abilities of the 3D pharmacophore models developed here, which may thus be further utilized as a 3D query tool in the virtual screening of new chemical entities with potent antioxidant activities.