Stabilities of multiply charged dimers and clusters of fullerenes

The authors find even-odd variations as functions of r (...+[C60]2(r+)([C60C70](r+)) electron-transfer collisions. This even-odd behavior is in sharp contrast to the smooth one for fullerene monomers and may be related to even-odd effects in dimer ionization energies in agreement with results from an electrostatic model. The kinetic energy releases for dimer dissociations [predominantly yielding intact fullerenes [C60]2(r+)-->C60(r1+)+C60(r2+) in the same (r1=r2) or nearby (r1=r2+/-1) charge states] are found to be low in comparison with the corresponding model results indicating that internal excitations of the separating (intact) fullerenes are important. Experimental appearance sizes for the heavier clusters of fullerenes [C60]n(r+) (n>3 and r=2-5) compare well with predictions from a new nearest-neighbor model assuming that r unit charges in [C60]n(r+) are localized to r C60 molecules such that the Coulomb energy of the system is minimized. The system is then taken to be stable if (i) two (singly) charged C60 are not nearest neighbors and (ii) the r C60(+) molecules have binding energies to their neutral nearest neighbors which are larger than the repulsive energies for the (r-1) C60(+)-C60(+) pairs. Essential ingredients in the nearest-neighbor model are cluster geometries and the present results on dimer stabilities.     

Molecular structure, conformation, and potential to internal rotation of 2,6- and 3,5-difluoronitrobenzene studied by gas-phase electron diffraction and quantum chemical calculations

3,5-Difluoronitrobenzene (3,5-DFNB) and 2,6-difluoronitrobenzene (2,6-DFNB) have been studied by gas-phase electron diffraction (GED), MP2 ab initio, and by B3LYP density functional calculations. Refinements of r h1 and r e static and r h1 dynamic GED models were carried out for both molecules. Equilibrium r e structures were determined using anharmonic vibrational corrections to the internuclear distances ( r e - r a) calculated from B3LYP/cc-pVTZ cubic force fields. 3,5-DFNB possesses a planar structure of C 2 v symmetry with the following r e values for bond lengths and bond angles: r(C-C) av = 1.378(4) A, r(C-N) = 1.489(6) A, r(N-O) = 1.217(2) A, r(C-F) = 1.347(5) A, angleC6-C1-C2 = 122.6(6) degrees , angleC1-C2-C3 = 117.3(3) degrees , angleC2-C3-C4 = 123.0(3) degrees , angleC3-C4-C5 = 116.9(6) degrees , angleC-C-N = 118.7(3) degrees , angleC-N-O = 117.3(4) degrees , angleO-N-O = 125.5(7) degrees , angleC-C-F = 118.6(7) degrees . The uncertainties in parentheses are three times the standard deviations. As in the case of nitrobenzene, the barrier to internal rotation of the nitro group in 3,5-DFNB, V 90 = 10 +/- 4 kJ/mol, is substantially lower than that predicted by quantum chemical calculations. The presence of substituents in the ortho positions force the nitro group to rotate about the C-N bond, out of the plane of the benzene ring. For 2,6-DFNB, a nonplanar structure of C 2 symmetry with a torsional angle of phi(C-N) = 53.8(14) degrees and the following r e values for structural parameters was determined by the GED analysis: r(C-C) av = 1.383(5) A, r(C-N) = 1.469(7) A, r(N-O) = 1.212(2) A, r(C-F) = 1.344(4) A, angleC6-C1-C2 = 118.7(5) degrees , angleC1-C2-C3 = 121.2(2) degrees , angleC2-C3-C4 = 119.0(2) degrees , angleC3-C4-C5 = 121.1(4) degrees , angleC-C-N = 120.6(2) degrees , angleC-N-O = 115.7(4) degrees , angleO-N-O = 128.6(7) degrees , angleC-C-F = 118.7(5) degrees . The refinement of a dynamic model led to barriers V 0 = 16.5 +/- 1.5 kJ/mol and V 90 = 2.2 +/- 0.5 kJ/mol, which are in good agreement with values predicted by B3LYP/6-311++G(d,p) and MP2/ cc-pVTZ calculations. The values of C-F bond lengths are similar in both molecules. This is in contrast to the drastic shortening of the C-F bond in the ortho position in 2-fluoronitrobenzene compared to the C-F bond length in the meta and para position in 3- and 4-fluoronitrobenzene observed in an earlier GED study.     

Molecular structure and conformation of chloronitromethane as determined by gas-phase electron diffraction and theoretical calculations

The molecular structure of chloronitromethane was studied in the gas phase at a nozzle-tip temperature of 373 K. The experimental data were interpreted using a dynamic model where the molecules are undergoing torsional motion governed by a potential function: V = V2/2x(1 - cos 2tau) + V4/2x(1 - cos 4tau) with V2 = 0.81(30) and V4 = 0.12(40) kcal/mol (tau is the dihedral angle between the C-Cl and N-O bond). The conformer with a zero degree dihedral angle is the most stable conformer. Comparison with results from HF/MP2/B3LYP 6-311G(d,p) calculations were made. The important geometrical parameter values (for the eclipsed form) obtained from least-squares refinements are the following: r(C-H) = 1.061(18)A, r(C-N) = 1.509 (5)A, r(N-O) = 1.223(1)A, r(C-Cl) = 1.742(2)A, angleClCN = 115.2(7) degrees, angleO4NC = 118.9(10) degrees, angleO5NC = 114.9(16) degrees, and angleClCH 115(4) degrees.     

A comparative study of quantitative structure-activity relationship methods based on gallic acid derivatives

By using hologram quantitative structure-activity relationship (HQSAR) and comparative molecular field analysis (CoMFA) methods, the relationships between the structures of 49 gallic acid derivatives and their analgesic activity have been investigated to yield statistically reliable models with considerable predictive power. The best HQSAR model was generated using atoms, bond and connectivity as fragment distinction parameters and fragment size 5-7 from a hologram length of 307 with 3 components. High conventional r2 (r2 = 0.825) and cross-validation r2 (r2(cv) = 0.726) values were obtained. CoMFA analyses varying lattice size and location, grid spacing, probe charges and using, Tripos standard and Indicator force field were performed. The best model was developed with 4 components using sp3-hybridized carbon atom with +1.0 charge as probe, grid spacing (2 A), lattice offset (1.0, 3.0, -2.5). The CoMFA model showed a conventional correlation coefficient r2 of 0.889 and across-validation r2(cv) equals to 0.633. The robustness and predictive ability of the HQSAR and CoMFA models have been validated by means of an external test set. The results indicate that both models possess high statistical quality in the prediction of analgesic potency of novel gallic acid analogs.     

有电荷转移络合物参与的共聚体系竞聚率

从有电荷转移络合物参与的共聚体系增长基元反应出,提出一种有ＣＴＣ参与的共聚体系竞聚率的求法,从中得到４个竞聚率,自由单体和ＣＴＣ的相对活性等信息。     

An extensive ecdysteroid CoMFA

The ecdysteroid agonist activity of 71 HPLC-purified ecdysteroids was measured in the Drosophila melanogaster BII tumorous blood cell line assay. The resultant log(ED50) values, spanning almost 6 orders of magnitude, were used to construct a comparative molecular field analysis (CoMFA) model in which conformations were selected by homology to the crystal structure of ecdysone. Model A was constructed by utilization of the region- focused electrostatic indicator field (q2=0.631, r2=0.903, 5 components, 4 outliers). Model B made use of region-focused electrostatic and steric indicator fields along with MlogP (q2=0.694, r2=0.892, 5 components, 4 outliers). The model and its underlying bioassay data support a pharmacophore hypothesis in which ecdysteroid binding is understood to be due principally to the summation of localized interactions from approximately six specific loci. This is in contrast to previous structure-activity relationship hypotheses which are formulated in terms of the presence or absence of essential functional groups, without which ecdysteroid receptor affinity would be completely absent. The present CoMFA model is utilized to predict the activities of heretofore unknown ecdysteroids.     

Three-dimensional quantitative structure-activity and structure-selectivity relationships of dihydrofolate reductase inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) modelling using comparative molecular similarity indices analysis (CoMSIA) was applied to a series of 406 structurally diverse dihydrofolate reductase (DHFR) inhibitors from Pneumocystis carinii (pc) and rat liver (rl). X-ray crystal structures of three inhibitors bound to pcDHFR were used for defining the alignment rule. For pcDHFR, a QSAR model containing 6 components was selected using leave-10%-out cross-validation (n= 240, q2 = 0.65), while a 4-component model was selected for rlDHFR (n= 237, q2 = 0.63); both include steric, electrostatic and hydrophobic contributions. The models were validated using a large test set, designed to maximise its diversity and to verify the predictive accuracy of models for extrapolation. The pcDHFR model has r2 = 0.60 and mean absolute error (MAE) = 0.57 for the test set after removing 4 outliers, and the rlDHFR model has r2 = 0.60 and MAE = 0.69 after removing 4 test set outliers. In addition, classification models predicting selectivity for pcDHFR over rlDHFR were developed using soft independent modelling by class analogy (SIMCA), with a selectivity ratio of 2 (IC50,rlDHFR/ IC50,pcDHFR) used for delimiting classes. A 5-component model including steric and electrostatic contributions has cross-validated and test set classification rates of 0.67 and 0.68 for selective inhibitors, and 0.85 and 0.72 for unselective inhibitors. The predictive accuracy of models, together with the identification of important contributions in QSAR and classification models, offer the possibility of designing potent selective inhibitors and estimating their activity prior to synthesis.     

