150 likes | 301 Views
BL23A BioSAXS check list. fitted solution envelope of native cyt . c.
E N D
BL23A BioSAXS check list fitted solution envelope of native cyt. c Acknowledge EMBL Hamburg BioSAXS team memebers’ effort in ATSAS software package. The data analysis software, ATSAS software package, are downloaded from http://www.embl-hamburg.de/biosaxs/atsas-online/download.php If one use some programs in the work, please cite the corresponding papers.
Bio-SAXS Check list 1. Check the radiation damage 2. Check aggregation by AutoRg program. Calculate theoretical I0, MW, and Vtotwith protein SAXS I0 estimtaion.xls. 3. Check inter-particle interaction with concentration dependence I0/c 4. Profile merging merge the low and high concentration curves. 5. Crysol Program Compare (fit) solution structure with PDB crystal structure. 6. Transform data to P(r) with GNOM program and data fitting. 7. Dammin/ GasborProgram An ensemble of dummy atom/residues model simulation with P(r) output by GNOM.
1. Check the radiation damage • Instruction: • Plot each frame and combined data. • Check the low-q region in the all profiles. • If all profiles are overlapped well at low-q region, no radiation damage occurs. • If some frames of them appear up-turn at low-q region, radiation damage may occur. • The radiation damage might result from long data requisition time or aggregation sample. • Resolving ❶1. combine the frames without up-turn and out a low-q reliable curve. • 2. combine all frames and output high-q reliable curve. • 3. fine-tune and merge the two curves and output into one profile. • Resolving ❷ 1. Run HPLC 1 2 well-overlapped
2. Check aggregation by AutoRg program. Instruction: Start program: Primus -> AutoRg; input file type: txt, dat. InPlot tab, see the raw data. In Guinier tab, check the aggregation and define the qmin before unreliable range. InInfo tab, read the sRg limits (Guinier region), I0, and Rgvalues. 2 1 4 3 Define qmin before unreliable range.
One can calculate theoretical MW, I0, and Vtotwith protein SAXS I0 estimtaion.xls. (spreadsheet edited by Dr. Jeng) Instruction: Fill in protein amino acids numbers in AQ column. Fill in protein concentration in B18 cell. Fill in the aggregation number in B(13) cell. Check the molecular weight calculated in B(17). Check the calculated I(0) at B(19). Input measured I(0) in B(27) to check the aggregation no. of oligomers. e.g. Cytochrome c (PDB ID: 1HRC), MW~12368 3 1 4 2 5 6
3. Check inter-particle interaction with concentration dependence I0/c. • Instruction: • Plot all different concentration curves. • Confirm if concentration dependence is distinguishable at low-q region. • e.g. There is concentration dependence behavior in below figure. The high conc. profile displays lower I0/c value due to the interference from inter-particle interaction. inter-particle interaction Attention! Need to measure low concentration sample as low as possible. Make sure that more than two I0/c curves are overlapped well.
4. Profile merging – merge the low and high concentration curves. Instruction: Start program: Primus Click “Tools” to select the low and high concentration data for merging. Click “plot” to plot the two curves. 2 * * *nBeg value First 45 points in #1 and #2 are cut because of beamstop. 3
According to the q-cut range (see ① and ②) , input the parameters into nBeg (begin point #) and nEnd (end point #) of Data Processing. (5) Click “Plot Scale ” and make sure that the curves are superimposable. ① remove concentration dependence ②remove divergence region 4 • q1-cut range • (high conc. Profile) q2-cut range (low conc. profile) Input q-cut region based on ① and ② The high-q endpoints of low conc. profile do not meet the other curve. 5 The scale multiplier is based on the intensity of the first file (#1)
(6) Fine tune the eEnd point to make sure that the endpoints are superimposed with the merged curve. 6 Each Merge-click can output one merged file. Filename of 1st profile-merging: Merge00.dat Filename of 2nd profile-merging: Merge01.dat • Attention! • Need to measure low concentration sample as low as possible. • Need to check concentration dependence before model simulation. • Make sure that all different concentration curves are superimposable after q-cut.
5. Compare (fit) solution structure with PDB crystal structure by Crysol Program (in ATSAS package). • Instruction: • Enter your option <0>: • Select the PDB file (protein data bank) • Maximum order of harmonics <15>: • Order of Fibonacci grid <17>: • Maximum s value <0.5>: • Number of points <51>: • Account for explicit hydrogens <no>: • Fit the experimental curve <yes>: • Enter data file (experimental dat file) • Subtract constant <no>: • Angular units in the input file <1>: • Electron density of the solvent <0.334>: • Plot the fit <yes>: • Another set of parameters <no>: • Press “” to terminate the program e.g. Cytochrome c (1HRC), MW~12368 2 9 Experimental SAXS data
(16) Plot the 1st (experimental scattering vector) and 2nd (theoretical intensity in solution) columns of fit file. (17) Open the log file to read the experimental and theoretical Rg values. fitted solution envelope of native cyt. c Refer to http://www.embl-hamburg.de/biosaxs/manual_crysol.html#un
6. Transform data to probability distribution function P(r) with GNOM program and data fitting • Instruction: • Printer type [ postscr ] : • Input data, first file : (select the file) • Output file [ gnom.out ] : filename.out • No of start points to skip [ 0 ] : • Input data, second file [ none ] : • No of end points to omit [ 0 ] : • Angular scale (1/2/3/4) [ 1 ] : • Plot input data (Y/N) [ Yes ] : • File containing expert parameters [ none ] : • Kernel already calculated (Y/N) [ No ] : • Type of system (0/1/2/3/4/5/6) [ 0 ] : • Zero condition at r=rmin(Y/N) [ Yes ] : • Zero condition at r=rmax (Y/N) [ Yes ] : • Rmaxfor evaluating p(r) : (input a proper value) • Number of points in real space [ 101 ] : • Kernel-storage file name [ kern.bin ] : • Experimental setup (0/1/2) [ 0 ] : • Initial ALPHA [ 0.0 ] : • Plot results (Y/N) [ Yes ] : • Press CR to continue probability distribution function P(r) In this case, Rmax (diameter of the object) = 40 Å
7. Dammin Program Dummy atom model simulation with P(r) output by GNOM program. • Mode: <[F]>ast, [S]low, [J]ag, [K]eep, [E]xpert < Fast >: • Log file name < .log >: filename.log • (PLEASE SELECT THE INPUT FILE NAME) • (3) Input data, GNOM output file name < .out >: filename • File to be opened: m.out • Project identificator: m • (4) Enter project description : • Random sequence initialized from : 172235 • ** Information read from the GNOM file ** • Data set title: Merge of: cytc-5.txt cytc-10.txt • Raw data file name: Merge01.dat • Maximum diameter of the particle : 40.00 • Solution at Alpha = 0.203E+01 • Rg : 0.132E+02 I(0) : 0.388E-01 • Radius of gyration read : 13.20 • Number of GNOM data points : 428 • (5) Angular units in the input file: • 4*pi*sin(theta)/lambda [1/angstrom] (1) • 4*pi*sin(theta)/lambda [1/nm ] (2) < 1 >: • Maximum s value [1/angstrom] : 0.3924 • Number of Shannon channels . : 4.996 • Portion of the curve to be fitted < 1.000 >: • Number of knots in the curve to fit : 20 • *** Warning: constant reduced to avoid oversubtraction • A constant was subtracted : 1.788e-4 • Maximum order of harmonics : 10 • (6) Initial DAM: type S for sphere [default], E for ellipsoid, C for cylinder, P for parallelepiped or start file name <.pdb>: (select one type) • (7) Symmetry: P1...19 or Pn2 (n=1,..,12) or P23 or P432 or PICO < P1 >: • (8) Expected particle shape: <P>rolate, <O>blate, or <U>nknown < unknown >: • Start running simulation Blue-colored text come from out file read by Dammin χ2=1.179 Dammin Output Files: log file: contain the same information as the screen output. -0.pdb file: contain the beads of the solvent inside the search volume -1.pdb: represent the modeled particle. fir file: Fit of the simulated scattering curve versus the experimental data. Columns in the output file are: 's', 'Iexp', 'Errexp'and 'Isim'.
7. GasborProgram An ensemble of dummy residuesmodel simulation with P(r) output by GNOM program. Computation mode (User or expert)……<User>: Log file name <.log>:filename .log ***Please Select The Input Filename*** File to be opened: m.out (3) Enter Project decsription: (4) Angular units in the input file: 4*pi*sin(theta)/lambda [1/angstrom] (1) 4*pi*sin(theta)/lambda [1/nm ] (2) < 1 >: (5) Portion of the curve to be fitted……<1.000>: (6) Initial DRM (CR for random)……< .pdb>: (7) Symmetry: P1…19 or Pn2 (n=1,..,12) or P23 or P432 or PIC0….< P1>: (8) number of residues in asymmetric part..< 62>: (9) Fibonacci grid order……<9>: (10) Expected particle shape: <P>rolate, <O>blate, or <U>nknown……<Unknown>: View along x Cytochrome c (1HRC) • Start running simulation View along y Displayed by VMD χ2=2.286 View along z GASBOR Output Files: log file: contain the same information as the screen output. fir file: fit to the raw experimental data pdb file: resulting model in PDB-like format that can be viewed with MASSHA (in ATSAS package) or VMD (molecular visualization program)
Appendix. On-line fit solution structure (experimental profile) with crystal structure (PDB or upload file). • Instruction: • Open the FoXS website http://modbase.compbio.ucsf.edu/foxs/index.html • Input the PDB ID or upload the PDB file. • Upload the experimental profile. • Click “Submit Form”. FoXS: Fast SAXS Profile Computation with Debye Formula