An algorithm to find distant repeats in a pair of protein sequences

Distant repeats between a pair of protein sequences can be exploited to study the various aspects of proteins such as structure-function relationship, disorders due to protein malfunction, evolutionary analysis, etc. An in-depth analysis of the distant repeats would facilitate to establish a stable evolutionary relation of the repeats with respect to their three-dimensional structure. To this effect, an algorithm has been devised to identify the distant repeats in a pair of protein sequences by essentially using the scores of PAM (Percent Accepted Mutation) matrices. The proposed algorithm will be of much use to researchers involved in the comparative study of various organisms based on the amino-acid repeats in protein sequences. (C) 2010 Elsevier B.V. All rights reserved.

DNA polymorphism and local variation in base-pair orientation: a theoretical rationale

Basepair stacking calculations have been carried out to understand the conformational polymorphism of DNA and its sequence dependence. The recently developed self-consistent parameter set, which is specially suitable for describing irregular DNA structures, has been used to describe the geometry of a basepair doublet. While for basepairs without any propeller, the favourable stacking patterns do not appear to have very strong features, much more noticeable sequence dependent stacking patterns emerge once a propeller is applied to the basepairs. The absolute minima for most sequences occurs for a doublet geometry close to the B-DNA fibre models. Hence in the B-DNA region, no strong sequence dependent features are found, but the range of doublet geometries observed in the crystal structures generally lie within the low energy contours, obtained from stacking energy calculations. The doublet geometry corresponding to the A-DNA fibre model is not energetically favourable for the purine-pyrimidine sequences, which prefer small roll angle values when the slide has a large negative value as in A-DNA. However positive roll with large negative slide is allowed for GG, GA, AG and the pyrimidine-purine steps. This is consistent with the observed geometries of various steps in A-DNA crystals. Thus the general features of the basepair doublets predicted from these theoretical studies agree very well with the results from crystal structure analysis. However, since most sequences show an overall preference for B-type doublet geometry, the B --> A transition for random sequence DNA cannot be explained on the basis of basepair stacking interactions.

Diastereomeric Drug Glucuronides: Enzymatic Glucuronidation and Analysis by Capillary Electrophoresis and Liquid Chromatography-Mass Spectrometry

Foreign compounds, such as drugs are metabolised in the body in numerous reactions. Metabolic reactions are divided into phase I (functionalisation) and phase II (conjugation) reactions. Uridine diphosphoglucuronosyltransferase enzymes (UGTs) are important catalysts of phase II metabolic system. They catalyse the transfer of glucuronic acid to small lipophilic molecules and convert them to hydrophilic and polar glucuronides that are readily excreted from the body. Liver is the main site of drug metabolism. Many drugs are racemic mixtures of two enantiomers. Glucuronidation of a racemic compound yields a pair of diastereomeric glucuronides. Stereoisomers are interesting substrates in glucuronidation studies since some UGTs display stereoselectivity. Diastereomeric glucuronides of O-desmethyltramadol (M1) and entacapone were selected as model compounds in this work. The investigations of the thesis deal with enzymatic glucuronidation and the development of analytical methods for drug metabolites, particularly diastereomeric glucuronides. The glucuronides were analysed from complex biological matrices, such as urine or from in vitro incubation matrices. Various pretreatment techniques were needed to purify, concentrate and isolate the analytes of interest. Analyses were carried out by liquid chromatography (LC) with ultraviolet (UV) or mass spectrometric (MS) detection or with capillary electromigration techniques. Commercial glucuronide standards were not available for the studies. Enzyme-assisted synthesis with rat liver microsomes was therefore used to produce M1 glucuronides as reference compounds. The glucuronides were isolated by LC/UV and ultra performance liquid chromatography (UPLC)/MS, while tandem mass spectrometry (MS/MS) and nuclear magnetic resonance (NMR) spectroscopy were employed in structural characterisation. The glucuronides were identified as phenolic O-glucuronides of M1. To identify the active UGT enzymes in (±)-M1 glucuronidation recombinant human UGTs and human tissue microsomes were incubated with (±)-M1. The study revealed that several UGTs can catalyse (±)-M1 glucuronidation. Glucuronidation in human liver microsomes like in rat liver microsomes is stereoselective. The results of the studies showed that UGT2B7, most probably, is the main UGT responsible for (±)-M1 glucuronidation in human liver. Large variation in stereoselectivity of UGTs toward (±)-M1 enantiomers was observed. Formation of M1 glucuronides was monitored with a fast and selective UPLC/MS method. Capillary electromigration techniques are known for their high resolution power. A method that relied on capillary electrophoresis (CE) with UV detection was developed for the separation of tramadol and its free and glucuronidated metabolites. The suitability of the method to identify tramadol metabolites in an authentic urine samples was tested. Unaltered tramadol and four of its main metabolites were detected in the electropherogram. A micellar electrokinetic chromatography (MEKC) /UV method was developed for the separation of the glucuronides of entacapone in human urine. The validated method was tested in the analysis of urine samples of patients. The glucuronides of entacapone could be quantified after oral entacapone dosing.

Conformational analysis and design of cross-strand disulfides in antiparallel beta-sheets

Cross-strand disulfides bridge two cysteines in a registered pair of antiparallel beta-strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross-strand disulfides. Seventy-six cross-strand disulfides were found of which 75 and 1 occurred at non-hydrogen-bonded (NHB) and hydrogen-bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive chi(1) value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1 degrees C in T-m. All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (Delta Delta G(0) = -3.3 to -6.7 kcal/mol). The data demonstrate that introduction of cross-strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity. Proteins 2011; 79: 244-260. (C) 2010 Wiley-Liss, Inc.

Free vibration analysis Of stochastic strings

The free vibration of strings with randomly varying mass and stiffness is considered. The joint probability density functions of the eigenvalues and eigenfunctions are characterized in terms of the solution of a pair of stochastic non-linear initial value problems. Analytical solutions of these equations based on the method of stochastic averaging are obtained. The effects of the mean and autocorrelation of the mass process are included in the analysis. Numerical results for the marginal probability density functions of eigenvalues and eigenfunctions are obtained and are found to compare well with Monte Carlo simulation results. The random eigenvalues, when normalized with respect to their corresponding deterministic values, are observed to tend to become first order stochastically stationary with respect to the mode count.

Germline BRCA1 mutation analysis in Indian breast/ovarian cancer families

Most of the predisposition to hereditary breast and ovarian cancer has been attributed to inherited defects in two tumor suppressor genes BRCA1 and BRCA2. To explore the contribution of BRCA1 mutations to hereditary breast cancer among Indian women, we examined the coding sequence of the BRCA1 gene in 14 breast cancer patients with a positive family history of breast and/or ovarian cancer. Mutation analysis was carried out using conformation sensitive gel electrophoresis (CSGE) followed by sequencing. Three mutations (21%) in the BRCA1 gene were identified. Two of them are novel mutations of which one is a missense mutation in exon 7 near the RING finger domain, while the other is a one base pair deletion in exon 11 which results in protein truncation. The third mutation, 185delAG, has been previously described in Ashkenazi Jewish families. To our knowledge this is the first report of a study of germline BRCA1 mutation analysis in familial breast cancer in India. Our data from 14 different families suggests a lower prevalence but definite involvement of germline mutations in the BRCA1 gene among Indian women with breast cancer and a family history of breast cancer.

Topological analysis of general linear networks

A new method of network analysis, a generalization in several different senses of existing methods and applicable to all networks for which a branch-admittance (or impedance) matrix can be formed, is presented. The treatment of network determinants is very general and essentially four terminal rather than three terminal, and leads to simple expressions based on trees of a simple graph associated with the network and matrix, and involving products of low-order, usually(2 times 2)determinants of tree-branch admittances, in addition to tree-branch products as in existing methods. By comparison with existing methods, the total number of trees and of tree pairs is usually considerably reduced, and this fact, together with an easy method of tree-pair sign determination which is also presented, makes the new method simpler in general. The method can be very easily adapted, by the use of infinite parameters, to accommodate ideal transformers, operational amplifiers, and other forms of network constraint; in fact, is thought to be applicable to all linear networks.

