Local sequential minimization of double stranded B-DNA using Monte Carlo annealing

A software algorithm has been developed to investigate the folding process in B-DNA structures in vacuum under a simple and accurate force field. This algorithm models linear double stranded B-DNA sequences based on a local, sequential minimization procedure. The original B-DNA structures were modeled using initial nucleotide structures taken from the Brookhaven database. The models contain information at the atomic level allowing one to investigate as accurately as possible the structure and characteristics of the resulting DNA structures. A variety of DNA sequences and sizes were investigated containing coding and non-coding, random and real, homogeneous or heterogeneous sequences in the range of 2 to 40 base pairs. The force field contains terms such as angle bend, Lennard-Jones, electrostatic interactions and hydrogen bonding which are set up using the Dreiding II force field and defined to account for the helical parameters such as twist, tilt and rise. A close comparison was made between this local minimization algorithm and a global one (previously published) in order to find out advantages and disadvantages of the different methods. From the comparison, this algorithm gives better and faster results than the previous method, allowing one to minimize larger DNA segments. DNA segments with a length of 40 bases need approximately 4 h, while 2.5 weeks are needed with the previous method. After each minimization the angles between phosphate–oxygen-carbon A₁, the oxygen–phosphate–oxygen A₂ and the average helical twists were calculated. From the generated fragments it was found that the bond angles are A₁=150°±2°and A₂=130°±10°, while the helical twist is 36.6°±2° in the A strand and A₁=150°±6° and A₂=130±6° with helical twist 39.6°±2° in the B strand for the DNA segment with the same sequence as the Dickerson dodecamer.Figure The final minimized DNA segment of the Dickerson dodecamer sequence represented by ball drawings and viewed (left) perpendicular and (right) down the helical axis     

Glutathione S-transferase in the insect Apis mellifera macedonica kinetic characteristics and effect of stress on the expression of GST isoenzymes in the adult worker bee

The glutathione S-transferase present in the adult worker bee Apis mellifera macedonica was purified and analyzed for its physicochemical and kinetic properties. The enzyme is heterodimeric with subunit molecular masses of 29 and 25 kDa, respectively. Two main isoenzymes with distinct kinetic properties are present, with isoelectric points of 7.40 for the alkaline and 4.58 for the acidic forms, respectively. The two enzymes are induced independently by factors such as insecticide treatments and environmental conditions, including low temperatures or starvation.     

BAK and NOXA Are Critical Determinants of Mitochondrial Apoptosis Induced by Bortezomib in Mesothelioma

Based on promising preclinical efficacy associated with the 20S proteasome inhibitor bortezomib in malignant pleural mesothelioma (MPM), two phase II clinical trials have been initiated (EORTC 08052 and ICORG 05–10). However, the potential mechanisms underlying resistance to this targeted drug in MPM are still unknown. Functional genetic analyses were conducted to determine the key mitochondrial apoptotic regulators required for bortezomib sensitivity and to establish how their dysregulation may confer resistance. The multidomain proapoptotic protein BAK, but not its orthologue BAX, was found to be essential for bortezomib-induced apoptosis in MPM cell lines. Immunohistochemistry was performed on tissues from the ICORG-05 phase II trial and a TMA of archived mesotheliomas. Loss of BAK was found in 39% of specimens and loss of both BAX/BAK in 37% of samples. However, MPM tissues from patients who failed to respond to bortezomib and MPM cell lines selected for resistance to bortezomib conserved BAK expression. In contrast, c-Myc dependent transactivation of NOXA was abrogated in the resistant cell lines. In summary, the block of mitochondrial apoptosis is a limiting factor for achieving efficacy of bortezomib in MPM, and the observed loss of BAK expression or NOXA transactivation may be relevant mechanisms of resistance in the clinic.     

The “clustered structure” of the purines/pyrimidines distribution in DNA distinguishes systematically between coding and non-coding sequences

A method allowing to measure the inhomogeneous distribution of purines/pyrimidines in nucleotide sequences is developed. We show that this measure relates to the coding or non-coding character of the considered sequence. Coding sequences present a near to the random Pu or Py distribution. This property is shared by both protein-coding DNA and functional RNA-coding DNA. Non-coding sequences present a highly clustered inhomogeneity. We propose the hypothesis, corroborated with appropriate computer simulations, that this is due to the action of various transposition events accumulated for long time periods.     

Autophagy Correlates with the Therapeutic Responsiveness of Malignant Pleural Mesothelioma in 3D Models

