Quantitatively Identifying the Roles of Interfacial Water and Solid Surface in Governing Peptide Adsorption

The impact of chemical properties of the solid-liquid-adsorbate interfaces on the entropy-enthalpy compensation involved in adsorption

Despite extensive studies on the thermodynamic mechanism governing molecular adsorption at the solid-water interface, a comprehensive understanding of the crucial role of interface properties in mediating the entropy-enthalpy compensation during adsorption is lacking, particularly at a quantitative level. Herein, we employed two types of surface models (hydroxyapatite and graphene) along with a series of amino acids to successfully elucidate how distinct interfacial features dictate the delicate balance between entropy and enthalpy variations. The adsorption of all amino acids on the hydroxyapatite surface is an enthalpy-dominated process, where the water-induced enthalpic component of the free energy and the surface-adsorbate electrostatic interaction term alternatively act as the driving force for adsorption in different regions of the surface. Although favorable interactions are observed between amino acids and the graphene surface, the entropy-enthalpy compensation exhibits dependence on the molecular size of the adsorbates. For small amino acids, favorable enthalpy changes predominantly determine their adsorption behavior; however, larger amino acids tend to bind more tightly with the graphene surface, which is thermodynamically dominated by the entropy variations despite the structural characteristics of amino acids. This study reveals specific entropy-enthalpy mechanisms underlying amino acid adsorption at the solid-liquid interface, providing guidance for surface design and synthesis of new biomolecules.

Physical insight into the entropy-driven ion association

The ion association is widely believed to be dominated by the favorable entropy change arising from the release of water molecules from ion hydration shells. However, no direct thermodynamic evidence exists to validate the reliability and suitability of this view. Herein, we employ complicated free energy calculations to rigorously split the free energy including its entropic and enthalpic components into the water-induced contributions and ion-ion interaction terms for several ion pairs from monatomic to polyatomic ions, spanning the size range from small kosmotropes to large chaotropes (Na⁺ , Cs⁺ , Ca²⁺ , F^- , I^- , CO₃ ^2- , and HPO₄ ^2- ). Our results successfully reveal that though ion associations are indeed determined by a delicate balance between the favorable entropy variation and the repulsive enthalpy change, the entropy gain dominated by the solvent occurs only for the monatomic ion pairing. The water-induced entropic contribution significantly goes against the ion pairing between polyatomic anion and cation, which is, alternatively, dominated by the favorable entropy from the ion-ion interaction term, due to the configurational arrangement of polyatomic anions involved in ion association. The structural and dynamic analysis demonstrates that the entropy penalty from the water phase is primarily ascribed to the enhanced stability of water molecules around the cation imposed by the incoming anion. Our study successfully provides a fundamental understanding of water-mediated ion associations and highlights disparate lengthscale dependencies of the dehydration thermodynamics on the specific types of ions.

Adsorption characteristics of peptides on ω-functionalized self-assembled monolayers: a molecular dynamics study

Molecular dynamics simulations were employed to investigate the adsorption behavior of a variety of amino-acid side-chain analogs (SCAs) and a β-hairpin (HP7) peptide on a series of liquid-like self-assembled monolayers (SAMs) with terminal functional groups of -OH, -OCH₃, -CH₃, and -CF₃. The relationships between the adsorption free energy of the SCAs and the interfacial properties of water on the SAMs were examined to determine the acute predictors of protein adsorption on the SAM surfaces. The structural changes of HP7 on the SAM surfaces were also investigated to understand the relationship between the surface nature and protein denaturation. It was found that the adsorption free energy of the SCAs was linearly related to the surface hydrophobicity, which was computed as the free energy of cavity formation near the SAM-water interfaces. In addition, the hydrophobic -CH₃ and -CF₃ SAMs produced substantial conformational changes in HP7 because of the strong hydrophobic attractions to the nonpolar side chains. The hydrophilic surface terminated by -OH also promoted structural changes in HP7 resulting from the formation of hydrogen bonds between the hydrophilic tail and HP7. Consequently, the moderate amphiphilic surface terminated by -OCH₃ avoided the denaturation of HP7 most efficiently, thus improving the biocompatibility of the surface. In conclusion, these results provide a deep understanding of protein adsorption for a wide range of polymeric surfaces, and they can potentially aid the design of appropriate biocompatible coatings for medical applications.

