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.

Computing the phase diagram of binary mixtures: A patchy particle case study

We investigate the phase behaviour of 2D mixtures of bi-functional and three-functional patchy particles and 3D mixtures of bi-functional and tetra-functional patchy particles by means of Monte Carlo simulations and Wertheim theory. We start by computing the critical points of the pure systems and then we investigate how the critical parameters change upon lowering the temperature. We extend the successive umbrella sampling method to mixtures to make it possible to extract information about the phase behaviour of the system at a fixed temperature for the whole range of densities and compositions of interest. (C) 2013 AIP Publishing LLC.

Simulations numériques de la dynamique des protéines : translation de ligands, flexibilité et dynamique des boucles

La flexibilité est une caractéristique intrinsèque des protéines qui doivent, dès le mo- ment de leur synthèse, passer d’un état de chaîne linéaire à un état de structure tridimen- sionnelle repliée et enzymatiquement active. Certaines protéines restent flexibles une fois repliées et subissent des changements de conformation de grande amplitude lors de leur cycle enzymatique. D’autres contiennent des segments si flexibles que leur structure ne peut être résolue par des méthodes expérimentales. Dans cette thèse, nous présentons notre application de méthodes in silico d’analyse de la flexibilité des protéines : • À l’aide des méthodes de dynamique moléculaire dirigée et d’échantillonnage pa- rapluie, nous avons caractérisé les trajectoires de liaison de l’inhibiteur Z-pro- prolinal à la protéine Prolyl oligopeptidase et identifié la trajectoire la plus pro- bable. Nos simulations ont aussi identifié un mode probable de recrutement des ligands utilisant une boucle flexible de 19 acides aminés à l’interface des deux domaines de la protéine. • En utilisant les méthodes de dynamique moléculaire traditionnelle et dirigée, nous avons examiné la stabilité de la protéine SAV1866 dans sa forme fermée insérée dans une membrane lipidique et étudié un des modes d’ouverture possibles par la séparation de ses domaines liant le nucléotide. • Nous avons adapté auproblème de la prédiction de la structure des longues boucles flexibles la méthode d’activation et de relaxation ART-nouveau précédemment uti- lisée dans l’étude du repliement et de l’agrégation de protéines. Appliqué au replie- ment de boucles de 8 à 20 acides aminés, la méthode démontre une dépendance quadratique du temps d’exécution sur la longueur des boucles, rendant possible l’étude de boucles encore plus longues.

Monte-Carlo-Simulationen zum Phasen- und Keimbildungsverhalten von Polymerlösungen

Das Phasenverhalten und die Grenzflächeneigenschaften vonPolymeren in superkritischer Lösung werden anhand einesvergröberten Kugel-Feder-Modells für das ReferenzsystemHexadekan-CO2 untersucht. Zur Bestimmung der Parameter imPotential setzt man die kritischen Punkte von Simulation undExperiment gleich. Wechselwirkungen zwischen beidenKomponenten werden durch eine modifizierteLorentz-Berthelot-Regel modelliert. Die Übereinstimmung mitden Experimenten ist sehr gut - insbesondere kann dasPhasendiagramm des Mischsystems inklusive kritischer Linienreproduziert werden. Ein Vergleich mit numerischenStörungsrechnungen (TPT1) liefert eine qualitativeÜbereinstimmung und Hinweise zur Verbesserung derverwendeten Zustandsgleichung. Aufbauend auf diesen Betrachtungen werden die Frühstadiender Keimbildung untersucht. Für das Lennard-Jones-Systemwird zum ersten Mal der Übergang vom homogenen Gas zu einemeinzelnen Tropfen im endlichen Volumen direkt nachgewiesenund quantifiziert. Die freie Energie von kleinen Clusternwird mit einem einfachen, klassischen Nukleationsmodellbestimmt und nach oben abgeschätzt. Die vorgestellten Untersuchungen wurden durch eineWeiterentwicklung des Umbrella-Sampling-Algorithmusermöglicht. Hierbei wird die Simulation in mehrereSimulationsfenster unterteilt, die nacheinander abgearbeitetwerden. Die Methode erlaubt eine Bestimmung derFreien-Energie-Landschaft an einer beliebigen Stelle desPhasendiagramms. Der Fehler ist kontrollierbar undunabhängig von der Art der Unterteilung.

Mixtures of colloidal rods and spheres in bulk and in confinement

