예제 #1
0
        /// <summary>  Generate canonical SMILES from the <code>molecule</code>. This method
        /// canonicaly lables the molecule but dose not perform any checks on the
        /// chemical validity of the molecule. Does not care about multiple molecules.
        /// IMPORTANT: A precomputed Set of All Rings (SAR) can be passed to this 
        /// SmilesGenerator in order to avoid recomputing it. Use setRings() to 
        /// assign the SAR.
        /// 
        /// </summary>
        /// <param name="molecule">                The molecule to evaluate
        /// </param>
        /// <param name="chiral">                  true=SMILES will be chiral, false=SMILES
        /// will not be chiral.
        /// </param>
        /// <param name="doubleBondConfiguration"> Description of Parameter
        /// </param>
        /// <returns>                          Description of the Returned Value
        /// </returns>
        /// <exception cref="CDKException">        At least one atom has no Point2D;
        /// coordinates are needed for creating the chiral smiles. This excpetion
        /// can only be thrown if chiral smiles is created, ignore it if you want a
        /// non-chiral smiles (createSMILES(AtomContainer) does not throw an
        /// exception).
        /// </exception>
        /// <seealso cref="org.openscience.cdk.graph.invariant.CanonicalLabeler.canonLabel(IAtomContainer)">
        /// </seealso>
        //UPGRADE_NOTE: Synchronized keyword was removed from method 'createSMILESWithoutCheckForMultipleMolecules'. Lock expression was added. "ms-help://MS.VSCC.v80/dv_commoner/local/redirect.htm?index='!DefaultContextWindowIndex'&keyword='jlca1027'"
        public virtual System.String createSMILESWithoutCheckForMultipleMolecules(IMolecule molecule, bool chiral, bool[] doubleBondConfiguration)
        {
            lock (this)
            {
                if (molecule.AtomCount == 0)
                {
                    return "";
                }
                canLabler.canonLabel(molecule);
                brokenBonds.Clear();
                ringMarker = 0;
                IAtom[] all = molecule.Atoms;
                IAtom start = null;
                for (int i = 0; i < all.Length; i++)
                {
                    IAtom atom = all[i];
                    if (chiral && atom.getPoint2d() == null)
                    {
                        throw new CDKException("Atom number " + i + " has no 2D coordinates, but 2D coordinates are needed for creating chiral smiles");
                    }
                    //System.out.println("Setting all VISITED flags to false");
                    atom.setFlag(CDKConstants.VISITED, false);
                    if ((long)((System.Int64)atom.getProperty("CanonicalLable")) == 1)
                    {
                        start = atom;
                    }
                }

                //detect aromaticity
                if (rings == null)
                {
                    if (ringFinder == null)
                    {
                        ringFinder = new AllRingsFinder();
                    }
                    rings = ringFinder.findAllRings(molecule);
                }
                HueckelAromaticityDetector.detectAromaticity(molecule, rings, false);
                if (chiral && rings.AtomContainerCount > 0)
                {
                    System.Collections.ArrayList v = RingPartitioner.partitionRings(rings);
                    //System.out.println("RingSystems: " + v.size());
                    for (int i = 0; i < v.Count; i++)
                    {
                        int counter = 0;
                        IAtomContainer allrings = RingSetManipulator.getAllInOneContainer((IRingSet)v[i]);
                        for (int k = 0; k < allrings.AtomCount; k++)
                        {
                            if (!BondTools.isStereo(molecule, allrings.getAtomAt(k)) && hasWedges(molecule, allrings.getAtomAt(k)) != null)
                            {
                                IBond bond = molecule.getBond(allrings.getAtomAt(k), hasWedges(molecule, allrings.getAtomAt(k)));
                                if (bond.Stereo == CDKConstants.STEREO_BOND_UP)
                                {
                                    allrings.getAtomAt(k).setProperty(RING_CONFIG, UP);
                                }
                                else
                                {
                                    allrings.getAtomAt(k).setProperty(RING_CONFIG, DOWN);
                                }
                                counter++;
                            }
                        }
                        if (counter == 1)
                        {
                            for (int k = 0; k < allrings.AtomCount; k++)
                            {
                                allrings.getAtomAt(k).setProperty(RING_CONFIG, UP);
                            }
                        }
                    }
                }

                System.Text.StringBuilder l = new System.Text.StringBuilder();
                createSMILES(start, l, molecule, chiral, doubleBondConfiguration);
                rings = null;
                return l.ToString();
            }
        }