Study of the thymine molecule: equilibrium structure from joint analysis of gas-phase electron diffraction and microwave data and assignment of vibrational spectra using results of ab initio calculations

Thymine is one of the nucleobases which forms the nucleic acid (NA) base pair with adenine in DNA. The study of molecular structure and dynamics of nucleobases can help to understand and explain some processes in biological systems and therefore it is of interest. Because the scattered intensities on the C, N, and O atoms as well as some bond lengths in thymine are close to each other the structural problem cannot been solved by the gas phase electron diffraction (GED) method alone. Therefore the rotational constants from microvawe (MW) studies and differences in the groups of N-C, C=O, N-H, and C-H bond lengths from MP2 (full)/cc-pVQZ calculations were used as supplementary data. The analysis of GED data was based on the C(s) molecular symmetry according to results of the structure optimizations at the MP2 (full) level using 6-311G (d,p), cc-pVTZ, and cc-pVQZ basis sets confirmed by vibrational frequency calculations with 6-311G (d,p) and cc-pVTZ basis sets. Mean-square amplitudes as well as harmonic and anharmonic vibrational corrections to the internuclear distances (r(e)-r(a)) and to the rotational constants (B(e)(k)-B(0)(k), where k = A, B, C) were calculated from the quadratic (MP2 (full)/cc-pVTZ) and cubic (MP2 (full)/6-311G (d,p)) force constants (the latter were used only for anharmonic corrections). The harmonic force field was scaled using published IR and Raman spectra of the parent and N1,N3-dideuterated species, which were for the first time completely assigned in the present work. The main equilibrium structural parameters of the thymine molecule determined from GED data supplemented by MW rotational constants and results of MP2 calculations are the following (bond lengths in Angstroms and bond angles in degrees with 3sigma in parentheses): r(e) (C5=C6) = 1.344 (16), r(e) (C5-C9) = 1.487 (8), r(e) (N1-C6) = 1.372 (3), r(e) (N1-C2) = 1.377 (3), r(e) (C2-N3) = 1.378 (3), r(e) (N3-C4) = 1.395 (3), r(e) (C2=O7) = 1.210 (1), r(e) (C4=O8) = 1.215 (1), angle e (N1-C6=C5) = 123.1 (5), angle e (C2-N1-C6) = 123.7 (5), angle e (N1-C2-N3) = 112.8 (5), angle e (C2-N3-C4) = 128.0 (5), angle e (N3-C4-C5) = 114.8 (5), angle e (C6=C5-C9) = 124.4 (9). The experimental structural parameters are in good agreement with those from MP2 (full) calculations with use of cc-pVTZ and cc-pVQZ basis sets.     

