def CalculateVSAEstate(mol, bins=None):
"""
#################################################################
MOE-type descriptors using Estate indices and surface
area contributions.
vsaBins=[4.78,5.00,5.410,5.740,6.00,6.07,6.45,7.00,11.0]
You can specify your own bins to compute some descriptors
Usage:
result=CalculateVSAEstate(mol)
Input: mol is a molecule object
Output: result is a dict form
#################################################################
"""
temp = EVSA.VSA_EState_(mol, bins, force=1)
res = {}
for i, j in enumerate(temp):
res['VSAEstate' + str(i)] = round(j, 3)
return res
def CalculateVSAEstate(mol, bins=None):
"""
MOE-type descriptors using Estate indices and surface
area contributions.
"""
temp = EVSA.VSA_EState_(mol, bins, force=1)
res = {}
for i, j in enumerate(temp):
res['VSAEstate' + str(i)] = round(j, 3)
return res
def CalculateVSAEstate(mol: Chem.Mol, bins: List[float] = None) -> dict:
"""Get MOE-type descriptors using SA and Estate indices contributions.
:param bins: interval boundaries used in the P_VSA calculation.
The default VSA bins are [4.78,5.00,5.410,5.740,6.00,6.07,6.45,7.00,11.0].
"""
temp = EVSA.VSA_EState_(mol, bins, force=1)
res = {}
for i, j in enumerate(temp):
res[f'VSAEstate{i}'] = round(j, 3)
return res
def compute_MOE_descriptors(self):
"""compute the MOE-type descriptors.
Ref:???
Returns:
MOE_dict: MOE dictionary, data type: float
"""
assert type(self.Molecule) == Chem.rdchem.Mol
MOE_dict = {}
SlogP_VSA_names = []
for i in range(1, 13):
SlogP_VSA_names.append('SlogP_VSA' + str(i))
MOE_dict.update(
dict(zip(SlogP_VSA_names, rdDesc.SlogP_VSA_(self.Molecule))))
SMR_VSA_names = []
for i in range(1, 11):
SMR_VSA_names.append('SMR_VSA' + str(i))
MOE_dict.update(
dict(zip(SMR_VSA_names, rdDesc.SMR_VSA_(self.Molecule))))
PEOE_VSA_names = []
for i in range(1, 15):
PEOE_VSA_names.append('PEOE_VSA' + str(i))
MOE_dict.update(
dict(zip(PEOE_VSA_names, rdDesc.PEOE_VSA_(self.Molecule))))
EState_VSA_names = []
for i in range(1, 12):
EState_VSA_names.append('EState_VSA' + str(i))
MOE_dict.update(
dict(zip(EState_VSA_names, EState.EState_VSA_(self.Molecule))))
VSA_EState_names = []
for i in range(1, 12):
VSA_EState_names.append('VSA_EState' + str(i))
MOE_dict.update(
dict(zip(VSA_EState_names, EState.VSA_EState_(self.Molecule))))
return MOE_dict
def CalculateVSAEstate(mol, bins=None):
temp = EVSA.VSA_EState_(mol, bins, force=1)
res = {}
for i, j in enumerate(temp):
res['VSAEstate' + str(i)] = round(j, 3)
return res