Insights into the conformation changes of SARS-CoV-2 spike receptor-binding domain on graphene

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been widely spread in the world, causing more than two million deaths and seriously threatening human life. Effective protection measures are important to prevent the infection and spreading of the virus. To explore the effects of graphene on the virus adsorption and its biological properties, the adsorption process of the receptor binding domain (RBD) of SARS-CoV-2 on graphene has been investigated by molecular dynamics simulations in this paper. The results show that RBD can be quickly adsorbed onto the surface of graphene due to π - π stacking and hydrophobic interactions. Residue PHE486 with benzene ring has stronger adsorption force and the maximum contact area with graphene. Graphene significantly affects the secondary structure of RBD area, especially on the three key sites of binding with human ACE2, GLY476, PHE486 and ASN487. The binding free energy of RBD and graphene shows that the adsorption is irreversible. Undoubtedly, these changes will inevitably affect the pathogenicity of the virus. Therefore, this study provides a theoretical basis for the application of graphene in the protection of SARS-CoV-2, and also provides a reference for the potential application of graphene in the biomedical field.

Spontaneous desorption of protein from self-assembled monolayer (SAM)-coated gold nanoparticles induced by high temperature

The nonspecific binding of proteins with nanomaterials (NMs) is a dynamic reversible process including both protein adsorption and desorption parts, which is crucial for controlled release of protein drug loaded by nanocarriers. The nonspecific binding of proteins is susceptible to high temperature, whereas its underlying mechanism still remains elusive. Here, the binding behavior of human serum albumin (HSA) with an amino-terminated self-assembled monolayer (SAM)-coated gold (111) surface was investigated by using molecular dynamics (MD) simulations. HSA binds to the SAM surface through salt bridges at 300 K. As the temperature increases to 350 K, HSA maintains its native structure, while the salt bridges largely diminish owing to the considerable lateral diffusion of HSA on the SAM. Moreover, the interfacial water located between HSA and the SAM gets increased and prevents the reformation of the salt bridges of HSA with the SAM, which reduces the binding affinity of HSA. And HSA eventually desorbs from the SAM. The depiction of thermally induced spontaneous protein desorption enriches our understanding of reversible binding behavior of protein with NMs, and may provide new insights into the controlled release of protein drugs delivered by using nanocarriers under the regulation of high temperature.

Lessons from a Challenging System: Accurate Adsorption Free Energies at the Amino Acid/ZnO Interface

We undertake steps to overcome four challenges that have hindered the understanding of ZnO/biomolecule interfaces at the atomic scale: parametrization of a classical force field, ZnO surface termination and amino acid protonation state in methanol, and convergence of enhanced sampling molecular dynamics simulations. We predict adsorption free energies for histidine, serine, cysteine, and tryptophan in remarkable agreement with experimental measurements obtained via a novel indicator-displacement assay. Adsorption is driven by direct surface/amino-acid interactions mediated by terminal hydroxyl groups and stabilized by strongly structured methanol solvation shells.

Simulation Insight into the Synergic Role of Citrate and Polyaspartic Peptide in Biomineralization

Hydroxyapatite (HAP) is one of the most important inorganic components in biological minerals such as bones and teeth. More than 90% of the total citrate is accumulated in human bones and other biomineralized tissues. In addition, mineralizing proteins are enriched in glutamate and aspartate residues, which are important for their mineral-regulating properties. However, how citrate ions (CITs) and/or acidic amino acids regulate the formation of HAP is still unclear. In this work, molecular dynamics simulations were performed to study how CIT regulates the adsorption behavior of polyaspartic acid (PASP) on the HAP surface in the calcium phosphate solution. The simulation results indicate that PASP can be used as an ion chelator to complex Ca²⁺ and can serve as templates for HAP mineralization by templating the distribution of Ca²⁺ on its surface, which are attributed to the -COO^- and α-helix structure. Most importantly, the orientation distributions of PASP in all systems are narrower with the help of CIT, thereby PASP can be adsorbed on the HAP surface stably with a "lying-down" orientation. This indicates that CIT can be used as a bridging agent to bond the acidic peptide to the HAP surface in biomineralization. Thus, the synergic role of CIT and the acidic peptide on the HAP surface were revealed in this work, which can provide new insights into the interfacial phenomena during the biomineralization.

Advancing Rational Control of Peptide-Surface Complexes

