Intrinsically Semi-disordered State and Its Role in Induced Folding and Protein Aggregation

Intrinsically disordered proteins (IDPs) refer to those proteins without fixed three-dimensional structures under physiological conditions. Although experiments suggest that the conformations of IDPs can vary from random coils, semi-compact globules, to compact globules with different contents of secondary structures, computational efforts to separate IDPs into different states are not yet successful. Recently, we developed a neural-network-based disorder prediction technique SPINE-D that was ranked as one of the top performing techniques for disorder prediction in the biannual meeting of critical assessment of structure prediction techniques (CASP 9, 2010). Here, we further analyze the results from SPINE-D prediction by defining a semi-disordered state that has about 50 % predicted probability to be disordered or ordered. This semi-disordered state is partially collapsed with intermediate levels of predicted solvent accessibility and secondary structure content. The relative difference in compositions between semi-disordered and fully disordered regions is highly correlated with amyloid aggregation propensity (a correlation coefficient of 0.86 if excluding four charged residues and proline, 0.73 if not). In addition, we observed that some semi-disordered regions participate in induced folding, and others play key roles in protein aggregation. More specifically, a semi-disordered region is amyloidogenic in fully unstructured proteins (such as alpha-synuclein and Sup35) but prone to local unfolding that exposes the hydrophobic core to aggregation in structured globular proteins (such as SOD1 and lysozyme). A transition from full disorder to semi-disorder at about 30–40 Qs is observed in the poly-Q (poly-glutamine) tract of huntingtin. The accuracy of using semi-disorder to predict binding-induced folding and aggregation is compared with several methods trained for the purpose. These results indicate the usefulness of three-state classification (order, semi-disorder, and full-disorder) in distinguishing nonfolding from induced-folding and aggregation-resistant from aggregation-prone IDPs and in locating weakly stable, locally unfolding, and potentially aggregation regions in structured proteins. A comparison with five representative disorder-prediction methods showed that SPINE-D is the only method with a clear separation of semi-disorder from ordered and fully disordered states.     

Identification of male-specific amh duplication,sexually differentially expressed genes and microRNAs at early embryonic development of Nile tilapia (Oreochromis niloticus)

Background

The probable influence of genes and the environment on sex determination in Nile tilapia suggests that it should be regarded as a complex trait. Detection of sex determination genes in tilapia has both scientific and commercial importance. The main objective was to detect genes and microRNAs that were differentially expressed by gender in early embryonic development.

Results

Artificial fertilization of Oreochromis niloticus XX females with either sex-reversed ΔXX males or genetically-modified YY ‘supermales’ resulted in all-female and all-male embryos, respectively. RNA of pools of all-female and all-male embryos at 2, 5 and 9 dpf were used as template for a custom Agilent eArray hybridization and next generation sequencing. Fifty-nine genes differentially expressed between genders were identified by a false discovery rate of p < 0.05. The most overexpressed genes were amh and tspan8 in males, and cr/20β-hsd, gpa33, rtn4ipl and zp3 in females (p < 1 × 10⁻⁹). Validation of gene expression using qPCR in embryos and gonads indicated copy number variation in tspan8, gpa33, cr/20β-hsd and amh. Sequencing of amh identified a male-specific duplication of this gene, denoted amhy, differing from the sequence of amh by a 233 bp deletion on exonVII, hence lacking the capability to encode the protein motif that binds to the transforming growth factor beta receptor (TGF-β domain). amh and amhy segregated in the mapping family in full concordance with SD-linked marker on LG23 signifying the QTL for SD. We discovered 831 microRNAs in tilapia embryos of which nine had sexually dimorphic expression patterns by a false discovery rate of p < 0.05. An up-regulated microRNA in males, pma-mir-4585, was characterized with all six predicted target genes including cr/20β-hsd, down-regulated in males.

Conclusions

This study reports the first discovery of sexually differentially expressed genes and microRNAs at a very early stage of tilapia embryonic development, i.e. from 2 dpf. Genes with sexually differential expression patterns are enriched for copy number variation. A novel male-specific duplication of amh, denoted amhy, lacking the TGF-β domain was identified and mapped to the QTL region on LG23 for SD, thus indicating its potential role in SD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-774) contains supplementary material, which is available to authorized users.     

A global machine learning based scoring function for protein structure prediction

