Discontinuous Molecular Dynamics Simulations of Biomolecule Interfacial Behavior: Study of Ovispirin-1 Adsorption on a Graphene Surface

Graph Clustering Analyses of Discontinuous Molecular Dynamics Simulations: Study of Lysozyme Adsorption on a Graphene Surface

Understanding the interfacial behaviors of biomolecules is crucial to applications in biomaterials and nanoparticle-based biosensing technologies. In this work, we utilized autoencoder-based graph clustering to analyze discontinuous molecular dynamics (DMD) simulations of lysozyme adsorption on a graphene surface. Our high-throughput DMD simulations integrated with a Go̅-like protein-surface interaction model makes it possible to explore protein adsorption at a large temporal scale with sufficient accuracy. The graph autoencoder extracts a low-dimensional feature vector from a contact map. The sequence of the extracted feature vectors is then clustered, and thus the evolution of the protein molecule structure in the absorption process is segmented into stages. Our study demonstrated that the residue-surface hydrophobic interactions and the π-π stacking interactions play key roles in the five-stage adsorption. Upon adsorption, the tertiary structure of lysozyme collapsed, and the secondary structure was also affected. The folding stages obtained by autoencoder-based graph clustering were consistent with detailed analyses of the protein structure. The combination of machine learning analysis and efficient DMD simulations developed in this work could be an important tool to study biomolecules' interfacial behaviors.

Understanding interactions between biomolecules and two-dimensional nanomaterials using in silico microscopes

Two-dimensional (2D) nanomaterials such as graphene are increasingly used in research and industry for various biomedical applications. Extensive experimental and theoretical studies have revealed that 2D nanomaterials are promising drug delivery vehicles, yet certain materials exhibit toxicity under biological conditions. So far, it is known that 2D nanomaterials possess strong adsorption propensities for biomolecules. To mitigate potential toxicity and retain favorable physical and chemical properties of 2D nanomaterials, it is necessary to explore the underlying mechanisms of interactions between biomolecules and nanomaterials for the subsequent design of biocompatible 2D nanomaterials for nanomedicine. The purpose of this review is to integrate experimental findings with theoretical observations and facilitate the study of 2D nanomaterial interaction with biomolecules at the molecular level. We discuss the current understanding and progress of 2D nanomaterial interaction with proteins, lipid membranes, and DNA based on molecular dynamics (MD) simulation. In this review, we focus on the 2D graphene nanosheet and briefly discuss other 2D nanomaterials. With the ever-growing computing power, we can image nanoscale processes using MD simulation that are otherwise not observable in experiment. We expect that molecular characterization of the complex behavior between 2D nanomaterials and biomolecules will help fulfill the goal of designing effective 2D nanomaterials as drug delivery platforms.

Multiscale Simulation of Protein Corona Formation on Silver Nanoparticles: Study of Ovispirin-1 Peptide Adsorption

The exposure of nanoparticles (NPs) to biofluids leads to the rapid coverage of proteins, named protein corona, which alters the NPs' chemicophysical and biological properties. Fundamental studies of the protein corona are thus critical to the increasing applications of NPs in nanotechnology and nanomedicines. The present work utilizes multiscale simulations of a model biological system, small ovispirin-1 peptides, and bare silver nanoparticles (AgNPs) to examine the NPs' size and surface hydrophilicity effects on formation dynamics and the structure of the peptide corona. Our simulations revealed the different adsorption dynamics of ovispirin-1 peptides on the NPs, including the direct adsorption of a single peptide and peptide aggregates and multistep adsorption, as well as an intermediate cycle of desorption and readsorption. Notably, the whole process of peptide adsorption on hydrophilic AgNP surfaces can be generalized as three stages: diffusion to the surface, initial landing via hydrophilic residues, and the final attachment. The decrease in AgNP's size leads to faster adsorption with more heterogeneous peptide interfacial dynamics, a denser and inhomogeneous peptide packing structure, and a wider distribution of adsorption orientations. Subsequent atomistic molecular dynamics simulations demonstrated that on the hydrophilic AgNP surfaces, adsorbed peptides display moderate changes in their secondary structure, resulting in further changes of corona composition, i.e., amino acid residue distribution on the surface.

Molecular Structure of the Surface-Immobilized Super Uranyl Binding Protein

Recently, a super uranyl binding protein (SUP) was developed, which exhibits excellent sensitivity/selectivity to bind uranyl ions. It can be immobilized onto a surface in sensing devices to detect uranyl ions. Here, sum frequency generation (SFG) vibrational spectroscopy was applied to probe the interfacial structures of surface-immobilized SUP. The collected SFG spectra were compared to the calculated orientation-dependent SUP SFG spectra using a one-excitonic Hamiltonian approach based on the SUP crystal structures to deduce the most likely surface-immobilized SUP orientation(s). Furthermore, discrete molecular dynamics (DMD) simulation was applied to refine the surface-immobilized SUP conformations and orientations. The immobilized SUP structures calculated from DMD simulations confirmed the SUP orientations obtained from SFG data analyzed based on the crystal structures and were then used for a new round of SFG orientation analysis to more accurately determine the interfacial orientations and conformations of immobilized SUP before and after uranyl ion binding, providing an in-depth understanding of molecular interactions between SUP and the surface and the effect of uranyl ion binding on the SUP interfacial structures. We believe that the developed method of combining SFG measurements, DMD simulation, and Hamiltonian data analysis approach is widely applicable to study biomolecules at solid/liquid interfaces.