http://pubs.acs.org/doi/abs/10.1021/acsami.5b05615 Studies reveal that biomolecules can form intriguing molecular structures with fascinating functionalities upon interaction with graphene. Then, interesting questions arise. How does silk fibroin interact with graphene? Does such interaction lead to an enhancement in its mechanical properties? In this study, using large-scale molecular dynamics simulations, we first examine the interaction of graphene with several typical peptide structures of silk fibroin extracted from different domains of silk fibroin, including pure amorphous (P1), pure crystalline (P2), a segment from N-terminal (P3), and a combined amorphous and crystalline segment (P4), aiming to reveal their structural modifications. Our study shows that graphene can have intriguing influences on the structures formed by the peptides with sequences representing different domains of silk fibroin. In general, for protein domains with stable structure and strong intramolecular interaction (e.g., β-sheets), graphene tends to compete with the intramolecular interactions and thus weaken the interchain interaction and reduce the contents of β-sheets. For the silk domains with random or less ordered secondary structures and weak intramolecular interactions, graphene tends to enhance the stability of peptide structures; in particular, it increases the contents of helical structures. Thereafter, tensile simulations were further performed on the representative peptides to investigate how such structure modifications affect their mechanical properties. It was found that the strength and resilience of the peptides are enhanced through their interaction with graphene. The present work reveals interesting insights into the interactions between silk peptides and graphene, and contributes in the efforts to enhance the mechanical properties of silk fibroin.

We consider the mesoscopic model of protein materials composed of protein crystals with given space group for understanding the mechanical properties of protein materials with respect to their structures. This preprint was accepted for publication at Journal of Computational Chemistry.

Abstract
Mechanical characterization of protein molecules has played a role on gaining insight into the biological functions of proteins, since some proteins perform the mechanical function. Here, we present the mesoscopic model of biological protein materials composed of protein crystals prescribed by Go potential for characterization of elastic behavior of protein materials. Specifically, we consider the representative volume element (RVE) containing the protein crystals represented by alpha-carbon atoms, prescribed by Go potential, with application of constant normal strain to RVE. The stress-strain relationship computed from virial stress theory provides the nonlinear elastic behavior of protein materials and their mechanical properties such as Young's modulus, quantitatively and/or qualitatively comparable to mechanical properties of biological protein materials obtained from experiments and/or atomistic simulations. Further, we discuss the role of native topology on the mechanical properties of protein crystals. It is shown that parallel strands (hydrogen bonds in parallel) enhances the mechanical resilience of protein materials.