4
$\begingroup$

I am interested in assigning 3D coordinates to (atoms in) some 10K molecules that I have, currently represented by SMILES. This is because as have been shown by many chemoinformatics papers, 3D and 2D structures together can better represent a molecule. This is also emphasized in large-scale molecular representation learning challenges, where 3D structures are estimated from DFT.

However, it seems time-costly to compute DFT with publicly available libraries (e.g. pyscf), or maybe I am missing something. And the faster approach used by RDKit or OpenBabel (force field-based, e.g. MMFF, UFF) generates less accurate 3D coordinates (yes, I did notice this post). I have been thinking of some ways: (1) find some possible conversions of my compounds to nominal identifiers (e.g. DrugBank ID) and extract 3D structures from those databases; (2) just run pyscf; (3) find some papers that predicts 3D structure from SMILES (does this kind of papers even exist?)

I'm not from a chemistry background but just started to look into interesting problems in computational chemistry. Would appreciate you experts' suggestions!

$\endgroup$
8
  • 3
    $\begingroup$ Welcome. For part 1, while I don't know DrugBank in particular, most of these databases use RDKit, Open Babel, or OMEGA (a commercial product) to generate coordinates using force fields. One challenge, if you're not from a chemistry background, is to know that for most molecules, there are many energetically-accessible conformers (e.g., you need multiple 3D geometries). $\endgroup$ Apr 28 at 16:51
  • 1
    $\begingroup$ You don't have to go full dft, you can try semiempiric models. But you have to start somewhere for that, too. Also the larger the molecules get, the more important Geoff's comment becomes. Here's a practical recommendation though: xtb. $\endgroup$ Apr 28 at 19:51
  • 1
    $\begingroup$ Even though I'm on board with this question (is on topic here), you'll might find a better audience at Matter Modeling. If you want, we can push it there; please do not cross post the question. $\endgroup$ Apr 28 at 19:55
  • $\begingroup$ What physical/chemical/biological property of the molecules / of the material is the property of interest the model should learn? Or are there multiple to account for? Perhaps there already are data compilations available (see an earlier post to narrow the question. Is the interest in the conformations of global energetic minima (per molecule), or equally for more shallow minima (which, depending on the property in mind, can be more interesting / more relevant)? $\endgroup$
    – Buttonwood
    May 3 at 14:18
  • $\begingroup$ @GeoffHutchison Thank you for the suggestions! I have sent an inquiry to DrugBank since they didn't provide any information re. the source of 3D structures. I am aware of conformers, and I'm happy to inject those equivariances into my encoder, so that as long as there is one valid (not sure how to define this chemically) 3D structure it should be fine. $\endgroup$
    – jasperhyp
    May 4 at 18:44

1 Answer 1

2
$\begingroup$

I'd say that the way to approach this problem greatly depends on what you want to do and why you need the 3D structure. What grade of precision do you need and what do you need to measure?

Molecules move all the time unless crystallized. Do you need the crystal structure or do you need a representative conformation that is often assumed in a solvent at room temperature? If you just need to represent the molecules in 3D in a conformation that "makes sense" for small molecules, a force field based method will generally work and produce conformations that work well enough; otherwise you can try some semiempirical methods, those are quite fast and cheap as well. DFT will allow you to predict with better precision the conformation such a molecule would assume in vacuum at 0K (unless you go through several passages to try and model something else). A third option could be using experimentally measured crystal structures, though I don't know how comfortable it is to go finding and cleaning up 10,000 of them. I would not use databases containing computed structures such as DrugBank as you'd have no control over the process nor the possibility to check the quality of the results.

You may obtain quite different conformations using the different methodologies. Your model may work on some molecules because they have a particular conformation and fail on others that would normally have a similar one, however also have another different one that comes out slightly favorable in the calculations. 3D QSAR is particularly complex because of this and due to the difficulties in encoding molecular information in a comparable way. In a three dimensional object it's much more difficult to look for similarities and dissimilarities.

Regarding DFT computation time, that greatly depends on the sampling method you use to search the conformational space. Performing a steepest descent/BFGS may give you results in a couple minutes when we're talking about small molecules. However that won't work when searching for a global minimum, I'd say that computing 10,000 structures using DFT and with no access to a large cluster where to parallelize the computations would take too much time.

Since you're new to chemistry; I'd advise you to stay away from 3D QSAR for the time being. It does provide its advantages, but carries along a large quantity of complexities which can just be skimmed over when working with 2D structures.

EDIT: I just noticed in the website you link they use models to predict the structure of a molecule as modeled by DFT. In that case you necessarily have to use DFT to obtain the structure. It appears as if they were suggesting to use the 3D structure as input to the model; but I doubt that's the case, they evidently meant that you should train a model on 3D data so that it will later be able to predict other 3D structures from the 2D information.

$\endgroup$
8
  • $\begingroup$ Thank you for the very insightful answer! I'm not sure how precise I need actually -- I am working on a protein-ligand docking problem and it seems natural to model 3D molecular structure in addition to 3D protein structure. I am pretty sure I need just a conformation since I can inject equivariances into the molecular encoder. $\endgroup$
    – jasperhyp
    May 4 at 18:58
  • $\begingroup$ "If you just need to represent the molecules in 3D in a conformation that "makes sense" for small molecules, a force field based method will generally work and produce conformations that work well enough" -- this sounds very encouraging! Perhaps I'll start with those structures. Thanks! Re. the drugbank comment, I am not sure if the 3d structures for drugs are sort of more experimentally-based since many drugs are pretty well-studied. Have sent an inquiry to them. Totally agree that 3D QSAR might be too hard, but I guess my use case is 3D inherently so that might make sense to include. $\endgroup$
    – jasperhyp
    May 4 at 19:00
  • $\begingroup$ For the website's use of model predictions -- I noticed this "3D molecular structures provided. We additionally provide 3D structures for training molecules. These structures are calculated by DFT and are obtained together with the HOMO-LUMO gap." Does thi mean they still used DFT? $\endgroup$
    – jasperhyp
    May 4 at 19:03
  • $\begingroup$ I'm not sure what you're trying to do with protein ligand docking; if you're predicting the conformation of the ligand it would make more sense to use a docked conformation during the training. You definitely can use 2D otherwise. I'm not sure what you mean by injecting equivariance; molecular conformations can change a lot depending on the conditions. In that website they build a model which predicts the minimum DFT conformation. They do indeed compute it with DFT for the training data, there's no way around it. Unfortunately, they do not explain HOW they obtain such structures... $\endgroup$
    – stanton63
    May 5 at 13:01
  • 1
    $\begingroup$ @jasperhyp if all conformations have the same representation; then you're not really providing any conformational information to your model, thus I don't see how that would improve the model. Having to go through the effort of optimizing the structure through DFT can take time, you should also confirm that the effort is worth it due to speed or efficiency when compared to other methods such as molecular docking or molecular dynamics. If your model requires significant amounts of preprocessing time, why not just use those methods? $\endgroup$
    – stanton63
    May 8 at 7:26

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.