Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Biophysical characterization of tubulin tyrosine ligase-like proteins on microtubule dynamics

Report for the scientific sojourn carried out at the Cell Biology and Biophysics Unit from the National Institutes of Health, from 2010 to 2012.

Functional dynamics of PDZ binding domains: a normal-mode analysis.

Postsynaptic density-95/disks large/zonula occludens-1 (PDZ) domains are relatively small (80-120 residues) protein binding modules central in the organization of receptor clusters and in the association of cellular proteins. Their main function is to bind C-terminals of selected proteins that are recognized through specific amino acids in their carboxyl end. Binding is associated with a deformation of the PDZ native structure and is responsible for dynamical changes in regions not in direct contact with the target. We investigate how this deformation is related to the harmonic dynamics of the PDZ structure and show that one low-frequency collective normal mode, characterized by the concerted movements of different secondary structures, is involved in the binding process. Our results suggest that even minimal structural changes are responsible for communication between distant regions of the protein, in agreement with recent NMR experiments. Thus, PDZ domains are a very clear example of how collective normal modes are able to characterize the relation between function and dynamics of proteins, and to provide indications on the precursors of binding/unbinding events.

Telomere length dynamics in human lymphocyte subpopulations measured by flow cytometry.

To measure the average length of telomere repeats at chromosome ends in individual cells we developed a flow cytometry method using fluorescence in situ hybridization (flow FISH) with labeled peptide nucleic acid (PNA) probes. Results of flow FISH measurements correlated with results of conventional telomere length measurements by Southern blot analysis (R = 0.9). Consistent differences in telomere length in CD8+ T-cell subsets were identified. Naive and memory CD4+ T lymphocytes in normal adults differed by around 2.5 kb in telomere length, in agreement with known replicative shortening of telomeres in lymphocytes in vivo. T-cell clones grown in vitro showed stabilization of telomere length after an initial decline and rare clones capable of growing beyond 100 population doublings showed variable telomere length. These results show that flow FISH can be used to measure specific nucleotide repeat sequences in single cells and indicate that the very large replicative potential of lymphocytes is only indirectly related to telomere length.

The role of computational fluid dynamics in the management of unruptured intracranial aneurysms: a clinicians' view

Objective: The importance of hemodynamics in the etiopathogenesis of intracranial aneurysms (IAs) is widely accepted.Computational fluid dynamics (CFD) is being used increasingly for hemodynamic predictions. However, alogn with thecontinuing development and validation of these tools, it is imperative to collect the opinion of the clinicians. Methods: A workshopon CFD was conducted during the European Society of Minimally Invasive Neurological Therapy (ESMINT) Teaching Course,Lisbon, Portugal. 36 delegates, mostly clinicians, performed supervised CFD analysis for an IA, using the @neuFuse softwaredeveloped within the European project @neurIST. Feedback on the workshop was collected and analyzed. The performancewas assessed on a scale of 1 to 4 and, compared with experts’ performance. Results: Current dilemmas in the management ofunruptured IAs remained the most important motivating factor to attend the workshop and majority of participants showedinterest in participating in a multicentric trial. The participants achieved an average score of 2.52 (range 0–4) which was 63% (range 0–100%) of an expert user. Conclusions: Although participants showed a manifest interest in CFD, there was a clear lack ofawareness concerning the role of hemodynamics in the etiopathogenesis of IAs and the use of CFD in this context. More effortstherefore are required to enhance understanding of the clinicians in the subject.

Leveraging explicitly disclosed location information to understand tourist dynamics: a case study