Understanding peptide-surface interactions is crucial for programming self-assembly of peptides at surfaces and in realizing their applications, such as biosensors and biomimetic materials. In this study, we developed insights into the dependence of a residue's interaction with a surface on its neighboring residue in a tripeptide using molecular dynamics simulations. This knowledge is integral for designing rational mutations to control peptide-surface complexes. Using graphene as our model surface, we estimated the free energy of adsorption (ΔA_ads) and extracted predominant conformations of 26 tripeptides with the motif LNR-CR-Gly, where LNR and CR are variable left-neighboring and central residues, respectively. We considered a combination of strongly adsorbing (Phe, Trp, and Arg) and weakly adsorbing (Ala, Val, Leu, Ser, and Thr) amino acids on graphene identified in a prior study to form the tripeptides. Our results indicate that ΔA_ads of a tripeptide cannot be estimated as the sum of ΔA_ads of each residue indicating that the residues in a tripeptide do not behave as independent entities. We observed that the contributions from the strongly adsorbing amino acids were dominant, which suggests that such residues could be used for strengthening peptide-graphene interactions irrespective of their neighboring residues. In contrast, the adsorption of weakly adsorbing central residues is dependent on their neighboring residues. Our structural analysis revealed that the dihedral angles of LNR are more correlated with that of CR in the adsorbed state than in bulk state. Together with ΔA_ads trends, this implies that different backbone structures of a given CR can be accessed for a similar ΔA_ads by varying the LNR. Therefore, incorporation of context effects in designing mutations can lead to desired peptide structure at surfaces. Our results also emphasize that these cooperative effects in ΔA_ads and structure are not easily predicted a priori. The collective results have applications in guiding rational mutagenesis techniques to control orientation of peptides at surfaces and in developing peptide structure prediction algorithms in adsorbed state from its sequence.

Entropy-Enthalpy Compensation in Peptide Adsorption on Solid Surfaces: Dependence on Surface Hydration

Although protein adsorption at the solid-water interface is of immense importance, understanding the crucial role of the water phase in mediating protein-surface interactions is lacking, particularly due to the lack of fundamental thermodynamic data. Herein, we have performed complicated free energy calculations and successfully extracted the entropy and enthalpy changes of molecular adsorption on solids. Using the gold and graphene as the surface models with distinct affinities to the water phase, we successfully unravel the sharply opposite manners of entropy-enthalpy compensation in driving water and tripeptide adsorptions on two surfaces. Though the thermodynamic features of water adsorption on surface are enthalpically dominated based on the positions of free energy barriers and minima, the favorable entropy term significantly decreases the free energy barrier and further stabilizes the adsorbate at the adsorption site on the graphene surface. For the peptide, the shape of the adsorption free energy profile is jointly determined by the enthalpy and entropy changes, which, however, alternatively act the driving force to promote the peptide adsorption on the Au surface and graphene surface. The distinct structural and dynamic properties of solid-liquid interfaces account for the special role of the interfacial water phase in regulating the competitive relationship between the entropy and enthalpy variations.

Structure-Activity Relationships of Hydroxyapatite-Binding Peptides

Elucidating the structure-activity relationships between biomolecules and hydroxyapatite (HAP) is essential to understand bone mineralization mechanisms, develop HAP-based implants, and design drug delivery vectors. Here, four peptides identified by phage display were selected as model HAP-binding peptides (HBPs) to examine the effects of primary amino acid sequence, phosphorylation of serine, presence of charged amino acid residues, and net charge of the peptide on (1) HAP-binding affinity, (2) secondary conformation, and (3) HAP nucleation and crystal growth. Binding affinities were determined by obtaining adsorption isotherms by mass depletion, and the conformations of the peptides in solution and bound states were observed by circular dichroism. Results showed that the magnitude of the net charge primarily controlled binding affinity, with little dependence on the other HBP features. The binding affinity and conformation results were in good agreement with our previous molecular dynamics simulation results, thus providing an excellent benchmark for the simulations. Transmission electron microscopy was used to explore the effect of these HBPs on calcium phosphate (Ca-PO₄) nucleation and growth. Results indicated that HBPs may inhibit nucleation of Ca-PO₄ nanoparticles and their phase transition to crystalline HAP, as well as control crystal growth rates in specific crystallographic directions, thus changing the classical needle-like morphology of inorganically grown HAP crystals to a biomimetic plate-like morphology.

Bioremediation of gaseous methyl tert-butyl ether by combination of sulfuric acid modified bagasse activated carbon-bone biochar beads and Acinetobacter indicus screened from petroleum contaminated soil

