Search NASA⌕ Search

SEARCH · Search NASA

Results for “Pore Network Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

In-situ Imaging of Pyrolyzing Aerospace Materials

Tracking morphological changes of materials during heating is crucial to understand its response in fire protection, biofuel production, thermal protection systems (TPS) for hypersonic flight. As materials are heated, they undergo physical and chemical changes due to water loss, stretching or shrinking, pyrolysis and chemical reactions in the ambient environment. The effects of these changes can have a profound impact on the material’s performance, indicated by changes in on the porosity and volume. While materials such as wood shrink as they pyrolyze and lose mass, others swell due to their inherent characteristics when exposed to heat [1]. This study focuses on experiments conducted at the Advanced Light Source (ALS) beamline 8.3.2, where in situ micro-computed tomography (µ-CT) is performed on materials as they are being pyrolyzed. Through in situ µ-CT, the change in total volume and porosity can be obtained in real-time, allowing for better understanding of the underlying thermophysical and chemical processes as a function of temperature. This study also focuses on the implementation of the Porous Microstructure Analysis software (PuMA) [2] to obtain thermal conductivity, permeability, and other properties of the material from the 3D tomographies. The information gained from these tomographies will supplement microscale model development of material morphological change and will aid macroscale modeling for high-temperature applications. For this study, Room Temperature Vulcanizing silicone (RTV) [3-5] is heated from room temperature to 1000°C using an infrared lamp heating system, and tomographies are continuously collected as the sample is heated. The tomographies are then segmented to obtain solid and void phases, from which estimates of pore size, porosity and total volume are extracted as a function of temperature. PuMA is deployed on the segmented tomographies to obtain thermal conductivity, permeability, and other properties as a function of temperature. Preliminary results show that RTV first intumesces (swells) as pyrolysis begins, due to build-up of pyrolysis gases in closed pores, and then shrinks significantly as more open pores are formed and the pyrolysis gases outgas. Pore network visualization of the tomographies using OpenPNM [6] showed the increase in pore connectivity with increase in temperature. Future work will focus on using PuMA to obtain macroscopic properties of RTV as a function of temperature.

Tomography↗
DOE OSTI2020

Beam Hardening Artifact Reduction in X-Ray CT Reconstruction of 3D Printed Metal Parts Leveraging Deep Learning and CAD Models

Nondestructive evaluation (NDE) of additively manufactured (AM) parts is important for understanding the impacts of various process parameters and qualifying the built part. X-ray computed tomography (XCT) has played a critical role in rapid NDE and characterization of AM parts. However, XCT of metal AM parts can be challenging because of artifacts produced by standard reconstruction algorithms as a result of a confounding effect called “beam hardening.” Beam hardening artifacts complicate the analysis of XCT images and adversely impact the process of detecting defects, such as pores and cracks, which is key to ensuring the quality of the parts being printed. In this work, we propose a novel framework based on using available computer-aided design (CAD) models for parts to be manufactured, accurate XCT simulations, and a deep-neural network to produce high-quality XCT reconstructions from data that are affected by noise and beam hardening. Using extensive experiments with simulated data sets, we demonstrate that our method can significantly improve the reconstruction quality, thereby enabling better detection of defects compared with the state of the art. We also present promising preliminary results of applying the deep networks trained using CAD models to experimental data obtained from XCT of an AM jet-engine turbine blade.

Ziabari, Amir↗
DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Electron tomography unravels new insights into fiber cell wall nanostructure; exploring 3D macromolecular biopolymeric nano-architecture of spruce fiber secondary walls

Lignocellulose biomass has a tremendous potential as renewable biomaterials for fostering the “bio-based society” and circular bioeconomy paradigm. It requires efficient use and breakdown of fiber cell walls containing mainly cellulose, hemicellulose and lignin biopolymers. Despite their great importance, there is an extensive debate on the true structure of fiber walls and knowledge on the macromolecular nano-organization is limited and remains elusive in 3D. We employed dual-axis electron tomography that allows visualization of previously unseen 3D macromolecular organization/biopolymeric nano-architecture of the secondary S2 layer of Norway spruce fiber wall. Unprecedented 3D nano-structural details with novel insights into cellulose microfibrils (~2 nm diameter), macrofibrils, nano-pore network and cell wall chemistry (volume %) across the S2 were explored and quantified including simulation of structure related permeability. Matrix polymer association with cellulose varied between microfibrils and macrofibrils with lignin directly associated with MFs. Simulated bio-nano-mechanical properties revealed stress distribution within the S2 and showed similar properties between the idealized 3D model and the native S2 (actual tomogram). Present work has great potential for significant advancements in lignocellulose research on nano-scale understanding of cell wall assembly/disassembly processes leading to more efficient industrial processes of functionalization, valorization and target modification technologies.

