Conformational analysis of perezone and dihydroperezone using vibrational circular dichroism

The vibrational circular dichroism (VCD) spectra of perezone and dihydroperezone measured from CDCl₃ solutions were quite similar, suggesting analogous conformations for both molecules. Their absolute configurations were confirmed by comparison of the experimental VCD spectrum of each compound with curves generated from theoretical calculations using density functional theory (DFT) at the B3LYP/DGDZVP level of theory taking into account their conformational mobility. Conformational analysis of the 8-(R) enantiomer showed 19 low energy conformers in a 2.4 kcal/mol energy range, while for 8-(R), with the saturated side alkyl chain, 34 conformers were considered in the first 2 kcal/mol. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field, all MMFF94 conformers were geometrically optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by reoptimization and calculations of their vibrational frequencies at the B3LYP/DGDZVP level. Good agreement between the theoretical 8-(R) enantiomers and experimental VCD curves were observed for both.     

Structure and absolute configuration of ginkgolide B characterized by IR‐ and VCD spectroscopy

Experimental and calculated (B3LYP/6‐31G(d)) vibrational circular dichroism (VCD) and IR spectra are compared, illustrating that the structure and absolute configuration of ginkgolide B (GB) may be characterized directly in solution. A conformational search for GB using MacroModel and subsequent DFT optimizations (B3LYP/6‐31G(d)) provides a structure for the lowest energy conformer which agrees well with the structure determined by X‐ray diffraction. In addition, a conformer at an energy of 7 kJ mol^?1 (B3LYP/6‐311+G(2d,2p)) with respect to the lowest energy conformer is predicted, displaying different intramolecular hydrogen bonding. Differences between measured and calculated IR and VCD spectra for GB at certain wavenumbers are rationalized in terms of interactions with solvent, intermolecular GB‐GB interactions, and the potential presence of more than one conformer. This is the first detailed investigation of the spectroscopic fingerprint region (850?1300 cm^?1) of the natural product GB employing infrared absorption and VCD spectroscopy. Chirality, 2010. © 2009 Wiley‐Liss, Inc.     

Determination of the absolute configurations of chiral organometallic complexes via density functional theory calculations of their vibrational circular dichroism spectra: The chiral chromium tricarbonyl complex of N-pivaloyl-tetrahydroquinoline

The racemate of the chiral tricarbonyl-η⁶-arene-chromium(0) complex, tricarbonyl-η⁶-N-pivaloyl-tetrahydroquinoline-chromium(0), 1, has been synthesized and resolved using chromatography on a (R,R)-Whelk-O1 column. The Absolute Configuration (AC) of 1 has been determined using vibrational circular dichroism (VCD) spectroscopy. The VCD spectrum of 1 has been predicted using the Stephens equation for vibrational rotational strengths, implemented using density functional theory (DFT) in the gaussian program. Using the B3PW91 functional and the 6-311++G (2d,2p) basis set, the predicted VCD spectrum of S-1 is in excellent agreement with the experimental VCD spectrum of (+)-1, leading unambiguously to the AC S-(+). It is concluded that VCD is a useful technique for determining the ACs of chiral organometallic complexes, given the use of optimum functionals and basis sets.     

Absolute configuration of eremophilanoids by vibrational circular dichroism

The absolute configurations (AC) of natural occurring 6-hydroxyeuryopsin (1), of its acetyl derivative 2, and of eremophilanolide 8 were confirmed by comparison of the experimental vibrational circular dichroism (VCD) spectra with theoretical curves generated from density functional theory (DFT) calculations. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field. All MMFF94 conformers were further optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by calculations of their vibrational frequencies at the B3LYP/6-31G(d,p); the VCD spectra of 2 and 8 were also calculated at the B3PW91/DGDZVP level of theory. Good agreement between theoretical and experimental VCD curves unambiguously verified the 4S,5R,6S absolute configuration for 1 and 2, and the 1S,4S,5R,6S,8S,10S configuration for 8.     

Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.     

Comparative theoretical study of the UV/Vis absorption spectra of styrylpyridine compounds using TD-DFT calculations

