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.     

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.     

Relative stability of the open and closed conformations of the active site loop of streptavidin

The eight-residue surface loop, 45-52 (Ser, Ala, Val, Gly, Asn, Ala, Glu, Ser), of the homotetrameric protein streptavidin has a "closed" conformation in the streptavidin-biotin complex, where the corresponding binding affinity is one of the strongest found in nature (ΔG ～ -18 kcal∕mol). However, in most of the crystal structures of apo (unbound) streptavidin, the loop conformation is "open" and typically exhibits partial disorder and high B-factors. Thus, it is plausible to assume that the loop structure is changed from open to closed upon binding of biotin, and the corresponding difference in free energy, ΔF = F(open) - F(closed) in the unbound protein, should therefore be considered in the total absolute free energy of binding. ΔF (which has generally been neglected) is calculated here using our "hypothetical scanning molecular-dynamics" (HSMD) method. We use a protein model in which only the atoms closest to the loop are considered (the "template") and they are fixed in the x-ray coordinates of the free protein; the x-ray conformation of the closed loop is attached to the same (unbound) template and both systems are capped with the same sphere of TIP3P water. Using the force field of the assisted model building with energy refinement (AMBER), we carry out two separate MD simulations (at temperature T = 300 K), starting from the open and closed conformations, where only the atoms of the loop and water are allowed to move (the template-water and template-loop interactions are considered). The absolute F(open) and F(closed) (of loop + water) are calculated from these trajectories, where the loop and water contributions are obtained by HSMD and a thermodynamic integration (TI) process, respectively. The combined HSMD-TI procedure leads to total (loop + water) ΔF = -27.1 ± 2.0 kcal∕mol, where the entropy TΔS constitutes 34% of ΔF, meaning that the effect of S is significant and should not be ignored. Also, ΔS is positive, in accord with the high flexibility of the open loop observed in crystal structures, while the energy ΔE is unexpectedly negative, thus also adding to the stability of the open loop. The loop and the 250 capped water molecules are the largest system studied thus far, which constitutes a test for the efficiency of HSMD-TI; this efficiency and technical issues related to the implementation of the method are also discussed. Finally, the result for ΔF is a prediction that will be considered in the calculation of the absolute free energy of binding of biotin to streptavidin, which constitutes our next project.     

Secondary electron emission control in X-ray photoelectron spectroscopy

Secondary electron emission (SEE) is a major player in surface charging during X-ray photoelectron spectroscopy (XPS); its characteristics and applicability as a current source for electrical measurements are studied. We employ sample biasing and a top retarding grid to control the photoelectron current, and further compare their I–V characteristics with direct spectroscopy of the secondary electrons. Using silica-coated gold substrates, the effect of sample work function on the emitted secondary electrons is shown and fine control over the surface potential gradients, in the range of 10–100 meV, is achieved. XPS-based chemically resolved electrical measurements (CREM) can thus be extended to the positive current regime.     

Energy level and band alignment for GaAs-alkylthiol monolayer-Hg junctions from electrical transport and photoemission experiments

A series of p- and n-GaAs-S-C(n)H(2n+1) || Hg junctions are prepared, and the electronic transport through them is measured. From current-voltage measurements, we find that, for n-GaAs, transport occurs by both thermionic emission and tunneling, with the former dominating at low forward bias and the latter dominating at higher forward bias. For p-GaAs, tunneling dominates at all bias voltages. By combining the analysis of the transport data with results from direct and inverse photoemission spectroscopy, we deduce an energy band diagram of the system, including the tunnel barrier and, with this barrier and within the Simmons tunneling model, extract an effective mass value of 1.5-1.6m(e) for the electronic carriers that cross the junctions. We find that transport is well-described by lowest unoccupied and highest occupied states at 1.3-1.4 eV above and 2.0-2.2 eV below the Fermi level. At the same time, the photoemission data indicate that there are continua of states from the conduction band minimum and the valence band maximum, the density of which varies with energy. On the basis of our results, it appears likely that, for both types of junctions, electrons are the main carrier type, although holes may contribute significantly to the transport in the p-GaAs system.