Ordering surfaces on the nanoscale: implications for protein adsorption

Monolayer-protected metal nanoparticles (MPMNs) are a newly discovered class of nanoparticles with an ordered, striped domain structure that can be readily manipulated by altering the ratio of the hydrophobic to hydrophilic ligands. This property makes them uniquely suited to systematic studies of the role of nanostructuring on biomolecule adsorption, a phenomenon of paramount importance in biomaterials design. In this work, we examine the interaction of the simple, globular protein cytochrome C (Cyt C) with MPMN surfaces using experimental protein assays and computational molecular dynamics simulations. Experimental assays revealed that adsorption of Cyt C generally increased with increasing surface polar ligand content, indicative of the dominance of hydrophilic interactions in Cyt C-MPMN binding. Protein-surface adsorption enthalpies calculated from computational simulations employing rigid-backbone coarse-grained Cyt C and MPMN models indicate a monotonic increase in adsorption enthalpy with respect to MPMN surface polarity. These results are in qualitative agreement with experimental results and suggest that Cyt C does not undergo significant structural disruption upon adsorption to MPMN surfaces. Coarse-grained and atomistic simulations furthermore elucidated the important role of lysine in facilitating Cyt C adsorption to MPMN surfaces. The amphipathic character of the lysine side chain enables it to form close contacts with both polar and nonpolar surface ligands simultaneously, rendering it especially important for interactions with surfaces composed of adjacent nanoscale chemical domains. The importance of these structural characteristics of lysine suggests that proteins may be engineered to specifically interact with nanomaterials by targeted incorporation of unnatural amino acids possessing dual affinity to differing chemical motifs.     

Adsorption of homopolypeptides on gold investigated using atomistic molecular dynamics

We investigate the role of dynamics on adsorption of peptides to gold surfaces using all-atom molecular dynamics simulations in explicit solvent. We choose six homopolypeptides [Ala(10), Ser(10), Thr(10), Arg(10), Lys(10), and Gln(10)], for which experimental surface coverages are not correlated with amino acid level affinities for gold, with the idea that dynamic properties may also play a role. To assess dynamics we determine both conformational movement and flexibility of the peptide within a given conformation. Low conformational movement indicates stability of a given conformation and leads to less adsorption than homopolypeptides with faster conformational movement. Likewise, low flexibility within a given conformation also leads to less adsorption. Neither amino acid affinities nor dynamic considerations alone predict surface coverage; rather both quantities must be considered in peptide adsorption to gold surfaces.     

Formation of stacking defects at surfaces: From atomistic simulations to density functional theory calculations

We have performed electronic structure calculations to study the evolution of the stacking fault energy at (111) surfaces of metals. We first apply an sp–d tight-binding model and then increase the accuracy on the electronic structure by using density functional theory (DFT) calculations. We show in this way the relative importance of sp–d hybridization both in the formation of defects at the surface of metals and in reconstruction phenomena as a function of band filling especially at the end of transition metal series. Comparing our results with atomistic simulations it is concluded that although atomistic calculations are powerful tools to investigate relaxation mechanisms at surfaces, a higher degree of accuracy on electronic structure is necessary to quantify the energy of some defects at surfaces like stacking faults. In particular long range interactions associated to less localized sp electrons are playing a rather important role in reconstruction phenomena for metals like platinum and gold. These results are backed up by DFT calculations applied to iridium, platinum and gold (111) surfaces.     

Simulation studies of the protein-water interface. I. Properties at the molecular resolution

We report molecular dynamics simulations of three globular proteins: ubiquitin, apo-calbindin D(9K), and the C-terminal SH2 domain of phospholipase C-gamma1 in explicit water. The proteins differ in their overall charge and fold type and were chosen to represent to some degree the structural variability found in medium-sized proteins. The length of each simulation was at least 15 ns, and larger than usual solvent boxes were used. We computed radial distribution functions, as well as orientational correlation functions about the surface residues. Two solvent shells could be clearly discerned about charged and polar amino acids. Near apolar amino acids the water density near such residues was almost devoid of structure. The mean residence time of water molecules was determined for water shells about the full protein, as well as for water layers about individual amino acids. In the dynamic properties, two solvent shells could be characterized as well. However, by comparison to simulations of pure water it could be shown that the influence of the protein reaches beyond 6 A, i.e., beyond the first two shells. In the first shell (r < or =3.5 A), the structural and dynamical properties of solvent waters varied considerably and depended primarily on the physicochemical properties of the closest amino acid side chain, with which the waters interact. By contrast, the solvent properties seem not to depend on the specifics of the protein studied (such as the net charge) or on the secondary structure element in which an amino acid is located. While differing considerably from the neat liquid, the properties of waters in the second solvation shell (3.5< r < or =6 A) are rather uniform; a direct influence from surface amino acids are already mostly shielded.     

