Testing statistical significance scores of sequence comparison methods with structure similarity

Background  

In the past years the Smith-Waterman sequence comparison algorithm has gained popularity due to improved implementations and rapidly increasing computing power. However, the quality and sensitivity of a database search is not only determined by the algorithm but also by the statistical significance testing for an alignment. The e-value is the most commonly used statistical validation method for sequence database searching. The CluSTr database and the Protein World database have been created using an alternative statistical significance test: a Z-score based on Monte-Carlo statistics. Several papers have described the superiority of the Z-score as compared to the e-value, using simulated data. We were interested if this could be validated when applied to existing, evolutionary related protein sequences.     

Identifying dynamical modules from genetic regulatory systems: applications to the segment polarity network

Background  

It is widely accepted that genetic regulatory systems are 'modular', in that the whole system is made up of smaller 'subsystems' corresponding to specific biological functions. Most attempts to identify modules in genetic regulatory systems have relied on the topology of the underlying network. However, it is the temporal activity (dynamics) of genes and proteins that corresponds to biological functions, and hence it is dynamics that we focus on here for identifying subsystems.     

Prodepth: Predict Residue Depth by Support Vector Regression Approach from Protein Sequences Only

Residue depth (RD) is a solvent exposure measure that complements the information provided by conventional accessible surface area (ASA) and describes to what extent a residue is buried in the protein structure space. Previous studies have established that RD is correlated with several protein properties, such as protein stability, residue conservation and amino acid types. Accurate prediction of RD has many potentially important applications in the field of structural bioinformatics, for example, facilitating the identification of functionally important residues, or residues in the folding nucleus, or enzyme active sites from sequence information. In this work, we introduce an efficient approach that uses support vector regression to quantify the relationship between RD and protein sequence. We systematically investigated eight different sequence encoding schemes including both local and global sequence characteristics and examined their respective prediction performances. For the objective evaluation of our approach, we used 5-fold cross-validation to assess the prediction accuracies and showed that the overall best performance could be achieved with a correlation coefficient (CC) of 0.71 between the observed and predicted RD values and a root mean square error (RMSE) of 1.74, after incorporating the relevant multiple sequence features. The results suggest that residue depth could be reliably predicted solely from protein primary sequences: local sequence environments are the major determinants, while global sequence features could influence the prediction performance marginally. We highlight two examples as a comparison in order to illustrate the applicability of this approach. We also discuss the potential implications of this new structural parameter in the field of protein structure prediction and homology modeling. This method might prove to be a powerful tool for sequence analysis.     

Hybrid approaches to molecular simulation

Molecular dynamics (MD) simulation is an established method for studying the conformational changes that are important for protein function. Recent advances in hardware and software have allowed MD simulations over the same timescales as experiment, improving the agreement between theory and experiment to a large extent. However, running such simulations are costly, in terms of resources, storage, and trajectory analysis. There is still a place for techniques that involve short MD simulations. In order to overcome the sampling paucity of short time-scales, hybrid methods that include some form of MD simulation can exploit certain features of the system of interest, often combining experimental information in surprising ways. Here, we review some recent hybrid approaches to the simulation of proteins.     

Telomere-related sequences at interstitial sites in the human genome

The ends (telomeres) of eukaryotic chromosomes are protected from degradation and from loss during DNA replication by buffers of simple tandem repetitive sequence. The nucleotide sequence of these telomeric arrays is fundamental to telomere function as a site for protein and ribonucleoprotein binding and varies only slightly in a wide range of organisms. We present evidence that arrays of this human telomeric sequence, TTAGGG, are present not only at the ends of human chromosomes but also at numerous interstitial sites. These interstitial loci share nucleotide sequence similarity outside the repetitive array, suggesting that they are related functionally or have evolved from a common progenitor locus.     

One-step purification of the serotonin transporter located at the human platelet plasma membrane.

