Euclidean distance between atoms using RDKit

Question

I'm trying to find the Euclidean distance between two atoms in the molecule with SMILES representation O=CC1OC12CC1OC12 using the rdkit package. Looking online, I have converged to the following code.

import numpy as np
from rdkit import Chem

mol = Chem.MolFromSmiles('O=CC1OC12CC1OC12')

conf = mol.GetConformer()
at1Coords = np.array(conf.GetAtomPosition(bond_i.GetBeginAtomIdx()))
at2Coords = np.array(conf.GetAtomPosition(bond_i.GetEndAtomIdx()))
print(np.linalg.norm(at2Coords - at1Coords))

But I get the error

conf = mol.GetConformer()  
ValueError: Bad Conformer Id

I was wondering how to fix the issue. Thanks!

Answer 1

The issue was that I had to embed the molecule to have a 3D conformation. So I added following line:

AllChem.EmbedMolecule(mol)

Feel free to answer with better suggestions and alternatives.