Computational investigations of ion conduction mechanisms using enhanced sampling methods

Proper ion channels’ functioning is a prerequisite for a normal cell and disorders involving ion channels, or channelopathies, underlie many human diseases. Long QT syndromes (LQTS) for example may arise from the malfunctioning of hERG channel, caused either by the binding of drugs or mutations in HERG gene. In the first part of this thesis I present a framework to investigate the mechanism of ion conduction through hERG channel. The free energy profile governing the elementary steps of ion translocation in the pore was computed by means of umbrella sampling simulations. Compared to previous studies, we detected a different dynamic behavior: according to our data hERG is more likely to mediate a conduction mechanism which has been referred to as “single-vacancy-like” by Roux and coworkers (2001), rather then a “knock-on” mechanism. The same protocol was applied to a model of hERG presenting the Gly628Ser mutation, found to be cause of congenital LQTS. The results provided interesting insights about the reason of the malfunctioning of the mutant channel. Since they have critical functions in viruses’ life cycle, viral ion channels, such as M2 proton channel, are considered attractive targets for antiviral therapy. A deep knowledge of the mechanisms that the virus employs to survive in the host cell is of primary importance in the identification of new antiviral strategies. In the second part of this thesis I shed light on the role that M2 plays in the control of electrical potential inside the virus, being the charge equilibration a condition required to allow proton influx. The ion conduction through M2 was simulated using metadynamics technique. Based on our results we suggest that a potential anion-mediated cation-proton exchange, as well as a direct anion-proton exchange could both contribute to explain the activity of the M2 channel.

A novel map-matching procedure for low-sampling GPS data with applications to traffic flow analysis

An extensive sample (2%) of private vehicles in Italy are equipped with a GPS device that periodically measures their position and dynamical state for insurance purposes. Having access to this type of data allows to develop theoretical and practical applications of great interest: the real-time reconstruction of traffic state in a certain region, the development of accurate models of vehicle dynamics, the study of the cognitive dynamics of drivers. In order for these applications to be possible, we first need to develop the ability to reconstruct the paths taken by vehicles on the road network from the raw GPS data. In fact, these data are affected by positioning errors and they are often very distanced from each other (~2 Km). For these reasons, the task of path identification is not straightforward. This thesis describes the approach we followed to reliably identify vehicle paths from this kind of low-sampling data. The problem of matching data with roads is solved with a bayesian approach of maximum likelihood. While the identification of the path taken between two consecutive GPS measures is performed with a specifically developed optimal routing algorithm, based on A* algorithm. The procedure was applied on an off-line urban data sample and proved to be robust and accurate. Future developments will extend the procedure to real-time execution and nation-wide coverage.

Quantifying the OH radical by global scale NMHC measurements from the NOAA - ESRL cooperative flask sampling network

Summary PhD Thesis Jan Pollmann: This thesis focuses on global scale measurements of light reactive non-methane hydrocarbon (NMHC), in the volatility range from ethane to toluene with a special focus on ethane, propane, isobutane, butane, isopentane and pentane. Even though they only occur at the ppt level (nmol mol-1) in the remote troposphere these species can yield insight into key atmospheric processes. An analytical method was developed and subsequently evaluated to analyze NMHC from the NOAA – ERSL cooperative air sampling network. Potential analytical interferences through other atmospheric trace gases (water vapor and ozone) were carefully examined. The analytical parameters accuracy and precision were analyzed in detail. It was proven that more than 90% of the data points meet the Global Atmospheric Watch (GAW) data quality objective. Trace gas measurements from 28 measurement stations were used to derive the global atmospheric distribution profile for 4 NMHC (ethane, propane, isobutane, butane). A close comparison of the derived ethane data with previously published reports showed that northern hemispheric ethane background mixing ratio declined by approximately 30% since 1990. No such change was observed for southern hemispheric ethane. The NMHC data and trace gas data supplied by NOAA ESRL were used to estimate local diurnal averaged hydroxyl radical (OH) mixing ratios by variability analysis. Comparison of the variability derived OH with directly measured OH and modeled OH mixing ratios were found in good agreement outside the tropics. Tropical OH was on average two times higher than predicted by the model. Variability analysis was used to assess the effect of chlorine radicals on atmospheric oxidation chemistry. It was found that Cl is probably not of significant relevance on a global scale.

