Search

Iron-sulfur (Fe-S) clusters are one of the earliest known metal complexes in biological molecules. Suf system is one of the Fe-S biogenesis pathways. SufA belongs to the Suf pathway. It is an A-type carrier protein that transfers Fe-S clusters from the scaffold to target proteins. Structural studies were performed for the Suf pathway protein, SufA, in order to explore the conformational changes that probably aid in the transfer of Fe-S clusters to target proteins. Three-dimensional (3D) structure of Plasmodium falciparum (Pf) SufA homodimer was obtained by homology modeling using 3D structure of Escherichia coli (Ec) SufA as template. Molecular dynamics (MD) simulation of Pf SufA and Ec SufA homodimers followed by trajectory and pocket analyses were carried out. A co-ordinated displacement of the homodimeric chains in the interfacial region, resembling a swinging trapeze-like movement was observed. Potential involvement of this swinging trapeze-like movement of the residues belonging to the interfacial region has been proposed as a probable mechanism that assists in the transfer of Fe-S cluster from SufA to apo proteins. This was substantiated by protein-protein interaction studies in Pf SufA by performing molecular docking of 3D conformations of Pf SufA obtained from MD trajectory at every 1 ns interval with Pf ferredoxin.Communicated by Ramaswamy H. Sarma.

Structural flexibility of the peptide linker connecting two domains is essential for the functioning of multi-domain complex. Nitric oxide synthase (NOS) isoforms contain the oxygenase and the reductase domains connected by calmodulin binding linker (CBL) region. Additionally, the endothelial NOS (eNOS) isoform contain an auto-inhibitory loop (AI) in the FMN reductase sub-domain which represses the inter-domain electron transfer process. Binding of Ca 2+ -Calmodulin complex on the CBL region relieves the AI loop repression and facilitates electron transfer from FMN in the reductase domain to the heme in the oxygenase domain. Few experimental studies have reported that in vitro mutation of Serine-615 (S615D) and Serine-633 (S633D) in the FMN reductase sub-domain to aspartic acid increased NO production and increased Ca 2+ sensitivity. To understand the role of AI loop in eNOS repression and activation in serine mutants (S615D and S633D), we modelled the FMN reductase sub-domain of human eNOS protein with and without the CBL region. Molecular dynamics simulations performed indicated that the mutant protein AI loop structure was stabilized by salt bridge formed between D615 and R602. It was also found that mutation increased the flexibility of C-terminal residues of eNOS CBL region. The hinge-like movement of the AI loop allowed rotation of the FMN sub-domain clockwise which may favour electron-transfer in the mutant protein. This study provides insight on mutation (S615D and S633D) induced changes in AI loop and increased flexibility of CBL region which may lead to the protein activation and may also facilitate Calmodulin binding at physiological Ca 2+ concentration. Graphical Abstract Mutation of amino acid residues contribute to structural changes at molecular level leading to alteration in protein dynamics and its function. Serine-615 and Serine-633 in the auto-inhibitory loop of human eNOS reductase model was mutated to aspartic acid in silico and molecular dynamics simulations of the protein showed that steric hindrance due to mutation altered the auto-inhibitory loop rearrangement and the FMN sub-domain movement favouring electron transfer.

