Learning composition-transferable coarse-grained models: Designing external potential ensembles to maximize thermodynamic information.

Achieving thermodynamic faithfulness and transferability across state points is an outstanding challenge in the bottom-up coarse graining of molecular models, with many efforts focusing on augmenting the form of coarse-grained interaction potentials to improve transferability. Here, we revisit the critical role of the simulation ensemble and the possibility that even simple models can be made more predictive through a smarter choice of ensemble. We highlight the efficacy of coarse graining from ensembles where variables conjugate to the thermodynamic quantities of interest are forced to respond to applied perturbations. For example, to learn activity coefficients, it is natural to coarse grain from ensembles with spatially varying external potentials applied to one species to force local composition variations and fluctuations. We apply this strategy to coarse grain both an atomistic model of water and methanol and a binary mixture of spheres interacting via Gaussian repulsions and demonstrate near-quantitative capture of activity coefficients across the whole composition range. Furthermore, the approach is able to do so without explicitly measuring and targeting activity coefficients during the coarse graining process; activity coefficients are only computed after-the-fact to assess accuracy. We hypothesize that ensembles with applied thermodynamic potentials are more "thermodynamically informative." We quantify this notion of informativeness using the Fisher information metric, which enables the systematic design of optimal bias potentials that promote the learning of thermodynamically faithful models. The Fisher information is related to variances of structural variables, highlighting the physical basis underlying the Fisher information's utility in improving coarse-grained models.

Conjugal Transfer, Whole-Genome Sequencing, and Plasmid Analysis of Four mcr-1-Bearing Isolates from U.S. Patients.

Four Enterobacteriaceae clinical isolates bearing mcr-1 gene-harboring plasmids were characterized. All isolates demonstrated the ability to transfer colistin resistance to Escherichia coli; plasmids were stable in conjugants after multiple passages on nonselective media. mcr-1 was located on an IncX4 (n = 3) or IncN (n = 1) plasmid. The IncN plasmid harbored 13 additional antimicrobial resistance genes. Results indicate that the mcr-1-bearing plasmids in this study were highly transferable in vitro and stable in the recipients.

E22G Pathogenic Mutation of ß-Amyloid (Aß) Enhances Misfolding of Aß40 by Unexpected Prion-like Cross Talk between Aß42 and Aß40.

Cross-seeding of misfolded amyloid proteins is postulated to induce cross-species infection of prion diseases. In sporadic Alzheimer's disease (AD), misfolding of 42-residue ß-amyloid (Aß) is widely considered to trigger amyloid plaque deposition. Despite increasing evidence that misfolded Aß mimics prions, interactions of misfolded 42-residue Aß42 with more abundant 40-residue Aß40 in AD are elusive. This study presents in vitro evidence that a heterozygous E22G pathogenic ("Arctic") mutation of Aß40 can enhance misfolding of Aß via cross-seeding from wild-type (WT) Aß42 fibril. Thioflavin T (ThT) fluorescence analysis suggested that misfolding of E22G Aß40 was enhanced by adding 5% (w/w) WT Aß42 fibril as "seed", whereas WT Aß40 was unaffected by Aß42 fibril seed. 13C SSNMR analysis revealed that such cross-seeding prompted formation of E22G Aß40 fibril that structurally mimics the seed Aß42 fibril, suggesting unexpected cross talk of Aß isoforms that potentially promotes early onset of AD. The SSNMR approach is likely applicable to elucidate structural details of heterogeneous amyloid fibrils produced in cross-seeding for amyloids linked to neurodegenerative diseases.

Multicenter Evaluation of the Modified Carbapenem Inactivation Method and the Carba NP for Detection of Carbapenemase-Producing Pseudomonas aeruginosa and Acinetobacter baumannii.

