Extracellular volatilomic alterations induced by hypoxia in breast cancer cells.

INTRODUCTION: The metabolic shift induced by hypoxia in cancer cells has not been explored at volatilomic level so far. The volatile organic metabolites (VOMs) constitute an important part of the metabolome and their investigation could provide us crucial aspects of hypoxia driven metabolic reconfiguration in cancer cells. OBJECTIVE: To identify the altered volatilomic response induced by hypoxia in metastatic/aggressive breast cancer (BC) cells. METHODS: BC cells were cultured under normoxic and hypoxic conditions and VOMs were extracted using HS-SPME approach and profiled by standard GC-MS system. Univariate and multivariate statistical approaches (p < 0.05, Log2 FC ≥ 0.58/≤ - 0.58, PC1 > 0.13/< - 0.13) were applied to select the VOMs differentially altered after hypoxic treatment. Metabolic pathway analysis was also carried out in order to identify altered metabolic pathways induced by the hypoxia in the selected BC cells. RESULTS: Overall, 20 VOMs were found to be significantly altered (p < 0.05, PC1 > 0.13/< - 0.13) upon hypoxic exposure to BC cells. Further, cell line specific volatilomic alterations were extracted by comparative metabolic analysis of aggressive (MDA-MB-231) vs. non-aggressive (MCF-7) cells incubated under hypoxia and normoxia. In this case, 15 and 12 VOMs each were found to be significantly altered in aggressive cells when exposed to hypoxic and normoxic condition respectively. Out of these, 9 VOMs were found to be uniquely associated with hypoxia, 6 were specific to normoxia and 6 were found common to both the conditions. Formic acid was identified as the most prominent molecule with higher abundance levels in aggressive as compared to non-aggressive cells in both conditions. Furthermore, metabolic pathway analyses revealed that fatty acid biosynthesis and nicotinate and nicotinamide metabolism were significantly altered in aggressive as compared to non-aggressive cells in normoxia and hypoxia respectively. CONCLUSIONS: Higher formate overflow was observed in aggressive cells compared to non-aggressive cells incubated under both the conditions, reinforcing its correlation with aggressive and invasive cancer type. Moreover, under hypoxia, aggressive cells preferred to be bioenergetically more efficient whereas, under normoxia, fatty acid biosynthesis was favoured when compared to non-aggressive cells.

Novel mutations of the arylsulphatase B (ARSB) gene in Indian patients with mucopolysaccharidosis type VI.

BACKGROUND & OBJECTIVES: Mucopolysaccharidosis type VI (MPS VI) is a rare, autosomal recessive lysosomal storage disorder caused by deficient enzymatic activity of N-acetyl galactosamine-4-sulphatase resulting from mutations in the arylsulphatase B (ARSB) gene. The ARSB gene is located on chromosome 5q11-q13 and is composed of eight exons. More than hundred ARSB mutations have been reported so far, but the mutation spectrum of MPS VI in India is still unknown. Hence, the aim of the present study was to identify the mutational spectrum in patients with MPS VI in India and to study the genotype-phenotype association and functional outcomes of these mutations. METHODS: Molecular characterization of the ARSB gene by Sanger sequencing was done for 15 patients (aged 15 months to 11 yr) who were enzymatically confirmed to have MPS VI. Age of onset, clinical progression and enzyme activity levels in each patient were studied to look for genotype-phenotype association. Haplotype analysis performed for unrelated patients with the recurring mutation W450C, was suggestive of a founder effect. Sequence and structural analyses of the ARSB protein using standard software were carried out to determine the impact of detected mutations on the function of the ARSB protein. RESULTS: A total of 12 mutations were identified, of which nine were novel mutations namely, p.D53N, p.L98R, p.Y103SfsX9, p.W353X, p.H393R, p.F166fsX18, p.I220fsX5, p.W450L, and p.W450C, and three were known mutations (p.D54N, p.A237D and p.S320R). The nine novel sequence variants were confirmed not to be polymorphic variants by performing sequencing in 50 unaffected individuals from the same ethnic population. INTERPRETATION & CONCLUSIONS: Nine novel mutations were identified in MPS VI cases from India in the present study. The study also provides some insights into the genotype-phenotype association in MPS VI.

