Thermodynamic and Structural Signatures of Arginine Self-Assembly Across Concentration Regimes

Arginine plays a critical role in biomolecular interactions due to its guanidinium side chain, which enables multivalent electrostatic and hydrogen bonding contacts. In this study, atomistic molecular dynamics simulations were conducted across a broad concentration range (26–605 mM) to investigate the thermodynamic and structural features of arginine self-assembly in aqueous solution. Key observables—including hydrogen bond count, radius of gyration, contact number, and isobaric heat capacity—were analyzed to characterize emergent behavior. A three-regime aggregation pattern (dilute, cooperative, and saturated) was identified and quantitatively modeled using the Hill equation, revealing a non-linear transition in clustering behavior. Spatial analyses were supplemented with trajectory-based clustering and radial distribution functions. The heat capacity peak observed near 360 mM was interpreted as a thermodynamic signature of hydration rearrangement. Trajectory analyses utilized both GROMACS tools and the MDAnalysis library. While force field limitations and single-replica sampling are acknowledged, the results offer mechanistic insight into how arginine concentration modulates molecular organization—informing the understanding of biomolecular condensates, protein–nucleic acid complexes, and the design of functional supramolecular systems. The findings are in strong agreement with experimental observations from small-angle X-ray scattering and differential scanning calorimetry. Overall, this work establishes a cohesive framework for understanding amino acid condensation and reveals arginine’s concentration-dependent behavior as a model for weak, reversible molecular association.