Wednesday, September 19, 2012

پست بیست و دوم: ادامه پست یازدهم پرسسینگ

ادامه پست یازدهم 

وقتی استریپهای زشت و مزاحمی در طیف داریم با تبدیل فاکتور 
c
  از یک به نیم درست میشه 

در واقع اگر نگاه کنیم این حرف "سی"* در مقابل ساین** یا كساین*** قرار گرفته و این فاکتور ساین رو به کساین تبدیل میکنه اگر نیم باشه و ساین میمونه با تاخیر ۱ دور یا ساین هست بدون تاخیر!
*c
**syn
***cos
c=0
c=0.5
c=1

Monday, July 23, 2012

پست بیست و یکم: همه کامندهای گرافیکی و محاسباتی pymol

*******گرافیکی*******

برای اینکه هر  ۲٠تا ساختار رو ببینیم کامند زیر اجرا میکنیم:

split_states filename


برای اینکه نحوه نمایش ساختار رو به مدلی که دوست داریم در بیاریم
لینک زیر خیلی کمک  میکنه:
http://www.pymolwiki.org/index.php/Cartoon
در واقع اینها مدلهای مختلف نمایش کارتونی هستند و بعد از اجرای کامنت اگر کارتون رو انتخاب کنم میبینم که مدلی که میخواستم رو نشون میده، من هی کامنت رو میزدم بعد میدیدم اتفاقی نیافتاد در حالیکه اصلا نمایشم کارتونی نبود که چیزی ببینم،بعد که کارتون رو انتخاب کردم درست شد


cartoon type, (selection)
type:
  1. skip
  2. automatic
  3. loop
  4. rectangle
  5. oval
  6. tube
  7. arrow
  8. dumbbell

برای تنظیم سایز نمایش کره ای اتمها:

set sphere_scale, 0.5


set stick_radius, 0.12


برای تنظیم قطر نمایش روبانی پروتئین می شه کامند زیر استفاده میکنیم:

set ribbon_width, 2

برای اینکه زمینه شفاف باشد


set ray_opaque_background, off

و برای ایجاد کیفیت بالا ازکامند  زیر استفاده میکنیم:
ray 2048, 1536




*******محاسباتی*******




برای اینکه دو تاساختار رو با هم الاین کنیم:

align pdbfile ////CA, pdbfile , object=alignment
یک روش دیگه برای الاین کردن دادن این کامنت هست:

 
set sphere_scale, 0.5
 align structure2 and resi 1-100, structure1 and resi 300-400
فرم کوتاه

align structure2 & i. 1-100, structure 1 & i. 300-400
یکی از کامنتهای خوب این هست که ساختارها رو با هم الاین کنیم. قبلا کامتش رو اینجا گذاشته بودم
 ولی به نظرم این وبسایت بهتر توضیح داده:

که توش میگیم که از چه شماره آمینو اسید تا چه شماره ای رو بگیره، و اینجوری 
RMSD
رو هم حساب میکنه.

اگر فقط بخواهیم زنجیره اصلی رو برامون حساب کنه و الاین کنه میتونیم کامنت زیر رو بدیم:


align structure2 & i. 1-100 & n. n+ca+c+o, structure1 & i. 300-400 & n. n+ca+c+o
اگر بخوایم خطوطی که اتمها رو هم به هم متصل میکنه ببینیم میتونیم آبجکت رو به آخر کامند اضافه کنیم


align structure2 & i. 1-100 & n. n+ca+c+o, structure1 & i. 300-400 & n. n+ca+c+o, object=matches


برای نشان دادن پروتونهایی که در فاصله خاصی از آمینو اسید مورد نظر قرار دارند:


select H_824, name h* and not resi 824 within 5 of resi 824 and name ha*

یک کار دیگه نشان دادن همه کانفورمرها هست، که برای این از کامند زیر استفاده میکنیم:


split_states pdbfilename


and save each state:

for i in cmd.get_object_list(): cmd.save("%s.pdb" %i, i)


امروز ۱۴( مارچ ۲۰۱۳) وقتی دنبال این بودم که بتونم پایمول رو از ترمینال اجرا کنم، یک کشف خیلی جالب کردم که تو این لینک کلی قابلیتها برای تعریف کارهای مختلف توسط پایمول وجود داره، مثلا وقتی تایپ کنی سفید، زمینه رو سفید کنه...
معنی خیلیهاش رو نمیدونم ولی همه رو دونه دونه اجرا میکنم بلکه گاهی به دردم بخوره... یا اینجا بهش دوباره نگاه کنم و به مرور یاد بگیرم:

کاربرد ماکرو در پایمول


 
سوپروایزرم گیر داده که از ماکرو برای اجرا کردن عملیات روی پروتیین استفاده کنم: من هم رفتم سراغ این لینک که .چیزی ازش سر در نمیارم!! 

چند تا ماکرو این زیر میگذارم امیدوارم که مفید باشه 
 توی درک ماکروها این جدول میتونه مفید باشه که از این لینک پیداش کردم :
Matching
Property
Selector
Short Form
Selector
Identifier
and Example
symbole.chemical-symbol-list
list of 1- or 2-letter chemical symbols from the periodic table
PyMOL> select polar, symbol o+n
namen.atom-name-list
list of up to 4-letter codes for atoms in proteins or nucleic acids
PyMOL> select carbons, name ca+cb+cg+cd
resnr.residue-name-list
list of 3-letter codes for amino acids
PyMOL> select aas, resn asp+glu+asn+gln 
or list of up to 2-letter codes for nucleic acids
PyMOL> select bases, resn a+g
resii.residue-identifier-list
list of up to 4-digit residue numbers
PyMOL> select mults10, resi 1+10+100+1000
residue-identifier-range
PyMOL> select nterm, resi 1-10
altaltalternate-conformation-identifier-list
list of single letters
PyMOL> select altconf, alt a+""
chainc.chain-identifier-list
list of single letters or sometimes numbers
PyMOL> select firstch, chain a
segis.segment-identifier-list
list of up to 4 letter identifiers
PyMOL> select ligand, segi lig 
flagf.flag-number
a single integer from 0 to 31
PyMOL> select f1, flag 0
  numeric_type  nt.type-number
a single integer
PyMOL> select type1, nt. 5
text_typett.type-string
a list of up to 4 letter codes
PyMOL> select subset, text_type HA+HC
ididexternal-index-number
a single integer
PyMOL> select idno, id 23
indexidx.internal-index-number
a single integer
PyMOL> select intid, index 11
sssssecondary-structure-type
list of single letters
PyMOL> select allstrs, ss h+s+l+""



RMSD: Root Mean Square Deviation is the square root of the mean of the square of the distances between the matched atoms.
RMSD = SQRT[{SUM(dii)2}/N]
where dii is the distance between the ith atom of structure 1 and the ith atom of structure 2 and N is the number of atoms matched in each structure.


Thursday, July 5, 2012

پست بیستم: چگونه یک فایل قابل اجرا را در لینوکس پیدا کنم؟



دو تا کامند که خیلی دیدم
which lpr?
whereis lpr?
اولی هدایت میکنه به اینجا: lpr: /usr/ucb/lpr /usr/man/man1/lpr.1
In this example, the query asked about the lpr command, which spools jobs to printers. The operating system returned two answers, and thus two paths. The first path is the location of the lpr executable, and the second path is the location of the lpr manual page.

/var/bsd/lpr
This means that when you enter lpr at the command line, the system is really executing /var/bsd/lpr.


حالا مشکلم اینه که اینارو که میزنم هیچ جوابی نمیگیرم و نمیدونم این یعنی اینکه من این فایل رو ندارم یا معنی دیگه ای میده؟

Friday, April 27, 2012

پست نوزدهم: پیک ۹۰ درجه خارج از فاز


گاهی اوقات کوتاه بودن زمان تاخیری d1 باعث میشود که به نظر برسد که پیک ۹۰ درجه خارج از فاز است، به اون نوع پیک که حاصل آن دو پیک مثبت و منفی کنار هم هست، دیسپرسیو میگویند و اغلب باید با ۹۰ درجه تغییر فاز درست شود وگرنه شاید مشکل از d1 باشد...
همچنین در بعد غیر مستقیم نیز این امکان وجود دارد که زمان اکوییزیشن کوتاه بوده باشد و بخش مهمی از  FID از دست رفته باشد. به این پدیده ناقص سازی FID گویند. 

Tuesday, January 10, 2012

درس هجدهم: تجزیه تحلیل فایلهای پردازش


