public static string RDKit_Smiles(string text,
object isomericSmiles = null, bool kekuleSmiles = false,
object rootedAtAtom = null, object canonical = null,
bool allBondsExplicit = false, bool allHsExplicit = false,
bool doRandom = false)
{
var _isomericSmiles = ExcelDnaUtility.ToBoolean(isomericSmiles, true);
var _rootedAtAtom = ExcelDnaUtility.ToInt32(rootedAtAtom, -1);
var _canonical = ExcelDnaUtility.ToBoolean(canonical, true);
if (_isomericSmiles &&
!kekuleSmiles &&
_rootedAtAtom == -1 &&
_canonical &&
!allBondsExplicit &&
!allHsExplicit &&
!doRandom)
{
return(RDKit_CalcDesc(text, "RDKit_Smiles", mol => RDKit.Chem.MolToSmiles(mol)));
}
{
ROMol mol;
if (kekuleSmiles)
{
var _mol = RWMol.MolFromSmiles(text);
RDKFuncs.Kekulize(_mol);
mol = _mol;
}
else
{
mol = Parse(text);
}
return(RDKit.Chem.MolToSmiles(mol, _isomericSmiles, kekuleSmiles, _rootedAtAtom, _canonical, allBondsExplicit, allHsExplicit, doRandom));
}
}
static void Main()
{
// ----- Object creation -----
Console.WriteLine("Creating some objects:");
ROMol m1 = RWMol.MolFromSmiles("c1ccccc1");
Console.WriteLine(" mol: " + m1 + " " + m1.getNumAtoms());
ROMol m2 = RWMol.MolFromSmiles("c1ccccn1");
Console.WriteLine(" smi: " + m1 + " " + m1.MolToSmiles());
Console.WriteLine(" smi2: " + m2 + " " + m2.MolToSmiles());
ExplicitBitVect fp1 = RDKFuncs.LayeredFingerprintMol(m1);
ExplicitBitVect fp2 = RDKFuncs.LayeredFingerprintMol(m2);
Console.WriteLine(" sim: " + RDKFuncs.TanimotoSimilarityEBV(fp1, fp2));
//rxnTest();
//smiTest();
morganTest();
Console.WriteLine("Goodbye");
}
/// <summary>
/// Tries to parses text as SMILES, InChI, or MolBlock and returns RDKit ROMol.
/// </summary>
/// <param name="molIdent">Text expression of molecule like SMILES, InChI, or MolBlock.</param>
/// <returns>RDKit ROMol.</returns>
public static ROMol Parse(string molIdent)
{
if (molIdent == null)
{
throw new ArgumentNullException(nameof(molIdent));
}
var mol = MolecularCache.GetOrAdd(molIdent, a_ident =>
{
string notationType = null;
ROMol a_mol;
a_mol = RWMol.MolFromSmiles(molIdent);
if (a_mol != null)
{
notationType = "SMILES";
goto L_Found;
}
using (var rv = new ExtraInchiReturnValues())
{
a_mol = RDKFuncs.InchiToMol(molIdent, rv);
if (a_mol != null)
{
notationType = "InChI";
goto L_Found;
}
}
a_mol = RWMol.MolFromMolBlock(molIdent);
if (a_mol != null)
{
RDKFuncs.assignStereochemistryFrom3D(a_mol);
notationType = "MolBlock";
goto L_Found;
}
L_Found:
if (a_mol == null)
{
a_mol = nullMol;
}
else
{
if (notationType != null)
{
a_mol.setProp("source", notationType);
}
}
return(a_mol);
});
if (object.ReferenceEquals(mol, nullMol))
{
return(null);
}
return(mol);
}
static void Main(string[] args)
{
var smiles = "c1ccccc1C";
var mol = RWMol.MolFromSmiles(smiles);
var re_smi = RDKFuncs.MolToSmiles(mol);
Console.WriteLine($"Hello toluene, {re_smi}.");
}
static void BulkMemoryLeakTest()
{
Console.WriteLine("Bulk memory leak test");
for (uint i = 0; i < 10000; ++i)
{
ROMol m4 = RWMol.MolFromSmiles("Clc1cccc(N2CCN(CCC3CCC(CC3)NC(=O)c3cccs3)CC2)c1Cl");
if ((i % 1000) == 0)
{
Console.WriteLine(" Done: " + i);
}
m4.Dispose();
//GC.Collect();
}
}
// static void rxnTest() {
// Console.WriteLine( "Reaction tests" );
// var rxn = ChemicalReaction.ReactionFromSmarts("[N:1][C:2].[OH][C:3]=[O:4]>>[C:2][N:1][C:3]=[O:4]");
// var amine = RWMol.MolFromSmiles("CCCN");
// var acid = RWMol.MolFromSmiles("C1CC1CC(=O)O");
// ROMol[] rs = {amine,acid};
// ROMol_Vect rv = new ROMol_Vect(rs);
// for(var i=0;i<100000;i++){
// var ps=rxn.runReactants(rv);
// if(i%100 == 0) {
// Console.WriteLine( "\t{0}", i );
// }
// }
// Console.WriteLine( "Goodbye" );
// }
// static void smiTest() {
// Console.WriteLine( "repeatedly from smiles" );
// for(var i=0;i<1000000;i++){
// ROMol m1=RDKFuncs.MolFromSmiles("c1ccccc1");
// if(i%1000 == 0) {
// Console.WriteLine( "\t{0}", i );
// }
// }
// Console.WriteLine( "Goodbye" );
// }
static void morganTest()
{
// ----- Object creation -----
Console.WriteLine("Creating some objects:");
ROMol m1 = RWMol.MolFromSmiles("c1ccccc1");
Console.WriteLine(" mol: " + m1 + " " + m1.getNumAtoms());
ROMol m2 = RWMol.MolFromSmiles("c1ccccn1");
var fp1 = RDKFuncs.MorganFingerprintMol(m1, 2);
var fp2 = RDKFuncs.MorganFingerprintMol(m2, 2);
Console.WriteLine(" sim: " + RDKFuncs.DiceSimilarity(fp1, fp2));
}
public void TestDrawMolecule()
{
var mol = RWMol.MolFromSmiles("c1ccc(C)c(C)c1C");
var drawer = new MolDraw2DCairo(300, 300, -1, -1, true);
drawer.drawOptions().addAtomIndices = true;
drawer.drawMolecule(mol);
drawer.finishDrawing();
var png1 = drawer.getImage().ToArray();
drawer.writeDrawingText("test.png");
byte[] png2 = File.ReadAllBytes("test.png");
Assert.Equal(png1.Length, png2.Length);
for (int i = 0; i < png1.Length; i++)
{
Assert.Equal(png1[i], png2[i]);
}
File.Delete("test.png");
}
static void Demo()
{
var toluene = RWMol.MolFromSmiles("Cc1ccccc1");
var mol1 = RWMol.MolFromMolFile(Path.Combine("Data", "input.mol"));
var stringWithMolData = new StreamReader(Path.Combine("Data", "input.mol")).ReadToEnd();
var mol2 = RWMol.MolFromMolBlock(stringWithMolData);
using (var suppl = new SDMolSupplier(Path.Combine("Data", "5ht3ligs.sdf")))
{
while (!suppl.atEnd())
{
var mol = suppl.next();
if (mol == null)
{
continue;
}
Console.WriteLine(mol.getAtoms().Count);
using (var maccs = RDKFuncs.MACCSFingerprintMol(mol))
{
Console.WriteLine(ToString(maccs));
}
}
}
using (var gzsuppl = new ForwardSDMolSupplier(new gzstream("Data/actives_5ht3.sdf.gz")))
{
while (!gzsuppl.atEnd())
{
var mol = gzsuppl.next();
if (mol == null)
{
continue;
}
Console.WriteLine(mol.getAtoms().Count);
using (var maccs = RDKFuncs.MACCSFingerprintMol(mol))
{
Console.WriteLine(ToString(maccs));
}
}
}
}
static void MakePicture(string smiles, string filename)
{
int width = 200;
int height = 200;
RWMol mol = null;
mol = RWMol.MolFromSmiles(smiles);
if (mol == null)
{
mol = RWMol.MolFromSmarts(smiles);
}
if (mol == null)
{
throw new Exception($"Cannot recognize: '{smiles}'");
}
RDKFuncs.prepareMolForDrawing(mol);
if (filename.EndsWith(".svg"))
{
var view = new MolDraw2DSVG(width, height);
view.drawMolecule(mol);
view.finishDrawing();
using (var w = new StreamWriter(filename))
{
w.Write(view.getDrawingText());
Console.WriteLine($"{filename} is drawn.");
}
}
else if (filename.EndsWith(".png"))
{
var view = new MolDraw2DCairo(width, height);
view.drawMolecule(mol, Path.GetFileNameWithoutExtension(filename));
view.finishDrawing();
view.writeDrawingText(filename);
}
else
{
throw new Exception($"Not supported: {filename}");
}
}
static void Main(string[] args)
{
// ----- Object creation -----
Console.WriteLine("Creating some objects:");
ROMol m1 = RWMol.MolFromSmiles("c1ccccc1");
Console.WriteLine(" mol: " + m1 + " " + m1.getNumAtoms());
ROMol m2 = RWMol.MolFromSmiles("c1ccccn1");
Console.WriteLine(" smi: " + m1 + " " + m1.MolToSmiles());
Console.WriteLine(" smi2: " + m2 + " " + m2.MolToSmiles());
ExplicitBitVect fp1 = RDKFuncs.LayeredFingerprintMol(m1);
ExplicitBitVect fp2 = RDKFuncs.LayeredFingerprintMol(m2);
Console.WriteLine(" sim: " + RDKFuncs.TanimotoSimilarityEBV(fp1, fp2));
//rxnTest();
//smiTest();
//morganTest();
ROMol m3 = RWMol.MolFromSmiles("c1ccccc1");
uint nAtoms = m3.getNumAtoms(true);
Console.WriteLine("Bulk memory leak test");
for (uint i = 0; i < 10000; ++i)
{
ROMol m4 = RWMol.MolFromSmiles("Clc1cccc(N2CCN(CCC3CCC(CC3)NC(=O)c3cccs3)CC2)c1Cl");
if ((i % 1000) == 0)
{
Console.WriteLine(" Done: " + i);
}
m4.Dispose();
//GC.Collect();
}
Console.WriteLine("Goodbye");
}
static void CreateSomeObjects()
{
// ----- Object creation -----
Console.WriteLine("Creating some objects:");
ROMol m1 = RWMol.MolFromSmiles("c1ccccc1");
Console.WriteLine(" mol: " + m1 + " " + m1.getNumAtoms());
ROMol m2 = RWMol.MolFromSmiles("c1ccccn1");
Console.WriteLine(" smi: " + m1 + " " + m1.MolToSmiles());
Console.WriteLine(" smi2: " + m2 + " " + m2.MolToSmiles());
ExplicitBitVect fp1 = RDKFuncs.LayeredFingerprintMol(m1);
ExplicitBitVect fp2 = RDKFuncs.LayeredFingerprintMol(m2);
Console.WriteLine(" sim: " + RDKFuncs.TanimotoSimilarityEBV(fp1, fp2));
const string SMILES_benzene = "c1ccccc1";
ROMol m3 = RWMol.MolFromSmiles(SMILES_benzene);
uint nAtoms = m3.getNumAtoms(true);
Console.WriteLine($"The number of atoms in {SMILES_benzene} is {nAtoms}.");
}
static void ErrorHandle()
{
var invalid = RWMol.MolFromSmiles("dfep3js2g");
}
public static RWMol MolFromSmiles(string smiles, bool sanitize = true)
{
return(RWMol.MolFromSmiles(smiles, 0, sanitize));
}