The tripyrrolic antibiotic prodigiosin causes diverse reactions on its targets like energy spilling, membrane leakage, loss of motility and phototoxicity. It has bacteriostatic, bactericidal, anti-fungal, anti-cancer and immunosuppressive properties. Most of the functions suggest the role of prodigiosin in membrane disruption but the exact mechanism remains unknown. A molecular dynamics study was performed to understand the interactions of prodigiosin with the membrane. It was seen that prodigiosin from the solvent enters the membrane immediately either individually or as small clusters. Prodigiosin clusters with more than eight molecules do not appear to enter the membrane. Upon entry, the molecules orient themselves along the membrane-water interface with the pyrrole rings interacting with lipid head groups and with water. This orientation is stabilised by hydrogen bonding and hydrophobic interactions. The presence of prodigiosin molecules in the membrane changes the local lipid architecture and reduces the solvent accessibility of the membrane. The membrane fluidity, thickness or area per lipid head are largely unaffected. This suggests that prodigiosin could cause most damage in the vicinity of a membrane protein and thus could also explain the reason for varied effects on the targets., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Mycolic acids are long chain alpha-alkyl beta-hydroxy fatty acids that are major constituents of the cell wall of Mycobacterium tuberculosis. M. tuberculosis produces three main types of mycolic acids, alpha mycolic acids and keto and methoxy mycolic acids. Cycloproponated mycolic acids make the cell wall less permeable, contribute to antibiotic resistance and host immunomodulation and protect from injury. Cyclopropanation is catalyzed by enzymes of the Cyclopropane Mycolic Acid Synthase (CMAS) family. In the current study, we addressed two CMAS enzymes, proximal alpha cyclopropane mycolic acid synthase (PcaA/CmaA3) and keto cyclopropane Mycolic acid synthase (CmaA2). All-atom Molecular Dynamics (MD) simulations were performed for these enzymes for a timeframe of 100ns each (in triplicate), using GROMACS. Based on the PDB structures of apo and holo states of related CMAS enzymes, we generated a framework which helped us correlate active or inactive states of the enzymes to different conformations sampled by the enzymes during MD simulations. Dynamics suggested that the free or unbound enzymes have intrinsic memory and sample different states of catalysis even in the absence of the substrate/cofactor. Additionally, we find that F200, P201 and W204 may have functional significance. MD simulation of CmaA2 was performed with the objective of gaining insights into the putative role of a loop insert. Analysis showed that acidic residues of this loop possibly play an important role during the active state by forming salt bridges. The insights gained in this study can potentially be utilized for design of effective inhibitors against CMAS enzymes., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Sulfolobus solfataricus β-glycosidase (SS-βGly) belongs to Glycosyl Hydrolase family1 (GH1) with broad substrate specificity. SS-βGly catalyzes both hydrolysis and transglycosylation reactions. SS-βGly is commonly used to synthesize variety of galacto-oligosaccharides. A comparison of SS-βGly with bacterial and eukaryotic homologs, using DALI search, revealed unique inserts. Free enzyme molecular dynamics (MD) simulation was performed under two different pH conditions (pH 6.5 and 2.5) at a constant temperature of 65 °C using GROMACS. A probable active-site loop (residues 331-364) in SS-βGly was identified. Dynamics of substrate binding cavity revealed that it was buried and inaccessible during most timeframes at pH 6.5 whereas open and accessible at pH 2.5. New cavities identified during both simulations may act as probable water channel or product egress path. Analyses of docked complexes of 3D structures obtained at every 1ns interval with compounds, involved in hydrolysis and tranglycosylation reactions, demonstrated that conformational states sampled by SS-βGly during free enzyme dynamics mimic the stages in enzyme catalysis thereby providing a mechanistic perspective. Current study revealed that conformational changes were conducive for hydrolysis at pH 6.5 and multiple cycles of transglycosylation at pH 2.5. Probable role of salt-bridge interactions in determining the type of reaction mechanism was also explored., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Thalidomide is an infamous teratogen and it is continuously being explored for its anticancer properties. Fibroblast growth factor receptors (FGFRs) are implicated in embryo development and cancer pathophysiology. With striking similarities observed between FGFR implicated conditions and thalidomide embryopathy, we hypothesized thalidomide targets FGFRs. We utilized three different cell lines and chicken embryo model to investigate the effects of thalidomide and analogs on FGFR expression. We performed molecular docking, KINOMEscan analysis, and kinase activity assays to study the drug-protein interactions. The expression of FGFR1 and FGFR2 was differentially regulated by all the three drugs in cells as well as in developing organs. Transcriptome analysis of thalidomide-treated chick embryo strongly suggests the modulation of FGFR signaling and key transcription factors. Corroboration with previous studies suggests that thalidomide might affect FGFR expression through the transcription factor, E2F1. At the protein level, molecular docking predicted all three analogs to interact with lysine residue at 517th and 508th positions of FGFR2 and FGFR3, respectively. This lysine coordinates the ATP binding site of FGFR, thus hinting at the possible perturbation of FGFR activity by thalidomide. Kinome analysis revealed that kinase activities of FGFR2 and FGFR3 (G697C) reduced by 31% and 65%, respectively, in the presence of 10 μM thalidomide. Further, we checked and confirmed that the analogs inhibited the FGFR2 kinase activity in a dose-dependent manner. This study suggests that FGFRs could be potential targets of thalidomide and the two analogs, and also endorses the link between the teratogenicity and antitumor activities of the drugs.

