Insight into the K channel's selectivity from binding of K+, Na+ and water to N-methylacetamide.

In potassium channels that conduct K+ selectively over Na+, which sites are occupied by K+ or water and the mechanism of selectivity are unresolved questions. The combination of the energetics and the constraints imposed by the protein structure yield the selective permeation and occupancy. To gain insight into the combination of structure and energetics, we performed density functional theory (DFT) calculations of multiple N-methyl acetamide (NMA) ligands binding to K+ and Na+, relative to hydrated K+ and Na+. NMA is an analogue of the amino acid backbone and provides the carbonyl binding to the ions that occurs in most binding sites of the K+ channel. Unconstrained optimal structures are obtained through geometry optimization calculations of the NMA ligand binding. The complexes formed by 8 NMA binding to the cations have the O atoms positioned in nearly identical locations as the O atoms in the selectivity filter. The transfer free energies between bulk water and K+ or Na+ bound to 8 NMA are almost identical, implying there is no selectivity by a single site. For water optimized with 8 NMA, binding is weak and O atoms are not positioned as in the K+ channel selectivity filter, suggesting that the ions are much more favored than water. Optimal structures of 8 NMA binding with two cations (K+ or Na+) are stable and have lower binding free energy than the optimal structures with just one cation. However, in the Na+ case, the optimal structure deforms and does not match the K+ channel; that is, two bound Na+ are destabilizing. In contrast, the two K+ structure is stabilized and the selectivity free energy favors K+. Overall, this study shows that binding site occupancy and the mechanism for K+ selectivity involves multiple K+ binding in multiple neighboring layers or sites of the K+ channel selectivity filter.

Limited effects of azithromycin on the oropharyngeal microbiome in children with CF and early pseudomonas infection.

BACKGROUND: Tobramycin inhalation solution (TIS) and chronic azithromycin (AZ) have known clinical benefits for children with CF, likely due to antimicrobial and anti-inflammatory activity. The effects of chronic AZ in combination with TIS on the airway microbiome have not been extensively investigated. Oropharyngeal swab samples were collected in the OPTIMIZE multicenter, randomized, placebo-controlled trial examining the addition of AZ to TIS in 198 children with CF and early P. aeruginosa infection. Bacterial small subunit rRNA gene community profiles were determined. The effects of TIS and AZ were assessed on oropharyngeal microbial diversity and composition to uncover whether effects on the bacterial community may be a mechanism of action related to the observed changes in clinical outcomes. RESULTS: Substantial changes in bacterial communities (total bacterial load, diversity and relative abundance of specific taxa) were observed by week 3 of TIS treatment for both the AZ and placebo groups. On average, these shifts were due to changes in non-traditional CF taxa that were not sustained at the later study visits (weeks 13 and 26). Bacterial community measures did not differ between the AZ and placebo groups. CONCLUSIONS: This study provides further evidence that the mechanism for AZ's effect on clinical outcomes is not due solely to action on airway microbial composition.

Building blocks of non-Euclidean ribbons: size-controlled self-assembly via discrete frustrated particles.

Geometric frustration offers a pathway to soft matter self-assembly with controllable finite sizes. While the understanding of frustration in soft matter assembly derives almost exclusively from continuum elastic descriptions, a current challenge is to understand the connection between microscopic physical properties of misfitting "building blocks" and emergent assembly behavior at the mesoscale. We present and analyze a particle-based description of what is arguably the best studied example for frustrated soft matter assembly, negative-curvature ribbon assembly, observed in both assemblies of chiral surfactants and shape-frustrated nanoparticles. Based on our particle model, known as saddle wedge monomers, we numerically test the connection between microscopic shape and interactions of the misfitting subunits and the emergent behavior at the supra-particle scale, specifically focussing on the propagation and relaxation of inter-particle strains, the emergent role of extrinsic shape on frustrated ribbons and the equilibrium regime of finite width selection. Beyond the intuitive role of shape misfit, we show that self-limitation is critically dependent on the finite range of cohesive interactions, with larger size finite assemblies requiring increasing short-range interparticle forces. Additionally, we demonstrate that non-linearities arising from discrete particle interactions alter self-limiting behavior due to both strain-softening in shape-flattened assembly and partial yielding of highly strained bonds, which in turn may give rise to states of hierarchical, multidomain assembly. Tracing the regimes of frustration-limited assembly to the specific microscopic features of misfitting particle shapes and interactions provides necessary guidance for translating the theory of size-programmable assembly into design of intentionally-frustrated colloidal particles.

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics.

Strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties.

Mucosal Microbiota Associated With Eosinophilic Esophagitis and Eosinophilic Gastritis.

OBJECTIVE: The aim of the study was to determine the mucosal microbiota associated with eosinophilic esophagitis (EoE) and eosinophilic gastritis (EoG) in a geographically diverse cohort of patients compared to controls. METHODS: We conducted a prospective study of individuals with eosinophilic gastrointestinal disease (EGID) in the Consortium of Eosinophilic Gastrointestinal Disease Researchers, including pediatric and adult tertiary care centers. Eligible individuals had clinical data, mucosal biopsies, and stool collected. Total bacterial load was determined from mucosal biopsy samples by quantitative polymerase chain reaction (PCR). Community composition was determined by small subunit rRNA gene amplicons. RESULTS: One hundred thirty-nine mucosal biopsies were evaluated corresponding to 93 EoE, 17 EoG, and 29 control specimens (18 esophageal) from 10 sites across the United States. Dominant community members across disease activity differed significantly. When comparing EoE and EoG with controls, the dominant taxa in individuals with EGIDs was increased ( Streptococcus in esophagus; Prevotella in stomach). Specific taxa were associated with active disease for both EoE ( Streptococcus , Gemella ) and EoG ( Leptotrichia ), although highly individualized communities likely impacted statistical testing. Alpha diversity metrics were similar across groups, but with high variability among individuals. Stool analyses did not correlate with bacterial communities found in mucosal biopsy samples and was similar in patients and controls. CONCLUSIONS: Dominant community members ( Streptococcus for EoE, Prevotella for EoG) were different in the mucosal biopsies but not stool of individuals with EGIDs compared to controls; taxa associated with EGIDs were highly variable across individuals. Further study is needed to determine if therapeutic interventions contribute to the observed community differences.

Binding of carboxylate and water to monovalent cations.

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. To gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li+, Na+ and K+. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents.

We develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

Carboxylate binding prefers two cations to one.

Almost all studies of specific ion binding by carboxylates (-COO-) have considered only a single cation, but clustering of ions and ligands is a common phenomenon. We apply density functional theory to investigate how variations in the number of acetate ligands in binding to two monovalent cations affects ion binding preferences. We study a series of monovalent (Li+, Na+, K+, Cs+) ions relevant to experimental work on many topics, including ion channels, battery storage, water purification and solar cells. We find that the preferred optimal structure has 3 acetates except for Cs+, which has 2 acetates. The optimal coordination of the cation by the carboxylate O atoms is 4 for both Na+ and K+, and 3 for Li+ and Cs+. There is a 4-fold coordination minimum just a few kcal mol-1 higher than the optimal 3-fold structure for Li+. For two cations, multiple minima occur in the vicinity of the lowest free energy state. We find that, for Li, Na and K, the preferred optimal structure with two cations is favored over a mixture of single cation complexes, providing a basis for understanding ionic cluster formation that is relevant for engineering proteins and other materials for rapid, selective ion transport.

Influence of Acid Blockade on the Aerodigestive Tract Microbiome in Children With Cystic Fibrosis.