Molecular analyses of novel ASAH1 mutations causing Farber lipogranulomatosis: analyses of exonic splicing enhancer inactivating mutation.

Farber lipogranulomatosis is a rare autosomal recessive lysosomal storage disorder caused by mutations in the ASAH1 gene. In the largest ever study, we identified and characterized ASAH1 mutations from 11 independent Farber disease (FD) families. A total of 13 different mutations were identified including 1 splice, 1 polypyrimidine tract (PPT) deletion and 11 missense mutations. Eleven mutations were exclusive to the Indian population. The IVS6+4A>G splice and IVS5-16delTTTTC PPT deletion mutations resulted in skipping of exon 6 precluding thereby the region responsible for cleavage of enzyme precursor. A missense mutation (p.V198A) resulted in skipping of exon 8 due to inactivation of an exonic splicing enhancer (ESE) element. This is the first report of mutations affecting PPT and ESE in the ASAH1 gene resulting in FD.

Distal arthrogryposis type 5D with a novel ECEL1 gene mutation.

Distal arthrogryposis syndromes (DAs) show wide clinical variability and overlapping clinical findings with the other DAs classified by Bamshad et al. [1996]. Most of the DAs are inherited as autosomal dominant disorders. DA type 5D is a subtype of DA type 5 inherited as autosomal recessive disorder, clinically characterized by congenital distal joint contractures, knee extension contractures, congenital hip dislocation, club foot, ptosis and other eye findings, furrowed tongue, and scoliosis. Here, we report on a family with clinical features of DA type 5D with novel mutations in the ECEL1 gene.

GALNS mutations in Indian patients with mucopolysaccharidosis IVA.

Mucopolysaccharidosis IV A (Morquio syndrome A, MPS IVA) is a lysosomal storage disease caused by the deficiency of N-acetylgalactosamine-6-sulfatase (GALNS). The mutation spectrum in this condition is yet to be determined in Indians. We aimed to analyze the mutations in the GALNS gene in Asian Indians with MPS IVA. All the exons and the adjacent intronic regions of the gene were amplified and sequenced in sixty-eight unrelated Indian families. We identified 136 mutant alleles comprising of 40 different mutations. We report twenty-two novel mutations that comprise of seventeen missense (p.Asn32Thr, p.Leu36Arg, p.Pro52Leu, p.Pro77Ser, p.Cys79Arg, p.His142Pro, p.Tyr191Asp, p.Asn204Thr, p.Gly188Ser, p.Phe216Ser, p.Trp230Cys, p.Ala291Ser, p.Gly317Arg, p.His329Pro, p.Arg386Ser, p.Glu450Gly, p.Cys501Ser), three splice-site variants (c.120+1G>C, c.1003-3C>G, c.1139+1G>A), one nonsense mutation (p.Gln414*) and one frameshift mutation (p.Pro420Leufs*440). Eighteen mutations have been reported earlier. Among these p.Ser287Leu (8.82%), p.Phe216Ser (7.35%), p.Asn32Thr (6.61%) and p.Ala291Ser (5.88%) were the most frequent mutations in Indian patients but were rare in the mutational profiles reported in other populations. These results indicate that the Indian patients may have a distinct mutation spectrum compared to those of other populations. Mutant alleles in exon 1, 7 and 8 accounted for 44.8% of the mutations, and sequencing of these exons initially may be a cost-effective approach in Asian Indian patients. This is the largest study on molecular analysis of patients with MPS IVA reported in the literature, and the first report from India.

Molecular genetic analysis of MSUD from India reveals mutations causing altered protein truncation affecting the C-termini of E1α and E1ß.

