Frequency-Dependent Structural Responses of Amyloid-β Fibrils to Terahertz Pulses: A Nonequilibrium Molecular Dynamics Study

The abnormal aggregation of amyloid proteins into fibrils is closely associated with several neurodegenerative diseases, highlighting the importance of developing effective strategies to modulate fibril structural stability. Here, we employed nonequilibrium molecular dynamics(NEMD) simulations to examine the effects of intense terahertz-to-mid-infrared pulses on three amyloid fibrils (2LMP, 2MXU, and 5OQV), focusing on the combined effects of pulse frequency and electric-field strength. Frequency scans from 0.5 to 60 THz in 0.5 THz steps revealed frequency-dependent losses of β-sheet content, with several common response patterns in these fibrils. Comparisons at 2, 23, 33, and 50 THz showed that pulse-induced changes in secondary structure differed from those in high-temperature controls. Hydrogen-bond analysis showed that total hydrogen-bond counts could partially recover even as interchain hydrogen bonds continued to decrease, with distinct responses among representative chain pairs. We therefore propose a two-threshold electric-field framework for fibril disruption: Eth1 describes the onset of structural perturbation, whereas Eth2 indicates persistent structural loss after the pulse. This framework offers a testable approach for characterizing amyloid structural responses to external fields, although the underlying energy-coupling mechanisms, quantitative thresholds, and long-term recovery require further validation. These findings may provide a theoretical basis for exploring external-field-based strategies to modulate protein aggregation.