Combination of sulfuric acid modified bagasse activated carbon-bone biochar beads and Acinetobacter indicus screened from petroleum contaminated soil was the best condition for gaseous methyl tert-butyl ether (MTBE) removal. It was found that H₂SO₄ modified bagasse AC in powder form had higher adsorption capacity (989.33 mg g^-1) than that in bead form (1.94 mg g^-1). In addition, bone biochar in powder form (3.51 mg g^-1) also had higher adsorption capacity than that in bead form (1.63 mg g^-1). This was the fact that material beads contained high moisture content that inhibited the penetration of gaseous MTBE into the material. And a mixed material of H₂SO₄ modified bagasse AC-bone biochar beads had the highest adsorption capacity (2.22 mg g^-1) compared to individual H₂SO₄ modified bagasse AC beads (1.94 mg g^-1) and bone biochar beads (1.63 mg g^-1) due to a mixed material had more rough surface and high surface area on its material. So, gaseous MTBE can penetrate through this material more easily. Although the maximum adsorption capacity of H₂SO₄ modified bagasse AC in powder form was the highest but microorganism cannot sustain and survive in this form for a long time. Therefore, the material beads were more suitable for microorganism to grow and degrade gaseous MTBE. Microorganism can degrade MTBE and caused no secondary wastes. Moreover, A. indicus was a novel strain for MTBE removal that has not been previously reported. Therefore, a combination of A. indicus-mixed material beads was a good choice for MTBE removal in a biofilter system.

Understanding the Thermodynamics of the Binding of PAMAM Dendrimers to Graphene: A Combined Analytical and Simulation Study

We investigate the thermodynamics of the binding of a poly(amidoamine) dendrimer to an uncharged graphene sheet as a function of the pH level using umbrella sampling simulations and a mean-field theory for generations three and four. We find that the dendrimer strongly binds to the graphene sheet ( O (100) kcal/mol) from our potential of mean force (PMF) calculations. In specific, we find that the dendrimer binds the most at neutral pH (∼7) and the least at low pH (∼4). We explain this nonmonotonic nature of the dendrimer's adsorption by studying the interactions contributing to the PMF, i.e., the dendrimer-graphene, dendrimer-water, and dendrimer-ion interactions. We also corroborate our PMF calculations with molecular mechanics generalized Born surface area analysis and free energies obtained from a mean-field theory of Flory-Huggins-Debye-Hückel type [ Muthukumar , M. , J. Chem. Phys. 2010 , 132 , 084901 ], including electrostatic interactions. We find that the van der Waals interactions between the dendrimer and the graphene alone cannot capture the accurate trends in the binding free energies (BEs) as a function of pH. The solvent and the counterions present in the system are also found to have a major influence on these trends. We demonstrate that the dendrimer-graphene and dendrimer-water interactions become favorable, whereas the dendrimer-ion interaction becomes unfavorable, as the dendrimer binds to graphene. These opposing effects lead to the observed nonmonotonicity in the BE trends. Our theoretical model also reproduces these trends in the subinteractions contributing to the PMF. To the best of our knowledge, this is a novel attempt where an equivalence between theory and simulations is made in the context of the dendrimer's adsorption.

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.

The assembly of silk fibroin and graphene-based nanomaterials with enhanced mechanical/conductive properties and their biomedical applications

The assembly of silk fibroin and graphene-based nanomaterials would present fantastic properties and functions via optimizing the interaction between each other, and can be processed into various formats to tailor specific biomedical applications.

Physical origin underlying the prenucleation-cluster-mediated nonclassical nucleation pathways for calcium phosphate

The involvement of prenucleation clusters (PNCs) in crystallization from a supersaturated solution has been recently admitted within the framework of nonclassical nucleation theory; however, little is known about PNCs, at the quantitative level, for their formation mechanism and stability, the new phase formed by them, as well as their impact on nucleation barriers. Herein, using the sophisticated free energy calculations with a cumulative simulation time of over 5 μs, we identify a thermodynamically favored pathway of the PNC-mediated nucleation for calcium phosphate, starting with the ion pair association in solution. We demonstrate that such an ion association occurs not only between cations and anions, but also for the polyatomic species with charges of the same sign, which, in fact, leads to PNC formation via the consecutive coordination of the phosphate ions to calcium. The free energy decomposition calculations illustrate that the water phase is capable to either hinder or promote ion association for the abovementioned processes, and its specific role is intricately related to the characteristics of the hydration shell around calcium ions. The favorable interactions between the highly charged species play a crucial role in stabilizing the PNC complexes and the aggregates formed by PNCs. Furthermore, our present work reveals that the uptake of an extra calcium ion is the first and mandatory step to trigger PNC aggregation into amorphous calcium phosphate (ACP) by eliminating the related free energy barriers. Our theoretical study successfully provides quantitative explanations to a large set of experimental data in the field, which is currently under intense discussions associated with the nonclassical nucleation mechanism. The combination of computational methods developed in our present work offers a feasible and general solution to quantitatively and systematically study ion associations and crystal nucleation/growth in an aqueous solution at the atomic level, which are normally inaccessible to most of the existing experimental acquisitions.