WO2015009016A1 - 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템 - Google Patents

분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템 Download PDF

Info

Publication number
WO2015009016A1
WO2015009016A1 PCT/KR2014/006364 KR2014006364W WO2015009016A1 WO 2015009016 A1 WO2015009016 A1 WO 2015009016A1 KR 2014006364 W KR2014006364 W KR 2014006364W WO 2015009016 A1 WO2015009016 A1 WO 2015009016A1
Authority
WO
WIPO (PCT)
Prior art keywords
block
molecular orbital
mbs
molecular
similarity
Prior art date
Application number
PCT/KR2014/006364
Other languages
English (en)
French (fr)
Inventor
이승엽
조혜성
Original Assignee
주식회사 엘지화학
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 주식회사 엘지화학 filed Critical 주식회사 엘지화학
Priority to US14/904,651 priority Critical patent/US10102349B2/en
Priority to JP2015533002A priority patent/JP5933130B2/ja
Publication of WO2015009016A1 publication Critical patent/WO2015009016A1/ko

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C07ORGANIC CHEMISTRY
    • C07BGENERAL METHODS OF ORGANIC CHEMISTRY; APPARATUS THEREFOR
    • C07B61/00Other general methods
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like

Definitions

  • the present invention relates to a method for evaluating molecular orbital similarity deviation and a system using the same. More particularly, the present invention relates to a method for evaluating molecular orbital similarity deviation and a system using the same. will be.
  • the electrochemical properties of the material are greatly influenced by the electron behavior. Because these electronic behaviors are very complex and difficult to interpret, it is difficult to evaluate them by direct measurements through experiments. Therefore, electron behavior is interpreted through calculations, not experiments, and the concept introduced for this is Molecular Orbital (MO).
  • MO Molecular Orbital
  • Molecular orbitals represent the probability that an electron can be present at each position in a molecule. Molecular orbitals can only be identified through calculation by quantum mechanical methods. Traditionally, the evaluation of molecular orbitals relies on a qualitative method of generating a calculated molecular orbital and then relying on visual judgment. However, qualitative methods are useful for roughly evaluating the overall properties of molecular orbitals, but it is difficult to objectively and accurately evaluate molecular orbitals like quantitative methods. Because qualitative methods are evaluated according to subjective criteria, very different evaluations may be made by different evaluators even for the same molecular orbital.
  • the present inventors have overcome the limitations of such a qualitative evaluation method and have developed various methods capable of quantitatively evaluating molecular orbitals for systematically and accurately evaluating molecular orbital information.
  • each subregion of the entire molecular structure of the molecule is formed into blocks, and then the entire molecular structure is generated from a combination of blocks, and (2) the molecular orbitals associated with each block are added to the sum of the total orbitals.
  • the AC2B method which introduces the concept of blocks, simplifies the pattern of complex molecular orbitals and enables intuitive and clear evaluation of molecular orbital properties.
  • the rearranged block spectra computed by AC2B intuitively represent single-molecule orbital properties but are difficult to use directly to compare different molecular orbital properties. This is because the block spectrum is a sequential arrangement of blocks according to positions occupied by each block and thus cannot be immediately used for quantitative comparison.
  • the present inventors intuitively and clearly evaluate the molecular orbital properties distributed in a complex and diverse pattern within the molecular structure by using a new concept of blocks, and use the calculated block spectra to different molecular orbitals.
  • the 2BS-score method is an evaluation method that quantitatively evaluates the similarity of molecular orbitals by evaluating the identity of sequences of block spectra calculated for different molecular orbitals.
  • the 2BS-score method is a method of quantitatively determining the overall molecular orbital similarity by comparing the same sequence in the same sequence in three steps in the block spectrum for the molecular orbitals to be compared. If the sequences of the two molecular orbital block spectra to be compared are exactly the same, the 2BS-score is 100%, and the larger the sequence difference, the smaller the value, indicating poor molecular orbital similarity.
  • the present inventors confirmed that the molecular orbital similarity can be accurately and quantitatively represented by the 2BS-score method.
  • the molecular orbital similarity can be quantitatively evaluated through the 2BS-score method, it is not known whether the similarity has a certain degree of variation. For example, assuming that the similarity between different molecular orbital pairs A0-A1 and A2-A3 is 2BS-score, it is determined that both cases have the same similarity if they have the same 2BS-score. The similarity between and A2-A3 is not exactly the same. This is because the similarity between A0-A1 and A2-A3 may be different even if the 2BS-score values are the same.
  • the similarity deviation can be accurately measured, it is expected that molecular orbital similarity can be more accurately evaluated, and thus it is very useful for understanding similarity characteristics.
  • the inventors have quantitatively evaluated the molecular orbital similarity calculated according to the 2BS-score method. Developed.
  • an object of the present invention is to provide a new method capable of quantitatively and accurately assessing similarity deviations of molecular orbitals.
  • the present invention is a) N x multi-block spectrum (MBS) of different block size by the method of steps i) to v) after selecting two target molecular orbitals to compare the molecular orbital similarity And combining N x MBS pairs of the same block size: i) calculating the molecular orbital distribution using quantum mechanics, ii) radial at the molecular center in the molecular structure of the molecule direction), iii) calculating the molecular orbital ratio (BX ( k )) ( k is the number of blocks and is a natural number from 1 to N) associated with each block, iv) the molecule Sequentially rearranging the blocks based on the size of the orbital ratio BX ( k ) to obtain a block spectrum, and repeating steps ii) to iv), N x having different block sizes for each target molecular orbital.
  • MBS multi-block spectrum
  • the present invention enables the comprehensive evaluation and understanding of the characteristics related to molecular orbital similarity by quantitatively and accurately evaluating the deviations included in the molecular orbital similarity. It can be usefully used for material development.
  • FIG. 1 is a flowchart illustrating a process for evaluating molecular orbital similarity deviation of the present invention.
  • NPB N, N'-Di [(1-naphthyl) -N, N'-diphenyl] -1,1 '-(biphenyl) -4,4'-diamine
  • FIG. 3 shows NPB (N, N'-Di [(1-naphthyl) -N, N'-diphenyl] -1,1 '-(biphenyl) -4,4'-diamine) having different molecular orbitals
  • T1 and T3 calculated by the quantum mechanical method in different electronic states
  • the molecular orbital pairs T1 and T3 with deviations of ⁇ 24.9 obtained by applying the SSDE-MO method are illustrated.
  • the method for evaluating molecular orbital similarity deviation of the present invention includes a) selecting two target molecular orbitals to compare molecular orbital similarities, and then N x multi-block spectrums (MBS) having different block sizes by the method of steps i) to v).
  • MSS multi-block spectrums
  • N x MBS pairs of the same block size i) calculating the molecular orbital distribution using quantum mechanics, ii) radial at the molecular center in the molecular structure of the molecule direction), iii) calculating the molecular orbital ratio (BX ( k )) ( k is the number of blocks and is a natural number from 1 to N) associated with each block, iv) the molecule
  • BX ( k ) the number of blocks and is a natural number from 1 to N
  • the molecule By rearranging the blocks sequentially based on the size of the orbital ratio (BX ( k )) to obtain a block spectrum, and v) repeating steps ii) to iv), the block size is different for each target molecular orbital.
  • N x MBS pairs of the same block size after obtaining N x multi-block spectrums b) performing a multi-step identity evaluation on the N x MBS pairs to calculate a TSS (m) value (m is 1 to N x as the number of MBSs) for each MBS pair; And c) calculating a standard deviation of TSS (m) values for each MBS pair and quantitatively evaluating molecular orbital similarity deviations.
  • the SSDE-MO method is to quantitatively calculate molecular orbital similarity deviations through a statistical estimation method.
  • the SSDE-MO method generates (1) N multi-level block spectra (MBS) with different block sequence sizes to calculate similarity deviations of different molecular orbitals, and (2) evaluates the identity of the generated MBS sequences.
  • the similarity of each MBS pair is calculated through (3) and the similarity deviation is statistically calculated using the calculated similarity data for each MBS pair.
  • the calculation process of SSDE-MO method is shown by the flowchart of FIG. Hereinafter, the calculation process of the SSDE-MO method will be described in detail.
  • N x multi-block spectrums having different block sizes are obtained by the method of steps i) to v), and then the same block size.
  • Step a) is a step of using the "Assembly of Consecutive Building Block” (AC2B) method of the molecular orbital characteristic analysis developed by the present inventor.
  • AC2B Consecutive Building Block
  • Molecular orbitals can be defined as mathematical simulations of the wave-like behavior of electrons in a molecule.
  • the region in which the molecular orbital is present can be obtained through quantum mechanical calculation.
  • the quantum mechanical calculation method is not particularly limited as long as it is a method using quantum mechanics, but preferably an orbital wave at each point calculated in the molecular structure of a material It can be used to calculate through the distribution of the electron density ( ⁇ 2 ), which is the square of the orbital wave function ( ⁇ ), it is also possible to use a single point energy calculation or geometry optimization calculation.
  • the inventors of the present invention calculated molecular orbitals using DMol3 of MATERIAL STUDIO developed by ACCELYS, based on Density Functional Theory (DFT), and used VISUALIZER of MATERIAL STUDIO package to generate molecular orbital pictures. It was.
  • DFT Density Functional Theory
  • the total molecular structure of the target molecule is composed of a combination of N sequential blocks generated on the basis of the molecular center according to the designated total block number N of molecules.
  • the RDM represents a mesh containing all constituent atoms in the molecular structure starting from the center of the molecule in the radial direction.
  • the method for obtaining the intramolecular center (x c , y c , z c ) is shown in Equations 1-1 to 1-3.
  • N Coord represents the total number of atomic coordinates constituting the molecule.
  • the total number N of blocks is not particularly limited but is preferably in the range of 3 to 100.
  • the entire molecular structure of the target molecule can be blocked based on the distance from the molecular center.
  • Step a) may include calculating a molecular orbital ratio BX ( k ) associated with each block.
  • the molecular orbital ratio (BX ( k )) associated with each of the blocks refers to the amount occupied by the molecular orbital associated with the kth block among the total molecular orbitals.
  • K is the number of blocks and a natural number from 1 to N.
  • the molecular orbital ratio (BX ( k )) associated with each of the blocks is calculated by quantum mechanical calculation of the molecular orbital (BMO ( k )) associated with each of the blocks, thereby calculating SUM, which is the total molecular orbital sum, It can be found by calculating the ratio (BX ( k )) of the molecular orbitals associated with each block to the total molecular orbital sum.
  • the step a) may include rearranging the combination of the blocks sequentially based on the size of the molecular orbital ratio BX ( k ) to obtain rearranged block spectra.
  • the "rearranged block spectrum” refers to a combination of sequential blocks obtained by rearranging an assembly of consecutive building blocks (AC2B) obtained in step a) sequentially based on BX ( k ).
  • Step a) includes repeating steps ii) to v) to obtain N x MBSs (MBs) having different block sizes.
  • the different molecular orbitals to compare this process are described as T1 and T2.
  • the default block size for creating a block spectrum (BS) in AC2B is five.
  • the total number N x of MBSs and the corresponding block size set N ( k ) must be determined.
  • N x must be greater than 1 and may have various values depending on the efficiency of the calculation and the purpose of use. For example, if N x is designated as 6, the corresponding block size set may be set as follows.
  • N x MBSs multi-block spectrums having different block sizes for each target molecular orbital can be obtained.
  • MMS pair means a pair consisting of block spectra of the same size combined from N x multi-block spectrum (MBS) obtained for each target molecule orbital.
  • a total of 6 pairs of MBSs having different block sizes may be combined.
  • the present invention includes b) calculating a TSS (m) value (m is 1 to N x as the number of MBS) for each MBS pair by performing a multi-step identity evaluation on the N x MBS pairs. Characterized in that.
  • the multi-stage identity evaluation of step b) may be a step of calculating similarity of the N x MBS pairs by a total of three stages of evaluation of the first to third stages of identity evaluation. In all procedures, identity evaluation is performed based on comparisons in the same sequence of MBS pairs. The following describes the three-step process of the multi-step in detail.
  • the primary identity is evaluated for the total N x MBS pairs by the following sequence comparison.
  • MBS1 ⁇ m, n ⁇ and MBS2 ⁇ m, n ⁇ represent the identity of the block at position n in MBS having a block size of N ⁇ m ⁇ generated in the T1 and T2 molecular orbitals, respectively.
  • SC ⁇ m, n ⁇ represents the degree of identity at the n position for the MBS pair having the size of N ⁇ m ⁇ . This assesses the primary identity at the MBS sequence position.
  • X is greater than 0 and less than 0.6
  • Y is greater than 0 and less than 0.5
  • Z is less than 0 and less than 0.3
  • X + Y + (N-2) ⁇ Z 1.0 is established between X, Y, and Z. .
  • the first step identity evaluation compares each block spectral sequence of the N x MBS pairs at the same sequence position and gives SC (m, n) values if the blocks are identical, but increases n (SC, m, n) It is characterized by giving a low value.
  • the third step identity evaluation sets a block sequence distance (d_BL ⁇ m, n ⁇ ) criterion at each same sequence position of the MBS pair , and the block sequence distance at each position is the block sequence distance d_BL ⁇ l ⁇ .
  • SC (m, n) value is given if the same as, but the SC (m, n) value is n It is characterized in that the value is given to decrease as it increases.
  • step b the similarity (TSS, total similarity score) for each MBS pair is calculated using Equation 1 using the result of the identity evaluation for each MBS pair through the multi-stage identity evaluation.
  • m is 1 to N X as the number of MBS .
  • the TSS ⁇ m ⁇ represents the identity for the MBS pair of the block size N ⁇ m ⁇ .
  • the TSS indicates a value of 100% and the lower the similarity, the smaller the value.
  • the present invention is characterized by including c) calculating the standard deviation of the TSS (m) value for each MBS pair and thereby quantitatively assessing molecular orbital similarity deviations.
  • the present invention also provides a molecular orbital similarity deviation evaluation system using the molecular orbital similarity deviation evaluation method.
  • the molecular orbital similarity deviation evaluation system comprises the steps of: a) selecting two target molecular orbitals to compare molecular orbital similarities, i) calculating the molecular orbital distribution using quantum mechanical calculations, ii) molecules within the molecular structure of the molecule Making N blocks in a radial direction from the center, iii) calculating a molecular orbital ratio (BX ( k )) ( k : natural number from 1 to N) associated with each block, iv) Sequentially rearranging the blocks based on the size of the molecular orbital ratio (BX ( k )) to obtain a block spectrum, and v) repeating steps ii) to iv) of N x MBSs having different block sizes ( Blocking module for obtaining N x MBS (multi-block spectrum) having a different block size through the step of obtaining a multi-block spectrum; b) combining N x MBS pairs of the same block size from each of the N x multi-block
  • a data input module configured to receive data by calculating a TSS (m) value for each MBS pair; And c) an evaluation module for calculating a standard deviation of TSS (m) values for each MBS pair and quantitatively evaluating molecular orbital similarity deviations.
  • module refers to one unit for processing a specific function or operation, which may be implemented by hardware or software or a combination of hardware and software.
  • NPB N, N'-Di [(1-naphthyl) -N, N'-diphenyl] -1,1 '-(biphenyl) -4,4'-diamine
  • TSS values for MBS pairs were obtained by applying the SSDE-MO method of the present invention to T1 and T2, which are molecular orbital pairs calculated by the quantum mechanical method in the state.
  • the final standard deviation was ⁇ 8.6. This is a very small value, indicating little variation in the similarity between T1 and T2 molecular orbitals.
  • molecular orbitals T1 and T2 were visualized using VISUALIZER, a program in the MATERIAL STUDIO package, and are shown in FIG. 2. As can be seen in Figure 2 it can be seen that the molecular orbitals are not the same but very similar. Therefore, even qualitatively confirming that the variation in the similarity between the T1 and T2 molecular orbital is not large. Through this, it can be confirmed that the SSDE-MO method of the present invention quantitatively accurately represents molecular orbital similarity deviations.

