Skip to main content

This package has functions for the conversion of amino acid sequences to physicochemical vectors and the subsequent analysis of those vector sequences.

Project description

PCDTW is a package that implements the conversion of amino acid sequences to physicochemical vectors and subsequently allows for alignment of the sequences based on those vectors, development of consensus vectors that can be used to search databases for similar physicochemical profiles, development of the DTW distance between two physicochemical vectors and a few other functions. The basis for this package can be found in three publications and should be consulted for further background [1–3].

To install PCDTW (Two Options): -Use ‘pip install PCDTW’ in a powershell prompt -Use ‘!pip install PCDTW’ in a jupyter notebook

To use PCDTW: Use ‘import PCDTW’

Citations

1)Dixson, J.D.; Vumma, L.; Azad, R.K. An Analysis of Combined Molecular Weight and Hydrophobicity Similarity between the Amino Acid Sequences of Spike Protein Receptor Binding Domains of Betacoronaviruses and Functionally Similar Sequences from Other Virus Families. Microorganisms 2024, 12.

2)Dixson, J.D.; Azad, R.K. Physicochemical Evaluation of Remote Homology in the Twilight Zone. Proteins Struct. Funct. Bioinforma. 2024, n/a, doi:https://doi.org/10.1002/prot.26742.

3)Dixson, J.D.; Azad, R.K. A Novel Predictor of ACE2-Binding Ability among Betacoronaviruses. Evol. Med. Public Heal. 2021, 9, 360–373, doi:10.1093/EMPH/EOAB032.

Usage:

  1. To convert an amino acid sequence to vector form using two physicochemical properties:

    PCDTW.PCDTWConvert(x, PCProp1='Mass', PCProp2='HydroPho', normalize=False)
    

    PCProp1/PCProp2 options:

    • 'HydroPho'
    • 'HydroPhIl'
    • 'Hbond'
    • 'SideVol'
    • 'Polarity'
    • 'Polarizability'
    • 'SASA'
    • 'NCI'
    • 'Mass'

    Normalization: If normalize is set to True then the individual physicochemical scalar values for each amino acid are absolute maximum normalized before converting the amino acid sequence to vector form.

  2. To align two amino acid sequences using DTW and two physicochemical properties:

    PCDTW.PCDTWAlign(inputseq1str, inputseq2str, PCProp1='Mass', PCProp2='HydroPho', Penalty=0, Window=3)
    
    • window = size of Sakoe-Chiba band
    • penalty = somewhat equivalent to mismatch penalty in standard dynamic programming based alignment

    Returns a dictionary containing the following values:

    • 'Seq1AlignedString'
    • 'Seq2AlignedString'
    • 'FullAlignment'
    • 'Identity'
    • 'ConsensusVector'

    Example to get the full alignment and identity:

    seq1 = "MSDSNQGNNQQNYQQYSQNGNQQQGNNRYQG"
    seq2 = "MMNNNGNQVSNLSNALRQVNIGNRNSNTTT"
    print(PCDTWAlign(seq1, seq2)['FullAlignment'])
    print(PCDTWAlign(seq1, seq2)['Identity'])
    
  3. To get the PCDTW distance between two sequences normalized to the number of amino acids in the alignment:

    Dist=PCDTW.PCDTWDist(Seq1, Seq2)
    print(Dist)
    

    Example to get the distance:

    seq1 = "MSDSNQGNNQQNYQQYSQNGNQQQGNNRYQG"
    seq2 = "MMNNNGNQVSNLSNALRQVNIGNRNSNTTT"
    print(PCDTWDist(seq1, seq2))
    
  4. To get synthetically evolved homologs for an input sequence:

    SynHomologs=PCDTW.PCEvolve(Seq='GALM', PCProp1='Mass', PCProp2='HydroPho', BaseName='ProtX')
    print(SynHomologs)
    

    PCProp1/PCProp2 options:

    • 'HydroPho'
    • 'HydroPhIl'
    • 'Hbond'
    • 'SideVol'
    • 'Polarity'
    • 'Polarizability'
    • 'SASA'
    • 'NCI'
    • 'Mass'
  5. To get a newick format tree using PCDTW that represents the physicochemical similarity of protein sequences:

    Newick=PCDTW.PCDTWTree(FastaFile='Your_File_Location.fasta',PCProp1='Mass', PCProp2='HydroPho')
    print(Newick)
    

    PCProp1/PCProp2 options:

    • 'HydroPho'
    • 'HydroPhIl'
    • 'Hbond'
    • 'SideVol'
    • 'Polarity'
    • 'Polarizability'
    • 'SASA'
    • 'NCI'
    • 'Mass'

    This function is derived from the original algorithm used in Dixson and Azad, 2021. Unlike the original algorithm the two physicochemical properties used can be set to any two from the nine included above. If the PCProps are not specified by the user then they default to mass and hydrophobicity. The hydrophobicity values used in this package vary slightly from those used in the original algorithm.

