Force fields for simulating the interaction of surfaces with biological molecules

Iron oxide magnetic aggregates: Aspects of synthesis, computational approaches and applications

Superparamagnetic magnetite nanoparticles have been central to numerous investigations in the past few decades for their use in many applications, such as drug delivery, medical diagnostics, magnetic separation, and material science. However, the properties of single magnetic nanoparticles are sometimes not sufficient to accomplish tasks where a strong magnetic response is required. In light of this, aggregated magnetite nanoparticles have been proposed as an alternative advanced material, which may expand and combine some of the advantages of single magnetic nanoparticles, including superparamagnetism, with an enhanced magnetic moment and increased colloidal stability. This review comprehensively discusses the current literature on aggregates made of magnetic iron oxide nanoparticles. This review is divided into three sections. First, the current synthetic strategies for magnetite nanoparticle aggregates are discussed, together with the influence of different stabilizers on the primary crystals and the final aggregate size and morphology. The second section is dedicated to computational approaches, such as density functional methods (which permit accurate predictions of electronic and magnetic properties and shed light on the behavior of surfactant molecules on iron oxide surfaces) and molecular dynamics simulations (which provide additional insight into the influence of ligands on the surface chemistry of iron oxide nanocrystals). The last section discusses current and possible future applications of iron oxide magnetic aggregates, including wastewater treatment, water purification, medical applications, and magnetic aggregates for materials displaying structural colors.

Computational Insights into the Adsorption of Ligands on Gold Nanosurfaces

We study the adsorption process of model peptides, nucleobases, and selected standard ligands on gold through the development of a computational protocol based on fully atomistic classical molecular dynamics (MD) simulations combined with umbrella sampling techniques. The specific features of the interface components, namely, the molecule, the metallic substrate, and the solvent, are taken into account through different combinations of force fields (FFs), which are found to strongly affect the results, especially changing absolute and relative adsorption free energies and trends. Overall, noncovalent interactions drive the process along the adsorption pathways. Our findings also show that a suitable choice of the FF combinations can shed light on the affinity, position, orientation, and dynamic fluctuations of the target molecule with respect to the surface. The proposed protocol may help the understanding of the adsorption process at the microscopic level and may drive the in-silico design of biosensors for detection purposes.

An ancestral molecular response to nanomaterial particulates

The varied transcriptomic response to nanoparticles has hampered the understanding of the mechanism of action. Here, by performing a meta-analysis of a large collection of transcriptomics data from various engineered nanoparticle exposure studies, we identify common patterns of gene regulation that impact the transcriptomic response. Analysis identifies deregulation of immune functions as a prominent response across different exposure studies. Looking at the promoter regions of these genes, a set of binding sites for zinc finger transcription factors C2H2, involved in cell stress responses, protein misfolding and chromatin remodelling and immunomodulation, is identified. The model can be used to explain the outcomes of mechanism of action and is observed across a range of species indicating this is a conserved part of the innate immune system.

Multiscale Modeling of Bio-Nano Interactions of Zero-Valent Silver Nanoparticles

Understanding the specifics of interaction between the protein and nanomaterial is crucial for designing efficient, safe, and selective nanoplatforms, such as biosensor or nanocarrier systems. Routing experimental screening for the most suitable complementary pair of biomolecule and nanomaterial used in such nanoplatforms might be a resource-intensive task. While a range of computational tools are available for prescreening libraries of proteins for their interactions with small molecular ligands, choices for high-throughput screening of protein libraries for binding affinities to new and existing nanomaterials are very limited. In the current work, we present the results of the systematic computational study of interaction of various biomolecules with pristine zero-valent noble metal nanoparticles, namely, AgNPs, by using the UnitedAtom multiscale approach. A set of blood plasma and dietary proteins for which the interaction with AgNPs was described experimentally were examined computationally to evaluate the performance of the UnitedAtom method. A set of interfacial descriptors (log P^NM, adsorption affinities, and adsorption affinity ranking), which can characterize the relative hydrophobicity/hydrophilicity/lipophilicity of the nanosized silver and its ability to form bio(eco)corona, was evaluated for future use in nano-QSAR/QSPR studies.

Adsorption of amino acids on graphene: assessment of current force fields

We compare the free energies of adsorption (ΔAads) and the structural preferences of amino acids on graphene obtained using the non-polarizable force fields-Amberff99SB-ILDN/TIP3P, CHARMM36/modified-TIP3P, OPLS-AA/M/TIP3P, and Amber03w/TIP4P/2005. The amino acid-graphene interactions are favorable irrespective of the force field. While the magnitudes of ΔAads differ between the force fields, the relative free energy of adsorption across amino acids is similar for the studied force fields. ΔAads positively correlates with amino acid-graphene and negatively correlates with graphene-water interaction energies. Using a combination of principal component analysis and density-based clustering technique, we grouped the structures observed in the graphene adsorbed state. The resulting population of clusters, and the conformation in each cluster indicate that the structures of the amino acid in the graphene adsorbed state vary across force fields. The differences in the conformations of amino acids are more severe in the graphene adsorbed state compared to the bulk state for all the force fields. Our findings suggest that the force fields studied will give qualitatively consistent relative strength of adsorption across proteins but different structural preferences in the graphene adsorbed state.

Poly(vinyl alcohol) as a water protecting agent for silver nanoparticles: the role of polymer size and structure

Chemical modification of silver nanoparticles (AgNPs) with a stabilizing agent, such as poly(vinyl alcohol) (PVA), plays an important role in shape-controlled seeded-growth and colloidal stability. However, theoretical aspects of the stabilizing mechanism of PVA are still poorly understood. To gain a better understanding of the role of PVA in water protecting effects for silver nanoparticles, we developed an atomistic model of a AgNP grafted with single-chain PVA of various lengths. Our model, designed for classical molecular dynamics (MD) simulations, approximates the AgNP as a quasi-spherical silver nanocrystal with 3.9 nm diameter and uses a united-atom representation for PVA with its polymer chain length varying from 220 up to 1540 repeating units. We found that PVA adsorbs onto the AgNP surface through multiple non-covalent interactions, among which non-covalent bonding of the hydroxyl groups plays a key role. The analysis of adsorption isotherms by using the Hill, Scatchard, and McGhee & von Hippel models exhibits evidence for positive binding cooperativity with the cooperativity parameter varying from 1.55 to 2.12. Our results indicate that the size of the PVA polymer rather than its structure plays a crucial role in providing water protecting effects for the AgNP core, varying from 40% up to 91%. The water-protecting efficiency was well approximated by the Langmuir-Freundlich equation, allowing us to predict that the saturated coverage of the nanoparticle of a given diameter of 3.9 nm should occur when the PVA molecular weight approaches 115 kDa, which corresponds to the number of vinyl alcohol monomers being equal to 3100 units.

Interplay between adsorbed peptide structure, trapped water, and surface hydrophobicity

Atomistic molecular dynamics simulations were used to study the influence of interfacial water on the orientation and conformation of a facewise amphipathic α-helical peptide adsorbed to hydrophilic and hydrophobic substrates. Water behavior beneath the peptide adsorbed to a hydrophilic surface was observed to vary with the height of the peptide above the surface. In general, the orientation of water close to the peptide (with the oxygen atom pointing up toward the peptide) was complementary to that observed near the hydrophilic surface in the absence of peptide. That is, no change in orientation of water trapped between the peptide and a hydrophilic surface is required as the peptide approaches the surface. The adsorption of the peptide to the hydrophilic surface was observed to be mediated by a layer of ordered water. Water was found to be largely excluded on adsorption to the hydrophobic surface. However, the small amount of water present was observed to be highly ordered. At the closest point of contact to the hydrophobic surface, the peptide was observed to make direct contact. These findings shed light on the fundamental driving forces of peptide adsorption to hydrophobic and hydrophilic surfaces in aqueous environments.

