Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

DNA-Interactive Properties of Crotamine, a Cell-Penetrating Polypeptide and a Potential Drug Carrier

Crotamine, a 42-residue polypeptide derived from the venom of the South American rattlesnake Crotalus durissus terrificus, has been shown to be a cell-penetrating protein that targets chromosomes, carries plasmid DNA into cells, and shows specificity for actively proliferating cells. Given this potential role as a nucleic acid-delivery vector, we have studied in detail the binding of crotamine to single- and double-stranded DNAs of different lengths and base compositions over a range of ionic conditions. Agarose gel electrophoresis and ultraviolet spectrophotometry analysis indicate that complexes of crotamine with long-chain DNAs readily aggregate and precipitate at low ionic strength. This aggregation, which may be important for cellular uptake of DNA, becomes less likely with shorter chain length. 25-mer oligonucleotides do not show any evidence of such aggregation, permitting the determination of affinities and size via fluorescence quenching experiments. The polypeptide binds non-cooperatively to DNA, covering about 5 nucleotide residues when it binds to single (ss) or (ds) double stranded molecules. The affinities of the protein for ss-vs. ds-DNA are comparable, and inversely proportional to salt levels. Analysis of the dependence of affinity on [NaCl] indicates that there are a maximum of,3 ionic interactions between the protein and DNA, with some of the binding affinity attributable to non-ionic interactions. Inspection of the three-dimensional structure of the protein suggests that residues 31 to 35, Arg-Trp-Arg-Trp-Lys, could serve as a potential DNA-binding site. A hexapeptide containing this sequence displayed a lower DNA binding affinity and salt dependence as compared to the full-length protein, likely indicative of a more suitable 3D structure and the presence of accessory binding sites in the native crotamine. Taken together, the data presented here describing crotamine-DNA interactions may lend support to the design of more effective nucleic acid drug delivery vehicles which take advantage of crotamine as a carrier with specificity for actively proliferating cells. Citation: Chen P-C, Hayashi MAF, Oliveira EB, Karpel RL (2012) DNA-Interactive Properties of Crotamine, a Cell-Penetrating Polypeptide and a Potential Drug Carrier. PLoS ONE 7(11): e48913. doi:10.1371/journal.pone.0048913

Production of human factor VIII-FL in 293T cells using the bicistronic MGMT(P140K)-retroviral vector

Hemophilia A is the most common X-linked bleeding disorder; it is caused by deficiency of coagulation factor VIII (FVIII). Replacement therapy with rFVIII produced from human cell line is a major goal for treating hemophilia patients. We prepared a full-length recombinant FVIII (FVIII-FL), using the pMFG-P140K retroviral vector. The IRES DNA fragment was cloned upstream to the P140K gene, providing a 9.34-kb bicistronic vector. FVIII-FL cDNA was then cloned upstream to IRES, resulting in a 16.6-kb construct. In parallel, an eGFP control vector was generated, resulting in a 10.1-kb construct. The 293T cells were transfected with these constructs, generating the 293T-FVIII-FL/P140K and 293T-eGFP/P140K cell lines. In 293T-FVIII-FL/P140K cells, FVIII and P140K mRNAs levels were 4,410 (+/- 931.7)- and 295,400 (+/- 75,769)-fold higher than in virgin cells. In 293T-eGFP/P140K cells, the eGFP and P140K mRNAs levels were 1,501,000 (+/- 493,700)- and 308,000 (+/- 139,300)-fold higher than in virgin cells. The amount of FVIII-FL was 0.2 IU/mL and 45 ng/mL FVIII cells or 4.4 IU/mu g protein. These data demonstrate the efficacy of the bicistronic retroviral vector expressing FVIII-FL and MGMT(P140K), showing that it could be used for producing the FVIII-FL protein in a human cell line.