The purpose of this study was to develop the modified carbapenem inactivation method (mCIM) for the detection of carbapenemase-producing Pseudomonas aeruginosa (CP-PA) and carbapenemase-producing Acinetobacter baumannii (CP-AB) and perform a multicenter evaluation of the mCIM and Carba NP tests for these nonfermenters. Thirty P. aeruginosa and 30 A. baumannii isolates previously characterized by whole-genome sequencing from the CDC-FDA Antibiotic Resistance Isolate Bank were evaluated, including CP isolates (Ambler class A, B, and D), non-carbapenemase-producing (non-CP) carbapenem-resistant isolates, and carbapenem-susceptible isolates. Initial comparison of a 1-µl versus 10-µl loop inoculum for the mCIM was performed by two testing sites and showed that 10 µl was required for reliable detection of carbapenemase production among P. aeruginosa and A. baumannii Ten testing sites then evaluated the mCIM using a 10-µl loop inoculum. Overall, the mean sensitivity and specificity of the mCIM for detection of CP-PA across all 10 sites were 98.0% (95% confidence interval [CI], 94.3 to 99.6; range, 86.7 to 100) and 95% (95% CI, 89.8 to 97.7; range, 93.3 to 100), whereas the mean sensitivity and specificity among CP-AB were 79.8% (95% CI, 74.0 to 84.9; range, 36.3 to 95.7) and 52.9% (95% CI, 40.6 to 64.9; range, 28.6 to 100), respectively. At three sites that evaluated the performance of the Carba NP test using the same set of isolates, the mean sensitivity and specificity of the Carba NP test were 97.8% (95% CI, 88.2 to 99.9; range, 93.3 to 100) and 97.8% (95% CI, 88.2 to 99.9; range, 93.3 to 100) for P. aeruginosa and 18.8% (95% CI, 10.4 to 30.1; range, 8.7 to 26.1) and 100% (95% CI, 83.9 to 100; range, 100) for A. baumannii Overall, we found both the mCIM and the Carba NP test to be accurate for detection of carbapenemase production among P. aeruginosa isolates and less reliable for use with A. baumannii isolates.

Cassandra: An open source Monte Carlo package for molecular simulation.

Cassandra is an open source atomistic Monte Carlo software package that is effective in simulating the thermodynamic properties of fluids and solids. The different features and algorithms used in Cassandra are described, along with implementation details and theoretical underpinnings to various methods used. Benchmark and example calculations are shown, and information on how users can obtain the package and contribute to it are provided. © 2017 Wiley Periodicals, Inc.

Discrete Fractional Component Monte Carlo Simulation Study of Dilute Nonionic Surfactants at the Air-Water Interface.

We present a newly developed Monte Carlo scheme to predict bulk surfactant concentrations and surface tensions at the air-water interface for various surfactant interfacial coverages. Since the concentration regimes of these systems of interest are typically very dilute (≪10-5 mol. frac.), Monte Carlo simulations with the use of insertion/deletion moves can provide the ability to overcome finite system size limitations that often prohibit the use of modern molecular simulation techniques. In performing these simulations, we use the discrete fractional component Monte Carlo (DFCMC) method in the Gibbs ensemble framework, which allows us to separate the bulk and air-water interface into two separate boxes and efficiently swap tetraethylene glycol surfactants C10E4 between boxes. Combining this move with preferential translations, volume biased insertions, and Wang-Landau biasing vastly enhances sampling and helps overcome the classical "insertion problem", often encountered in non-lattice Monte Carlo simulations. We demonstrate that this methodology is both consistent with the original molecular thermodynamic theory (MTT) of Blankschtein and co-workers, as well as their recently modified theory (MD/MTT), which incorporates the results of surfactant infinite dilution transfer free energies and surface tension calculations obtained from molecular dynamics simulations.

Molecular Mechanism of Ionic-Liquid-Induced Membrane Disruption: Morphological Changes to Bilayers, Multilayers, and Vesicles.

The application of ionic liquids (ILs) in many industrially relevant processes provides an urgent need to better understand their molecular interactions with biological systems. A detailed understanding of the cytotoxicity mechanism of ILs can be helpful in facilitating the molecular design of nontoxic ILs. Using coarse-grained molecular dynamics (MD) simulations, we investigate the effects of imidazolium-based ILs on several lipid bilayer morphologies. Our results demonstrate that the asymmetric insertion of IL cations into one side of a lipid bilayer leaflet enhances the leaflet strain, which upon reaching a critical value triggers a morphological disruption in the bilayer. Consistently, the bending modulus of the bilayer is reduced by 1 to 2 orders of magnitude relative to that of an IL-free planar bilayer prior to the disruption event. Our results suggest that ILs that can easily insert into the lipid bilayer without diffusing across or inducing lipid flip-flop can be more disruptive to a lipid biomembrane.

Capturing a reactive state of amyloid aggregates: NMR-based characterization of copper-bound Alzheimer disease amyloid ß-fibrils in a redox cycle.

