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.

Proteins in silico-modeling and sampling

Proteins are linear chain molecules made out of amino acids. Only when they fold to their native states, they become functional. This dissertation aims to model the solvent (environment) effect and to develop & implement enhanced sampling methods that enable a reliable study of the protein folding problem in silico. We have developed an enhanced solvation model based on the solution to the Poisson-Boltzmann equation in order to describe the solvent effect. Following the quantum mechanical Polarizable Continuum Model (PCM), we decomposed net solvation free energy into three physical terms– Polarization, Dispersion and Cavitation. All the terms were implemented, analyzed and parametrized individually to obtain a high level of accuracy. In order to describe the thermodynamics of proteins, their conformational space needs to be sampled thoroughly. Simulations of proteins are hampered by slow relaxation due to their rugged free-energy landscape, with the barriers between minima being higher than the thermal energy at physiological temperatures. In order to overcome this problem a number of approaches have been proposed of which replica exchange method (REM) is the most popular. In this dissertation we describe a new variant of canonical replica exchange method in the context of molecular dynamic simulation. The advantage of this new method is the easily tunable high acceptance rate for the replica exchange. We call our method Microcanonical Replica Exchange Molecular Dynamic (MREMD). We have described the theoretical frame work, comment on its actual implementation, and its application to Trp-cage mini-protein in implicit solvent. We have been able to correctly predict the folding thermodynamics of this protein using our approach.

PLUMED 2: New feathers for an old bird

Enhancing sampling and analyzing simulations are central issues in molecular simulation. Recently, we introduced PLUMED, an open-source plug-in that provides some of the most popular molecular dynamics (MD) codes with implementations of a variety of different enhanced sampling algorithms and collective variables (CVs). The rapid changes in this field, in particular new directions in enhanced sampling and dimensionality reduction together with new hardware, require a code that is more flexible and more efficient. We therefore present PLUMED 2 here a,complete rewrite of the code in an object-oriented programming language (C++). This new version introduces greater flexibility and greater modularity, which both extends its core capabilities and makes it far easier to add new methods and CVs. It also has a simpler interface with the MD engines and provides a single software library containing both tools and core facilities. Ultimately, the new code better serves the ever-growing community of users and contributors in coping with the new challenges arising in the field.

Program summary

Program title: PLUMED 2

Catalogue identifier: AEEE_v2_0

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEEE_v2_0.html

Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland

Licensing provisions: Yes

No. of lines in distributed program, including test data, etc.: 700646

No. of bytes in distributed program, including test data, etc.: 6618136

Distribution format: tar.gz

Programming language: ANSI-C++.

Computer: Any computer capable of running an executable produced by a C++ compiler.

Operating system: Linux operating system, Unix OSs.

Has the code been vectorized or parallelized?: Yes, parallelized using MPI.

RAM: Depends on the number of atoms, the method chosen and the collective variables used.

Classification: 3, 7.7, 23. Catalogue identifier of previous version: AEEE_v1_0.

Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 1961.

External routines: GNU libmatheval, Lapack, Bias, MPI. (C) 2013 Elsevier B.V. All rights reserved.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

Development of models and methods to simulate peptide-assisted nucleation and growth of calcium-minerals