A 68-kDa glycoprotein bearing the biological activity of the plasma membrane serotonin (5-hydroxytryptamine, 5-HT) transporter has been purified from human blood platelets, a classical cell model for the study of 5-HT uptake. After treatment of the whole platelet population or its plasma membrane fraction by sulfhydryl-dependent bacterial protein toxins or by digitonin, purification was reproducibly obtained by a one-step affinity chromatography using two different columns with 5-HT or 6-fluorotryptamine as ligands and elution by 5-HT or Na(+)-free buffer. The purified fraction migrated on sodium dodecyl sulfate-polyacrylamide gel electrophoresis as a single band with an apparent molecular mass of 68 kDa and exhibited an apparent isoelectric point of 5.6-6.2. Two sialic acid residues were detected in the purified material. The purified glycoprotein bound the 5-HT uptake blocker [3H]paroxetine with a Kd (0.25 nM) similar to the one observed for intact human platelets. It also bound [3H] 5-HT but neither [3H]hydroxytetrabenazine nor [3H] ouabain, the respective markers of the granular monoamine transporter and of the Na+,K(+)-ATPase associated to the plasma membrane 5-HT transporter. 5-HT derivatives and 5-HT uptake inhibitors exhibited similar Ki values for 5-HT uptake and paroxetine binding in intact human platelets and in the purified glycoprotein. Under laser UV irradiation, 40% of this purified glycoprotein could be labeled by either [3H]paroxetine or [3H]cyanoimipramine. No labeling was detected with either [3H] gamma-aminobutyric acid or [3H]GBR 12783, the respective markers of gamma-aminobutyric acid and dopamine carriers. The purified 68-kDa protein is therefore likely to correspond at least to the binding domain of the 5-HT transporter located at the human platelet plasma membrane.     

Study of bacteriophage T4-encoded Dam DNA (adenine-N6)-methyltransferase binding with substrates by rapid laser UV cross-linking

DNA methyltransferases of the Dam family (including bacteriophage T4-encoded Dam DNA (adenine-N(6))-methyltransferase (T4Dam)) catalyze methyl group transfer from S-adenosyl-L-methionine (AdoMet), producing S-adenosyl-L-homocysteine (AdoHcy) and methylated adenine residues in palindromic GATC sequences. In this study, we describe the application of direct (i.e. no exogenous cross-linking reagents) laser UV cross-linking as a universal non-perturbing approach for studying the characteristics of T4Dam binding with substrates in the equilibrium and transient modes of interaction. UV irradiation of the enzyme.substrate complexes using an Nd(3+):yttrium aluminum garnet laser at 266 nm resulted in up to 3 and >15% yields of direct T4Dam cross-linking to DNA and AdoMet, respectively. Consequently, we were able to measure equilibrium constants and dissociation rates for enzyme.substrate complexes. In particular, we demonstrate that both reaction substrates, specific DNA and AdoMet (or product AdoHcy), stabilized the ternary complex. The improved substrate affinity for the enzyme in the ternary complex significantly reduced dissociation rates (up to 2 orders of magnitude). Several of the parameters obtained (such as dissociation rate constants for the binary T4Dam.AdoMet complex and for enzyme complexes with a nonfluorescent hemimethylated DNA duplex) were previously inaccessible by other means. However, where possible, the results of laser UV cross-linking were compared with those of fluorescence analysis. Our study suggests that rapid laser UV cross-linking efficiently complements standard DNA methyltransferase-related tools and is a method of choice to probe enzyme-substrate interactions in cases in which data cannot be acquired by other means.     

Modelling of Thyroid Peroxidase Reveals Insights into Its Enzyme Function and Autoantigenicity

Thyroid peroxidase (TPO) catalyses the biosynthesis of thyroid hormones and is a major autoantigen in Hashimoto’s disease—the most common organ-specific autoimmune disease. Epitope mapping studies have shown that the autoimmune response to TPO is directed mainly at two surface regions on the molecule: immunodominant regions A and B (IDR-A, and IDR-B). TPO has been a major target for structural studies for over 20 years; however, to date, the structure of TPO remains to be determined. We have used a molecular modelling approach to investigate plausible modes of TPO structure and dimer organisation. Sequence features of the C-terminus are consistent with a coiled-coil dimerization motif that most likely anchors the TPO dimer in the apical membrane of thyroid follicular cells. Two contrasting models of TPO were produced, differing in the orientation and exposure of their active sites relative to the membrane. Both models are equally plausible based upon the known enzymatic function of TPO. The “trans” model places IDR-B on the membrane-facing side of the myeloperoxidase (MPO)-like domain, potentially hindering access of autoantibodies, necessitating considerable conformational change, and perhaps even dissociation of the dimer into monomers. IDR-A spans MPO- and CCP-like domains and is relatively fragmented compared to IDR-B, therefore most likely requiring domain rearrangements in order to coalesce into one compact epitope. Less epitope fragmentation and higher solvent accessibility of the “cis” model favours it slightly over the “trans” model. Here, IDR-B clusters towards the surface of the MPO-like domain facing the thyroid follicular lumen preventing steric hindrance of autoantibodies. However, conformational rearrangements may still be necessary to allow full engagement with autoantibodies, with IDR-B on both models being close to the dimer interface. Taken together, the modelling highlights the need to consider the oligomeric state of TPO, its conformational properties, and its proximity to the membrane, when interpreting epitope-mapping data.