Skip to main content
zenodoopen

rbLEC - restricted backbone Local Euler Characteristic - from CATH database

<p>-----------------------------------------------------------------------------------------------------------------------------------</p> <p><strong>Author: Rodrigo A. Moreira (C) 2023<br> https://orcid.org/0000-0002-7605-8722<br> LICENSE: CC BY-NC-ND 4.0 (https://creativecommons.org/licenses/by-nc-nd/4.0/)</strong></p> <p>----------------------------------------------------------------------------------------------------------------------------------</p> <p><strong>rbLEC - Local Euler Charactersitics - from CATH database</strong></p> <p>----------------------------------------------------------------------------------------------------------------------------------</p> <p><strong>A. rbLEC NETWORK</strong></p> <p>&nbsp;&nbsp;&nbsp; [I] The networks for each PDB[1] structure is defined by the PDB atoms N,CA,C of each residue as nodes of a graph G.<br> &nbsp;&nbsp; &nbsp;[II] An edge of G is set if the distance between two atom in [I] is greater than 2.0 Angstrons.<br> &nbsp;&nbsp; &nbsp;[III] The graph G is defined in the files with extensions &quot;.network_backboneRE_heavy_gt2&quot;</p> <p>Equation (1) [6,7]<br> &nbsp;&nbsp; &nbsp;\begin{equation}<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;\chi = \sum_{k=1}^{N} \kappa_k = \sum_{k=1}^{N} \underbrace{ \left(1 + \sum_{l=1}^{\infty} (-1)^{l} \frac{v_{l-1}}{l+1} \right)_{k}}_{\kappa_k}<br> &nbsp;&nbsp; &nbsp;\end{equation}</p> <p>Equation (2)<br> &nbsp;&nbsp; &nbsp;\begin{equation}<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;LEC = \sum_{m \in R} \kappa_m = \kappa_{N} + \kappa_{CA} + \kappa_{C}<br> &nbsp;&nbsp; &nbsp;\end{equation}</p> <p><strong>B. FILENAME EXTENSIONS</strong></p> <p>&nbsp;<strong> B.1</strong> Basic files</p> <p>&quot;.fixed&quot;<br> &nbsp;&nbsp; &nbsp;PDB file after use of pdbfixer[2] in structures from CATH database.</p> <p>&quot;.dssp&quot;<br> &nbsp;&nbsp; &nbsp;Output of DSSP[3] software</p> <p>&quot;.stride&quot;<br> &nbsp;&nbsp; &nbsp;Output of STRIDE[4] software</p> <p>&nbsp; <strong>B.2</strong> Data files</p> <p>&quot;.network_backboneRE_heavy_gt2&quot; - Generate by <strong>D.2</strong> below.<br> &nbsp;&nbsp;&nbsp; Describe the network graph, as described in <strong>A.</strong> above.</p> <p>&quot;.knill_curvature&quot; - Generate by <strong>D.1</strong> below.<br> &nbsp;&nbsp; &nbsp;Contain the filtration of kappas for each vertice of the network.</p> <p>&quot;.residues_curvature&quot;&nbsp; - Generate by <strong>D.1</strong> below.<br> &nbsp;&nbsp; &nbsp;They are the filtration of LEC, Equation (2) above, for each residue, namely summation of 3 kappas from respective &#39;.knill_curvature&#39;, correspoings to PDB atoms N,CA and C, describe in <strong>A.</strong> above.</p> <p>&quot;.label&quot;&nbsp; - Generated by <strong>D.3 </strong>below<br> &nbsp;&nbsp; &nbsp;Extra file for easier assesment of structures. They have the same information about LEC as described in respective &quot;.residue_curvature&quot; file extensions, but merge also the information from &quot;.dssp&quot; and &quot;.stride&quot; classes as well as residue name and residue ID for each molecule.<br> &nbsp;&nbsp; &nbsp;Format of columns:<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;cutoff resname resid DSSP_class STRIDE_class LEC</p> <p><strong>C. FOLDERS</strong></p> <p>&nbsp;&nbsp;&nbsp; CATH_FIXED (after uncompress cath_fixed.tar.xz, approximately 13GB)<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;contains the fixed PDBs and LECs from CATH[5] database</p> <p><strong>D. SOFTWARE</strong><br> &nbsp;&nbsp;&nbsp; <strong>D.1</strong> lec.py:&nbsp; compute the kappas in Equation (1) above.<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;Example usage:<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;$ python3 lec.py CATH_FIXED/2x0qA02/2x0qA02<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;It will create the files with extension &quot;.kappas&quot; and &quot;.relec&quot;, which reproduces the respectively the files with extension &quot;<strong>.knill_curvature</strong>&quot; and &quot;<strong>.residue_curvature</strong>&quot;.</p> <p>&nbsp;&nbsp;&nbsp;&nbsp;<strong> D.2</strong> pdb2network.lua: creates rbLEC network file (number of nodes and edges list) from PDB to be used as input by lec.py.<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;Example usage:<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;$ lua pdb2rbLEC.lua CATH_FIXED/2x0qA02/2x0qA02.fixed<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; Output reproduces the file CATH_FIXED/2x0qA02/2x0qA02.pdb.<strong>network_backboneRE_heavy_gt2</strong></p> <p>&nbsp;&nbsp;&nbsp;&nbsp; <strong>D.3</strong> label.lua: create files with extension &#39;*.label&#39; from files &#39;*.pdb.stride&#39;, &#39;*.pdb.dssp&#39; and &#39;*.pdb.network_backboneRE_heavy_gt2.residues_curvature.<br> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Example usage:<br> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; $ lua label.lua CATH_FIXED/2x0qA02/2x0qA02.pdb<br> &nbsp; &nbsp; &nbsp; &nbsp; Output reproduces the file CATH_FIXED/2x0qA02/2x0qA02.pdb.<strong>network_backboneRE_heavy_gt2.residues_curvature.label</strong></p> <p><strong>REFERENCES</strong><br> [1] Herman, H., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T., Weissig, H., Shindyalov, I., &amp; Bourne, P. (2000). The protein data bank. Nucleic acids research, 28, 235&ndash;42.<br> [2] Eastman, P., Swails, J., Chodera, J., McGibbon, R., Zhao, Y., Beauchamp, K., Wang, L.P., Simmonett, A., Harrigan, M., Stern, C., &amp; others (2017). OpenMM 7: Rapid development of high performance algorithms for molecular dynamics. PLoS computational biology, 13(7), e1005659.<br> [3] Kabsch, W., &amp; Sander, C. (1983). Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers: Original Research on Biomolecules, 22(12), 2577&ndash;2637.<br> [4] Frishman, D., &amp; Argos, P. (1995). Knowledge-based protein secondary structure assignment. Proteins: Structure, Function, and Bioinformatics, 23(4), 566&ndash;579.<br> [5] Knudsen, M., &amp; Wiuf, C. (2010). The CATH database. Human genomics, 4(3), 1&ndash;6.<br> [6] Levitt, N. (1992). The Euler characteristic is the unique locally determined numerical homotopy invariant of finite complexes. Discrete &amp; computational geometry, 7, 59&ndash;67.<br> [7] Knill, O. (2011). A graph theoretical Gauss-Bonnet-Chern theorem. arXiv preprint arXiv:1111.5395.</p> <p>&nbsp;</p> <p>&nbsp;</p>

ShareScore

40/100

Overall dataset sharing score

Score breakdown

These five areas show where the dataset supports — or may limit — practical reuse.

Stewardship
8
Harmonization
8
Access
16
Reuse readiness
8
Engagement
0

Topics