This study examined absorption properties of 2-styrylpyridine, trans-2-(m-cyanostyryl)pyridine, trans-2-[3-methyl-(m-cyanostyryl)]pyridine, and trans-4-(m-cyanostyryl)pyridine compounds based on theoretical UV/Vis spectra, with comparisons between time-dependent density functional theory (TD-DFT) using B3LYP, PBE0, and LC-ωPBE functionals. Basis sets 6–31G(d), 6–31G(d,p), 6–31+G(d,p), and 6–311+G(d,p) were tested to compare molecular orbital energy values, gap energies, and maxima absorption wavelengths. UV/Vis spectra were calculated from fully optimized geometry in B3LYP/6–311+G(d,p) in gas phase and using the IEFPCM model. B3LYP/6–311+G(d,p) provided the most stable form, a planar structure with parameters close to 2-styrylpyridine X-ray data. Isomeric structures were evaluated by full geometry optimization using the same theory level. Similar energetic values were found: ～4.5 kJ?mol^?1 for 2-styrylpyridine and ～1 kJ?mol^?1 for derivative compound isomers. The 2-styrylpyridine isomeric structure differed at the pyridine group N-atom position; structures considered for the other compounds had the cyano group attached to the phenyl ring m-position equivalent. The energy difference was almost negligible between m-cyano-substituted molecules, but high energy barriers existed for cyano-substituted phenyl ring torsion. TD-DFT appeared to be robust and accurate approach. The B3LYP functional with the 6–31G(d) basis set produced the most reliable λ_max values, with mean errors of 0.5 and 12 nm respect to experimental values, in gas and solution, respectively. The present data describes effects on the λ_max changes in the UV/Vis absorption spectra of the electron acceptor cyano substituent on the phenyl ring, the electron donor methyl substituent, and the N-atom position on the electron acceptor pyridine ring, causing slight changes respect to the 2-styrylpyridine title compound.     

Synthesis,characterization and absolute configurations of methyl ladderanoates

Methyl esters of [5]-ladderanoic acid and [3]-ladderanoic acid were prepared by esterification of the acids isolated from biomass at a wastewater treatment plant. Optical rotations at six different wavelengths (633, 589, 546, 436, 405 and 365 nm) and vibrational circular dichroism (VCD) spectra in the 1800–900 cm⁻¹ region were measured in CDCl₃ solvent and compared with quantum chemical (QC) predictions using B3LYP functional and 6-311++G(2d,2p) basis set with polarizing continuum model representing the solvent. QC predictions gave negative optical rotations at all six wavelengths for (R)-methyl [5]-ladderanoate and positive optical rotations for (R)-methyl [3]-ladderanoate, the same signs as previously reported for the corresponding acids. The crystal structure of (−)-methyl [5]-ladderanoate independently confirmed (R) configuration. The QC-predicted VCD spectra using Boltzmann population weighted spectra of individual conformers did not provide satisfactory quantitative agreement with the experimental VCD spectra. An improved quantitative agreement for VCD spectra could be obtained when conformer populations were optimized to maximize the similarity between experimental and predicted VCD spectra, but more improvements in VCD predictions are needed.     

Vibrational Circular Dichroism Discrimination of Diastereomeric Cedranol Acetates

The reliability of vibrational circular dichroism (VCD) spectroscopy to discriminate four diastereomeric cedranol acetates 1 , 2 , 3 , 4 by means of their absolute configuration is examined. The usage of CompareVOA software to quantify comparisons of the measured infrared (IR) and VCD spectra with the corresponding simulated spectra at the B3LYP/DGDZVP and B3PW91/DGDZVP levels of theory for each diastereomer enabled the B3PW91 functional to be qualified as superior to the B3LYP functional for vibrational calculations of 1 , 2 , 3 , 4 . Analogously, a set of quantitative VCD spectra cross‐comparisons of 1 , 2 , 3 , 4 unambiguously distinguished the diastereomers using B3PW91 and failed using B3LYP. Remarkably, quantitative IR spectra cross‐comparisons of 1 , 2 , 3 , 4 using B3PW91 or B3LYP functionals demonstrated that the achiral spectroscopic IR technique is not able to distinguish cedranol acetate diastereomers. VCD comparisons using anisotropy g‐factor values of bands in the 1550–950 cm^‐1 region of the spectra were of aid to facilitate visual spectra matching for each diastereomer. Chirality 25:939‐951, 2013. © 2013 Wiley Periodicals, Inc.     

Absolute configuration of tropane alkaloids from Schizanthus species by vibrational circular dichroism

