skip to main content


Title: Development of a PointNet for Detecting Morphologies of Self-Assembled Block Oligomers in Atomistic Simulations
ABSTRACT: Molecular simulations with atomistic or coarse- 6 grained force fields are a powerful approach for understanding and 7 predicting the self-assembly phase behavior of complex molecules. 8 Amphiphiles, block oligomers, and block polymers can form 9 mesophases with different ordered morphologies describing the 10 spatial distribution of the blocks, but entirely amorphous nature for 11 local packing and chain conformation. Screening block oligomer 12 chemistry and architecture through molecular simulations to find 13 promising candidates for functional materials is aided by effective 14 and straightforward morphology identification techniques. Captur- 15 ing 3-dimensional periodic structures, such as ordered network 16 morphologies, is hampered by the requirement that the number of 17 molecules in the simulated system and the shape of the periodic simulation box need to be commensurate with those of the resulting 18 network phase. Common strategies for structure identification include structure factors and order parameters, but these fail to 19 identify imperfect structures in simulations with incorrect system sizes. Building upon pioneering work by DeFever et al. [Chem. Sci. 20 2019, 10, 7503−7515] who implemented a PointNet (i.e., a neural network designed for computer vision applications using point 21 clouds) to detect local structure in simulations of single-bead particles and water molecules, we present a PointNet for detection of 22 nonlocal ordered morphologies of complex block oligomers. Our PointNet was trained using atomic coordinates from molecular 23 dynamics simulation trajectories and synthetic point clouds for ordered network morphologies that were absent from previous 24 simulations. In contrast to prior work on simple molecules, we observe that large point clouds with 1000 or more points are needed 25 for the more complex block oligomers. The trained PointNet model achieves an accuracy as high as 0.99 for globally ordered 26 morphologies formed by linear diblock, linear triblock, and 3-arm and 4-arm star-block oligomers, and it also allows for the discovery 27 of emerging ordered patterns from nonequilibrium systems.  more » « less
Award ID(s):
2011401
NSF-PAR ID:
10228125
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The journal of physical chemistry
ISSN:
1520-6106
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Identifying local structure in molecular simulations is of utmost importance. The most common existing approach to identify local structure is to calculate some geometrical quantity referred to as an order parameter. In simple cases order parameters are physically intuitive and trivial to develop ( e.g. , ion-pair distance), however in most cases, order parameter development becomes a much more difficult endeavor ( e.g. , crystal structure identification). Using ideas from computer vision, we adapt a specific type of neural network called a PointNet to identify local structural environments in molecular simulations. A primary challenge in applying machine learning techniques to simulation is selecting the appropriate input features. This challenge is system-specific and requires significant human input and intuition. In contrast, our approach is a generic framework that requires no system-specific feature engineering and operates on the raw output of the simulations, i.e. , atomic positions. We demonstrate the method on crystal structure identification in Lennard-Jones (four different phases), water (eight different phases), and mesophase (six different phases) systems. The method achieves as high as 99.5% accuracy in crystal structure identification. The method is applicable to heterogeneous nucleation and it can even predict the crystal phases of atoms near external interfaces. We demonstrate the versatility of our approach by using our method to identify surface hydrophobicity based solely upon positions and orientations of surrounding water molecules. Our results suggest the approach will be broadly applicable to many types of local structure in simulations. 
    more » « less
  2. The PHF6 (Val-Gln-Ile-Val-Tyr-Lys) motif, found in all isoforms of the microtubule-associated protein tau, forms an integral part of ordered cores of amyloid fibrils formed in tauopathies and is thought to play a fundamental role in tau aggregation. Because PHF6 as an isolated hexapeptide assembles into ordered fibrils on its own, it is investigated as a minimal model for insight into the initial stages of aggregation of larger tau fragments. Even for this small peptide, however, the large length and time scales associated with fibrillization pose challenges for simulation studies of its dynamic assembly, equilibrium configurational landscape, and phase behavior. Here, we develop an accurate, bottom-up coarse-grained model of PHF6 for large-scale simulations of its aggregation, which we use to uncover molecular interactions and thermodynamic driving forces governing its assembly. The model, not trained on any explicit information about fibrillar structure, predicts coexistence of formed fibrils with monomers in solution, and we calculate a putative equilibrium phase diagram in concentration-temperature space. We also characterize the configurational and free energetic landscape of PHF6 oligomers. Importantly, we demonstrate with a model of heparin that this widely studied cofactor enhances the aggregation propensity of PHF6 by ordering monomers during nucleation and remaining associated with growing fibrils, consistent with experimentally characterized heparin–tau interactions. Overall, this effort provides detailed molecular insight into PHF6 aggregation thermodynamics and pathways and, furthermore, demonstrates the potential of modern multiscale modeling techniques to produce predictive models of amyloidogenic peptides simultaneously capturing sequence-specific effects and emergent aggregate structures.

     
    more » « less
  3. Self-assembling dendrimers have facilitated the discovery of periodic and quasiperiodic arrays of supramolecular architectures and the diverse functions derived from them. Examples are liquid quasicrystals and their approximants plus helical columns and spheres, including some that disregard chirality. The same periodic and quasiperiodic arrays were subsequently found in block copolymers, surfactants, lipids, glycolipids, and other complex molecules. Here we report the discovery of lamellar and hexagonal periodic arrays on the surface of vesicles generated from sequence-defined bicomponent monodisperse oligomers containing lipid and glycolipid mimics. These vesicles, known as glycodendrimersomes, act as cell-membrane mimics with hierarchical morphologies resembling bicomponent rafts. These nanosegregated morphologies diminish sugar–sugar interactions enabling stronger binding to sugar-binding proteins than densely packed arrangements of sugars. Importantly, this provides a mechanism to encode the reactivity of sugars via their interaction with sugar-binding proteins. The observed sugar phase-separated hierarchical arrays with lamellar and hexagonal morphologies that encode biological recognition are among the most complex architectures yet discovered in soft matter. The enhanced reactivity of the sugar displays likely has applications in material science and nanomedicine, with potential to evolve into related technologies. 
    more » « less
  4. Context. The 1°.3 (G1.3) and 1°.6 (G1.6) cloud complexes in the central molecular zone (CMZ) of our Galaxy have been proposed to possibly reside at the intersection region of the X1 and X2 orbits for several reasons. This includes the detection of co-spatial low- and high-velocity clouds, high velocity dispersion, high fractional molecular abundances of shock-tracing molecules, and kinetic temperatures that are higher than for usual CMZ clouds. Aims. By investigating the morphology and deriving physical properties as well as chemical composition, we want to find the origin of the turbulent gas and, in particular, whether evidence of an interaction between clouds can be identified. Methods. We mapped both cloud complexes in molecular lines in the frequency range from 85 to 117 GHz with the IRAM 30 m telescope. The APEX 12m telescope was used to observe higher frequency transitions between 210 and 475 GHz from selected molecules that are emitted from higher energy levels. We performed non-local thermodynamic equilibrium (non-LTE) modelling of the emission of an ensemble of CH 3 CN lines to derive kinetic temperatures and H 2 volume densities. These were used as starting points for non-LTE modelling of other molecules, for which column densities and abundances were determined and compared with values found for other sources in the CMZ. Results. The kinematic structure of G1.3 reveals an ‘emission bridge’ at intermediate velocities (~150 km s −1 ) connecting low-velocity (~100 km s −1 ) and high-velocity (~180 km s −1 ) gas and an overall fluffy shell-like structure. These may represent observational evidence of cloud-cloud interactions. Low- and high-velocity gas components in G1.6 do not show this type of evidence of an interaction, suggesting that they are spatially separated. We selected three positions in each cloud complex for further analysis. Each position reveals several gas components at various peak velocities and of various line widths. We derived kinetic temperatures of 60–100 K and H 2 volume densities of 10 4 –10 5 cm −3 in both complexes. Molecular abundances relative to H 2 suggest a similar chemistry of the two clouds, which is moreover similar to that of other GC clouds and, especially, agrees well with that of G+0.693 and G−0.11. Conclusions. We conclude that G1.3 may indeed exhibit signs of cloud-cloud interactions. In particular, we propose an interaction of gas that is accreted from the near-side dust lane to the CMZ, with gas pre-existing at this location. Low- and high-velocity components in G1.6 are rather coincidentally observed along the same line of sight. They may be associated with either overshot decelerated gas from the far-side dust line or actual CMZ gas and high-velocity gas moving on a dust lane. These scenarios would be in agreement with numerical simulations. 
    more » « less
  5. Context. Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims. We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods. We present new observations of the dust polarized emission at 850 μ m toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with the James Clerk Maxwell Telescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity ( PI ), the polarization fraction ( PF ), and the plane-of-the-sky B-field angle ( χ B_POS ) toward the whole region, along the 10 pc-long ridge and along the sub-filaments connected to the ridge and the hubs. We derived the power spectra of the intensity and χ B POS along the ridge crest and compared them with the results obtained from simulated filaments. Results. The observations span ~3 orders of magnitude in Stokes I and PI and ~2 orders of magnitude in PF (from ~0.2 to ~ 20%). A large scatter in PI and PF is observed for a given value of I . Our analyses show a complex B-field structure when observed over the whole region (~ 10 pc); however, at smaller scales (~1 pc), χ B POS varies coherently along the crests of the filament network. The observed power spectrum of χ B POS can be well represented with a power law function with a slope of − 1.33 ± 0.23, which is ~20% shallower than that of I . We find that this result is compatible with the properties of simulated filaments and may indicate the physical processes at play in the formation and evolution of star-forming filaments. Along the sub-filaments, χ B POS rotates frombeing mostly perpendicular or randomly oriented with respect to the crests to mostly parallel as the sub-filaments merge with the ridge and hubs. This variation of the B-field structure along the sub-filaments may be tracing local velocity flows of infalling matter in the ridge and hubs. Our analysis also suggests a variation in the energy balance along the crests of these sub-filaments, from magnetically critical or supercritical at their far ends to magnetically subcritical near the ridge and hubs. We also detect an increase in PF toward the high-column density ( N H 2 ≳ 10 23  cm −2 ) star cluster-forming hubs. These latter large PF values may be explained by the increase in grain alignment efficiency due to stellar radiation from the newborn stars, combined with an ordered B-field structure. Conclusions. These observational results reveal for the first time the characteristics of the small-scale (down to ~ 0.1 pc) B-field structure of a 10 pc-long hub-filament system. Our analyses show variations in the polarization properties along the sub-filaments that may be tracing the evolution of their physical properties during their interaction with the ridge and hubs. We also detect an impact of feedback from young high-mass stars on the local B-field structure and the polarization properties, which could put constraints on possible models for dust grain alignment and provide important hints as to the interplay between the star formation activity and interstellar B-fields. 
    more » « less