We present a knowledge‐based function to score protein decoys based on their similarity to native structure. A set of features is constructed to describe the structure and sequence of the entire protein chain. Furthermore, a qualitative relationship is established between the calculated features and the underlying electromagnetic interaction that dominates this scale. The features we use are associated with residue–residue distances, residue–solvent distances, pairwise knowledge‐based potentials and a four‐body potential. In addition, we introduce a new target to be predicted, the fitness score, which measures the similarity of a model to the native structure. This new approach enables us to obtain information both from decoys and from native structures. It is also devoid of previous problems associated with knowledge‐based potentials. These features were obtained for a large set of native and decoy structures and a back‐propagating neural network was trained to predict the fitness score. Overall this new scoring potential proved to be superior to the knowledge‐based scoring functions used as its inputs. In particular, in the latest CASP (CASP10) experiment our method was ranked third for all targets, and second for freely modeled hard targets among about 200 groups for top model prediction. Ours was the only method ranked in the top three for all targets and for hard targets. This shows that initial results from the novel approach are able to capture details that were missed by a broad spectrum of protein structure prediction approaches. Source codes and executable from this work are freely available at http://mathmed.org /#Software and http://mamiris.com/ . Proteins 2014; 82:752–759. © 2013 Wiley Periodicals, Inc.     

From network structure to network reorganization: implications for adult neurogenesis

Networks can be dynamical systems that undergo functional and structural reorganization. One example of such a process is adult hippocampal neurogenesis, in which new cells are continuously born and incorporate into the existing network of the dentate gyrus region of the hippocampus. Many of these introduced cells mature and become indistinguishable from established neurons, joining the existing network. Activity in the network environment is known to promote birth, survival and incorporation of new cells. However, after epileptogenic injury, changes to the connectivity structure around the neurogenic niche are known to correlate with aberrant neurogenesis. The possible role of network-level changes in the development of epilepsy is not well understood. In this paper, we use a computational model to investigate how the structural and functional outcomes of network reorganization, driven by addition of new cells during neurogenesis, depend on the original network structure. We find that there is a stable network topology that allows the network to incorporate new neurons in a manner that enhances activity of the persistently active region, but maintains global network properties. In networks having other connectivity structures, new cells can greatly alter the distribution of firing activity and destroy the initial activity patterns. We thus find that new cells are able to provide focused enhancement of network only for small-world networks with sufficient inhibition. Network-level deviations from this topology, such as those caused by epileptogenic injury, can set the network down a path that develops toward pathological dynamics and aberrant structural integration of new cells.     

Genome holography: deciphering function-form motifs from gene expression data

Background

DNA chips allow simultaneous measurements of genome-wide response of thousands of genes, i.e. system level monitoring of the gene-network activity. Advanced analysis methods have been developed to extract meaningful information from the vast amount of raw gene-expression data obtained from the microarray measurements. These methods usually aimed to distinguish between groups of subjects (e.g., cancer patients vs. healthy subjects) or identifying marker genes that help to distinguish between those groups. We assumed that motifs related to the internal structure of operons and gene-networks regulation are also embedded in microarray and can be deciphered by using proper analysis.

Methodology/Principal Findings

The analysis presented here is based on investigating the gene-gene correlations. We analyze a database of gene expression of Bacillus subtilis exposed to sub-lethal levels of 37 different antibiotics. Using unsupervised analysis (dendrogram) of the matrix of normalized gene-gene correlations, we identified the operons as they form distinct clusters of genes in the sorted correlation matrix. Applying dimension-reduction algorithm (Principal Component Analysis, PCA) to the matrices of normalized correlations reveals functional motifs. The genes are placed in a reduced 3-dimensional space of the three leading PCA eigen-vectors according to their corresponding eigen-values. We found that the organization of the genes in the reduced PCA space recovers motifs of the operon internal structure, such as the order of the genes along the genome, gene separation by non-coding segments, and translational start and end regions. In addition to the intra-operon structure, it is also possible to predict inter-operon relationships, operons sharing functional regulation factors, and more. In particular, we demonstrate the above in the context of the competence and sporulation pathways.

Conclusions/Significance

We demonstrated that by analyzing gene-gene correlation from gene-expression data it is possible to identify operons and to predict unknown internal structure of operons and gene-networks regulation.     

Allometric relationships in young seedlings of faba bean (Vicia faba L.) following removal of certain root types

Sink-source relationships and allometric ratios were studied in young seedlings of faba bean (Vicia faba L.) following pruning of some root types. The plants were grown in an aeroponic system allowing an easy access to each part of the root system, throughout the experiment, without disturbing the others. Root, leaf and stem growth as well as their mineral content were determined in one group of undisturbed plants (CTRL) and in four groups of plants treated as follows: TAP – the distal-free portion of the taproot was removed; HALF – half the laterals were removed; ALL – all lateral roots were removed, and TAP+HALF – both the distal part of the taproot and half of the laterals were removed. Removal of all the lateral roots (ALL) induced the development of a longer taproot but severely arrested shoot growth. Phosphorus, potassium and calcium contents were lower in the plants of the ALL treatment. However, the content of Mg was practically unaffected. The effect of the HALF treatment was hardly noticeable but the effects of TAP+HALF treatment were cumulative. The allometric relationships between the surface area of the roots and that of the leaves were restored within the experimental period, apparently due to reduction in shoot growth. Removal of the distal parts of the taproot did not cause an increase in shoot growth. This indicates that the strength of the sinks (mostly of lateral roots) rather than that of the source determines these relationships.     

