20.109(F21):M2D6

From Course Wiki
Revision as of 18:21, 2 September 2021 by Noreen Lyell (Talk | contribs)

Jump to: navigation, search
20.109(F21): Laboratory Fundamentals of Biological Engineering
Drawing provided by Marissa A., 20.109 student in Sp21 term.  Schematic generated using BioRender.

Fall 2021 schedule        FYI        Assignments        Homework        Class data        Communication        Accessibility

       Module 1: Genomic instability                          Module 2: Drug discovery       


Introduction

Protocols

Part 3: Examine binding shifts

You will receive two Excel sheets containing raw data from each well of a 384 well plate over the specified range of temperatures. The Excel sheet with "Melt Curve RFU Results" in the file name will contain raw fluorescence intensity data, while the other sheet with "Melt Curve Derivative Results" in its name will have the values for the first derivative of the melt curve. The teaching faculty will inform you which wells correspond to which conditions for your group.

One basic way to determine the "melting temperature," or Tm of the protein is to determine temperature at the inflection point of the melting curve. This inflection point would occur at the maximum value of the first derivative. The BioRad CFX machine we use actually exports the negative of the first derivative in the Excel file, so we will find the minimum value in the first derivative Excel file, and take the corresponding temperature to be the Tm of FKBP12 in each condition.

  1. Open the Excel file corresponding to the first derivative data
  2. Column B should contain temperature information in Celsius.
  3. At a row on the bottom of column C, type in the following command: =INDEX($B$FirstRow:$B$LastRow, MATCH(MIN(CFirstRow:CLastRow),CFirstRow:CLastRow,0)), where FirstRow corresponds to the row number of the first row containing data, and LastRow contains the row number of the last row containing data.
  4. Press enter, and double check that the listed temperature occurs at the minimum value of the first derivative.
  5. Then, drag the bottom left corner of the cell across all relevant columns to apply the formula to those columns of interest.
  6. Plot the columns relevant to your data set by making a scatter plot ("straight marked scatter"), having the temperature (values in column B) on the x-axis, and the first derivative values on the y-axis.
  7. Double check by eye that the values you calculated to be the melting temperatures correspond to the minimum values on the curves. (See example plot in the introduction section of this wiki page)
  8. Next, you may also check to see what the melting curves look like in terms of raw fluorescence by plotting fluorescence intensity vs. temperature in the "Melt Curve RFU Results" file. Again, validate the results you found by eye to see if the Tms correspond to the inflection point of the raw fluorescence melt curves.
  9. Check to see if the Tm of the control protein shifted when its ligand was added. Quantify the shift.
  10. Check to see if the Tm of FKBP12 shifted when Rapamycin or other compounds were added. Quantify the shifts.
  11. By varying the concentration of Rapamycin, you will be able to determine an apparent dissociation constant of Rapamycin and FKBP12. Here is the reference for finding the apparent dissociation constant that was mentioned in lecture: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4692391/
    • You may use this MATLAB script (ApparentKd.m) to help fit the Tm vs. Rapamycin concentration curve to a single binding site model and find a value for the apparent KD. The function uses a nonlinear regression.
      • Based on the given article linked above, the single binding site model is as follows, where the fit parameters include Tm_min (minimum Tm at no ligand concentration), Tm_max (maximum Tm at infinite ligand concentration), and KD is the apparent KD value. For our experiment, the concentration of FKBP12 (written as [FKBP12]) was 8.5 uM and the concentration of Rapamycin (written as [Rap]) was variable.
      • SingleSiteBindingModelEquation.png
    • Create an array of Rapamycin concentrations by typing RapConc = [A, B, C, D, E, etc.] at the MATLAB command prompt, where A, B, C are the various Rapamycin concentrations in units of uM.
    • The 10 concentrations of Rapamycin are: 20uM, 10uM, 5uM, 1uM, 0.1uM, 0.05uM, 0.01uM, 0.005uM, 0.001uM, and 0.0001uM
    • Create an array of Tms by typing Tm = [A2, B2, C2, D2, E2, et c.] at the MATLAB command prompt, where A2, B2, C2 are Tms corresponding to concentrations A, B, and C in the previous array.
    • Making sure the MATLAB function ApparentKd.m is in your current working directory, type in ApparentKd(RapConc, Tm) at the command prompt and press enter
    • The function performs a nonlinear regression of your data with a single binding site model, and will return an apparent KD value in units of uM from the best fit.
    • If this regression does not work well, you may use alternative methods to estimate an apparent KD value
      • You can find the EC50 value by fitting with the following formula in this matlab function EC50.m
      • Making sure the MATLAB function EC50.m is in your current working directory, type in EC50(RapConc, Tm) at the command prompt and press enter
      • This file performs nonlinear regression using the following equation, where the fit parameters include Tm_min (minimum Tm at no ligand concentration, Tm_max (maximum Tm at infinite ligand concentration, and EC50 is the EC50 or apparent KD value, and a Hill coefficient which should not be meaningful in this context.
      • EC50equation.png

Reagents

Introduction

Protocols

Reagent list

Navigation links

Next day: Complete data analysis

Previous day: Prepare for secondary assay to test putative small molecule binders