Molecular Dynamics Simulations Cyclotide Kalata B1 Interactions with Lipid Bilayers

Cyclotides are exceptionally stable plant peptides whose biological activity is widely attributed to interactions with lipid membranes, yet the molecular mechanisms underlying these interactions remain incompletely resolved. Here, we employ microsecond-scale (1 μs) all-atom molecular dynamics simulations to investigate the membrane association of the cyclotide kalata B1 with phospholipid bilayers of distinct headgroup composition, including POPC, POPE, and POPG. This extended timescale enables full bilayer equilibration and allows observation of slower peptide-induced membrane responses that are not accessible in shorter simulations. Across all systems, kalata B1 rapidly adsorbs to the membrane surface and remains predominantly surface-associated throughout the simulations, while the cyclic cystine knot motif remains structurally intact, confirming the exceptional robustness of the cyclotide fold during membrane engagement. Lipid-dependent differences arise primarily from variations in peptide orientation, conformational flexibility, and interfacial dynamics rather than deep bilayer insertion or pore formation. Zwitterionic POPC membranes favor compact, upright peptide configurations, whereas POPE and POPG bilayers promote enhanced lateral spreading and dynamic reorganization driven by hydrogen bonding and electrostatic interactions, respectively. Leaflet-resolved analyses of lipid contacts, membrane thickness, and area per lipid reveal localized, asymmetric perturbations confined to the peptide-exposed leaflet, with no evidence of sustained bilayer thinning or global destabilization. Together, these results support an interfacial, headgroup-dependent mechanism of cyclotide membrane activity and reconcile previous experimental observations. This work provides molecular-level insight into lipid selectivity and early-stage cyclotide-membrane interactions that may inform future design of cyclotide-based bioactive agents.