Strong patterns in homooligomer tracts occurrences in non-coding and in potential regulatory sites in eukaryotic genomes

Previous studies of the dinucleotides flanking both the 5' and 3' ends of homooligomer tracts have shown that some flanks are consistently preferred over others (1,2). In the first preferred group, the homooligomer tracts are flanked by the same nucleotide and/or the complementary nucleotides, e.g.,ATAn,TTAn,CCGn, where n = 2-5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group An/Tn are flanked by C and/or G; Gn/Cn are flanked by A/T, e.g.,CGAn,TnGG,GnAT). The frequencies of runs flanked by A or T, and G or C ("mixed"group) are as expected. Here we seek the origin of this effect and its relevance to protein-DNA interactions. Surprisingly, within the first group, runs flanked by their complements with a pyrimidine-purine junction (e.g.,TTAn,CnGG) are greatly preferred. The frequencies of their purine-pyrimidine junction mirror-images is just as expected. This effect, as well as additional ones enumerated below, is seen universally in eukaryotes and in prokaryotes, although it is stronger in the former. Detailed analysis of regulatory regions shows these strong trends, particularly in GC sequences. The potential relationship to DNA conformation and DNA-protein interaction is discussed.     

Effect of intravenously-administered putative and potential antagonists of ethanol on sleep time in ethanol-narcotized mice

Groups of male CD-1 mice (n = 12/group) were injected intraperitoneally (IP) with 5 g ethanol/kg of body weight. After loss of righting reflex, they were given vehicle or one of 2-3 doses of reputed or potential antagonists of ethanol intravenously (IV). Sleep time was measured from loss to return of righting reflex. Mean sleep time (MST) was increased significantly (P less than 0.05) by a large dose of dl-amphetamine (24 mg/kg) and by 4-aminopyridine (1, 5 mg/kg). Significant (P less than 0.01) increases were also produced by small and large doses of aminophylline (25, 100 mg/kg) and by yohimbine (1, 5 mg/kg). MST was not altered significantly by small and medium doses of dl-amphetamine (6, 12 mg/kg), a medium dose of aminophylline (50 mg/kg), or by any doses of naloxone, thyrotropin-releasing hormone, propranolol, physostigmine, doxapram, or Ro 15-4513. When Ro 15-4513 was given IP 15 minutes before ethanol (n = 6/group), onset and duration of narcosis were not altered. None of the compounds tested was an effective IV antidote for deep ethanol narcosis because of drug side effects, toxicity, prolongation of MST, or insufficient shortening of MST.     

Potentials 'R'Us web-server for protein energy estimations with coarse-grained knowledge-based potentials

Background  

Knowledge-based potentials have been widely used in the last 20 years for fold recognition, protein structure prediction from amino acid sequence, ligand binding, protein design, and many other purposes. However generally these are not readily accessible online.     

Structural interpretation of protein-protein interaction network

Background

Currently a huge amount of protein-protein interaction data is available from high throughput experimental methods. In a large network of protein-protein interactions, groups of proteins can be identified as functional clusters having related functions where a single protein can occur in multiple clusters. However experimental methods are error-prone and thus the interactions in a functional cluster may include false positives or there may be unreported interactions. Therefore correctly identifying a functional cluster of proteins requires the knowledge of whether any two proteins in a cluster interact, whether an interaction can exclude other interactions, or how strong the affinity between two interacting proteins is.

Methods

In the present work the yeast protein-protein interaction network is clustered using a spectral clustering method proposed by us in 2006 and the individual clusters are investigated for functional relationships among the member proteins. 3D structural models of the proteins in one cluster have been built – the protein structures are retrieved from the Protein Data Bank or predicted using a comparative modeling approach. A rigid body protein docking method (Cluspro) is used to predict the protein-protein interaction complexes. Binding sites of the docked complexes are characterized by their buried surface areas in the docked complexes, as a measure of the strength of an interaction.

Results

The clustering method yields functionally coherent clusters. Some of the interactions in a cluster exclude other interactions because of shared binding sites. New interactions among the interacting proteins are uncovered, and thus higher order protein complexes in the cluster are proposed. Also the relative stability of each of the protein complexes in the cluster is reported.

Conclusions