Maple Syrup Urine Disease is a rare metabolic disorder caused by reduced/absent activity of the branched chain α-Ketoacid dehydrogenase enzyme complex. Mutations in BCKDHA, BCKDHB, and DBT, that encode important subunits of the enzyme complex namely E1α, E1ß, and E2, are the primary cause for the disease. We have performed the first molecular genetic analysis of MSUD from India on nine patients exhibiting classical MSUD symptoms. BCKDHA and BCKDHB mutations were identified in four and five patients, respectively including seven novel mutations namely the BCKDHA c.1249delC, c.1312T>C, and c.1561T>A and the BCKDHB c.401T>A, c.548G>A, c.964A>G, and c.1065delT. The BCKDHB c.970C>T (p.R324X) mutation was shown to trigger nonsense mediated decay-based degradation of the transcript. Seven of the total 11 mutations resulted in perturbations in the E1α or E1ß C-termini either through altered termination or through an amino acid change; these are expected to result in disruption of E1 enzyme complex assembly. Our study has therefore revealed that BCKDHA and BCKDHB mutations might be primarily responsible for MSUD in the Indian population.

A study on mutational dynamics of simple sequence repeats in relation to mismatch repair system in prokaryotic genomes.

Mutational bias toward expansion or contraction of simple sequence repeats (SSRs) is referred to as directionality of SSR evolution. In this communication, we report the mutational bias exhibited by mononucleotide SSRs occurring in the non-coding regions of several prokaryotic genomes. Our investigations revealed that the strains or species lacking mismatch repair (MMR) system generally show higher number of polymorphic SSRs than those species/strains having MMR system. An exception to this observation was seen in the mycobacterial genomes that are MMR deficient where only a few SSR tracts were seen with mutations. This low incidence of SSR mutations even in the MMR-deficient background could be attributed to the high fidelity of the DNA polymerases as a consequence of high generation time of the mycobacteria. MMR system-deficient species generally did not show any bias toward mononucleotide SSR expansions or contractions indicating a neutral evolution of SSRs in these species. The MMR-proficient species in which the observed mutations correspond to secondary mutations showed bias toward contraction of polymononucleotide tracts, perhaps, indicating low efficiency of MMR system to repair SSR-induced slippage errors on template strands. This bias toward deletion in the mononucleotide SSR tracts might be a probable reason behind scarcity for long poly A|T and G|C tracts in prokaryotic systems which are mostly MMR proficient. In conclusion, our study clearly demonstrates mutational dynamics of SSRs in relation to the presence/absence of MMR system in the prokaryotic system.

A founder ectodysplasin A receptor (EDAR) mutation results in a high frequency of the autosomal recessive form of hypohidrotic ectodermal dysplasia in India.

BACKGROUND: Hypohidrotic/anhidrotic ectodermal dysplasia (HED) is a rare Mendelian disorder affecting ectodermal tissues. The disease is primarily caused by inactivation of any one of three genes, namely ectodysplasin A1 (EDA-A1), which encodes a ligand belonging to the tumour necrosis factor (TNF) superfamily; ectodysplasin A receptor (EDAR), encoding the EDA-A1 receptor and ectodysplasin A receptor-associated death domain (EDARADD), encoding an adaptor protein. X-linked recessive (EDA-A1), the predominant form of HED, as well as autosomal recessive and dominant (EDAR and EDARADD) inheritance patterns have been identified in affected families. OBJECTIVES: To determine the common genes causing HED in India. METHODS: We performed mutation analysis on 26 HED families from India (including 30 patients). In addition, we carried out sequence and structural analysis of missense/nonsense and insertion/deletion mutations. RESULTS: Among the 26 families analysed, disease-causing EDAR mutations were identified in 12 (46%) while EDA-A1 mutations were detected in 11 (42%). Four novel mutations in EDAR and five in EDA-A1 were identified. More importantly, a possible founder EDAR mutation, namely c.1144G>A, was identified in five independent families, thus accounting for about one-fifth of affected families in whom mutation was detected. A majority of EDA-A1 mutations localized to the TNF-like domain while the location of EDAR mutations was more widespread. CONCLUSIONS: This is the first report of a founder EDAR mutation and of a significantly high frequency of autosomal recessive HED.

Support Vector Machine-based classification of protein folds using the structural properties of amino acid residues and amino acid residue pairs.

