Coarse-grained peptide models: conformational sampling, peptide association and dynamical properties for peptides

In dieser Arbeit wird ein vergröbertes (engl. coarse-grained, CG) Simulationsmodell für Peptide in wässriger Lösung entwickelt. In einem CG Verfahren reduziert man die Anzahl der Freiheitsgrade des Systems, so dass manrngrössere Systeme auf längeren Zeitskalen untersuchen kann. Die Wechselwirkungspotentiale des CG Modells sind so aufgebaut, dass die Peptid Konformationen eines höher aufgelösten (atomistischen) Modells reproduziert werden.rnIn dieser Arbeit wird der Einfluss unterschiedlicher bindender Wechsel-rnwirkungspotentiale in der CG Simulation untersucht, insbesondere daraufhin,rnin wie weit das Konformationsgleichgewicht der atomistischen Simulation reproduziert werden kann. Im CG Verfahren verliert man per Konstruktionrnmikroskopische strukturelle Details des Peptids, zum Beispiel, Korrelationen zwischen Freiheitsgraden entlang der Peptidkette. In der Dissertationrnwird gezeigt, dass diese “verlorenen” Eigenschaften in einem Rückabbildungsverfahren wiederhergestellt werden können, in dem die atomistischen Freiheitsgrade wieder in die CG-Strukturen eingefügt werden. Dies gelingt, solange die Konformationen des CG Modells grundsätzlich gut mit der atomistischen Ebene übereinstimmen. Die erwähnten Korrelationen spielen einerngrosse Rolle bei der Bildung von Sekundärstrukturen und sind somit vonrnentscheidender Bedeutung für ein realistisches Ensemble von Peptidkonformationen. Es wird gezeigt, dass für eine gute Übereinstimmung zwischen CG und atomistischen Kettenkonformationen spezielle bindende Wechselwirkungen wie zum Beispiel 1-5 Bindungs- und 1,3,5-Winkelpotentiale erforderlich sind. Die intramolekularen Parameter (d.h. Bindungen, Winkel, Torsionen), die für kurze Oligopeptide parametrisiert wurden, sind übertragbarrnauf längere Peptidsequenzen. Allerdings können diese gebundenen Wechselwirkungen nur in Kombination mit solchen nichtbindenden Wechselwirkungspotentialen kombiniert werden, die bei der Parametrisierung verwendet werden, sind also zum Beispiel nicht ohne weiteres mit einem andere Wasser-Modell kombinierbar. Da die Energielandschaft in CG-Simulationen glatter ist als im atomistischen Modell, gibt es eine Beschleunigung in der Dynamik. Diese Beschleunigung ist unterschiedlich für verschiedene dynamische Prozesse, zum Beispiel für verschiedene Arten von Bewegungen (Rotation und Translation). Dies ist ein wichtiger Aspekt bei der Untersuchung der Kinetik von Strukturbildungsprozessen, zum Beispiel Peptid Aggregation.rn

Charakterisierung einer pflanzlichen Variante des Cytoskelettproteins gamma-Tubulin durch Strukturmodellierung, Überexpression und Analyse der Interaktionspartner