When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. Monte Carlo simulations were performed to study the phase diagram of such rod-polymer mixtures. The colloidal rods were modelled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while overlap of polymers and rods is forbidden. In this thesis the emphasis was on the depletion effects caused by the addition of spheres on the isotropic phase of rod-like particles. Although most of the present experimental studies consider systems close to or beyond the isotropic-nematic transition, the isotropic phase with depletion interactions turns out to be a not less interesting topic. First, the percolation problem was studied in canonical simulations of a system of hard rods and soft spheres, where the amount of depletant was kept low to prevent phase separation of the mixture. The lowering of the percolation threshold seen in experiment is confirmed to be due to the depletion interactions. The local changes in the structure of the fluid of rods, which were measured in the simulations, indicated that the depletion forces enhance local alignment and aggregation of the rods. Then, the phase diagram of isotropic-isotropic demixing of short spherocylinders was calculated using grand canonical ensemble simulations with successive umbrella sampling. Finite size scaling analysis allowed to estimate the location of the critical point. Also, estimates for the interfacial tension between the coexisting isotropic phases and analyses of its power-law behaviour on approach of the critical point are presented. The obtained phase diagram was compared to the predictions of the free volume theory. After an analysis of the bulk, the phase behaviour in confinement was studied. The critical point of gas-liquid demixing is shifted to higher concentrations of rods and smaller concentrations of spheres due to the formation of an orientationally ordered surface film. If the separation between the walls becomes very small, the critical point is shifted back to smaller concentrations of rods because the surface film breaks up. A method to calculate the contact angle of the liquid-gas interface with the wall is introduced and the wetting behaviour on the approach to the critical point is analysed.

Parallelization of the algorithm WHAM with NVIDIA CUDA.

The aim of my thesis is to parallelize the Weighting Histogram Analysis Method (WHAM), which is a popular algorithm used to calculate the Free Energy of a molucular system in Molecular Dynamics simulations. WHAM works in post processing in cooperation with another algorithm called Umbrella Sampling. Umbrella Sampling has the purpose to add a biasing in the potential energy of the system in order to force the system to sample a specific region in the configurational space. Several N independent simulations are performed in order to sample all the region of interest. Subsequently, the WHAM algorithm is used to estimate the original system energy starting from the N atomic trajectories. The parallelization of WHAM has been performed through CUDA, a language that allows to work in GPUs of NVIDIA graphic cards, which have a parallel achitecture. The parallel implementation may sensibly speed up the WHAM execution compared to previous serial CPU imlementations. However, the WHAM CPU code presents some temporal criticalities to very high numbers of interactions. The algorithm has been written in C++ and executed in UNIX systems provided with NVIDIA graphic cards. The results were satisfying obtaining an increase of performances when the model was executed on graphics cards with compute capability greater. Nonetheless, the GPUs used to test the algorithm is quite old and not designated for scientific calculations. It is likely that a further performance increase will be obtained if the algorithm would be executed in clusters of GPU at high level of computational efficiency. The thesis is organized in the following way: I will first describe the mathematical formulation of Umbrella Sampling and WHAM algorithm with their apllications in the study of ionic channels and in Molecular Docking (Chapter 1); then, I will present the CUDA architectures used to implement the model (Chapter 2); and finally, the results obtained on model systems will be presented (Chapter 3).

Equilibrium fluid-solid coexistence of hard spheres

We present a tethered Monte Carlo simulation of the crystallization of hard spheres. Our method boosts the traditional umbrella sampling to the point of making practical the study of constrained Gibbs’ free energies depending on several crystalline order parameters. We obtain high-accuracy estimates of the fluid-crystal coexistence pressure for up to 2916 particles (enough to accommodate fluid-solid interfaces). We are able to extrapolate to infinite volume the coexistence pressure [p_(co) = 11.5727(10)k_(B)T/σ^(3)] and the interfacial free energy [γ_({100}) = 0.636(11)k_(B)T/σ^(2)].

Topics in equilibrium and nonequilibrium thermodynamics: computing crystalline free energies and engineering Maxwell’s demon.

This dissertation covers two separate topics in statistical physics. The first part of the dissertation focuses on computational methods of obtaining the free energies (or partition functions) of crystalline solids. We describe a method to compute the Helmholtz free energy of a crystalline solid by direct evaluation of the partition function. In the many-dimensional conformation space of all possible arrangements of N particles inside a periodic box, the energy landscape consists of localized islands corresponding to different solid phases. Calculating the partition function for a specific phase involves integrating over the corresponding island. Introducing a natural order parameter that quantifies the net displacement of particles from lattices sites, we write the partition function in terms of a one-dimensional integral along the order parameter, and evaluate this integral using umbrella sampling. We validate the method by computing free energies of both face-centered cubic (FCC) and hexagonal close-packed (HCP) hard sphere crystals with a precision of $10^{-5}k_BT$ per particle. In developing the numerical method, we find several scaling properties of crystalline solids in the thermodynamic limit. Using these scaling properties, we derive an explicit asymptotic formula for the free energy per particle in the thermodynamic limit. In addition, we describe several changes of coordinates that can be used to separate internal degrees of freedom from external, translational degrees of freedom. The second part of the dissertation focuses on engineering idealized physical devices that work as Maxwell's demon. We describe two autonomous mechanical devices that extract energy from a single heat bath and convert it into work, while writing information onto memory registers. Additionally, both devices can operate as Landauer's eraser, namely they can erase information from a memory register, while energy is dissipated into the heat bath. The phase diagrams and the efficiencies of the two models are solved and analyzed. These two models provide concrete physical illustrations of the thermodynamic consequences of information processing.

