RESUMEN
Block copolymer electrolytes (BCE) such as polystyrene-block-poly(ethylene oxide) (SEO) blended with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and composed of mechanically robust insulating and rubbery conducting nanodomains are promising solid-state electrolytes for Li batteries. Here, we compare ionic solvation, association, distribution, and conductivity in SEO-LiTFSI BCEs and their homopolymer PEO-LiTFSI analogs toward a fundamental understanding of the maximum in conductivity and transport mechanisms as a function of salt concentration. Ionic conductivity measurements reveal that SEO-LiTFSI and PEO-LiTFSI exhibit similar behaviors up to a Li/EO ratio of 1/12, where roughly half of the available solvation sites in the system are filled, and conductivity is maximized. As the Li/EO ratios increase to 1/5 the conductivity, of the PEO-LiTFSI drops nearly 3-fold, while the conductivity of SEO-LiTFSI remains constant. FTIR spectroscopy reveals that additional Li cations in the homopolymer electrolyte are complexed by additional EO units when the Li/EO ratio exceeds 1/12, while in the BCE, the proportion of complexed and uncomplexed EO units remains constant; Raman spectroscopy data at the same concentrations show that Li cations in the SEO-LiTFSI samples tend to coordinate more to their counteranions. Atomistic-scale molecular dynamics simulations corroborate these results and further show that associated ions tend to segregate to the SEO-LiTFSI domain interfaces. The opportunity for "excess" salt to be sequestered at BCE interfaces results in the retention of an optimum ratio of uncompleted and complexed PEO solvation sites in the middle of the conductive nanodomains of the BCE and maximized conductivity over a broad range of salt concentrations.
RESUMEN
Force-field development has undergone a revolution in the past decade with the proliferation of quantum chemistry based parametrizations and the introduction of machine learning approximations of the atomistic potential energy surface. Nevertheless, transferable force fields with broad coverage of organic chemical space remain necessary for applications in materials and chemical discovery where throughput, consistency, and computational cost are paramount. Here, we introduce a force-field development framework called Topology Automated Force-Field Interactions (TAFFI) for developing transferable force fields of varying complexity against an extensible database of quantum chemistry calculations. TAFFI formalizes the concept of atom typing and makes it the basis for generating systematic training data that maintains a one-to-one correspondence with force-field terms. This feature makes TAFFI arbitrarily extensible to new chemistries while maintaining internal consistency and transferability. As a demonstration of TAFFI, we have developed a fixed-charge force-field, TAFFI-gen, from scratch that includes coverage for common organic functional groups that is comparable to established transferable force fields. The performance of TAFFI-gen was benchmarked against OPLS and GAFF for reproducing several experimental properties of 87 organic liquids. The consistent performance of these force fields, despite their distinct origins, validates the TAFFI framework while also providing evidence of the representability limitations of fixed-charge force fields.
Asunto(s)
Aprendizaje Automático , Compuestos OrgánicosRESUMEN
A configurational sampling algorithm based on nested layerings of Markov chains (Layered Nested Markov Chain Monte Carlo or L-NMCMC) is presented for simulations of systems characterized by rugged free energy landscapes. The layerings are generated using a set of auxiliary potential energy surfaces. The implementation of the method is demonstrated in the context of a rugged, two-dimensional potential energy surface. The versatility of the algorithm is next demonstrated on a simple, many-body system, namely, a canonical Lennard-Jones fluid in the liquid state. In that example, different layering schemes and auxiliary potentials are used, including variable cutoff distances and excluded-volume tempering. In addition to calculating a variety of properties of the system, it is also shown that L-NMCMC, when combined with a free-energy perturbation formalism, provides a straightforward means to construct approximate free-energy surfaces at no additional computational cost using the sampling distributions of each auxiliary Markov chain. The proposed L-NMCMC scheme is general in that it could be complementary to any number of methods that rely on sampling from a target distribution or methods that exploit a hierarchy of time scales and/or length scales through decomposition of the potential energy.
RESUMEN
Molecular simulation has emerged as an essential tool for modern-day research, but obtaining proper results and making reliable conclusions from simulations requires adequate sampling of the system under consideration. To this end, a variety of methods exist in the literature that can enhance sampling considerably, and increasingly sophisticated, effective algorithms continue to be developed at a rapid pace. Implementation of these techniques, however, can be challenging for experts and non-experts alike. There is a clear need for software that provides rapid, reliable, and easy access to a wide range of advanced sampling methods and that facilitates implementation of new techniques as they emerge. Here we present SSAGES, a publicly available Software Suite for Advanced General Ensemble Simulations designed to interface with multiple widely used molecular dynamics simulations packages. SSAGES allows facile application of a variety of enhanced sampling techniques-including adaptive biasing force, string methods, and forward flux sampling-that extract meaningful free energy and transition path data from all-atom and coarse-grained simulations. A noteworthy feature of SSAGES is a user-friendly framework that facilitates further development and implementation of new methods and collective variables. In this work, the use of SSAGES is illustrated in the context of simple representative applications involving distinct methods and different collective variables that are available in the current release of the suite. The code may be found at: https://github.com/MICCoM/SSAGES-public.
RESUMEN
We combine path-integral Monte Carlo methods with high-quality potential energy surfaces to compute equilibrium isotope effects in a variety of systems relevant to 'clumped' isotope analysis and isotope geochemistry, including CO2, N2O, methane, and propane. Through a systematic study of heavy-atom isotope-exchange reactions, we quantify and analyze errors that arise in the widely used Urey model for predicting equilibrium constants of isotope-exchange reactions using reduced partition function ratios. These results illustrate that the Urey model relies on a nontrivial cancellation of errors that can shift the apparent equilibrium temperature by as much as 35 K for a given distribution of isotopologues. The calculations reported here provide the same level of precision as the best existing analytical instrumentation, resolving the relative enrichment of certain isotopologues to as little as 0.01. These findings demonstrate path-integral methods to be a rigorous and viable alternative to more approximate methods for heavy-atom geochemical applications.
Asunto(s)
Dióxido de Carbono/química , Metano/química , Óxido Nitroso/química , Propano/química , Isótopos/química , Método de MontecarloRESUMEN
Polymers, with the capacity to tunably alter properties and response based on manipulation of their chemical characteristics, are attractive components in biomaterials. Nevertheless, their potential as functional materials is also inhibited by their complexity, which complicates rational or brute-force design and realization. In recent years, machine learning has emerged as a useful tool for facilitating materials design via efficient modeling of structure-property relationships in the chemical domain of interest. In this Spotlight, we discuss the emergence of data-driven design of polymers that can be deployed in biomaterials with particular emphasis on complex copolymer systems. We outline recent developments, as well as our own contributions and takeaways, related to high-throughput data generation for polymer systems, methods for surrogate modeling by machine learning, and paradigms for property optimization and design. Throughout this discussion, we highlight key aspects of successful strategies and other considerations that will be relevant to the future design of polymer-based biomaterials with target properties.
Asunto(s)
Materiales Biocompatibles , Polímeros , Polímeros/química , Materiales Biocompatibles/química , Aprendizaje Automático , Simulación por ComputadorRESUMEN
We introduce a lattice framework that incorporates elements of Flory-Huggins solution theory and the q-state Potts model to study the phase behavior of polymer solutions and single-chain conformational characteristics. Without empirically introducing temperature-dependent interaction parameters, standard Flory-Huggins theory describes systems that are either homogeneous across temperatures or exhibit upper critical solution temperatures. The proposed Flory-Huggins-Potts framework extends these capabilities by predicting lower critical solution temperatures, miscibility loops, and hourglass-shaped spinodal curves. We particularly show that including orientation-dependent interactions, specifically between monomer segments and solvent particles, is alone sufficient to observe such phase behavior. Signatures of emergent phase behavior are found in single-chain Monte Carlo simulations, which display heating- and cooling-induced coil-globule transitions linked to energy fluctuations. The framework also capably describes a range of experimental systems. Importantly, and by contrast to many prior theoretical approaches, the framework does not employ any temperature- or composition-dependent parameters. This work provides new insights regarding the microscopic physics that underpin complex thermoresponsive behavior in polymers.
RESUMEN
Phase-separated biomolecular condensates exhibit a wide range of dynamic properties, which depend on the sequences of the constituent proteins and RNAs. However, it is unclear to what extent condensate dynamics can be tuned without also changing the thermodynamic properties that govern phase separation. Using coarse-grained simulations of intrinsically disordered proteins, we show that the dynamics and thermodynamics of homopolymer condensates are strongly correlated, with increased condensate stability being coincident with low mobilities and high viscosities. We then apply an "active learning" strategy to identify heteropolymer sequences that break this correlation. This data-driven approach and accompanying analysis reveal how heterogeneous amino acid compositions and nonuniform sequence patterning map to a range of independently tunable dynamic and thermodynamic properties of biomolecular condensates. Our results highlight key molecular determinants governing the physical properties of biomolecular condensates and establish design rules for the development of stimuli-responsive biomaterials.
Asunto(s)
Proteínas Intrínsecamente Desordenadas , Aprendizaje Basado en Problemas , Termodinámica , Aminoácidos , Materiales BiocompatiblesRESUMEN
Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.
RESUMEN
Single-chain nanoparticles (SCNPs) are intriguing materials inspired by proteins that consist of a single precursor polymer chain that has collapsed into a stable structure. In many prospective applications, such as catalysis, the utility of a single-chain nanoparticle will intricately depend on the formation of a mostly specific structure or morphology. However, it is not generally well understood how to reliably control the morphology of single-chain nanoparticles. To address this knowledge gap, we simulate the formation of 7680 distinct single-chain nanoparticles from precursor chains that span a wide range of, in principle, tunable patterning characteristics of cross-linking moieties. Using a combination of molecular simulation and machine learning analyses, we show how the overall fraction of functionalization and blockiness of cross-linking moieties biases the formation of certain local and global morphological characteristics. Importantly, we illustrate and quantify the dispersity of morphologies that arise due to the stochastic nature of collapse from a well-defined sequence as well as from the ensemble of sequences that correspond to a given specification of precursor parameters. Moreover, we also examine the efficacy of precise sequence control in achieving morphological outcomes in different regimes of precursor parameters. Overall, this work critically assesses how precursor chains might be feasibly tailored to achieve given SCNP morphologies and provides a platform to pursue future sequence-based design.
RESUMEN
The wettability of a polymer surfaceârelated to its hydrophobicity or tendency to repel waterâcan be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating the relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon-66, and poly(vinyl alcohol). We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientationâeven when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. These results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.
RESUMEN
Polymer-protein hybrids are intriguing materials that can bolster protein stability in non-native environments, thereby enhancing their utility in diverse medicinal, commercial, and industrial applications. One stabilization strategy involves designing synthetic random copolymers with compositions attuned to the protein surface, but rational design is complicated by the vast chemical and composition space. Here, a strategy is reported to design protein-stabilizing copolymers based on active machine learning, facilitated by automated material synthesis and characterization platforms. The versatility and robustness of the approach is demonstrated by the successful identification of copolymers that preserve, or even enhance, the activity of three chemically distinct enzymes following exposure to thermal denaturing conditions. Although systematic screening results in mixed success, active learning appropriately identifies unique and effective copolymer chemistries for the stabilization of each enzyme. Overall, this work broadens the capabilities to design fit-for-purpose synthetic copolymers that promote or otherwise manipulate protein activity, with extensions toward the design of robust polymer-protein hybrid materials.
Asunto(s)
Polímeros , Procedimientos Quirúrgicos Robotizados , Aprendizaje Automático , Polímeros/química , Proteínas/químicaRESUMEN
Among the many molecules that contribute to glial scarring, chondroitin sulfate proteoglycans (CSPGs) are known to be potent inhibitors of neuronal regeneration. Chondroitinase ABC (ChABC), a bacterial lyase, degrades the glycosaminoglycan (GAG) side chains of CSPGs and promotes tissue regeneration. However, ChABC is thermally unstable and loses all activity within a few hours at 37 °C under dilute conditions. To overcome this limitation, the discovery of a diverse set of tailor-made random copolymers that complex and stabilize ChABC at physiological temperature is reported. The copolymer designs, which are based on chain length and composition of the copolymers, are identified using an active machine learning paradigm, which involves iterative copolymer synthesis, testing for ChABC thermostability upon copolymer complexation, Gaussian process regression modeling, and Bayesian optimization. Copolymers are synthesized by automated PET-RAFT and thermostability of ChABC is assessed by retained enzyme activity (REA) after 24 h at 37 °C. Significant improvements in REA in three iterations of active learning are demonstrated while identifying exceptionally high-performing copolymers. Most remarkably, one designed copolymer promotes residual ChABC activity near 30%, even after one week and notably outperforms other common stabilization methods for ChABC. Together, these results highlight a promising pathway toward sustained tissue regeneration.
Asunto(s)
Condroitina ABC Liasa , Traumatismos de la Médula Espinal , Axones/metabolismo , Teorema de Bayes , Condroitina ABC Liasa/química , Condroitina ABC Liasa/metabolismo , Condroitina ABC Liasa/farmacología , Humanos , Regeneración NerviosaAsunto(s)
Estética , Cara/anatomía & histología , Incisivo , Adolescente , Puntos Anatómicos de Referencia , Cefalometría , Femenino , Humanos , Masculino , Maxilar , Ortodoncia , Valores de ReferenciaRESUMEN
The design of new functional polymers depends on the successful navigation of their structure-function landscapes. Advances in combinatorial polymer chemistry and machine learning provide exciting opportunities for the engineering of fit-for-purpose polymeric materials.
RESUMEN
Polymers are uniquely suited for drug delivery and biomaterial applications due to tunable structural parameters such as length, composition, architecture, and valency. To facilitate designs, researchers may explore combinatorial libraries in a high throughput fashion to correlate structure to function. However, traditional polymerization reactions including controlled living radical polymerization (CLRP) and ring-opening polymerization (ROP) require inert reaction conditions and extensive expertise to implement. With the advent of air-tolerance and automation, several polymerization techniques are now compatible with well plates and can be carried out at the benchtop, making high throughput synthesis and high throughput screening (HTS) possible. To avoid HTS pitfalls often described as "fishing expeditions," it is crucial to employ intelligent and big data approaches to maximize experimental efficiency. This is where the disruptive technologies of machine learning (ML) and artificial intelligence (AI) will likely play a role. In fact, ML and AI are already impacting small molecule drug discovery and showing signs of emerging in drug delivery. In this review, we present state-of-the-art research in drug delivery, gene delivery, antimicrobial polymers, and bioactive polymers alongside data-driven developments in drug design and organic synthesis. From this insight, important lessons are revealed for the polymer therapeutics community including the value of a closed loop design-build-test-learn workflow. This is an exciting time as researchers will gain the ability to fully explore the polymer structural landscape and establish quantitative structure-property relationships (QSPRs) with biological significance.
Asunto(s)
Automatización , Diseño de Fármacos , Polímeros/uso terapéutico , Animales , Ensayos Analíticos de Alto Rendimiento , HumanosRESUMEN
The chemical design of polymers with target structural and/or functional properties represents a grand challenge in materials science. While data-driven design approaches are promising, success with polymers has been limited, largely due to limitations in data availability. Here, we demonstrate the targeted sequence design of single-chain structure in polymers by combining coarse-grained modeling, machine learning, and model optimization. Nearly 2000 unique coarse-grained polymers are simulated to construct and analyze machine learning models. We find that deep neural networks inexpensively and reliably predict structural properties with limited sequence information as input. By coupling trained ML models with sequential model-based optimization, polymer sequences are proposed to exhibit globular, swollen, or rod-like behaviors, which are verified by explicit simulations. This work highlights the promising integration of coarse-grained modeling with data-driven design and represents a necessary and crucial step toward more complex polymer design efforts.
RESUMEN
Knowledge of intrinsic properties is of central importance for materials design and assessing suitability for specific applications. Self-assembling block copolymer electrolytes (BCEs) are of great interest for applications in solid-state energy storage devices. A fundamental understanding of ion transport properties, however, is hindered by the difficulty in deconvoluting extrinsic factors, such as defects, from intrinsic factors, such as the presence of interfaces between the domains. Here, we quantify the intrinsic ion transport properties of a model BCE system consisting of poly(styrene-block-ethylene oxide) (SEO) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt using a generalizable strategy of depositing thin films on interdigitated electrodes and self-assembling fully connected parallel lamellar structures throughout the films. Comparison between conductivity in homopolymer poly(ethylene oxide) (PEO)-LiTFSI electrolytes and the analogous conducting material in SEO over a range of salt concentrations (r, molar ratio of lithium ion to ethylene oxide repeat units) and temperatures reveals that between 20% and 50% of the PEO in SEO is inactive. Using mean-field theory calculations of the domain structure and monomer concentration profiles at domain interfaces-both of which vary substantially with salt concentration-the fraction of inactive PEO in the SEO, as derived from conductivity measurements, can be quantitatively reconciled with the fraction of PEO that is mixed with greater than a few volume percent of polystyrene. Despite the detrimental interfacial effects for ion transport in BCEs, the intrinsic conductivity of the SEO studied here (ca. 10-3 S/cm at 90 °C, r = 0.085) is an order of magnitude higher than reported values from bulk samples of similar molecular weight SEO (ca. 10-4 S/cm at 90 °C, r = 0.085). Overall, this work provides motivation and methods for pursuing improved BCE chemical design, interfacial engineering, and processing.
RESUMEN
A novel methodology is introduced here to generate coarse-grained (CG) representations of molecular models for simulations. The proposed strategy relies on basic graph-theoretic principles and is referred to as graph-based coarse-graining (GBCG). It treats a given system as a molecular graph and derives a corresponding CG representation by using edge contractions to combine nodes in the graph, which correspond to atoms in the molecule, into CG sites. A key element of this methodology is that the nodes are combined according to well-defined protocols that rank-order nodes based on the underlying chemical connectivity. By iteratively performing these operations, successively coarser representations of the original atomic system can be produced to yield a systematic set of CG mappings with hierarchical resolution in an automated fashion. These capabilities are demonstrated in the context of several test systems, including toluene, pentadecane, a polysaccharide dimer, and a rhodopsin protein. In these examples, GBCG yields multiple, intuitive structures that naturally preserve the chemical topology of the system. Importantly, these representations are rendered from algorithmic implementation rather than an arbitrary ansatz, which, until now, has been the conventional approach for defining CG mapping schemes. Overall, the results presented here indicate that GBCG is efficient, robust, and unambiguous in its application, making it a valuable tool for future CG modeling.
RESUMEN
Computational studies aimed at understanding conformationally dependent electronic structure in soft materials require a combination of classical and quantum-mechanical simulations, for which the sampling of conformational space can be particularly demanding. Coarse-grained (CG) models provide a means of accessing relevant time scales, but CG configurations must be back-mapped into atomistic representations to perform quantum-chemical calculations, which is computationally intensive and inconsistent with the spatial resolution of the CG models. A machine learning approach, denoted as artificial neural network electronic coarse graining (ANN-ECG), is presented here in which the conformationally dependent electronic structure of a molecule is mapped directly to CG pseudo-atom configurations. By averaging over decimated degrees of freedom, ANN-ECG accelerates simulations by eliminating backmapping and repeated quantum-chemical calculations. The approach is accurate, consistent with the CG spatial resolution, and can be used to identify computationally optimal CG resolutions.