Dynamic aerofracture of dense granular packings

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 A deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNA(Phe), pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses.

We have applied all atoms discrete molecular dynamics (DMD) based on a quasi-physical potential to study the flexibility of an extended set of proteins for which atomistic MD simulations are available. The method uses pure physical potentials supplemented by information on secondary structure and despite its simplicity is able to reproduce with good accuracy the dynamics of proteins in solution. The method presents a clear improvement with respect to coarse-grained methods based on structure potentials and opens the possibility to explore dynamics of proteins out from the equilibrium and to trace conformational changes induced by interaction of proteins with both small and macromolecular ligands.

α-Synuclein (α-syn) is the major component of the intraneuronal inclusions called Lewy bodies, which are the pathological hallmark of Parkinson's disease. α-Syn is capable of self-assembly into many different species, such as soluble oligomers and fibrils. Even though attempts to resolve the structures of the protein have been made, detailed understanding about the structures and their relationship with the different aggregation steps is lacking, which is of interest to provide insights into the pathogenic mechanism of Parkinson's disease. Here we report the structural flexibility of α-syn monomers and dimers in an aqueous solution environment as probed by single-molecule time-lapse high-speed AFM. In addition, we present the molecular basis for the structural transitions using discrete molecular dynamics (DMD) simulations. α-Syn monomers assume a globular conformation, which is capable of forming tail-like protrusions over dozens of seconds. Importantly, a globular monomer can adopt fully extended conformations. Dimers, on the other hand, are less dynamic and show a dumbbell conformation that experiences morphological changes over time. DMD simulations revealed that the α-syn monomer consists of several tightly packed small helices. The tail-like protrusions are also helical with a small β-sheet, acting as a "hinge". Monomers within dimers have a large interfacial interaction area and are stabilized by interactions in the non-amyloid central (NAC) regions. Furthermore, the dimer NAC-region of each α-syn monomer forms a β-rich segment. Moreover, NAC-regions are located in the hydrophobic core of the dimer.

Visualization of Gravity Drainage in a Fractured System During Gas Injection Using Glass Micromodel

Computational insights into the aggregation mechanism of human calcitonin.

The abnormal aggregation of human calcitonin (hCT) hormone peptides impairs their physiological function, leading to harmful immune responses and cytotoxicity, which limits their clinical utility. Interestingly, a representative hCT analog incorporating Y12L and N17H substitutions (DM-hCT) has shown reduced aggregation tendencies while maintaining bioactivity. But the molecular mechanism of Y12L and N17H substitutions on the conformational dynamics of hCT remains unclear. Here, we systematically investigated the folding and self-assembly dynamics of hCT and DM-hCT using atomistic discrete molecular dynamics (DMD) simulations. Our findings revealed that hCT monomers predominantly adopted unstructured conformations with dynamic helices. Oligomerization of hCT resulted in the formation of β-sheet-rich aggregates and β-barrel intermediates. The Y12L and N17H substitutions enhanced helical conformations and suppressed β-sheet formation in both monomers and oligomers. These substitutions stabilized the dynamic helices and disrupted aromatic interactions responsible for β-sheet formation at residue 12. Notably, DM-hCT assemblies still exhibited β-sheets in phenylalanine-rich and C-terminal hydrophobic regions, suggesting that future optimizations should focus on these areas. Our simulations provide insights into the molecular mechanisms underlying hCT aggregation and the amyloid-resistant effects of Y12L and N17H substitutions. These findings have valuable implications for the development of clinical hCT analogs.

The mechanical behavior of nanoparticle assemblies depends strongly on particle size, yet the underlying mechanisms remain insufficiently understood. In present study, we employ a scheme combining discrete element method (DEM) and molecular dynamics (MD) simulations to examine size-dependent strength and deformation in disordered copper nanoparticle assemblies. Granular packings generated by DEM were transformed into atomic models and subjected to uniaxial compression in MD simulations. Assemblies composed of nanoparticles with radius smaller than ~2.5 nm fully densify during relaxation, forming nanopolycrystalline solids, whereas larger particles preserve porous architectures. This structural divergence governs subsequent deformation. Small-particle assemblies deform through grain boundary migration and grain growth, exhibiting an inverse Hall–Petch-type strength dependence. In contrast, large-particle assemblies deform primarily via interparticle contact evolution and densification, with strength conforming to a Gibson–Ashby-type prediction. A scaling law captures the strength variation across size range in this regime. These results establish the competition between surface energy-driven densification and contact-dominated deformation as the controlling factor in the mechanical response of nanoparticle assemblies, providing guidance for designing nanoparticle-based materials with tailored mechanical performance.

Multiscale Modeling of Nucleosome Dynamics