The absolute configuration (AC) of 6β-hydroxy-3α-senecioyloxytropane (1), 3α-hydroxy-6β-tigloyloxytropane (2), 3α-hydroxy-6β-senecioyloxytropane (3), and 3α-hydroxy-6β-angeloyloxytropane (4) was assigned as (1R,3R,5S,6R) using density functional theory (DFT) calculations at the B3LYP/DGDZVP level of theory in combination with experimental vibrational circular dichroism (VCD) measurements and comparison with the spectra of similar tropanes. The AC of 1 followed from a sample isolated from Schizanthus grahamii, while those of the mixture of 2 and 3, isolated from the same source, were determined by comparing the VCD measurement to a weighted calculation of the individual VCD spectra according to a 69:31 ratio of 2:3 determined by ¹H NMR signal integration. In turn, Schizanthus pinnatus provided a 7:3 mixture of 1:4 whose AC was determined using the experimental VCD absorptions in the 1150-950 cm⁻¹ spectral region which were compared with those observed for 1-3 and with those described for other 3α,6β-tropanediol derivatives.     

Density functional theory investigation of electrophilic addition reaction of bromine to tricyclo[4.2.2.22,5]dodeca-1,5-diene1

The geometry and the electronic structure of tricyclo[4.2.2.2^2,5]dodeca-1,5-diene (TCDD) molecule were investigated by DFT/B3LYP and /B3PW91 methods using the 6-311G(d,p) and 6-311++G(d,p) basis sets. The double bonds of TCDD molecule are syn-pyramidalized. The structure of π-orbitals and their mutual interactions for TCDD molecule were investigated. Potential energy surface (PES) of the TCDD-Br₂ system was studied by B3LYP/6-311++G(d,p) method and the configurations [molecular charge-transfer (CT) complex, transition states (TS1 and TS2), intermediate (INT) and product (P)] corresponding to the stationary points (minima or saddle points) were determined. Initially, a molecular CT-complex forms between Br₂ and TCDD. With a barrier of 22.336 kcal mol^-1 the CT-complex can be activated to an intermediate (INT) with energy 15.154 kcal mol^-1 higher than that of the CT-complex. The intermediate (INT) then transforms easily (barrier 5.442 kcal mol^-1) into the final, N-type product. The total bromination is slightly exothermic. Accompanying the breaking of Br-Br bond, C1-Br, C5-Br and C2-C6 bonds are formed, and C1 = C2 and C5 = C6 double bonds transform into single bonds. The direction of the reaction is determined by the direction of intramolecular skeletal rearrangement that is realized by the formation of C2-C6 bond.

An assessment to evaluate the validity of different methods for the description of some corrosion inhibitors

New research and development efforts using computational chemistry in studying an assessment of the validity of different quantum chemical methods to describe the molecular and electronic structures of some corrosion inhibitors were introduced. The standard and the highly accurate CCSD method with 6-311++G(d,p), ab initio calculations using the HF/6-31G++(d,p) and MP2 with 6-311G(d,p), 6-31++G(d,p), and 6-311++G(2df,p) methods as well as DFT method at the B3LYP, BP86, B3LYP*, M06L, and M062x/6-31G++(d,p) basis set level were performed on some triazole derivatives and sulfur containing compounds used as corrosion inhibitors. Quantum chemical parameters, such as the energy of the highest occupied molecular orbital energy (E_HOMO), the energy of the lowest unoccupied molecular orbital energy (E_LUMO), energy gap (ΔE), dipole moment (μ), sum of total negative charges (TNC), chemical potential (Pi), electronegativity (χ), hardness (η), softness (σ), local softness (s), Fukui functions (f ⁺,f ^?), electrophilicity (ω), the total energy change (?E_T) and the solvation energy (S.E), were calculated. Furthermore, the accuracy and the applicability of these methods were estimated relative to the highest accuracy and standard CCSD with 6-311++G(d,p) method. Good correlations between the quantum chemical parameters and the corresponding inhibition efficiency (IE%) were found.     

Zwitterionic conformers of pyrrolysine and their interactions with metal ions—a theoretical study

