Wednesday, February 12, 2014

Translational fidelity / transcriptional fidelity, a collection


http://www.yeastgenome.org/go/GO:1990145/overview


EFB1 essential gene
SUI1 essential gene
EFT1 YOR133W non-essential gene, null has very short lifespan
EFT2 is a paralog of EFT1, non-essential. 

Michael Sovaneau's 1979 PNAS paper on tRNA aminocylation proofread.

High fidelity translation in nude mouse, Nature paper 2013.


Yeast, SUP45 (essential gene)
http://www.yeastgenome.org/reference/S000043151/overview
Overexpression SUP45 leads to decreased vegetative growth

"SUP45 TCH tetO"


http://www.yeastgenome.org/reference/S000151121/overview


 http://www.ncbi.nlm.nih.gov/pubmed/24082110

Proc Natl Acad Sci U S A. 2013 Sep 30. [Epub ahead of print]

Naked mole-rat has increased translational fidelity compared with the mouse, as well as a unique 28S ribosomal RNA cleavage.

Source

Department of Biology, University of Rochester, Rochester, NY 14627.

Abstract

The naked mole-rat (Heterocephalus glaber) is a subterranean eusocial rodent with a markedly long lifespan and resistance to tumorigenesis. Multiple data implicate modulation of protein translation in longevity. Here we report that 28S ribosomal RNA (rRNA) of the naked mole-rat is processed into two smaller fragments of unequal size. The two breakpoints are located in the 28S rRNA divergent region 6 and excise a fragment of 263 nt. The excised fragment is unique to the naked mole-rat rRNA and does not show homology to other genomic regions. Because this hidden break site could alter ribosome structure, we investigated whether translation rate and amino acid incorporation fidelity were altered. We report that naked mole-rat fibroblasts have significantly increased translational fidelity despite having comparable translation rates with mouse fibroblasts. Although we cannot directly test whether the unique 28S rRNA structure contributes to the increased fidelity of translation, we speculate that it may change the folding or dynamics of the large ribosomal subunit, altering the rate of GTP hydrolysis and/or interaction of the large subunit with tRNA during accommodation, thus affecting the fidelity of protein synthesis. In summary, our results show that naked mole-rat cells produce fewer aberrant proteins, supporting the hypothesis that the more stable proteome of the naked mole-rat contributes to its longevity.


H3K36 methylation promotes longevity by enhancing transcriptional fidelity.
Epigenetic mechanisms, including histone post-translational modifications, control longevity in diverse organisms. Relatedly, loss of proper transcriptional regulation on a global scale is an emerging phenomenon of shortened life span, but the specific mechanisms linking these observations remain to be uncovered. Here, we describe a life span screen in Saccharomyces cerevisiae that is designed to identify amino acid residues of histones that regulate yeast replicative aging. Our results reveal that lack of sustained histone H3K36 methylation is commensurate with increased cryptic transcription in a subset of genes in old cells and with shorter life span. In contrast, deletion of the K36me2/3 demethylase Rph1 increases H3K36me3 within these genes, suppresses cryptic transcript initiation, and extends life span. We show that this aging phenomenon is conserved, as cryptic transcription also increases in old worms. We propose that epigenetic misregulation in aging cells leads to loss of transcriptional precision that is detrimental to life span, and, importantly, this acceleration in aging can be reversed by restoring transcriptional fidelity.

Educational problems


Student cheating problems

http://www.nytimes.com/2012/09/08/education/studies-show-more-students-cheat-even-high-achievers.html?_r=0

http://www.nytimes.com/2013/08/04/education/edlife/educators-study-online-dishonesty.html

http://www.nytimes.com/2013/09/17/education/students-accused-of-cheating-return-awkwardly-to-a-changed-harvard.html

http://education-portal.com/articles/75_to_98_Percent_of_College_Students_Have_Cheated.html

http://abcnews.go.com/Primetime/story?id=132376

Honor code
http://www.nytimes.com/2014/04/13/education/edlife/the-fading-honor-code.html?WT.mc_id=D-NYT-MKTG-MOD-70311-06-10-HD&WT.mc_ev=click&WT.mc_c=${CAMPAIGN_ID}

ms02 R, batch run

# ms02-2014Feb12.R
#permuate merged yeast PPI+GIN

#2014 Feb 12, re-name function to ms02_singlerun
#2014 Jan 31, fixed a bug that inserted "NA" into new network. The bug seems to be caused by spliting the
# arrays. I rewrote the spliting portion.



#require(igraph)
rm(list=ls())
debug = 0
setwd("~/projects/0.ginppi.reliability.simulation/ms02GINPPI")
#set.seed(2014)

