Protein Folding Simulations: Confinement, External Fields and Sequence Design

The present Thesis looks at the problem of protein folding using Monte Carlo and Langevin simulations, three topics in protein folding have been studied: 1) the effect of confining potential barriers, 2) the effect of a static external field and 3) the design of amino acid sequences which fold in a short time and which have a stable native state (global minimum). Regarding the first topic, we studied the confinement of a small protein of 16 amino acids known as 1NJ0 (PDB code) which has a beta-sheet structure as a native state. The confinement of proteins occurs frequently in the cell environment. Some molecules called Chaperones, present in the cytoplasm, capture the unfolded proteins in their interior and avoid the formation of aggregates and misfolded proteins. This mechanism of confinement mediated by Chaperones is not yet well understood. In the present work we considered two kinds of potential barriers which try to mimic the confinement induced by a Chaperon molecule. The first kind of potential was a purely repulsive barrier whose only effect is to create a cavity where the protein folds up correctly. The second kind of potential was a barrier which includes both attractive and repulsive effects. We performed Wang-Landau simulations to calculate the thermodynamical properties of 1NJ0. From the free energy landscape plot we found that 1NJ0 has two intermediate states in the bulk (without confinement) which are clearly separated from the native and the unfolded states. For the case of the purely repulsive barrier we found that the intermediate states get closer to each other in the free energy landscape plot and eventually they collapse into a single intermediate state. The unfolded state is more compact, compared to that in the bulk, as the size of the barrier decreases. For an attractive barrier modifications of the states (native, unfolded and intermediates) are observed depending on the degree of attraction between the protein and the walls of the barrier. The strength of the attraction is measured by the parameter $\epsilon$. A purely repulsive barrier is obtained for $\epsilon=0$ and a purely attractive barrier for $\epsilon=1$. The states are changed slightly for magnitudes of the attraction up to $\epsilon=0.4$. The disappearance of the intermediate states of 1NJ0 is already observed for $\epsilon =0.6$. A very high attractive barrier ($\epsilon \sim 1.0$) produces a completely denatured state. In the second topic of this Thesis we dealt with the interaction of a protein with an external electric field. We demonstrated by means of computer simulations, specifically by using the Wang-Landau algorithm, that the folded, unfolded, and intermediate states can be modified by means of a field. We have found that an external field can induce several modifications in the thermodynamics of these states: for relatively low magnitudes of the field ($<2.06 \times 10^8$ V/m) no major changes in the states are observed. However, for higher magnitudes than ($6.19 \times 10^8$ V/m) one observes the appearance of a new native state which exhibits a helix-like structure. In contrast, the original native state is a $\beta$-sheet structure. In the new native state all the dipoles in the backbone structure are aligned parallel to the field. The design of amino acid sequences constitutes the third topic of the present work. We have tested the Rate of Convergence criterion proposed by D. Gridnev and M. Garcia ({\it work unpublished}). We applied it to the study of off-lattice models. The Rate of Convergence criterion is used to decide if a certain sequence will fold up correctly within a relatively short time. Before the present work, the common way to decide if a certain sequence was a good/bad folder was by performing the whole dynamics until the sequence got its native state (if it existed), or by studying the curvature of the potential energy surface. There are some difficulties in the last two approaches. In the first approach, performing the complete dynamics for hundreds of sequences is a rather challenging task because of the CPU time needed. In the second approach, calculating the curvature of the potential energy surface is possible only for very smooth surfaces. The Rate of Convergence criterion seems to avoid the previous difficulties. With this criterion one does not need to perform the complete dynamics to find the good and bad sequences. Also, the criterion does not depend on the kind of force field used and therefore it can be used even for very rugged energy surfaces.

Soil, water and nutrient loss under conventional and organic vegetable production managed in small farms