In den vergangenen Jahren wurden einige bislang unbekannte Phänomene experimentell beobachtet, wie etwa die Existenz unterschiedlicher Prä-Nukleations-Strukturen. Diese haben zu einem neuen Verständnis von Prozessen, die auf molekularer Ebene während der Nukleation und dem Wachstum von Kristallen auftreten, beigetragen. Die Auswirkungen solcher Prä-Nukleations-Strukturen auf den Prozess der Biomineralisation sind noch nicht hinreichend verstanden. Die Mechanismen, mittels derer biomolekulare Modifikatoren, wie Peptide, mit Prä-Nukleations-Strukturen interagieren und somit den Nukleationsprozess von Mineralen beeinflussen könnten, sind vielfältig. Molekulare Simulationen sind zur Analyse der Formation von Prä-Nukleations-Strukturen in Anwesenheit von Modifikatoren gut geeignet. Die vorliegende Arbeit beschreibt einen Ansatz zur Analyse der Interaktion von Peptiden mit den in Lösung befindlichen Bestandteilen der entstehenden Kristalle mit Hilfe von Molekular-Dynamik Simulationen.rnUm informative Simulationen zu ermöglichen, wurde in einem ersten Schritt die Qualität bestehender Kraftfelder im Hinblick auf die Beschreibung von mit Calciumionen interagierenden Oligoglutamaten in wässrigen Lösungen untersucht. Es zeigte sich, dass große Unstimmigkeiten zwischen etablierten Kraftfeldern bestehen, und dass keines der untersuchten Kraftfelder eine realistische Beschreibung der Ionen-Paarung dieser komplexen Ionen widerspiegelte. Daher wurde eine Strategie zur Optimierung bestehender biomolekularer Kraftfelder in dieser Hinsicht entwickelt. Relativ geringe Veränderungen der auf die Ionen–Peptid van-der-Waals-Wechselwirkungen bezogenen Parameter reichten aus, um ein verlässliches Modell für das untersuchte System zu erzielen. rnDas umfassende Sampling des Phasenraumes der Systeme stellt aufgrund der zahlreichen Freiheitsgrade und der starken Interaktionen zwischen Calciumionen und Glutamat in Lösung eine besondere Herausforderung dar. Daher wurde die Methode der Biasing Potential Replica Exchange Molekular-Dynamik Simulationen im Hinblick auf das Sampling von Oligoglutamaten justiert und es erfolgte die Simulation von Peptiden verschiedener Kettenlängen in Anwesenheit von Calciumionen. Mit Hilfe der Sketch-Map Analyse konnten im Rahmen der Simulationen zahlreiche stabile Ionen-Peptid-Komplexe identifiziert werden, welche die Formation von Prä-Nukleations-Strukturen beeinflussen könnten. Abhängig von der Kettenlänge des Peptids weisen diese Komplexe charakteristische Abstände zwischen den Calciumionen auf. Diese ähneln einigen Abständen zwischen den Calciumionen in jenen Phasen von Calcium-Oxalat Kristallen, die in Anwesenheit von Oligoglutamaten gewachsen sind. Die Analogie der Abstände zwischen Calciumionen in gelösten Ionen-Peptid-Komplexen und in Calcium-Oxalat Kristallen könnte auf die Bedeutung von Ionen-Peptid-Komplexen im Prozess der Nukleation und des Wachstums von Biomineralen hindeuten und stellt einen möglichen Erklärungsansatz für die Fähigkeit von Oligoglutamaten zur Beeinflussung der Phase des sich formierenden Kristalls dar, die experimentell beobachtet wurde.

Structural flexibility and oxygen diffusion pathways in monomeric fluorescent proteins

Fluorescent proteins are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. We used molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variants, and its GFP homolog citrine. We employed implicit ligand sampling and locally enhanced sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. The pathway contains several oxygen hosting pockets, which were identified by the amino acid residues that form the pocket. We calculated the free-energy of an oxygen molecule at points along the path. We also investigated an RFP variant known to be significantly less photostable than mCherry and find much easier oxygen access in this variant. We showed that oxygen pathways can be blocked or altered, and barrel fluctuations can be reduced by strategic amino acid substitutions. The results provide a better understanding of the mechanism of molecular oxygen access into the fully folded mCherry protein barrel and provide insight into the photobleaching process in these proteins.

Enhanced Sensitivity by Nonuniform Sampling Enables Multidimensional MAS NMR Spectroscopy of Protein Assemblies

We report dramatic sensitivity enhancements in multidimensional MAS NMR spectra by the use of nonuniform sampling (NUS) and introduce maximum entropy interpolation (MINT) processing that assures the linearity between the time and frequency domains of the NUS acquired data sets. A systematic analysis of sensitivity and resolution in 2D and 3D NUS spectra reveals that with NUS, at least 1.5- to 2-fold sensitivity enhancement can be attained in each indirect dimension without compromising the spectral resolution. These enhancements are similar to or higher than those attained by the newest-generation commercial cryogenic probes. We explore the benefits of this NUS/MaxEnt approach in proteins and protein assemblies using 1-73-(U-C-13,N-15)/74-108-(U-N-15) Escherichia coil thioredoxin reassembly. We demonstrate that in thioredoxin reassembly, NUS permits acquisition of high-quality 3D-NCACX spectra, which are inaccessible with conventional sampling due to prohibitively long experiment times. Of critical importance, issues that hinder NUS-based SNR enhancement in 3D-NMR of liquids are mitigated in the study of solid samples in which theoretical enhancements on the order of 3-4 fold are accessible by compounding the NUS-based SNR enhancement of each indirect dimension. NUS/MINT is anticipated to be widely applicable and advantageous for multidimensional heteronuclear MAS NMR spectroscopy of proteins, protein assemblies, and other biological systems.

Enhanced conformational sampling in Monte Carlo simulations of proteins: application to a constrained peptide.