An analysis of substitution, deletion and insertion mutations in cancer genes

Cancer-associated mutations in cancer genes constitute a diverse set of mutations associated with the disease. To gain insight into features of the set, substitution, deletion and insertion mutations were analysed at the nucleotide level, from the COSMIC database. The most frequent substitutions were c -> t, g -> a, g -> t, and the most frequent codon changes were to termination codons. Deletions more than insertions, FS (frameshift) indels more than I-F (in-frame) ones, and single-nucleotide indels, were frequent. FS indels cause loss of significant fractions of proteins. The 5'-cut in FS deletions, and 5'-ligation in FS insertions, often occur between pairs of identical bases. Interestingly, the cut-site and 3'-ligation in insertions, and 3'-cut and join-pair in deletions, were each found to be the same significantly often (p < 0.001). It is suggested that these features aid the incorporation of indel mutations. Tumor suppressors undergo larger numbers of mutations, especially disruptive ones, over the entire protein length, to inactivate two alleles. Proto-oncogenes undergo fewer, less-disruptive mutations, in selected protein regions, to activate a single allele. Finally, catalogues, in ranked order, of genes mutated in each cancer, and cancers in which each gene is mutated, were created. The study highlights the nucleotide level preferences and disruptive nature of cancer mutations.

NMR Analysis of Cross Strand Aromatic Interactions in an 8 Residue Hairpin and a 14 Residue Three Stranded beta-Sheet Peptide

Cross strand aromatic interactions between a facing pair of phenylalanine residues in antiparallel beta-sheet structures have been probed using two structurally defined model peptides. The octapeptide Boc-(LFVPPLFV)-P-D-P-L-OMe (peptide 1) favors the beta-hairpin conformation nucleated by the type II' beta-turn formed by the (D)Pro-(L)Pro segment, placing Phe2 and Phe7 side chains in proximity. Two centrally positioned (D)Pro-(L)Pro segments facilitate the three stranded beta-sheet formation in the 14 residue peptide Boc-LFV(D)P(L)PLFVA(D)P(L)PLFV-OMe (peptide 2) in which the Phe2/Phe7 orientations are similar to that in the octapeptide. The anticipated folded conformations of peptides 1 and 2 are established by the delineation of intramolecularly hydrogen bonded NH groups and by the observation of specific cross strand NOEs. The observation of ring current shifted aromatic protons is a diagnostic of close approach of the Phe2 and Phe7 side chains. Specific assignment of aromatic proton resonances using HSQC and HSQC-TOCSY methods allow an analysis of interproton NOEs between the spatially proximate aromatic rings. This approach facilitates specific assignments in systems containing multiple aromatic rings in spectra at natural abundance. Evidence is presented for a dynamic process which invokes a correlated conformational change about the C-alpha-C-beta(chi(1)) bond for the pair of interacting Phe residues. NMR results suggest that aromatic ring orientations observed in crystals are maintained in solution. Anomalous temperature dependence of ring current induced proton chemical shifts suggests that solvophobic effects may facilitate aromatic ring clustering in apolar solvents.

Design and analysis of a peer-to-peer based distributed clock synchronisation protocol for wireless sensor networks

Clock synchronisation is an important requirement for various applications in wireless sensor networks (WSNs). Most of the existing clock synchronisation protocols for WSNs use some hierarchical structure that introduces an extra overhead due to the dynamic nature of WSNs. Besides, it is difficult to integrate these clock synchronisation protocols with sleep scheduling scheme, which is a major technique to conserve energy. In this paper, we propose a fully distributed peer-to-peer based clock synchronisation protocol, named Distributed Clock Synchronisation Protocol (DCSP), using a novel technique of pullback for complete sensor networks. The pullback technique ensures that synchronisation phases of any pair of clocks always overlap. We have derived an exact expression for a bound on maximum synchronisation error in the DCSP protocol, and simulation study verifies that it is indeed less than the computed upper bound. Experimental study using a few TelosB motes also verifies that the pullback occurs as predicted.

Improved LSB stegananalysis based on analysis of adjacent pixel pairs

