ReaxFF MD Simulations of Peptide-Grafted Gold Nanoparticles

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.

ReaxFF Molecular Dynamics Simulations of Large Gold Nanocrystals

A systematic study of gold nanocrystals is carried out using molecular dynamics simulations with reactive force fields. The nanocrystal size is varied between 2 and 10 nm with methane and butane thiolate as ligands. The reactive force fields allow investigation of the formation of staples. The simulations explain several experimental observations such as the number of staples per thiolate of about 40% and the occupation of the top adsorption sites on the facets. They also show that the frequency of staples is increased on the edges, which leads to a desorption of gold atoms from the nanocrystal edges. In contrast to previous nonreactive simulations, no difference between the distances of the ligands on the nanocrystal edges and facets is observed. Except for the 2 nm particles, the nanocrystal size and the alkane chain length of the ligands have only a small influence on the nanocrystal properties. The occupation of adsorption sites and staple frequencies are very slowly converging properties, taking more than ns.

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.

Peptide-Conjugated Nano Delivery Systems for Therapy and Diagnosis of Cancer

Peptides are strings of approximately 2-50 amino acids, which have gained huge attention for theranostic applications in cancer research due to their various advantages including better biosafety, customizability, convenient process of synthesis, targeting ability via recognizing biological receptors on cancer cells, and better ability to penetrate cell membranes. The conjugation of peptides to the various nano delivery systems (NDS) has been found to provide an added benefit toward targeted delivery for cancer therapy. Moreover, the simultaneous delivery of peptide-conjugated NDS and nano probes has shown potential for the diagnosis of the malignant progression of cancer. In this review, various barriers hindering the targeting capacity of NDS are addressed, and various approaches for conjugating peptides and NDS have been discussed. Moreover, major peptide-based functionalized NDS targeting cancer-specific receptors have been considered, including the conjugation of peptides with extracellular vesicles, which are biological nanovesicles with promising ability for therapy and the diagnosis of cancer.

Highly Heterogeneous Polarization and Solvation of Gold Nanoparticles in Aqueous Electrolytes

The performance of gold nanoparticles (NPs) in applications depends critically on the structure of the NP-solvent interface, at which the electrostatic surface polarization is one of the key characteristics that affects hydration, ionic adsorption, and electrochemical reactions. Here, we demonstrate significant effects of explicit metal polarizability on the solvation and electrostatic properties of bare gold NPs in aqueous electrolyte solutions of sodium salts of various anions (Cl^-, BF₄^-, PF₆^-, nitrophenolate, and 3- and 4-valent hexacyanoferrate), using classical molecular dynamics simulations with a polarizable core-shell model for the gold atoms. We find considerable spatial heterogeneity of the polarization and electrostatic potentials on the NP surface, mediated by a highly facet-dependent structuring of the interfacial water molecules. Moreover, ion-specific, facet-dependent ion adsorption leads to considerable alterations of the interfacial polarization. Compared to nonpolarizable NPs, surface polarization modifies water local dipole densities only slightly but has substantial effects on the electrostatic surface potentials and leads to significant lateral redistributions of ions on the NP surface. Besides, interfacial polarization effects cancel out in the far field for monovalent ions but not for polyvalent ions, as anticipated from continuum "image-charge" concepts. Far-field effective Debye-Hückel surface potentials change accordingly in a valence-specific fashion. Hence, the explicit charge response of metal NPs is crucial for the accurate description and interpretation of interfacial electrostatics (e.g., for charge transfer and interfacial polarization in catalysis and electrochemistry).

The Peptide Functionalized Inorganic Nanoparticles for Cancer-Related Bioanalytical and Biomedical Applications

In order to improve their bioapplications, inorganic nanoparticles (NPs) are usually functionalized with specific biomolecules. Peptides with short amino acid sequences have attracted great attention in the NP functionalization since they are easy to be synthesized on a large scale by the automatic synthesizer and can integrate various functionalities including specific biorecognition and therapeutic function into one sequence. Conjugation of peptides with NPs can generate novel theranostic/drug delivery nanosystems with active tumor targeting ability and efficient nanosensing platforms for sensitive detection of various analytes, such as heavy metallic ions and biomarkers. Massive studies demonstrate that applications of the peptide-NP bioconjugates can help to achieve the precise diagnosis and therapy of diseases. In particular, the peptide-NP bioconjugates show tremendous potential for development of effective anti-tumor nanomedicines. This review provides an overview of the effects of properties of peptide functionalized NPs on precise diagnostics and therapy of cancers through summarizing the recent publications on the applications of peptide-NP bioconjugates for biomarkers (antigens and enzymes) and carcinogens (e.g., heavy metallic ions) detection, drug delivery, and imaging-guided therapy. The current challenges and future prospects of the subject are also discussed.

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide-Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide-nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure-function relationship, for technological and nanomedical applications.

