An in silico exploration of the neutral network in protein sequence space

Designating amino-acid sequences that fold into a common main-chain structure as "neutral sequences" for the structure, regardless of their function or stability, we investigated the distribution of neutral sequences in protein sequence space. For four distinct target structures (alpha, beta,alpha/beta and alpha+beta types) with the same chain length of 108, we generated the respective neutral sequences by using the inverse folding technique with a knowledge-based potential function. We assumed that neutral sequences for a protein structure have Z scores higher than or equal to fixed thresholds, where thresholds are defined as the Z score for the corresponding native sequence (case 1) or much greater Z score (case 2). An exploring walk simulation suggested that the neutral sequences mapped into the sequence space were connected with each other through straight neutral paths and formed an inherent neutral network over the sequence space. Through another exploring walk simulation, we investigated contiguous regions between or among the neutral networks for the distinct protein structures and obtained the following results. The closest approach distance between the two neutral networks ranged from 5 to 29 on the Hamming distance scale, showing a linear increase against the threshold values. The sequences located at the "interchange" regions between the two neutral networks have intermediate sequence-profile-scores for both corresponding structures. Introducing a "ball" in the sequence space that contains at least one neutral sequence for each of the four structures, we found that the minimal radius of the ball that is centered at an arbitrary position ranged from 35 to 50, while the minimal radius of the ball that is centered at a certain special position ranged from 20 to 30, in the Hamming distance scale. The relatively small Hamming distances (5-30) may support an evolution mechanism by transferring from a network for a structure to another network for a more beneficial structure via the interchange regions.     

An inverse correlation between expression of a preprocathepsin B-related protein with cysteine-rich sequences and steroid 11beta -hydroxylase in adrenocortical cells

A cDNA encoding a secretory protein hitherto unknown was cloned from mouse adrenocortical cells by subtractive hybridization between the cells without and with expressing steroid 11beta-hydroxylase (Cyp11b-1), a marker for the functional differentiation of cells in the zonae fasciculata reticularis (zFR). The deduced protein consisting of 466 amino acids contained a secretory signal, epidermal growth factor-like repeats, and a proteolytically inactive cathepsin B-related sequence. The amino acid sequence was 89% identical with that of human tubulointerstitial nephritis antigen-related protein. Among the mouse organs examined, adrenal glands prominently expressed its mRNA. The mRNA and its encoded protein were detected in the outer adrenocortical zones that do not express Cyp11b-1, i.e. the zona glomerulosa and the undifferentiated cell zone, while being undetectable in zFR that express Cyp11b-1. The new protein was designated as adrenocortical zonation factor 1 (AZ-1). Clonal lines with different levels of AZ-1 expression were established from Y-1 adrenocortical cells that originally express Cyp11b-1 but little AZ-1. Analyses of the clonal lines revealed that Cyp11b-1 is detected in the clonal lines maintaining little AZ-1 expression and becomes undetectable in those expressing AZ-1. On the other hand, irrespective of the AZ-1 expression, all clones expressed cholesterol side-chain cleavage enzyme, which occurs throughout the cortical zones. These results demonstrated that adrenocortical cells expressing AZ-1 do not express Cyp11b-1, whereas those with little AZ-1 express this zFR marker in vitro and in vivo, implying a putative role of AZ-1 in determining the zonal differentiation of adrenocortical cells.     

Quantitative proteomic analysis of Myc oncoprotein function

This study applies a new quantitative proteomics technology to the analysis of the function of the Myc oncoprotein in mammalian cells. Employing isotope-coded affinity tag (ICAT) reagent labeling and tandem mass spectrometry, the global pattern of protein expression in rat myc-null cells was compared with that of myc-plus cells (myc-null cells in which myc has been introduced) to generate a differential protein expression catalog. Expression differences among many functionally related proteins were identified, including reduction of proteases, induction of protein synthesis pathways and upregulation of anabolic enzymes in myc-plus cells, which are predicted to lead to increased cell mass (cell growth). In addition, reduction in the levels of adhesion molecules, actin network proteins and Rho pathway proteins were observed in myc-plus cells, leading to reduced focal adhesions and actin stress fibers as well as altered morphology. These effects are dependent on the highly conserved Myc Box II region. Our results reveal a novel cytoskeletal function for Myc and indicate the feasibility of quantitative whole-proteome analysis in mammalian cells.     

An efficient ligation method in the making of an in vitro virus for in vitro protein evolution

The “in vitro virus” is a molecular construct to perform evolutionary protein engineering. The “virion(=viral particle)”(mRNA-peptide fusion), is made by bonding a nascent protein with its coding mRNA via puromycin in a test tube for in vitro translation. In this work, the puromycin-linker was attached to mRNA using the Y-ligation, which was a method of two single-strands ligation at the end of a double-stranded stem to make a stem-loop structure. This reaction gave a yield of about 95%. We compared the Y-ligation with two other ligation reactions and showed that the Y-ligation gave the best productivity. An efficient amplification of the in vitro virus with this “viral genome” was demonstrated. Published: October 28, 2002     

Negative correlations between resistance to three organophosphate insecticides and productivity within a natural population of Drosophila melanogaster (Diptera: Drosophilidae)