Malignant pleural mesothelioma is a highly chemoresistant solid tumor. We have studied this apoptotic resistance using in vitro and ex vivo three-dimensional models, which acquire a high level of chemoresistance that can be reduced by PI3K/mTOR inhibitors. Here, we investigate the activity of GDC-0980, a novel dual PI3K/mTOR inhibitor, which has been proposed to be effective in mesothelioma. In this work, we aimed to identify mechanisms and markers of efficacy for GDC-0980 by utilizing 3D models of mesothelioma, both in vitro multicellular spheroids and ex vivo tumor fragment spheroids grown from patient tumor samples. We found that a subset of mesothelioma spheroids is sensitive to GDC-0980 alone and to its combination with chemotherapy. Unexpectedly, this sensitivity did not correlate with the activation of the Akt/mTOR pathway. Instead, sensitivity to GDC-0980 correlated with the presence of constitutive ATG13 puncta, a feature of autophagy, a cellular program that supports cells under stress. In tumor fragment spheroids grown from 21 tumors, we also found a subset (n = 11) that was sensitive to GDC-0980, a sensitivity that also correlated with the presence of ATG13 puncta. Interference with autophagy by siRNA of ATG7, an essential autophagic protein, increased the response to chemotherapy, but only in the sensitive multicellular spheroids. In the spheroids resistant to GDC-0980, autophagy appeared to play no role. In summary, we show that GDC-0980 is effective in mesothelioma 3D models that display ATG13 puncta, and that blockade of autophagy increases their response to chemotherapy. For the first time, we show a role for autophagy in the response to chemotherapy of 3D models of mesothelioma and propose ATG13 as a potential biomarker of the therapeutic responsiveness of mesothelioma.     

Glutathione S-transferase in the developmental stages of the insect Apis mellifera macedonica

We investigated the pattern of glutathione S-transferase (GST) activity in the course of the development of Apis mellifera macedonica. GST activity is present in all developmental stages of A. mellifera macedonica. The highest activity towards the substrate 1-chloro-2,4-dinitrobenzene (CDNB) is found in the adult stage and the lowest in the egg. The kinetic characteristics of the whole enzyme change as the insect develops. Significant changes are observed in substrate specificity, inhibitor sensitivity and V(max). The number of isoenzymes and their rate of expression vary as the insect develops. However, two main isoenzymes are present in all developmental stages, one in the alkaline area and the other in the acidic. While in the larval stage the acidic isoenzyme is expressed at a slightly higher rate (52.2% over 47.8% for the alkaline isoenzyme), in the adult stage, the rate is reversed dramatically (13.24% and 84.2%, respectively).     

Evolutionary Constraints Favor a Biophysical Model Explaining Hox Gene Collinearity

The Hox gene collinearity enigma has often been approached using models based on biomolecular mechanisms. The biophysical model is an alternative approach based on the hypothesis that collinearity is caused by physical forces pulling the Hox genes from a territory where they are inactive to a distinct spatial domain where they are activated in a step by step manner. Such Hox gene translocations have recently been observed in support of the biophysical model. Genetic engineering experiments, performed on embryonic mice, gave rise to several unexpected mutant expressions that the biomolecular models cannot predict. On the contrary, the biophysical model offers convincing explanation. Evolutionary constraints consolidate the Hox clusters and as a result, denser and well organized clusters may create more efficient physical forces and a more emphatic manifestation of gene collinearity. This is demonstrated by stochastic modeling with white noise perturbing the expression of Hox genes. As study cases the genomes of mouse and amphioxus are used. The results support the working hypothesis that vertebrates have adopted their comparably more compact Hox clustering as a tool needed to develop more complex body structures. Several experiments are proposed in order to test further the physical forces hypothesis.     

Alu and LINE1 distributions in the human chromosomes: evidence of global genomic organization expressed in the form of power laws

Spatial distribution and clustering of repetitive elements are extensively studied during the last years, as well as their colocalization with other genomic components. Here we investigate the large-scale features of Alu and LINE1 spatial arrangement in the human genome by studying the size distribution of interrepeat distances. In most cases, we have found power-law size distributions extending in several orders of magnitude. We have also studied the correlations of the extent of the power law (linear region in double-logarithmic scale) and of the corresponding exponent (slope) with other genomic properties. A model has been formulated to explain the formation of the observed power laws. According to the model, 2 kinds of events occur repetitively in evolutionary time: random insertion of several types of intruding sequences and occasional loss of repeats belonging to the initial population due to "elimination" events. This simple mechanism is shown to reproduce the observed power-law size distributions and is compatible with our present knowledge on the dynamics of repeat proliferation in the genome.