Isomolar semigrand ensemble molecular dynamics: development and application to liquid-liquid equilibria

An extended system molecular dynamics method for the isomolar semigrand ensemble (fixed number of particles, pressure, temperature, and fugacity fraction) is developed and applied to the calculation of liquid-liquid equilibria (LLE) for two Lennard-Jones mixtures. The method utilizes an extended system variable to dynamically control the fugacity fraction xi of the mixture by gradually transforming the identity of particles in the system. Two approaches are used to compute coexistence points. The first approach uses multiple-histogram reweighting techniques to determine the coexistence xi and compositions of each phase at temperatures near the upper critical solution temperature. The second approach, useful for cases in which there is no critical solution temperature, is based on principles of small system thermodynamics. In this case a coexistence point is found by running N-P-T-xi simulations at a common temperature and pressure and varying the fugacity fraction to map out the difference in chemical potential between the two species A and B (mu(A)-mu(B)) as a function of composition. Once this curve is known the equal-distance/equal-area criterion is used to determine the coexistence point. Both approaches give results that are comparable to those of previous Monte Carlo (MC) simulations. By formulating this approach in a molecular dynamics framework, it should be easier to compute the LLE of complex molecules whose intramolecular degrees of freedom are often difficult to properly sample with MC techniques.     

Nucleic acid reactivity: Challenges for next‐generation semiempirical quantum models

Semiempirical quantum models are routinely used to study mechanisms of RNA catalysis and phosphoryl transfer reactions using combined quantum mechanical (QM)/molecular mechanical methods. Herein, we provide a broad assessment of the performance of existing semiempirical quantum models to describe nucleic acid structure and reactivity to quantify their limitations and guide the development of next‐generation quantum models with improved accuracy. Neglect of diatomic differential overlap and self‐consistent density‐functional tight‐binding semiempirical models are evaluated against high‐level QM benchmark calculations for seven biologically important datasets. The datasets include: proton affinities, polarizabilities, nucleobase dimer interactions, dimethyl phosphate anion, nucleoside sugar and glycosidic torsion conformations, and RNA phosphoryl transfer model reactions. As an additional baseline, comparisons are made with several commonly used density‐functional models, including M062X and B3LYP (in some cases with dispersion corrections). The results show that, among the semiempirical models examined, the AM1/d‐PhoT model is the most robust at predicting proton affinities. AM1/d‐PhoT and DFTB3‐3ob/OPhyd reproduce the MP2 potential energy surfaces of 6 associative RNA phosphoryl transfer model reactions reasonably well. Further, a recently developed linear‐scaling “modified divide‐and‐conquer” model exhibits the most accurate results for binding energies of both hydrogen bonded and stacked nucleobase dimers. The semiempirical models considered here are shown to underestimate the isotropic polarizabilities of neutral molecules by approximately 30%. The semiempirical models also fail to adequately describe torsion profiles for the dimethyl phosphate anion, the nucleoside sugar ring puckers, and the rotations about the nucleoside glycosidic bond. The modeling of pentavalent phosphorus, particularly with thio substitutions often used experimentally as mechanistic probes, was problematic for all of the models considered. Analysis of the strengths and weakness of the models suggests that the creation of robust next‐generation models should emphasize the improvement of relative conformational energies and barriers, and nonbonded interactions. © 2015 Wiley Periodicals, Inc.     

Structural and rate studies of the 1,2-additions of lithium phenylacetylide to lithiated quinazolinones: influence of mixed aggregates on the reaction mechanism

The 1,2-addition of lithium phenylacetylide (PhCCLi) to quinazolinones was investigated using a combination of structural and rate studies. (6)Li, (13)C, and (19)F NMR spectroscopies show that deprotonation of quinazolinones and phenylacetylene in THF/pentane solutions with lithium hexamethyldisilazide affords a mixture of lithium quinazolinide/PhCCLi mixed dimer and mixed tetramer along with PhCCLi dimer. Although the mixed tetramer dominates at high mixed aggregate concentrations and low temperatures used for the structural studies, the mixed dimer is the dominant form at the low total mixed aggregate concentrations, high THF concentrations, and ambient temperatures used to investigate the 1,2-addition. Monitoring the reaction rates using (19)F NMR spectroscopy revealed a first-order dependence on mixed dimer, a zeroth-order dependence on THF, and a half-order dependence on the PhCCLi concentration. The rate law is consistent with the addition of a disolvated PhCCLi monomer to the mixed dimer. Investigation of the 1,2-addition of PhCCLi to an O-protected quinazolinone implicates reaction via trisolvated PhCCLi monomers.     