#permute.pairs.wo.selfpairs = function( inpairs,  ncycles=10, debug=1 ) { 
ms02_singlerun = function( inpairs,  ncycles=10, indebug=0 ) { # Renamed, 2014 Feb 12
  if (ncycles >= 1 ) {
    if(indebug>0) {
      print(paste('ncycles=', ncycles))
    }
    longids = c(as.character(inpairs[,1]), as.character(inpairs[,2]) )
    longids = sample(longids)
    len = length(inpairs[,1])
    newpairs = data.frame( cbind( longids[1:len], longids[(len+1): (2*len)]) )
    names(newpairs) = c('id1', 'id2')
    newpairs$id1 = as.character( newpairs$id1)
    newpairs$id2 = as.character( newpairs$id2)   
    newpairs$selfpairs = ifelse( newpairs$id1 == newpairs$id2, 1, 0 )
    self.tb = newpairs[ newpairs$selfpairs==1, ]
    nonself.tb = newpairs[newpairs$selfpairs==0, ]
    if(indebug>0) {
      print(paste("===selfpairs===="),NULL)
      print(self.tb)
      print(paste("================="),NULL)
    }
    if( length(self.tb[,1])>=1 ) {
      if ( ncycles == 0) {
        #return (c(NA,NA, NA) );
        print(paste("ncycles reached zero, ncycles"),ncycles)
        print(paste("Abort!"),NULL)
        stop;
      } else {
        ncycles = ncycles - 1
        splitPos = round( length(self.tb[,1]) * sqrt(ncycles) ) + 5  #2014Jan31 change
        splitPos = min( splitPos, (length(nonself.tb[,1])-1 ) )
        selectedpairs = rbind(self.tb,  nonself.tb[1: splitPos, ] )
        restpairs = nonself.tb[ (splitPos + 1): length(nonself.tb[,1]), ]
        #return( rbind(restpairs, permute.pairs.wo.selfpairs(selectedpairs, ncycles)))
        return( rbind(restpairs, ms02_singlerun(selectedpairs, ncycles)))  #2014 Feb 12
      }
    } else { 
      return (newpairs)
    }
  } else {
    return( c(NA,NA,NA ))
  }
}

#net = read.table("repeat.tab")
#write.table(pairs, "merged_PPIGIN_2014Jan20.tab", quote=F, row.names=F, col.names=F, sep='\t')
net = read.table( "merged_PPIGIN_2014Jan20.tab", header=F, sep="\t", colClass = c("character", "character") )
head(net)
if(debug==9) {
  #net = read.table('pair.tab',header=F)
 net = net[1:90000,]
}

for( i in 1:100) {
 net.ms02 = ms02_singlerun( net, indebug=0 )
 cmnd = paste( "mkdir dipgin.ms02.output/", i, sep="")
 system( cmnd )
 outputname = paste( 'dipgin.ms02.output/', i, '/', "ms02_",i,".tab", sep="")
 write.csv(net.ms02, outputname)
}


#do they have the same degree?
#t1 = table(c(net[,1],net[,2]))
#t2 = table(c(net.ms02[,1],net.ms02[,2]))
#comp <- t1 == t2
#table(comp)
#tf = comp[comp==F]; tf
#t1[names(tf)[1]]
#t1[names(tf)]
#t2[names(tf)]




DMG problem in Snow Leopard

RStudio-0.98.501.dmg file is not automatically mounted on a snow leopard laptop. Open it with "Disk Utility", then click "Mount". It worked.


VPN Spelman



http://princess.spelman.edu/accountreq.nsf/vpninstructions

sudo sh *.sh

install 'Cisco' folder to Application, which contains Cisco Anyconnect client.


Tuesday, February 11, 2014

Youtube screen ratio is 16:9

The 16:9 ratio for youtube means 1280x720 screen.

Bio125, Feb 11, Tue, 2014, ApE, identify mutations, gel electrophoresis

I manually generated mutations and give students for analysis.  I made two mistakes and were found out by studetns.

I posted a youtube tutorial and asked student finish the ApE assignment in class. This seems to push student to focus on the work. At 1:30pm, the college was closed and some student left. Others stay and kept working on the ApE exercieses.

Gel were poured by lab instructor. Student loaded their samples, and gels were run. I found that ApaI and Xab I gave two similar-sized bands for pMSH2 and one sized-band for pRS413. When the plasmid are overloaded, the one 4.5kb smear band is similar to the 4.4 and 4.9 kb bands. We decided to run the gel longer to get better resolution.  (Next time we should pick different RE, or we can ask student to pick better RE. 0.8% gel might be used). 



Because of the HA tag, locations can be different between MSH2 and  HA-MSH2, (info based on D. Wang). 
1.     G704S: AG41_G693S

2.     R542P: AG29_R542P

3.     H669Y: AG33_H658Y

4.     G704D: AG91_G693D

5.     M718I: AG371_M707I

6.     H669R: AG372_H658R

7.     A618V: AG421_A618V

8.     D621N: AG422_D621N

Monday, February 10, 2014

GAS 2014 annual meeting prep, math ramp research


Friday and Saturday, 28-29 March, 2014 at Georgia Regents University, Augusta, GA. 

meeting registration site
http://www.gru.edu/trends/gas/accommodations.php
http://www.gaacademy.org/annualmeeting10.htm

hotel:
Hilton Garden Inn Augusta Meeting rate ($89) but not available anymore. 
1065 Stevens Creek Road, Augusta, Georgia, 30907