DYNAMICS OF SINGLE-MOLECULE ROTATIONS ON SURFACES DEPEND ON SYMMETRY, INTERACTIONS AND MOLECULAR SIZES

Rotating surface-mounted molecules have attracted attention of many research groups as a way to develop new nanoscale devices and materials. However, mechanisms of motion of these rotors at the single-molecule level are still not well understood. Theoretical and experimental studies on thioether molecular rotors on gold surfaces suggest that the size of the molecules, their flexibility and steric repulsions with the surface are important for dynamics of the system. A complex combination of these factors leads to the observation that the rotation speeds have not been hindered by increasing the length of the alkyl chains. However, experiments on diferrocene derivatives indicated that a significant increase in the rotational barriers for longer molecules. We present here a comprehensive theoretical study that combines molecular dynamics simulations and simple models to investigate what factors influence single-molecule rotations on the surfaces. Our results suggest that rotational dynamics is determined by the size and by the symmetry of the molecules and surfaces, and by interactions with surfaces. Our theoretical predictions are in excellent agreement with current experimental observations.     

Reversible oriented immobilization of histidine-tagged proteins on gold surfaces using a chelator thioalkane

The reversible oriented immobilization of proteins on solid surfaces is a prerequisite for the investigation of molecular interactions at interfaces or the construction of supramolecular assemblies. We demonstrate a generally applicable method using a synthetic chelator thioalkane which can self-assemble on a gold surface via its thiol group. It exposes its nitrilotriacetic acid group which serves as a chelator for transition metal ions. Reversible binding of a F_ab fragment modified with a C-terminal hexahistidine extension was monitored in situ using surface plasmon resonance. The directed immobilization of proteins on surfaces opens new ways for structural investigations of proteins and the development of biosensors.     

Atomistic Simulations of Dopamine Diffusion Dynamics on a Pristine Graphene Surface**

Carbon microelectrodes enable in vivo detection of neurotransmitters, and new electrodes aim to optimize the carbon surface. However, atomistic detail on the diffusion and orientation of neurotransmitters near these surfaces is lacking. Here, we employ molecular dynamics simulations to investigate the surface diffusion of dopamine (DA), its oxidation product dopamine-o-quinone (DOQ), and their protonated forms on the pristine basal plane of flat graphene. We find that all DA species rapidly adsorb to the surface and remain adsorbed, even without a holding potential or graphene surface defects. We also find that the diffusivities of the adsorbed and the fully solvated DA are similar and that the protonated species diffuse more slowly on the surface than their corresponding neutral forms, while the oxidized species diffuse more rapidly. Structurally, we find that the underlying graphene lattice has little influence over the molecular adsorbate's lateral position, and the vertical placement of the amine group on dopamine is highly dependent upon its charge. Finally, we find that solvation has a large effect on surface diffusivities. These first results from molecular dynamics simulations of dopamine at the aqueous-graphene interface show that dopamine diffuses rapidly on the surface, even without an applied potential, and provide a basis for future simulations of neurotransmitter structure and dynamics on advanced carbon materials electrodes.     

Preparation of metal nanoband microelectrode on poly(dimethylsiloxane) for chip-based amperometric detection

We proposed herein a novel approach for fabricating nanoband microelectrodes for electrochemical detection on an electrophoresis microchip. The metal films were first obtained via region-selective electroless deposition of gold or copper films on PDMS substrates by selective region plasma oxidation through shadow masking. Both metal films show uniform surfaces with the thickness at the level of 100 nm. By casting another PDMS layer on the metal films, the cross section of the sandwich structures can be used as nanoband microelectrodes, which can be renewed just by cutting. These nanoband microelectrodes are successfully used as electrochemical detectors in microchip electrophoresis for the detection of amino acids, proteins and neurotransmitter molecules. Moreover, integrating an Au-Cu double-metal detector with a double-channel electrophoresis system, we can easily distinguish electroactive amino acids from that of non-electroactive amino acids.     

Free energy surfaces for the interaction of D-glucose with planar aromatic groups in aqueous solution

