WAP to find out the percentage of amino acid residues from a protein data bank file ( PDB ID: 1XYZ.pdb)

Sure, here is a Python script that can be used to calculate the percentage of amino acid residues from a PDB file. This script assumes that you have the Biopython library installed. If not, you can install it using pip:

pip install biopython

Here is the Python script:

from Bio.PDB import PDBParser

def count_amino_acids(pdb_file):
    parser = PDBParser()
    structure = parser.get_structure('protein', pdb_file)

    amino_acids = ['ALA', 'ARG', 'ASN', 'ASP', 'CYS', 'GLU', 'GLN', 'GLY', 'HIS', 'ILE', 'LEU', 'LYS', 'MET', 'PHE', 'PRO', 'SER', 'THR', 'TRP', 'TYR', 'VAL']
    aa_count = {aa: 0 for aa in amino_acids}

    for model in structure:
        for chain in model:
            for residue in chain:
                if residue.get_resname() in aa_count:
                    aa_count[residue.get_resname()] += 1

    total_residues = sum(aa_count.values())
    percentages = {aa: (count / total_residues) * 100 for aa, count in aa_count.items()}

    return percentages

pdb_file = '1XYZ.pdb'  # replace with your pdb file path
percentages = count_amino_acids(pdb_file)

for aa, percentage in percentages.items():
    print(f'{aa}: {percentage}%')

This script first defines a list of the 20 standard amino acids. It then creates a dictionary to keep track of the count of each amino acid in the PDB file. It uses the Biopython PDBParser to parse the PDB file and iterates over each residue in each chain in each model in the structure. If the residue is an amino acid (i.e., its name is in the list of