Portal hypertension: a comparison between portal-venous pressure measurement and ARFI measurement of liver and spleen stiffness

PURPOSE. Portal pressure is measured invasively as Hepatic Venous Pressure Gradient (HVPG) in the angiography room. Liver stiffness measured by Fibroscan was shown to correlate with HVPG values below 12 mmHg. This is not surprising, since in cirrhosis the increase of portal pressure is not directly linked with liver fibrosis and consequently to liver stiffness. We hypothesized that, given the spleen’s privileged location upstream to the whole portal system, splenic stiffness could provide relevant information about portal pressure. Aim of the study was to assess the relationship between liver and spleen stiffness measured by Virtual Touch™ (ARFI) and HVPG in cirrhotic patients. METHODS. 40 consecutive patients (30 males, mean age 62y, mean BMI=26, mean Child-Pugh A6, mean platelet count=92.000/mmc, 19 HCV+, 7 with ascites) underwent to ARFI stiffness measurement (10 valid measurements in right liver lobe both surface and centre, left lobe and 20 in the spleen) and HPVG, blindly to each other. Median ARFI values of 10 samplings on every liver area and of 20 samplings on spleen were calculated. RESULTS. Stiffness could be easily measured in all patients with ARFI, resulting a mean of 2,61±0,76, 2,5±0,62 and 2,55±0,66 m/sec in the liver areas and 3.3±0,5 m/s in the spleen. Median HPVG was 14 mmHg (range 5-27); 28 patients showed values ≥10 mmHg. A positive significant correlation was found between spleen stiffness and HPVG values (r=0.744, p<0.001). No significant correlation was found between all liver stiffness and HVPG (p>0,05). AUROC was calculated to test spleen stiffness ability in discriminating patients with HVPG ≥10. AUROC = 0.911 was obtained, with sensitivity of 69% and specificity of 91% at a cut-off of 3.26 m/s. CONCLUSION. Spleen stiffness measurement with ARFI correlates with HVPG in patients with cirrhosis, with a potential of identifying patients with clinically significant portal hypertension.

Local Earthquake Tomography by trans-dimensional Monte Carlo Sampling

In this study a new, fully non-linear, approach to Local Earthquake Tomography is presented. Local Earthquakes Tomography (LET) is a non-linear inversion problem that allows the joint determination of earthquakes parameters and velocity structure from arrival times of waves generated by local sources. Since the early developments of seismic tomography several inversion methods have been developed to solve this problem in a linearized way. In the framework of Monte Carlo sampling, we developed a new code based on the Reversible Jump Markov Chain Monte Carlo sampling method (Rj-McMc). It is a trans-dimensional approach in which the number of unknowns, and thus the model parameterization, is treated as one of the unknowns. I show that our new code allows overcoming major limitations of linearized tomography, opening a new perspective in seismic imaging. Synthetic tests demonstrate that our algorithm is able to produce a robust and reliable tomography without the need to make subjective a-priori assumptions about starting models and parameterization. Moreover it provides a more accurate estimate of uncertainties about the model parameters. Therefore, it is very suitable for investigating the velocity structure in regions that lack of accurate a-priori information. Synthetic tests also reveal that the lack of any regularization constraints allows extracting more information from the observed data and that the velocity structure can be detected also in regions where the density of rays is low and standard linearized codes fails. I also present high-resolution Vp and Vp/Vs models in two widespread investigated regions: the Parkfield segment of the San Andreas Fault (California, USA) and the area around the Alto Tiberina fault (Umbria-Marche, Italy). In both the cases, the models obtained with our code show a substantial improvement in the data fit, if compared with the models obtained from the same data set with the linearized inversion codes.

