CN102708269A - Method for predicting inhibiting concentration of inhibitor of cytosolic phospholipase A2alpha based on support vector machine - Google Patents

Method for predicting inhibiting concentration of inhibitor of cytosolic phospholipase A2alpha based on support vector machine Download PDF

Info

Publication number
CN102708269A
CN102708269A CN201110325636XA CN201110325636A CN102708269A CN 102708269 A CN102708269 A CN 102708269A CN 201110325636X A CN201110325636X A CN 201110325636XA CN 201110325636 A CN201110325636 A CN 201110325636A CN 102708269 A CN102708269 A CN 102708269A
Authority
CN
China
Prior art keywords
descriptor
inhibitor
value
inhibitor molecules
cytosolic phospholipase
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
CN201110325636XA
Other languages
Chinese (zh)
Other versions
CN102708269B (en
Inventor
卢小泉
姬东琴
周喜斌
陈晶
史海材
刘冬
李亚亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwest Normal University
Original Assignee
Northwest Normal University
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 Northwest Normal University filed Critical Northwest Normal University
Priority to CN201110325636.XA priority Critical patent/CN102708269B/en
Publication of CN102708269A publication Critical patent/CN102708269A/en
Application granted granted Critical
Publication of CN102708269B publication Critical patent/CN102708269B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

The invention relates to a method for predicting the inhibiting concentration of an inhibitor of cytosolic phospholipase A2alpha based on a support vector machine, and belongs to the cross field of chemometrics and cheminformatics. By establishing a relational model between the molecular structure and the inhibiting concentration of one inhibitor through the support vector machine, the inhibiting concentration corresponding to the inhibitor can be predicted as long as the molecular structure of the inhibitor is known. The prediction accuracy rate of the method can achieve over 90 percent, the model is stable, the risk of inhibitor development in a later period can be reduced, and the cost of research and development can be reduced.

Description

Based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor
Technical field
The present invention relates to a kind of based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor belongs to the crossing domain of Chemical Measurement and Chemoinformatics.
Background technology
Cytosolic phospholipase A 2α is cytosolic phospholipase A 2One of three hypotypes (being respectively α, β, γ), it can limit arachidonic generation, in addition cytosolic phospholipase A 2The activation of α has caused the generation of a large amount of grease mediators, for example leukotrienes, prostaglandin, platelet activating factor etc.Owing to suppress the endochylema phospholipase A 2The activation of α has great significance, so cytosolic phospholipase A 2The exploitation of alpha inhibitor has caused drug research person's interest.The inhibition concentration experimental data test of inhibitor molecules is difficult, expense is higher; Limited obtaining of great amount of samples; And inhibition concentration is to estimate the leading indicator of inhibitor effect, therefore at the suppressant initial stage of development, accurately predicts the inhibition concentration of inhibitor molecules; Can reduce the risk of later stage suppressant exploitation, reduce R&D costs.
SVMs is based on a kind of new machine learning method of Statistical Learning Theory, and purpose is according to the estimation of given training sample to dependence between certain system's input and output, makes it make prediction as far as possible accurately to the unknown output.
Summary of the invention
The objective of the invention is to solve cytosolic phospholipase A in the prior art 2The problem of the inhibition concentration experimental data test difficulty of alpha inhibitor provides a kind of and accurately predicts cytosolic phospholipase A based on SVMs 2The method of the inhibition concentration of alpha inhibitor.
The object of the invention realizes through following technical scheme,
A kind of based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that, comprises the steps:
1) foundation of sample set: collect cytosolic phospholipase A 2The molecular structure of alpha inhibitor;
2) structure of inhibitor molecules descriptor set: input cytosolic phospholipase A 2The molecular structure of alpha inhibitor calculates the molecule descriptor value corresponding with it, and this molecule descriptor contains several components;
3) simplify the inhibitor molecules descriptor set;
4) scale again of inhibitor molecules descriptor set: the inhibitor molecules descriptor set after will simplifying is mapped to [1; + 1] interval, the mapping formula is:
Figure DEST_PATH_IMAGE001
Wherein, xBe the original value of inhibitor molecules descriptor, x Pre Be again the value after the scale, x Max With x MinThe maximal value and the minimum value of the corresponding inhibitor molecules descriptor of difference, y Max With y MinInterval maximal value of difference correspondence mappings and minimum value are promptly+1 with-1;
5) will pass through step 2) to 4) sample set after handling at random be divided into training set and test set, utilize training set data, adopt 10 to roll over cross validation method at random, the supporting vector machine model parameter is optimized;
6) with the described training set of step 5) with optimize after the SVMs parameter that obtains set up the relational model of inhibitor molecules structure and inhibition concentration;
7) with the described test set data of step 5) input step 6) relational model set up, the inhibition concentration of prediction suppressant.
Further, step 2) conformation of the molecular structure of said suppressant is in the minimum energy state.
Further, the calculating of molecule descriptor is to adopt online drug molecule descriptor computation software MODEL to accomplish step 2).
Further, the described simplification process of step 3) is:
(a) deletion and the little inhibitor molecules descriptor of inhibitor molecules structural dependence;
(b) deleting descriptor value again is 0 inhibitor molecules descriptor, the descriptor that deletion all equates for the pairing molecule descriptor value of all suppressant;
(c) use stepwise regression method that remaining inhibitor molecules descriptor is screened again.
Further, in step 5), sample set is that the ratio random division in 4:1 is training set and test set.
Further, the supporting vector machine model parameter optimisation procedure described in the step 5) is:
Capacity factor measure C is set is fixed as 100; The maximal value of the ε value variation range of ε insensitive loss function is 1, and minimum value is-1, and change step is 0.01; The maximal value of the value variation range of kernel function parameter γ is 1, and minimum value-1, change step are 0.01; Kernel function K selects the radially basic kernel function of Gauss for use;
With training set at random be divided into 10 groups; 9 groups of relational models that are used for setting up inhibitor molecules structure and inhibition concentration wherein; Remaining one group is used for verifying this model; Successively each group is carried out one-time authentication, with the mean value of resultant 10 results' in checking back accuracy rate as estimation accurately;
The value of the ε value of pairing capacity factor measure C, ε insensitive loss function, kernel function parameter γ was the optimal value of supporting vector machine model parameter when accuracy rate was the highest.
The present invention compared with prior art has the following advantages: the present invention effectively utilizes support vector machine method to set up model prediction cytosolic phospholipase A 2The inhibition concentration of alpha inhibitor; It is based on a kind of new machine learning method of Statistical Learning Theory; It is a theoretical foundation with the structural risk minimization principle; Have the Complex Nonlinear System of approaching, stronger study generalization ability and good recurrence performance and the unique advantage of the outside huge data of computing machine fast processing thereof, forecast quality and forecasting efficiency are greatly improved.At the initial stage of suppressant exploitation, through computing machine input inhibitor molecules structure,, accomplish prediction to its inhibition concentration based on support vector machine method, can reduce the risk of later stage suppressant exploitation, reduce R&D costs.Through the model of being set up is estimated, calculating training set and test set correlation coefficient r is 0.9170,0.9667, cross validation coefficient Q 2Be 0.9901,0.9943, predictablity rate of the present invention can reach more than 90%, and model is more stable.
Description of drawings
Fig. 1 utilizes the present invention to realize based on SVM prediction cytosolic phospholipase A 2The method flow diagram of the inhibition concentration of alpha inhibitor.
Fig. 2 based on support vector machine method to the endochylema phospholipase A 2The inhibition concentration of the alpha inhibitor figure that predicts the outcome.
Wherein, the IC among Fig. 2 50Be meant the 503nhibiting concentration of suppressant, pIC 50=-logIC 50
Embodiment
Below in conjunction with embodiment and accompanying drawing the present invention is further specified.
Based on SVM prediction endochylema phospholipase A 2The method of the inhibition concentration of alpha inhibitor, concrete steps are:
1, the foundation of sample set: collected 49 kinds of cytosolic phospholipase A 2Alpha inhibitor molecule and to cytosolic phospholipase A 2The inhibition concentration of α.Utilize the computer software Gaussian view 3.09 every kind of cytosolic phospholipase A that draws 2The structure of alpha inhibitor molecule is carried out conformation optimization with 7 pairs of above-mentioned molecular structures of computer software Hyperchem Release again, seeks out every kind of cytosolic phospholipase A 2The minimum energy state of alpha inhibitor molecule, the semiempirical MP3 method of selecting for use software Hyperchem to be provided when conformation is optimized.
2, the structure of inhibitor molecules descriptor set: each the inhibitor molecules sample in the sample set that step 1 is set up all will be by a descriptor vector representation; This descriptor vector has included the structural information of suppressant; This descriptor vector comprises 3778 kinds of molecule descriptors: specifically can be divided into electronics descriptor, ingredient descriptor, topological index descriptor, physicochemical property descriptor, geometrical molecular descriptor, quantum chemistry descriptor six big classes, the molecule descriptor of all samples has constituted the inhibitor molecules descriptor set.Concrete calculating is with the cytosolic phospholipase A through optimizing 2The conformation of alpha inhibitor molecule imports the software for calculation MODEL (referring to http://jing.cz3.nus.edu.sg/cgi-bin/model/model.cgi) of online drug molecule descriptor; Thereby calculate the molecule descriptor value of above-mentioned every kind of suppressant automatically; Obtain 49 * 3778 molecule descriptor values altogether, this online drug molecule descriptor computation software discloses the computing method (referring to http://jing.cz3.nus.edu.sg/model/) of described various molecule descriptor values.Wherein, Described electronics descriptor comprises that dipole moment, electric density connect index, topological electric charge index; Described ingredient descriptor comprises number of rings, atomicity, hydrogen bond number, atomic weight, hydrogen bond donor, the hydrogen bond receptor of drug molecule; Described topological index descriptor comprises that Schultz Topological Index, Gutman Topological Index, Balaban molecule connection index, Wiener chemical bond index, CHI molecule connect index, kappa shape index, Hosoya Molecular Graphs index, Zagreb Molecular Graphs index, Moreau-Broto topology auto-correlation descriptor, Moran topology auto-correlation descriptor; Described physicochemical property descriptor comprises profit partition factor, polarizability; Described geometrical molecular descriptor comprises principal moments of inertia, drug molecule volume, drug molecule surface area, and described quantum chemistry descriptor comprises high occupied orbital, lowest unoccupied molecular orbital, gross energy.
3, simplify the inhibitor molecules descriptor set: a) deletion and cytosolic phospholipase A 2The inhibitor molecules descriptor that alpha inhibitor molecular structure correlativity is little, b) deleting descriptor value again is 0 inhibitor molecules descriptor, the descriptor that deletion all equates for the pairing molecule descriptor value of all suppressant; After the deletion through a, b two steps, also surplus 1430 the drug molecule descriptors of each sample; C) use stepwise regression method that remaining molecule descriptor is screened; Each sample is finally represented by the descriptor vector that contains 5 drug molecule descriptors; The purpose of simplifying the inhibitor molecules descriptor set be deletion do not have researching value the inhibitor molecules descriptor to reduce redundance, improve the accuracy rate and the stability of prediction.
4, the scale again of inhibitor molecules descriptor set: the inhibitor molecules descriptor set after will simplifying is mapped to [1; + 1] interval, the mapping formula is:
Figure 129971DEST_PATH_IMAGE001
Wherein, xBe the original value of inhibitor molecules descriptor, x Pre Be again the value after the scale, x Max With x MinThe maximal value and the minimum value of the corresponding inhibitor molecules descriptor of difference, y Max With y MinInterval maximal value of difference correspondence mappings and minimum value are promptly+1 with-1.
5, will pass through after step 2 is handled to 4 sample set at random be divided into training set and test set; Utilize training set data; Adopt 10 to roll over cross validation method at random; The supporting vector machine model parameter is optimized, and described parameter comprises ε, the kernel function parameter γ of capacity factor measure C, ε insensitive loss function, and concrete optimizing process is following:
5.1, with 49 kinds of cytosolic phospholipase A in the sample set 2The alpha inhibitor molecule at random be divided into training set and test set, wherein training set comprises 39 samples, test set comprises 10 samples.
5.2, capacity factor measure C be set be fixed as 100; The maximal value of the ε value variation range of ε insensitive loss function is 1, and minimum value is-1, and change step is 0.01; The maximal value of the value variation range of kernel function parameter γ is 1, and minimum value-1, change step are 0.01; Kernel function K selects the radially basic kernel function of Gauss for use.
5.3, with training set at random be divided into 10 groups; Utilize 9 groups of relational models that are used for setting up inhibitor molecules structure and inhibition concentration wherein; Remaining one group is used for verifying this model; And successively each group is carried out one-time authentication, with the mean value of resultant 10 results' in checking back accuracy rate as estimation accurately.
5.4 the value of the ε value of pairing capacity factor measure C, ε insensitive loss function, kernel function parameter γ was an optimal value when accuracy rate was the highest, optimized the C, ε, the γ value that obtain and was respectively 100,0.25,0.03;
6, adopt the described training set of step 5 and optimize supported vector machine model parameter (C, ε, γ), set up the relational model of inhibitor molecules structure and inhibition concentration.
7, the relational model of step 5 described test set data substitution step 6 being set up obtains the predicted value of the inhibition concentration of corresponding inhibitor molecules.Model to being set up is estimated, and calculating training set and test set correlation coefficient r is 0.9170,0.9667, cross validation coefficient Q 2Be 0.9901,0.9943, root-mean-square error RMSE is 0.0085.