The interaction of redox-active copper ions with misfolded amyloid ß (Aß) is linked to production of reactive oxygen species (ROS), which has been associated with oxidative stress and neuronal damages in Alzheimer disease. Despite intensive studies, it is still not conclusive how the interaction of Cu(+)/Cu(2+) with Aß aggregates leads to ROS production even at the in vitro level. In this study, we examined the interaction between Cu(+)/Cu(2+) and Aß fibrils by solid-state NMR (SSNMR) and other spectroscopic methods. Our photometric studies confirmed the production of ~60 µM hydrogen peroxide (H2O2) from a solution of 20 µM Cu(2+) ions in complex with Aß(1-40) in fibrils ([Cu(2+)]/[Aß] = 0.4) within 2 h of incubation after addition of biological reducing agent ascorbate at the physiological concentration (~1 mM). Furthermore, SSNMR (1)H T1 measurements demonstrated that during ROS production the conversion of paramagnetic Cu(2+) into diamagnetic Cu(+) occurs while the reactive Cu(+) ions remain bound to the amyloid fibrils. The results also suggest that O2 is required for rapid recycling of Cu(+) bound to Aß back to Cu(2+), which allows for continuous production of H2O2. Both (13)C and (15)N SSNMR results show that Cu(+) coordinates to Aß(1-40) fibrils primarily through the side chain Nδ of both His-13 and His-14, suggesting major rearrangements from the Cu(2+) coordination via Nε in the redox cycle. (13)C SSNMR chemical shift analysis suggests that the overall Aß conformations are largely unaffected by Cu(+) binding. These results present crucial site-specific evidence of how the full-length Aß in amyloid fibrils offers catalytic Cu(+) centers.

Amphiphilic interactions of ionic liquids with lipid biomembranes: a molecular simulation study.

Current bottlenecks in the large-scale commercial use of many ionic liquids (ILs) include their high costs, low biodegradability, and often unknown toxicities. As a proactive effort to better understand the molecular mechanisms of ionic liquid toxicities, the work herein presents a comprehensive molecular simulation study on the interactions of 1-n-alkyl-3-methylimidazolium-based ILs with a phosphatidylcholine (PC) lipid bilayer. We explore the effects of increasing alkyl chain length (n = 4, 8, and 12) in the cation and anion hydrophobicity on the interactions with the lipid bilayer. Bulk atomistic molecular dynamics (MD) simulations performed at millimolar (mM) IL concentrations show spontaneous insertion of cations into the lipid bilayer regardless of the alkyl chain length and a favorable orientational preference once a cation is inserted. Cations also exhibit the ability to "flip" inside the lipid bilayer (as is common for amphiphiles) if partially inserted with an unfavorable orientation. Moreover, structural analysis of the lipid bilayer show that cationic insertion induces roughening of the bilayer surface, which may be a precursor to bilayer disruption. To overcome the limitation in the timescale of our simulations, free energies for a single IL cation and anion insertion have been determined based on potential of mean force calculations. These results show a decrease in free energy in response to both short and long alkyl chain IL cation insertion, and likewise for a single hydrophobic anion insertion, but an increase in free energy for the insertion of a hydrophilic chloride anion. Both bulk MD simulations and free energy calculations suggest that toxicity mechanisms toward biological systems are likely caused by ILs behaving as ionic surfactants. [Yoo et al., Soft Matter, 2014].

Predicting surfactant phase behavior with a molecularly informed field theory.

HYPOTHESIS: The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations. SIMULATION WORK: We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements. FINDINGS: The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.

Upper lumbar spine far lateral disc herniations masquerading as peripheral nerve sheath tumors: illustrative cases.

BACKGROUND: Migratory disc herniations can mimic neoplasms clinically and on imaging. Far lateral lumbar disc herniations usually compress the exiting nerve root and can be challenging to distinguish from a nerve sheath tumor due to the proximity of the nerve and characteristics on magnetic resonance imaging (MRI). These lesions can occasionally present in the upper lumbar spine region at the L1-2 and L2-3 levels. OBSERVATIONS: The authors describe 2 extraforaminal lesions in the far lateral space at the L1-2 and L2-3 levels, respectively. On MRI, both lesions tracked along the corresponding exiting nerve roots with avid postcontrast rim enhancement and edema in the adjacent muscle tissue. Thus, they were initially concerning for peripheral nerve sheath tumors. One patient underwent fluorodeoxyglucose positron emission tomography-computed tomography (FDG PET-CT) screening and demonstrated moderate FDG uptake on PET-CT scan. In both cases, intraoperative and postoperative pathology revealed fibrocartilage disc fragments. LESSONS: Differential diagnosis for lumbar far lateral lesions that are peripherally enhancing on MRI should include migratory disc herniation, regardless of the level of the disc herniations. Accurate preoperative diagnosis can aid in decision making for management, surgical approach, and resection.

