Skip to content

New plugin: ProtVar#818

Merged
jamie-m-a merged 13 commits intoEnsembl:mainfrom
nakib103:protvar
Feb 4, 2026
Merged

New plugin: ProtVar#818
jamie-m-a merged 13 commits intoEnsembl:mainfrom
nakib103:protvar

Conversation

@nakib103
Copy link
Copy Markdown
Contributor

@nakib103 nakib103 commented Dec 17, 2025

ENSVAR-6378

Test

Example variant

7 38269535 rs1482920289 T C
22 39964438 rs934530393 A G
17 42845309 rs2055585941 C G

Check the JIRA ticket for example data file for these variant.

@nakib103 nakib103 marked this pull request as draft December 17, 2025 16:12
@nakib103 nakib103 marked this pull request as ready for review January 5, 2026 09:43
@nakib103 nakib103 requested a review from jamie-m-a January 5, 2026 09:44
Comment thread ProtVar.pm
$pos-- if $translation_seq =~ /^X/;

# we need position of peptide in the ALL_AAS array
my $peptide_number = ( first_index { $_ eq $peptide } @ALL_AAS );
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is going to return -1 if the AA is X or * or an ambiguous code. Maybe they don't exist in the data? Maybe safer to handle anyway as it'll probably break the next function that uses $peptide_number

Comment thread ProtVar.pm Outdated
=head1 DESCRIPTION

An Ensembl VEP plugin that retrieves data from ProtVar resource providing contexualised information for
missense variation such destabilization of protein structure, overlapping protein pocket,
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
missense variation such destabilization of protein structure, overlapping protein pocket,
missense variation such as destabilization of protein structure, overlapping protein pocket,

Comment thread ProtVar.pm Outdated
$header{ProtVar_pocket} .= "MpLDDT - Mean pLDDT score of all the residues from AlphaFold2 model used to form the pocket, ";
$header{ProtVar_pocket} .= "energy - Pocket energy per volume (in kcal/mol unit), ";
$header{ProtVar_pocket} .= "burriedness - Pocket burriedness (0.0 means entirely exposed and 1.0 means entierly burried), ";
$header{ProtVar_pocket} .= "RoG - Pocket compactness in terms of radius of gyration (in Angstrom unit), ";
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
$header{ProtVar_pocket} .= "RoG - Pocket compactness in terms of radius of gyration (in Angstrom unit), ";
$header{ProtVar_pocket} .= "RoG - Pocket compactness in terms of radius of gyration (in Angstroms), ";

Comment thread ProtVar.pm Outdated
my %header;

if ( defined $self->{stability} ) {
$header{ProtVar_stability} = "Impact on protein stability in ddG enegry difference.";
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
$header{ProtVar_stability} = "Impact on protein stability in ddG enegry difference.";
$header{ProtVar_stability} = "Impact on protein stability in ddG energy difference (kcal/mol). Values over 2 are considered likely to be destabilising";

Copy link
Copy Markdown
Contributor

@jamie-m-a jamie-m-a left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

LGTM

@jamie-m-a jamie-m-a merged commit 676a40c into Ensembl:main Feb 4, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants