Semiempirical Study of the Bergman Reaction: Towards a Computationally Efficient and Accurate Method for Modeling Enediyne Anticancer Antibiotics

Over the past 7 years, the enediyne anticancer antibiotics have been widely studied due to their DNA cleaving ability. The focus of these antibiotics, represented by kedarcidin chromophore, neocarzinostatin chromophore, calicheamicin, esperamicin A, and dynemicin A, is on the enediyne moiety contained within each of these antibiotics. In its inactive form, the moiety is benign to its environment. Upon suitable activation, the system undergoes a Bergman cycloaromatization proceeding through a 1,4-dehydrobenzene diradical intermediate. It is this diradical intermediate that is thought to cleave double-stranded dna through hydrogen atom abstraction. Semiempirical, semiempiricalci, Hartree–Fock ab initio, and mp2 electron correlation methods have been used to investigate the inactive hex-3-ene-1,5-diyne reactant, the 1,4-dehydrobenzene diradical, and a transition state structure of the Bergman reaction. Geometries calculated with different basis sets and by semiempirical methods have been used for single-point calculations using electron correlation methods. These results are compared with the best experimental and theoretical results reported in the literature. Implications of these results for computational studies of the enediyne anticancer antibiotics are discussed.

Rotational constants and structure of para-difluorobenzene determined by femtosecond Raman coherence spectroscopy: A new transient type

Femtosecond Raman rotational coherence spectroscopy (RCS) detected by degenerate four-wave mixing is a background-free method that allows to determine accurate gas-phase rotational constants of non-polar molecules. Raman RCS has so far mostly been applied to the regular coherence patterns of symmetric-top molecules, while its application to nonpolar asymmetric tops has been hampered by the large number of RCS transient types, the resulting variability of the RCS patterns, and the 10³–10⁴ times larger computational effort to simulate and fit rotational Raman RCS transients. We present the rotational Raman RCS spectra of the nonpolar asymmetric top 1,4-difluorobenzene (para-difluorobenzene, p-DFB) measured in a pulsed Ar supersonic jet and in a gas cell over delay times up to ~2.5 ns. p-DFB exhibits rotational Raman transitions with ΔJ = 0, 1, 2 and ΔK = 0, 2, leading to the observation of J −, K −, A −, and C–type transients, as well as a novel transient (S–type) that has not been characterized so far. The jet and gas cell RCS measurements were fully analyzed and yield the ground-state (v = 0) rotational constants Aₒ = 5637.68(20) MHz, Bₒ = 1428.23(37) MHz, and Cₒ = 1138.90(48) MHz (1σ uncertainties). Combining the Aₒ, Bₒ, and Cₒ constants with coupled-cluster with single-, double- and perturbatively corrected triple-excitation calculations using large basis sets allows to determine the semi-experimental equilibrium bond lengths rₑ(C₁–C₂) = 1.3849(4) Å, rₑ(C₂–C³) = 1.3917(4) Å, rₑ(C–F) = 1.3422(3) Å, and rₑ(C₂–H₂) = 1.0791(5) Å.

Gaussian operator bases for correlated fermions

We formulate a general multi-mode Gaussian operator basis for fermions, to enable a positive phase-space representation of correlated Fermi states. The Gaussian basis extends existing bosonic phase-space methods to Fermi systems and thus allows first-principles dynamical or equilibrium calculations in quantum many-body Fermi systems. We prove the completeness of the basis and derive differential forms for products with one- and two-body operators. Because the basis satisfies fermionic superselection rules, the resulting phase space involves only c-numbers, without requiring anticommuting Grassmann variables. Furthermore, because of the overcompleteness of the basis, the phase-space distribution can always be chosen positive. This has important consequences for the sign problem in fermion physics.

Robust automatic mapping algorithms in a network monitoring scenario