Profilin is one of the actin-binding proteins that regulate dynamics of actin polymerization. It plays a key role in cell motility and invasion. It also interacts with several other proteins notably through its poly-L-proline (PLP) binding site. Profilin in apicomplexa is characterized by a unique mini-domain consisting of a large β-hairpin extension and an acidic loop which is relatively longer in Plasmodium species. Profilin is essential for the invasive blood stages of Plasmodium falciparum. In the current study, unbound profilins from Plasmodium falciparum (Pf), Toxoplasma gondii (Tg), and Homo sapiens (Hs) were subjected to molecular dynamics (MD) simulations for a timeframe of 100 ns each to understand the conformational dynamics of these proteins. It was found that the β-hairpin of profilins from Pf and Tg shows a hinge-like movement. This movement in Pf profilin may possibly be driven by the loss of a salt-bridge within profilin. The impact of this conformational change on actin binding was assessed by docking three dimensional (3D) structures of profilin from Pf and Tg with their corresponding actins using ClusPro2.0. The stability of docked Pf profilin-actin complex was assessed through a 50 ns MD simulation. As Hs profilin I does not have the apicomplexa specific mini-domain, MD simulation was performed for this protein and its dynamics was compared to that of profilins from Pf and Tg. Using an immunoinformatics approach, potential epitope regions were predicted for Pf profilin. This has a potential application in the design of vaccines as they mapped to its unique mini-domain., (© 2017 Wiley Periodicals, Inc.)

Hormone sensitive lipases (HSLs) play an important role in the survival of M. tuberculosis during dormancy. They help in the utilization of fatty acids from stored lipids. The objective of the current study was to identify all HSLs from the proteome of M. tuberculosis H37Rv. We have developed a novel HSL identification pipeline, based on amino acid sequence homology, presence of conserved motifs and other sequence features deciphered from known HSL dataset. Through this pipeline, we identified 10 proteins as putative HSLs in M. tuberculosis. We have annotated a lipase LipT, as putative p-nitrobenzyl esterase and also identified a new motif "PGG" which is a possible characteristic motif of a subfamily of HSLs., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Prevotella is part of the oral bacterial community implicated in periodontitis. Pan genome analyses of eight oral Prevotella species, P. dentalis, P. enoeca, P. fusca, P. melaninogenica, P. denticola, P. intermedia 17, P. intermedia 17-2 and P. sp. oral taxon 299 are presented in this study. Analysis of the Prevotella pan genome revealed features such as secretion systems, resistance to oxidative stress and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems that enable the bacteria to adapt to the oral environment. We identified the presence of type VI secretion system (T6SS) in P. fusca and P. intermedia strains. For some VgrG and Hcp proteins which were not part of the core T6SS loci, we used gene neighborhood analysis and identified putative effector proteins and putative polyimmunity loci in P. fusca and polymorphic toxin systems in P. intermedia strains. Earlier studies have identified the presence of Por secretion system (PorSS) in P. gingivalis, P. melaninogenica and P. intermedia. We noted the presence of their homologs in six other oral Prevotella studied here. We suggest that in Prevotella, PorSS is used to secrete cysteine proteases such as interpain and C-terminal domain containing proteins with a "Por_secre_tail" domain. We identified subtype I-B CRISPR-Cas system in P. enoeca. Putative CRISPR-Cas system subtypes for 37 oral Prevotella and 30 non-oral Prevotella species were also predicted. Further, we performed a BLASTp search of the Prevotella proteins which are also conserved in the red-complex pathogens, against the human proteome to identify potential broad-spectrum drug targets. In summary, the use of a pan genome approach enabled identification of secretion systems and defense mechanisms in Prevotella that confer adaptation to the oral cavity.

Antigen85 (Ag85) proteins of Mycobacterium tuberculosis are mycolyl transferases that aid in cell wall biosynthesis. MPT51 (Ag85D) is closely related to Ag85 proteins. We have performed a comparative molecular dynamics (MD) simulation study of Ag85 proteins (Ag85A, Ag85B, and Ag85C) and MPT51. We observe that helix α5, β7-α9 loop, and N-terminal region of helix α9 of Ag85 proteins are mobile, suggestive of lid like movement over the active site. Further, in Ag85B, we observe the proposed scooting mode of the hydrophobic gating residue Phe232. Our simulations also show a similar scooting mode for Phe232 of Ag85A and Trp158 of Ag85C. We also found aromatic residue clusters at the ends of the hydrophobic channel of Ag85 proteins, which may have functional significance. Although MPT51 lacks the tunnel, it has the aromatic clusters. The aromatic cluster region has the ability to bind trehalose. From an immunoinformatics study, a promiscuous linear epitope was identified in MPT51 which could be useful in subunit vaccine studies. Recent studies have shown that a mycobacterial protein HupB, interacts with Ag85 proteins and has a regulatory role in cell wall biogenesis, with implications in growth rate and latency. We performed molecular docking studies of HupB protein with Ag85 proteins and predicted potential sites of interaction in Ag85 proteins. The insights gained through the current study can potentially pave way for newer therapeutic interventions. Graphical Abstract Dynamics of antigen85 proteins and MPT51 from Mycobacterium tuberculosis.