Multidimensional potentials of mean force for the interactions in aqueous solution of both anomers of D-glucopyranose with two planar aromatic molecules, indole and para-methyl-phenol, have been calculated using molecular dynamics simulations with umbrella sampling and were subsequently used to estimate binding free energies. Indole and para-methyl-phenol serve as models for the side chains of the amino acids tryptophan and tyrosine, respectively. In all cases, a weak affinity between the glucose molecules and the flat aromatic surfaces was found. The global minimum for these interactions was found to be for the case when the pseudoplanar face of β-D-glucopyranose is stacked against the planar surfaces of the aromatic residues. The calculated binding free energies are in good agreement with both experiment and previous simulations. The multidimensional free energy maps suggest a mechanism that could lend kinetic stability to the complexes formed by sugars bound to sugar-binding proteins.     

Modeling van der Waals interactions between proteins and inorganic surfaces from time-dependent density functional theory calculations

In this work we show how the ab initio determination of van der Waals coefficients within time-dependent density functional theory can be used to build efficient and accurate atomistic models that describe the long-range interactions of proteins with other proteins and of proteins with semi-conducting surfaces. The model parameters are fitted so that they reproduce the ab initio van der Waals coefficients of amino acids and dipeptides. We then assess the quality of our results by comparing ab initio van der Waals coefficients for larger peptides with the coefficients yielded by the models. The different sets of parameters can be easily incorporated in current empirical force field methods, thus providing an essential ingredient for molecular dynamics simulations of proteins close to surfaces.     

Self-assembly of organic molecules at metal surfaces

Self-assembly represents a promising strategy for surface functionalisation as well as creating nanostructures with well-controlled, tailor-made properties and functionality. Molecular self-assembly at solid surfaces is governed by the subtle interplay between molecule–molecule and molecule–substrate interactions that can be tuned by varying molecular building blocks, surface chemistry and structure as well as substrate temperature.In this review, basic principles behind molecular self-assembly of organic molecules on metal surfaces will be discussed. Controlling these formation principles allows for creating a wide variety of different molecular surface structures ranging from well-defined clusters, quasi one-dimensional rows to ordered, two-dimensional overlayers. An impressive number of studies exist, demonstrating the ability of molecular self-assembly to create these different structural motifs in a predictable manner by tuning the molecular building blocks as well as the metallic substrate.Here, the multitude of different surface structures of the natural amino acid cysteine on two different gold surfaces observed with scanning tunnelling microscopy will be reviewed. Cysteine on Au(110)-(1×2) represents a model system illustrating the formation of all the above mentioned structural motifs without changing the molecular building blocks or the substrate surface. The only parameters in this system are substrate temperature and molecular coverage, controlling both the molecular adsorption state (physisorption versus chemisorption) and molecular surface mobility. By tuning the adsorption state and the molecular mobility, distinctly different molecular structures are formed, exemplifying the variety of structural motifs that can be achieved by molecular self-assembly.     

Multiplexed Discrimination of Single Amino Acid Residues in Polypeptides in a Single SERS Hot Spot

The SERS-based detection of protein sequences with single-residue sensitivity suffers from signal dominance of aromatic amino acid residues and backbones, impeding detection of non-aromatic amino acid residues. Herein, we trap a gold nanoparticle in a plasmonic nanohole to generate a single SERS hot spot for single-molecule detection of 2 similar polypeptides (vasopressin and oxytocin) and 10 distinct amino acids that constitute the 2 polypeptides. Significantly, both aromatic and non-aromatic amino acids are detected and discriminated at the single-molecule level either at individual amino acid molecules or within the polypeptide chains. Correlated with molecular dynamics simulations, our results suggest that the signal dominance due to large spatial occupancy of aromatic rings of the polypeptide sidechains on gold surfaces can be overcome by the high localization of the single hot spot. The superior spectral and spatial discriminative power of our approach can be applied to single-protein analysis, fingerprinting, and sequencing.     

Molecular Dynamics Simulations of a β-Hairpin Fragment of Protein G by Means of Atom-Bond Electronegativity Equalization Method Fused into Molecular Mechanics (ABEEMσπ/MM)