Sampling methods for low-frequency electromagnetic imaging

For the detection of hidden objects by low-frequency electromagnetic imaging the Linear Sampling Method works remarkably well despite the fact that the rigorous mathematical justification is still incomplete. In this work, we give an explanation for this good performance by showing that in the low-frequency limit the measurement operator fulfills the assumptions for the fully justified variant of the Linear Sampling Method, the so-called Factorization Method. We also show how the method has to be modified in the physically relevant case of electromagnetic imaging with divergence-free currents. We present numerical results to illustrate our findings, and to show that similar performance can be expected for the case of conducting objects and layered backgrounds.

A sampling method for detecting buried objects using electromagnetic scattering

We consider a simple (but fully three-dimensional) mathematical model for the electromagnetic exploration of buried, perfect electrically conducting objects within the soil underground. Moving an electric device parallel to the ground at constant height in order to generate a magnetic field, we measure the induced magnetic field within the device, and factor the underlying mathematics into a product of three operations which correspond to the primary excitation, some kind of reflection on the surface of the buried object(s) and the corresponding secondary excitation, respectively. Using this factorization we are able to give a justification of the so-called sampling method from inverse scattering theory for this particular set-up.

Visualisierung vergrabener Objekte: Ein Sampling-Verfahren für die Maxwell-Gleichungen

Wir untersuchen die numerische Lösung des inversen Streuproblems der Rekonstruktion der Form, Position und Anzahl endlich vieler perfekt leitender Objekte durch Nahfeldmessungen zeitharmonischer elektromagnetischer Wellen mit Hilfe von Metalldetektoren. Wir nehmen an, dass sich die Objekte gänzlich im unteren Halbraum eines unbeschränkten zweischichtigen Hintergrundmediums befinden. Wir nehmen weiter an, dass der obere Halbraum mit Luft und der untere Halbraum mit Erde gefüllt ist. Wir betrachten zuerst die physikalischen Grundlagen elektromagnetischer Wellen, aus denen wir zunächst ein vereinfachtes mathematisches Modell ableiten, in welchem wir direkt das elektromagnetische Feld messen. Dieses Modell erweitern wir dann um die Messung des elektromagnetischen Feldes von Sendespulen mit Hilfe von Empfangsspulen. Für das vereinfachte Modell entwickeln wir, unter Verwendung der Theorie des zugehörigen direkten Streuproblems, ein nichtiteratives Verfahren, das auf der Idee der sogenannten Faktorisierungsmethode beruht. Dieses Verfahren übertragen wir dann auf das erweiterte Modell. Wir geben einen Implementierungsvorschlag der Rekonstruktionsmethode und demonstrieren an einer Reihe numerischer Experimente die Anwendbarkeit des Verfahrens. Weiterhin untersuchen wir mehrere Abwandlungen der Methode zur Verbesserung der Rekonstruktionen und zur Verringerung der Rechenzeit.

Validation of a new method to assess the long-term Hypothalamic-Pituitary-Adrenal Activity using hair glucocorticoids as biomarkers. - Studies on pigs and laboratory Sprague-Dawley Rats

The evaluation of chronic activity of the hypothalamic-pituitary-adrenal (HPA) axis is critical for determining the impact of chronic stressful situations. The potential use of hair glucocorticoids as a non-invasive, retrospective, biomarker of long term HPA activity is of great interest, and it is gaining acceptance in humans and animals. However, there are still no studies in literature examining hair cortisol concentration in pigs and corticosterone concentration in laboratory rodents. Therefore, we developed and validated, for the first time, a method for measuring hair glucocorticoids concentration in commercial sows and in Sprague-Dawley rats. Our preliminary data demonstrated: 1) a validated and specific washing protocol and extraction assay method with a good sensitivity in both species; 2) the effect of the reproductive phase, housing conditions and seasonality on hair cortisol concentration in sows; 3) similar hair corticosterone concentration in male and female rats; 4) elevated hair corticosterone concentration in response to chronic stress manipulations and chronic ACTH administration, demonstrating that hair provides a good direct index of HPA activity over long periods than other indirect parameters, such adrenal or thymus weight. From these results we believe that this new non-invasive tool needs to be applied to better characterize the overall impact in livestock animals and in laboratory rodents of chronic stressful situations that negatively affect animals welfare. Nevertheless, further studies are needed to improve this methodology and maybe to develop animal models for chronic stress of high interest and translational value in human medicine.

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