Not only in silico drug discovery: Molecular modeling towards in silico drug delivery formulations

The use of methods at molecular scale for the discovery of new potential active ligands, as well as previously unknown binding sites for target proteins, is now an established reality. Literature offers many successful stories of active compounds developed starting from insights obtained in silico and approved by Food and Drug Administration (FDA). One of the most famous examples is raltegravir, a HIV integrase inhibitor, which was developed after the discovery of a previously unknown transient binding area thanks to molecular dynamics simulations. Molecular simulations have the potential to also improve the design and engineering of drug delivery devices, which are still largely based on fundamental conservation equations. Although they can highlight the dominant release mechanism and quantitatively link the release rate to design parameters (size, drug loading, et cetera), their spatial resolution does not allow to fully capture how phenomena at molecular scale influence system behavior. In this scenario, the "computational microscope" offered by simulations at atomic scale can shed light on the impact of molecular interactions on crucial parameters such as release rate and the response of the drug delivery device to external stimuli, providing insights that are difficult or impossible to obtain experimentally. Moreover, the new paradigm brought by nanomedicine further underlined the importance of such computational microscope to study the interactions between nanoparticles and biological components with an unprecedented level of detail. Such knowledge is a fundamental pillar to perform device engineering and to achieve efficient and safe formulations. After a brief theoretical background, this review aims at discussing the potential of molecular simulations for the rational design of drug delivery systems.

ReaxFF/AMBER-A Framework for Hybrid Reactive/Nonreactive Force Field Molecular Dynamics Simulations

Combined quantum mechanical/molecular mechanical (QM/MM) models using semiempirical and ab initio methods have been extensively reported on over the past few decades. These methods have been shown to be capable of providing unique insights into a range of problems, but they are still limited to relatively short time scales, especially QM/MM models using ab initio methods. An intermediate approach between a QM based model and classical mechanics could help fill this time-scale gap and facilitate the study of a range of interesting problems. Reactive force fields represent the intermediate approach explored in this paper. A widely used reactive model is ReaxFF, which has largely been applied to materials science problems and is generally used as a stand-alone (i.e., the full system is modeled using ReaxFF). We report a hybrid ReaxFF/AMBER molecular dynamics (MD) tool, which introduces ReaxFF capabilities to capture bond breaking and formation within the AMBER MD software package. This tool enables us to study local reactive events in large systems at a fraction of the computational costs of QM/MM models. We describe the implementation of ReaxFF/AMBER, validate this implementation using a benzene molecule solvated in water, and compare its performance against a range of similar approaches. To illustrate the predictive capabilities of ReaxFF/AMBER, we carried out a Claisen rearrangement study in aqueous solution. In a first for ReaxFF, we were able to use AMBER's potential of mean force (PMF) capabilities to perform a PMF study on this organic reaction. The ability to capture local reaction events in large systems using combined ReaxFF/AMBER opens up a range of problems that can be tackled using this model to address both chemical and biological processes.

Ion-Specific Adsorption on Bare Gold (Au) Nanoparticles in Aqueous Solutions: Double-Layer Structure and Surface Potentials

We study the solvation and electrostatic properties of bare gold (Au) nanoparticles (NPs) of 1-2 nm in size in aqueous electrolyte solutions of sodium salts of various anions with large physicochemical diversity (Cl^-, BF₄^-, PF₆^-, Nip^- (nitrophenolate), 3- and 4-valent hexacyanoferrate (HCF)) using nonpolarizable, classical molecular dynamics computer simulations. We find a substantial facet selectivity in the adsorption structure and spatial distribution of the ions at the AuNPs: while sodium and some of the anions (e.g., Cl^-, HCF^3-) adsorb more at the "edgy" (100) and (110) facets of the NPs, where the water hydration structure is more disordered, other ions (e.g., BF₄^-, PF₆^-, Nip^-) prefer to adsorb strongly on the extended and rather flat (111) facets. In particular, Nip^-, which features an aromatic ring in its chemical structure, adsorbs strongly and perturbs the first water monolayer structure on the NP (111) facets substantially. Moreover, we calculate adsorptions, radially resolved electrostatic potentials as well as the far-field effective electrostatic surface charges and potentials by mapping the long-range decay of the calculated electrostatic potential distribution onto the standard Debye-Hückel form. We show how the extrapolation of these values to other ionic strengths can be performed by an analytical Adsorption-Grahame relation between the effective surface charge and potential. We find for all salts negative effective surface potentials in the range from -10 mV for NaCl down to about -80 mV for NaNip, consistent with typical experimental ranges for the zeta potential. We discuss how these values depend on the surface definition and compare them to the explicitly calculated electrostatic potentials near the NP surface, which are highly oscillatory in the ±0.5 V range.