There are some controversial opinions about the origin of folding β‐hairpin stability in aqueous solution. In this study, the structural and dynamic behavior of a 16‐residue β‐hairpin from B1 domain of protein G has been investigated at 280, 300, 350 and 450 K using molecular dynamics (MD) simulations by means of Atom‐Bond Electronegativity Equalization Method Fused into Molecular Mechanics i.e., ABEEMδπ/MM and the explicit ABEEM‐7P water solvent model. In addition, a 300 K simulation of one mutant having the aromatic residues substituted with alanines has been performed. The hydrophobic surface area, hydrophilic surface area and some structural properties have been used to measure the role of the hydrophobic interactions. It is found that the aromatic residues substituted with alanines have shown an evident destabilization of the structure and unfolding started after 1.5 ns. It is also found that the number of the main chain hydrogen bonds have different distributions through three different simulations. All above demonstrate that the hydrophobic interactions and the main chain hydrogen bonds play an important role in the stability of the folding structure of β‐hairpin in solution. Furthermore, through the structural analyses of the β‐hairpin structures from four temperature simulations and the comparison with other MD simulations of β‐hairpin peptides, the new ABEEMδπ force field can reproduce the structural data in good agreement with the experimental data.     

Structural and Dynamic Properties of Monoclonal Antibodies Immobilized on CNTs: A Computational Study

Due to the widespread application of carbon nanotube (CNT)‐based materials in nanomedicine, it is nowadays of paramount importance to unravel at the atomistic level of detail the structural properties of such bioconjugates in order to rationalize and predict the effect exerted by the graphitic framework on the bio‐active counterpart. In this paper, we report for the first time all‐atom explicit solvent molecular dynamics (MD) simulations investigating the structural and dynamic properties of a noncovalent bioconjugate in which the monoclonal Cetuximab antibody (Ctx) is adsorbed on a CNT surface. Upon selection of the three most representative adsorption modes as obtained by docking studies, force‐field MD and DFT simulations unambiguously showed that hydrophobic interactions mainly govern the adsorption of the protein on the graphitic surface. Two main adsorption poses have been predicted: a pose‐fab (p‐fab) and pose‐fc (p‐fc) (fab = fragment antigen binding region; fc = fragment crystallizable region), the former being favored with small‐diameter tubes (≤40 Å). In all the predicted poses, the secondary structure of Ctx is largely unaffected by the presence of the graphitic surface and, consistently with previous literature studies, our simulations reveal that positively charged amino acidic residues, such as Lys and Arg, predominantly contribute to the stabilization of the CNT?Ctx complex acting like surfactants. The predicted structural models are consistent with the experimental data, for which the immobilization of the antibody on CNTs does not disrupt the structural and recognition properties of the Ctx, consequently supporting the reliability of the used bioconjugation strategy for engineering stable and responsive hybrid nanomaterials for therapeutic applications. Moreover, a remarkable structural similarity of Ctx with antibodies of different isotypes suggests that in principle the CNT framework can interact in the same manner with all antibodies currently used in clinical applications.     

Synthesis,experimental and in silico studies of N‐fluorenylmethoxycarbonyl‐O‐tert‐butyl‐N‐methyltyrosine,coupled with CSD data: a survey of interactions in the crystal structures of Fmoc–amino acids

Recently, fluorenylmethoxycarbonyl (Fmoc) amino acids (e.g. Fmoc–tyrosine or Fmoc–phenylalanine) have attracted growing interest in biomedical research and industry, with special emphasis directed towards the design and development of novel effective hydrogelators, biomaterials or therapeutics. With this in mind, a systematic knowledge of the structural and supramolecular features in recognition of those properties is essential. This work is the first comprehensive summary of noncovalent interactions combined with a library of supramolecular synthon patterns in all crystal structures of amino acids with the Fmoc moiety reported so far. Moreover, a new Fmoc‐protected amino acid, namely, 2‐{[(9H‐fluoren‐9‐ylmethoxy)carbonyl](methyl)amino}‐3‐{4‐[(2‐hydroxypropan‐2‐yl)oxy]phenyl}propanoic acid or N‐fluorenylmethoxycarbonyl‐O‐tert‐butyl‐N‐methyltyrosine, Fmoc‐N‐Me‐Tyr(t‐Bu)‐OH, C₂₉H₃₁NO₅, was successfully synthesized and the structure of its unsolvated form was determined by single‐crystal X‐ray diffraction. The structural, conformational and energy landscape was investigated in detail by combined experimental and in silico approaches, and further compared to N‐Fmoc‐phenylalanine [Draper et al. (2015). CrystEngComm, 42 , 8047–8057]. Geometries were optimized by the density functional theory (DFT) method either in vacuo or in solutio. The polarizable conductor calculation model was exploited for the evaluation of the hydration effect. Hirshfeld surface analysis revealed that H…H, C…H/H…C and O…H/H…O interactions constitute the major contributions to the total Hirshfeld surface area in all the investigated systems. The molecular electrostatic potentials mapped over the surfaces identified the electrostatic complementarities in the crystal packing. The prediction of weak hydrogen‐bonded patterns via Full Interaction Maps was computed. Supramolecular motifs formed via C—H…O, C—H…π, (fluorenyl)C—H…Cl(I), C—Br…π(fluorenyl) and C—I…π(fluorenyl) interactions are observed. Basic synthons, in combination with the Long‐Range Synthon Aufbau Modules, further supported by energy‐framework calculations, are discussed. Furthermore, the relevance of Fmoc‐based supramolecular hydrogen‐bonding patterns in biocomplexes are emphasized, for the first time.     

