In Silico Discovery and Validation of Neuropeptide-Y-Binding Peptides for Sensors.

Wearable sensors for human health, performance, and state monitoring, which have a linear response to the binding of biomarkers found in sweat, saliva, or urine, are of current interest for many applications. A critical part of any device is a biological recognition element (BRE) that is able to bind a biomarker at the surface of a sensor with a high affinity and selectivity to produce a measurable signal response. In this study, we discover and compare 12-mer peptides that bind to neuropeptide Y (NPY), a stress and human health biomarker, using independent and complimentary experimental and computational approaches. The affinities of the NPY-binding peptides discovered by both methods are equivalent and below the micromolar level, which makes them suitable for application in sensors. The in silico design protocol for peptide-based BREs is low cost, highly efficient, and simple, suggesting its utility for discovering peptide binders to a variety of biomarker targets.

Advancing Peptide-Based Biorecognition Elements for Biosensors Using in-Silico Evolution.

Sensors for human health and performance monitoring require biological recognition elements (BREs) at device interfaces for the detection of key molecular biomarkers that are measurable biological state indicators. BREs, including peptides, antibodies, and nucleic acids, bind to biomarkers in the vicinity of the sensor surface to create a signal proportional to the biomarker concentration. The discovery of BREs with the required sensitivity and selectivity to bind biomarkers at low concentrations remains a fundamental challenge. In this study, we describe an in-silico approach to evolve higher sensitivity peptide-based BREs for the detection of cardiac event marker protein troponin I (cTnI) from a previously identified BRE as the parental affinity peptide. The P2 affinity peptide, evolved using our in-silico method, was found to have ∼16-fold higher affinity compared to the parent BRE and ∼10 fM (0.23 pg/mL) limit of detection. The approach described here can be applied towards designing BREs for other biomarkers for human health monitoring.

Peptide interactions with zigzag edges in graphene.

Recognition and manipulation of graphene edges enable the control of physical properties of graphene-based devices. Recently, the authors have identified a peptide that preferentially binds to graphene edges from a combinatorial peptide library. In this study, the authors examine the functional basis for the edge binding peptide using experimental and computational methods. The effect of amino acid substitution, sequence context, and solution pH value on the binding of the peptide to graphene has been investigated. The N-terminus glutamic acid residue plays a key role in recognizing and binding to graphene edges. The protonation, substitution, and positional context of the glutamic acid residue impact graphene edge-binding. Our findings provide insights into the binding mechanisms and the design of peptides for recognizing and functionalizing graphene edges.

In silico carbon molecular beam epitaxial growth of graphene on the h-BN substrate: carbon source effect on van der Waals epitaxy.

Against the presumption that hexagonal boron-nitride (h-BN) should provide an ideal substrate for van der Waals (vdW) epitaxy to grow high quality graphene films, carbon molecular beam epitaxy (CMBE) techniques using solid carbon sublimation have reported relatively poor quality of the graphene. In this article, the CMBE growth of graphene on the h-BN substrate is numerically studied in order to identify the effect of the carbon source on the quality of the graphene film. The carbon molecular beam generated by the sublimation of solid carbon source materials such as graphite and glassy carbon is mostly composed of atomic carbon, carbon dimers and carbon trimers. Therefore, the graphene film growth becomes a complex process involving various deposition characteristics of a multitude of carbon entities. Based on the study of surface adsorption and film growth characteristics of these three major carbon entities comprising graphite vapour, we report that carbon trimers convey strong traits of vdW epitaxy prone to high quality graphene growth, while atomic carbon deposition is a surface-reaction limited process accompanied by strong chemisorption. The vdW epitaxial behaviour of carbon trimers is found to be substantial enough to nucleate and develop into graphene like planar films within a nanosecond of high flux growth simulation, while reactive atomic carbons tend to impair the structural integrity of the crystalline h-BN substrate upon deposition to form an amorphous interface between the substrate and the growing carbon film. The content of reactive atomic carbons in the molecular beam is suspected to be the primary cause of low quality graphene reported in the literature. A possible optimization of the molecular beam composition towards the synthesis of better quality graphene films is suggested.

Biotic-Abiotic Interactions: Factors that Influence Peptide-Graphene Interactions.

