Amino acid interaction preferences in helical membrane proteins

Membrane proteins are involved in a number of important biological functions. Yet, they are poorly understood from the structure and folding point of view. The external environment being drastically different from that of globular proteins, the intra-protein interactions in membrane proteins are also expected to be different. Hence, statistical potentials representing the features of inter-residue interactions based exclusively on the structures of membrane proteins are much needed. Currently, a reasonable number of structures are available, making it possible to undertake such an analysis on membrane proteins. In this study we have examined the inter-residue interaction propensities of amino acids in the membrane spanning regions of the alpha-helical membrane (HM) proteins. Recently we have shown that valuable information can be obtained on globular proteins by the evaluation of the pair-wise interactions of amino acids by classifying them into different structural environments, based on factors such as the secondary structure or the number of contacts that a residue can make. Here we have explored the possible ways of classifying the intra-protein environment of HM proteins and have developed scoring functions based on different classification schemes. On evaluation of different schemes, we find that the scheme which classifies amino acids to different intra-contact environment is the most promising one. Based on this classification scheme, we also redefine the hydrophobicity scale of amino acids in HM proteins.

Identity, energetics, dynamics and environment of interfacial water molecules in a micellar solution

The structure and energetics of interfacial water molecules in the aqueous micelle of cesium perfluorooctanoate have been investigated, using large-scale atomistic molecular dynamics simulations, with the primary objective of classifying them. The simulations show that the water molecules at the interface fall into two broad classes: bound and free, present in a ratio of 9:1. The bound water molecules can be further categorized on the basis of the number of hydrogen bonds (one or two) that they form with the surfactant headgroups. The hydrogen bonds of the doubly hydrogen-bonded species are found to be, on the average, slightly weaker than those in the singly bonded species. The environment around interfacial water molecules is more ordered than that in the bulk. The surface water molecules have substantially lower potential energy, because of interaction with the micelle. In particular, both forms of bound water have energies that are lower by �2.5-4.0 kcal/ mol. Entropy is found to play an important role in determining the relative concentration of the species.

NETGEM: Network Embedded Temporal GEnerative Model for gene expression data

Background: Temporal analysis of gene expression data has been limited to identifying genes whose expression varies with time and/or correlation between genes that have similar temporal profiles. Often, the methods do not consider the underlying network constraints that connect the genes. It is becoming increasingly evident that interactions change substantially with time. Thus far, there is no systematic method to relate the temporal changes in gene expression to the dynamics of interactions between them. Information on interaction dynamics would open up possibilities for discovering new mechanisms of regulation by providing valuable insight into identifying time-sensitive interactions as well as permit studies on the effect of a genetic perturbation. Results: We present NETGEM, a tractable model rooted in Markov dynamics, for analyzing the dynamics of the interactions between proteins based on the dynamics of the expression changes of the genes that encode them. The model treats the interaction strengths as random variables which are modulated by suitable priors. This approach is necessitated by the extremely small sample size of the datasets, relative to the number of interactions. The model is amenable to a linear time algorithm for efficient inference. Using temporal gene expression data, NETGEM was successful in identifying (i) temporal interactions and determining their strength, (ii) functional categories of the actively interacting partners and (iii) dynamics of interactions in perturbed networks. Conclusions: NETGEM represents an optimal trade-off between model complexity and data requirement. It was able to deduce actively interacting genes and functional categories from temporal gene expression data. It permits inference by incorporating the information available in perturbed networks. Given that the inputs to NETGEM are only the network and the temporal variation of the nodes, this algorithm promises to have widespread applications, beyond biological systems. The source code for NETGEM is available from https://github.com/vjethava/NETGEM

Interaction of G-Quadruplexes with Nonintercalating Duplex-DNA Minor Groove Binding Ligands

