GyroNovo: 基于误差引导的片段插补与质量感知注意力的从头肽测序
摘要
GyroNovo通过基于解码器误差自适应地插补缺失片段,并利用旋转嵌入整合质量差异,改进了从头肽测序,在基准测试中达到了最先进的性能。
arXiv:2609.30542v1 Announce Type: new
Abstract: De novo peptide sequencing from tandem mass spectra is essential for identifying peptides without relying on reference databases. Despite advances in deep learning, accurate sequencing remains challenging because experimental spectra are often sparse, noisy, and incomplete, leaving informative b- and y-ion fragments unobserved. Existing methods attempt to recover this missing evidence via latent-space imputation before autoregressive decoding. However, they typically treat imputation as a fixed reconstruction task, without considering which missing fragments are most relevant to decoder errors. Moreover, existing peak representations do not explicitly model mass differences between peaks, despite their fundamental importance.
We introduce GyroNovo, a framework with two main contributions. First, we use decoder errors observed during training to adapt the imputation objective, prioritizing fragments associated with frequent decoding errors. We further use the decoder error distribution to construct easy and hard augmented views of each spectrum, enabling the decoder to learn under varying degrees of spectral corruption and missing-fragment severity. Second, we introduce a mass-aware inductive bias into self-attention by using rotary embeddings to encode pairwise mass differences between spectral peaks. Together, these components align missing-fragment recovery with decoder behavior while explicitly incorporating the mass relationships that underlie peptide fragmentation. At inference time, GyroNovo retains a standard encoder-imputer-decoder architecture and requires neither additional inputs nor auxiliary search procedures. Experiments on NovoBench show gains of about 9 percentage points in peptide-level precision and 7 percentage points in amino-acid-level precision over the state-of-the-art baseline. Code: https://github.com/UBC-NLP/gyronovo.
查看缓存全文
缓存时间: 2026/09/29 09:39
# GyroNovo: Error-Guided Fragment Imputation with Mass-Aware Attention for De Novo Peptide Sequencing
Source: [https://arxiv.org/html/2609.30542](https://arxiv.org/html/2609.30542)
###### Abstract
De novopeptide sequencing from tandem mass spectra is essential for identifying peptides without relying on reference databases\. Despite recent advances in deep learning, accurate sequencing remains challenging because experimental spectra are often sparse, noisy, and incomplete, causing informative b\- and y\-ion fragments to be unobserved\. Existing methods have attempted to recover this missing evidence through latent\-space imputation before autoregressive decoding\. However, they typically formulate imputation as a fixed reconstruction task, without considering which missing fragments are most relevant to the decoder’s current sequencing errors\. Moreover, existing peak representations do not explicitly model mass differences between peaks, despite their fundamental importance inde novopeptide sequencing\.
We introduce GyroNovo, a framework with two main contributions\. First, we use decoder errors observed during training to adapt the imputation objective, assigning greater emphasis to fragments associated with frequent decoding errors\. We further use the decoder error distribution to construct easy and hard augmented views of each peptide spectrum, enabling the decoder to learn under varying degrees of spectral corruption and missing\-fragment severity\. Second, we introduce a mass\-aware inductive bias into self\-attention by using rotary embeddings to encode pairwise mass differences between spectral peaks\. Together, these components align missing\-fragment recovery with the decoder’s sequencing behavior while explicitly incorporating the mass relationships that underlie peptide fragmentation\. At inference time, GyroNovo retains a standard encoder–imputer–decoder architecture and requires neither additional inputs nor auxiliary search procedures\. Experiments on NovoBench demonstrate improvements of approximately 9 percentage points in peptide\-level precision and 7 percentage points in amino\-acid\-level precision over the previous state\-of\-the\-art baseline\. Our code is publicly available athttps://github\.com/UBC\-NLP/gyronovo\.
1The University of British Columbia
2Canada Research Chair in NLP and ML
\{abdellah\.elmekki, muhammad\.mageed\}@ubc\.ca, laks@cs\.ubc\.ca
## Introduction
De novopeptide sequencing is a central task in computational proteomics: given a tandem mass spectrum, the goal is to infer the underlying amino\-acid sequence without relying on a protein database\. Accuratede novosequencing enables discovery\-oriented workflows, including antibody and antigen characterization\([Tran et al\. 2019](https://arxiv.org/html/2609.30542#bib.bib20)\), supports the analysis of organisms or samples with incomplete or unavailable reference databases\([Yilmaz et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib1)\), and helps identify novel or modified peptides that are difficult to recover using conventional database search\([HE et al\. 2013](https://arxiv.org/html/2609.30542#bib.bib21);[Chick et al\. 2015](https://arxiv.org/html/2609.30542#bib.bib22)\)\. Recent neural approaches, including influential Transformer\-based sequence decoders\([Yilmaz et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib1)\)and other systems systematically evaluated in NovoBench\([Zhou et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib4)\), have substantially advanced the field by learning rich representations of mass spectra and decoding peptide sequences directly from observed peaks\. Despite these advances, practical spectra remain difficult because fragmentation is incomplete, ion intensities are irregular, noise peaks are common, and post\-translational modifications \(PTMs\) can obscure otherwise useful mass evidence\([McDonnell et al\. 2022](https://arxiv.org/html/2609.30542#bib.bib13)\)\.
A major source of sequencing error is missing fragmentation evidence\([Mao et al\. 2023](https://arxiv.org/html/2609.30542#bib.bib16);[Du et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib2)\)\. Autoregressive decoders predict each amino acid from a partial sequence and a spectrum representation, even when evidence for the relevant cleavage is weak or absent\. LIPNovo addresses this limitation through latent imputation: learnable peak queries infer representations of theoretical fragments, and the peptide decoder jointly uses observed and imputed evidence\([Du et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib2)\)\. This design directly targets missing fragments and provides a strong foundation for robustde novosequencing\. However, LIPNovo does not explicitly prioritize the missing fragments most useful for correcting current decoding errors\. Its imputation supervision is determined mainly by theoretical peak targets and bipartite matching, which encourages plausible fragment reconstruction but treats easy and difficult residue decisions similarly\. As a result, the model may emphasize valid peaks with limited downstream value while underweighting fragments whose absence causes sequencing errors\.
Current methods also provide minor explicit structure for relationships among observed peaks\. Tandem mass spectra are organized by continuous mass differences rather than discrete sequence positions: residue transitions, fragmentation ladders, complementary ions, and neutral losses all induce characteristic mass separations\. Standard transformer encoders represent each peak using absolute m/z and intensity, leaving self\-attention to infer these relative relationships implicitly\. Introducing a relative\-mass inductive bias could therefore better align the encoder with the geometry of tandem mass spectra\.
We propose GyroNovo, a method that addresses both limitations through decoder\-guided latent imputation and continuous mass\-aware rotary attention\. During training, token\-level decoder errors are mapped to nearby theoretical fragment cleavages, assigning greater imputation weight to evidence associated with difficult residue decisions\. The same feedback guides peptide\-conditioned spectral augmentation, emphasizing informative spectra and cleavage\-local regions while constructing easier and harder views that preserve observed peaks and selectively modify theoretical fragments\. We also introduce continuous mass\-based rotary attention in the spectrum encoder\. Instead of using peak\-list indices, rotation angles are derived from each peak’s continuous m/z coordinate at multiple wavelengths\. Attention therefore depends on relative mass displacement, enabling multi\-resolution modeling of mass relationships\.
In summary, this work makes the following contributions: \(i\) a decoder\-conditioned latent imputation objective that projects token\-level sequencing errors onto cleavage\-local theoretical fragment targets; \(ii\) an observed\-anchor, peptide\-conditioned strategy that produces realistic easy and hard spectral views with both imputation and decoding supervision; \(iii\) a continuous mass\-based rotary self\-attention mechanism that introduces relative mass differences directly into the spectrum encoder; and \(iv\) an empirical NovoBench evaluation demonstrating state\-of\-the\-artde novopeptide sequencing, with approximately 9\- and 7\-percentage\-point gains in peptide\- and amino\-acid\-level precision, respectively\.
## Related Work
### De NovoPeptide Sequencing
Classicalde novosequencing methods formulate peptide identification as path finding or dynamic programming over spectrum graphs, using handcrafted or probabilistic scoring functions\([Ma et al\. 2003](https://arxiv.org/html/2609.30542#bib.bib11);[Frank and Pevzner 2005](https://arxiv.org/html/2609.30542#bib.bib14);[Ma 2015](https://arxiv.org/html/2609.30542#bib.bib15)\)\. DeepNovo introduced neural next\-residue prediction, followed by architectures that model spectra as unordered peak sets or graphs, including PointNovo, GraphNovo, and Spectralis\([Tran et al\. 2017](https://arxiv.org/html/2609.30542#bib.bib5);[Xu et al\. 2023](https://arxiv.org/html/2609.30542#bib.bib8);[Mao et al\. 2023](https://arxiv.org/html/2609.30542#bib.bib16);[Klaproth\-Andrade et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib17)\)\. Nevertheless, sequencing accuracy remains highly sensitive to noise and missing fragment ions\([Zhou et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib4)\)\.
Deep learning models for mass spectrometry have employed convolutional, recurrent, point\-set, and graph architectures to learn representations directly from irregular MS/MS peaks\. Most methods encode the observed spectrum and immediately decode the peptide\. LIPNovo instead imputes latent representations of theoretical fragment peaks before prediction, directly targeting incomplete fragmentation\([Du et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib2)\)\.
### Transformer\-basedDe NovoPeptide Sequencing
Transformers provide a natural encoder\-decoder formulation for translating spectra into peptide sequences\([Vaswani et al\. 2017](https://arxiv.org/html/2609.30542#bib.bib18)\)\. Casanovo introduced direct autoregressive decoding from peak sequences, while subsequent methods such as ContraNovo, InstaNovo,π\\pi\-PrimeNovo, AdaNovo, andπ\\pi\-HelixNovo improved scalability, robustness, or decoding\([Yilmaz et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib1);[Jin et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib23);[Eloff et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib6);[Zhang et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib19);[Xia et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib9);[Yang et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib7)\)\. Conventional peak embeddings primarily encode absolute mass\-to\-charge values, requiring chemically meaningful mass differences to be learned implicitly\. Recent methods incorporate pairwise mass differences through attention biases\([Lapin et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib12)\)\.
## Proposed Method
Figure 1:Overview of GyroNovo ’s error\-guided imputation and difficulty\-aware multi\-view training framework\.### Problem Formulation
We consider the problem ofde novopeptide sequencing, whose objective is to recover the amino acid sequence of a peptide from an observed tandem mass spectrum and its precursor information\. An MS/MS spectrum is represented as a variable\-size set of peaksx=\{\(mi,Ii\)\}i=1N,x=\\\{\(m\_\{i\},I\_\{i\}\)\\\}\_\{i=1\}^\{N\},wheremim\_\{i\}andIiI\_\{i\}denote the mass\-to\-charge ratio and intensity of theii\-th peak, respectively\. Each spectrum is accompanied by precursor metadatap=\(M,z\)p=\(M,z\), whereMMis the precursor mass andzzis the precursor charge\. The target peptide is represented asy=\(y1,…,yL\),y=\(y\_\{1\},\\ldots,y\_\{L\}\),whereyly\_\{l\}is an amino acid from a vocabulary containing both unmodified and modified residues, together with a stop token\. The sequencing model factorizes the conditional probability of the peptide autoregressively:
P\(y∣x,p\)=∏l=1LP\(yl∣y<l,x,p\)\.P\(y\\mid x,p\)=\\prod\_\{l=1\}^\{L\}P\(y\_\{l\}\\mid y\_\{<l\},x,p\)\.Unlike database\-search methods,de novosequencing performs this prediction without access to a protein sequence database\.
A central challenge is that experimental spectra provide only a partial and noisy observation of the underlying fragmentation process\. Informativebb\- andyy\-ions may be absent, have low intensity, or overlap with unrelated noise peaks\. Consequently, the decoder must often infer amino acids from incomplete fragment\-ion ladders\.
A common neural formulation, exemplified by Casanovo\([Yilmaz et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib1)\), uses a Transformer encoder to represent the observed spectrum and an autoregressive Transformer decoder to generate the peptide sequence\. LIPNovo\([Du et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib2)\)extends this formulation with a latent imputation module that predicts representations of missing fragment ions before decoding\. In this work, we retain the latent imputation backbone of LIPNovo and introduce GyroNovo which has two key innovations: \(i\) decoder\-conditioned supervision that concentrates imputation and augmentation on fragments associated with difficult amino acids, and \(ii\) continuous mass rotary attention that makes spectrum self\-attention directly sensitive to relative mass differences\.
### Latent Imputation Backbone
We follow LIPNovo\([Du et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib2)\)for the latent imputation backbone\. Each observed peak is first mapped to a peak embedding by combining a continuous encoding of its m/z value with a learned projection of its intensity\. A spectrum encoderEθE\_\{\\theta\}then produces contextualized peak representations:
H=Eθ\(x,p\)=\[hcls,h1,…,hN\],H=E\_\{\\theta\}\(x,p\)=\[h\_\{\\mathrm\{cls\}\},h\_\{1\},\\ldots,h\_\{N\}\],wherehclsh\_\{\\mathrm\{cls\}\}is a learned global spectrum representation andhih\_\{i\}is the contextual representation of theii\-th peak\.
An imputation moduleGθG\_\{\\theta\}usesQQlearned queries to cross\-attend to the encoded spectrum:
\{\(u^j,r^j\)\}j=1Q=Gθ\(H\)\.\\\{\(\\hat\{u\}\_\{j\},\\hat\{r\}\_\{j\}\)\\\}\_\{j=1\}^\{Q\}=G\_\{\\theta\}\(H\)\.Here,u^j\\hat\{u\}\_\{j\}is a candidate latent fragment representation andr^j\\hat\{r\}\_\{j\}is its confidence logit\.
During training, the ground\-truth peptide is used to construct a theoretical spectrum containing itsbb\- andyy\-ions\. The theoretical spectrum is processed by the same spectrum encoder, yielding target fragment representations
U⋆=\{uk⋆\}k=1K\.U^\{\\star\}=\\\{u\_\{k\}^\{\\star\}\\\}\_\{k=1\}^\{K\}\.The predicted queries and theoretical fragments are aligned through bipartite matching\. For a training sampleii, the matching is obtained by minimizing
πi=argmin∑\(j,k\)∈ππ\[‖u^ij−uik⋆‖22\+1−σ\(r^ij\)\],\\pi\_\{i\}=\\arg\\min\_\{\\pi\}\\sum\_\{\(j,k\)\\in\\pi\}\\left\[\\\|\\hat\{u\}\_\{ij\}\-u\_\{ik\}^\{\\star\}\\\|\_\{2\}^\{2\}\+1\-\\sigma\(\\hat\{r\}\_\{ij\}\)\\right\],whereσ\(⋅\)\\sigma\(\\cdot\)is the sigmoid function\. Matched queries are trained to reconstruct theoretical fragment representations, while the confidence head is trained to distinguish matched queries from unmatched ones\.
The decoder is then fed with the confident imputed fragments and the observed peak representations:
H~=\[hcls;\{u^j:σ\(r^j\)\>τ\};h1,…,hN\],\\widetilde\{H\}=\\left\[h\_\{\\mathrm\{cls\}\};\\\{\\hat\{u\}\_\{j\}:\\sigma\(\\hat\{r\}\_\{j\}\)\>\\tau\\\};h\_\{1\},\\ldots,h\_\{N\}\\right\],whereτ\\tauis a confidence threshold\. It then generates the peptide sequence autoregressively\. The theoretical spectrum is required only to define training targets; at inference time, missing fragment representations are predicted directly from the observed spectrum\.
We now describe our contributions\. Figure[1](https://arxiv.org/html/2609.30542#Sx3.F1)provides an overview of all the contributions\. Below, we describe each component in detail\.
### Error\-Guided Fragment Imputation Supervision
The standard imputation objective assigns equal importance to all theoretical fragments\. However, accurately reconstructing every missing fragment is not equally useful for peptide sequencing\. In particular, fragments adjacent to amino acids that are already decoded correctly provide less useful training feedback than fragments associated with persistent decoder errors\.
We therefore use the current decoder errors to adapt the representation\-learning task in two complementary ways\. First, we increase the imputation weight of theoretical fragments associated with difficult amino acids\. Second, we construct an easy and a hard augmented spectrum view: the easy view exposes and strengthens difficult fragment evidence, whereas the hard view removes or weakens the same evidence to improve robustness\.
#### Imputation Loss Re\-Weighting\.
We first decode the peptide from the observed\-spectrum branch and compute the cross\-entropy loss at each valid amino acid position:
ℓil=−logPθ\(yil∣yi,<l,H~i,pi\)\.\\ell\_\{il\}=\-\\log P\_\{\\theta\}\\left\(y\_\{il\}\\mid y\_\{i,<l\},\\widetilde\{H\}\_\{i\},p\_\{i\}\\right\)\.Because raw cross\-entropy values are unbounded, we map them to a bounded token\-error score:
eil=1−exp\(−ℓil\)\.e\_\{il\}=1\-\\exp\\left\(\-\\ell\_\{il\}\\right\)\.
Each theoretical fragmentkkis annotated with its ion family and cleavage index\. These annotations determine the amino acid positions adjacent to the corresponding cleavage\. Let𝒩i\(k\)\\mathcal\{N\}\_\{i\}\(k\)denote these positions together with a small local neighborhood\.111In this work, we set the neighborhood size to 1\.The weight assigned to fragmentkkis
wik=clip\(1\+maxl∈𝒩i\(k\)eil,wmin,wmax\),w\_\{ik\}=\\operatorname\{clip\}\\left\(1\+\\max\_\{l\\in\\mathcal\{N\}\_\{i\}\(k\)\}e\_\{il\},\\,w\_\{\\min\},w\_\{\\max\}\\right\),wherewmin=1\.0w\_\{\\min\}=1\.0andwmax=2\.0w\_\{\\max\}=2\.0are hyperparameters\. Thus, a theoretical fragment receives a larger weight when it is close to an amino acid that the decoder currently predicts poorly\.
The decoder\-conditioned imputation loss is
ℒimp=∑\(j,k\)∈πiwik‖u^ij−uik⋆‖22∑\(j,k\)∈πiwik\+ℒconf,\\mathcal\{L\}\_\{\\mathrm\{imp\}\}=\\frac\{\\displaystyle\\sum\_\{\(j,k\)\\in\\pi\_\{i\}\}w\_\{ik\}\\\|\\hat\{u\}\_\{ij\}\-u\_\{ik\}^\{\\star\}\\\|\_\{2\}^\{2\}\}\{\\displaystyle\\sum\_\{\(j,k\)\\in\\pi\_\{i\}\}w\_\{ik\}\}\+\\mathcal\{L\}\_\{\\mathrm\{conf\}\},whereπi\\pi\_\{i\}is the bipartite matching for sampleii\. The confidence lossℒconf\\mathcal\{L\}\_\{\\mathrm\{conf\}\}is a binary cross\-entropy objective in which matched queries are positive and unmatched queries are negative\. The weightswikw\_\{ik\}are also applied to matched positive queries, encouraging the confidence head to retain fragments associated with difficult residues\.
#### Observed Mass Spectrum Augmentations\.
Additionally, we use decoder errors to construct two peptide\-conditioned augmented views of each observed spectrum\. The conditioning signal is obtained from a teacher\-forced forward pass through the observed\-spectrum branch\. For sampleii, letℓil\\ell\_\{il\}denote the cross\-entropy loss for ground\-truth amino acidll, and letLiL\_\{i\}denote the peptide length, excluding padding and the end\-of\-sequence token\. We define the sample\-level decoding difficulty as
ℓ¯i=1Li∑l=1Liℓil,e¯i=1−exp\(−ℓ¯i/τ\),\\bar\{\\ell\}\_\{i\}=\\frac\{1\}\{L\_\{i\}\}\\sum\_\{l=1\}^\{L\_\{i\}\}\\ell\_\{il\},\\qquad\\bar\{e\}\_\{i\}=1\-\\exp\(\-\\bar\{\\ell\}\_\{i\}/\\tau\),and compute the augmentation scale
si=clip\(1\+αe¯i,smin,smax\),s\_\{i\}=\\operatorname\{clip\}\\left\(1\+\\alpha\\bar\{e\}\_\{i\},\\,s\_\{\\min\},\\,s\_\{\\max\}\\right\),whereτ=1\\tau=1,α=1\\alpha=1,smin=1s\_\{\\min\}=1, andsmax=1\.5s\_\{\\max\}=1\.5\. This scale modulates a linear augmentation curriculum, with a lower base corruption strength for the easy view than for the hard view\. Consequently, samples with larger decoding errors receive stronger background perturbations, subject to operator\-specific upper bounds\.
To localize the augmentation, we additionally convert amino\-acid losses into residue\-level errors,eil=1−exp\(−ℓil/τ\)e\_\{il\}=1\-\\exp\(\-\\ell\_\{il\}/\\tau\), and associate each theoretical fragment with the residues surrounding its cleavage site\. The fragment weightwikw\_\{ik\}is determined by the largest residue\-level error in this local neighborhood\. These weights are normalized within each spectrum to obtain focus scoresfik∈\[0,1\]f\_\{ik\}\\in\[0,1\], which control where the two views preserve or remove evidence\. All conditioning quantities are detached from the computation graph\.
##### View construction\.
Both views are constructed from the observed spectrum\. We match each theoretical singly chargedbb\- oryy\-ion to its nearest observed peak when their absolutem/zm/zdifference is at most0\.50\.5Da\. This identifies matched observed peaks, unannotated observed peaks, and missing theoretical fragments\. The perturbation design is informed by error analysis on a nine\-species validation set separate from the test set, while the augmentation hyperparameters are selected by random search\.
The easy view increases the retention probability and intensity of observed peaks associated with high decoder error, and increases the insertion probability of corresponding missing or dropped theoretical fragments\. The hard view decreases these probabilities and attenuates the same observed peaks\. Both views also incorporate stochastic peak removal, ion\-family\-specific dropout, contiguous fragmentation\-ladder gaps, reduced visibility of modification\-containing fragments, and insertion of structured noise peaks\. Newly inserted theoretical peaks receive intensity andm/zm/zjitter\. Thus, the easy view supplies additional evidence near difficult residues, whereas the hard view encourages reconstruction and decoding when that evidence is incomplete\.
Both augmented views pass through the shared encoder, imputation module, and peptide decoder\. Each receives imputation supervision from the theoretical fragment representations and sequence supervision from the ground\-truth peptide; the easy and hard views contribute equally to the augmented\-view losses\. These augmentations are used only during training\. Full construction details and hyperparameters are provided in Appendix[A](https://arxiv.org/html/2609.30542#A1)\.
### Continuous Mass Rotary Attention
The baseline spectrum encoder represents each peak independently using its absolutem/zm/zand intensity\. Thus, although the query and key vectors contain absolute mass information, vanilla self\-attention does not explicitly expose the displacement between two peaks\. Such displacements are informative for peptide fragmentation: Fragments derived from the same precursor exhibit residue\-dependentm/zm/zincrements, and consecutive ions form approximately regular fragmentation ladders\.
To introduce this inductive bias, we retain the baseline peak embeddings but replace vanilla spectrum self\-attention with a continuous\-m/zm/zadaptation of Rotary Position Embedding \(RoPE\)\([Su et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib3)\)\. Unlike conventional RoPE, whose rotation angles are determined by discrete token indices, our rotation angles are determined directly by the continuousm/zm/zcoordinate of each peak\.
Letdhd\_\{h\}be the dimension of one attention head and assume it is even\. We partition each query and key vector intodh/2d\_\{h\}/2two\-dimensional pairs\. For pairrr, we use the fixed wavelength
λr=λmin\(λmaxλmin\)rdh/2−1,r=0,…,dh2−1,\\lambda\_\{r\}=\\lambda\_\{\\min\}\\left\(\\frac\{\\lambda\_\{\\max\}\}\{\\lambda\_\{\\min\}\}\\right\)^\{\\frac\{r\}\{d\_\{h\}/2\-1\}\},\\qquad r=0,\\ldots,\\frac\{d\_\{h\}\}\{2\}\-1,so that the wavelengths are logarithmically spaced betweenλmin\\lambda\_\{\\min\}= 1 andλmax\\lambda\_\{\\max\}= 10,000\. A peak atm/zm/zvaluemminduces the phase
ϕr\(m\)=2πmλr\\phi\_\{r\}\(m\)=\\frac\{2\\pi m\}\{\\lambda\_\{r\}\}and the corresponding rotation
Rr\(m\)=\[cosϕr\(m\)−sinϕr\(m\)sinϕr\(m\)cosϕr\(m\)\]\.R\_\{r\}\(m\)=\\begin\{bmatrix\}\\cos\\phi\_\{r\}\(m\)&\-\\sin\\phi\_\{r\}\(m\)\\\\ \\sin\\phi\_\{r\}\(m\)&\\cos\\phi\_\{r\}\(m\)\\end\{bmatrix\}\.LetR\(m\)=diag\(R0\(m\),…,Rdh/2−1\(m\)\)R\(m\)=\\operatorname\{diag\}\(R\_\{0\}\(m\),\\ldots,R\_\{d\_\{h\}/2\-1\}\(m\)\)\. We rotate only the query and key vectors, leaving the value vectors unchanged:
q~i=R\(mi\)qi,k~j=R\(mj\)kj\.\\widetilde\{q\}\_\{i\}=R\(m\_\{i\}\)q\_\{i\},\\qquad\\widetilde\{k\}\_\{j\}=R\(m\_\{j\}\)k\_\{j\}\.The resulting attention logit is
aij\\displaystyle a\_\{ij\}=q~i⊤k~jdh\\displaystyle=\\frac\{\\widetilde\{q\}\_\{i\}^\{\\top\}\\widetilde\{k\}\_\{j\}\}\{\\sqrt\{d\_\{h\}\}\}=qi⊤R\(mi\)⊤R\(mj\)kjdh\\displaystyle=\\frac\{q\_\{i\}^\{\\top\}R\(m\_\{i\}\)^\{\\top\}R\(m\_\{j\}\)k\_\{j\}\}\{\\sqrt\{d\_\{h\}\}\}=qi⊤R\(mj−mi\)kjdh\.\\displaystyle=\\frac\{q\_\{i\}^\{\\top\}R\(m\_\{j\}\-m\_\{i\}\)k\_\{j\}\}\{\\sqrt\{d\_\{h\}\}\}\.The final equality follows from the composition property of planar rotations\. Consequently, the rotary component of the attention score depends on the relative displacementΔmij=mj−mi\\Delta m\_\{ij\}=m\_\{j\}\-m\_\{i\}, rather than on the two coordinates only independently\. The same learned query\-key interaction can therefore recognize a characteristic fragmentation gap regardless of where it occurs in the spectrum\.
### Training Objective and Inference
The complete objective combines peptide decoding from the observed spectrum, peptide decoding from the theoretical spectrum, imputation on the observed spectrum, and imputation and decoding on the augmented views:
ℒ\\displaystyle\\mathcal\{L\}=ℒdec\(x\)\+ℒdec\(x⋆\)\+ℒimp\(x\)\\displaystyle=\\mathcal\{L\}\_\{\\mathrm\{dec\}\}\(x\)\+\\mathcal\{L\}\_\{\\mathrm\{dec\}\}\(x^\{\\star\}\)\+\\mathcal\{L\}\_\{\\mathrm\{imp\}\}\(x\)\(1\)\+ℒimp\(x~\)\+ℒdec\(x~\),\\displaystyle\+\\mathcal\{L\}\_\{\\mathrm\{imp\}\}\(\\widetilde\{x\}\)\+\\mathcal\{L\}\_\{\\mathrm\{dec\}\}\(\\widetilde\{x\}\),wherex⋆x^\{\\star\}is the theoretical spectrum andx~\\widetilde\{x\}denotes the set of training\-time easy and hard augmented views\. The augmented\-view losses are computed as weighted averages over the two views\.
At inference time, the model requires only the observed spectrum and precursor metadata\. The theoretical spectrum, decoder\-conditioned weights, and augmented views are discarded\. The observed spectrum is encoded using continuous mass \(m/z\) rotary attention, the learned imputation queries predict candidate missing\-fragment representations, and low\-confidence candidates are removed using the thresholdτ\\tau\(we set it to 0\.8 following[Du et al\. \(2025\)](https://arxiv.org/html/2609.30542#bib.bib2)\)\. Finally, the peptide decoder attends jointly to the observed and imputed spectrum representations and generates the peptide sequence using beam search\.
## Experiments
### Datasets
Consistent with the benchmark established by[Zhou et al\. \(2024\)](https://arxiv.org/html/2609.30542#bib.bib4), we evaluate our method on three datasets: Nine\-species and Seven\-species\([Tran et al\. 2017](https://arxiv.org/html/2609.30542#bib.bib5)\), and HC\-PT\([Eloff et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib6)\)\. The Nine\-species dataset, which has been widely adopted in prior research, includes mass spectra collected from nine species\. Spectra from yeast are reserved for testing, while data from the remaining eight species are used for model training and validation\. Similarly, the Seven\-species dataset contains mass spectra from seven species, with yeast serving as the held\-out test species and the other six used for training and validation\. The HC\-PT dataset comprises spectra from human\-derived peptides, including synthetic tryptic peptides covering canonical human proteins and their isoforms, peptides produced using alternative proteases, and human leukocyte antigen–associated peptides\.
### Baselines
We compare our proposed GyroNovo against strong baselines, including DeepNovo\([Tran et al\. 2017](https://arxiv.org/html/2609.30542#bib.bib5)\), PointNovo\([Xu et al\. 2023](https://arxiv.org/html/2609.30542#bib.bib8)\), CasaNovo\([Yilmaz et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib1)\), AdaNovo\([Xia et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib9)\), LIPNovo\([Du et al\. 2025](https://arxiv.org/html/2609.30542#bib.bib2)\), and the recent concurrent work LIPNovo\+\([Du et al\. 2026](https://arxiv.org/html/2609.30542#bib.bib10)\)\. Unless otherwise noted, the reported results are taken from NovoBench\. Because LIPNovo is our most direct competitor, we report both its NovoBench results and our reproduced results to ensure a fair comparison\. We refer to the reproduced version as the “baseline” throughout the paper and in the result tables\.
### Evaluation Metrics
We assessde novopeptide sequencing at three levels of granularity: amino acid, PTM, and complete peptide\. An amino acid prediction is considered correct if its mass differs from that of the corresponding reference residue by less than0\.10\.1Da and its prefix or suffix mass differs from that of the reference sequence by less than0\.50\.5Da\([Tran et al\. 2017](https://arxiv.org/html/2609.30542#bib.bib5)\)\. We report precision and recall at both the amino acid and PTM levels\.
For peptide\-level evaluation, a prediction is regarded as correct only when the entire predicted sequence exactly matches the reference sequence\. We report the peptide\-level precision\. Finally, we report the area under the precision–recall curve \(AUC\) to characterize performance across confidence\-score thresholds\.
### Implementation Details
To maintain a controlled comparison with prior work, we adopt the hyperparameter configuration reported by[Du et al\. \(2025\)](https://arxiv.org/html/2609.30542#bib.bib2)\. Our model follows the approach of[Du et al\. \(2025\)](https://arxiv.org/html/2609.30542#bib.bib2), which builds on an implementation of Casanovo\([Yilmaz et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib1)\)by incorporating complementary spectra\([Yang et al\. 2024](https://arxiv.org/html/2609.30542#bib.bib7)\)as an additional input modality and introducing an imputation module\.
The encoder and decoder each contain nine Transformer layers, with a hidden dimension of 512, a feed\-forward dimension of 1024, and eight attention heads\. The imputation module consists of three layers and uses 100 peak queries\. Input spectra are truncated or padded to a maximum of 150 peaks, and the maximum peptide sequence length is set to 100 residues\. The decoder vocabulary contains tokens for the 20 standard amino acids, three post\-translational modifications, and an end\-of\-sequence token\.
We train once the model using distributed data parallelism on four NVIDIA A100 GPUs with a global batch size of 32\. Training is performed for at most 30 epochs using a learning rate of5×10−45\\times 10^\{\-4\}, with early stopping applied when the validation performance does not improve for five consecutive epochs\. During inference, peptide sequences are generated using beam search with a beam width of five\.
## Results and Analysis
### Main Results
Table 1:Amino\-acid\- and peptide\-level performance on the Nine\-Species, Seven\-Species, and HC\-PT datasets\. Best results are shown in bold\.The empirical results of our experiments are summarized in Tables[1](https://arxiv.org/html/2609.30542#Sx5.T1)and[2](https://arxiv.org/html/2609.30542#Sx5.T2)\. Overall, GyroNovo consistently outperforms the baseline methods across all datasets and evaluation metrics\. Table[1](https://arxiv.org/html/2609.30542#Sx5.T1)presents the amino acid\- and peptide\-level results\. At the amino acid level, GyroNovo achieves the highest precision and recall on all datasets and exceeds the strongest baseline \(LIPNovo\) by an average of 6\.4 percentage points in precision\. Similarly, at the peptide level, GyroNovo outperforms all baseline methods across all metrics and surpasses the baseline by an average of 8\.7 percentage points in precision\. These improvements are consistent across the three evaluated datasets\. Table[2](https://arxiv.org/html/2609.30542#Sx5.T2)presents the PTM\-level results and shows that GyroNovo outperforms the strongest baseline by an average of 5\.6 percentage points in precision\.
In addition to the preceding comparisons, we compare GyroNovo with the recently released,concurrentLIPNovo\+\([Du et al\. 2026](https://arxiv.org/html/2609.30542#bib.bib10)\)\. As shown in Tables[1](https://arxiv.org/html/2609.30542#Sx5.T1)and[2](https://arxiv.org/html/2609.30542#Sx5.T2), GyroNovo consistently outperforms LIPNovo\+ across all metrics on all three datasets\. In particular, GyroNovo improves amino\-acid\-level, peptide\-level, and PTM\-level precision by 1\.20, 2\.93, and 2\.10 average percentage points, respectively\.
Table 2:PTM\-level performance on the Nine\-Species, Seven\-Species, and HC\-PT datasets\. Best results are shown in bold\.
### Ablation Study and Analysis
#### Component Ablation\.
Table[3](https://arxiv.org/html/2609.30542#Sx5.T3)presents a cumulative ablation study of GyroNovo on the nine\-species dataset, where each component is added on top of all components included in the preceding rows\. Imputation\-loss reweighting provides only marginal changes, improving amino\-acid precision and recall by0\.40\.4and0\.20\.2percentage points, respectively, and peptide precision by0\.30\.3percentage points\. Adding augmented views yields the largest incremental gains, improving amino\-acid precision and recall by3\.13\.1percentage points each, and peptide precision by5\.35\.3percentage points\. Mass Rotary Attention further improves these metrics by1\.91\.9,1\.91\.9, and3\.53\.5percentage points, respectively\. Overall, the complete model improves upon the LIPNovo baseline by5\.45\.4and5\.25\.2percentage points in amino\-acid precision and recall, respectively, and by9\.19\.1percentage points in peptide precision\.
#### Mechanistic Analysis of Mass Rotary Attention\.
To assess whether Mass Rotary Attention improves access to chemically related peaks, we formulate spectrum self\-attention as a neighbor\-retrieval task\. For each spectrum, we match observed peaks to theoretical fragments of the ground\-truth peptide, and label pairs corresponding to adjacentbb\-ions oryy\-ions as true ladder neighbors\. For every eligible query peak, encoder layer, and attention head, we rank the remaining observed peaks using the query–key attention logits and compute Hit@KK, indicating whether at least one true ladder neighbors appear among theKKhighest\-ranked peaks\. We average scores across queries, heads, and layers within each spectrum and then macro\-average across spectra\. We compare the vanilla encoder with Mass Rotary Attention and include two controls: a within\-spectrum intervention that permutes only the mass assignments supplied to the rotary operation while preserving the peak embeddings and intensities, and an exact random\-ranking baseline accounting for the number of candidate and positive peaks\.
As shown in Figure[2](https://arxiv.org/html/2609.30542#Sx5.F2), Mass Rotary Attention consistently improves adjacent ion\-ladder neighbor retrieval across all values ofKK\. In particular, Hit@1 increases from approximately 13% with vanilla attention to 19% with Mass Rotary Attention, while Hit@20 improves from approximately 58% to 65%\. Shuffling the mass assignments used by the rotary operation substantially reduces retrieval performance, bringing it closer to random ranking and below the vanilla model\. These results show that the improvement arises from correctly encoding relative mass relationships, enabling self\-attention to prioritize peptide\-consistent peak neighbors\.
Figure 2:Adjacent ion\-ladder neighbor retrieval performance of vanilla and mass\-aware rotary attention on Nine\-Species dataset\.Table 3:Cumulative ablation study of GyroNovo on the Nine\-Species dataset\. Components are added sequentially to the reproduced baseline\. Best results are shown in bold\.Table 4:Training and inference efficiency of the baseline, GyroNovo\-Pairwise, and our GyroNovo\. Throughput and peak memory usage are normalized to the baseline; higher throughput and lower memory usage are better\.
#### Time and Memory Overhead\.
We compare the computational overhead of GyroNovo and LIPNovo in terms of throughput and peak memory usage\. We also include GyroNovo\-Pairwise in this comparison, a variant that replaces rotary mass embeddings with an additive attention bias computed from a precomputed pairwise mass\-difference matrix, following[Lapin et al\. \(2025\)](https://arxiv.org/html/2609.30542#bib.bib12)\. This variant allows us to directly evaluate whether Mass Rotary Attention represents mass relationships more efficiently than explicitly constructing and storing a pairwise mass\-difference matrix\. Table[4](https://arxiv.org/html/2609.30542#Sx5.T4)summarizes the results relative to LIPNovo\. During training, GyroNovo achieves\(0\.30×\)\(0\.30\\times\)the baseline throughput while using\(2\.08×\)\(2\.08\\times\)the baseline memory, compared with\(0\.27×\)\(0\.27\\times\)throughput and\(4\.02×\)\(4\.02\\times\)memory for GyroNovo\-Pairwise\. At inference, both variants maintain baseline throughput; however, GyroNovo uses only\(1\.01×\)\(1\.01\\times\)the baseline memory, whereas GyroNovo\-Pairwise requires\(3\.86×\)\(3\.86\\times\)\. Finally, these results demonstrate that rotary mass embeddings capture mass relationships substantially more memory\-efficiently than a precomputed pairwise mass\-difference matrix, particularly during inference\.
Both LIPNovo and GyroNovo converged at epoch 13 under the same early\-stopping protocol, despite a maximum training budget of 30 epochs\. Therefore, the observed improvements are not attributable to GyroNovo receiving more optimization epochs, although its multi\-view training still entails greater computation per epoch\.
## Conclusion
In this paper, we propose GyroNovo, which addresses two major limitations of existing Transformer\-basedde novopeptide sequencing methods: uniformly supervising the imputation of missing fragment ions and insufficiently modeling relative mass relationships between spectral peaks\. First, our decoder\-guided learning strategy uses current amino\-acid prediction difficulty to prioritize fragments that are most relevant to unresolved residues\. The same signal guides complementary easy and hard spectrum augmentations, providing clearer evidence for difficult residues while improving robustness to incomplete fragmentation ladders\. Second, continuous mass \(m/z\) rotary attention explicitly incorporates relative mass differences into spectrum self\-attention without introducing additional learned positional parameters\.
These components retain the latent\-imputation framework and require no theoretical spectrum or augmentation procedure at inference time\. Our results demonstrate that aligning fragment representation learning with decoder uncertainty, together with an appropriate relative\-mass inductive bias, provides a more effective approach to learning from sparse and noisy tandem mass spectra\. We hope this work inspires further investigation into task\-aware spectrum representation learning and physically informed attention mechanisms for database\-free peptide sequencing\.
## Acknowledgments
This work was supported by the UBC AI and Health Network at the University of British Columbia\. Muhammad Abdul\-Mageed also acknowledges support from Canada Research Chairs \(CRC\), Canada Foundation for Innovation \(CFI; 37771\), Digital Research Alliance of Canada,222https://alliancecan\.caand UBC ARC\-Sockeye\.
## References
- Chicket al\.\(2015\)J\. M\. Chick, D\. Kolippakkam, D\. P\. Nusinow, B\. Zhai, R\. Rad, E\. L\. Huttlin, and S\. P\. GygiA mass\-tolerant database search identifies a large proportion of unassigned spectra in shotgun proteomics as modified peptides\.Nature Biotechnology33\(7\),pp\. 743–749\.External Links:ISSN 1546\-1696,[Document](https://dx.doi.org/10.1038/nbt.3267),[Link](https://doi.org/10.1038/nbt.3267)Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p1.1)\.
- Duet al\.\(2026\)Y\. Du, Q\. Niu, C\. Yang, N\. Yu, and S\. WangLIPNovo\+: self\-reflective latent imputation for robust de novo peptide sequencing\.SSRN Electronic Journal\.Note:Posted June 6, 2026External Links:[Document](https://dx.doi.org/10.2139/ssrn.6890054),[Link](https://ssrn.com/abstract=6890054)Cited by:[Baselines](https://arxiv.org/html/2609.30542#Sx4.SSx2.p1.1),[Main Results](https://arxiv.org/html/2609.30542#Sx5.SSx1.p2.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.11.1)\.
- Duet al\.\(2025\)Y\. Du, C\. Yang, N\. Yu, W\. Lin, Q\. Zhao, and S\. WangLatent imputation before prediction: a new computational paradigm for de novo peptide sequencing\.InInternational Conference on Machine Learning,Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p2.1),[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p2.1),[Problem Formulation](https://arxiv.org/html/2609.30542#Sx3.SSx1.p3.1),[Latent Imputation Backbone](https://arxiv.org/html/2609.30542#Sx3.SSx2.p1.1),[Training Objective and Inference](https://arxiv.org/html/2609.30542#Sx3.SSx5.p2.1),[Baselines](https://arxiv.org/html/2609.30542#Sx4.SSx2.p1.1),[Implementation Details](https://arxiv.org/html/2609.30542#Sx4.SSx4.p1.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.10.1)\.
- Eloffet al\.\(2025\)K\. Eloff, K\. Kalogeropoulos, A\. Mabona, O\. Morell, R\. Catzel, E\. Rivera\-de\-Torre, J\. Berg Jespersen, W\. Williams, S\. P\. B\. van Beljouw, M\. J\. Skwark, A\. H\. Laustsen, S\. J\. J\. Brouns, A\. Ljungars, E\. M\. Schoof, J\. Van Goey, U\. auf dem Keller, K\. Beguir, N\. Lopez Carranza, and T\. P\. JenkinsInstaNovo enables diffusion\-powered de novo peptide sequencing in large\-scale proteomics experiments\.Nature Machine Intelligence7\(4\),pp\. 565–579\.External Links:ISSN 2522\-5839,[Document](https://dx.doi.org/10.1038/s42256-025-01019-5),[Link](https://doi.org/10.1038/s42256-025-01019-5)Cited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1),[Datasets](https://arxiv.org/html/2609.30542#Sx4.SSx1.p1.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.9.1)\.
- Frank and Pevzner \(2005\)A\. Frank and P\. PevznerPepNovo: de novo peptide sequencing via probabilistic network modeling\.Analytical Chemistry77\(4\),pp\. 964–973\.External Links:ISSN 0003\-2700,[Document](https://dx.doi.org/10.1021/ac048788h),[Link](https://doi.org/10.1021/ac048788h),https://pubs\.acs\.org/ancham/article\-pdf/77/4/964/37852098/ac048788h\.pdfCited by:[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1)\.
- HEet al\.\(2013\)L\. HE, X\. HAN, and B\. MADE novo sequencing with limited number of post\-translational modifications per peptide\.Journal of Bioinformatics and Computational Biology11\(04\),pp\. 1350007\.Note:PMID: 23859271External Links:[Document](https://dx.doi.org/10.1142/S0219720013500078),[Link](https://doi.org/10.1142/S0219720013500078),https://doi\.org/10\.1142/S0219720013500078Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p1.1)\.
- Jinet al\.\(2024\)Z\. Jin, S\. Xu, X\. Zhang, T\. Ling, N\. Dong, W\. Ouyang, Z\. Gao, C\. Chang, and S\. SunContraNovo: a contrastive learning approach to enhance de novo peptide sequencing\.Proceedings of the AAAI Conference on Artificial Intelligence38\(1\),pp\. 144–152\.External Links:[Document](https://dx.doi.org/10.1609/aaai.v38i1.27765),[Link](https://doi.org/10.1609/aaai.v38i1.27765)Cited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1)\.
- Klaproth\-Andradeet al\.\(2024\)D\. Klaproth\-Andrade, J\. Hingerl, Y\. Bruns, N\. H\. Smith, J\. Träuble, M\. Wilhelm, and J\. GagneurDeep learning\-driven fragment ion series classification enables highly precise and sensitive de novo peptide sequencing\.Nature Communications15\(1\),pp\. 151\.External Links:ISSN 2041\-1723,[Document](https://dx.doi.org/10.1038/s41467-023-44323-7),[Link](https://doi.org/10.1038/s41467-023-44323-7)Cited by:[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1)\.
- Lapinet al\.\(2025\)J\. Lapin, A\. Nilsson, M\. Wilhelm, and L\. KällPairwise attention: leveraging mass differences to enhance de novo sequencing of mass spectra\.Journal of Proteome Research24\(7\),pp\. 3722–3730\.External Links:ISSN 1535\-3893,[Document](https://dx.doi.org/10.1021/acs.jproteome.5c00063),[Link](https://doi.org/10.1021/acs.jproteome.5c00063),https://pubs\.acs\.org/jprobs/article\-pdf/24/7/3722/42334399/pr5c00063\.pdfCited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1),[Time and Memory Overhead\.](https://arxiv.org/html/2609.30542#Sx5.SSx2.SSSx3.p1.1)\.
- Maet al\.\(2003\)B\. Ma, K\. Zhang, C\. Hendrie, C\. Liang, M\. Li, A\. Doherty\-Kirby, and G\. LajoiePEAKS: powerful software for peptide de novo sequencing by tandem mass spectrometry\.Rapid Communications in Mass Spectrometry17\(20\),pp\. 2337–2342\.External Links:[Document](https://dx.doi.org/https%3A//doi.org/10.1002/rcm.1196),[Link](https://analyticalsciencejournals.onlinelibrary.wiley.com/doi/abs/10.1002/rcm.1196),https://analyticalsciencejournals\.onlinelibrary\.wiley\.com/doi/pdf/10\.1002/rcm\.1196Cited by:[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.4.1)\.
- Ma \(2015\)B\. MaNovor: real\-time peptide de novo sequencing software\.Journal of the American Society for Mass Spectrometry26\(11\),pp\. 1885–1894\.External Links:[Document](https://dx.doi.org/10.1007/s13361-015-1204-0),ISSN 1044\-0305Cited by:[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1)\.
- Maoet al\.\(2023\)Z\. Mao, R\. Zhang, L\. Xin, and M\. LiMitigating the missing\-fragmentation problem in de novo peptide sequencing with a two\-stage graph\-based deep learning model\.Nature Machine Intelligence5\(11\),pp\. 1250–1260\.External Links:ISSN 2522\-5839,[Document](https://dx.doi.org/10.1038/s42256-023-00738-x),[Link](https://doi.org/10.1038/s42256-023-00738-x)Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p2.1),[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1)\.
- McDonnellet al\.\(2022\)K\. McDonnell, E\. Howley, and F\. AbramThe impact of noise and missing fragmentation cleavages on <i\>de novo</i\> peptide identification algorithms\.Computational and Structural Biotechnology Journal20\(\),pp\. 1402–1412\.External Links:[Document](https://dx.doi.org/10.1016/j.csbj.2022.03.008),[Link](https://spj.science.org/doi/abs/10.1016/j.csbj.2022.03.008),https://spj\.science\.org/doi/pdf/10\.1016/j\.csbj\.2022\.03\.008Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p1.1)\.
- Qiaoet al\.\(2021\)R\. Qiao, N\. H\. Tran, L\. Xin, X\. Chen, M\. Li, B\. Shan, and A\. GhodsiComputationally instrument\-resolution\-independent de novo peptide sequencing for high\-resolution devices\.Nature Machine Intelligence3\(5\),pp\. 420–425\.External Links:ISSN 2522\-5839,[Document](https://dx.doi.org/10.1038/s42256-021-00304-3),[Link](https://doi.org/10.1038/s42256-021-00304-3)Cited by:[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.6.1)\.
- Suet al\.\(2024\)J\. Su, M\. Ahmed, Y\. Lu, S\. Pan, W\. Bo, and Y\. LiuRoFormer: enhanced transformer with rotary position embedding\.Neurocomputing568,pp\. 127063\.External Links:ISSN 0925\-2312,[Document](https://dx.doi.org/https%3A//doi.org/10.1016/j.neucom.2023.127063),[Link](https://www.sciencedirect.com/science/article/pii/S0925231223011864)Cited by:[Continuous Mass Rotary Attention](https://arxiv.org/html/2609.30542#Sx3.SSx4.p2.1)\.
- Tranet al\.\(2019\)N\. H\. Tran, R\. Qiao, L\. Xin, X\. Chen, C\. Liu, X\. Zhang, B\. Shan, A\. Ghodsi, and M\. LiDeep learning enables de novo peptide sequencing from data\-independent\-acquisition mass spectrometry\.Nature Methods16\(1\),pp\. 63–66\.External Links:ISSN 1548\-7105,[Document](https://dx.doi.org/10.1038/s41592-018-0260-3),[Link](https://doi.org/10.1038/s41592-018-0260-3)Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p1.1)\.
- Tranet al\.\(2017\)N\. H\. Tran, X\. Zhang, L\. Xin, B\. Shan, and M\. LiDe novo peptide sequencing by deep learning\.Proceedings of the National Academy of Sciences114\(31\),pp\. 8247–8252\.External Links:[Document](https://dx.doi.org/10.1073/pnas.1705691114),[Link](https://www.pnas.org/doi/abs/10.1073/pnas.1705691114),https://www\.pnas\.org/doi/pdf/10\.1073/pnas\.1705691114Cited by:[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1),[Datasets](https://arxiv.org/html/2609.30542#Sx4.SSx1.p1.1),[Baselines](https://arxiv.org/html/2609.30542#Sx4.SSx2.p1.1),[Evaluation Metrics](https://arxiv.org/html/2609.30542#Sx4.SSx3.p1.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.5.1)\.
- Vaswaniet al\.\(2017\)A\. Vaswani, N\. Shazeer, N\. Parmar, J\. Uszkoreit, L\. Jones, A\. N\. Gomez, Ł\. Kaiser, and I\. PolosukhinAttention is all you need\.InProceedings of the 31st International Conference on Neural Information Processing Systems,NIPS’17,Red Hook, NY, USA,pp\. 6000–6010\.External Links:ISBN 9781510860964Cited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1)\.
- Xiaet al\.\(2024\)J\. Xia, S\. Chen, J\. Zhou, X\. Shan, W\. Du, Z\. Gao, C\. Tan, B\. Hu, J\. Zheng, and S\. Z\. LiAdaNovo: towards robust de novo peptide sequencing in proteomics against data biases\.InThe Thirty\-eighth Annual Conference on Neural Information Processing Systems,External Links:[Link](https://openreview.net/forum?id=0zfUiSX5si)Cited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1),[Baselines](https://arxiv.org/html/2609.30542#Sx4.SSx2.p1.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.8.1)\.
- Xuet al\.\(2023\)X\. Xu, C\. Yang, Q\. He, K\. Shu, Y\. Xinpu, Z\. Chen, Y\. Zhu, and T\. ChenPGPointNovo: an efficient neural network\-based tool for parallel de novo peptide sequencing\.Bioinformatics Advances3\(1\),pp\. vbad057\.External Links:ISSN 2635\-0041,[Document](https://dx.doi.org/10.1093/bioadv/vbad057),[Link](https://doi.org/10.1093/bioadv/vbad057),https://academic\.oup\.com/bioinformaticsadvances/article\-pdf/3/1/vbad057/50146875/vbad057\.pdfCited by:[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1),[Baselines](https://arxiv.org/html/2609.30542#Sx4.SSx2.p1.1)\.
- Yanget al\.\(2024\)T\. Yang, T\. Ling, B\. Sun, Z\. Liang, F\. Xu, X\. Huang, L\. Xie, Y\. He, L\. Li, F\. He, Y\. Wang, and C\. ChangIntroducingπ\\pi\-helixnovo for practical large\-scale de novo peptide sequencing\.Briefings in Bioinformatics25\(2\),pp\. bbae021\.External Links:ISSN 1477\-4054,[Document](https://dx.doi.org/10.1093/bib/bbae021),[Link](https://doi.org/10.1093/bib/bbae021),https://academic\.oup\.com/bib/article\-pdf/25/2/bbae021/57169420/bbae021\.pdfCited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1),[Implementation Details](https://arxiv.org/html/2609.30542#Sx4.SSx4.p1.1)\.
- Yilmazet al\.\(2024\)M\. Yilmaz, W\. E\. Fondrie, W\. Bittremieux, C\. F\. Melendez, R\. Nelson, V\. Ananth, S\. Oh, and W\. S\. NobleSequence\-to\-sequence translation from mass spectra to peptides with a transformer model\.Nature Communications15\(1\),pp\. 6427\.External Links:ISSN 2041\-1723,[Document](https://dx.doi.org/10.1038/s41467-024-49731-x),[Link](https://doi.org/10.1038/s41467-024-49731-x)Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p1.1),[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1),[Problem Formulation](https://arxiv.org/html/2609.30542#Sx3.SSx1.p3.1),[Baselines](https://arxiv.org/html/2609.30542#Sx4.SSx2.p1.1),[Implementation Details](https://arxiv.org/html/2609.30542#Sx4.SSx4.p1.1),[Table 1](https://arxiv.org/html/2609.30542#Sx5.T1.1.1.1.1.1.1.1.7.1)\.
- Zhanget al\.\(2025\)X\. Zhang, T\. Ling, Z\. Jin, S\. Xu, Z\. Gao, B\. Sun, Z\. Qiu, J\. Wei, N\. Dong, G\. Wang, G\. Wang, L\. Li, M\. Abdul\-Mageed, L\. V\. S\. Lakshmanan, F\. He, W\. Ouyang, C\. Chang, and S\. Sunπ\\pi\-PrimeNovo: an accurate and efficient non\-autoregressive deep learning model for de novo peptide sequencing\.Nature Communications16\(1\),pp\. 267\.External Links:ISSN 2041\-1723,[Document](https://dx.doi.org/10.1038/s41467-024-55021-3),[Link](https://doi.org/10.1038/s41467-024-55021-3)Cited by:[Transformer\-basedDe NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx2.p1.1)\.
- Zhouet al\.\(2024\)J\. Zhou, S\. Chen, J\. Xia, S\. Liu, T\. Ling, W\. Du, Y\. Liu, J\. Yin, and S\. Z\. LiNovoBench: benchmarking deep learning\-based de novo sequencing methods in proteomics\.InThe Thirty\-eight Conference on Neural Information Processing Systems Datasets and Benchmarks Track,External Links:[Link](https://openreview.net/forum?id=RQlbMrA5XL)Cited by:[Introduction](https://arxiv.org/html/2609.30542#Sx1.p1.1),[De NovoPeptide Sequencing](https://arxiv.org/html/2609.30542#Sx2.SSx1.p1.1),[Datasets](https://arxiv.org/html/2609.30542#Sx4.SSx1.p1.1)\.
## Appendix ADecoder\-Conditioned Spectrum Augmentation
##### Decoder\-derived conditioning\.
Conditioning uses the current decoder’s teacher\-forced predictions on the observed\-spectrum branch\. The token losses used for conditioning are ordinary cross\-entropy losses without label smoothing; padding and end\-of\-sequence tokens are excluded\. The sample\-level scalesis\_\{i\}controls overall corruption strength, while fragment\-level weights determine its spatial distribution along the fragmentation ladders\.
For theoretical fragmentkk, let𝒩ik\\mathcal\{N\}\_\{ik\}contain the two residues adjacent to its cleavage site and their immediate neighbors, restricted to valid peptide positions\. This corresponds to a neighborhood radius of one around each cleavage\-adjacent residue\. The mapping accounts for the ion family and the decoder’s sequence direction\. We define
wik=clip\(1\+βmaxl∈𝒩ikeil,wmin,wmax\),w\_\{ik\}=\\operatorname\{clip\}\\left\(1\+\\beta\\max\_\{l\\in\\mathcal\{N\}\_\{ik\}\}e\_\{il\},\\,w\_\{\\min\},\\,w\_\{\\max\}\\right\),whereβ=1\\beta=1,wmin=1w\_\{\\min\}=1, andwmax=2w\_\{\\max\}=2\. Taking the maximum ensures that a single difficult residue can emphasize nearby fragments without being diluted by easier neighboring residues\.
The augmentation focus is
fik=\{wik−1maxk′\(wik′−1\),maxk′\(wik′−1\)\>10−8,0,otherwise\.f\_\{ik\}=\\begin\{cases\}\\dfrac\{w\_\{ik\}\-1\}\{\\max\_\{k^\{\\prime\}\}\(w\_\{ik^\{\\prime\}\}\-1\)\},&\\max\_\{k^\{\\prime\}\}\(w\_\{ik^\{\\prime\}\}\-1\)\>10^\{\-8\},\\\\\[6\.0pt\] 0,&\\text\{otherwise\}\.\\end\{cases\}Accordingly, the most emphasized fragment in a spectrum has focus one whenever the conditioning signal is nonzero\. This normalization makes focus a relative measure within each spectrum; the sample scalesis\_\{i\}separately captures absolute decoding difficulty\. Decoder\-derived weights are computed for the theoretical targets retained by the imputation branch, which contains at most 100 targets in our configuration\. Additional theoretical fragments receive the neutral weightwik=1w\_\{ik\}=1\.
##### Curriculum and view\-specific strength\.
Lett∈\[0,1\]t\\in\[0,1\]denote training progress, computed as the zero\-based current epoch divided by the configured maximum number of epochs\. We represent corruption strength by
𝐝=\(dmiss,db,dy,dgap,σI,σm,λN,dPTM\),\\mathbf\{d\}=\(d\_\{\\mathrm\{miss\}\},d\_\{b\},d\_\{y\},d\_\{\\mathrm\{gap\}\},\\sigma\_\{I\},\\sigma\_\{m\},\\lambda\_\{N\},d\_\{\\mathrm\{PTM\}\}\),corresponding to general peak missingness,bb\-ion dropout,yy\-ion dropout, ladder\-gap probability, intensity jitter,m/zm/zjitter, structured\-noise count, and PTM visibility reduction\. The curriculum interpolates linearly from the zero vector to
𝐝max=\(0\.5,0\.5,0\.5,0\.8,0\.5,0\.2,24,0\.8\)\.\\mathbf\{d\}\_\{\\max\}=\(0\.5,\\,0\.5,\\,0\.5,\\,0\.8,\\,0\.5,\\,0\.2,\\,24,\\,0\.8\)\.For viewvv, the effective strength is
𝐝i\(v\)=clip\[0,𝐝max\]\(siavt\(𝐦i⊙𝐝max\)\),\\mathbf\{d\}\_\{i\}^\{\(v\)\}=\\operatorname\{clip\}\_\{\[0,\\mathbf\{d\}\_\{\\max\}\]\}\\left\(s\_\{i\}\\,a\_\{v\}\\,t\\,\\bigl\(\\mathbf\{m\}\_\{i\}\\odot\\mathbf\{d\}\_\{\\max\}\\bigr\)\\right\),whereaeasy=0\.5a\_\{\\mathrm\{easy\}\}=0\.5,ahard=1a\_\{\\mathrm\{hard\}\}=1, and clipping is componentwise\. The multiplier𝐦i\\mathbf\{m\}\_\{i\}equals one except for its PTM component, which is1\.21\.2for peptides containing modified residues\. These entries are operator strengths: depending on the operator, they act as probability penalties, Gaussian standard deviations, or Poisson count parameters\.
The fixed retention and insertion probabilities and decoder\-focus adjustments remain active even whent=0t=0\. The curriculum therefore increases additional corruption over training rather than delaying all augmentation until later epochs\.
##### Matching and observed\-peak retention\.
Each theoretical fragment is independently assigned to its nearest observed peak if their absolutem/zm/zdifference is at most0\.50\.5; matching is not constrained to be one\-to\-one\. When multiple fragments match the same observed peak, the peak receives their maximum focus score\. Unannotated observed peaks receive focus zero\.
For observed peakjj, letFijF\_\{ij\}denote its focus, and letBijB\_\{ij\},YijY\_\{ij\}, andPijP\_\{ij\}indicate its assignedbb\-ion,yy\-ion, and modification\-containing annotations\. These indicators are zero for unannotated peaks\. Before applying ladder gaps, the retention probability is
pijkeep,v=clip\[0,1\]\[\\displaystyle p\_\{ij\}^\{\\mathrm\{keep\},v\}=\\operatorname\{clip\}\_\{\[0,1\]\}\\Bigl\[pij0−0\.20dmiss\\displaystyle p\_\{ij\}^\{0\}\-0\.20d\_\{\\mathrm\{miss\}\}−0\.20Bijdb−0\.20Yijdy\\displaystyle\-0\.20B\_\{ij\}d\_\{b\}\-0\.20Y\_\{ij\}d\_\{y\}−0\.10PijdPTM\+ηv0\.40Fij\],\\displaystyle\-0\.10P\_\{ij\}d\_\{\\mathrm\{PTM\}\}\+\\eta\_\{v\}\\,0\.40F\_\{ij\}\\Bigr\],
wherepij0=0\.90p\_\{ij\}^\{0\}=0\.90for matched peaks and0\.700\.70for unannotated peaks, andηeasy=\+1\\eta\_\{\\mathrm\{easy\}\}=\+1andηhard=−1\\eta\_\{\\mathrm\{hard\}\}=\-1\. Here and below, operator strengths are those of the corresponding sample and view\.
Retained observed peaks preserve their measuredm/zm/zvalues\. Their intensities are adjusted as
Iij\(v\)=Iij\(1\+ηv0\.20Fij\)\.I\_\{ij\}^\{\(v\)\}=I\_\{ij\}\\left\(1\+\\eta\_\{v\}\\,0\.20F\_\{ij\}\\right\)\.Thus, maximally focused observed peaks are strengthened by20%20\\%in the easy view and attenuated by20%20\\%in the hard view\. Unannotated observed peaks retain their original intensities\. No additional independent focus\-based forced\-drop stage is used\.
##### Insertion of missing or dropped fragments\.
A theoretical fragment is eligible for insertion if it has no observed match or its matched observed peak was removed\. Its insertion probability is
pikadd,v=clip\[0,0\.90\]\[\\displaystyle p\_\{ik\}^\{\\mathrm\{add\},v\}=\\operatorname\{clip\}\_\{\[0,0\.90\]\}\\Bigl\[0\.30\\displaystyle 0\.30\+ηv0\.30fik−0\.15dmiss\\displaystyle\+\\eta\_\{v\}\\,0\.30f\_\{ik\}\-0\.15d\_\{\\mathrm\{miss\}\}−0\.15PikdPTM\],\\displaystyle\-0\.15P\_\{ik\}d\_\{\\mathrm\{PTM\}\}\\Bigr\],
subject to the ladder\-gap restrictions below\. The insertion operation can therefore restore both originally missing evidence and evidence removed earlier in view construction\.
For a fragment of ordinal indexr∈\{1,…,Li−1\}r\\in\\\{1,\\ldots,L\_\{i\}\-1\\\}, the baseline synthetic intensity is
I^ik\\displaystyle\\widehat\{I\}\_\{ik\}=Ri\[0\.75\+0\.35\(1−2\|rLi−1−12\|\)\],\\displaystyle=R\_\{i\}\\biggl\[0\.75\+0\.35\\biggl\(1\-2\\,\\biggl\|\\frac\{r\}\{L\_\{i\}\-1\}\-\\frac\{1\}\{2\}\\biggr\|\\biggr\)\\biggr\],whereRi=max\(1,maxjIij\)\.\\displaystyle\\quad\\text\{where \}R\_\{i\}=\\max\\Bigl\(1,\\,\\max\_\{j\}I\_\{ij\}\\Bigr\)\.The same baseline profile is used forbb\- andyy\-ions, with no additional PTM\-specific intensity boost\. Inserted fragment intensities and locations are perturbed according to
Iikadd=I^ikmax\(0\.05,1\+ϵI\),mikadd=miktheory\+ϵm,I\_\{ik\}^\{\\mathrm\{add\}\}=\\widehat\{I\}\_\{ik\}\\max\(0\.05,1\+\\epsilon\_\{I\}\),\\qquad m\_\{ik\}^\{\\mathrm\{add\}\}=m\_\{ik\}^\{\\mathrm\{theory\}\}\+\\epsilon\_\{m\},whereϵI∼𝒩\(0,σI2\)\\epsilon\_\{I\}\\sim\\mathcal\{N\}\(0,\\sigma\_\{I\}^\{2\}\)andϵm∼𝒩\(0,σm2\)\\epsilon\_\{m\}\\sim\\mathcal\{N\}\(0,\\sigma\_\{m\}^\{2\}\)\.
##### Contiguous ladder gaps\.
A ladder\-gap event is sampled with probabilitydgapd\_\{\\mathrm\{gap\}\}\. When activated, its run length is sampled uniformly from
\{1,…,min\(Li−1,1\+round\(3dgap\)\)\},\\left\\\{1,\\ldots,\\min\\left\(L\_\{i\}\-1,\\,1\+\\operatorname\{round\}\(3d\_\{\\mathrm\{gap\}\}\)\\right\)\\right\\\},and a valid starting index is chosen uniformly\. Under the configured maximumdgap=0\.8d\_\{\\mathrm\{gap\}\}=0\.8, this produces runs of at most three consecutive fragment indices\. The selected ordinal indices are applied to both ion ladders, groupingbrb\_\{r\}andyry\_\{r\}by fragment index\.
In the hard view, observed peaks assigned to the selected indices are removed and theoretical insertion at those indices is disabled\. In the easy view, matched peaks with focus at least0\.50\.5are exempt from this additional ladder\-gap removal, although they remain subject to the ordinary retention draw\. Insertion within an easy\-view gap is permitted only for fragments with focus strictly greater than0\.50\.5\.
##### Structured noise and final spectrum assembly\.
In addition to retaining a subset of unannotated observed peaks, we sample
Ni\(v\)∼Poisson\(0\.35λN\)N\_\{i\}^\{\(v\)\}\\sim\\operatorname\{Poisson\}\(0\.35\\lambda\_\{N\}\)candidate structured\-noise peaks\. Each candidate is generated from a sampled theoretical parent fragment by subtracting an offset of either18\.01056518\.010565or17\.02654917\.026549, corresponding to water\- or ammonia\-loss masses, and adding Gaussianm/zm/zjitter with standard deviation0\.010\.01\. Its intensity is sampled uniformly between0\.080\.08and0\.350\.35times the parent’s baseline synthetic intensity and capped at0\.35Ri0\.35R\_\{i\}\. Noise candidates belowm/z=50m/z=50are discarded\.
Retained observed peaks, inserted theoretical fragments, and structured\-noise peaks are combined\. Invalid or nonpositive entries are removed, the 150 most intense peaks are retained if necessary, and the resulting spectrum is sorted bym/zm/z\. If no peak survives, the strongest baseline theoretical fragment is retained to ensure a nonempty input\.
##### Supervision and hyperparameter selection\.
Both views retain the original precursor information and ground\-truth peptide label\. They are processed by the same encoder, imputation module, and decoder as the observed\-spectrum branch\. The fragment weightswikw\_\{ik\}also emphasize difficult targets in the imputation objective\. Easy\- and hard\-view imputation losses are averaged with equal weights, as are their sequence losses\. Conditioning signals are recomputed from the current observed\-spectrum predictions during training and are treated as constants when constructing the augmented views\.
The augmentation hyperparameters reported above were selected by random search\. Error analysis on the nine\-species validation set, held separate from the test set, motivated the perturbation types; their numerical strengths and probability coefficients are tunable augmentation parameters\.相似文章
@pushmeet:今天,我们的团队发布了新系统 AlphaProtein Novo 的预印本,用于生成式从头酶设计。在 @GoogleDeepMind……
Google DeepMind 与其合作者发布了 AlphaProtein Novo,这是一个用于生成式从头酶设计的机器学习流程,能够产生自然界中不存在的酶,并在基准反应中达到最先进活性水平,其中包括用于合成哌啶的酶以及用于降解环境毒素 DEHP 的酶。
ProHiFlo:具有功能引导的分层流匹配用于从头蛋白质生成
介绍了ProHiFlo,一种用于从头蛋白质生成的分层流匹配框架,具有从粗到细的生成、功能引导和SE(3)等变架构,以4倍更少的采样步骤实现了最先进的性能。
AgForce 实现抗原条件生成式抗体设计
本文识别了现有抗体设计方法中的三种失败模式(抗原盲区、词汇崩溃、收敛到边缘分布),并提出 AgForce,一种使用图神经网络和混合密度网络的新型编码器-解码器架构,在 Chimera-Bench 基准测试上实现了最先进的结合质量和序列恢复。
引导蛋白质语言模型中的离群流形塌缩
本文识别了引导蛋白质语言模型中一种名为“离群流形塌缩”的故障模式,即生成序列出现退化,并提出马氏过滤(Mahalanobis filtering)作为一种无需训练的事后处理方法,以提升序列质量和结构合理性。
DeepRHP:一种用于设计随机杂聚物作为蛋白质模拟物的混合变分自编码器
DeepRHP是一种混合变分自编码器,能够指导随机杂聚物作为蛋白质模拟物的设计,并通过稳定非天然环境中的膜蛋白(如Aquaporin Z)得到了验证。