The Coupled Electronic-Ionic Monte Carlo Simulation Method

Quantum Monte Carlo (QMC) methods such as Variational Monte Carlo, Diffusion Monte Carlo or Path Integral Monte Carlo are the most accurate and general methods for computing total electronic energies. We will review methods we have developed to perform QMC for the electrons coupled to a classical Monte Carlo simulation of the ions. In this method, one estimates the Born-Oppenheimer energy E(Z) where Z represents the ionic degrees of freedom. That estimate of the energy is used in a Metropolis simulation of the ionic degrees of freedom. Important aspects of this method are how to deal with the noise, which QMC method and which trial function to use, how to deal with generalized boundary conditions on the wave function so as to reduce the finite size effects. We discuss some advantages of the CEIMC method concerning how the quantum effects of the ionic degrees of freedom can be included and how the boundary conditions can be integrated over. Using these methods, we have performed simulations of liquid H2 and metallic H on a parallel computer.Comment: 27 pages, 10 figure

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Cellular Levels and Binding of c-di-GMP Control Subcellular Localization and Activity of the Vibrio cholerae Transcriptional Regulator VpsT

The second messenger, cyclic diguanylate (c-di-GMP), regulates diverse cellular processes in bacteria. C-di-GMP is produced by diguanylate cyclases (DGCs), degraded by phosphodiesterases (PDEs), and receptors couple c-di-GMP production to cellular responses. In many bacteria, including Vibrio cholerae, multiple DGCs and PDEs contribute to c-di-GMP signaling, and it is currently unclear whether the compartmentalization of c-di-GMP signaling components is required to mediate c-di-GMP signal transduction. In this study we show that the transcriptional regulator, VpsT, requires c-di-GMP binding for subcellular localization and activity. Only the additive deletion of five DGCs markedly decreases the localization of VpsT, while single deletions of each DGC do not impact VpsT localization. Moreover, mutations in residues required for c-di-GMP binding, c-di-GMP-stabilized dimerization and DNA binding of VpsT abrogate wild type localization and activity. VpsT does not co-localize or interact with DGCs suggesting that c-di-GMP from these DGCs diffuses to VpsT, supporting a model in which c-di-GMP acts at a distance. Furthermore, VpsT localization in a heterologous host, Escherichia coli, requires a catalytically active DGC and is enhanced by the presence of VpsT-target sequences. Our data show that c-di-GMP signaling can be executed through an additive cellular c-di-GMP level from multiple DGCs affecting the localization and activity of a c-di-GMP receptor and furthers our understanding of the mechanisms of second messenger signaling

Concentration Dependent Ion Selectivity in VDAC: A Molecular Dynamics Simulation Study

The voltage-dependent anion channel (VDAC) forms the major pore in the outer mitochondrial membrane. Its high conducting open state features a moderate anion selectivity. There is some evidence indicating that the electrophysiological properties of VDAC vary with the salt concentration. Using a theoretical approach the molecular basis for this concentration dependence was investigated. Molecular dynamics simulations and continuum electrostatic calculations performed on the mouse VDAC1 isoform clearly demonstrate that the distribution of fixed charges in the channel creates an electric field, which determines the anion preference of VDAC at low salt concentration. Increasing the salt concentration in the bulk results in a higher concentration of ions in the VDAC wide pore. This event induces a large electrostatic screening of the charged residues promoting a less anion selective channel. Residues that are responsible for the electrostatic pattern of the channel were identified using the molecular dynamics trajectories. Some of these residues are found to be conserved suggesting that ion permeation between different VDAC species occurs through a common mechanism. This inference is buttressed by electrophysiological experiments performed on bean VDAC32 protein akin to mouse VDAC

Bafilomycin A1 activates respiration of neuronal cells via uncoupling associated with flickering depolarization of mitochondria

Bafilomycin A1 (Baf) induces an elevation of cytosolic Ca2+ and acidification in neuronal cells via inhibition of the V-ATPase. Also, Baf uncouples mitochondria in differentiated PC12 (dPC12), dSH-SY5Y cells and cerebellar granule neurons, and markedly elevates their respiration. This respiratory response in dPC12 is accompanied by morphological changes in the mitochondria and decreases the mitochondrial pH, Ca2+ and ΔΨm. The response to Baf is regulated by cytosolic Ca2+ fluxes from the endoplasmic reticulum. Inhibition of permeability transition pore opening increases the depolarizing effect of Baf on the ΔΨm. Baf induces stochastic flickering of the ΔΨm with a period of 20 ± 10 s. Under conditions of suppressed ATP production by glycolysis, oxidative phosphorylation impaired by Baf does not provide cells with sufficient ATP levels. Cells treated with Baf become more susceptible to excitation with KCl. Such mitochondrial uncoupling may play a role in a number of (patho)physiological conditions induced by Baf

Complex c-di-GMP Signaling Networks Mediate Transition between Virulence Properties and Biofilm Formation in Salmonella enterica Serovar Typhimurium