Experimental findings suggest that oligomeric forms of the amyloid beta protein (Abeta) play a critical role in Alzheimer's disease. Thus, elucidating their structure and the mechanisms of their formation is critical for developing therapeutic agents. We use discrete molecular dynamics simulations and a four-bead protein model to study oligomerization of two predominant alloforms, Abeta40 and Abeta42, at the atomic level. The four-bead model incorporates backbone hydrogen-bond interactions and amino acid-specific interactions mediated through hydrophobic and hydrophilic elements of the side chains. During the simulations we observe monomer folding and aggregation of monomers into oligomers of variable sizes. Abeta40 forms significantly more dimers than Abeta42, whereas pentamers are significantly more abundant in Abeta42 relative to Abeta40. Structure analysis reveals a turn centered at Gly-37-Gly-38 that is present in a folded Abeta42 monomer but not in a folded Abeta40 monomer and is associated with the first contacts that form during monomer folding. Our results suggest that this turn plays an important role in Abeta42 pentamer formation. Abeta pentamers have a globular structure comprising hydrophobic residues within the pentamer's core and hydrophilic N-terminal residues at the surface of the pentamer. The N termini of Abeta40 pentamers are more spatially restricted than Abeta42 pentamers. Abeta40 pentamers form a beta-strand structure involving Ala-2-Phe-4, which is absent in Abeta42 pentamers. These structural differences imply a different degree of hydrophobic core exposure between pentamers of the two alloforms, with the hydrophobic core of the Abeta42 pentamer being more exposed and thus more prone to form larger oligomers.

Lipid aggregation into fluid bilayers is an essential process for sustaining life. Simplified models of lipid structure, which allow for long time scales or large length scales not obtainable with all-atom simulations, have recently been developed and show promise for describing lipid dynamics in biological systems. Here, we describe two simplified models, a reduced-lipid model and a bola-lipid model for thermophilic bacterial membranes, developed for use with the rapid discrete molecular dynamics simulation method. In the reduced-lipid model, we represent the lipid chain by a series of three beads interacting through pairwise discrete potentials that model hydrophobic attractions between hydrocarbon tails in implicit solvent. Our phase diagram recapitulates those produced by continuous potential models with similar coarse-grained lipid representations. We also find that phase transition temperatures for our reduced-lipid model are dependent upon the flexibility of the lipid chain, giving an insight into archaebacterial membrane stability and prompting development of a bola-lipid model specific for archaebacteria lipids. With both the reduced-lipid and bola-lipid model, we find that the reduced flexibility inherent in archaebacteria lipids yields more stable bilayers as manifested by increased phase transition temperatures. The results of these studies provide a simulation methodology for lipid molecules in biological systems and show that discrete molecular dynamics is applicable to lipid aggregation and dynamics.

This paper presents a three-layer magnetic functional film structure composed of a substrate, a dielectric layer, and a two-dimensional structural layer, with its electromagnetic wave absorption characteristics analyzed through numerical simulation. The results demonstrate that: (1) When the substrate and the two-dimensional layer are made of the same magnetic material, increasing the dielectric layer thickness ( d ) causes the absorption peak frequency ( f max ) to first blueshift and then redshift, while the full width at half maximum (FWHM) initially decreases and then increases. This phenomenon originates from multiple reflection interference effects; (2) When the substrate is fixed as a perfect electric conductor (PEC), an increase in the permeability ( μ ) of the two-dimensional layer leads to a red-shift in the peak frequency ( f max ) and a broadening of the bandwidth—both consistent with skin effect theory; (3) When the two-dimensional structural layer is fixed as PEC, the variation pattern of substrate permeability shows similar characteristics, which aligns with the reflection phase properties of magnetic materials. The study illuminates the synergistic mechanism by which dielectric layer thickness ( d ) and magnetic parameters regulate wave absorption performance. Through optimization, nearly perfect absorption (≈ 100%) can be achieved, demonstrating significant potential for applications in microwave stealth and electromagnetic compatibility.≈

A conductive heat transfer model for particle flows over immersed surfaces

Powder discharge from a hopper-standpipe system modelled with CPFD

The N-methyl-D-aspartate (NMDA)-sensitive glutamate receptor (NMDAR) helps assemble downstream signaling pathways through protein interactions within the postsynaptic density (PSD), which are mediated by its intracellular C-terminal domain (CTD). The most abundant NMDAR subunits in the brain are GluN2A and GluN2B, which are associated with a developmental switch in NMDAR composition. Previously, we used single molecule fluorescence resonance energy transfer (smFRET) to show that the GluN2B CTD contained an intrinsically disordered region with slow, hop-like conformational dynamics. The CTD from GluN2B also undergoes liquid–liquid phase separation (LLPS) with synaptic proteins. Here, we extend these observations to the GluN2A CTD. Sequence analysis showed that both subunits contain a form of intrinsic disorder classified as weak polyampholytes. However, only GluN2B contained matched patterning of arginine and aromatic residues, which are linked to LLPS. To examine the conformational distribution, we used discrete molecular dynamics (DMD), which revealed that GluN2A favors extended disordered states containing secondary structures while GluN2B favors disordered globular states. In contrast to GluN2B, smFRET measurements found that GluN2A lacked slow conformational dynamics. Thus, simulation and experiments found differences in the form of disorder. To understand how this affects protein interactions, we compared the ability of these two NMDAR isoforms to undergo LLPS. We found that GluN2B readily formed condensates with PSD-95 and SynGAP, while GluN2A failed to support LLPS and instead showed a propensity for colloidal aggregation. That GluN2A fails to support this same condensate formation suggests a developmental switch in LLPS propensity.

From Disordered Amyloid β-Proteins to Soluble Oligomers and Protofibrils Using Discrete Molecular Dynamics