Anharmonic vibrational spectroscopy calculations with electronic structure potentials: comparison of MP2 and DFT for organic molecules

Density functional theory (DFT) technique is the most commonly used approach when it comes to computation of vibrational spectra of molecular species. In this study, we compare anharmonic spectra of several organic molecules such as allene, propyne, glycine, and imidazole, computed from ab initio MP2 potentials and DFT potentials based on commonly used BLYP and B3LYP functionals. Anharmonic spectra are obtained using the direct vibrational self-consistent field (VSCF) method and its correlation-corrected extension (CC-VSCF). The results of computations are compared with available experimental data. It is shown that the most accurate vibrational frequencies are obtained with the MP2 method, followed by the DFT/B3LYP method, while DFT/BLYP results are often unsatisfactory. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

Crystallographic analysis of metal-ion binding to human ubiquitin

The metal-binding ability of human ubiquitin (hUb) towards a selection of biologically relevant metal ions and complexes has been probed. Different techniques have been used to obtain crystals suitable for crystallographic analysis. In the first type of experiments, crystals of hUb have been soaked in solutions containing copper(II) acetate and two metallodrugs, Zeise salt (K[PtCl(3)(η(2)-C(2)H(4))]·H(2)O) and cisplatin (cis-[PtCl(2)(NH(3))(2)]). The Zeise salt is used in a test for hepatitis, whereas cisplatin is one of the most powerful anticancer drugs in clinical use. The Zeise salt readily reacts with hUb crystals to afford an adduct with three platinum residues per protein molecule, Pt(3)-hUb. In contrast, copper(II) acetate and cisplatin were found to be unreactive for contact times up to one hour and to cause degradation of the hUb crystals for longer times. In the second type of experiments, hUb was cocrystallized with a solution of copper(II) or zinc(II) acetate or cisplatin. Zinc(II) acetate gives, at low metal-to-protein molar ratios (8:1), crystals containing one metal ion per three molecules of protein, Zn-hUb(3) (already reported in previous work), whereas at high metal-to-protein ratios (70:1) gives crystals containing three Zn(II) ions per protein molecule, Zn(3)-hUb. In contrast, once again, copper(II) acetate and cisplatin, even at low metal-to-protein ratios, do not give crystalline material. In the soaking experiment, the Zeise anion leads to simultaneous platination of His68, Met1, and Lys6. Present and previous results of cocrystallization experiments performed with Zn(II) and other Group 12 metal ions allow a comprehensive understanding of the metal-ion binding properties of hUb with His68 as the main anchoring site, followed by Met1 and carboxylic groups of Glu16, Glu18, Glu64, Asp21, and Asp32, to be reached. In the case of platinum, Lys6 can also be a binding site. The amount of bound metal ion, with respect to that of the protein, appears to be a relevant parameter influencing crystal packing.     

Gas transport in TiO2 nanoparticle-filled poly(1-trimethylsilyl-1-propyne)

Titanium dioxide (TiO₂) nanoparticles were dispersed via solution processing in poly(1-trimethylsilyl-1-propyne) (PTMSP) to form nanocomposite films. Nanoparticle dispersion was investigated using atomic force microscopy and transmission electron microscopy. At low-particle loadings, nanoparticles were dispersed individually and in nanoscale aggregates. At high-particle loadings, some nanoparticles formed micron-sized aggregates. The gas transport and density exhibited a strong dependence on nanoparticle loading. At low-TiO₂ loadings, the composite density was similar to or slightly higher than that predicted by a two-phase additive model. However, at particle loadings exceeding approximately 7 nominal vol.%, the density was markedly lower than predicted, suggesting that the particles induced the creation of void space within the nanocomposite. For example, when the TiO₂ nominal volume fraction was 0.35, the polymer/particle composite density was 40% lower than expected based on a two-phase additive model for density. At low-nanoparticle loading, light gas permeability was lower than that of the unfilled polymer. At higher nanoparticle loadings, light gas permeability (i.e., CO₂, N₂, and CH₄) increased to more than four times higher than in unfilled PTMSP. At most, selectivity changed only slightly with particle loading.     

Conformational transitions of glycine induced by vibrational excitation of the O-H stretch