Neuropeptide Y modulates steroid production of human adrenal H295R cells through Y1 receptors

Neuropeptide Y (NPY) is abundantly expressed in the nervous system and acts on target cells through NPY receptors. The human adrenal cortex and adrenal tumors express NPY receptor subtype Y1, but its function is unknown. We studied Y1-mediated signaling, steroidogenesis and cell proliferation in human adrenal NCI-H295R cells. Radioactive ligand binding studies showed that H295R cells express Y1 receptor specifically. NPY treatment of H295R cells stimulated the MEK/ERK1/2 pathway, confirming that H295R cells express functional Y1 receptors. Studies of the effect of NPY and related peptide PYY on adrenal steroidogenesis revealed a decrease in 11-deoxycortisol production. RIA measurements of cortisol from cell culture medium confirmed this finding. Co-treatment with the Y1 antagonist BIBP2336 reversed the inhibitory effect of NPY on cortisol production proving specificity of this effect. At mRNA level, NPY decreased HSD3B2 and CYP21A2 expression. However NPY revealed no effect on cell proliferation. Our data show that NPY can directly regulate human adrenal cortisol production.

Control of Steroid 21-oic Acid Synthesis by Peroxisome Proliferator-activated Receptor and Role of the Hypothalamic-Pituitary-Adrenal Axis

A previous study identified the peroxisome proliferator-activated receptor alpha (PPARalpha) activation biomarkers 21-steroid carboxylic acids 11beta-hydroxy-3,20-dioxopregn-4-en-21-oic acid (HDOPA) and 11beta,20-dihydroxy-3-oxo-pregn-4-en-21-oic acid (DHOPA). In the present study, the molecular mechanism and the metabolic pathway of their production were determined. The PPARalpha-specific time-dependent increases in HDOPA and 20alpha-DHOPA paralleled the development of adrenal cortex hyperplasia, hypercortisolism, and spleen atrophy, which was attenuated in adrenalectomized mice. Wy-14,643 activation of PPARalpha induced hepatic FGF21, which caused increased neuropeptide Y and agouti-related protein mRNAs in the hypothalamus, stimulation of the agouti-related protein/neuropeptide Y neurons, and activation of the hypothalamic-pituitary-adrenal (HPA) axis, resulting in increased adrenal cortex hyperplasia and corticosterone production, revealing a link between PPARalpha and the HPA axis in controlling energy homeostasis and immune regulation. Corticosterone was demonstrated as the precursor of 21-carboxylic acids both in vivo and in vitro. Under PPARalpha activation, the classic reductive metabolic pathway of corticosterone was suppressed, whereas an alternative oxidative pathway was uncovered that leads to the sequential oxidation on carbon 21 resulting in HDOPA. The latter was then reduced to the end product 20alpha-DHOPA. Hepatic cytochromes P450, aldehyde dehydrogenase (ALDH3A2), and 21-hydroxysteroid dehydrogenase (AKR1C18) were found to be involved in this pathway. Activation of PPARalpha resulted in the induction of Aldh3a2 and Akr1c18, both of which were confirmed as target genes through introduction of promoter luciferase reporter constructs into mouse livers in vivo. This study underscores the power of mass spectrometry-based metabolomics combined with genomic and physiologic analyses in identifying downstream metabolic biomarkers and the corresponding upstream molecular mechanisms.

Ease of Application of Medical Compression-stocking Systems for the Treatment of Venous Ulcers