Landscapes

  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Complex Calculations (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

본 발명은 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계; b) 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m)(m는 MBS의 번호로서 1 내지 Nx 이다) 값을 계산하는 단계; 및 c) 상기 Nx개의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 단계를 포함하는 분자 오비탈 유사성 편차 평가 방법에 관한 것이다.

Description

분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템
본 발명은 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템에 관한 것으로서, 보다 자세하게는 정량적으로 분자 오비탈(molecular orbital)의 유사성 편차를 정량적으로 정확하게 평가할 수 있는 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템에 관한 것이다.
물질의 전기화학적 특성은 전자 거동에 큰 영향을 받는다. 이러한 전자 거동은 매우 복잡하고 해석하기 어렵기 때문에 이를 실험을 통한 직접적인 측정으로 평가하는 것은 어렵다. 따라서 전자 거동은 실험이 아닌 계산을 통해 해석되고 있고 이를 위해서 도입된 개념이 바로 분자 오비탈(Molecular Orbital, MO)이다.
분자 오비탈은 분자 내 각 위치에서 전자가 존재할 수 있는 확률을 나타낸다. 분자 오비탈은 양자역학 방법에 의한 계산을 통해서만 파악될 수 있다. 전통적으로 분자 오비탈의 평가는 계산된 분자 오비탈을 그림으로 생성한 후 이를 시각적인 판단에 의존하는 정성적인 방법에 의존한다. 그러나 정성적인 방법은 분자 오비탈의 전체적인 특성을 대략적으로 평가하는 데는 유용하지만 정량적인 방법과 같이 분자 오비탈을 객관적이고 정확하게 평가하기는 어렵다. 왜냐하면 정성적인 방법은 주관적인 기준에 따라서 평가하기 때문에 동일한 분자 오비탈에 대해서도 평가자에 따라 매우 다른 평가가 나올 수 있기 때문이다.
따라서, 본 발명자는 이와 같은 정성적인 평가 방법이 가지는 한계를 극복하고 분자 오비탈 정보를 체계적으로 정확하게 평가하여 이를 이용하기 위한 분자 오비탈을 정량적으로 평가할 수 있는 여러 방법을 개발하게 되었다.
예를 들어, 전체 분자 구조에서 분자 오비탈이 분포할 수 있는 패턴의 수는 거의 무한대인데, 이와 같이 복잡한 패턴을 가지고 분포하는 분자 오비탈을 직관적이고 정확하게 평가하기 위해 블록 개념을 도입해 이를 이용하는 새로운 방법인 AC2B(Assembly of Consecutive Building Block)법을 개발했다.
상기 AC2B법은 (1) 분자의 전체 분자 구조의 각각의 세부 영역을 블록으로 만든 후 전체 분자 구조를 블록의 조합으로 생성한 후, (2) 각각의 블록에 연관된 분자 오비탈이 전체 오비탈의 합에 대해 차지하는 비율을 계산해 이를 바탕으로 순차적으로 블록을 재배열하여 얻은 재배열된 순차적 블록 스펙트럼을 얻음으로써 분자 오비탈 특성을 해석하는 방법을 의미한다. 재배열된 순차적 블록 스펙트럼에서는 가장 앞에 위치한 블록에 가장 많은 분자 오비탈이 분포하고 가장 마지막에 위치한 블록에 가장 적은 분자 오비탈이 분포한다. 새롭게 블록이란 개념을 도입한 AC2B법을 이용하여 복잡한 분자 오비탈의 패턴을 단순화시켜 분자 오비탈 특성을 직관적이고 명확하게 평가할 수 있었다.
하지만 AC2B를 통해 계산되는 재배열된 블록 스펙트럼은 단일 분자 오비탈 특성을 직관적으로 나타내주지만 서로 다른 분자 오비탈 특성을 비교하는데 직접적으로 사용하기 어렵다. 이는 상기 블록 스펙트럼은 각 블록이 차지하는 위치에 따른 순차적인 블록의 배열이기 때문에 이를 정량적인 비교에 바로 이용할 수 없기 때문이다.
따라서, 본 발명자는 분자 구조 내에서 복잡하고 다양한 패턴으로 분포하고 있는 분자 오비탈 특성을 블록이라는 새로운 개념을 도입해 계산된 블록 스펙트럼을 이용해 직관적이고 명확하게 평가하고, 블록 스펙트럼을 이용해 서로 다른 분자 오비탈에 대한 유사성을 정량적으로 정확하게 비교할 수 있는 방법인 2BS-score법을 개발하였다.
2BS-score법은 서로 다른 분자 오비탈에 대해 계산된 블록 스펙트럼의 서열에 대한 동일성 평가를 통해 분자 오비탈의 유사성을 정량적으로 평가하는 평가법이다. 구체적으로 설명하면, 2BS-score법은 비교하려는 분자 오비탈에 대한 블록 스펙트럼에서 같은 서열에서의 동일성을 3단계에 걸쳐 비교하여 전체 분자 오비탈 유사도를 정량적으로 측정하는 방법이다. 비교하려는 2개의 분자 오비탈 블록 스펙트럼의 서열이 정확하게 동일한 경우 2BS-score는 100%이고, 서열 차이가 클수록 점점 작은 값을 가져 분자 오비탈 유사성이 좋지 않다는 것을 나타낸다. 본 발명자는 상기 2BS-score법에 따라 분자 오비탈 유사성을 정량적으로 정확하게 나타낼 수 있음을 확인하였다.
그러나, 상기 2BS-score법을 통해 분자 오비탈 유사성을 정량적으로 평가할 수 있지만 어느 정도의 편차를 가지는 유사성인지 여부는 알 수 없었다. 예를 들면 서로 다른 분자 오비탈 pair인 A0-A1과 A2-A3의 유사성을 2BS-score로 평가한다고 했을 때, 계산 결과 두 경우 모두 동일한 2BS-score를 가지면 동일한 유사성을 가진다고 판정되는데, 실제 A0-A1과 A2-A3사이의 유사성은 정확하게 동일하다고 볼 수 없다. 왜냐하면 2BS-score값이 동일하다고 해도 A0-A1과 A2-A3가 가지는 유사성은 편차가 존재할 수 있기 때문이다.
따라서, 이러한 유사성 편차를 정확하게 측정할 수 있다면 분자 오비탈 유사성을 보다 정확하게 평가할 수 있어 유사성 특성을 이해하는데 매우 유용할 것으로 기대되었다. 본 발명자는 이러한 필요성에 따라 2BS-score법에 따라 계산된 분자 오비탈 유사성에 대한 정량적 평가를 통해 분자 오비탈 유사성에 대한 새로운 특성, 즉 분자 오비탈 유사성에 대한 편차(similarity deviation)를 통계적으로 평가하는 방법을 개발하게 되었다.
상기 종래 기술의 문제점을 해결하기 위해 본 발명은 분자 오비탈의 유사성 편차를 정량적으로 정확하게 평가할 수 있는 새로운 방법을 제공하는 것을 목적으로 한다.
상기 목적을 달성하기 위하여, 본 발명은 a) 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 하기 i) 내지 v) 단계의 방법에 의해 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계: i) 양자역학 계산법을 이용하여 상기 분자 오비탈 분포를 계산하는 단계, ii) 상기 분자의 분자 구조 내의 분자 중심에서 방사방향(radial direction)으로 N개의 블록을 만드는 단계, iii) 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))(k는 블록의 번호이고 1 내지 N의 자연수이다)을 계산하는 단계, iv) 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록을 순차적으로 재배열하여 블록 스펙트럼을 얻는 단계, 및 v) 상기 ii)~iv) 단계를 반복하여 대상 분자 오비탈 별로 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계; b) 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m)(m는 MBS의 번호로서 1 내지 Nx 이다) 값을 계산하는 단계; 및 c) 상기 Nx개의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 단계를 포함하는 분자 오비탈 유사성 편차 평가 방법을 제공한다.
본 발명은 분자 오비탈 유사성에 포함되어 있는 편차를 정량적으로 정확하게 평가함으로써 분자 오비탈 유사성에 관련된 특성을 포괄적으로 평가하고 이해하게 할 수 있어, OLED(organic light-emitting diode) 또는 태양전지 등과 같은 분야의 신규 물질 개발에 유용하게 이용될 수 있다.
도 1은 본 발명의 분자 오비탈 유사성 편차 평가 과정을 나타낸 Flow chart 차트이다.
도 2는 서로 다른 분자 오비탈을 갖는 NPB(N,N'-Di[(1-naphthyl)-N,N'-diphenyl]-1,1'-(biphenyl)-4,4'-diamine)에 대해 서로 다른 전자 상태에서 양자역학 방법으로 계산된 분자 오비탈 T1과 T2에 대해 각각 SSDE-MO법을 적용하여 얻은 편차가 ±8.6인 분자 오비탈 pair T1과 T2를 그림으로 나타낸 것이다.
도 3은 서로 다른 분자 오비탈을 갖는 NPB(N,N'-Di[(1-naphthyl)-N,N'-diphenyl]-1,1'-(biphenyl)-4,4'-diamine)에 대해 서로 다른 전자 상태에서 양자역학 방법으로 계산된 분자 오비탈 T1과 T3에 대해 각각 SSDE-MO법을 적용하여 얻은 편차가 ±24.9인 분자 오비탈 pair T1과 T3를 그림으로 나타낸 것이다.
이하, 본 발명을 상세하게 설명한다.
본 발명의 분자 오비탈 유사성 편차 평가 방법은 a) 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 하기 i) 내지 v) 단계의 방법에 의해 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계: i) 양자역학 계산법을 이용하여 상기 분자 오비탈 분포를 계산하는 단계, ii) 상기 분자의 분자 구조 내의 분자 중심에서 방사방향(radial direction)으로 N개의 블록을 만드는 단계, iii) 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))(k는 블록의 번호이고 1 내지 N의 자연수이다)을 계산하는 단계, iv) 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록을 순차적으로 재배열하여 블록 스펙트럼을 얻는 단계, 및 v) 상기 ii)~iv) 단계를 반복하여 대상 분자 오비탈 별로 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계; b) 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m)(m는 MBS의 번호로서 1 내지 Nx 이다) 값을 계산하는 단계; 및 c) 상기 각각의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 단계를 포함하는 것을 특징으로 한다.
본 발명자는 상기 분자 오비탈 유사성 편차 평가 방법을 "SDE-MO(Statistical Similarity Deviation Estimation-Molecular Orbital)" 법 이라고 명명하였다. 상기 SSDE-MO법은 통계적인 추정 방법을 통해 분자 오비탈 유사성 편차를 정량적으로 계산하는 것이다. 상기 SSDE-MO법은 서로 다른 분자 오비탈의 유사성 편차를 계산하기 위해 (1) 블록 서열 크기가 다른 N개의 MBS(Multi-level Block Spectrum)를 생성하고 (2) 생성된 MBS 서열에 대한 동일성 평가를 통해 각각의 MBS pair에 대한 유사성을 계산하고 (3) 계산된 각각의 MBS pair에 대한 유사성 데이터를 이용해 유사성 편차를 통계적으로 계산하는 과정을 수행한다. SSDE-MO법의 계산 과정을 도 1의 Flow chart로 나타냈다. 이하 SSDE-MO 법의 계산 과정을 자세히 설명한다.
본 발명은 a) 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 하기 i) 내지 v) 단계의 방법에 의해 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계: i) 양자역학 계산법을 이용하여 상기 분자 오비탈 분포를 계산하는 단계, ii) 상기 분자의 분자 구조 내의 분자 중심에서 방사방향(radial direction)으로 N개의 블록을 만드는 단계, iii) 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))(k는 블록의 번호이고 1 내지 N의 자연수이다)을 계산하는 단계, iv) 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록을 순차적으로 재배열하여 블록 스펙트럼을 얻는 단계, 및 v) 상기 ii)~iv) 단계를 반복하여 대상 분자 오비탈 별로 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계를 포함하는 것을 특징으로 한다.
상기 a) 단계는 본 발명자가 개발한 분자 오비탈 특성 해석 방법인 "AC2B(Assembly of Consecutive Building Block)" 법을 이용하는 단계이다.
분자 오비탈은 분자 내에서의 전자의 파동적(wave-like) 거동을 나타내는 수학적인 모사로 정의할 수 있다. 분자 오비탈이 존재하는 영역은 양자역학 계산을 통해 구할 수 있으며, 상기 양자역학 계산 방법으로는 양자역학을 이용한 방법이라면 특별한 제한은 없으나, 바람직하게는 물질의 분자 구조에서 계산되는 각 지점에서의 오비탈 파동 함수(orbital wave function, ψ)의 제곱인 전자 밀도(ψ2)의 분포를 통하여 계산하는 것을 사용할 수 있고, 단일지점 에너지(single point energy) 계산 또는 기하학적 최적화(geometry optimization) 계산을 이용할 수도 있다. 구체적으로, 본 발명의 발명자들은 DFT(Density Functional Theory)에 근간을 둔 ACCELYS 사에서 개발한 MATERIAL STUDIO의 DMol3를 이용하여 분자 오비탈을 계산하였고, 분자 오비탈 그림을 생성하기 위해 MATERIAL STUDIO 패키지의 VISUALIZER를 이용하였다.
본 발명에서, 상기 대상 분자의 전체 분자 구조는, 지정된 분자 전체 블록 개수(N)에 따라 분자 중심을 기준으로 생성된 N 개의 순차적인 블록의 조합으로 구성된다.
이를 위해, 분자 중심(r=0.0)을 시작점으로 해서 방사방향(radial direction)으로 분자 전체를 포함하는 크기가 가장 큰 RDM(radially discrete mesh)을 계산하고 이때의 크기를 r=1.0으로 지정한다. 상기 RDM은 분자의 중심으로부터 출발해서 방사방향(radial direction)으로 분자 구조 내의 모든 구성 원자(atom)를 포함하는 메쉬(mesh)를 나타낸다. 상기 RDM에 의한 분자 구조 계산에 있어서, 분자 내 중심(xc, yc, zc)을 구하는 방법은 하기 수학식 1-1 내지 1-3과 같다.
(수학식 1-1)
Figure PCTKR2014006364-appb-I000001
(수학식 1-2)
Figure PCTKR2014006364-appb-I000002
(수학식 1-3)
Figure PCTKR2014006364-appb-I000003
상기 수학식 1-1 내지 1-3에서 NCoord는 분자를 구성하는 원자 좌표의 총 개수를 나타낸다.
블록의 총 개수 N 값은 특별한 제한은 없으나 바람직하게는 3 내지 100의 범위를 갖는다.
상기와 같이 RDM에 의한 분자 구조 계산에 의해, 상기 대상 분자의 전체 분자 구조는 분자 중심으로부터의 거리를 기준으로 블록화될 수 있다.
상기 a)단계는 각각의 블록에 연관된 분자 오비탈 비율(BX(k))을 계산하는 단계를 포함할 수 있다.
상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))은 전체 분자 오비탈 중에서 k번째 블록에 연관되어 있는 분자 오비탈이 차지하고 있는 양을 의미한다. 상기 k는 블록의 번호이고 1 내지 N의 자연수이다. 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))은 상기 각각의 블록에 연관된 분자 오비탈(BMO(k))을 양자역학 계산을 통해 계산하고 이를 통해서 전체 분자 오비탈 합인 SUM을 계산한 다음, 전체 분자 오비탈 합에 대한 각각의 블록에 연관된 분자 오비탈의 비율(BX(k))을 계산함으로써 구할 수 있다.
또한, 상기 a) 단계는 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록의 조합을 순차적으로 재배열하여 재배열된 블록 스펙트럼을 얻는 단계를 포함할 수 있다. 상기 "재배열된 블록 스펙트럼"은 상기 a) 단계에서 얻은 순차적 블록의 조합(AC2B, Assembly of Consecutive Building Block)을 BX(k)를 기준으로 순차적으로 재배열하여 얻은 순차적 블록의 조합을 의미한다.
상기 a) 단계는 상기 ii)~v) 단계를 반복하여 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻는 단계를 포함하는 것을 특징으로 한다. 분자 오비탈 유사성 편차를 통계적으로 추정하기 위해서는 다양한 상황에서 분자 오비탈 유사성 변화에 대한 정보가 필요하다. 이와 같은 데이터를 얻기 위해 본 발명에서는 상기 AC2B법을 이용하여 비교하려는 분자 오비탈에 대해 블록 크기(=전체 블록 개수)가 다른 Nx개의 MBS를 생성한다. 이 과정을 비교하려는 서로 다른 분자 오비탈을 T1과 T2라고 하여 설명한다. AC2B에서 BS(Block Spectrum)를 생성하기 위해 기본으로 지정된 블록 크기는 5개이다. MBS를 생성하기 위해서는 MBS의 전체 개수 Nx와 그에 따른 블록 크기 집합인 N(k)를 정해야 한다. Nx는 1보다 커야 하고 계산의 효율성을 고려하고 사용 목적에 따라 다양한 값을 가질 수 있다. 예를 들어 Nx를 6으로 지정한 경우 이에 따른 블록 크기 집합은 다음과 같이 설정될 수 있다.
N={5,6,7,8,9,10}
상기 블록 크기에 따라 상기 ii)~iv) 단계를 반복하면 대상 분자 오비탈 별로 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻을 수 있다.
상기 Nx개의 MBS(multi-block spectrum)로부터 동일한 블록 크기의 Nx개의 MBS pair를 조합한다.
본 발명에서 "MBS pair"는 대상 분자 오비탈 별로 얻은 Nx개의 MBS(multi-block spectrum)로부터 조합한 동일한 크기의 블록 스펙트럼으로 이루어진 쌍을 의미한다.
상기의 예에서 서로 다른 블록 크기를 가지는 총 6 pair(쌍)의 MBS를 조합할 수 있다.
본 발명은, b) 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m)(m는 MBS의 번호로서 1 내지 Nx 이다) 값을 계산하는 단계를 포함하는 것을 특징으로 한다.
상기 b) 단계의 다단계 동일성 평가는 제1 단계 내지 제3 단계 동일성 평가의 총 3단계 평가로 상기 Nx개의 MBS pair의 유사성을 계산하는 단계일 수 있다. 모든 과정에서 동일성 평가는 MBS pair의 동일 서열에서의 비교를 바탕으로 수행된다. 아래에 Multi-step의 3단계 과정에 대해 자세히 설명했다.
<제1 단계 동일성 평가>
총 Nx개의 MBS pair에 대해 아래와 같은 서열 비교를 통해 1차적 동일성을 평가한다.
Figure PCTKR2014006364-appb-I000004
상기 MBS1{m,n}과 MBS2{m,n}은 각각 T1과 T2 분자 오비탈에 생성된 N{m}의 블록 크기를 가지는 MBS에서 n 위치에 있는 블록의 동일성(identity)을 나타낸다.
상기 SC{m,n}은 N{m}크기를 가지는 MBS pair에 대해 n 위치에서 동일성 정도를 나타낸다. 이를 통해 MBS 서열 위치에서의 일차적 동일성을 평가한다.
상기 X는 0초과 내지 0.6이하이고, Y는 0초과 내지 0.5이하이고, Z는 0 내지 0.3이하이며, X, Y 및 Z 간에는 X+Y+(N-2)×Z=1.0의 관계가 성립한다.
상기 제1 단계 동일성 평가는 상기 Nx개의 MBS pair의 각각의 블록 스펙트럼 서열을 동일한 서열 위치에서 비교하여 블록이 동일하면 SC(m, n) 값을 부여하되, n 값이 증가할수록 SC(m, n) 값을 낮게 부여하는 것을 특징으로 한다.
상기 제1 단계 동일성 평가는 m이 1인 경우 n=1 위치에서 블록이 동일하면 SC(m, n)= X, ℓ=2 위치에서 블록이 동일하면 SC(m, n)= Y, ℓ>2 위치에서 블록이 동일하면 SC(m, n)= Z인 것으로 평가되고, 상기 X는 0초과 0.6이하이고, Y는 0초과 0.5이하이고, Z는 0 내지 0.3이하이며, X, Y 및 Z 간에는 X+Y+(N-2)×Z=1의 관계가 성립하는 것을 특징으로 한다.
이를 통해 MBS 서열 위치에서의 일차적 동일성을 평가한다.
<제2 단계 동일성 평가>
MBS pair에 대한 2차적 동일성 평가를 다음과 같이 수행한다. MBS 서열에서 가장 높은 중요도를 가지는 것은 1번째와 2번째 위치에 배열되어 있는 블록이다. 그 이유는 1번째와 2번째에 분자 오비탈이 가장 많이 분포하기 때문이다. 이와 같은 점을 고려한 2차적 동일성 평가를 아래와 같이 수행한다.
Figure PCTKR2014006364-appb-I000005
상기 제2 단계 동일성 평가는 i) 상기 Nx개의 MBS pair의 각각의 블록 스펙트럼에서 제1 블록 스펙트럼의 n=1 위치와 제2 블록 스펙트럼의 n=2 위치에서의 블록이 동일하고, 제1 블록 스펙트럼의 n=2 위치와 제1 블록 스펙트럼의 n=1 위치에서의 블록이 동일하면 SC(m, n) 값을 부여하는 단계; 및 ii) 제1 블록 스펙트럼의 n=2 위치와 제2 블록 스펙트럼의 n=3 위치에서의 블록이 동일하고, 제1 블록 스펙트럼의 n=3 위치와 제2 블록 스펙트럼의 n=2 위치에서의 블록이 동일하면 SC(m, n) 값을 부여하는 단계를 포함하고, 상기 SC(m, n) 값은 상기 i) 단계에서 부여된 값이 상기 ii) 단계에서 부여된 값보다 높은 것을 특징으로 한다.
<제3 단계 동일성 평가>
Figure PCTKR2014006364-appb-I000006
MBS pair의 각각의 동일한 서열 위치에서 블록 거리(d_BL{m,n})을 계산한다. 즉, d_BL{m,n}은 N{m} 크기를 가지는 블록 스펙트럼에서 n 위치에서의 블록 사이의 교차 거리를 나타낸다. 예를 들면 N{1} 블록 크기를 가지는 MBS pair에서 MBS1{1,1} = BL5 이고 MBS2{1,1} = BL4 이면 d_BL{1,1} = 5 - 4 = 1이다. 이를 이용해 3차적 동일성 평가를 아래와 같은 알고리즘으로 수행한다.
상기 제3 단계 동일성 평가는 상기 MBS pair의 각각의 동일한 서열 위치에서 블록서열거리(d_BL{m,n}) 기준을 설정하고, 각 위치에서의 블록서열거리가 상기 블록서열거리(d_BL{ℓ}) 기준과 동일하면 SC(m, n) 값을 부여하되, 상기 SC(m, n) 값은 n 값이 증가함에 따라 감소되도록 부여되는 것을 특징으로 한다.
상기 b) 단계에서 다단계 동일성 평가를 통해 각각의 MBS pair에 대해 동일성 평가 결과를 이용해서 하기 수학식 1을 이용해 각각의 MBS pair에 대한 유사성(TSS, total similarity score)을 계산한다.
[수학식 1]
Figure PCTKR2014006364-appb-I000007
상기 식에서 m은 MBS의 번호로서 1 내지 NX 이다.
상기 TSS{m}은 블록 크기 N{m}의 MBS pair에 대한 동일성을 나타낸다. MBS pair가 서로 정확하게 동일한 경우는 TSS는 100%값을 나타내고 유사도가 떨어질수록 100%보다 작은 값을 나타낸다.
본 발명은, c) 상기 각각의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 단계를 포함하는 것을 특징으로 한다.
Nx개의 TSS에 대해서 표준 편차(Standard Deviation)를 계산한다. 상기 표준 편차는 전체 MBS의 편차이기 때문에 결국 분자 오비탈 유사성 편차를 나타내게 된다. 상기 표준 편차(Standard Deviation)는 0보다 큰 값을 나타낸다. 모든 TSS(m)가 동일한 값을 가져서 TSS(m) 값의 표준편차(Standard Deviation)가 0인 경우 분자 오비탈 유사성 편차가 없는 것이고 표준편차가 0보다 큰 값을 나타낼수록 분자 오비탈 유사성 편차가 크다는 것을 의미한다. 이와 같은 표준편차(Standard Deviation) 값을 통해 정량적으로 분자 오비탈 유사성 편차를 평가할 수 있다.
또한, 본 발명은 상기 분자 오비탈 유사성 편차 평가 방법을 이용한 분자 오비탈 유사성 편차 평가 시스템을 제공한다.
상기 분자 오비탈 유사성 편차 평가 시스템은, a) 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 i) 양자역학 계산법을 이용하여 상기 분자 오비탈 분포를 계산하는 단계, ii) 상기 분자의 분자 구조 내의 분자 중심에서 방사방향(radial direction)으로 N개의 블록을 만드는 단계, iii) 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))(k : 1 내지 N의 자연수)을 계산하는 단계, iv) 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록을 순차적으로 재배열하여 블록 스펙트럼을 얻는 단계, 및 v) 상기 ii)~iv) 단계를 반복하여 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻는 단계를 통하여 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻는 블록화 모듈; b) 상기 2개의 대상 분자 오비탈의 각각의 Nx개의 MBS(multi-block spectrum) 로부터 동일한 블록 크기의 Nx개의 MBS pair를 조합한 후, 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m) 값을 계산하여 데이터를 입력받는 데이터 입력 모듈; 및 c) 상기 각각의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 평가 모듈을 포함하는 것을 특징으로 한다.
본 발명에서 모듈(module)이란 용어는 특정한 기능이나 동작을 처리하는 하나의 단위를 의미하며, 이는 하드웨어나 소프트웨어 또는 하드웨어 및 소프트웨어의 결합으로 구현할 수 있다.
이하 본 발명을 실시예에 기초하여 더욱 상세하게 설명하지만, 하기에 개시되는 본 발명의 실시 형태는 어디까지 예시로서, 본 발명의 범위는 이들의 실시 형태에 한정되지 않는다. 본 발명의 범위는 특허청구범위에 표시되었고, 더욱이 특허 청구범위 기록과 균등한 의미 및 범위 내에서의 모든 변경을 함유하고 있다.
실시예
서로 다른 분자 오비탈을 갖는 NPB(N,N'-Di[(1-naphthyl)-N,N'-diphenyl]-1,1'-(biphenyl)-4,4'-diamine)에 대해 서로 다른 전자 상태에서 양자역학 방법으로 계산된 분자 오비탈 pair인 T1과 T2에 대해 본 발명의 SSDE-MO법을 적용해 MBS pair에 대한 TSS값을 구하였다. 그 결과 최종적인 표준 편차는 ±8.6이었다. 이는 매우 작은 값으로 T1과 T2 분자 오비탈 사이의 유사도에는 편차가 별로 없다는 것을 나타낸다. 정성적인 평가를 위해 MATERIAL STUDIO 패키지의 프로그램인 VISUALIZER를 사용하여 분자 오비탈 T1과 T2을 시각화하여 도 2에 나타냈다. 도 2에서도 확인할 수 있듯이 분자 오비탈이 동일하지는 않지만 유사성이 매우 크다는 것을 알 수 있다. 따라서 정성적으로 확인해 보아도 T1과 T2 분자 오비탈의 유사성에 대한 편차는 크지 않다는 것을 알 수 있다. 이를 통해서 본 발명의 SSDE-MO법은 분자 오비탈 유사도 편차를 정량적으로 정확하게 나타낸다는 것을 확인할 수 있다.
서로 다른 분자 오비탈을 갖는 NPB(N,N'-Di[(1-naphthyl)-N,N'-diphenyl]-1,1'-(biphenyl)-4,4'-diamine)에 대해 서로 다른 전자 상태에서 양자역학 방법으로 계산된 분자 오비탈 pair인 T1와 T3에 대해서 SSDE-MO법을 적용한 결과 표준 편차는 ±24.9이었다. 상기 T1과 T3 분자 오비탈 사이의 유사도는 상대적으로 큰 편차를 나타냈다. 이와 같은 정량적인 SSDE-MO 평가는 도 3과 같이 그림을 통한 정성적인 평가와 정확하게 일치한다.
이를 통해 본 발명의 SSDE-MO법은 분자 오비탈 유사성에 대한 편차를 정량적으로 정확하게 평가한다는 것을 확인했다.

