Skip to content

Instantly share code, notes, and snippets.

View JonasR's full-sized avatar
⛰️

JonasR

⛰️
  • AbbVie
  • Germany, Rosenheim area (remote)
  • 10:19 (UTC +02:00)
View GitHub Profile
#!/usr/bin/env python
from __future__ import division
import pymongo as pym
import datetime as dt
import codecs
import re
import sys
import matplotlib.pyplot as plt
import numpy as np
@JonasR
JonasR / gist:2919934
Created June 12, 2012 20:29
helloworld
#!/bim/bash
echo "This even allows versioning and forking!"
@JonasR
JonasR / gist:2944226
Created June 17, 2012 11:16
P06865 sequence
>P06865
MTSSRLWFSLLLAAAFAGRATALWPWPQNFQTSDQRYVLYPNNFQFQYDVSSAAQPGCSV
LDEAFQRYRDLLFGSGSWPRPYLTGKRHTLEKNVLVVSVVTPGCNQLPTLESVENYTLTI
NDDQCLLLSETVWGALRGLETFSQLVWKSAEGTFFINKTEIEDFPRFPHRGLLLDTSRHY
LPLSSILDTLDVMAYNKLNVFHWHLVDDPSFPYESFTFPELMRKGSYNPVTHIYTAQDVK
EVIEYARLRGIRVLAEFDTPGHTLSWGPGIPGLLTPCYSGSEPSGTFGPVNPSLNNTYEF
MSTFFLEVSSVFPDFYLHLGGDEVDFTCWKSNPEIQDFMRKKGFGEDFKQLESFYIQTLL
DIVSSYGKGYVVWQEVFDNKVKIQPDTIIQVWREDIPVNYMKELELVTKAGFRALLSAPW
YLNRISYGPDWKDFYIVEPLAFEGTPEQKALVIGGEACMWGEYVDNTNLVPRLWPRAGAV
AERLWSNKLTSDLTFAYERLSHFRCELLRRGVQAQPLNVGFCEQEFEQT
@JonasR
JonasR / gist:2944236
Created June 17, 2012 11:20
SIFT mutation input file
M1V
L39R
C58Y
L127R
R170W
R178H
S210F
D258H
L451V
E482K
@JonasR
JonasR / gist:2944224
Created June 17, 2012 11:15
PolyPhen 2 Batch file
P06865 1 M V
P06865 39 L R
P06865 58 C Y
P06865 127 L R
P06865 170 R W
P06865 178 R H
P06865 210 S F
P06865 258 D H
P06865 451 L V
P06865 482 E K
@JonasR
JonasR / prepareFiles.pml
Created June 21, 2012 16:02
Prepare 2gjx to be used in structure based mutation analysis
###Create base file
fetch 2gjx
select chainA, chain A
select rest, not chainA
remove rest
select firstDomain, resi 1-166
select secondDomain, not firstDomain
remove firstDomain
zoom secondDomain
@JonasR
JonasR / allmuts_*.sh
Created June 21, 2012 16:18
Preparing sequences for SCWRL and then calling SCWRL on them
#!/bin/bash
for i in R178H R178C P183L D207E S293I F434L L451V E482K L484Q E506D; do
perl prepForSQRWL2.pl -i 2gjxA_sequence -mut $i > SCWRL_$i.sequence
done
for i in R178H R178C P183L D207E S293I F434L L451V E482K L484Q E506D; do
/opt/SS12-Practical/scwrl4/Scwrl4 -s SCWRL_$i.sequence -i ../2gjxA_secondD.pdb -o SCWRL_$i.pdb &>> SCWRL.log
done
@JonasR
JonasR / prepForSQRWL2.pl
Created June 21, 2012 16:19
Lowercase sequence and create one sequence per mutation
#!/usr/bin/perl
use Getopt::Long;
our $opt_f = ();
sub usage {
printf "Usage: $0 [flags] file...\n";
printf "flags:\n";
@JonasR
JonasR / run.txt
Created June 22, 2012 00:55
FoldX runfile
<TITLE>FOLDX_runscript;
<JOBSTART>#;
<PDBS>#;
<BATCH>list.txt;
<COMMANDS>FOLDX_commandfile;
<BuildModel>#,individual_list.txt;
<END>#;
<OPTIONS>FOLDX_optionfile;
<Temperature>298;
<R>#;