Wie alle Eukaryoten besitzen auch höhere Pflanzen ein mikrotubuläres Cytoskelett. Einige Funktionen dieses Cytoskeletts sind relativ stark konserviert, andere dagegen scheinen sehr pflanzenspezifisch zu sein. Dies betrifft insbesondere charakteristische mikrotubuläre Netzwerke, die bei der Neubildung und der Verstärkung der Zellwände wichtige Rollen übernehmen. Wie der Aufbau dieser Netzwerke kontrolliert wird, ist bisher relativ unklar. Typische Mikrotubuli organisierende Zentren (MTOC), insbesondere Centrosomen oder Spindelpolkörper, sind bei höheren Pflanzen nicht beobachtet worden. Von pilzlichen und tierischen Organismen weiß man, dass gamma-Tubulin (gTUB) mit seinen assoziierten Proteinen in den MTOC bei der Nukleation von Mikrotubuli eine Schlüsselfunktion hat. Dieses Mitglied der Tubulin-Superfamilie wird aber auch in Pflanzen gefunden, dessen genaue Funktion bisher unbekannt ist. Zu Beginn der Arbeit wurden mittels in silico Berechnungen Strukturmodelle des pflanzlichen gTUBs aus Nicotiana tabacum erarbeitet, da die Struktur, die zu einem Verständnis der pflanzlichen Wachstumsregulation beitragen könnte, bisher unbekannt ist. Auf Grundlage der bioinformatischen Daten konnte für weitere Studien eine notwendige gTUB-Deletionsmutante entwickelt werden. Für Röntgendiffraktionsstudien und gTUB-Interaktionspartneranalysen war die Verfügbarkeit verhältnismäßig großer Proteinmengen notwendig. Die Expression der gTUB-Volllängensequenz in gelöster und aktiver Form stellte einen immanent wichtigen Zwischenschritt dar. Das Escherichia coli T7/lacO-Expressionssystem lieferte, trotz vielversprechender Erfolge in der Vergangenheit, kein gelöstes rekombinantes gTUB. So wurden zwar verhältnismäßig hohe Expressionsraten erzielt, aber das rekombinante gTUB lag quantitativ als Inclusion bodies vor. Eine Variationen der Expressionsparameter sowie umfangreiche Versuche mittels verschiedenster Konstrukte sowie potentiell die Löslichkeit erhöhenden Tags gTUB in gelöster Form in E. coli zu exprimieren blieben erfolglos. Eine Denaturierung der Inclusion bodies und Rückfaltung wurde aufgrund der wohl bei der Tubulinfaltung notwendigen komplexeren Chaperone sowie thermodynamischer Überlegungen ausgeschlossen. Die höher evolvierte Chaperonausstattung war ein Hauptgrund für die Verwendung der eukaryotischen Hefe-Expressionssysteme K. lactis und des S. cerevisiae-Stammes FGY217 zur gTUB-Expression. So konnten nach der Selektion nur transgene Hefe-Zellen dokumentiert werden, die die gTUB-Expressionskassette nachweislich an der vorgesehenen Zielposition in ihrem Genom integrierten, aber keine dokumentierbare Expression zeigten. Die wahrscheinlichste Begründung hierfür ist, dass ein erhöhter intrazellulärer gTUB-Titer mit dem Zellwachstum und der Zellteilung dieser eukaryotischen Organismen interferierte und durch Rückkopplungen die rekombinante gTUB-CDS aus N. tabacum ausgeschaltet wurde. Der Versuch einer transienten gTUB-Überexpression in differenzierten Blattgeweben höherer Pflanzen war eine logische Konsequenz aus den vorherigen Ergebnissen und lieferte, wenn auch nicht die für eine Proteinkristallisation notwendigen Mengen, gelöstes gTUB. Bestrebungen einer stabilen Transfektion von A. thaliana oder BY-2-Zellkulturen mit einer gTUB-CDS lieferten keine transgenen Organismen, was starke Interferenzen der rekombinanten gTUB-CDS in den Zellen vermuten lies. Transfektionsversuche mit nur GFP tragenden Konstrukten ergaben hingegen eine hohe Anzahl an transgenen Organismen, die auch verhältnismäßig starke Expressionsraten zeigten. Die erzielten Proteinmengen bei der transienten gTUB-Überexpression in N. benthamiana Blattgeweben, in Co-Expression mit dem Posttransriptional Gene Silencing-Suppressorprotein p19, waren für einen Pull-Down sowie eine massenspektroskopische Analyse der Interaktionspartner ausreichend und ergaben Befunde. Eine abschließende Auswertung des erarbeiteten massenspektroskopischen Datensatzes wird jedoch erst dann möglich sein, wenn das Tabak-Proteom vollständig sequenziert ist. Die Erweiterung der bestehenden pflanzlichen Vergleichsdatenbanken um das bisher bekannte Tabak-Proteom vervielfachte die Anzahl der in dieser Studie identifizierten gTUB-Interaktionspartner. Interaktionen mit dem TCP1-Chaperon untermauern die Hypothese der zur Faltung pflanzlichen gTUBs notwendigen Chaperone. Beobachtete gTUB-Degradationsmuster in Verbindung mit Interaktionen des 26S-Proteasoms deuten auf eine Gegenregulationen bei erhöhtem gTUB-Titer auf Proteinebene hin. Da Blattgewebe selbst nur noch über eine sehr geringe und inhomogene Teilungsaktivität verfügen ist diese Regulation hoch spannend. Auch konnte durch Co-Expression des PTGS-Suppressorproteins p19 gezeigt werden, dass bei der gTUB-Expression eine Regulation auf RNA-Ebene erfolgt.

Evaluation of 3D cell culture systems for host-pathogen interaction studies