Claims (13)

  1. a) 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 하기 i) 내지 v) 단계의 방법에 의해 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계: i) 양자역학 계산법을 이용하여 상기 분자 오비탈 분포를 계산하는 단계, ii) 상기 분자의 분자 구조 내의 분자 중심에서 방사방향(radial direction)으로 N개의 블록을 만드는 단계, iii) 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))(k는 블록의 번호이고 1 내지 N의 자연수이다)을 계산하는 단계, iv) 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록을 순차적으로 재배열하여 블록 스펙트럼을 얻는 단계, 및 v) 상기 ii)~iv) 단계를 반복하여 대상 분자 오비탈 별로 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻은 후, 동일한 블록 크기의 Nx개의 MBS pair를 조합하는 단계;
    b) 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m)(m는 MBS의 번호로서 1 내지 Nx 이다) 값을 계산하는 단계; 및
    c) 상기 Nx개의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 단계를 포함하는 분자 오비탈 유사성 편차 평가 방법.
  2. 청구항 1에 있어서, 상기 i) 단계의 양자역학 계산법은 물질의 분자 구조에서 계산되는 각 지점에서의 오비탈 파동 함수(orbital wave function, ψ)의 제곱인 전자 밀도(ψ2)의 분포를 통하여 계산하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  3. 청구항 1에 있어서, 상기 i) 단계의 양자역학 계산법은 단일지점 에너지(single point energy) 계산 또는 기하학적 최적화(geometry optimization) 계산을 이용하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  4. 청구항 1에 있어서, 상기 iii) 단계의 각각의 블록에 연관된 단계의 분자 오비탈 비율(BX(k))은 상기 각각의 블록에 연관된 분자 오비탈(BMO(k))을 양자역학 계산을 통해 계산하고 이를 통해서 전체 분자 오비탈 합인 SUM을 계산한 다음, 전체 분자 오비탈 합에 대한 각각의 블록에 연관된 분자 오비탈의 비율(BX(k))을 계산함으로써 구하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  5. 청구항 1에 있어서, 상기 a) 단계는 RDM계산방법을 이용하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  6. 청구항 1에 있어서, 상기 b) 단계의 다단계 동일성 평가는 제1 단계 내지 제3 단계 동일성 평가의 총 3단계 평가로 상기 Nx개의 MBS pair의 유사성을 계산하는 단계인 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  7. 청구항 6에 있어서, 상기 제1 단계 동일성 평가는 상기 Nx개의 MBS pair의 각각의 블록 스펙트럼 서열을 동일한 서열 위치에서 비교하여 블록이 동일하면 SC(m, n) 값(m는 MBS의 번호로서 1 내지 Nx 이며, n은 서열번호로서 1 내지 N(m)이다), n 값이 증가할수록 SC(m, n) 값을 낮게 부여하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  8. 청구항 6에 있어서, 상기 제1 단계 동일성 평가는 m이 1인 경우 n=1 위치에서 블록이 동일하면 SC(m, n)= X, n=2 위치에서 블록이 동일하면 SC(m, n)= Y, n>2 위치에서 블록이 동일하면 SC(m, n)= Z인 것으로 평가되고,
    상기 X는 0초과 0.6이하이고, Y는 0초과 0.5이하이고, Z는 0 내지 0.3이하이며, X, Y 및 Z 간에는 X+Y+(N-2)×Z=1.0의 관계가 성립하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  9. 청구항 6에 있어서, 상기 제2 단계 동일성 평가는 i) 상기 Nx개의 MBS pair의 각각의 블록 스펙트럼에서 제1 블록 스펙트럼의 n=1 위치와 제2 블록 스펙트럼의 n=2 위치에서의 블록이 동일하고, 제1 블록 스펙트럼의 n=2 위치와 제1 블록 스펙트럼의 n=1 위치에서의 블록이 동일하면 SC(m, n) 값을 부여하는 단계; 및
    ii) 제1 블록 스펙트럼의 n=2 위치와 제2 블록 스펙트럼의 n=3 위치에서의 블록이 동일하고, 제1 블록 스펙트럼의 n=3 위치와 제2 블록 스펙트럼의 n=2 위치에서의 블록이 동일하면 SC(m, n) 값을 부여하는 단계를 포함하고,
    상기 SC(m, n) 값은 상기 i) 단계에서 부여된 값이 상기 ii) 단계에서 부여된 값보다 높은 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  10. 청구항 6에 있어서, 상기 제3 단계 동일성 평가는 상기 MBS pair의 각각의 동일한 서열 위치에서 블록서열거리(d_BL{m,n}) 기준을 설정하고, 각 위치에서의 블록서열거리가 상기 블록서열거리(d_BL{ℓ}) 기준과 동일하면 SC(m, n) 값을 부여하되, 상기 SC(m, n) 값은 n 값이 증가함에 따라 감소되도록 부여되는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  11. 청구항 1에 있어서, 상기 b) 단계의 TSS(m) 값은 하기 수학식 1에 따라 계산되는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
    [수학식 1]
    Figure PCTKR2014006364-appb-I000008
    상기 식에서 m은 MBS의 번호로서 1 내지 Nx 이다.
  12. 청구항 1에 있어서, 상기 c) 단계에서 상기 표준 편차가 0인 경우 분자 오비탈 유사성 편차가 없는 것이고 0보다 큰 값을 나타낼수록 분자 오비탈 유사성 편차가 큰 것으로 분자 오비탈 유사성 편차를 평가하는 것을 특징으로 하는 분자 오비탈 유사성 편차 평가 방법.
  13. a) 분자 오비탈 유사성을 비교할 2개의 대상 분자 오비탈을 선택한 후 i) 양자역학 계산법을 이용하여 상기 분자 오비탈 분포를 계산하는 단계, ii) 상기 분자의 분자 구조 내의 분자 중심에서 방사방향(radial direction)으로 N개의 블록을 만드는 단계, iii) 상기 각각의 블록에 연관된 분자 오비탈 비율(BX(k))(k : 1 내지 N의 자연수)을 계산하는 단계, iv) 상기 분자 오비탈 비율(BX(k))의 크기를 기준으로 상기 블록을 순차적으로 재배열하여 블록 스펙트럼을 얻는 단계, 및 v) 상기 ii)~iv) 단계를 반복하여 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻는 단계를 통하여 블록 크기가 다른 Nx개의 MBS(multi-block spectrum)를 얻는 블록화 모듈;
    b) 상기 2개의 대상 분자 오비탈의 각각의 Nx개의 MBS(multi-block spectrum) 로부터 동일한 블록 크기의 Nx개의 MBS pair를 조합한 후, 상기 Nx개의 MBS pair에 대하여 다단계 동일성 평가를 수행하여 각각의 MBS pair에 대한 TSS(m) 값을 계산하여 데이터를 입력받는 데이터 입력 모듈; 및
    c) 상기 각각의 MBS pair에 대한 TSS(m)값의 표준 편차를 계산하고 이를 통해 정량적으로 분자 오비탈 유사성 편차를 평가하는 평가 모듈을 포함하는 분자 오비탈 유사성 편차 평가 시스템.