MOTIVATION: Fold recognition is a key step in the protein structure discovery process, especially when traditional sequence comparison methods fail to yield convincing structural homologies. Although many methods have been developed for protein fold recognition, their accuracies remain low. This can be attributed to insufficient exploitation of fold discriminatory features. RESULTS: We have developed a new method for protein fold recognition using structural information of amino acid residues and amino acid residue pairs. Since protein fold recognition can be treated as a protein fold classification problem, we have developed a Support Vector Machine (SVM) based classifier approach that uses secondary structural state and solvent accessibility state frequencies of amino acids and amino acid pairs as feature vectors. Among the individual properties examined secondary structural state frequencies of amino acids gave an overall accuracy of 65.2% for fold discrimination, which is better than the accuracy by any method reported so far in the literature. Combination of secondary structural state frequencies with solvent accessibility state frequencies of amino acids and amino acid pairs further improved the fold discrimination accuracy to more than 70%, which is approximately 8% higher than the best available method. In this study we have also tested, for the first time, an all-together multi-class method known as Crammer and Singer method for protein fold classification. Our studies reveal that the three multi-class classification methods, namely one versus all, one versus one and Crammer and Singer method, yield similar predictions. AVAILABILITY: Dataset and stand-alone program are available upon request.

A non-invasive approach to explore the discriminatory potential of the urinary volatilome of invasive ductal carcinoma of the breast.

Worldwide, breast invasive ductal carcinoma (IDC) accounts for the majority of the reported cases of this form of cancer. IDC effective management, as for any form of cancer, would greatly benefit from early diagnosis. This, however, due to various socio-economic reasons, is very far for the reality in developing countries like India, where cancer diagnosis is often carried out at late stages when disease management is troublesome. With the present work, we aim to evaluate a simple analytical methodology to identify a set of volatile organic compounds (VOCs) in urine samples, as a biosignature for IDC. Using solid-phase microextraction followed by gas chromatography/mass spectrometry, a panel of 14 urinary VOCs was found to discriminate IDC (n = 65) from a healthy control (HC) group (n = 70) through multivariate statistical treatments. Furthermore, metabolic pathway analysis revealed various dysregulated pathways involved in IDC patients hinting that their detailed investigations could lead to novel mechanistic insights into the disease pathophysiology. In addition, we validated the expression pattern of five of these VOCs namely 2-ethyl-1-hexanol, isolongifolenone, furan, dodecanoic acid, 2-methoxy-phenol in another external cohort of 59 urinary samples (IDC = 32 and HC = 27) and found their expression pattern to be consistent with the primary sample set. To our knowledge, this is the first study exploring breast IDC volatome alterations in Indian patients.

Volatilomic insight of head and neck cancer via the effects observed on saliva metabolites.

Head and neck cancer (HNC) is a heterogeneous malignant disease with distinct global distribution. Metabolic adaptations of HNC are significantly gaining clinical interests nowadays. Here, we investigated effects of HNC on differential expression of volatile metabolites in human saliva. We applied headspace solid phase microextraction coupled with gas chromatography-mass spectrometry analysis of saliva samples collected from 59 human subjects (HNC - 32, Control - 27). We identified and quantified 48 volatile organic metabolites (VOMs) and observed profound effects of HNC on these metabolites. These effects were VOM specific and significantly differed in the biologically comparable healthy controls. HNC induced changes in salivary VOM composition were well attributed to in vivo metabolic effects. A panel of 15 VOMs with variable importance in projection (VIP) score >1, false discovery rate (FDR) corrected p-value < 0.05 and log2 fold change (log2 FC) value of ≥0.58/≤-0.58 were regarded as discriminatory metabolites of pathophysiological importance. Afterwards, receiver operator characteristic curve (ROC) projected certain VOMs viz., 1,4-dichlorobenzene, 1,2-decanediol, 2,5-bis1,1-dimethylethylphenol and E-3-decen-2-ol with profound metabolic effects of HNC and highest class segregation potential. Moreover, metabolic pathways analysis portrayed several dysregulated pathways in HNC, which enhanced our basic understanding on salivary VOM changes. Our observations could redefine several known/already investigated systemic phenomenons (e.g. biochemical pathways). These findings will inspire further research in this direction and may open unconventional avenues for non-invasive monitoring of HNC and its therapy in the future.

