import mdtraj as md
import numpy as np
def protein_surface_area(pdb_file):
"""
1. This function loads a pdb
2. Calculates the Solvent Accessible Surface Area (sasa) per residue
3. Returns it
"""
[insert]
return sasa
import mdtraj as md
import numpy as np
def protein_surface_area(pdb_file):
"""
1. This function loads a pdb
2. Calculates the Solvent Accessible Surface Area (sasa) per residue
3. Returns it
"""
traj = md.load(pdb_file)
sasa = md.shrake_rupley(traj, mode='residue')
sasa = np.array(sasa)
return sasa
Go to Google Form