In recent years, the large deployment of mobile devices has led to a massiveincrease in the volume of records of where people have been and when they were there.The analysis of these spatio-temporal data can supply high-level human behaviorinformation valuable to urban planners, local authorities, and designer of location-basedservices. In this paper, we describe our approach to collect and analyze the history ofphysical presence of tourists from the digital footprints they publicly disclose on the web.Our work takes place in the Province of Florence in Italy, where the insights on thevisitors’ flows and on the nationalities of the tourists who do not sleep in town has beenlimited to information from survey-based hotel and museums frequentation. In fact, mostlocal authorities in the world must face this dearth of data on tourist dynamics. In thiscase study, we used a corpus of geographically referenced photos taken in the provinceby 4280 photographers over a period of 2 years. Based on the disclosure of the locationof the photos, we design geovisualizations to reveal the tourist concentration and spatiotemporalflows. Our initial results provide insights on the density of tourists, the points ofinterests they visit as well as the most common trajectories they follow.

Bilinear models for spatio-temporal point distribution analysis: application to extrapolation of left ventricular, biventricular and whole heart cardiac dynamics

In this work we describe the usage of bilinear statistical models as a means of factoring the shape variability into two components attributed to inter-subject variation and to the intrinsic dynamics of the human heart. We show that it is feasible to reconstruct the shape of the heart at discrete points in the cardiac cycle. Provided we are given a small number of shape instances representing the same heart atdifferent points in the same cycle, we can use the bilinearmodel to establish this. Using a temporal and a spatial alignment step in the preprocessing of the shapes, around half of the reconstruction errors were on the order of the axial image resolution of 2 mm, and over 90% was within 3.5 mm. From this, weconclude that the dynamics were indeed separated from theinter-subject variability in our dataset.

Investigation and monitoring of rock slope instabilities in Norway by terrestrial laser scanning

Long-range Terrestrial Laser Scanning (TLS) is widely used in studies on rock slope instabilities. TLS point clouds allow the creation of high-resolution digital elevation models for detailed mapping of landslide morphologies and the measurement of the orientation of main discontinuities. Multi-temporal TLS datasets enable the quantification of slope displacements and rockfall volumes. We present three case studies using TLS for the investigation and monitoring of rock slope instabilities in Norway: 1) the analysis of 3D displacement of the Oksfjellet rock slope failure (Troms, northern Norway); 2) the detection and quantification of rockfalls along the sliding surfaces and at the front of the Kvitfjellet rock slope instability (Møre og Romsdal, western Norway); 3) the analysis of discontinuities and rotational movements of an unstable block at Stampa (Sogn og Fjordane, western Norway). These case studies highlight the possibilities but also limitations of TLS in investigating and monitoring unstable rock slopes.

Dynamics of diversity within the Mossos d’Esquadra

This chapter offers a case-study of diversity and diversity policies withinthe Mossos d’Esquadra, the police force of the Catalan autonomous communityin Spain. The case is described in a comprehensive way (includingpolicies in all relevant policy areas: recruitment, retention, and promotion)and at the same time analyzed with a new analytical framework (includingthe definition of diversity, the motivation for diversity within the organisation,and the facilitation of diversity within the organisation with policies).The goal of the chapter is twofold. First, offer a deeper understanding ofthe dynamics of diversity within this police force. Second, demonstrate theacademic potential of this new analytical framework.

The dynamics of the Swiss electricity market : three essays

Executive Summary Electricity is crucial for modern societies, thus it is important to understand the behaviour of electricity markets in order to be prepared to face the consequences of policy changes. The Swiss electricity market is now in a transition stage from a public monopoly to a liberalised market and it is undergoing an "emergent" liberalisation - i.e. liberalisation taking place without proper regulation. The withdrawal of nuclear capacity is also being debated. These two possible changes directly affect the mechanisms for capacity expansion. Thus, in this thesis we concentrate on understanding the dynamics of capacity expansion in the Swiss electricity market. A conceptual model to help understand the dynamics of capacity expansion in the Swiss electricity market is developed an explained in the first essay. We identify a potential risk of imports dependence. In the second essay a System Dynamics model, based on the conceptual model, is developed to evaluate the consequences of three scenarios: a nuclear phase-out, the implementation of a policy for avoiding imports dependence, and the combination of both. We conclude that the Swiss market is not well prepared to face unexpected changes of supply and demand, and we identify a risk of imports dependence, mainly in the case of a nuclear phase-out. The third essay focus on the opportunity cost of hydro-storage power generation, one of the main generation sources in Switzerland. We use and extended version of our model to test different policies for assigning an opportunity cost to hydro-storage power generation. We conclude that the preferred policies are different for different market participants and depend on market structure.