A Monte Carlo simulation method for globular proteins, called extended-scaled-collective-variable (ESCV) Monte Carlo, is proposed. This method combines two Monte Carlo algorithms known as entropy-sampling and scaled-collective-variable algorithms. Entropy-sampling Monte Carlo is able to sample a large configurational space even in a disordered system that has a large number of potential barriers. In contrast, scaled-collective-variable Monte Carlo provides an efficient sampling for a system whose dynamics is highly cooperative. Because a globular protein is a disordered system whose dynamics is characterized by collective motions, a combination of these two algorithms could provide an optimal Monte Carlo simulation for a globular protein. As a test case, we have carried out an ESCV Monte Carlo simulation for a cell adhesive Arg-Gly-Asp-containing peptide, Lys-Arg-Cys-Arg-Gly-Asp-Cys-Met-Asp, and determined the conformational distribution at 300 K. The peptide contains a disulfide bridge between the two cysteine residues. This bond mimics the strong geometrical constraints that result from a protein's globular nature and give rise to highly cooperative dynamics. Computation results show that the ESCV Monte Carlo was not trapped at any local minimum and that the canonical distribution was correctly determined.

Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification

In recent years, concern has arisen over the effects of increasing carbon dioxide (CO2) in the earth's atmosphere due to the burning of fossil fuels. One way to mitigate increase in atmospheric CO2 concentration and climate change is carbon sequestration to forest vegeta-tion through photosynthesis. Comparable regional scale estimates for the carbon balance of forests are therefore needed for scientific and political purposes. The aim of the present dissertation was to improve methods for quantifying and verifying inventory-based carbon pool estimates of the boreal forests in the mineral soils. Ongoing forest inventories provide a data based on statistically sounded sampling for estimating the level of carbon stocks and stock changes, but improved modelling tools and comparison of methods are still needed. In this dissertation, the entire inventory-based large-scale forest carbon stock assessment method was presented together with some separate methods for enhancing and comparing it. The enhancement methods presented here include ways to quantify the biomass of understorey vegetation as well as to estimate the litter production of needles and branches. In addition, the optical remote sensing method illustrated in this dis-sertation can be used to compare with independent data. The forest inventory-based large-scale carbon stock assessment method demonstrated here provided reliable carbon estimates when compared with independent data. Future ac-tivity to improve the accuracy of this method could consist of reducing the uncertainties regarding belowground biomass and litter production as well as the soil compartment. The methods developed will serve the needs for UNFCCC reporting and the reporting under the Kyoto Protocol. This method is principally intended for analysts or planners interested in quantifying carbon over extensive forest areas.

Remote sensing of leaf area index: enhanced retrieval from close-range and remotely sensed optical observations

A wide range of models used in agriculture, ecology, carbon cycling, climate and other related studies require information on the amount of leaf material present in a given environment to correctly represent radiation, heat, momentum, water, and various gas exchanges with the overlying atmosphere or the underlying soil. Leaf area index (LAI) thus often features as a critical land surface variable in parameterisations of global and regional climate models, e.g., radiation uptake, precipitation interception, energy conversion, gas exchange and momentum, as all areas are substantially determined by the vegetation surface. Optical wavelengths of remote sensing are the common electromagnetic regions used for LAI estimations and generally for vegetation studies. The main purpose of this dissertation was to enhance the determination of LAI using close-range remote sensing (hemispherical photography), airborne remote sensing (high resolution colour and colour infrared imagery), and satellite remote sensing (high resolution SPOT 5 HRG imagery) optical observations. The commonly used light extinction models are applied at all levels of optical observations. For the sake of comparative analysis, LAI was further determined using statistical relationships between spectral vegetation index (SVI) and ground based LAI. The study areas of this dissertation focus on two regions, one located in Taita Hills, South-East Kenya characterised by tropical cloud forest and exotic plantations, and the other in Gatineau Park, Southern Quebec, Canada dominated by temperate hardwood forest. The sampling procedure of sky map of gap fraction and size from hemispherical photographs was proven to be one of the most crucial steps in the accurate determination of LAI. LAI and clumping index estimates were significantly affected by the variation of the size of sky segments for given zenith angle ranges. On sloping ground, gap fraction and size distributions present strong upslope/downslope asymmetry of foliage elements, and thus the correction and the sensitivity analysis for both LAI and clumping index computations were demonstrated. Several SVIs can be used for LAI mapping using empirical regression analysis provided that the sensitivities of SVIs at varying ranges of LAI are large enough. Large scale LAI inversion algorithms were demonstrated and were proven to be a considerably efficient alternative approach for LAI mapping. LAI can be estimated nonparametrically from the information contained solely in the remotely sensed dataset given that the upper-end (saturated SVI) value is accurately determined. However, further study is still required to devise a methodology as well as instrumentation to retrieve on-ground green leaf area index . Subsequently, the large scale LAI inversion algorithms presented in this work can be precisely validated. Finally, based on literature review and this dissertation, potential future research prospects and directions were recommended.

Microporation Techniques for Enhanced Delivery of Therapeutic Agents