Temporal variation of the phytoplankton community at short sampling intervals in the Mundaú reservoir, Northeastern Brazil

The aim of this study was to determine how abiotic factors drive the phytoplankton community in a water supply reservoir within short sampling intervals. Samples were collected at the subsurface (0.1 m) and bottom of limnetic (8 m) and littoral (2 m) zones in both the dry and rainy seasons. The following abiotic variables were analyzed: water temperature, dissolved oxygen, electrical conductivity, total dissolved solids, turbidity, pH, total nitrogen, nitrite, nitrate, total phosphorus, total dissolved phosphorus and orthophosphate. Phytoplankton biomass was determined from biovolume values. The role abiotic variables play in the dynamics of phytoplankton species was determined by means of Canonical Correspondence Analysis. Algae biomass ranged from 1.17×10(4) to 9.21×10(4) µg.L-1; cyanobacteria had biomass values ranging from 1.07×10(4) to 8.21×10(4) µg.L-1. High availability of phosphorous, nitrogen limitation, alkaline pH and thermal stability all favored cyanobacteria blooms, particularly during the dry season. Temperature, pH, total phosphorous and turbidity were key factors in characterizing the phytoplankton community between sampling times and stations. Of the species studied, Cylindrospermopsis raciborskii populations were dominant in the phytoplankton in both the dry and rainy seasons. We conclude that the phytoplankton was strongly influenced by abiotic variables, particularly in relation to seasonal distribution patterns.

Precision of distances and ordering of microsatellite markers in consensus linkage maps of chromosomes 1, 3 and 4 from two reciprocal chicken populations using bootstrap sampling

Some factors complicate comparisons between linkage maps from different studies. This problem can be resolved if measures of precision, such as confidence intervals and frequency distributions, are associated with markers. We examined the precision of distances and ordering of microsatellite markers in the consensus linkage maps of chromosomes 1, 3 and 4 from two F 2 reciprocal Brazilian chicken populations, using bootstrap sampling. Single and consensus maps were constructed. The consensus map was compared with the International Consensus Linkage Map and with the whole genome sequence. Some loci showed segregation distortion and missing data, but this did not affect the analyses negatively. Several inversions and position shifts were detected, based on 95% confidence intervals and frequency distributions of loci. Some discrepancies in distances between loci and in ordering were due to chance, whereas others could be attributed to other effects, including reciprocal crosses, sampling error of the founder animals from the two populations, F(2) population structure, number of and distance between microsatellite markers, number of informative meioses, loci segregation patterns, and sex. In the Brazilian consensus GGA1, locus LEI1038 was in a position closer to the true genome sequence than in the International Consensus Map, whereas for GGA3 and GGA4, no such differences were found. Extending these analyses to the remaining chromosomes should facilitate comparisons and the integration of several available genetic maps, allowing meta-analyses for map construction and quantitative trait loci (QTL) mapping. The precision of the estimates of QTL positions and their effects would be increased with such information.

Phylogenetic relationships within the speciose family Characidae (Teleostei: Ostariophysi: Characiformes) based on multilocus analysis and extensive ingroup sampling

Background: With nearly 1,100 species, the fish family Characidae represents more than half of the species of Characiformes, and is a key component of Neotropical freshwater ecosystems. The composition, phylogeny, and classification of Characidae is currently uncertain, despite significant efforts based on analysis of morphological and molecular data. No consensus about the monophyly of this group or its position within the order Characiformes has been reached, challenged by the fact that many key studies to date have non-overlapping taxonomic representation and focus only on subsets of this diversity. Results: In the present study we propose a new definition of the family Characidae and a hypothesis of relationships for the Characiformes based on phylogenetic analysis of DNA sequences of two mitochondrial and three nuclear genes (4,680 base pairs). The sequences were obtained from 211 samples representing 166 genera distributed among all 18 recognized families in the order Characiformes, all 14 recognized subfamilies in the Characidae, plus 56 of the genera so far considered incertae sedis in the Characidae. The phylogeny obtained is robust, with most lineages significantly supported by posterior probabilities in Bayesian analysis, and high bootstrap values from maximum likelihood and parsimony analyses. Conclusion: A monophyletic assemblage strongly supported in all our phylogenetic analysis is herein defined as the Characidae and includes the characiform species lacking a supraorbital bone and with a derived position of the emergence of the hyoid artery from the anterior ceratohyal. To recognize this and several other monophyletic groups within characiforms we propose changes in the limits of several families to facilitate future studies in the Characiformes and particularly the Characidae. This work presents a new phylogenetic framework for a speciose and morphologically diverse group of freshwater fishes of significant ecological and evolutionary importance across the Neotropics and portions of Africa.

