A method to construct polyelectrolyte multilayers film containing gold nanoparticles

Gold nanoparticles in polyelectrolyte multilayers film can be easily prepared by repeating immersion of a substrate in poly(diallyl dimethylammonium) chloride (PDDA)-AuCl₄⁻ complexes solution followed by reduction Au³⁺ through heating. UV-vis spectroscopy, cyclic voltammetry (CV) and tapping-mode atomic force microscopy (AFM) are used to confirm the successful construction of the polyelectrolyte multilayers film and the formation of gold nanoparticles. The multilayers film shows electrocatalytic activity to dioxygen reduction.     

Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis

QM/MM methods have been developed as a computationally feasible solution to QM simulation of chemical processes, such as enzyme-catalyzed reactions, within a more approximate MM representation of the condensed-phase environment. However, there has been no independent method for checking the quality of this representation, especially for highly nonisotropic protein environments such as those surrounding enzyme active sites. Hence, the validity of QM/MM methods is largely untested. Here we use the possibility of performing all-QM calculations at the semiempirical PM3 level with a linear-scaling method (MOZYME) to assess the performance of a QM/MM method (PM3/AMBER94 force field). Using two model pathways for the hydride-ion transfer reaction of the enzyme dihydrofolate reductase studied previously (Titmuss et al., Chem Phys Lett 2000, 320, 169-176), we have analyzed the reaction energy contributions (QM, QM/MM, and MM) from the QM/MM results and compared them with analogous-region components calculated via an energy partitioning scheme implemented into MOZYME. This analysis further divided the MOZYME components into Coulomb, resonance and exchange energy terms. For the model in which the MM coordinates are kept fixed during the reaction, we find that the MOZYME and QM/MM total energy profiles agree very well, but that there are significant differences in the energy components. Most significantly there is a large change (approximately 16 kcal/mol) in the MOZYME MM component due to polarization of the MM region surrounding the active site, and which arises mostly from MM atoms close to (<10 A) the active-site QM region, which is not modelled explicitly by our QM/MM method. However, for the model where the MM coordinates are allowed to vary during the reaction, we find large differences in the MOZYME and QM/MM total energy profiles, with a discrepancy of 52 kcal/mol between the relative reaction (product-reactant) energies. This is largely due to a difference in the MM energies of 58 kcal/mol, of which we can attribute approximately 40 kcal/mol to geometry effects in the MM region and the remainder, as before, to MM region polarization. Contrary to the fixed-geometry model, there is no correlation of the MM energy changes with distance from the QM region, nor are they contributed by only a few residues. Overall, the results suggest that merely extending the size of the QM region in the QM/MM calculation is not a universal solution to the MOZYME- and QM/MM-method differences. They also suggest that attaching physical significance to MOZYME Coulomb, resonance and exchange components is problematic. Although we conclude that it would be possible to reparameterize the QM/MM force field to reproduce MOZYME energies, a better way to account for both the effects of the protein environment and known deficiencies in semiempirical methods would be to parameterize the force field based on data from DFT or ab initio QM linear-scaling calculations. Such a force field could be used efficiently in MD simulations to calculate free energies.     

Single polymer molecules adsorbed to mica and the oppositely charged polymer/surfactant complexes formed at the air–water interface visualized by atomic force microscopy

In the present study, the structure and morphology of single sodium poly(styrenesulfonate) (PSS) molecules adsorbed to mica surface from the natural aqueous solution is investigated using atomic force microscopy technique. Results show that single PSS molecules are observed which show a morphology of wormlike coils. Meanwhile, single sodium poly(styrenesulfonate) (PSS)/Hexadecyltrimethylammonium bromide (CTAB) complexes deposited on mica from the air–water interface are also observed. However, the PSS/CTA⁺ complexes show different conformations by appearing in the morphology of circular patches. Experimental data are in fair agreement with the theoretical analysis.     

A quantitative method for dual‐pass electrostatic force microscopy phase measurements

Electrostatic force microscopy (EFM) has become a powerful tool for investigating charges on surfaces. The use of phase measurement in EFM is a direct and fast way to detect electrostatic force gradients, but only qualitatively. With the dual‐pass scheme, the phase signal at lifted height is often assumed to exclude any influences from the topography, but it does not. We report the collection of both topography and phase data by EFM on charged, micron‐sized metal wires. In order to quantify the electrostatic force, a cone model and finite element analysis are provided to integrate the force gradient from the phase signal. Copyright © 2007 John Wiley & Sons, Ltd.     

A case study in conformation design: learning by doing

Efforts are described to design simple, fully flexible but conformationally preorganised omega-hydroxy-nonanoic acids that could serve as the conformation controlling unit in analogues of the potent protein-kinase C activator aplysiatoxin. Such analogues are macrodilactones incorporating the designed omega-hydroxy-nonanoic acid and 3,4-dihydroxy-pentanoic acid, which contains the pharmacophoric groups. The design process (replacement of CH(2) groups by an oxygen atom, annelation of a six-membered ring and placement of alkyl substituents) of the omega-hydroxy-nonanoic acids was monitored by force-field calculations. In the end of this process simple analogues of aplysiatoxin are proposed in which the proper disposition of the pharmacophoric groups is secured by a conformationally flexible but preorganised template structure as part of the macrodilactone ring.     