Mycobacterium tuberculosis (Mtb) is the causative organism of tuberculosis. Extensively drug resistant strains and latency have posed formidable challenges in the treatment of tuberculosis. The current study addresses an alpha/beta hydrolase fold bearing enzyme, epoxide hydrolase Rv1938 from Mtb. Epoxide hydrolases are involved in detoxification processes, catabolism and regulation of signaling molecules. Using GROMACS, a 100ns Molecular Dynamics (MD) simulation was performed for Rv1938. Cavities were identified within the protein at various time frames of the simulation and their volumes were computed. During MD simulation, in addition to the substrate binding cavity, opening of two new cavities located behind the active site was observed. These cavities may be similar to the backdoor proposed for acetylcholinesterase. Structural superimposition of epoxide hydrolase from Mtb with the epoxide hydrolase of Agrobacterium radiobacter1 AD1 (Ephy) indicates that cavity1 in Mtb lies at an identical position to that of the water tunnel in Ephy. Further, docking of the substrate and an inhibitor with protein structures obtained from MD simulation at various time frames was also performed. The potential role of these cavities is discussed., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Rv0045c is an esterase involved in lipid metabolism of Mycobacterium tuberculosis. It belongs to the α/β hydrolase family. In the current study, we performed sequence- and structure-based analysis of Rv0045c followed by molecular dynamics (MD) simulation for 100 ns to investigate conformational changes in the enzyme. Sequence analysis revealed that this enzyme is possibly a hormone-sensitive lipase. Further, through structural analysis, a putative catalytic tetrad containing "Ser-Asp-Ser-His" and residues involved in the formation of an oxyanion hole were identified. MD simulation of Rv0045c revealed a conformational transition from an open to a closed state. The active site pocket was found to be gated by four loops. The potential role of the cap domain and the mobile histidine is discussed. From the simulation, we see that the conformational changes mimic the different stages in the reaction mechanism of Rv0045c. These results support the hypothesis that free enzyme simulation encompasses all the conformations necessary for the different stages of catalysis. Our findings add to the growing knowledge of an important family of esterases in Mycobacterium tuberculosis.

Adhesion of uropathogenic E. coli (UPEC) to uroepithelial cell receptors is facilitated through the lectin domain of FimH adhesin. In the current study, Molecular Dynamics (MD) simulations were performed for the lectin domain of FimH from UPEC J96. The high affinity state lectin domain was found to be stable and rigid during the simulations. Further, based on conserved subsequences around one of the disulfide forming cysteines, two sequence motifs were designed. An immunoinformatics approach was utilized to identify linear and discontinuous epitopes for the lectin domain of FimH. We propose that the accessibility of predicted epitopes should also be assessed in a dynamic aqueous environment to evaluate the potential of vaccine candidates. Since MD simulation data enables assessing the accessibility in a dynamic environment, we evaluated the accessibility of the top ranked discontinuous and linear epitopes using structures obtained at every nanosecond (ns) in the 1-20 ns MD simulation timeframe. Knowledge gained in this study has a potential utility in the design of vaccine candidates for Urinary Tract Infection (UTI)., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Chemotherapeutic resistance due to over expression of Inhibitor of Apoptosis Proteins (IAPs) XIAP, survivin and livin has been observed in various cancers. In the current study, Molecular Dynamics (MD) simulations were carried out for all three IAPs and a common ligand binding scaffold was identified. Further, a novel sequence based motif specific to these IAPs was designed. SMAC is an endogenous inhibitor of IAPs. Screening of ChemBank for compounds similar to lead SMAC-non-peptidomimetics yielded a cemadotin related compound NCIMech_000654. Cemadotin is a derivative of natural anti-tumor peptide dolastatin-15; hence these compounds were docked against all three IAPs. Based on our analysis, we propose that NCIMech_000654/dolastatin-15/cemadotin derivatives may be investigated for their potential in inhibiting XIAP, survivin and livin., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