The enzyme telomerase synthesizes the G-rich DNA strands of the telomere and its activity is often associated with cancer. The telomerase may be therefore responsible for the ability of a cancer cell-to escape apoptosis. The G-rich DNA sequences often adopt tetra-stranded structure, known as the G-quadruplex DNA (G4-DNA). The stabilization of the telomeric DNA into the G4-DNA structures by small molecules has been the focus of many researchers for the design and development of new anticancer agents. The compounds which stabilize the G-quadruplex in the telomere inhibit the telomerase activity. Besides telomeres, the G4-DNA forming sequences are present in the genomic regions of biological significance including the transcriptional regulatory and promoter regions of several oncogenes. Inducing a G-quadruplex structure within the G-rich promoter sequences is a potential way of achieving selective gene regulation. Several G-quadruplex stabilizing ligands are known. Minor groove binding ligands (MGBLs) interact with the double-helical DNA through the minor grooves sequence-specifically and interfere with several DNA associated processes. These MGBLs when suitably modified switch their preference sometimes from the duplex DNA to G4-DNA and stabilize the G4-DNA as well. Herein, we focus on the recent advances in understanding the G-quadruplex structures, particularly made by the human telomeric ends, and review the results of various investigations of the interaction of designed organic ligands with the G-quadruplex DNA while highlighting the importance of MGBL-G-quadruplex interactions.

Nonspecific Interaction between DNA and Protein allows for Cooperativity: A Case Study with Mycobacterium DNA Binding Protein

Different DNA-binding proteins have different interaction modes with DNA. Sequence-specific DNA protein interaction has been mostly associated with regulatory processes inside a cell, and as such extensive studies have been made. Adequate data is also available on nonspecific DNA protein interaction, as an intermediate to protein's search for its cognate partner. Multidomain nonspecific DNA protein interaction involving physical sequestering of DNA has often been implicated to regulate gene expression indirectly. However, data available on this type of interaction is limited. One such interaction is the binding of DNA with mycobacterium DNA binding proteins. We have used the Langmuir-Blodgett technique to evaluate for the first time the kinetics and thermodynamics of Mycobacterium smegmatis Dps 1 binding to DNA. By immobilizing one of the interacting partners, we have shown that, when a kinetic bottleneck is applied, the binding mechanism showed cooperative binding (n = 2.72) at lower temperatures, but the degree of cooperativity gradually reduces (n = 1.38) as the temperature was increased We have also compared the kinetics and thermodynamics of sequence-specific and nonspecific DNA protein interactions under the same set of conditions.

Interaction of nucleic acids with carbon nanotubes and dendrimers

Nucleic acid interaction with nanoscale objects like carbon nanotubes (CNTs) and dendrimers is of fundamental interest because of their potential application in CNT separation, gene therapy and antisense therapy. Combining nucleic acids with CNTs and dendrimers also opens the door towards controllable self-assembly to generate various supra-molecular and nano-structures with desired morphologies. The interaction between these nanoscale objects also serve as a model system for studying DNA compaction, which is a fundamental process in chromatin organization. By using fully atomistic simulations, here we report various aspects of the interactions and binding modes of DNA and small interfering RNA (siRNA) with CNTs, graphene and dendrimers. Our results give a microscopic picture and mechanism of the adsorption of single- and double-strand DNA (ssDNA and dsDNA) on CNT and graphene. The nucleic acid-CNT interaction is dominated by the dispersive van der Waals (vdW) interaction. In contrast, the complexation of DNA (both ssDNA and dsDNA) and siRNA with various generations of poly-amido-amine (PAMAM) dendrimers is governed by electrostatic interactions. Our results reveal that both the DNA and siRNA form stable complex with the PAMAM dendrimer at a physiological pH when the dendrimer is positively charged due to the protonation of the primary amines. The size and binding energy of the complex increase with increase in dendrimer generation. We also give a summary of the current status in these fields and discuss future prospects.

A Comparative Kinetic and Thermodynamic Perspective of the sigma-Competition Model in Escherichia coli