ELASTIC BEHAVIOR OF PROTEIN-LIKE SINGLE CHAIN

The conformational properties and elastic behaviors of protein-like single chains in the process of tensile elongation were investigated by means of Monte Carlo method. The sequences of protein-like single chains contain two types of residues: hydrophobic (H) and hydrophilic (P). The average conformations and thermodynamics statistical properties of protein-like single chains with various elongation ratio λ were calculated. It was found that the mean-square end-to-end distance r increases with elongation ratio,λ. The tensor eigenvalues ratio of : decreases with elongation ratio λ for short (HP)x protein-like polymers, however, the ratio of : increases with elongation ratioλ,especially for long (H)x sequence. Average energy per bond increases with elongation ratioλ, especially for(H)x protein-like single chains. Helmholtz free energy per bond also increases with elongation ratioλ. Elastic force (f), energy contribution to force (fU) and entropy contribution to force (fs) for different protein-like single chains were also calculated.These investigations may provide some insights into elastic behaviors of proteins.     

Flexible matching of test ligands to a 3D pharmacophore using a molecular superposition force field: Comparison of predicted and experimental conformations of inhibitors of three enzymes

Summary A computer procedure TFIT, which uses a molecular superposition force field to flexibly match test compounds to a 3D pharmacophore, was evaluated to find out whether it could reliably predict the bioactive conformations of flexible ligands. The program superposition force field optimizes the overlap of those atoms of the test ligand and template that are of similar chemical type, by applying an attractive force between atoms of the test ligand and template which are close together and of similar type (hydrogen bonding, charge, hydrophobicity). A procedure involving Monte Carlo torsion perturbations, followed by torsional energy minimization, is used to find conformations of the test ligand which cominimize the internal energy of the ligand and the superposition energy of ligand and template. The procedure was tested by applying it to a series of flexible ligands for which the bioactive conformation was known experimentally. The 15 molecules tested were inhibitors of thermolysin, HIV-1 protease or endothiapepsin for which X-ray structures of the bioactive conformation were available. For each enzyme, one of the molecules served as a template and the others, after being conformationally randomized, were fitted. The fitted conformation was then compared to the known binding geometry. The matching procedure was successful in predicting the bioactive conformations of many of the structures tested. Significant deviation from experimental results was found only for parts of molecules where it was readily apparent that the template did not contain sufficient information to accurately determine the bioactive conformation.     

Interconversion behavior of the C-H bond in the CH4+ radical cation: ab initio molecular dynamics study

Thermal motion of CH4+ is investigated by performing an ab initio molecular dynamics method with the second-order M?ller-Plesset (MP2)/6-311G** force field. In the trajectories obtained at 400 K, we have observed rapid interconversion behavior of the geometrical parameters of CH4+ with the frequency of 0.6/ps, where the C-H pair forming the small angle around 55 degrees is switched to another pair on subpicosecond time scale. The switching patterns are found to be classified into the following two types. Type 1: one C-H of the small angled C-H pair is switched to one C-H of the other C-H pair. Type 2: the small angled C-H pair is switched to the other C-H pair, which has been newly observed in the present ab initio MD calculation. The four C-H bonds of CH4+ are characterized by the long and short C-H bonds in a time region of the trajectories, and also for the time-evolution of C-H bonds such interconversion behavior is observed. The switching patterns of the geometrical parameters are compared with those in the interconversion scheme between six equivalent C2v symmetry structures of CH4+ [Paddon-Row, M. N. et al., J Am Chem Soc 1985, 107, 7696]. We have also investigated the electronic energy fluctuation due to thermal motion of CH4+. The standard deviation of total electronic energy at 400 K is evaluated to be 1.2 kcal/mol.     

Thermodynamics of liquid Fe-Si and Fe-Ge alloys

The enthalpies of mixing of liquid binary Fe-Ge (1765±5 K) and Fe-Si (1750±5 K) alloys were determined using a high-temperature isoperibolic calorimeter. The thermodynamic properties of Fe-Ge melts were also studied by electromotive force method in the temperature range of 1250-1580 K. The comparison of our measurement results with literature data has been performed. The extreme negative values of integral enthalpy of mixing and alternating-sign deviations from Raoult's low for germanium can be explained by the influence of binary clusters formation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Radical addition reactions: factors determining the transition-state geometry

Interatomic distances in the reaction centers of the addition reactions of (i) H^· to the C=C, C=O, N≡C, and C≡C bonds, (ii) ^·CH₃ radical to the C=C, C=O, and C≡C bonds, and (iii) alkyl, aminyl, and alkoxyl radicals to olefin C=C bonds were determined using a new semiempirical method for calculating transition-state geometries of radical reactions. For all reactions of the type X^· + Y=Z → X— Y—Z^· the r ^# _X...Y distance in the transition state is a linear function of the enthalpy of reaction. Parameters of this dependence were determined for seventeen classes of radical addition reactions. The bond elongation, Δr ^# _X...Y, in the transition state decreases as the triplet repulsion, electronegativity difference between the atoms X and Y in the reaction center, and the force constant of the attacked multiple bond increase. __________ Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 894–902, April, 2005.