دستکاری FID با اهداف زیر صورت میگیرد:
۱. پیش بینی خطی (LP) برای تصحیح نقاط بد ( انتخابی)
۲. شیفت فرکانس FID ( انتخابی)
۳. تفریق حلال
۴. تابع وزن دهی ( تقریبا همیشه)
۵. پیش بینی خطی برای گسترش اطلاعات دومین زمان ( انتخابی)
۶. پرکردن با صفر (همیشه) procXY.com

تفریق حلال -
nmrPipe -fn SOL
برای تفریق حلال یک تابع کانولوشن دومین زمان به کار میبرد
nmrPipe -fn POLY -time
.یک عبارت چند جمله ای را ازاطلاعات دومین زمان یا فرکانس تفریق می کند. در دومین زمان، این میتواند برای تفریق حلال استفاده بشود
مثال:
طیف یک بعدی بدون تفریق حلال

طیف یک بعدی با تفریق حلال

طیف دو بعدی بدون تفریق حلال

طیف دو بعدی با تفریق حلال
SP تابع وزن دهی -
تابع وزن دهی برای به صفررساندن
FID
استفاده می شود، که مانع از ازدیاد سیگنالهای
FID
در پایان دومین زمان، پس از تبدیل فوریر که پیکهای جانبی در جوانب پیک اصلی ایجاد میشود، میشود
معمولا تابع وزنی مورد استفاده در
NMR
برای سیگنال به نویز ماکزیمم (ماکزیمم حساسیت و مینیمم تفکیک) یک تابع سینوسی با ۹۰ درجه شیفت (تابع کسینوسی ) هست. برای افزایش قدرت تفکیک(کاهش حساسیت) باید افست را به مثلا ۴۵ درجه کاهش داد

پر کردن با صفر -
پر کردن با صفر فرایند افزودن صفر به پایان
FID
برای بالابردن ظاهر پیک پس از تبدیل فوریه است

پیش بینی خطی -
پیش بینی خطی برای بالا بردن نسبت سیگنال به نویز و قدرت تفکیک با افزایش تعداد نقاط دیتاانجام میشود. به طور معمول ، نباید بیشتر از تعداد واقعی نقاط دیتا پیش بینی صورت بگیرد.

تبدیل فوریه -
سیگنال را از دومین زمان به دومین فرکانس منتقل میکند و به هر یک از سیگنالها اجازه میدهد که از سینوس کمپلکس بازیافت شوند! واقع که عجب ترجمه ای
کسی که اینجا سر نزده تو این مدت خودم هم که میفهمم منظورم چیه
اگر شما نفهمیدید تو همون لینک اول این مطلب برید اصل مطلب رو به زبون انگلیسی میتونید ببینید..

تنظیم فاز -
تنظیم فاز یک طیف بسیار مهم هست
دو راه برای تنظیم فاز وجود دارد. اول، فاز کل طیف، مستقل از دومین فرکانس قابل تنظیم هست، وقتی وقتی تنظیم فاز درجه صفر انجام میشود، طیف اغلب خوب تنظیم فاز شده است. به هر حال، به خاطر بعضی نواقص در توالی پالس، گاهی لازم است که از پارامتر درجه اول فاز استفاده شود که فاز را بسته به دومین فرکانس تصحیح میکند.




#!/bin/csh

#
# Processing of a 2D Plane from 3D States HN-Detected data:
# "POLY -time" subtracts solvent by polynomial fitting.
# "EXT -left" extracts left half of HN dimension.
nmrPipe -in fid/test001.fid \
| nmrPipe -fn POLY -time \
| nmrPipe -fn SP -off 0.5 -end 0.98 -pow 2 -c 0.5 \
| nmrPipe -fn ZF -size 2048 \
| nmrPipe -fn FT \
| nmrPipe -fn PS -p0 -95 -p1 0.0 -di \
| nmrPipe -fn EXT -x1 5ppm -xn 12ppm -sw -verb \
| nmrPipe -fn TP \
| nmrPipe -fn LP -fb -ord 8 -pred 16 -after -verb \
| nmrPipe -fn SP -off 0.5 -end 0.98 -pow 2 -c 0.5 \
| nmrPipe -fn ZF -size 128 \
| nmrPipe -fn FT -alt \
| nmrPipe -fn REV \
| nmrPipe -fn PS -p0 0.0 -p1 0.0 -di \
| nmrPipe -fn TP \
-verb -ov -out ft/XY001.ft2