Bilateral similarity function: A novel and universal method for similarity analysis of biological sequences

Bilateral similarity function is designed for analyzing the similarities of biological sequences such as DNA, RNA secondary structure or protein in this paper. The defined function can perform comprehensive comparison between sequences remarkably well, both in terms of the Hamming distance of two compared sequences and the corresponding location difference. Compared with the existing methods for similarity analysis, the examination of similarities/dissimilarities illustrates that the proposed method with the computational complexity of O(N) is effective for these three kinds of biological sequences, and bears the universality for them.     

Molecular modeling of cardiac glycoside binding by the human sequence monoclonal antibody 1B3

The amino acid sequences of the heavy- and light-chain variable regions of the high-affinity human sequence antidigoxin monoclonal antibody 1B3 (mAb 1B3) were determined, and a structural model for the mAb's variable region was developed by homology modeling techniques. The structural model provided the basis for computationally docking digoxin and eight related cardiac glycosides into the putative binding site of mAb 1B3. Analysis of the consensus binding mode obtained for digoxin showed that the cardenolide moiety of digoxin is deeply embedded in a predominantly hydrophobic, narrow cavity, whereas the terminal, gamma-carbohydrate group is solvent-exposed. The docking results indicated that the primary driving forces for digoxin binding by mAb 1B3 are hydrophobic interactions with the digoxin steroid ring system and hydrogen bonds with the digitoxose groups. The binding model accounts for the experimentally observed variations in mAb 1B3 binding affinity for various structural analogs of digoxin used previously to develop a 3D structure-activity relationship model of drug binding (Farr CD, Tabet MR, Ball WJ Jr, Fishwild DM, Wang X, Nair AC, Welsh WJ. Three-dimensional quantitative structure-activity relationship analysis of ligand binding to human sequence antidigoxin monoclonal antibodies using comparative molecular field analysis. J Med Chem 2002;45:3257-3270). In particular, the hydrogen bond pattern is consistent with the unique sensitivity of mAb 1B3's binding affinity to the number of sugar residues present in a cardiac glycoside. The hydrophobic environment about the steroid moiety of digoxin is compatible with the mAb's reduced affinity for ligands that possess hydrophilic hydroxyl and acetyl group modifications in this region. The model also indicated that most of the amino acid residues in contact with the ligand reside in or about the three complementarity determining regions (CDRs) of the heavy chain and the third CDR of the light chain. A comparison of the 1B3 binding model with the crystal structures of two murine antidigoxin mAbs revealed similar binding patterns used by the three mAbs, such as a high frequency of occurrence of aromatic, hydrophobic residues in the CDRs and a dominant role of the heavy chain CDR3 in antigen binding.     

U7 snRNAs： A Computational Survey

U7 small nuclear RNA (snRNA) sequences have been described only for a handful of animal species in the past. Here we describe a computational search for func- tional U7 snRNA genes throughout vertebrates including the upstream sequence elements characteristic for snRNAs transcribed by polymerase Ⅱ. Based on the results of this search, we discuss the high variability of U7 snRNAs in both se- quence and structure, and report on an attempt to find U7 snRNA sequences in basal deuterostomes and non-drosophilids insect genomes based on a combination of sequence, structure, and promoter features. Due to the extremely short se- quence and the high variability in both sequence and structure, no unambiguous candidates were found. These results cast doubt on putative U7 homologs in even more distant organisms that are reported in the most recent release of the Rfam database.     

谷胱甘肽硫转移酶(GST)的固定化及酶学特性研究

对谷胱甘肽硫转移酶的固定化、游离酶和固定化酶的酶学特性进行了研究,通过试验,确定谷胱甘肽硫转移酶的最佳固定化条件为先用2％壳聚糖吸附酶,然后再加戊二醛交联,交联用戊二醛浓度为1．2％,交联时间6h;游离酶的最适温度为45—55℃,最适pH值为6．5-7．0：固定化酶的最适温度为45-50℃,最适pH值为7．0;游离酶和固定化酶的最适酶促反应时间为30min。     

Importance of anisotropy in detachment rates for force production and cargo transport by a team of motor proteins

Many cellular processes are driven by collective forces generated by a team consisting of multiple molecular motor proteins. One aspect that has received less attention is the detachment rate of molecular motors under mechanical force/load. While detachment rate of kinesin motors measured under backward force increases rapidly for forces beyond stall‐force; this scenario is just reversed for non‐yeast dynein motors where detachment rate from microtubule decreases, exhibiting a catch‐bond type behavior. It has been shown recently that yeast dynein responds anisotropically to applied load, i.e. detachment rates are different under forward and backward pulling. Here, we use computational modeling to show that these anisotropic detachment rates might help yeast dynein motors to improve their collective force generation in the absence of catch‐bond behavior. We further show that the travel distance of cargos would be longer if detachment rates are anisotropic. Our results suggest that anisotropic detachment rates could be an alternative strategy for motors to improve the transport properties and force production by the team.     

GGIP: Structure and sequence‐based GPCR–GPCR interaction pair predictor

G Protein‐Coupled Receptors (GPCRs) are important pharmaceutical targets. More than 30% of currently marketed pharmaceutical medicines target GPCRs. Numerous studies have reported that GPCRs function not only as monomers but also as homo‐ or hetero‐dimers or higher‐order molecular complexes. Many GPCRs exert a wide variety of molecular functions by forming specific combinations of GPCR subtypes. In addition, some GPCRs are reportedly associated with diseases. GPCR oligomerization is now recognized as an important event in various biological phenomena, and many researchers are investigating this subject. We have developed a support vector machine (SVM)‐based method to predict interacting pairs for GPCR oligomerization, by integrating the structure and sequence information of GPCRs. The performance of our method was evaluated by the Receiver Operating Characteristic (ROC) curve. The corresponding area under the curve was 0.938. As far as we know, this is the only prediction method for interacting pairs among GPCRs. Our method could accelerate the analyses of these interactions, and contribute to the elucidation of the global structures of the GPCR networks in membranes. Proteins 2016; 84:1224–1233. © 2016 Wiley Periodicals, Inc.     

A facile method for controlling the reaction equilibrium of sphingolipid ceramide N-deacylase for lyso-glycosphingolipid production

Lyso-glycosphingolipids (lyso-GSLs), the N-deacylated forms of glycosphingolipids (GSLs), are important synthetic intermediates for the preparation of GSL analogs. Although lyso-GSLs can be produced by hydrolyzing natural GSLs using sphingolipid ceramide N-deacylase (SCDase), the yield for this reaction is usually low because SCDase also catalyzes the reverse reaction, ultimately establishing an equilibrium between hydrolysis and synthesis. In the present study, we developed an efficient method for controlling the reaction equilibrium by introducing divalent metal cation and detergent in the enzymatic reaction system. In the presence of both Ca²⁺ and taurodeoxycholate hydrate, the generated fatty acids were precipitated by the formation of insoluble stearate salts and pushing the reaction equilibrium toward hydrolysis. The yield of GM1 hydrolysis can be achieved as high as 96%, with an improvement up to 45% compared with the nonoptimized condition. In preparative scale, 75 mg of lyso-GM1 was obtained from 100 mg of GM1 with a 90% yield, which is the highest reported yield to date. The method can also be used for the efficient hydrolysis of a variety of GSLs and sphingomyelin. Thus, this method should serve as a facile, easily scalable, and general tool for lyso-GSL production to facilitate further GSL research.     

A Refined Open State of the Glycine Receptor Obtained via Molecular Dynamics Simulations

1. Download : Download high-res image (389KB)
2. Download : Download full-size image

     

Persistent degenerative changes in thymic organ function revealed by an inducible model of organ regrowth

The thymus is the most rapidly aging tissue in the body, with progressive atrophy beginning as early as birth and not later than adolescence. Latent regenerative potential exists in the atrophic thymus, because certain stimuli can induce quantitative regrowth, but qualitative function of T lymphocytes produced by the regenerated organ has not been fully assessed. Using a genome-wide computational approach, we show that accelerated thymic aging is primarily a function of stromal cells, and that while overall cellularity of the thymus can be restored, many other aspects of thymic function cannot. Medullary islet complexity and tissue-restricted antigen expression decrease with age, representing potential mechanisms for age-related increases in autoimmune disease, but neither of these is restored by induced regrowth, suggesting that new T cells produced by the regrown thymus will probably include more autoreactive cells. Global analysis of stromal gene expression profiles implicates widespread changes in Wnt signaling as the most significant hallmark of degeneration, changes that once again persist even at peak regrowth. Consistent with the permanent nature of age-related molecular changes in stromal cells, induced thymic regrowth is not durable, with the regrown organ returning to an atrophic state within 2 weeks of reaching peak size. Our findings indicate that while quantitative regrowth of the thymus is achievable, the changes associated with aging persist, including potential negative implications for autoimmunity.