Importance of van der Waals Interactions in QM/MM Simulations

The importance of accurately treating van der Waals interactions between the quantum mechanical (QM) and molecular mechanical (MM) atoms in hybrid QM/MM simulations has been investigated systematically. First, a set of van der Waals (vdW) parameters was optimized for an approximate density functional method, the self-consistent charge-tight binding density functional (SCC-DFTB) approach, based on small hydrogen-bonding clusters. The sensitivity of condensed phase observables to the SCC-DFTB vdW parameters was then quantitatively investigated by SCC-DFTB/MM simulations of several model systems using the optimized set and two sets of extreme vdW parameters selected from the CHARMM22 forcefield. The model systems include a model FAD molecule in solution and a solvated enediolate, and the properties studied include the radial distribution functions of water molecules around the solute (model FAD and enediolate), the reduction potential of the model FAD and the potential of mean force for an intramolecular proton transfer in the enediolate. Although there are noticeable differences between parameter sets for gas-phase clusters and solvent structures around the solute, thermodynamic quantities in the condensed phase (e.g., reduction potential and potential of mean force) were found to be less sensitive to the numerical values of vdW parameters. The differences between SCC-DFTB/MM results with the three vdW parameter sets for SCC-DFTB atoms were explained in terms of the effects of the parameter set on solvation. The current study has made it clear that efforts in improving the reliability of QM/MM methods for energetical properties in the condensed phase should focus on components other than van der Waals interactions between QM and MM atoms.     

Temperature-sensitive polypeptide brushes-coated mesoporous silica nanoparticles for dual-responsive drug release

A multifunctional nanohybrid based on mesoporous silica nanoparticle and biocompatible polypeptide was fabricated for targeted and dual-responsive therapy of tumor cells.     

ICP-AES determination of trace elements after preconcentrated with p-dimethylaminobenzaldehyde-modified nanometer SiO2 from sample solution

A new method that utilizes p-dimethylaminobenzaldehyde-modified nanometer SiO₂ (SiO₂-p-DMABD) as a solid phase extractant has been developed for simultaneous preconcentration of trace Cr(III), Cu(II), Fe(III) and Pb(II) prior to the measurement by inductively coupled plasma atomic emission spectrometry (ICP-AES). The preconcentration conditions of analytes were investigated, including the pH value, the shaking time, the mass of sorbent, the sample flow rate and volume, the elution condition and the interfering ions. The adsorption capacity of nanometer SiO₂-p-DMABD was found to be (mg g^− 1) Cr(III): 6.2, Cu(II): 18.6, Fe(III): 4.7 and Pb(II): 6.0 at pH 4. The adsorbed metals were quantitatively eluted with 4 mL of 1.0 mol L^− 1 HCl. According to the definition of IUPAC, the detection limits (3σ) of this method for Cr(III), Cu(II), Fe(III) and Pb(II) were 0.79, 1.27, 0.40 and 1.79 ng mL^− 1, respectively. The proposed method achieved satisfied results when it was applied to the determination of trace Cr(III), Cu(II), Fe(III) and Pb(II) in biological and water samples.     

Synthesis, structure, and fluorescence of the novel cadmium(II)-trimesate coordination polymers with different coordination architectures