A total of 16 pyrrolysine conformers in their zwitterionic forms are studied in gas and simulated aqueous phase using a polarizable continuum model (PCM). These conformers are selected on the basis of our study on the intrinsic conformational properties of non-ionic pyrrolysine molecule in gas phase [Das and Mandal (2013) J Mol Model 19:1695?1704]. In aqueous phase, the stable zwitterionic pyrrolysine conformers are characterized by full geometry optimization and vibrational frequency calculations using B3LYP/6-311++G(d,p) level of theory. Single point calculations are also carried out at MP2/6-311++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. The calculated relative energy range of the conformers at B3LYP/6-311++G(d,p) level is 5.19 kcal mol^?1 whereas the same obtained by single point calculations at MP2/6-311++G(d,p) level is 4.58 kcal mol^?1. A thorough analysis reveals that four types of intramolecular H-bonds are present in the conformers; all of which play key roles in determining the energetics and in imparting the observed conformations to the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of the H-bonds. This study also points out that conformers with diverse structural motifs may differ in their thermodynamical stability by a narrow range of relative energy. The effects of metal coordination on the relative stability order and structural features of the conformers are examined by complexing five zwitterionic conformers of pyrrolysine with Cu⁺² through their carboxylate groups. The interaction enthalpies and Gibbs energies, rotational constants, vibrational frequencies and dipole moments of the metal complexes calculated at B3LYP level are also reported. The zwitterionic conformers of pyrrolysine are not stable in gas phase; after geometry optimization they are converted to the non-ionic forms.     

Studies on molecular structure and tautomerism of a vitamin B<Subscript>6</Subscript> analog with density functional theory

This work presents a computational study on the molecular structure and tautomeric equilibria of a novel Schiff base L derived from pyridoxal (PL) and o-phenylenediamine by using the density functional method B3LYP with basis sets 6-31 G(d,p), 6-31++G(d,p), 6-311 G(d,p) and 6-311++G(d,p). The optimized geometrical parameters obtained by B3LYP/6-31 G(d,p) method showed the best agreement with the experimental values. Tautomeric stability study of L inferred that the enolimine form is more stable than its ketoenamine form in both gas phase and solution. However, protonation of the pyridoxal nitrogen atom (LH) have accelerated the formation of ketoenamine form, and therefore, both ketoenamine and enolimine forms could be present in acidic media.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Structure and electronic properties of azadirachtin

We performed a combined DFT and Monte Carlo ¹³C NMR chemical-shift study of azadirachtin A, a triterpenoid that acts as a natural insect antifeedant. A conformational search using a Monte Carlo technique based on the RM1 semiempirical method was carried out in order to establish its preferred structure. The B3LYP/6-311++G(d,p), wB97XD/6-311++G(d,p), M06/6-311++G(d,p), M06-2X/6-311++G(d,p), and CAM-B3LYP/6-311++G(d,p) levels of theory were used to predict NMR chemical shifts. A Monte Carlo population-weighted average spectrum was produced based on the predicted Boltzmann contributions. In general, good agreement between experimental and theoretical data was obtained using both methods, and the ¹³C NMR chemical shifts were predicted highly accurately. The geometry was optimized at the semiempirical level and used to calculate the NMR chemical shifts at the DFT level, and these shifts showed only minor deviations from those obtained following structural optimization at the DFT level, and incurred a much lower computational cost. The theoretical ultraviolet spectrum showed a maximum absorption peak that was mainly contributed by the tiglate group.     

Probing the antioxidant potential of phloretin and phlorizin through a computational investigation

The structures and energetics of two dihydrochalcones (phloretin and its glycoside phlorizin) were examined with density functional theory, using the B3LYP, M06-2X, and LC-ω PBE functionals with both the 6-311G(d,p) and 6-311 + G(d,p) basis sets. Properties connected to antioxidant activity, i.e., bond dissociation enthalpies (BDEs) for OH groups and ionization potentials (IPs), were computed in a variety of environments including the gas-phase, n-hexane, ethanol, methanol, and water. The smallest BDEs among the four OH groups for phloretin (three for phlorizin) were determined (using B3LYP/6-311 + G(d,p) in water) to be 79.36 kcal/mol for phloretin and 79.98 kcal/mol for phlorizin while the IPs (at the same level of theory) were obtained as 139.48 and 138.98 kcal/mol, respectively. By comparing with known antioxidants, these values for the BDEs indicate both phloretin and phlorizin show promise for antioxidant activity. In addition, the presence of the sugar moiety has a moderate (0-6 kcal/mol depending on functional) effect on the BDEs for all OH groups. Interestingly, the BDEs suggest that (depending on the functional chosen) the sugar moiety can lead to an increase, decrease, or no change in the antioxidant activity. Therefore, further experimental tests are encouraged to understand the substituent effect on the BDEs for phloretin and to help determine the most appropriate functional to probe BDEs for dihydrochalcones.     