The Lol system in Escherichia coli is involved in localization of lipoproteins and hence is essential for growth of the organism. LolA is a periplasmic chaperone that binds to outer-membrane specific lipoproteins and transports them from inner membrane to outer membrane through LolB. The hydrophobic lipid-binding cavity of LolA consists of α-helices which act as a lid in regulating the transfer of lipoproteins from LolA to LolB. The current study aims to investigate the structural changes observed in LolA during the transition from open to closed conformation in the absence of lipoprotein. Molecular dynamics (MD) simulations were carried out for two LolA crystal structures; LolA(R43L), and in silico mutated MsL43R for a simulation time of 50 ns in water environment. We have performed an in silico point mutation of leucine to arginine in MsL43R to evaluate the importance of arginine to induce structural changes and impact the stability of protein structure. A complete dynamic analysis of open to closed conformation reveals the existence of two distinct levels; closing of lid and closing of entrance of hydrophobic cavity. Our analysis reveals that the structural flexibility of LolA is an important factor for its role as a periplasmic chaperone., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Protein-Protein Interactions (PPI) are vital to many cellular processes. The availability of high-throughput protein interaction data has provided us with an opportunity to assess domain associations in interacting proteins using computational approaches. High throughput PPI data, wherein the interaction status of every protein in the dataset has been experimentally tested against all the other proteins in the dataset contains information not only on protein interactions but also on proteins which do not interact with each other. We call such datasets "all against all" datasets. In the current study, using these datasets and the Pfam domain composition of the proteins in the sets, we have developed a matrix based method for predicting PPI. We infer positive and negative Domain-Domain Associations (DDA) by our method. We have generated more than a million domain association values which can be utilized for predicting new PPI. The performance of the algorithm was evaluated against a test set and the sensitivity and specificity was found to be 68.1% and 65.3%, respectively. The overall prediction accuracy of the algorithm with individual test sets from IntAct, DIP, 3did, iPfam databases and a literature curated set from Saccharomyces cerevisiae was found to be around 70%. The insights gained in the study have a potential application in providing leads for experimental interaction studies and understanding host pathogen interactions amongst others., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

With the current knowledge of the multitude of microbes that inhabit the human body, it is increasingly clear that they constitute an integral component of the host. The gut microbiota community is principally involved in the metabolism of dietary constituents such as carbohydrates which account for majority of the energy intake from diet. Diet has gained an important role in shaping the composition of gut inhabitants. The quantity and type of food consumed is recognized as a causal factor for metabolic disorders such as obesity and diabetes. Analysis of host-microbe interactions can thus contribute to the understanding of such metabolic disorders. In this study, data from Kyoto Encyclopedia of Genes and Genomes and Carbohydrate Active EnZYmes Database was utilized as a starting point. Enzyme information from the host Homo sapiens coupled with details of the three predominant phyla of gut bacteria, namely Firmicutes, Bacteroidetes and Actinobacteria were used in the creation of a comprehensive metabolic network, which we refer to as 'meta-metabolome'. This 'meta-metabolome' provides a perspective of the degree to which microbes influence carbohydrate metabolism, in conjunction with host specific enzymes. Analysis of reactions in the network reveals the amplification of monosaccharide content brought about by microbial enzyme activity. The framework outlined in this study provides a holistic approach to assess host-microbe symbiosis. It also provides us with a means of analyzing how diet can be modulated to provide beneficial effects to the host or how probiotics can potentially be used to relieve certain metabolic disorders., (Copyright © 2012 Elsevier Inc. All rights reserved.)

