First-Principles Modeling of Protein/Surface Interactions. Polyglycine Secondary Structure Adsorption on the TiO2 (101) Anatase Surface Adopting a Full Periodic Approach

Computational modeling of protein/surface systems is challenging since the conformational variations of the protein and its interactions with the surface need to be considered at once. Adoption of first-principles methods to this purpose is overwhelming and computationally extremely expensive so that, in many cases, dramatically simplified systems (e.g., small peptides or amino acids) are used at the expenses of modeling nonrealistic systems. In this work, we propose a cost-effective strategy for the modeling of peptide/surface interactions at a full quantum mechanical level, taking the adsorption of polyglycine on the TiO₂ (101) anatase surface as a test case. Our approach is based on applying the periodic boundary conditions for both the surface model and the polyglycine peptide, giving rise to full periodic polyglycine/TiO₂ surface systems. By proceeding this way, the considered complexes are modeled with a drastically reduced number of atoms compared with the finite-analogous systems, modeling the polypeptide structures at the same time in a realistic way. Within our modeling approach, full periodic density functional theory calculations (including implicit solvation effects) and ab initio molecular dynamics (AIMD) simulations at the PBE-D2* theory level have been carried out to investigate the adsorption and relative stability of the different polyglycine structures (i.e., extended primary, β-sheet, and α-helix) on the TiO₂ surface. It has been found that, upon adsorption, secondary structures become partially denatured because the peptide C═O groups form Ti-O═C dative bonds. AIMD simulations have been fundamental to identify these phenomena because thermal and entropic effects are of paramount importance. Irrespective of the simulated environments (gas phase and implicit solvent), adsorption of the α-helix is more favorable than that of the β-sheet because in the former, more Ti-O═C bonds are formed and the adsorbed secondary structure results less distorted with respect to the isolated state. Under the implicit water solvent, additionally, adsorbed β-sheet structures weaken with respect to their isolated states as the H-bonds between the strands are longer due to solvation effects. Accordingly, the results indicate that the preferred conformation upon adsorption is the α-helix over the β-sheet.

Environment-controlled water adsorption at hydroxyapatite/collagen interfaces

Water contributes to the structure of bone by coupling hydroxyapatite to collagen over the hierarchical levels of tissue organization. Bone water exists in two states, bound or mobile, each accomplishing different roles. Although many experimental studies show that the amount of bound water correlates with bone strength, a molecular understanding of the interactions between hydroxyapatite, collagen and water is missing. In this work, we unveil the water adsorption properties of bone tissues at the nanoscale using advanced density functional theory methods. We demonstrate that environmental factors such as collagen conformation or degree of confinement, rather than the surface itself, dictate the adsorption mode, strength and density of water on hydroxyapatite. While the results derived in this paper come from a simplified model of bone tissues, they are consistent with experimental observations and constitute a reasonable starting point for more realistic models of bone tissues. For example, we show that environmental changes expected in aging bone lead to reduced water adsorption capabilities, which is consistent with weaker bones at the macroscale. Our findings provide a new interpretation of molecular interactions in bone tissues with the potential to impact bone repair strategies.

Amino acid adsorption on anatase (101) surface at vacuum and aqueous solution: a density functional study

The adsorption of 20 amino acids (AAs) on the (101) surface of anatase titanium dioxide (TiO₂) has been investigated under the scheme of density functional theory. Through the analysis of adsorption geometries, amino group and side chains of AAs have been identified as the major side to adsorb on TiO₂, while the carboxyl group prefers to stay outside to avoid the repulsion between negatively charged oxygen from TiO₂ and AAs. On the surface, two-coordinated oxygen is the major site to stabilize AAs through O-H interactions. The above conclusion does not change when it is in the aqueous solution based on the calculations with AAs surrounded by explicit water molecules. The above knowledge is helpful in predicting how AAs and even peptides adsorb on inorganic materials. Graphical abstract The adsorption of 20 amino acids (AAs) on the (101) surface of anatase titanium dioxide (TiO₂) has been investigated under the scheme of density functional theory.

The effect of water on the binding of glycosaminoglycan saccharides to hydroxyapatite surfaces: a molecular dynamics study

Classical molecular dynamics (MD) simulations have been employed to study the interaction of the saccharides glucuronic acid (GlcA) and N-acetylgalactosamine (GalNAc) with the (0001) and (011̄0) surfaces of the mineral hydroxyapatite (HAP). GlcA and GalNAc are the two constituent monosaccharides of the glycosaminoglycan chondroitin sulfate, which is commonly found in bone and cartilage and has been implicated in the modulation of the hydroxyapatite biomineralization process. MD simulations of the mineral surfaces and the saccharides in the presence of solvent water allowed the calculation of the adsorption energies of the saccharides on the HAP surfaces. The calculations show that GalNAc interacts with HAP principally through the sulfate and the carbonyl of acetyl amine groups, whereas the GlcA interacts primarily through the carboxylate functional groups. The mode and strength of the interaction depends on the orientation of the saccharide with respect to the surface and the level of disruption of the layer of water competing with the saccharide for adsorption sites on the HAP surface, suggesting that chondroitin 4-sulfate binds to the layer of solvent water rather than to HAP.