Development of a Broth Microdilution Method To Characterize Chlorhexidine MICs among Bacteria Collected from 2005 to 2019 at Three U.S. Sites.

Chlorhexidine bathing to prevent transmission of multidrug-resistant organisms has been adopted by many U.S. hospitals, but increasing chlorhexidine use has raised concerns about possible emergence of resistance. We sought to establish a broth microdilution method for determining chlorhexidine MICs and then used the method to evaluate chlorhexidine MICs for bacteria that can cause health care-associated infections. We adapted a broth microdilution method for determining chlorhexidine MICs, poured panels, established quality control ranges, and tested Staphylococcus aureus, Escherichia coli, Klebsiella pneumoniae, and Enterobacter cloacae complex isolates collected at three U.S. sites. Chlorhexidine MICs were determined for 535 isolates including 129 S. aureus, 156 E. coli, 142 K. pneumoniae, and 108 E. cloacae complex isolates. The respective MIC distributions for each species ranged from 1 to 8 mg/L (MIC50 = 2 mg/L and MIC90 = 4 mg/L), 1 to 64 mg/L (MIC50 = 2 mg/L and MIC90 = 4 mg/L), 4 to 64 mg/L (MIC50 = 16 mg/L and MIC90 = 32 mg/L), and 1 to >64 mg/L (MIC50 = 16 mg/L and MIC90 = 64 mg/L). We successfully adapted a broth microdilution procedure that several laboratories were able to use to determine the chlorhexidine MICs of bacterial isolates. This method could be used to investigate whether chlorhexidine MICs are increasing. IMPORTANCE Chlorhexidine bathing to prevent transmission of multidrug-resistant organisms and reduce health care-associated infections has been adopted by many hospitals. There is concern about the possible unintended consequences of using this agent widely. One possible unintended consequence is decreased susceptibility to chlorhexidine, but there are not readily available methods to perform this evaluation. We developed a method for chlorhexidine MIC testing that can be used to evaluate for possible unintended consequences.

Involvement of a putative ATP-Binding Cassette (ABC) Involved in manganese transport in virulence of Listeria monocytogenes.

Listeria monocytogenes is a foodborne pathogen and the causative agent of listeriosis, a disease associated with high fatality (20-30%) and hospitalization rates (>95%). ATP-Binding Cassette (ABC) transporters have been demonstrated to be involved in the general stress response. In previous studies, in-frame deletion mutants of the ABC transporter genes, LMOf2365_1875 and LMOf2365_1877, were constructed and analyzed; however, additional work is needed to investigate the virulence potential of these deletion mutants. In this study, two in vitro methods and one in vivo model were used to investigate the virulence potential of in-frame deletion mutants of ABC transporter genes. First, the invasion efficiency in host cells was measured using the HT-29 human cell line. Second, cell-to-cell spread activity was measured using a plaque forming assay. Lastly, virulence potential of the mutants was tested in the Galleria mellonella wax moth model. Our results demonstrated that the deletion mutant, ⊿LMOf2365_1875, displayed decreased invasion and cell-to-cell spread efficiency in comparison to the wild-type, LMOf2365, indicating that LMOf2365_1875 may be required for virulence. Furthermore, the reduced virulence of these mutants was confirmed using the Galleria mellonella wax moth model. In addition, the expression levels of 15 virulence and stress-related genes were analyzed by RT-PCR assays using stationary phase cells. Our results showed that virulence-related gene expression levels from the deletion mutants were elevated (15/15 genes from ⊿LMOf2365_1877 and 7/15 genes from ⊿LMOf2365_1875) compared to the wild type LMOf2365, suggesting that ABC transporters may negatively regulate virulence gene expression under specific conditions. The expression level of the stress-related gene, clpE, also was increased in both deletion mutants, indicating the involvement of ABC transporters in the stress response. Taken together, our findings suggest that ABC transporters may be used as potential targets to develop new therapeutic strategies to control L. monocytogenes.

Bias free multiobjective active learning for materials design and discovery.