Comparative analysis of two sampling techniques for pollen gathered by Nannotrigona testaceicornis Lepeletier (Apidae, Meliponini)

Pollen counts from samples taken from storage pots throughout one year (from October to September) were adjusted by Tasei's volumetric correction coefficient for the determination of pollen sources exploited by two colonies of Nannotrigona testaceicornis in Sao Paulo, Brazil. The results obtained by this sampling technique for seven months (December to June) were compared with those from corbicula load samples taken within the same period. This species visited a large variety of plant species, but few of them were frequently used. As a rule, pollen sources that appeared at frequencies greater than 1% were found with both sampling methods and significant positive correlations (Spearman correlation coefficient) were found between their values. The pollen load sample data showed that N. testaceicornis gathered pollen throughout the external activity period.

Method for cadmium and lead longitudinal profiles determination in hair by solid sampling graphite furnace atomic absorption spectrometry

This paper describes methods for the direct determination of Cd and Pb in hair segments (c.a. 5 mm similar to 80 mu g) by solid sampling graphite furnace atomic absorption spectrometry, becoming possible longitudinal profiles in a single strand of hair. To distinguish endogenous and exogenous content. strands of hair were washed by using two different procedures: IAEA protocol (acetone + water + acetone) and the combination of IAEA protocol with HCl washing (acetone + water + acetone + 0.1 mol l(-1) HCl). The concentration of Cd and Pb increased from the root Until the tip of hair washed according to IAEA protocol. However, when the strand of hair was washed using the combination of IAEA protocol and 0.1 mol l(-1) HCl, Cd concentrations decreased in all segments, and Pb concentrations decreased drastically near to the root (5 to 12 mm) and was systematically higher ill the end. The proposed method showed to be useful to assess the temporal variation to Cd and Pb exposure and call be Used for toxicological and environmental investigations. The limits of detection were 2.8 ng g(-1) for Cd and 40 ng g(-1) for Pb. The characteristic masses based oil integrated absorbance were 2.4 pg for Cd and 22 pg for Pb.

Simultaneous determination of Cr, Fe, Ni and V in crude oil by emulsion sampling graphite furnace atomic absorption spectrometry

In this work a simple and reliable method for the simultaneous determination of Cr, Fe, Ni and V in crude oil, using emulsion sampling graphite furnace atomic absorption spectrometry is proposed. Under the best conditions, sample masses around 50 mg were weighed in polypropylene tubes and emulsified in a mixture of 0.5% (v v(-1)) hexane + 6% (m v(-1)) Triton X-100 (R). Considering the compromised conditions, the pyrolysis an atomization temperatures for the simultaneous determination of Cr, Fe, Ni and V were 1400 degrees C and 2500 degrees C, respectively. Aliquots of 20 mu L of reference solution and sample emulsion were co-injected into the graphite tube with 10 mu L of 1.0 g L(-1) Mg(NO(3))(2) as chemical modifier. The detection limits (n = 10, 3 sigma) and characteristic masses were, respectively: 0.07 mu g g(-1) and 19 pg for Cr; 2.15 mu g g(-1) and 31 pg for Fe; 1.25 mu g g(-1) and 44 pg for Ni; and 1.15 mu g g(-1) and 149 pg for V. The reliability of the proposed method was checked by fuel oil Standard Reference Material (SRMTriton X-100 (R) 1634c - NIST) analysis. The concentrations found presented no statistical differences compared to the certified values at 95% confidence level.

A sampling device to quantify offspring release of sessile marine invertebrates

Quantifying the rate of propagule release is of most importance to estimate reproductive output of natural populations, but simple methods to obtain such data are seldom reported. We designed and tested an inexpensive apparatus capable of reliably measure the release of gametes, eggs or larvae of sessile marine invertebrates in vertical walls. A population of the acom barnacle Chthamalus bisinuatus was sampled with this trap over 68d to obtain a time series of naupliar release. An apparent semilunar trend is shown, indicating the effectiveness of this sampling method.