DNA methylation in Dictyostelium discoideum

DNA methyltransferases of type Dnmt2 are a highly conserved protein family with enigmatic function. The aim of this work was to characterize DnmA, the Dnmt2 methyltransferase in Dictyostelium discoideum, and further to investigate its implication in DNA methylation and transcriptional gene silencing. The genome of the social amoeba Dictyostelium encodes DnmA as the sole DNA methyltransferase. The enzyme bears all ten characteristic DNA methyltransferase motifs in its catalytic domain. The DnmA mRNA was found by RT-PCR to be expressed during vegetative growth and down regulated during development. Investigations using fluorescence microscopy showed that both DnmA-myc and DnmA-GFP fusions predominantly localised to the nucleus. The function of DnmA remained initially unclear, but later experiment revealed that the enzyme is an active DNA methyltransferase responsible for all DNA (cytosine) methylation in Dictyostelium. Neither in gel retardation assays, nor by the yeast two hybrid system, clues on the functionality of DnmA could be obtained. However, immunological detection of the methylation mark with an α - 5mC antibody gave initial evidence that the DNA of Dictyostelium was methylated. Furthermore, addition of 5-aza-cytidine as demethylating agent to the Dictyostelium medium and subsequent in vitro incubation of the DNA isolated from these cells with recombinant DnmA showed that the enzyme binds slightly better to this target DNA. In order to investigate further the function of the protein, a gene knock-out for dnmA was generated. The gene was successfully disrupted by homologous recombination, the knock-out strain, however, did not show any obvious phenotype under normal laboratory conditions. To identify specific target sequences for DNA methylation, a microarray analysis was carried out. Setting a threshold of at least 1.5 fold for differences in the strength of gene expression, several such genes in the knock-out strain were chosen for further investigation. Among the up-regulated genes were the ESTs representing the gag and the RT genes respectively of the retrotransposon skipper. In addition Northern blot analysis confirmed the up-regulation of skipper in the DnmA knock-out strain. Bisufite treatment and sequencing of specific DNA stretches from skipper revealed that DnmA is responsible for methylation of mostly asymmetric cytosines. Together with skipper, DIRS-1 retrotransposon was found later also to be methylated but was not present on the microarray. Furthermore, skipper transcription was also up-regulated in strains that had genes disrupted encoding components of the RNA interference pathway. In contrast, DIRS 1 expression was not affected by a loss of DnmA but was strongly increased in the strain that had the RNA directed RNA polymerase gene rrpC disrupted. Strains generated by propagating the usual wild type Ax2 and the DnmA knock-out cells over 16 rounds in development were analyzed for transposon activity. Northern blot analysis revealed activation for skipper expression, but not for DIRS-1. A large number of siRNAs were found to be correspondent to the DIRS-1 sequence, suggesting concerted regulation of DIRS-1 expression by RNAi and DNA methylation. In contrast, no siRNAs corresponding to the standard skipper element were found. The data show that DNA methylation plays a crucial role in epigenetic gene regulation in Dictyostelium and that different, partially overlapping mechanisms control transposon silencing for skipper and DIRS-1. To elucidate the mechanism of targeting the protein to particular genes in the Dictyostelium genome, some more genes which were up-regulated in the DnmA knock-out strain were analyzed by bisulfite sequencing. The chosen genes are involved in the multidrug response in other species, but their function in Dictyostelium is uncertain. Bisulfite data showed that two of these genes were methylated at asymmetrical C-residues in the wild type, but not in DnmA knock-out cells. This suggested that DNA methylation in Dictyostelium is involved not only in transposon regulation but also in transcriptional silencing of specific genes.

Ground state correlation energy of the Be-sequence for Z = 4-20 in MCDF approximation

The ground state (J = 0) electronic correlation energy of the 4-electron Be-sequence is calculated in the Multi-Configuration Dirac-Fock approximation for Z = 4-20. The 4 electrons were distributed over the configurations arising from the 1s, 2s, 2p, 3s, 3p and 3d orbitals. Theoretical values obtained here are in good agreement with experimental correlation energies.

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.