Automatically generating maps of a measured variable of interest can be problematic. In this work we focus on the monitoring network context where observations are collected and reported by a network of sensors, and are then transformed into interpolated maps for use in decision making. Using traditional geostatistical methods, estimating the covariance structure of data collected in an emergency situation can be difficult. Variogram determination, whether by method-of-moment estimators or by maximum likelihood, is very sensitive to extreme values. Even when a monitoring network is in a routine mode of operation, sensors can sporadically malfunction and report extreme values. If this extreme data destabilises the model, causing the covariance structure of the observed data to be incorrectly estimated, the generated maps will be of little value, and the uncertainty estimates in particular will be misleading. Marchant and Lark [2007] propose a REML estimator for the covariance, which is shown to work on small data sets with a manual selection of the damping parameter in the robust likelihood. We show how this can be extended to allow treatment of large data sets together with an automated approach to all parameter estimation. The projected process kriging framework of Ingram et al. [2007] is extended to allow the use of robust likelihood functions, including the two component Gaussian and the Huber function. We show how our algorithm is further refined to reduce the computational complexity while at the same time minimising any loss of information. To show the benefits of this method, we use data collected from radiation monitoring networks across Europe. We compare our results to those obtained from traditional kriging methodologies and include comparisons with Box-Cox transformations of the data. We discuss the issue of whether to treat or ignore extreme values, making the distinction between the robust methods which ignore outliers and transformation methods which treat them as part of the (transformed) process. Using a case study, based on an extreme radiological events over a large area, we show how radiation data collected from monitoring networks can be analysed automatically and then used to generate reliable maps to inform decision making. We show the limitations of the methods and discuss potential extensions to remedy these.

Modern Computational Chemistry Methods for Prediction of Ground- and Excited-State Properties in Open-Shell Systems

The main goal of the research presented in this work is to provide some important insights about computational modeling of open-shell species. Such projects are: the investigation of the size-extensivity error in Equation-of-Motion Coupled Cluster methods, the analysis of the Long-Range corrected scheme in predicting UV-Vis spectra of Cu(II) complexes with the 4-imidazole acetate and its ethylated derivative, and the exploration of the importance of choosing a proper basis set for the description of systems such as the lithium monoxide anion. The most significant findings of this research are: (i) The contribution of the left operator to the size-extensivity error of the CR-EOMCC(2,3) approach, (ii) The cause of d-d shifts when varying the range-separation parameter and the amount of the exact exchange arising from the imbalanced treatment of localized vs. delocalized orbitals via the "tuned" CAM-B3LYP* functional, (iii) The proper acidity trend of the first-row hydrides and their lithiated analogs that may be reversed if the basis sets are not correctly selected.

Multicomponent kinetic spectrophotometric determination of pefloxacin and norfloxacin in pharmaceutical preparations and human plasma samples with the aid of chemometrics

A spectrophotometric method for the simultaneous determination of the important pharmaceuticals, pefloxacin and its structurally similar metabolite, norfloxacin, is described for the first time. The analysis is based on the monitoring of a kinetic spectrophotometric reaction of the two analytes with potassium permanganate as the oxidant. The measurement of the reaction process followed the absorbance decrease of potassium permanganate at 526 nm, and the accompanying increase of the product, potassium manganate, at 608 nm. It was essential to use multivariate calibrations to overcome severe spectral overlaps and similarities in reaction kinetics. Calibration curves for the individual analytes showed linear relationships over the concentration ranges of 1.0–11.5 mg L−1 at 526 and 608 nm for pefloxacin, and 0.15–1.8 mg L−1 at 526 and 608 nm for norfloxacin. Various multivariate calibration models were applied, at the two analytical wavelengths, for the simultaneous prediction of the two analytes including classical least squares (CLS), principal component regression (PCR), partial least squares (PLS), radial basis function-artificial neural network (RBF-ANN) and principal component-radial basis function-artificial neural network (PC-RBF-ANN). PLS and PC-RBF-ANN calibrations with the data collected at 526 nm, were the preferred methods—%RPET not, vert, similar 5, and LODs for pefloxacin and norfloxacin of 0.36 and 0.06 mg L−1, respectively. Then, the proposed method was applied successfully for the simultaneous determination of pefloxacin and norfloxacin present in pharmaceutical and human plasma samples. The results compared well with those from the alternative analysis by HPLC.