Claims (6)

1. one kind based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that, comprises the steps:
1) foundation of sample set: collect cytosolic phospholipase A 2The molecular structure of alpha inhibitor;
2) structure of inhibitor molecules descriptor set: input cytosolic phospholipase A 2The molecular structure of alpha inhibitor calculates the molecule descriptor value corresponding with it, and this molecule descriptor contains several components;
3) simplify the inhibitor molecules descriptor set;
4) scale again of inhibitor molecules descriptor set: the inhibitor molecules descriptor set after will simplifying is mapped to [1; + 1] interval, the mapping formula is:
Wherein, xBe the original value of inhibitor molecules descriptor, x Pre Be again the value after the scale, x Max With x MinThe maximal value and the minimum value of the corresponding inhibitor molecules descriptor of difference, y Max With y MinInterval maximal value of difference correspondence mappings and minimum value are promptly+1 with-1;
5) will pass through step 2) to 4) sample set after handling at random be divided into training set and test set, utilize training set data, adopt 10 to roll over cross validation method at random, the supporting vector machine model parameter is optimized;
6) with the described training set of step 5) with optimize after the SVMs parameter that obtains set up the relational model of inhibitor molecules structure and inhibition concentration;
7) with the described test set data of step 5) input step 6) relational model set up, the inhibition concentration of prediction suppressant.
2. according to claim 1 based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that: step 2) conformation of the molecular structure of said suppressant is in the minimum energy state.
3. according to claim 1 based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that: step 2) in the calculating of molecule descriptor be to adopt online drug molecule descriptor computation software MODEL to accomplish.
4. according to claim 1 based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that, the described simplification process of step 3) is:
(a) deletion and the little inhibitor molecules descriptor of inhibitor molecules structural dependence;
(b) deleting descriptor value again is 0 inhibitor molecules descriptor, the descriptor that deletion all equates for the pairing molecule descriptor value of all suppressant;
(c) use stepwise regression method that remaining inhibitor molecules descriptor is screened again.
5. according to claim 1 based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that: in step 5), sample set is that the ratio random division in 4:1 is training set and test set.
6. according to claim 1 based on SVM prediction cytosolic phospholipase A 2The method of the inhibition concentration of alpha inhibitor is characterized in that, the supporting vector machine model parameter optimisation procedure described in the step 5) is:
Capacity factor measure C is set is fixed as 100; The maximal value of the ε value variation range of ε insensitive loss function is 1, and minimum value is-1, and change step is 0.01; The maximal value of the value variation range of kernel function parameter γ is 1, and minimum value-1, change step are 0.01; Kernel function K selects the radially basic kernel function of Gauss for use;
With training set at random be divided into 10 groups; 9 groups of relational models that are used for setting up inhibitor molecules structure and inhibition concentration wherein; Remaining one group is used for verifying this model; Successively each group is carried out one-time authentication, with the mean value of resultant 10 results' in checking back accuracy rate as estimation accurately;
The value of the ε value of pairing capacity factor measure C, ε insensitive loss function, kernel function parameter γ was the optimal value of supporting vector machine model parameter when accuracy rate was the highest.
CN201110325636.XA 2011-10-24 2011-10-24 Based on SVM prediction cytosolic phospholipase A2The method of the inhibition concentration of alpha inhibitor Expired - Fee Related CN102708269B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110325636.XA CN102708269B (en) 2011-10-24 2011-10-24 Based on SVM prediction cytosolic phospholipase A2The method of the inhibition concentration of alpha inhibitor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110325636.XA CN102708269B (en) 2011-10-24 2011-10-24 Based on SVM prediction cytosolic phospholipase A2The method of the inhibition concentration of alpha inhibitor

