Using Variable Rate Models to Identify Genes Under Selection in Sequence Pairs: Their Validity and Limitations for EST Sequences

Using likelihood-based variable selection models, we determined if positive selection was acting on 523 EST sequence pairs from two lineages of sunflower and lettuce. Variable rate models are generally not used for comparisons of sequence pairs due to the limited information and the inaccuracy of estimates of specific substitution rates. However, previous studies have shown that the likelihood ratio test (LRT) is reliable for detecting positive selection, even with low numbers of sequences. These analyses identified 56 genes that show a signature of selection, of which 75% were not identified by simpler models that average selection across codons. Subsequent mapping studies in sunflower show four of five of the positively selected genes identified by these methods mapped to domestication QTLs. We discuss the validity and limitations of using variable rate models for comparisons of sequence pairs, as well as the limitations of using ESTs for identification of positively selected genes. [Reviewing Editor: Dr. Rasmus Nielsen]     

Effect of succinylation of antibodies on their conformation and interaction with the antigen

Using succinic anhydride, a succinylated derivative of anti-urease IgG having 49 ± 6% modification was prepared and its physicochemical and immunological properties were studied. IgG undergoes substantial changes in its native conformation on succinylation, which was mainly attributed to electrostatic destabilization of the native protein conformation. The modified IgG exhibited a decrease in its cross-reactivity with urease. This decrease is attributed to the conformational change in IgG upon succinylation and/or is due to the disruption of the lysine residues in the antigen-binding site of IgG upon succinylation, which may be involved in binding the antigen. IgG was able to bind to the specific antigen although its conformation was partially modified. Therefore, partial modification of the conformation of the antigen-binding site of IgG is permissible in order to bind to the antigen. Published in Russian in Biokhimiya, 2006, Vol. 71, No. 12, pp. 1642–1647.     

A DIFFERENTIAL GEOMETRIC APPROACH TO PROTEIN STRUTURE ALIGNMENT

A new method of comparing protein structures is deseribed, based on the differential-geometric representaion of protein conformation with the information of global dissimilarity and residue relatedness.The utility of this method for comparing closely and distantly related homologous proteins is demostrated by multiple alignment of globins, serine proteinases and aspartic proteinases.     

Sequence instability in the long terminal repeats of avian spleen necrosis virus and reticuloendotheliosis virus

Summary Sequence divergence between the 3 long terminal repeats (LTR) of avian reticuloendotheliosis virus (REV), deletion variant proviral clone 2-20-4, and spleen necrosis virus (SNV)—proviral clones 14-44, 60, and 70—was found to involve two classes of base substitutions: low-frequency interspersed and high-frequency clustered substitutions. Clones 2-20-4 and 14-44 have diverged 4.4% owing to low-frequency substitutions. In contrast, two high-frequency substitution segments have diverged by 30% and 29%, respectively. Clustered substitutions appear to be located either within or next to tandem repeats, suggesting their introduction concomitant with sequence deletions and duplications commonly associated with such repeats. A new 19-bp tandem repeat is found in clone 2-20-4. Its sequence could have evolved from the 26-bp repeats found in the SNV clones.     

Folding/Unfolding Kinetics of Apomyoglobin

The equilibrium and kinetic folding/unfolding of apomyoglobin (ApoMb) were studied at pH 6.2, 11 °C by recording tryptophan fluorescence. The equilibrium unfolding of ApoMb in the presence of urea was shown to involve accumulation of an intermediate state, which had a higher fluorescence intensity as compared with the native and unfolded states. The folding proceeded through two kinetic phases, a rapid transition from the unfolded to the intermediate state and a slow transition from the intermediate to the native state. The accumulation of the kinetic intermediate state was observed in a wide range of urea concentrations. The intermediate was detected even in the region corresponding to the unfolding limb of the chevron plot. Urea concentration dependence was obtained for the observed folding/unfolding rate. The shape of the dependence was compared with that of two-state proteins characterized by a direct transition from the unfolded to the native state.     

Determination of the Phosphorylated Metabolites of Gemcitabine and of Difluorodeoxyuridine by LCMSMS

Gemcitabine is an established chemotherapy agent in several solid tumors. Its mechanism of action has been theoretically established and this is supported with strong experimental evidence. However, certain aspects of the resistance mechanism for this agent remain elusive. We present a method of analysis using tandem liquid chromatography and mass spectrometry that provides a broader, yet more focused view of the action of gemcitabine and its primary metabolite, difluorodeoxyuridine in relation to the (deoxy) nucleoside and (deoxy) nucleotide pools in tumor cell lines.

Alcoholic cytosole extracts were incubated with alkaline phosphatase reducing the nucleotide pools to their respective nucleosides. Determination of the nucleoside content by a sensitive LCMSMS method before and after incubation enables the calculation of the total amount of phosphorylation of each (deoxy) nucleoside in the cell. Incubation with clinically relevant levels of gemcitabine (dFdC) or difluorodeoxyuridine (dFdU) for 24 hours enabled the determination of the changes in the (deoxy) nucleotide pools in relation to chemotherapeutic and toxicological effects. Confirmation of the presence of dFdC phosphorylation is presented as well as direct evidence of dFdU phosphorylation after both dFdC and dFdU treatment. Differences in the nucleotide pools are presented after dFdC and dFdU incubation, indicating that dFdU might have more chemotherapeutic properties than previously believed.     