Three novel complexes, Cd3tma2*13H2O (1), Cd3tma2*dabco*2H2O (2), and Cd3Htma3*8H2O (3) (tma = trimesate), of cadmium(II)-trimesate coordination polymers are obtained from hydrothermal reaction. 1 (C18H32O25Cd3) crystallizes in the monoclinic C2/c space group [a = 18.985(2) A, b = 7.3872(6) A, c = 20.432(2) A, = 97.1660(10), and Z = 4]. 2 (C24H22N2O14Cd3) crystallizes in the monoclinic P2(1)/c space group [a = 10.1323(2) A, b = 19.5669(5) A, c = 13.15880(10) A, = 108.9810(10), and Z = 4]. 3 (C27H28O26Cd3) belongs to the trigonal P31c space group [a = 15.7547(3) A, b = 15.7547(3) A, c = 7.93160(10) A, and Z = 2]. The Cd(II) centers in the three complexes are bridged by tma ligands in the coordination fashion of unidentate, bridging unidentate, bidentate, chelating bis-bidentate, chelating/bridging bis-bidentate, or chelating/bridging bidentate to form the T-shaped molecular bilayer motif for 1, chicken-wire-like motif for 2, and honeycomb-like porous structure for 3, respectively, in which the T-shaped molecular bilayer motif and chicken-wire-like motif are further interlinked in interdigitating or alternating fashion to construct the different coordination architectures. These three complexes exhibit strong fluorescent emission bands at 355 nm (lambda(ex) = 220 nm) for 1, 437 nm (lambda(ex) = 365 nm) for 2, and 353 nm (lambda(ex) = 218 nm) for 3 in the solid state at room temperature.     

Preparation of multihollow P(St-MAA) particles by sequence soap-free/soap emulsion polymerization and followed by stepwise alkali/acid posttreatment

The effects of ionic emulsifier, sodium dodecylbenzene sulfate (SDBS), on the formation of the multihollow structures in sub-micron sized polymer particles produced by alkali/acid posttreatment were investigated. The original latex particles with narrow size distribution were synthesized by a new sequence emulsifier-free/emulsifier emulsion copolymerization of styrene (St) and methacrylic acid (MAA). Results indicated that the pore size decreased and the pore number increased with the increase of SDBS amount, and the morphology of the posttreated latex particles was also significantly influenced by the introducing time of SDBS in the preparation of the original latex particles, and a suitable introducing time was 3 h of polymerization.     

Stannous oxalate as a novel catalyst for the synthesis of polytrimethylene terephthalate

Stannous oxalate was prepared efficiently and characterized by XRD and FT-IR. It exhibited higher catalytic activity and had profitable effect than tetrabutyl titanate and stannous octoate for the synthesis of polytrimethylene terephthalate (PTT) via esterification-route. Over this catalyst, the degree of esterification of pure terephthalic acid was up to 94.4% at 260°C after 1.5 h, while the intrinsic viscosity and content of terminal carboxyl groups of the corresponding PTT polyester, polymerized at 260°C, 60 Pa for 2 h, was 0.8950 dL/g and 15 mol/t, respectively. Stannous oxalate was a promising catalyst for the synthesis of PTT polyester.     

CO2选择性氧化甲烷制C2烃催化剂La2O3/ZnO的氧化还原和CO2吸附性能

Ｌａ２Ｏ３／ＺｎＯ催化剂体系在以二氧化碳作为氧化剂的甲烷氧化偶联反应中具有很高的Ｃ２烃选择性和稳定性．采用ＣＯ２－ＴＰＤ－ＭＳ和ＴＰＲ技术考察了Ｌａ２Ｏ３／ＺｎＯ对ＣＯ２的吸附性质及其氧化还原行为．结果表明：（１）Ｌａ２Ｏ３／ＺｎＯ催化剂体系存在着强、弱两种碱中心，其中弱碱中心数量随样品中Ｌａ２Ｏ３含量增加而减少，强碱中心强度随样品中Ｌａ２Ｏ３含量增加而增强．（２）由于组分相互作用，高温下，Ｌａ２Ｏ３／ＺｎＯ易产生晶格氧空位，使之对ＣＯ２的吸附增强，吸附后的ＣＯ２与晶格氧作用形成立方晶型Ｌａ２Ｏ２ＣＯ３．（３）Ｌａ２Ｏ３／ＺｎＯ表面的Ｌａ３＋和Ｚｎ２＋可以部分被还原，由于组分间的相互作用，使得二者的还原都较单一组分存在时更难．（４）Ｈ２－ＣＯ２－Ｈ２氧化还原循环实验表明，Ｌａ２Ｏ３／ＺｎＯ表面被部分还原后，ＣＯ２可以将部分被还原的表面再氧化．在此基础上对Ｌａ２Ｏ３／ＺｎＯ催化剂上甲烷与ＣＯ２转化为Ｃ２烃的机制也进行了讨论．     