Publications (2)

Publication Number Publication Date
CN102708269A true CN102708269A (en) 2012-10-03
CN102708269B CN102708269B (en) 2016-05-25

Family

ID=46901032

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110325636.XA Expired - Fee Related CN102708269B (en) 2011-10-24 2011-10-24 Based on SVM prediction cytosolic phospholipase A2The method of the inhibition concentration of alpha inhibitor

Country Status (1)

Country Link
CN (1) CN102708269B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103258244A (en) * 2013-04-28 2013-08-21 西北师范大学 Method for predicting inhibiting concentration of pyridazine HCV NS5B polymerase inhibitor based on particle swarm optimization support vector machine
CN104866710A (en) * 2015-05-08 2015-08-26 西北师范大学 Method for predicting inhibition concentration of cytochrome P450 enzyme CYP1A2 inhibitor by utilizing simplified partial least squares
CN109524064A (en) * 2018-11-12 2019-03-26 云南省烟草农业科学研究院 A kind of virtual screening method of polyphenol oxidase enzyme inhibitor
CN109927675A (en) * 2019-04-09 2019-06-25 深圳创维汽车智能有限公司 A kind of rain brush control method, device, equipment and storage medium
CN114783506A (en) * 2022-03-17 2022-07-22 大连理工大学 Method for predicting half-inhibitory concentration of inhibitor on coronavirus main protease
WO2022166129A1 (en) * 2021-02-08 2022-08-11 江西煌上煌集团食品股份有限公司 Screening method for heterologous competitive antigen for use in improvement of immunodetection sensitivity

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101339180A (en) * 2008-08-14 2009-01-07 南京工业大学 Organic compound combustion and explosion characteristic prediction method based on support vector machine
CN101630346A (en) * 2009-06-26 2010-01-20 上海大学 Method based on support vector machine for on-line prediction of interaction of protein and nucleic acid

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101339180A (en) * 2008-08-14 2009-01-07 南京工业大学 Organic compound combustion and explosion characteristic prediction method based on support vector machine
CN101630346A (en) * 2009-06-26 2010-01-20 上海大学 Method based on support vector machine for on-line prediction of interaction of protein and nucleic acid

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈俏等: "支持向量机应用于大气污染物浓度预测", 《计算机技术与发展》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103258244A (en) * 2013-04-28 2013-08-21 西北师范大学 Method for predicting inhibiting concentration of pyridazine HCV NS5B polymerase inhibitor based on particle swarm optimization support vector machine
CN104866710A (en) * 2015-05-08 2015-08-26 西北师范大学 Method for predicting inhibition concentration of cytochrome P450 enzyme CYP1A2 inhibitor by utilizing simplified partial least squares
CN104866710B (en) * 2015-05-08 2017-11-10 西北师范大学 The method for predicting Cytochrome P450 1A2 inhibitor inhibition concentrations
CN109524064A (en) * 2018-11-12 2019-03-26 云南省烟草农业科学研究院 A kind of virtual screening method of polyphenol oxidase enzyme inhibitor
CN109524064B (en) * 2018-11-12 2020-10-20 云南省烟草农业科学研究院 Virtual screening method of polyphenol oxidase inhibitor
CN109927675A (en) * 2019-04-09 2019-06-25 深圳创维汽车智能有限公司 A kind of rain brush control method, device, equipment and storage medium
CN109927675B (en) * 2019-04-09 2022-02-08 深圳创维汽车智能有限公司 Windshield wiper control method, device, equipment and storage medium
WO2022166129A1 (en) * 2021-02-08 2022-08-11 江西煌上煌集团食品股份有限公司 Screening method for heterologous competitive antigen for use in improvement of immunodetection sensitivity
CN114783506A (en) * 2022-03-17 2022-07-22 大连理工大学 Method for predicting half-inhibitory concentration of inhibitor on coronavirus main protease