One of the challenges faced by Mycobacterium tuberculosis (M. tuberculosis) in dormancy is hypoxia. DosR/DevR of M. tuberculosis is a two component dormancy survival response regulator which induces the expression of 48 genes. In this study, we have used DosR regulon proteins of M. tuberculosis H37Rv as the query set and performed a comprehensive homology search against the non-redundant database. Homologs were found in environmental mycobacteria, environmental bacteria and archaebacteria. Analysis of genomic context of DosR regulon revealed that they are distributed as nine blocks in the genome of M. tuberculosis with many transposases and integrases in their vicinity. Further, we classified DosR regulon proteins into eight functional categories. One of the hypothetical proteins Rv1998c could probably be a methylisocitrate lyase or a phosphonomutase. Another hypothetical protein, Rv0572 was found only in mycobacteria. Insights gained in this study can potentially aid in the development of novel therapeutic interventions., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Context: Stress is a state of mental or emotional strain or tension, which can lead to underperformance and adverse clinical conditions. Adaptogens are herbs that help in combating stress. Ayurvedic classical texts, animal studies and clinical studies describe Ashwagandha as a safe and effective adaptogen., Aims: The aim of the study was to evaluate the safety and efficacy of a high-concentration full-spectrum extract of Ashwagandha roots in reducing stress and anxiety and in improving the general well-being of adults who were under stress., Settings and Design: Single center, prospective, double-blind, randomized, placebo-controlled trial., Materials and Methods: A total of 64 subjects with a history of chronic stress were enrolled into the study after performing relevant clinical examinations and laboratory tests. These included a measurement of serum cortisol, and assessing their scores on standard stress-assessment questionnaires. They were randomized to either the placebo control group or the study drug treatment group, and were asked to take one capsule twice a day for a period of 60 days. In the study drug treatment group, each capsule contained 300 mg of high-concentration full-spectrum extract from the root of the Ashwagandha plant. During the treatment period (on Day 15, Day 30 and Day 45), a follow-up telephone call was made to all subjects to check for treatment compliance and to note any adverse reactions. Final safety and efficacy assessments were done on Day 60., Statistical Analysis: t-test, Mann-Whitney test., Results: The treatment group that was given the high-concentration full-spectrum Ashwagandha root extract exhibited a significant reduction (P<0.0001) in scores on all the stress-assessment scales on Day 60, relative to the placebo group. The serum cortisol levels were substantially reduced (P=0.0006) in the Ashwagandha group, relative to the placebo group. The adverse effects were mild in nature and were comparable in both the groups. No serious adverse events were reported., Conclusion: The findings of this study suggest that a high-concentration full-spectrum Ashwagandha root extract safely and effectively improves an individual's resistance towards stress and thereby improves self-assessed quality of life.