Traditional cell culture models have limitations in extrapolating functional mechanisms that underlie strategies of microbial virulence. Indeed during the infection the pathogens adapt to different tissue-specific environmental factors. The development of in vitro models resembling human tissue physiology might allow the replacement of inaccurate or aberrant animal models. Three-dimensional (3D) cell culture systems are more reliable and more predictive models that can be used for the meaningful dissection of host–pathogen interactions. The lung and gut mucosae often represent the first site of exposure to pathogens and provide a physical barrier against their entry. Within this context, the tracheobronchial and small intestine tract were modelled by tissue engineering approach. The main work was focused on the development and the extensive characterization of a human organotypic airway model, based on a mechanically supported co-culture of normal primary cells. The regained morphological features, the retrieved environmental factors and the presence of specific epithelial subsets resembled the native tissue organization. In addition, the respiratory model enabled the modular insertion of interesting cell types, such as innate immune cells or multipotent stromal cells, showing a functional ability to release pertinent cytokines differentially. Furthermore this model responded imitating known events occurring during the infection by Non-typeable H. influenzae. Epithelial organoid models, mimicking the small intestine tract, were used for a different explorative analysis of tissue-toxicity. Further experiments led to detection of a cell population targeted by C. difficile Toxin A and suggested a role in the impairment of the epithelial homeostasis by the bacterial virulence machinery. The described cell-centered strategy can afford critical insights in the evaluation of the host defence and pathogenic mechanisms. The application of these two models may provide an informing step that more coherently defines relevant molecular interactions happening during the infection.

Convergence results for stochastic particle systems with social interaction

We consider stochastic individual-based models for social behaviour of groups of animals. In these models the trajectory of each animal is given by a stochastic differential equation with interaction. The social interaction is contained in the drift term of the SDE. We consider a global aggregation force and a short-range repulsion force. The repulsion range and strength gets rescaled with the number of animals N. We show that for N tending to infinity stochastic fluctuations disappear and a smoothed version of the empirical process converges uniformly towards the solution of a nonlinear, nonlocal partial differential equation of advection-reaction-diffusion type. The rescaling of the repulsion in the individual-based model implies that the corresponding term in the limit equation is local while the aggregation term is non-local. Moreover, we discuss the effect of a predator on the system and derive an analogous convergence result. The predator acts as an repulsive force. Different laws of motion for the predator are considered.

Modelling and simulating in systems biology: an approach based on multi-agent systems

Systems Biology is an innovative way of doing biology recently raised in bio-informatics contexts, characterised by the study of biological systems as complex systems with a strong focus on the system level and on the interaction dimension. In other words, the objective is to understand biological systems as a whole, putting on the foreground not only the study of the individual parts as standalone parts, but also of their interaction and of the global properties that emerge at the system level by means of the interaction among the parts. This thesis focuses on the adoption of multi-agent systems (MAS) as a suitable paradigm for Systems Biology, for developing models and simulation of complex biological systems. Multi-agent system have been recently introduced in informatics context as a suitabe paradigm for modelling and engineering complex systems. Roughly speaking, a MAS can be conceived as a set of autonomous and interacting entities, called agents, situated in some kind of nvironment, where they fruitfully interact and coordinate so as to obtain a coherent global system behaviour. The claim of this work is that the general properties of MAS make them an effective approach for modelling and building simulations of complex biological systems, following the methodological principles identified by Systems Biology. In particular, the thesis focuses on cell populations as biological systems. In order to support the claim, the thesis introduces and describes (i) a MAS-based model conceived for modelling the dynamics of systems of cells interacting inside cell environment called niches. (ii) a computational tool, developed for implementing the models and executing the simulations. The tool is meant to work as a kind of virtual laboratory, on top of which kinds of virtual experiments can be performed, characterised by the definition and execution of specific models implemented as MASs, so as to support the validation, falsification and improvement of the models through the observation and analysis of the simulations. A hematopoietic stem cell system is taken as reference case study for formulating a specific model and executing virtual experiments.

Noxious stimulation response index: a novel anesthetic state index based on hypnotic-opioid interaction

The noxious stimulation response index (NSRI) is a novel anesthetic depth index ranging between 100 and 0, computed from hypnotic and opioid effect-site concentrations using a hierarchical interaction model. The authors validated the NSRI on previously published data.

Video-based quantification of body movement during social interaction indicates the severity of negative symptoms in patients with schizophrenia