POP1 might be recruiting its type-Ia interface for NLRP3-mediated PYD-PYD interaction: Insights from MD simulation

Inflammasomes are multiprotein caspase-activating complexes that enhance the maturation and release of proinflammatory cytokines (IL-1β and IL-18) in response to the invading pathogen and/or host-derived cellular stress. These are assembled by the sensory proteins (viz NLRC4, NLRP1, NLRP3, and AIM-2), adaptor protein (ASC), and effector molecule procaspase-1. In NLRP3-mediated inflammasome activation, ASC acts as a mediator between NLRP3 and procaspase-1 for the transmission of signals. A series of homotypic protein-protein interactions (NLRP3^PYD :ASC^PYD and ASC^CARD :CASP1^CARD ) propagates the downstream signaling for the production of proinflammatory cytokines. Pyrin-only protein 1 (POP1) is known to act as the regulator of inflammasome. It modulates the ASC-mediated inflammasome assembly by interacting with pyrin domain (PYD) of ASC. However, despite similar electrostatic surface potential, the interaction of POP1 with NLRP3^PYD is obscured till date. Herein, to explore the possible PYD-PYD interactions between NLRP3^PYD and POP1, a combined approach of protein-protein docking and molecular dynamics simulation was adapted. The current study revealed that POP1's type-Ia interface and type-Ib interface of NLRP3^PYD might be crucial for 1:1 PYD-PYD interaction. In addition to type-I mode of interaction, we also observed type-II and type-III interaction modes in two different dynamically stable heterotrimeric complexes (POP1-NLRP3-NLRP3 and POP1-NLRP3-POP1). The inter-residual/atomic distance calculation exposed several critical residues that possibly govern the said interaction, which need further investigation. Overall, the findings of this study will shed new light on hitherto concealed molecular mechanisms underlying NLRP3-mediated inflammasome, which will have strong future therapeutic implications.

Atomistic Modelling of Confined Polypropylene Chains between Ferric Oxide Substrates at Melt Temperature

The interactions and conformational characteristics of confined molten polypropylene (PP) chains between ferric oxide (Fe₂O₃) substrates were investigated by molecular dynamics (MD) simulations. A comparative analysis of the adsorbed amount shows strong adsorption of the chains on the high-energy surface of Fe₂O₃. Local structures formed in the polymer film were studied utilizing density profiles, orientation of bonds, and end-to-end distance of chains. At interfacial regions, the backbone carbon-carbon bonds of the chains preferably orient in the direction parallel to the surface while the carbon-carbon bonds with the side groups show a slight tendency to orient normal to the surface. Based on the conformation tensor data, the chains are compressed in the normal direction to the substrates in the interfacial regions while they tend to flatten in parallel planes with respect to the surfaces. The orientation of the bonds as well as the overall flattening of the chains in planes parallel to the solid surfaces are almost identical to that of the unconfined PP chains. Also, the local pressure tensor is anisotropic closer to the solid surfaces of Fe₂O₃ indicating the influence of the confinement on the buildup imbalance of normal and tangential pressures.

Integrated multiscale biomaterials experiment and modelling: a perspective

Advances in multiscale models and computational power have enabled a broad toolset to predict how molecules, cells, tissues and organs behave and develop. A key theme in biological systems is the emergence of macroscale behaviour from collective behaviours across a range of length and timescales, and a key element of these models is therefore hierarchical simulation. However, this predictive capacity has far outstripped our ability to validate predictions experimentally, particularly when multiple hierarchical levels are involved. The state of the art represents careful integration of multiscale experiment and modelling, and yields not only validation, but also insights into deformation and relaxation mechanisms across scales. We present here a sampling of key results that highlight both challenges and opportunities for integrated multiscale experiment and modelling in biological systems.