Although the methods used are computationally expensive and require human intervention and judgment, they can identify the interactions that could occur together or ones that are mutually exclusive. In addition indirect interactions through another intermediate protein can be identified. These theoretical predictions might be useful for crystallographers to select targets for the X-ray crystallographic determination of protein complexes.

     

GFP-tagged PrP supports compromised prion replication in transgenic mice

The ability of green fluorescent protein (GFP)-prion protein (PrP) fusions to support prion propagation has not been demonstrated. Here, we show that while transgenic mice expressing PrP tagged at its amino terminus with enhanced GFP, referred to as EGFPrP-N, supported prion replication, disease onset was prolonged, the brains of diseased mice did not exhibit typical disease neuropathology and disease-associated EGFPrP-N lacked the full spectrum of biochemical properties normally associated with PrP(Sc). Co-expression of wild-type PrP and EGFPrP-N substantially reduced prion incubation times and resulted in accumulation of protease-resistant EGFPrP(Sc)-N in the brains of transgenic mice as well as chronically infected cultured cells, suggesting that wild-type PrP rescued a compromised amino terminal function in EGFPrP-N. While our results show that EGFPrP(C)-N adopts a conformation necessary for the production of infectious prions, the synergistic interaction of wild-type and EGFPrP-N underscores the importance of the amino terminus in modulating prion pathogenesis.     

A comparison of the effects of prenatal exposure of CD‐1 mice to three imidazolium‐based ionic liquids

BACKGROUND: Ionic liquids (ILs; salts with melting points below 100°C) exhibit wide liquid ranges, non‐flammability, and thermal stability among other properties. These unique salts are best known as “green” alternatives to traditional volatile organic solvents, which are utilized in both academia and industry. Our current study compares the developmental toxicity potential of three representative ionic liquids, with various chain lengths: 1‐ethyl‐3‐methylimidazolium chloride ([C₂mim]Cl), 1‐butyl‐3‐methylimidazolium chloride ([C₄mim]Cl), and 1‐decyl‐3methylimidazolium chloride ([C₁₀mim]Cl). METHODS: From gestation days (GD) 6‐16, mated CD‐1 mice were orally dosed with one of the following: 1,000, 2,000, or 3,000 mg/kg/day [C₂mim]Cl; 113, 169, or 225 mg/kg/day [C₄mim]Cl; 50, 75, or 100 mg/kg/day [C₁₀mim]Cl; or the vehicle only. Dams were sacrificed on GD 17, and their litters were examined for adverse effects. RESULTS: Fetal weight was significantly decreased in the two highest dosage groups exposed to [C₄mim]Cl and [C₁₀mim]Cl in comparison with their controls, but the [C₂mim]Cl treated groups were not affected. An apparent teratogenic effect was associated with both [C₄mim]Cl and [C₁₀mim]Cl, as the offspring exhibited certain uncommon morphological defects. However, the incidences of malformations were low and no correlation between incidence and dosage could be made. No morphological defects were observed in any of the [C₂mim]Cl‐treated groups, despite maternal morbidity at the highest dosage level. CONCLUSIONS: This study indicates that [C₄mim]Cl and [C₁₀mim]Cl may have adverse effects on development at high maternal exposures and strongly supports the supposition that the toxicity of imidazolium‐based ILs is influenced by alkyl chain length. Birth Defects Res (Part B) 89:233–238, 2010. © 2010 Wiley‐Liss, Inc.     

Distributions of amino acids suggest that certain residue types more effectively determine protein secondary structure

Exponential growth in the number of available protein sequences is unmatched by the slower growth in the number of structures. As a result, the development of efficient and fast protein secondary structure prediction methods is essential for the broad comprehension of protein structures. Computational methods that can efficiently determine secondary structure can in turn facilitate protein tertiary structure prediction, since most methods rely initially on secondary structure predictions. Recently, we have developed a fast learning optimized prediction methodology (FLOPRED) for predicting protein secondary structure (Saraswathi et al. in JMM 18:4275, 2012). Data are generated by using knowledge-based potentials combined with structure information from the CATH database. A neural network-based extreme learning machine (ELM) and advanced particle swarm optimization (PSO) are used with this data to obtain better and faster convergence to more accurate secondary structure predicted results. A five-fold cross-validated testing accuracy of 83.8 % and a segment overlap (SOV) score of 78.3 % are obtained in this study. Secondary structure predictions and their accuracy are usually presented for three secondary structure elements: α-helix, β-strand and coil but rarely have the results been analyzed with respect to their constituent amino acids. In this paper, we use the results obtained with FLOPRED to provide detailed behaviors for different amino acid types in the secondary structure prediction. We investigate the influence of the composition, physico-chemical properties and position specific occurrence preferences of amino acids within secondary structure elements. In addition, we identify the correlation between these properties and prediction accuracy. The present detailed results suggest several important ways that secondary structure predictions can be improved in the future that might lead to improved protein design and engineering.     