Transcription is the most fundamental step in gene expression in any living organism. Various environmental cues help in the maturation of core RNA polymerase (RNAP; alpha(2)beta beta'omega) with different sigma-factors, leading to the directed recruitment of RNAP to different promoter DNA sequences. Thus it is essential to determine the sigma-factors that affect the preferential partitioning of core RNAP among various a-actors, and the role of sigma-switching in transcriptional gene regulation. Further, the macromolecular assembly of holo RNAP takes place in an extremely crowded environment within a cell, and thus far the kinetics and thermodynamics of this molecular recognition process have not been well addressed. In this study we used a site-directed bioaffinity immobilization method to evaluate the relative binding affinities of three different Escherichia coli sigma-factors to the same core RNAP with variations in temperature and ionic strength while emulating the crowded cellular milieu. Our data indicate that the interaction of core RNAP-sigma is susceptible to changes in external stimuli such as osmolytic and thermal stress, and the degree of susceptibility varies among different sigma-factors. This allows for a reversible sigma-switching from housekeeping factors to alternate sigma-factors when the organism senses a change in its physiological conditions.

MODELING AND SIMULATION OF HYDRODYNAMIC INTERACTION OF DNA IN A MICRO-FLUIDIC CHANNEL

The motion of DNA (in the bulk solution) and the non-Newtonian effective fluid behavior are considered separately and self-consistently with the fluid motion satisfying the no-slip boundary condition on the surface of the confining geometry in the presence of channel pressure gradients. A different approach has been developed to model DNA in the micro-channel. In this study the DNA is assumed as an elastic chain with its characteristic Young's modulus, Poisson's ratio and density. The force which results from the fluid dynamic pressure, viscous forces and electromotive forces is applied to the elastic chain in a coupled manner. The velocity fields in the micro-channel are influenced by the transport properties. Simulations are carried out for the DNAs attached to the micro-fluidic wall. Numerical solutions based on a coupled multiphysics finite element scheme are presented. The modeling scheme is derived based on mass conservation including biomolecular mass, momentum balance including stress due to Coulomb force field and DNA-fluid interaction, and charge transport associated to DNA and other ionic complexes in the fluid. Variation in the velocity field for the non-Newtonian flow and the deformation of the DNA strand which results from the fluid-structure interaction are first studied considering a single DNA strand. Motion of the effective center of mass is analyzed considering various straight and coil geometries. Effects of DNA statistical parameters (geometry and spatial distribution of DNAs along the channel) on the effective flow behavior are analyzed. In particular, the dynamics of different DNA physical properties such as radius of gyration, end-to-end length etc. which are obtained from various different models (Kratky-Porod, Gaussian bead-spring etc.) are correlated to the nature of interaction and physical properties under the same background fluid environment.

Quorum sensing and pathogenesis: role of small signalling molecules in bacterial persistence

The pathogenesis of Mycobacterium tuberculosis is associated with its ability to survive inside the human host and the bacteria use a variety of mechanism to evade the host's defence. A clearer understanding of the host pathogen interaction is needed to follow the pathogenicity and virulence. Recent advances in the study of inter and intra-cellular communication in bacteria had prompted us to study the role of quorum sensing in bacterial survival and pathogenicity. The cell cell communication in bacteria (quorum sensing) is mediated through the exchange of small molecules called as autoinducers that allow bacteria to modulate their gene expression in response to change in cell-population density. It is a coordinated response that confers multicellularity to a bacterial population in response to stress from external environment. Quorum sensing molecules are the global regulators and regulate a wide range of physiological processes including biofilm formation, motility, cell differentiation, long-term survival and many others. Many bacterial pathogens require quorum sensing to produce the virulence factors in response to host pathogen interaction. Here, we summarize our current understanding on small molecule signalling and their role in the bacterial persistence. New discoveries in these areas have enriched our knowledge on intracellular signalling and their role in the long-term survival of mycobacteria under nutrient starvation.

DNA STRUCTURAL FEATURES AND ARCHITECTURE OF PROMOTER REGIONS PLAY A ROLE IN GENE RESPONSIVENESS OF S. CEREVISIAE

Gene expression is the most fundamental biological process, which is essential for phenotypic variation. It is regulated by various external (environment and evolution) and internal (genetic) factors. The level of gene expression depends on promoter architecture, along with other external factors. Presence of sequence motifs, such as transcription factor binding sites (TFBSs) and TATA-box, or DNA methylation in vertebrates has been implicated in the regulation of expression of some genes in eukaryotes, but a large number of genes lack these sequences. On the other hand, several experimental and computational studies have shown that promoter sequences possess some special structural properties, such as low stability, less bendability, low nucleosome occupancy, and more curvature, which are prevalent across all organisms. These structural features may play role in transcription initiation and regulation of gene expression. We have studied the relationship between the structural features of promoter DNA, promoter directionality and gene expression variability in S. cerevisiae. This relationship has been analyzed for seven different measures of gene expression variability, along with two different regulatory effect measures. We find that a few of the variability measures of gene expression are linked to DNA structural properties, nucleosome occupancy, TATA-box presence, and bidirectionality of promoter regions. Interestingly, gene responsiveness is most intimately correlated with DNA structural features and promoter architecture.

Efficient Dendrimer-DNA Complexation and Gene Delivery Vector Properties of Nitrogen-Core Poly(propyl ether imine) Dendrimer in Mammalian Cells

Dendrimers as vectors for gene delivery were established, primarily by utilizing few prominent dendrimer types so far. We report herein studies of DNA complexation efficacies and gene delivery vector properties of a nitrogen-core poly(propyl ether imine) (PETIM) dendrimer, constituted with 22 tertiary amine internal branches and 24 primary amines at the periphery. The interaction of the dendrimer with pEGFPDNA was evaluated through UV-vis, circular dichroism (CD) spectral studies, ethidium bromide fluorescence emission quenching, thermal melting, and gel retardation assays, from which most changes to DNA structure during complexation was found to occur at a weight ratio of dendrimer:DNA similar to 2:1. The zeta potential measurements further confirmed this stoichiometry at electroneutrality. The structure of a DNA oligomer upon dendrimer complexation was simulated through molecular modeling and the simulation showed that the dendrimer enfolded DNA oligomer along both major and minor grooves, without causing DNA deformation, in 1:1 and 2:1 dendrimer-to-DNA complexes. Atomic force microscopy (AFM) studies on dendrimer-pEGFP DNA complex showed an increase in the average z-height as a result of dendrimers decorating the DNA, without causing a distortion of the DNA structure. Cytotoxicity studies involving five different mammalian cell lines, using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide] (MTT) assay, reveal the dendrimer toxicity profile (IC50) values of similar to 400-1000 mu g mL(-1), depending on the cell line tested. Quantitative estimation, using luciferase assay, showed that the gene transfection was at least 100 times higher when compared to poly(ethylene imine) branched polymer, having similar number of cationic sites as the dendrimer. The present study establishes the physicochemical behavior of new nitrogen-core PETIM dendrimer-DNA complexes, their lower toxicities, and efficient gene delivery vector properties.

