QSPR Modeling of Gibbs Free Energy of Organic Compounds by Weighting of Nearest Neighboring Codes

We examine the encoding of chemical structure of organic compounds by Labeled Hydrogen-Filled Graphs (LHFGs). Quantitative Structure-Property Relationships (QSPR) for a representative set of 150 organic molecules have been derived by means of the optimization of correlation weights of local invariants of the LHFGs. We have tested as local invariants Morgan extended connectivity of zero- and first order, numbers of path of length 2 (P2) and valence shells of distance of 2 (S2) associated with each atom in the molecular structure, and the Nearest Neighboring Codes (NNC). The best statistical characteristics for the Gibbs free energy has been obtained for the NNC weighting. Statistical parameters corresponding to this model are the following n = 100, r² = 0.9974, s = 5.136 kJ/mol, F = 38319 (training set); n = 50, r² = 0.9990, s = 3.405 kJ/mol, F = 48717 (test set). Some possible further developments are pointed out.     

Observations of trends for singular values and eigenvectors from library searching ultraviolet-visible spectra with spectral subtraction

A promising procedure for computerized library searching and identification of ultraviolet-visible spectra for one- and two-component mixtures is evaluated. The procedure is based on singular value decomposition and utilizes band position and shape. Searching more than one library spectrum at a time is possible. Library spectra which have similar spectra (collinear) with the sample and each other are indicated by two measures. First, a large condition index must be obtained, and second, two or more large variance-decomposition proportions in the same row need to be associated with the large condition index. The searching procedure has a significant degree of differentiation between the actual sample and potential candidates and is compared with the dot product, Euclidean distance and the correlation coefficient.     

A proposal for the molecular basis of μ and δ opiate receptor differentiation based on modeling of two types of cyclic enkephalins and a narcotic alkaloid

Summary The molecular basis underlying the divergent receptor selectivity of two cyclic opioid peptides Tyr-c[N -d-Orn²-Gly-Phe-Leu-] (c-ORN) and [d-Pen², l-Cys⁵]-enkephalinamide (c-PEN) was investigated using a molecular modeling approach. Ring closure and conformational searching procedures were used to determine low-energy cyclic backbone conformers. Following reinsertion of amino acid side chains, the narcotic alkaloid 7-[(1R)-1-methyl-1-hydroxy-3-phenylpropyl]-6,14-endoethenotetrahydro oripavine (PEO) was used as a flexible template for bimolecular superpositions with each of the determined peptide ring conformers using the coplanarity and cocentricity of the phenolic rings as the minimum constraint. A vector space of PEO, accounting for all possible orientations for the C₂₁-aromatic ring of PEO served as a geometrical locus for the aromatic ring of the Phe⁴ residue in the opioid peptides. Although a vast number of polypeptide conformations satisfied the criteria of the opiate pharmacophore, they could be grouped into three classes differing in magnitude and sign of the torsional angle values of the tyrosyl side chain. Only class III conformers for both c-ORN and c-PEN, having tyramine dihedral angles ₁ =–150° ± 30° and ₂=–155° ± 20°, had significant structural and conformational properties that were mutually compatible while respecting the PEO vector space. Comparison of these properties in the context of the divergent receptor selectivity of the studied opioid peptides suggests that the increased distortion of the peptide backbone in the closure region of c-PEN together with the pendant ,-dimethyl group, combine to generate a steric volume which is absent in c-ORN and that may be incompatible with a restrictive topography of the receptor. The nature and stereo-chemistry of substituents adjacent to the closure region of the peptides could also modulate receptor selection by interacting with a charged () or neutral () subsite.     

Grouped coordinate minimization using Newton's method for inexact minimization in one vector coordinate

Letf(x,y) be a function of the vector variablesx R ⁿ andy R ^m. The grouped (variable) coordinate minimization (GCM) method for minimizingf consists of alternating exact minimizations in either of the two vector variables, while holding the other fixed at the most recent value. This scheme is known to be locally,q-linearly convergent, and is most useful in certain types of statistical and pattern recognition problems where the necessary coordinate minimizers are available explicitly. In some important cases, the exact minimizer in one of the vector variables is not explicitly available, so that an iterative technique such as Newton's method must be employed. The main result proved here shows that a single iteration of Newton's method solves the coordinate minimization problem sufficiently well to preserve the overall rate of convergence of the GCM sequence.The authors are indebted to Professor R. A. Tapia for his help in improving this paper.     

A long-time tail for random walk in random scenery

Consider a simple random walk on ^d whose sites are colored black or white independently with probabilityq, resp. 1–q. Walk and coloring are independent. Letn _k be the number of steps by the walk between itskth and (k+1) th visits to a black site (i.e., the length of itskth white run), and let _k =E(n _k )–q ^–1. Our main result is a proof that (*) lim _k k ^d/2 _k = (1 –q)q ^{d/2 – 2}(d/2) ^d/2. Since it is known thatq ^{– 1} _k =E(n ₁ n _{k + 1} B) –E(n ₁ B)E(n _{k + 1} B), withB the event that the origin is black, (*) exhibits a long-time tail in the run length autocorrelation function. Numerical calculations of _k (1k100) ind=1, 2, and 3 show that there is an oscillatory behavior of _k for smallk. This damps exponentially fast, following which the power law sets in fairly rapidly. We prove that if the coloring is not independent, but is convex in the sense of FKG, then the decay of _k cannot be faster than (*).     