Molecular structure,vibrational spectra and HOMO,LUMO analysis of 4-piperidone by density functional theory and ab initio Hartree–Fock calculations

Vibrational frequencies and geometrical parameters of 4-piperidone (4-PID) in the ground state have been calculated by using the Hartree–Fock (HF) and density functional methods (B3LYP) with 6-311++G(d,p) and 6-311+G(3df,2p) basis sets. These methods are proposed as a tool to be applied in the structural characterisation of 4-PID (C₅H₉NO). The title molecule has C _s point group symmetry, thus providing useful support in the interpretation of experimental IR and Raman data. The DFT-B3LYP/6-311+G(3df,2p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 4-PID. The calculated highest occupied molecular orbital and lowest unoccupied molecular orbital energies show that charge transfer occurs within the molecule. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title molecule have been constructed.     

A computational investigation on the antioxidant potential of myricetin 3,4′-di-<Emphasis Type="Italic">O</Emphasis>-<Emphasis Type="Italic">α</Emphasis>-<Emphasis Type="Italic">L</Emphasis>-rhamnopyranoside

In this work, we present a computational study on the antioxidant potential of myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside (Compound M). A density functional theory (DFT) approach with the B3LYP and LC-ωPBE functionals and with both the 6-311G(d,p) and 6-311+G(d,p) basis sets was used. The focus of the investigation was on the structural and energetic parameters including both bond dissociation enthalpies (BDEs) and ionization potentials (IPs), which provide information on the potential antioxidant activity. The properties computed were compared with BDEs and IPs available in the literature for myricetin, a compound well known for presenting antioxidant activity (and the parent molecule of the compound of interest in the present work). Myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside presented the lowest BDE to be 79.13 kcal/mol (as determined using B3LYP/6-311G(d,p) in water) while myricetin has a quite similar value (within 3.4 kcal/mol). IPs computed in the gas phase [B3LYP/6-311G(d,p)] are 157.18 and 161.4 kcal/mol for myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside and myricetin, respectively. As the values of BDEs are considerably lower than the ones probed for IPs (in the gas phase or in any given solvent environment), the hydrogen atom transfer mechanism is preferred over the single electron transfer mechanism. The BDEs obtained suggest that myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside can present antioxidant potential as good as the parent molecule myricetin (a well-known antioxidant). Therefore, experimental tests on the antioxidant activity of Compound M are encouraged.     

Racemization of the gastrointestinal antisecretory chiral drug esomeprazole magnesium via the pyramidal inversion mechanism: A theoretical study

The pyramidal inversion mechanisms of the 6‐methoxy and the 5‐methoxy tautomers of (S)‐omeprazole were studied, employing ab initio and DFT methods. The conformational space of the model molecule (S)‐2‐[(3‐methyl‐2‐pyridinyl)methyl]sulfinyl‐1H‐benzimidazole was calculated, with respect to rotations around single bonds, at the B3LYP/6‐311G(d,p) level. All of the resulting conformations were used as starting points for full optimizations of (S)‐omeprazole, at B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), B3LYP/6‐311++G(d,p), B3LYP/6‐311G(2df,2pd), MP2/6‐31G(d), and MP2/6‐311G(d,p) levels. Four distinct pathways were found for enantiomerization via the pyramidal inversion mechanism for each of the tautomers of (S)‐omeprazole. Each transition state, in which the sulfur, the oxygen and the two carbon atoms connected directly to the sulfur are in one plane, connects two diastereomeric minima. The enantiomerization is completed by free rotation around the sulfur–methylene bond, and around the methylene–pyridine ring bond. The effective Gibbs' free energy barrier for racemization ΔG of the two tautomers of (S)‐omeprazole are 39.8 kcal/mol (5‐methoxy tautomer) and 40.0 kcal/mol (6‐methoxy tautomer), indicating that the enantiomers of omeprazole are stable at room temperature (in the gas phase). The 5‐methoxy tautomer of (S)‐omeprazole was found to be slightly more stable than the 6‐methoxy tautomer, in the gas phase. The energy barrier (ΔG^?) for the(S,M) (S,P) diastereomerization of (S)‐omeprazole due to the rotation around the pyridine chiral axis was very low, 5.8 kcal/mole at B3LYP/6‐311G(d,p). Chirality 2010. © 2010 Wiley‐Liss, Inc.