Agricultural systems with conventional tillage and intensive use of agrochemicals, especially those on high slopes and with shallow soils, have the potential to release pollutants. This study aimed at evaluating the soil, water and nutrient lost via agricultural runoff in large plots (small catchments) under conventional and organic farming of vegetables as well as under forest (control) system in a Cambisol in the Campestre catchment. Samples of runoff were collected biweekly for one year through a Coshocton wheel. The soil and water losses from the conventional farming were 218 and 6 times higher, respectively, than forest. Under organic farming the soil and water losses were 12 and 4 times higher, respectively, than forest. However the soil losses (0.5 to 114 kg ha^(−1) year^(−1)) are considered low in agronomy but environmentally represent a potential source of surface water contamination by runoff associated pollutants. The concentrations and losses of all forms of phosphorus (P) were higher in the conventional system (9.5, 0.9 and 0.3 mg L^(−1) of total P for conventional, organic and forest systems, respectively), while the organic system had the highest concentrations and losses of soluble nitrogen (4.7, 38.6 and 0.4 mg L^(−1) of NO_3-N, respectively). The percentage of bioavailable P was proportionally higher in the organic system (91% of total P lost was as bioavailable P), indicating greater potential for pollution in the short term.

Loss of wobble uridine modification in tRNA anticodons interferes with TOR pathway signaling

Previous work in yeast has suggested that modification of tRNAs, in particular uridine bases in the anticodon wobble position (U34), is linked to TOR (target of rapamycin) signaling. Hence, U34 modification mutants were found to be hypersensitive to TOR inhibition by rapamycin. To study whether this involves inappropriate TOR signaling, we examined interaction between mutations in TOR pathway genes (tip41Δ, sap190Δ, ppm1Δ, rrd1Δ) and U34 modification defects (elp3Δ, kti12Δ, urm1Δ, ncs2Δ) and found the rapamycin hypersensitivity in the latter is epistatic to drug resistance of the former. Epistasis, however, is abolished in tandem with a gln3Δ deletion, which inactivates transcription factor Gln3 required for TOR-sensitive activation of NCR (nitrogen catabolite repression) genes. In line with nuclear import of Gln3 being under control of TOR and dephosphorylation by the Sit4 phosphatase, we identify novel TOR-sensitive sit4 mutations that confer rapamycin resistance and importantly, mislocalise Gln3 when TOR is inhibited. This is similar to gln3Δ cells, which abolish the rapamycin hypersensitivity of U34 modification mutants, and suggests TOR deregulation due to tRNA undermodification operates through Gln3. In line with this, loss of U34 modifications (elp3Δ, urm1Δ) enhances nuclear import of and NCR gene activation (MEP2, GAP1) by Gln3 when TOR activity is low. Strikingly, this stimulatory effect onto Gln3 is suppressed by overexpression of tRNAs that usually carry the U34 modifications. Collectively, our data suggest that proper TOR signaling requires intact tRNA modifications and that loss of U34 modifications impinges on the TORsensitive NCR branch via Gln3 misregulation.

Loss of Anticodon Wobble Uridine Modifications Affects tRNALys Function and Protein Levels in Saccharomyces cerevisiae

In eukaryotes, wobble uridines in the anticodons of tRNALysUUU, tRNAGluUUC and tRNAGlnUUG are modified to 5-methoxy-carbonyl-methyl-2-thio-uridine (mcm5s2U). While mutations in subunits of the Elongator complex (Elp1-Elp6), which disable mcm5 side chain formation, or removal of components of the thiolation pathway (Ncs2/Ncs6, Urm1, Uba4) are individually tolerated, the combination of both modification defects has been reported to have lethal effects on Saccharomyces cerevisiae. Contrary to such absolute requirement of mcm5s2U for viability, we demonstrate here that in the S. cerevisiae S288C-derived background, both pathways can be simultaneously inactivated, resulting in combined loss of tRNA anticodon modifications (mcm5U and s2U) without a lethal effect. However, an elp3 disruption strain displays synthetic sick interaction and synergistic temperature sensitivity when combined with either uba4 or urm1 mutations, suggesting major translational defects in the absence of mcm5s2U modifications. Consistent with this notion, we find cellular protein levels drastically decreased in an elp3uba4 double mutant and show that this effect as well as growth phenotypes can be partially rescued by excess of tRNALysUUU. These results may indicate a global translational or protein homeostasis defect in cells simultaneously lacking mcm5 and s2 wobble uridine modification that could account for growth impairment and mainly originates from tRNALysUUU hypomodification and malfunction.