BACKGROUND: Acid blockade is commonly prescribed in patients with cystic fibrosis (CF). Growing concerns, however, exist about its possible role in the pathophysiology of pulmonary infections. We aimed to investigate if acid blockade alters esophageal and respiratory microbiota leading to dysbiosis and inflammation. METHODS: We performed a cross sectional study of children with CF who were either prescribed acid blockade or not. Samples from the gastrointestinal and respiratory tracts were obtained and microbiome analyzed. Mixed effect models were used to compare outcomes between cohorts and across sampling sites. A random subject intercept was included to account for the multiple sampling sites per individual. RESULTS: A cohort of 25 individuals, 44% girls with median age of 13.8 years [IQR 11.2--14.8] were enrolled. Alpha diversity, total bacterial load, and beta diversity were similar across anatomic compartments, across the upper gastrointestinal tract, and in respiratory samples. Similar alpha diversity, total bacterial load, and beta diversity results were also observed when comparing individuals on versus those off acid blockade. IL-8 was elevated in the distal versus proximal esophagus in the whole cohort (P < 0.01). IL-8 concentrations were similar in the distal esophagus in patients on and off acid blockade, but significantly greater in the proximal esophagus of subjects on treatment (P < 0.01). CONCLUSIONS: On the basis of these data, acid blockade use does not appear to influence the microbiome of the aerodigestive tract in children with cystic fibrosis suggesting a complex interplay between these medications and the bacterial composition of the esophagus and lung.

Simulation of polymerization induced phase separation in model thermosets.

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (Tg) and chain length and thus have the potential for phase separation. One crosslinker has a high Tg corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low Tg. A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

Single-stranded nucleic acid elasticity arises from internal electrostatic tension.

Understanding of the conformational ensemble of flexible polyelectrolytes, such as single-stranded nucleic acids (ssNAs), is complicated by the interplay of chain backbone entropy and salt-dependent electrostatic repulsions. Molecular elasticity measurements are sensitive probes of the statistical conformation of polymers and have elucidated ssNA conformation at low force, where electrostatic repulsion leads to a strong excluded volume effect, and at high force, where details of the backbone structure become important. Here, we report measurements of ssDNA and ssRNA elasticity in the intermediate-force regime, corresponding to 5- to 100-pN forces and 50-85% extension. These data are explained by a modified wormlike chain model incorporating an internal electrostatic tension. Fits to the elastic data show that the internal tension decreases with salt, from [Formula: see text]5 pN under 5 mM ionic strength to near zero at 1 M. This decrease is quantitatively described by an analytical model of electrostatic screening that ascribes to the polymer an effective charge density that is independent of force and salt. Our results thus connect microscopic chain physics to elasticity and structure at intermediate scales and provide a framework for understanding flexible polyelectrolyte elasticity across a broad range of relative extensions.

Self-assembled highly ordered acid layers in precisely sulfonated polyethylene produce efficient proton transport.

Recent advances in polymer synthesis have allowed remarkable control over chain microstructure and conformation. Capitalizing on such developments, here we create well-controlled chain folding in sulfonated polyethylene, leading to highly uniform hydrated acid layers of subnanometre thickness with high proton conductivity. The linear polyethylene contains sulfonic acid groups pendant to precisely every twenty-first carbon atom that induce tight chain folds to form the hydrated layers, while the methylene segments crystallize. The proton conductivity is on par with Nafion 117, the benchmark for fuel cell membranes. We demonstrate that well-controlled hairpin chain folding can be utilized for proton conductivity within a crystalline polymer structure, and we project that this structure could be adapted for ion transport. This layered polyethylene-based structure is an innovative and versatile design paradigm for functional polymer membranes, opening doors to efficient and selective transport of other ions and small molecules on appropriate selection of functional groups.

Polyelectrolyte Conformation Controlled by a Trivalent-Rich Ion Jacket.