Sheraton Augusta Hotel  , prepaid rate $107 per night, refundable rate $119
1069 Stevens Creek Road, Augusta, GA 30907
Sheraton Augusta Hotel
1069 Stevens Creek Road
Augusta, Georgia 30907 United States
Phone: (706) 396-1000 Fax: (706) 396-1001



Parking?

-------------------------
Prepare of a 10 minutes talk.

Research talk:

10 slides.

  1. background, gene network, essential and non-essential,
  2. background, what is aging, our quantitative definition
  3. specific aim: how reliable is the yeast gene interaction network.
  4. method: how to generate random networks.
  5. method: evaluate the failure rates of original and random networks
  6. method: how network failure is simulated and evaluated
  7. result: histogram of failure rate, average lifespan of random networks, with original network
  8. result: Gompertz versus Weibull (biological versus machine aging)
  9. summary and conclusions
  10. future directions:

acknowledge, funding support

things to do:
->flow chart to generate 100 random networks and simulate their failures.
-> run the two R codes on network aging simulation and network permutation.




BIO125, FOA plate media and western blot media

=> For FOA plates, the media is -HIS-TRP-LEU-THR+URA 

-LEU -THR for URA3's LEU2 promoter to be active.  The presence of Leucine and Theorine lead to 10 fold drop of LEU2 expression.

-HIS to keep pMSH2

-TRP to keep pSH44 reporter plasmid

+URA for FOA assay, because we are selecting for mutant URA3.
So the FOA media is HIS-TRP-LEU-THR + URA


However, before the FOA plate, yeast cells should be grown in - URA to ensure intact dinucleotide repeat in pSH44.
So the regular media is -URA-HIS-TRP (-LEU-THR)

See Gammie04CBE


=> For MSH2 expression:

Gammie07,

Three type of pMSH2 was mentioned. pMSH2, pGBD-MSH2, and pGAL-MSH2.

The pMSH2 is under endogenous promoter. So, the media should be -His to keep pMSH2.

The pGBD-MSH2 was designed to have a constitutive expression of MSH2 protein with HA tag. This plasmid seems to be used for two-hybrid.


 The pGAL-MSH2 contains a GAL10 promoter.


References
http://www.ncbi.nlm.nih.gov/pmc/articles/PMC294742/pdf/jbacter00572-0372.pdf


self-pairing test, ms02, R, network permutation

test network: 



Expectation: When ncycles run to 1, the last row should be "NA"

Summary: the code 'ms02' passed the test. No extra row or self-pairs are found.

2014April 7. 'ms02' did not check id1-id2 versus id2-id1. 



exercises on rate, R code

MM = function(x, Km, Vmax) {
 ret = Vmax /(1 + (Km/x)^10);
}
x = seq(-0.1,90, by=0.1);
Km = 40;
Vmax = 45; 
y = MM( x,Km,Vmax);
plot( y ~ x,  type='l', lwd=4, ylab='concentration', xlab='time'); #sigmoid shape transport

#get the rate

dy = NA;
for ( i in 2:length(y)) {
  dy[i]= y[i] - y[i-1]
}
plot( dy ~ x, type='l', lwd=4, ylab='rate', xlab='time')

dy2 = 0.25 - dy

plot( dy2 ~ x, type='l', lwd=4, ylab='rate', xlab='time')

KM2 = 40

Vmax2 = 0.25
y2 = MM( x,KM2,Vmax2);
plot( y2 ~ x,  type='l', lwd=4, ylab='rate', xlab='time'); #sigmoid shape transport

y2b = 0.25 - y2

plot( y2b ~ x,  type='l', lwd=4, ylab='rate', xlab='time'); #sigmoid shape transport

y3 = c(rep(0,40),rep(0.25,40))

x3 = seq(1,80,by=1)
plot( y3 ~ x3, type='l', lwd=4, ylab='rate',xlab='time')

y3b = 0.25 - y3

plot( y3b ~ x3, type='l', lwd=4, ylab='rate',xlab='time')

Thursday, February 6, 2014

online alignment tools

t-coffee
clustalw
muschle

Bio125, Xba1 Apa1, pRS413, pMSH2

1hour, went through protocol for XbaI and Apa1 cutting.

Use ApE to do XbaI and ApaI mapping. Dr. Kioko ordered FastDigest enzyme and buffer. Both enzyme work at 37C. The 37C incubation did not finish because the heat block somehow could not reach 37C.

Explained why 6bp cutters have fewer sites than 4bp cutters. Most students use plus instead of multiplication to do the calculation. I used 3 colored play dough again to enumerate the 9 possibility for 2 positions. I initially tried 2 colored play dough, but ran into the problem of 2+2 and 2x2 yield the same result.

Kioko mixed XbaI and ApaI. Student took 1ul of the mixture.




Bio233, Feb 6 2014

Played music science video on glycolysis and Kreb cycle.

Emphasized electron transport chain, went through the calculation of 38 ATP.

Emphasized negative feedbacks.

Summarized key learning objectives at the end.