text: import rdkit
  from rdkit import Chem
  
  mol = rdkit.Chem.rdmolfiles.MolFromPDBFile('bene.pdb')
  SMILES_pdb = Chem.MolToSmiles(mol)
  
  from rdkit import Chem
  from rdkit.Chem.Draw import IPythonConsole
  from rdkit.Chem import AllChem
  from rdkit.Chem import rdFreeSASA
  
  mol1 = Chem.MolFromSmiles(SMILES_pdb)
  mol2 = Chem.MolFromSmiles('Oc1ccncc1')
  
  #add H atoms
  hmol1 = Chem.AddHs(mol1)
  hmol2 = Chem.AddHs(mol2)
  
  AllChem.EmbedMolecule(hmol1)
  AllChem.EmbedMolecule(hmol2)
  
  radii1 = rdFreeSASA.classifyAtoms(hmol1)
  radii2 = rdFreeSASA.classifyAtoms(hmol2)
  #Here we can calculate the sasa of the micromolecule
  rdFreeSASA.CalcSASA(hmol1, radii1)
  rdFreeSASA.CalcSASA(hmol2, radii2)
  
  atoms1 = hmol1.GetAtoms()
  atoms2 = hmol2.GetAtoms()
  
  for i in range(len(atoms1)):
      print(atoms1[i].GetSymbol(), atoms1[i].GetProp('SASAClassName'), atoms1[i].GetProp("SASA"))
      
  sum(float(a.GetProp("SASA")) for a in atoms1)
  
  for i in range(len(atoms2)):
      print(atoms2[i].GetSymbol(), atoms2[i].GetProp('SASAClassName'), atoms2[i].GetProp("SASA"))
       
  sum(float(a.GetProp("SASA")) for a in atoms2)