We propose a simple, reliable method based on probability of transitions and distribution of adjacent pixel pairs for steganalysis on digital images in spatial domain subjected to Least Significant Bit replacement steganography. Our method is sensitive to the statistics of underlying cover image and is a variant of Sample Pair Method. We use the new method to estimate length of hidden message reliably. The novelty of our method is that it detects from the statistics of the underlying image, which is invariant with embedding, whether the results it calculate are reliable or not. To our knowledge, no steganalytic method so far predicts from the properties of the stego image, whether its results are accurate or not.

Static analysis for checking data format compatibility of programs

Large software systems are developed by composing multiple programs. If the programs manip-ulate and exchange complex data, such as network packets or files, it is essential to establish that they follow compatible data formats. Most of the complexity of data formats is associated with the headers. In this paper, we address compatibility of programs operating over headers of network packets, files, images, etc. As format specifications are rarely available, we infer the format associated with headers by a program as a set of guarded layouts. In terms of these formats, we define and check compatibility of (a) producer-consumer programs and (b) different versions of producer (or consumer) programs. A compatible producer-consumer pair is free of type mismatches and logical incompatibilities such as the consumer rejecting valid outputs gen-erated by the producer. A backward compatible producer (resp. consumer) is guaranteed to be compatible with consumers (resp. producers) that were compatible with its older version. With our prototype tool, we identified 5 known bugs and 1 potential bug in (a) sender-receiver modules of Linux network drivers of 3 vendors and (b) different versions of a TIFF image library.

Game theoretic analysis of collusions in nonneutral networks

This paper studies the impact of exclusive contracts between a content provider (CP) and an internet service provider (ISP) in a nonneutral network. We consider a simple linear demand function for the CPs. We studywhen an exclusive contract is benefcial to the colluding pair and evaluate its impact on the noncolluding players at equilibrium. For the case of two CPs and one ISP we show that collusion may not always be benefcial. We derive an explicit condition in terms of the advertisement revenues of the CPs that tells when a collusion is proftable to the colluding entities.

Curvature based mobility analysis and form closure of smooth planar curves with multiple contacts

This paper presents a simple second-order, curvature based mobility analysis of planar curves in contact. The underlying theory deals with penetration and separation of curves with multiple contacts, based on relative configuration of osculating circles at points of contact for a second-order rotation about each point of the plane. Geometric and analytical treatment of mobility analysis is presented for generic as well as special contact geometries. For objects with a single contact, partitioning of the plane into four types of mobility regions has been shown. Using point based composition operations based on dual-number matrices, analysis has been extended to computationally handle multiple contacts scenario. A novel color coded directed line has been proposed to capture the contact scenario. Multiple contacts mobility is obtained through intersection of the mobility half-spaces. It is derived that mobility region comprises a pair of unbounded or a single bounded convex polygon. The theory has been used for analysis and synthesis of form closure configurations, revolute and prismatic kinematic pairs. (C) 2013 Elsevier Ltd. All rights reserved.

PROTEIN DISULFIDE ANALYSIS AND DESIGN

The significant contribution of naturally occurring disulfide bonds to protein stability has encouraged development of methods to engineer non-native disulfides in proteins. These have yielded mixed results. We summarize applications of the program MODIP for disulfide engineering. The program predicts sites in proteins where disulfides can be stably introduced. The program has also been used as an aid in conformational analysis of naturally occurring disulfides in a-helices, antiparallel and parallel beta-strands. Disulfides in a-helices occur only at N-termini, where the first cysteine residue is the N-cap residue of the helix. The disulfide occurs as a CXXC motif and can possess redox activity. In antiparallel beta-strands, disulfides occur exclusively at non-hydrogen bonded (NHB) registered pairs of antiparallel beta-sheets with only 1 known natural example occurring at a hydrogen bonded (HB) registered pair. Conformational analysis suggests that disulfides between HB residue pairs are under torsional strain. A similar analysis to characterize disulfides in parallel beta-strands was carried out. We observed that only 9 instances of cross-strand disulfides exist in a non-redundant dataset. Stereochemical analysis shows that while tbe chi(square) angles are similar to those of other disulfides, the chi(1) and chi(2) angles show more variation and that one of tbe strands is generally an edge strand.