Optimized Multiple Circular Sequence Alignment for Cyclic Peptide Motif Discovery

Abstract

Head-to-tail (H2T) cyclized peptides are an increasingly important modality in drug discovery, combining high target affinity and selectivity with metabolic stability. Because their underlying chemistry is still that of a linear amino-acid chain, their linear sequence representation is the native input format of main-stream sequence generative models now driving de novo peptide design (Slough et al., 2018; Rettie et al., 2025a;b). Discovering the conserved motifs responsible for a family’s function across a library of such candidates requires a multiple sequence alignment (MSA). Because a cyclic peptide can be linearised at any residue, the alignment must additionally solve for the unknown rotation of each sequence, which is the multiple circular sequence alignment (MCSA) problem. However, leading MCSA heuristics (e.g. Ayad & Pissis, 2017) were tuned for the genomic regime (a few tens of long sequences) and become prohibitively slow on the cyclic peptide library regime (hundreds to thousands of shorter sequences). We close this gap by identifying quality-preserving optimisation opportunities, notably the library-scale preset tailored to short-sequence inputs (algorithmic details in Appendix A), and by adding an orthogonal multi-core and SIMD backend for further performance tuning, which gives near-linear thread scaling on the pairwise-comparison stage. We validate the pipeline on a library of 1,000 H2T cyclized peptides of length 18 targeting the oncoprotein Mouse double minute 2 human homolog (MDM2) produced by an internal peptide-design engine. In this practical setup, the optimised MCSA recovers the underlying positional motif of MDM2 binders at the same fidelity as the original MCSA implementation while running over 650× faster. Our optimised MCSA tool thus enables library-scale cyclic peptide sequence alignment and is publicly available at https://github.com/IVB-Generative-Biology/mars-turbo.