OBJECTIVE: To evaluate the ease of application of two-piece, graduated, compression systems for the treatment of venous ulcers. METHODS: Four kits used to provide limb compression in the management of venous ulcers were evaluated. These have been proven to be non-inferior to various types of bandages in clinical trials. The interface pressure exerted above the ankle by the under-stocking and the complete compression system and the force required to pull the over-stocking off were assessed in vitro. Ease of application of the four kits was evaluated in four sessions by five nurses who put stockings on their own legs in a blinded manner. They expressed their assessment of the stockings using a series of visual analogue scales (VASs). RESULTS: The Sigvaris Ulcer X((R)) kit provided a mean interface pressure of 46 mmHg and required a force in the range of 60-90 N to remove it. The Mediven((R)) ulcer kit exerted the same pressure but required force in the range of 150-190 N to remove it. Two kits (SurePress((R)) Comfort and VenoTrain((R)) Ulcertec) exerted a mean pressure of only 25 mmHg and needed a force in the range of 100-160 N to remove them. Nurses judged the Ulcer X and SurePress kits easiest to apply. Application of the VenoTrain kit was found slightly more difficult. The Mediven kit was judged to be difficult to use. CONCLUSIONS: Comparison of ease of application of compression-stocking kits in normal legs revealed marked differences between them. Only one system exerted a high pressure and was easy to apply. Direct comparison of these compression kits in leg-ulcer patients is required to assess whether our laboratory findings correlate with patient compliance and ulcer healing.

Lack of prophylaxis before the onset of acute venous thromboembolism among hospitalized cancer patients: the SWIss Venous ThromboEmbolism Registry (SWIVTER)

Venous thromboembolism (VTE) prophylaxis remains underutilized, particularly in cancer patients. We explored clinical predictors of prophylaxis in hospitalized cancer patients before the onset of acute VTE.

Comparison of low-strength compression stockings with bandages for the treatment of recalcitrant venous ulcers

OBJECTIVE: To compare the proportion and rate of healing, pain, and quality of life of low-strength medical compression stockings (MCS) with traditional bandages applied for the treatment of recalcitrant venous leg ulcers. METHODS: A single-center, randomized, open-label study was performed with consecutive patients. Sigvaris prototype MCS providing 15 mm Hg-25 mm Hg at the ankle were compared with multi-layer short-stretch bandages. In both groups, pads were placed above incompetent perforating veins in the ulcer area. The initial static pressure between the dressing-covered ulcer and the pad was 29 mm Hg and 49 mm Hg with MCS and bandages, respectively. Dynamic pressure measurements showed no difference. Compression was maintained day and night and changed every week. The primary endpoint was healing within 90 days. Secondary endpoints were healing within 180 days, time to healing, pain (weekly Likert scales), and monthly quality of life (ChronIc Venous Insufficiency Quality of Life [CIVIQ] questionnaire). RESULTS: Of 74 patients screened, 60 fulfilled the selection criteria and 55 completed the study; 28 in the MCS and 27 in the bandage group. Ulcers were recurrent (48%), long lasting (mean, 27 months), and large (mean, 13 cm2). All but one patient had deep venous reflux and/or incompetent perforating veins in addition to trunk varices. Characteristics of patients and ulcers were evenly distributed (exception: more edema in the MCS group; P = .019). Healing within 90 days was observed in 36% with MCS and in 48% with bandages (P = .350). Healing within 180 days was documented in 50% with MCS and in 67% with bandages (P = .210). Time to healing was identical. Pain scored 44 and 46 initially (on a scale in which 100 referred to maximum and 0 to no pain) and decreased within the first week to 20 and 28 in the MCS and bandage groups, respectively (P < .001 vs .010). Quality of life showed no difference between the treatment groups. In both groups, pain at 90 days had decreased by half, independent of completion of healing. Physical, social, and psychic impairment improved significantly in patients with healed ulcers only. CONCLUSION: Our study illustrates the difficulty of bringing large and long-standing venous ulcers to heal. The effect of compression with MCS was not different from that of compression with bandages. Both treatments alleviated pain promptly. Quality of life was improved only in patients whose ulcers had healed.