Scanning the genome for essential hypertension loci

1. Essential hypertension occurs in people with an underlying genetic predisposition who subject themselves to adverse environmental influences. The number of genes involved is unknown, as is the extent to which each contributes to final blood pressure and the severity of the disease. 2. In the past, studies of potential candidate genes have been performed by association (case-control) analysis of unrelated individuals or linkage (pedigree or sibpair) analysis of families. These studies have resulted in several positive findings but, as one may expect, also an enormous number of negative results. 3. In order to uncover the major genetic loci for essential hypertension, it is proposed that scanning the genome systematically in 100- 200 affected sibships should prove successful. 4. This involves genotyping sets of hypertensive sibships to determine their complement of several hundred microsatellite polymorphisms. Those that are highly informative, by having a high heterozygosity, are most suitable. Also, the markers need to be spaced sufficiently evenly across the genome so as to ensure adequate coverage. 5. Tests are performed to determine increased segregation of alleles of each marker with hypertension. The analytical tools involve specialized statistical programs that can detect such differences. Non- parametric multipoint analysis is an appropriate approach. 6. In this way, loci for essential hypertension are beginning to emerge.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

Studies on the relative biopotencies and intestinal absorption of different apo-beta-carotenoids in rats and chickens.

1. The biopotencies relative to beta-carotene of several apocarotenoids, such as 8'-, 10'- and 12'-apo-beta-carotenal and methyl 8'-apo-beta-carotenoate, were investigated in rats, on a molar basis, by both curative-growth assay and liver-storage tests. 2. In the curative-growth assays, on a molar basis the biopotencies of 8'-, 10'- and 12'-apo-beta-carotenal and methyl 8'-apo-beta-carotenoate were 72, 78, 72 and 53% respectively, whereas on a weight basis the corresponding values were 93, 111, 111 and 63%, with respect to beta-carotene taken as 100%. In terms of yield of vitamin A, these values were much lower in the liver-storage tests. 3. When 8'-apo-beta-carotenal was fed, the unchanged aldehyde together with small amounts of the corresponding alcohol and larger proportions of the acid rapidly appeared in the tissues of both rats and chickens. The 8'-apocarotenol, 8'-apocarotenoic acid and its methyl ester were absorbed unchanged. The free acid disappeared most rapidly from the tissues, but its methyl ester persisted in the tissues longest. 4. On the basis of these observations it is suggested that most of an apocarotenal is oxidized to the corresponding acid, which, in turn, is mostly degraded to retinoic acid, with small proportions of it being attacked by the dioxygenase system giving retinal.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

AB-initio study of alpha-silyl and beta-silyl substituent effects on vinyl radicals

The geometries of alpha- and beta-silyl substituted vinyl radicals and of alpha,beta-disilylvinyl radical have been optimised with the STO-3G and the STO-3G* basis sets. The relative stabilities of various conformers have been determined at the UMP2/6-31G* level. The stabilisation of vinyl radicals through alpha-silyl substitution is larger than that due to corresponding alkyl groups. The presence of an alpha-silyl group also leads to a tendency towards linearisation of the vinyl radical centre and a corresponding reduction in the inversion barrier. In marked contrast, the beta-silyl effect is negligible. The geometric, conformational and energetic consequences are insignificant. Overall, the silyl substituent effect at vinyl radicals is very different from that computed earlier for the vinyl cations, but qualitatively similar to that found in carbanions.

SCF and electron correlation studies on structure, force constants and vibrational spectra of borane monoammoniate complex

