public static void TestAnalysisEngineTests()
{
List <DigestionMotif> motifs = new List <DigestionMotif> {
new DigestionMotif("K", null, 1, null)
};
Protease protease = new Protease("Custom Protease5", CleavageSpecificity.Full, null, null, motifs);
ProteaseDictionary.Dictionary.Add(protease.Name, protease);
CommonParameters CommonParameters = new CommonParameters(
digestionParams: new DigestionParams(
protease: protease.Name,
maxMissedCleavages: 0,
minPeptideLength: 1,
maxModificationIsoforms: 1042),
scoreCutoff: 1,
productMassTolerance: new PpmTolerance(10));
var fsp = new List <(string fileName, CommonParameters fileSpecificParameters)>();
fsp.Add(("", CommonParameters));
List <Modification> localizeableModifications = new List <Modification>();
List <Modification> variableModifications = new List <Modification>();
List <Modification> fixedModifications = new List <Modification>();
Dictionary <Modification, ushort> modsDictionary = new Dictionary <Modification, ushort>();
foreach (var mod in fixedModifications)
{
modsDictionary.Add(mod, 0);
}
int i = 1;
foreach (var mod in variableModifications)
{
modsDictionary.Add(mod, (ushort)i);
i++;
}
foreach (var mod in localizeableModifications)
{
modsDictionary.Add(mod, (ushort)i);
i++;
}
var proteinList = new List <Protein> {
new Protein("MNNNKQQQ", "accession")
};
var modPep = proteinList.First().Digest(CommonParameters.DigestionParams, fixedModifications, variableModifications).Last();
HashSet <PeptideWithSetModifications> value1 = new HashSet <PeptideWithSetModifications> {
modPep
};
PeptideWithSetModifications compactPeptide1 = value1.First();
Assert.AreEqual("QQQ", value1.First().BaseSequence);
var modPep2 = proteinList.First().Digest(CommonParameters.DigestionParams, fixedModifications, variableModifications).First();
HashSet <PeptideWithSetModifications> value2 = new HashSet <PeptideWithSetModifications> {
modPep2
};
PeptideWithSetModifications compactPeptide2 = value2.First();
Assert.AreEqual("MNNNK", value2.First().BaseSequence);
var modPep3 = proteinList.First().Digest(CommonParameters.DigestionParams, fixedModifications, variableModifications).ToList()[1];
HashSet <PeptideWithSetModifications> value3 = new HashSet <PeptideWithSetModifications> {
modPep3
};
PeptideWithSetModifications compactPeptide3 = value3.First();
Assert.AreEqual("NNNK", value3.First().BaseSequence);
Ms2ScanWithSpecificMass scanA = new Ms2ScanWithSpecificMass(new MsDataScan(new MzSpectrum(new double[] { 1 }, new double[] { 1 }, false), 2, 1, true, Polarity.Positive, double.NaN, null, null, MZAnalyzerType.Orbitrap, double.NaN, null, null, "scan=1", double.NaN, null, null, double.NaN, null, DissociationType.AnyActivationType, 1, null), 1, 1, null, new CommonParameters());
Ms2ScanWithSpecificMass scanB = new Ms2ScanWithSpecificMass(new MsDataScan(new MzSpectrum(new double[] { 1 }, new double[] { 1 }, false), 3, 1, true, Polarity.Positive, double.NaN, null, null, MZAnalyzerType.Orbitrap, double.NaN, null, null, "scan=2", double.NaN, null, null, double.NaN, null, DissociationType.AnyActivationType, 1, null), 2 + 132.040, 1, null, new CommonParameters());
Ms2ScanWithSpecificMass scanC = new Ms2ScanWithSpecificMass(new MsDataScan(new MzSpectrum(new double[] { 1 }, new double[] { 1 }, false), 4, 1, true, Polarity.Positive, double.NaN, null, null, MZAnalyzerType.Orbitrap, double.NaN, null, null, "scan=3", double.NaN, null, null, double.NaN, null, DissociationType.AnyActivationType, 1, null), 3, 1, null, new CommonParameters());
PeptideSpectralMatch matchA = new PeptideSpectralMatch(compactPeptide1, 0, 0, 0, scanA, CommonParameters, new List <MatchedFragmentIon>());
PeptideSpectralMatch matchB = new PeptideSpectralMatch(compactPeptide2, 0, 0, 0, scanB, CommonParameters, new List <MatchedFragmentIon>());
PeptideSpectralMatch matchC = new PeptideSpectralMatch(compactPeptide3, 0, 0, 0, scanC, CommonParameters, new List <MatchedFragmentIon>());
var newPsms = new List <PeptideSpectralMatch> {
matchA, matchB, matchC
};
MsDataFile myMsDataFile = new TestDataFile(new List <PeptideWithSetModifications> {
value1.First(), value2.First(), value3.First()
});
var searchMode = new SinglePpmAroundZeroSearchMode(5);
Action <List <PeptideSpectralMatch>, string, List <string> > action2 = (List <PeptideSpectralMatch> l, string s, List <string> sdf) => {; };
Tolerance DeconvolutionMassTolerance = new PpmTolerance(5);
var arrayOfMs2ScansSortedByMass = MetaMorpheusTask.GetMs2Scans(myMsDataFile, null, new CommonParameters()).OrderBy(b => b.PrecursorMass).ToArray();
Action <BinTreeStructure, string> action1 = (BinTreeStructure l, string s) =>
{
Assert.AreEqual(1, l.FinalBins.Count);
};
FdrAnalysisEngine engine = new FdrAnalysisEngine(newPsms, searchMode.NumNotches, CommonParameters, fsp, new List <string> {
"ff"
});
engine.Run();
}
public static void TestModificationAnalysis()
{
IScan scan = new ThisTestScan();
ModificationMotif.TryGetMotif("N", out ModificationMotif motif1);
Modification mod1 = new Modification(_originalId: "mod1", _modificationType: "myModType", _target: motif1, _locationRestriction: "Anywhere.", _monoisotopicMass: 10);
ModificationMotif.TryGetMotif("L", out ModificationMotif motif2);
Modification mod2 = new Modification(_originalId: "mod2", _modificationType: "myModType", _target: motif2, _locationRestriction: "Anywhere.", _monoisotopicMass: 10);
IDictionary <int, List <Modification> > oneBasedModifications = new Dictionary <int, List <Modification> >
{
{ 2, new List <Modification> {
mod1
} },
{ 5, new List <Modification> {
mod2
} },
{ 7, new List <Modification> {
mod1
} },
};
Protein protein1 = new Protein("MNLDLDNDL", "prot1", oneBasedModifications: oneBasedModifications);
Dictionary <int, Modification> allModsOneIsNterminus1 = new Dictionary <int, Modification>
{
{ 2, mod1 },
};
PeptideWithSetModifications pwsm1 = new PeptideWithSetModifications(protein1, new DigestionParams(), 2, 9, CleavageSpecificity.Unknown, null, 0, allModsOneIsNterminus1, 0);
Dictionary <int, Modification> allModsOneIsNterminus2 = new Dictionary <int, Modification>
{
{ 2, mod1 },
{ 7, mod1 },
};
PeptideWithSetModifications pwsm2 = new PeptideWithSetModifications(protein1, new DigestionParams(), 2, 9, CleavageSpecificity.Unknown, null, 0, allModsOneIsNterminus2, 0);
Dictionary <int, Modification> allModsOneIsNterminus3 = new Dictionary <int, Modification>
{
{ 7, mod1 },
};
PeptideWithSetModifications pwsm3 = new PeptideWithSetModifications(protein1, new DigestionParams(), 2, 9, CleavageSpecificity.Unknown, null, 0, allModsOneIsNterminus3, 0);
Dictionary <int, Modification> allModsOneIsNterminus4 = new Dictionary <int, Modification>
{
{ 8, mod1 },
};
PeptideWithSetModifications pwsm4 = new PeptideWithSetModifications(protein1, new DigestionParams(), 1, 9, CleavageSpecificity.Unknown, null, 0, allModsOneIsNterminus4, 0);
CommonParameters CommonParameters = new CommonParameters(
digestionParams: new DigestionParams(
maxMissedCleavages: 0,
minPeptideLength: 1,
maxModificationIsoforms: int.MaxValue),
scoreCutoff: 1);
var newPsms = new List <PeptideSpectralMatch>
{
new PeptideSpectralMatch(pwsm1, 0, 10, 0, scan, CommonParameters.DigestionParams, new List <MatchedFragmentIon>()),
new PeptideSpectralMatch(pwsm1, 0, 10, 0, scan, CommonParameters.DigestionParams, new List <MatchedFragmentIon>()),
new PeptideSpectralMatch(pwsm2, 0, 10, 0, scan, CommonParameters.DigestionParams, new List <MatchedFragmentIon>()),
new PeptideSpectralMatch(pwsm3, 0, 10, 0, scan, CommonParameters.DigestionParams, new List <MatchedFragmentIon>()),
new PeptideSpectralMatch(pwsm4, 0, 10, 0, scan, CommonParameters.DigestionParams, new List <MatchedFragmentIon>()),
};
foreach (var psm in newPsms)
{
psm.ResolveAllAmbiguities();
}
MassDiffAcceptor searchMode = new SinglePpmAroundZeroSearchMode(5);
List <Protein> proteinList = new List <Protein> {
protein1
};
FdrAnalysisEngine fdrAnalysisEngine = new FdrAnalysisEngine(newPsms, searchMode.NumNotches, CommonParameters, new List <string>());
fdrAnalysisEngine.Run();
ModificationAnalysisEngine modificationAnalysisEngine = new ModificationAnalysisEngine(newPsms, new CommonParameters(), new List <string>());
var res = (ModificationAnalysisResults)modificationAnalysisEngine.Run();
Assert.AreEqual(2, res.CountOfEachModSeenOnProteins.Count());
Assert.AreEqual(2, res.CountOfEachModSeenOnProteins[mod1.IdWithMotif]);
Assert.AreEqual(1, res.CountOfEachModSeenOnProteins[mod2.IdWithMotif]);
Assert.AreEqual(1, res.CountOfModsSeenAndLocalized.Count());
Assert.AreEqual(2, res.CountOfModsSeenAndLocalized[mod1.IdWithMotif]);
Assert.AreEqual(0, res.CountOfAmbiguousButLocalizedModsSeen.Count());
Assert.AreEqual(0, res.CountOfUnlocalizedMods.Count());
Assert.AreEqual(0, res.CountOfUnlocalizedFormulas.Count());
}
public static void FdrTestMethod()
{
MassDiffAcceptor searchModes = new DotMassDiffAcceptor(null, new List <double> {
0, 1.0029
}, new PpmTolerance(5));
List <string> nestedIds = new List <string>();
Protein p = new Protein("MNKNNKNNNKNNNNK", null);
DigestionParams digestionParams = new DigestionParams();
var digested = p.Digest(digestionParams, new List <Modification>(), new List <Modification>()).ToList();
PeptideWithSetModifications pep1 = digested[0];
PeptideWithSetModifications pep2 = digested[1];
PeptideWithSetModifications pep3 = digested[2];
PeptideWithSetModifications pep4 = digested[3];
TestDataFile t = new TestDataFile(new List <PeptideWithSetModifications> {
pep1, pep2, pep3
});
MsDataScan mzLibScan1 = t.GetOneBasedScan(2);
Ms2ScanWithSpecificMass scan1 = new Ms2ScanWithSpecificMass(mzLibScan1, pep1.MonoisotopicMass.ToMz(1), 1, null, new CommonParameters());
PeptideSpectralMatch psm1 = new PeptideSpectralMatch(pep1, 0, 3, 0, scan1, digestionParams, new List <MatchedFragmentIon>());
MsDataScan mzLibScan2 = t.GetOneBasedScan(4);
Ms2ScanWithSpecificMass scan2 = new Ms2ScanWithSpecificMass(mzLibScan2, pep2.MonoisotopicMass.ToMz(1), 1, null, new CommonParameters());
PeptideSpectralMatch psm2 = new PeptideSpectralMatch(pep2, 1, 2, 1, scan2, digestionParams, new List <MatchedFragmentIon>());
MsDataScan mzLibScan3 = t.GetOneBasedScan(6);
Ms2ScanWithSpecificMass scan3 = new Ms2ScanWithSpecificMass(mzLibScan3, pep3.MonoisotopicMass.ToMz(1), 1, null, new CommonParameters());
PeptideSpectralMatch psm3 = new PeptideSpectralMatch(pep3, 0, 1, 2, scan3, digestionParams, new List <MatchedFragmentIon>());
psm3.AddOrReplace(pep4, 1, 1, true, new List <MatchedFragmentIon>(), 0);
var newPsms = new List <PeptideSpectralMatch> {
psm1, psm2, psm3
};
foreach (PeptideSpectralMatch psm in newPsms)
{
psm.ResolveAllAmbiguities();
}
FdrAnalysisEngine fdr = new FdrAnalysisEngine(newPsms, searchModes.NumNotches, new CommonParameters(), nestedIds);
fdr.Run();
Assert.AreEqual(2, searchModes.NumNotches);
Assert.AreEqual(0, newPsms[0].FdrInfo.CumulativeDecoyNotch);
Assert.AreEqual(1, newPsms[0].FdrInfo.CumulativeTargetNotch);
Assert.AreEqual(0, newPsms[1].FdrInfo.CumulativeDecoyNotch);
Assert.AreEqual(1, newPsms[1].FdrInfo.CumulativeTargetNotch);
Assert.AreEqual(0, newPsms[2].FdrInfo.CumulativeDecoyNotch);
Assert.AreEqual(1, newPsms[2].FdrInfo.CumulativeTargetNotch);
Assert.AreEqual(0, newPsms[0].FdrInfo.CumulativeDecoy);
Assert.AreEqual(1, newPsms[0].FdrInfo.CumulativeTarget);
Assert.AreEqual(0, newPsms[1].FdrInfo.CumulativeDecoy);
Assert.AreEqual(2, newPsms[1].FdrInfo.CumulativeTarget);
Assert.AreEqual(0, newPsms[2].FdrInfo.CumulativeDecoy);
Assert.AreEqual(3, newPsms[2].FdrInfo.CumulativeTarget);
}
public static void TestModificationAnalysisWithNonLocalizedPtms()
{
IScan scan = new ThisTestScan();
ModificationMotif.TryGetMotif("N", out ModificationMotif motif1);
Modification mod1 = new Modification(_originalId: "mod1", _modificationType: "mt", _target: motif1, _locationRestriction: "Anywhere.", _monoisotopicMass: 10, _neutralLosses: new Dictionary <DissociationType, List <double> > {
{ MassSpectrometry.DissociationType.AnyActivationType, new List <double> {
10
} }
});
IDictionary <int, List <Modification> > oneBasedModifications = new Dictionary <int, List <Modification> >
{
{ 2, new List <Modification> {
mod1
} },
{ 7, new List <Modification> {
mod1
} },
};
Protein protein1 = new Protein("MNLDLDNDL", "prot1", oneBasedModifications: oneBasedModifications);
Dictionary <int, Modification> allModsOneIsNterminus1 = new Dictionary <int, Modification>
{
{ 2, mod1 },
};
PeptideWithSetModifications pwsm1 = new PeptideWithSetModifications(protein1, new DigestionParams(), 2, 9, CleavageSpecificity.Unknown, null, 0, allModsOneIsNterminus1, 0);
Dictionary <int, Modification> allModsOneIsNterminus3 = new Dictionary <int, Modification>
{
{ 7, mod1 },
};
PeptideWithSetModifications pwsm2 = new PeptideWithSetModifications(protein1, new DigestionParams(), 2, 9, CleavageSpecificity.Unknown, null, 0, allModsOneIsNterminus3, 0);
CommonParameters CommonParameters = new CommonParameters(digestionParams: new DigestionParams(maxMissedCleavages: 0, minPeptideLength: 1), scoreCutoff: 1);
PeptideSpectralMatch myPsm = new PeptideSpectralMatch(pwsm1, 0, 10, 0, scan, new DigestionParams(), new List <MatchedFragmentIon>());
myPsm.AddOrReplace(pwsm2, 10, 0, true, new List <MatchedFragmentIon>(), 0);
myPsm.ResolveAllAmbiguities();
MassDiffAcceptor searchMode = new SinglePpmAroundZeroSearchMode(5);
List <Protein> proteinList = new List <Protein> {
protein1
};
FdrAnalysisEngine fdrAnalysisEngine = new FdrAnalysisEngine(new List <PeptideSpectralMatch> {
myPsm
}, searchMode.NumNotches, CommonParameters, new List <string>());
fdrAnalysisEngine.Run();
ModificationAnalysisEngine modificationAnalysisEngine = new ModificationAnalysisEngine(new List <PeptideSpectralMatch> {
myPsm
}, new CommonParameters(), new List <string>());
var res = (ModificationAnalysisResults)modificationAnalysisEngine.Run();
Assert.AreEqual(1, res.CountOfEachModSeenOnProteins.Count());
Assert.AreEqual(2, res.CountOfEachModSeenOnProteins[mod1.IdWithMotif]);
Assert.AreEqual(0, res.CountOfModsSeenAndLocalized.Count());
Assert.AreEqual(0, res.CountOfAmbiguousButLocalizedModsSeen.Count);
Assert.AreEqual(1, res.CountOfUnlocalizedMods[mod1.IdWithMotif]); // Saw it, but not sure where!
Assert.AreEqual(0, res.CountOfUnlocalizedFormulas.Count());
}
public static void TestAnalysisEngineTests()
{
CommonParameters CommonParameters = new CommonParameters
{
DigestionParams = new DigestionParams
{
Protease = new Protease("Custom Protease", new List <string> {
"K"
}, new List <string>(), TerminusType.C, CleavageSpecificity.Full, null, null, null),
MinPeptideLength = null,
MaxMissedCleavages = 0,
MaxModificationIsoforms = 1042,
},
ConserveMemory = false,
ScoreCutoff = 1,
ProductMassTolerance = new PpmTolerance(10),
};
List <ModificationWithMass> localizeableModifications = new List <ModificationWithMass>();
List <ModificationWithMass> variableModifications = new List <ModificationWithMass>();
List <ModificationWithMass> fixedModifications = new List <ModificationWithMass>();
Dictionary <ModificationWithMass, ushort> modsDictionary = new Dictionary <ModificationWithMass, ushort>();
foreach (var mod in fixedModifications)
{
modsDictionary.Add(mod, 0);
}
int i = 1;
foreach (var mod in variableModifications)
{
modsDictionary.Add(mod, (ushort)i);
i++;
}
foreach (var mod in localizeableModifications)
{
modsDictionary.Add(mod, (ushort)i);
i++;
}
var proteinList = new List <Protein> {
new Protein("MNNNKQQQ", "accession")
};
var modPep = proteinList.First().Digest(CommonParameters.DigestionParams, fixedModifications, variableModifications).Last();
HashSet <PeptideWithSetModifications> value1 = new HashSet <PeptideWithSetModifications> {
modPep
};
CompactPeptide compactPeptide1 = new CompactPeptide(value1.First(), TerminusType.None);
Assert.AreEqual("QQQ", value1.First().BaseSequence);
var modPep2 = proteinList.First().Digest(CommonParameters.DigestionParams, fixedModifications, variableModifications).First();
HashSet <PeptideWithSetModifications> value2 = new HashSet <PeptideWithSetModifications> {
modPep2
};
CompactPeptide compactPeptide2 = new CompactPeptide(value2.First(), TerminusType.None);
Assert.AreEqual("MNNNK", value2.First().BaseSequence);
var modPep3 = proteinList.First().Digest(CommonParameters.DigestionParams, fixedModifications, variableModifications).ToList()[1];
HashSet <PeptideWithSetModifications> value3 = new HashSet <PeptideWithSetModifications> {
modPep3
};
CompactPeptide compactPeptide3 = new CompactPeptide(value3.First(), TerminusType.None);
Assert.AreEqual("NNNK", value3.First().BaseSequence);
//newPsms[0] = new List<PsmParent>[] { new List<PsmParent>{ new PsmModern(compactPeptide1, null, 1, 1, 2, 2, 1,1, 1, 1, 3,0) },
// new List<PsmParent>{ new PsmModern(compactPeptide2, null, 2,2+132.040,3,3,2,2,2,2,2,0) },
// new List<PsmParent>{ new PsmModern(compactPeptide3, null, 3, 3, 4, 3, 3, 3, 3, 3, 3, 0)} };
Ms2ScanWithSpecificMass scanA = new Ms2ScanWithSpecificMass(new MzmlScanWithPrecursor(2, new MzmlMzSpectrum(new double[] { 1 }, new double[] { 1 }, false), 1, true, Polarity.Positive, double.NaN, null, null, MZAnalyzerType.Orbitrap, double.NaN, double.NaN, null, null, double.NaN, null, DissociationType.AnyActivationType, 1, null, null, "scan=1"), 1, 1, null);
Ms2ScanWithSpecificMass scanB = new Ms2ScanWithSpecificMass(new MzmlScanWithPrecursor(3, new MzmlMzSpectrum(new double[] { 1 }, new double[] { 1 }, false), 1, true, Polarity.Positive, double.NaN, null, null, MZAnalyzerType.Orbitrap, double.NaN, double.NaN, null, null, double.NaN, null, DissociationType.AnyActivationType, 1, null, null, "scan=2"), 2 + 132.040, 1, null);
Ms2ScanWithSpecificMass scanC = new Ms2ScanWithSpecificMass(new MzmlScanWithPrecursor(4, new MzmlMzSpectrum(new double[] { 1 }, new double[] { 1 }, false), 1, true, Polarity.Positive, double.NaN, null, null, MZAnalyzerType.Orbitrap, double.NaN, double.NaN, null, null, double.NaN, null, DissociationType.AnyActivationType, 1, null, null, "scan=3"), 3, 1, null);
PeptideSpectralMatch matchA = new PeptideSpectralMatch(compactPeptide1, 0, 0, 0, scanA);
PeptideSpectralMatch matchB = new PeptideSpectralMatch(compactPeptide2, 0, 0, 0, scanB);
PeptideSpectralMatch matchC = new PeptideSpectralMatch(compactPeptide3, 0, 0, 0, scanC);
var newPsms = new List <PeptideSpectralMatch> {
matchA, matchB, matchC
};
IMsDataFile <IMsDataScan <IMzSpectrum <IMzPeak> > > myMsDataFile = new TestDataFile(new List <PeptideWithSetModifications> {
value1.First(), value2.First(), value3.First()
});
var searchMode = new SinglePpmAroundZeroSearchMode(5);
Action <List <PeptideSpectralMatch>, string, List <string> > action2 = (List <PeptideSpectralMatch> l, string s, List <string> sdf) => {; };
bool DoPrecursorDeconvolution = true;
bool UseProvidedPrecursorInfo = true;
double DeconvolutionIntensityRatio = 4;
int DeconvolutionMaxAssumedChargeState = 10;
Tolerance DeconvolutionMassTolerance = new PpmTolerance(5);
var arrayOfMs2ScansSortedByMass = MetaMorpheusTask.GetMs2Scans(myMsDataFile, null, DoPrecursorDeconvolution, UseProvidedPrecursorInfo, DeconvolutionIntensityRatio, DeconvolutionMaxAssumedChargeState, DeconvolutionMassTolerance).OrderBy(b => b.PrecursorMass).ToArray();
Action <BinTreeStructure, string> action1 = (BinTreeStructure l, string s) =>
{
Assert.AreEqual(1, l.FinalBins.Count);
};
SequencesToActualProteinPeptidesEngine sequencesToActualProteinPeptidesEngine = new SequencesToActualProteinPeptidesEngine(newPsms, proteinList, fixedModifications, variableModifications, new List <ProductType> {
ProductType.B, ProductType.Y
}, new List <IDigestionParams> {
CommonParameters.DigestionParams
}, CommonParameters.ReportAllAmbiguity, new List <string>());
var res = (SequencesToActualProteinPeptidesEngineResults)sequencesToActualProteinPeptidesEngine.Run();
var compactPeptideToProteinPeptideMatching = res.CompactPeptideToProteinPeptideMatching;
foreach (var huh in newPsms)
{
if (huh != null)
{
huh.MatchToProteinLinkedPeptides(compactPeptideToProteinPeptideMatching);
}
}
FdrAnalysisEngine engine = new FdrAnalysisEngine(newPsms, searchMode.NumNotches, false, new List <string> {
"ff"
});
engine.Run();
}
public static void FdrTestMethod()
{
MassDiffAcceptor searchModes = new DotMassDiffAcceptor(null, new List <double> {
0, 1.0029
}, new PpmTolerance(5));
List <string> nestedIds = new List <string>();
Protein p = new Protein("MNKNNKNNNKNNNNK", null);
DigestionParams digestionParams = new DigestionParams();
var digested = p.Digest(digestionParams, new List <ModificationWithMass>(), new List <ModificationWithMass>()).ToList();
PeptideWithSetModifications pep1 = digested[0];
PeptideWithSetModifications pep2 = digested[1];
PeptideWithSetModifications pep3 = digested[2];
PeptideWithSetModifications pep4 = digested[3];
TestDataFile t = new TestDataFile(new List <PeptideWithSetModifications> {
pep1, pep2, pep3
});
CompactPeptide peptide1 = new CompactPeptide(pep1, TerminusType.None);
IMsDataScanWithPrecursor <IMzSpectrum <IMzPeak> > mzLibScan1 = t.GetOneBasedScan(2) as IMsDataScanWithPrecursor <IMzSpectrum <IMzPeak> >;
Ms2ScanWithSpecificMass scan1 = new Ms2ScanWithSpecificMass(mzLibScan1, peptide1.MonoisotopicMassIncludingFixedMods.ToMz(1), 1, null);
PeptideSpectralMatch psm1 = new PeptideSpectralMatch(peptide1, 0, 3, 0, scan1);
CompactPeptide peptide2 = new CompactPeptide(pep2, TerminusType.None);
IMsDataScanWithPrecursor <IMzSpectrum <IMzPeak> > mzLibScan2 = t.GetOneBasedScan(4) as IMsDataScanWithPrecursor <IMzSpectrum <IMzPeak> >;
Ms2ScanWithSpecificMass scan2 = new Ms2ScanWithSpecificMass(mzLibScan2, peptide2.MonoisotopicMassIncludingFixedMods.ToMz(1), 1, null);
PeptideSpectralMatch psm2 = new PeptideSpectralMatch(peptide2, 1, 2, 1, scan2);
CompactPeptide peptide3 = new CompactPeptide(pep3, TerminusType.None);
IMsDataScanWithPrecursor <IMzSpectrum <IMzPeak> > mzLibScan3 = t.GetOneBasedScan(6) as IMsDataScanWithPrecursor <IMzSpectrum <IMzPeak> >;
Ms2ScanWithSpecificMass scan3 = new Ms2ScanWithSpecificMass(mzLibScan3, peptide3.MonoisotopicMassIncludingFixedMods.ToMz(1), 1, null);
PeptideSpectralMatch psm3 = new PeptideSpectralMatch(peptide3, 0, 1, 2, scan3);
CompactPeptide peptide4 = new CompactPeptide(pep4, TerminusType.None);
psm3.AddOrReplace(peptide4, 1, 1, true);
Dictionary <CompactPeptideBase, HashSet <PeptideWithSetModifications> > matching = new Dictionary <CompactPeptideBase, HashSet <PeptideWithSetModifications> >
{
{
peptide1, new HashSet <PeptideWithSetModifications> {
pep1
}
},
{
peptide2, new HashSet <PeptideWithSetModifications> {
pep2
}
},
{
peptide3, new HashSet <PeptideWithSetModifications> {
pep3
}
},
{
peptide4, new HashSet <PeptideWithSetModifications> {
pep4
}
},
};
psm1.MatchToProteinLinkedPeptides(matching);
psm2.MatchToProteinLinkedPeptides(matching);
psm3.MatchToProteinLinkedPeptides(matching);
var newPsms = new List <PeptideSpectralMatch> {
psm1, psm2, psm3
};
FdrAnalysisEngine fdr = new FdrAnalysisEngine(newPsms, searchModes.NumNotches, true, nestedIds);
fdr.Run();
Assert.AreEqual(2, searchModes.NumNotches);
Assert.AreEqual(0, newPsms[0].FdrInfo.CumulativeDecoyNotch);
Assert.AreEqual(1, newPsms[0].FdrInfo.CumulativeTargetNotch);
Assert.AreEqual(0, newPsms[1].FdrInfo.CumulativeDecoyNotch);
Assert.AreEqual(1, newPsms[1].FdrInfo.CumulativeTargetNotch);
Assert.AreEqual(0, newPsms[2].FdrInfo.CumulativeDecoyNotch);
Assert.AreEqual(1, newPsms[2].FdrInfo.CumulativeTargetNotch);
Assert.AreEqual(0, newPsms[0].FdrInfo.CumulativeDecoy);
Assert.AreEqual(1, newPsms[0].FdrInfo.CumulativeTarget);
Assert.AreEqual(0, newPsms[1].FdrInfo.CumulativeDecoy);
Assert.AreEqual(2, newPsms[1].FdrInfo.CumulativeTarget);
Assert.AreEqual(0, newPsms[2].FdrInfo.CumulativeDecoy);
Assert.AreEqual(3, newPsms[2].FdrInfo.CumulativeTarget);
}
public static void TestModificationAnalysis()
{
IScan scan = new ThisTestScan();
ModificationMotif.TryGetMotif("N", out ModificationMotif motif1);
ModificationWithMass mod1 = new ModificationWithMass("mod1", null, motif1, TerminusLocalization.Any, 10);
ModificationMotif.TryGetMotif("L", out ModificationMotif motif2);
ModificationWithMass mod2 = new ModificationWithMass("mod2", null, motif2, TerminusLocalization.Any, 10);
IDictionary <int, List <Modification> > oneBasedModifications = new Dictionary <int, List <Modification> >
{
{ 2, new List <Modification> {
mod1
} },
{ 5, new List <Modification> {
mod2
} },
{ 7, new List <Modification> {
mod1
} },
};
Protein protein1 = new Protein("MNLDLDNDL", "prot1", oneBasedModifications: oneBasedModifications);
Dictionary <int, ModificationWithMass> allModsOneIsNterminus1 = new Dictionary <int, ModificationWithMass>
{
{ 2, mod1 },
};
PeptideWithSetModifications pwsm1 = new PeptideWithSetModifications(0, protein1, 2, 9, allModsOneIsNterminus1);
CompactPeptideBase pep1 = new CompactPeptide(pwsm1, TerminusType.None);
Dictionary <int, ModificationWithMass> allModsOneIsNterminus2 = new Dictionary <int, ModificationWithMass>
{
{ 2, mod1 },
{ 7, mod1 },
};
PeptideWithSetModifications pwsm2 = new PeptideWithSetModifications(0, protein1, 2, 9, allModsOneIsNterminus2);
CompactPeptideBase pep2 = new CompactPeptide(pwsm2, TerminusType.None);
Dictionary <int, ModificationWithMass> allModsOneIsNterminus3 = new Dictionary <int, ModificationWithMass>
{
{ 7, mod1 },
};
PeptideWithSetModifications pwsm3 = new PeptideWithSetModifications(0, protein1, 2, 9, allModsOneIsNterminus3);
CompactPeptideBase pep3 = new CompactPeptide(pwsm3, TerminusType.None);
Dictionary <int, ModificationWithMass> allModsOneIsNterminus4 = new Dictionary <int, ModificationWithMass>
{
{ 8, mod1 },
};
PeptideWithSetModifications pwsm4 = new PeptideWithSetModifications(0, protein1, 1, 9, allModsOneIsNterminus4);
CompactPeptideBase pep4 = new CompactPeptide(pwsm4, TerminusType.None);
CommonParameters CommonParameters = new CommonParameters(
digestionParams: new DigestionParams(
maxMissedCleavages: 0,
minPeptideLength: 1,
maxModificationIsoforms: int.MaxValue),
scoreCutoff: 1);
var newPsms = new List <PeptideSpectralMatch>
{
new PeptideSpectralMatch(pep1, 0, 10, 0, scan, CommonParameters.DigestionParams),
new PeptideSpectralMatch(pep1, 0, 10, 0, scan, CommonParameters.DigestionParams),
new PeptideSpectralMatch(pep2, 0, 10, 0, scan, CommonParameters.DigestionParams),
new PeptideSpectralMatch(pep3, 0, 10, 0, scan, CommonParameters.DigestionParams),
new PeptideSpectralMatch(pep4, 0, 10, 0, scan, CommonParameters.DigestionParams),
};
MassDiffAcceptor searchMode = new SinglePpmAroundZeroSearchMode(5);
List <Protein> proteinList = new List <Protein> {
protein1
};
SequencesToActualProteinPeptidesEngine sequencesToActualProteinPeptidesEngine = new SequencesToActualProteinPeptidesEngine
(newPsms, proteinList, new List <ModificationWithMass>(), new List <ModificationWithMass>(), new List <ProductType>
{
ProductType.B, ProductType.Y
}, new List <DigestionParams> {
CommonParameters.DigestionParams
}, CommonParameters.ReportAllAmbiguity, CommonParameters, new List <string>());
var nice = (SequencesToActualProteinPeptidesEngineResults)sequencesToActualProteinPeptidesEngine.Run();
foreach (var psm in newPsms)
{
psm.MatchToProteinLinkedPeptides(nice.CompactPeptideToProteinPeptideMatching);
}
FdrAnalysisEngine fdrAnalysisEngine = new FdrAnalysisEngine(newPsms, searchMode.NumNotches, CommonParameters, new List <string>());
fdrAnalysisEngine.Run();
ModificationAnalysisEngine modificationAnalysisEngine = new ModificationAnalysisEngine(newPsms, new CommonParameters(), new List <string>());
var res = (ModificationAnalysisResults)modificationAnalysisEngine.Run();
Assert.AreEqual(2, res.AllModsOnProteins.Count());
Assert.AreEqual(2, res.AllModsOnProteins[mod1.id]);
Assert.AreEqual(1, res.AllModsOnProteins[mod2.id]);
Assert.AreEqual(1, res.ModsSeenAndLocalized.Count());
Assert.AreEqual(2, res.ModsSeenAndLocalized[mod1.id]);
Assert.AreEqual(0, res.AmbiguousButLocalizedModsSeen.Count());
Assert.AreEqual(0, res.UnlocalizedMods.Count());
Assert.AreEqual(0, res.UnlocalizedFormulas.Count());
}
public static void TestModificationAnalysisWithNonLocalizedPtms()
{
IScan scan = new ThisTestScan();
ModificationMotif.TryGetMotif("N", out ModificationMotif motif1);
ModificationWithMass mod1 = new ModificationWithMass("mod1", "mt", motif1, TerminusLocalization.Any, 10, neutralLosses: new List <double> {
10
});
IDictionary <int, List <Modification> > oneBasedModifications = new Dictionary <int, List <Modification> >
{
{ 2, new List <Modification> {
mod1
} },
{ 7, new List <Modification> {
mod1
} },
};
Protein protein1 = new Protein("MNLDLDNDL", "prot1", oneBasedModifications: oneBasedModifications);
Dictionary <int, ModificationWithMass> allModsOneIsNterminus1 = new Dictionary <int, ModificationWithMass>
{
{ 2, mod1 },
};
PeptideWithSetModifications pwsm1 = new PeptideWithSetModifications(0, protein1, 2, 9, allModsOneIsNterminus1);
CompactPeptideBase pep1 = new CompactPeptide(pwsm1, TerminusType.None);
Dictionary <int, ModificationWithMass> allModsOneIsNterminus3 = new Dictionary <int, ModificationWithMass>
{
{ 7, mod1 },
};
PeptideWithSetModifications pwsm3 = new PeptideWithSetModifications(0, protein1, 2, 9, allModsOneIsNterminus3);
CompactPeptideBase pep3 = new CompactPeptide(pwsm3, TerminusType.None);
var newPsms = new List <PeptideSpectralMatch>
{
new PeptideSpectralMatch(pep1, 0, 10, 0, scan),
new PeptideSpectralMatch(pep3, 0, 10, 0, scan),
};
MassDiffAcceptor searchMode = new SinglePpmAroundZeroSearchMode(5);
List <Protein> proteinList = new List <Protein> {
protein1
};
CommonParameters CommonParameters = new CommonParameters
{
DigestionParams = new DigestionParams
{
MinPeptideLength = null,
MaxMissedCleavages = 0,
MaxModificationIsoforms = int.MaxValue
},
ConserveMemory = false,
ScoreCutoff = 1,
};
SequencesToActualProteinPeptidesEngine sequencesToActualProteinPeptidesEngine = new SequencesToActualProteinPeptidesEngine(newPsms, proteinList, new List <ModificationWithMass>(), new List <ModificationWithMass>(), new List <ProductType> {
ProductType.B, ProductType.Y
}, new List <IDigestionParams> {
CommonParameters.DigestionParams
}, CommonParameters.ReportAllAmbiguity, new List <string>());
var nice = (SequencesToActualProteinPeptidesEngineResults)sequencesToActualProteinPeptidesEngine.Run();
foreach (var psm in newPsms)
{
psm.MatchToProteinLinkedPeptides(nice.CompactPeptideToProteinPeptideMatching);
}
Assert.AreEqual(2, nice.CompactPeptideToProteinPeptideMatching[pep1].Count);
FdrAnalysisEngine fdrAnalysisEngine = new FdrAnalysisEngine(newPsms, searchMode.NumNotches, false, new List <string>());
fdrAnalysisEngine.Run();
ModificationAnalysisEngine modificationAnalysisEngine = new ModificationAnalysisEngine(newPsms, new List <string>());
var res = (ModificationAnalysisResults)modificationAnalysisEngine.Run();
Assert.AreEqual(1, res.AllModsOnProteins.Count());
Assert.AreEqual(2, res.AllModsOnProteins[mod1.id]);
Assert.AreEqual(0, res.ModsSeenAndLocalized.Count());
Assert.AreEqual(0, res.AmbiguousButLocalizedModsSeen.Count);
Assert.AreEqual(1, res.UnlocalizedMods[mod1.id]); // Saw it, but not sure where!
Assert.AreEqual(0, res.UnlocalizedFormulas.Count());
}