Vibrational energy flow and conformational transitions following excitation of the OH stretching mode of the most stable conformer of glycine are studied by classical trajectories. "On the fly" simulations with the PM3 semiempirical electronic structure method for the potential surface are used. Initial conditions are selected to correspond to the ν=1 excitation of the OH stretch. The main findings are: (1) An an equilibrium-like ratio is established between the populations of the 3 lowest-lying conformers after about 10 picoseconds. (2) There is a high probability throughout the 150 ps of the simulations for finding the molecule in geometries far from the equilibrium structures of the lowest-energy conformers. (3) Energy from the initial excited OH (ν=1) stretch flows preferentially to 5 other vibrational modes, including the bending motion of the H atom. (4) RRK theory yields conformational transition rates that deviate substantially from the classical trajectory results. Possible implication of these results for vibrational energy flow and conformational transitions in small biological molecules are discussed.     

Small subgraphs of random regular graphs

The main aim of this short paper is to answer the following question. Given a fixed graph H, for which values of the degree d does a random d-regular graph on n vertices contain a copy of H with probability close to one?     

On the Value of a Random Minimum Weight Steiner Tree

Consider a complete graph on n vertices with edge weights chosen randomly and independently from an exponential distribution with parameter 1. Fix k vertices and consider the minimum weight Steiner tree which contains these vertices. We prove that with high probability the weight of this tree is (1+o(1))(k-1)(log n-log k)/n when k =o(n) and n.* Research supported in part by NSF grant DSM9971788 Research supported in part by NSF grants DMS-0106589, CCR-9987845 and by the State of New Jersey. Part of this research was done while visiting IBM T. J. Watson Research Center.     

A New Lower Bound For A Ramsey-Type Problem

Let 3 ≤ r < s be fixed integers and let G be a graph on n vertices not containing a complete graph on s vertices. The main aim of this paper is to provide a new lower bound on the size of the maximum subset of G without a copy of complete graph K_r. Our results substantially improve previous bounds of Krivelevich and Bollobás and Hind. * Research supported in part by NSF grants DMS-0106589, CCR-9987845 and by the State of New Jersey. Part of this research was done while visiting Microsoft Research.     

Optimal pricing for a heterogeneous portfolio for a given risk factor and convex distance measure

Consider a portfolio containing heterogeneous risks. The premiums of the policyholders might not cover the amount of the payments which an insurance company pays the policyholders. When setting the premium, this risk has to be taken into consideration. On the other hand the premium that the insured pays has to be fair. This fairness is measured by a function of the difference between the risk and the premium paid—we call this function a distance function. For a given small probability of insolvency, we find the premium for each class, such that the distance function is minimized. Next we formulate and solve the dual problem, which is minimizing the insolvency probability under the constraint that the distance function does not exceed a given level. This paper generalizes a previous paper [Zaks, Y., Frostig, E., Levikson, B., 2006. Optimal pricing of a heterogeneous portfolio for a given risk level. Astin Bull. 36 (1), 161–185] where only a square distance function was considered.     

Evaluation of extraction methods for quantification of aqueous fullerenes in urine

There is a growing concern about the human and environmental health effects of fullerenes (e.g., C(60)) due to their increasing application in research, medicine, and industry. Toxicological and pharmacokinetic research requires standard methods for extraction and detection of fullerenes from biological matrices such as urine. The present study validates the use of liquid-liquid extraction (LLE) and solid-phase extraction (SPE) methods in conjunction with liquid chromatography-mass spectrometry (LC-MS) for the quantitative determination of C(60) in human and synthetic urine as compared with ultrapure water. Glacial acetic acid, which is necessary to prevent emulsions during LLE, inhibited C(60) detection by LC-MS, but this could be mitigated with evaporation. Aqueous C(60) aggregates (nC(60)) were spiked at 180 μg/L into the components of a synthetic urine recipe to determine their individual impacts on extraction and detection. Urea, creatinine, and a complex protein (i.e., gelatin) were found to impair SPE, leading to a low recovery rate of 43 ± 4% for C(60) spiked into human urine. In contrast, C(60) was consistently recovered from synthetic matrices using LLE, and recovery in human urine was 80 ± 6%. These results suggest that LLE combined with LC-MS is suitable for studying the clearance of fullerenes from the body. LLE is a robust technique that holds promise for extracting C(60) from other complex biological matrices (e.g., blood, sweat, amniotic fluid) in toxicological studies, enabling a better understanding of the behavior of fullerenes in human and animal systems and facilitating a more comprehensive risk evaluation of fullerenes.