QSAR prediction of antagonistic activity of PCBs towards human PXR by using heuristic method and best subset modeling

Polychlorinated biphenyls(PCBs) can antagonize human pregnane X receptor(hPXR) activation.Such chemicals could pose a serious threat to the reproductive and developmental ability of humans.The quantitative structure activity relationship(QSAR) provides a promising method for the estimation of PCBs’ antagonistic activity.In this investigation,a QSAR model was developed by using heuristic method and best subset modeling(r2 = 0.873,q2LOO=0.742).The built model was validated externally by splitting the original data set into training and prediction sets.The results of the model derived are as follows:r2 = 0.907,q2LOO=0.709,r2pred=0.676,suggesting developed QSAR model had good robustness and predictive ability.The applicability domain(AD) of the model was assessed by Williams plot.The antagonistic activity(?logKi) of 108 PCBs,which are unavailable by experiment at present,was predicted within the applicability domain of the model.The critical structural features related to the activity of PCBs were identified.     

Quantitative determination of melatonin in human plasma and cerebrospinal fluid with high-performance liquid chromatography and fluorescence detection

A validated new and precise reversed-phase high-performance liquid chromatographic method for the determination of melatonin in human plasma and cerebrospinal fluid, with 5-fluorotryptamine as internal standard, is described. Liquid-liquid extraction with dichloromethane was performed under alkaline conditions. After evaporation of the organic solvent, the extract was dissolved in eluent and chromatographed on a base-deactivated octadecyl column, using an eluent composed of 650 mL potassium dihydrogenphosphate solution (0.07 mol/L water), adjusted to a pH of 3.0 with a 43% phosphoric acid solution, mixed with 350 mL methanol. Fluorescence detection at an excitation wavelength of 224 nm and an emission wavelength of 348 nm was used for quantitation. Melatonin and 5-fluorotryptamine chromatographed with retention times of 5.3 and 9. 3 min, respectively. Mean recoveries of 96% (n = 10) and 95% (n = 5) were found for melatonin in plasma and cerebrospinal fluid respectively. 5-Fluorotryptamine was found to have a mean recovery of 90% (n = 10) and 82% (n = 5) in plasma and cerebrospinal fluid, respectively. The repeatability coefficients of variation for both melatonin and 5-fluorotryptamine in plasma were 4-5% [five different samples (r = 5) on two consecutive days (n = 2)], with reproducibility coefficients of 1.6-7% (n = 2, r = 5) and 0.9-4% (n = 2, r = 5) for melatonin and internal standard, respectively. In cerebrospinal fluid the repeatability coefficient of variation of the extraction procedure was 5% (n = 1, r = 5) for melatonin and 7% (n = 1, r = 5) for 5-fluorotryptamine. The correlation coefficients of the calibration curves were 0.9998 (n = 2) in plasma at a concentration range of 0.108-25.9 ng/mL and 0.9994 (n = 2) at a concentration range of 0.108-25.9 ng/mL in cerebrospinal fluid. The limit of detection was determined at 8 pg/mL which enables to measure melatonin concentrations at physiological concentrations reached during daytime.     

