Electronic properties of a graphene device with peptide adsorption: insight from simulation.

In this work, to explain doping behavior of single-layer graphene upon HSSYWYAFNNKT (P1) and HSSAAAAFNNKT (P1-3A) adsorption in field-effect transistors (GFETs), we applied a combined computational approach, whereby peptide adsorption was modeled by molecular dynamics simulations, and the lowest energy configuration was confirmed by density functional theory calculations. On the basis of the resulting structures of the hybrid materials, electronic structure and transport calculations were investigated. We demonstrate that π-π stacking of the aromatic residues and proximate peptide backbone to the graphene surface in P1 have a role in the p-doping. These results are consistent with our experimental observation of the GFET's p-doping even after a 24-h annealing procedure. Upon substitution of three of the aromatic residues to Ala in (P1-3A), a considerable decrease from p-doping is observed experimentally, demonstrating n-doping as compared to the nonadsorbed device, yet not explained based on the atomistic MD simulation structures. To gain a qualitative understanding of P1-3A's adsorption over a longer simulation time, which may differ from aromatic amino acid residues' swift anchoring on the surface, we analyzed equilibrated coarse-grain simulations performed for 500 ns. Desorption of the Ala residues from the surface was shown computationally, which could in turn affect charge transfer, yet a full explanation of the mechanism of n-doping will require elucidation of differences between various aromatic residues as dependent on peptide composition, and inclusion of effects of the substrate and environment, to be considered in future work.

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.

Halogenated aliphatic toxicity QSARs employing metabolite descriptors.

The toxic effects from exposure to halogenated hydrocarbons (HAs), which are produced in large amounts and used in a variety of applications, are well-known. Previously, QSARs for the toxicity of a series of HAs in vitro have been studied extensively. In this work, using a composite toxicity metric calculated from a set of five in vitro hepatotoxicity endpoints determined for 20 HAs, we find that QSARs derived using quantum descriptors calculated from the neutral HA species are statistically similar to QSARs calculated from HA metabolites. In most cases, QSARs derived using descriptors calculated from both neutral HAs and metabolites are statistically superior to those derived using either neutral-HA descriptors or metabolite descriptors. However, to properly utilize metabolite descriptors, multiple QSARs, each of which utilizes a set of HAs that form unique metabolites, must be derived and toxicity values calculated therefrom must be averaged. These average toxicity values agree better with experiment than those calculated from the neutral-HA QSARs.

Risk assessment of high-energy chemicals by in vitro toxicity screening and quantitative structure-activity relationships.

Hydrazine propellants pose a substantial operational concern to the U.S. Air Force and to the aerospace industry because of their toxicity. In our continuing efforts to develop methods for the prediction of the toxicological response to such materials, we have measured in vitro toxicity endpoints for a series of high-energy chemicals (HECs) that were recently proposed as propellants. The HECs considered are structurally diverse and can be classified into four chemical types (hydrazine-based, amino-based, triazoles, and a quaternary ammonium salt), although most are hydrazine derivatives. We measured the following endpoints in primary cultures of isolated rat hepatocytes: mitochondrial function (MTT), lactate dehydrogenase leakage (LDH), generation of reactive oxygen species (ROS), and total glutathione content (GSH). In several instances, effective concentrations (EC) were indeterminate, and only lower limits to the measured endpoints could be ascertained. Using molecular descriptors calculated with a semiempirical molecular orbital method, quantitative structure-activity relationships (QSARs) were derived for MTT (EC25) and for GSH (EC50). Correlation coefficients for 2- and 3-parameter QSARs of about 0.9 enable us to predict the toxicity for similar compounds. Furthermore, except in one case, predicted EC values for the uncertain endpoints were consistent with experiment. Descriptors comprising the QSARs for MTT were consistent with the biophysical mechanism of toxic response found experimentally for hydrazine derivatives. Application of our derived QSARs will assist in predicting toxicity for newly proposed propellants.