You should use quantum chemistry codes such as Gaussian to perform single point energy task based on the structure in .pdb file, the .wfn file can be outputted if you use proper keywords (for example, use "out=wfn" in Gaussian input file. See beginning of Chapter 4 of Multiwfn manual for detail)
How do you converter a pdb file or xsd file to wfn file?