Effects of basis set and electron correlation on the equilibrium geometry, force constants and vibrational spectra of BH3NH3 have been studied. A series of basis sets ranging from double zeta to triple zeta including polarization and diffuse functions have been utilized. All the SCF based calculations overestimate the dative B-N bond distance and considerable improvement occurs when the treatment for electron correlation is introduced. Detailed vibrational analysis for BH3NH3 has been carried out. The mean absolute percentage deviation of the ab initio predicted vibration frequencies of (BH3NH3)-B-11 from the experiment is about 10% for the SCF based calculations and the MP2 method shows better agreement, the overall deviation being 5-6%. The ground state effective force constants of BH3NH3 were obtained using RECOVES procedure. The RECOVES sets of force constants are found to be highly satisfactory for the prediction of the vibrational frequencies of different isotopomers of BH3NH3. The mean absolute percentage deviation of the calculated frequencies of different isotopomers from the experiment is much less than 1%. The RECOVES-MP2/augDZP set of force constants was found to be the best set among the different sets for this molecule. Theoretical infrared intensities are in fair agreement with the observed spectral features.

Ab initio studies of solvent effect on conformational equilibria and vibrational spectra of dipropionamide

Systematic ab initio molecular orbital studies of the conformational equilibria and vibrational spectra of dipropionamide using the basis sets 6-31g(d) and 6-31++G(d,p) have been carried out. The vibrational spectra of dipropionamide have been satisfactorily interpreted taking into account the agreement between the calculated frequencies, infrared and Raman band intensities and the shifts in the spectra of deuterated molecules with those observed. The previous assignments of most of the vibrational bands are well confirmed, a few bands need reassignment, however. The solvent effects were investigated by self-consistent reaction field theory using dipole and self-consistent isodensity polarized continuum model methods. The introduction of a dielectric medium has only a marginal effect on the conformational equilibria and vibrational spectra. However, the calculated changes in geometry and vibrational spectra on going from the gas phase to the solution phase are in accord with the increasing weight of the dipolar resonance structure in polar solvents. (C) 2002 Elsevier Science B.V. All rights reserved.

Thermal decomposition of haloethanols: single pulse shock tube and ab initio studies

This paper reports single pulse shock tube and ab initio studies on thermal decomposition of 2-fluoro and 2-chloroethanol at T=1000–1200 K. Both molecules have HX (X = F/Cl) and H2O molecular elimination channels. The CH3CHO formed by HX elimination is chemically active and undergoes secondary decomposition resulting in the formation of CH4, C2H6, and C2H4. A detailed kinetic simulation indicates that the formation of C2H4 could not be quantitatively explained as arising exclusively from secondary CH3CHO decomposition. Contributions from primary radical processes need to be considered to explain C2H4 quantitatively. Ab initio calculations on HX and H2O elimination reactions from the haloethanols at HF, MP2, and DFT levels with various basis sets up to 6/311++G**are reported. It is pointed out that due to strong correlations between A and Eα, comparison of these two parameters between experimental and theoretical results could be misleading.

Prospects of zero Schottky barrier height in a graphene-inserted MoS2-metal interface

A low Schottky barrier height (SBH) at source/drain contact is essential for achieving high drive current in atomic layer MoS(2-)channel-based field effect transistors. Approaches such as choosing metals with appropriate work functions and chemical doping are employed previously to improve the carrier injection from the contact electrodes to the channel and to mitigate the SBH between the MoS2 and metal. Recent experiments demonstrate significant SBH reduction when graphene layer is inserted between metal slab (Ti and Ni) and MoS2. However, the physical or chemical origin of this phenomenon is not yet clearly understood. In this work, density functional theory simulations are performed, employing pseudopotentials with very high basis sets to get insights of the charge transfer between metal and monolayer MoS2 through the inserted graphene layer. Our atomistic simulations on 16 different interfaces involving five different metals (Ti, Ag, Ru, Au, and Pt) reveal that (i) such a decrease in SBH is not consistent among various metals, rather an increase in SBH is observed in case of Au and Pt; (ii) unlike MoS2-metal interface, the projected dispersion of MoS2 remains preserved in any MoS2-graphene- metal system with shift in the bands on the energy axis. (iii) A proper choice of metal (e.g., Ru) may exhibit ohmic nature in a graphene-inserted MoS2-metal contact. These understandings would provide a direction in developing high-performance transistors involving heteroatomic layers as contact electrodes. (c) 2016 AIP Publishing LLC.