Mutations in LPIN1 cause recurrent acute myoglobinuria in childhood

Recurrent episodes of life-threatening myoglobinuria in childhood are caused by inborn errors of glycogenolysis, mitochondrial fatty acid beta-oxidation, and oxidative phosphorylation. Nonetheless, approximately half of the patients do not suffer from a defect in any of these pathways. Using homozygosity mapping, we identified six deleterious mutations in the LPIN1 gene in patients who presented at 2-7 years of age with recurrent, massive rhabdomyolysis. The LPIN1 gene encodes the muscle-specific phosphatidic acid phosphatase, a key enzyme in triglyceride and membrane phospholipid biosynthesis. Of six individuals who developed statin-induced myopathy, one was a carrier for Glu769Gly, a pathogenic mutation in the LPIN1 gene. Analysis of phospholipid content disclosed accumulation of phosphatidic acid and lysophospholipids in muscle tissue of the more severe genotype. Mutations in the LPIN1 gene cause recurrent rhabdomyolysis in childhood, and a carrier state may predispose for statin-induced myopathy.     

Management of synchronized network activity by highly active neurons

Increasing evidence supports the idea that spontaneous brain activity may have an important functional role. Cultured neuronal networks provide a suitable model system to search for the mechanisms by which neuronal spontaneous activity is maintained and regulated. This activity is marked by synchronized bursting events (SBEs)--short time windows (hundreds of milliseconds) of rapid neuronal firing separated by long quiescent periods (seconds). However, there exists a special subset of rapidly firing neurons whose activity also persists between SBEs. It has been proposed that these highly active (HA) neurons play an important role in the management (i.e. establishment, maintenance and regulation) of the synchronized network activity. Here, we studied the dynamical properties and the functional role of HA neurons in homogeneous and engineered networks, during early network development, upon recovery from chemical inhibition and in response to electrical stimulations. We found that their sequences of inter-spike intervals (ISI) exhibit long time correlations and a unimodal distribution. During the network's development and under intense inhibition, the observed activity follows a transition period during which mostly HA neurons are active. Studying networks with engineered geometry, we found that HA neurons are precursors (the first to fire) of the spontaneous SBEs and are more responsive to electrical stimulations.     

Is xylem cavitation resistance a relevant criterion for screening drought resistance among Prunus species?

Fruit trees are likely to suffer from the effects of severe drought in the future; however, sound criteria for evaluating the species' ability to survive these extreme conditions are largely missing. Here, we evaluated the feasibility of using xylem cavitation resistance as a tool for screening Prunus species for drought resistance. Ten different Prunus species were selected to cover a large range of water requirements, from hydrophilic to xerophilic types. Shoot cavitation resistance was evaluated with the new Cavitron technique. At this inter-specific level, cavitation resistance was related to species drought resistance, with xerophilic species being less vulnerable to cavitation. The Cavitron technique enabled species characterization that required a short time and small amounts of plant material. This technique could be used to evaluate the drought resistance of a limited number of fruit tree genotypes. Genotype screening on a larger scale, however, would likely require another approach. Out of a number of anatomical traits tested, a significant correlation was found between cavitation resistance and inter-vessel wall thickness across species. This anatomical trait is, therefore, suggested as a possible alternative to direct cavitation estimates and could be included as a screening criterion in breeding programs for drought resistance of Prunus genotypes.     

The effect of phosphorus availability on the carbon economy of contrasting common bean (Phaseolus vulgaris L.) genotypes

A common response to low phosphorus availability is increased relative biomass allocation to roots. The resulting increase in root:shoot ratio presumably enhances phosphorus acquisition, but may also reduce growth rates by diverting carbon to the production of heterotrophic rather than photosynthetic tissues. To assess the importance of increased carbon allocation to roots for the adaptation of plants to low P availability, carbon budgets were constructed for four common bean genotypes with contrasting adaptation to low phosphorus availability in the field ("phosphorus efficiency"). Solid-phase-buffered silica sand provided low (1 microM), medium (10 microM), and high (30 microM) phosphorus availability. Compared to the high phosphorus treatment, plant growth was reduced by 20% by medium phosphorus availability and by more than 90% by low phosphorus availability. Low phosphorus plants utilized a significantly larger fraction of their daytime net carbon assimilation on root respiration (c. 40%) compared to medium and high phosphorus plants (c. 20%). No significant difference was found among genotypes in this respect. Genotypes also had similar rates of P absorption per unit root weight and plant growth per unit of P absorbed. However, P-efficient genotypes allocated a larger fraction of their biomass to root growth, especially under low P conditions. Efficient genotypes had lower rates of root respiration than inefficient genotypes, which enabled them to maintain greater root biomass allocation than inefficient genotypes without increasing overall root carbon costs.