Microwave spectra, ab initio calculations, and r0 structural parameters for (methylamino)thiophosphoryl difluoride

The microwave spectra of (methylamino)thiophosphoryl difluoride, CH(3)NHP(=S)F(2), and two deuterated species, CH(3)NDP(=S)F(2) and CD(3)NHP(=S)F(2), have been investigated in the region from 26.5 to 39.0 GHz. The rotational constants of the ground vibrational state have been determined and have been shown to be only consistent with the trans conformer (CH(3) group antiperiplanar to the P=S bond) with C(s) symmetry. The a-type R branch transitions have been assigned for the trans conformer for the three isotopomers on the basis of the rigid rotor model. Near-trans and near-cis forms without molecular planes of symmetry are predicted by all ab initio calculations with the near-trans form being more stable. However, the double-well potentials governing the interchange between the two enantiomeric near-trans as well as the two near-cis forms are too shallow to accommodate the zero-point energies of the nu(24) asymmetric torsion. Thus, the trans conformation with C(s) symmetry may be more accurate in explaining the microwave experimental data. The "adjusted" r(0) structural parameters have been obtained by systematically adjusting the ab initio MP2(full)/6-311+G(d,p) structure of the trans conformer with C(s) symmetry to fit the microwave rotational constants. The determined heavy atom distances are r(C-N) = 1.459(5), r(P-N) = 1.621(5), r(P=S) = 1.879(5), and r(P-F) = 1.550(5) A, and the heavy atom angles are angleCNP = 124.7(5) degrees , angleNPS = 118.3(5) degrees , angleNPF = 103.2(5) degrees , angleFPS = 117.0(5) degrees , and angleFPF = 94.6(5) degrees . The adjusted r(0) parameters have also been obtained for aminodifluorophosphine, H(2)NPF(2), with a slightly pyramidal -PNH(2) moiety. The results indicate that the previously reported short distance of 0.981(5) A for the N-H(o)(outer) bond from the microwave study is too short, and the adjusted r(0) value of 1.007(3) A is obtained from the combined data. Adjusted r(0) parameters are also reported for (dimethylamino)difluorophosphine, (CH(3))(2)NPF(2), with C(s) symmetry with the PNC(2) portion of the molecule being planar. The previously reported C-H distances from the electron diffraction study are too long, and the anglePNC(i) and angleC(o)NC(i) angles are also found to be in error. These results provide a reasonable explanation why the microwave and electron diffraction results differ for the structures of these latter two molecules.     

