Calculation of the Entropy of Lattice Polymer Models from Monte Carlo Trajectories

While lattice models are used extensively for macromolecules (synthetic polymers proteins, etc.), calculation of the absolute entropy, S, and the free energy, F, from a given Monte Carlo (MC) trajectory is not straightforward. Recently, we have developed the hypothetical scanning MC (HSMC) method for calculating S and F of fluids. Here we extend HSMC to self-avoiding walks on a square lattice and discuss its wide applicability to complex polymer lattice models. HSMC is independent of existing techniques and thus constitutes an independent research tool; it provides rigorous upper and lower bounds for F, which can be obtained from a very small sample and even from a single chain conformation.     

Where is the sodium in self-assembled monolayers of single-stranded DNA?

Monolayers of single-stranded DNA (ssDNA) immobilized on surfaces form the basis of a number of important biotechnology applications, including DNA microarrays and biosensors. The organization of ssDNA as layer on a solid substrate allows one to investigate various properties of the DNA in a controlled manner and to use DNA for analytical applications as well as for exploring futuristic schemes for molecular electronics. It is commonly assumed that the adsorbed DNA layer contains some structural water and the cations. Here we show, based on XPS studies, that when monolayers of ssDNA are formed from sodium phosphate buffer and washed thoroughly, no Na+ signal is detected. A finite concentration of ions is observed when the DNA is made from a solution of Mg2+ ions, but it is still only a fifth of what it would be if all the phosphate ions were fully neutralized by the metal cations.     

Solution-processed anodes from layer-structure materials for high-efficiency polymer light-emitting diodes

The development of low-cost, large-area electronic applications requires the deposition of active materials in simple and inexpensive techniques at room temperature, properties usually associated with polymer films. In this study, we demonstrate the integration of solution-processed inorganic films in light-emitting diodes. The layered transition metal dichalcogenide (LTMDC) films are deposited through Li intercalation and exfoliation in aqueous solution and partially oxidized in an oxygen plasma generator. The chemical composition and thickness of the LTMDC and corresponding transition metal oxide (TMO) films are investigated by X-ray photoelectron spectroscopy. The morphology and topography of the films are studied by atomic force microscopy. X-ray powder diffraction is used to determine the orientation of the LTMDC film. Finally, the LTMDC and their corresponding oxides are utilized as hole-injecting and electron-blocking materials in polymer light-emitting diodes with the general structure ITO/LTMDC/TMO/polyfluorene/Ca/Al. Efficient hole injection and electron blocking by the inorganic layers result in outstanding device performance and high efficiency.     

Polymer collapse, protein folding, and the percolation threshold

We study the transition of polymers in the dilute regime from a swollen shape at high temperatures to their low-temperature structures. The polymers are modeled by a single self-avoiding walk (SAW) on a lattice for which l of the monomers (the H monomers) are self-attracting, i.e., if two nonbonded H monomers become nearest neighbors on the lattice they gain energy of interaction (epsilon = -/epsilon/); the second type of monomers, denoted P, are neutral. This HP model was suggested by Lau and Dill (Macromolecules 1989, 22, 3986-3997) to study protein folding, where H and P are the hydrophobic and polar amino acid residues, respectively. The model is simulated on the square and simple cubic (SC) lattices using the scanning method. We show that the ground state and the sharpness of the transition depend on the lattice, the fraction g of the H monomers, as well as on their arrangement along the chain. In particular, if the H monomers are distributed at random and g is larger than the site percolation threshold of the lattice, a collapsed transition is very likely to occur. This conclusion, drawn for the lattice models, is also applicable to proteins where an effective lattice with coordination number between that of the SC lattice and the body centered cubic lattice is defined. Thus, the average fraction of hydrophobic amino acid residues in globular proteins is found to be close to the percolation threshold of the effective lattice.     

Optimal vibration control of flexible spacecraft during a minimum-time maneuver

This paper is concerned with the simultaneous maneuver and vibration control of a flexible spacecraft. The problem is solved by means of a perturbation approach whereby the slewing of the spacecraft regarded as rigid represents the zero-order problem and the control of vibration, as well as of perturbations from the rigid-body maneuver, represents the first-order problem. The zero-order control is to be carried out in minimum time, which implies bang-bang control. On the other hand, the first-order control is a time-dependent linear quadratic regulator including integral feedback and prescribed convergence rate.This research was sponsored by USAF/ASD and AFOSR Research Grant F33615-86-C-3233 monitored by Drs. A. K. Amos and V. B. Venkayya, whose support is greatly appreciated.     

Calculation of the Absolute Free Energy of Binding and Related Entropies with the HSMD-TI Method: The FKBP12-L8 Complex

The hypothetical scanning molecular dynamics (HSMD) method is used here for calculating the absolute free energy of binding, ΔA(0) of the complex of the protein FKBP12 with the ligand SB2 (also denoted L8) - a system that has been studied previously for comparing the performance of different methods. Our preliminary study suggests that considering long-range electrostatics is imperative even for a hydrophobic ligand such as L8. Therefore the system is modeled by the AMBER force field using Particle Mesh Ewald (PME). HSMD consists of three stages applied to both the ligand-solvent and ligand-protein systems. (1) A small set of system configurations (frames) is extracted from an MD trajectory. (2) The entropy of the ligand in each frame is calculated by a reconstruction procedure. (3) The contribution of water and protein to ΔA(0) is calculated for each frame by gradually increasing the ligand-environment interactions from zero to their full value using thermodynamic integration (TI). Unlike the conventional methods, the structure of the ligand is kept fixed during TI, and HSMD is thus free from the end-point problem encountered with the double annihilation method (DAM); therefore, the need for applying restraints is avoided. Furthermore, unlike the conventional methods, the entropy of the ligand and water is obtained directly as a byproduct of the simulation. In this paper, in addition to the difference in the internal entropies of the ligand in the two environments, we calculate for the first time the external entropy of the ligand, which provides a measure for the size of the active site. We obtain ΔA(0) = -10.7 ±1.0 as compared to the experimental values -10.9 and -10.6 kcal/mol. However, a protein/water system treated by periodic boundary conditions grows significantly with increasing protein size and the computation of ΔA(0) would become expensive by all methods. Therefore, we also apply HSMD to FKBP12-L8 described by the GSBP/SSBP model of Roux's group (implemented in the software CHARMM) where only part of the protein and water around the active site are considered and long-range electrostatic effects are taken into account. For comparison this model was also treated by the double decoupling method (DDM). The two methods have led to comparable results for ΔA(0) which are somewhat lower than the experimental value. The ligand was found to be more confined in the active site described by GSBP/SSBP than by PME where its entropy in solvent is larger than in the active site by 1.7 and by 5.5 kcal/mol, respectively.     

Exploring the free energy surfaces of clusters using reconnaissance metadynamics

A new approach is proposed for exploring the low-energy structures of small to medium-sized aggregates of atoms and molecules. This approach uses the recently proposed reconnaissance metadynamics method [G. A. Tribello, M. Ceriotti, and M. Parrinello. Proc. Natl. Acad. Sci. U.S.A. 107(41), 17509 (2010)] in tandem with collective variables that describe the average structure of the coordination sphere around the atoms/molecules. We demonstrate this method on both Lennard-Jones and water clusters and show how it is able to quickly find the global minimum in the potential energy surface, while exploring the finite temperature free energy surface.