PCT/KR2014/006364 2013-07-15 2014-07-15 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템 WO2015009016A1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US14/904,651 US10102349B2 (en) 2013-07-15 2014-07-15 Molecular orbital similarity deviation evaluation method, and system using same
JP2015533002A JP5933130B2 (ja) 2013-07-15 2014-07-15 分子オービタルの類似性偏差評価方法およびこれを用いたシステム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020130082930A KR101586382B1 (ko) 2013-07-15 2013-07-15 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템
KR10-2013-0082930 2013-07-15

Publications (1)

Publication Number Publication Date
WO2015009016A1 true WO2015009016A1 (ko) 2015-01-22

Family

ID=52346405

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/KR2014/006364 WO2015009016A1 (ko) 2013-07-15 2014-07-15 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템

Country Status (4)

Country Link
US (1) US10102349B2 (ko)
JP (1) JP5933130B2 (ko)
KR (1) KR101586382B1 (ko)
WO (1) WO2015009016A1 (ko)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101586388B1 (ko) * 2013-07-18 2016-01-18 주식회사 엘지화학 분자 오비탈 분포에 대한 정량적 비교 분석 방법 및 이를 이용한 시스템
US20160378910A1 (en) * 2015-06-25 2016-12-29 Florida Institute of Technology, Inc. Molecular active center identification using tunneling barriers and/or associated measurements for sub-molecular qsar

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06332880A (ja) * 1993-05-26 1994-12-02 Nec Corp 一般原子軌道基底に対する分子軌道初期値生成装置
US20080059549A1 (en) * 2004-09-27 2008-03-06 Japan Science And Technology Agency Molecular Orbital Computing Device for Elongation Method
US7570055B1 (en) * 2006-03-21 2009-08-04 The University Of Vermont And State Agricultural College Molecular identification and electron resonance system and method
US20110253991A1 (en) * 2008-12-26 2011-10-20 Nippon Steel Chemical Co., Ltd. Organic electroluminescence device
US20120185513A1 (en) * 2002-11-07 2012-07-19 International Business Machines Corporation Method and apparatus for storing and processing molecular information

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2005029385A1 (ja) * 2003-09-22 2006-11-30 日本電気株式会社 分子シミュレーション方法及び装置
JP2006023156A (ja) * 2004-07-07 2006-01-26 Toshiba Corp 化学物質毒性の評価方法および評価装置
WO2008085804A2 (en) * 2007-01-03 2008-07-17 Blacklight Power, Inc. System and method of computing and rendering the nature of molecules, molecular ions, compounds and materials
US20100082306A1 (en) * 2008-01-02 2010-04-01 Mills Randell L System and method of computing the nature of atoms and molecules using classical physical laws
JP2011173821A (ja) 2010-02-24 2011-09-08 Sumitomo Chemical Co Ltd 化学物質の活性度の予測方法
US8706427B2 (en) * 2010-02-26 2014-04-22 The Board Of Trustees Of The Leland Stanford Junior University Method for rapidly approximating similarities
JP5998580B2 (ja) * 2012-03-29 2016-09-28 富士通株式会社 決定プログラム、決定装置、および決定方法
JP5357320B1 (ja) * 2012-11-16 2013-12-04 学校法人立教学院 フラグメントモデル作成装置、フラグメントモデル作成システム、フラグメントモデル作成方法、及び、フラグメントモデル作成プログラム
JP2016510310A (ja) * 2012-12-18 2016-04-07 ジェネンテック, インコーポレイテッド 分子の生理活性予測
KR101472417B1 (ko) * 2013-07-09 2014-12-12 주식회사 엘지화학 순차적 블록 구성을 통한 분자 오비탈 특성 해석 방법 및 이를 이용한 시스템
KR101586386B1 (ko) * 2013-07-18 2016-01-18 주식회사 엘지화학 배타적 분자 오비탈 분포를 갖는 분자 오비탈 라이브러리 및 이를 이용한 분자 오비탈 분포 영역 평가 방법 및 이를 이용한 시스템
KR101655912B1 (ko) * 2013-07-18 2016-09-08 주식회사 엘지화학 전하 상태에 따른 분자 오비탈 분포 특성의 정량적 비교 분석 방법 및 이를 이용한 시스템
KR101586388B1 (ko) * 2013-07-18 2016-01-18 주식회사 엘지화학 분자 오비탈 분포에 대한 정량적 비교 분석 방법 및 이를 이용한 시스템
US20150051882A1 (en) * 2013-08-16 2015-02-19 Technology S.G., Lp Artificially Simulating Emissions of a Chemical Compound
US20150142398A1 (en) * 2013-11-20 2015-05-21 California Institute Of Technology Methods for a multi-scale description of the electronic structure of molecular systems and materials and related applications
US20160378910A1 (en) * 2015-06-25 2016-12-29 Florida Institute of Technology, Inc. Molecular active center identification using tunneling barriers and/or associated measurements for sub-molecular qsar

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06332880A (ja) * 1993-05-26 1994-12-02 Nec Corp 一般原子軌道基底に対する分子軌道初期値生成装置
US20120185513A1 (en) * 2002-11-07 2012-07-19 International Business Machines Corporation Method and apparatus for storing and processing molecular information
US20080059549A1 (en) * 2004-09-27 2008-03-06 Japan Science And Technology Agency Molecular Orbital Computing Device for Elongation Method
US7570055B1 (en) * 2006-03-21 2009-08-04 The University Of Vermont And State Agricultural College Molecular identification and electron resonance system and method
US20110253991A1 (en) * 2008-12-26 2011-10-20 Nippon Steel Chemical Co., Ltd. Organic electroluminescence device