Upon Salmonella enterica serovar Typhimurium infection of the gut, an early line of defense is the gastrointestinal epithelium which senses the pathogen and intrusion along the epithelial barrier is one of the first events towards disease. Recently, we showed that high intracellular amounts of the secondary messenger c-di-GMP in S. typhimurium inhibited invasion and abolished induction of a pro-inflammatory immune response in the colonic epithelial cell line HT-29 suggesting regulation of transition between biofilm formation and virulence by c-di-GMP in the intestine. Here we show that highly complex c-di-GMP signaling networks consisting of distinct groups of c-di-GMP synthesizing and degrading proteins modulate the virulence phenotypes invasion, IL-8 production and in vivo colonization in the streptomycin-treated mouse model implying a spatial and timely modulation of virulence properties in S. typhimurium by c-di-GMP signaling. Inhibition of the invasion and IL-8 induction phenotype by c-di-GMP (partially) requires the major biofilm activator CsgD and/or BcsA, the synthase for the extracellular matrix component cellulose. Inhibition of the invasion phenotype is associated with inhibition of secretion of the type three secretion system effector protein SipA, which requires c-di-GMP metabolizing proteins, but not their catalytic activity. Our findings show that c-di-GMP signaling is at least equally important in the regulation of Salmonella-host interaction as in the regulation of biofilm formation at ambient temperature

A Phylogenetic Analysis of the Globins in Fungi

BACKGROUND: ALL GLOBINS BELONG TO ONE OF THREE FAMILIES: the F (flavohemoglobin) and S (sensor) families that exhibit the canonical 3/3 α-helical fold, and the T (truncated 3/3 fold) globins characterized by a shortened 2/2 α-helical fold. All eukaryote 3/3 hemoglobins are related to the bacterial single domain F globins. It is known that Fungi contain flavohemoglobins and single domain S globins. Our aims are to provide a census of fungal globins and to examine their relationships to bacterial globins. RESULTS: Examination of 165 genomes revealed that globins are present in >90% of Ascomycota and ∼60% of Basidiomycota genomes. The S globins occur in Blastocladiomycota and Chytridiomycota in addition to the phyla that have FHbs. Unexpectedly, group 1 T globins were found in one Blastocladiomycota and one Chytridiomycota genome. Phylogenetic analyses were carried out on the fungal globins, alone and aligned with representative bacterial globins. The Saccharomycetes and Sordariomycetes with two FHbs form two widely divergent clusters separated by the remaining fungal sequences. One of the Saccharomycete groups represents a new subfamily of FHbs, comprising a previously unknown N-terminal and a FHb missing the C-terminal moiety of its reductase domain. The two Saccharomycete groups also form two clusters in the presence of bacterial FHbs; the surrounding bacterial sequences are dominated by Proteobacteria and Bacilli (Firmicutes). The remaining fungal FHbs cluster with Proteobacteria and Actinobacteria. The Sgbs cluster separately from their bacterial counterparts, except for the intercalation of two Planctomycetes and a Proteobacterium between the Fungi incertae sedis and the Blastocladiomycota and Chytridiomycota. CONCLUSION: Our results are compatible with a model of globin evolution put forward earlier, which proposed that eukaryote F, S and T globins originated via horizontal gene transfer of their bacterial counterparts to the eukaryote ancestor, resulting from the endosymbiotic events responsible for the origin of mitochondria and chloroplasts

Thermal Transport in Micro- and Nanoscale Systems

Small-scale (micro-/nanoscale) heat transfer has broad and exciting range of applications. Heat transfer at small scale quite naturally is influenced – sometimes dramatically – with high surface area-to-volume ratios. This in effect means that heat transfer in small-scale devices and systems is influenced by surface treatment and surface morphology. Importantly, interfacial dynamic effects are at least non-negligible, and there is a strong potential to engineer the performance of such devices using the progress in micro- and nanomanufacturing technologies. With this motivation, the emphasis here is on heat conduction and convection. The chapter starts with a broad introduction to Boltzmann transport equation which captures the physics of small-scale heat transport, while also outlining the differences between small-scale transport and classical macroscale heat transport. Among applications, examples are thermoelectric and thermal interface materials where micro- and nanofabrication have led to impressive figure of merits and thermal management performance. Basic of phonon transport and its manipulation through nanostructuring materials are discussed in detail. Small-scale single-phase convection and the crucial role it has played in developing the thermal management solutions for the next generation of electronics and energy-harvesting devices are discussed as the next topic. Features of microcooling platforms and physics of optimized thermal transport using microchannel manifold heat sinks are discussed in detail along with a discussion of how such systems also facilitate use of low-grade, waste heat from data centers and photovoltaic modules. Phase change process and their control using surface micro-/nanostructure are discussed next. Among the feature considered, the first are microscale heat pipes where capillary effects play an important role. Next the role of nanostructures in controlling nucleation and mobility of the discrete phase in two-phase processes, such as boiling, condensation, and icing is explained in great detail. Special emphasis is placed on the limitations of current surface and device manufacture technologies while also outlining the potential ways to overcome them. Lastly, the chapter is concluded with a summary and perspective on future trends and, more importantly, the opportunities for new research and applications in this exciting field