Multiplexed Discrimination of Single Amino Acid Residues in Polypeptides in a Single SERS Hot Spot

The SERS‐based detection of protein sequences with single‐residue sensitivity suffers from signal dominance of aromatic amino acid residues and backbones, impeding detection of non‐aromatic amino acid residues. Herein, we trap a gold nanoparticle in a plasmonic nanohole to generate a single SERS hot spot for single‐molecule detection of 2 similar polypeptides (vasopressin and oxytocin) and 10 distinct amino acids that constitute the 2 polypeptides. Significantly, both aromatic and non‐aromatic amino acids are detected and discriminated at the single‐molecule level either at individual amino acid molecules or within the polypeptide chains. Correlated with molecular dynamics simulations, our results suggest that the signal dominance due to large spatial occupancy of aromatic rings of the polypeptide sidechains on gold surfaces can be overcome by the high localization of the single hot spot. The superior spectral and spatial discriminative power of our approach can be applied to single‐protein analysis, fingerprinting, and sequencing.     

Nanostructure‐Directed Physisorption vs Chemisorption at Semiconductor Interfaces: The Inverse of the HSAB Concept

A concept, complementary to that of hard and soft acid–base interactions (HSAB‐dominant chemisorption) and consistent with dominant physisorption to a semiconductor interface, is presented. We create a matrix of sensitivities and interactions with several basic gases. The concept, based on the reversible interaction of hard‐acid surfaces with soft bases, hard‐base surfaces with soft acids, or vice versa, corresponds 1) to the inverse of the HSAB concept and 2) to the selection of a combination of semiconductor interface and analyte materials, which can be used to direct a physisorbed vs chemisorbed interaction. The technology, implemented on nanopore coated porous silicon micropores, results in the coupling of acid–base chemistry with the depletion or enhancement of majority carriers in an extrinsic semiconductor. Using the inverse‐HSAB (IHSAB) concept, significant and predictable changes in interface sensitivity for a variety of gases can be implemented. Nanostructured metal oxide particle depositions provide selectivity and complement a highly efficient electrical contact to a porous silicon nanopore covered microporous interface. The application of small quantities (much less than a monolayer) of nanostructured metals, metal oxides, and catalysts which focus the physisorbtive and chemisorbtive interactions of the interface, can be made to create a range of notably higher sensitivities for reversible physisorption. This is exemplified by an approach to reversible, sensitive, and selective interface responses. Nanostructured metal oxides developed from electroless gold (Au_xO), tin (SnO₂), copper (Cu_xO), and nickel (NiO) depositions, nanoalumina, and nanotitania are used to demonstrate the IHSAB concept and provide for the detection of gases, including NH₃, PH₃, CO, NO, and H₂S, in an array‐based format to the sub‐ppm level.     

Role of hydration force in the self-assembly of collagens and amyloid steric zipper filaments

In protein self-assembly, types of surfaces determine the force between them. Yet the extent to which the surrounding water contributes to this force remains as a fundamental question. Here we study three self-assembling filament systems that respectively have hydrated (collagen), dry nonpolar, and dry polar (amyloid) interfaces. Using molecular dynamics simulations, we calculate and compare local hydration maps and hydration forces. We find that the primary hydration shells are formed all over the surface, regardless of the types of the underlying amino acids. The weakly oscillating hydration force arises from coalescence and depletion of hydration shells as two filaments approach, whereas local water diffusion, orientation, or hydrogen-bonding events have no direct effect. Hydration forces between hydrated, polar, and nonpolar interfaces differ in the amplitude and phase of the oscillation relative to the equilibrium surface separation. Therefore, water-mediated interactions between these protein surfaces, ranging in character from "hydrophobic" to "hydrophilic", have a common molecular origin based on the robustly formed hydration shells, which is likely applicable to a broad range of biomolecular assemblies whose interfacial geometry is similar in length scale to those of the present study.