A fundamental understanding of catechol and water adsorption on a hydrophilic silica surface: exploring the underwater adhesion mechanism of mussels on an atomic scale

Mussels have a remarkable ability to bond to solid surfaces under water. From a microscopic perspective, the first step of this process is the adsorption of dopa molecules to the solid surface. In fact, it is the catechol part of the dopa molecule that is interacting with the surface. These molecules are able to make reversible bonds to a wide range of materials, even underwater. Previous experimental and theoretical efforts have produced only a limited understanding of the mechanism and quantitative details of the competitive adsorption of catechol and water on hydrophilic silica surfaces. In this work, we uncover the nature of this competitive absorption by atomic scale modeling of water and catechol adsorbed at the geminal (001) silica surface using density functional theory calculations. We find that catechol molecules displace preadsorbed water molecules and bond directly on the silica surface. Using molecular dynamics simulations, we observe this process in detail. We also calculate the interaction force as a function of distance, and observe a maximum of 0.5 nN of attraction. The catechol has a binding energy of 23 kcal/mol onto the silica surface with adsorbed water molecules.

Co-adsorption of water and glycine on Cu{110}

In this study, we use density functional theory (DFT) to investigate the surface co-adsorption of glycine with water on Cu{110}. Our results show that, under UHV conditions and for a wide range of temperatures, a pure glycine monolayer is more stable than either mixed gly-water phases or pure water (ice) monolayers, but for a high water pressure half-dissociated water layers can appear on the surface at low and medium temperatures.

Does adsorption at hydroxyapatite surfaces induce peptide folding? Insights from large-scale B3LYP calculations

Large-scale periodic quantum mechanical calculations (509 atoms, 7852 atomic orbitals) based on the hybrid B3LYP functional focused on the peptide folding induced by the adsorption on the (001) and (010) hydroxyapatite (HA) surfaces give interesting insights on the role of specific interactions between surface sites and the peptide, which stabilize the helix conformation over the "native" random coil ones for in silico designed model peptides. The two peptides were derived from the 12-Gly oligomer, with one (P1, C-tGGKGGGGGGEGGN-t) and two (P2, C-tGGKGGKEGGEGGN-t) glutamic acid (E) and lysine (K) residue mutations. The most stable gas-phase "native" conformation for both peptides resulted in a random coil (RC) structure, with the helix (H) conformation being ≈100 kJ mol(-1) higher in free energy. The two peptide conformations interact with the HA (001) and (010) surfaces by C═O groups via Ca(2+) ions, by hydrogen bond between NH(2) groups and the basic PO(4)(3-) groups and by a relevant fraction due to dispersion forces. Peptide adsorption was studied on the dry (001) surface, the wet one envisaging 2 H(2)O per surface Ca(2+) and, on the latter, also considering the adsorption of microsolvated peptides with 4 H(2)O molecules located at sites responsible of the interaction with the surface. The P1 mutant does prefer to be adsorbed as a random coil by ≈160 kJ/mol, whereas the reverse is computed for P2, preferring the helix conformation by ≈50 kJ/mol. Adsorption as helix of both P1 and P2 mutants brings about proton transfer toward the HA surfaces with a large charge transfer component to the interaction energy.

Ab initio modelling of protein-biomaterial interactions: influence of amino acid polar side chains on adsorption at hydroxyapatite surfaces

The adsorption from the gas phase of five different amino acids (AAs), namely Gly, Ser, Lys, Gln and Glu, on three surface models of hexagonal hydroxyapatite (HA) has been studied at B3LYP level with Gaussian type basis set within a periodic approach. The AA adsorption was simulated on the (001) and (010) stoichiometric surfaces, the latter both in its pristine and water-reacted form. Low/high AA coverage has been studied by doubling the HA unit cell size. The AAs have been docked to the HA surfaces following the electrostatic complementarity between the electrostatic potentials of AA and the bare HA. Gly adsorbs as a zwitterion at the (001) surface, whereas at the (010) ones, the proton of the COOH group is transferred to the surface resulting in an HA(+)/Gly(-) ion pair. For the other AAs, the common COOH-CH-NH(2) moiety behaves like in Gly, while the specific side-chain functionalities adsorb at the HA surfaces by maximizing electrostatic and H-bond interactions. The interactions between the side chains and the HA surface impart a higher stability compared with the Gly case, with Glu being the strongest adsorbate owing to its high Ca affinity and H-bond donor propensity. For AAs of large size, the adsorption is more favourable in conditions of low coverage as repulsion between adjacent AAs is avoided. For all considered AAs, the strongest interaction is always established on the (010) faces rather than on the (001) one, therefore suggesting an easier growth along the c-direction of HA crystals from AA solutions.