Understanding the factors that influence the interaction between biomolecules and abiotic surfaces is of utmost interest in biosensing and biomedical research. Through phage display technology, several peptides have been identified as specific binders to abiotic material surfaces, such as gold, graphene, silver, and so forth. Using graphene-peptide as our model abiotic-biotic pair, we investigate the effect of graphene quality, number of layers, and the underlying support substrate effect on graphene-peptide interactions using both experiments and computation. Our results indicate that graphene quality plays a significant role in graphene-peptide interactions. The graphene-biomolecule interaction appears to show no significant dependency on the number of graphene layers or the underlying support substrate.

Thermal anisotropy in nano-crystalline MoS2 thin films.

In this work, we grow thin MoS2 films (50-150 nm) uniformly over large areas (>1 cm(2)) with strong basal plane (002) or edge plane (100) orientations to characterize thermal anisotropy. Measurement results are correlated with molecular dynamics simulations of thermal transport for perfect and defective MoS2 crystals. The correlation between predicted (simulations) and measured (experimental) thermal conductivity are attributed to factors such as crystalline domain orientation and size, thereby demonstrating the importance of thermal boundary scattering in limiting thermal conductivity in nano-crystalline MoS2 thin films. Furthermore, we demonstrate that the cross-plane thermal conductivity of the films is strongly impacted by exposure to ambient humidity.

Conformational response to solvent interaction and temperature of a protein (Histone h3.1) by a multi-grained monte carlo simulation.

Interaction with the solvent plays a critical role in modulating the structure and dynamics of a protein. Because of the heterogeneity of the interaction strength, it is difficult to identify multi-scale structural response. Using a coarse-grained Monte Carlo approach, we study the structure and dynamics of a protein (H3.1) in effective solvent media. The structural response is examined as a function of the solvent-residue interaction strength (based on hydropathy index) in a range of temperatures (spanning low to high) involving a knowledge-based (Miyazawa-Jernigan(MJ)) residue-residue interaction. The protein relaxes rapidly from an initial random configuration into a quasi-static structure at low temperatures while it continues to diffuse at high temperatures with fluctuating conformation. The radius of gyration (Rg ) of the protein responds non-monotonically to solvent interaction, i.e., on increasing the residue-solvent interaction strength (fs ), the increase in Rg (fs ≤fsc ) is followed by decay (fs ≥fsc ) with a maximum at a characteristic value (fsc ) of the interaction. Raising the temperature leads to wider spread of the distribution of the radius of gyration with higher magnitude of fsc . The effect of solvent on the multi-scale (λ: residue to Rg ) structures of the protein is examined by analyzing the structure factor (S( q ),|q| = 2π/λ is the wave vector of wavelength, λ) in detail. Random-coil to globular transition with temperature of unsolvated protein (H3.1) is dramatically altered by the solvent at low temperature while a systematic change in structure and scale is observed on increasing the temperature. The interaction energy profile of the residues is not sufficient to predict its mobility in the solvent. Fine-grain representation of protein with two-node and three-node residue enhances the structural resolution; results of the fine-grained simulations are consistent with the finding described above of the coarse-grained description with one-node residue.

A hierarchical coarse-grained (all-atom-to-all-residue) computer simulation approach: self-assembly of peptides.

A hierarchical computational approach (all-atom residue to all-residue peptide) is introduced to study self-organizing structures of peptides as a function of temperature. A simulated residue-residue interaction involving all-atom description, analogous to knowledge-based analysis (with different input), is used as an input to a phenomenological coarse-grained interaction for large scales computer simulations. A set of short peptides P1 ((1)H (2)S (3)S (4)Y (5)W (6)Y (7)A (8)F (9)N (10)N (11)K (12)T) is considered as an example to illustrate the utility. We find that peptides assemble rather fast into globular aggregates at low temperatures and disperse as random-coil at high temperatures. The specificity of the mass distribution of the self-assembly depends on the temperature and spatial lengths which are identified from the scaling of the structure factor. Analysis of energy and mobility profiles, gyration radius of peptide, and radial distribution function of the assembly provide insight into the multi-scale (intra- and inter-chain) characteristics. Thermal response of the global assembly with the simulated residue-residue interaction is consistent with that of the knowledge-based analysis despite expected quantitative differences.

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.