Mammalian gastric lipases are stable and active under acidic conditions and also in the duodenal lumen. There has been considerable interest in acid stable lipases owing to their potential application in the treatment of pancreatic exocrine insufficiency. In order to gain insights into the domain movements of these enzymes, molecular dynamics simulations of human gastric lipase was performed at an acidic pH and under neutral conditions. For comparative studies, simulation of dog gastric lipase was also performed at an acidic pH. Analyses show, that in addition to the lid region, there is another region of high mobility in these lipases. The potential role of this novel region is discussed., (Copyright © 2010 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

Background and Purpose: Nitric oxide (NO) promotes angiogenesis by activating endothelial cells. Thalidomide arrests angiogenesis by interacting with the NO pathway, but its putative targets are not known. Here, we have attempted to identify these targets., Experimental Approach: Cell-based angiogenesis assays (wound healing of monolayers and tube formation in ECV304, EAhy926 and bovine arterial endothelial cells), along with ex vivo and in vivo angiogenesis assays, were used to explore interactions between thalidomide and NO. We also carried out in silico homology modelling and docking studies to elucidate possible molecular interactions of thalidomide and soluble guanylyl cyclase (sGC)., Key Results: Thalidomide inhibited pro-angiogenic functions in endothelial cell cultures, whereas 8-bromo-cGMP, sildenafil (a phosphodiesterase inhibitor) or a NO donor [sodium nitroprusside (SNP)] increased these functions. The inhibitory effects of thalidomide were reversed by adding 8-bromo-cGMP or sildenafil, but not by SNP. Immunoassays showed a concentration-dependent decrease of cGMP in endothelial cells with thalidomide, without affecting the expression level of sGC protein. These results suggested that thalidomide inhibited the activity of sGC. Molecular modelling and docking experiments revealed that thalidomide could interact with the catalytic domain of sGC, which would explain the inhibitory effects of thalidomide on NO-dependent angiogenesis., Conclusion and Implications: Our results showed that thalidomide interacted with sGC, suppressing cGMP levels in endothelial cells, thus exerting its anti-angiogenic effects. These results could lead to the formulation of thalidomide-based drugs to curb angiogenesis by targeting sGC.

Molecular Dynamics simulations were performed on aminoglycoside N-acetyltransferase ((AAC) (2')-Ic) of Mycobacterium tuberculosis, to gain insight into enzyme flexibility and conformation around the ligand and acetyl-CoA binding sites in nano scale. Simulations were performed in water and methanol to further study the effect of solvation on the enzyme. AAC(2')-Ic of M. tuberculosis consists of negatively charged residues of aspartic and glutamic acids, which are believed to form a docking platform for the positively charged aminoglycosides. The acetyl-CoA binding site involves a wide groove surrounded by helices alphal, alpha3, alpha4 and strand beta4. Simulation results revealed the enzyme to be stable in water and the enzyme was found to be flexible around the docking platform (ligand binding site) with Asp 35, Asp 40 and Asp 179 exhibiting maximum movement. An interesting observation was the flexibility and conformational change at beta4, hydrophobic pocket of the acetyl-CoA binding site a prerequisite for catalysis to occur with Leu 95 exhibiting maximum displacement with a root mean square fluctuation of 0.18 nm. The enzyme was found to be very unstable in methanol with root mean square deviation not stabilizing even at 10 ns.

Metal ion binding domains are found in proteins that mediate transport, buffering or detoxification of metal ions. In this study, we have performed an in silico analysis of metal binding proteins and have identified putative metal binding motifs for the ions of cadmium, cobalt, zinc, arsenic, mercury, magnesium, manganese, molybdenum and nickel. A pattern search against the UniProtKB/Swiss-Prot and UniProtKB/TrEMBL databases yielded true positives in each case showing the high-specificity of the motifs. Motifs were also validated against PDB structures and site directed mutagenesis studies.

BRCA1 and BRCA2 are the two major breast cancer susceptibility genes. We tested members of a family with multiple cases of breast cancer, for mutations in the BRCA1 gene. Analysis of the PCR amplicons of all the exons (22) of the BRCA1 gene using conformation sensitive gel electrophoresis (CSGE) revealed a heteroduplex band pattern in exon 2 of the proband (III-3) in this family. The amplicon was further sequenced to assess the nature of the mutation, which revealed a deletion of AG nucleotides at the 185th position (185delAG). The two base pair deletion introduces a stop codon at the 39th amino acid residue. A similar analysis was carried out on other extended family members to evaluate their allelic status. We detected the same deletion in 7 of the 19 members tested. Two of them are males. Haplotype analysis suggested an independent origin for this mutation. Our study highlights the importance of testing hereditary cases of breast/ovarian cancer for BRCA1 mutations in extended families in order to identify high-risk individuals at a pre-clinical stage and provide genetic counseling.

A novel methodology to predict the local conformational changes in a protein as a consequence of missense mutations is proposed. A pentapeptide at the locus of mutation plays the dominant role and it is analyzed in terms of tripeptides. A measure for spatial and temporal fluctuations in a pentapeptide is devised and validated. The method does not involve any prior knowledge of structural templates from sequence homology studies. Structural deformations can be predicted with about 70-80% reliability in any protein. Disease causing mutations and benign mutations have been addressed. In particular, p53, retinoblastoma protein and lipoprotein lipase are studied in detail.

The emergence of multidrug resistant varieties of Mycobacterium tuberculosis has led to a search for novel drug targets. We have performed an insilico comparative analysis of metabolic pathways of the host Homo sapiens and the pathogen M. tuberculosis. Enzymes from the biochemical pathways of M. tuberculosis from the KEGG metabolic pathway database were compared with proteins from the host H. sapiens, by performing a BLASTp search against the non-redundant database restricted to the H. sapiens subset. The e-value threshold cutoff was set to 0.005. Enzymes, which do not show similarity to any of the host proteins, below this threshold, were filtered out as potential drug targets. We have identified six pathways unique to the pathogen M. tuberculosis when compared to the host H. sapiens. Potential drug targets from these pathways could be useful for the discovery of broad spectrum drugs. Potential drug targets were also identified from pathways related to lipid metabolism, carbohydrate metabolism, amino acid metabolism, energy metabolism, vitamin and cofactor biosynthetic pathways and nucleotide metabolism. Of the 185 distinct targets identified from these pathways, many are in various stages of progress at the TB Structural Genomics Consortium. However, 67 of our targets are new and can be considered for rational drug design. As a case study, we have built a homology model of one of the potential drug targets MurD ligase using WHAT IF software. The model could be further explored for insilico docking studies with suitable inhibitors. The study was successful in listing out potential drug targets from the M. tuberculosis proteome involved in vital aspects of the pathogen's metabolism, persistence, virulence and cell wall biosynthesis. This systematic evaluation of metabolic pathways of host and pathogen through reliable and conventional bioinformatic methods can be extended to other pathogens of clinical interest.

Pseudomonas aeruginosa lipase is a 29-kDa protein that, following the determination of its crystal structure, was postulated to have a lid that stretched between residues 125 and 148. In this paper, using molecular dynamics simulations, we propose that there exists, in addition to the above-mentioned lid, a novel second lid in this lipase. We further show that the second lid, covering residues 210-222, acts as a triggering lid for the movement of the first. We also investigate the role of hydrophobicity in the movement of the lids and show that two residues, Phe214 and Ala217, play important roles in lid movement. To our knowledge, this is the first time that a double-lid movement of the type described in our manuscript has been presented to the scientific community. This work also elucidates the interplay of hydrophobic interactions in the dynamics, and hence the function, of an enzyme.

Purpose: To study some functional candidate genes in cataract families of Indian descent., Methods: Nine Indian families, clinically documented to have congenital/childhood cataracts, were screened for mutations in candidate genes such as CRYG (A-->D), CRYBB2, and GJA8 by PCR analyses and sequencing. Genomic DNA samples of either probands or any representative affected member of each family were PCR amplified and sequenced commercially. Documentation of single nucleotide polymorphisms (SNPs) and candidate mutations was done through BLAST SEARCH (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi?)., Results: Several single nucleotide polymorphisms in CRYG, CRYBB2, and GJA8 genes were observed. Because they do not co-segregate with the phenotype, they were excluded as candidates for the cataract formation in these patients. However, a substitution (W151C in exon 6 of CRYBB2) was identified as the most likely causative mutation underlying the phenotype of central nuclear cataract in all affected members of family C176. Protein structural interpretations demonstrated that no major structural alterations could be predicted and that even the hydrogen bonds to the neighboring Leu166 were unchanged. Surprisingly, hydropathy analysis of the mutant betaB2-crystallin featuring the amino acids at position 147 to 155, further increased the hydrophobicity, which might impair the solubility of the mutant protein. Finally, the Cys residue at position 151 might possibly be involved in intramolecular disulphide bridges with other cysteines during translation, possibly leading to dramatic structural changes., Conclusions: Exon 6 of CRYBB2 appears to be a critical region susceptible for mutations leading to lens opacity.

Background: Proteins having similar functions from different sources can be identified by the occurrence in their sequences, a conserved cluster of amino acids referred to as pattern, motif, signature or fingerprint. The wide usage of protein sequence analysis in par with the growth of databases signifies the importance of using patterns or signatures to retrieve out related sequences. Blue copper proteins are found in the electron transport chain of prokaryotes and eukaryotes. The signatures already existing in the databases like the type 1 copper blue, multiple copper oxidase, cyt b/b6, photosystem 1 psaA&B, psaG&K, and reiske iron sulphur protein are not specified signatures for blue copper proteins as the name itself suggests. Most profile and motif databases strive to classify protein sequences into a broad spectrum of protein families. This work describes the signatures designed based on the copper metal binding motifs in blue copper proteins. The common feature in all blue copper proteins is a trigonal planar arrangement of two nitrogen ligands [each from histidine] and one sulphur containing thiolate ligand [from cysteine], with strong interactions between the copper center and these ligands., Results: Sequences that share such conserved motifs are crucial to the structure or function of the protein and this could provide a signature of family membership. The blue copper proteins chosen for the study were plantacyanin, plastocyanin, cucumber basic protein, stellacyanin, dicyanin, umecyanin, uclacyanin, cusacyanin, rusticyanin, sulfocyanin, halocyanin, azurin, pseudoazurin, amicyanin and nitrite reductase which were identified in both eukaryotes and prokaryotes. ClustalW analysis of the protein sequences of each of the blue copper proteins was the basis for designing protein signatures or peptides. The protein signatures and peptides identified in this study were designed involving the active site region involving the amino acids bound to the copper atom. It was highly specific for each kind of blue copper protein and the false picks were minimized. The set of signatures designed specifically for the BCP's was entirely different from the existing broad spectrum signatures as mentioned in the background section., Conclusions: These signatures can be very useful for the annotation of uncharacterized proteins and highly specific to retrieve blue copper protein sequences of interest from the non redundant databases containing a large deposition of protein sequences.

Background: An efficient building block for protein structure prediction can be tripeptides. 8000 different tripeptides from a dataset of 1220 high resolution (