Synthesis of 2‐Dihydrooxadiazolinyl Chromones

Reactions of the substituted 2‐formyl chromones with aroylhydrazines gave corresponding 2‐(aroylhydrazonomethylidyne) chromones. Then 2‐(3′‐acetyl‐5′‐aryl‐2′,3′‐dihydro‐1′,3′,4′‐oxadiazol‐2′‐yl) chromones were prepared by these 2‐(aroylhydrazonomethylidyne) chromones under refluxing with Ac₂O. All target compounds were characterized through elemental analysis and IR, ¹H NMR, MS.     

Development and certification of a coal fly ash certified reference material for selected polycyclic aromatic hydrocarbons

The development and certification of a coal fly ash certified reference material (CRM) for polycyclic aromatic hydrocarbons (PAH) is described; this is the first natural matrix CRM for organic environmental analysis in China. The homogeneity and stability of this material have been tested by HPLC. The concentrations of several PAH were determined by use of two independent, different methods--solvent extraction-HPLC analysis with UV detection coupled with fluorescence detection (FLD) and solvent extraction, isolation with a silica column, and GC analysis with flame ionization detection (FID). Five certified values were determined: phenanthrene 7.1 +/- 2.6 microg g(-1), anthracene 2.0 +/- 0.8 microg g(-1), fluoranthene 7.4 +/- 1.9 microg g(-1), pyrene 7 +/- 2 microg g(-1), and benzo[a]pyrene 1.3 +/- 0.3 microg g(-1). Reference values for several other PAH are also suggested.     

pKa calculations in solution and proteins with QM/MM free energy perturbation simulations: a quantitative test of QM/MM protocols

The accuracy of biological simulations depends, in large part, on the treatment of electrostatics. Due to the availability of accurate experimental values, calculation of pKa provides stringent evaluation of computational methods. The generalized solvent boundary potential (GSBP) and Ewald summation electrostatic treatments were recently implemented for combined quantum mechanical and molecular mechanics (QM/MM) simulations by our group. These approaches were tested by calculating pKa shifts due to differences in electronic structure and electrostatic environment; the shifts were determined for a series of small molecules in solution, using various electrostatic treatments, and two residues (His 31, Lys 102) in the M102K T4-lysozyme mutant with large pKa shifts, using the GSBP approach. The calculations utilized a free energy perturbation scheme with the QM/MM potential function involving the self-consistent charge density functional tight binding (SCC-DFTB) and CHARMM as the QM and MM methods, respectively. The study of small molecules demonstrated that inconsistent electrostatic models produced results that were difficult to correct in a robust manner; by contrast, extended electrostatics, GSBP, and Ewald simulations produced consistent results once a bulk solvation contribution was carefully chosen. In addition to the electrostatic treatment, the pKa shifts were also sensitive to the level of the QM method and the scheme of treating QM/MM Coulombic interactions; however, simple perturbative corrections based on SCC-DFTB/CHARMM trajectories and higher level single point energy calculations were found to give satisfactory results. Combining all factors gave a root-mean-square difference of 0.7 pKa units for the relative pKa values of the small molecules compared to experiment. For the residues in the lysozyme, an accurate pKa shift was obtained for His 31 with multiple nanosecond simulations. For Lys 102, however, the pKa shift was estimated to be too large, even after more than 10 nanosecond simulations for each lambda window; the difficulty was due to the significant, but slow, reorganization of the protein and water structure when Lys 102 was protonated. The simulations support that Lys 102 is deprotonated in the X-ray structure and the protein is highly destabilized when this residue is protonated.