Biomolecule-biomaterial interaction: a DFT-D study of glycine adsorption and self-assembly on hydroxylated Cr2O3 surfaces

The adsorption of glycine, the building block of amino acids, on hydroxylated (0001)-Cr2O3 model surfaces, representing the stainless steel passive film surface, was modeled by means of the GGA + U method. The roles of glycine coverage and surface termination (hydroxylated Cr- and O-terminated surfaces) on the adsorption mode and self-assembly properties were explored. The hydroxylated Cr-terminated Cr2O3 surface, which presents two types of (H)OH groups exhibiting different acidic character, is more reactive than the hydroxylated O-terminated surface, where one single type of OH group is present, for all adsorption modes and coverages considered. Outer sphere adsorption occurs in the zwitterion form, stabilized at low coverage through H-bond formation with coadsorbed water molecules, and at the monolayer coverage by glycine self-assembling. The OH substitution by glycinate is favored on the hydroxylated Cr-terminated surface and not on the O-terminated one. The inclusion of dispersion forces does not change the observed tendencies. An atomistic thermodynamics approach suggests that outer sphere adsorption is thermodynamically favored over inner sphere adsorption in the whole domain of glycine concentration. The obtained SAM's free energies of formation are rationalized in a model considering the balance between sublimation and solvation free energies, and extrapolated to other amino acids, to predict the SAMs formation above hydroxylated surfaces. It is found that hydrophobic AA tend to self-assemble at the surface, whereas hydrophilic ones do not.

Affinity of hydroxyapatite (001) and (010) surfaces to formic and alendronic acids: a quantum-mechanical and infrared study

The affinity of the (001) and of the water reacted (010)WR hydroxyapatite surfaces towards formic and alendronic acids is studied with density functional theory (PBE functional) using periodic boundary conditions based on Gaussian basis set. Structures, energetic of the adsorption and vibrational features of the adsorbates are computed in order to understand at the atomic level both the cariogenic processes (for the formic acid) and the features of anti-osteoporosis drugs (for the alendronic acid). For both molecules the interaction energy is very high on an absolute scale, and for all examined cases, it is higher on the (010)WR HA surface than on the (001) one. For the latter, a number of cases by which the acidic proton of the adsorbate is transferred to the HA surface are also characterized. For the formic acid case, experimental infrared spectra are also measured and the position and nature of the C=O stretching bands have been found to be in excellent agreement with the quantum mechanical simulations. For alendronic acid IR experiments are still not available and the present predicted infrared spectra will be useful as a guide to interpret future experimental studies.

Structural evolution of protein-biofilms: Simulations and experiments

The control of biofilm formation is a challenging goal that has not been reached yet in many aspects. One unsolved question is the role of van der Waals forces and another is the importance of mutual interactions between the adsorbing and the adsorbed biomolecules ("critical crowding"). In this study, a combined experimental and theoretical approach is presented, which fundamentally probes both aspects. On three model proteins-lysozyme, α-amylase, and bovine serum albumin-the adsorption kinetics is studied experimentally. Composite substrates are used enabling a separation of the short- and the long-range forces. Although usually neglected, experimental evidence is given for the influence of van der Waals forces on the protein adsorption as revealed by in situ ellipsometry. The three proteins were chosen for their different conformational stabilities in order to investigate the influence of conformational changes on the adsorption kinetics. Monte Carlo simulations are used to develop a model for these experimental results by assuming an internal degree of freedom to represent conformational changes. The simulations also provide data on the distribution of adsorption sites. By in situ atomic force microscopy we can also test this distribution experimentally, which opens the possibility to, e.g., investigate the interactions between adsorbed proteins.

A density functional theory study of the interaction of collagen peptides with hydroxyapatite surfaces

Density functional theory calculations were applied to investigate the binding of four peptide strands, which are important in the collagen protein, to the bone and tooth mineral hydroxyapatite: amphiphilic PRO-HYP-GLY and HYP-PRO-GLY, and hydrophobic PRO-LYS-GLY and PRO-HYL-GLY. The particular peptide sequences are chosen for their different functional groups, containing (i) hydrophobic; (ii) uncharged polar; and (iii) charged polar side groups, thus allowing direct comparison of the general effect of these carboxylic acid and amine functional groups, as well as hydroxylation and charge, on their interactions with two major hydroxyapatite surfaces, (0001) and (0110). The calculated results are consistent with experiments, confirming that the terminal carboxyl groups and amine groups mainly contribute to the adsorption of the peptides to the hydroxyapatite surfaces and primarily to the (0110) surface rather than the dominant (0001) plane. Of the side groups in the tripeptide motifs representing the collagen protein, the -OH and positively charged -NH(3)(+) groups in particular bind strongly to the surfaces, and their presence should therefore promote hydroxyapatite growth.