A 16-Gene Signature Distinguishes Anaplastic Astrocytoma from Glioblastoma

Anaplastic astrocytoma (AA; Grade III) and glioblastoma (GBM; Grade IV) are diffusely infiltrating tumors and are called malignant astrocytomas. The treatment regimen and prognosis are distinctly different between anaplastic astrocytoma and glioblastoma patients. Although histopathology based current grading system is well accepted and largely reproducible, intratumoral histologic variations often lead to difficulties in classification of malignant astrocytoma samples. In order to obtain a more robust molecular classifier, we analysed RT-qPCR expression data of 175 differentially regulated genes across astrocytoma using Prediction Analysis of Microarrays (PAM) and found the most discriminatory 16-gene expression signature for the classification of anaplastic astrocytoma and glioblastoma. The 16-gene signature obtained in the training set was validated in the test set with diagnostic accuracy of 89%. Additionally, validation of the 16-gene signature in multiple independent cohorts revealed that the signature predicted anaplastic astrocytoma and glioblastoma samples with accuracy rates of 99%, 88%, and 92% in TCGA, GSE1993 and GSE4422 datasets, respectively. The protein-protein interaction network and pathway analysis suggested that the 16-genes of the signature identified epithelial-mesenchymal transition (EMT) pathway as the most differentially regulated pathway in glioblastoma compared to anaplastic astrocytoma. In addition to identifying 16 gene classification signature, we also demonstrated that genes involved in epithelial-mesenchymal transition may play an important role in distinguishing glioblastoma from anaplastic astrocytoma.

Lactoylglutathione lyase, a critical enzyme in methylglyoxal detoxification, contributes to survival of Salmonella in the nutrient rich environment

Glyoxalase I which is synonymously known as lactoylglutathione lyase is a critical enzyme in methylglyoxal (MG) detoxification. We assessed the STM3117 encoded lactoylglutathione lyase (Lgl) of Salmonella Typhimurium, which is known to function as a virulence factor, due in part to its ability to detoxify methylglyoxal. We found that STM3117 encoded Lgl isomerises the hemithioacetal adduct of MG and glutathione (GSH) into S-lactoylglutathione. Lgl was observed to be an outer membrane bound protein with maximum expression at the exponential growth phase. The deletion mutant of S. Typhimurium (lgl) exhibited a notable growth inhibition coupled with oxidative DNA damage and membrane disruptions, in accordance with the growth arrest phenomenon associated with typical glyoxalase I deletion. However, growth in glucose minimal medium did not result in any inhibition. Endogenous expression of recombinant Lgl in serovar Typhi led to an increased resistance and growth in presence of external MG. Being a metalloprotein, Lgl was found to get activated maximally by Co2+ ion followed by Ni2+, while Zn2+ did not activate the enzyme and this could be attributed to the geometry of the particular protein-metal complex attained in the catalytically active state. Our results offer an insight on the pivotal role of the virulence associated and horizontally acquired STM3117 gene in non-typhoidal serovars with direct correlation of its activity in lending survival advantage to Salmonella spp.