Also Published As

Publication number Publication date
KR101586382B1 (ko) 2016-01-18
JP5933130B2 (ja) 2016-06-08
JP2016502158A (ja) 2016-01-21
KR20150008677A (ko) 2015-01-23
US10102349B2 (en) 2018-10-16
US20160171182A1 (en) 2016-06-16

Similar Documents

Publication Publication Date Title
Alneberg et al. CONCOCT: clustering contigs on coverage and composition
US20190139624A1 (en) Identifying ancestral relationships using a continuous stream of input
WO2014069769A1 (ko) 리드 전체를 고려한 염기 서열 정렬 시스템 및 방법
Bazin et al. panRGP: a pangenome-based method to predict genomic islands and explore their diversity
Wang et al. MBBC: an efficient approach for metagenomic binning based on clustering
Schumer et al. Versatile simulations of admixture and accurate local ancestry inference with mixnmatch and ancestryinfer
WO2015009048A1 (ko) 분자 오비탈 분포에 대한 정량적 비교 분석 방법 및 이를 이용한 시스템
WO2015009016A1 (ko) 분자 오비탈 유사성 편차 평가 방법 및 이를 이용한 시스템
WO2015005668A1 (ko) 순차적 블록 구성을 통한 분자 오비탈 특성 해석 방법 및 이를 이용한 시스템
Kokosiński et al. On sum coloring of graphs with parallel genetic algorithms
WO2015009049A1 (ko) 전하 상태에 따른 분자 오비탈 분포 특성의 정량적 비교 분석 방법 및 이를 이용한 시스템
WO2012144683A1 (ko) 유망 기술 라이프 사이클을 이용한 유망 단계 판단 장치 및 방법
Gu et al. An evolutionary algorithm based on decomposition for multimodal optimization problems
Chuang et al. SNP-SNP interaction using gauss chaotic map particle swarm optimization to detect susceptibility to breast cancer
Elishco et al. The capacity of some Pólya string models
CN110322928B (zh) Dna甲基化谱检测方法
Liu et al. Unsupervised binning of metagenomic assembled contigs using improved fuzzy C-Means method
WO2014069766A1 (ko) 염기 서열 재조합 시스템 및 방법
WO2016187898A1 (zh) 代谢物ms/ms质谱计算机仿真方法
Hwang et al. Moving clusters within a memetic algorithm for graph partitioning
Fertin et al. Finding exact and maximum occurrences of protein complexes in protein-protein interaction graphs
WO2014069770A1 (ko) 중복을 고려한 염기 서열 재조합 시스템 및 방법
Karamouz et al. Development of an algorithm for vulnerability zoning of water distribution network
CN115985400B (zh) 一种宏基因组多重比对序列重分配的方法及应用
Kazem et al. An evolutionary approach for partitioning weighted module dependency graphs

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14825624

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015533002

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14904651

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14825624

Country of ref document: EP

Kind code of ref document: A1