Modeling of thermal transport in pillared-graphene architectures.

Carbon nanotubes (CNT) and graphene are considered as potential future candidates for many nano/microscale integrated devices due to their superior thermal properties. Both systems, however, exhibit significant anisotropy in their thermal conduction, limiting their performance as three-dimensional thermal transport materials. From thermal management perspective, one way to tailor this anisotropy is to consider designing alternative carbon-based architectures. This paper investigates the thermal transport in one such novel architecture-a pillared-graphene (PG) network nanostructure which combines graphene sheets and carbon nanotubes to create a three-dimensional network. Nonequilibrium molecular dynamics simulations have been carried out using the AIREBO potential to calculate the thermal conductivity of pillared-graphene structures along parallel (in-plane) as well as perpendicular (out-of-plane) directions with respect to the graphene plane. The resulting thermal conductivity values for PG systems are discussed and compared with simulated values for pure CNT and graphite. Our results show that in these PG structures, the thermal transport is governed by the minimum interpillar distance and the CNT-pillar length. This is primarily attributed to scattering of phonons occurring at the CNT-graphene junctions in these nanostructures. We foresee that such architecture could potentially be used as a template for designing future structurally stable microscale systems with tailorable in-plane and out-of-plane thermal transport.

Nature of molecular interactions of peptides with gold, palladium, and Pd-Au bimetal surfaces in aqueous solution.

We investigated molecular interactions involved in the selective binding of several short peptides derived from phage-display techniques (8-12 amino acids, excluding Cys) to surfaces of Au, Pd, and Pd-Au bimetal. The quantitative analysis of changes in energy and conformation upon adsorption on even {111} and {100} surfaces was carried out by molecular dynamics simulation using an efficient computational screening technique, including 1000 explicit water molecules and physically meaningful peptide concentrations at pH = 7. Changes in chain conformation from the solution to the adsorbed state over the course of multiple nanoseconds suggest that the peptides preferably interact with vacant sites of the face-centered cubic lattice above the metal surface. Residues that contribute to binding are in direct contact with the metal surfaces, and less-binding residues are separated from the surface by one or two water layers. The strength of adsorption ranges from 0 to -100 kcal/(mol peptide) and scales with the surface energy of the metal (Pd surfaces are more attractive than Au surfaces), the affinity of individual residues versus the affinity of water, and conformation aspects, as well as polarization and charge transfer at the metal interface (only qualitatively considered here). A hexagonal spacing of approximately 1.6 A between available lattice sites on the {111} surfaces accounts for the characteristic adsorption of aromatic side groups and various other residues (including Tyr, Phe, Asp, His, Arg, Asn, Ser), and a quadratic spacing of approximately 2.8 A between available lattice sites on the {100} surface accounts for a significantly lower affinity to all peptides in favor of mobile water molecules. The combination of these factors suggests a "soft epitaxy" mechanism of binding. On a bimetallic Pd-Au {111} surface, binding patterns are similar, and the polarity of the bimetal junction can modify the binding energy by approximately 10 kcal/mol. The results are semiquantitatively supported by experimental measurements of the affinity of peptides and small molecules to metal surfaces as well as results from quantum-mechanical calculations on small peptide and surface fragments. Interfaces were modeled using the consistent valence force field extended for Lennard-Jones parameters for fcc metals which accurately reproduce surface and interface energies [Heinz, H.; Vaia, R. A.; Farmer, B. L.; Naik, R. R. J. Phys. Chem. C 2008, 112, 17281-17290].

Computational study of the absorption spectra of green fluorescent protein mutants.

In this work, we present a theoretical study of the relationship between molecular structure and the red-shift in absorption spectra of S65G and S65T green fluorescent protein (GFP) mutants. To identify the effects of the protein environment, we combined results from molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics calculations to obtain structural properties, and applied time-dependent density functional theory to calculate the excitation energies. By using results from the MD simulations, we were able to provide a systematic analysis of the structural details that may effect the red-shift in the absorption spectra when taking into account temperature effects. Furthermore, a detailed study of hydrogen bonding during the MD simulations demonstrated differences between S65G and S65T, for example, regarding hydrogen bonding with Glu222. An analysis of the absorption spectra for different forms of the chromophore emphasized the dominance of the anionic forms in solution for the S65G and S65T GFP mutants.

Molecular modeling of green fluorescent protein: structural effects of chromophore deprotonation.

Molecular dynamics (MD) simulations were carried out to study the conformational rearrangement induced by deprotonation of the fluorescent chromophore in GFP, as well as the associated changes in the hydrogen-bonding network. For both the structures with either a neutral or an anionic chromophore, it was found that the beta-barrel was stable and rigid, and the conformation of the chromophore was consistent with the available x-ray structure. The conformational change in Thr203 due to deprotonation was also found to be consistent with the three-state isomerization model. Although GFP is highly fluorescent, denatured-GFP is nonfluorescent, indicating that the environment of the protein plays an important role in its fluorescence behavior. Our MD simulations, which explore the effect of the protein shell on the conformation of the chromophore, find the flexibility of the central chromophore to be significantly restricted due to the rigid nature of the protein shell. The hydrogen-bonding between the chromophore and neighboring residues was also shown to contribute to the chromophore rigidity. In addition to the MD studies, quantum mechanics/molecular mechanics (QM/MM) ONIOM calculations were carried out to investigate the effect of the beta-barrel on the internal rotation in the chromophore. Along with providing quantitative values for torsional rotation barriers about the bridging bond in the chromophore, the ONIOM calculations also validate our MD force field parameters.