Early detection of microcirculatory perfusion changes with a high resolution, real time laser Doppler imaging camera--frostbite case study.

A 41-year-old male presented with severe frostbite that was monitored clinically and with a new laser Doppler imaging (LDI) camera that records arbitrary microcirculatory perfusion units (1-256 arbitrary perfusion units (APU's)). LDI monitoring detected perfusion differences in hand and foot not seen visually. On day 4-5 after injury, LDI showed that while fingers did not experience any significant perfusion change (average of 31±25 APUs on day 5), the patient's left big toe did (from 17±29 APUs day 4 to 103±55 APUs day 5). These changes in regional perfusion were not detectable by visual examination. On day 53 postinjury, all fingers with reduced perfusion by LDI were amputated, while the toe could be salvaged. This case clearly demonstrates that insufficient microcirculatory perfusion can be identified using LDI in ways which visual examination alone does not permit, allowing prognosis of clinical outcomes. Such information may also be used to develop improved treatment approaches.

Establishment of a reliable laser-Induced Choroidal Neovascularization Animal Model

Purpose: Pathologic choroidal neovascularizations (CNV) are implicated in the wet form of age-related macular degeneration (ARMD). Abnormal vessel growth is also observed in disease when hypoxia and/or inflammation occur. Our goal is to establish a standard protocol of laser-induced CNV in mice that have different levels of pigmentation to identify the most reliable animal model.Methods: CNV was induced by 4 burns around the optic disk, using a green argon laser (100μm diameter spot size; 0,05 sec. duration) in C57/Bl6, DBA/1 and Balb/c to ascertain the efficacy of the method in function of retina pigmentation. Five different intensities were tested and Bruch's membrane disruption was identified by the appearance of a bubble at the site of photocoagulation. Fluorescein angiographies (FA) were undertaken 14 days post lesion and CNV area was quantified by immunohistochemistry on cryosections.Results: CNV retina area was related to spot intensity after laser injury. While 180mW and 200mW do not induce reliable CNV (respectively 27.85±0.35% and 29±1.67% of the retina surface), 260mW is required to induce 51,07±8.52% of CNV in C57/Bl6 mice. For the DBA/1 strain, less pigmented, 200mW was sufficient to induce 49.35±3.9% of CNV, indicating that lower intensity are required to induce CNV. Furthermore, an intensity of 180mW induced greater CNV (35.55±6.01%) than in C57/Bl6 mice. Nevertheless, laser did not induce reproducible 50% CNV in Balb/c albino mice for all intensities tested. Isolectin-B4 and GFAP stainings revealed neovessel formation and photoreceptor (PR) degeneration at the impact site. The presence of glia was observed throughout all the retinal layers and angiograms showed fluorescein leakage in pigmented mice.Conclusions: The establishment of a standard protocol to induce CNV and subsequent PR degeneration is of prime importance for the use of the laser-induced CNV model and will allow to evaluate the therapeutic potency of agents to prevent CNV and retinal degeneration.

Financing Constraints, Irreversibility, and Investment Dynamics

We develop a model of an industry with many heterogeneous firms that face both financing constraints and irreversibility constraints. The financing constraint implies that firms cannot borrow unless the debt is secured by collateral; the irreversibility constraint that they can only sell their fixed capital by selling their business. We use this model to examine the cyclical behavior of aggregate fixed investment, variable capital investment, and output in the presence of persistent idiosyncratic and aggregate shocks. Our model yields three main results. First, the effect of the irreversibility constraint on fixed capital investment is reinforced by the financing constraint. Second, the effect of the financing constraint on variable capital investment is reinforced by the irreversibility constraint. Finally, the interaction between the two constraints is key for explaining why input inventories and material deliveries of US manufacturing firms are so volatile and procyclical, and also why they are highly asymmetrical over the business cycle.