SilkSatDb: a microsatellite database of the silkworm, Bombyx mori.

The SilkSatDb (silkmoth microsatellite database) (http://www.cdfd.org.in/silksatdb) is a relational database of microsatellites extracted from the available expressed sequence tags and whole genome shotgun sequences of the silkmoth, Bombyx mori. The database has been rendered with a simple and robust web-based search facility, developed using PHP. The SilkSatDb also stores information on primers developed and validated in the laboratory. Users can retrieve information on the microsatellite and the protocols used, along with informative figures and polymorphism status of those microsatellites. In addition, the interface is coupled with Autoprimer, a primer-designing program, using which users can design primers for the loci of interest.

Escherichia coli RNA polymerase mutations located near the upstream edge of an RNA:DNA hybrid and the beginning of the RNA-exit channel are defective for transcription antitermination by the N protein from lambdoid phage H-19B.

Transcription antitermination is an important mechanism that can control regulation of gene expression. The N protein of lambdoid phages modifies the transcription elongation complex (EC) and helps it to overcome downstream terminators. In this modified EC, the C-terminal domain of N makes specific interactions with RNA polymerase (RNAP). The interacting surface of RNAP for N is unknown. Here, we report five mutations in the beta (G1045D) and beta' (P251S, P254L, R270C and G336S) subunits of RNAP that are specifically defective for antitermination by N protein of the lambdoid phage, H-19B. A mutation in the C-terminal domain of N, L108F, suppresses the defect of beta'-P254L. Purified mutant holoenzymes exhibit less processive antitermination. The amino acid substitutions in the mutant RNAPs cluster very close to the RNA:DNA hybrid at the beginning of the RNA-exit channel of the EC. We suggest that the action of H-19B N is exerted through the region defined by these amino acids. Wild-type N stabilizes the EC at terminator sites and in this modified EC a part of the terminator hairpin may form but appears to be unstable. We propose that the action of N close to the active center alters the RNAP-nucleic acid interactions around the RNA:DNA hybrid, which impairs proper folding of the terminator hairpin or stabilizes the weak RNA:DNA hybrid, or both.

Interaction and localization diversities of global and local hubs in human protein-protein interaction networks.

Hubs, the highly connected nodes in protein-protein interaction networks (PPINs), are associated with several characteristic properties and are known to perform vital roles in cells. We defined two classes of hubs, global (housekeeping) and local (tissue-specific) hubs. These two categories of hubs are distinct from each other with respect to their abundance, structure and function. However, how distinct are the spatial expression pattern and other characteristics of their interacting partners is still not known. Our investigations revealed that the partners of the local hubs compared with those of global hubs are conserved across the tissues in which they are expressed. Partners of local hubs show diverse subcellular localizations as compared with the partners of global hubs. We examined the nature of interacting domains in both categories of hubs and found that they are promiscuous in global hubs but not so in local hubs. Deletion of some of the local and global hubs has an impact on the characteristic path length of the network indicating that those hubs are inter-modular in nature. Our present study has, therefore, shed further light on the characteristic features of the local and global hubs in human PPIN. This knowledge of different topological aspects of hubs with regard to their types and subtypes is essential as it helps in better understanding of roles of hub proteins in various cellular processes under various conditions including those caused by host-pathogen interactions and therefore useful in prioritizing targets for drug design and repositioning.

Spectrum of mutations in Glutaryl-CoA dehydrogenase gene in glutaric aciduria type I--Study from South India.

BACKGROUND: Glutaric aciduria type I is an autosomal recessive organic acid disorder. The primary defect is the deficiency of Glutaryl-CoA dehydrogenase (EC number 1.3.99.7) enzyme that is involved in the catabolic pathways of the amino acids l-lysine, l-hydroxylysine, and l-tryptophan. It is a treatable neuro-metabolic disorder. Early diagnosis and treatment helps in preventing brain damage. METHODS: The Glutaryl-CoA dehydrogenase gene (GCDH) gene was sequenced to identify disease causing mutations by direct sequencing of all the exons in twelve patients who were biochemically confirmed with GA I. RESULTS: We identified eleven mutations of which nine are homozygous mutations, one heterozygous and two synonymous mutations. Among the eleven mutations, four mutations p.Q162R, p.P286S, p.W225X in two families and p.V410M are novel. A milder clinical presentation is observed in those families who are either heterozygous or with a benign synonymous SNP. Multiple sequence alignment (MSA) of GCDH with its homologues revealed that the observed novel mutations are not tolerated by protein structure and function. CONCLUSIONS: The present study indicates genetic heterogeneity in GCDH gene mutations among South Indian population. Genetic analysis is useful in prenatal diagnosis and prevention. Mutation analysis is a useful tool in the absence of non-availability of enzyme assay in GA I.

Stereochemical punctuation marks in protein structures: glycine and proline containing helix stop signals.

An analysis on the nature of alpha-helix stop signals has been carried out, using a dataset of 1057 helices identified from 250 high resolution (

Recurrent and novel GLB1 mutations in India.

GM1 gangliosidosis is a lysosomal storage disorder caused by mutations in the GLB1 gene, leading to the deficiency of the enzyme ß-d-galactosidase. In this study, we report molecular findings in 50 Asian Indian families with GM1 gangliosidosis. We sequenced all the exons and flanking intronic sequences of GLB1 gene. We identified 33 different mutations (20 novel and 13 previously reported). The novel mutations include 12 missense (p.M1?, p.E129Q, p.G134R, p.L236P, p.G262E, p.L297F, p.Y331C, p.G414V, p.K493N, p.L514P, p.P597L, p.T600I), four splicing (c.246-2A>G, c.397-2A>G, c.552+1G>T, c.956-2A>G), three indels (p.R22Qfs*8, p.L24Cfs*47, p.I489Qfs*4) and one nonsense mutation (p.Q452*). Most common mutations identified in this study were c.75+2InsT (14%) and p.L337P (10%). Known mutations accounted for 67% of allele frequency in our cohort of patients, suggesting that these mutations in GLB1 are recurrent across different populations. Twenty three mutations were localized in the TIM barrel domain, ß-domain 1 and ß-domain 2. In silico sequence and structure analysis of GLB1 reveal that all the novel mutations affect the function and structure of the protein. We hereby report on the largest series of patients with GM1 gangliosidosis and the first from India.

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins.

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.

Influence of solvent molecules on the stereochemical code of glycyl residues in proteins.

The Ramachandran steric map and energy diagrams of the glycyl residue are symmetric. A plot of (phi,psi) angles of glycyl residues in 250 nonhomologous and high-resolution protein structures is also largely symmetric. However, there is a clear aberration in the symmetry. Although there is a cluster of points corresponding to the right-handed alpha-helical region, the "equivalent" cluster is clearly shifted to in and around the (phi,psi) values of (90 degrees, 0 degrees ) instead of being centered at the left-handed alpha-helical region of (60 degrees, 40 degrees ). This lack of symmetry exists even in the (phi,psi) distribution of residues from non-alpha-helical regions in proteins. Here we provide an explanation for this observation. An analysis of glycyl conformations in small peptide structures and in "coil" proteins, which are largely devoid of helical and sheet regions, shows that glycyl residues prefer to adopt conformations around (+/-90 degrees, 0 degrees ) instead of right- and left-handed alpha-helical regions. By using theoretical calculations, such conformations are shown to have highest solvent accessibility in a system of two-linked peptide units with glycyl residue at the central C(alpha) atom. This finding is consistent with the observations from 250 nonhomologous protein structures where glycyl residues with conformations close to (+/-90 degrees, 0 degrees ) are seen to have high solvent accessibility. Analysis of a subset of nonhomologous structures with very high resolution (1.5 A or better) shows that water molecules are indeed present at distances suitable for hydrogen bond interaction with glycyl residues possessing conformations close to (+/-90 degrees, 0 degrees ). It is suggested that water molecules play a key role in determining and stabilizing these conformations of glycyl residues and explain the aberration in the symmetry of glycyl conformations in proteins.