UTF-8Modeling hydrophobic effectBGSU CHEM 4450/5450 General Biochemistry I Final Project
Over the past two weeks you've explored properties that affect miscibility of two liquids and considered the difference between properties that affect individual particles vs. those that arise from collections of particles. The Final Project of this semester for our evening sessions is focused on investigating what gives rise to the tendency for soluble proteins to fold with their hydrophobic amino acids on the interior and 'hidden' from the aqueous solvent. You can work towards a better understanding of this phenomenon by modifying the following model to recapitulate a similar behavior ('understanding' in this case means knowing what is required to cause the behavior).
Below, there are two compartments with the following characteristics:
- each has a polymer of 24 linked, green particles (all of the same particle type, 'Pl'). Note: any change you make to the particles of one polymer will by definition affect the other polymer too!
- each polymer is 'pinned' to the center of the compartments (small white 'x' symbol)
- the radius of gyration (Rg) for each polymer is plotted on the graph at the right as a measure of how compact or folded each polymer is. Rg is a calculated value that essentially measures the average distance between all monomers in a polymer. Extended polymers have larger Rg values, whereas folded polymers have smaller Rg values.
The Challenge: Make changes to the compartments (adding solvent particles, changing particle properties, etc.) to control how folded the polymers are. Your goal will be to make the polymer in the left compartment as compact/folded as possible and the other as extended as possible. In other words, create conditions in the two compartments that maximize the difference between their radii of gyration.
Rules:
- the properties and composition of the polymers must remain identical (i.e. any change to one polymer must be made for the other as well).
- any other particles you choose to add must behave like liquids
Note:
- to select an individual particle in the polymer, hold down the <Shift> key and left-click with your mouse (double-left-click to modify the polymer particle properties).
- you are asked to output the radius of gyration data to Excel for questions below. To do this, click the "List data sets" button (leftmost button) above the graph and double click on one of the channels. Channel 0 corresponds to the left polymer, channel 1 to the right (see pink labels on polymers). Paste into MS Excel. There are three columns:
- column 1: 'i' (the index of the data point measured)
- column 2: the time point (i.e. measured in thousands of femtoseconds)
- column 3: the radius of gyration for the given polymer at that index/time point
Answer the following questions in the document below. Send me a screenshot of your responses by e-mail. Then collaborate with your group to write up a short report as instructed in the Google Document linked here: http://tinyurl.com/obvo2xo
JENNIFER GEBES
1.) The goal of this exercise was to design solvent conditions that would cause two identical polymers to behave differently - one folded and the other as extended as possible. Describe in complete sentences what you did to solve the challenge (one folded and one extended polymer, ‘liquid’ solvent, hydrophobic polymers)? Be sure to describe whether your strategy worked and how what evidence you used to come to that conclusion. For both polymers, I labeled it with several small negative charges. To get it to fold, I changed the interactions between the solution particles to have a smaller electrostatic potential (0.05) and the interactions between te different molecules to just about double that. To unfold, I did almost the opposite, making the solution have a larger potential (0.2) and the interaction with the solution and polymer to about half that.
2.) Let your model run for enough time for the two polymers to achieve their different behaviors. Output the Rg vs. time data and plot it in Excel. Paste the overlaid plots of Rg vs. time to support your argument. See attached excel spreadsheet
3.) According to your model, what is the most important factor controlling how folded/extended the polymers are? It really seemed to be the interactions between the molecules and the solution, as well as the interaction of the solution to itself.
4.) A similar phenomenon occurs in nature where proteins fold with hydrophobic surfaces sequestered on their interior away from aqueous solvent. Does your polymer exhibit similar behavior (i.e. is your polymer ‘hydrophobic’ and does your solvent have properties that resemble aqueous solvent? For the most part my folded molecule shows hydrophobic properties, as some of the inner monomers are completely sequestered, not many, but a few. And my solution for the folded protein is very liquid in character, but the solution on the unfolded was a little less, ok, but the interactions were high enough they had a nice mobility, but clung more to the polymer.
10 0 25524truetrue18truetrue161616160 153 0true1616161616255 0 0true161616255 0 0true16255 0 0true16161616true161616161616161org.concord.mw2d.activity.AtomContainerJ. Gebes BGSU_biochem-FinalProject_2013-1209$0.mmlorg.concord.modeler.PageXYGraphorg.concord.mw2d.models.MolecularModel0Time (fs)Channel 0Channel 1350380Time (fs)Radius of gyration (Angstrom)20060-333334-1Raised Bevel11org.concord.modeler.PageSliderorg.concord.mw2d.models.MolecularModel05.00.05.050Temperature (K) (exaggerated)Set temperature40060falsetrue{value="1.0", label="300"}{value="2.0", label="310"}{value="3.0", label="320"}{value="4.0", label="330"}{value="5.0", label="340"}{value="0.0", label="290"}Execute MW script255 0 255true16161616161616161616161616true14true14