Random walk in an inhomogeneous medium with local impurities

In spite of Sinai's result that the decay of the velocity autocorrelation function for a random walk on ^d (d=2) can drastically change if local impurities are present, it is shown that local impurities can not abolish weak convergence to the Brownian motion if d2.     

On the convergence of the LJ search method

The convergence of the Luus-Jaakola search method for unconstrained optimization problems is established.Notation E _n Euclideann-space - f Gradient off(x) - ² f Hessian matrix - (·) ^T Transpose of (·) - I Index set {1, 2, ...,n} - [x _i1 ^*(j) ] Point around which search is made in the (j + 1)th iteration, i.e., [x _1l ^*(j) ,x _2l ^*(j) ,...,x _n1 ^*(j) ] - r _i ⁽ⁱ⁾ Range ofx _il ^*(i) in the (j + 1)th iteration - l ₁ min_i {r _i ⁽⁰⁾ } - l ₂ min_i {r _i ⁽⁰⁾ } - A _j Region of search in thejth iteration, i.e., {x E _n:x _il ^*(j-1) –0.5r _i ^(j-1) x _ix _il ^*(j-1) +0.5r _i ^(j-1) ,i I} - S _j Closed sphere with center origin and radius _j - Reduction factor in each iteration - 1– - (·) Gamma function Many discussions with Dr. S. N. Iyer, Professor of Electrical Engineering, College of Engineering, Trivandrum, India, are gratefully acknowledged. The author has great pleasure to thank Dr. K. Surendran, Professor, Department of Electrical Engineering, P.S.G. College of Technology, Coimbatore, India, for suggesting this work.     

Properties of the random search in global optimization

From theorems which we prove about the behavior of gaps in a set ofN uniformly random points on the interval [0, 1], we determine properties of the random search procedure in one-dimensional global optimization. In particular, we show that the uniform grid search is better than the random search when the optimum is chosen using the deterministic strategy, that a significant proportion of large gaps are contained in the uniformly random search, and that the error in the determination of the point at which the optimum occurs, assuming that it is unique, will on the average be twice as large using the uniformly random search compared with the uniform grid. In addition, some of the properties of the largest gap are verified numerically, and some extensions to higher dimensions are discussed. The latter show that not all of the conclusions derived concerning the inadequacies of the one-dimensional random search extend to higher dimensions, and thaton average the random search is better than the uniform grid for dimensions greater than 6.This paper is based on work started in the Statistics Department of Princeton University when the first author was visiting as a Research Associate. Part of this research was supported by the Office of Naval Research, Contract No. 0014-67-A-0151-0017, and by the US Army Research Office—Durham, Contract No. DA-31-124-ARO-D-215.2 The authors wish to thank B. Omodei for his careful work in preparing the programs for the results of Figs. 1–2 and Table 1. The computations were performed on the IBM 360/50 of the Australian National University's Computer Centre. Thanks are also due to R. Miles for suggestions regarding the extension of the results to multidimensional regions, and to P. A. P. Moran and R. Brent for suggestions regarding the evaluation of the integral ₀ ¹ ... ^0/1 (x ₁ ² + ... +x _n ^/2 )^1/2 dx ₁ ...dx _n.     

变系数计量经济学联立模型的局部线性工具向量估计及其性质

联立方程计量经济学模型在经济政策制定、经济结构分析和经济预测方面起着重要作用。本文首次提出了我国宏观经济的一类变系数联立模型,并建立了函数系数的局部线性工具向量估计,同时在时间点列固定设计、经济变量随机设计条件下,研究了估计量的大样本性质。与我国宏观经济经典线性联立模型相比,变系数联立模型拟合效果更优。另外它也有助于克服我国宏观经济数据不多而造成的非参数方法应用困难的现实情况。     

Millimeter-Wave Spectroscopy of Methylfuran Isomers: Local vs. Global Treatments of the Internal Rotation

Methylfurans are methylated aromatic heterocyclic volatile organic compounds and primary or secondary pollutants in the atmosphere due to their capability to form secondary organic aerosols in presence of atmospheric oxidants. There is therefore a significant interest to monitor these molecules in the gas phase. High resolution spectroscopic studies of methylated furan compounds are generally limited to pure rotational spectroscopy in the vibrational ground state. This lack of results might be explained by the difficulties arisen from the internal rotation of the methyl group inducing non-trivial patterns in the rotational spectra. In this study, we discuss the benefits to assign the mm-wave rotational-torsional spectra of methylfuran with the global approach of the BELGI-Cs code compared to local approaches such as XIAM and ERHAM. The global approach reproduces the observed rotational lines of 2-methylfuran and 3-methylfuran in the mm-wave region at the experimental accuracy for the ground vt=0 and the first torsional vt=1 states with a unique set of molecular parameters. In addition, the V3 and V6 parameters describing the internal rotation potential barrier may be determined with a high degree of accuracy with the global approach. Finally, a discussion with other heterocyclic compounds enables the study of the influence of the electronic environment on the hindered rotation of the methyl group.