1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
//
// Copyright (C) 2025 Greg Landrum and other RDKit contributors
// @@ All Rights Reserved @@
// This file is part of the RDKit.
// The contents are covered by the terms of the BSD license
// which is included in the file license.txt, found at the root
// of the RDKit source tree.
//
#include <string>
#include "RDGeneral/test.h"
#include <catch2/catch_all.hpp>
#include <GraphMol/RDKitBase.h>
#include <GraphMol/Chirality.h>
#include <GraphMol/FileParsers/FileParsers.h>
#include <GraphMol/FileParsers/MolSupplier.h>
#include <GraphMol/SmilesParse/SmilesWrite.h>
#include <GraphMol/test_fixtures.h>
using namespace RDKit;
TEST_CASE("explicit tests") {
// making some of the other atropisomer tests explicitly test what we expect
std::string rdbase = getenv("RDBASE");
REQUIRE(rdbase != "");
SECTION("basics") {
std::string file =
rdbase +
"/Code/GraphMol/FileParsers/test_data/atropisomers/AtropManyChirals.sdf";
auto mol = v2::FileParsers::MolFromMolFile(file);
REQUIRE(mol);
CHECK(mol->getBondWithIdx(7)->getStereo() ==
Bond::BondStereo::STEREOATROPCCW);
}
SECTION("macrocycle") {
std::string file =
rdbase +
"/Code/GraphMol/FileParsers/test_data/atropisomers/macrocycle-5-meta-Cl-ortho-hash.sdf";
auto mol = v2::FileParsers::MolFromMolFile(file);
REQUIRE(mol);
CHECK(mol->getBondWithIdx(15)->getStereo() ==
Bond::BondStereo::STEREOATROPCW);
}
SECTION("with stereo groups and other chiral centers") {
REQUIRE(rdbase != "");
std::string file =
rdbase +
"/Code/GraphMol/FileParsers/test_data/atropisomers/AtropManyChiralsEnhanced.sdf";
auto mol = v2::FileParsers::MolFromMolFile(file);
REQUIRE(mol);
CHECK(mol->getBondWithIdx(7)->getStereo() ==
Bond::BondStereo::STEREOATROPCCW);
CHECK(mol->getAtomWithIdx(6)->getChiralTag() ==
Atom::ChiralType::CHI_UNSPECIFIED);
CHECK(mol->getAtomWithIdx(9)->getChiralTag() !=
Atom::ChiralType::CHI_UNSPECIFIED);
CHECK(mol->getAtomWithIdx(13)->getChiralTag() !=
Atom::ChiralType::CHI_UNSPECIFIED);
CHECK(mol->getAtomWithIdx(16)->getChiralTag() !=
Atom::ChiralType::CHI_UNSPECIFIED);
CHECK(mol->getAtomWithIdx(10)->getChiralTag() !=
Atom::ChiralType::CHI_UNSPECIFIED);
REQUIRE(mol->getStereoGroups().size() == 2);
{
const auto &sg = mol->getStereoGroups()[0];
CHECK(sg.getAtoms().size() == 2);
CHECK(sg.getBonds().size() == 1);
}
{
const auto &sg = mol->getStereoGroups()[1];
CHECK(sg.getAtoms().size() == 2);
CHECK(sg.getBonds().size() == 0);
}
}
}
TEST_CASE("atropisomerism and dependent chirality") {
SECTION("chirality allows stereoisomers") {
{
auto m = "C=C(C)C1=CC([C@H](C)O)=CC([C@@H](C)O)=C1 |wD:3.14|"_smiles;
REQUIRE(m);
CHECK(m->getBondWithIdx(2)->getStereo() ==
Bond::BondStereo::STEREOATROPCW);
}
#if 0
// FIX: this test does not work, this is Github #8340
{
// the two stereocenters are identical, we should not have atropisomerism
auto m = "C=C(C)C1=CC([C@H](C)O)=CC([C@H](C)O)=C1 |wD:3.14|"_smiles;
REQUIRE(m);
CHECK(m->getBondWithIdx(2)->getStereo() == Bond::BondStereo::STEREONONE);
}
#endif
}
SECTION("atropisomers allow chirality") {
#if 0
// FIX: this test does not work
{
// here the atropisomeric bonds have opposite stereochemistry, so the
// central atom should be chiral, this is Github #8341
auto m =
"C=C(C1C=C(C)C=CC=1)[C@H](C)C(=C)C1C=C(C)C=CC=1 |wD:2.2,wD:13.14|"_smiles;
REQUIRE(m);
CHECK(m->getBondWithIdx(1)->getStereo() ==
Bond::BondStereo::STEREOATROPCCW);
CHECK(m->getBondWithIdx(12)->getStereo() ==
Bond::BondStereo::STEREOATROPCCW);
CHECK(m->getAtomWithIdx(9)->getChiralTag() !=
Atom::ChiralType::CHI_UNSPECIFIED);
}
#endif
{
// here the atropisomeric bonds have the same stereochemistry, so the
// central atom should not be chiral
auto m =
"C=C(C1C=C(C)C=CC=1)[C@H](C)C(=C)C1C=C(C)C=CC=1 |wD:2.2,wU:13.14|"_smiles;
REQUIRE(m);
CHECK(m->getBondWithIdx(1)->getStereo() ==
Bond::BondStereo::STEREOATROPCCW);
CHECK(m->getBondWithIdx(12)->getStereo() ==
Bond::BondStereo::STEREOATROPCW);
CHECK(m->getAtomWithIdx(9)->getChiralTag() ==
Atom::ChiralType::CHI_UNSPECIFIED);
}
}
}
|