The design rules for materials are clear for applications with a single objective. For most applications, however, there are often multiple, sometimes competing objectives where there is no single best material and the design rules change to finding the set of Pareto optimal materials. In this work, we leverage an active learning algorithm that directly uses the Pareto dominance relation to compute the set of Pareto optimal materials with desirable accuracy. We apply our algorithm to de novo polymer design with a prohibitively large search space. Using molecular simulations, we compute key descriptors for dispersant applications and drastically reduce the number of materials that need to be evaluated to reconstruct the Pareto front with a desired confidence. This work showcases how simulation and machine learning techniques can be coupled to discover materials within a design space that would be intractable using conventional screening approaches.

Molecularly Informed Field Theories from Bottom-up Coarse-Graining.

Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.

PAX6 is expressed in pancreatic cancer and actively participates in cancer progression through activation of the MET tyrosine kinase receptor gene.

Tumors of the exocrine pancreas have a poor prognosis. Several proteins are overexpressed in this cancer type, including the MET tyrosine kinase receptor and the transcription factor PAX6. In this report, we find that PAX6(5a), an alternately spliced variant form of PAX6, is expressed in pancreatic carcinoma cell lines at higher levels than the canonical PAX6 protein. Both protein forms of PAX6 bind directly to an enhancer element in the MET promoter and activate the expression of the MET gene. In addition, inhibition of PAX6 transcripts leads to a decline in cell growth and survival, differentiation, and a concurrent reduction of MET protein expression. These data support a model for a neoplastic pathway, where expression of a transcription factor from development activates the MET receptor, a protein that has been directly linked to protumorigenic processes of resisting apoptosis, tumor growth, invasion, and metastasis.

LMOf2365_0442 Encoding for a Fructose Specific PTS Permease IIA May Be Required for Virulence in L. monocytogenes Strain F2365.

Listeria monocytogenes is a foodborne pathogen that causes listeriosis, which is a major public health concern due to the high fatality rate. LMOf2365_0442, 0443, and 0444 encode for fructose-specific EIIABC components of phosphotransferase transport system (PTS) permease that is responsible for sugar transport. In previous studies, in-frame deletion mutants of a putative fructose-specific PTS permease (LMOf2365_0442, 0443, and 0444) were constructed and analyzed. However, the virulence potential of these deletion mutants has not been studied. In this study, two in vitro methods were used to analyze the virulence potential of these L. monocytogenes deletion mutants. First, invasion assays were used to measure the invasion efficiencies to host cells using the human HT-29 cell line. Second, plaque forming assays were used to measure cell-to-cell spread in host cells. Our results showed that the deletion mutant ΔLMOf2365_0442 had reduced invasion and cell-to-cell spread efficiencies in human cell line compared to the parental strain LMOf2365, indicating that LMOf2365_0442 encoding for a fructose specific PTS permease IIA may be required for virulence in L. monocytogenes strain F2365. In addition, the gene expression levels of 15 virulence and stress-related genes were analyzed in the stationary phase cells of the deletion mutants using RT-PCR assays. Virulence-related gene expression levels were elevated in the deletion mutants ΔLMOf2365_0442-0444 compared to the wild type parental strain LMOf2365, indicating the down-regulation of virulence genes by this PTS permease in L. monocytogenes. Finally, stress-related gene clpC expression levels were also increased in all of the deletion mutants, suggesting the involvement of this PTS permease in stress response. Furthermore, these deletion mutants displayed the same pressure tolerance and the same capacity for biofilm formation compared to the wild-type parental strain LMOf2365. In summary, our findings suggest that the LMOf2365_0442 gene can be used as a potential target to develop inhibitors for new therapeutic and pathogen control strategies for public health.

Molecular mechanisms of ionic liquid cytotoxicity probed by an integrated experimental and computational approach.

Ionic liquids (ILs) are salts that remain liquid down to low temperatures, and sometimes well below room temperature. ILs have been called "green solvents" because of their extraordinarily low vapor pressure and excellent solvation power, but ecotoxicology studies have shown that some ILs exhibit greater toxicity than traditional solvents. A fundamental understanding of the molecular mechanisms responsible for IL toxicity remains elusive. Here we show that one mode of IL toxicity on unicellular organisms is driven by swelling of the cell membrane. Cytotoxicity assays, confocal laser scanning microscopy, and molecular simulations reveal that IL cations nucleate morphological defects in the microbial cell membrane at concentrations near the half maximal effective concentration (EC50) of several microorganisms. Cytotoxicity increases with increasing alkyl chain length of the cation due to the ability of the longer alkyl chain to more easily embed in, and ultimately disrupt, the cell membrane.