mittinatten / freesasa

C-library for calculating Solvent Accessible Surface Areas
http://freesasa.github.io/
MIT License
103 stars 37 forks source link

Wrong Number of Atoms were calculated #74

Closed maxsw2 closed 2 years ago

maxsw2 commented 2 years ago

Hi, So I was using FreeSASA to calculate the SASA of my polymer system that has voids. So on PDB it says, I have 32998 atoms. However, on FreeSASA output, I get a much lower total of atoms. I have attached my PDB file below, and here is my output for FreeSASA.

PARAMETERS algorithm : Shrake & Rupley probe-radius : 1.400 threads : 2 testpoints : 100

INPUT source : traj.pdb chains : A model : 0 atoms : 13710

traj.txt

RESULTS (A^2) Total : 84985.85 Apolar : 0.00 Polar : 0.00 Unknown : 84985.85 CHAIN A : [84985.85]

mittinatten commented 2 years ago

Hi, FreeSASA ignores hydrogens by default, could that be an explanation? You can include hydrogens, but the atomic radii are calibrated for the hydrogens being excluded, so then you might want to use another set of radii.

Simon

On Fri, 18 Mar 2022 at 17:58, maxsw2 @.***> wrote:

Hi, So I was using FreeSASA to calculate the SASA of my polymer system that has voids. So on PDB it says, I have 32998 atoms. However, on FreeSASA output, I get a much lower total of atoms. I have attached my PDB file below, and here is my output for FreeSASA.

PARAMETERS algorithm : Shrake & Rupley probe-radius : 1.400 threads : 2 testpoints : 100

INPUT source : traj.pdb chains : A model : 0 atoms : 13710

traj.txt https://github.com/mittinatten/freesasa/files/8306381/traj.txt

RESULTS (A^2) Total : 84985.85 Apolar : 0.00 Polar : 0.00 Unknown : 84985.85 CHAIN A : [84985.85]

— Reply to this email directly, view it on GitHub https://github.com/mittinatten/freesasa/issues/74, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAVV4HEDX5BWG4B37CIHFKDVASY2NANCNFSM5RCMHRGA . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you are subscribed to this thread.Message ID: @.***>

mittinatten commented 2 years ago

Closing due to inactivity, feel free to reopen if you have more questions.