Line 115: | Line 115: | ||
<p style="text-align:center; font-size:20px;"><font face="verdana">At chemical equilibrium, the total gibbs free energy of the system will be minimized. The equation for the Gibbs Free Energy of the aptamer-protein interaction is given by: | <p style="text-align:center; font-size:20px;"><font face="verdana">At chemical equilibrium, the total gibbs free energy of the system will be minimized. The equation for the Gibbs Free Energy of the aptamer-protein interaction is given by: | ||
<i>G =(Go+RTln(Q))</i>, where <i>G</i> denotes the gibbs free energy of the system, <i>Go</i> denotes the free energy change of the reaction under standard conditions (standard gibbs free energy of formation), <i>R</i> denotes the ideal gas constant (8.314 <i>Jk-1 mol-1</i>), <i>T</i> denotes the temperature in Kelvin, and <i>Q</i> denotes the reaction quotient, which is given by <i>Q=[product(s)][reactant(s)]</i>, at the given state. Because aptamer interactions are often defined in terms of dissociation constants (<i>kd</i>), rather than by the standard gibbs free energy of formation (Go), we will define Go in terms of <i>kd</i> to calculate the gibbs free energy change. Because <i>Go = -RTln(1/kd)</i> at equilibrium, we can use <i>kd</i> to find <i>Go</i>.</p> | <i>G =(Go+RTln(Q))</i>, where <i>G</i> denotes the gibbs free energy of the system, <i>Go</i> denotes the free energy change of the reaction under standard conditions (standard gibbs free energy of formation), <i>R</i> denotes the ideal gas constant (8.314 <i>Jk-1 mol-1</i>), <i>T</i> denotes the temperature in Kelvin, and <i>Q</i> denotes the reaction quotient, which is given by <i>Q=[product(s)][reactant(s)]</i>, at the given state. Because aptamer interactions are often defined in terms of dissociation constants (<i>kd</i>), rather than by the standard gibbs free energy of formation (Go), we will define Go in terms of <i>kd</i> to calculate the gibbs free energy change. Because <i>Go = -RTln(1/kd)</i> at equilibrium, we can use <i>kd</i> to find <i>Go</i>.</p> | ||
+ | |||
+ | <p style="text-align:center; font-size:20px;"><font face="verdana">Restrictions will be placed on this equation before the minimum gibbs free energy of the system is calculated. Our main restrictions are the mass balance equations. In other words, the concentration of aptamer one fed into the reaction (as reactants) should equal the amount that comes out of the reaction in the form of complexes. The same logic can be applied to the other three reactants (aptamer 2, connector, and protein), to provide a total of four restrictions (mass balance equations). Our final restrictions involve the fact that all concentrations have to be greater than or equal to zero (check “Code for Mathematica” for details).</p> | ||
+ | |||
+ | <p style="text-align:center; font-size:20px;"><font face="verdana">Once our restrictions are in place, we will employ a built-in function on mathematica to find the minimum gibbs free energy of the system. The syntax for the minimization function is given by: Minimize[{f,cons},{x,y,…}]. We will find the minimum gibbs free energy of the system for several different concentrations of the reactants, in order to find the concentration that optimizes the aptamer-protein interaction.</p> | ||
</div> | </div> | ||
Revision as of 21:46, 18 October 2016