Nuclear resonance vibrational spectra of five-coordinate imidazole-ligated iron(II) porphyrinates

Nuclear resonance vibrational spectra have been obtained for six five-coordinate imidazole-ligated iron(II) porphyrinates, [Fe(Por)(L)] (Por = tetraphenylporphyrinate, octaethylporphyrinate, tetratolylporphyrinate, or protoporphyrinate IX and L = 2-methylimidazole or 1,2-dimethylimidazole). Measurements have been made on both powder and oriented crystal samples. The spectra are dominated by strong signals around 200-300 cm(-1). Although the in-plane and out-of-plane vibrations are seriously overlapped, oriented crystal spectra allow their deconvolution. Thus, oriented crystal experimental data, along with density functional theory (DFT) calculations, enable the assignment of key vibrations in the spectra. Molecular dynamics are also discussed. The nature of the Fe-N(Im) vibrations has been elaborated further than was possible from resonance Raman studies. Our study suggests that the Fe motions are coupled with the porphyrin core and peripheral groups motions. Both peripheral groups and their conformations have significant influence on the vibrational spectra (position and shape).     

Fluence Rate-Dependent Photobleaching of Intratumorally Administered Pc 4 Does not Predict Tumor Growth Delay

We examined effects of fluence rate on the photobleaching of the photosensitizer Pc 4 during photodynamic therapy (PDT) and the relationship between photobleaching and tumor response to PDT. BALB/c mice with intradermal EMT6 tumors were given 0.03 mg kg^?1 Pc 4 by intratumor injection and irradiated at 667 nm with an irradiance of 50 or 150 mW cm^?2 to a fluence of 100 J cm^?2. While no cures were attained, significant tumor growth delay was demonstrated at both irradiances compared with drug‐only controls. There was no significant difference in tumor responses to these two irradiances (P = 0.857). Fluorescence spectroscopy was used to monitor the bleaching of Pc 4 during irradiation, with more rapid bleaching with respect to fluence shown at the higher irradiance. No significant correlation was found between fluorescence photobleaching and tumor regrowth for the data interpreted as a whole. Within each treatment group, weak associations between photobleaching and outcome were observed. In the 50 mW cm^?2 group, enhanced photobleaching was associated with prolonged growth delay (P = 0.188), while at 150 mW cm^?2 this trend was reversed (P = 0.308). Thus, it appears that Pc 4 photobleaching is not a strong predictor of individual tumor response to Pc 4‐PDT under these treatment conditions.     

Sodium periodate as a primary oxidant for water-oxidation catalysts

Sodium periodate was characterized as a primary chemical oxidant for the catalytic evolution of oxygen at neutral pH using a variety of water-oxidation catalysts. The visible spectra of solutions formed from Cp*Ir(bpy)SO(4) during oxygen-evolution catalysis were measured. NMR spectroscopy suggests that the catalyst remains molecular after several turnovers with sodium periodate. Two of our [Cp*Ir(bis-NHC)][PF(6)](2) complexes, along with other literature catalysts, such as the manganese terpyridyl dimer, Hill's cobalt polyoxometallate, and Meyer's blue dimer, were also tested for activity. Sodium periodate was found to function only for water-oxidation catalysts with low overpotentials. This specificity is attributed to the relatively low oxidizing capability of sodium periodate solutions relative to solutions of other common primary oxidants. Studying oxygen-evolution catalysis by using sodium periodate as a primary oxidant may, therefore, provide preliminary evidence that a given catalyst has a low overpotential.     

Modelling the circularly and linearly polarised spectrum of [Ni(en)3]2+

The relative intensity and band shapes of the low energy spin-allowed transitions in the linearly polarised and circular dichroism spectrum of [Ni(en)(3)](2+) have been calculated using a time-dependent density functional theory approach. The effect of the trigonal ligand-field is minimal and no splitting of the bands is predicted by the simulations or observed experimentally. The 'd-d' transitions of the [Ni(en)(3)](2+) ion are electric dipole allowed but gain much of their intensity through Herzberg-Teller vibronic coupling. Its CD spectrum is dominated by the low energy band, which gains its rotatory strength through the magnetic dipole-allowed character of the parent octahedral transition and the electric dipole character due to the trigonal field. The simulation of the spectrum incorporates the contribution from all inducing vibrational modes with significant involvement of the {NiN(6)} unit. Vibrations which are centred on the chelate rings are not important in generating intensity, reflecting the localised d-d' character of the transitions. Simulated linearly polarised and circular dichroism spectra of such an open-shell system are presented for the first time and predict the essential elements of the experimental spectra.     

Low-temperature tunneling and rotational dynamics of the ammonium cations in (NH4)2B12H12

Low-temperature neutron scattering spectra of diammonium dodecahydro-closo-dodecaborate [(NH(4))(2)B(12)H(12)] reveal two NH(4)(+) rotational tunneling peaks (e.g., 18.5 μeV and 37 μeV at 4 K), consistent with the tetrahedral symmetry and environment of the cations. The tunneling peaks persist between 4 K and 40 K. An estimate was made for the tunnel splitting of the first NH(4)(+) librational state from a fit of the observed ground-state tunnel splitting as a function of temperature. At temperatures of 50 K-70 K, classical neutron quasi-elastic scattering appears to dominate the spectra and is attributed to NH(4)(+) cation jump reorientation about the four C(3) axes defined by the N-H bonds. A reorientational activation energy of 8.1 ± 0.6 meV (0.79 ± 0.06 kJ/mol) is determined from the behavior of the quasi-elastic linewidths in this temperature regime. This activation energy is in accord with a change in NH(4)(+) dynamical behavior above 70 K. A low-temperature inelastic neutron scattering feature at 7.8 meV is assigned to a NH(4)(+) librational mode. At increased temperatures, this feature drops in intensity, having shifted entirely to higher energies by 200 K, suggesting the onset of quasi-free NH(4)(+) rotation. This is consistent with neutron-diffraction-based model refinements, which derive very large thermal ellipsoids for the ammonium-ion hydrogen atoms at room temperature in the direction of reorientation.     

Continuous photochemical generation of catalytically active [CpRu]+ complexes from CpRu(η6-C6H6)PF6

Inter- and intramolecular ene-yne coupling reactions catalyzed by a species generated in situ via photolysis of CpRu(η(6)-C(6)H(6))PF(6)--an inexpensive, readily available, and shelf-stable complex--have been demonstrated under conditions of continuous flow. Importantly, the catalyst can be recovered quantitatively at the end of the reaction. Various functional groups are tolerated by the reaction, which affords skipped diene products in high yields.     

Monitoring cytoplasmic protein complexes with blue native gel electrophoresis and stable isotope labelling with amino acids in cell culture: analysis of changes in the 20S proteasome

Analysis of protein complexes is of increasing interest in the field of proteomics. A challenge is to develop methods for monitoring changes in the quantity and subunit composition of protein complexes on a proteome-wide scale. Here, we describe the combination of 1-D blue native polyacrylamide gel electrophoresis (BN-PAGE) with stable isotope labelling of amino acids in cell culture (SILAC) and tandem mass spectrometry (MS/MS). Cleared lysates from normal and perturbed samples, one incorporating heavy stable isotopes and the other light isotopes, are co-separated by blue native PAGE and then analysed and quantitated with MS/MS and appropriate software. This permits the analysis of cytoplasmic complexes. To demonstrate this technique, we explored how the 20S proteasome changes when the Pre9/α3 subunit, the only non-essential subunit of this complex, was deleted. Our results showed that ΔPre9/α3 cells can form the 20S proteasome complex, although with reduced efficiency. This involves an increase in expression of the α4 subunit. Our findings suggest this technique as an approach for the study of quantitative and qualitative differences in protein complexes, from cleared cell lysates.