Structural and gene expression analyses of uptake hydrogenases and other proteins involved in nitrogenase protection in Frankia

The actinorhizal bacterium Frankia expresses nitrogenase and can therefore convert molecular nitrogen into ammonia and the by-product hydrogen. However, nitrogenase is inhibited by oxygen. Consequently, Frankia and its actinorhizal hosts have developed various mechanisms for excluding oxygen from their nitrogen-containing compartments. These include the expression of oxygen-scavenging uptake hydrogenases, the formation of hopanoid-rich vesicles, enclosed by multi-layered hopanoid structures, the lignification of hyphal cell walls, and the production of haemoglobins in the symbiotic nodule. In this work, we analysed the expression and structure of the so-called uptake hydrogenase (Hup), which catalyses the in vivo dissociation of hydrogen to recycle the energy locked up in this ‘waste’ product. Two uptake hydrogenase syntons have been identified in Frankia: synton 1 is expressed under free-living conditions while synton 2 is expressed during symbiosis. We used qPCR to determine synton 1 hup gene expression in two Frankia strains under aerobic and anaerobic conditions. We also predicted the 3D structures of the Hup protein subunits based on multiple sequence alignments and remote homology modelling. Finally, we performed BLAST searches of genome and protein databases to identify genes that may contribute to the protection of nitrogenase against oxygen in the two Frankia strains. Our results show that in Frankia strain ACN14a, the expression patterns of the large (HupL1) and small (HupS1) uptake hydrogenase subunits depend on the abundance of oxygen in the external environment. Structural models of the membrane-bound hydrogenase subunits of ACN14a showed that both subunits resemble the structures of known [NiFe] hydrogenases (Volbeda et al. 1995), but contain fewer cysteine residues than the uptake hydrogenase of the Frankia DC12 and Eu1c strains. Moreover, we show that all of the investigated Frankia strains have two squalene hopane cyclase genes (shc1 and shc2). The only exceptions were CcI3 and the symbiont of Datisca glomerata, which possess shc1 but not shc2. Four truncated haemoglobin genes were identified in Frankia ACN14a and Eu1f, three in CcI3, two in EANpec1 and one in the Datisca glomerata symbiont (Dg).     

A Model of the Strain-induced B-Z Transition

Abstract

The B-to-Z transition in supercoiled circular DNA is modeled as a strain-induced nonlinear excitation process. Using a model, in which DNA is regarded as a chain of units with a bistable energy function along the twisting coordinate together with a harmonic inter-unit interaction, we show that a Z region and the accompanying two B-Z junctions of finite width appear naturally as a solution of nonlinear equations, when the strain exceeds a critical value. We examine the B-Z transition behaviour as a function of twist under various situations. We also analyse available experimental results on B-Z transition in supercoiled plasmid with G-C insertions by this mechanistic model in order to estimate the magnitude of model parameters. The energy barrier of the B-Z transition is estimated to be of the order of 1 kcal/mole per base pair. The analysis shows that if the length of the insertion is less than a certain value, the entire insertion converts to Z form at a transition point, but if the insertion is much longer, the B-Z transition exhibits a different behavior, in which part of the insertion flips to Z form and the Z region expands linearly upon changing linking number.     

Loop Folds in Proteins and Evolutionary Conservation of Folding Nuclei

Abstract

We show that loops of close contacts involving hydrophobic residues are important in protein folding. Contrary to Berezovsky and Trifonov (J. Biomol. Struct. Dyn. 20, 5–6, 2002) the loops important in protein folding usually are much larger in size than 23–31 residues, being instead comparable to the size of the protein for single domain proteins. Additionally what is important are not single loop contacts, but a highly interconnected network of such loop contacts, which provides extra stability to a protein fold and which leads to their conservation in evolution.