The configuration of charged polymers is heavily dependent on interactions with surrounding salt ions, typically manifesting as a sensitivity to the bulk ionic strength. Here, we use single-molecule mechanical measurements to show that a charged polysaccharide, hyaluronic acid, shows a surprising regime of insensitivity to ionic strength in the presence of trivalent ions. Using simulations and theory, we propose that this is caused by the formation of a "jacket" of ions, tightly associated with the polymer, whose charge (and thus effect on configuration) is robust against changes in solution composition.

Diverse balances of tubulin interactions and shape change drive and interrupt microtubule depolymerization.

Microtubules are stiff biopolymers that self-assemble via the addition of GTP-tubulin (αß-dimer bound to GTP), but hydrolysis of GTP- to GDP-tubulin within the tubules destabilizes them toward catastrophically-fast depolymerization. The molecular mechanisms and features of the individual tubulin proteins that drive such behavior are still not well-understood. Using molecular dynamics simulations of whole microtubules built from a coarse-grained model of tubulin, we demonstrate how conformational shape changes (i.e., deformations) in subunits that frustrate tubulin-tubulin binding within microtubules drive depolymerization of stiff tubules via unpeeling "ram's horns" consistent with experiments. We calculate the sensitivity of these behaviors to the length scales and strengths of binding attractions and varying degrees of binding frustration driven by subunit shape change, and demonstrate that the dynamic instability and mechanical properties of microtubules can be produced based on either balanced or imbalanced strengths of lateral and vertical binding attractions. Finally, we show how catastrophic depolymerization can be interrupted by small regions of the microtubule containing undeformed dimers, corresponding to incomplete lattice hydrolysis. The results demonstrate a mechanism by which microtubule rescue can occur.

Catastrophic depolymerization of microtubules driven by subunit shape change.

Microtubules exhibit a dynamic instability between growth and catastrophic depolymerization. GTP-tubulin (αß-dimer bound to GTP) self-assembles, but dephosphorylation of GTP- to GDP-tubulin within the tubule results in destabilization. While the mechanical basis for destabilization is not fully understood, one hypothesis is that dephosphorylation causes tubulin to change shape, frustrating bonds and generating stress. To test this idea, we perform molecular dynamics simulations of microtubules built from coarse-grained models of tubulin, incorporating a small compression of α-subunits associated with dephosphorylation in experiments. We find that this shape change induces depolymerization of otherwise stable systems via unpeeling "ram's horns" characteristic of microtubules. Depolymerization can be averted by caps with uncompressed α-subunits, i.e., GTP-rich end regions. Thus, the shape change is sufficient to yield microtubule behavior.

Importance of corners in fracture of highly crosslinked polymeric adhesives.

Very large molecular dynamics simulations with open ends between two solid adherends have been performed treating tensile deformation of coarse-grained, highly crosslinked polymer networks modeling epoxy systems. The open boundary and the presence of corners dramatically alter the fracture behavior. In contrast to systems with periodic boundaries, the failure strain decreases with increasing system size until a critical size is reached. This decrease greatly reduces the difference in the crack initiation strains between simulation and experiment. In the open geometry, the sides of the polymer network contract inward forming wedge shaped corners. The stress and strain are concentrated in the corners where the shear component is present and large. The nonuniformity of the strain results in accumulation of bond breaking in the corners and crack initiation there. Moreover, the corner strain is system size dependent, which results in a system size dependence of the failure strain.

The effect of chain stiffness and salt on the elastic response of a polyelectrolyte.

We present simulations of the force-extension curves of strong polyelectrolytes with varying intrinsic stiffness as well as specifically treating hyaluronic acid, a polyelectrolyte of intermediate stiffness. Whereas fully flexible polyelectrolytes show a high-force regime where extension increases nearly logarithmically with force, we find that the addition of even a small amount of stiffness alters the short-range structure and removes this logarithmic elastic regime. This further confirms that the logarithmic regime is a consequence of the short-ranged "wrinkles" in the flexible chain. As the stiffness increases, the force-extension curves tend toward and reach the wormlike chain behavior. Using the screened Coulomb potential and a simple bead-spring model, the simulations are able to reproduce the hyaluronic acid experimental force-extension curves for salt concentrations ranging from 1 to 500 mM. Furthermore, the simulation data can be scaled to a universal curve like the experimental data. The scaling analysis is consistent with the interpretation that, in the low-salt limit, the hyaluronic acid chain stiffness scales with salt with an exponent of -0.7, rather than either of the two main theoretical predictions of -0.5 and -1. Furthermore, given the conditions of the simulation, we conclude that this exponent value is not due to counterion condensation effects, as had previously been suggested.