Rrp1B gene polymorphism (1307T > C) in metastatic progression of breast cancer

Rrp1B (ribosomal RNA processing1 homolog B) is a novel candidate metastasis modifier gene in breast cancer. Functional gene assays demonstrated that a physical and functional interaction existing between Rrp1b and metastasis modifier gene SIPA1 causes reduction in the tumor growth and metastatic potential. Ectopic expression of Rrp1B modulates various metastasis predictive extra cellular matrix (ECM) genes associated with tumor suppression. The aim of this study is to determine the functional significance of single nucleotide polymorphism (SNP) in human Rrp1B gene (1307 T > C; rs9306160) with breast cancer development and progression. The study consists of 493 breast cancer cases recruited from Nizam's Institute of Medical Sciences, Hyderabad, and 558 age-matched healthy female controls from rural and urban areas. Genomic DNA was isolated by non-enzymatic method. Genotyping was done by amplification refractory mutation system (ARMS-PCR) method. Genotypes were reconfirmed by sequencing and results were analyzed statistically. We have performed Insilco analysis to know the RNA secondary structure by using online tool m fold. The TT genotype and T allele frequencies of Rrp1B1307 T > C polymorphism were significantly elevated in breast cancer (chi (2); p = < 0.008) cases compared to controls under different genetic models. The presence of T allele had conferred 1.75-fold risk for breast cancer development (OR = 1.75; 95 % CI = 1.15-2.67). The frequency of TT genotype of Rrp1b 1307T > C polymorphism was significantly elevated in obese patients (chi (2); p = 0.008) and patients with advanced disease (chi (2); p = 0.01) and with increased tumor size (chi (2); p = 0.01). Moreover, elevated frequency of T allele was also associated with positive lymph node status (chi (2); p = 0.04) and Her2 negative receptor status (chi (2); p = 0.006). Presence of Rrp1b1307TT genotype and T allele confer strong risk for breast cancer development and progression.

Co-liposomes of redox-active alkyl-ferrocene modified low MW branched PEI and DOPE for efficacious gene delivery in serum

Herein, we present six new lipopolymers based on low molecular weight, branched polyethylenimine (BPEI 800 Da) which are hydrophobically modified using ferrocene terminated alkyl tails of variable lengths. The effects of degree of grafting, spacer length and the redox state of ferrocene in the lipopolymers on the self assembly properties were investigated in detail by TEM, AFM, DLS and zeta potential measurements. The assemblies displayed an oxidation induced increase in the size of the aggregates. The co-liposomes comprising the lipopolymer and a helper lipid, 1,2-dioleoyl phosphatidyl ethanolamine (DOPE), showed excellent gene (pDNA) delivery capability in a serum containing environment in two cancer cell lines (HeLa and U251 cells). Optimized formulations showed remarkably higher transfection activity than BPEI (25 kDa) and were also significantly better than a commercial transfection reagent, Lipofectamine 2000 as evidenced from both the luciferase activity and GFP expression analysis. Oxidation of ferrocene in the lipopolymers led to drastically reduced levels of gene transfection which was substantiated by reduced cellular internalization of fluorescently labelled pDNA as detected using confocal microscopy and flow cytometry. Moreover, the transfection inactive oxidized lipopolyplexes could be turned transfection active by exposure to ascorbic acid (AA) in cell culture medium during transfection. Endocytosis inhibition experiments showed that gene expression mediated by reduced formulations involved both clathrin and caveolae mediated pathways while the oxidized formulations were routed via the caveolae. Cytotoxicity assays revealed no obvious toxicity for the lipopolyplexes in the range of optimized transfection levels in both the cell lines studied. Overall, we have exploited the redox activity of ferrocene in branched PEI-based efficient polymeric gene carriers whose differential transfection activities could be harnessed for spatial or temporal cellular transfections.