Reversible intramolecular hydrogen transfer between cysteine thiyl radicals and glycine and alanine in model peptides: absolute rate constants derived from pulse radiolysis and laser flash photolysis

The intramolecular reaction of cysteine thiyl radicals with peptide and protein alphaC-H bonds represents a potential mechanism for irreversible protein oxidation. Here, we have measured absolute rate constants for these reversible hydrogen transfer reactions by means of pulse radiolysis and laser flash photolysis of model peptides. For N-Ac-CysGly6 and N-Ac-CysGly2AspGly3, Cys thiyl radicals abstract hydrogen atoms from Gly with k(f) = (1.0-1.1 x 10(5) s(-1), generating carbon-centered radicals, while the reverse reaction proceeds with k(r) = (8.0-8.9) x 10(5) s(-1). The forward reaction shows a normal kinetic isotope effect of k(H)/k(D) = 6.9, while the reverse reaction shows a significantly higher normal kinetic isotope effect of 17.6, suggesting a contribution of tunneling. For N-Ac-CysAla2AspAla3, cysteine thiyl radicals abstract hydrogen atoms from Ala with k(f) = (0.9-1.0) x 10(4) s(-1), while the reverse reaction proceeds with k(r) = 1.0 x 10(5) s(-1). The order of reactivity, Gly > Ala, is in accord with previous studies on intermolecular reactions of thiyl radicals with these amino acids. The fact that k(f) < k(r) suggests some secondary structure of the model peptides, which prevents the adoption of extended conformations, for which calculations of homolytic bond dissociation energies would have predicted k(f) > k(r). Despite k(f) < k(r), model calculations show that intramolecular hydrogen abstraction by Cys thiyl radicals can lead to significant oxidation of other amino acids in the presence of physiologic oxygen concentrations.     

Implementation of a genetically tuned neural platform in optimizing fluorescence from receptor-ligand binding interactions on microchips

This paper describes the use of a genetically tuned neural network platform to optimize the fluorescence realized upon binding 5-carboxyfluorescein-D-Ala-D-Ala-D-Ala (5-FAM-(D-Ala)(3) ) (1) to the antibiotic teicoplanin from Actinoplanes teichomyceticus electrostatically attached to a microfluidic channel originally modified with 3-aminopropyltriethoxysilane. Here, three parameters: (i) the length of time teicoplanin was in the microchannel; (ii) the length of time 1 was in the microchannel, thereby, in equilibrium with teicoplanin, and; (iii) the amount of time buffer was flushed through the microchannel to wash out any unbound 1 remaining in the channel, are examined at a constant concentration of 1, with neural network methodology applied to optimize fluorescence. Optimal neural structure provided a best fit model, both for the training set (r(2) = 0.985) and testing set (r(2) = 0.967) data. Simulated results were experimentally validated demonstrating efficiency of the neural network approach and proved superior to the use of multiple linear regression and neural networks using standard back propagation.