Perhaps the greatest barrier to development of the field of transmembrane drug delivery is that only a limited number of drugs are amenable to administration by this route. The highly lipophilic nature and barrier function of the uppermost layer of the skin, the stratum corneum, for example, restricts the permeation of hydrophilic, high molecular weight and charged compounds into the systemic circulation. Other membranes in the human body can also present significant barriers to drug permeation. In order to successfully deliver hydrophilic drugs, and macromolecular agents of interest, including peptides, DNA and small interfering RNA, many research groups and pharmaceutical companies Worldwide are focusing on the use of microporation methods and devices. Whilst there are a variety of microporation techniques, including the use of laser, thermal ablation, electroporation, radiofrequency, ultrasound, high pressure jets, and microneedle technology, they share the common goal of enhancing the permeability of a biological membrane through the creation of transient aqueous transport pathways of micron dimensions across that membrane. Once created, these micropores are orders of magnitude larger than molecular dimensions and, therefore, should readily permit the transport of hydrophilic macromolecules. Additionally, microporation devices also enable minimally-invasive sampling and monitoring of biological fluids. This review deals with the innovations relating to microporation-based methods and devices for drug delivery and minimally invasive monitoring, as disclosed in recent patent literature. © 2010 Bentham Science Publishers Ltd.

Enhanced proton flux in the MeV range by defocused laser irradiation

Thin Al foils (50 nm and 6 mu m) were irradiated at intensities of up to 2x10(19) W cm(-2) using high contrast (10(8)) laser pulses. Ion emission from the rear of the targets was measured using a scintillator-based Thomson parabola and beam sampling 'footprint' monitor. The variation of the ion spectra and beam profile with focal spot size was systematically studied. The results show that while the maximum proton energy is achieved around tight focus for both target thicknesses, as the spot size increases the ion flux at lower energies is seen to peak at significantly increased spot sizes. Measurements of the proton footprint, however, show that the off-axis proton flux is highest at tight focus, indicating that a previously identified proton deflection mechanism may alter the on-axis spectrum. One-dimensional particle-in-cell modelling of the experiment supports our hypothesis that the observed change in spectra with focal spot size is due to the competition of two effects: decrease in laser intensity and an increase in proton emission area.

Pilot study of the utility and acceptability of tampon sampling for the diagnosis of Neisseria gonorrhoeae and Chlamydia trachomatis infections by duplex realtime polymerase chain reaction in United Kingdom sex workers

We aimed to evaluate the acceptability of self-collected tampon samples for the screening of female sex workers for sexually transmitted infections. We recruited 65 sex workers, and 63 agreed to provide tampon samples. The tampon samples were processed by realtime polymerase chain reaction (PCR) targeting Neisseria gonorrhoeae and Chlamydia trachomatis. Urethral and endocervical swabs were also obtained from 61 of 63 participants and tested using culture (N. gonorrhoeae) and the BD ProbeTec strand displacement amplification (SDA) (C. trachomatis) assay. Tampon sampling was preferred by 95% of the women and all favoured being tested away from genitourinary medicine clinics; the most common reasons cited were avoidance of embarrassment (40%) and convenience (30%). Besides near-universal acceptability of tampon sampling, the tampon sampling-PCR approach described in this study appeared to have enhanced sensitivity compared with conventional testing, suggesting the possibility of a residual hidden burden of N. gonorrhoeae and/or C. trachomatis genital infections in UK female sex workers.

Radiation dose rates and exposure associated to Technologically Enhanced Naturally Occurring Radioactive Materials (TENORM) at locations of a Portuguese coal-fired power plant

Coal contains trace quantities of natural radionuclides such as Th-232, U-235, U-238, as well as their radioactive decay products and 40K. These radionuclides can be released as fly ash in atmospheric emissions from coal-fired power plants, dispersed into the environment and deposited on the surrounding top soils. Therefore, the natural radiation background level is enhanced and consequently increase the total dose for the nearby population. A radiation monitoring programme was used to assess the external dose contribution to the natural radiation background, potentially resulting from the dispersion of coal ash in past atmospheric emissions. Radiation measurements were carried out by gamma spectrometry in the vicinity of a Portuguese coal-fired power plant. The radiation monitoring was achieved both on and off site, being the boundary delimited by a 20 km circle centered in the stacks of the coal plant. The measured radionuclides concentrations for the uranium and thorium series ranged from 7.7 to 41.3 Bq/kg for Ra-226 and from 4.7 to 71.6 Bq/kg for Th-232, while K-40 concentrations ranged from 62.3 to 795.1 Bq/kg. The highest values were registered near the power plant and at distances between 6 and 20 km from the stacks, mainly in the prevailing wind direction. The absorbed dose rates were calculated for each sampling location: 13.97-84.00 ηGy/h, while measurements from previous studies carried out in 1993 registered values in the range of 16.6-77.6 ηGy/h. The highest values were registered at locations in the prevailing wind direction (NW-SE). This study has been primarily done to assess the radiation dose rates and exposure to the nearby population in the surroundings of a coal-fired power plant. The results suggest an enhancement or at least an influence in the background radiation due to the coal plant past activities.