Variation in structure of a protein (H2AX) with knowledge-based interactions.

The structure of a protein (H2AX) as a function of temperature is examined by three knowledge-based phenomenological interactions, MJ (Miyazawa and Jernigan), BT (Betancourt and Thirumalai), and BFKV (Bastolla et al.) to identify similarities and differences in results. Data from the BT and BFKV residue-residue interactions verify finding with the MJ interaction, i.e., the radius of gyration (Rg ) of H2AX depends non-monotonically on temperature. The increase in Rg is followed by a decay on raising the temperature with a maximum at a characteristic value, Tc , which depends on the knowledge-based contact matrix, TcBFKV ≤ TcMJ ≤ TcBT . The range (ΔT) of non-monotonic thermal response and its decay pattern with the temperature are sensitive to interaction. A rather narrow temperature range of ΔTMJ ≈ 0.015-0.022 with the MJ interaction expands and shifts up to ΔTBT ≈ 0.018-0.30 at higher temperatures with the BT interaction and shifts down with the BFKV interaction to ΔTBFKV ≈ 0.011-0.018. The scaling of the structure factor with the wave vector reveals that the structure of the protein undergoes a transformation from a random coil at high temperature to a globular conformation at low temperatures.

The effect of single wall carbon nanotube metallicity on genomic DNA-mediated chirality enrichment.

Achieving highly enriched single wall carbon nanotubes (SWNTs) is one of the major hurdles today because their chirality-dependent properties must be uniform and predictable for use in nanoscale electronics. Due to the unique wrapping and groove-binding mechanism, DNA has been demonstrated as a highly specific SWNT dispersion and fractionation agent, with its enrichment capabilities depending on the DNA sequence and length as well as the nanotube properties. Salmon genomic DNA (SaDNA) offers an inexpensive and scalable alternative to synthetic DNA. In this study, SaDNA enrichment capabilities were tested on SWNT separation with varying degrees of metallicity that were formulated from mixtures of commercial metallic (met-) and semiconducting (sem-) abundant SWNTs. The results herein demonstrate that the degree of metallicity of the SWNT sample has a significant effect on the SaDNA enrichment capabilities, and this effect is modeled based on deconvolution of the near-infrared (NIR) absorption spectra and verified with photoluminescence emission (PLE) measurements. Using molecular dynamics and circular dichroism, the preferential SaDNA mediated separation of the (6, 5) sem-tube is shown to be largely influenced by the presence of met-SWNTs.

Random coil to globular thermal response of a protein (H3.1) with three knowledge-based coarse-grained potentials.

The effect of temperature on the conformation of a histone (H3.1) is studied by a coarse-grained Monte Carlo simulation based on three knowledge-based contact potentials (MJ, BT, BFKV). Despite unique energy and mobility profiles of its residues, the histone H3.1 undergoes a systematic (possibly continuous) structural transition from a random coil to a globular conformation on reducing the temperature. The range over which such a systematic response in variation of the radius of gyration (R(g)) with the temperature (T) occurs, however, depends on the potential, i.e. ΔT(MJ) ≈ 0.013-0.020, ΔT(BT) ≈ 0.018-0.026, and ΔT(BFKV) ≈ 0.006-0.013 (in reduced unit). Unlike MJ and BT potentials, results from the BFKV potential show an anomaly where the magnitude of R(g) decreases on raising the temperature in a range ΔT(A) ≈ 0.015-0.018 before reaching its steady-state random coil configuration. Scaling of the structure factor, S(q) ∝ q(-1/ν), with the wave vector, q=2π/λ, and the wavelength, λ, reveals a systematic change in the effective dimension (D(e)∼1/ν) of the histone with all potentials (MJ, BT, BFKV): D(e)∼3 in the globular structure with D(e)∼2 for the random coil. Reproducibility of the general yet unique (monotonic) structural transition of the protein H3.1 with the temperature (in contrast to non-monotonic structural response of a similar but different protein H2AX) with three interaction sets shows that the knowledge-based contact potential is viable tool to investigate structural response of proteins. Caution should be exercise with the quantitative comparisons due to differences in transition regimes with these interactions.

A novel nano-configuration for thermoelectrics: helicity induced thermal conductivity reduction in nanowires.

