Structural and elastic properties of DNA and chromatin

We investigate a chain consisting of two coupled worm-like chains withconstant distance between the strands. The effects due todouble-strandedness of the chain are studied. In a previous analyticalstudy of this system an intrinsic twist-stretch coupling and atendency of kinking is predicted. Even though a local twist structureis observed the predicted features are not recovered. A new model for DNA at the base-pair level is presented. Thebase-pairs are treated as flat rigid ellipsoids and thesugar-phosphate backbones are represented as stiff harmonic springs.The base-pair stacking interaction is modeled by a variant of theGay-Berne potential. It is shown by systematic coarse-graininghow the elastic constants of a worm-like chain are related to thelocal fluctuations of the base-pair step parameters. Even though a lotof microscopic details of the base-pair geometry is neglected themodel can be optimized to obtain a B-DNA conformation as ground stateand reasonable elastic properties. Moreover the model allows tosimulate much larger length scales than it is possible with atomisticsimulations due to the simplification of the force-field and inparticular due to the possibility of non-local Monte-Carlo moves. Asa first application the behavior under stretching is investigated. Inagreement with micromanipulation experiments on single DNA moleculesone observes a force-plateau in the force-extension curvescorresponding to an overstretching transition from B-DNA to aso-called S-DNA state. The model suggests a structure for S-DNA withhighly inclined base-pairs in order to enable at least partialbase-pair stacking. Finally a simple model for chromatin is introduced to study itsstructural and elastic properties. The underlying geometry of themodeled fiber is based on a crossed-linker model. The chromatosomesare treated as disk-like objects. Excluded volume and short rangenucleosomal interaction are taken into account by a variant of theGay-Berne potential. It is found that the bending rigidity and thestretching modulus of the fiber increase with more compact fibers. Fora reasonable parameterization of the fiber for physiologicalconditions and sufficiently high attraction between the nucleosomes aforce-extension curve is found similar to stretching experiments onsingle chromatin fibers. For very small stretching forces a kinkedfiber forming a loop is observed. If larger forces are applied theloop formation is stretched out and a decondensation of the fibertakes place.

Simulations of amphiphilic peptides at the air-water interface

Amphiphile Peptide, Pro-Glu-(Phe-Glu)n-Pro, Pro-Asp-(Phe-Asp)n-Pro, und Phe-Glu-(Phe-Glu)n-Phe, können so aus n alternierenden Sequenzen von hydrophoben und hydrophilen Aminosäuren konstruiert werden, dass sie sich in Monolagen an der Luft-Wasser Grenzfläche anordnen. In biologischen Systemen können Strukturen an der organisch-wässrigen Grenzfläche als Matrix für die Kristallisation von Hydroxyapatit dienen, ein Vorgang der für die Behandlung von Osteoporose verwendet werden kann. In der vorliegenden Arbeit wurden Computersimulationenrneingesetzt, um die Strukturen und die zugrunde liegenden Wechselwirkungen welche die Aggregation der Peptide auf mikroskopischer Ebene steuern, zu untersuchen. Atomistische Molekulardynamik-Simulationen von einzelnen Peptidsträngen zeigen, dass sie sich leicht an der Luft-Wasser Grenzfläche anordnen und die Fähigkeit haben, sich in β-Schleifen zu falten, selbst für relativ kurze Peptidlängen (n = 2). Seltene Ereignisse wie diese (i.e. Konformationsänderungen) erfordern den Einsatz fortgeschrittener Sampling-Techniken. Hier wurde “Replica Exchange” Molekulardynamik verwendet um den Einfluss der Peptidsequenzen zu untersuchen. Die Simulationsergebnisse zeigten, dass Peptide mit kürzeren azidischen Seitenketten (Asp vs. Glu) gestrecktere Konformationen aufwiesen als die mit längeren Seitenketten, die in der Lage waren die Prolin-Termini zu erreichen. Darüber hinaus zeigte sich, dass die Prolin-Termini (Pro vs. Phe) notwendig sind, um eine 2D-Ordnung innerhalb derrnAggregate zu erhalten. Das Peptid Pro-Asp-(Phe-Asp)n-Pro, das beide dieser Eigenschaften enthält, zeigt das geordnetste Verhalten, eine geringe Verdrehung der Hauptkette, und ist in der Lage die gebildeten Aggregate durch Wasserstoffbrücken zwischen den sauren Seitenketten zu stabilisieren. Somit ist dieses Peptid am besten zur Aggregation geeignet. Dies wurde auch durch die Beurteilung der Stabilität von experimentnah-aufgesetzten Peptidaggregaten, sowie der Neigung einzelner Peptide zur Selbstorganisation von anfänglich ungeordneten Konfigurationen unterstützt. Da atomistische Simulationen nur auf kleine Systemgrößen und relativ kurze Zeitskalen begrenzt sind, wird ein vergröbertes Modell entwickelt damit die Selbstorganisation auf einem größeren Maßstab studiert werden kann. Da die Selbstorganisation an der Grenzfläche vonrnInteresse ist, wurden existierenden Vergröberungsmethoden erweitert, um nicht-gebundene Potentiale für inhomogene Systeme zu bestimmen. Die entwickelte Methode ist analog zur iterativen Boltzmann Inversion, bildet aber das Update für das Interaktionspotential basierend auf der radialen Verteilungsfunktion in einer Slab-Geometrie und den Breiten des Slabs und der Grenzfläche. Somit kann ein Kompromiss zwischen der lokalen Flüssigketsstruktur und den thermodynamischen Eigenschaften der Grenzfläche erreicht werden. Die neue Methode wurde für einen Wasser- und einen Methanol-Slab im Vakuum demonstriert, sowie für ein einzelnes Benzolmolekül an der Vakuum-Wasser Grenzfläche, eine Anwendung die von besonderer Bedeutung in der Biologie ist, in der oft das thermodynamische/Grenzflächenpolymerisations-Verhalten zusätzlich der strukturellen Eigenschaften des Systems erhalten werden müssen. Daraufrnbasierend wurde ein vergröbertes Modell über einen Fragment-Ansatz parametrisiert und die Affinität des Peptids zur Vakuum-Wasser Grenzfläche getestet. Obwohl die einzelnen Fragmente sowohl die Struktur als auch die Wahrscheinlichkeitsverteilungen an der Grenzfläche reproduzierten, diffundierte das Peptid als Ganzes von der Grenzfläche weg. Jedoch führte eine Reparametrisierung der nicht-gebundenen Wechselwirkungen für eines der Fragmente der Hauptkette in einem Trimer dazu, dass das Peptid an der Grenzfläche blieb. Dies deutet darauf hin, dass die Kettenkonnektivität eine wichtige Rolle im Verhalten des Petpids an der Grenzfläche spielt.