Chain Folding Produces a Multilayered Morphology in a Precise Polymer: Simulations and Experiments.

Precise control over polymer architecture unlocks the potential for engineered self-assembled crystal structures with useful features on the nanometer length scale. Here we elucidate the structure of the ordered phase of a semicrystalline, functional polyethylene having a precise linear architecture, namely, pendant carboxylic acid groups precisely every 21st backbone carbon atom. By comparing the results of atomistic molecular dynamics simulations with experimental X-ray scattering and Raman spectroscopy data, we find that the polymer chains are folded in a hairpin manner near each carboxylic acid group, giving rise to multiple embedded layers of functional groups that have an interlayer distance of 2.5 nm. This is in contrast to other precise polyethylenes, where the chains are mostly trans within the crystals. Such layers could act as two-dimensional pathways for ionic or molecular transport given an appropriate choice of functional group.

Airway microbiota across age and disease spectrum in cystic fibrosis.

Our objectives were to characterise the microbiota in cystic fibrosis (CF) bronchoalveolar lavage fluid (BALF), and determine its relationship to inflammation and disease status.BALF from paediatric and adult CF patients and paediatric disease controls undergoing clinically indicated bronchoscopy was analysed for total bacterial load and for microbiota by 16S rDNA sequencing.We examined 191 BALF samples (146 CF and 45 disease controls) from 13 CF centres. In CF patients aged <2 years, nontraditional taxa (e.gStreptococcus, Prevotella and Veillonella) constituted ∼50% of the microbiota, whereas in CF patients aged ≥6 years, traditional CF taxa (e.gPseudomonas, Staphylococcus and Stenotrophomonas) predominated. Sequencing detected a dominant taxon not traditionally associated with CF (e.gStreptococcus or Prevotella) in 20% of CF BALF and identified bacteria in 24% of culture-negative BALF. Microbial diversity and relative abundance of Streptococcus, Prevotella and Veillonella were inversely associated with airway inflammation. Microbiota communities were distinct in CF compared with disease controls, but did not differ based on pulmonary exacerbation status in CF.The CF microbiota detected in BALF differs with age. In CF patients aged <2 years, Streptococcus predominates, whereas classic CF pathogens predominate in most older children and adults.

Longitudinal and Source-to-Tap New Orleans, LA, U.S.A. Drinking Water Microbiology.

The two municipal drinking water systems of New Orleans, LA, U.S.A. were sampled to compare the microbiology of independent systems that treat the same surface water from the Mississippi River. To better understand temporal trends and sources of microbiology delivered to taps, these treatment plants and distribution systems were subjected to source-to-tap sampling over four years. Both plants employ traditional treatment by chloramination, applied during or after settling, followed by filtration before distribution in a warm, low water age system. Longitudinal samples indicated microbiology to have stability both spatially and temporally, and between treatment plants and distribution systems. Disinfection had the greatest impact on microbial composition, which was further refined by filtration and influenced by distribution and premise plumbing. Actinobacteria spp. exhibited trends with treatment. In particular, Mycobacterium spp., very low in finished waters, occurred idiosyncratically at high levels in some tap waters, indicating distribution and/or premise plumbing as main contributors of mycobacteria. Legionella spp., another genus containing potential opportunistic pathogens, also occurred ubiquitously. Source water microbiology was most divergent from tap water, and each step of treatment brought samples more closely similar to tap waters.