In schizophrenia, nonverbal behavior, including body movement, is of theoretical and clinical importance. Although reduced nonverbal expressiveness is a major component of the negative symptoms encountered in schizophrenia, few studies have objectively assessed body movement during social interaction. In the present study, 378 brief, videotaped role-play scenes involving 27 stabilized outpatients diagnosed with paranoid-type schizophrenia were analyzed using Motion Energy Analysis (MEA). This method enables the objective measuring of body movement in conjunction with ordinary video recordings. Correlations between movement parameters (percentage of time in movement, movement speed) and symptom ratings from independent PANSS interviews were calculated. Movement parameters proved to be highly reliable. In keeping with predictions, reduced movement and movement speed correlated with negative symptoms. Accordingly, in patients who exhibited noticeable movement for less than 20% of the observation time, prominent negative symptoms were highly probable. As a control measure, the percentage of movement exhibited by the patients during role-play scenes was compared to that of their normal interactants. Patients with negative symptoms differed from normal interactants by showing significantly reduced head and body movement. Two specific positive symptoms were possibly related to movement parameters: suspiciousness tended to correlate with reduced head movement, and the expression of unusual thought content tended to relate to increased movement. Overall, a close and theoretically meaningful association between the objective movement parameters and the symptom profiles was found. MEA appears to be an objective, reliable and valid method for quantifying nonverbal behavior, an aspect which may furnish new insights into the processes related to reduced expressiveness in schizophrenia.

Development of Dynamic Constraint Models for a Model Based Transient Calibration Process

Model based calibration has gained popularity in recent years as a method to optimize increasingly complex engine systems. However virtually all model based techniques are applied to steady state calibration. Transient calibration is by and large an emerging technology. An important piece of any transient calibration process is the ability to constrain the optimizer to treat the problem as a dynamic one and not as a quasi-static process. The optimized air-handling parameters corresponding to any instant of time must be achievable in a transient sense; this in turn depends on the trajectory of the same parameters over previous time instances. In this work dynamic constraint models have been proposed to translate commanded to actually achieved air-handling parameters. These models enable the optimization to be realistic in a transient sense. The air handling system has been treated as a linear second order system with PD control. Parameters for this second order system have been extracted from real transient data. The model has been shown to be the best choice relative to a list of appropriate candidates such as neural networks and first order models. The selected second order model was used in conjunction with transient emission models to predict emissions over the FTP cycle. It has been shown that emission predictions based on air-handing parameters predicted by the dynamic constraint model do not differ significantly from corresponding emissions based on measured air-handling parameters.

Transmission of Chlamydia trachomatis through sexual partnerships: a comparison between three individual-based models and empirical data

Chlamydia trachomatis is the most common bacterial sexually transmitted infection (STI) in many developed countries. The highest prevalence rates are found among young adults who have frequent partner change rates. Three published individual-based models have incorporated a detailed description of age-specific sexual behaviour in order to quantify the transmission of C. trachomatis in the population and to assess the impact of screening interventions. Owing to varying assumptions about sexual partnership formation and dissolution and the great uncertainty about critical parameters, such models show conflicting results about the impact of preventive interventions. Here, we perform a detailed evaluation of these models by comparing the partnership formation and dissolution dynamics with data from Natsal 2000, a population-based probability sample survey of sexual attitudes and lifestyles in Britain. The data also allow us to describe the dispersion of C. trachomatis infections as a function of sexual behaviour, using the Gini coefficient. We suggest that the Gini coefficient is a useful measure for calibrating infectious disease models that include risk structure and highlight the need to estimate this measure for other STIs.

[(99m)Tc]Demomedin C, a radioligand based on human gastrin releasing peptide(18-27): synthesis and preclinical evaluation in gastrin releasing peptide receptor-expressing models

The synthesis and preclinical evaluation of [(99m)Tc]Demomedin C in GRPR-expressing models are reported. Demomedin C resulted by coupling a Boc-protected N(4)-chelator to neuromedin C (human GRP(18-27)), which, after (99m)Tc-labeling, afforded [(99m)Tc]Demomedin C. Demomedin C showed high affinity and selectivity for the GRPR during receptor autoradiography on human cancer samples (IC(50) in nM: GRPR, 1.4 ± 0.2; NMBR, 106 ± 18; and BB(3)R, >1000). It triggered GRPR internalization in HEK-GRPR cells and Ca(2+) release in PC-3 cells (EC(50) = 1.3 nM). [(99m)Tc]Demomedin C rapidly and specifically internalized at 37 °C in PC-3 cells and was stable in mouse plasma. [(99m)Tc]Demomedin C efficiently and specifically localized in human PC-3 implants in mice (9.84 ± 0.81%ID/g at 1 h pi; 6.36 ± 0.85%ID/g at 4 h pi, and 0.41 ± 0.07%ID/g at 4 h pi block). Thus, human GRP-based radioligands, such as [(99m)Tc]Demomedin C, can successfully target GRPR-expressing human tumors in vivo while displaying attractive biological features--e.g. higher GRPR-selectivity--vs their frog-homologues.