1
|
Gu M, Wang H, Chen G, Yu A, Dang W, Ling X. Experimental and numerical study on explosion resistance of polyurea-coated shelter in petrochemical industry. Sci Rep 2024; 14:20643. [PMID: 39232177 PMCID: PMC11375104 DOI: 10.1038/s41598-024-71339-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2024] [Accepted: 08/27/2024] [Indexed: 09/06/2024] Open
Abstract
To reduce the number of casualties in explosion accidents, blast-resistant shelters can be used to protect personnel in high-risk areas of petrochemical processing plants. In this work, the deformation behaviours of uncoated and polyurea-coated blast-resistant plates were studied through gas explosion tests. An ANSYS/LS-DYNA model of a polyurea-coated shelter was established, and the dynamic responses of the shelter under various explosion loads were analysed. A series of fuel-air explosion tests were carried out to investigate the explosion resistance of the full-scale shelter. The results showed that compared with the uncoated blast-resistant plate, the deformation of the polyurea-coated blast-resistant plate was significantly reduced. The overall deformation of the shelter was the central depression of the wall and the inward bending of the frame. The damage effect of a typical high-overpressure, low-duration load was greater than that of typical low-overpressure, long-duration load. The shelter remained intact under three repeated explosive loads, with cracks appearing on the inner wall but no collapse or debris splashing. The shock wave attenuation rate of the shelter reached over 90%, which could significantly reduce the number of indoor casualties.
Collapse
Affiliation(s)
- Meng Gu
- State Key Laboratory of Chemical Safety, Qingdao, 266104, China
- SINOPEC Research Institute of Safety Engineering Co., Ltd., Qingdao, 266104, China
- College of Mechanical and Electronic Engineering, China University of Petroleum (East China), Qingdao, 266580, China
| | - Haozhe Wang
- State Key Laboratory of Chemical Safety, Qingdao, 266104, China.
- SINOPEC Research Institute of Safety Engineering Co., Ltd., Qingdao, 266104, China.
| | - Guoxin Chen
- State Key Laboratory of Chemical Safety, Qingdao, 266104, China
- SINOPEC Research Institute of Safety Engineering Co., Ltd., Qingdao, 266104, China
| | - Anfeng Yu
- State Key Laboratory of Chemical Safety, Qingdao, 266104, China
- SINOPEC Research Institute of Safety Engineering Co., Ltd., Qingdao, 266104, China
| | - Wenyi Dang
- State Key Laboratory of Chemical Safety, Qingdao, 266104, China
- SINOPEC Research Institute of Safety Engineering Co., Ltd., Qingdao, 266104, China
| | - Xiaodong Ling
- State Key Laboratory of Chemical Safety, Qingdao, 266104, China
- SINOPEC Research Institute of Safety Engineering Co., Ltd., Qingdao, 266104, China
| |
Collapse
|
2
|
Zhang XZ, Shi R, Lu ZY, Qian HJ. Chemically Specific Systematic Coarse-Grained Polymer Model with Both Consistently Structural and Dynamical Properties. JACS AU 2024; 4:1018-1030. [PMID: 38559727 PMCID: PMC10976574 DOI: 10.1021/jacsau.3c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 02/23/2024] [Accepted: 02/28/2024] [Indexed: 04/04/2024]
Abstract
The coarse-grained (CG) model serves as a powerful tool for the simulation of polymer systems; its reliability depends on the accurate representation of both structural and dynamical properties. However, strong correlations between structural and dynamical properties on different scales and also a strong memory effect, enforced by chain connectivity between monomers in polymer systems, render developing a chemically specific systematic CG model a formidable task. In this study, we report a systematic CG approach that combines the iterative Boltzmann inversion (IBI) method and the generalized Langevin equation (GLE) dynamics. Structural properties are ensured by using conservative CG potentials derived from the IBI method. To retrieve the correct dynamical properties in the system, we demonstrate that using a combination of a Rouse-type delta function and a time-dependent short-time kernel in the GLE simulation is practically efficient. The former can be used to adjust the long-time diffusion dynamics, and the latter can be reconstructed from an iterative procedure according to the velocity autocorrelation function (ACF) from all-atomistic (AA) simulations. Taking the polystyrene as an example, we show that not only structural properties of radial distribution function, intramolecular bond, and angle distributions can be reproduced but also dynamical properties of mean-square displacement, velocity ACF, and force ACF resulted from our CG model have quantitative agreement with the reference AA model. In addition, reasonable agreements are observed in other collective properties between our GLE-CG model and the AA simulations as well.
Collapse
Affiliation(s)
| | | | - Zhong-Yuan Lu
- State Key Laboratory of Supramolecular
Structure and Materials, Institute of Theoretical Chemistry, College
of Chemistry, Jilin University, Changchun 130021, China
| | - Hu-Jun Qian
- State Key Laboratory of Supramolecular
Structure and Materials, Institute of Theoretical Chemistry, College
of Chemistry, Jilin University, Changchun 130021, China
| |
Collapse
|
3
|
DeLuca M, Sensale S, Lin PA, Arya G. Prediction and Control in DNA Nanotechnology. ACS APPLIED BIO MATERIALS 2024; 7:626-645. [PMID: 36880799 DOI: 10.1021/acsabm.2c01045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2023]
Abstract
DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.
Collapse
Affiliation(s)
- Marcello DeLuca
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Sebastian Sensale
- Department of Physics, Cleveland State University, Cleveland, Ohio 44115, United States
| | - Po-An Lin
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Gaurav Arya
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
4
|
Zhang D, Wang Y, Safaripour M, Bellido-Aguilar DA, Van Donselaar KR, Webster DC, Croll AB, Xia W. Energy renormalization for temperature transferable coarse-graining of silicone polymer. Phys Chem Chem Phys 2024; 26:4541-4554. [PMID: 38241021 DOI: 10.1039/d3cp05969c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
The bottom-up prediction of thermodynamic and mechanical behaviors of polymeric materials based on molecular dynamics (MD) simulation is of critical importance in polymer physics. Although the atomistically informed coarse-grained (CG) model can access greater spatiotemporal scales and retain essential chemical specificity, the temperature-transferable CG model is still a big challenge and hinders widespread application of this technique. Herein, we use a silicone polymer, i.e., polydimethylsiloxane (PDMS), having an incredibly low chain rigidity as a model system, combined with an energy-renormalization (ER) approach, to systematically develop a temperature-transferable CG model. Specifically, by introducing temperature-dependent ER factors to renormalize the effective distance and cohesive energy parameters, the developed CG model faithfully preserved the dynamics, mechanical and conformational behaviors compared with the target all-atomistic (AA) model from glassy to melt regimes, which was further validated by experimental data. With the developed CG model featuring tremendously improved computational efficiency, we systematically explored the influences of cohesive interaction strength and temperature on the dynamical heterogeneity and mechanical response of polymers, where we observed consistent trends with other linear polymers with varying chain rigidity and monomeric structures. This study serves as an extension of our proposed ER approach of developing temperature transferable CG models with diverse segmental structures, highlighting the critical role of cohesive interaction strength on CG modeling of polymer dynamics and thermomechanical behaviors.
Collapse
Affiliation(s)
- Dawei Zhang
- Department of Civil, Construction, and Environmental Engineering, North Dakota State University, Fargo, ND 58108, USA
| | - Yang Wang
- Zernike Institute for Advanced Materials, University of Groningen, 9747 AG, Groningen, The Netherlands
| | - Maryam Safaripour
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, ND 58108, USA
| | - Daniel A Bellido-Aguilar
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, ND 58108, USA
| | | | - Dean C Webster
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, ND 58108, USA
| | - Andrew B Croll
- Department of Physics, North Dakota State University, Fargo, ND 58108, USA
| | - Wenjie Xia
- Department of Aerospace Engineering, Iowa State University, Ames, IA 50011, USA.
| |
Collapse
|
5
|
Zong T, Liu X, Zhang X, Yang Q. Efficient characterization of double-cross-linked networks in hydrogels using data-inspired coarse-grained molecular dynamics model. J Chem Phys 2024; 160:024115. [PMID: 38197443 DOI: 10.1063/5.0180847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Accepted: 12/14/2023] [Indexed: 01/11/2024] Open
Abstract
The network structure within polymers significantly influences their mechanical properties, including their strength, toughness, and fatigue resistance. All-atom molecular dynamics (AAMD) simulations offer a method to investigate the energy dissipation mechanism within polymers during deformation and fracture; Such an approach is, however, computationally inefficient when used to analyze polymers with complex network structures, such as the common chemically double-networked hydrogels. Alternatively, coarse-grained molecular dynamics (CGMD) models, which reduce the computational degrees of freedom by concentrating a set of adjacent atoms into a coarse-grained bead, can be employed. In CGMD simulations, a coarse-grained force field (CGFF) is a critical factor affecting the simulation accuracy. In this paper, we proposed a data-based method for predicting the CGFF parameters to improve the simulation efficiency of complex cross-linked network in polymers. Here, we utilized a typical chemically double-networked hydrogel as an example. An artificial neural network was selected, and it was trained with the tensile stress-strain data from the CGMD simulations using different CGFF parameters. The CGMD simulations using the predicted CGFF parameters show good agreement with the AAMD simulations and are almost fifty times faster. The data-inspired CGMD model presented here broadens the applicability of molecular dynamics simulations to cross-linked polymers and has the potential to provide insights that will aid the design of polymers with desirable mechanical properties.
Collapse
Affiliation(s)
- Ting Zong
- Beijing University of Technology, Beijing 100124, China
| | - Xia Liu
- Beijing University of Technology, Beijing 100124, China
| | - Xingyu Zhang
- Beijing University of Technology, Beijing 100124, China
| | | |
Collapse
|
6
|
Wu Z, Wu JW, Michaudel Q, Jayaraman A. Investigating the Hydrogen Bond-Induced Self-Assembly of Polysulfamides Using Molecular Simulations and Experiments. Macromolecules 2023; 56:5033-5049. [PMID: 38362140 PMCID: PMC10865372 DOI: 10.1021/acs.macromol.3c01093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2023] [Revised: 06/08/2023] [Indexed: 02/17/2024]
Abstract
In this paper, we present a synergistic, experimental, and computational study of the self-assembly of N,N'-disubstituted polysulfamides driven by hydrogen bonds (H-bonds) between the H-bonding donor and acceptor groups present in repeating sulfamides as a function of the structural design of the polysulfamide backbone. We developed a coarse-grained (CG) polysulfamide model that captures the directionality of H-bonds between the sulfamide groups and used this model in molecular dynamics (MD) simulations to study the self-assembly of these polymers in implicit solvent. The CGMD approach was validated by reproducing experimentally observed trends in the extent of crystallinity for three polysulfamides synthesized with aliphatic and/or aromatic repeating units. After validation of our CGMD approach, we computationally predicted the effect of repeat unit bulkiness, length, and uniformity of segment lengths in the polymers on the extent of orientational and positional order among the self-assembled polysulfamide chains, providing key design principles for tuning the extent of crystallinity in polysulfamides in experiments. Those computational predictions were then experimentally tested through the synthesis and characterization of polysulfamide architectures.
Collapse
Affiliation(s)
- Zijie Wu
- Department
of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy St., Newark, Delaware 19716, United States
| | - Jiun Wei Wu
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
| | - Quentin Michaudel
- Department
of Chemistry, Texas A&M University, College Station, Texas 77843, United States
- Department
of Materials Science & Engineering, Texas A&M University, College Station, Texas 77843, United States
| | - Arthi Jayaraman
- Department
of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy St., Newark, Delaware 19716, United States
- Department
of Materials Science and Engineering, University
of Delaware, 201 DuPont Hall, Newark, Delaware 19716, United States
| |
Collapse
|
7
|
Ryu WH, Voth GA. Coarse-Graining of Imaginary Time Feynman Path Integrals: Inclusion of Intramolecular Interactions and Bottom-up Force-Matching. J Phys Chem A 2022; 126:6004-6019. [PMID: 36007243 PMCID: PMC9466601 DOI: 10.1021/acs.jpca.2c04349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Revised: 08/16/2022] [Indexed: 11/29/2022]
Abstract
Feynman's imaginary time path integral formalism of quantum statistical mechanics and the corresponding quantum-classical isomorphism provide a tangible way of incorporating nuclear quantum effect (NQE) in the simulation of condensed matter systems using well-developed classical simulation techniques. Our previous work has presented the many-body coarse-graining of path integral (CG-PI) theory that builds an isomorphism between the quantum partition function of N distinguishable particles and the classical partition function of 2N pseudoparticles. In this present work, we develop a generalized version of the many-body CG-PI theory that incorporates many-body interactions in the force field. Based on the new derivation, we provide a numerical CG-PI (n-CG-PI) modeling strategy parametrized from the underlying path integral molecular dynamics (PIMD) trajectories using force matching and Boltzmann inversion. The n-CG-PI models for two liquid systems are shown to capture well both the intramolecular and intermolecular structural correlations of the reference PIMD simulations. The generalized derivation of the many-body CG-PI theory and the n-CG-PI model presented in this work extend the scope of the CG-PI formalism by generalizing the previously limited theory to incorporate force fields of realistic molecular systems.
Collapse
Affiliation(s)
- Won Hee Ryu
- Department of Chemistry,
James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry,
James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
8
|
Polyurea for Blast and Impact Protection: A Review. Polymers (Basel) 2022; 14:polym14132670. [PMID: 35808715 PMCID: PMC9269495 DOI: 10.3390/polym14132670] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2022] [Revised: 06/26/2022] [Accepted: 06/28/2022] [Indexed: 12/10/2022] Open
Abstract
Polyurea has attracted extensive attention from researchers and engineers in the field of blast and impact protection due to its excellent quasi-static mechanical properties and dynamic mechanical properties. Its mechanical properties and energy absorption capacity have been tuned by means of formulation optimization, molecular dynamics (MD) simulation and the addition of reinforcing materials. Owing to the special molecular structure of polyurea, the mechanism of polyurea protection against blasts and impacts is the simultaneous effect of multiple properties. For different substrates and structures, polyurea needs to provide different performance characteristics, including adhesion, hardness, breaking elongation, etc., depending on the characteristics of the load to which it is subjected. The current article reviews relevant publications in the field of polyurea blast and impact protection, including material optimization, protection mechanisms and applications in blast and impact protection.
Collapse
|
9
|
Zheng T, Li T, Shi J, Wu T, Zhuang Z, Xu J, Guo B. Molecular Insight into the Toughness of Polyureas: A Hybrid All-Atom/Coarse-Grained Molecular Dynamics Study. Macromolecules 2022. [DOI: 10.1021/acs.macromol.1c02453] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tianze Zheng
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Ting Li
- Aerospace Research Institute of Materials and Processing Technology, Beijing 100076, China
| | - Jiaxin Shi
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Tianyu Wu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Zhuo Zhuang
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, China
| | - Jun Xu
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Baohua Guo
- Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
10
|
Dhamankar S, Webb MA. Chemically specific coarse‐graining of polymers: Methods and prospects. JOURNAL OF POLYMER SCIENCE 2021. [DOI: 10.1002/pol.20210555] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Satyen Dhamankar
- Department of Chemical and Biological Engineering Princeton University Princeton New Jersey USA
| | - Michael A. Webb
- Department of Chemical and Biological Engineering Princeton University Princeton New Jersey USA
| |
Collapse
|
11
|
Zheng T, Zhang Y, Shi J, Xu J, Guo B. Revealing the role of hydrogen bonding in polyurea with multiscale simulations. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1967346] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Tianze Zheng
- Advanced Materials Laboratory of Ministry of Education (MOE), Department of Chemical Engineering, Tsinghua University, Beijing, People’s Republic of China
| | - Yao Zhang
- Advanced Materials Laboratory of Ministry of Education (MOE), Department of Chemical Engineering, Tsinghua University, Beijing, People’s Republic of China
| | - Jiaxin Shi
- Advanced Materials Laboratory of Ministry of Education (MOE), Department of Chemical Engineering, Tsinghua University, Beijing, People’s Republic of China
| | - Jun Xu
- Advanced Materials Laboratory of Ministry of Education (MOE), Department of Chemical Engineering, Tsinghua University, Beijing, People’s Republic of China
| | - Baohua Guo
- Advanced Materials Laboratory of Ministry of Education (MOE), Department of Chemical Engineering, Tsinghua University, Beijing, People’s Republic of China
| |
Collapse
|
12
|
Wu C, Li K, Ning X, Zhang L. An Enhanced Scheme for Multiscale Modeling of Thermomechanical Properties of Polymer Bulks. J Phys Chem B 2021; 125:8612-8626. [PMID: 34291641 DOI: 10.1021/acs.jpcb.1c02663] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
While multiscale modeling significantly enhances the capability of molecular simulations of polymer systems, it is well realized that the systematically derived coarse-grained (CG) models generally underestimate the thermomechanical properties. In this work, a charge-based mapping scheme has been adopted to include explicit electrostatic interactions and benchmarked against two typical polymers, atactic poly(methyl methacrylate) (PMMA) and polystyrene (PS). The CG potentials are parameterized against the oligomer bulks of nine monomers per chain to match the essential structural features and the two basic pressure-volume-temperature (PVT) properties, which are obtained from the all-atomistic (AA) molecular dynamics (MD) simulations at a single elevated temperature. The so-parameterized CG potentials are extended with the MD method to simulate the two polymer bulks of one hundred monomers per chain over a wide temperature range. Without any scaling, all the simulated results, including mass densities and bulk moduli at room temperature, thermal expansion coefficients at rubbery and glassy states, and glass transition temperatures (Tg), compare well with the corresponding experimental data. The proposed scheme not only contributes to realistically simulating various thermomechanical properties of both apolar and polar polymers but also allows for directly simulating their electrical properties.
Collapse
Affiliation(s)
- Chaofu Wu
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| | - Kewen Li
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| | - Xutao Ning
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| | - Lei Zhang
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials, School of Materials and Environmental Engineering, Hunan University of Humanities, Science and Technology, Loudi, Hunan 417000, P. R. China
| |
Collapse
|
13
|
Zhang Z, Krajniak J, Ganesan V. A Multiscale Simulation Study of Influence of Morphology on Ion Transport in Block Copolymeric Ionic Liquids. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00025] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Zidan Zhang
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Jakub Krajniak
- Independent researcher, os. Kosmonautow 13/56, 61-631 Poznan, Poland
| | - Venkat Ganesan
- McKetta Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
14
|
Mousavi AA, Arash B, Rolfes R. Optimization assisted coarse-grained modeling of agglomerated nanoparticle reinforced thermosetting polymers. POLYMER 2021. [DOI: 10.1016/j.polymer.2021.123741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
15
|
Johnson LC, Phelan FR. Dynamically consistent coarse-grain simulation model of chemically specific polymer melts via friction parameterization. J Chem Phys 2021; 154:084114. [PMID: 33639746 PMCID: PMC10075510 DOI: 10.1063/5.0034910] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Coarse-grained (CG) models of polymers involve grouping many atoms in an all-atom (AA) representation into single sites to reduce computational effort yet retain the hierarchy of length and time scales inherent to macromolecules. Parameterization of such models is often via "bottom-up" methods, which preserve chemical specificity but suffer from artificially accelerated dynamics with respect to the AA model from which they were derived. Here, we study the combination of a bottom-up CG model with a dissipative potential as a means to obtain a chemically specific and dynamically correct model. We generate the conservative part of the force-field using the iterative Boltzmann inversion (IBI) method, which seeks to recover the AA structure. This is augmented with the dissipative Langevin thermostat, which introduces a single parameterizable friction factor to correct the unphysically fast dynamics of the IBI-generated force-field. We study this approach for linear polystyrene oligomer melts for three separate systems with 11, 21, and 41 monomers per chain and a mapping of one monomer per CG site. To parameterize the friction factor, target values are extracted from the AA dynamics using translational monomer diffusion, translational chain diffusion, and rotational chain motion to test the consistency of the parameterization across different modes of motion. We find that the value of the friction parameter needed to bring the CG dynamics in line with AA target values varies based on the mode of parameterization with short-time monomer translational dynamics requiring the highest values, long-time chain translational dynamics requiring the lowest values, and rotational dynamics falling in between. The friction ranges most widely for the shortest chains, and the span narrows with increasing chain length. For longer chains, a practical working value of the friction parameter may be derived from the rotational dynamics, owing to the contribution of multiple relaxation modes to chain rotation and a lack of sensitivity of the translational dynamics at these intermediate levels of friction. A study of equilibrium chain structure reveals that all chains studied are non-Gaussian. However, longer chains better approximate ideal chain dimensions than more rod-like shorter chains and thus are most closely described by a single friction parameter. We also find that the separability of the conservative and dissipative potentials is preserved.
Collapse
Affiliation(s)
- Lilian C Johnson
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Frederick R Phelan
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
16
|
Effectiveness of coarse graining degree and speedup on the dynamic properties of homopolymer. J Mol Model 2021; 27:55. [PMID: 33511476 DOI: 10.1007/s00894-020-04661-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Accepted: 12/20/2020] [Indexed: 10/22/2022]
Abstract
Evaluation of effective coarse graining (CG) degree and reasonable speedup relative to all-atomistic (AA) model was conducted to provide a basis for building appropriate larger-scale model. The reproducibility of atomistic conformation and temperature transferability both act as the analysis criteria to resolve the maximum acceptable CG degree. Taking short- and long time spans into account simultaneously in the estimation of computational speedup, a dynamic scaling factor is accessible in fitting mean squared displacement ratio of CG to AA as an exponential function. Computing loss in parallel running is an indispensable component in acceleration, which was also added in the evaluation. Subsequently, a quantified prediction of CG speedup arises as a multiplication of dynamic scaling factor, computing loss, time step, and the square of reduction in the number of degrees of freedom. Polyethylene oxide was adopted as a reference system to execute the direct Boltzmann inversion and iterative Boltzmann inversion. Bonded and non-bonded potentials were calculated in CG models with 1~4 monomers per bead. The effective CG degree was determined as two at the most with a speedup of four orders magnitude over AA in this study. Determination of effectiveness CG degree and the corresponding speedup prediction provide available tools in larger spatiotemporal-scale calculations.
Collapse
|
17
|
Eastmond T, Hu J, Alizadeh V, Hrubiak R, Oswald J, Amirkhizi A, Peralta P. Probing High-Pressure Structural Evolution in Polyurea with In Situ Energy-Dispersive X-ray Diffraction and Molecular Dynamics Simulations. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c02266] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tyler Eastmond
- Mechanical and Aerospace Engineering, Arizona State University, 501 E Tyler Mall, Tempe, Arizona 85287, United States
| | - Jing Hu
- Mechanical and Aerospace Engineering, Arizona State University, 501 E Tyler Mall, Tempe, Arizona 85287, United States
| | - Vahidreza Alizadeh
- Department of Mechanical Engineering, University of Massachusetts Lowell, Dandeneau Hall
219, 1 University Ave, Lowell, Massachusetts 01854, United States
| | - Rostislav Hrubiak
- High Pressure Collaborative Access Team, X-ray Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Jay Oswald
- Mechanical and Aerospace Engineering, Arizona State University, 501 E Tyler Mall, Tempe, Arizona 85287, United States
| | - Alireza Amirkhizi
- Department of Mechanical Engineering, University of Massachusetts Lowell, Dandeneau Hall
219, 1 University Ave, Lowell, Massachusetts 01854, United States
| | - Pedro Peralta
- Mechanical and Aerospace Engineering, Arizona State University, 501 E Tyler Mall, Tempe, Arizona 85287, United States
| |
Collapse
|
18
|
Yoshimoto Y, Sugiyama S, Shimada S, Kaneko T, Takagi S, Kinefuchi I. Molecular Insights into the Mechanical Properties of Polymer–Fullerene Bulk Heterojunctions for Organic Photovoltaic Applications. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c02278] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yuta Yoshimoto
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Sou Sugiyama
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Shuntaro Shimada
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Toshihiro Kaneko
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Shu Takagi
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Ikuya Kinefuchi
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| |
Collapse
|
19
|
Park S, Moon J, Kim B, Cho M. Multi-scale coarse-grained molecular dynamics simulation to investigate the thermo-mechanical behavior of shape-memory polyurethane copolymers. POLYMER 2021. [DOI: 10.1016/j.polymer.2020.123228] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
20
|
Ugarte La Torre D, Takada S. Coarse-grained implicit solvent lipid force field with a compatible resolution to the Cα protein representation. J Chem Phys 2020; 153:205101. [PMID: 33261497 DOI: 10.1063/5.0026342] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Biological membranes have been prominent targets for coarse-grained (CG) molecular dynamics simulations. While minimal CG lipid models with three beads per lipid and quantitative CG lipid models with >10 beads per lipid have been well studied, in between them, CG lipid models with a compatible resolution to residue-level CG protein models are much less developed. Here, we extended a previously developed three-bead lipid model into a five-bead model and parameterized it for two phospholipids, POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine). The developed model, iSoLF, reproduced the area per lipid, hydrophobic thickness, and phase behaviors of the target phospholipid bilayer membranes at the physiological temperature. The model POPC and DPPC membranes were in liquid and gel phases, respectively, in accordance with experiments. We further examined the spontaneous formation of a membrane bilayer, the temperature dependence of physical properties, the vesicle dynamics, and the POPC/DPPC two-component membrane dynamics of the CG lipid model, showing some promise. Once combined with standard Cα protein models, the iSoLF model will be a powerful tool to simulate large biological membrane systems made of lipids and proteins.
Collapse
Affiliation(s)
- Diego Ugarte La Torre
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| |
Collapse
|
21
|
Kapoor U, Kulshreshtha A, Jayaraman A. Development of Coarse-Grained Models for Poly(4-vinylphenol) and Poly(2-vinylpyridine): Polymer Chemistries with Hydrogen Bonding. Polymers (Basel) 2020; 12:E2764. [PMID: 33238611 PMCID: PMC7709027 DOI: 10.3390/polym12112764] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2020] [Revised: 11/06/2020] [Accepted: 11/09/2020] [Indexed: 11/16/2022] Open
Abstract
In this paper, we identify the modifications needed in a recently developed generic coarse-grained (CG) model that captured directional interactions in polymers to specifically represent two exemplary hydrogen bonding polymer chemistries-poly(4-vinylphenol) and poly(2-vinylpyridine). We use atomistically observed monomer-level structures (e.g., bond, angle and torsion distribution) and chain structures (e.g., end-to-end distance distribution and persistence length) of poly(4-vinylphenol) and poly(2-vinylpyridine) in an explicitly represented good solvent (tetrahydrofuran) to identify the appropriate modifications in the generic CG model in implicit solvent. For both chemistries, the modified CG model is developed based on atomistic simulations of a single 24-mer chain. This modified CG model is then used to simulate longer (36-mer) and shorter (18-mer and 12-mer) chain lengths and compared against the corresponding atomistic simulation results. We find that with one to two simple modifications (e.g., incorporating intra-chain attraction, torsional constraint) to the generic CG model, we are able to reproduce atomistically observed bond, angle and torsion distributions, persistence length, and end-to-end distance distribution for chain lengths ranging from 12 to 36 monomers. We also show that this modified CG model, meant to reproduce atomistic structure, does not reproduce atomistically observed chain relaxation and hydrogen bond dynamics, as expected. Simulations with the modified CG model have significantly faster chain relaxation than atomistic simulations and slower decorrelation of formed hydrogen bonds than in atomistic simulations, with no apparent dependence on chain length.
Collapse
Affiliation(s)
- Utkarsh Kapoor
- Department of Chemical and Biomolecular Engineering, Colburn Laboratory, University of Delaware, 150 Academy Street, Newark, DE 19716, USA; (U.K.); (A.K.)
| | - Arjita Kulshreshtha
- Department of Chemical and Biomolecular Engineering, Colburn Laboratory, University of Delaware, 150 Academy Street, Newark, DE 19716, USA; (U.K.); (A.K.)
| | - Arthi Jayaraman
- Department of Chemical and Biomolecular Engineering, Colburn Laboratory, University of Delaware, 150 Academy Street, Newark, DE 19716, USA; (U.K.); (A.K.)
- Department of Materials Science and Engineering, University of Delaware, Newark, DE 19716, USA
| |
Collapse
|
22
|
Wang KW, Wang Y, Hall CK. Development of a coarse-grained lipid model, LIME 2.0, for DSPE using multistate iterative Boltzmann inversion and discontinuous molecular dynamics simulations. FLUID PHASE EQUILIBRIA 2020; 521:112704. [PMID: 37982069 PMCID: PMC10655612 DOI: 10.1016/j.fluid.2020.112704] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
We suggest an improved version of the intermediate resolution implicit solvent model for lipids, LIME, that was previously developed for use with discontinuous molecular dynamics (DMD) simulations. LIME gets its geometrical and the energy parameters between bonded and nonbonded pairs of coarse-grained (CG) sites from atomistic simulations. The improved model, LIME 2.0, uses multiple square wells rather than the single square well used in original LIME to obtain intermolecular interactions that more faithfully mimic those from atomistic simulations. The multi-state iterative Boltzmann inversion (MS-IBI) scheme is used to determine the interaction parameters. This means that a single set of interaction parameters between coarse-grained sites can be used to represent the lipid bilayers at different temperatures. The physical properties of CG DSPE lipid bilayer are calculated using CG simulations and compared to atomistic simulations results to verify the improved model. The phase transition temperature of the lipid bilayer is measured accurately and the lipid translocation phenomenon, " flip-flop" is observed through CG simulation. These results suggest that CG parameterization using multiple square-well and the MS-IBI scheme is well suited to the study of lipid bilayers cross a range of temperatures with DMD simulations.
Collapse
Affiliation(s)
- Kye Won Wang
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina, 27695
| | - Yiming Wang
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina, 27695
| | - Carol K. Hall
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina, 27695
| |
Collapse
|
23
|
Joshi SY, Deshmukh SA. A review of advancements in coarse-grained molecular dynamics simulations. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1828583] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Soumil Y. Joshi
- Department of Chemical Engineering, Virginia Tech, Blacksburg, VA, USA
| | | |
Collapse
|
24
|
Wang S, Ma Z, Pan W. Data-driven coarse-grained modeling of polymers in solution with structural and dynamic properties conserved. SOFT MATTER 2020; 16:8330-8344. [PMID: 32785383 DOI: 10.1039/d0sm01019g] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present data-driven coarse-grained (CG) modeling for polymers in solution, which conserves the dynamic as well as structural properties of the underlying atomistic system. The CG modeling is built upon the framework of the generalized Langevin equation (GLE). The key is to determine each term in the GLE by directly linking it to atomistic data. In particular, we propose a two-stage Gaussian process-based Bayesian optimization method to infer the non-Markovian memory kernel from the data of the velocity autocorrelation function (VACF). Considering that the long-time behaviors of the VACF and memory kernel for polymer solutions can exhibit hydrodynamic scaling (algebraic decay with time), we further develop an active learning method to determine the emergence of hydrodynamic scaling, which can accelerate the inference process of the memory kernel. The proposed methods do not rely on how the mean force or CG potential in the GLE is constructed. Thus, we also compare two methods for constructing the CG potential: a deep learning method and the iterative Boltzmann inversion method. With the memory kernel and CG potential determined, the GLE is mapped onto an extended Markovian process to circumvent the expensive cost of directly solving the GLE. The accuracy and computational efficiency of the proposed CG modeling are assessed in a model star-polymer solution system at three representative concentrations. By comparing with the reference atomistic simulation results, we demonstrate that the proposed CG modeling can robustly and accurately reproduce the dynamic and structural properties of polymers in solution.
Collapse
Affiliation(s)
- Shu Wang
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| | - Zhan Ma
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| | - Wenxiao Pan
- Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| |
Collapse
|
25
|
Li Q, Guo Y, Tong J, He H, Zhang X, Huo F. Development of a coarse-grained force field model of polymeric 1-vinyl-3-ethylimidazolium tetrafluoroborate ionic liquids. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137656] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
26
|
|
27
|
|
28
|
Jayaraman A. 100th Anniversary of Macromolecular Science Viewpoint: Modeling and Simulation of Macromolecules with Hydrogen Bonds: Challenges, Successes, and Opportunities. ACS Macro Lett 2020; 9:656-665. [PMID: 35648569 DOI: 10.1021/acsmacrolett.0c00134] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Macromolecular materials with directional interactions such as hydrogen bonds exhibit numerous attractive features in terms of structure, thermodynamics, and dynamics. Besides enabling precise tuning of desirable geometries in the assembled state (e.g., programmable coordination numbers depending on the valency of the directional interaction), mixing in a blend/composite through stabilization via hydrogen bonds between the various components, hydrogen bonds can also impart responsiveness to external stimuli (e.g., temperature, pH). In biomacromolecules (e.g., proteins, DNA, polysaccharides), hydrogen bonds play a key role in stabilizing secondary and tertiary structures, which in turn define the function of these macromolecules. In this Viewpoint, I present the challenges, successes, and opportunities for molecular modeling and simulations to conduct fundamental and application-focused research on macromolecular materials with hydrogen bonding interactions. The past successes and limitations of atomistic simulations are discussed first, followed by highlights from recent developments in coarse-grained modeling and their use in studies of (synthetic and biologically relevant) macromolecular materials. Model development focused on polynucleotides (e.g., DNA, RNA, etc.), polypeptides, polysaccharides, and synthetic polymers at experimentally relevant conditions are highlighted. This viewpoint ends with potential future directions for macromolecular modeling and simulations with other types of directional interactions beyond hydrogen bonding.
Collapse
Affiliation(s)
- Arthi Jayaraman
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
- Department of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
29
|
Ransom TC, Gamache RM, Mason BP, Ladouceur HD, Roland CM. DEFORMATION INDEX APPLIED TO IMPACT. RUBBER CHEMISTRY AND TECHNOLOGY 2020. [DOI: 10.5254/rct.20.80362] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
ABSTRACT
Almost three decades ago, S. Futamura devised the deformation index concept for determining the control parameter for the viscoelastic response of deformed elastomers. We have extended this concept to impact mitigation, wherein material hardness and energy dissipation typically both affect the behavior. Laboratory impact tests were carried out on a series of compounds to deduce the deformation index pertinent to the rubber component. We then analyzed ballistic experiments, wherein material failure is associated with more complex conditions. The utility and limitations of this approach are discussed.
Collapse
Affiliation(s)
- T. C. Ransom
- Naval Research Laboratory, Chemistry Division, Code 6105, Washington, D.C. 20375-5342
| | - R. M. Gamache
- Naval Postgraduate School, Department of Physics, Monterey, CA 93943-5216
| | - B. P. Mason
- Naval Postgraduate School, Department of Physics, Monterey, CA 93943-5216
| | - H. D. Ladouceur
- Naval Research Laboratory, Chemistry Division, Code 6105, Washington, D.C. 20375-5342
| | - C. M. Roland
- Naval Research Laboratory, Chemistry Division, Code 6105, Washington, D.C. 20375-5342
| |
Collapse
|
30
|
Shahidi N, Chazirakis A, Harmandaris V, Doxastakis M. Coarse-graining of polyisoprene melts using inverse Monte Carlo and local density potentials. J Chem Phys 2020; 152:124902. [PMID: 32241142 DOI: 10.1063/1.5143245] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Bottom-up coarse-graining of polymers is commonly performed by matching structural order parameters such as distribution of bond lengths, bending and dihedral angles, and pair distribution functions. In this study, we introduce the distribution of nearest-neighbors as an additional order parameter in the concept of local density potentials. We describe how the inverse-Monte Carlo method provides a framework for forcefield development that is capable of overcoming challenges associated with the parameterization of interaction terms in polymer systems. The technique is applied on polyisoprene melts as a prototype system. We demonstrate that while different forcefields can be developed that perform equally in terms of matching target distributions, the inclusion of nearest-neighbors provides a straightforward route to match both thermodynamic and conformational properties. We find that several temperature state points can also be addressed, provided that the forcefield is refined accordingly. Finally, we examine both the single-particle and the collective dynamics of the coarse-grain models, demonstrating that all forcefields present a similar acceleration relative to the atomistic systems.
Collapse
Affiliation(s)
- Nobahar Shahidi
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Antonis Chazirakis
- Department of Mathematics and Applied Mathematics, University of Crete, Heraklion GR-71110, Greece
| | - Vagelis Harmandaris
- Department of Mathematics and Applied Mathematics, University of Crete, Heraklion GR-71110, Greece
| | - Manolis Doxastakis
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| |
Collapse
|
31
|
Bartnik K, Barth A, Pilo-Pais M, Crevenna AH, Liedl T, Lamb DC. A DNA Origami Platform for Single-Pair Förster Resonance Energy Transfer Investigation of DNA-DNA Interactions and Ligation. J Am Chem Soc 2020; 142:815-825. [PMID: 31800234 DOI: 10.1021/jacs.9b09093] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
DNA double-strand breaks (DSBs) pose an everyday threat to the conservation of genetic information and therefore life itself. Several pathways have evolved to repair these cytotoxic lesions by rejoining broken ends, among them the nonhomologous end-joining mechanism that utilizes a DNA ligase. Here, we use a custom-designed DNA origami nanostructure as a model system to specifically mimic a DNA DSB, enabling us to study the end-joining of two fluorescently labeled DNA with the T4 DNA ligase on the single-molecule level. The ligation reaction is monitored by Förster resonance energy transfer (FRET) experiments both in solution and on surface-anchored origamis. Due to the modularity of DNA nanotechnology, DNA double strands with different complementary overhang lengths can be studied using the same DNA origami design. We show that the T4 DNA ligase repairs sticky ends more efficiently than blunt ends and that the ligation efficiency is influenced by both DNA sequence and the incubation conditions. Before ligation, dynamic fluctuations of the FRET signal are observed due to transient binding of the sticky overhangs. After ligation, the FRET signal becomes static. Thus, we can directly monitor the ligation reaction through the transition from dynamic to static FRET signals. Finally, we revert the ligation process using a restriction enzyme digestion and religate the resulting blunt ends. The here-presented DNA origami platform is thus suited to study complex multistep reactions occurring over several cycles of enzymatic treatment.
Collapse
Affiliation(s)
- Kira Bartnik
- Department of Chemistry, Center for Nanoscience (CeNS), Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CIPSM) , Ludwig-Maximilians-Universität München , 81377 Munich , Germany
| | - Anders Barth
- Department of Chemistry, Center for Nanoscience (CeNS), Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CIPSM) , Ludwig-Maximilians-Universität München , 81377 Munich , Germany
| | - Mauricio Pilo-Pais
- Department of Physics and Center for Nanoscience (CeNS) , Ludwig-Maximilians-Universität , 80539 Munich , Germany
| | - Alvaro H Crevenna
- Department of Chemistry, Center for Nanoscience (CeNS), Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CIPSM) , Ludwig-Maximilians-Universität München , 81377 Munich , Germany
| | - Tim Liedl
- Department of Physics and Center for Nanoscience (CeNS) , Ludwig-Maximilians-Universität , 80539 Munich , Germany
| | - Don C Lamb
- Department of Chemistry, Center for Nanoscience (CeNS), Nanosystems Initiative Munich (NIM) and Center for Integrated Protein Science Munich (CIPSM) , Ludwig-Maximilians-Universität München , 81377 Munich , Germany
| |
Collapse
|
32
|
Liu M, Oswald J. Coarse–grained molecular modeling of the microphase structure of polyurea elastomer. POLYMER 2019. [DOI: 10.1016/j.polymer.2019.04.039] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
33
|
Hu C, Lu T, Guo H. Developing a Transferable Coarse-Grained Model for the Prediction of Thermodynamic, Structural, and Mechanical Properties of Polyimides at Different Thermodynamic State Points. J Chem Inf Model 2019; 59:2009-2025. [DOI: 10.1021/acs.jcim.8b00887] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Chenchen Hu
- Beijing National Laboratory for Molecular Sciences, Joint Laboratory of Polymer Sciences and Materials, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Teng Lu
- Beijing National Laboratory for Molecular Sciences, Joint Laboratory of Polymer Sciences and Materials, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Hongxia Guo
- Beijing National Laboratory for Molecular Sciences, Joint Laboratory of Polymer Sciences and Materials, State Key Laboratory of Polymer Physics and Chemistry, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
34
|
Guenza MG, Dinpajooh M, McCarty J, Lyubimov IY. Accuracy, Transferability, and Efficiency of Coarse-Grained Models of Molecular Liquids. J Phys Chem B 2018; 122:10257-10278. [DOI: 10.1021/acs.jpcb.8b06687] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- M. G. Guenza
- Department of Chemistry and Biochemistry and Institute of Theoretical Science, University of Oregon, Eugene, Oregon 97403, United States
| | - M. Dinpajooh
- Department of Chemistry and Biochemistry and Institute of Theoretical Science, University of Oregon, Eugene, Oregon 97403, United States
| | - J. McCarty
- Department of Chemistry and Biochemistry and Institute of Theoretical Science, University of Oregon, Eugene, Oregon 97403, United States
| | - I. Y. Lyubimov
- Department of Chemistry and Biochemistry and Institute of Theoretical Science, University of Oregon, Eugene, Oregon 97403, United States
| |
Collapse
|
35
|
Chakraborty D, Hori N, Thirumalai D. Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA. J Chem Theory Comput 2018; 14:3763-3779. [PMID: 29870236 PMCID: PMC6423546 DOI: 10.1021/acs.jctc.8b00091] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.
Collapse
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Naoto Hori
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - D. Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
36
|
Zadok I, Srebnik S. Coarse-Grained Simulation of Protein-Imprinted Hydrogels. J Phys Chem B 2018; 122:7091-7101. [DOI: 10.1021/acs.jpcb.8b03774] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Israel Zadok
- Department of Chemical Engineering, Technion—Israel Institute of Technology, Haifa 32000, Israel
| | - Simcha Srebnik
- Department of Chemical Engineering, Technion—Israel Institute of Technology, Haifa 32000, Israel
| |
Collapse
|
37
|
Yoshimoto Y, Li Z, Kinefuchi I, Karniadakis GE. Construction of non-Markovian coarse-grained models employing the Mori–Zwanzig formalism and iterative Boltzmann inversion. J Chem Phys 2017; 147:244110. [DOI: 10.1063/1.5009041] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yuta Yoshimoto
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Zhen Li
- Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA
| | - Ikuya Kinefuchi
- Department of Mechanical Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - George Em Karniadakis
- Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
38
|
Wu C. Multiscale Modeling Scheme for Simulating Polymeric Melts: Application to Poly(Ethylene Oxide). MACROMOL THEOR SIMUL 2017. [DOI: 10.1002/mats.201700066] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- C. Wu
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials; Hunan University of Humanities Science & Technology; Dixing Road 487, Louxing District Loudi 417000 Hunan Province P. R. China
- College of Materials and Environment Engineering; Hunan University of Humanities Science & Technology; Dixing Road 487, Louxing District Loudi 417000 Hunan Province P. R. China
| |
Collapse
|
39
|
Kacar G. Dissipative particle dynamics parameterization and simulations to predict negative volume excess and structure of PEG and water mixtures. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.10.052] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
40
|
Dinpajooh M, Guenza MG. On the Density Dependence of the Integral Equation Coarse-Graining Effective Potential. J Phys Chem B 2017; 122:3426-3440. [DOI: 10.1021/acs.jpcb.7b10494] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mohammadhasan Dinpajooh
- Department of Chemistry and Biochemistry, and Institute of Theoretical Science, University of Oregon, Eugene, Oregon 97403, United States
| | - Marina G. Guenza
- Department of Chemistry and Biochemistry, and Institute of Theoretical Science, University of Oregon, Eugene, Oregon 97403, United States
| |
Collapse
|
41
|
Ransom TC, Ahart M, Hemley RJ, Roland CM. Vitrification and Density Scaling of Polyurea at Pressures up to 6 GPa. Macromolecules 2017. [DOI: 10.1021/acs.macromol.7b01676] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Timothy C. Ransom
- Chemistry
Division, Naval Research Laboratory, Code 6105, Washington, D.C. 20375-53452, United States
| | - Muhtar Ahart
- Geophysical
Laboratory, Carnegie Institution of Washington, Washington, D.C. 20015, United States
| | - Russell J. Hemley
- Department
of Civil and Environmental Engineering, The George Washington University, Washington, D.C. 20052, United States
| | - C. Michael Roland
- Chemistry
Division, Naval Research Laboratory, Code 6105, Washington, D.C. 20375-53452, United States
| |
Collapse
|
42
|
Kuo AT, Okazaki S, Shinoda W. Transferable coarse-grained model for perfluorosulfonic acid polymer membranes. J Chem Phys 2017; 147:094904. [DOI: 10.1063/1.4986287] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Shi Z, Castro CE, Arya G. Conformational Dynamics of Mechanically Compliant DNA Nanostructures from Coarse-Grained Molecular Dynamics Simulations. ACS NANO 2017; 11:4617-4630. [PMID: 28423273 DOI: 10.1021/acsnano.7b00242] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Structural DNA nanotechnology, the assembly of rigid 3D structures of complex yet precise geometries, has recently been used to design dynamic, mechanically compliant nanostructures with tunable equilibrium conformations and conformational distributions. Here we use coarse-grained molecular dynamics simulations to provide insights into the conformational dynamics of a set of mechanically compliant DNA nanostructures-DNA hinges that use single-stranded DNA "springs" to tune the equilibrium conformation of a layered double-stranded DNA "joint" connecting two stiff "arms" constructed from DNA helix bundles. The simulations reproduce the experimentally measured equilibrium angles between hinge arms for a range of hinge designs. The hinges are found to be structurally stable, except for some fraying of the open ends of the DNA helices comprising the hinge arms and some loss of base-pairing interactions in the joint regions coinciding with the crossover junctions, especially in hinges designed to exhibit a small bending angle that exhibit large local stresses resulting in strong kinks in their joints. Principal component analysis reveals that while the hinge dynamics are dominated by bending motion, some twisting and sliding of hinge arms relative to each other also exists. Forced deformation of the hinges reveals distinct bending mechanisms for hinges with short, inextensible springs versus those with longer, more extensible springs. Lastly, we introduce an approach for rapidly predicting equilibrium hinge angles from individual force-deformation behaviors of its single- and double-stranded DNA components. Taken together, these results demonstrate that coarse-grained modeling is a promising approach for designing, predicting, and studying the dynamics of compliant DNA nanostructures, where conformational fluctuations become important, multiple deformation mechanisms exist, and continuum approaches may not yield accurate properties.
Collapse
Affiliation(s)
- Ze Shi
- Department of NanoEngineering, University of California, San Diego , La Jolla, California 92093, United States
| | - Carlos E Castro
- Department of Mechanical and Aerospace Engineering, The Ohio State University , Columbus, Ohio 43210, United States
| | - Gaurav Arya
- Department of NanoEngineering, University of California, San Diego , La Jolla, California 92093, United States
| |
Collapse
|
44
|
Root SE, Savagatrup S, Printz AD, Rodriquez D, Lipomi DJ. Mechanical Properties of Organic Semiconductors for Stretchable, Highly Flexible, and Mechanically Robust Electronics. Chem Rev 2017; 117:6467-6499. [DOI: 10.1021/acs.chemrev.7b00003] [Citation(s) in RCA: 465] [Impact Index Per Article: 66.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Samuel E. Root
- Department of NanoEngineering, University of California, San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Suchol Savagatrup
- Department of NanoEngineering, University of California, San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Adam D. Printz
- Department of NanoEngineering, University of California, San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Daniel Rodriquez
- Department of NanoEngineering, University of California, San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Darren J. Lipomi
- Department of NanoEngineering, University of California, San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| |
Collapse
|
45
|
Li Z, Bian X, Yang X, Karniadakis GE. A comparative study of coarse-graining methods for polymeric fluids: Mori-Zwanzig vs. iterative Boltzmann inversion vs. stochastic parametric optimization. J Chem Phys 2017; 145:044102. [PMID: 27475343 DOI: 10.1063/1.4959121] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We construct effective coarse-grained (CG) models for polymeric fluids by employing two coarse-graining strategies. The first one is a forward-coarse-graining procedure by the Mori-Zwanzig (MZ) projection while the other one applies a reverse-coarse-graining procedure, such as the iterative Boltzmann inversion (IBI) and the stochastic parametric optimization (SPO). More specifically, we perform molecular dynamics (MD) simulations of star polymer melts to provide the atomistic fields to be coarse-grained. Each molecule of a star polymer with internal degrees of freedom is coarsened into a single CG particle and the effective interactions between CG particles can be either evaluated directly from microscopic dynamics based on the MZ formalism, or obtained by the reverse methods, i.e., IBI and SPO. The forward procedure has no free parameters to tune and recovers the MD system faithfully. For the reverse procedure, we find that the parameters in CG models cannot be selected arbitrarily. If the free parameters are properly defined, the reverse CG procedure also yields an accurate effective potential. Moreover, we explain how an aggressive coarse-graining procedure introduces the many-body effect, which makes the pairwise potential invalid for the same system at densities away from the training point. From this work, general guidelines for coarse-graining of polymeric fluids can be drawn.
Collapse
Affiliation(s)
- Zhen Li
- Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA
| | - Xin Bian
- Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA
| | - Xiu Yang
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - George Em Karniadakis
- Division of Applied Mathematics, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
46
|
All-atom molecular dynamics simulations of multiphase segregated polyurea under quasistatic, adiabatic, uniaxial compression. POLYMER 2016. [DOI: 10.1016/j.polymer.2016.10.053] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
47
|
Agrawal V, Peralta P, Li Y, Oswald J. A pressure-transferable coarse-grained potential for modeling the shock Hugoniot of polyethylene. J Chem Phys 2016; 145:104903. [DOI: 10.1063/1.4962255] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Vipin Agrawal
- School for the Engineering of Matter, Transport, and Energy, Arizona State University, Tempe, Arizona 85287, USA
| | - Pedro Peralta
- School for the Engineering of Matter, Transport, and Energy, Arizona State University, Tempe, Arizona 85287, USA
| | - Yiyang Li
- School for the Engineering of Matter, Transport, and Energy, Arizona State University, Tempe, Arizona 85287, USA
| | - Jay Oswald
- School for the Engineering of Matter, Transport, and Energy, Arizona State University, Tempe, Arizona 85287, USA
| |
Collapse
|
48
|
Uddin MS, Ju J. Multiscale modeling of a natural rubber: Bridging a coarse-grained molecular model to the rubber network theory. POLYMER 2016. [DOI: 10.1016/j.polymer.2016.08.037] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
49
|
Root SE, Savagatrup S, Pais CJ, Arya G, Lipomi DJ. Predicting the Mechanical Properties of Organic Semiconductors Using Coarse-Grained Molecular Dynamics Simulations. Macromolecules 2016. [DOI: 10.1021/acs.macromol.6b00204] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Samuel E. Root
- Department of NanoEngineering, University of California San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Suchol Savagatrup
- Department of NanoEngineering, University of California San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Christopher J. Pais
- Department of NanoEngineering, University of California San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Gaurav Arya
- Department of NanoEngineering, University of California San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| | - Darren J. Lipomi
- Department of NanoEngineering, University of California San Diego, 9500 Gilman Drive, Mail Code 0448, La Jolla, California 92093-0448, United States
| |
Collapse
|
50
|
Agrawal V, Holzworth K, Nantasetphong W, Amirkhizi AV, Oswald J, Nemat‐Nasser S. Prediction of viscoelastic properties with coarse‐grained molecular dynamics and experimental validation for a benchmark polyurea system. ACTA ACUST UNITED AC 2016. [DOI: 10.1002/polb.23976] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Vipin Agrawal
- School for the Engineering of Matter, Transport and EnergyArizona State UniversityTempe Arizona 85287
| | - Kristin Holzworth
- Department of Mechanical and Aerospace EngineeringCenter of Excellence for Advanced Materials, University of CaliforniaSan DiegoLa Jolla California92093‐0416
| | - Wiroj Nantasetphong
- Department of Mechanical and Aerospace EngineeringCenter of Excellence for Advanced Materials, University of CaliforniaSan DiegoLa Jolla California92093‐0416
| | | | - Jay Oswald
- School for the Engineering of Matter, Transport and EnergyArizona State UniversityTempe Arizona 85287
| | - Sia Nemat‐Nasser
- Department of Mechanical and Aerospace EngineeringCenter of Excellence for Advanced Materials, University of CaliforniaSan DiegoLa Jolla California92093‐0416
| |
Collapse
|