How to extract the atomic coordinates of a specific amino acid from a pdb file
1 view (last 30 days)
Show older comments
Hello, Is it possible in Matlab bioinformatics tool, to extract the coordinates of a specific amino acid like lysine and store it in a Matrix.. There will be many lysine residues in a protein. How to store the atomic coordinates of each and every lysine residue from a protein in matrices Lysine1_coordinates, Lysine 2_coordinates and so on from a pdb file using Matlab?
0 Comments
Accepted Answer
Arthur Goldsipe
on 16 Feb 2023
What is your definition for the location of lysine amino acid? The reason I ask is that a PDB file contain information about the coordinates of specific atoms, not amino acids. As I'm sure you know from the other questions you've already posted, you can use getpdb to access this information in MATLAB. This information includes whether specific atoms are part of a lysine (from the resName field) and which amino acide in the seuqence (from the resSeq field), as well as the coordinates in the X, Y, and Z fields.
My assessment is that this question is a bit too general right now for you to get much help. I suggest you to attempt to sketch out what you want to do in more detail in a combination of MATLAB code and pseudocode. Then, perhaps we can offer more specific advice on whatever parts you're stuck on.
3 Comments
Arthur Goldsipe
on 22 Feb 2023
The simplest way I can think of to fix your code is to do the following:
Lys = {};
for i = 1:5
Lys{i} = LysXYZ((22*i-21):(22*i),:);
end
Note that you have to write Lys{i} instead of Lys(i). That's because I'm creating Lys as a cell array. You can read more about them here.
More Answers (0)
See Also
Categories
Find more on Nucleotide Sequence Analysis in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!