In this article, we propose a novel helical nano-configuration towards the designing of high ZT thermoelectric materials. Non-equilibrium molecular dynamics (NEMD) simulations for 'model' bi-component nanowires indicate that a significant reduction in thermal conductivity, similar to that of flat superlattice nanostructures, can be achieved using a helical geometric configuration. The reduction is attributed to a plethora of transmissive and reflective phonon scattering events resulting from the steady alteration of phonon propagating direction that emerges from the continuous rotation of the helical interface. We also show that increasing the relative mass ratio of the two components lowers the phonon energy transmission at the interface due to differences in vibrational frequency spectra, thereby relatively 'easing' the phonon energy propagation along the helical pathway. While the proposed mechanisms result in a reduced lattice thermal conductivity, the continuous nature of the bi-component nanowire would not be expected to significantly reduce its electrical counterpart, as often occurs in superlattice/alloy nanostructures. Hence, we postulate that the helical configuration of atomic arrangement provides an attractive and general framework for improved thermoelectric material assemblies independent of the specific chemical composition.

Thermal rectification in three-dimensional asymmetric nanostructure.

Previously, thermal rectification has been reported in several low-dimensional shape-asymmetric nanomaterials. In this Letter, we demonstrate that a three-dimensional crystalline material with an asymmetric shape also displays as strong thermal rectification as low-dimensional materials do. The observed rectification is attributed to the stronger temperature dependence of vibration density of states in the narrower region of the asymmetric material, resulting from the small number of atomic degrees of freedom directly interacting with the thermostat. We also demonstrate that the often reported "device shape asymmetry" is not a sufficient condition for thermal rectification. Specifically, the size asymmetry in boundary thermal contacts is equally important toward determining the magnitude of thermal rectification. When the boundary thermal contacts retain the same size asymmetry as the nanomaterial, the overall system displays notable thermal rectification, in accordance with existing literature. However, when the wider region of the asymmetric nanomaterial is partially thermostatted by a smaller sized contact, thermal rectification decreases dramatically and even changes direction.

Structure of a peptide adsorbed on graphene and graphite.

Noncovalent functionalization of graphene using peptides is a promising method for producing novel sensors with high sensitivity and selectivity. Here we perform atomic force microscopy, Raman spectroscopy, infrared spectroscopy, and molecular dynamics simulations to investigate peptide-binding behavior to graphene and graphite. We studied a dodecamer peptide identified with phage display to possess affinity for graphite. Optical spectroscopy reveals that the peptide forms secondary structures both in powder form and in an aqueous medium. The dominant structure in the powder form is α-helix, which undergoes a transition to a distorted helical structure in aqueous solution. The peptide forms a complex reticular structure upon adsorption on graphene and graphite, having a helical conformation different from α-helix due to its interaction with the surface. Our observation is consistent with our molecular dynamics calculations, and our study paves the way for rational functionalization of graphene using biomolecules with defined structures and, therefore, functionalities.

Conformational temperature-dependent behavior of a histone H2AX: a coarse-grained Monte Carlo approach via knowledge-based interaction potentials.

Histone proteins are not only important due to their vital role in cellular processes such as DNA compaction, replication and repair but also show intriguing structural properties that might be exploited for bioengineering purposes such as the development of nano-materials. Based on their biological and technological implications, it is interesting to investigate the structural properties of proteins as a function of temperature. In this work, we study the spatial response dynamics of the histone H2AX, consisting of 143 residues, by a coarse-grained bond fluctuating model for a broad range of normalized temperatures. A knowledge-based interaction matrix is used as input for the residue-residue Lennard-Jones potential.We find a variety of equilibrium structures including global globular configurations at low normalized temperature (T* = 0.014), combination of segmental globules and elongated chains (T* = 0.016,0.017), predominantly elongated chains (T* = 0.019,0.020), as well as universal SAW conformations at high normalized temperature (T* ≥ 0.023). The radius of gyration of the protein exhibits a non-monotonic temperature dependence with a maximum at a characteristic temperature (T(c)* = 0.019) where a crossover occurs from a positive (stretching at T* ≤ T(c)*) to negative (contraction at T* ≥ T(c)*) thermal response on increasing T*.