Also Published As

Publication number Publication date
CN102708269B (en) 2016-05-25

Similar Documents

Publication Publication Date Title
Jiang et al. Dl-traff: Survey and benchmark of deep learning models for urban traffic prediction
CN102708269A (en) Method for predicting inhibiting concentration of inhibitor of cytosolic phospholipase A2alpha based on support vector machine
Torabi et al. A hybrid machine learning approach for daily prediction of solar radiation
CN108090652A (en) The method for building the electricity transaction index system based on big data technology
Northrop et al. Efficient simulation and reformulation of lithium-ion battery models for enabling electric transportation
Tang et al. Parallelization strategies for rapid and robust evolutionary multiobjective optimization in water resources applications
Qin et al. Machine learning K-means clustering algorithm for interpolative separable density fitting to accelerate hybrid functional calculations with numerical atomic orbitals
Jardak et al. Parallel processing of data from very large-scale wireless sensor networks
Zhang et al. Study on spiral source models revealing groundwater transfusion effects on pile foundation ground heat exchangers
Molaeimanesh et al. Agglomerate modeling of cathode catalyst layer of a PEM fuel cell by the lattice Boltzmann method
Kiss et al. Efficient handling of Gaussian charge distributions: An application to polarizable molecular models
CN102542131A (en) Method for predicting medicament pharmacokinetic properties and toxicity on basis of genetic algorithm and artificial neural network
CN104112167B (en) Can power generating wind resource distribution acquisition methods
Psaltis et al. Comparing charge transport predictions for a ternary electrolyte using the Maxwell–Stefan and Nernst–Planck equations
Zhang et al. Improved particle swarm optimization-extreme learning machine modeling strategies for the accurate lithium-ion battery state of health estimation and high-adaptability remaining useful life prediction
CN104199870A (en) Method for building LS-SVM prediction model based on chaotic search
CN107807961A (en) A kind of power distribution network big data multidomain treat-ment method based on Spark computing engines
Feng et al. Molecular design and performance improvement in organic solar cells guided by high‐throughput screening and machine learning
Szarek et al. Most probable distance between the nucleus and HOMO electron: the latent meaning of atomic radius from the product of chemical hardness and polarizability
CN109447512B (en) Large power grid reliability assessment method based on uniform design
CN101839946A (en) Resistor-capacitor decoupling soft sensing method based on magnitude-phase characteristics
CN104751253B (en) Distribution power flow Forecasting Methodology based on B- spline Basis bottom developed curve cluster
Li et al. Research on the method of traffic organization and optimization based on dynamic traffic flow model
Liang et al. HSMA: An O (N) electrostatics package implemented in LAMMPS
CN116316617B (en) Multi-station intelligent fusion new energy generation power region prediction method and system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160525

Termination date: 20161024