To investigate the relationship between resistance to organophosphate insecticides and fitness components, we first measured resistance to three organophosphates, malathion, prothiophos, and fenitrothion, and productivity, a measure of fitness components, for each of the isofemale lines from the same natural population of Drosophila melanogaster (Meigen). Pearson correlation coefficients indicated that positive correlations among resistance to the organophosphates and negative correlations between resistance to each of the organophosphates and the productivity existed within the natural population. We further investigated the genetic basis of the correlations among resistance to the organophosphates and the productivity, by using chromosome-substituted lines between a resistant and a susceptible inbred line established from the same natural population. Chromosomal analyses indicated that the third chromosome from the resistant line exhibited not only significant, positive effects on resistance to all of the organophosphates tested but also a significant negative effect on the productivity, suggesting positive genetic correlations between resistance to each organophosphate and negative genetic correlations between resistance to each organophosphate and the productivity. In addition, a significant negative effect on the productivity was also detected from the second chromosome, which did not exhibit significant major effects on resistance to the organophosphates. This suggests that fitness components of resistant lines could be also affected by factors independent of insecticide resistance. The dynamics of genetic variation in resistance to the organophosphates within the natural population of D. melanogaster are discussed from the standpoint of negative genetic correlations between resistance to the organophosphates and the productivity.     

Molecular diversity of skin mucus lectins in fish

Among lectins in the skin mucus of fish, primary structures of four different types of lectin have been determined. Congerin from the conger eel Conger myriaster and AJL-1 from the Japanese eel Anguilla japonica were identified as galectin, characterized by its specific binding to β-galactoside. Eel has additionally a unique lectin, AJL-2, which has a highly conserved sequence of C-type lectins but displays Ca²⁺-independent activity. This is rational because the lectin exerts its function on the cutaneous surface, which is exposed to a Ca²⁺ scarce environment when the eel is in fresh water. The third type lectin is pufflectin, a mannose specific lectin in the skin mucus of pufferfish Takifugu rubripes. This lectin showed no sequence similarity with any known animal lectins but, surprisingly, shares sequence homology with mannose-binding lectins of monocotyledonous plants. The fourth lectin was found in the ponyfish Leiognathus nuchalis and exhibits homology with rhamnose-binding lectins known in eggs of some fish species. These lectins, except ponyfish lectin, showed agglutination of certain bacteria. In addition, pufflectin was found to bind to a parasitic trematode, Heterobothrium okamotoi. Taken together, these results demonstrate that skin mucus lectins in fish have wide molecular diversity.     

The comparison of intrinsic rates of increase among chromosome-substituted lines resistant and susceptible to organophosphate insecticides in Drosophila melanogaster

To investigate the genetic basis of the seasonal fluctuations in resistance to three organophosphates, observed within a natural population of Drosophila melanogaster (Meigen), we compared the intrinsic rate of increase, generation time and net reproduction rate among chromosome substitution lines derived from a resistant and a susceptible line, obtained from this natural population. There was significant variation among substituted lines; lines possessing the third chromosome from the resistant line, which confers resistance to the three organophosphates, generally showed lower mean values of these fitness measures. Chromosomal analyses also indicated significant negative contributions of the third chromosome from the resistant line. However, significant positive contributions of the interactions among chromosomes from the resistant line to these fitness measures were also detected. We further conducted a local stability analysis, in which each chromosome-substituted line was assumed to be introduced at a low frequency into the initial susceptible population. It was demonstrated that the resistance factor(s) on the third chromosome tend to decrease in their frequency under both density-independent and juvenile density-regulated conditions. Based on these results, a possible explanation for the seasonal fluctuations in resistance to the three organophosphates observed in the natural population was proposed.     

A longitudinal study on hand and wrist skeletal maturation in chimpanzees ( Pan troglodytes), with emphasis on growth in linear dimensions

Skeletal maturation in the chimpanzee hand and wrist (the RUS system; radius, ulna, and short bones) was studied both longitudinally and cross-sectionally. Maturity states were evaluated in each of the 13 bones of the RUS system based on the TW2 method (Tanner and Whitehouse method), and the RUS score was calculated by the summation of scores for these bones. Individual variation was examined by means of residual curves and pseudo-velocity curves of RUS score and anterior trunk length (ATL). Norms of the age change pattern in RUS skeletal maturation and the growth of ATL were determined for each sex, and the relationships among ATL growth and skeletal and reproductive maturation were examined. We found a fairly good relationship between ATL growth and RUS skeletal maturation. Comparison of growth and development between humans and chimpanzees showed that growth characteristics are coupled with each other at puberty in male chimpanzees and in both sexes of humans. Although nutritional condition influenced ATL growth in infancy, it had no effect on the RUS maturational process. Social relationships appeared to influence both ATL growth and RUS maturation. Analyses on relationships between RUS skeletal maturation, ATL growth, and reproductive maturation, showed that RUS skeletal maturation is a good indicator of "physiological age".     

Extraordinarily low bone mineral density in an old female chimpanzee (<Emphasis Type="Italic">Pan troglodytes schweinfurthii</Emphasis>) from the Mahale Mountains National Park

We examined bone mineral density (BMD) of the femoral neck and lumbar vertebrae of four chimpanzee skeletons from Mahale Mountains National Park, Tanzania, and four captive ones, with a dual energy X-ray absorptiometer. The BMD of Wansombo, an old female chimpanzee from Mahale , was remarkably lower than the mean of the other six younger adult female chimpanzees and categorized as osteoporosis. Posture, locomotion, and trunk-sacral anatomy of chimpanzees may have prevented fractures in Wansombo, whose BMD was below human osteoporosis criteria. Electronic Publication