tDNA^Trp识别位碱基的研究

识别位碱基Ｇ７３,为ｔＲＮＡ＾Ｔｒｐ的主要个性元件之一,ｔＤＮＡ＾Ｔｒｐ－ＮＣＣｒＡ的氨酰化活力测定及圆二色谱均表明在相同的ｐＨ条件下,识别位碱基的改变将影响到ｔＲＮＡ３′端的构象,而具有同一识别位碱基的ｔＤＮＡ＾Ｔｒｐ在不同ｐＨ条件下,亦显示出不同的构象及氨酰化活力,ｔＤＮＡ＾Ｔｒｐ的研究表明,是ｔＲＮＡ的构象而不是碱基本身决定了ｔＲＮＡ与其相应合成酶之间的精确识别。     

A multiscale investigation of mechanical properties of bio-inspired scaffolds

Abstract

Finding a structural design which allows the scaffold to have a high porosity and large pore size while retaining high strength is essential. Here, a bio-inspired scaffold is designed based on the observed geometrical pattern of the apatite atomic crystal structure, and mechanical properties are compared with other common scaffold geometries. The bio-inspired scaffold design is proven superior using a multiscale computational approach, which combines density functional theory and finite element analysis to predict the stress reaction and substitution effects on the scaffolds. This study provides insight into better scaffold design using bio-inspired structures and the effect of substitutions.     

Structural characterization of unfolded states of apomyoglobin using residual dipolar couplings

The conformational propensities of unfolded states of apomyoglobin have been investigated by measurement of residual dipolar couplings between (15)N and (1)H in backbone amide groups. Weak alignment of apomyoglobin in acid and urea-unfolded states was induced with both stretched and compressed polyacrylamide gels. In 8 M urea solution at pH 2.3, conditions under which apomyoglobin contains no detectable secondary or tertiary structure, significant residual dipolar couplings of uniform sign were observed for all residues. At pH 2.3 in the absence of urea, a change in the magnitude and/or sign of the residual dipolar couplings occurs in local regions of the polypeptide where there is a high propensity for helical secondary structure. These results are interpreted on the basis of the statistical properties of the unfolded polypeptide chain, viewed as a polymer of statistical segments. For a folded protein, the magnitude and sign of the residual dipolar couplings depend on the orientation of each bond vector relative to the alignment tensor of the entire molecule, which reorients as a single entity. For unfolded proteins, there is no global alignment tensor; instead, residual dipolar couplings are attributed to alignment of the statistical segments or of transient elements of secondary structure. For apomyoglobin in 8 M urea, the backbone is highly extended, with phi and psi dihedral angles favoring the beta or P(II) regions. Each statistical segment has a highly anisotropic shape, with the N-H bond vectors approximately perpendicular to the long axis, and becomes weakly aligned in the anisotropic environment of the strained acrylamide gels. Local regions of enhanced flexibility or chain compaction are characterized by a decrease in the magnitude of the residual dipolar couplings. The formation of a small population of helical structure in the acid-denatured state of apomyoglobin leads to a change in sign of the residual dipolar couplings in local regions of the polypeptide; the population of helix estimated from the residual dipolar couplings is in excellent agreement with that determined from chemical shifts. The alignment model described here for apomyoglobin can also explain the pattern of residual dipolar couplings reported previously for denatured states of staphylococcal nuclease and other proteins. In conjunction with other NMR experiments, residual dipolar couplings can provide valuable insights into the dynamic conformational propensities of unfolded and partly folded states of proteins and thereby help to chart the upper reaches of the folding landscape.     

Variation in synonymous substitution rates among mammalian genes and the correlation between synonymous and nonsynonymous divergences

Using mammalian gene sequences, the variances in the numbers of synonymous and nonsynonymous substitutions among genes were estimated together with the correlation coefficient between the two. The expected correlation coefficient can be obtained under the neutral theory using these estimated values of the variances. The expected coefficient is found to often be one-half to two-thirds of the observed value. Possible causes for the disagreement were discussed, such as correlated selective constraints on the two types of substitutions and excess doublet mutations. The variance of mutation rate and that of selective constraint were also estimated. The results show that the coefficient of variation of the former is 0.2–0.3, whereas that of the latter is 0.7–0.9.Correspondence to: T. Ohta     

TM-MOTIF: an alignment viewer to annotate predicted transmembrane helices and conserved motifs in aligned set of sequences