Influence of the shape of nanostructured metal surfaces on adsorption of single peptide molecules in aqueous solution.

Self-assembly and function of biologically modified metal nanostructures depend on surface-selective adsorption; however, the influence of the shape of metal surfaces on peptide adsorption mechanisms has been poorly understood. The adsorption of single peptide molecules in aqueous solution (Tyr(12) , Ser(12) , A3, Flg-Na(3) ) is investigated on even {111} surfaces, stepped surfaces, and a 2 nm cuboctahedral nanoparticle of gold using molecular dynamics simulation with the CHARMM-METAL force field. Strong and selective adsorption is found on even surfaces and the inner edges of stepped surfaces (-20 to -60 kcal/mol peptide) in contrast to weaker and less selective adsorption on small nanoparticles (-15 to -25 kcal/mol peptide). Binding and selectivity appear to be controlled by the size of surface features and the extent of co-ordination of epitaxial sites by polarizable atoms (N, O, C) along the peptide chain. The adsorption energy of a single peptide equals a fraction of the sum of the adsorption energies of individual amino acids that is characteristic of surface shape, epitaxial pattern, and conformation constraints (often ß-strand and random coil). The proposed adsorption mechanism is supported and critically evaluated by earlier sequence data from phage display, dissociation constants of small proteins as a function of nanoparticle size, and observed shapes of peptide-stabilized nanoparticles. Understanding the interaction of single peptides with shaped metal surfaces is a key step towards control over self-organization of multiple peptides on shaped metal surfaces and the assembly of superstructures from nanostructures.

Importance of interfaces in governing thermal transport in composite materials: modeling and experimental perspectives.

Thermal management in polymeric composite materials has become increasingly critical in the air-vehicle industry because of the increasing thermal load in small-scale composite devices extensively used in electronics and aerospace systems. The thermal transport phenomenon in these small-scale heterogeneous systems is essentially controlled by the interface thermal resistance because of the large surface-to-volume ratio. In this review article, several modeling strategies are discussed for different length scales, complemented by our experimental efforts to tailor the thermal transport properties of polymeric composite materials. Progress in the molecular modeling of thermal transport in thermosets is reviewed along with a discussion on the interface thermal resistance between functionalized carbon nanotube and epoxy resin systems. For the thermal transport in fiber-reinforced composites, various micromechanics-based analytical and numerical modeling schemes are reviewed in predicting the transverse thermal conductivity. Numerical schemes used to realize and scale the interface thermal resistance and the finite mean free path of the energy carrier in the mesoscale are discussed in the frame of the lattice Boltzmann-Peierls-Callaway equation. Finally, guided by modeling, complementary experimental efforts are discussed for exfoliated graphite and vertically aligned nanotubes based composites toward improving their effective thermal conductivity by tailoring interface thermal resistance.

Single mode phonon energy transmission in functionalized carbon nanotubes.

Although the carbon nanotube (CNT) features superior thermal properties in its pristine form, the chemical functionalization often required for many applications of CNT inevitably degrades the structural integrity and affects the transport of energy carriers. In this article, the effect of the side wall functionalization on the phonon energy transmission along the symmetry axis of CNT is studied using the phonon wave packet method. Three different functional groups are studied: methyl (-CH(3)), vinyl (-C(2)H(3)), and carboxyl (-COOH). We find that, near Γ point of the Brillouin zone, acoustic phonons show ideal transmission, while the transmission of the optical phonons is strongly suppressed. A positive correlation between the energy transmission coefficient and the phonon group velocity is observed for both acoustic and optical phonon modes. On comparing the transmission due to functional groups with equivalent point mass defects on CNT, we find that the chemistry of the functional group, rather than its molecular mass, has a dominant role in determining phonon scattering, hence the transmission, at the defect sites.

Preferential binding of peptides to graphene edges and planes.

Peptides identified from combinatorial peptide libraries have been shown to bind to a variety of abiotic surfaces. Biotic-abiotic interactions can be exploited to create hybrid materials with interesting electronic, optical, or catalytic properties. Here we show that peptides identified from a combinatorial phage display peptide library assemble preferentially to the edge or planar surface of graphene and can affect the electronic properties of graphene. Molecular dynamics simulations and experiments provide insight into the mechanism of peptide binding to the graphene edge.