Dependency Citations:

Bio.Phylo:

Talevich, E., Invergo, B.M., Cock, P.J.A., & Chapman, B.A. (2012). Bio.Phylo: A unified toolkit for processing, analyzing, and visualizing phylogenetic trees in Biopython. BMC Bioinformatics, 13, 209

Biopython:

Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics, 25(11), 1422–1423. https://doi.org/10.1093/bioinformatics/btp163

dtaidistance:

Wannes Meert, Kilian Hendrickx, Toon Van Craenendonck, Pieter Robberechts, Hendrik Blockeel, & Jesse Davis. (2022). DTAIDistance (Version v2). Zenodo. http://doi.org/10.5281/zenodo.5901139

Matplotlib:

Hunter, J. D. (2007). Matplotlib: A 2D graphics environment. Computing in Science & Engineering, 9(3), 90–95.

numpy:

Harris, C.R., Millman, K.J., van der Walt, S.J. et al. (2020). Array programming with NumPy. Nature 585, 357–362. DOI: 10.1038/s41586-020-2649-2.

pandas:

McKinney, W. (2010). Data Structures for Statistical Computing in Python. Proceedings of the 9th Python in Science Conference (SciPy 2010).

SciPy:

Pauli Virtanen, Ralf Gommers, Travis E. Oliphant, Matt Haberland, Tyler Reddy, David Cournapeau, Evgeni Burovski, Pearu Peterson, Warren Weckesser, Jonathan Bright, Stéfan J. van der Walt, Matthew Brett, Joshua Wilson, K. Jarrod Millman, Nikolay Mayorov, Andrew R. J. Nelson, Eric Jones, Robert Kern, Eric Larson, CJ Carey, İlhan Polat, Yu Feng, Eric W. Moore, Jake VanderPlas, Denis Laxalde, Josef Perktold, Robert Cimrman, Ian Henriksen, E.A. Quintero, Charles R Harris, Anne M. Archibald, Antônio H. Ribeiro, Fabian Pedregosa, Paul van Mulbregt, and SciPy 1.0 Contributors. (2020) SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nature Methods, 17(3), 261-272.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

pcdtw-0.3.1.tar.gz (9.3 kB view details)

Uploaded Source

Built Distribution

pcdtw-0.3.1-py3-none-any.whl (10.3 kB view details)

Uploaded Python 3

File details

Details for the file pcdtw-0.3.1.tar.gz.

File metadata

  • Download URL: pcdtw-0.3.1.tar.gz
  • Upload date:
  • Size: 9.3 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.1.0 CPython/3.12.7

File hashes

Hashes for pcdtw-0.3.1.tar.gz
Algorithm Hash digest
SHA256 97ee57ddea48e1db90e48fdfb589fb1b26758ca8522ad15c73ff2931c0c00a4f
MD5 a6b1353a4b28b3d96376b93a64d4bcfa
BLAKE2b-256 57fe5dafcd43c2a3fbd04f93fbc5d5393732a2683a44ac597de33a1bff15f5ef

See more details on using hashes here.

File details

Details for the file pcdtw-0.3.1-py3-none-any.whl.

File metadata

  • Download URL: pcdtw-0.3.1-py3-none-any.whl
  • Upload date:
  • Size: 10.3 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/6.1.0 CPython/3.12.7

File hashes

Hashes for pcdtw-0.3.1-py3-none-any.whl
Algorithm Hash digest
SHA256 8d0deab1de2f7120a479dfac8a45edb0dd2f7cdd14c74fe3b7a8bbf0480ea92f
MD5 3100274094dc713627d37874c3928077
BLAKE2b-256 a2c7065c4b6fffabadd22ba330875084d5d33ceed2e3d7ff01a083a39f03baa4

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page