3-D reconstruction↗
NASA NTRS1997

Electrical signal transmission and gap junction regulation in a bone cell network: a cable model for an osteon

A cable model is formulated to estimate the spatial distribution of intracellular electric potential and current, from the cement line to the lumen of an osteon, as the frequency of the loading and the conductance of the gap junction are altered. The model predicts that the characteristic diffusion time for the spread of current along the membrane of the osteocytic processes, 0.03 sec, is nearly the same as the predicted pore pressure relaxation time in Zeng et al. (Annals of Biomedical Engineering. 1994) for the draining of the bone fluid into the osteonal canal. This approximate equality of characteristic times causes the cable to behave as a high-pass, low-pass filter cascade with a maximum in the spectral response for the intracellular potential at approximately 30 Hz. This behavior could be related to the experiments of Rubin and McLeod (Osteoporosis, Academic Press, 1996) which show that live bone appears to be selectively responsive to mechanical loading in a specific frequency range (15-30 Hz) for several species.

NASA Discipline Musculoskeletal↗
DOE OSTI2021

The Application of Convolutional Neural Networks (CNNs) to Recognize Defects in 3D-Printed Parts

Cracks and pores are two common defects in metallic additive manufacturing (AM) parts. In this paper, deep learning-based image analysis is performed for defect (cracks and pores) classification/detection based on SEM images of metallic AM parts. Three different levels of complexities, namely, defect classification, defect detection and defect image segmentation, are successfully achieved using a simple CNN model, the YOLOv4 model and the Detectron2 object detection library, respectively. The tuned CNN model can classify any single defect as either a crack or pore at almost 100% accuracy. The other two models can identify more than 90% of the cracks and pores in the testing images. In addition to the application of static image analysis, defect detection is also successfully applied on a video which mimics the AM process control images. The trained Detectron2 model can identify almost all the pores and cracks that exist in the original video. This study lays a foundation for future in situ process monitoring of the 3D printing process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reconstructing porous media using generative flow networks

One area of intense scientific interest for the study of sandstones, carbonates, and shale at the pore scale is the use of limited image and petrophysical data to generate multiple realizations of a rock’s pore structure. Such images aid efforts to quantify uncertainty in petrophysical properties, including porosity–permeability transforms. We develop and evaluate a deep learning-based method to synthesize porous media volumes using a so-called generative flow model trained on x-ray microscope images of rock texture and pore structure. These models are optimized on a log-likelihood objective and they synthesize large and realistic three-dimensional images. Further, we demonstrate the rapid generation of sandstone image volumes that display realism as gauged by quantitative comparison of topological features using Minkowski functionals of porosity, specific surface area, and the Euler–Poincaré characteristic (i.e., pore connectivity). We also evaluate the single-phase permeability using Navier–Stokes and lattice Boltzmann methods and show that transport properties of the generated samples match measured trends.

58 GEOSCIENCES↗
DOE OSTI2025

A percolating path to green iron

About 1.9 gigatons of steel is produced every year, emitting 8% (3.6 gigatons) of global CO 2 in the process. More than 50% of the CO 2 emissions come from a single step of steel production, known as ironmaking. Hydrogen- based direct reduction (HyDR) of iron oxide to iron has emerged as an emission-free ironmaking alternative. However, multiple physical and chemical phenomena ranging from nanometers to meters inside HyDR reactors alter the microstructure and pore networks in iron oxide pellets, in ways that resist gaseous transport of H 2 /H 2 O, slow reaction rates, and disrupt continuous reactor operation. Using synchrotron nano X-ray computed tomography and percolation theory, we quantify the evolution of pores in iron oxide pellets and demonstrate how nanoscale pore connectivity influences micro- and macroscale flow properties such as permeability, diffusivity, and tortuosity. Our modeling framework connects disparate scales and offers opportunities to accelerate HyDR.

hydrogen↗
DOE OSTI2023

Machine Learning Application to Assess Occurrence and Saturations of Methane Hydrate in Marine Deposits Offshore India

Artificial Neural Networks (ANN) were used to assess methane hydrate occurrence and saturation in marine sediments offshore India. The ANN analysis classifies the gas hydrate occurrence into three types: methane hydrate in pore space, methane hydrate in fractures, or no methane hydrate. Further, predicted saturation characterizes the volume of gas hydrate with respect to the available void volume. Log data collected at six wells, which were drilled during the India National Gas Hydrate Program Expedition 02 (NGHP-02), provided a combination of well log measurements that were used as input for machine learning (ML) models. Well log measurements included density, porosity, electrical resistivity, natural gamma radiation, and acoustic wave velocity. Combinations of well logs used in the ML models provide good overall balanced accuracy (0.79 to 0.86) for the prediction of the gas hydrate occurrence and good accuracy (0.68 to 0.92) for methane hydrate saturation prediction in the marine accumulations against reference data. The accuracy scores indicate that the ML models can successfully predict reservoir characteristics for marine methane hydrate deposits. In conclusion, the results indicate that the ML models can either augment physics-driven methods for assessing the occurrence and saturation of methane hydrate deposits or serve as an independent predictive tool for those characteristics.

58 GEOSCIENCES↗
DOE OSTI2020

Multi-Scale Microfluidics for Transport in Shale Fabric

We develop a microfluidic experimental platform to study solute transport in multi-scale fracture networks with a disparity of spatial scales ranging between two and five orders of magnitude. Using the experimental scaling relationship observed in Marcellus shales between fracture aperture and frequency, the microfluidic design of the fracture network spans all length scales from the micron (1 μ) to the dm (10 dm). This intentional `tyranny of scales’ in the design, a determining feature of shale fabric, introduces unique complexities during microchip fabrication, microfluidic flow-through experiments, imaging, data acquisition and interpretation. Here, we establish best practices to achieve a reliable experimental protocol, critical for reproducible studies involving multi-scale physical micromodels spanning from the Darcy- to the pore-scale (dm to μm). With this protocol, two fracture networks are created: a macrofracture network with fracture apertures between 5 and 500 μm and a microfracture network with fracture apertures between 1 and 500 μm. The latter includes the addition of 1 μm ‘microfractures’, at a bearing of 55°, to the backbone of the former. Comparative analysis of the breakthrough curves measured at corresponding locations along primary, secondary and tertiary fractures in both models allows one to assess the scale and the conditions at which microfractures may impact passive transport.

58 GEOSCIENCES↗
DOE OSTI2022

Deep learning-based quasi-continuum theory for structure of confined fluids

Predicting the structural properties of water and simple fluids confined in nanometer scale pores and channels is essential in, for example, energy storage and biomolecular systems. Classical continuum theories fail to accurately capture the interfacial structure of fluids. In this work, we develop a deep learning-based quasi-continuum theory (DL-QT) to predict the concentration and potential profiles of a Lennard-Jones (LJ) fluid and water confined in a nanochannel. The deep learning model is built based on a convolutional encoder–decoder network (CED) and is applied for high-dimensional surrogate modeling to relate the fluid properties to the fluid–fluid potential. The CED model is then combined with the interatomic potential-based continuum theory to determine the concentration profiles of a confined LJ fluid and confined water. Further, we show that the DL-QT model exhibits robust predictive performance for a confined LJ fluid under various thermodynamic states and for water confined in a nanochannel of different widths. The DL-QT model seamlessly connects molecular physics at the nanoscale with continuum theory by using a deep learning model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Elucidating the Atomic Structures of the Gel Layer Formed during Aluminoborosilicate Glass Dissolution: An Integrated Experimental and Simulation Study

The altered glasses produced during aqueous dissolution of silicate and borosilicate glasses are among the most complex structures to understand at the atomic level due to their amorphous nature, random porosity and various levels of hydration. In this study, we gained insights of the complex atomic structures of altered aluminoborosilicate glasses by combining a range of experimental and computational approaches. The altered glasses were prepared by dissolution of three glasses with varying level of alumina in acid for 7 days. A comprehensive set of experimental [elemental analysis, high-energy X-ray diffraction, 29 Si and 27 Al solid-state nuclear magnetic resonance (NMR), and O 1s X-ray photoelectron spectroscopy (XPS)] and modeling (molecular dynamics (MD) simulations using non-reactive and reactive force fields) were used to study the atomic structures of these altered glasses. Elemental analysis showed that most of the B in the pristine glasses was leached into the solution and was not contained in the altered glass. The 29 Si and 27 Al solid-state NMR spectra revealed that the altered glasses have more polymerized silicate network as compared to those in the pristine glasses due to the reformation of linkages among Si and Al oxygen polyhedral in the altered glasses. The bridging and non-bridging (or hydroxyl O) atoms in the altered glasses were also quantified from their O 1s XPS spectra. Atomic structure models of the altered glasses were constructed by using MD simulations using the reactive force field based on the compositional information obtained from experiments. Various pore structures were generated using the charge-scaling (CS) method by using different initial densities and CS temperatures; the best CS parameters of each alternated glass were then determined by comparing with the experimental structure factors obtained from high energy X-ray diffraction. Pore and the atomic structures and vibrational properties around these pore surfaces were analyzed. Furthermore, these results from this comprehensive study thus provides a realistic insight into the pore morphology, atomic structure, and vibrational properties of altered glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗
DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Neural network–based pore flow field prediction in porous media using super resolution

Direct pore-scale simulations of fluid flow through porous media are computationally expensive to perform for realistic systems. Previous works have demonstrated using the geometry of the microstructure of porous media to predict the velocity fields therein based on neural networks. However, such trained neural networks do not perform well for unseen porous media with a large degree of heterogeneity. In this study we propose that incorporating a coarse velocity field in the input of neural networks is an effective way to improve the prediction performance. The coarse velocity field can be simulated with a low computational cost and provides global information to regularize the ill-posedness of the learning problem, which is usually caused by the use of local geometries due to the computational resource constraints. We show that incorporating the coarse-mesh velocity field significantly improves the prediction accuracy of the fine-mesh velocity field by comparison to the prediction that relies on geometric information alone, especially for the porous medium with a large interior vuggy pore space. We also show the flexibility of training the network in using coarse velocity fields with various resolutions. The results suggest that even using coarse velocity field with a very low resolution, the predictions are still enhanced and close to the ground truths. The feasibility of the method is further demonstrated by testing the trained network on real rocks. This study highlights the merits of incorporating a coarse-mesh velocity field into the input for neural networks, which provides global, physics-based information for the model, thereby improving the model's generalization capability.

42 ENGINEERING↗
DOE OSTI2023

One thousand soils for molecular understanding of belowground carbon cycling

While significant progress has been made in understanding global carbon (C) cycling, the mechanisms regulating belowground C fluxes and storage are still uncertain. New molecular technologies have the power to elucidate these processes, yet we have no widespread standardized implementation of molecular techniques. To address this gap, we introduce the Molecular Observation Network (MONet), a decadal vision from the Environmental Molecular Sciences Laboratory (EMSL), to develop a national network for understanding the molecular composition, physical structure, and hydraulic and biological properties of soil and water. These data are essential for advancing the next generation of multiscale Earth systems models. In this paper, we discuss the 1000 Soils Pilot for MONet, including a description of standardized sampling materials and protocols and a use case to highlight the utility of molecular-level and microstructural measurements for assessing the impacts of wildfire on soil. While the 1000 Soils Pilot generated a plethora of data, we focus on assessments of soil organic matter (SOM) chemistry via Fourier-transform ion cyclotron resonance-mass spectrometry and microstructural properties via X-ray computed tomography to highlight the effects of recent fire history in forested ecosystems on belowground C cycling. We observed decreases in soil respiration, microbial biomass, and potential enzyme activity in soils with high frequency burns. Additionally, the nominal oxidation state of carbon in SOM increased with burn frequency in surface soils. This results in a quantifiable shift in the molecular signature of SOM and shows that wildfire may result in oxidation of SOM and structural changes to soil pore networks that persist into deeper soils.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

The M2 Proton Channel of Influenza Virus: How Does It Work?

The transport of protons across membranes is an essential process for both bioenergetics of modem cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATR), synthesized from adenosine diphosphate. ATR, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this study was: how the same process can be accomplished with the aid of similar, but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC), which is a good model of the biological membranes focusing cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M2 protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M2 protein is 97 amino acids in length, but a fragment 25 amino acids long, which contains a transmembrane domain of 19 amino acids flanked by 3 amino acids on each side, is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules, and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional, important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel, in which the histidine residues were in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel looses its structural integrity and one helix moves away from the remaining three. This result indicates that such a mechanism of proton transport is unlikely. The simulations revealed that translocation along a network of water molecules in the channel and tautomerization of the histidine residues in the M2 proteins in the most likely mechanism of proton transport. The results not only explain how a remarkably simple protein system can efficiently aid in the formation of proton gradients across cell walls, but also suggest how this system can be genetically re-engineered to become a directional, reversible proton pump. Such a pump can provide energy to laboratory-built models of simple cellular systems. If they were successfully constructed it would greatly advance our understanding of the beginnings of life and find important applications in medicine and pharmacology.

Pohorille, Andrew↗