Multiple sequence alignments become biologically meaningful only if conserved and functionally important residues and secondary structural elements preserved can be identified at equivalent positions. This is particularly important for transmembrane proteins like G-protein coupled receptors (GPCRs) with seven transmembrane helices. TM-MOTIF is a software package and an effective alignment viewer to identify and display conserved motifs and amino acid substitutions (AAS) at each position of the aligned set of homologous sequences of GPCRs. The key feature of the package is to display the predicted membrane topology for seven transmembrane helices in seven colours (VIBGYOR colouring scheme) and to map the identified motifs on its respective helices /loop regions. It is an interactive package which provides options to the user to submit query or pre-aligned set of GPCR sequences to align with a reference sequence, like rhodopsin, whose structure has been solved experimentally. It also provides the possibility to identify the nearest homologue from the available inbuilt GPCR or Olfactory Receptor cluster dataset whose association is already known for its receptor type. AVAILABILITY: The database is available for free at mini@ncbs.res.in.     

Insights from the analysis of conserved motifs and permitted amino acid exchanges in the human, the fly and the worm GPCR clusters

G-protein coupled receptors (GPCRs) belong to biologically important and functionally diverse and largest super family of membrane proteins. GPCRs retain a characteristic membrane topology of seven alpha helices with three intracellular, three extracellular loops and flanking N' and C' terminal residues. Subtle differences do exist in the helix boundaries (TM-domain), loop lengths, sequence features such as conserved motifs, and substituting amino acid patterns and their physiochemical properties amongst these sequences (clusters) at intra-genomic and inter-genomic level (please re-phrase into 2 statements for clarity). In the current study, we employ prediction of helix boundaries and scores derived from amino acid substitution exchange matrices to identify the conserved amino acid residues (motifs) as consensus in aligned set of homologous GPCR sequences. Co-clustered GPCRs from human and other genomes, organized as 32 clusters, were employed to study the amino acid conservation patterns and species-specific or cluster-specific motifs. Critical analysis on sequence composition and properties provide clues to connect functional relevance within and across genome for vast practical applications such as design of mutations and understanding of disease-causing genetic abnormalities.     

Secondary structure of peptidese: 13. 15N and 13C n.m.r. spectroscopic investigation of the helix stability of poly(l-alanine) in acidic solutions

¹⁵N-enriched poly(l-alanines) of various molecular weights were prepared from $\text{l}\text{-alanine}\text{-N-}\text{carboxyanhydride}$ (l-Ala-NCA) and their helix/coil equilibrium in trifluoroacetic acid (TFA) investigated by means of 40.5 MHz ¹⁵N nuclear magneic resonance (n.m.r.), 22.3 MHz ¹³C n.m.r. and circular dichroism (c.d.) spectra. The ¹⁵N n.m.r. spectra exhibit at least three peaks, and the dependence of their intensities on molecular weight, molecular weight distribution and temperature, as well as dynamic nuclear Overhauser effect (NOE) measurements, indicate that the high-field peak represents the helix fraction. All three spectroscopic methods agree that a helix→coil transition takes place with decreasing concentration. Furthermore, poly(l-alanines) containing d-alanine or glycine in various mole ratios were synthesizsed by copolymerizations of N-carboxyanhydrides (NCAs). The ¹⁵N n.m.r. spetra demonstrate that one d-Ala unit per 100 l-Ala units suffices to affect significantly the helix/coil equilibrium in TFA. In other words, the helix content under equilibrium conditions is highly sensitive to racemization. Furthermore, ¹³ C n.m.r. cross-polarization/magic angle spinning (CP/MAS) spectra demonstrate that the presence of d-Ala units also affects the α-helix content in the solid state.     

Histidyl residues at the active site of the Na/succinate cotransporter in rabbit renal brush borders

Summary Mono-, dicarboxylic acid-, andd-glucose transport were measured in brush border vesicles from renal cortex after treatment with reagents known to modify terminal amino, lysyl, -amino, guanidino, serine/threonine, histidyl, tyrosyl, tryptophanyl and carboxylic residues. All three sodium-coupled cotransport systems proved to possess sulfhydryl (and maybe tryptophanyl sulfhydryl, disulfide, thioether and tyrosyl) residues but not at the substrate site or at the allosteric cavity for the Na coion. Histidyl groups seem to be located in the active site of the dicarboxylic transporter in that the simultaneous presence of Na and succinate protects the transporter against the histidyl specific reagent diethylpyrocarbonate. Lithium, which specifically competes for sodium sites in the dicarboxylic acid transporter, substantially blocked the protective effect of Na and succinate. Hydroxylamine specifically reversed the covalent binding of diethylpyrocarbonate to the succinate binding site. The pH dependence of the Na/succinate cotransport is consistent with an involvement of histidyl and sulfhydryl residues. We conclude that a histidyl residue is at, or is close to, the active site of the dicarboxylate transporter in renal brush border membranes.     

银桦次生木质部导管分子观察研究

对银桦 (Grevillearobusta)次生木质部导管分子进行了观察研究。在银桦次生木质部导管分子中存在着许多不同的样式 ,分别对其进行了描述 ,并从导管分子个体发育与系统发育的角度进行了讨论。银桦导管分子中还存在着特殊的样式 ,如具三个单穿孔的导管以及导管分子侧壁上存在明显的皱褶 ,对其也分别进行了讨论